傅里叶变换一步步详细推导

前言

在本地Markdown下写的,传上来后公式不能对齐😠​
在大学的时候接触过几次傅里叶变换的知识,但是从来没真正懂过.也在网上找过很多资料,看过很多视频,但是,这些内容要么举些简单的例子说说直观上的理解,要么就是直接堆出公式没有任何推导.

直到一个巧合在B站上看到这样一个视频才真正搞懂,非常感谢这位UP主DR_CAN.这篇博客也主要是对视频中的推导模仿一遍.同时记录下笔记方便复习.

另外,记录当时一条印象很深的弹幕: 根本就没有人会学不懂数学,只是缺乏推导过程.

详细过程

1. 推导中的重要数学公式

三角函数系： { 1 , s i n ( x ) , c o s ( x ) , s i n ( 2 x ) , c o s ( 2 x ) , s i n ( 3 x ) , c o s ( 3 x ) , . . . , s i n ( n x ) , c o s ( n x ) } \{1,sin(x),cos(x),sin(2x),cos(2x),sin(3x),cos(3x),...,sin(nx),cos(nx)\} {1,sin(x),cos(x),sin(2x),cos(2x),sin(3x),cos(3x),...,sin(nx),cos(nx)}

在这个三角函数系里有个重要的性质,也就是傅里叶变换的推导精髓所在：

在三角函数系里任取两个不同的函数$f(x),g(x)$(1也可以看作一个函数).都有一下公式成立:
$${\int }_{-\pi }^{\pi }f(x)g(x)dx=0$$
比如说:${\int }_{-\pi }^{\pi }sin(4x)cos(7x)dx=0$,${\int }_{-\pi }^{\pi }1\cdot sin(9x)dx=0$,${\int }_{-\pi }^{\pi }sin(x)cos(x)dx=0$ , … …

如果取得$f(x),g(x)$两个函数是相同的,那么
$${\int }_{-\pi }^{\pi }1\cdot 1dx=2\pi ,{\int }_{-\pi }^{\pi }sin(nx)sin(nx)dx=\pi ,{\int }_{-\pi }^{\pi }cos(nx)cos(nx)dx=\pi (n=1,2,3...)$$
(这就不证明了,积分算一下就可以得到)

结论就是：任取三角函数系的不同函数相乘在一个周期的积分都为0.只有相同的函数才不为0.

欧拉公式：
$$e^{i\theta }=cos(\theta )+isin(\theta )$$
简单地理解,欧拉公式就是复数的另一种表示形式.$e^{i\theta }$就看作复数.

2. 傅里叶级数

先讨论周期函数$f(t)$,周期为$T$. 也就是$f(t)=f(t+T)$.

傅里叶级数本质就是用这无穷多个三角函数来表示或者无限近似函数$f(t)$.公式如下
$$f(t)=a_0+\sum _{n=1}^{\infty }[a_{n}cos(nt)+b_{n}sin(nt)]$$
在上面的三角函数系里除了$1,sin(x),cos(x)$,其他函数其实都是在$sin(x),cos(x)$的乘上一个系数.这里$sin(x),cos(x)$的$T=\frac{2\pi }{1}$.我们令$w=\frac{2\pi }{T}$. 称$w$为基频率.在上面的三角函数系里基频率就是1.如果把最开始的$sin(x),cos(x)$换成$sin(1.5x),cos(1.5x)$.那么整个三角函数系就变了啊.基频率也就变了.所以我们还需要一个变量$w$来描述三角函数系长什么样子.上面公式没有$w$只是因为碰巧$w=1$而已.所以傅里叶级数完整公式如下.(太啰嗦了… …)
$$f(t)=a_0+\sum _{n=1}^{\infty }[a_{n}cos(nwt)+b_{n}sin(nwt)]$$
写到这里其实就差不多了,但是欧拉说这里我必须要插一脚.这个公式太难看了.

根据欧拉公式有
$$\begin{gathered} e^{inwt}=cos(nwt)+isin(nwt) \\ {e}^{-inwt}=cos(nwt)-isin(nwt) \end{gathered}$$
那么
$$cos(nwt)=\frac{e^{inwt}+e^{-inwt}}{2},sin(nwt)=\frac{e^{inwt}-e^{-inwt}}{2i}$$
带入上式
$$\begin{gathered} f(t)=a_0+\sum _{n=1}^{\infty }[a_n\frac{e^{inwt}+e^{-inwt}}{2}+b_n\frac{e^{inwt}-e^{-inwt}}{2i}] \\ f(t)=a_0+\sum _{n=1}^{\infty }(\frac{a_n-i{b}_n}{2})e^{inwt}+\sum _{n=1}^{\infty }(\frac{a_n+i{b}_n}{2})e^{-inwt} \end{gathered}$$
还不够好看啊,好像也不能再化简了吧.(但是接下来这个操作很帅,注意细节)
$$\begin{gathered} f(t)=\sum _{n=0}^0{a}_n{e}^{inwt}+\sum _{n=1}^{\infty }(\frac{a_n-i{b}_n}{2})e^{inwt}+\sum _{n=-1}^{-\infty }(\frac{a_{-n}+i{b}_{-n}}{2})e^{inwt} \\ f(t)=\sum _{n=-\infty }^{\infty }{c}_n{e}^{inwt} \end{gathered}$$
其中
$$\begin{gathered} 当n=0时,c_n=a_0 \\ 当n=1,2,3...时,c_n=\frac{a_n-i{b}_n}{2} \\ 当n=-1,-2,-3...时,c_n=\frac{a_{-n}+i{b}_{-n}}{2} \end{gathered}$$

