Levinson-Duin算法仿真与第七章学习参考.pdfVIP

文档给出的三个仿真作为第七章学习参考,与PPt内容结合

Levinson-Durbin算法

1.算法程序

function[a,k,Po]=Durbin(r,M)Po(2)=P;

a=zeros(M+1,1);%Recursion

k=zeros(M,1);form=2:M

Po=zeros(M+1,1);

a=(r(2:m))*flipud((a(1:m-1)))+r(m+1);

%Initializationk(m)=-a/P;

P=r(1);a(1:m)=[(a(1:m-1))0]...

Po(1)=P;+[(flipud(a(1:m-1)))1]*k(m);

a=r(2);P=P+a*k(m);

k(1)=-a/P;Po(m+1)=P;

a(1)=k(1);end

P=P+a*k(1);

2.例7.4.2P357

clear;clc;disp(输出-k);

M=3;%fori=2:M

r=[3;2;1;0.5];%disp(k(i));

[a,k,Po]=Durbin(r,M);end

disp(输出模型参数-a);disp(输出MMSE-P);

fori=1:Mfori=1:M+1

disp(a(i