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()

 

SRE实战 互联网时代守护先锋，助力企业售后服务体系运筹帷幄！一键直达领取阿里云限量特价优惠。
扫码关注我们
微信号：SRE实战
拒绝背锅 运筹帷幄