计算机科学与技术系数值计算方法第3章2矩阵的三角分解法矩阵求逆-ppt课件

,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,3.2 矩阵的三角分解法,我们知道对矩阵进行一次初等变换，就相当于用相应的初等矩阵去左乘原来的矩阵。因此我们这个观点来考察Gauss消元法用矩阵乘法来表示，即可得到求解线性方程组的另一种直接法：矩阵的三角分解。,3.2 矩阵的三角分解法我们知道对矩阵进行一次初等变换，就,1,3.2.1 Gauss消元法的矩阵形式,3.2.1 Gauss消元法的矩阵形式,2,计算机科学与技术系数值计算方法第3章2矩阵的三角分解法矩阵求逆-ppt课件,3,计算机科学与技术系数值计算方法第3章2矩阵的三角分解法矩阵求逆-ppt课件,4,计算机科学与技术系数值计算方法第3章2矩阵的三角分解法矩阵求逆-ppt课件,5,计算机科学与技术系数值计算方法第3章2矩阵的三角分解法矩阵求逆-ppt课件,6,3.2.2 Doolittle分解,3.2.2 Doolittle分解,7,计算机科学与技术系数值计算方法第3章2矩阵的三角分解法矩阵求逆-ppt课件,8,Doolittle分解,若矩阵A有分解：A=LU，其中L为单位下三角阵，U为上三角阵，则称该分解为Doolittle分解，可以证明，当A的各阶顺序主子式均不为零时，Doolittle分解可以实现并且唯一。,Doolittle分解若矩阵A有分解：A=LU，其中L为单位,9,A的各阶顺序主子式均不为零，即,A的各阶顺序主子式均不为零，即,10,Doolittle分解,Doolittle分解,11,Doolittle分解,Doolittle分解,12,Doolittle分解,Doolittle分解,13,Doolittle分解,Doolittle分解,14,Doolittle分解,Doolittle分解,15,Doolittle分解,Doolittle分解,16,例题,例题,17,例题,例题,18,例题,例题,19,例题,例题,20,例题,例题,21,Doolittle分解,Doolittle分解,22,3.2.3 对称矩阵的Cholesky分解,在应用数学中，线性方程组大多数的系数矩阵为对称正定这一性质，因此利用对称正定矩阵的三角分解式求解对称正定方程组的一种有效方法，且分解过程无需选主元，有良好的数值稳定性。,3.2.3 对称矩阵的Cholesky分解在应用数学中，线性,23,对称矩阵的Cholesky分解,A对称：A,T,=A,A正定：A的各阶顺序主子式均大于零。即,对称矩阵的Cholesky分解 A对称：AT=A,24,对称矩阵的Cholesky分解,由Doolittle分解，A有唯一分解,对称矩阵的Cholesky分解由Doolittle分解，A有,25,对称矩阵的Cholesky分解,定理3.2.4 设,A,为对称正定矩阵，则存在唯一分解,A,=,LDL,T,，其中,L,为单位下三角阵，,D,=diag(,d,1,d,2,d,n,)且,d,i,0(,i,=1,n,),对称矩阵的Cholesky分解定理3.2.4 设A为对称正定,26,对称矩阵的Cholesky分解,证明：,对称矩阵的Cholesky分解证明：,27,对称矩阵的Cholesky分解,对称矩阵的Cholesky分解,28,对称矩阵的Cholesky分解,对称矩阵的Cholesky分解,29,对称矩阵的Cholesky分解,推论：设A为对称正定矩阵，则存在唯一分解,其中L为具有主对角元素为正数的下三角矩阵。,对称矩阵的Cholesky分解推论：设A为对称正定矩阵，则存,30,对称矩阵的Cholesky分解,证明：,对称矩阵的Cholesky分解证明：,31,Cholesky分解的求法,Cholesky分解的求法,32,Cholesky分解的求法,Cholesky分解的求法,33,Cholesky分解的求法,Cholesky分解的求法,34,Cholesky分解法,Cholesky分解法缺点及优点,优点：可以减少存储单元。,缺点：存在开方运算，可能会出现根号下负数。,Cholesky分解法Cholesky分解法缺点及优点,35,改进Cholesky分解法,改进的cholesky分解,A,=,LDL,T,改进Cholesky分解法改进的cholesky分解A=LD,36,改进的cholesky分解,改进的cholesky分解,37,改进的cholesky分解,改进的cholesky分解,38,计算机科学与技术系数值计算方法第3章2矩阵的三角分解法矩阵求逆-ppt课件,39,改进的cholesky分解算法,改进的cholesky分解算法,40,改进的cholesky分解算法,改进的cholesky分解算法,41,例题,例题,42,例题,例题,43,例题,例题,44,例题,例题,45,A=LDL,T,分解，既适合于解对称正定方程组，也适合求解A为对称，而各阶顺序主子式不为零的方程组,而对A=LL,T,只适合于对称正定方程组,A=LDLT分解，既适合于解对称正定方程组，也适合求解A为对,46,3.2.4 三对角方程组求解的追赶法,3.2.4 三对角方程组求解的追赶法,47,三对角方程组求解的追赶法,三对角方程组求解的追赶法,48,三对角方程组求解的追赶法,三对角方程组求解的追赶法,49,三对角方程组求解的追赶法,三对角方程组求解的追赶法,50,三对角方程组求解的追赶法,其计算工作量为5n-4次乘除法。工作量小，其实现的条件为,q,i,不为零。有以下定理可得证三对角矩阵求解的充分性条件。,三对角方程组求解的追赶法其计算工作量为5n-4次乘除法。工作,51,解三对角矩阵线性方程组的追赶法程序框图,解三对角矩阵线性方程组的追赶法程序框图,52,3.3 矩阵求逆,3.3 矩阵求逆,53,矩阵求逆,矩阵求逆,54,矩阵求逆,矩阵求逆,55,为使求逆过程不断提高求解精度，因此增加选主元工作，最常用的是选列主元求逆。因此增加一个数组,Z,(,n,)，记录选主元的交换号，最后在消元工作完成后，根据,Z,(,n,)对,A,中的元素进行相应的列交换，得到,A,-1,GaussJordan原地求逆法,为使求逆过程不断提高求解精度，因此增加选主元工作，最常用的是,56,算法（原地求逆法）,算法（原地求逆法）,57,计算机科学与技术系数值计算方法第3章2矩阵的三角分解法矩阵求逆-ppt课件,58,例题,例题,59,例题,例题,60,例题,例题,61,例题,例题,62,例题,例题,63,例题,例题,64,