leinlin

3,459 人赞同了该文章

这是一个傅里叶变化系列的公式推导及其编程应用，公式上有什么不对的，大家可以随时在评论区给我留言，我一定积极修改，不误人子弟。

目前系列文章有：

傅里叶系列（一）傅里叶级数的推导3459 赞同

傅里叶系列（二）傅里叶变换的推导608 赞同 · 87 评论

傅里叶系列（三）离散傅里叶变换（DFT）184 赞同 · 23 

首先，隆重推出傅里叶级数的公式，不过这个东西属于“文物”级别的，诞生于19世纪初，因为傅里叶他老人家生于1768年，死于1830年。

　　但傅里叶级数在数论、组合数学、信号处理、概率论、统计学、密码学、声学、光学等领域都有着广泛的应用，这不由得让人肃然起敬。一打开《信号与系统》、《锁相环原理》等书籍，动不动就跳出一个“傅里叶级数”或“傅里叶变换”，弄一长串公式，让人云山雾罩。

　　如下就是傅里叶级数的公式：

(1)f(t)=a02+a1cos(ωt)+b1sin(ωt)+a2cos(2ωt)+b2sin(2ωt)+...=a02+∑n=1∞[ancos(nωt)+bnsin(nωt)]

其中：

(2)an=2T∫t0t0+Tf(t)cos(nωt)dt(3)bn=2T∫t0t0+Tf(t)sin(nωt)dt

不客气地说，这个公式可以说是像“臭婆娘的裹脚布——又臭又长”，而且来历相当蹊跷，不知那个傅里叶什么时候灵光乍现，把一个周期函数f(t)硬生生地写成这么一大堆东西。单看那个{1}式，就是把周期函数f(t)描述成一个常数系数a0、及1倍ω的sin和cos函数、2倍ω的sin和cos函数等、到n倍ω的sin和cos函数等一系列式子的和，且每项都有不同的系数，即An和Bn，至于这些系数，需要用积分来解得，即{2}{3}{4}式，不过为了积分方便，积分区间一般设为[-π, π]，也相当一个周期T的宽度。

　　能否从数学的角度推导出此公式，以使傅里叶级数来得明白些，让我等能了解它的前世今生呢？下面来详细解释一下此公式的得出过程：

１、把一个周期函数表示成三角级数：

　　首先，周期函数是客观世界中周期运动的数学表述，如物体挂在弹簧上作简谐振动、单摆振动、无线电电子振荡器的电子振荡等，大多可以表述为：

(4)f(t)=Asin(ωt+ψ)

　　这里 t 表示时间， A 表示振幅， ω 为角频率， ψ 为初相（与考察时设置原点位置有关，可以理解为一个常量）。

　　然而，世界上许多周期信号并非正弦函数那么简单，如方波、三角波等。傅里叶就想，能否用一系列的三角函数Ansin(nωt+ψn) 之和来表示那个较复杂的周期函数 f(t) 呢？因为正弦函数sin可以说是最简单的周期函数了。于是，傅里叶写出下式：（关于傅里叶推导纯属猜想、已经有知乎的网友指出是解热方程和弦振动导出的，有机会找找相关资料）

(5)f(t)=A0+∑n=1∞Ansin(nωt+ψn)

　　这里，t是变量，其他都是常数。与上面最简单的正弦周期函数相比，5式中多了一个n，且n从1到无穷大。这里f(t)是已知函数，也就是需要分解的原周期函数。从公式5来看，傅里叶是想把一个周期函数表示成许多正弦函数的线性叠加，这许许多多的正弦函数有着不同的幅度分量（即式中An）、有不同的周期或说是频率（是原周期函数的整数倍，即n）、有不同的初相角（即 ψn ），当然还有一项常数项（即 A0 ）。要命的是，这个n是从1到无穷大，也就是是一个无穷级数。

