Games101 の第 6 講ではフーリエ変換に関する内容が含まれており、授業後に少し時間をかけてこのことを理解しました(大学で学んだことはありますが、ほとんど忘れてしまいました)。そこで、学びの心得を記録するためにこの記事を書きました。
参考資料として、私は DR_CAN 先生の『純干貨数学推導』シリーズを強くお勧めします!この記事の大部分は彼の講義ノートに基づいています。
この記事では、以下のいくつかの側面から順に推導します。
三角関数の直交性
周期が 2 π \pi π の関数のフーリエ級数展開
周期が T の関数のフーリエ級数展開
フーリエ級数の複素形式
非周期関数のフーリエ級数展開
注意が必要なのは、私の知識が不足しているため、フーリエ級数展開の公式の証明を提供することはできません。G=mg がなぜ成り立つのかを証明できないのと同じように、真理をそのまま使う方が簡単かもしれません。しかし、理解できないことが気になる読者には、王繼良先生が書いた記事 をお勧めします。この文章では、フーリエが熱伝導方程式からフーリエ級数をどのように抽出するかの過程が詳しく説明されています。
フーリエ級数展開 : f ( x ) = a 0 2 + Σ n = 1 + ∞ ( a n cos n x + b n sin n x ) f(x) = \frac{a_0}{2} + \varSigma_{n=1}^{+\infty}(a_{n}\cos nx + b_{n}\sin nx) f ( x ) = 2 a 0 + Σ n = 1 + ∞ ( a n cos n x + b n sin n x )
三角関数の直交性#
もし関数 f (x) と g (x) が存在し、両者が ∫ a b f ( x ) g ( x ) d x = 0 \int_{a}^{b} f(x)g(x) dx = 0 ∫ a b f ( x ) g ( x ) d x = 0 を満たすなら、関数 f (x) と g (x) は区間 [a, b] で直交すると呼びます。
三角関数 0 ( sin 0 ) , 1 ( cos 0 ) , sin x , cos x , sin 2 x , cos 2 x . . . . . sin n x , cos n x 0(\sin 0), 1(\cos 0), \sin x, \cos x, \sin 2x, \cos 2x..... \sin nx, \cos nx 0 ( sin 0 ) , 1 ( cos 0 ) , sin x , cos x , sin 2 x , cos 2 x ..... sin n x , cos n x は区間 [- π \pi π , π \pi π ] 上で直交します。
Copy { ∫ − π π cos ( n x ) sin ( m x ) d x = 0 , ∫ − π π cos ( n x ) cos ( m x ) d x = 0 , n ≠ m ∫ − π π sin ( n x ) sin ( m x ) d x = 0 n ≠ m \begin{cases}
\int_{-\pi}^{\pi} \cos (nx) \sin (mx) dx = 0, \\
\int_{-\pi}^{\pi} \cos (nx) \cos (mx) dx = 0, & n \neq m \\
\int_{-\pi}^{\pi} \sin (nx) \sin (mx) dx = 0 & n \neq m
\end{cases} ⎩ ⎨ ⎧ ∫ − π π cos ( n x ) sin ( m x ) d x = 0 , ∫ − π π cos ( n x ) cos ( m x ) d x = 0 , ∫ − π π sin ( n x ) sin ( m x ) d x = 0 n = m n = m
上記の積分は三角関数の積和公式を用いて証明できますが、ここではその計算過程の一例を示します。
Copy ∫ − π π cos ( n x ) sin ( m x ) d x = 1 2 ∫ − π π sin ( n + m ) x d x − 1 2 ∫ − π π sin ( n − m ) x d x = 1 2 [ − 1 n + m cos ( n + m ) x ∣ − π π + 1 n − m cos ( n − m ) x ∣ − π π ] = 0 \int_{-\pi}^{\pi} \cos (nx) \sin (mx) dx
= \frac{1}{2} \int_{-\pi}^{\pi} \sin(n+m)x dx - \frac{1}{2} \int_{-\pi}^{\pi} \sin(n-m)x dx \\
= \frac{1}{2} [-\frac{1}{n+m} \cos (n+m)x |_{-\pi}^{\pi} + \frac{1}{n-m} \cos (n-m)x |_{-\pi}^{\pi}] \\
= 0 ∫ − π π cos ( n x ) sin ( m x ) d x = 2 1 ∫ − π π sin ( n + m ) x d x − 2 1 ∫ − π π sin ( n − m ) x d x = 2 1 [ − n + m 1 cos ( n + m ) x ∣ − π π + n − m 1 cos ( n − m ) x ∣ − π π ] = 0
では、n = m の場合はどうなるでしょうか?sin と cos の組み合わせでは影響はありませんが、cos と cos または sin と sin の場合は、同じ三角関数を掛け合わせて積分することになります。以下に cos m x \cos mx cos m x と cos m x \cos mx cos m x の組み合わせの解法を示します。
Copy ∫ − π π cos ( m x ) cos ( m x ) d x = ∫ − π π 1 + cos ( 2 m ) x 2 d x = 1 2 ( x + 1 2 m sin 2 m x ) ∣ − π π = π \int_{-\pi}^{\pi} \cos (mx) \cos (mx) dx = \int_{-\pi}^{\pi} \frac{1 + \cos (2m)x}{2} dx \\
= \frac{1}{2} (x + \frac{1}{2m} \sin{2mx})|_{-\pi}^{\pi}
= \pi ∫ − π π cos ( m x ) cos ( m x ) d x = ∫ − π π 2 1 + cos ( 2 m ) x d x = 2 1 ( x + 2 m 1 sin 2 m x ) ∣ − π π = π
ここで導出した結果を覚えておいてください。次のセクションでこれらの内容が使われます。
周期が 2 π \pi π の関数のフーリエ級数展開#
最初にフーリエ級数の関連公式を直接示しましたが、次に三角関数の直交性を利用して公式内の係数を求める方法を紹介します。
a 0 a_{0} a 0 の求め方#
フーリエ級数展開公式の両辺に − π -\pi − π から π \pi π までの積分を求めます。
Copy ∫ − π π f ( x ) d x = ∫ − π π ( a 0 2 + Σ n = 1 + ∞ ( a n cos n x + b n sin n x ) ) d x = ∫ − π π a 0 2 d x + a n ∫ − π π Σ n = 1 + ∞ 1 cos n x d x + b n ∫ − π π Σ n = 1 + ∞ 1 sin n x d x \int_{-\pi}^{\pi} f(x) dx \\
= \int_{-\pi}^{\pi} (\frac{a_0}{2} + \varSigma_{n=1}^{+\infty}(a_{n}\cos nx + b_{n}\sin nx)) dx \\
= \int_{-\pi}^{\pi} \frac{a_0}{2} dx + a_{n} \int_{-\pi}^{\pi} \varSigma_{n=1}^{+\infty} 1 \cos nx dx + b_{n} \int_{-\pi}^{\pi} \varSigma_{n=1}^{+\infty} 1 \sin nx dx ∫ − π π f ( x ) d x = ∫ − π π ( 2 a 0 + Σ n = 1 + ∞ ( a n cos n x + b n sin n x )) d x = ∫ − π π 2 a 0 d x + a n ∫ − π π Σ n = 1 + ∞ 1 cos n x d x + b n ∫ − π π Σ n = 1 + ∞ 1 sin n x d x
両辺を積分した後、1、 cos n x \cos nx cos n x と 1、 sin n x \sin nx sin n x の組み合わせに注目すると、三角関数の直交性によりこれらの積分結果は必ず 0 になることがわかります。したがって、以下の等式を得ることができます。
Copy ∫ − π π f ( x ) d x = ∫ − π π a 0 2 d x = a 0 2 ∫ − π π 1 d x = a 0 2 x ∣ − π π = a 0 π \int_{-\pi}^{\pi} f(x) dx = \int_{-\pi}^{\pi} \frac{a_0}{2} dx \\
= \frac{a_0}{2} \int_{-\pi}^{\pi} 1 dx \\
= \frac{a_0}{2} x|_{-\pi}^{\pi} \\
= a_{0} \pi ∫ − π π f ( x ) d x = ∫ − π π 2 a 0 d x = 2 a 0 ∫ − π π 1 d x = 2 a 0 x ∣ − π π = a 0 π
したがって、a 0 a_{0} a 0 は次のように求められます。
Copy a 0 = 1 2 ∫ − π π f ( x ) d x a_{0} = \frac{1}{2}\int_{-\pi}^{\pi} f(x) dx a 0 = 2 1 ∫ − π π f ( x ) d x
a n a_{n} a n の求め方#
フーリエ級数展開の両辺に cos m x \cos mx cos m x を掛けて積分します。
Copy ∫ − π π f ( x ) cos m x d x = ∫ − π π ( a 0 2 cos m x + Σ n = 1 + ∞ ( a n cos n x cos m x + b n sin n x cos m x ) ) d x = a 0 2 ∫ − π π 1 cos m x d x + a n ∫ − π π Σ n = 1 + ∞ cos n x cos m x d x + b n ∫ − π π Σ n = 1 + ∞ sin n x cos m x d x \int_{-\pi}^{\pi} f(x) \cos mx dx \\
= \int_{-\pi}^{\pi} (\frac{a_0}{2} \cos mx + \varSigma_{n=1}^{+\infty}(a_{n}\cos nx \cos mx + b_{n}\sin nx \cos mx)) dx \\
= \frac{a_0}{2} \int_{-\pi}^{\pi} 1 \cos mx dx + a_{n} \int_{-\pi}^{\pi} \varSigma_{n=1}^{+\infty} \cos nx \cos mx dx + b_{n} \int_{-\pi}^{\pi} \varSigma_{n=1}^{+\infty} \sin nx \cos mx dx ∫ − π π f ( x ) cos m x d x = ∫ − π π ( 2 a 0 cos m x + Σ n = 1 + ∞ ( a n cos n x cos m x + b n sin n x cos m x )) d x = 2 a 0 ∫ − π π 1 cos m x d x + a n ∫ − π π Σ n = 1 + ∞ cos n x cos m x d x + b n ∫ − π π Σ n = 1 + ∞ sin n x cos m x d x
この時、三角関数の直交性により a 0 2 ∫ − π π 1 cos m x d x \frac{a_0}{2} \int_{-\pi}^{\pi} 1 \cos mx dx 2 a 0 ∫ − π π 1 cos m x d x と b n ∫ − π π Σ n = 1 + ∞ sin n x cos m x d x b_{n} \int_{-\pi}^{\pi} \varSigma_{n=1}^{+\infty} \sin nx \cos mx dx b n ∫ − π π Σ n = 1 + ∞ sin n x cos m x d x を直接削除できます。最終的にどちらも 0 になるからです。しかし、中央の a n ∫ − π π Σ n = 1 + ∞ cos n x cos m x d x a_{n} \int_{-\pi}^{\pi} \varSigma_{n=1}^{+\infty} \cos nx \cos mx dx a n ∫ − π π Σ n = 1 + ∞ cos n x cos m x d x はそのまま削除できません。なぜなら、m = n m = n m = n の場合、その積分結果は π \pi π になるからです。
したがって、次のようになります。
Copy ∫ − π π f ( x ) cos m x d x = a n π \int_{-\pi}^{\pi} f(x) \cos mx dx = a_{n} \pi ∫ − π π f ( x ) cos m x d x = a n π
変形すると
Copy a n = 1 π ∫ − π π f ( x ) cos m x d x a_{n} = \frac{1}{\pi} \int_{-\pi}^{\pi} f(x) \cos mx dx a n = π 1 ∫ − π π f ( x ) cos m x d x
b n b_{n} b n の求め方#
フーリエ級数展開の両辺に sin m x \sin mx sin m x を掛けて積分します。
Copy ∫ − π π f ( x ) sin m x d x = ∫ − π π ( a 0 2 sin m x + Σ n = 1 + ∞ ( a n cos n x sin m x + b n sin n x sin m x ) ) d x \int_{-\pi}^{\pi} f(x) \sin mx dx = \int_{-\pi}^{\pi} (\frac{a_0}{2} \sin mx + \varSigma_{n=1}^{+\infty}(a_{n}\cos nx \sin mx + b_{n}\sin nx \sin mx)) dx ∫ − π π f ( x ) sin m x d x = ∫ − π π ( 2 a 0 sin m x + Σ n = 1 + ∞ ( a n cos n x sin m x + b n sin n x sin m x )) d x
証明過程は省略しますが、上記の考え方と同じで、次のようになります。
Copy b n = 1 π ∫ − π π f ( x ) sin m x d x b_{n} = \frac{1}{\pi} \int_{-\pi}^{\pi} f(x) \sin mx dx b n = π 1 ∫ − π π f ( x ) sin m x d x
周期が T の関数のフーリエ級数展開#
周期が T の関数の展開方法は少し巧妙で、全体的にはパラメータの変換を通じて式を強制的に π \pi π の形に変えるというものです。
T = 2L とすると、f (t) = f (t + 2L) となります。
x = π L \frac{\pi}{L} L π t とすると、t = L x π \frac{Lx}{\pi} π Lx となり、f (t) = f ( L x π \frac{Lx}{\pi} π Lx ) となります。
g (x) を f ( L π \frac{L}{\pi} π L (x)) と見なすと、上表の対応関係から、g (x) のグラフは f (t) に対応します。例えば、f (2L) = g (2 π \pi π ) という具合です。g (x) の周期は 2 π \pi π ですが、前のセクションで周期が 2 π \pi π の関数のフーリエ展開の係数をすでに得ています。
Copy { a 0 = 1 2 ∫ − π π g ( x ) d x a n = 1 π ∫ − π π g ( x ) cos n x d x b n = 1 π ∫ − π π g ( x ) sin n x d x \begin{cases}
a_{0} = \frac{1}{2}\int_{-\pi}^{\pi} g(x) dx \\
a_{n} = \frac{1}{\pi} \int_{-\pi}^{\pi} g(x) \cos nx dx \\
b_{n} = \frac{1}{\pi} \int_{-\pi}^{\pi} g(x) \sin nx dx
\end{cases} ⎩ ⎨ ⎧ a 0 = 2 1 ∫ − π π g ( x ) d x a n = π 1 ∫ − π π g ( x ) cos n x d x b n = π 1 ∫ − π π g ( x ) sin n x d x
次に、x = π L \frac{\pi}{L} L π t を代入します。
Copy { cos n x = cos n π L t sin n x = sin n π L t ∫ − π π 1 d x = ∫ − L L 1 d π t L g ( x ) = f ( t ) \begin{cases}
\cos nx = \cos \frac{n \pi}{L}t \\
\sin nx = \sin \frac{n \pi}{L}t \\
\int_{-\pi}^{\pi} 1 dx = \int_{-L}^{L} 1 d\frac{\pi t}{L} \\
g(x) = f(t)
\end{cases} ⎩ ⎨ ⎧ cos n x = cos L nπ t sin n x = sin L nπ t ∫ − π π 1 d x = ∫ − L L 1 d L π t g ( x ) = f ( t )
これにより、周期が T のフーリエ級数展開が得られます。
Copy f ( t ) = a 0 2 + Σ n = 1 + ∞ ( a n cos n π L t + b n sin n π L t ) f(t) = \frac{a_0}{2} + \varSigma_{n=1}^{+\infty}(a_{n}\cos \frac{n \pi}{L}t + b_{n}\sin \frac{n \pi}{L}t) f ( t ) = 2 a 0 + Σ n = 1 + ∞ ( a n cos L nπ t + b n sin L nπ t )
ここで(計算過程は省略、解法は前のセクションと同じです)
Copy { a 0 = 1 L ∫ − L L f ( t ) d t a n = 1 L ∫ − L L f ( t ) cos n π L t d t b n = 1 L ∫ − L L f ( t ) sin n π L t d t \begin{cases}
a_{0} = \frac{1}{L}\int_{-L}^{L} f(t) dt \\
a_{n} = \frac{1}{L} \int_{-L}^{L} f(t) \cos \frac{n \pi}{L} t dt \\
b_{n} = \frac{1}{L} \int_{-L}^{L} f(t) \sin \frac{n \pi}{L} t dt
\end{cases} ⎩ ⎨ ⎧ a 0 = L 1 ∫ − L L f ( t ) d t a n = L 1 ∫ − L L f ( t ) cos L nπ t d t b n = L 1 ∫ − L L f ( t ) sin L nπ t d t
L を T に置き換えると、次のようになります。
Copy f ( t ) = a 0 2 + Σ n = 1 + ∞ ( a n cos 2 n π T t + b n sin 2 n π T t ) f(t) = \frac{a_0}{2} + \varSigma_{n=1}^{+\infty}(a_{n}\cos \frac{2n \pi}{T}t + b_{n}\sin \frac{2n \pi}{T}t) f ( t ) = 2 a 0 + Σ n = 1 + ∞ ( a n cos T 2 nπ t + b n sin T 2 nπ t )
Copy { a 0 = 2 T ∫ 0 T f ( t ) d t a n = 2 T ∫ 0 T f ( t ) cos 2 n π T t d t b n = 2 T ∫ 0 T f ( t ) sin 2 n π T t d t \begin{cases}
a_{0} = \frac{2}{T}\int_{0}^{T} f(t) dt \\
a_{n} = \frac{2}{T} \int_{0}^{T} f(t) \cos \frac{2 n \pi}{T} t dt \\
b_{n} = \frac{2}{T} \int_{0}^{T} f(t) \sin \frac{2 n \pi}{T} t dt
\end{cases} ⎩ ⎨ ⎧ a 0 = T 2 ∫ 0 T f ( t ) d t a n = T 2 ∫ 0 T f ( t ) cos T 2 nπ t d t b n = T 2 ∫ 0 T f ( t ) sin T 2 nπ t d t
フーリエ級数の複素形式#
このセクションではオイラーの公式を使用します。オイラーの公式の証明については、後で別の記事を書きます(証明にはテイラー展開を用いる必要がありますが、オイラー自身はこの方法で得たわけではありません)。
Copy e i θ = cos θ + i sin θ e^{i\theta} = \cos \theta + i \sin \theta e i θ = cos θ + i sin θ
Copy { e i θ = cos θ + i sin θ e i − θ = cos − θ + i sin − θ = cos θ − i sin θ \begin{cases}
e^{i\theta} = \cos \theta + i \sin \theta \\
e^{i-\theta} = \cos -\theta + i \sin -\theta = \cos \theta - i \sin \theta
\end{cases} { e i θ = cos θ + i sin θ e i − θ = cos − θ + i sin − θ = cos θ − i sin θ
上式から得られるのは
Copy { sin θ = e i θ − e − i θ 2 i = − i e i θ − e − i θ 2 cos θ = e i θ + e − i θ 2 \begin{cases}
\sin \theta = \frac{e^{i \theta} - e^{-i \theta}}{2i} = -i \frac{e^{i \theta} - e^{-i \theta}}{2} \\
\cos \theta = \frac{e^{i \theta} + e^{- i \theta}}{2}
\end{cases} { sin θ = 2 i e i θ − e − i θ = − i 2 e i θ − e − i θ cos θ = 2 e i θ + e − i θ
x = 2 π t T x = \frac{2 \pi t}{T} x = T 2 π t と置き、この式をフーリエ級数に代入すると
Copy f ( x ) = a 0 2 + Σ n = 1 + ∞ ( a n e i n x + e − i n x 2 − i b n e i n x − e − i n x 2 ) = a 0 2 + Σ n = 1 + ∞ a n − i b n 2 e i n x + Σ n = 1 + ∞ a n + i b n 2 e − i n x f(x) = \frac{a_0}{2} + \varSigma_{n=1}^{+\infty}(a_{n} \frac{e^{i nx} + e^{- i nx}}{2} -i b_{n} \frac{e^{i nx} - e^{-i nx}}{2}) \\
= \frac{a_0}{2} + \varSigma_{n=1}^{+\infty} \frac{a_{n} - i b_{n}}{2} e^{i nx} + \varSigma_{n=1}^{+\infty} \frac{a_{n}+i b_{n}}{2} e^{-i nx} f ( x ) = 2 a 0 + Σ n = 1 + ∞ ( a n 2 e in x + e − in x − i b n 2 e in x − e − in x ) = 2 a 0 + Σ n = 1 + ∞ 2 a n − i b n e in x + Σ n = 1 + ∞ 2 a n + i b n e − in x
次の変換は少し難解ですが、理解できると「なるほど、そういうことか」と感じるでしょう。
a 0 2 \frac{a_{0}}{2} 2 a 0 を Σ n = 0 0 a 0 2 e i n x \varSigma_{n=0}^{0} \frac{a_{0}}{2} e^{i nx} Σ n = 0 0 2 a 0 e in x と見なすことができます。なぜなら、n = 0 のとき、 e i n x e^{i nx} e in x も 1 になるため、成立します。
Σ n = 1 + ∞ e − i n x \varSigma_{n=1}^{+\infty} e^{-i nx} Σ n = 1 + ∞ e − in x の変形も面白く、n を -1 から - ∞ \infty ∞ と見なすと、式は Σ n = − 1 − ∞ e i n x \varSigma_{n=-1}^{-\infty} e^{i nx} Σ n = − 1 − ∞ e in x になります。
したがって、次のようになります。
Copy f ( x ) = Σ n = − ∞ + ∞ C n e i n x f(x) = \varSigma_{n = -\infty}^{+\infty} C_{n} e^{i nx} f ( x ) = Σ n = − ∞ + ∞ C n e in x
x を置き換えると
Copy f ( t ) = Σ n = − ∞ + ∞ C n e i 2 n π T t f(t) = \varSigma_{n = -\infty}^{+\infty} C_{n} e^{i \frac{2 n \pi}{T} t} f ( t ) = Σ n = − ∞ + ∞ C n e i T 2 nπ t
ここで
Copy C n = { a 0 2 , n = 0 a n − i b n 2 n = 1 , 2 , 3... a − n + i b − n 2 n = − 1 , − 2 , − 3... Cn =
\begin{cases}
\frac{a_{0}}{2}, & n = 0 \\
\frac{a_{n} - i b_{n}}{2} & n = 1,2,3... \\
\frac{a_{-n}+i b_{-n}}{2} & n = -1, -2, -3...
\end{cases} C n = ⎩ ⎨ ⎧ 2 a 0 , 2 a n − i b n 2 a − n + i b − n n = 0 n = 1 , 2 , 3... n = − 1 , − 2 , − 3...
前のセクションで周期 T 関数のフーリエ級数展開を求めたので、
Copy { a 0 = 2 T ∫ 0 T f ( t ) d t a n = 2 T ∫ 0 T f ( t ) cos 2 n π T t d t b n = 2 T ∫ 0 T f ( t ) sin 2 n π T t d t \begin{cases}
a_{0} = \frac{2}{T}\int_{0}^{T} f(t) dt \\
a_{n} = \frac{2}{T} \int_{0}^{T} f(t) \cos \frac{2 n \pi}{T} t dt \\
b_{n} = \frac{2}{T} \int_{0}^{T} f(t) \sin \frac{2 n \pi}{T} t dt
\end{cases} ⎩ ⎨ ⎧ a 0 = T 2 ∫ 0 T f ( t ) d t a n = T 2 ∫ 0 T f ( t ) cos T 2 nπ t d t b n = T 2 ∫ 0 T f ( t ) sin T 2 nπ t d t
代入すると
Copy C n = { 1 T ∫ 0 T f ( t ) d t , n = 0 1 T ∫ 0 T f ( t ) ( cos 2 n π T t − i sin 2 n π T t ) d t n > 0 1 T ∫ 0 T f ( t ) ( cos 2 n π T t − i sin 2 n π T t ) d t n < 0 Cn =
\begin{cases}
\frac{1}{T}\int_{0}^{T} f(t) dt, & n = 0 \\
\frac{1}{T} \int_{0}^{T} f(t) (\cos \frac{2n\pi}{T} t - i \sin \frac{2n\pi}{T} t) dt & n > 0 \\
\frac{1}{T} \int_{0}^{T} f(t) (\cos \frac{2n\pi}{T} t - i \sin \frac{2n\pi}{T} t) dt & n < 0
\end{cases} C n = ⎩ ⎨ ⎧ T 1 ∫ 0 T f ( t ) d t , T 1 ∫ 0 T f ( t ) ( cos T 2 nπ t − i sin T 2 nπ t ) d t T 1 ∫ 0 T f ( t ) ( cos T 2 nπ t − i sin T 2 nπ t ) d t n = 0 n > 0 n < 0
ただし、表現が少し煩雑なので、オイラーの公式を使って簡略化します。
Copy cos 2 n π T − i sin 2 n π T = cos − 2 n π T + i sin − 2 n π T = e − i 2 n π T \cos \frac{2n\pi}{T} - i \sin \frac{2n\pi}{T} = \cos -\frac{2n\pi}{T} + i \sin -\frac{2n\pi}{T} = e^{-i \frac{2n\pi}{T}} cos T 2 nπ − i sin T 2 nπ = cos − T 2 nπ + i sin − T 2 nπ = e − i T 2 nπ
Copy C n = { 1 T ∫ 0 T f ( t ) d t , n = 0 1 T ∫ 0 T f ( t ) e − i 2 n π T t d t n > 0 1 T ∫ 0 T f ( t ) e − i 2 n π T t d t n < 0 Cn =
\begin{cases}
\frac{1}{T}\int_{0}^{T} f(t) dt, & n = 0 \\
\frac{1}{T} \int_{0}^{T} f(t) e^{-i \frac{2n\pi}{T} t} dt & n > 0 \\
\frac{1}{T} \int_{0}^{T} f(t) e^{-i \frac{2n\pi}{T} t} dt & n < 0
\end{cases} C n = ⎩ ⎨ ⎧ T 1 ∫ 0 T f ( t ) d t , T 1 ∫ 0 T f ( t ) e − i T 2 nπ t d t T 1 ∫ 0 T f ( t ) e − i T 2 nπ t d t n = 0 n > 0 n < 0
注意してほしいのは、n > 0 と n <0 の結果が同じであることです。また、n = 0 の場合も、n> 0 と n < 0 の表現と同じであることがわかります(n=0 を代入すると e − i 2 n π T e^{-i \frac{ 2n \pi }{T}} e − i T 2 nπ の結果は 1 になります)。
したがって、C n C_{n} C n の表現は次のように統一できます。
Copy C n = 1 T ∫ 0 T f ( t ) e − i 2 n π T t d t C_{n} = \frac{1}{T} \int_{0}^{T} f(t) e^{-i \frac{ 2n \pi }{T} t} dt C n = T 1 ∫ 0 T f ( t ) e − i T 2 nπ t d t
この時、ω 0 \omega_{0} ω 0 を 2 π T \frac{ 2 \pi }{T} T 2 π (工学上の角速度)とすると、
Copy C n = 1 T ∫ 0 T f ( t ) e − i n ω 0 t d t C_{n} = \frac{1}{T} \int_{0}^{T} f(t) e^{-i n \omega_{0} t} dt C n = T 1 ∫ 0 T f ( t ) e − in ω 0 t d t
これにより、フーリエ級数の複素形式が得られます。
Copy { f ( t ) = Σ n = − ∞ + ∞ C n e i n ω 0 t C n = 1 T ∫ 0 T f ( t ) e − i n ω 0 t d t \begin{cases}
f(t) = \varSigma_{n = -\infty}^{+\infty} C_{n} e^{i n \omega_{0} t} \\
C_{n} = \frac{1}{T} \int_{0}^{T} f(t) e^{-i n \omega_{0} t} dt
\end{cases} { f ( t ) = Σ n = − ∞ + ∞ C n e in ω 0 t C n = T 1 ∫ 0 T f ( t ) e − in ω 0 t d t
次に、小さな変換を行います。
積分の上下限を [-T/2, T/2] に置き換え、元の式に代入すると(周期関数の定積分の上下限が 1 周期の差であれば、得られる結果は同じです)
Copy { f ( t ) = Σ n = − ∞ + ∞ C n e i n ω 0 t C n = 1 T ∫ − T 2 T 2 f ( t ) e − i n ω 0 t d t \begin{cases}
f(t) = \varSigma_{n = -\infty}^{+\infty} C_{n} e^{i n \omega_{0} t} \\
C_{n} = \frac{1}{T} \int_{-\frac{T}{2}}^{\frac{T}{2}} f(t) e^{-i n \omega_{0} t} dt
\end{cases} { f ( t ) = Σ n = − ∞ + ∞ C n e in ω 0 t C n = T 1 ∫ − 2 T 2 T f ( t ) e − in ω 0 t d t
非周期関数のフーリエ級数展開#
非周期関数のフーリエ級数展開の処理も非常に巧妙です。大まかな考え方は、非周期関数を 周期が無限 の周期関数と見なし、その後周期関数の展開公式を用いて変換するというものです。
1 T \frac{1}{T} T 1 を f f f とし、T = ∞ \infty ∞ であるため、f f f は無限小に近づきます。したがって、
Copy f ( t ) = Σ n = − ∞ + ∞ ( ∫ − ∞ ∞ f ( t ) e − i n ω 0 t d t ) e i n ω 0 t f f(t) = \varSigma_{n = -\infty}^{+\infty} (\int_{-\infty}^{\infty} f(t) e^{-i n \omega_{0} t} dt) e^{i n \omega_{0} t} f f ( t ) = Σ n = − ∞ + ∞ ( ∫ − ∞ ∞ f ( t ) e − in ω 0 t d t ) e in ω 0 t f
見覚えのある形ですね(定積分の定義)。f f f が無限に近づくため、ω 0 = 2 π T = 2 π f \omega_{0} = \frac{ 2 \pi }{T} = 2 \pi f ω 0 = T 2 π = 2 π f に基づいて 和を積分に変える ことができます。
Copy f ( t ) = ∫ − ∞ ∞ ( ∫ − ∞ ∞ f ( t ) e − i 2 π f t d t ) e i 2 π f t d f f(t) = \int_{-\infty}^{\infty} (\int_{-\infty}^{\infty} f(t) e^{-i 2 \pi f t} dt) e^{i 2 \pi f t} df f ( t ) = ∫ − ∞ ∞ ( ∫ − ∞ ∞ f ( t ) e − i 2 π f t d t ) e i 2 π f t df
ただし、皆さんはおそらく次の形式の方が馴染みがあるでしょう。
Copy f ( t ) = 1 2 π ∫ − ∞ ∞ ( ∫ − ∞ ∞ f ( t ) e − i ω t d t ) e i ω t d ω f(t) = \frac{1}{ 2 \pi } \int_{-\infty}^{\infty} (\int_{-\infty}^{\infty} f(t) e^{-i \omega t} dt) e^{i \omega t} d \omega f ( t ) = 2 π 1 ∫ − ∞ ∞ ( ∫ − ∞ ∞ f ( t ) e − iω t d t ) e iω t d ω
二つの表現方法は実際には同じで、単に一つは周波数から出発し、もう一つは角速度から出発するということです。
ここで
Copy F ( ω ) = ∫ − ∞ ∞ f ( t ) e − i ω t d t F(\omega) = \int_{-\infty}^{\infty} f(t) e^{-i \omega t} dt F ( ω ) = ∫ − ∞ ∞ f ( t ) e − iω t d t
はフーリエ変換です。
また、
Copy f ( t ) = 1 2 π ∫ − ∞ ∞ F ( ω ) e i ω t d ω f(t) = \frac{1}{ 2 \pi } \int_{-\infty}^{\infty} F(\omega) e^{i \omega t} d \omega f ( t ) = 2 π 1 ∫ − ∞ ∞ F ( ω ) e iω t d ω
はフーリエ逆変換です。
参考資料#
純干貨数学推導_フーリエ級数とフーリエ変換_Part5_フーリエ級数からフーリエ変換を導出
IoT の最前線実践 - フーリエ級数
Bilibili 初公開!ミジンコでも理解できる【フーリエ変換】の解説、清華大学の李永楽先生がフーリエ変換の理解方法を教えます。美顔と声変換の原理を見分ける...
フーリエ変換の視覚的表示