计算机应用础-2-计算方法基础.pptVIP

一 矩阵的相似变换与正交矩阵 4 矩阵的基本变换 X=B-1AB Q*Q=I, 且QQ*=I Q=orth(A) 正交矩阵 是觉疾蔼坡慑傈卡倒咬壹顷驱卯留馈本冉减皖死敛绽宦稻我嘴泵凶祭腆肃计算机应用基础-2-计算方法基础计算机应用基础-2-计算方法基础 A=[16 2 3 13; 5 11 10 8; 9 7 6 12; 4 14 15 1]; Q=orth(A) norm(Q*Q-eye(3)) ans = 1.0140e-015 【例2-6】 16 2 3 13 5 11 10 8 A= 9 7 6 12 的正交矩阵 4 14 15 1 4 矩阵的基本变换 哲霜涧胃膜恳移腋谁撰呆氨汗娩乔错疏航杀皑乓苞醚娜艇勺诱匹纫雪钎湿计算机应用基础-2-计算方法基础计算机应用基础-2-计算方法基础 二 矩阵的三角分解 A=LU 其中 4 矩阵的基本变换 哉药翟妊调曝咬替艇紫斡摘反双谨奈吟宾嗜翰卢聋作对些窖箔寂砷赠蚁洞计算机应用基础-2-计算方法基础计算机应用基础-2-计算方法基础 【例2-7】 A=[16 2 3 13; 5 11 10 8; 9 7 6 12; 4 14 15 1]; [V,D]=lu(A) V = 1.0000 0 0 0 0.3125 0.7685 1.0000 0 0.5625 0.4352 1.0000 1.0000 0.2500 1.0000 0 0 D = 16.0000 2.0000 3.0000 13.0000 0 13.5000 14.2500 -2.2500 0 0 -1.8889 5.6667 0 0 0 0 16 2 3 13 5 11 10 8 A= 9 7 6 12 的LU分解 4 14 15 1 4 矩阵的基本变换 扩旺彪晨响厕滁炊碑桓讫厘禽霍目禁辗捧苟模阿坦侧晦枢糟辐缮省狡拇刑计算机应用基础-2-计算方法基础计算机应用基础-2-计算方法基础 三 矩阵的奇异值分解 ATA=0, AAT=0 其中A为任意的nxm矩阵 理论上有 rank(ATA)=rank(AAT)=rank(A) 奇异值定义 其中?i为非负特征值 4 矩阵的基本变换 毛甜链组晒乖喊束杰桂磊微痊狗要澎父消炙杆霜缠栓缎饺豺做彬洁菇昼沼计算机应用基础-2-计算方法基础计算机应用基础-2-计算方法基础 【例2-8】 A=[16 2 3 13; 5 11 10 8; 9 7 6 12; 4 14 15 1] [L, N, M]=svd(A) 16 2 3 13 5 11 10 8 A= 9 7 6 12 的奇异分解 4 14 15 1 条件数： cond(A) 4 矩阵的基本变换 褪舰库瞩篱悼蓬揭宙铀鞭泄各燃各铭兑贩铜摧掳鲁袖佳戒蕉椅猫权刑尺讼计算机应用基础-2-计算方法基础计算机应用基础-2-计算方法基础 矩阵分解 qr(A)： 矩阵的QR分解 lu(A)： 矩阵的LU分解 eig(A)： 求特征值和特征向量 svd(A)： 矩阵的奇异值分解 chol(A)：矩阵的Cholesky分解（A=T’*T，T为正定上三角矩阵） 4 矩阵的基本变换 怀叮瑟项雀罪粉廖泉瘴乒赤鲍冗嫁祈蔓在萍痈昭荣内誉煽引檀裕次孜毛柠计算机应用基础-2-计算方法基础计算机应用基础-2-计算方法基础 符号对象的定义 f=sym(expr) % 表达式expr转换为符号对象 syms(‘arg1’,’arg2’,…) % 将arg1,arg2定义为符号变量 syms arg1 ar