从傅里叶级数到傅里叶变换

---

傅里叶级数

为什么要有傅里叶级数

傅里叶级数（Fourier Series）是用一系列正弦波（Sinusoid）来描述任何周期函数的一种方法。图1中的三条曲线分别是周期为1秒的方波，正弦波和三角波。由于正弦和余弦只有相位差，故统称正弦波。

图1. 周期为1秒的方波，正弦波，三角波

在介绍傅里叶级数之前，让我们先来回顾一下级数的概念。级数是用一个无穷数列的加和来逼急一个数。函数项级数则是用一个函数列的加和来逼近一个函数。设$u_1(x), u_2(x), u_3(x), \cdots, u_n(x), \cdots$为定义在$(a, b)$内的函数序列。则

$$\sum_{n=1}^\infty u(x)=u_1(x)+u_2(x)+u_3(x)+\cdots+u_n(x)+\cdots$$
称为定义在 $(a, b)$内的函数项级数。为什么要把一个看似简单的函数分解成一大堆函数的和呢？因为有些函数直接研究起来比较困难，以某种形式的级数进行展开，对里面的每一项单独研究，会变得更简单，也使得计算更加容易。级数有千千万万种，如泰勒级数，等比级数，调和级数等等。但是有一种由正弦函数组合而成的级数，显得尤为重要。这就是傅里叶级数。为什么傅里叶技术格外重要呢？这要归功于正弦函数优秀的性质。我们将函数展开成级数是为了获得更加简便和易于计算的形式。而当正弦波输入一个系统时，输出仍然是一个正弦波，只有振幅、相位和频率会发生变化，而不像其他的级数会使函数形式本身发生改变。这使得傅里叶级数在分析函数时具有了巨大的优势。此外，由于通信系统中电磁场与电磁波，以及诸多物理原理都与正弦信号有关，所以造就了傅里叶级数如此重要的地位。

傅里叶级数是怎么来的

傅里叶级数的得出

假如有两个周期函数（Periodic Function），它们的频率分别为$f_1$和$f_2$，那么他们的叠加还是一个周期函数吗？频率又是多少呢？显然，两个不同频率的周期函数叠加仍然是一个周期函数，叠加后函数的周期是两个原函数周期的最小公倍数。因此，当一组频率为$1Hz, 2Hz, 3Hz, \cdots, nHz, \cdots$的周期函数叠加时，叠加后的函数频率必然为$1Hz$。然而，如果采用了诸如$1.1Hz, 2.5Hz, 3,12435Hz$之类频率的级数项，则输出频率将陷入混乱，所以这里只选取如$1Hz, 2Hz, 3Hz, \cdots, nHz, \cdots$的频率作为级数项。$1Hz$可以作为基本频率，改写作$fHz$， 则级数项将变为$fHz, 2fHz, 3fHz, \cdots, nfHz, \cdots$。回想图1中周期为1秒的方波函数，我们可以将它表示成

$$f(t)=\sum_{n=0}^\infty b_n \cdot \sin(2\pi\cdot n\cdot t) \quad n\in N$$
然而，上面我们所表示的函数恰好是一个周期为1秒的奇函数。如果用上面的公式来逼近一个偶函数则无法实现。所以，若 $f(t)$是一个周期为1秒的偶函数，则
$$f(t)=\sum_{n=0}^\infty a_n \cdot \cos(2\pi\cdot n\cdot t) \quad n\in N$$
因此，当 $f(t)$是一个周期为 $T$，频率为 $f$的一般函数，既有奇函数成分也有偶函数成分，此外，作为奇函数或偶函数对称点可能相对原点产生位移，易知这个位移不会影响函数的形状，可以用一个常数来表示，为了后续计算方便，这个位移记作 $\frac{a_0}{2}$，则有
$$f(t)=\frac{a_0}{2} + \sum_{n=0}^\infty \left[a_n \cos(2 \pi f n t)＋b_n \sin(2 \pi f n t)\right] \quad n\in N$$
至此，我们已经得到了傅里叶级数的完整表达形式。

傅里叶级数中参数的确定与函数的正交性

那么如何确定上面公式中的$b_n$呢？在这之前，然我们来谈谈什么是函数的正交性。学过线性代数的同学都知道，两个向量的正交是通过内积为零来定义的。而内积则是将向量的对应项相乘再求和来得到的。假设我们有一个任意长度的向量，每两个元素之间的距离无限小，那么我们就可以把这样两个向量看作两个连续的函数。类比内积的概念，两个函数正交也就是将两个函数赋予相同的自变量，再相乘，再做积分，如果积分等于零，则说明这两个函数在积分域上是正交的。

我们高兴的发现，不同频率的三角函数具有如下的正交性。其中$n, m \in N$

∫T2−T2sin(2πfnt)sin(2πfmt)dt=⎧⎩⎨⎪⎪mcos(πTfm)sin(πTfn)−ncos(πTfn)sin(πTfm)πf(m2−n2)=0T2−sin(2πTfn)4πfn=T2n≠mn=m∫T2−T2cos(2πfnt)cos(2πfmt)dt=⎧⎩⎨⎪⎪mcos(πTfn)sin(πTfm)−ncos(πTfm)sin(πTfn)πf(m2−n2)=0T2