矩形窗的DFT

Question

我正在尝试计算矩形窗口的DFT的大小和角度。下面的代码得到了我期望的大小，但阶段看起来是错误的。有人能帮我找出我做错了什么吗？谢谢!

这是我得到的：

​

​

这是我正在读的书(理查德·莱昂斯著)的内容：

​

​

fs = 256  # sample freq - number of samples per second
k  = 11
window = np.array([0] * 64 + [1] * k + [0] * 64)
freq = fft(window)

plt.stem(fftfreq(len(freq), 1/fs), np.abs(freq), 'r')
plt.show()
plt.stem(fftfreq(len(freq), 1/fs), np.angle(freq), 'r')
plt.xlim([-30, 30])
plt.show()

Answer 1

https://stackoverflow.com/a/60764461/89806解决了我的问题。在调用fft之前，我只需要添加fftshift(window)。