MATLAB矩阵及其运算21变量和数据操作22MATLAB矩阵.pptVIP

2.4.3 矩阵的逆与伪逆1．矩阵的逆对于一个方阵A，如果存在一个与其同阶的方阵B，使得：A·B=B·A=I (I为单位矩阵)则称B为A的逆矩阵，当然，A也是B的逆矩阵。 inv(A)函数求方阵A的逆矩阵。 用求逆矩阵的方法解线性方程组：Ax=b其解为：x=A-1b A = [1,-1,1;5,-4,3;2,1,1] A = 1 -1 1 5 -4 3 2 1 1 B = inv(A) B = -1.4000 0.4000 0.2000 0.2000 -0.2000 0.4000 2.6000 -0.6000 0.2000 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.0000 0.0000 -0.0000 1.0000 2．矩阵的伪逆如果矩阵A不是一个方阵，或者A是一个非满秩的方阵时，矩阵A没有逆矩阵，但可以找到一个与A的转置矩阵A‘同型的矩阵B，使得：A·B·A=AB·A·B=B此时称矩阵B为矩阵A的伪逆，也称为广义逆矩阵。pinv(A)求一个矩阵伪逆。 3.4.1 方阵的行列式值把一个方阵看作一个行列式，并对其按行列式的规则求值。 det(A)求方阵A所对应的行列式的值。 A = [1 2 3;3 2 1; 3 1 2] A = 1 2 3 3 2 1 3 1 2 det(A) ans = -12 2.4.5 矩阵的秩与迹1．矩阵的秩矩阵线性无关的行数与列数称为矩阵的秩。 rank(A)求矩阵秩 。2．矩阵的迹矩阵的迹等于矩阵的对角线元素之和，也等于矩阵的特征值之和。 trace(A)求矩阵的迹。 1．向量的3种常用范数及其计算函数在MATLAB中，求向量范数的函数为：(1) norm(V)或norm(V,2)：计算向量V的2—范数。(2) norm(V,1)：计算向量V的1—范数。(3) norm(V,inf)：计算向量V的∞—范数。 2．矩阵的范数及其计算函数MATLAB提供了求3种矩阵范数的函数，其函数调用格式与求向量的范数的函数完全相同。 2.4.7 矩阵的条件数在MATLAB中，计算矩阵A的3种条件数的函数是：(1) cond(A,1) 计算A的1—范数下的条件数。(2) cond(A)或cond(A,2) 计算A的2—范数数下的条件数。(3) cond(A,inf) 计算A的 ∞—范数下的条件数。 2.4.8 矩阵的特征值与特征向量对于n阶方阵A，求数λ 和向量ζ，使得等式Aζ = λζ 成立，则满足等式的λ 成为A的特征值， 向量ζ为A的特征向量 (1) E=eig(A)：求矩阵A的全部特征值，构成向量E。(2) [V,D]=eig(A)：求矩阵A的全部特征值，构成对角阵D，并求A的特征向量构成V的列向量。 (3) [V,D]=eig(A,‘nobalance’)：与第2种格式类似，但第2种格式中先对A作相似变换后求矩阵A的特征值和特征向量，而格式3直接求矩阵A的特征值和特征向量。 A = 1.0000 1.0000 0.5000 1.0000 1.0000 0.2500 0.5000 0.2500 2.0000 [V,D] = eig(A) V = 0.7212 0.4443 0.5315 -0.6863 0.5621 0.4615 -0.0937 -0.6976 0.7103 D = -0.0166 0 0 0 1.4801 0 0 0 2.5365 3.6 矩阵的超越函数 超越函数直接作用于矩阵，输入参数A必须是方阵1．矩阵平方根sqrtmsqrtm(A)计算矩阵A的平方根。2．矩阵对数logmlogm(A)计算矩阵A的自然对数。 3．矩阵指数expm、expm1、expm2、expm3expm(A)、expm1(A)、expm2(A)、expm3(A)的功能都求矩阵指数eA。 4．普通矩阵函