仅作为记录,大量请跳过。
信号长度为1s时:
对原信号做傅里叶变换(不做频率轴x轴的处理,仅通过python图像的放大工具观察)
待无限放大后:
发现400Hz和400.2Hz的频谱无法分开
直接上代码:
import numpy as np
import matplotlib.pyplot as plt
fs=2000
t=np.arange(0,1,1/fs)
a=np.cos(2*np.pi*400*t)+np.cos(2*np.pi*400.2*t);data=a
N=len(data)
f=fs*np.arange(0,fs/2)/fs
fy=np.fft.fft(data,N)
fy=np.abs(fy)
# fy=fy[0:len(f)]
# fy=20*np.log(fy)
plt.figure()
plt.plot(fy)
plt.show()
print('ok')
10s的信号:
以上面同样的放大过程,可看到
频谱中有两条,即400Hz和400.2Hz的谱
直接上代码:
import numpy as np
import matplotlib.pyplot as plt
fs=2000
t=np.arange(0,10,1/fs)
a=np.cos(2*np.pi*400*t)+np.cos(2*np.pi*400.2*t);data=a
N=len(data)
f=fs*np.arange(0,fs/2)/fs
fy=np.fft.fft(data,N)
fy=np.abs(fy)
# fy=fy[0:len(f)]
# fy=20*np.log(fy)
plt.figure()
plt.plot(fy)
plt.show()
print('ok')