数值分析计算实习题 列主元高斯消去法解线性方程组.docVIP

数值分析计算实习题 第5 章 解线性方程组的直接方法 【选题】 列主元高斯消去法解线性方程组。 书上的计算实习题1、2、3都要求用列主元高斯消去法解线性方程组，所以考虑写一个普适的程序来实现。 对于线性方程组Ax=b，程序允许用户从文件读入矩阵数据或直接在屏幕输入数据。 文件输入格式要求： （1）第一行为一个整数n（2= n = 100）,表示矩阵阶数。 （2）第2 ~ n+1行为矩阵A各行列的值。 （3）第n+2 ~ n+n+2行为矩阵b各行的值。 屏幕输入：按提示输入各个数据。 输出：A、b、det(A)、列主元高斯消去计算过程、解向量x。 【算法说明】 设有线性方程组Ax=b，其中设A为非奇异矩阵。方程组的增广矩阵为 第1步（k=1）：首先在A的第一列中选取绝对值最大的元素，作为第一步的主元素: 然后交换（A，b）的第1行与第l行元素，再进行消元计算。 设列主元素消去法已经完成第1步到第k-1步的按列选主元，交换两行，消元计算得到与原方程组等价的方程组 A(k)x=b(k) 第k步计算如下： 对于k=1，2，…，n-1 （1）按列选主元：即确定t使 （2）如果t≠k，则交换[A，b]第t行与第k行元素。 （3）消元计算 消元乘数mik满足: （4）回代求解 【程序】 /* 【普适列主元消去法解线性方程组】 对于线性