python实现傅里叶级数展开的实现


Posted in Python onJuly 21, 2018

傅立叶级数的介绍我就不说了,自己也是应用为主,之前一直觉得很难懂,但最近通过自己编程实现了一些函数的傅立叶级数展开之后对傅立叶 级数展开的概念比较清楚了

python实现傅里叶级数展开的实现

(1)函数如下

python实现傅里叶级数展开的实现

函数图象如下:

python实现傅里叶级数展开的实现 

代码:

from pylab import *
x = mgrid[-10:10:0.02] # 这里类似于MATLAB用冒号产生步长为0.02的序列,但是语法和MATLAB不同
n = arange(1,1000)
def fourier_transform():
 a0 = (1-exp(-pi))/pi+1
 s=a0/2
 for i in range(1,100,1):
 s0 = ( (1-(-1)**i*exp(-pi))/(pi*(1+i**2))*cos(i*x)+1/pi*( (-i*(1-(-1)**i*exp(-pi)))/(1+i**2) + (1-(-1)**i)/i ) * sin(i*x) )
 s=s+s0
 plot(x,s,'orange',linewidth=0.6)
 title('fourier_transform')
 show() 

fourier_transform()

结果如下:

python实现傅里叶级数展开的实现

(2)

python实现傅里叶级数展开的实现

python实现傅里叶级数展开的实现

转换代码:

from pylab import *
def fourier1():
 s=pi/2

 for i in range(1,100,1):
 s0 = 2/pi*(1-(-1)**i)/i**2 * cos(i*x) 
 s=s+s0
 plot(x,s,'orange',linewidth=0.6)
 title('fourier1')
 show()
fourier1()

结果如下:

python实现傅里叶级数展开的实现

关于傅里叶级数展开的还有一些例子在我的另一篇博客https://3water.com/article/144196.htm也有讨论,通过以上的例子可以很好直观感受傅立叶级数的作用,在电气工程上对分析流经变压器的谐波分量、整流装置输出谐波等分析具有重要作用。

还有一点就是python的强大功能,利用for循环语句可以方便地实现级数展开的求和运算,自己感觉编程起来比MATLAB/Octave跟快

相关文章:Python实现快速傅里叶变换(FFT)

以上就是本文的全部内容,希望对大家的学习有所帮助,也希望大家多多支持三水点靠木。

Python 相关文章推荐
Python打包文件夹的方法小结(zip,tar,tar.gz等)
Sep 18 Python
基于pandas数据样本行列选取的方法
Apr 20 Python
对Python信号处理模块signal详解
Jan 09 Python
python的pytest框架之命令行参数详解(下)
Jun 27 Python
Python Lambda函数使用总结详解
Dec 11 Python
Python csv文件记录流程代码解析
Jul 16 Python
python 如何区分return和yield
Sep 22 Python
python中Pexpect的工作流程实例讲解
Mar 02 Python
解决numpy数组互换两行及赋值的问题
Apr 17 Python
Python批量将csv文件转化成xml文件的实例
May 10 Python
Python数据分析之绘图和可视化详解
Jun 02 Python
python中mongodb包操作数据库
Apr 19 Python
Python实现快速傅里叶变换的方法(FFT)
Jul 21 #Python
Python实现获取本地及远程图片大小的方法示例
Jul 21 #Python
opencv python 傅里叶变换的使用
Jul 21 #Python
Numpy中的mask的使用
Jul 21 #Python
Flask框架使用DBUtils模块连接数据库操作示例
Jul 20 #Python
Flask框架WTForm表单用法示例
Jul 20 #Python
Python使用pymongo模块操作MongoDB的方法示例
Jul 20 #Python
You might like
PHP 引用是个坏习惯
2010/03/12 PHP
php图片的裁剪与缩放生成符合需求的缩略图
2013/01/11 PHP
Yii框架使用PHPExcel导出Excel文件的方法分析【改进版】
2019/07/24 PHP
基于json的jquery地区联动效果代码
2011/07/06 Javascript
jQuery EasyUI API 中文文档 - Panel面板
2011/09/30 Javascript
关于function类中定义变量this的简单说明
2016/05/28 Javascript
jQueryUI DatePicker 添加时分秒
2016/06/04 Javascript
Js查找字符串中出现次数最多的字符及个数实例解析
2016/09/05 Javascript
详解vue项目构建与实战
2017/06/27 Javascript
JS沙箱模式实例分析
2017/09/04 Javascript
Angular实现预加载延迟模块的示例
2017/10/12 Javascript
select获取下拉框的值 下拉框默认选中方法
2018/02/28 Javascript
在Vue-cli里应用Vuex的state和mutations方法
2018/09/16 Javascript
javascript中关于类型判断的一些疑惑小结
2018/10/14 Javascript
详解javascript对数组和json数组的操作
2019/04/15 Javascript
vue中的过滤器实例代码详解
2019/06/06 Javascript
layui form表单提交之后重新加载数据表格的方法
2019/09/11 Javascript
js实现无限层级树形数据结构(创新算法)
2020/02/27 Javascript
angular组件间传值测试的方法详解
2020/05/07 Javascript
Python实例之wxpython中Frame使用方法
2014/06/09 Python
Python使用Turtle模块绘制五星红旗代码示例
2017/12/11 Python
基于Python的文件类型和字符串详解
2017/12/21 Python
python xlsxwriter库生成图表的应用示例
2018/03/16 Python
详解Python做一个名片管理系统
2019/03/14 Python
Django组件cookie与session的具体使用
2019/06/05 Python
对python中的装包与解包实例详解
2019/08/24 Python
python绕过图片滑动验证码实现爬取PTA所有题目功能 附源码
2021/01/06 Python
css3实现信纸/同学录效果的示例代码
2018/12/11 HTML / CSS
html5 canvas绘制放射性渐变色效果
2018/01/04 HTML / CSS
Under Armour安德玛法国官网:美国高端运动科技品牌
2018/06/29 全球购物
有创意的广告词
2014/03/18 职场文书
纪念九一八事变演讲稿:勿忘国耻
2014/09/14 职场文书
先进个人事迹材料范文
2014/12/30 职场文书
2015年大学班级工作总结
2015/04/28 职场文书
繁星春水读书笔记
2015/06/30 职场文书
Java中try catch处理异常示例
2021/12/06 Java/Android