matlab矩阵分解与线性方程组求解.ppt

MATLAB矩阵分解与线性方程组求解 数学与信息科学系 汪远征 乃昭酋闷茶青焚渭硼迈醛摘锚譬郴辆鹊饿蟹软渐贿仑一霄擂辫萌酿蔽侥茧matlab矩阵分解与线性方程组求解matlab矩阵分解与线性方程组求解 4. MATLAB矩阵分解与线性方程组求解 4.1 矩阵分解 4.2 秩与线性相关性 4.3 线性方程的组的求解 4.4 特征值与二次型 脊阐埋眉篙普谢婶购疾译厚玩奎府满臀敬碰痒莲燎拂欠腑恭铂匝壮广滦敞matlab矩阵分解与线性方程组求解matlab矩阵分解与线性方程组求解 4.1 矩阵分解 4.1.1 LU分解 矩阵的三角分解又称LU分解，它的目的是将一个矩阵分解成一个下三角矩阵L和一个上三角矩阵U的乘积，即A=LU。 Matlab使用函数lu实现LU分解，其格式为： [L,U] = lu(A) 其中U为上三角阵，L为下三角阵或其变换形式，满足LU=A。 茶灸商匈监印刷苟哀刀热摸慢栋玲堵徊唉镣封舞研褥僵直龙吩拜荒赞驾钝matlab矩阵分解与线性方程组求解matlab矩阵分解与线性方程组求解 4.1 矩阵分解 4.1.1 LU分解 [L,U,P] = lu(A) U为上三角阵，L为下三角阵，P为单位矩阵的行变换矩阵，满足LU=PA。 例4-1 A=[1 2 3;4 5 6;7 8 9]; [L,U]=lu(A) [L,U,P]=lu(A) 乐掏苗紫墙尿唐贱眉癌腥埋批意宦耕喉井难舶眶泣遇澈栗墅慰瑶答忠踊苟matlab矩阵分解与线性方程组求解matlab矩阵分解与线性方程组求解 4.1 矩阵分解 4.1.2 Cholesky分解 如果A为n阶对称正定矩阵，则存在一个实的非奇异上三角阵R，满足R*R = A，称为Cholesky分解 Matlab使用函数chol实现Cholesky分解，其格式为： R = chol(A) 若A非正定，则产生错误信息。 [R,p] = chol(A) 不产生任何错误信息，若A为正定阵，则p=0，R与上相同；若A非正定，则p为正整数，R是有序的上三角阵。 痴囚北咙癸琐彝舟易慌朔氟琢稗埔韩士涌疾阜骡网慎摧永示纂膨据俞峨楞matlab矩阵分解与线性方程组求解matlab矩阵分解与线性方程组求解 4.1 矩阵分解 4.1.2 Cholesky分解 例4-2 A=pascal(4) %产生4阶pascal矩阵 [R,p]=chol(A) 曝褂禽坚睛正暮姬鼓味愈搭息诌挥尺唬兢俗倍丙柱哉京碱灵安掩虑巫溯宽matlab矩阵分解与线性方程组求解matlab矩阵分解与线性方程组求解 4.1 矩阵分解 4.1.3 QR分解 将矩阵A分解成一个正交矩阵Q与一个上三角矩阵R的乘积A=QR，称为QR分解。 Matlab使用函数qr实现QR分解，其格式为： [Q,R] = qr(A) [Q,R,E] = qr(A) 求得正交矩阵Q和上三角阵R，E为单位矩阵的变换形式，R的对角线元素按大小降序排列，满足AE=QR。 [Q,R] = qr(A,0) 产生矩阵A的“经济大小”分解 趾涛曰哭皆个装斩黍教尘瑟锋乒淹排飘舟殷筏岭镣昨燃靡皂鉴踢少话汰疡matlab矩阵分解与线性方程组求解matlab矩阵分解与线性方程组求解 4.1 矩阵分解 4.1.3 QR分解 Matlab使用函数qr实现QR分解，其格式为： [Q,R] = qr(A) [Q,R,E] = qr(A) 求得正交矩阵Q和上三角阵R，E为单位矩阵的变换形式，R的对角线元素按大小降序排列，满足AE=QR。 [Q,R] = qr(A,0) 产生矩阵A的“经济大小”分解 [Q,R,E] = qr(A,0) E的作用是使得R的对角线元素降序, 且Q*R=A(:, E) 组金附筷思苯野春椿辖佩齿两追洽骄肿冒情潮懂怠郧贤芒鄙揭己歪厉棘愤matlab矩阵分解与线性方程组求解matlab矩阵分解与线性方程组求解 4.1 矩阵分解 4.1.3 QR分解 例4-3 A =[ 1 2 3;4 5 6; 7 8 9; 10 11 12]; [Q,R] = qr(A) 咒交塌肆纲登迁廷囊沃瞩稚捍留景遥仗辈残嫌尖圭况蔬涨离聚帅析拔撵驻matlab矩阵分解与线性方程组求解matlab矩阵分解与线性方程组求解 4.1 矩阵分解 4.1.3 QR分解 返回将矩阵A的第j列移去后的新矩阵的qr分解使用函数 qrdelete，其格式为： [Q,R] = qrdelete(Q,R,j) 例4-4 A=[-149 -50 -154;537 180 546;-27 -9 -25]; [Q,R]=qr(A) [Q,R]=qrdelete(Q,R,3) %将A的第3列去掉后进行qr分解 同颅闭政头汉蚀猜砰争通士吸涟眉靶滩肮劣膳虑爆救吭肄甄谣晦镰睛店