谈谈傅里叶级数和傅里叶变换

谈谈傅里叶级数和傅里叶变换

前言

​ 作为工科大学生，相信很多人都避免不了要学习傅里叶变换和傅里叶级数，我是属于每次用到都查询资料代入公式，所以对于它的理解只停留在了科普层面，没有真正理解其中奥秘，所以决心用一篇文章让自己记忆深刻。

​ 首先举个例子，从简单讲一讲这个让人又爱又恨的傅里叶级数。我们想象一个能够将污水里面各种成分都分别过滤的净水器，这里的污水就是我们平时用到的信号，而净水器则是我们的傅里叶变换。

​ 如果你生活在农村，那么下面这张图片应该更能说明问题，这个木制的谷风车就是我们的傅里叶变换或者傅里叶级数展开，而这里信号的各种成分是相对于信号的不同频率而言的。

---

一、正交函数集

​ 要学习傅里叶级数，首先就要搞清楚什么叫做正交函数集，因为我们在对周期信号进行傅里叶级数展开时，我们通常是将它变成正交函数的和式。那么我们先来复习一下向量的正交。

下面我们令两个向量的坐标分别为a（0，1，1），b（1，0，0）

那么这两个向量的点积计算过程为 0×1+1×0+1×0=0。

我们继续假设有这样两个向量相互正交

a ⃗ = { a 1 , a 2 , ⋯   , a n } \vec{a}=\lbrace a_1,a_2,\cdots,a_n\rbrace a ={a1​,a2​,⋯,an​}, b ⃗ = { b 1 , b 2 , ⋯   , b n } \vec{b}=\lbrace b_1,b_2,\cdots,b_n\rbrace b ={b1​,b2​,⋯,bn​},则有

$\sum _{n=1}^{\infty }{a}_n\cdot b_n=0$

这样一个式子看起来是不是很像两个数列，我们知道数列就是数轴上离散的点，把这些离散的点连接起来，就得到了正交函数的定义，同时也是我们在后面要用到的一条重要性质。

${\int }_{-\infty }^{+\infty }f(x)\cdot g(x)dx=0$

---

总结

通常，我们对周期函数进行傅里叶展开时我们使用的正交函数集为三角函数集，即：

{ 0 , 1 , s i n x , c o s x , s i n 2 x , c o s 2 x , ⋯   , s i n ( n x ) , c o s ( n x ) } \lbrace0,1,sinx,cosx,sin2x,cos2x,\cdots,sin(nx),cos(nx)\rbrace {0,1,sinx,cosx,sin2x,cos2x,⋯,sin(nx),cos(nx)}

至于三角函数的正交性，可以用积化和差公式进行推导，此处不做详细证明。

二、周期函数的傅里叶级数

1、周期为$2\pi$的情况

​ 按照前面所说，我们把周期函数展开成正交函数的和，而正交函数集我们采用三角函数的话，我们就可以把一个周期函数表示成为另外一种形式，如下：
$$f(t)=\sum _{n=0}^{\infty }{a}_{n}cos(nt)+\sum _{n=0}^{\infty }{b}_{n}sin(nt)$$
多数教科书是下面这种形式
$$f(t)=\frac{a_0}{2}+\sum _{n=1}^{\infty }{a}_{n}cos(nt)+\sum _{n=1}^{\infty }{b}_{n}sin(nt)$$
接下来我们来分析这两种形式的区别以及这其中的三个参数的求法。

求解$a_0$

$f(t)=a_{0}cos0t+b_{0}sin0t+\sum _{n=1}^{\infty }{a}_{n}cos(nt)+\sum _{n=1}^{\infty }{b}_{n}sin(nt)$

$f(t)=a_0+\sum _{n=1}^{\infty }{a}_{n}cosnt+\sum _{n=1}^{\infty }{b}_{n}sin(nt)$ （2-1）

我们对该式的左右两边同时积分，得

${\int }_{-\pi }^{\pi }f(t)dt={\int }_{-\pi }^{\pi }{a}_{0}dt+\sum _{n=1}^{\infty }{a}_n{\int }_{-\pi }^{\pi }cos(nt)dt+\sum _{n=1}^{\infty }{b}_n{\int }_{-\pi }^{\pi }sin(nt)dt$

${\int }_{-\pi }^{\pi }f(t)dt=2\pi a_0$

