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

定义1：傅里叶变换定义式

$$F(\omega )={\int }_{-\infty }^{+\infty }f(t)e^{-j\omega t}dt\text{\hspace{1em}(1)}$$

定义2：傅里叶级数
对于再$[-T/2,T/2]$上，周期为$T$的函数$f_T(t)$，有
$$f_T(t)\sim \frac{a_0}{2}+\sum _{n=1}^{+\infty }[a_{n}cos(n\omega t)+b_{n}sin(n\omega t)]\text{\hspace{1em}(2)}$$
其中 $$\omega =2\pi /T$$
$$a_n=\frac{2}{T}{\int }_{-T/2}^{T/2}{f}_T(t)cos(n\omega t)dt$$
$$b_n=\frac{2}{T}{\int }_{-T/2}^{T/2}{f}_T(t)sin(n\omega t)dt$$
显然，
$$a_0=\frac{2}{T}{\int }_{-T/2}^{T/2}{f}_T(t)dt$$

这些东西都是怎么来的呢？ 下面我们分解一下

引用：https://blog.csdn.net/constantin_ouc/article/details/78881709

1.前提： 我们基于一个认识，即函数$f(t)$可以分解为若干个余弦函数的加权和的形式。即，
$$f(t)=A_{0}cos(k_{0}t+{\psi }_0)+A_{1}cos(k_{1}t+{\psi }_1)+\cdots +A_{n}cos(k_{n}t+{\psi }_n)\text{\hspace{1em}(1.1)}$$
其中，$A_i$为余弦波的振幅，${\psi }_i$为对应余弦分量$A_{1}cos(k_{1}t+{\psi }_1)$的相位角，$i=0,1,\dots$

