线性方程组矩阵求解算法.docVIP

线性方程组的矩阵求解算法 摘要 线性方程组的矩阵求解算法,只需在约当消元法的基础上,再对方程组的增广矩阵的行最简形进行行(列)删除和增加行,交换行等运算即可得到方程组的解,并且这种方法既可求解有唯一解的方程组.因而算法简单,易于实现. 关键词 线性方程组;解向量;解法;约当消元法 矩阵求解算法 设有线性方程组,其增广矩阵,算法的步骤如下: 第一步:利用约当消元法,把增广矩阵化为行最简形,设行最简形为.若则方程组无解;否则设并执行以下步骤; 第二步:删除中的所有零行和每一行第一个非零元素(这个非零元素一定是1)所在的列,得到矩阵并记录每行的第一个非零元所在的列标,放在一维数组中,如第行的第一个非零元在第列,则; 第三步:构造矩阵,其中 第四步:对矩阵中的行作交换运算:把中的第行(即从第行开始直到第一行)依次与其下一行交换,使之成为第行,交换运算结果后的矩阵记为,则中的前个维列向量即为方程组的一个基础解系,最后一列向量即为方程组的一个特解; 第五步:写出方程组的通解. 算法证明 先证一个特殊情形,增广矩阵的行最简形矩阵的左上角为一阶的单位矩阵,即第行的第一个非零元的列标为,即,所以设为 则 由上述算法可得为 由于,故从得到时,中的行不需交换位置,