矩阵分析课件chapt4 矩阵分解例题详解.docVIP

第4章 矩阵分解 (I)高斯消去法 假设矩阵A 的顺序主子式?0 (i=1,…,n), 则我们可以进行以下的顺序消元过程 1．消元过程 等价于用初等矩阵分别 左乘和，即 (1) 其中，, 我们称为消元因子，为主元素； 消元过程的一个重要性质是：消元过程不改 变矩阵的主子矩阵的行列式(主子式)的值。 例 ,顺序主子式为，1，5，?10 ,顺序主子式为，1，5，?10 ，顺序主子式为，1，5，?10 引理：约化的主元素?0的充要条件是 矩阵A 的顺序主子式?0 (i=1,…,k); 推论：若矩阵A 的顺序主子式?0 (i=1,…,k)，则 ，； 由此有若A对称正定或严格对角占优，而 它们的顺序主子矩阵也是对称正定或严格 对角占优，从而顺序主子式不为0，顺序高斯 消去过程可进行； 2.回代过程： 设， , 用高斯消去法解线性方程Ax=b. 增广矩阵为 ， 因此，问题的解为 3. 数值稳定性 1)选列主元 ； 2)选全主元； 3)高斯若当(Gauss-Jordan)消去法，求矩阵的逆; , 求A?1. 增广矩阵为 从而 4．高斯顺序消元法解方程的计算量 1)乘除次数： 2)加减次数： 3)求矩阵的逆的计算量为o() (II) 顺序消元过程与矩阵的三角分解 (1) (2) 若 则有 ，从而 ， (3) 由 有 故有 A=LU，其中. 这时L为单位下三角矩阵。 矩阵的三角分解 A=(a1,a2,…,an)=PR=(p1,p2,…,pn)R (a) LU分解(Doolittle分解) (1) 存在唯一的条件;(顺序主子式不为0) (2) 公式推导;(矩阵乘法) 定义4.1 如果方阵A可分解成一个下三角 矩阵L和一个上三角矩阵U的乘积，则称 A可作三角分解. 如果方阵A可分解成 A=LDU,其中L为一个单位下三角矩阵， D为对角矩阵，U是一个单位上三角矩阵。 推论：若矩阵A 的顺序主子式?k?0 (k=1,…,n-1)，则A可唯一分解为 A=LDU, 其中L为一个单位下三角矩阵， D为对角矩阵，U是一个单位上三角矩阵。 dk=?k/?k?1 (?0=1)。 定理4.2 设A是n阶非奇异矩阵，则存在 置换矩阵P使得PA=LDU, 其中L为一个 单位下三角矩阵，D为对角矩阵，U是一个 单位上三角矩阵。 定义4.2 设A存在唯一的LDU分解.若把 A=LDU中的D和U结合起来，并且用 U’表示，则得到唯一的LU分解 A=LU’ 称为Doolittle分解；若把A=LDU中的L和D 结合成L’，就得到A=L’U 称为Crout分解。 LU分解的公式： lik=aik-(li1u1k+…+li,k-1uk-1,k) ukj= [akj-(lk1u1j+…+lk,k-1uk-1,j)]/lkk Crout分解类似。 (b)平方根法 (1)对称矩阵的三角分解定理； (2)对称正定矩阵的三角分解(Cholesky分解) 递推公式 gii=(aii-)1/2 gij=[aij-(gi1gj1+gi2gj2+…+gi,j-1gj,j-1)]/gii, ij gij=0 ij 四、分块矩阵的拟LU分解和拟LDU分解 若A11可逆，作 则 从而det(A)=det(A11)?det() 同样若A22可逆，可得类似结果。 矩阵求逆引理 (A+BC)?1=A?1? A?1B(I+CA?1B)?1CA?1 证明：求方程(A+BC)x = b, 令y=Cx, 则有 Ax+By=b ?Cx+y=0 写成矩阵为 ， 利用高斯消去法有 因此可得 x = (A?1? A?1B(I+CA?1B)?1CA?1)b, 而由(A+BC)x = b可得x = (A+BC)?1b 由于b的任意性可得 (A+BC)?1=A?1? A?1B(I+CA?1B)?1CA?1 从而命题得证。 推论：(A+BD?1C)?1=A?1? A?1B(D+CA?1B)?1CA?1 例= = 由于 利用矩阵求逆引理有 矩阵的QR分解 Givens变换与Householder变换 Givens变换 定义：设实数c与s满足 c2+s2=1, 称 为Givens 矩阵，也可记作Tij=Tij(c,s). 由Givens矩阵确定 线性变换称为Givens 变换。 性质1. Givens矩阵是正交矩阵，且有 [Tij(c,s)]?1=[Tij(c,s)]T=Tij(c,?s) det[Tij(c,s)]=1 性质2 设x=(?1,?2,?,?n)T,y=Tijx=(?1,?2,?,?n)T,则有