前言
本人水平有限,文章仅作为自己的学习记录,且多用不严谨的、通俗易懂的语言去介绍。如有错误,希望大佬们多多批评指正。
一、历史小故事
傅里叶(1768-1830)是法国数学物理学家,1807年他在法国科学学会上发表了一篇论文,用正弦曲线来描述温度分布,并提出了一个论断:任何连续周期的信号,都可以表示为一组适当加权的正弦曲线之和。没错,就是现在我们学习的傅里叶级数!
审查这篇论文的人包括了拉格朗日和拉普拉斯,拉普拉斯及其他审查者同意发表该论文,然而拉格朗日坚决反对。最终法国科学学会迫于拉格朗日的威望,拒绝了傅里叶的论文,直到拉格朗日死后15年,这篇论文才被发表。
参考文献:《通信之道》杨学志
二、初识傅里叶级数概念
2.1什么是频率
如图2.1.1,我画了一个连续、周期(T=3/2)的方波。我们可以用傅里叶级数研究它。
频率f概念的基础来源于周期信号,时间是无穷的。f=1/T
但是我们经常接触到的信号并不是这样周期,时间无限的。例如我要研究一段持续时间为2s的来自微信聊天框的语音怎么办?别急,这样的信号叫连续、非周期的信号,我们可以用傅里叶变换研究它。
2.2为什么用正弦信号(三角函数)
在傅里叶级数的公式中,是将一个信号表示为很多个不同频率的正弦信号(sin和cos统称为正弦)之和,为什么要费这么大力气分解为这么多正弦信号?答曰:三角函数作为输入f(t)进入线性系统,出来之后的输出y(t),仍然是个三角函数,而且频率根本不会变!
初学者可能还不知道什么是傅里叶级数公式,什么是线性系统,也不知道经过线性系统后f不变意味着什么,别急,先知道这个事儿,有兴趣的话还可以先看看我证明这个结论!
如图,是一个不严谨的证明,但是可以说明这个特点,输入是一个正弦信号,输出也是一个正弦信号,且w=2*pi*f不变,只是幅度和相位有一定改变。
PS.初中学的正弦函数是y=sin(x),现在将x替换为wt或者2*pi*f*t,则叫正弦信号。
参考文献:《通信之道》杨学志
2.3用方波引入傅里叶级数
有了频率f和正弦信号的概念了之后,现在我们要正式进入傅里叶级数!我的介绍并不是从源头介绍傅里叶级数怎么来的,我只是介绍,这个傅里叶级数我们要怎么直观地看懂。
图2.1.1连续、周期的方波我用Matlab将其画出,注意我用Matlab只画了一个周期。
看到红色的曲线没有?它就是一系列的余弦信号叠加而成的,是不是最后叠加了之后有点像方波了?我现在把余弦信号加多一点试试,看看会不会拟合的更好!
显然,已经拟合的像一个方波了!不过在方波的t=正负0.5出,还有一点“突起”,这个叫吉布斯Gibbs现象。当用来叠加的余弦信号越多,这个误差会越小。
至此,我们已经初步的认识到了傅里叶级数传递给我们的信息:任何连续周期信号,都可以表示为一组适当加权的正弦曲线之和。
三、傅里叶级数公式学习
3.1三角形式的傅里叶级数公式
如图为傅里叶级数公式,an,bn为傅里叶级数系数。连续、周期信号f(t)可以被表达为一个直流a0,和无数个不同频率的正弦信号之和。
还是拿图2.1.1的方波为例子,这个方波是一个偶函数,在用公式计算bn的时候,因为积分的对称性质,bn全部等于0;在计算an的时候,系数如下图所示:
可以看到这个方波可以被分解为【直流】+【1、2、4、5、7、8、10、11、13、14......次谐波(余弦)】谐波越多,拟合的越像原信号。
3.1.1用实际生活的例子理解傅里叶级数
现实中的信号,例如乐器发出的声音,都可以看做连续时间信号,都可以在时域中用时间轴和幅度轴表示出来。对其作傅里叶级数展开,可以看到不同的谐波成分。在这里分享一个up主搬运的视频。真的超级棒!https://www.bilibili.com/video/av927069444/?vd_source=736edacaecc16480e6b3922dcdee5363https://www.bilibili.com/video/av927069444/?vd_source=736edacaecc16480e6b3922dcdee5363
3.1.2对傅里叶级数公式的一点理解
截止目前,梳理一下思路:我们先认识了频率、连续周期函数、三角函数的概念,然后通过一个方波展现了一个连续周期信号确实可以被分解为一系列正弦信号。最后给出了傅里叶级数公式,告诉我们具体由哪些频率的余弦、哪些频率的正弦组成,以及对应的加权是多少。
现在,我来按照自己的逻辑浅述一下我对傅里叶级数的一些感受和理解。
①本科教材提到过三角函数的正交完备性、正交函数集的概念,如下图所示。
正弦函数和余弦函数一定是正交的。
不同频率的正弦函数是正交的。
不同频率的余弦函数是正交的。
②那么也就是说{sin0t、cos0t、sin1t、cos1t、sin2t、cos2t、.......sinnt、cosnt}可以构成一个空间,在这个空间里的每一个信号彼此之间都是正交的。
③线性代数中有向量空间、正交基的概念,我们先看看二维空间,也就是我们接触的最多的平面直角坐标系。
可以用红色或绿色任意一组正交基表示二维空间中的任意向量。
④现在将空间扩展至无穷维度,傅里叶级数中的每一个频率成分都是一个维度。
{sin0t、cos0t、sin1t、cos1t、sin2t、cos2t、.......sinnt、cosnt}组成三角函数正交基,于是乎,可以表示无穷维空间中的任意信号。一旦f(t)被表示(即用傅里叶级数展开),f(t)的频率成分也就被一一分解出来。
⑤相位谱和幅度谱
如图f(t)用傅里叶级数展开,w=2*pi/T,w为基频,显然f(t)包含直流、基频(w)、5次谐波分量(5w)。
以频率w或f为横轴,以幅度A为纵轴,即画出幅度谱。完成了时域到频域的转换。
以频率w或f为横轴,以相位为纵轴,即画出相位谱。
参考文献:《通信之道》杨学志,鲜枣课堂《深入浅出数字信号处理》
3.2复指数形式的傅里叶级数公式
3.2.1由三角形式傅里叶级数推导复指数形式
这里的推导我直接上书上的图片了,写的太太太太清楚啦!!
参考文献:《通信之道》杨学志
3.2.2为什么出现了负频率?
经过欧拉公式推导后我们得到了cn,仍然以频率为横轴,幅度为纵轴,可以画出幅度谱。
注意现在n的范围从(-∞,+∞),此时幅度谱和三角形式的幅度谱相比,变成了双边谱。
三角形式的单边谱,直流不变,其它位置的幅度减半,再关于y轴对称,形成的双边谱即为cn。
关于cn双边谱的例子我就不举例了,因为其实我在学习过程中,对三角形式的理解更多,虽然傅里叶级数肯定是用复指数形式,但是我的理解只是利用欧拉公式做了一个变换更加简洁而已。如果您想对复指数形式的傅里叶级数有更多了解,可以参考up主,DR_CAN.
3.2.3关于复指数的思考
这里埋个坑,下一节写傅里叶变换,也会写关于复指数信号的内容。