傅里叶级数理论详讲&实例应用

过冷水诚挚邀请你加入Matlab仿真秀官方交流群进行Matlab学习、问题咨询、 Matlab相关资料下载,群号:927550334


久久不更,更则必是大招,本期推文有理论有案例,有兴趣的读者仔细阅读。
过冷水之前有和大家讲傅里叶级数,并给出以一个函数用傅里叶级数近似的案例。本期就进一步详讲傅里叶级数。傅里叶级数展开时基底函数取1,cosx、sinx,cos2x、sin2x.....cosnx、sinnx,傅里叶级数一般情况下表示为:

a0、an、bn是展开系数。假定一个周期为2π的函数f(x 2π)=f(x)现在计算其系数。这就需要一点灵活的数学思想来解决问题。我们对上式两边在[0,2π]范围内积分。可得:

我们继续采用两边积分的方法求系数an、bn,对泰勒级数两边同时乘以cos(mx),然后在[0,2π]范围积分,如下:

则当m=n时:同理推bn:

许多实际问题中,函数f(x)是一个定义在有限区间0<x<X上的任意函数。对这样的函数如何展开傅里叶级数?大师告诉我们依旧可以用正弦函数、余弦函数展开:

现在求展开系数的表达式:

讲完傅里叶级数的理论,我们来看看实际傅里叶级数实际应用。过冷水在学习过程中遇到如下一幅图,需要对该函数进行积分求值:

数值图像求积分的问题过冷水的往期推文数值计算——MATLAB数值积分原理详讲、数值优化—三种复杂函数数值积分方法实例演示都有讲。过冷水本期的重点不是求数值积分的方法讨论,就想弄清楚如何很好的给出很好拟合该线段的函数表达式。这个函数图像算是比较复杂的,很少接触到所示的图像函数式,在此只能给出一个较为接近的拟合公式:


函数拟合代码很简单,后期会有详细推文讲解函数拟合问题,本期先先给出代码。整体来看还算拟合的好,部分拟合能够看出来有明显的错差,该方法并没有给出具体的函数解析式,类似于差值法,可以用于原图像的光滑处理。过冷水的实际问题是用具体函数表达式表示出第一个谷底区间的函数表达式,不需要给出整个图像区间的函数解析式,样问题就变得容易了,实际拟合一下图像:

使用函数拟合工具箱直接给出多项式拟合、高斯拟合、傅里叶拟合的函数图像,由拟合图像可知,傅里叶拟合和高斯拟合明显较好,使用起来比较简单。再一次告诫大家,在拟合图像的时候只知道多项式拟合就显得太不成熟了。

在此给出傅里叶拟合的代码:

y4即为想获得函数解析式 。过冷水在学习的过程中有接触到较为自己编辑傅里叶级数公式代码见从泰勒级数说傅里叶级数。

过冷水发表于 仿真秀 平台原创文章,未经授权禁止私自转载,如需转载请需要和作者沟通表明授权声明,未授权文章皆视为侵权行为,必将追责。如果您希望加入Matlab仿真秀官方交流群进行Matlab学习、问题咨询、 Matlab相关资料下载均可加群:927550334。
精品回顾
matlab绘制农夫过河动态图
分子动力学的原子空间运动轨迹演示编程
过冷水带你用matlab制作演示动画
python批量移动文件&重命名代码分享
过冷水和你分享 matlab读取存储各种文件的方法 文末有独家金曲分享
