傅里叶级数到傅里叶变换

1. 欧拉公式

欧拉公式是把复指函数和三角函数联系起来的函数,建立了三角函数和指数函数的关联:
e i x = c o s x + i s i n x e^{ix} = cosx + isinx eix=cosx+isinx
它是根据泰勒公式观察而来,证明如下:
s i n x = ∑ k = 0 ∞ ( − 1 ) k ⋅ x 2 k + 1 ( 2 k + 1 ) ! = x − x 3 3 ! + x 5 5 ! − x 7 7 ! … sinx=\sum_{k=0}^{\infty}(-1)^k \cdot \frac{x^{2k+1}}{(2k+1)!}=x - \frac{x^3}{3!}+\frac{x^5}{5!}-\frac{x^7}{7!} \dots sinx=k=0(1)k(2k+1)!x2k+1=x3!x3+5!x57!x7
c o s x = ∑ k = 0 ∞ ( − 1 ) k ⋅ x 2 k ( 2 k ) ! = 1 − x 2 2 ! + x 4 4 ! − x 6 6 ! … cosx=\sum_{k=0}^{\infty}(-1)^k \cdot \frac{x^{2k}}{(2k)!}=1 - \frac{x^2}{2!}+\frac{x^4}{4!}-\frac{x^6}{6!} \dots cosx=k=0(1)k(2k)!x2k=12!x2+4!x46!x6
e x = ∑ k = 0 ∞ ( − 1 ) k ⋅ x k ( k ) ! = 1 + x + x 2 2 ! + x 3 3 ! + x 4 4 ! … e^x=\sum_{k=0}^{\infty}(-1)^k \cdot \frac{x^{k}}{(k)!}=1 +x + \frac{x^2}{2!}+\frac{x^3}{3!}+\frac{x^4}{4!} \dots ex=k=0(1)k(k)!xk=1+x+2!x2+3!x3+4!x4
e i x = 1 + i x − x 2 2 ! − x 3 3 ! i + x 4 4 ! ⋯ e^{ix}=1 + ix -\frac{x^2}{2!}-\frac{x^3}{3!}i+\frac{x^4}{4!}\cdots eix=1+ix2!x23!x3i+4!x4
根据观察即可得到欧拉公式。我们还可以根据欧拉公式可以得到以下的结论:

