您好,登錄后才能下訂單哦!
小編給大家分享一下如何使用python實現(xiàn)傅里葉級數(shù)展開,希望大家閱讀完這篇文章之后都有所收獲,下面讓我們一起去探討吧!
傅立葉級數(shù)的介紹我就不說了,自己也是應(yīng)用為主,之前一直覺得很難懂,但最近通過自己編程實現(xiàn)了一些函數(shù)的傅立葉級數(shù)展開之后對傅立葉 級數(shù)展開的概念比較清楚了
(1)函數(shù)如下
函數(shù)圖象如下:
代碼:
from pylab import * x = mgrid[-10:10:0.02] # 這里類似于MATLAB用冒號產(chǎn)生步長為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()
結(jié)果如下:
(2)
轉(zhuǎn)換代碼:
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()
結(jié)果如下:
看完了這篇文章,相信你對“如何使用python實現(xiàn)傅里葉級數(shù)展開”有了一定的了解,如果想了解更多相關(guān)知識,歡迎關(guān)注億速云行業(yè)資訊頻道,感謝各位的閱讀!
免責(zé)聲明:本站發(fā)布的內(nèi)容(圖片、視頻和文字)以原創(chuàng)、轉(zhuǎn)載和分享為主,文章觀點不代表本網(wǎng)站立場,如果涉及侵權(quán)請聯(lián)系站長郵箱:is@yisu.com進行舉報,并提供相關(guān)證據(jù),一經(jīng)查實,將立刻刪除涉嫌侵權(quán)內(nèi)容。