资源描述:
《高斯—约旦算法用来求高阶矩阵的逆》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、clc;clearall;x=[1,3,6,7,5;1,2,3,4,11;5,3,7,9,8;8,10,13,11,19;1,12,23,45,22];N=length(x);maxrow=zeros(1,N);maxcolumn=zeros(1,N);fork=1:Nmax=x(k,k);maxrow(k)=k;maxcolumn(k)=k;forrow=k:Nforcolumn=k:Nifabs(x(row,column))>maxmax=x(row,column);maxrow(k)=row;maxcolu
2、mn(k)=column;endendend%½»»»ÐÐtemp=x(k,:);x(k,:)=x(maxrow(k),:);x(maxrow(k),:)=temp;%½»»»ÁÐtemp=x(:,k);x(:,k)=x(:,maxcolumn(k));x(:,maxcolumn(k))=temp;x(k,k)=1/x(k,k);forcolumn=1:Nifcolumn~=kx(k,column)=x(k,column)*x(k,k);endendforrow=1:Nifrow~=kforcolumn=1:Ni
3、fcolumn~=kx(row,column)=x(row,column)-x(row,k)*x(k,column);endendendendforrow=1:Nifrow~=kx(row,k)=-x(row,k)*x(k,k);endendendfork=N:-1:1temp=x(:,k);x(:,k)=x(:,maxrow(k));x(:,maxrow(k))=temp;temp=x(k,:);x(k,:)=x(maxcolumn(k),:);x(maxcolumn(k),:)=temp;endx