如何在 Python 中用 Matplotlib 绘制幅度频谱?
可以采取以下步骤绘制幅度频谱:-
- 设置图形大小,并调整子图间和周围的边距。
- 获取随机种子值。
- 初始化采样间隔 dt 并找出采样频率。
- 为 t 创建随机数据点。
- 若要生成噪声,请使用 numpy 获取 nse、r、cnse 和 s
- 使用 subplots() 方法创建一个图形及其子图。
- 设置绘图标题。
- 绘制幅度频谱。
- 若要显示图形,请使用 show() 方法。
示例
import matplotlib.pyplot as plt import numpy as np plt.rcParams["figure.figsize"] = [7.50, 3.50] plt.rcParams["figure.autolayout"] = True np.random.seed(0) dt = 0.01 # sampling interval Fs = 1 / dt # sampling frequency t = np.arange(0, 10, dt) # generate noise: nse = np.random.randn(len(t)) r = np.exp(-t / 0.05) cnse = np.convolve(nse, r) * dt cnse = cnse[:len(t)] s = 0.1 * np.sin(4 * np.pi * t) + cnse fig, axs = plt.subplots() axs.set_title("Magnitude Spectrum") axs.magnitude_spectrum(s, Fs=Fs, color='C1') plt.show()
输出
广告