线性方程组及MATLAB应用.pdfVIP

实用文档 数值实验 线性方程组与 MATLAB应用 王 1．实验目的：理解矩阵的范数与条件数。 实验内容：已知矩阵 1 1 1 1 1 1 1 1 A 求 A 1 ， A 2 ， A 和 cond 2 ( A ) 。 1 1 1 1 1 1 1 1 解：编写了一个 M文件来求矩阵 A 的范数与条件数： test3_1.m 如下： A=[1 1 1 1;-1 1 -1 1;-1 -1 1 1;1 -1 -1 1]; norm(A,1) norm(A,2) norm(A,inf) cond(A,2) 计算结果依次是： 4 2 4 1.0000 2．实验目的：研究高斯消去法的数值稳定性（出现小主元）。 实验内容：设方程组 Ax b ，其中两个矩阵如下，分别对以上两个方程组 15 0.3 10 59.14 3 1 59.17 5.291 6.130 1 2 46.78 （1） A ， b 1 1 11.2 9 5 2 1 1 2 1 1 2 10 7 0 1 8 3 2.0999999999 9999 6 2 5.9000000000 0001 （2 ） A ， b 2 2 5 1 5 1 5 0 1 0 2 1 （1）计算矩阵的条件数，判断系数矩阵是良态的还是病态的？ 解： 本题编写了一个 test3_21 的 M 文件如下： A1=[0.3*1e-15 59.14 3 1;5.291 -6.130 -1 2;11.2 9 5 2;1 2 1 1]; A2=[10 -7 0 1;-3 2.099999999999999 6 2;5 -1 5 -1;0 1 0 2]; cond(A1) cond(A2) 求得两个矩阵的条件数分别为 68.4296 和 8.9939 ，易知这矩阵 A1 的条件数远远大于 1，而 矩阵 A2 的条件数刚大于 1，故这 , 矩阵 A1 为病态矩阵，矩阵 A2 为良态矩阵。 （2 ）用列主元消去法求得 L 和 U及解向量 x ,x R4 ；