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