{ e − i x = c o s x − i s i n x ( 1 ) e i x = c o s x + i s i n x ( 2 ) \left\{ \begin{aligned} &e^{-ix} & = cosx-isinx \qquad (1)\\ &e^{ix} & = cosx+isinx \qquad (2)\\ \end{aligned} \right. { eixeix=cosxisinx(1)=cosx+isinx(2)
式(1)减去式(2)可以得到 s i n x = e i x − e − i x 2 i sinx=\frac{e^{ix}-e^{-ix}}{2i} sinx=2ieixeix,式(1)加上式(2)可以得到 c o s x = e i x + e − i x 2 cosx=\frac{e^{ix}+e^{-ix}}{2} cosx=2eix+eix
x = π x=\pi x=π时,可以得到欧拉恒等式:
e i π + 1 = 0 e^{i\pi}+1=0 eiπ+1=0

2. 三角函数系

定义三角函数系:
R = { 0 , 1 , s i n x , c o s x , s i n 2 x , c o s 2 x , ⋯   , s i n n x , c o s n x } \mathbb{R}=\{0, 1, sinx, cosx, sin2x, cos2x, \cdots, sinnx, cosnx\} R={ 0,1,sinx,cosx,sin2x,cos2x,,sinnx,cosnx}
在三角函数系中有这样的正交关系:
{ ∫ − π π sin ⁡ n x cos ⁡ m x = 0 ∫ − π π cos ⁡ n x cos ⁡ m x = 0 , n ≠ m ∫ − π π sin ⁡ n x sin ⁡ m x = 0 , n ≠ m \left\{ \begin{aligned} &\int_{-\pi}^{\pi}\sin{nx}\cos{mx}=0\\ &\int_{-\pi}^{\pi}\cos{nx}\cos{mx}=0, \quad n\neq m\\ &\int_{-\pi}^{\pi}\sin{nx}\sin{mx}=0, \quad n\neq m\\ \end{aligned} \right. ππsinnxcosmx=0ππcosnxcosmx=0,n=mππsinnxsinmx=0,n=m

3. 傅里叶级数

3.1 周期为 2 π 2\pi 2π的函数展开

定义 f ( x ) = f ( x + 2 π ) f(x)=f(x+2\pi) f(x)=f(x+2π),那么其傅里叶级数展开为:
f ( x ) = ∑ 0 ∞ a n cos ⁡ n x + ∑ 0 ∞ b n sin ⁡ n x = a 0 + ∑ 1 ∞ a n cos ⁡ n x + ∑ 1 ∞ b n sin ⁡ n x ( 1 ) \begin{aligned} f(x) &= \sum_{0}^{\infty}a_n\cos{nx}+\sum_{0}^{\infty}b_n\sin{nx} \\ &=a_0+\sum_{1}^{\infty}a_n\cos{nx}+\sum_{1}^{\infty}b_n\sin{nx} \qquad (1) \end{aligned} f(x)=0ancosnx+0bnsinnx=a0+1ancosnx+1bnsinnx(1)

  1. 计算 a 0 a_0 a0
    只需对式(1)两边做积分, 并且根据三角函数系的正交可得:
    ∫ − π π f ( x ) d x = 2 π a 0 + a n ∫ − π π ∑ n = 1 ∞ cos ⁡ n x d x + b n ∫ − π π ∑ n = 1 ∞ sin ⁡ n x d x = 2 π a 0 \begin{aligned} \int_{-\pi}^{\pi}f(x)dx &=2\pi a_0 + a_n\int_{-\pi}^{\pi}\sum_{n=1}^{\infty}\cos{nx}dx+b_n\int_{-\pi}^{\pi}\sum_{n=1}^{\infty}\sin{nx}dx\\ &=2\pi a_0 \end{aligned} ππf(x)dx=2πa0+anππn=1cosnxdx+bnππn=1sinnxdx=2πa0
    那么, a 0 = 1 2 π ∫ − π π f ( x ) d x a_0=\frac{1}{2\pi}\int_{-\pi}^{\pi}f(x)dx a0=2π1ππf(x)dx
    为了公式的统一,通常我们使 a 0 = 1 π ∫ − π π f ( x ) d x a_0=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x)dx a0=π1ππf(x)dx,因此式(1)可以转化为:
    f ( x ) = a 0 2 + ∑ 1 ∞ a n cos ⁡ n x + ∑ 1 ∞ b n sin ⁡ n x ( 2 ) , f(x)=\frac{a_0}{2}+\sum_{1}^{\infty}a_n\cos{nx}+\sum_{1}^{\infty}b_n\sin{nx} \qquad (2), f(x)=2a0+1ancosnx+1bnsinnx(2),
    这也是最为常见的表达形式。

  2. 计算 a n a_n an
    式(2)乘以 cos ⁡ m x \cos{mx} cosmx再做积分:
    ∫ − π π f ( x ) cos ⁡ m x d x = ∫ − π π a 0 2 cos ⁡ m x d x + ∫ − π π ∑ n = 1 ∞ a n cos ⁡ n x cos ⁡ m x d x + ∫ − π π ∑ n = 1 ∞ b n sin ⁡ n x cos ⁡ m x d x = ∫ − π π ∑ n = 1 ∞ a n cos ⁡ n x cos ⁡ m x d x \begin{aligned} \int_{-\pi}^{\pi}f(x) \cos{mx}dx \\ &= \int_{-\pi}^{\pi}\frac{a_0}{2}\cos{mx}dx+\int_{-\pi}^{\pi}\sum_{n=1}^{\infty}a_n\cos{nx}\cos{mx}dx+\int_{-\pi}^{\pi}\sum_{n=1}^{\infty}b_n\sin{nx}\cos{mx}dx\\ &=\int_{-\pi}^{\pi}\sum_{n=1}^{\infty}a_n\cos{nx}\cos{mx}dx \end{aligned} ππf(x)cosmxdx=ππ2a0cosmxdx+ππn=1ancosnxcosmxdx+ππn=1bnsinnxcosmxdx=ππn=1ancosnxcosmxdx
    根据正交性可知,当 n = m n = m n=m时:
    ∫ − π π ∑ n = 1 ∞ a n cos ⁡ n x cos ⁡ m x d x = ∫ − π π a n cos ⁡ 2 n x d x = a n π \int_{-\pi}^{\pi}\sum_{n=1}^{\infty}a_n\cos{nx}\cos{mx}dx=\int_{-\pi}^{\pi}a_n\cos^2{nx}dx=a_n\pi ππn=1ancosnxcosmxdx=ππancos2nxdx=anπ
    因此,可得:
    a n = 1 π ∫ − π π f ( x ) cos ⁡ n x d x a_n=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x) \cos{nx}dx an=π1ππf(x)cosnxdx

  3. 计算 b n b_n bn
    同理,式(2)乘以 sin ⁡ m x \sin{mx} sinmx再做积分:
    ∫ − π π f ( x ) sin ⁡ m x d x = ∫ − π π a 0 2 sin ⁡ m x d x + ∫ − π π ∑ n = 1 ∞ a n cos ⁡ n x sin ⁡ m x d x + ∫ − π π ∑ n = 1 ∞ b n sin ⁡ n x sin ⁡ m x d x = ∫ − π π ∑ n = 1 ∞ b n sin ⁡ n x sin ⁡ m x d x \begin{aligned} \int_{-\pi}^{\pi}f(x) \sin{mx}dx \\ &= \int_{-\pi}^{\pi}\frac{a_0}{2}\sin{mx}dx+\int_{-\pi}^{\pi}\sum_{n=1}^{\infty}a_n\cos{nx}\sin{mx}dx+\int_{-\pi}^{\pi}\sum_{n=1}^{\infty}b_n\sin{nx}\sin{mx}dx\\ &=\int_{-\pi}^{\pi}\sum_{n=1}^{\infty}b_n\sin{nx}\sin{mx}dx \end{aligned} ππf(x)sinmxdx=ππ2a0sinmxdx+ππn=1ancosnxsinmxdx+ππn=1bnsinnxsinmxdx=ππn=1bnsinnxsinmxdx
    根据正交性可知,当 n = m n = m n=m时:
    ∫ − π π ∑ n = 1 ∞ a n sin ⁡ n x sin ⁡ m x d x = ∫ − π π a n sin ⁡ 2 n x d x = b n π \int_{-\pi}^{\pi}\sum_{n=1}^{\infty}a_n\sin{nx}\sin{mx}dx=\int_{-\pi}^{\pi}a_n\sin^2{nx}dx=b_n\pi ππn=1ansinnxsinmxdx=ππansin2nxdx=bnπ
    因此,可得:
    b n = 1 π ∫ − π π f ( x ) sin ⁡ n x d x b_n=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x) \sin{nx}dx bn=π1ππf(x)sinnxdx

综上所述,我们可以得到 T = 2 π T=2\pi T=2π的函数傅里叶级数展开为:
f ( x ) = a 0 2 + ∑ 1 ∞ a n cos ⁡ n x + ∑ 1 ∞ b n sin ⁡ n x f(x)=\frac{a_0}{2}+\sum_{1}^{\infty}a_n\cos{nx}+\sum_{1}^{\infty}b_n\sin{nx} f(x)=2a0+1ancosnx+1bnsinnx
其中:
{ a 0 = 1 π ∫ − π π f ( x ) d x a n = 1 π ∫ − π π f ( x ) cos ⁡ n x d x b n = 1 π ∫ − π π f ( x ) sin ⁡ n x d x \left\{ \begin{aligned} &a_0=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x)dx\\ &a_n=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x) \cos{nx}dx\\ &b_n=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x) \sin{nx}dx\\ \end{aligned} \right. a0=π1ππf(x)dxan=π1ππf(x)cosnxdxbn=π1ππf(x)sinnxdx

3.2 周期为 2 L 2L 2L的函数展开

定义 f ( t ) = f ( t + 2 L ) f(t)=f(t+2L) f(t)=f(t+2L),对其换元,令 x = π t L x=\frac{\pi t}{L} x=Lπt:
f ( t ) = f ( L x π ) = g ( x ) f(t)=f(\frac{Lx}{\pi})=g(x) f(t)=f(πLx)=g(x)
因此,可以将 x = π t L x=\frac{\pi t}{L} x=Lπt代入上一节推导得出的展开式:
f ( t ) = a 0 2 + ∑ 1 ∞ a n cos ⁡ n π t L + ∑ 1 ∞ b n sin ⁡ n π t L f(t)=\frac{a_0}{2}+\sum_{1}^{\infty}a_n\cos{\frac{n\pi t}{L}}+\sum_{1}^{\infty}b_n\sin{\frac{n\pi t}{L}} f(t)=2a0+1ancosLnπt+1bnsinLnπt
其中:
{ a 0 = 1 L ∫ − L L f ( t ) d t a n = 1 L ∫ − L L f ( t ) cos ⁡ n π t L d t b n = 1 L ∫ − L L f ( t ) sin ⁡ n π t L d t \left\{ \begin{aligned} &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 t}{L}}dt\\ &b_n=\frac{1}{L}\int_{-L}^{L}f(t) \sin{\frac{n\pi t}{L}}dt\\ \end{aligned} \right. a0=L1LLf(t)dtan=L1LLf(t)cosLnπtdtbn=L1LLf(t)sinLnπtdt

3.3 总结

通常, t t t从0开始, T = 2 L T=2L T=2L ω = π L = 2 π T \omega=\frac{\pi}{L}=\frac{2\pi}{T} ω=Lπ=T2π:
f ( t ) = a 0 2 + ∑ 1 ∞ a n cos ⁡ n ω t + ∑ 1 ∞ b n sin ⁡ n ω t f(t)=\frac{a_0}{2}+\sum_{1}^{\infty}a_n\cos{ {n\omega t}}+\sum_{1}^{\infty}b_n\sin{ {n\omega t}} f(t)=2a0+1ancosnωt+1bnsinnωt

{ a 0 = 2 T ∫ 0 T f ( t ) d t a n = 2 T ∫ 0 T f ( t ) cos ⁡ n ω t d t b n = 2 T ∫ 0 T f ( t ) sin ⁡ n ω t d t \left\{ \begin{aligned} &a_0=\frac{2}{T}\int_{0}^{T}f(t)dt\\ &a_n=\frac{2}{T}\int_{0}^{T}f(t) \cos{ {n\omega t}}dt\\ &b_n=\frac{2}{T}\int_{0}^{T}f(t) \sin{ {n\omega t}}dt\\ \end{aligned} \right. a0=T20Tf(t)dtan=T20Tf(t)cosnωtdtbn=T20Tf(t)sinnωtdt

4. 傅里叶级数的复数形式

根据欧拉公式:
{ sin ⁡ θ = − i 2 ( e i θ − e − i θ ) cos ⁡ θ = 1 2 ( e i θ + e − i θ ) \left\{ \begin{aligned} &\sin{\theta}=- \frac{i}{2}(e^{i \theta} - e^{-i \theta})\\ &\cos{\theta}=\frac{1}{2}(e^{i \theta} + e^{-i \theta})\\ \end{aligned} \right. sinθ=2i(eiθeiθ)cosθ=21(eiθ+eiθ)
因此可以将上一章最后得出的傅里叶级数写成复数形式:
f ( t ) = a 0 2 + ∑ n = 1 ∞ a n − i b n 2 e i n ω t + ∑ n = 1 ∞ a n + i b n 2 e − i n ω t f(t)=\frac{a_0}{2}+\sum_{n=1}^{\infty}\frac{a_n-ib_n}{2}e^{in\omega t}+\sum_{n=1}^{\infty}\frac{a_n+ib_n}{2}e^{-in\omega t} f(t)=2a0+n=12anibneinωt+n=12an+ibneinωt
接下来有意思的来了,我们可以将第一项写为求和的形式:
a 0 2 = ∑ n = 0 0 a 0 2 e i n ω t \frac{a_0}{2}=\sum_{n=0}^{0}\frac{a_0}{2}e^{in\omega t} 2a0=n=002a0einωt
再令第三项的 n = − n n=-n n=n:
∑ n = 1 ∞ a n + i b n 2 e − i n ω t = ∑ n = − ∞ − 1 a − n + i b − n 2 e i n ω t \sum_{n=1}^{\infty}\frac{a_n+ib_n}{2}e^{-in\omega t}=\sum_{n=-\infty}^{-1}\frac{a_{-n}+ib_{-n}}{2}e^{in\omega t} n=12an+ibneinωt=n=12an+ibneinωt
我们惊喜地发现这三项竟然可以合并起来了:
f ( t ) = ∑ − ∞ ∞ C n e i n ω t f(t)=\sum_{-\infty}^{\infty}C_ne^{in\omega t} f(t)=Cneinωt
系数 C n C_n Cn有三种情况:
{ C n = a 0 2 , n = 0 C n = a n − i b n 2 , n > 0 C n = a − n + i b − n 2 , n < 0 \left\{ \begin{aligned} &C_n=\frac{a_0}{2}, \qquad n=0\\ &C_n=\frac{a_n-ib_n}{2}, \qquad n \gt 0\\ &C_n=\frac{a_{-n}+ib_{-n}}{2}, \qquad n \lt 0\\ \end{aligned} \right. Cn=2a0,n=0Cn=2anibn,n>0Cn=2an+ibn,n<0
接着,我们可以将 a 0 a_0 a0, a n a_n an b n b_n bn代入可得:
{ C n = 1 T ∫ 0 T f ( t ) d t , n = 0 C n = 1 T ∫ 0 T f ( t ) e − i n ω t d t , n > 0 C n = 1 T ∫ 0 T f ( t ) e − i n ω t d t , n < 0 \left\{ \begin{aligned} &C_n=\frac{1}{T}\int_{0}^{T}f(t)dt, \qquad n=0\\ &C_n=\frac{1}{T}\int_{0}^{T}f(t)e^{-in\omega t}dt, \qquad n \gt 0\\ &C_n=\frac{1}{T}\int_{0}^{T}f(t)e^{-in\omega t}dt, \qquad n \lt 0\\ \end{aligned} \right. Cn=T10Tf(t)dt,n=0Cn=T10Tf(t)einωtdt,n>0Cn=T10Tf(t)einωtdt,n<0
我们可以发现三个式子本质上都是一样的,因此我们可以得出最后的结论,:
f ( t ) = ∑ − ∞ ∞ C n e i n ω t f(t)=\sum_{-\infty}^{\infty}C_ne^{in\omega t} f(t)=Cneinωt
C n = 1 T ∫ 0 T f ( t ) e − i n ω t d t C_n=\frac{1}{T}\int_{0}^{T}f(t)e^{-in\omega t}dt Cn=T10Tf(t)einωtdt

5. 傅里叶变换

假设一个函数的周期趋于无穷( T → ∞ T \to \infty T),此时可以视为 lim ⁡ T → ∞ f T ( t ) = f ( t ) \lim \limits_{T \to \infty}f_T(t)=f(t) TlimfT(t)=f(t),一个周期函数就成了一个周期无限长的普通函数。
我们重新写一下傅里叶级数:
f T ( t ) = ∑ − ∞ ∞ C n e i n ω 0 t f_T(t)=\sum_{-\infty}^{\infty}C_ne^{in\omega_0 t} fT(t)=Cneinω0t
C n = 1 T ∫ − T 2 T 2 f T ( t ) e − i n ω 0 t d t C_n=\frac{1}{T}\int_{-\frac{T}{2}}^{\frac{T}{2}}f_T(t)e^{-in\omega_0 t}dt Cn=T12T2TfT(t)einω0tdt

频率之间的差值 △ ω = ( n + 1 ) ω 0 − n ω 0 = ω 0 = 2 π T \triangle \omega =(n+1)\omega_0-n \omega_0=\omega_0=\frac{2\pi}{T} ω=(n+1)ω0nω0=ω0=T2π,如果周期趋于无穷时,频率和频率之间的距离无限接近,因此从离散成为了连续,即 n ω 0 → ω n\omega_0 \to \omega nω0ω
此外,由于变成了连续值,那么求和就变成了求积分,因此将傅里叶级数改写为:
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}dte^{i \omega t}d\omega f(t)=2π1f(t)eiωtdteiωtdω,
其中我们将 F ( ω ) = ∫ − ∞ ∞ f ( t ) e − i ω t d t F(\omega)=\int_{-\infty}^{\infty}f(t)e^{-i \omega t}dt F(ω)=f(t)eiωtdt称为傅里叶变换,而 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π1F(ω)eiωtdω称为傅里叶逆变换。

猜你喜欢

转载自blog.csdn.net/loki2018/article/details/123949773