虽然很简洁,但是本质上都是一样的,只是公式的恒等变形.🤒

3. 周期函数傅里叶变换

分析上面的公式,$t$是变量,其实sin,cos里面都是固定的.不知道的只有$a_0,a_n,b_n$.如果知道了$a_0,a_n,b_n$,那函数$f(t)$就写成了三角函数之和的形式了.下面就来求$a_0,a_n,b_n$.

先求$a_0$,对上式左右两边同时乘上$1$,也就是同时乘上$cos(0t)$,再从$0$到$T$积分
$$\begin{gathered} {\int }_0^{T}f(t)dt={\int }_0^T{a}_{0}dt+{\int }_0^T\sum _{n=1}^{\infty }[a_{n}cos(nwt)+b_{n}sin(nwt)]dt \\ {\int }_0^{T}f(t)cos(0t)dt={\int }_0^T{a}_{0}cos(0t)dt+{\int }_0^T\sum _{n=1}^{\infty }[a_{n}cos(nwt)cos(0t)+b_{n}sin(nwt)cos(0t)]dt \end{gathered}$$
写到这里应该就可以知道三角函数系的结论能干什么了吧😃

上式最右边那项里面$cos(nwt)cos(0t)$,对它积分不就等于0嘛.$sin(nwt)cos(0t)$,对它积分也等于0.因为他们在三角函数系是不同的函数啊.因此
$$\begin{gathered} {\int }_0^{T}f(t)dt={\int }_0^T{a}_{0}cos(0t)dt+0 \\ {\int }_0^{T}f(t)dt={\int }_0^T{a}_{0}dt=T\cdot a_0 \\ {a}_0=\frac{1}{T}{\int }_0^{T}f(t)dt \end{gathered}$$
接着算$a_n$,这里$n=1,2,3...$ 对傅里叶级数公式左右两边同时乘上$cos(mwt)$,再从$0$到$T$积分
$$\begin{gathered} f(t)=a_0+\sum _{n=1}^{\infty }[a_{n}cos(nwt)+b_{n}sin(nwt)] \\ {\int }_0^{T}f(t)cos(mwt)dt={\int }_0^T{a}_{0}cos(mwt)dt+{\int }_0^{T}cos(mwt)\sum _{n=1}^{\infty }[a_{n}cos(nwt)+b_{n}sin(nwt)]dt \\ {\int }_0^{T}f(t)cos(mwt)dt={\int }_0^T{a}_{0}cos(0t)cos(mwt)dt+{\int }_0^{T}cos(mwt)\sum _{n=1}^{\infty }[a_{n}cos(nwt)+b_{n}sin(nwt)]dt \end{gathered}$$
看到了吧,三角函数系的结论又来了,对$cos(0t)cos(mwt),sin(nwt)cos(mwt)$的积分又都是0了.只有$cos(nwt)cos(nwt)$的积分才不等于0,对于其他$m\ne n,cos(mwt)cos(nwt)$的积分也是0,因此求和符号也就没了
$$\begin{gathered} {\int }_0^{T}f(t)cos(mwt)dt=0+{\int }_0^T\sum _{n=1}^{\infty }{a}_{n}cos(nwt)cos(mwt)dt+{\int }_0^T\sum _{n=1}^{\infty }{b}_{n}sin(nwt)cos(mwt)dt \\ {\int }_0^{T}f(t)cos(mwt)dt={\int }_0^T{a}_{m}cos(mwt)cos(mwt)dt+0 \\ {\int }_0^{T}f(t)cos(mwt)dt={\int }_0^T{a}_{m}co{s}^2(mwt)dt=a_m{\int }_0^T\frac{1+cos(2mwt)}{2}dt=\frac{T}{2}{a}_m \\ {a}_m=\frac{2}{T}{\int }_0^{T}f(t)cos(mwt)dt(m=1,2,3...) \\ {a}_n=\frac{2}{T}{\int }_0^{T}f(t)cos(nwt)dt(n=1,2,3...) \end{gathered}$$
求得了$a_n$,来算$b_n$. 其实很明显能观察到的相似性,我们根据$a_n$的规律直接得出$b_n$
$$b_n=\frac{2}{T}{\int }_0^{T}f(t)sin(nwt)dt(n=1,2,3...)$$

现在把$a_0,a_n,b_n$,带入上面求$c_n$的公式

