numpy中的fft和scipy中的fft,fftshift以及fftfreq
阿新 • • 發佈:2019-01-26
numpy中有一個fft的庫,scipy中也有一個fftpack的庫,各自都有fft函式,兩者的用法基本是一致的: 在畫頻譜圖的時候,要給出橫座標的數字頻率,這裡可以用fftfreq給出,對於fftfreq的說明如下:scipy.fftpack.fftfreq(n, d=1.0)第一個引數n是FFT的點數,一般取FFT之後的資料的長度(size)第二個引數d是取樣週期,其倒數就是取樣頻率Fs,即d=1/Fs需要說明的是,DFT變換中,頻率的解析度為Fs/n=1/d*nfftfreq得到的結果為各個數字頻率 k*Fs/n = k/d*n最終fftfreq輸出的數字頻率的結果為
舉例:
舉例如下:
fftshift也可以應用在數字頻率上,即
最後,可以利用freq_shift作為橫座標,xfn_shift的幅值mag_shift作為縱座標,作圖如下: