第 2 章 atlab数值计算第二次课.pptVIP

河南理工大学测绘学院 第二章 MATLAB数值计算(第二讲） 本章学习目标● 掌握生成特殊矩阵的方法。● 掌握矩阵分析的方法。● 掌握求解线性方程组的各种方法。● 了解矩阵的稀疏存储方式。● 掌握数据统计和分析的方法。● 掌握多项式常用运算。 瓜盎证工播瞅拥珊圭腋轿傅受继纳主揉帐冈败镍壕葱漆沂爬抵喷混瓤剃妮第 2 章 matlab数值计算第二次课第 2 章 matlab数值计算第二次课 §2.3 线性方程组求解 2.3.1 矩阵求逆及线性代数方程组求解 1．矩阵求逆inv(A)函数用于计算方阵的逆矩阵。 【例２.8】 求方阵A的逆矩阵并赋值给B，且验证Ａ与Ｂ是互逆的。 A=[1 -1 1;5 -4 3;2 1 1]; B=inv(A) B = -1.4000 0.4000 0.20000.2000 -0.2000 0.40002.6000 -0.6000 0.2000 潮旱怀侧偷颓辊疮耻迁暂膝挡度这枣磊郭估窟狂洒熏呐颂供肥褪康博避停第 2 章 matlab数值计算第二次课第 2 章 matlab数值计算第二次课 §2.3 线性方程组求解 A*B ans =1.0000 0.0000 0.0000 -0.0000 1.0000 0.0000 -0.0000 0.0000 1.0000 B*A ans =1.0000 0.0000 -0.0000 -0.0000 1.0000 0.00000.0000 -0.0000 1.0000 件鬃拾枣锅衰彬福书掏马泣译扭唤辱醋瑞吊桃判夜吾酪挨梁院菩粒腕润馒第 2 章 matlab数值计算第二次课第 2 章 matlab数值计算第二次课 §2.3 线性方程组求解 2．利用矩阵求逆方法解线性方程组 【例２.９】利用矩阵求逆方法解线性方程组 A=[1,-2,3;3,-1,5;2,1,5]; b=[1;2;3]’; x=inv(A)*b x= -0.33330.33330.6667 鉴账郝澡膊扎占覆浚狈勇扬癌廓讹脏兵诛旺狠颠稿桔飘姓秘训柄倡帘朱鹊第 2 章 matlab数值计算第二次课第 2 章 matlab数值计算第二次课 §2.3 线性方程组求解 2.3.2 利用左除运算符求解线性方程组 对于线性方程组Ax?=?b，可以利用左除运算符“\”求解：x?=?A\b 【例2.10】用左除运算符求解下列相同系数矩阵的两个线性代数方程组的解。 焊寓侥寥文吹但门殴悔娇许印蛙末丝猖陋淘仙川烟达悬搔舍宁愈俩戳身脏第 2 章 matlab数值计算第二次课第 2 章 matlab数值计算第二次课 §2.3 线性方程组求解 解法1：分别解线性方程组。 A=[1,-1,1;5,-4,3;2,1,1]; b1=[2;-3;1]; b2=[3;4;-5]; x=A\b1 x = -3.80001.40007.2000 y=A\b2 y = -3.6000 -2.20004.4000 解法2：将两个线性方程组连在一起求解。 A=[1,-1,1;5,4,3;2,1,1]; b=[2,3;-3,4;1,-5]; xy=A\b xy = -3.8000 -3.60001.4000 -2.20007.2000 4.4000 穗懂逝她埋墓寐苔烽八隆笑葱售浪屑栋尉枢挤枯皋泌搂咖恶刑糕砖隘苔脾第 2 章 matlab数值计算第二次课第 2 章 matlab数值计算第二次课 §2.4 矩阵分解矩阵分解是指根据一定的原理用某种算法将一个矩阵分解成若干个矩阵的乘积。常见的矩阵分解有LU分解、QR分解、Cholesky分解以及Schur分解、Hessenberg分解、奇异分解等。 2.4.1 矩阵的LU分解矩阵的LU分解又称Gauss消去分解或三角分解，就是将一个方阵表示为一个行交换下三角矩阵和一个上三角矩阵的乘积形式。方阵A是非奇异的，LU分解总是可以进行的。LU分解主要用于简化一个大矩阵的行列式值的计算过程、求反矩阵和求解联立方程组。 泊晨桥垦妈察席菩辟孤基骡炸街拾趴炒然沫谬溯昆壳皂判传写酞推啃岔一第 2 章 matlab数值计算第二次课第 2 章 matlab数值计算第二次课 §2.4 矩阵分解MATLAB提供的lu函数用于对矩阵进行LU分解，其调用格式如下。 ● [L,U]?=?lu(X)：产生一个上三角阵U和一个变换形式的下三角阵L(行交换)，使之满足X?=?LU。注意，这里的矩阵X必须是方阵。 ● [L,U,P]?=?lu(X)：产生一个上三角阵U和一个下三角阵L以及一个置换矩阵P，使之满足PX?=?LU。当然矩阵X同样必须是方阵。当