%用Burg法計(jì)算AR模型的參數(shù)！
%下面是一個(gè)例子：大家如果遇見(jiàn)其他的題的話，
%可以直接把濾波器的階次K和數(shù)據(jù)矢量改變一下，就可以求出參數(shù)！
%非常方便！
K=2; %有待輸入
x=[1 2 3 4 5];%有待輸入

N=max(size(x));
r=[];
%初始化。
ef(1,:)=x;
eb(1,:)=x;
sigma(1)=sum(x.^2)/N;
a=zeros(K+1,K+1);
a(:,1)=1;
for p=2:K+1;
    r(p-1)=2*(ef(p-1,[p:N])*eb(p-1,[p-1:N-1])')/(ef(p-1,[p:N])*ef(p-1,[p:N])'+eb(p-1,[p-1:N-1])*eb(p-1,[p-1:N-1])');
    sigma(p)=(1-r(p-1).^2)*sigma(p-1);
     
    for i=2:p-1
        a(p,i)=a(p-1,i)-r(p-1)*a(p-1,p-i+1);        
    end
    a(p,p)=-r(p-1); 
    for n=2:N
        ef(p,n)=ef(p-1,n)-r(p-1)*eb(p-1,n-1);
        eb(p,n)=-r(p-1)*ef(p-1,n)+eb(p-1,n-1);
    end    
end
r
a=a([2:K+1],[2:K+1])
sigma