傅里叶变换的一些总结

傅里叶变换的一些总结

1.三角函数的正交性

三角函数系：{1(cos(0x)),sin(x),cos(x),sin(2x),cos(2x),sin(3x),cos(3x),...,sin(nx),cos(nx)}\{ 1(cos(0x)), sin(x),cos(x),sin(2x),cos(2x),sin(3x),cos(3x),...,sin(nx),cos(nx) \}{1(cos(0x)),sin(x),cos(x),sin(2x),cos(2x),sin(3x),cos(3x),...,sin(nx),cos(nx)}​​

从三角函数系中任选两个函数，他们的乘积在$[-\pi ,\pi ]$​​​上的积分具有一定的正交性。（正交，理解：两个二维向量，内积为0，图像上表示为垂直关系，例子：向量a={a1,a2,...,an}a = \{a_1,a_2,...,a_n\}a={a1​,a2​,...,an​}​​​和向量b={b1,b2,...,bn}b = \{b_1,b_2,...,b_n\}b={b1​,b2​,...,bn​}​​​的内积为0，即$a_1{b}_1+a_2{b}_2+...+a_n{b}_n=0$​​​）。
$$\begin{gathered} {\int }_{-\pi }^{\pi }sin(nx)cos(mx)=0 \\ {\int }_{-\pi }^{\pi }sin(nx)sin(mx)=0,m\ne n \\ {\int }_{-\pi }^{\pi }cos(nx)cos(mx)=0,m\ne n \end{gathered}$$
证明，利用三角函数积化和差化简。

2.周期为$2\pi$的连续函数

$f(t)=f(t+2\pi )$​，展开为傅里叶级数为：
$$\begin{gathered} f(t)=a_0+a_{1}cos(t)+b_{1}sin(t)+a_{2}cos(2t)+b_{2}sin(2t)+...+a_{\infty }cos(\infty t)+b_{\infty }sin(\infty t) \\ =a_0/2+\sum _{n=1}^{\infty }{a}_{n}cos(nt)+\sum _{n=1}^{\infty }{b}_{n}sin(nt)\text{\hspace{1em}(1)} \end{gathered}$$

$$\begin{gathered} a_0=\frac{1}{\pi }{\int }_{-\pi }^{\pi }f(t)dt \\ {a}_n=\frac{1}{\pi }{\int }_{-\pi }^{\pi }f(t)cos(nt)dt \\ {b}_n=\frac{1}{\pi }{\int }_{-\pi }^{\pi }f(t)sin(nt)dt\text{\hspace{1em}(2)} \end{gathered}$$

求解$a_0$，对等式$(1)$两边求积分${\int }_{-\pi }^{\pi }dt$，${\int }_{-\pi }^{\pi }{\sum }_{n=1}^{\infty }{b}_{n}sin(nt)dt=0$ ,${\int }_{-\pi }^{\pi }{\sum }_{n=1}^{\infty }{a}_{n}cos(nt)dt=0$

求解$a_n$​，对等式$(1)$​两边求积分${\int }_{-\pi }^{\pi }cos(kx)dt$​，化简。

求解$b_n$，对等式$(1)$两边求积分${\int }_{-\pi }^{\pi }sin(kx)dt$​，化简。

3.周期为2l的连续函数

$f(t)=f(t+2l)$，其傅里叶级数为：
$$f(t)=a_0/2+\sum _{n=1}^{\infty }{a}_{n}cos(n\omega t)+\sum _{n=1}^{\infty }{b}_{n}sin(n\omega t)\text{\hspace{1em}(3)}$$
利用周期为$2\pi$的傅里叶级数进行换元（换元法），令$x=\pi t/l,t=lx/\pi$​。​

构造函数:$f(t)=f(lx/\pi )=g(x)$​，根据(1)和(2)，$g(x)$可以分解为傅里叶级数。将x替换为t，进行化简。其中，周期$T=2l,\omega =\frac{2\pi }{T}$
$$\begin{gathered} a_0=\frac{2}{T}{\int }_{-l}^{l}f(t)dt \\ {a}_n=\frac{2}{T}{\int }_{-l}^{l}cos(n\omega t)f(t)dt \\ {b}_n=\frac{2}{T}{\int }_{-l}^{l}sin(n\omega t)f(t)dt\text{\hspace{1em}(4)} \end{gathered}$$
4.连续周期傅里叶级数的指数形式

欧拉公式：$e^{i\theta }=cos(\theta )+isin(\theta )$，

正弦函数可以使用指数表示：$cos(\theta )=\frac{e^{i\theta }+e^{-i\theta }}{2},sin(\theta )=-i(\frac{e^{i\theta }-e^{-i\theta }}{2})$​

