AR模型功率谱估计burg算法

clear all;
u=[1018266353172092154125856838231024831321311189067604741211664714344543484228108201512143546413024167428173650626771482881357122138103866337241115406298124966664543921742355949677594447301673774];
L=length（u）;
p=5;
K=zeros（1p）;
var=zeros（1p）;
e=zeros（pL）;
b=zeros（pL）;
f=0;
a=zeros（pp）;
e（1:）=u;
b（1:）=u;
ejw=zeros（pp）;
for i=1:L;
    f=f+u（i）*u（i）;
end
var（1）=f*2; 
for k=1:p
    fz=0;fm=0;
    for n=k:L-1
        fz=fz+e（kn+1）*b（kn）;
        fm=fm+e（kn+1）^2+b（kn）^2;
    end
    K（k）=2*fz/fm;
    a（kk）=-K（k）;
    for i=1:k
        if i            a（ki）=a（k-1i）-K（k）*a（k-1k-i）;
        end
    end
    for i=1:L
        if（i==1）