这里强调一下，傅里叶级数中对不同频率的波有一个要求就是给定一个初始的频率 ω0 ,之后的角频率必须是 ω0 的整数倍， 这就是DTF（离散傅里叶变化）中的角频率取值的原则。

应该说，傅里叶是一个天才，想得那么复杂。一般人不太会把一个简单的周期函数弄成这么一个复杂的表示式。但傅里叶认为，式子右边一大堆的函数，其实都是最简单的正弦函数，有利于后续的分析和计算。当然，这个式能否成立，关键是级数中的每一项都有一个未知系数，如A0、An等，如果能把这些系数求出来，那么5式就可以成立。当然在5式中，唯一已知的就是原周期函数f(t),那么只需用已知函数f(t)来表达出各项系数，上式就可以成立，也能计算了。

因为 ψn 是个常数， An 也是常数。解过常微分方程的人都知道，方程中的常数能整合到一起就整合到一起。

于是乎，傅里叶首先对式5作如下变形：

Ansin(nωt+ψn)=Ansinψncos(nωt)+Ancosψnsin(nωt)

这个变化并不陌生，源自于三角公式：

sin(α±β)=sinαcosβ±cosαsinβ

式中，蓝色项即为我们需要合并的常数项，

记：an=An⋅sinψn ， bn=An⋅cosψn

这样，公式{5}就可以写成如下公式{6}的形式：

(6)f(t)=A0+∑n=1∞[ancos(nωt)+bnsin(nωt)]

到了这一步我们只要解出 A0 、 an、 bn 的值即可。

2、麦克劳林公式中的待定系数法：

这里为解出 A0 、 an、 bn值奠定下思路：

泰勒级数即为任意一个函数都可以用一个多项式来逼近，记为：

f(x)=A+Bx+Cx2+Dx3+...备注：2,3等数字是X的幂，看原文2023.6.7

那么，麦克劳林令：

f′(x)=B+2Cx+3Dx2f″(x)=2C+6Dx........

在每个等式中令x=0,然后使用待定系数法就可以解出A,B,C...的值

A=f(0)

B=f′(0)

C=f″(0)/2

D=f‴(0)/(1∗2∗3)

即：

N=fn(x)/n!        此处fn，n是f的n次导数，看原文2023.6.7

而众所周知三角函数在一个周期内的积分为0，如图

我们只要对(6)左右进行积分后即可求出 A0 的值，然后依次代入即可解出 an 、 bn 使用 f(t) 表达的公式。

3、三角函数的正交性：

　　这是为下一步傅里叶级数展开时所用积分的准备知识。一个三角函数系：1，cosx , sinx , cos2x , sin2x , … , cosnx , sinnx , … 如果这一堆函数（包括常数1）中任何两个不同函数的乘积在区间[-π, π]上的积分等于零，就说三角函数系在区间[-π, π]上正交，即有如下式子：

∫−ππcos nxdx=0(n=1,2,3,...)∫−ππsin nxdx=0(n=1,2,3,...)∫−ππsin kx⋅cos nxdx=0(k,n=1,2,3,...;k≠n)∫−ππcos kx⋅cos nxdx=0(k,n=1,2,3,...;k≠n)∫−ππsin kx⋅sin nxdx=0(k,n=1,2,3,...;k≠n)

　　以上各式在区间[-π, π]的定积分均为0，第１第２式可视为三角函数cos和sin与１相乘的积分；第3-5式则为sin和cos的不同组合相乘的积分式。除了这5个式子外，不可能再有其他的组合了。注意，第4第5两个式中，k不能等于n，否则就不属于“三角函数系中任意两个不同函数”的定义了，变成同一函数的平方了。但第3式中，k与n可以相等，相等时也是二个不同函数。下面通过计算第4式的定积分来验证其正确性，第4式中二函数相乘可以写成：

cos kx⋅cos nx=12[cos(k+n)x+cos(k−n)x]sin kx⋅sin nx=−12[cos(k+n)x−cos(k−n)x]sin kx⋅cos nx=12[sin(k+n)x+sin(k−n)x]