那么，对于数轴$t$上不同的值，$t_1,t_2,\dots ,t_m$有
$$\{\begin{array}{l}f(t_1)=A_{0}cos(k_0{t}_1+{\psi }_0)+A_{1}cos(k_1{t}_1+{\psi }_1)+\cdots +A_{n}cos(k_n{t}_1+{\psi }_n)\\ f(t_2)=A_{0}cos(k_0{t}_2+{\psi }_0)+A_{1}cos(k_1{t}_2+{\psi }_1)+\cdots +A_{n}cos(k_n{t}_2+{\psi }_n)\\ \dots \\ f(t_m)=A_{0}cos(k_0{t}_m+{\psi }_0)+A_{1}cos(k_1{t}_m+{\psi }_1)+\cdots +A_{n}cos(k_n{t}_m+{\psi }_n)\end{array}\text{\hspace{1em}(1.2)}$$

令$$A=\left[\begin{array}{ccccccc}cos(k_0{t}_1+{\psi }_0)& & cos(k_1{t}_1+{\psi }_1)& & \dots & & cos(k_n{t}_1+{\psi }_n)\\ cos(k_0{t}_2+{\psi }_0)& & cos(k_1{t}_2+{\psi }_1)& & \dots & & cos(k_n{t}_1+{\psi }_n)\\ \dots & & \dots & & \dots & & \dots \\ cos(k_0{t}_m+{\psi }_0)& & cos(k_1{t}_m+{\psi }_1)& & \dots & & cos(k_n{t}_m+{\psi }_n)\end{array}\right]$$,
$$X={\left[\begin{array}{ccccccc}{a}_0& & a_1& & \dots & & a_n\end{array}\right]}^{\mathrm{T}}$$,
$$B={\left[\begin{array}{ccccccc}f(t_1)& & f(t_2)& & \dots & & f(t_m)\end{array}\right]}^{\mathrm{T}}$$
得到$$B=AX$$
即$$X=A^{-1}B$$可以得到系数向量$X$。

在(1.1)中，$f(t)=A_{0}cos(k_{0}t+{\psi }_0)+A_{1}cos(k_{1}t+{\psi }_1)+\cdots +A_{n}cos(k_{n}t+{\psi }_n)$的角频率 { k i } \{k_i\} {ki​}为随意选取的值。显然，不同的角频率 { k i } \{k_i\} {ki​}可以得到不同的余弦分量系数 { a i } \{a_i\} {ai​}。

为了便于处理，我们规定角频率 { k i } \{k_i\} {ki​}的选取。令$$k_i=i,(i=0,1,\dots )$$

我们看到$f(t)={\sum }_{i=0}^n{A}_{i}cos(k_{i}t+{\psi }_i)$中，相位角${\psi }_i$ 太讨厌了！ 不利于分析！我们把相位角处理掉！
$$\begin{gathered} f(t)=\sum _{i=0}^n{A}_{i}cos(k_{i}t+{\psi }_i) \\ =\sum _{i=0}^n{A}_i(cos(k_{i}t)cos({\psi }_i)-sin(k_{i}t)sin({\psi }_i)) \\ =\sum _{i=0}^n{A}_{i}cos({\psi }_i)cos(k_{i}t)-A_{i}sin({\psi }_i)sin(k_{i}t) \end{gathered}$$
令$a_i=-A_{i}sin({\psi }_i),b_i=A_{i}cos({\psi }_i)$，则上式有

$$f(t)=\sum _{i=0}^n[a_{i}sin(k_{i}t)+b_{i}cos(k_{i}t)]\text{\hspace{1em}(1.3)}$$
以上为傅里叶级数的展开形式。
\\ \\
这样做有什么好处呢？
答：
好处1(周期性)：对于(1.3)中$f(t)={\sum }_{i=0}^n[a_{i}sin(k_{i}t)+b_{i}cos(k_{i}t)]$而言，记$T_i$为期周期，则
$$T_i=\frac{2\pi }{k_i}$$
由于$k_i=i,(i=0,1,\dots )$，使得 { T i } \{T_i\} {Ti​}有公共的周期 $T=2\pi$，这使得$f(t)$成为周期为$2\pi$的函数。

好处2：正弦函数、余弦函数在一个周期内的积分等于0
$${\int }_{-\frac{T}{2}}^{\frac{T}{2}}{a}_{i}cos(k_{i}t)dt={\int }_0^T{a}_{i}cos(k_{i}t)dt=0\text{\hspace{1em}(1.4)}$$

\\
好处说完了，但是(1.3)中$f(t)={\sum }_{i=0}^n[a_{i}sin(k_{i}t)+b_{i}cos(k_{i}t)]$中的系数$a_i,b_i$怎么求呢？
鉴于$f(t)$有周期$T=2\pi$，我们设计一个积分形式

(1.5) ∫ − π π f ( t ) c o s ( k i t ) d t = ∫ − π π ∑ j = 0 n [ a j s i n ( k j t ) + b j c o s ( k j t ) ] c o s ( k i t ) d t = ∫ − π π ∑ j = 0 n [ a j s i n ( k j t ) c o s ( k i t ) + b j c o s ( k j t ) c o s ( k i t ) ] d t = ∫ − π π { b 0 c o s ( k i t ) + ∑ j = 1 n [ a j s i n ( k j t ) c o s ( k i t ) + b j c o s ( k j t ) c o s ( k i t ) ] } d t = 2 π b 0 c o s ( k i t ) + ∫ − π π ∑ j = 1 n [ a j s i n ( k j t ) c o s ( k i t ) ] d t + ∫ − π π ∑ j = 1 n [ b j c o s ( k j t ) c o s ( k i t ) ] d t \int_{-\pi}^{\pi}f(t)cos(k_it)dt=\int_{-\pi}^{\pi}\sum_{j=0}^{n}[a_jsin(k_jt)+b_jcos(k_jt)]cos(k_it)dt\\ =\int_{-\pi}^{\pi}\sum_{j=0}^{n}[a_jsin(k_jt)cos(k_it)+b_jcos(k_jt)cos(k_it)]dt\\ =\int_{-\pi}^{\pi} \{b_0cos(k_it)+\sum_{j=1}^{n}[a_jsin(k_jt)cos(k_it)+b_jcos(k_jt)cos(k_it)]\}dt\\ =2\pi b_0cos(k_it)+\int_{-\pi}^{\pi}\sum_{j=1}^{n}[a_jsin(k_jt)cos(k_it)]dt+\int_{-\pi}^{\pi}\sum_{j=1}^{n}[b_jcos(k_jt)cos(k_it)]dt \tag{1.5}\\ ∫−ππ​f(t)cos(ki​t)dt=∫−ππ​j=0∑n​[aj​sin(kj​t)+bj​cos(kj​t)]cos(ki​t)dt=∫−ππ​j=0∑n​[aj​sin(kj​t)cos(ki​t)+bj​cos(kj​t)cos(ki​t)]dt=∫−ππ​{b0​cos(ki​t)+j=1∑n​[aj​sin(kj​t)cos(ki​t)+bj​cos(kj​t)cos(ki​t)]}dt=2πb0​cos(ki​t)+∫−ππ​j=1∑n​[aj​sin(kj​t)cos(ki​t)]dt+∫−ππ​j=1∑n​[bj​cos(kj​t)cos(ki​t)]dt(1.5)
上式中第一项，正、余弦函数在一个周期内的积分为0，故第一项为0；第二项为奇函数积分，为0；第三项，根据三角函数的正交性，当$k_j̸=k_i$时，为0；仅当$k_j=k_i$时不为0。

故，第三项为
$$\begin{gathered} {\int }_{-\pi }^{\pi }\sum _{j=1}^n[b_{j}cos(k_{j}t)cos(k_{i}t)]dt \\ ={\int }_{-\pi }^{\pi }[b_{i}co{s}^2(k_{i}t)]dt \\ =\frac{b_i}{2}{\int }_{-\pi }^{\pi }(1-cos(2k_{i}t))dt（半角公式） \\ =\pi b_i-0=\pi b_i\text{\hspace{1em}(1.6)} \end{gathered}$$

结合(1.5)、(1.6)，${\int }_{-\pi }^{\pi }f(t)cos(k_{i}t)dt=\pi b_i$，可得
$$b_i=\frac{1}{\pi }{\int }_{-\pi }^{\pi }f(t)cos(k_{i}t)dt\text{\hspace{1em}(1.7)}$$

同理，与(1.5)一致，取${\int }_{-\pi }^{\pi }f(t)sin(k_{i}t)dt$，可得
$$a_i=\frac{1}{\pi }{\int }_{-\pi }^{\pi }f(t)sin(k_{i}t)dt\text{\hspace{1em}(1.8)}$$

以上就是傅里叶级数的由来

下面我们要推出傅里叶变换。

$$傅里叶级数\stackrel{e^{ix}=cosx+isinx}{⟶}傅里叶变换$$

设$j^2=-1$。
有
$$\begin{gathered} cos\theta =\frac{e^{j\theta }+e^{-j\theta }}{2} \\ sin\theta =\frac{e^{j\theta }-e^{-j\theta }}{2}=-j\frac{e^{j\theta }-e^{-j\theta }}{2}\text{\hspace{1em}(1.9)} \end{gathered}$$
（1.3）写成
$$\begin{gathered} f_T(x)=\frac{a_0}{2}+\sum _{n=1}^{\infty }(a_{n}sinn\omega x+b_{n}cosn\omega x) \\ =\frac{a_0}{2}+\sum _{n=1}^{\infty }(-j{a}_n\frac{e^{jn\omega x}-e^{-jn\omega x}}{2}+b_n\frac{e^{jn\omega x}+e^{-jn\omega x}}{2}) \\ =\frac{a_0}{2}+\sum _{n=1}^{\infty }(\frac{b_n-j{a}_n}{2}{e}^{jn\omega x}+\frac{b_n-j{a}_n}{2}{e}^{-jn\omega x}) \end{gathered}$$

令
$$c_0=\frac{a_0}{2}=\frac{1}{T}{\int }_{-\frac{T}{2}}^{\frac{T}{2}}{f}_T(x)dx$$

$$\begin{gathered} c_n=\frac{a_n-j{b}_n}{2}=\frac{1}{T}[{\int }_{-\frac{T}{2}}^{\frac{T}{2}}{f}_T(x)cosn\omega xdx-j{\int }_{-\frac{T}{2}}^{\frac{T}{2}}{f}_T(x)sinn\omega xdx] \\ =\frac{1}{T}{\int }_{-\frac{T}{2}}^{\frac{T}{2}}{f}_T(x)(cosn\omega x-jsinn\omega x)dx \\ =\frac{1}{T}{\int }_{-\frac{T}{2}}^{\frac{T}{2}}{f}_T(x)e^{-jn\omega x}dx \end{gathered}$$
其中$n=\pm 1,\pm 2,\dots ,$
若令$w_n=n\omega$，则
$$\begin{gathered} f_T(x)=c_0+\sum _{n=1}^{\infty }(c_n{e}^{j{w}_{n}x}+c_{-n}{e}^{-j{w}_{n}x})=\sum _{n=-\infty }^{\infty }{c}_n{e}^{j{w}_{n}x}=\sum _{n=-\infty }^{\infty }[\frac{1}{T}{\int }_{-\frac{T}{2}}^{\frac{T}{2}}{f}_T(x)e^{-j{w}_{n}x}dx]e^{j{w}_{n}x} \\ =\frac{1}{T}\sum _{n=-\infty }^{\infty }[{\int }_{-\frac{T}{2}}^{\frac{T}{2}}{f}_T(x)e^{-j{w}_{n}x}dx]e^{j{w}_{n}x}\text{\hspace{1em}(1.10)} \end{gathered}$$
我们来解读下(1.10)

$$\begin{gathered} f_T(x)=\frac{1}{T}\sum _{n=-\infty }^{\infty }[{\int }_{-\frac{T}{2}}^{\frac{T}{2}}{f}_T(x)e^{-j{w}_{n}x}dx]e^{j{w}_{n}x} \\ =\frac{1}{T}\sum _{n=-\infty }^{\infty }[F(w_n)]e^{j{w}_{n}x} \end{gathered}$$
上式表明，周期函数$f_T(x)$由幅值$F(w_n)$和频率$e^{j{w}_{n}x}$的累加和组成。