小编给大家分享一下如何使用python实现傅里叶级数展开,希望大家阅读完这篇文章之后都有所收获,下面让我们一起去探讨吧!
傅立叶级数的介绍我就不说了,自己也是应用为主,之前一直觉得很难懂,但最近通过自己编程实现了一些函数的傅立叶级数展开之后对傅立叶 级数展开的概念比较清楚了
(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()
结果如下:
看完了这篇文章,相信你对“如何使用python实现傅里叶级数展开”有了一定的了解,如果想了解更多相关知识,欢迎关注天达云行业资讯频道,感谢各位的阅读!