当 k≠n 时，有：

∫−ππcos kx⋅cos nxdx=12∫−ππ[cos(k+n)x+cos(k−n)x]dx=12[sin(k+n)xk+n + sin(k−n)xk−n]|−ππ=12[0+0]=0

可见在指定[-π, π]的区间里，该式的定积分为0。其他式也可逐一验证。

4、函数展开成傅里叶级数：

　　先把傅里叶级数表示为下式，即⑥式：

(6)f(t)=A0+∑n=1∞[ancos(nωt)+bnsin(nωt)]

　　对⑥式从[-π, π]积分，得：

∫−ππf(t)=∫−ππA0+∫−ππ∑n=1∞[ancos(nωt)+bnsin(nωt)]=∫−ππA0+0=A0|−ππ=(π−(−π))A0=2πA0

解出： A0=12π∫−ππf(t)

这就求得了第一个系数 A0 的表达式，即最上边傅里叶级数公式里的(2)式。接下来再求 an 和 bn 的表达式。用 cos(kωt) 乘{6}式的二边得：

∑n=1∞[ancos(nωt)+bnsin(nωt)

f(t)⋅cos(kωt)=A0cos(kωt)+∑n=1∞[ancos(nωt)⋅cos(kωt)+bnsin(nωt)⋅cos(kωt)]

然后对上式从 −π 到 π 逐项积分：

∫−ππf(t)⋅cos(kωt)dt=A0∫−ππcos(kωt)dt+∑n=1∞[an∫−ππcos(nωt)⋅cos(kωt)dt+bn∫−ππsin(nωt)⋅cos(kωt)dt]

　　根据三角函数系的正交性，红色积分为0，蓝色项中仅当 k=n 时积分不为0，其余项积分为0，所以有：

半角公式∫−ππcos(kωt)⋅f(t)=an∑n=1∞∫−ππcos(kωt)⋅cos(nωt)dt=an∫−ππcos2(nωt)dt=an2∫−ππ(1+cos2nωt)dt(半角公式)=an2(∫−ππ1dt+∫−ππcos2nωt dt)=an2⋅2π=anπ

解得：

an=1π∫−ππcos(nωt)⋅f(t)dt(k=n)

同理用 sin(kωt) 乘{6}式的二边得：

bn=1π∫−ππsin(nωt)⋅f(t)dt(k=n)

我们发现 A0 的分母为 2π 而 、an、bn 为 π ,为了统一分母我们令 a0=2A0 有：

a0=2T∫−ππf(t)dt

(6)变形为：

f(t)=a02+∑n=1∞[ancos(nωt)+bnsin(nωt)]

推导的时候假设 T=2π ,代入即可得到(2)、(3)、(4)

(2)an=2T∫t0t0+Tf(t)cos(nωt)dt(3)bn=2T∫t0t0+Tf(t)sin(nωt)dt

至此，已经求得傅里叶级数中各系数的表达式，当然这里有个条件： ∫t0t0+Tf(t)dt 积分存在，这里涉及到勒贝格可积的问题，因为离散傅里叶变化涉及到周期内有无限个可去间断点的问题，狄利克雷条件仅仅是个充分条件，一个函数有傅里叶级数但是它也存在无限个间断点以及极大值极小值比如方波信号。

至于勒贝格可积有空另开篇文章进行证明。

综上，傅里叶级数的产生过程可以分为以下三步：

1、设想可以把一个周期函数f(t)通过最简单的一系列正弦函数来表示，即5式；

2、通过变形后用三角级数（含sin和cos）来表示；

3、通过积分，把各未知系数用f(t)的积分式来表达；

4、最后得到的4个表达式就是傅里叶级数公式。

没错一切都是 命运石之门的选择

傅里叶系列（一）傅里叶级数的推导 - 知乎