资源预览内容
第1页 / 共5页
第2页 / 共5页
第3页 / 共5页
第4页 / 共5页
第5页 / 共5页
亲,该文档总共5页全部预览完了,如果喜欢就下载吧!
资源描述
对矩阵分解及其应用矩阵分解是指将一个矩阵表示为结构简单或具有特殊性质若干矩阵之积或之和,大体分为三角分解、QR 分解、满秩分解和奇异值分解。矩阵的分解是很重要的一部分内容,在线性代数中时常用来解决各种复杂的问题,在各个不同的专业领域也有重要的作用。秩亏网平差是测量数据处理中的一个难点,不仅表现在原理方面,更表现在计算方面,而应用矩阵分解来得到未知数的估计数大大简化了求解过程和难度。1. 矩阵的三角分解如果方阵 A 可表示为一个下三角矩阵 L 和一个上三角矩阵 U 之积,即A=LU,则称 A 可作三角分解。矩阵三角分解是以 Gauss 消去法为根据导出的,因此矩阵可以进行三角分解的条件也与之相同,即矩阵 A 的前 n-1 个顺序主子式都不为 0,即 .所以在对矩阵 A 进行三角分解的着手的第一步应该是判0断是否满足这个前提条件,否则怎么分解都没有意义。矩阵的三角分解不是唯一的,但是在一定的前提下,A=LDU 的分解可以是唯一的,其中 D 是对角矩阵。矩阵还有其他不同的三角分解,比如 Doolittle 分解和 Crout 分解,它们用待定系数法来解求 A 的三角分解,当矩阵阶数较大的时候有其各自的优点,使算法更加简单方便。矩阵的三角分解可以用来解线性方程组 Ax=b。由于 A=LU,所以 Ax=b 可以变换成 LU x=b,即有如下方程组:=先由 依次递推求得 , ,, ,再由方程 依次递推求得 ,= 1 2 = , , . 1 1必须指出的是,当可逆矩阵 A 不满足 时,应该用置换矩阵 P 左乘 A 以便0使 PA 的 n 个顺序主子式全不为零,此时有:=这样,应用矩阵的三角分解,线性方程组的解求就可以简单很多了。2. 矩阵的 QR 分解矩阵的 QR 分解是指,如果实非奇异矩阵 A 可以表示为 A=QR,其中 Q 为正交矩阵,R 为实非奇异上三角矩阵。QR 分解的实际算法各种各样,有 Schmidt 正交方法、Givens 方法和 Householder 方法,而且各有优点和不足。2.1Schmidt 正交方法的 QR 分解Schmidt 正交方法解求 QR 分解原理很简单,容易理解。步骤主要有:1)把A 写成 m 个列向量 a=( , , ) ,并进行 Schmidt 正交化得1 2 =( , , , ) ;2) 单位化,并令 Q=( , , ) ,1 2 1 2 R=diag( , , )K,其中 a= K;3)A=QR. 这种方法来进行 QR 分1 2 解,过程相对较为复杂,尤其是计算量大,尤其是阶数逐渐变大时,就显得更加不方便。2.2Givens 方法的 QR 分解Givens 方法求 QR 分解是利用旋转初等矩阵,即 Givens 矩阵 Tij(c,s)来得到的,T ij(c,s)是正交矩阵,并且 det(Tij(c,s)=1。T ij(c,s)的第 i 行第 i 列和第 j 行第 j 列为 cos ,第 i 行第 j 列和第 j 行第 i 列分别为 sin 和-sin ,其他的都为 0.任何 n 阶实非奇异矩阵 A 可通过左连乘 Tij(c,s)矩阵(乘积为T)化为上三角矩阵 R,另 Q=T-,就有 A=QR。该方法最主要的是在把矩阵化为列向量的基础上找出 c 和 s,然后由此把矩阵的一步步向上三角矩阵靠近。Givens 方法相对 Schmidt 正交方法明显的原理要复杂得多,但是却计算量小得多,矩阵 Tij(c,s)固有的性质很特别可以使其在很多方面的应用更加灵活。2.3Householder 方法的 QR 分解Householder 方法分解矩阵是利用反射矩阵,即 Householder 矩阵,其中 u 是单位列向量,H 是正交矩阵, 。可以证明,=-2 =-1两个 H 矩阵的乘积就是 Givens 矩阵,并且任何实非奇异矩阵 A 可通过连乘Householder 矩阵(乘积为 S)化为上三角矩阵 R,则 A= QR。这种方法首要的就是寻找合适的单位列向量去构成矩阵 H,过程和 Givens 方法基本相似,但是计算量要小一些。矩阵的 QR 分解可以用来解决线性最小二乘法的问题,也可以用来降低矩阵求逆的代价。矩阵的求逆是件不小的工程,尤其是阶数慢慢变大的情况时,而用先把矩阵 QR 分解成正交矩阵和上三角矩阵,就容易多了,况且正交矩阵的转置就是逆,这一点是其他的矩阵分解无法比拟的。在解求线性方程组中,如果系数矩阵的阶数比较大,可以利用 QR 分解来使计算简单化。另外,QR 分解考虑的是 n 阶矩阵,其他的矩阵是不能用这种方法进行分解,由于 QR 分解的这一前提条件,使得下面提到的满秩矩阵分解和奇异值分解就有了其特殊的意义。3. 满秩分解满秩分解也称最大秩分解,前面的 QR 分解是面对 n 阶矩阵的,而满秩分解可以处理长方阵。满秩分解是指,把秩为 r 的 mxn 矩阵 A 分解成 A=FG,其中 F是秩为 r 的 mxr 阶矩阵,G 是秩为 r 的 rxn 阶矩阵。满秩矩阵的解求可以通过初等变换法,但是必须经过多次求逆,所以就利用 Hermite 行标准形来完成。把矩阵 A 经过变换成为 Hermite 行标准形 B,B 的 j1,j2,,j r列为单位矩阵Im的前 r 列,另 A 的第 j1,j2,,j r列为矩阵 F,B 的前 r 行为矩阵 G,则有A=FG。在广义逆中,满秩分解有很多的应用。在证明 A1的存在性时就需要用到 Hermite 行标准形来得到“对于任一的矩阵,总是存在非奇异矩阵 Q 和置换矩阵 P,使 ”,之后才能构造 来证明 A1是存在=( 0 0 0) =( 0 0 )的。用矩阵的满秩分解还能构造 A+,若矩阵 A 有满秩分解,即 A= FG,则可以证明有 。+=()14. 奇异值分解矩阵的奇异值分解是线性代数中一种重要的矩阵分解,在最优化问题、特征值问题、最小二乘问题和广义逆问题及统计学问题中都有重要的应用。对秩为 r 的 mxn 阶矩阵 A 进行奇异值分解的步骤是:1)求得 AHA 的特征值,及对应的特征向量并正交单位化,得矩阵 V,使得1, 2, ;2)将 V 的前 r 列作为 ,令()=2 00 0=(1, 2,) 1,再扩张 U1成 m 阶的矩阵 U;3)那么 。从计算过1=11 =00 0程中可以看出,矩阵的奇异值分解解求是由矩阵的特征值开始的,因此这种分解自然和特征值的问题有莫大联系的。在广义逆问题中,矩阵的奇异值分解的作用一样不可代替。在证明A1,2,3的存在性时,首先就需要用奇异分解来得到一个结论:r(A HA)= r(AAH)= r(AH)= r(A),由此得到的 AH可以由 AHA 表示,再去证明 A1,2,3应该满足的条件就方便得多了。另外,在构造 A+的过程中也有应用,若 A 有奇异值分解 ,则有可以得到 。 +=(M 0 0 0) +=(M1 0 0 0)5. 奇异值分解应用于秩亏网平差在经典平差中,都是以已知的起算数据为基础,将控制网固定在已知数据上,比如水准网必须至少知道已知网中某一点的高程,平面网至少要已知一个点的坐标、一条边的边长和一条边的方位角。此时,误差方程的系数矩阵 B 总是列满秩的,由此得出的法方程系数阵 是个对称的满秩方阵,即=,法方程有唯一解。当网中没有必要的起算数据时(引起秩亏的原()=()因) ,网中所有点均为待定点,就为自由网,B 为列亏矩阵,秩亏数为 d(必要的起算数据个数),误差方程为:=组成的法方程为: =0若是按照直接解法用如下的方程组来解求 x 的解:=-=0=(a)可以得到 ,即该方程组有解但不唯一,虽然满足最小二乘准则,但|=0有 无穷多组解,无法求得 的唯一解,这是与经典平差的根本区别。 为了求得唯一解,必须增加新的约束条件。秩亏自由网平差就是在满足最小二乘 和最小范数 的条件下,求参数一组最佳估值的平= =差方法,也就是通过对如下的方程组来解求 的唯一解:(b)=-=0=这是个复杂的方程组,如果按部就班按照正常求解的方法是很困难的,下面我们把矩阵的奇异值分解融合进来。我们首先根据前面矩阵奇异分解的步骤求得矩阵 B 的奇异值分解:,在此基础上令矩阵 。通过矩阵理论的学习我=00 0 =1 00 0们知道,我们可以通过如下的方式来验证 G 就是 B 的广义逆:(1) = 00 0 1 00 000 0=00 0=(2) =1 00 0 00 01 00 0=1 00 0=(3) ()=(00 01 00 0)=(4) ()=(1 00 000 0)=我们知道,对于不相容方程组 ,使得 为极小范数最小二乘的充= =要条件是 为 的广义逆。而我们已经得到了 就是 的广义逆,那么就说明 G 是 满足该方程式的极小范数最小二乘解。也就是说,我们得到未知参数的估值。通过这种方式,我们求解方程组(b)就简单多了,=1 00 0矩阵的奇异分解令问题很容易的简单化了。6.结论矩阵的分解还有很多的应用,比如可以用来求矩阵的秩,对于阶数偏大的矩阵,即使用初等变换的方法,也是计算量很大的,而把矩阵分解后可以使计算简单。再如,在线性代数中求矩阵的 n 次幂是很常见的,若是一板一眼的进行矩阵相乘,当 n 较大时计算量可想而知,况且,当 n 逐渐增大或是非纯数据间的运算的情况下,根本就没有计算的可能,此时,矩阵分解方法的应用可以令问题变得简单而易懂。判断矩阵的正定性需要不断的计算行列式,计算量大而复杂,矩阵分解可以使之更简单直接。矩阵的分解作用很广泛,在不同的领域都发挥着其独特的作用,只要应用得好,肯定可以使原有的问题简单而易于理解。我们知道,矩阵理论就其理论来说,对于除了数学本专业的人而言,意义是不大的。纯理论的学习是枯燥而乏味的,只有和是具体问题的结合才会显出它的强大生命力。单看一个定理还是推论,我们会觉得它是简单而几乎没有意义的,甚至不知道怎么去理解它以及存在的意义,当运用到实际的领域,一方面我们可以更好的了解相关的知识,重要的是解决了具体的问题。这应该就是学习的乐趣所在。在测量平差的秩亏网平差中,解求未知数的估计值时候和奇异值分解结合起来,不仅可以使得运算更加简单化,并且得到的结果更利于理解
网站客服QQ:2055934822
金锄头文库版权所有
经营许可证:蜀ICP备13022795号 | 川公网安备 51140202000112号