从傅里叶级数到傅里叶变换(连续、离散)

傅里叶级数

任何周期函数都可以用正弦函数和余弦函数构成的无穷级数来表示:

\displaystyle f(t)=\frac{a_0}{2}+\sum_{n=1}^\infty[a_n\cos(nw_0t)+b_n\sin(nw_0t)]

由于正弦函数和余弦函数构成正交函数系,进行一些推导:

\displaystyle \int_{-\frac{\pi}{w_0}}^{\frac{\pi}{w_0}}f(t)dt=\int_{-\frac{\pi}{w_0}}^{\frac{\pi}{w_0}}\frac{a_0}{2}dt=\frac{\pi}{w_0}a_0\\{}\\ \Rightarrow a_0=\frac{w_0}{\pi}\int_{-\frac{\pi}{w_0}}^{\frac{\pi}{w_0}}f(t)dt=\frac{2}{T_0}\int_{-\frac{T_0}{2}}^{\frac{T_0}{2}}f(t)dt\\{}\\ \int_{-\frac{\pi}{w_0}}^{\frac{\pi}{w_0}}f(t)\cos(nw_0t)dt=\int_{-\frac{\pi}{w_0}}^{\frac{\pi}{w_0}}a_n\cos^2(nw_0t)dt=\frac{\pi}{w_0}a_n\\{}\\ \Rightarrow a_n=\frac{w_0}{\pi}\int_{-\frac{\pi}{w_0}}^{\frac{\pi}{w_0}}f(t)\cos(nw_0t)dt=\frac{2}{T_0}\int_{-\frac{T_0}{2}}^{\frac{T_0}{2}}f(t)\cos(nw_0t)dt\\{}\\ \int_{-\frac{\pi}{w_0}}^{\frac{\pi}{w_0}}f(t)\sin(nw_0t)dt=\int_{-\frac{\pi}{w_0}}^{\frac{\pi}{w_0}}b_n\sin^2(nw_0t)dt=\frac{\pi}{w_0}b_n\\{}\\ \Rightarrow b_n=\frac{w_0}{\pi}\int_{-\frac{\pi}{w_0}}^{\frac{\pi}{w_0}}f(t)\sin(nw_0t)dt=\frac{2}{T_0}\int_{-\frac{T_0}{2}}^{\frac{T_0}{2}}f(t)\sin(nw_0t)dt

周期函数可以通过上述推导表示成三角级数形式,然后使用欧拉公式将三角级数化为了复指数形式,于是引出傅里叶级数:
欧拉公式

\displaystyle e^{jx}=\cos(x)+j\sin(x)

使用欧拉公式得:

\displaystyle \cos(nw_0t)=\frac{e^{jnw_0t}+e^{-jnw_0t}}{2}\\{}\\ \sin(nw_0t)=\frac{e^{jnw_0t}-e^{-jnw_0t}}{2j}

代入到f(t)

\displaystyle f(t)=\frac{a_0}{2}+\sum_{n=1}^{\infty}\left(\frac{a_n-jb_n}{2}e^{jnw_0t}+\frac{a_n+jb_n}{2}e^{-jnw_0t}\right)\\{}\\ F_n=\frac{a_n-jb_n}{2}\Rightarrow F_{-n}=\frac{a_n+jb_n}{2},F_0=\frac{a_0}{2}\\{}\\ \Rightarrow f(t)=F_0+\sum_{n=1}^\infty\left(F_ne^{jnw_0t}+F_{-n}e^{-jnw_0t}\right)=\sum_{n=-\infty}^{\infty}F_ne^{hnw_0t}

a_n,b_n代到F_n中:

\displaystyle F_n=\frac{1}{T_0}\int_{-\frac{T_0}{2}}^{\frac{T_0}{2}}f(t)e^{-jnw_0t}dt

连续傅里叶变换

傅里叶级数解决的是周期函数的表达问题,其频谱图是离散的,而傅里叶变换解决的则是非周期函数的表达问题,其频谱图是连续的。
对于非周期函数,有T_0\to\infty,w_0\to dw,nw_0\to w
将式子\displaystyle F_n=\frac{1}{T_0}\int_{-\frac{T_0}{2}}^{\frac{T_0}{2}}f(t)e^{-jnw_0t}dt两边同乘T_0取极限:

\displaystyle \lim_{T_0\to\infty}F_nT_0=\lim_{T_0\to\infty}\int_{-\frac{T_0}{2}}^{\frac{T_0}{2}}f(t)e^{-jnw_0tdt}=\int_{-\infty}^\infty f(t)e^{-jwt}dt

F(w)=\lim_{T_0\to\infty}F_nT_0得到傅里叶正变换

\displaystyle F(w)=\int_{-\infty}^\infty f(t)e^{-jwt}dt

将式子F(w)=\lim_{T_0\to\infty}F_nT_0调整一下:

\displaystyle \lim_{T_0\to\infty}F_n=\lim_{T_0\to\infty}\frac{F(w)}{T_0}=\lim_{T_0\to\infty}\frac{F(w)w_0}{2\pi}=\frac{F(w)dw}{2\pi}

将式子\displaystyle f(t)=\sum_{n=-\infty}^\infty F_ne^{jnw_0t}取极限得到傅里叶逆变换

\displaystyle f(t)=\lim_{T_0\to\infty}\sum_{n=-\infty}^\infty F_ne^{jnw_0t}=\frac{1}{2\pi}\int_{-\infty}^\infty F(w)e^{jwt}dw

离散傅里叶变换

离散傅里叶变换得到的频谱也是离散的,因此其公式可以由傅里叶级数公式推导出。
w_0=\frac{2\pi}{T_0}代入到F_n

\displaystyle F_n=\frac{1}{T_0}\int_{-\frac{T_0}{2}}^{\frac{T_0}{2}}f(t)e^{-jnw_0t}dt=\frac{1}{T_0}\int_{-\frac{T_0}{2}}^{\frac{T_0}{2}}f(t)e^{-j2\pi\frac{nt}{T_0}}dt

转换成离散情况:

\displaystyle F(u)=\frac{1}{N}\sum_{x=0}^Nf(x)e^{-j2\pi\frac{ux}{N}},u=0,1,\cdots,N-1

一般离散傅里叶变换把\frac{1}{N}移到逆变换上,所以离散傅里叶正变换最终表达式为:

\displaystyle F(u)=\sum_{x=0}^Nf(x)e^{-j2\pi\frac{ux}{N}},u=0,1,\cdots,N-1

离散傅里叶逆变换:

\displaystyle f(x)=\frac{1}{N}\sum_{u=0}^NF(u)e^{j2\pi\frac{ux}{N}},x=0,1,\cdots,N-1

ai
本作品采用《CC 协议》,转载必须注明作者和本文链接

不要试图用百米冲刺的方法完成马拉松比赛。

讨论数量: 0
(= ̄ω ̄=)··· 暂无内容!

请勿发布不友善或者负能量的内容。与人为善,比聪明更重要!