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生成随机mac地址的方法
Mar 16 Python
python统计文本字符串里单词出现频率的方法
May 26 Python
Python实现冒泡排序的简单应用示例
Dec 11 Python
Python网络编程之TCP套接字简单用法示例
Apr 09 Python
TensorFlow Session会话控制&Variable变量详解
Jul 30 Python
python中map的基本用法示例
Sep 10 Python
Python3 max()函数基础用法
Feb 19 Python
Python3.5面向对象程序设计之类的继承和多态详解
Apr 24 Python
python可视化篇之流式数据监控的实现
Aug 07 Python
python的Jenkins接口调用方式
May 12 Python
通过Django Admin+HttpRunner1.5.6实现简易接口测试平台
Nov 11 Python
python 基于opencv去除图片阴影
Jan 26 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定时执行任务的3种方法详解
2015/12/21 PHP
PHP中的使用curl发送请求(GET请求和POST请求)
2017/02/08 PHP
JavaScript Event学习第五章 高级事件注册模型
2010/02/07 Javascript
javascript动态向网页中添加表格实现代码
2014/02/19 Javascript
Javascript 函数parseInt()转换时出现bug问题
2014/05/20 Javascript
JavaScript实现自动变换表格边框颜色
2015/05/08 Javascript
浅谈jQuery构造函数分析
2015/05/11 Javascript
手机端页面rem宽度自适应脚本
2015/05/20 Javascript
JS实现转动随机数抽奖特效代码
2020/04/16 Javascript
轻松学习jQuery插件EasyUI EasyUI创建菜单与按钮
2015/11/30 Javascript
JavaScript编写检测用户所使用的浏览器的代码示例
2016/05/05 Javascript
JS中去掉array中重复元素的方法
2017/05/26 Javascript
vue语法之拼接字符串的示例代码
2017/10/25 Javascript
Vue2.0 http请求以及loading展示实例
2018/03/06 Javascript
在 Angular6 中使用 HTTP 请求服务端数据的步骤详解
2018/08/06 Javascript
微信小程序实现获取用户信息并存入数据库操作示例
2019/05/07 Javascript
vue element-ui读取pdf文件的方法
2019/11/26 Javascript
vue props 一次传多个值实例
2020/07/22 Javascript
[29:10]Ti4 冒泡赛第二天 NEWBEE vs Titan 3
2014/07/15 DOTA
剖析Python的Twisted框架的核心特性
2016/05/25 Python
python机器学习库常用汇总
2017/11/15 Python
Python制作词云图代码实例
2019/09/09 Python
Python使用tkinter模块实现推箱子游戏
2019/10/08 Python
python列表推导式入门学习解析
2019/12/02 Python
python读取yaml文件后修改写入本地实例
2020/04/27 Python
jupyter使用自动补全和切换默认浏览器的方法
2020/11/18 Python
一些Solaris面试题
2013/03/22 面试题
护士专业推荐信
2013/11/02 职场文书
表彰大会主持词
2014/03/26 职场文书
文化活动实施方案
2014/03/28 职场文书
感恩寄语大全
2014/04/11 职场文书
高中社区服务活动报告
2015/02/05 职场文书
2015年护士长个人工作总结
2015/04/24 职场文书
2015年除四害工作总结
2015/07/23 职场文书
2016入党积极分子党课学习心得体会
2015/10/09 职场文书
Promise面试题详解之控制并发
2021/05/14 面试题