数值分析实验报告之矩阵的分解.docVIP

矩阵分解实验报告 学院：数学与计算机学院 专业： 班级： 姓名 学号 实验组 实验时间 指导教师 成绩 实验项目名称 矩阵分解 实验目的 利用高斯消去法将一个非奇异矩阵分解为一个上三角矩阵和一个下三角矩阵。 实验要求 掌握矩阵分解的数学原理和算法。 实验原理 对于线性方程组中的矩阵，将分解为一个上三角矩阵和一个下三角矩阵的乘积，即，为单位下三角矩阵，为上三角矩阵。 ，即分解为和两个上（下）三角矩阵。 实验仪器 1、计算机一台。 2、安装IIS。 实验步骤 Step1：选择数据文件； Step2：将文件数据读取到内存； Step3: 初始化上三角矩阵和下三角矩阵，即， ,; Step4: ； Step5: 计算的第行元素; Step6: 判断是否为0，若为0，则显示原矩阵为非奇异矩阵，不可分解，过程结束。 Step7: 的第列元素； Step8: 判断是否大于，若是，转到Step9，否则，令，返回Step5。 Step9: 将得到的矩阵结果保存至文件，过程结束。 实 验 总 结 本次试验也是针对求解线性方程组解的问题。将一个非奇异矩阵分解之后，很快就可以将方程组的解求出来了。在实验的过程中，我首先复习了一遍矩阵分解的数学原理，然后开始在自己的理解之下进行算法的实现。因为理解了实验所需的数学原理，因此很快将算法实现了。我首先测试了课本上的一个三阶的矩阵，发现求解的结果是错误的，仔细查找才发现是初始化的时候由于手误将循环变量写错了，改过来之后，结果便对了，接着又从网上找到了两组测试数据，结果均正确。然后才测试老师所给的800阶的矩阵，发现算法还需要一段运行时间才能将结果呈现出来。本次试验使我对矩阵分解更加熟悉了，对用矩阵分解来解线性方程组有了更进一步的认识。在实验课之后，我对矩阵分解之后再利用分解的结果求解线性方程组的解进行了实现。 实 验 流 程 图