import numpy as np from scipy.fftpack import fft,ifft from matplotlib.pylab import plt# Fs = 1000 #采樣頻率 # T = 1/Fs #采樣周期，只相鄰兩數據點的時間間隔 # L = 1500 #信號長度Fs = 5000000000 #采樣頻率 T = 1/Fs #采樣周期，只相鄰兩數據點的時間間隔 L = 10000#信號長度nnn = 500000t = np.arange(L)*T #間隔 T * 個數 npp = np.piX = np.sin(2*np.pi *nnn*t) # X = S + np.random.rand(L)plt.plot(t[:nnn], X[:nnn]) plt.xlabel("Time(s)") plt.ylabel("Amplitude") plt.title("Signol with random noise") plt.show()Y = fft(X) p2 = np.abs(Y) # 雙側頻譜 p1 = p2[:int(L/2)] f = np.arange(int(L/2))*Fs/L; # plt.plot(f,2*p1/L) plt.plot(f, p1) plt.title('Single-Sided Amplitude Spectrum of X(t)') plt.xlabel('f (Hz)') plt.ylabel('|P1(f)|') plt.show()

總結

以上是生活随笔為你收集整理的正弦波及FFT变换py程序的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。