【正文】
答案,并與使用inv(A)所得結(jié)果和運(yùn)行時(shí)間進(jìn)行比較。②令i=1,2,3……n,找出第i列中從i到n的絕對(duì)值最大元素C(k,i),如果k~=I,則兩行調(diào)換位置,并且把換位信息傳給記錄矩陣flag,同時(shí)對(duì)第i行進(jìn)行化一處理,對(duì)第i+1行到n行通過(guò)行列變換使得C[:,i]=0,最后得到C為前n*n是上三角的單位陣;③通過(guò)行列變換把矩陣C的前n*n部分化成單位陣;④把C的后n*n部分賦給B,則B就是所求的逆矩陣。B=zeros(n,n)。for i=1:n for j=1:n C(i,j)=A(i,j)。 else C(i,k)=1。%記錄行順序向量for i=1:n1 t=find(abs(C(i:n,i))==max(abs(C(i:n,i))))。 flag(i)=t。C(i,:)=C(t,:)。 %換位 end if C(i,i)==0 error(39。)。 for j=i+1:n C(j,:)=C(j,:)C(i,:)*C(j,i)。%因?yàn)檠h(huán)結(jié)構(gòu)問(wèn)題,所以C[n,:]并沒(méi)有化一,所以在最后才把最后一行化一for i=n:1:2 for j=1:i1