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）
?????