嗨,我有一个快速傅立叶变换,这是相当嘈杂。如何应用到我的代码汉明窗口,使它不那么嘈杂。看看我的代码:
plt.subplot(212)
plt.title('Fast Fourier Transform')
plt.ylabel('Power [a.u.]')
plt.xlabel('Frequency Hz')
fft1 = (Bx[51:-14])
fft2 = (By[1:-14])
for dataset in [fft1]:
dataset = np.asarray(dataset)
psd = np.abs(np.fft.fft(dataset))**2.5
freq = np.fft.fftfreq(dataset.size, float(300)/dataset.size)
plt.semilogy(freq[freq>0], psd[freq>0]/dataset.size**2, color='r')
for dataset2 in [fft2]:
dataset2 = np.asarray(dataset2)
psd2 = np.abs(np.fft.fft(dataset2))**2.5
freq2 = np.fft.fftfreq(dataset2.size, float(300)/dataset2.size)
plt.semilogy(freq2[freq2>0], psd2[freq2>0]/dataset2.size**2, color='b')什么是plt.show()

我需要的是:

我已经看过(https://docs.scipy.org/doc/scipy-0.13.0/reference/generated/scipy.signal.hamming.html)和这个(https://docs.scipy.org/doc/numpy-1.13.0/reference/generated/numpy.hamming.html),但仍然不知道如何将它应用到我的代码中。有什么想法吗?正如我所说的,你在第二张图片中看到了我需要的东西。也许Blackman窗口也可以很好地应用,但仍然不知道如何添加它。
应用以下内容:
freqs, psd = scipy.signal.welch(dataset, fs=300, window='hamming')给我的,它看起来不像我想要的图表。

https://stackoverflow.com/questions/47613183
复制相似问题