$a_0=\frac{1}{2\pi }{\int }_{-\pi }^{\pi }f(t)dt$

接下来我们令，$a_0=\frac{a_0}{2}$，这时我们就能够得到教科书上面的形式，且

$a_0=\frac{1}{\pi }{\int }_{-\pi }^{\pi }f(t)dt$ （2-2）

求解$a_n$

我们对式2-1的左右两边同时乘以$cos(nt)$,再对式子左右两边同时积分，得

${\int }_{-\pi }^{+\pi }f(t)\cdot cos(nt)dt=a_0{\int }_{-\pi }^{+\pi }cos(nt)dt+\sum _{n=1}^{\infty }{a}_n{\int }_{-\pi }^{+\pi }co{s}^2(nt)dt+\sum _{n=1}^{\infty }{b}_n{\int }_{-\pi }^{\pi }sin(nt)\cdot cos(nt)dt$

$\because a_0{\int }_{-\pi }^{\pi }cos(nt)dt=0$,$\sum _{n=1}^{\infty }{b}_n{\int }_{-\pi }^{\pi }sin(nt)\cdot cos(nt)dt=0$

$\therefore {\int }_{-\pi }^{+\pi }f(t)\cdot cos(nt)dt=\sum _{n=1}^{\infty }{a}_n{\int }_{-\pi }^{\pi }co{s}^2(nt)dt=\pi \sum _{n=1}^{\infty }{a}_n$

$\therefore a_n=\frac{1}{\pi }{\int }_{-\pi }^{\pi }f(t)\cdot cos(nt)dt$ （2-3）

求解$b_n$

我们对式2-1的左右两边同时乘以$sin(nt)$,再对式子两边积分，同理可得

$b_n=\frac{1}{\pi }{\int }_{-\pi }^{\pi }f(t)\cdot sin(nt)dt$ （2-4）

---

2、周期为T时$f(t)=f(t+T)$

​ 到了这一步，我想不到一个直接求取任意周期函数的傅里叶展开式的方法，所以我们来试试把任意周期的函数转换到周期为2Π的函数，令：$x=\frac{2\pi }{T}t$，则，
$$f(t)=f(\frac{T}{2\pi }x)$$

$$f(t+T)=f[\frac{T}{2\pi }(x+2\pi )]$$

从以上两个式子我们能够得出，当自变量$t$以$T$为周期变化时，自变量$x$则以$2\pi$为周期在变化，所以令：
$$f(\frac{T}{2\pi }x)=g(x)$$
则函数$g(x)$就是一个以$2\pi$为周期的周期函数，也就是我们熟悉的朋友了，接下来我们把这个熟悉的朋友解剖一下。
$$g(x)=\frac{a_0}{2}+\sum _{n=1}^{\infty }{a}_{n}cos(nx)+\sum _{n=1}^{\infty }{b}_{n}sin(nx)$$
其中

$a_0=\frac{1}{\pi }{\int }_{-\pi }^{\pi }g(x)dx$ （2-5）

$a_n=\frac{1}{\pi }{\int }_{-\pi }^{\pi }g(x)\cdot cos(nx)dx$ （2-6）

$b_n=\frac{1}{\pi }{\int }_{-\pi }^{\pi }g(x)\cdot sin(nx)dx$ （2-7）

为了得到周期为$T$的周期函数的这三个参数，就需要把上面三个式子中的自变量$x$转换成自变量$t$，将$x=\frac{2\pi }{T}t$分别代入式**（2-5），（2-6），（2-7）**中，可得：
$$a_0=\frac{2}{T}{\int }_{-\frac{T}{2}}^{\frac{T}{2}}f(t)dt$$

$$a_n=\frac{2}{T}{\int }_{-\frac{T}{2}}^{\frac{T}{2}}f(t)\cdot cos(\frac{2n\pi }{T}t)dt$$

$$b_n=\frac{2}{T}{\int }_{-\frac{T}{2}}^{\frac{T}{2}}f(t)\cdot sin(\frac{2n\pi }{T}t)dt$$

这样我们就得到了任意周期函数的傅里叶展开式求取方法。在工程领域，我们一般把$\frac{2\pi }{T}$记作$\omega$，并且不存在负时间轴，所以上面三式可以改写成

$$a_n=\frac{2}{T}{\int }_0^{T}f(t)\cdot cos(n\omega t)dt$$

