python实现傅里叶级数展开的实现
2018-07-21 14:27
851 查看
傅立叶级数的介绍我就不说了,自己也是应用为主,之前一直觉得很难懂,但最近通过自己编程实现了一些函数的傅立叶级数展开之后对傅立叶 级数展开的概念比较清楚了
(1)函数如下
函数图象如下:
代码:
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()
结果如下:
(2)
转换代码:
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()
结果如下:
关于傅里叶级数展开的还有一些例子在我的另一篇博客https://www.geek-share.com/detail/2743154807.html也有讨论,通过以上的例子可以很好直观感受傅立叶级数的作用,在电气工程上对分析流经变压器的谐波分量、整流装置输出谐波等分析具有重要作用。
还有一点就是python的强大功能,利用for循环语句可以方便地实现级数展开的求和运算,自己感觉编程起来比MATLAB/Octave跟快
相关文章:Python实现快速傅里叶变换(FFT)
以上就是本文的全部内容,希望对大家的学习有所帮助,也希望大家多多支持脚本之家。
相关文章推荐
- c语言实现傅里叶级数展开
- c语言实现傅里叶级数展开
- Matlab 周期方波信号傅里叶级数展开
- 将大数组里面的小数组平行展开的实现(Making a flat list out of list of lists in Python)
- python实现傅立叶级数展开
- MachingLearning中的距离和相似性计算以及python实现
- python web py入门(13)- 实现发贴到论坛
- python实现简易数据库之三——join多表连接和group by分组
- Unicode 编码中文 + 未编码中文符号 转换成UTF-8 中文 python实现
- rm -rf python 实现 v0.1
- 两个栈实现队列和两个队列实现栈-Python实现-复杂度分析
- Python用threading实现多线程详解
- 机器学习实战学习1:python实现KNN
- Python中使用property实现类的特性
- bandit算法原理及Python实现
- 基于Python实现的爬虫源码(1)
- python实现多进程_multiprocessing
- python使用fcntl模块实现程序加锁功能示例
- Python实现发送邮件
- 简单scheme计算器-python实现