1.题目
书本592页14.12 利用习题14.11的实验数据(数据文件test.mat在https://github.com/jianjake/-test):
(1) 用自相关求解AR模型系数以估计其功率谱,模型阶次p=8,p=11,p=14,自己可调整。
(2) 用Burg方法重复(1)
2. 分析
自相关估计法(Yule-Walker AR方法),描述自回归序列参数与其协方差函数之间关系的方程,通过计算p+1阶自相关系数的Yule-Walker方程计算模型参数,先定义自相关函数为
对于真实信号,由自相关函数的偶对称性质可得:
因此上式可以写为:
接下来便可以使用Levinson-Durbin迭代算法计算AR模型的p+1个参数:
aryule进行参数估计,pyulear求出估计功率谱。
BurgBurg算法是首先令前后项预测误差功率之和为最小估计出反射系数ǩ(m ),