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基础之函数用法实例详解
Sep 10 Python
python动态性强类型用法实例
May 09 Python
Python实现excel转sqlite的方法
Jul 17 Python
mac系统安装Python3初体验
Jan 02 Python
Python实现多属性排序的方法
Dec 05 Python
Python之列表实现栈的工作功能
Jan 28 Python
Django框架验证码用法实例分析
May 10 Python
Django框架封装外部函数示例
May 28 Python
Python手绘可视化工具cutecharts使用实例
Dec 05 Python
python 插入日期数据到Oracle实例
Mar 02 Python
python中关于数据类型的学习笔记
Jul 19 Python
如何用Matlab和Python读取Netcdf文件
Feb 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 进度条实现代码
2009/03/10 PHP
解决laravel 出现ajax请求419(unknown status)的问题
2019/09/03 PHP
设为首页和收藏的Javascript代码(亲测兼容IE,Firefox,chrome等浏览器)
2013/11/18 Javascript
jQuery实现鼠标滑向当前图片高亮显示并且其它图片变灰的方法
2015/07/27 Javascript
详解Nodejs之npm&package.json
2017/06/15 NodeJs
node使用UEditor富文本编辑器的方法实例
2017/07/11 Javascript
mui back 返回刷新页面的实例
2017/12/06 Javascript
通过jquery的ajax请求本地的json文件方法
2018/08/08 jQuery
angularjs1.5 组件内用函数向外传值的实例
2018/09/30 Javascript
ztree加载完成后显示勾选节点的实现代码
2018/10/22 Javascript
javascript实现blob加密视频源地址的方法
2019/08/08 Javascript
基于JS抓取某高校附近共享单车位置 使用web方式展示位置变化代码实例
2019/08/27 Javascript
vue $router和$route的区别详解
2020/12/02 Vue.js
[01:06:07]2014 DOTA2国际邀请赛中国区预选赛5.21 DT VS CIS
2014/05/22 DOTA
Python自定义主从分布式架构实例分析
2016/09/19 Python
Python 编码规范(Google Python Style Guide)
2018/05/05 Python
TensorFlow实现卷积神经网络
2018/05/24 Python
selenium设置proxy、headers的方法(phantomjs、Chrome、Firefox)
2018/11/29 Python
Python实现Word表格转成Excel表格的示例代码
2020/04/16 Python
python模拟点击在ios中实现的实例讲解
2020/11/26 Python
有影响力的品牌之家:Our Social Collective
2019/06/08 全球购物
关于Assembly命名空间的三个面试题
2015/07/23 面试题
城市规划毕业生求职信
2013/10/10 职场文书
工商管理专业应届生求职信
2013/11/04 职场文书
岗位职责定义及内容
2013/11/08 职场文书
绿化先进工作者事迹材料
2014/01/30 职场文书
第一批党的群众路线教育实践活动工作总结
2014/03/03 职场文书
大学生安全责任书
2014/07/25 职场文书
放弃继承权公证书
2015/01/23 职场文书
行政主管岗位职责范本
2015/04/09 职场文书
民间借贷纠纷答辩状
2015/08/03 职场文书
nginx 反向代理之 proxy_pass的实现
2021/03/31 Servers
Python List remove()实例用法详解
2021/08/02 Python
MySQL七大JOIN的具体使用
2022/02/28 MySQL
十大经典日本动漫排行榜 海贼王第三,犬夜叉仅第八
2022/03/18 日漫
CSS子盒子水平和垂直居中的五种方法
2022/07/23 HTML / CSS