在命令行中绘制噪声频谱

要绘制输入-输出模型的频谱或时间序列模型的输出频谱，请使用 spectrum。若要自定义此类图，或通过编程方式为这些图启用置信域视图，请改用 spectrumplot。

要判断所估计的噪声模型是否足够好，可以将估计的噪声模型 H 的输出频谱与 v(t) 的估计输出频谱进行比较。要计算 v(t)（该量代表系统中的实际噪声项），请使用以下命令：

ysimulated = sim(m,data);
v = ymeasured-ysimulated;

ymeasured 即 data.y。v 是噪声项 v(t)，如 What Does a Noise Spectrum Plot Show? 所述，它对应于仿真响应 ysimulated 与实际响应 ymeasured 之间的差值。

要计算实际噪声的频率响应模型，请使用 spa：

V = spa(v);

该工具箱使用以下方程来计算实际噪声的频谱：

$Φ_{v} (ω) = \sum_{τ = - \infty}^{\infty} R_{v} (τ) e^{- i ω τ}$

协方差函数 $R_{v}$ 可用 E（表示数学期望）表示，具体如下：

$R_{v} (τ) = E v (t) v (t - τ)$

若要将参数化噪声模型 H 与实际噪声 v(t) 的（非参数化）频率响应估计值进行比较，请使用 spectrum：

spectrum(V,m)

如果噪声谱的参数估计值与非参数估计值不同，那么您可能需要采用更高阶的噪声模型。