傅里叶级数和傅里叶变换

最近在做音频处理相关的内容，所以打算对傅里叶的相关内容做一个小结。大学本科时学过复变函数，不过现在已经忘的差不多了，然后毕业后也经常会用到这个东西，无论是一维的信号还是二维的图像，所以就打算对此做一个理论推导。另外，在下文推导中全部假设所有函数都满足可进行傅里叶变换的条件假设，对于其他一些特殊函数，不在本文讨论范围之内。

基本定理

三角函数的正交性

对于三角函数组合 $1,sinx,cosx,sin2x,cos2x,...,sinnx,cosnx,...$中任意两个函数的乘积在$[-\pi ,\pi ]$上的积分都为0，即$$\underset{-\pi }{\overset{\pi }{\int }}sin(x)dx=0$$ $$\underset{-\pi }{\overset{\pi }{\int }}cos(x)dx=0$$ $$\underset{-\pi }{\overset{\pi }{\int }}sin(mx)cos(nx)dx=0$$$$\underset{-\pi }{\overset{\pi }{\int }}sin(mx)sin(nx)dx=0$$ $$\underset{-\pi }{\overset{\pi }{\int }}cos(mx)cos(nx)dx=0$$

欧拉公式

$$e^{ix}=cos(x)+isin(x)$$

傅里叶级数

三角形式

对于一个周期为 $T$ 的函数，可以用一系列三角函数的叠加组成，令$\omega =\frac{2\pi }{T}$，则
$$f(t)=\frac{b_0}{2}+\sum _{n=1}^{+\infty }[a_{n}sin(n\omega t)+b_{n}cos(n\omega t)]$$
对上式在某个周期上进行积分可得$$\underset{t_0}{\overset{t_0+T}{\int }}f(t)dt=\underset{t_0}{\overset{t_0+T}{\int }}\frac{b_0}{2}dt+\underset{t_0}{\overset{t_0+T}{\int }}\sum _{n=1}^{+\infty }[a_n\sin (n\omega t)+b_n\cos (n\omega t)]dt$$$$=\frac{b_0}{2}T+\frac{2\pi }{T}\underset{x_0}{\overset{x_0+2\pi }{\int }}\sum _{n=1}^{+\infty }[a_n\sin (nx)+b_n\cos (nx)]dx=b_{0}T$$得到，$b_0=\frac{2}{T}\underset{t_0}{\overset{t_0+T}{\int }}f(t)dt$。

对上式两边同时乘以$sin(m\omega t)$，再在某个周期上进行积分可得$$\underset{t_0}{\overset{t_0+T}{\int }}f(t)\sin (m\omega t)dt=a_0\underset{t_0}{\overset{t_0+T}{\int }}\sin (m\omega t)dt+\underset{t_0}{\overset{t_0+T}{\int }}\sum _{n=1}^{+\infty }[a_n\sin (n\omega t)\sin (m\omega t)+b_n\cos (n\omega t)\sin (m\omega t)]dt$$$$=a_m\underset{t_0}{\overset{t_0+T}{\int }}{\sin }^2(m\omega t)dt=\frac{a_m}{2}\underset{t_0}{\overset{t_0+T}{\int }}[1-\cos (2m\omega t)]dt=\frac{a_m}{2}T$$
得到$a_m=\frac{2}{T}\underset{t_0}{\overset{t_0+T}{\int }}f(t)\sin (m\omega t)dt$

对上式两边同时乘以$cos(m\omega t)$，再在某个周期上进行积分可得$$\underset{t_0}{\overset{t_0+T}{\int }}f(t)\cos (m\omega t)dt=a_0\underset{t_0}{\overset{t_0+T}{\int }}\cos (m\omega t)dt+\underset{t_0}{\overset{t_0+T}{\int }}\sum _{n=1}^{+\infty }[a_n\sin (n\omega t)\cos (m\omega t)+b_n\cos (n\omega t)\cos (m\omega t)]dt$$$$=b_m\underset{t_0}{\overset{t_0+T}{\int }}{\cos }^2(m\omega t)dt=\frac{b_m}{2}\underset{t_0}{\overset{t_0+T}{\int }}[\cos (2m\omega t)+1]dt=\frac{b_m}{2}T$$
得到$b_m=\frac{2}{T}\underset{t_0}{\overset{t_0+T}{\int }}f(t)\cos (m\omega t)dt$

综上，可以得到傅里叶级数及其系数的形式，如下
$$f(t)=\frac{b_0}{2}+\sum _{n=1}^{+\infty }[a_{n}sin(n\omega t)+b_{n}cos(n\omega t)]$$$$a_n=\frac{2}{T}\underset{t_0}{\overset{t_0+T}{\int }}f(t)\sin (n\omega t)dt$$ $$b_n=\frac{2}{T}\underset{t_0}{\overset{t_0+T}{\int }}f(t)\cos (n\omega t)dt$$

