傅里叶分析公式推导（最简单的傅里叶级数和傅里叶变换）

终于有机会可以把傅里叶分析推导一遍了。其实我对傅里叶变换一直停留在认识层面，今天就要好好梳理一下，为什么这么多人要用它来处理信号，它到底有什么魔力。

好，那我们就从傅里叶级数开始吧。

一、周期信号$x(t)$

傅里叶级数公式： $$x(t)=\sum _{k=-\infty }^{+\infty }{a}_k{e}^{jk{w}_{0}t}（1）$$
为什么上述公式需要用$e^{jk{w}_{0}t}$而不用其他呢？
$e^{jk{w}_{0}t}=\cos (k{w}_{0}t)+j\sin (k{w}_{0}t)$ 是由不同频率的正弦函数构成的。在自然界中，正弦是最普遍的现象，且易于表达和计算，所以首选是正弦函数。**傅里叶分析的基本思想是想将所有任意复杂的函数，都通用的表示为一组谐波相关的复指数。**对$e^{jk{w}_{0}t}$在一个周期内进行积分，可以得到
$${\int }_{T_0}{e}^{jk{w}_{0}t}dt={\int }_{T_0}(\cos (k{w}_{0}t)+j\sin (k{w}_{0}t))dt=\{\begin{array}{l}{T}_{0}k=0\\ 0k\ne 0\end{array}$$
积分结果非常的简单。
所以，在傅里叶级数中，只需要求出系数$a_k$就可以了。
对于周期函数，
$$\begin{array}{l}{\int }_{T_0}x(t)e^{-jn{w}_{0}t}dt={\int }_{T_0}{e}^{-jn{w}_{0}t}\sum _{k=-\infty }^{+\infty }{a}_k{e}^{jk{w}_{0}t}dt\\ =\sum _{k=-\infty }^{+\infty }{a}_k{\int }_{T_0}{e}^{j(k-n)w_{0}t}dt\\ =\{\begin{array}{l}{a}_k{T}_{0}k=n\\ 0k\ne n\end{array}\end{array}$$
即$\Rightarrow a_n=\frac{1}{T_0}{\int }_{T_0}x(t)e^{-jn{w}_{0}t}dt$

由于计算机无法处理无限长的函数，所以要进行截断。
现要处理非周期函数，利用上面周期函数得到的结论进行进一步的推导。

定义部分时间函数：
$x_N(t)=\sum _{k=-N}^N{a}_k{e}^{jk{w}_{0}t}$。
当$N\to \infty$时，$x_N(t)$是否和$x(t)$一致。这里引出狄利赫里条件，即
①在一个周期内，间断点数目有限，且为第一类间断点；
②在一个周期内，极大值极小值数目有限；
③在一个周期内，信号能量有限。
所以说，信号在每个连续点上都可以和$x(t)$一致。在不连续点处，会有误差，这是吉布斯效应，如下图所示。

二、非周期信号$x(t)$

以$T_0$为周期，将它复制成周期信号$\tilde{x}(t)$，$\tilde{x}(t)=x(t),|t|<\frac{T_0}{2}$。随着$T_0$趋近于无穷，周期信号会变成非周期信号，即$\tilde{x}(t)=x(t)$，两者可以相互转换，为下面的推导打下了基础。
对于周期信号$\tilde{x}(t)$来说，根据公式（1）可以得到
$$\tilde{x}(t)=\sum _{k=-\infty }^{+\infty }{a}_k{e}^{jk{w}_{0}t}（2）$$
那么，$$a_k=\frac{1}{T_0}{\int }_{-T_0/2}^{T_0/2}\tilde{x}(t)e^{-jk{w}_{0}t}dt$$
即，当$T_0\to \infty$时，$$a_k=\frac{1}{T_0}{\int }_{-\infty }^{\infty }x(t)e^{-jk{w}_{0}t}dt$$

定义傅里叶变换：$$X(w)={\int }_{-\infty }^{\infty }x(t)e^{-jwt}dt\Rightarrow X(w)=T_0{a}_k(w=k{w}_0)$$
那么，$a_k=\frac{1}{T_0}X(w)$。将该式带入式（2）得，
$$\begin{array}{l}\tilde{x}(t)=\sum _{k=-\infty }^{+\infty }{a}_k{e}^{jk{w}_{0}t}\\ =\sum _{k=-\infty }^{+\infty }\frac{1}{T_0}X(k{w}_0)e^{jk{w}_{0}t}\\ =\sum _{k=-\infty }^{+\infty }\frac{w_0}{2\pi }X(k{w}_0)e^{jk{w}_{0}t}\end{array}$$

当$T_0\to \infty \Rightarrow w_0\to 0,w_0\to dw,\sum \to \int ,\tilde{x}(t)\to x(t)$($w$变得很小，频域就变成函数了)
所以，$$x(t)=\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }X(w)e^{jwt}dw$$

综上所述，
傅里叶级数：
$$a_k=\frac{1}{T_0}{\int }_{T_0}x(t)e^{-jk{w}_{0}t}dt$$
$$x(t)=\sum _{k=-\infty }^{+\infty }{a}_k{e}^{jk{w}_{0}t}$$

傅里叶变换：
$$X(w)={\int }_{-\infty }^{\infty }x(t)e^{-jwt}dt$$
$$x(t)=\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }X(w)e^{jwt}dw$$

（PS：什么是频域？它是一种坐标系，只不过其空间中只有正弦函数。是构造出来的非真实空间。整正弦波是对频域的描述。因为任何时域波形都可以用正弦波表示。）