3.Z0=4,P(0)=946,P(0)’=1808;Z1=Z0-P(0)/P(0)’=3.4768;
代码:
a=[2 -8 1 -5 -7 3];
b=0;
af=a(end);
z0=4;
m=length(a)-1;
for k=m:-1:1
b=af+z0*b;
B(k)=b;
af=a(k)+z0*af;
A(k)=af;
end
B0=B(1);
A0=A(1);
Z1=z0-A0/B0
4.U=3,V=1,detU=-1.6518;detV=17.0324;
代码:
%%
a=[2 -8 1 -5 -7 3];
u=3;v=1;
b(6)=a(6);
c(6)=0;
c(5)=a(4);
n=length(a);
b(5)=a(5)+u*b(4);
for k=n-2:-1:1
b(k)=a(k)+u*b(k+1)+v*b(k+2);
c(k)=b(k+1)+u*c(k+1)+v*c(k+2);
end
J=c(1)*c(3)-c(2).^2
u=u+(c(2)*b(2)-c(3)*b(1))./J
v=v+(c(2)*b(1)-c(1)*b(2))./J