Python 3.x 为什么matplotlib震级谱函数似乎显示错误的震级？

我创建了一个由两个正弦波组成的1秒长的音频样本，然后使用matplotlibs幅值谱函数绘制频谱，结果似乎是错误的。在一秒钟的音频样本中，这两个波的振幅完全相同，但幅度却大不相同。这对我来说似乎很奇怪，所以我也使用了numpys函数来绘制DFT，其大小与我认为应该的完全相同。结果图如下图所示。有人知道为什么会这样吗？我的代码有没有做错什么？任何帮助都将是非常感激的

最简单的工作示例：

import matplotlib.pyplot as plt
import numpy as np

sr = 20000
freq1 = 200
freq2 = 100
duration = 1

x = np.linspace(0, duration, sr * duration)
y = np.concatenate([0.5*np.sin(freq1 * 2 * np.pi * x[:10000]) + 0.5*np.sin(freq2 * 2 * np.pi * x[:10000]), np.sin(freq1 * 2 * np.pi * x[10000:15000]), np.sin(freq2 * 2 * np.pi * x[15000:20000])])

fig, ax = plt.subplots(3, 1, figsize=(12, 10))

ax[0].plot(x, y)
ax[0].axis(xmin=0, xmax=1)
ax[0].set_xlabel('Time [s]')
ax[0].set_ylabel('Amplitude [-]')

ax[1].magnitude_spectrum(y, Fs=sr, color='C1')
ax[1].axis(xmin=0, xmax=500)
ax[1].set_xlabel('Frequency [Hz]')
ax[1].set_ylabel('Magnitude [-]')

ax[2].plot(np.fft.rfftfreq(sr, d=1/sr), np.abs(np.fft.rfft(y, norm='ortho'))/100)
ax[2].axis(xmin=0, xmax=500)
ax[2].set_xlabel('Frequency [Hz]')
ax[2].set_ylabel('Magnitude [-]')

plt.tight_layout()
plt.show()

---

我认为这与matplotlib中使用的窗口有关。默认情况下，它使用汉宁窗口，因此将“窗口类型”更改为“窗口无”。在这两种情况下，缩放的方式也不同。通过执行以下更改，您将看到它们都匹配

from matplotlib import mlab
ax[1].magnitude_spectrum(y, Fs=sr, color='C1', window=mlab.window_none)

ax[2].plot(np.fft.rfftfreq(sr, d=1/sr), np.abs(np.fft.rfft(y))/sr)

---

结果

这正是我想要的，你甚至解释了我不理解的比例，完美。非常感谢你！