傅里叶(Fourier)展开

傅里叶级数的定义如下:

设函数f(x)f(x)f(x)在区间[0,2π][0,2\pi][0,2π]上绝对可积,令
{an=1π∫02πf(x)cos⁡nxdx(n=0,1,2,⋯)bn=1π∫02πf(x)sin⁡nxdx(n=1,2,⋯)\left\{\begin{array}{ll} a_{n}=\frac{1}{\pi} \int_{0}^{2 \pi} f(x) \cos n x \mathrm{~d} x & (n=0,1,2, \cdots) \\ b_{n}=\frac{1}{\pi} \int_{0}^{2 \pi} f(x) \sin n x \mathrm{~d} x & (n=1,2, \cdots) \end{array}\right. {an​=π1​∫02π​f(x)cosnx dxbn​=π1​∫02π​f(x)sinnx dx​(n=0,1,2,⋯)(n=1,2,⋯)​
以an,bna_n,b_nan​,bn​为系数作三角级数
a02+∑n=1∞(ancosnx+bnsinnx)\frac{a_0}{2}+\sum_{n=1}^{\infty}(a_ncosnx+b_nsinnx) 2a0​​+n=1∑∞​(an​cosnx+bn​sinnx)
称为f(x)f(x)f(x)的傅里叶级数,an,bna_n,b_nan​,bn​称为傅里叶系数。

具体的傅里叶级数推导过程可以参考下面这篇文章:

https://zhuanlan.zhihu.com/p/41455378

根据定义可以编写计算傅里叶系数的函数:

function [a0,an,bn] = Fourierzpi(f)
%本函数旨在计算[0,2pi]上的傅里叶系数
syms x n
a0 = int(f,0,2*pi)/pi;
an = int(f*cos(n*x),0,2*pi)/pi;
bn = int(f*sin(n*x),0,2*pi)/pi;
end

利用该函数计算sinx+exx2\frac{sinx+e^x}{x^2}x2sinx+ex​在[0,2pi]上的傅里叶系数

syms x
f = (sin(x)+exp(x))/x^2;
[a0,an,bn] = Fourierzpi(f)a0 =Infan =int((cos(n*x)*(exp(x) + sin(x)))/x^2, x, 0, 2*pi)/pibn =int((sin(n*x)*(exp(x) + sin(x)))/x^2, x, 0, 2*pi)/pi

如果想要计算不同区间的傅里叶变化,只需要改变函数中范围,这样就可以计算不同区间,不同函数的傅里叶系数。

MATLAB之傅里叶展开(五)相关推荐

  1. 自适应滤波器在matlab仿真的程序_电气信息类专业课程之matlab系统仿真 第五章 BPSK通信系统(3)...

    继续讲解! 上一篇文章提出了那么多问题, 不知道大家是否能回答上来啊? 如果回答不了, 有没有去问问度娘呢? 程序写完了,回想了一下辅导2017届学生毕业设计的过程,那是真累.我要求他们在现有程序基础 ...

  2. 数学建模常用算法汇总及python,MATLAB实现(五) —— 拟合

    拟合 比较重要的就是2.1和2.3 2.2可以浅看一下, 自己敲着试一试 就拟合部分来说, MATLAB比python强大很多, 自带cftool工具包, 并且有很多快速的函数, 个人建议使用MATL ...

  3. MATLAB学习第五章矩阵

    MATLAB学习第五章 5.1矩阵 矩阵元素修改 矩阵运算 5.1矩阵 eye(n):创建n×n的矩阵 eye(m,n):创建m×n的矩阵 eye(size(A)):创建与A维数相同的单位阵 ones ...

  4. matlab 将图片立体化,MATLAB基础(五)——绘制立体图

    五.绘制立体图 基本XYZ立体绘图命令 在科学目视表示(Scientific visualization)中,三度空间的立体图是一个非常重要的技巧.本文将介绍MATLAB基本XYZ三度空间的各项绘图命 ...

  5. matlab simulink_简单五步实现 MATLAB/Simulink 锂电池建模

    锂电池是目前在各个能源密集型行业中用途广泛,例如新能源汽车.电力微网.航空航天等.电池模型的建立对研究电池的特性.SOC(state-of-charge)估计.SOH(state-of-health) ...

  6. matlab能画五维吗,进化算法之粒子群算法和Matlab实现(多维)

    前面一篇文章介绍了遗传算法,这里再介绍一种进化算法,称为粒子群算法.同遗传算法类似,粒子群算法也是仿照了自然界的生物现象得到的.这种现象就是鸟群在某个未知空间内寻找食物这一思想. 鸟群通过自身经验和种 ...

  7. 一个五位数取前三位matlab,【有五个小球,分别是1,2,3,4,5号,有放回的从中取三次,每次取一个,...-前三后五取一颗-数学-关偈邓同学...

    概述:本道作业题是关偈邓同学的课后练习,分享的知识点是前三后五取一颗,指导老师为习老师,涉及到的知识点涵盖:[有五个小球,分别是1,2,3,4,5号,有放回的从中取三次,每次取一个,...-前三后五取 ...

  8. matlab激光加工,科协五分钟|用Matlab实现激光扫描测距绘制面型

    科协五分钟|用Matlab实现激光扫描测距绘制面型 这次给大家介绍以下如何用Matlab控制步进电机和激光器进行激光扫描,并且将采集的点转化为图像.PS:如果你是科创大佬,或是上过专业认知实践课程,请 ...

  9. matlab篮球队需要五名队员,MATLAB应用与数学欣赏.doc

    MATLAB应用与数学欣赏 概 述:MATLAB具有弱小的计算功用和绘图功用,熟练运用MATLAB已成为大学生.硕士生.博士生必需把握 的基本技艺.经过MATLAB绘图.计算和数学建模综合运用等为学习 ...

最新文章

  1. python判断语句if...elif...else、if嵌套、三目运算符
  2. 内功重修之数据结构----数组
  3. ARM微处理器的体系结构
  4. IOC操作Bean管理XML方式(P名称空间注入)
  5. 优化搜索排序结果从而“ 提升CTR、CVR业务指标”
  6. 【百度地图API】多家地图API内存消耗对比测验(带源码)
  7. 力扣--125验证回文串/680验证回文字符串II
  8. HDOJ水题集合11:桶排序, 折半搜索
  9. 苹果设备型号代码 device model id / device codes(更新至iPhone 13 / iPhone SE3 / iPad Air 5代
  10. 一般期刊发表要求如何掌握
  11. ABAP报表添加点击跳转事务的方法
  12. 老师的经典口头禅,这一句最扎心
  13. Windows11中无法找到Office 2019 application的解决方案
  14. 回归基础:通过亮度恒定性和运动平滑度对光流进行无监督学习
  15. 用CSS实现设置网页文字为只读,不可复制的方法
  16. 抖音测试小程序功能开发
  17. android view.isshown,android View.isShown() 和 getVisibility() 的区别
  18. 后端 php mysql_搭建并部署PHP + MySQL后端环境
  19. 超时锁定计算机,Win10电脑设置锁定屏幕超时怎么办
  20. 高等数学:第三章 微分中值定理与导数的应用(8)曲率

热门文章

  1. 6-1 简单快速排序分数 10作者 唐艳琴单位 中国人民解放军陆军工程大学本题要求实现一个函数,可快速查找给定x(保证是整个数据中其值存在,如果x有多个,查找第一个x)在整个数据中的排名(数据
  2. Android 高德地图中路线规划绘制界面线路
  3. 理工男,我感觉 win10 挺好用的,一点不卡,为什么有那么多人买Macbook 呢?
  4. 3月16日---3月20日课程表
  5. 恢复被病毒强制隐藏的文件夹
  6. python线程和c++线程的区别_Python、线程、吉尔和C++
  7. MacBook Pro 时间机器备份(完美解决连接移动硬盘无反应)
  8. C语言中关于向函数中传入结构指针的易错点及解决方案
  9. 清华大学计算机学院武永卫,清华大学计算机科学与技术系导师简介:徐明伟
  10. 基于tkinter库txt文件操作的Python界面化图书管理系统