$$b_n=\frac{2}{T}{\int }_0^{T}f(t)\cdot sin(n\omega t)dt$$

其展开式就可以写成：
$$f(t)=\frac{a_0}{2}+\sum _{n=1}^{\infty }{a}_{n}cos(n\omega t)+\sum _{n=0}^{\infty }{b}_{n}sin(n\omega t)$$

---

3、傅里叶级数的复数形式

​ 上述方法其实已经可以解决周期函数的傅里叶级数问题，为什么又要引出一个它的复数形式呢？你可能会想，这不是在给自己找麻烦嘛。然而并不是，我们用它的复数形式可以使傅里叶级数的求解变得更加简单。接下来我们就来看看这个充满巧合的复数形式的傅里叶级数。

​ 首先让我们来认识一个将复数、三角函数，$e$指数联系在一起的神奇公式（欧拉公式），我们接下来的推导过程会使用到它。至于它的证明过程，我们在此处不做赘述。
$$e^{j\theta }=cos\theta +jsin\theta$$
我们先对欧拉公式进行一次改写
$$e^{-j\theta }=cos\theta -jsin\theta$$
看到这两个式子我们是不是就能想到用它来表示三角函数，从而将复数引入到傅里叶展开式中呢，让我们来试一下

$cos\theta =\frac{1}{2}[e^{j\theta }+e^{-j\theta }]$ (2-8)

$sin\theta =-\frac{1}{2}j[e^{j\theta }-e^{-j\theta }]$ (2-9)

将式**(2-8)** **(2-9)**代入到傅里叶展开式$f(t)=\frac{a_0}{2}+\sum _{n=1}^{\infty }{a}_{n}cos(n\omega t)+\sum _{n=0}^{\infty }{b}_{n}sin(n\omega t)$，则

$f(t)=\frac{a_0}{2}+\frac{1}{2}\sum _{n=1}^{\infty }{a}_n(e^{jn\omega t}+e^{-jn\omega t})-\frac{1}{2}j\sum _{n=0}^{\infty }{b}_n(e^{jn\omega t}-e^{-jn\omega t})$

$f(t)=\frac{a_0}{2}+\sum _{n=1}^{\infty }\frac{a_n-j{b}_n}{2}{e}^{jn\omega t}+\sum _{n=1}^{\infty }\frac{a_n+j{b}_n}{2}{e}^{-jn\omega t}$

把这个式子的最后一项进行改写可得，

$f(t)=\sum _{n=0}^0\frac{a_0}{2}\cdot e^{jn\omega t}+\sum _{n=1}^{\infty }\frac{a_n-j{b}_n}{2}{e}^{jn\omega t}+\sum _{n=-\infty }^{-1}\frac{a_{-n}+j{b}_{-n}}{2}{e}^{jn\omega t}$

提出公因式，系数用$c_n$来代替，且累加的区间是整个数轴，那么

$f(t)=\sum _{-\infty }^{\infty }{c}_n\cdot e^{jn\omega t}$ (2-10)

这就是我们得傅里叶级数展开式得复数形式，接下来我们来讨论表达式中系数的求解方法。

• 当$n=0$时

  $c_0=\frac{a_0}{2}=\frac{1}{T}{\int }_0^{T}f(t)dt$

• 当$n>0$时

  $c_n=\frac{a_n-j{b}_n}{2}$

  $c_n=\frac{1}{2}[\frac{2}{T}{\int }_0^{T}f(t)\cdot cos(n\omega t)dt-j\frac{2}{T}{\int }_0^{T}f(t)\cdot sin(n\omega t)dt]$

  $c_n=\frac{1}{T}{\int }_0^{T}f(t)\cdot [cos(n\omega t)-jsin(n\omega t)]dt$

  $c_n=\frac{1}{T}{\int }_0^{T}f(t)\cdot [cos(-n\omega t)+jsin(-n\omega t)]dt$

  根据欧拉公式可得，

  $c_n=\frac{1}{T}{\int }_0^{T}f(t)\cdot e^{-jn\omega t}dt$

• 当$n<0$时

  $c_n=\frac{a_{-n}+j{b}_{-n}}{2}$

  $c_n=\frac{1}{2}[\frac{2}{T}{\int }_0^{T}f(t)\cdot cos(-n\omega t)dt+j\frac{2}{T}{\int }_0^{T}f(t)\cdot sin(-n\omega t)dt]$

  $c_n=\frac{1}{T}{\int }_0^{T}f(t)\cdot [cos(-n\omega t)+jsin(-n\omega t)]dt$

  $c_n=\frac{1}{T}{\int }_0^{T}f(t)\cdot e^{-jn\omega t}dt$ (2-11)

