第二讲MATLAB的数值计算.pptVIP

矩阵与常数的运算 常数与此矩阵的各元素之间进行运算。 注意：进行数除时，常数通常只能做除数。 矩阵的逆运算 函数 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 矩阵的行列式运算 函数 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 矩阵的幂运算 与数字的幂运算形式相同，用“^”算符。 矩阵的指数运算 常用函数 expm expm1 expm2 expm3 矩阵的对数运算 函数 logm 矩阵的开方运算 函数 sqrtm 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 矩阵的基本函数运算 特征值函数 函数 [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 奇异值函数 函数 svd svds 矩阵非零特征值的平方根称作矩阵的奇异值。 矩阵的奇异值大小通常决定矩阵的性态，如果矩阵的奇异值变化特别大，则矩阵中某个参数的微小变化将严重影响到原矩阵的参数，这样的矩阵又称为病态矩阵，也称奇异矩阵。 例：A=[1,1;5*eps,0;0,5*eps]; sqrt(eig(A’*A)) 矩阵翻转 函数 fliplr flipud rot90 a = 7 3 -2 3 4 -1 -2 -1 3 fliplr(a) ans = -2 3 7 -1 4 3 3 -1 -2 flipud(a) ans = -2 -1 3 3 4 -1 7 3 -2 rot90(a) ans = -2 -1 3 3 4 -1 7 3 -2 条件数函数 判断矩阵的“病态”程度。 函数cond 计算矩阵的条件数的值 condest计算矩阵的1范数条件数的估计值 rcond 计算矩阵的条件数的倒数值 矩阵的条件数cond{A}=最大奇异值/最小奇异值，条件数越