线性代数方程组的直接解法详解演示文稿.pptVIP

例3.5 用矩阵的初等行变换求解解方程组 解: 用矩阵的初等行变换求解,对增广矩阵 (下面带下划线元素为主元素) 当前第31页\共有72页\编于星期二\12点 当前第32页\共有72页\编于星期二\12点 3.3 矩阵三角分解法 3.3.1 矩阵三角分解原理 应用高斯消去法解n阶线性方程组Ax=b, 经过n步消元之后, 得出一个等价的上三角型方程组 A(n) x=b(n), 对上三角形方程组用逐步回代就可以求出解来。上述过程可通过矩阵分解来实现。 将非奇异阵A分解成一个下三角阵L和一个上三角阵U的乘积 A=LU 称为对矩阵A的三角分解，又称LU分解 当前第33页\共有72页\编于星期二\12点 其中 当前第34页\共有72页\编于星期二\12点 方程组Ax=b的系数矩阵A经过顺序消元逐步化为上三角型A(n),相当于用一系列初等变换左乘A的结果。事实上，第1列消元将A(1)=A化为A(2)，若令： 则根据距阵左乘有L1A(1)=A(2) 当前第35页\共有72页\编于星期二\12点 第2列消元将A(2)化为A(3)，若令： 经计算可知 L2A(2)=A(3),依此类推,一般有LkA(k)=A(k+1) 当前第36页\共有72页\编于星期二\12点 mi1= a(1) i1/ a(1) 11 i=2,3,……n 于是矩阵 经过消元化为上三角阵 的过程可表示为 上述矩阵 是一类初等矩阵, 它们都是单位下三角阵，且其逆矩阵也是单位下三角阵,只需将 改为 ,就得到 。即 当前第37页\共有72页\编于星期二\12点 于是有 其中 当前第38页\共有72页\编于星期二\12点 L为由乘数构成的单位下三角阵，U为上三角阵，由此可见，在 的条件下，高斯消去法实质上是将方程组的系数矩阵A分解为两个三角矩阵的乘积A=LU。这种把非奇异矩阵A分解成一个下三角矩阵L和一个上三角矩阵U的乘积称为矩阵的三角分解，又称LU分解。 显然，如果 ,由行列式的性质知，方程组系数矩阵A的前n-1个顺序主子矩阵 非奇异，即顺序主子式不等于零，即 当前第39页\共有72页\编于星期二\12点 其中 （A的主子阵） 反之,可用归纳法证明,如果A的顺序主子式 则 于是得到下述定理： 当前第40页\共有72页\编于星期二\12点 定理3.5 设 。如果A顺序各阶主子式, ,则A可惟一地分解成 一个单位下三角阵L和一个非奇异的上三角阵U的乘积。 证：由于A各阶主子式不为零,则消元过程能进行到底, 前面已证明将方程组的系数矩阵A用初等变换的方法分解成两个三角矩阵的乘积A=LU的过程。 现仅证明分解的惟一性。 设A有两种LU分解 其中 为单位下三角阵， 为上三角阵 ∵ A的行列式 均为非奇异矩阵,有 上式两边左边同乘 ，右边同乘 得 上式左边为单位下三角阵,右边为上三角阵,故应为单位阵,即 惟一性得证。 当前第41页\共有72页\编于星期二\12点 把A分解成一个单位上三角阵L和一个下三角阵U的乘积称为杜利特尔（Doolittle）分解。其中 当前第42页\共有72页\编于星期二\12点 若把A分解成一个下三角阵L和一个单位上三角阵U的乘积称为克洛特分解Crout） 其中 当前第43页\共有72页\编于星期二\12点 3.3.2 用三角分解法解方程组 求解线性方程组Ax=b时,先对非奇异矩阵A进行LU分解使A=LU，那么方程组就化为 LU x=b 从而使问题转化为求解两个简单的的三角方程组 L y=b 求解 y U x=y 求解 x 这就是求解线性方程组的三角分解法的基本思想。下面只介绍杜利特尔（Doolittle）分解法。设A=LU为 当前第44页\共有72页\编于星期二\12点 由矩阵乘法规则 由此可得U的第1行元素和L的第1列元素 当前第45页\共有72页\编于星期二\12点 再确定U的第k行元素与L的第k列元素,对于k=2,3, …,n计