MATLAB,Python,Scilab,Julia比較 第5章 その59【マクローリン展開⑤】

MATLAB,Python,Scilab,Julia比較 第5章 その59【マクローリン展開⑤】 数値計算
MATLAB,Python,Scilab,Julia比較 第5章 その59【マクローリン展開⑤】

バックナンバーはこちら。
https://www.simulationroom999.com/blog/compare-matlabpythonscilabjulia5-backnumber/

はじめに

前回は指数関数のマクローリン展開の説明。
今回はcos関数のマクローリン展開について。

登場人物

博識フクロウのフクさん

指差しフクロウ

イラストACにて公開の「kino_k」さんのイラストを使用しています。
https://www.ac-illust.com/main/profile.php?id=iKciwKA9&area=1

エンジニア歴8年の太郎くん

技術者太郎

イラストACにて公開の「しのみ」さんのイラストを使用しています。
https://www.ac-illust.com/main/profile.php?id=uCKphAW2&area=1

【再掲】複素フーリエ級数に至る道

太郎くん
太郎くん

まずは複素フーリエ級数に至る道を再掲。

  • テイラー級数
  • マクローリン級数
  • 指数関数のマクローリン展開
  • cos(x)のマクローリン展開
  • sin(x)のマクローリン展開
  • オイラーの公式
  • 複素フーリエ級数
フクさん
フクさん

今回は、cos(x)のマクローリン展開。

cos(x)のマクローリン展開

太郎くん
太郎くん

前回は指数関数だったけど、今回はcos関数か。
感覚的には指数関数より楽そうなイメージがある。

フクさん
フクさん

指数関数よりもちょっと面倒かなー。
まぁ、手順は一緒なんだけど。

太郎くん
太郎くん

え゛

フクさん
フクさん

指数関数の時と同じく、
\(\cos(x)\)を微分しまくってみよう。

\(
\begin{eqnarray}
f(x)&=&\cos(x)\\
f^{\prime}(x)&=&-\sin(x)\\
f^{\prime\prime}(x)&=&-\cos(x)\\
f^{\prime\prime\prime}(x)&=&\sin(x)\\
f^{\prime\prime\prime\prime}(x)&=&\cos(x)\dots4階微分で\cos(x)に戻る\\
\end{eqnarray}
\)

太郎くん
太郎くん

何か、指数関数の時よりややこしいことになってる・・。

フクさん
フクさん

ただ、4階微分で戻るという特性は重要だ。
一旦整理するとこうなる。

\(
f(x)^n=\cases{
\cos(x)\dots(n=0,4,8,\dots)\\
-\sin(x)\dots(n=1,5,9,\dots)\\
-\cos(x)\dots(n=2,6,10,\dots)\\
\sin(x)\dots(n=3,7,11,\dots)\\
}
\)

太郎くん
太郎くん

一応、パターンみないなのはあるけど、
これを前回みたいな一つの式にまとめるってのはむつかしそう・・・。

フクさん
フクさん

そうでもない。
さらに、原点の時に限定すると、

\(
f(0)^n=\cases{
\cos(0)=1\dots(n=0,4,8,\dots)\\
-\sin(0)=0\dots(n=1,5,9,\dots)\\
-\cos(0)-1\dots(n=2,6,10,\dots)\\
\sin(0)=0\dots(n=3,7,11,\dots)\\
}
\)

太郎くん
太郎くん

1,0,-1しか出てこないのか。

フクさん
フクさん

nが奇数の時は0になり、偶数の時は符号反転ってことだな。
式として整理すると以下になる。

\(
\begin{eqnarray}
\displaystyle f(x)&=&f(0)+\frac{f^\prime(0)}{1!}x+\frac{f^\prime\prime(0)}{2!}x^2+\dots\\
\displaystyle &=&f(0)+\sum_{n=1}^\infty\frac{f^n(0)}{n!}x^n\\
\displaystyle \cos(x)&=&1-x+\frac{x^2}{2!}-\frac{x^4}{4!}+\frac{x^6}{6!}+\dots\\
\displaystyle &=& \sum_{n=0}^\infty\frac{(-1)^n x^{2n}}{(2n)!}
\end{eqnarray}
\)

太郎くん
太郎くん

0の時は消しちゃって良いから、1と-1の時だけ感がればシンプルになるのか。

フクさん
フクさん

そうそう。

cos関数のマクローリン展開の式を元にプロット

フクさん
フクさん

nを徐々に増やした場合のプロットも見せておこう。

cos関数をマクローリン展開
太郎くん
太郎くん

徐々に一致する範囲が広がっていく感じか。
無限にやれば、cos関数と同一にはなりそう。

まとめ

フクさん
フクさん

まとめだよ。

  • cos関数をマクローリン展開。
  • とりあえず微分しまくると4階微分の周期が見える。
  • これを元にマクローリン展開。
  • cos関数をマクローリン展開したプロットも出してみた。

バックナンバーはこちら。

マンガでわかるフーリエ解析

Amazon.co.jp

手を動かしてまなぶ フーリエ解析・ラプラス変換

Amazon.co.jp

物理数学 量子力学のためのフーリエ解析・特殊関数

物理数学 量子力学のためのフーリエ解析・特殊関数 | 柴田 尚和, 是常 隆 | 数学 | Kindleストア | Amazon
Amazonで柴田 尚和, 是常 隆の物理数学 量子力学のためのフーリエ解析・特殊関数。アマゾンならポイント還元本が多数。一度購入いただいた電子書籍は、KindleおよびFire端末、スマートフォンやタブレットなど、様々な端末でもお楽しみい...

単位が取れるフーリエ解析ノート

https://amzn.to/3V83fIl

今日から使えるフーリエ変換 普及版 式の意味を理解し、使いこなす

https://amzn.to/3ysbfvf

コメント

タイトルとURLをコピーしました