python中如何使用scipy.fftpack模塊計(jì)算快速傅里葉變換?
在編程中,快速傅里葉變換是工程中非常有價(jià)值的一類算法,它可以將時(shí)域和頻域的信號(hào)相互轉(zhuǎn)化,而在python中SciPy提供fftpack模塊,可讓用戶計(jì)算快速傅立葉變換,本文介紹python中使用scipy.fftpack模塊計(jì)算快速傅里葉變換的過(guò)程。
1、使用scipy.fftpack.fftfreq()函數(shù)將生成取樣頻率
因?yàn)楣β式Y(jié)果是對(duì)稱的,僅僅需要使用譜的正值部分來(lái)找出頻率:
In[48]:pidxs=np.where(sample_freq>0)
In[49]:freqs=sample_freq[pidxs]
In[50]:power=np.abs(sig_fft)[pidxs]
2、使用scipy.fftpack.fft()將計(jì)算快速傅里葉變換
>>>a=np.arange(1,5)
>>>a
array([1,2,3,4])
>>>b=fft(a)
>>>b
array([10.+0.j,-2.+2.j,-2.+0.j,-2.-2.j])
>>>b[1]
(-2+2j)
>>>b[1].real-2.0
>>>b[1].imag2.0
以上就是python中使用scipy.fftpack模塊計(jì)算快速傅里葉變換的過(guò)程,希望能對(duì)你有所幫助喲~更多Python學(xué)習(xí)教程請(qǐng)關(guān)注IT培訓(xùn)機(jī)構(gòu):千鋒教育。