LU和QR分解法解线性方程组.docVIP

LU和QR法解线性方程组 问题描述 求解方程组==, 要求： 编写用三角（LU）分解法求解线性方程组； 编写用正交三角（QR）分解法求解线性方程组。 问题分析 求解线性方程组Ax=b，其实质就是把它的系数矩阵A通过各种变换成一个下三角或上三角矩阵，从而简化方程组的求解。因此，在求解线性方程组的过程中，把系数矩阵A变换成上三角或下三角矩阵显得尤为重要，然而矩阵A的变换通常有两种分解方法：LU分解法和QR分解法。 LU分解法： 将A分解为一个下三角矩阵L和一个上三角矩阵U，即：A=LU, 其中 L=， U= QR分解法： 将A分解为一个正交矩阵Q和一个上三角矩阵R,即：A=QR 三、实验原理 LU分解法 解Ax=b 的问题就等价于要求解两个三角形方程组　⑴ Ly=b,求y;??? Ux=y,求x.设A为非奇异矩阵,且有分解式A=LU, L为单位下三角阵,U为上三角阵直接三角分解法解Ax=b（要求A的所有顺序主子式都不为零）的计算公式, ,i=2,3,…，n.计算U的第r行，L的第r列元素(i=2,3,n)：?②? ? ， i=r,r+1,…,n; ③ , i=r+1,…,n,且r≠n.求解Ly=b，Ux=y的计算公式； ⑤ QR分解法 实验步骤 1、LU分解法 1将矩阵A保存进计算机中，再定义2个空矩阵L，U以便保