在实现维纳滤波器和预测器的时候,需要计算数据的自相关矩阵的逆。但是当数据量比较大的时候,计算矩阵的逆花费的代价比较大,所以需要使用Levinson-Durbin算法来实现系数的求解。
一、数据模型
k阶前项维纳预测器:

对上述模型进行一下调整,可以得到

所以a(0)=1,

。
注意:在Levinson-Durbin算法中,求解的是a(i),不是w(i)。
二、Levinson-Durbin迭代算法的实现步骤
以下是迭代算法的步骤,其中m是预测器的阶数。





初始值:
其中r(i)是数据的自相关矩阵。