列主元消去法解方程组实验报告.docVIP

实验名称： 列主元消去法解方程组 引言 我们知道，高斯消去法是一个古老的解线性方程组的方法。而在用高斯消去法解Ax=b时，其中设A为非奇异矩阵，可能出现的情况，这时必须进行带行交换的高斯消去法。但在实际计算中即使但其绝对值很小时，用作除数，会导致中间结果矩阵元素数量级严重增长和舍入误差的扩散，使得最后的结果不可靠。因此，小主元可能导致计算的失败，我们应该避免采用绝对值很小的主元素。为此，我们在高斯消去法的每一步应该在系数矩阵或消元后的低阶矩阵中选取绝对值最大的元素作为主元素，保持乘数，以便减少计算过程中舍入误差对计算解的影响。 一种方式是完全主元消去法，这种消去法是在每次选主元时，选择为主元素。这种方法是解低阶稠密矩阵方程组的有效方法，但这种方法在选取主元时要花费一定的计算机时间。实际计算中我们常采用部分选主元的的消去法。列主元消去法即在每次选主元时，仅依次按列选取绝对值最大的元素作为主元素，且仅交换两行，再进行消元计算。 实验目的和要求 运用matlab编写一个.m文件，要求用列主元消去法求解方程组（实现PA=LU）： 要求输出以下内容： 计算解x； L,U； 整形数组IP（i）（i=1,2，…，n-1）（记录主行信息） 算法原理与流程图 算法原理 设有线性方程组Ax=b，A为非奇异矩阵。方程组的增广矩阵为 第1步（k=1）：A的第一列中选取绝