第3章:线性方程组求解 3/5 for(j=Kj<Nj++)A[[-=X*AK[j; 4利用约当消去法求矩阵的逆 对于MxM可逆矩阵A我们可以构造矩阵T=(AD,那么用约当消 去法把T的左边M列化为单位矩阵后其效果相当于用A1左乘 T从而得到A1(AD)=(A)所以T的后面的M列就是我们所要 求的A1 4实现选主元的源代码 j=K x=Math.abs(A[KIkd) for(=K+1;i<M++) i if(Math. abs(au[k)<x) continue x=Math.abs(Au][k]) 」=第 3 章：线性方程组求解 3/5 for(j=K;j<N;j++)A[i][j] -=x*A[K][j]; } } 4.利用约当消去法求矩阵的逆 对于 M×M 可逆矩阵A,我们可以构造矩阵 T=(A I),那么用约当消 去法把 T 的左边 M 列化为单位矩阵后，其效果相当于用 A-1左乘 T,从而得到 A-1 (A I)=(I A-1 ),所以 T 的后面的 M 列就是我们所要 求的 A-1 4 实现选主元的源代码 j=K; x=Math.abs(A[K][K]); for(i=K+1;i<M;i++) { if(Math.abs(A[i][K])<x) continue; x=Math.abs(A[i][K]); j=i; }