import numpy as np import matplotlib.pyplot as plt from scipy.fftpack import fft fs=100 #采样频率 N=128 #数据点数 n=np.arange(0,N) t=n/fs #时间序列 pi=3.14 x=0.5*np.sin(2*pi*15*t)+2*np.sin(2*pi*40*t) y=np.abs(fft(x)) #fft变换后的振幅 f=n*fs/N #频率序列 print(t.shape) plt.subplot(211) plt.plot(n,x) plt.subplot(212) plt.plot(n,y) plt.show()