理学院数值计算方法实验.docVIP

实验二 实验名称 线性方程组的数值解法 姓 名 张见 学 号班 级 信计（2） 指导教师 张昆 实验日期 2010-11-12 成 绩 实验目的 掌握Gauss消去法及Gauss列主元消去法，能用这两种方法求解方程组； 掌握Jacobi和G-S迭代法，能应用Jacobi和G-S迭代法求解方程组； 掌握相应数值算法的程序编写； 理解迭代法收敛的充要条件，会判断迭代法的收敛性。 实验题目 分别用列主元消去法与不选主元消去法求解，分析算法对结果的影响： a) b) 给定矩阵A=hilb(n)与向量b , 求的LU分解； 利用的LU分解求解方程组：；；x=b； 利用的LU分解求，n值自己给定。 实验原理 1、 当系数矩阵A的秩增广矩阵B的秩时，方程组无解向量； 当系数矩阵A的秩=增广矩阵B的秩n时，方程组有无数解向量； 当系数矩阵A的秩=增广矩阵B的秩=n时，方程组有唯一解向量； ①列主元消元法 当方程组有唯一解向量时， 增广矩阵B的第i()次消元后为，在 的第i+1列中与所在的行互换，然后把下的元素化为零。 最后增广矩阵化为上三角矩阵，再通过回代求得方程组的解向量 ②不选主元消去法 当方程组有唯一解向量时， 用高斯消元法求解方程组，即把方程组的增光矩阵化为上三角矩阵，再通过回代求得方程组的解向量 2.设A为n阶矩阵，如果A的顺序主子式，则A可分解为一个单位下三角矩阵L和一个上三角矩阵U的乘积，且这种分解是唯一的。即：A=LU 其中 首先，从第1行得出U的第1行的元素 再从第1列算出L的第1 列元素 其次，从第2行算出U的第2行元素 再从第2 列算出L的第2 列元素 一般的，设已经给出U的第1 行到第k-1行元素与L的第1列到第k-1列元素，则U的第k行元素为 L的第k列元素为 综上所述，可得用直接三角分解求L ,U 计算公式为： 求解Ax=bLUx=b,求得方程组的解向量 实验内容 1. ①列主元消元法 当方程组有唯一解向量时， 增广矩阵B的第i()次消元后为，在 的第i+1列中与所在的行互换，然后把下的元素化为零。 最后增广矩阵化为上三角矩阵，再通过回代求得方程组的解向量 ②不选主元消去法 当方程组有唯一解向量时， 用高斯消元法求解方程组，即把方程组的增光矩阵化为上三角矩阵，再通过回代求得方程组的解向量 2.a)求A的LU分解 设A为n阶矩阵，如果A的顺序主子式，则A可分解为一个单位下三角矩阵L和一个上三角矩阵U的乘积，且这种分解是唯一的。即：A=LU 其中 首先，从第1行得出U的第1行的元素 再从第1列算出L的第1 列元素 其次，从第2行算出U的第2行元素 再从第2 列算出L的第2 列元素 一般的，设已经给出U的第1 行到第k-1行元素与L的第1列到第k-1列元素，则U的第k行元素为 L的第k列元素为 综上所述，可得用直接三角分解求L ,U 计算公式为： 求解Ax=bLUx=b,求得方程组的解向量 b) 利用A的LU分解求解下列方程组： (1) Ax=b Ax=bLUx=b,求得方程组的解向量 (2) (3) c)利用A的LU分解求 实验结果 1) a) (1) 列主元消元法 e=7.1054273576 (2) 不选主元消元法 e=2.836931460819089 (1) 列主元消元法 e= 1.01753620972 (2) 不选主元消元法 e= 7.961155371425798 2) a) 求A的LU分解 , b) 利用A的LU分解求解下列方程组： ① Ax=b e= 6.301164130741664 KA =1.61889985892434 ② ； e= 2.8107995