傅里叶级数与傅里叶变换的推导

文章首发于:https://zhuanlan.zhihu.com/p/342982724

主要参考自 DR_CAN 的B站教程 纯干货傅里叶变换,某些地方加了一点自己的理解。。

0x01 三角函数系的正交性与傅里叶级数

我们首先看一下这个定理: 组成三角级数的函数系
0 ( sin 0 x ) , 1 ( cos 0 x ) , sin x , cos x , sin 2 x , cos 2 x , … , sin n x , cos n x \begin{aligned} 0(\text{sin}0x),1(\text{cos}0x),\text{sin}x,\text{cos}x,\text{sin}2x,\text{cos}2x,\dots,\text{sin}nx,\text{cos}nx \end{aligned} 0(sin0x),1(cos0x),sinx,cosx,sin2x,cos2x,,sinnx,cosnx
[ − π , π ] [-\pi,\pi] [π,π] 上正交,即其中任意两个不同的函数之积在 [ − π , π ] [-\pi,\pi] [π,π] 上的积分等于0,证明见附录1。

我们可以先记住这个结论,就是从上面的那组三角函数中随便选两个不同的三角函数,将他们的乘起来,然后在 [ − π , π ] [-\pi,\pi] [π,π] 下做积分,结果一定是0。

并且在三角函数系中两个相同的函数的乘积在 [ − π , π ] [-\pi,\pi] [π,π] 上的积分不等于0,且有:
∫ − π π cos n x ⋅ cos m x dx = ∫ − π π cos 2 n x dx = ∫ − π π 1 2 ( 1 + cos 2 n x ) dx = 1 2 [ ∫ − π π 1 dx + ∫ − π π cos 2 n x dx ] = π \begin{aligned} \int_{-\pi}^{\pi} \text{cos}nx·\text{cos}mx\text{dx} &= \int_{-\pi}^{\pi} \text{cos}^2nx\text{dx} \\ &=\int_{-\pi}^{\pi} \frac{1}{2}\left(1 +\text{cos}2nx\right)\text{dx} \\ &=\frac{1}{2} \left[\int_{-\pi}^{\pi}1\text{dx} + \int_{-\pi}^{\pi}\text{cos}2nx\text{dx}\right] \\ &=\pi \\ \end{aligned} ππcosnxcosmxdx=ππcos2nxdx=ππ21(1+cos2nx)dx=21[ππ1dx+ππcos2nxdx]=π
∫ − π π sin 2 n x dx = π ∫ − π π 1 ⋅ 1 dx = 2 π \begin{aligned} \int_{-\pi}^{\pi} \text{sin}^2nx\text{dx} =\pi \\ \int_{-\pi}^{\pi} 1·1\text{dx}=2\pi \end{aligned} ππsin2nxdx=πππ11dx=2π

我们知道简单的周期运动可写作: y = A sin ( ω t + φ ) y=A\text{sin}(\omega t+\varphi) y=Asin(ωt+φ) ,其中 A A A 为振幅, ω \omega ω 为角频率, φ \varphi φ 为初相。下图为 sin ⁡ ( x + π ) \sin\left(x+\pi\right) sin(x+π)

在这里插入图片描述

复杂的周期运动 f ( t ) f(t) f(t) 可以看作简单周期运动的叠加。

下图为: sin ⁡ ( x + π ) + 1 3 sin ⁡ ( 3 x + π ) + 1 5 sin ⁡ ( 5 x + π ) + 1 7 sin ⁡ ( 7 x + π ) + 1 9 sin ⁡ ( 9 x + π ) \sin\left(x+\pi\right)+\frac{1}{3}\sin\left(3x+\pi\right)+\frac{1}{5}\sin\left(5x+\pi\right)+\frac{1}{7}\sin\left(7x+\pi\right)+\frac{1}{9}\sin\left(9x+\pi\right) sin(x+π)+31sin(3x+π)+51sin(5x+π)+71sin(7x+π)+91sin(9x+π)

在这里插入图片描述

当累加的足够多,也就形成了一个周期性的方波。

因此周期函数 f ( t ) f(t) f(t) 的表达式可写作:
f ( t ) = ∑ n = 0 ∞ A n sin ( n ω t + φ n ) = A 0 sin φ 0 + ∑ n = 1 ∞ A n sin ( n ω t + φ n ) = A 0 sin φ 0 + ∑ n = 1 ∞ [ A n sin n ω t cos φ n + A n cos n ω t sin φ n ] \begin{aligned} f(t)&=\sum_{n=0}^{\infty}A_n\text{sin}{(n\omega t+\varphi_n)} \\ &=A_0\text{sin}{\varphi_0}+\sum_{n=1}^{\infty}A_n\text{sin}{(n\omega t+\varphi_n)} \\ &=A_0\text{sin}{\varphi_0}+\sum_{n=1}^{\infty}\left[A_n\text{sin}{n\omega t}\text{cos}{\varphi_n} + A_n\text{cos}{n\omega t}\text{sin}{\varphi_n}\right] \end{aligned} f(t)=n=0Ansin(nωt+φn)=A0sinφ0+n=1Ansin(nωt+φn)=A0sinφ0+n=1[Ansinnωtcosφn+Ancosnωtsinφn]
我们令 a n = A n sin φ n a_n=A_n\text{sin}{\varphi_n} an=Ansinφn b n = A n cos φ n b_n=A_n\text{cos}{\varphi_n} bn=Ancosφn ω t = x \omega t=x ωt=x(这里我们**假设周期为 2 π 2\pi 2π,**即 ω = 1 \omega=1 ω=1),代换后的周期函数表达式如下:
f ( x ) = A 0 sin φ 0 + ∑ n = 1 ∞ A n sin ( n ω t + φ n ) = a 0 + ∑ n = 1 ∞ ( a n cos n x + b n sin n x ) (1) \begin{aligned} f(x)&=A_0\text{sin}{\varphi_0}+\sum_{n=1}^{\infty}A_n\text{sin}{(n\omega t+\varphi_n)} \\ &=a_0+\sum_{n=1}^{\infty}(a_n\text{cos}{nx}+b_n\text{sin}{nx}) \\ \end{aligned} \tag{1} f(x)=A0sinφ0+n=1Ansin(nωt+φn)=a0+n=1(ancosnx+bnsinnx)(1)
实际上,(1)式只是一个猜想式,我们需要找到 a 0 a_0 a0 a n a_n an b n b_n bn 式该表达式成立。

