你能给我一些建议，用scipy在python中绘制DFT图

Question

很抱歉问你们这个， 但这是我第一次使用python和FFT。 我想使用abnormal.py整数数组绘制如下图，但到目前为止，网络搜索并不是很有用。

所以，看看内容并给我一些建议

（abnormal.npy包含3865个正整数值的数组，它们代表每天某些事件的数量）

以下是我想要绘制的图表

以下是我正在尝试的代码

import numpy as np
from scipy import fftpack
import pylab as pl
import numpy.fft as fft

data = np.load('C:\\Users\\python\\Desktop\\abnormal.npy')
freq = fft.fftfreq(len(data))
power = np.abs(freq)**2

pl.figure()
pl.plot(freq, power)
pl.xlabel('Frequency [Hz]')
pl.ylabel('plower')
pl.show()