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

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

验报告材料

本页仅作为文档页封面，使用时可以删除

This document is for reference only-rar21year.March

数值分析实验报告

一、 实验 题目：

考虑线性方程组b Ax =，n n R A ⨯∈，n R b ∈，编制一个能自动选取主元，又能手动选取主元的求解线性代数方程组的Gauss 消去过程。

（1）取矩阵⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢

⎢⎢⎣⎡=6816816816 A ，⎥⎥

⎥

⎥⎥

⎥⎦⎤

⎢⎢⎢⎢⎢⎢⎣⎡=1415157 b ，则方程有解()T x 1,,1,1*⋯=。取10=n 计算矩阵的条件数。分别用顺序Gauss 消元、列主元Gauss 消元和完全选主元Gauss 消元方法求解，结果如何？

（2）现选择程序中手动选取主元的功能，每步消去过程都选取模最小或按模尽可能小的元素作为主元进行消元，观察并记录计算结果，若每步消去过程总选取按模最大的元素作为主元，结果又如何分析实验的结果。

（3）取矩阵阶数n=20或者更大，重复上述实验过程，观察记录并分析不同的问题及消去过程中选择不同的主元时计算结果的差异，说明主元素的选取在消去过程中的作用。

（4）选取其他你感兴趣的问题或者随机生成的矩阵，计算其条件数，重复上述实验，观察记录并分析实验的结果。

1. 算法介绍

首先，分析各种算法消去过程的计算公式， 顺序高斯消去法：

第k 步消去中，设增广矩阵B 中的元素()

0k kk a ≠（若等于零则可以判定系数

矩阵为奇异矩阵，停止计算），则对k 行以下各行计算

()()

,1,2,

,k ik

ik k kk

a l i k k n a =

=++，分别用ik l -乘以增广矩阵B 的第k 行并加到第

1,2,

,k k n ++行，则可将增广矩阵B 中第k 列中()

k

kk

a 以下的元素消为零；重复此方法，从第1步进行到第n-1步，则可以得到最终的增广矩阵，即

()()(),n n n B A b ⎡⎤=⎣

⎦； 列主元高斯消去法：

第k 步消去中，在增广矩阵B 中的子方阵()()()()k k

kk

kn

k k nk

nn a a a a ⎡⎤⎢

⎥⎢

⎥⎢⎥⎣⎦

中，选取()

k k i k a 使得()(k)

max k k i k ik k i n

a a ≤≤=，当k i k ≠时，对B 中第k 行与第k i 行交换，然后按照和顺序消

去法相同的步骤进行。重复此方法，从第1步进行第n-1步，就可以得到最终的增广矩阵，即(

)

()()111,n n n B A b ⎡⎤=⎣

⎦； 完全主元高斯消去法：

第k 步消去中，在增广矩阵B 中对应的子方阵()()()()k k

kk

kn

k k nk

nn a a a a ⎡⎤⎢

⎥⎢

⎥⎢⎥⎣⎦

中，选取()k k k i j a 使得()(k)

max k k k i j ij k i n

k j n

a a ≤≤≤≤=，若k i k ≠或k j k ≠，则对B 中第k 行与第k i 行、第k 列与第

k j 列交换，然后按照和顺序消去法相同的步骤进行即可。重复此方法，从第1

步进行到第n-1步，就可以得到最终的增广矩阵，即()

()()222,n n n B

A b ⎡⎤=⎣

⎦；

接下来，分析回代过程求解的公式，容易看出，对上述任一种消元法，均有以下计算公式：

()()()

()()1;

;1,2,,1

n

n

n n nn

n

n n k kj j

j k k n kk

b x a b a x x k n n a

=+=-=

=--∑

2. 实验程序的设计

一、输入实验要求及初始条件；

二、计算系数矩阵A 的条件数及方程组的理论解； 三、对各不同方法编程计算，并输出最终计算结果。

3. 计算结果及分析

（1）

先计算系数矩阵的条件数，结果如下，

12()2557.500000,()1727.556025,()2557.50000cond A cond A cond A ∞===

可知系数矩阵的条件数较大，故此问题属于病态问题， b 或A 的扰动都可能引起解的较大误差；

采用顺序高斯消去法，计算结果为：

最终解为x=, , , , , , , , , T

使用无穷范数衡量误差，得到*X X ε∞

=-=，可以发现，采用顺序高斯消

元法求得的解与精确解之间误差较小。通过进一步观察，可以发现，按照顺序高斯消去法计算时，其选取的主元值和矩阵中其他元素大小相近，因此顺序高斯消去法方式并没有对结果造成特别大的影响。 若采用列主元高斯消元法，则结果为：

最终解为x=, , , , , , , , , T

同样使用无穷范数衡量误差，有*X X ε∞

=-=0；

若使用完全主元高斯消元法，则结果为

最终解x=, , , , , , , , , T

同样使用无穷范数衡量误差，有*X X ε∞

=-=0；

（2）

若每步都选取模最小或尽可能小的元素为主元，则计算结果为

最终解x= T