浅谈对矩阵理论、数值分析及matlab应用的理解 7.docVIP

附录1 Ax=b 利用高斯消元法 clear format rat A=randn(5) b=randn(5,1) [m,n]=size(A); for i=1:(m-1) numb=int2str(i); for j=(i+1):m A(j,:)=A(j,:)-A(i,:)*A(j,i)/A(i,i); end A end x(m)=A(m,n)/A(m,m); for i=(m-1):-1:1 x(i)=(A(i,n)-A(i,i+1:m)*x(i+1:m))/A(i,i); end x 2、利用高斯列主元消元法 clear format rat A=randn(5) b=randn(5,1) [m,n]=size(A); for i=1:(m-1) numb=int2str(i); disp([第,numb,次选列主元后的增广矩阵]) temp=max(abs(A(i:m,i))); [a,b]=find(abs(A(i:m,i))==temp); tempo=A(a(1)+i-1,:); A(a(1)+i-1,:)=A(i,:); A(i,:)=tempo disp([第,