依此类推，经过k-1步后，可得新矩阵A(): .(1) (1) 01k A(k)= () () akk akn (2.1) a限 恩 第k步：消第k列， () 设Q手0，开算：资三大+L候次将4的的第i行>)减去第至 的1：倍，将新得到的矩阵记为Ak+1),矩阵元素的更新公式为 ag+=a9-au周，6+=b-lb,i,j=k+1,k+2,n. (2.2) http://math.ecnu.edu.cn/-jypan 9/52 依此类推, 经过 k − 1 步后, 可得新矩阵 A(k) : A (k) =           a (1) 11 · · · a (1) 1k · · · a (1) 1n b (1) 1 . . . . . . . . . . . . a (k) kk · · · a (k) kn b (k) k . . . . . . . . . . . . a (k) nk · · · a (k) nn b (k) n           , (2.1) 第 k 步: 消第 k 列. 设 a (k) kk ̸= 0 , 计算 lik = a (k) ik a (k) kk , i = k + 1, . . . , n. 依次将 A(k) 的第 i 行 (i > k) 减去第 k 行 的 lik 倍, 将新得到的矩阵记为 A(k+1) , 矩阵元素的更新公式为 a (k+1) ij = a (k) ij − lika (k) kj , b(k+1) i = b (k) i − likb (k) k , i, j = k + 1, k + 2, . . . , n. (2.2) http://math.ecnu.edu.cn/~jypan 9/52