matlab 的 AR 模型参数谱估计,算是信号里一个挺基础但蛮有用的工具了。主要是用yule-walker方程配上levinson-durbin算法搞定参数估计,整个过程也不复杂,尤其用 Matlab 现成的函数,基本上几行代码就能跑起来。

AR 建模的套路挺适合做功率谱估计的,像你要某段时间序列的频率成分,用这个方法还挺方便的。yule-walker那套思路本身就比较稳,加上levinson递推,效率也不错,是你不想自己推矩阵的时候,直接调用aryule这些函数,能省不少事。

想再深入了解的话,可以看看这几个链接:

哦对了,如果你还在手动写 AR 方程的矩阵推导,不如直接上 Matlab 的内置函数,省时省力。功率谱估计常用在金融、通信、甚至生物信号里,用 Matlab 搞挺稳的。

如果你是刚接触谱估计,可以先跑跑 Matlab 的aryulepyulear,感受下不同阶数对结果的影响,理解起来也更直观些。