当$n=0$时
$$c_n=a_0=\frac{1}{T}{\int }_0^{T}f(t)dt=\frac{1}{T}{\int }_0^{T}f(t)e^{-iw\cdot 0}dt$$
当$n=1,2,3...$时
$$\begin{gathered} c_n=\frac{a_n-i{b}_n}{2}=\frac{\frac{2}{T}{\int }_0^{T}f(t)cos(nwt)dt-i\cdot \frac{2}{T}{\int }_0^{T}f(t)sin(nwt)dt}{2} \\ =\frac{1}{T}{\int }_0^{T}f(t)[cos(nwt)-isin(nwt)]dt \\ =\frac{1}{T}{\int }_0^{T}f(t)e^{-inwt}dt \end{gathered}$$
当$n=-1,-2,-3...$时
$$\begin{gathered} c_n=\frac{a_{-n}+i{b}_{-n}}{2}=\frac{\frac{2}{T}{\int }_0^{T}f(t)cos(-nwt)dt+i\cdot \frac{2}{T}{\int }_0^{T}f(t)sin(-nwt)dt}{2} \\ =\frac{1}{T}{\int }_0^{T}f(t)[cos(nwt)-isin(nwt)]dt \\ =\frac{1}{T}{\int }_0^{T}f(t)e^{-inwt}dt \end{gathered}$$
我们惊奇地发现对于任意整数$n$,$c_n$的表达式都是一样的.😂

总结上面的公式
$$f(t)=\sum _{n=-\infty }^{\infty }{c}_n{e}^{inwt}其中,c_n=\frac{1}{T}{\int }_0^{T}f(t)e^{-inwt}dt$$
其中$c_n$是可以通过上面公式算出来的复数.而$e$的指数$inwt$,对于每一个不同的函数$f(t)$其实都是一样的.那么决定函数的不同其实就是$c_n$.我们知道一系列这样的$(t,f(t))$对就可以画出函数图像,也就是说这样的$(t,f(t))$能够唯一确定函数.但是现在我们的函数由这样的$(n,c_n)$对来确定了.其实也就是说函数表示发生了转化–从时域到频域的转化.函数值是不变的,它们等号连接的嘛.只是表诉的方式改变了.如下图

(上图右边纵轴$c_n$其实是个复数,可以理解为应该有两个维度,一个实部,一个虚部,但是这里为了简单画图,就把它画成了实数.但其实它是个复数.)

4. 任意函数傅里叶变换

对于一个非周期函数,我们可以假设它是个周期函数,它的周期$T\to +\infty$.就可以用上面的公式求出了.

$w$是基频率.$w=\frac{2\pi }{T}$.当$T\to +\infty$时,$w\to 0$.$w=(n+1)w-nw=\Delta w$.所以$\Delta w\to 0$.

从上图可以看到,随着$\Delta w\to 0$,$nw$会从一个离散的量变为一个连续的变量.

对于上小节中$c_n$可以写作
$$c_n=\frac{1}{T}{\int }_{-\frac{T}{2}}^{\frac{T}{2}}f(t)e^{-inwt}dt$$
那么,对于任意函数的变换为
$$\begin{gathered} f(t)=\underset{T\to \infty }{\lim }\sum _{n=-\infty }^{\infty }\frac{1}{T}{\int }_{-\frac{T}{2}}^{\frac{T}{2}}f(t)e^{-inwt}dt{e}^{inwt} \\ =\underset{\Delta w\to 0}{\lim }\sum _{n=-\infty }^{\infty }\frac{\Delta w}{2\pi }{\int }_{-\infty }^{+\infty }f(t)e^{-inwt}dt{e}^{inwt} \end{gathered}$$
当$\Delta w\to 0$时,$(n+1)w$和$nw$就几乎相等.令$W=nw$,($w$是一经确定就不变的量,$n$才是离散的变量),那么就可以把$W$看作连续变化的量.所以上式求和就可写作求积分
$$\begin{gathered} f(t)={\int }_{-\infty }^{+\infty }\frac{1}{2\pi }({\int }_{-\infty }^{+\infty }f(t)e^{-iWt}dt)e^{iWt}dW \\ =\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }({\int }_{-\infty }^{+\infty }f(t)e^{-iWt}dt)e^{iWt}dW \end{gathered}$$
其中,${\int }_{-\infty }^{+\infty }f(t)e^{-iWt}dt$的积分变量是$t$,也就是说它是关于$W$的函数.令
$$F(W)={\int }_{-\infty }^{+\infty }f(t)e^{-iWt}dt$$
那么函数$F(W)$就是$f(t)$的傅里叶变换.把$F(W)$带入$f(t)$得
$$f(t)=\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }F(W)e^{iWt}dW$$
$f(t)$就是傅里叶变换的逆变换.

函数.令
$$F(W)={\int }_{-\infty }^{+\infty }f(t)e^{-iWt}dt$$
那么函数$F(W)$就是$f(t)$的傅里叶变换.把$F(W)$带入$f(t)$得
$$f(t)=\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }F(W)e^{iWt}dW$$
$f(t)$就是傅里叶变换的逆变换.