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图像处理之镜像实现方法
May 30 Python
解决Python安装后pip不能用的问题
Jun 12 Python
Python使用Selenium模块模拟浏览器抓取斗鱼直播间信息示例
Jul 18 Python
Python中创建二维数组
Oct 17 Python
nginx黑名单和django限速,最简单的防恶意请求方法分享
Aug 09 Python
Python的条件锁与事件共享详解
Sep 12 Python
PyTorch学习:动态图和静态图的例子
Jan 06 Python
Python3读写Excel文件(使用xlrd,xlsxwriter,openpyxl3种方式读写实例与优劣)
Feb 13 Python
python 回溯法模板详解
Feb 26 Python
Anaconda+Pycharm环境下的PyTorch配置方法
Mar 13 Python
基于Pyinstaller打包Python程序并压缩文件大小
May 28 Python
keras实现VGG16 CIFAR10数据集方式
Jul 07 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
磨咖啡豆的密诀
2021/03/03 冲泡冲煮
织梦模板标记简介
2007/03/11 PHP
php调用淘宝开放API实现根据卖家昵称获取卖家店铺ID的方法
2015/07/29 PHP
PHP中FTP相关函数小结
2016/07/15 PHP
PHP获取表单数据与HTML嵌入PHP脚本的实现
2017/02/09 PHP
详解PHP神奇又有用的Trait
2019/03/25 PHP
TP框架实现上传一张图片和批量上传图片的方法分析
2020/04/23 PHP
通过pjax实现无刷新翻页(兼容新版jquery)
2014/01/31 Javascript
JavaScript使用shift方法移除素组第一个元素实例分析
2015/04/06 Javascript
javascript结合CSS实现苹果开关按钮特效
2015/04/07 Javascript
jquery实现表单验证并阻止非法提交
2015/07/09 Javascript
纯javascript实现分页(两种方法)
2015/08/26 Javascript
微信小程序 触控事件详细介绍
2016/10/17 Javascript
真正好用的js验证上传文件大小的简单方法
2016/10/27 Javascript
如何在Angular应用中创建包含组件方法示例
2019/03/23 Javascript
小试小程序云开发(小结)
2019/06/06 Javascript
JS简单表单验证功能完整示例
2020/01/26 Javascript
Vue实现计算器计算效果
2020/08/17 Javascript
[55:44]完美世界DOTA2联赛决赛 FTD vs Phoenix 第二场 11.08
2020/11/11 DOTA
详解JavaScript编程中的window与window.screen对象
2015/10/26 Python
Python的Flask框架中配置多个子域名的方法讲解
2016/06/07 Python
python搭建虚拟环境的步骤详解
2016/09/27 Python
pandas.dataframe按行索引表达式选取方法
2018/10/30 Python
Python自定义函数计算给定日期是该年第几天的方法示例
2019/05/30 Python
PyCharm+Qt Designer+PyUIC安装配置教程详解
2019/06/13 Python
使用Python代码实现Linux中的ls遍历目录命令的实例代码
2019/09/07 Python
CSS3动画animation实现云彩向左滚动
2014/05/09 HTML / CSS
魔声耳机官方网站:Monster是世界第一品牌的高性能耳机
2016/10/26 全球购物
沙龙级头发造型工具:FOXYBAE
2018/07/01 全球购物
机关道德讲堂实施方案
2014/03/15 职场文书
邹越感恩父母演讲稿
2014/08/28 职场文书
大学生村官个人对照检查材料(群众路线)
2014/09/26 职场文书
银行资信证明
2015/06/17 职场文书
mysql字符串截取函数小结
2021/04/05 MySQL
Django实现聊天机器人
2021/05/31 Python
Python中生成随机数据安全性、多功能性、用途和速度方面进行比较
2022/04/14 Python