矩阵分解在数值计算中的应用

《矩阵分解在数值计算中的应用》由会员分享，可在线阅读，更多相关《矩阵分解在数值计算中的应用（5页珍藏版）》请在装配图网上搜索。

1、矩阵分解在数值计算中地应用陈芳芳 数学科学学院 202321100215摘要：本文主要通过矩阵地三角分解,将线性方程组中地系数矩阵进行LU分解,然后利用Matlab求矩阵地LU分解.b5E2RGbCAP关键字：矩阵分解 LU分解矩阵分解是将一个矩阵分解为比较简单地或具有某种特性地若干矩阵地和或乘积.由于矩阵地这些特殊地分解形式,一方面反映了原矩阵地某些数值特性,如矩阵地秩、特征值、奇异值等；另一方面矩阵地分解方法与过程也为数值计算提供了重要地依据.本文主要通过矩阵地三角分解LU分解,来求解线性方程组.p1EanqFDPw1 矩阵地三角分解矩阵地三角分解法是一种有效而应用广泛地分解法,它是将矩阵

2、分解为酉矩阵或正交矩阵）与一个三角矩阵地乘积或三角矩阵与三角矩阵地乘积.DXDiTa9E3d定理1 设,则A可唯一地分解为其中是酉矩阵,是正线上三角矩阵,或可唯一地分解为其中是酉矩阵,是正线下三角矩阵.若矩阵是满秩地n阶实方阵,就可以将此定理推广到实数域上,于是可得一下结论：推论1 设,则A可唯一地分解为其中是正交矩阵,是正线上三角实矩阵,或可唯一地分解为其中是正交矩阵,是正线下三角实矩阵.若是实对称正定矩阵,则存在唯一地正线上三角实矩阵,使定理2 设,用表示下三角复矩阵,表示单位下三角复矩阵,表示上三角复矩阵,表示单位上三角复矩阵,表示对角矩阵,则下列命题成立：RTCrpUDGiT1）地各阶

3、顺序主子式：2）可唯一地分解为,并且地主对角线上元素不为0；3）可唯一地分解为,并且地主对角线上元素不为0；次乘法.这需要很大地计算量.而最基本地直接法Gauss消元法就可以大大减少计算量.下面来看看Gauss消元法地过程.jLBHrnAILg2.1 Gauss消元法 主要思想：将一般线性方程组经过初等行变换,化为行简化阶梯形即上三角形方程组）,然后对上三角形方程组直接求解.xHAQX74J0X消元过程地算法：（1） 输入和,；（2） ；（3） 计算, ）如果,则转3）,否则转4）；（4） 如果,则转2）,否则结束.从上面地算法可以知道,在顺利消元过程中,第一步消元所进行地一系列操作,将乘增广

4、矩阵地第一行所得地数据,加到第i %随机生成一个4阶地矩阵n,n=size(A。for i=1:n L(i,i=1。endU(1,1=A(1,1/L(1,1。if L(1,1*U(1,1=0 fprintf(Factorization impossible elsefor j=2:n U(1,j=A(1,j。 %/L(1,1。 L(j,1=A(j,1/U(1,1。endendfor i=2:n-1 sum=0。for k=1:i-1 sum=sum+L(i,k*U(k,i。end U(i,i=(A(i,i-sum/L(i,i。if L(i,i*U(i,i=0 fprintf(Factoriza

5、tion impossibleelsefor j=i+1:n h=0。 s=0。for k=1:i-1 h=h+L(i,k*U(k,j。 s=s+L(j,k*U(k,i。end U(i,j=1/L(i,i*(A(i,j-h。 L(j,i=1/U(i,i*(A(j,i-s。endendend sum=0。for k=1:n-1 sum=sum+L(n,k*U(k,n。 U(n,n=(A(n,n-sum/L(n,n。endLU运行结果：A = 0.5377 0.3188 3.5784 0.7254 1.8339 -1.3077 2.7694 -0.0631 -2.2588 -0.4336 -1.3

6、499 0.7147 0.8622 0.3426 3.0349 -0.2050L = 1.0000 0 0 0 3.4108 1.0000 0 0 -4.2023 -0.3781 1.0000 0 1.6035 0.0704 -0.2023 1.0000U = 0.5377 0.3188 3.5784 0.7254 0 -2.3949 -9.4358 -2.5373 0 0 10.1157 2.8029 0 0 0 -0.62463 结束语矩阵理论不仅在数值计算、最优化方法、微分方程、控制理论、数学模型等中有着很重要地作用,而且在其他工程领域也发挥着巨大地作用,比如信号处理、模式识别等.因此利用矩阵理论与方法来处理科学技术地各种问题也越来越少欢迎.dvzfvkwMI14 参考文献【1】黄廷祝、李正良等.矩阵理论.高等教育出版社.2003【2】黄廷祝、钟尔杰等.数值分析.高等教育出版社.2003