【入門】フーリエの積分公式①【数値計算】

【入門】フーリエの積分公式①【数値計算】 数値計算
【入門】フーリエの積分公式①【数値計算】

角周波数へ

ここで、\(n\)を角周波数\(\omega\)に変更する。

\(n\)の段階では、周期\(2L\)に於ける波の数って感じだったのだが、
これを1周期を\(2\pi\)中の波の数にする感じ。
変換の式は簡単で、以下になる。

\(
\displaystyle \omega_n=n\frac{2\pi}{2L}=\frac{n\pi}{L}
\)

\(2L\)だったものを\(2\pi\)として読み替えるわけだから、こういう式になる。

これにより\(\displaystyle\frac{n\pi}{L}\)を\(\omega_n\)に置き換えられる。

\(
\displaystyle f(x)=\sum_{n=-\infty}^\infty \Big\{\frac{1}{2L}\int_{-L}^Lf(t)^{-i\omega_n t} \,dt\Big\}e^{i\omega_n x}
\)

ちょっぴりシンプルになった。

ωの刻みをΔωにする

ここで、角周波数\(\omega\)の刻みである\(\Delta\omega\)について考える。
端的に言うと、\(\omega_n-\omega_{n-1}\)を定義するってことにある。

定義方法は普通に導出すればOKだ。

\(
\begin{eqnarray}
\displaystyle \Delta\omega&=&\omega_n-\omega_{n-1}=\frac{n\pi}{L}-\frac{(n-1)\pi}{L}=\frac{\pi}{L}\\
L&=&\frac{\pi}{\Delta\omega}
\end{eqnarray}
\)

本当にそのまんま。

これを代入すると以下になる。

\(
\displaystyle f(x)=\sum_{n=-\infty}^\infty \Big\{\frac{1}{2\pi}\int_{-L}^Lf(t)^{-i\omega_n t} \,dt\Big\}e^{i\omega_n x}\Delta\omega
\)

どんどん式が変化しては行くけど、どこに向かってるかはわからないかもしれない。
方向性としては、
なんとか離散的な係数算出を連続的な周波数算出に持っていきたい
ってところになる。

まとめ

  • フーリエに積分公式は複素フーリエ級数と複素フーリエ係数から導出する。
  • 周期2Lの波の数を示すnを周期2πに於ける波の数である角周波数ωに変換。
  • 角周波数ωの刻みであるΔωについて説明。
  • Δωを定義することで、離散的な係数算出が連続的な角周波数算出に近づけていっている。

MATLAB、Python、Scilab、Julia比較ページはこちら

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

Amazon.co.jp

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

手を動かしてまなぶ フーリエ解析・ラプラス変換 | 山根 英司 | 数学 | Kindleストア | Amazon
Amazonで山根 英司の手を動かしてまなぶ フーリエ解析・ラプラス変換。アマゾンならポイント還元本が多数。一度購入いただいた電子書籍は、KindleおよびFire端末、スマートフォンやタブレットなど、様々な端末でもお楽しみいただけます。

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

Amazon.co.jp

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

Amazon.co.jp

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

Amazon.co.jp

コメント

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