levinson_durbin算法matlab程序

levinson_durbin算法matlab程序

function [ak]=levinson_durbin（RN）
%Levinson-Durbin递归算法

E=zeros（1N+1）;  %E（i）表示第i-1阶预测滤波器模型误差
k=zeros（1N）;    %k（j）表示反射系数
a=zeros（NN）;    %a（ji）表示第j阶滤波器的第i个系数
E（1）=R（1）;       %初始化参数

for j=1:1:N
    temp=0;      %初始化参数
    for i=1:1:j-1
        temp=temp+a（j-1i）*R（j-i+1）;  %更新参数
    end
    if E（j）~=0    %判断预测滤波器模型误差是否为0
                  %若不为0则更新反射系数，否者将反射系数置1.
        k（j）=-（R（j+1）+temp）/E（j）;
    else
        k（j）=1;
    end
    a（jj）=k（j）;   %第j阶模型的第j位的系数   
    for i=1:1:j-1  %更新预测滤波器模型的系数
        a（ji）=a（j-1i）+k（j）*a（j-1j-i）;
    end
    E（j+1）=E（j）*（1-abs（k（j））^2）;  %第j+1阶模型的均方预测误差
end
a=a（N:）;          %预测滤波器模型的系数

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        794  2015-05-10 00:41  levinson_durbin.m

----------- ---------  ---------- -----  ----

                  794                    1