清华大学数值分析实验报告.docx

清华大学数值分析实验报告

数值分析实验报告一、实验3.1题目：考虑线性方程组，，编制一个能自动选取主元，又能手动选取主元的求解线性代数方程组的Gauss消去过程。（1）取矩阵，则方程有解。取计算矩阵的条件数。分别用顺序Gauss消元、列主元Gauss消元和完全选主元Gauss消元方法求解，结果如何？（2）现选择程序中手动选取主元的功能，每步消去过程都选取模最小或按模尽可能小的元素作为主元进行消元，观察并记录计算结果，若每步消去过程总选取按模最大的元素作为主元，结果又如何？分析实验的结果。（3）取矩阵阶数n=20或者更大，重复上述实验过程，观察记录并分析不同的问题及消去过程中选择不同的主元时计算结果的差异，说明主元素的选取在消去过程中的作用。（4）选取其他你感兴趣的问题或者随机生成的矩阵，计算其条件数，重复上述实验，观察记录并分析实验的结果。1.算法介绍首先，分析各种算法消去过程的计算公式，顺序高斯消去法：第k步消去中，设增广矩阵中的元素（若等于零则可以判定系数矩阵为奇异矩阵，停止计算），则对k行以下各行计算，分别用乘以增广矩阵的第行并加到第行，则可将增广矩阵中第列中以下的元素消为零；重复此方法，从第1步进行到第n-1步，则可以得到最终的增广矩阵，即；列主元高斯消去法：第k步消去中，在增广矩阵中的子方阵中，选取使得，当时，对中第行与第行交换，然后按照和顺序消去法相同的步骤进行。重复此方法，从第1步进行第n-1步，就可以得到最终的增广矩阵，即；完全主元高斯消去法：第k步消去中，在增广矩阵中对应的子方阵中，选取使得，若或，则对中第行与第行、第列与第列交换，然后按照和顺序消去法相同的步骤进行即可。重复此方法，从第1步进行到第n-1步，就可以得到最终的增广矩阵，即；接下来，分析回代过程求解的公式，容易看出，对上述任一种消元法，均有以下计算公式：2.实验程序的设计一、输入实验要求及初始条件；二、计算系数矩阵A的条件数及方程组的理论解；三、对各不同方法编程计算，并输出最终计算结果。3.计算结果及分析（1）先计算系数矩阵的条件数，结果如下，可知系数矩阵的条件数较大，故此问题属于病态问题，b或A的扰动都可能引起解的较大误差；采用顺序高斯消去法，计算结果为：最终解为x=(1.***,1.***,1.***,1.***,0.***,1.***,0.***,1.***,0.***,1.***)T使用无穷范数衡量误差，得到=2.842***1e-14，可以发现，采用顺序高斯消元法求得的解与精确解之间误差较小。通过进一步观察，可以发现，按照顺序高斯消去法计算时，其选取的主元值和矩阵中其他元素大小相近，因此顺序高斯消去法方式并没有对结果造成特别大的影响。若采用列主元高斯消元法，则结果为：最终解为x=(1.***,1.***,1.***,1.***,1.***,1.***,1.***,1.***,1.***,1.***)T同样使用无穷范数衡量误差，有=0；若使用完全主元高斯消元法，则结果为最终解x=(1.***,1.***,1.***,1.***,1.***,1.***,1.***,1.***,1.***,1.***)T同样使用无穷范数衡量误差，有=0；（2）若每步都选取模最小或尽可能小的元素为主元，则计算结果为最终解x=(1.***1.***1.***1.***0.***1.***0.***1.***0.***1.***)T使用无穷范数衡量误差，有为2.842***1e-14；而完全主元消去法的误差为=0。从（1）和（2）的实验结果可以发现，列主元消去法和完全主元消去法都得到了精确解，而顺序高斯消去法和以模尽量小的元素为主元的消去法没有得到精确解。在后两种消去法中，由于程序计算时的舍入误差，对最终结果产生了一定的影响，但由于方程组的维度较低，并且元素之间相差不大，所以误差仍比较小。为进一步分析，计算上述4种方法每步选取的主元数值，并列表进行比较，结果如下：第n次消元顺序列主元完全主元模最小6.***6.***4.***4.***4.***4.***4.***33334.***33334.***4.***4.***4.***4.0***0634.0***0634.***4.***4.0039***4.0039***4.***0.0***4690.0***4694.***从上表可以发现，对这个方程组而言，顺序高斯消去选取的主元恰好事模尽量小的元素，而由于列主元和完全主元选取的元素为8，与4在数量级上差别小，所以计算过程中的累积误差也较小，最终4种方法的输出结果均较为精确。在这里，具体解释一下顺序法与模最小法的计算结果完全一致的原因。该矩阵在消元过程中，每次选取主元的一列只有两个非零元素，对角线上的元素为4左右，而其正下方的元素为8，该列其余位置的元素均为0。