将正弦函数带入公式（4），得：
$$\begin{gathered} f(t)=\frac{1}{T}\sum _{n=-\infty }^{\infty }{C}_n{e}^{in\omega t} \\ Cn={\int }_0^{T}f(t)e^{-n\omega t}dt\text{\hspace{1em}(5)} \end{gathered}$$
其中，$C_n$为傅里叶系数，$\omega =\frac{2\pi }{T}$​为基角频率，$T$为周期。

5.连续函数的傅里叶变换

公式(5)中，函数的周期为T，若$T->\infty$​​​​，则函数并不是周期函数，$\omega =\frac{2\pi }{T}$​,$T-\infty ,\omega -0,\omega =n_{x+1}\omega -n_x\omega =\Delta \omega$​​​​。

感受现象，若为非周期函数，函数的周期趋近于$\infty$，而$\omega 、\Delta \omega$趋近于0，离散的频谱变成了连续的频谱。

根据积分的基本定义，黎曼和，即${\int }_a^{b}f(x)dx=(f(b)-f(a))\ast \Delta x$​，对$f(t)$和$C_n$进行变换：
$$\begin{gathered} Cn={\int }_0^{T}f(t)e^{-n\omega t}dt \\ ={\int }_{-\infty }^{\infty }f(t)e^{-n\omega t} \\ f(t)=\frac{1}{T}\sum _{n=-\infty }^{\infty }{C}_n{e}^{in\omega t} \\ =\frac{\Delta \omega }{2\pi }\sum _{n=-\infty }^{\infty }{C}_n{e}^{in\omega t} \\ =\frac{1}{2\pi }{\int }_{-\infty }^{\infty }{C}_n{e}^{in\omega t}dn\omega \\ =\frac{1}{2\pi }{\int }_{-\infty }^{\infty }F({\omega }_x)e^{i{\omega }_{x}t}d{\omega }_x\text{\hspace{1em}(6)} \end{gathered}$$
其中记$n\omega ={\omega }_x,Cn=F({\omega }_x)$​。

总结，连续函数的傅里叶变换为
$$\begin{gathered} F({\omega }_x)={\int }_{-\infty }^{\infty }f(t)e^{-i{\omega }_{x}t}d(t) \\ f(t)=\frac{1}{2\pi }{\int }_{-\infty }^{\infty }F({\omega }_x)e^{i{\omega }_{x}t}d{\omega }_x\text{\hspace{1em}(7)} \end{gathered}$$
6.周期序列的离散傅里叶变换

$f_N(k)$​表示一个周期为N的离散时间序列，$f_N(k)=f_N(k+lN)$​，对于连续周期信号，可以分解为一系列角频率为$n\omega (n=0,1,-1,2,-2,...,)$​的虚指数之和$e^{i\omega nt}$​，同样周期为N的序列同样可以展开为多个虚指数之和$e^{i\omega nk}=e^{in\frac{2\pi }{N}k}$​，这些虚指数序列需要满足：$e^{in\frac{2\pi }{N}k}=e^{i(n+lN)\frac{2\pi }{N}k}$​​​​ , (理由，离散序列的频域是周期的，周期同样为N)，因此周期序列$f_N(k)$的傅里叶级数展开式可以使用有限项表示，取第一个周期$n=0,1,2,...,N-1$​，则其展开式可以写为：
$$\begin{gathered} f_N(k)=\sum _{k=0}^{N-1}{C}_n{e}^{jn\omega k} \\ {C}_n=\frac{1}{N}\sum _{k=0}^{N-1}{f}_N(k)e^{-jn\omega k}=\frac{1}{N}{F}_n(n\omega ) \\ {F}_n(n\omega )=\sum _{k=0}^{N-1}{f}_N(k)e^{-jn\omega k}\text{\hspace{1em}(8)} \end{gathered}$$
公式8中，$\omega =\frac{2\pi }{N}$​​​是一个常数，离散傅里叶系数变换以及逆变换。

7.有限长序列的离散傅里叶变换 （周期）

为了引用周期序列的知识，将有序序列$f(n)$延扩为周期序列$f_N(k)$.

对于周期序列$f_n(k)$，第一个周期$k=0,...,N-1$定义为主值区间，可以将$f(k)$看作$f_N(k)$​的主值区间序列。根据公式8，周期序列的傅里叶变换都只限于主值区间，因而将种种变换引申到相同的有限长序列，定义为有限长序列的离散傅里叶变换​。
$$\begin{gathered} F(n)=DFT[f(k)]=\sum _{k=0}^{N-1}f(k)e^{-i\frac{2\pi }{N}kn} \\ f(k)=iDFT[F(n)]=\frac{1}{N}\sum _{k=0}^{N-1}f(n)e^{i\frac{2\pi }{N}kn} \end{gathered}$$

参考：
https://www.bilibili.com/video/BV1Et411R78v/?spm_id_from=333.788.recommend_more_video.1 （连续傅里叶变换的推导）
http://www.elecfans.com/soft/49/52/2020/202002131164102.html（离散傅里叶变换的由来）