指数形式

根据欧拉公式可以得到
$$cos(x)=\frac{e^{ix}+e^{-ix}}{2}$$ $$sin(x)=\frac{e^{ix}-e^{-ix}}{2i}$$ 代入傅里叶级数-三角形式中可以得到
$$f(t)=\frac{b_0}{2}+\sum _{n=1}^{+\infty }[a_n\frac{e^{in\omega t}-e^{-in\omega t}}{2i}+b_n\frac{e^{in\omega t}+e^{-in\omega t}}{2}]$$$$=\frac{b_0}{2}+\sum _{n=1}^{+\infty }[a_n\frac{e^{in\omega t}-e^{-in\omega t}}{2i}+b_n\frac{e^{in\omega t}+e^{-in\omega t}}{2}]$$$$=\frac{b_0}{2}+\sum _{n=1}^{+\infty }[\frac{b_n-i{a}_n}{2}{e}^{in\omega t}+\frac{b_n+i{a}_n}{2}{e}^{-in\omega t}]$$
令${\omega }_n=n\omega$ 则
$$f(t)=\sum _{n=-\infty }^{+\infty }{c}_n{e}^{i{\omega }_{n}t}$$当$n=0$，$$c_0=\frac{b_0}{2}=\frac{1}{T}\underset{t_0}{\overset{t_0+T}{\int }}f(t)dt$$当$n>0$, $$c_n=\frac{b_n-i{a}_n}{2}=\frac{1}{T}\underset{t_0}{\overset{t_0+T}{\int }}f(t)[cos({\omega }_{n}t)-isin({\omega }_{n}t)]dt=\frac{1}{T}\underset{t_0}{\overset{t_0+T}{\int }}f(t)e^{-i{\omega }_{n}t}dt$$ 当$n<0$, $$c_{-n}=\frac{b_{-n}+i{a}_{-n}}{2}=\frac{1}{T}\underset{t_0}{\overset{t_0+T}{\int }}f(t)[cos({\omega }_{n}t)-isin({\omega }_{n}t)]dt=\frac{1}{T}\underset{t_0}{\overset{t_0+T}{\int }}f(t)e^{-i{\omega }_{n}t}dt$$
总上所述，可以得到傅里叶级数的指数形式$$f(t)=\frac{1}{T}\sum _{n=-\infty }^{+\infty }\underset{t_0}{\overset{t_0+T}{\int }}f(t)e^{-i{\omega }_{n}t}dt{e}^{i{\omega }_{n}t}$$

傅里叶变换

傅里叶积分定理

对于周期函数有上述的傅里叶级数形式，但对于非周期函数，则无法使用傅里叶级数。但非周期函数可以认为是周期趋于无穷大的周期函数，满足傅里叶积分定理。所谓傅里叶积分定理，是指满足一定条件的函数$f(x)$ 满足
$$f(t)=\frac{1}{2\pi }\underset{-\infty }{\overset{\infty }{\int }}F(\omega )e^{i\omega x}d\omega$$$$F(\omega )=\underset{-\infty }{\overset{\infty }{\int }}f(t)e^{-i\omega t}dt$$称$F(\omega )为f(t)的傅里叶变换函数，f(t)为F(\omega )$的原函数，那么接下来，我们来证明这个傅里叶积分定理。首先，我们回顾一下黎曼积分的具体形式，设$f(t)$为$[a,b]$上的可积函数，则$$\underset{a}{\overset{b}{\int }}f(t)dt=\underset{N\to +\infty }{\lim }\sum _{i=1}^{N}f(a+\frac{i(b-a)}{N})\frac{b-a}{N}$$那么，对于非周期函数，可以用周期区域无穷大的周期函数进行近似得到，则$$f(t)=\underset{T\to \infty }{\lim }\frac{1}{T}\sum _{n=-\infty }^{+\infty }\underset{t_0}{\overset{t_0+T}{\int }}f(t)e^{-i{\omega }_{n}t}dt{e}^{i{\omega }_{n}t}$$ $$=\underset{N\to +\infty }{\lim }\frac{1}{T}\sum _{n=-\infty }^{+\infty }F({\omega }_n)e^{i{\omega }_{n}t}$$ $$=\underset{N\to +\infty }{\lim }\frac{1}{2\pi }\sum _{n=-\infty }^{+\infty }F({\omega }_n)e^{i{\omega }_{n}t}\frac{2\pi }{T}$$$$=\frac{1}{2\pi }\underset{-\infty }{\overset{+\infty }{\int }}F(\omega )e^{i\omega t}d\omega$$而$$F(\omega )=\underset{-\infty }{\overset{\infty }{\int }}f(t)e^{-i\omega t}dt$$以上就是傅里叶级数和傅里叶变换推导的全过程。