【計算力学熱流体2級でよく出る物理数学】テーラー展開

こんにちは（@t_kun_kamakiri)

本日は「計算力学技術者試験の熱流体2級」で出てくるテーラー展開の解説をします。
テーラー展開は数値計算において微分を離散化する際には、必ず必要になってくる知識です。
理系学部の初学年の解析学で習うですが、社会人になって必ずしも皆が習っているわけではありませんので、これから計算力学技術者の熱流体を勉強する方を対象に問題を通して理解をしていければと思います。

こんな方が対象

計算力学技術者試験の熱流体2級を勉強している方
テーラー展開は大学で習ったけど忘れてしまったのでサクッと復習したい方

※微分や偏微分などは理解している方を対象にしています。

では、問題です。

問題1

$\sin x$を$x=0$のまわりでマクローリン展開しなさい。
第3項まで残しなさい。

問題2

$\cos x$を$x=\pi$のまわりでテーラー展開しなさい。
第2項まで残しなさい。

問題1の解説は本記事に書いていますが、問題2はご自身で解いてください。
解答はまとめに書いています。

テーラー展開
1. テーラー展開の次数
マクローリン展開
問題1の解説
まとめ
計算力学技術者のための問題アプリ
お勧めの参考書

テーラー展開

$x=a$のまわりでの関数の形を知りたい場合、テーラー展開を使います。

テーラー展開は以下の自身の関数の微分を使って、1次曲線、2次曲線、3次曲線・・・など多項式で表現したものです。

\begin{align*}f(x)=\sum_{n=0}^{\infty}\frac{f^{(n)}(x)}{n!}|_{x=a}(x-a)^{n}\tag{1}\end{align*}

と書くことができます。

Σの和があるとわかりにくいですが、展開しておくと多項式近似であることがわかります。

\begin{align*}f(x)&=f(a)+\frac{f^{(1)}(x)}{1!}|_{x=a}(x-a)\\
&+\frac{f^{(2)}(x)}{2!}|_{x=a}(x-a)^{2}\\
&+\cdot\cdot\cdot+\frac{f^{(n)}(x)}{n!}_{x=a}(x-a)^{n}+\\
&\cdot\cdot\cdot\end{align*}

※$f^{(n)}(x)=\frac{d^{n}f(x)}{dx^{n}}$：n階微分

念のため注意ですが$\frac{f^{(1)}(x)}{1!}|_{x=a}$は$x$で微分してから$x=a$を代入するって意味です。

$x=a$を代入してから微分してはいけませんよ（0にしかなりませんから）。

テーラー展開の次数

テーラー展開した際に何次の項まで残すかによって、数値的な離散化の精度を表すことがあります。

マクローリン展開

任意の$x=a$近傍での展開はテーラー展開と呼びますが、$x=0$近傍での展開を特にマクローリン展開と呼びます。

例えば「3次の項までマクローリン展開しなさい」と問われたら、以下のようにします。

\begin{align*}
f(x)=f(0)+\frac{1}{1!}\frac{df}{dx}|_{x=0}x+\frac{1}{2!}\frac{d^2f}{dx^2}|_{x=0}x^{2}+\frac{1}{3!}\frac{d^3f}{dx^3}|_{x=0}x^{3}
\end{align*}

マクローリン展開しなさいということなので、$x=0$の近傍で展開するということです。
わざわざ「$x=0$近傍でマクローリン展開をしなさい」と書かれないので注意が必要です。

問題1の解説

問題1

$\sin x$を$x=0$のまわりでマクローリン展開しなさい。
第3項まで残しなさい。

「$x=0$のまわりで」と書いているので、マクローリン展開自体が$x=0$のまわりで展開することなので、「$x=0$のまわりで」と表記されていない場合もあります。

「3次の項までマクローリン展開」は、以下のようにします。

\begin{align*}
f(x)=f(0)+\frac{1}{1!}\frac{df}{dx}|_{x=0}x+\frac{1}{2!}\frac{d^2f}{dx^2}|_{x=0}x^{2}+\frac{1}{3!}\frac{d^3f}{dx^3}|_{x=0}x^{3}\tag{1}
\end{align*}

ここで、$f(x)=\sin x$と考えれば良いです。

テーラー展開やマクローリン展開は、大学の解析学など1年生の段階で習うものですが、初めて見る人にとっては少々インパクトのある難しいイメージを持ちますが、(1)を見てわかる通りひたすら微分すればいいだけです。

では、微分してみます。