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

. . . 数值分析实验报告 一、 实验 3.1 题目： 考虑线性程组 Ax b ， A Rn n ， b R n ，编制一个能自动选取主元，又能手动 选取主元的求解线性代数程组的 Gauss消去过程。 6 1 7 8 6 1 15 * T A b x 1,1, ,1 （1）取矩阵 ， ，则程有解 。取 n 10 8 6 1 15 8 6 14 计算矩阵的条件数。分别用顺序 Gauss 消元、列主元 Gauss 消元和完全选主元 Gauss消元法求解，结果如？ （2 ）现选择程序中手动选取主元的功能，每步消去过程都选取模最小或按模尽 可能小的元素作为主元进行消元， 观察并记录计算结果， 若每步消去过程总选取 按模最大的元素作为主元，结果又如？分析实验的结果。 （3 ）取矩阵阶数 n=20 或者更大，重复上述实验过程，观察记录并分析不同的 问题及消去过程中选择不同的主元时计算结果的差异， 说明主元素的选取在消去 过程中的作用。 （4 ）选取其他你感兴趣的问题或者随机生成的矩阵，计算其条件数，重复上述 实验，观察记录并分析实验的结果。 z . . . z . . . 1. 算法介绍 首先，分析各种算法消去过程的计算公式， 顺序高斯消去法： k 第 k 步消去中，设增广矩阵 B 中的元素 akk 0 （若等于零则可以判定系数 k aik 矩阵为奇异矩阵，停止计算），则对 k 行以下各行计算 lik k ,i k 1,k 2, , n ， a