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类的专用方法实例分析
Jan 09 Python
python网络编程之文件下载实例分析
May 20 Python
python类和继承用法实例
Jul 07 Python
Python切片索引用法示例
May 15 Python
在Python中使用Neo4j的方法
Mar 14 Python
Python爬虫 批量爬取下载抖音视频代码实例
Aug 16 Python
python3中利用filter函数输出小于某个数的所有回文数实例
Nov 24 Python
基于MSELoss()与CrossEntropyLoss()的区别详解
Jan 02 Python
python中JWT用户认证的实现
May 18 Python
Python爬虫定时计划任务的几种常见方法(推荐)
Jan 15 Python
python实现计算图形面积
Feb 22 Python
python前后端自定义分页器
Apr 13 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
Win2003下APACHE+PHP5+MYSQL4+PHPMYADMIN 的简易安装配置
2006/11/18 PHP
PHP 手机归属地查询 api
2010/02/08 PHP
php中将指针移动到数据集初始位置的实现代码[mysql_data_seek]
2012/11/01 PHP
Yii不依赖Model的表单生成器用法实例
2014/12/04 PHP
分享一个漂亮的php验证码类
2016/09/29 PHP
如何用JavaScript定义一个类
2014/09/12 Javascript
JavaScript Sort 的一个错误用法示例
2015/03/20 Javascript
JavaScript中的parse()方法使用简介
2015/06/12 Javascript
深入学习jQuery Validate表单验证
2016/01/18 Javascript
JS判断指定dom元素是否在屏幕内的方法实例
2017/01/23 Javascript
EasyUI Datebox 日期验证之开始日期小于结束时间
2017/05/19 Javascript
Node.JS循环删除非空文件夹及子目录下的所有文件
2018/03/12 Javascript
JavaScript树的深度优先遍历和广度优先遍历算法示例
2018/07/30 Javascript
Vue props 单向数据流的实现
2018/11/06 Javascript
详解如何用webpack4从零开始构建react开发环境
2019/01/27 Javascript
JavaScript跳出循环的三种方法(break, return, continue)
2019/07/30 Javascript
使用 Element UI Table 的 slot-scope方法
2019/10/10 Javascript
详解搭建一个vue-cli的移动端H5开发模板
2020/01/17 Javascript
js利用拖放实现添加删除
2020/08/27 Javascript
python实现猜数字游戏(无重复数字)示例分享
2014/03/29 Python
python获取当前计算机cpu数量的方法
2015/04/18 Python
Python基于staticmethod装饰器标示静态方法
2020/10/17 Python
Python环境使用OpenCV检测人脸实现教程
2020/10/19 Python
在Pycharm中安装Pandas库方法(简单易懂)
2021/02/20 Python
CSS3 清除浮动的方法示例
2018/06/01 HTML / CSS
印尼最大的网上书店:Gramedia.com
2018/09/13 全球购物
速比涛英国官网:Speedo英国
2019/07/15 全球购物
致跳高运动员加油稿
2014/02/12 职场文书
珍惜资源保护环境的建议书
2014/05/14 职场文书
医院深入开展党的群众路线教育实践活动实施方案
2014/08/27 职场文书
店铺转让协议书(2014版)
2014/09/23 职场文书
2015年社区消防安全工作总结
2015/10/14 职场文书
电力培训学习心得体会
2016/01/11 职场文书
Vue3.0 手写放大镜效果
2021/07/25 Vue.js
HDFS免重启挂载新磁盘
2022/04/06 Servers
Java数据结构之堆(优先队列)
2022/05/20 Java/Android