《矩阵分析应用》编程作业.docx

矩阵分析与应用——Programming（编程） 1 题目测试矩阵 1.1编程要求：1）不限编程语言，程序为可执行文件，例如 .m文件等，不能是word或者txt文档；2）程序可以实现任意矩阵的LU分解，并附上简单实例；3）同时提交程序的说明文档，对程序简单说明； 1.2 测试矩阵 程序编写后，将从讲义的例子和习题中，选三个矩阵A、B、C进行测试，其中A、B为3*3矩阵，可进行LU分解，C为4*4矩阵，主对角元素有0，需进行PLU分解： 矩阵A： L： U： 矩阵B： L： U： 矩阵C： P L U 事先在MATLAB工作区中输入需要测试的矩阵A、B、C，如下图： 2 LU分解 2.1编程思路： 输入需要计算的矩阵： 第一步： 保存输入矩阵A的行数，用于循环。因为LU分解的前提条件是主子矩阵非奇异，非奇异可逆矩阵为方阵??满秩，所以矩阵A的行数=列数，即r=c。 第二步： 判断消元过程中主元是否为0，若为0，则结束程序，若不为0，进行第三步。 第三步： 高斯消元法，消元结果即上三角矩阵U，消元过程中的乘数作为下三角矩阵L对应的元素。 第四步： 输出计算结果，L、U矩阵。 2.2 程序编写 2.3 程序运行结果 矩阵A测试： 矩阵B测试： 矩阵C测试： A、B矩阵运行结果符合预期值，矩阵C不能进行LU分解。故程序符合要求。  3 PLU分解 3.1编程思路： 和LU分解相比，PLU分解需要记录交换矩阵P。 第一步： 循环计算，决定主元行。选取对角线上的主元和下方的元素进行绝对值比较，记录最大值所在行的位置； 第二步： 把最大值所在行交换到主元行； 第三步： 循环计算，高斯消元，步骤同LU分解。 第四步： 输出计算结果。 3.2 程序编写 3.3 程序运行结果 矩阵A： 矩阵B： 矩阵C： 虽然矩阵A、B的对角线元素不为0，但依然可以做PLU分解；矩阵C的计算结果与预期结果一致。故程序编写符合要求。