线性代数方程组的直接法.pptVIP

矩阵运算描述列主元素消去法！*先换行再消元！【定理4】（列主元素的三角分解）如果Ａ为非奇异矩阵，则存在排列矩阵P，使得PA=LU，其中L是单位下三角阵，U是上三角阵。第30页，共80页，星期日，2025年，2月5日【例5】用列主元素消去法解方程组我们用4位浮点数进行计算。解：消元：舍去或者说被“吃”已知精确解：*第31页，共80页，星期日，2025年，2月5日回代计算解：高斯选主元消去法的步骤：注：该解若取两位有效数字，则与精确解完全相同。优点：数值稳定修正方法：消元，回代。高斯-约当（Gauss-Jordam）消去法。缺点：既消元，又回代**第32页，共80页，星期日，2025年，2月5日3、高斯—约当（Gauss–Jordan)消去法高斯消去法始终是消去对角线下方的元素，Gauss–Jordan消去法要消去对角线下方和上方的元素。假设G--J消去法已完成第1步～第k-1步，得到与原方程组等价的方程组，其中第k步计算：（1）按列选主元：即确定ik，使（2）换行：交换的第k行和第ik行元素（3）消元计算：*第33页，共80页，星期日，2025年，2月5日（4）计算主行（主元素所在行）计算解：上述过程结束后，有**第34页，共80页，星期日，2025年，2月5日说明：在解方程组，一般不用高斯-约当消去法。因为计算量太大，但是在解多个方程组而它们的系数矩阵相同时，用此方法，即求系数矩阵的逆矩阵A-1时，比较合适。设为非奇异矩阵。如果用列主元G-J消去法将（A,I）化为（I,T），则。不用回代，将A化为单位矩阵，则解为常数项列。【定理5】（高斯—约当法求逆矩阵）优点：缺点：计算量较大，大约是次乘除法。注：该方法与高等代数中求逆矩阵方法的不同之处是有选主元，实际上选主元就是交换两行的位置，仍是初等变换，在一般的求逆矩阵方法中也有交换两行元素。**第35页，共80页，星期日，2025年，2月5日【例6】用列主元素G-J消去法求的逆矩阵A-1解：??**第36页，共80页，星期日，2025年，2月5日所以G-J列主元求逆算法见P176算法3！**第37页，共80页，星期日，2025年，2月5日**作业：P19912第38页，共80页，星期日，2025年，2月5日**§3Gauss消去法的变形直接三角分解法直接从矩阵A的元素得到计算Ｌ，U元素的递推公式，而不需要任何中间步骤，称之为直接三角分解法。这样，AX=b?L(UX)=b?LY=b，UX=Y。不选主元的三角分解(Doolittle分解法)第39页，共80页，星期日，2025年，2月5日通过比较法直接导出L和U的计算公式。思路*通过n步可以直接计算定出L，U的元素，其中第r步确定U的第r行和L的第r列元素。第1步：设第r-1步已确定U的第r-1行和L的第r-1列元素。第40页，共80页，星期日，2025年，2月5日*第r步：故：故：第41页，共80页，星期日，2025年，2月5日LU分解求解线性方程组:**第42页，共80页，星期日，2025年，2月5日直接三角分解法解AX=b的计算公式:对于r=2,3,…,n，(2)计算U的第r行元素(3)计算L的第r列元素(r?n)(1)**第43页，共80页，星期日，2025年，2月5日(4)(5)**第44页，共80页，星期日，2025年，2月5日【例7】用直接三角分解法解解：用分解计算公式得：求解：**第45页，共80页，星期日，2025年，2月5日**选主元的三角分解选主元素三角分解算法见P179!第46页，共80页，星期日，2025年，2月5日2.平方根法(对称正定矩阵而言)矩阵的LDR分解【定理6】若n阶矩阵A的所有顺序主子式均不等于零，则矩阵A存在唯一的分解式A=LDR，其中L和R分别是n阶单位下三角阵和单位上三角阵，D是n阶对角元素不为零的对角阵，上述分解也称为A的LDR分解。**第47页，共80页，星期日，2025年，2月5日平方根法如果A为对称矩阵，且A的所有顺序主子式均不等于零，