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的setuptools框架下生成egg的教程
Apr 13 Python
Python的Flask框架中实现分页功能的教程
Apr 20 Python
python操作ssh实现服务器日志下载的方法
Jun 03 Python
Python中条件判断语句的简单使用方法
Aug 21 Python
通过Python爬虫代理IP快速增加博客阅读量
Dec 14 Python
python GUI实例学习
Nov 21 Python
利用python将xml文件解析成html文件的实现方法
Dec 22 Python
Python文件打开方式实例详解【a、a+、r+、w+区别】
Mar 30 Python
这可能是最好玩的python GUI入门实例(推荐)
Jul 19 Python
详解Python3 中的字符串格式化语法
Jan 15 Python
用python发送微信消息
Dec 21 Python
pandas进行数据输入和输出的方法详解
Mar 23 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预定义常量
2006/12/25 PHP
PHP Directory 函数的详解
2013/03/07 PHP
使用php get_headers 判断URL是否有效的解决办法
2013/04/27 PHP
解析thinkphp中的导入文件标签
2013/06/20 PHP
PHP中提问频率最高的11个面试题和答案
2014/09/02 PHP
php实现约瑟夫问题的方法小结
2015/03/23 PHP
PHP strripos函数用法总结
2019/02/11 PHP
javascript之ESC(第二类混淆)
2007/05/06 Javascript
js限制textarea每行输入字符串长度的代码
2012/10/31 Javascript
JavaScrip常见的一些算法总结
2015/12/28 Javascript
JavaScript中将数组进行合并的基本方法讲解
2016/03/07 Javascript
javascript实现粘贴qq截图功能(clipboardData)
2016/05/29 Javascript
jQuery表单对象属性过滤选择器实例详解
2016/09/13 Javascript
前端常见跨域解决方案(全)
2017/09/19 Javascript
vue2.0.js的多级联动选择器实现方法
2018/02/09 Javascript
AngularJS使用Filter自定义过滤器控制ng-repeat去除重复功能示例
2018/04/21 Javascript
讲解vue-router之什么是嵌套路由
2018/05/28 Javascript
vue给组件传递不同的值方法
2018/09/29 Javascript
JavaScript之解构赋值的理解
2019/01/30 Javascript
微信小程序错误this.setData报错及解决过程
2019/09/18 Javascript
JS函数参数的传递与同名参数实例分析
2020/03/16 Javascript
JS端基于download.js实现图片、视频时直接下载而不是打开预览
2020/05/09 Javascript
[46:49]完美世界DOTA2联赛PWL S3 access vs Rebirth 第二场 12.19
2020/12/24 DOTA
Python将xml和xsl转换为html的方法
2015/03/10 Python
Python中二维列表如何获取子区域元素的组成
2017/01/19 Python
你眼中的Python大牛 应该都有这份书单
2017/10/31 Python
python3学习笔记之多进程分布式小例子
2018/02/13 Python
Python文件读写常见用法总结
2019/02/22 Python
基于Python执行dos命令并获取输出的结果
2019/12/30 Python
Canvas实现放大镜效果完整案例分析(附代码)
2020/11/26 HTML / CSS
教师查摆问题及整改措施
2014/10/11 职场文书
综治维稳工作汇报
2014/10/27 职场文书
社会实践活动报告
2015/02/05 职场文书
2016情人节宣传语
2015/07/14 职场文书
《比尾巴》教学反思
2016/02/24 职场文书
SQL实现LeetCode(196.删除重复邮箱)
2021/08/07 MySQL