python实现傅里叶级数展开的实现

站长资源 2024/11/27 佚名

13 1538 13

DDR爱好者之家 Design By 杰米

傅立叶级数的介绍我就不说了，自己也是应用为主，之前一直觉得很难懂，但最近通过自己编程实现了一些函数的傅立叶级数展开之后对傅立叶级数展开的概念比较清楚了

(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.jb51.net/article/144196.htm也有讨论，通过以上的例子可以很好直观感受傅立叶级数的作用，在电气工程上对分析流经变压器的谐波分量、整流装置输出谐波等分析具有重要作用。

还有一点就是python的强大功能，利用for循环语句可以方便地实现级数展开的求和运算，自己感觉编程起来比MATLAB/Octave跟快

相关文章：Python实现快速傅里叶变换（FFT）

以上就是本文的全部内容，希望对大家的学习有所帮助，也希望大家多多支持。

python,傅里叶级数展开

DDR爱好者之家 Design By 杰米

广告合作：本站广告合作请联系QQ：858582 申请时备注：广告合作（否则不回）
免责声明：本站资源来自互联网收集,仅供用于学习和交流,请遵循相关法律法规,本站一切资源不代表本站立场,如有侵权、后门、不妥请联系本站删除！

DDR爱好者之家 Design By 杰米

评论“python实现傅里叶级数展开的实现”

暂无评论...

python实现傅里叶级数展开的实现

Flask框架配置与调试操作示例

python实现时间o(1)的最小栈的实例代码

评论“python实现傅里叶级数展开的实现”

更新日志

友情链接