Excel VBA 数学教室ではアフィリエイトプログラムを利用して商品を紹介しています。

フーリエ変換(非周期関数の展開式)

フーリエ変換

周期 2L をもつ関数 f(x) は離散係数 cn を用いて
 f(x)=n=cnexp(inπxL)cn=12LLLf(t)exp(inπtL)
フーリエ級数展開できます。ここで L の極限をとれば、非周期関数の展開式も得られるはずです。そこで係数 cnf(x) の展開式に代入し、
 ωn=nπL,Δω=ωnωn1=πL
とおくと、
 f(x)=12πn=[LLf(t)eiωntdt]eiωnxΔω
という表式が得られます。L とすれば総和の記号が積分に変わって
 f(x)=12πdωf(t)eiω(tx)dt
となります。これをフーリエ積分公式とよびます。そして
 F(ω)=f(t)eiωtdt
とおけば、
 f(x)=12πF(ω)eiωxdω
と表すことができます。変数 t をあらためて x に置き直すと、

(A)F(ω)=f(x)eiωxdx(B)f(x)=12πF(ω)eiωxdω

となります。つまり f(x) から F(ω) を求めることができて、またその逆に F(ω) から f(x) を求めることができます。そこで F(ω)f(x)フーリエ変換f(x)F(ω)フーリエ逆変換とよびます。変換の対称性を考慮して係数を 1/2π とした

(C)F(ω)=12πf(x)eiωxdx(D)f(x)=12πF(ω)eiωxdω

という表式もよく用いられます。より簡略的に
 F(ω)=F[f(x)]f(x)=F1[F(ω)]
と表すこともあります。

単一矩形波のフーリエ変換

例として x∣≤d の範囲にだけ存在する波(単一矩形波)
 f(x)={12d(x∣≤d)0(x∣>d)
をフーリエ変換してみます。
 
Excel単一矩形波のフーリエ変換

定義にしたがって F(ω) を計算してみると
 F(ω)=f(x)eiωxdx=sinωdωd
となります。ω の変化に対する F(ω) を図示すると次のようになります。
 
Excel連続スペクトル

原点でピーク 1 をとり左右に減衰しています。これは f(x)ω=0 , すなわち y=1 という波(というより直線)を成分として最も多く含んでおり、ω の絶対値が大きくなるほど eiωx の寄与が小さくなっていくことを示しています。このように、F(ω) はフーリエ級数の係数 cn と本質的に全く同じものであり、関数 f(x) の中に、ある角振動数 ω の波 eiωx をどれだけ含んでいるかを示すものです。そうした意味で F(ω) のことを連続スペクトルとよぶことがあります。 

フーリエ変換の基本性質

以下にフーリエ変換の基本的性質を載せておきます。いずれも定義にしたがって簡単に導けるので証明は省略します。

f(x)g(x) のフーリエ変換をそれぞれ
F[f(x)]=F(ω),F[g(x)]=G(ω)とするとき、以下の式が成り立つ。

① 共役関数の変換
(E)F[f(x)]=F(ω)

② 重ね合わせの原理
(F)F[af(x)+bg(x)]=aF(ω)+bF(ω)

③ 尺度変換
(G)F[f(ax)]=1aF(ωa)(a0)

④ 対称性
(H)F[F(ω)]=2πf(x)

⑤ 周波数シフト
(I)F[f(t)eiωt]=F(ωω0)

⑥ 時間軸移動
(J)F[f(tt0)]=F(ω)eiωt0

工学・物理学等への応用で分かりやすいように、⑤ と ⑥ では変数 x を時間を表す変数 t に書き換えてあります。

フーリエ変換の微分

f(x) のフーリエ変換
 f(x)=12πF(ω)eiωxdω
x について微分すると
 f(x)=12πF(ω)(iω)eiωxdω
となります。これは f(x) のフーリエ変換が
 F[f(x)]=iωF(ω)
となることを表しています。n 回微分すると
 (K)F[f(n)(x)]=(iω)nF(ω)
となります。

エクセルや数学に関するコメントをお寄せください