首先求 a n a_n an的表达式,可将(1)式两边同时乘以 c o s m x cosmx cosmx,然后两边同时在 [ − π , π ] [-\pi,\pi] [π,π] 做积分,得到:
∫ − π π f ( x ) cos m x dx = ∫ − π π a 0 cos m x dx + ∫ − π π ∑ n = 1 ∞ a n cos n x cos m x dx + ∫ − π π ∑ n = 1 ∞ b n sin n x cos m x dx = a 0 ∫ − π π cos m x dx + a n ∫ − π π ∑ n = 1 ∞ cos n x cos m x dx + b n ∫ − π π ∑ n = 1 ∞ sin n x cos m x dx (2) \begin{aligned} \int_{-\pi}^{\pi}f(x)\text{cos}mx\text{dx}&=\int_{-\pi}^{\pi}a_0\text{cos}mx\text{dx}+\int_{-\pi}^{\pi}\sum_{n=1}^{\infty}a_n\text{cos}{nx}\text{cos}mx\text{dx}+\int_{-\pi}^{\pi}\sum_{n=1}^{\infty}b_n\text{sin}{nx}\text{cos}mx\text{dx} \\ &=a_0\int_{-\pi}^{\pi}\text{cos}mx\text{dx}+a_n\int_{-\pi}^{\pi}\sum_{n=1}^{\infty}\text{cos}{nx}\text{cos}mx\text{dx}+b_n\int_{-\pi}^{\pi}\sum_{n=1}^{\infty}\text{sin}{nx}\text{cos}mx \text{dx} \end{aligned} \tag{2} ππf(x)cosmxdx=ππa0cosmxdx+ππn=1ancosnxcosmxdx+ππn=1bnsinnxcosmxdx=a0ππcosmxdx+anππn=1cosnxcosmxdx+bnππn=1sinnxcosmxdx(2)
由前面证明的三角函数正交性可得:
∫ − π π cos m x dx = ∫ − π π cos m x sin 0 x dx = 0 ∫ − π π ∑ n = 1 ∞ sin n x cos m x dx = 0 ∫ − π π ∑ n = 1 ∞ cos n x cos m x dx = 0 ( n ≠ m ) ∫ − π π ∑ n = 1 ∞ cos n x cos m x dx = π ( n = m ) \begin{aligned} &\int_{-\pi}^{\pi}\text{cos}mx\text{dx}=\int_{-\pi}^{\pi}\text{cos}mx\text{sin}0x\text{dx}=0 \\ &\int_{-\pi}^{\pi}\sum_{n=1}^{\infty}\text{sin}{nx}\text{cos}mx \text{dx}=0 \\ &\int_{-\pi}^{\pi}\sum_{n=1}^{\infty}\text{cos}{nx}\text{cos}mx\text{dx}=0(n\neq m) \\ &\int_{-\pi}^{\pi}\sum_{n=1}^{\infty}\text{cos}{nx}\text{cos}mx\text{dx}=\pi(n=m) \end{aligned} ππcosmxdx=ππcosmxsin0xdx=0ππn=1sinnxcosmxdx=0ππn=1cosnxcosmxdx=0(n=m)ππn=1cosnxcosmxdx=π(n=m)
因此 (2) 式可以写作如下形式 ( n = m ) (n=m) (n=m)
∫ − π π f ( x ) cos m x dx = ∫ − π π f ( x ) cos n x dx = a n ∫ − π π ∑ n = 1 ∞ cos n x cos m x dx = a n π (3) \begin{aligned} \int_{-\pi}^{\pi}f(x)\text{cos}mx\text{dx}&=\int_{-\pi}^{\pi}f(x)\text{cos}nx\text{dx} \\ &=a_n\int_{-\pi}^{\pi}\sum_{n=1}^{\infty}\text{cos}{nx}\text{cos}mx\text{dx}\\ &=a_n\pi \\ \end{aligned} \tag{3} ππf(x)cosmxdx=ππf(x)cosnxdx=anππn=1cosnxcosmxdx=anπ(3)
可以求得:
a n = 1 π ∫ − π π f ( x ) cos n x dx (4) \begin{aligned} a_n=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x)\text{cos}nx\text{dx} \end{aligned} \tag{4} an=π1ππf(x)cosnxdx(4)
同理,我们将(1)式两边同时乘以 s i n m x sinmx sinmx,然后两边同时在 [ − π , π ] [-\pi,\pi] [π,π] 做积分,可以求得:
b n = 1 π ∫ − π π f ( x ) sin n x dx (5) \begin{aligned} b_n=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x)\text{sin}nx\text{dx} \end{aligned} \tag{5} bn=π1ππf(x)sinnxdx(5)
同理,对(1)式两边同时在 [ − π , π ] [-\pi,\pi] [π,π] 做积分,可以求得:
a 0 = 1 2 π ∫ − π π f ( x ) dx (6) \begin{aligned} a_0=\frac{1}{2\pi}\int_{-\pi}^{\pi}f(x)\text{dx} \end{aligned} \tag{6} a0=2π1ππf(x)dx(6)
对比(4)式和(6)式,不难看出当(4)式的 n = 0 n=0 n=0 时,它和 a 0 a_0 a0 只差了一个 1 2 \frac{1}{2} 21,因此为了将它们统一在一个表达式中,通常将(1)式写作如下的形式:
f ( x ) = a 0 2 + ∑ n = 1 ∞ ( a n cos n x + b n sin n x ) (7) \begin{aligned} f(x)&=\frac{a_0}{2}+\sum_{n=1}^{\infty}(a_n\text{cos}{nx}+b_n\text{sin}{nx}) \\ \end{aligned} \tag{7} f(x)=2a0+n=1(ancosnx+bnsinnx)(7)
最终得到的 (7) 式,就是三角级数的表示形式。当 a n a_n an b n b_n bn 由上面的 (4)(5) 式确定时,(7) 式即为 f ( x ) f(x) f(x)傅里叶级数,称此时的 a n a_n an b n b_n bn傅里叶系数。(这里 f ( x ) f(x) f(x) 周期为 2 π 2\pi 2π

0x02 傅里叶级数的三角形式

上面我们推导傅里叶级数时,假设了周期为 2 π 2\pi 2π ,得到了如下形式的傅里叶级数:
f ( x ) = a 0 2 + ∑ n = 1 ∞ ( a n cos n x + b n sin n x ) { a n = 1 π ∫ − π π f ( x ) cos n x dx b n = 1 π ∫ − π π f ( x ) sin n x dx (8) \begin{aligned} f(x)&=\frac{a_0}{2}+\sum_{n=1}^{\infty}(a_n\text{cos}{nx}+b_n\text{sin}{nx}) \\ &\left\{\begin{array}{ll} a_n=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x)\text{cos}nx\text{dx} \\ b_n=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x)\text{sin}nx\text{dx} \end{array}\right. \end{aligned} \tag {8} f(x)=2a0+n=1(ancosnx+bnsinnx){an=π1ππf(x)cosnxdxbn=π1ππf(x)sinnxdx(8)
现在可以进一步推广至周期为 T T T,即 f ( t ) = f ( t + T ) f(t)=f(t+T) f(t)=f(t+T),我们进行换元 x = 2 π T t x=\frac{2\pi}{T}t x=T2πt t = T 2 π x t=\frac{T}{2\pi}x t=2πTx,然后令 f ( t ) = f ( 2 L π x ) : = g ( x ) f(t)=f(\frac{2L}{\pi}x):= g(x) f(t)=f(π2Lx):=g(x),则有 g ( x ) = g ( x + 2 π ) g(x)=g(x+2\pi) g(x)=g(x+2π)

所以可以将 g ( x ) g(x) g(x) 用 (8) 式傅里叶级数展开:
g ( x ) = a 0 2 + ∑ n = 1 ∞ ( a n cos n x + b n sin n x ) { a n = 1 π ∫ − π π g ( x ) cos n x dx b n = 1 π ∫ − π π g ( x ) sin n x dx (9) \begin{aligned} g(x)&=\frac{a_0}{2}+\sum_{n=1}^{\infty}(a_n\text{cos}{nx}+b_n\text{sin}{nx}) \\ &\left\{\begin{array}{ll} a_n=\frac{1}{\pi}\int_{-\pi}^{\pi}g(x)\text{cos}nx\text{dx} \\ b_n=\frac{1}{\pi}\int_{-\pi}^{\pi}g(x)\text{sin}nx\text{dx} \end{array}\right. \end{aligned} \tag {9} g(x)=2a0+n=1(ancosnx+bnsinnx){an=π1ππg(x)cosnxdxbn=π1ππg(x)sinnxdx(9)
然后对 (9) 式再进行一次换元 x = 2 π T t x=\frac{2\pi}{T}t x=T2πt,得到:
f ( t ) = a 0 2 + ∑ n = 1 ∞ ( a n cos n 2 π L t + b n sin n 2 π L t ) = a 0 2 + ∑ n = 1 ∞ ( a n cos n ω t + b n sin n ω t ) (10) \begin{aligned} f(t)&=\frac{a_0}{2}+\sum_{n=1}^{\infty}(a_n\text{cos}{n\frac{2\pi}{L}t}+b_n\text{sin}{n\frac{2\pi}{L}t}) \\ &=\frac{a_0}{2}+\sum_{n=1}^{\infty}(a_n\text{cos}{n \omega t}+b_n\text{sin}{n \omega t}) \end{aligned} \tag{10} f(t)=2a0+n=1(ancosnL2πt+bnsinnL2πt)=2a0+n=1(ancosnωt+bnsinnωt)(10)
同理对 a n a_n an b n b_n bn 进行换元得到:
a n = 1 π ∫ − π π g ( x ) cos n x dx = 1 π ∫ − T 2 T 2 f ( t ) cos n 2 π T t d 2 π T t = 1 π 2 π T ∫ − T 2 T 2 f ( t ) cos n 2 π T t dt = 2 T ∫ 0 T f ( t ) cos n ω t dt (11) \begin{aligned} a_n&=\frac{1}{\pi}\int_{-\pi}^{\pi}g(x)\text{cos}nx\text{dx} \\ &= \frac{1}{\pi}\int_{-\frac{T}{2}}^{\frac{T}{2}}f(t)\text{cos}{n\frac{2\pi}{T}t}\text{d}\frac{2\pi}{T}t \\ &=\frac{1}{\pi}\frac{2\pi}{T}\int_{-\frac{T}{2}}^{\frac{T}{2}}f(t)\text{cos}{n\frac{2\pi}{T}t}\text{dt} \\ &=\frac{2}{T}\int_{0}^{T}f(t)\text{cos}{n \omega t}\text{dt} \end{aligned} \tag {11} an=π1ππg(x)cosnxdx=π12T2Tf(t)cosnT2πtdT2πt=π1T2π2T2Tf(t)cosnT2πtdt=T20Tf(t)cosnωtdt(11)

b n = 2 T ∫ 0 T f ( t ) sin n ω t dt (12) \begin{aligned} b_n&=\frac{2}{T}\int_{0}^{T}f(t)\text{sin}{n \omega t}\text{dt} \end{aligned} \tag {12} bn=T20Tf(t)sinnωtdt(12)

综上,也就说明了任何周期函数都能用傅里叶级数来表示为三角函数的形式。对于周期函数 f ( t ) f(t) f(t),周期为 T T T ,角频率为 ω 0 = 2 π f 0 = 2 π T \omega_0=2\pi f_0=\frac{2\pi}{T} ω0=2πf0=T2π,其三角形式的傅里叶级数如下:
f ( t ) = a 0 2 + ∑ n = 1 ∞ ( a n cos n ω 0 t + b n sin n ω 0 t ) { a n = 2 T ∫ 0 T f ( t ) cos n ω 0 t dt b n = 2 T ∫ 0 T f ( t ) sin n ω 0 t dt (13) \begin{aligned} f(t)&=\frac{a_0}{2}+\sum_{n=1}^{\infty}(a_n\text{cos}{n \omega_0 t}+b_n\text{sin}{n \omega_0 t}) \\ &\left\{\begin{array}{ll} a_n=\frac{2}{T}\int_{0}^{T}f(t)\text{cos}{n \omega_0 t}\text{dt} \\ b_n=\frac{2}{T}\int_{0}^{T}f(t)\text{sin}{n \omega_0 t}\text{dt} \end{array}\right. \end{aligned} \tag {13} f(t)=2a0+n=1(ancosnω0t+bnsinnω0t){an=T20Tf(t)cosnω0tdtbn=T20Tf(t)sinnω0tdt(13)

0x03 傅里叶级数的复指数形式

上一节我们推导了三角形式的傅里叶级数,其物理含义明确;但是实际操作上,指数形式的傅里叶级会更加便于数学处理,而且很容易与后面介绍的傅里叶变换统一起来。因此这一节我们就用欧拉公式将 (13) 式的傅里叶级数表示为复指数的形式。(欧拉公式的证明可见附录2

首先,由欧拉公式 e i θ = cos θ + i sin θ e^{i\theta}=\text{cos}{\theta}+i\text{sin}{\theta} eiθ=cosθ+isinθ ,分别将 θ \theta θ − θ -\theta θ 代入其中再消元,得到三角函数的复指数形式如下:
cos θ = 1 2 ( e i θ + e − i θ ) sin θ = − 1 2 i ( e i θ − e − i θ ) (14) \begin{aligned} \text{cos}{\theta} = \frac{1}{2}\left(e^{i\theta}+e^{-i\theta}\right) \\ \text{sin}{\theta} = -\frac{1}{2}i\left(e^{i\theta}-e^{-i\theta}\right) \end{aligned} \tag {14} cosθ=21(eiθ+eiθ)sinθ=21i(eiθeiθ)(14)
在这里插入图片描述

将 (14) 式的代入 (13) 式的傅里叶级数中得到:

f ( t ) = a 0 2 + ∑ n = 1 ∞ [ a n 1 2 ( e i n ω t + e − i n ω t ) − b n 1 2 i ( e i n ω t − e − i n ω t ) ] = a 0 2 + ∑ n = 1 ∞ ( a n − i b n 2 e i n ω t + a n + i b n 2 e − i n ω 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 = ∑ n = 0 0 a 0 2 e 0 ⋅ i ω t + ∑ n = 1 ∞ a n − i b n 2 e i n ω t + ∑ n = − 1 − ∞ a − n + i b − n 2 e i n ω t = ∑ n = − ∞ ∞ C n e i n ω t w h e r e   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 , … (15) \begin{aligned} f(t)&=\frac{a_0}{2}+\sum_{n=1}^{\infty}\left[a_n\frac{1}{2}\left(e^{in \omega t}+e^{-in \omega t}\right)-b_n\frac{1}{2}i\left(e^{in \omega t}-e^{-in \omega t}\right)\right] \\ &=\frac{a_0}{2}+\sum_{n=1}^{\infty}\left(\frac{a_n-ib_n}{2}e^{in \omega t}+\frac{a_n+ib_n}{2}e^{-in \omega t}\right) \\ &=\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} \\ &=\sum_{n=0}^{0}\frac{a_0}{2}e^{0·i \omega t}+\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} \\ &=\sum_{n=-\infty}^{\infty} C_n e^{in \omega t} \\ \\ where \space C_n=&\left\{\begin{array}{ll} \frac{a_0}{2},n=0 \\ \frac{a_n-ib_n}{2},n=1,2,3,\dots \\ \frac{a_{-n}+ib_{-n}}{2},n=-1,-2,-3,\dots \end{array}\right. \end{aligned} \tag {15} f(t)where Cn==2a0+n=1[an21(einωt+einωt)bn21i(einωteinωt)]=2a0+n=1(2anibneinωt+2an+ibneinωt)=2a0+n=12anibneinωt+n=12an+ibneinωt=n=002a0e0iωt+n=12anibneinωt+n=12an+ibneinωt=n=Cneinωt2a0,n=02anibn,n=1,2,3,2an+ibn,n=1,2,3,(15)

上面第 4 个等号第三项将 n n n 换元为 − n -n n,并将第一项写成单个累和的形式( e 0 ⋅ i ω t = 1 e^{0·i \omega t}=1 e0iωt=1);最后一个等号,可以看到三项的累和已经从 − ∞ -\infty 加到了 + ∞ +\infty +,我们将 e i n ω t e^{in \omega t} einωt 提出来, 剩余项项用 C n C_n Cn 表示。

进一步计算 C n C_n Cn 的形式,将 (13) 式的傅里叶系数代入 C n C_n Cn 中得到:

(1)当 n = 0 n=0 n=0 时:
C n = a 0 2 = 1 2 2 T ∫ 0 T f ( t ) dt = 1 T ∫ 0 T f ( t ) dt = 1 T ∫ 0 T f ( t ) e − i n ω t dt (16) \begin{aligned} C_n=\frac{a_0}{2}&=\frac{1}{2}\frac{2}{T}\int_{0}^{T}f(t)\text{dt} \\ &=\frac{1}{T}\int_{0}^{T}f(t)\text{dt} \\ &=\frac{1}{T}\int_{0}^{T}f(t)e^{-in \omega t}\text{dt} \end{aligned} \tag{16} Cn=2a0=21T20Tf(t)dt=T10Tf(t)dt=T10Tf(t)einωtdt(16)
上面的最后一个等号,因为在此情况下 n = 0 n=0 n=0,所以加上 e − i n ω t e^{-in \omega t} einωt 不改变结果。

(2)当 n = 1 , 2 , 3 , … n=1,2,3,\dots n=1,2,3, 时:
C n = a n − i b n 2 = 1 2 ( 2 T ∫ 0 T f ( t ) cos n ω t dt − i 2 T ∫ 0 T f ( t ) sin n ω t dt ) = 1 T ∫ 0 T f ( t ) ( cos n ω t − i sin n ω t ) dt = 1 T ∫ 0 T f ( t ) e − i n ω t dt (17) \begin{aligned} C_n&=\frac{a_n-ib_n}{2}\\ &=\frac{1}{2}\left(\frac{2}{T}\int_{0}^{T}f(t)\text{cos}{n \omega t}\text{dt}-i\frac{2}{T}\int_{0}^{T}f(t)\text{sin}{n \omega t}\text{dt}\right) \\ &=\frac{1}{T}\int_{0}^{T}f(t) \left(\text{cos}{n \omega t}-i\text{sin}{n \omega t}\right)\text{dt} \\ &=\frac{1}{T}\int_{0}^{T}f(t)e^{-in \omega t}\text{dt} \\ \end{aligned} \tag{17} Cn=2anibn=21(T20Tf(t)cosnωtdtiT20Tf(t)sinnωtdt)=T10Tf(t)(cosnωtisinnωt)dt=T10Tf(t)einωtdt(17)
上面第3个等式可由欧拉公式得到: cos n ω t − i sin n ω t = cos ( − n ω t ) + i sin ( − n ω t ) = e − i n ω t \text{cos}{n \omega t}-i\text{sin}{n \omega t}=\text{cos}{(-n \omega t)}+i\text{sin}{(-n \omega t)}=e^{-in \omega t} cosnωtisinnωt=cos(nωt)+isin(nωt)=einωt

(3)当 n = − 1 , − 2 , − 3 , … n=-1,-2,-3,\dots n=1,2,3, 时:
C n = a − n + i b − n 2 = 1 2 ( 2 T ∫ 0 T f ( t ) cos ( − n ω t ) dt + i 2 T ∫ 0 T f ( t ) sin ( − n ω t ) dt ) = 1 2 ( 2 T ∫ 0 T f ( t ) cos n ω t dt − i 2 T ∫ 0 T f ( t ) sin n ω t dt ) = 1 T ∫ 0 T f ( t ) ( cos n ω t − i sin n ω t ) dt = 1 T ∫ 0 T f ( t ) e − i n ω t dt (18) \begin{aligned} C_n&=\frac{a_{-n}+ib_{-n}}{2}\\ &=\frac{1}{2}\left(\frac{2}{T}\int_{0}^{T}f(t)\text{cos}{(-n \omega t)}\text{dt}+i\frac{2}{T}\int_{0}^{T}f(t)\text{sin}{(-n \omega t)}\text{dt}\right) \\ &=\frac{1}{2}\left(\frac{2}{T}\int_{0}^{T}f(t)\text{cos}{n \omega t}\text{dt}-i\frac{2}{T}\int_{0}^{T}f(t)\text{sin}{n \omega t}\text{dt}\right) \\ &=\frac{1}{T}\int_{0}^{T}f(t) \left(\text{cos}{n \omega t}-i\text{sin}{n \omega t}\right)\text{dt} \\ &=\frac{1}{T}\int_{0}^{T}f(t)e^{-in \omega t}\text{dt} \\ \end{aligned} \tag{18} Cn=2an+ibn=21(T20Tf(t)cos(nωt)dt+iT20Tf(t)sin(nωt)dt)=21(T20Tf(t)cosnωtdtiT20Tf(t)sinnωtdt)=T10Tf(t)(cosnωtisinnωt)dt=T10Tf(t)einωtdt(18)
观察 (16)(17)(18) 式,可以看到三个情况的最终结果竟然统一了!也就是对于 $ \forall n=0,±1,±2,±3,\dots$,有:
C n = 1 T ∫ 0 T f ( t ) e − i n ω t dt (19) \begin{aligned} C_n=\frac{1}{T}\int_{0}^{T}f(t)e^{-in \omega t}\text{dt} \end{aligned} \tag{19} Cn=T10Tf(t)einωtdt(19)
综上,对于周期函数 f ( t ) f(t) f(t),周期为 T T T,角频率为 w 0 = 2 π f 0 = 2 π T w_0=2\pi f_0=\frac{2\pi}{T} w0=2πf0=T2π,其复指数形式的傅里叶级数为:
f ( t ) = ∑ n = − ∞ ∞ C n e i n ω 0 t , C n = 1 T ∫ 0 T f ( t ) e − i n ω 0 t dt , n = 0 , ± 1 , ± 2 , ± 3 , … (20) \begin{aligned} &f(t)=\sum_{n=-\infty}^{\infty} C_n e^{in \omega_0 t}, \\ &C_n=\frac{1}{T}\int_{0}^{T}f(t)e^{-in \omega_0 t}\text{dt},n=0,±1,±2,±3,\dots \end{aligned} \tag {20} f(t)=n=Cneinω0t,Cn=T10Tf(t)einω0tdt,n=0,±1,±2,±3,(20)

0x04 傅里叶变换

傅里叶变换实际上就是对傅里叶级数的进一步扩展,即不再要求 f ( t ) f(t) f(t) 为周期函数,换句话也就是设周期 T T T 趋于 ∞ \infty

现在假设有一个周期为 T T T 的函数 f T ( t ) = f T ( t + T ) f_T(t)=f_T(t+T) fT(t)=fT(t+T),我们用复指数形式的傅里叶级数将其展开得到:
f T ( t ) = ∑ n = − ∞ ∞ C n e i n ω 0 t , C n = 1 T ∫ 0 T f T ( t ) e − i n ω 0 t dt , n = 0 , ± 1 , ± 2 , ± 3 , … (21) \begin{aligned} &f_T(t)=\sum_{n=-\infty}^{\infty} C_n e^{in \omega_0 t}, \\ &C_n=\frac{1}{T}\int_{0}^{T}f_T(t)e^{-in \omega_0 t}\text{dt},n=0,±1,±2,±3,\dots \end{aligned} \tag {21} fT(t)=n=Cneinω0t,Cn=T10TfT(t)einω0tdt,n=0,±1,±2,±3,(21)
其中 ω 0 = 2 π T \omega_0=\frac{2 \pi}{T} ω0=T2π 称为基频率。可以看到,实际上是由 C n C_n Cn 区分不同函数的,因为对于任意的傅里叶变换都是 C n C_n Cn 乘对应的复指数,再累和。假设 f T ( t ) f_T(t) fT(t)时域图像如下:
在这里插入图片描述
我们画出它在复平面上的频域图像如下:

20210108210120596

如果觉得上面这张三维复平面下的频域图不好理解,可以参考下面这张图(来自@李泽光一个相关回答里的图,图中 X n X_n Xn 与这里的 C n C_n Cn 含义类似):
在这里插入图片描述

所谓傅里叶级数就是将时域下的周期连续函数(上图),转换为频域上的非周期离散函数(下图),也就是说对于频域图像做运算: … , C − 2 e − 2 i ω 0 t + C − 1 e − i ω 0 t + C 0 e 0 + C 1 e i ω 0 t + C 2 e 2 i ω 0 t , … \dots,C_{-2}e^{-2i \omega_0t}+C_{-1}e^{-i \omega_0t}+C_{0}e^{0}+C_{1}e^{ i\omega_0t}+C_{2}e^{2i \omega_0t},\dots ,C2e2iω0t+C1eiω0t+C0e0+C1eiω0t+C2e2iω0t, 也就能得到时域下的函数图像。

而傅里叶变换要做的就是将时域下的非周期连续函数转换到频域上的非周期连续函数

img

前面说了非周期函数可以将其看作周期无限大( T → ∞ T→\infty T),也就是 lim ⁡ T → ∞ f T ( t ) = f ( t ) \lim_{T→\infty}f_T(t)=f(t) limTfT(t)=f(t)

我们令频域图像上的间隔为 Δ ω = ( n + 1 ) ω 0 − n ω 0 = ω 0 = 2 π T \Delta \omega=(n+1)\omega_0-n\omega_0=\omega_0=\frac{2\pi}{T} Δω=(n+1)ω0nω0=ω0=T2π,当 T → + ∞ T→+\infty T+ 时, Δ → − ∞ \Delta→-\infty Δ。从图像上看也就意味着每个 C i C_i Ci 之间的距离越来越来小,极限情况频域图像每个 C i C_i Ci 就挨在一起了,也就变成了连续的形式;并且离散情况下的横轴 n ω 0 n\omega_0 nω0,在连续的情况下就可以用一个连续的变量 ω \omega ω 表示。

进一步推导,将 (19) 式中的 C n C_n Cn 代入 f T ( t ) f_T(t) fT(t) ,且 T = 2 π Δ ω T=\frac{2\pi}{\Delta \omega} T=Δω2π
f ( t ) = ∑ n = − ∞ ∞ 1 T ∫ 0 T f ( t ) e − i n ω 0 t dt e i n ω 0 t = ∑ n = − ∞ ∞ Δ ω 2 π ∫ 0 T f T ( t ) e − i n ω 0 t dt e i n ω 0 t (22) \begin{aligned} f(t)&=\sum_{n=-\infty}^{\infty} \frac{1}{T}\int_{0}^{T}f(t)e^{-in \omega_0 t}\text{dt} e^{in \omega_0 t} \\ &=\sum_{n=-\infty}^{\infty} \frac{\Delta \omega}{2\pi}\int_{0}^{T}f_T(t)e^{-in \omega_0 t}\text{dt} e^{in \omega_0 t} \end{aligned} \tag {22} f(t)=n=T10Tf(t)einω0tdteinω0t=n=2πΔω0TfT(t)einω0tdteinω0t(22)
对于 (22) 式,当我们令 T → ∞ T→\infty T 时,相应的有 ∫ 0 T d t → ∫ − ∞ + ∞ \int_{0}^{T}dt→\int_{-\infty}^{+\infty} 0Tdt+ n ω 0 → ω n\omega_0→\omega nω0ω ∑ n = − ∞ ∞ Δ ω → ∫ − ∞ + ∞ d ω \sum_{n=-\infty}^{\infty}\Delta\omega→\int_{-\infty}^{+\infty}d\omega n=Δω+dω,得到:
f ( t ) = 1 2 π ∫ − ∞ ∞ ∫ − ∞ ∞ f ( t ) e − i ω t dt e i ω t d ω (23) \begin{aligned} f(t)&=\frac{1}{2\pi}\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}f(t)e^{-i \omega t}\text{dt} e^{i \omega t}\text{d}\omega \end{aligned} \tag {23} f(t)=2π1f(t)eiωtdteiωtdω(23)
最后,我们令上面 (23) 式的中间部分为为 F ( ω ) F(\omega) F(ω) 如下:
F ( w ) = ∫ − ∞ ∞ f ( t ) e − i ω t dt (24) \begin{aligned} F(w)&=\int_{-\infty}^{\infty}f(t)e^{-i \omega t}\text{dt} \end{aligned} \tag {24} F(w)=f(t)eiωtdt(24)
得到的 (24) 式即为傅里叶变换,而原来的 (23) 式即为傅里叶逆变换,可写成如下形式:
f ( t ) = 1 2 π ∫ − ∞ ∞ F ( ω ) e i ω t d ω (25) \begin{aligned} f(t)&=\frac{1}{2\pi}\int_{-\infty}^{\infty}F(\omega) e^{i \omega t}\text{d}\omega \end{aligned} \tag {25} f(t)=2π1F(ω)eiωtdω(25)

附录1:三角函数正交性证明

n ≠ m n≠m n=m时,由积化和差及牛顿莱布尼茨公式得:
∫ − π π cos n x ⋅ cos m x dx = ∫ − π π 1 2 [ cos ( n x − m x ) + cos ( n x + m x ) ] d x = 1 2 [ ∫ − π π cos ( n − m ) x dx + ∫ − π π cos ( n + m ) x dx ] = 1 2 [ 1 n − m sin ( n − m ) x ∣ − π π + 1 n + m sin ( n + m ) x ∣ − π π ] = 0 \begin{aligned} \int_{-\pi}^{\pi} \text{cos}nx·\text{cos}mx\text{dx} &= \int_{-\pi}^{\pi} \frac{1}{2}\left[\text{cos}(nx-mx) + \text{cos}(nx+mx)\right]dx \\ &=\frac{1}{2} \left[\int_{-\pi}^{\pi}\text{cos}(n-m)x\text{dx} + \int_{-\pi}^{\pi}\text{cos}(n+m)x\text{dx}\right] \\ &=\frac{1}{2} \left[\frac{1}{n-m}\text{sin}(n-m)x\bigg|_{-\pi}^{\pi}+ \frac{1}{n+m}\text{sin}(n+m)x\bigg|_{-\pi}^{\pi} \right] \\ &=0 \end{aligned} ππcosnxcosmxdx=ππ21[cos(nxmx)+cos(nx+mx)]dx=21[ππcos(nm)xdx+ππcos(n+m)xdx]=21[nm1sin(nm)xππ+n+m1sin(n+m)xππ]=0

∫ − π π cos n x ⋅ sin m x dx = ∫ − π π 1 2 [ sin ( n x + m x ) − sin ( n x − m x ) ] d x = 1 2 [ ∫ − π π sin ( n + m ) x dx + ∫ − π π sin ( n − m ) x dx ] = 0 \begin{aligned} \int_{-\pi}^{\pi} \text{cos}nx·\text{sin}mx\text{dx} &= \int_{-\pi}^{\pi} \frac{1}{2}\left[\text{sin}(nx+mx) - \text{sin}(nx-mx)\right]dx \\ &=\frac{1}{2} \left[\int_{-\pi}^{\pi}\text{sin}(n+m)x\text{dx} + \int_{-\pi}^{\pi}\text{sin}(n-m)x\text{dx}\right] \\ &=0 \end{aligned} ππcosnxsinmxdx=ππ21[sin(nx+mx)sin(nxmx)]dx=21[ππsin(n+m)xdx+ππsin(nm)xdx]=0

同理:
∫ − π π sin n x ⋅ cos m x dx = 0 ∫ − π π sin n x ⋅ sin m x dx = 0 \begin{aligned} \int_{-\pi}^{\pi} \text{sin}nx·\text{cos}mx\text{dx} &= 0 \\ \int_{-\pi}^{\pi} \text{sin}nx·\text{sin}mx\text{dx} &= 0 \end{aligned} ππsinnxcosmxdxππsinnxsinmxdx=0=0
n = m n=m n=m 时,
∫ − π π cos n x ⋅ sin m x dx = ∫ − π π sin n x ⋅ cos m x dx = ∫ − π π sin n x ⋅ cos n x dx = ∫ − π π 1 2 [ sin 2 n x + s i n 0 x ] dx = 1 2 ∫ − π π sin 2 n x dx = 0 \begin{aligned} \int_{-\pi}^{\pi} \text{cos}nx·\text{sin}mx\text{dx} &= \int_{-\pi}^{\pi} \text{sin}nx·\text{cos}mx\text{dx} \\ &=\int_{-\pi}^{\pi} \text{sin}nx·\text{cos}nx\text{dx} \\ &=\int_{-\pi}^{\pi} \frac{1}{2}\left[\text{sin}2nx+sin0x\right]\text{dx} \\ &=\frac{1}{2}\int_{-\pi}^{\pi}\text{sin}2nx\text{dx} \\ &=0 \end{aligned} ππcosnxsinmxdx=ππsinnxcosmxdx=ππsinnxcosnxdx=ππ21[sin2nx+sin0x]dx=21ππsin2nxdx=0

附录2:欧拉公式的证明

欧拉公式: e i θ = cos θ + i sin θ e^{i\theta}=\text{cos}{\theta}+i\text{sin}{\theta} eiθ=cosθ+isinθ

证明:

首先,令
f ( θ ) = e i θ cos θ + i sin θ \begin{aligned} f(\theta)=\frac{e^{i\theta}}{\text{cos}{\theta}+i\text{sin}{\theta}} \end{aligned} f(θ)=cosθ+isinθeiθ
对上式进行求导得到:
f ′ ( θ ) = i e i θ ( cos θ + i sin θ ) − e i θ ( − sin θ + i cos θ ) ( cos θ + i sin θ ) 2 = i e i θ cos θ − e i θ sin θ + e i θ sin θ − i e i θ cos θ ( cos θ + i sin θ ) 2 = 0 \begin{aligned} f'(\theta)&=\frac{ie^{i\theta}(\text{cos}{\theta}+i\text{sin}{\theta})-e^{i\theta}(-\text{sin}{\theta}+i\text{cos}{\theta})}{(\text{cos}{\theta}+i\text{sin}{\theta})^2} \\ &=\frac{ie^{i\theta}\text{cos}{\theta}-e^{i\theta}\text{sin}{\theta}+e^{i\theta}\text{sin}{\theta}-ie^{i\theta}\text{cos}{\theta}}{(\text{cos}{\theta}+i\text{sin}{\theta})^2} \\ &=0 \end{aligned} f(θ)=(cosθ+isinθ)2ieiθ(cosθ+isinθ)eiθ(sinθ+icosθ)=(cosθ+isinθ)2ieiθcosθeiθsinθ+eiθsinθieiθcosθ=0
因为 f ′ ( θ ) = 0 f'(\theta)=0 f(θ)=0,说明 f ( θ ) f(\theta) f(θ) 为一个常数,不妨设 θ = 0 \theta=0 θ=0,则:
f ( θ ) = f ( 0 ) = e i 0 cos 0 + i sin 0 = 1 \begin{aligned} f(\theta)=f(0)=\frac{e^{i0}}{\text{cos}{0}+i\text{sin}{0}}=1 \end{aligned} f(θ)=f(0)=cos0+isin0ei0=1
所以得证:
f ( θ ) = e i θ cos θ + i sin θ = 1 e i θ = cos θ + i sin θ \begin{aligned} &f(\theta)=\frac{e^{i\theta}}{\text{cos}{\theta}+i\text{sin}{\theta}}=1 \\ &e^{i\theta}=\text{cos}{\theta}+i\text{sin}{\theta} \end{aligned} f(θ)=cosθ+isinθeiθ=1eiθ=cosθ+isinθ

  • 0
    点赞
  • 6
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值