傅里叶级数与傅里叶变换

傅里叶级数与傅里叶变换

三角函数的正交性

三角函数系定义： { 0 , 1 , sin ⁡ x , cos ⁡ x , sin ⁡ 2 x , cos ⁡ 2 x , . . . , sin ⁡ n x , cos ⁡ n x , . . . } \{0,1,\sin{x},\cos{x},\sin{2x},\cos{2x},...,\sin{nx},\cos{nx},...\} {0,1,sinx,cosx,sin2x,cos2x,...,sinnx,cosnx,...}

频率不同的三角函数相乘在一个周期内($-\pi$ 到 $\pi$)的积分必定为 0：
$${\int }_{-\pi }^{\pi }\sin (mx)\sin (nx)dx=\frac{1}{2}{\int }_{-\pi }^{\pi }\cos [(m-n)x]-\cos [(m+n)x]dx=0,wheren\ne m$$

$${\int }_{-\pi }^{\pi }\cos (mx)\cos (nx)dx=\frac{1}{2}{\int }_{-\pi }^{\pi }\cos [(m-n)x]+\cos [(m+n)x]dx=0,wheren\ne m$$

$${\int }_{-\pi }^{\pi }\sin (mx)\cos (nx)dx=\frac{1}{2}{\int }_{-\pi }^{\pi }\sin [(m-n)x]+\sin [(m+n)x]dx=0$$

因此，在三角函数系中，取 2 个不同的函数，其乘积在 $[-\pi ,\pi ]$ 之间的定积分为 0。

傅里叶级数

周期为 2$\pi$ 的函数展开为傅里叶级数

将一个函数展开成傅里叶级数：
$$f(x)=\sum _{n=0}^{\infty }{a}_n\cos (nx)+\sum _{n=0}^{\infty }{b}_n\sin (nx)=\frac{a_0}{2}+\sum _{n=1}^{\infty }{a}_n\cos (nx)+\sum _{n=1}^{\infty }{b}_n\sin (nx)$$
求 $a_0$：

等式左右两边分别对 $-\pi$ 到 $\pi$ 求积分：
$${\int }_{-\pi }^{\pi }f(x)dx=2\pi a_0+{\int }_{-\pi }^{\pi }\sum _{n=1}^{\infty }{a}_n\cos (nx)dx+{\int }_{-\pi }^{\pi }\sum _{n=1}^{\infty }{b}_n\sin (nx)dx$$
根据三角函数的正交性：
$${\int }_{-\pi }^{\pi }f(x)dx=2\pi a_0\to \frac{a_0}{2}=\frac{1}{\pi }{\int }_{-\pi }^{\pi }f(x)dx$$
求 $a_n$：等式左右两边同时乘以 $\cos (mx)$ 再对 $-\pi$ 到 $\pi$ 求积分：
$${\int }_{-\pi }^{\pi }f(x)\cos (mx)dx={\int }_{-\pi }^{\pi }\frac{a_0}{2}\cos (mx)dx+{\int }_{-\pi }^{\pi }\sum _{n=1}^{\infty }{a}_n\cos (nx)\cos (mx)dx+{\int }_{-\pi }^{\pi }\sum _{n=1}^{\infty }{b}_n\sin (nx)\cos (mx)dx$$
根据三角函数的正交性：${\int }_{-\pi }^{\pi }f(x)\cos (mx)dx={\int }_{-\pi }^{\pi }{\sum }_{n=1}^{\infty }{a}_n\cos (nx)\cos (mx)dx$

仅当 $n=m$ 时：
$${\int }_{-\pi }^{\pi }f(x)\cos (mx)dx={\int }_{-\pi }^{\pi }{a}_n\cos (nx)\cos (nx)dx=a_n{\int }_{-\pi }^{\pi }{\cos }^2(nx)dx=a_n\pi$$
因此：
$$a_n=\frac{1}{\pi }{\int }_{-\pi }^{\pi }f(x)\cos (nx)dx$$
求 $b_n$：等式左右两边同时乘以 $\sin (mx)$ 再对 $-\pi$ 到 $\pi$ 求积分：
$$b_n=\frac{1}{\pi }{\int }_{-\pi }^{\pi }f(x)\sin (nx)dx$$

周期为 2$L$ 的函数展开为傅里叶级数

函数 $f(t)$ 的周期为 2$L$，利用换元的方法：
$$\begin{gathered} x=\frac{\pi }{L}t \\ t=\frac{L}{\pi } \\ f(t)=f(\frac{L}{\pi }x)\triangleq g(x) \end{gathered}$$
函数 $g(x)$ 即为周期 2$\pi$ 的函数。因此：
$$\begin{gathered} f(t)=\frac{a_0}{2}+\sum _{n=1}^{\infty }[a_n\cos (\frac{n\pi }{L}t)+b_n\sin (\frac{n\pi }{L}t)] \\ {a}_0=\frac{1}{L}{\int }_{-L}^{L}f(t)dt \\ {a}_n=\frac{1}{L}{\int }_{-L}^{L}f(t)\cos (\frac{n\pi }{L}t)dt \\ {b}_n=\frac{1}{L}{\int }_{-L}^{L}f(t)\sin (\frac{n\pi }{L}t)dt \end{gathered}$$
通常来说，$t$ 从 0 开始，周期为 $T=2L$，$\omega =\frac{\pi }{L}=\frac{2\pi }{T}$：
$$\begin{gathered} f(t)=\frac{a_0}{2}+\sum _{n=1}^{\infty }[a_n\cos (n\omega t)+b_n\sin (n\omega t)] \\ {a}_0=\frac{2}{T}{\int }_0^{T}f(t)dt \\ {a}_n=\frac{2}{T}{\int }_0^{T}f(t)\cos (n\omega t)dt \\ {b}_n=\frac{2}{T}{\int }_0^{T}f(t)\sin (n\omega t)dt \\ \text{\hspace{1em}(1)} \end{gathered}$$

