matlab第5讲byXR.pptVIP

矩阵的超越函数 矩阵对数logm，计算矩阵A的自然对数矩阵。 A=[4,9;1,5]; L=logm(A) L = 1.0639 2.4308 0.2701 1.3340 矩阵的超越函数 矩阵指数expm，计算矩阵指数 。 B=expm(L) B = 4.0000 9.0000 1.0000 5.0000 Ps：expm函数与logm函数是互逆的。 矩阵的超越函数 通用矩阵函数funm，funm(A,’fun’)对方阵A计算由fun定义的函数的矩阵函数值。 例如：fun是exp时，funm(A,’exp’)与expm(A)计算结果是一致的。 矩阵的超越函数 A=[2,-1;1,0] A = 2 -1 1 0 funm(A,exp) ans = 5.4366 -2.7183 2.7183 0.0000 expm(A) ans = 5.4366 -2.7183 2.7183 0.0000 由m个方程n个未知量构成的线性方程组 * 对于后两种情形的解法在学高等代数时已经介绍过 * null是用来求齐次线性方程组的基础解系的，加上r则求出的是一组最小正整数解，如果不加，则求出的是解空间的规范正交基。 * 3.4 矩阵求值 3.4.1 方阵的行列式把一个方阵看作一个行列式，并对其按行列式的规则求值，这个值就称为所对应的行列式的值。在MATLAB中，求方阵A所对应的行列式的值的函数是det(A)。 A=rand(5) A = 0.8147 0.0975 0.1576 0.1419 0.6557 0.9058 0.2785 0.9706 0.4218 0.0357 0.1270 0.5469 0.9572 0.9157 0.8491 0.9134 0.9575 0.4854 0.7922 0.9340 0.6324 0.9649 0.8003 0.9595 0.6787 B=det(A) B = -0.0250 3.4.2 矩阵的秩与迹1．矩阵的秩矩阵线性无关的行数与列数称为矩阵的秩。在MATLAB中，求矩阵秩的函数是rank(A)。 A=[2,2,-1,1;4,3,-1,2;8,5,-3,4;3,3,-2,2]; r=rank(A) r = 4 求矩阵 的秩。 解 可见R（B）= 2 ， 所以R（A）= 2 。 例 2．矩阵的迹矩阵的迹等于矩阵的对角线元素之和，也等于矩阵的特征值之和。在MATLAB中，求矩阵的迹的函数是trace(A)。 A=[2,2,3;4,5,-6;7,8,9]; trace(A) ans = 16 3.4.3 向量和矩阵的范数矩阵或向量的范数用来度量矩阵或向量在某种意义下的长度。范数有多种方法定义，其定义不同，范数值也就不同。 1．向量的3种常用范数及其计算函数在MATLAB中，求向量范数的函数为：(1) norm(V)或norm(V,2)：计算向量V的2—范数。(2) norm(V,1)：计算向量V的1—范数。(3) norm(V,inf)：计算向量V的∞—范数。 常用的向量范数: 范数意义下的单位向量: X=[x1, x2]T 1 -1 1 ||X||1 = 1 1 1 -1 -1 ||X||2 = 1 -1 1 1 -1 -1 ||X||∞ = 1 向量V的2—范数 公式： V=[-1,1/2,1]; v1=norm(V,1) v1 = 2.5000 向量V的1—范数 公式： V=[-1,1/2,1]; v2=norm(V) v2 = 1.5000 向量V的∞—范数 公式： V=[-1,1/2,1]; vinf=norm(V,inf) vinf = 1 2．矩阵的范数及其计算函数MATLAB提供了求3种矩阵范数的函数，其函数调用格式与求向量的范数的函数完全相同。 设A是一个m*n的矩阵，V是一个含有n个元素的列向量，定义： 称为A从属于向量的范数。 A=[17,0,1,0,15;23,5,7,14,16;4,0,13,0,22;10,12,19,21,3;11,18,25,2,19]; a1=norm(A,1) a1 = 75 其中λ1为A’A最大特征值 A=[17,0,1,0,