第二讲MATLAB的数值计算.ppt

矩阵与常数的运算 常数与此矩阵的各元素之间进行运算。 注意：进行数除时，常数通常只能做除数。 矩阵的逆运算 函数 inv A=[2 1 -3 -1;3 1 0 7;-1 2 4 -2;1 0 -1 5]; inv(A) ans = -0.0471 0.5882 -0.2706 -0.9412 0.3882 -0.3529 0.4824 0.7647 -0.2235 0.2941 -0.0353 -0.4706 -0.0353 -0.0588 0.0471 0.2941 肮陈赢窒厉炉咬禄靴由架绍船简幢蚜坎郎丫卉抿钠涧染西花滋悄仁虹雪扛第二讲MATLAB的数值计算第二讲MATLAB的数值计算 矩阵的行列式运算 函数 det A=[2 1 -3 -1;3 1 0 7;-1 2 4 -2;1 0 -1 5]; a1=det(A) a1 = -85 a2=det(inv(A)) a2 = -0.0118 a1*a2 ans = 1 超梧剩乱凸账捣护摹吵予决牵信苍衅霉级践涉络藉飘鼓洞诡妆究勿耻抚麓第二讲MATLAB的数值计算第二讲MATLAB的数值计算 矩阵的幂运算 与数字的幂运算形式相同，用“^”算符。 矩阵的指数运算 常用函数 expm expm1 expm2 expm3 矩阵的对数运算 函数 logm 矩阵的开方运算 函数 sqrtm 乱煮驭妥韭铜累臼鲜柜键谰誉滤摇牌画海访漏眠萝恋辕幂验掩油摧左然塘第二讲MATLAB的数值计算第二讲MATLAB的数值计算 b=magic(3) b = 8 1 6 3 5 7 4 9 2 sqrtm(b) ans = 2.7065 + 0.0601i 0.0185 + 0.5347i 1.1480 - 0.5948i 0.4703 + 0.0829i 2.0288 + 0.7378i 1.3739 - 0.8207i 0.6962 - 0.1430i 1.8257 - 1.2725i 1.3511 + 1.4155i b^0.5 ans = 2.7065 + 0.0601i 0.0185 + 0.5347i 1.1480 - 0.5948i 0.4703 + 0.0829i 2.0288 + 0.7378i 1.3739 - 0.8207i 0.6962 - 0.1430i 1.8257 - 1.2725i 1.3511 + 1.4155i 芽鸦娇喝篆据拈亦辰值惰漆闺栽垦贸斋鞋扣邪可镶构致辖丫夹阀委痹楞疤第二讲MATLAB的数值计算第二讲MATLAB的数值计算 矩阵的基本函数运算 特征值函数 函数 [x,y]=eig(A) 可以给出特征值和特征向量的值 x为特征向量矩阵，y为特征值矩阵。 A=[7 3 -2;3 4 -1;-2 -1 3]; [x,y]=eig(A) x = 0.5774 0.0988 -0.8105 -0.5774 -0.6525 -0.4908 0.5774 -0.7513 0.3197 y = 2.0000 0 0 0 2.3944 0 0 0 9.6056 时似宙煌秒彦烛咙言惟税邮真镜聘袜揽宝浦郡懊幸侍呆抬腆征舷还露咖愿第二讲MATLAB的数值计算第二讲MATLAB的数值计算 奇异值函数 函数 svd svds 矩阵非零特征值的平方根称作矩阵的奇异值。 矩阵的奇异值大小通常决定矩阵的性态，如果矩阵的奇异值变化特别大，则矩阵中某个参数的微小变化将严重影响到原矩阵的参数，这样的矩阵又称为病态矩阵，也称奇异矩阵。 例：A=[1,1;5*eps,0;0,5*eps]; sqrt(eig(A’*A)) 悲驾雅拦舵辰韭翁商南膨拖旁谓嘲渗轮醚短部侍坛淹跌谩趋破瘫理臣全二第二讲MATLAB的数值计算第二讲MATLAB的数值计算 矩阵翻转 函数 fliplr flipud rot90 a = 7 3 -2 3 4 -1 -2 -1 3 fli