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入门篇之数字
Oct 20 Python
使用Python编写类UNIX系统的命令行工具的教程
Apr 15 Python
实例讲解Python中函数的调用与定义
Mar 14 Python
基于Python 的进程管理工具supervisor使用指南
Sep 18 Python
Python+matplotlib实现填充螺旋实例
Jan 15 Python
Python使用googletrans报错的解决方法
Sep 25 Python
在Python中输入一个以空格为间隔的数组方法
Nov 13 Python
pandas 数据索引与选取的实现方法
Jun 21 Python
Django框架基础模板标签与filter使用方法详解
Jul 23 Python
手写一个python迭代器过程详解
Aug 27 Python
python3 requests库实现多图片爬取教程
Dec 18 Python
python识别验证码的思路及解决方案
Sep 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
php中文字符截取防乱码
2008/03/28 PHP
php单件模式结合命令链模式使用说明
2008/09/07 PHP
最常用的8款PHP调试工具
2014/07/06 PHP
Laravel中七个非常有用但很少人知道的Carbon方法
2017/09/21 PHP
php新建文件的方法实例
2019/09/26 PHP
js Array操作的最简短最容易理解方法
2013/12/09 Javascript
Javascript实现禁止输入中文或英文的例子
2014/12/09 Javascript
js中flexible.js实现淘宝弹性布局方案
2020/06/23 Javascript
下雪了 javascript实现雪花飞舞
2020/08/02 Javascript
JS简单实现DIV相对于浏览器固定位置不变的方法
2016/06/17 Javascript
Javascript中的 “&” 和 “|” 详解
2017/02/02 Javascript
Angular.JS内置服务$http对数据库的增删改使用教程
2017/05/07 Javascript
Mint UI 基于 Vue.js 移动端组件库
2017/11/07 Javascript
vue将时间戳转换成自定义时间格式的方法
2018/03/02 Javascript
[10:21]2018DOTA2国际邀请赛寻真——Winstrike
2018/08/11 DOTA
[45:52]完美世界DOTA2联赛PWL S3 Forest vs INK ICE 第二场 12.09
2020/12/12 DOTA
python抓取网页中的图片示例
2014/02/28 Python
深入Python函数编程的一些特性
2015/04/13 Python
Python3学习urllib的使用方法示例
2017/11/29 Python
用matplotlib画等高线图详解
2017/12/14 Python
python3获取当前文件的上一级目录实例
2018/04/26 Python
python3之模块psutil系统性能信息使用
2018/05/30 Python
使用OpCode绕过Python沙箱的方法详解
2019/09/03 Python
tensorboard 可以显示graph,却不能显示scalar的解决方式
2020/02/15 Python
Django 拼接两个queryset 或是两个不可以相加的对象实例
2020/03/28 Python
Python包资源下载路径报404解决方案
2020/11/05 Python
历史教育专业个人求职信
2013/12/13 职场文书
工厂采购员岗位职责
2014/04/08 职场文书
技术比武方案
2014/05/19 职场文书
学校地质灾害防治方案
2014/06/10 职场文书
不听老师话的万能检讨书
2014/10/04 职场文书
原告离婚代理词
2015/05/23 职场文书
二审代理词范文
2015/05/25 职场文书
教你快速开启Apache SkyWalking的自监控
2021/04/25 Servers
Python实现生成bmp图像的方法
2021/06/13 Python
PostgreSQL数据库创建并使用视图以及子查询
2022/04/11 PostgreSQL