zuisu.m

数值线性代数的Matlab应用程序包 共13个程序函数

function [xn,n]=zuisu(A,b,d)
[m,n]=size(A);
x(:,1)=zeros(m,1);
r(:,1)=b;
n=2;
while n>0
   t1(n)=(r(:,n-1)'*r(:,n-1))/(r(:,n-1)'*A*r(:,n-1));
   x(:,n)=x(:,n-1)+t1(n)*r(:,n-1);
   r(:,n)=r(:,n-1)-t1(n)*A*r(:,n-1);
   if r(:,n)'*r(:,n)<d
       break;    
   end
   n=n+1;
end
xn=x(:,n);
n=n-1