levinson_durbin算法_levinsdn_durbin算法c++代码

// An highlighted block
function [a_p,E_V]=Levinson_Durbin_Algo (x,p)
%输入信号的互相关函数计算算
    N=length (x) ;
    for i=1:N
        Rx ( i)=x(1:N-i+1)*(x( i :N))'/N ;%计算输入信号的自相关函数值数值
       end
%levinson算法迭代过程
    a(1)=1;
    a(2)=-Rx(2)/Rx(1);
%自适应滤波器的前两个系数个系数

    for k=1:p-1
        E_V(k+1)=Rx(1)+a(2:k+1)*Rx(2:k+1)';
%预测误差更新
        D(k+1)=a(1:k+1)*(fliplr(Rx(2:k+2)));%为扩大方程矩阵中的Dk的更新新
        gama(k+2 )=D(k+1)/E_V(k+1);%反射系数gama的更新的更新
        E_V(k+2)-E_V(k+1)*(1-(gama(k+2))^2);
%预测误差的更新预测系数的更新
        a_yuce (1)=1;
%由k阶系数预测第k+l阶系数
        for q=1:k
            a_yuce(q+1)=a(q+1)-gama(k+2)*a(k-q+2);
        end
    a_yuce (k+2)=-gama ( k+2) ;
    a=a_yuce ;
%准备下一次迭代
    end
    a_p=a;
%算法迭代完成输出模型系数
    E_p=E_V;
%输出系统预测误差
end



1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36