傅里叶级数的复数形式

欧拉公式：$e^{i\theta }=\cos \theta +i\sin \theta$

由此引申 2 个推导：
$$\begin{gathered} \cos \theta =\frac{1}{2}(e^{i\theta }+e^{-i\theta }) \\ \sin \theta =-\frac{1}{2}i(e^{i\theta }-e^{-i\theta }) \end{gathered}$$
将其带入公式 $(1)$：
$$\begin{gathered} f(t)=\frac{a_0}{2}+\sum _{n=1}^{\infty }[a_n\frac{e^{in\omega t}+e^{-in\omega t}}{2}+b_n\frac{-i(e^{in\omega t}+e^{-in\omega t})}{2}] \\ =\frac{a_0}{2}+\sum _{n=1}^{\infty }{e}^{in\omega t}\frac{a_n-i{b}_n}{2}+\sum _{n=1}^{\infty }{e}^{-in\omega t}\frac{a_n+i{b}_n}{2} \\ =\frac{a_0}{2}+\sum _{n=1}^{\infty }{e}^{in\omega t}\frac{a_n-i{b}_n}{2}+\sum _{n=-\infty }^{-1}{e}^{in\omega t}\frac{a_{-n}+i{b}_{-n}}{2} \\ =\sum _{n=-\infty }^{\infty }{c}_n{e}^{in\omega t} \end{gathered}$$
其中，
$$\begin{gathered} c_n=\frac{a_0}{2},wheren=0 \\ {c}_n=\frac{a_n-i{b}_n}{2},wheren>0 \\ {c}_n=\frac{a_{-n}+i{b}_{-n}}{2},wheren<0 \end{gathered}$$
即：
$$\begin{gathered} c_n=\frac{a_0}{2}=\frac{1}{T}{\int }_0^{T}f(t)dt,wheren=0 \\ {c}_n=\frac{a_n-i{b}_n}{2}=\frac{1}{T}{\int }_0^{T}f(t)e^{-in\omega t}dt,wheren>0 \\ {c}_n=\frac{a_{-n}+i{b}_{-n}}{2}=\frac{1}{T}{\int }_0^{T}f(t)e^{-in\omega t}dt,wheren<0 \end{gathered}$$
因此：
$$\begin{gathered} f(t)=\sum _{n=-\infty }^{\infty }{c}_n{e}^{in\omega t} \\ where{c}_n=\frac{1}{T}{\int }_0^{T}f(t)e^{-in\omega t}dt\text{\hspace{1em}(2)} \end{gathered}$$

傅里叶变换

为了表示的方便，在这里，定义 ${\omega }_0=\frac{2\pi }{T}$，而 $\omega$ 表示连续的角频率。

上面的公式 $(2)$ 可以写成：
$$f(t)=\sum _{n=-\infty }^{\infty }\frac{{\omega }_0}{2\pi }{\int }_{-\frac{T}{2}}^{\frac{T}{2}}f(t)e^{-in{\omega }_{0}t}dt{e}^{in{\omega }_{0}t}\text{\hspace{1em}(3)}$$

当一个函数为非周期函数时：$T\to \infty$，则：
$$\begin{gathered} {\int }_{-\frac{T}{2}}^{\frac{T}{2}}dt\to {\int }_{-\infty }^{\infty }dt \\ n{\omega }_0\to \omega \\ \sum _{n=-\infty }^{\infty }{\omega }_0\to {\int }_{-\infty }^{\infty }d\omega \end{gathered}$$
因此，非周期函数的傅里叶级数可以写成：
$$f(t)=\frac{1}{2\pi }{\int }_{-\infty }^{\infty }{\int }_{-\infty }^{\infty }f(t)e^{-i\omega t}dt{e}^{i\omega t}d\omega$$
其中，红色部分就称为傅里叶变换：
$$F(\omega )={\int }_{-\infty }^{\infty }f(t)e^{-i\omega t}dt\text{\hspace{1em}(4)}$$
剩下的就是傅里叶逆变换：
$$f(t)=\frac{1}{2\pi }{\int }_{-\infty }^{\infty }F(\omega )e^{i\omega t}d\omega$$

离散傅里叶变换 (DFT)

对于一个长度为 N 的序列 $x[n]$，其可以理解成周期为 N，公式 $(3)$ 可以表示成：
$$X[k]=\sum _{n=0}^{N-1}x[n]e^{-in\frac{2\pi }{N}k}$$
其中，$k=0,1,2,...,N-1$。