这几步讨论我们就可以看出复数形式傅里叶级数的巧合之处，即不论n取任何整数，我们都只需要求取一个$c_n$就能将一个周期函数分解成若干复指数函数相加的形式，并且$c_n$的求取方法都相同，式**(2-10)(2-11)**分别是它的表述形式和参数求取方法。

---

三、傅里叶变换

​ 以上的内容我们就完全地解决了周期函数的傅里叶级数展开问题，但我们的结论还是存在局限性，因为在生活中我们碰到的各种信号往往都不具有周期性，所以到这里我们还得说一句“革命尚未成功，同志仍需努力”。

​ 接下来就让我们一起思考一下怎么样解决这个非周期的问题。我们试想一下，一个非周期的函数我们是否能够把它看成一个永远也走不出它的第一个周期的周期函数呢？相信你已经有了这个解决方案。我们依然从周期函数的情况入手，它的展开式为：

$f(t)=\sum _{-\infty }^{\infty }{c}_n\cdot e^{jn{\omega }_{0}t}$ (3-1)

$c_n=\frac{1}{T}{\int }_0^{T}f(t)\cdot e^{-jn{\omega }_{0}t}dt$ (3-2)

我们先把目光集中在式**(3-2)**上面，它是一个关于$t$的定积分，所以在把这个积分结构计算出来以后，变量$t$是能够被消去的，所以对于$c_n$来说，最终能够影响它的变量就变成了$n{\omega }_0$，也即$c_n$是一个关于$n{\omega }_0$的序列，我们可以来随意作出一个它的图像

这里红色的竖线就是我们的$c_n$，至于这里使用了一个三维直角坐标系的原因就是$c_n$为一个复数，详情请参考$c_n$的求取过程。下面就是令人恍然大明白的公式推导了，我们先将**(3-2)代入(3-1)**中，可得：
$$f(t)=\sum _{-\infty }^{\infty }\frac{1}{T}{\int }_0^{T}f(t)\cdot e^{-jn{\omega }_{0}t}dt\cdot e^{jn{\omega }_{0}t}$$
根据之前的分析我们知道，一个非周期函数可以看作一个周期无穷大的周期函数，所以当上面这个用于周期函数的式子我们就可以拿来使用，但是在使用它之前，我们要对它其中的一些参数进行改写。

1. 关于$\frac{1}{T}$：

   $\because {\omega }_0=\frac{2\pi }{T}$

   $\therefore \frac{1}{T}=\frac{2\pi }{{\omega }_0}$

2. 关于${\omega }_0$：

   $\because {\omega }_0=\frac{2\pi }{T}$

   $\therefore T⟶\infty$时，${\omega }_0⟶0$

   $\therefore {\omega }_0=(n+1){\omega }_0-n{\omega }_0=\Delta \omega =d\omega$

3. $n{\omega }_0$的连续化：

   $\because {\omega }_0⟶0$

   $\therefore$图中的每一个离散点就变成了连续函数，我们就把$n{\omega }_0$写成$\omega$

4. 积分区间：

   $\because T⟶\infty$

   $\therefore$ 我们在区间$(-\infty ,+\infty )$对其积分

所以上面的式子我们就能够写成符合非周期函数的形式，即
$$f(t)=\frac{1}{2\pi }{\int }_{-\infty }^{\infty }{\int }_{-\infty }^{\infty }f(t)\cdot e^{-j\omega t}dt\cdot e^{j\omega t}d\omega$$
显然${\int }_{-\infty }^{\infty }f(t)\cdot e^{-j\omega t}dt$是一个关于$\omega$的函数，所以我们令

$F(\omega )={\int }_{-\infty }^{\infty }f(t)\cdot e^{-j\omega t}dt$ (3-3)

那么

$f(t)=\frac{1}{2\pi }{\int }_{-\infty }^{\infty }F(\omega )\cdot e^{j\omega t}d\omega$ (3-4)

式**(3-3)** **(3-4)**分别就是我们的傅立叶变换和傅里叶逆变换的定义式，至此我们就解决了非周期函数的傅里叶变换问题。