在Python中,使用numpy庫進行傅里葉變換非常簡單
import numpy as np
# 創建一個簡單的信號
time = np.linspace(0, 1, 1000)
signal = np.sin(2 * np.pi * 5 * time)
# 對信號進行傅里葉變換
fourier_transform = np.fft.fft(signal)
# 計算頻率軸
freqs = np.fft.fftfreq(len(signal), d=1/len(signal))
# 輸出結果
print("Fourier Transform:", fourier_transform)
print("Frequency Axis:", freqs)
在這個例子中,我們首先導入了numpy庫并創建了一個簡單的正弦波信號。然后,我們使用np.fft.fft()
函數對信號進行傅里葉變換。最后,我們計算了頻率軸并使用print()
函數輸出了結果。
注意:np.fft.fft()
函數返回的結果是一個復數數組,它包含了信號的實部和虛部。如果你只關心實部,可以使用np.real()
函數提取出來。同樣,如果你只關心虛部,可以使用np.imag()
函數提取出來。