你的理论传递函数G (z) = 1 / (1 + z z ^ ^ 2 + 0.5 * 4)。的频率响应函数的值是z = exp (j * w * Ts)。Ts = 2000。w是你想计算频谱的频率(说,linspace(0,π/ Ts 1000))
计算时间序列模型的光谱th_fb等,使用freqresp或预示命令。然而,请注意,基于“增大化现实”技术的预计模型平稳时间序列(或确定性情况下的脉冲响应)。所以我将应用AR diff (diff (Y))的各种模型;您可能需要设置noisevariance 1后估计进行比较。
使用FFT结果:
| G (w)
| =
fft (Y)
。/ | fft (x) |这不会不好运气因为各种各样的原因,除了非周期输入信号的性质。