MATLAB,Python,Scilab,Julia比較 第5章 その41【フーリエ係数⑤】

MATLAB,Python,Scilab,Julia比較 第5章 その41【フーリエ係数⑤】 数値計算
MATLAB,Python,Scilab,Julia比較 第5章 その41【フーリエ係数⑤】

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

はじめに

フーリエ係数に至る道。
フーリエ係数を求める式の一般化について。
今回はsin成分とa0

登場人物

博識フクロウのフクさん

指差しフクロウ

イラスト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

【再掲】フーリエ係数に至る道

太郎くん
太郎くん

まずは、フーリエ係数に至る道を再掲。

  • 偶関数
  • 奇関数
  • 関数の内積
  • 三角関数の加法定理
  • 三角関数の積和公式
  • 重要な極限値
  • 三角関数の直交性
  • フーリエ係数
フクさん
フクさん

今回はフーリエ係数の話の続き。
フーリエ係数を求める式の一般化について。
今回はsin成分とa0

フーリエ係数を求める式の一般化(sin関数)

太郎くん
太郎くん

前回は、フーリエ係数の\(a_n\)を求める式の一般化だったね。

フクさん
フクさん

今回は、\(b_n\)の方になるな。
やり方は\(a_n\)の時とほぼ一緒だ。

\(
\begin{eqnarray}
\displaystyle f(x)\cdot\sin(nx)&=&\frac{a_0}{2}+\sum_{n=0}^\infty (a_n\cos(nx)+b_n\sin(nx))\cdot\sin(nx)\\
\displaystyle &=&\frac{a_0}{2}\cdot\sin(nx)+\dots+{\color{red}b_n\sin(nx)\cdot\sin(nx)}+\dots\\
&=&b_n\sin(nx)\cdot\sin(nx)\\
&=&b_n\pi\\
\displaystyle\therefore b_n&=&\frac{f(x)\cdot\sin(nx)}{\pi}=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x)\sin(nx)dx
\end{eqnarray}
\)

太郎くん
太郎くん

そうか。
\(b_n\)の方は、cos関数じゃなくてsin関数を使えば良いのか。
最終的な式もcosがsinになった程度でほぼ一緒だね。

a0を求める

フクさん
フクさん

ついでに\(a_0\)を求める式も導出しよう。

太郎くん
太郎くん

あれ?
\(a_n\)の方はやったから\(a_0\)も求められるんじゃない?

フクさん
フクさん

まぁ、それでも求められなくはないのだが、ここでは\(a_0\)を特別扱いしておいた方が良い。
以下の感じで求める。

\(
\begin{eqnarray}
\displaystyle f(x)\cdot 1&=&\frac{a_0}{2}+\sum_{n=0}^\infty (a_n\cos(nx)+b_n\sin(nx))\cdot 1\\
\displaystyle &=&{\color{red}\frac{a_0}{2}\cdot 1}+a_1\cdot\cos(x)\cdot 1+\dots\\
\displaystyle &=&\frac{a_0}{2}\cdot 1=\int_{-\pi}^{\pi}\frac{a_0}{2}dx=\frac{a_0}{2}2\pi=a_0\pi\\
\displaystyle\therefore a_0&=&\frac{f(x)\cdot 1}{\pi}=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x)dx
\end{eqnarray}
\)

太郎くん
太郎くん

1で畳み込みをするってどういう状況だ???

フクさん
フクさん

常に1の関数との畳み込みだな。
まぁ、定積分に於いては矩形波をイメージすると認識しやすいかな。

太郎くん
太郎くん

それで、\(-\pi\sim\pi\)で定積分すると\(2\pi\)になるのか。

フクさん
フクさん

そうそう。

まとめ

フクさん
フクさん

まとめだよ。

  • フーリエ係数のbnを求める式の一般化。
  • ついでにa0を求める式も一般化。
  • 常に1のような定数関数は畳み込み積分に於いては矩形波をイメージすると認識しやすい。

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

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

https://amzn.to/4azyPDj

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

https://amzn.to/4bDH3vu

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

https://amzn.to/4bsygfU

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

https://amzn.to/3V83fIl

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

https://amzn.to/3ysbfvf

コメント

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