第2章节 矩阵及其运算.pptVIP

fix朝零方向取整，如fix(-1.3)=-1; fix(1.3)=1;floor朝负无穷方向取整，如floor(-1.3)=-2; floor(1.3)=1;ceil朝正无穷方向取整，如ceil(-1.3)=-1; ceil(1.3)=2;round四舍五入到最近的整数，如round(-1.3)=-1;round(-1.52)=-2;round(1.3)=1;round(1.52)=2 rem(3,2)=1 rem(-3,-2)=-1 rem(3,-2)=1 rem(-3,2)=-1 mod(3,2)=1 mod(-3,-2)=-1 mod(3,-2)=-1 mod(-3,2)=1 由此可以看出，rem和mod是有符号区别的！ 当除数与被除数的符号相同时，rem和mod的结果是完全相同的；当除数与被除数的符号不相同时，结果不同。 (1) 生成5阶随机方阵A。A=fix((90-10+1)*rand(5)+10) (2) 判断A的元素是否可以被3整除。 P=rem(A,3)==0其中，rem(A,3)是矩阵A的每个元素除以3的余数矩阵。此时，0被扩展为与A同维数的零矩阵，P是进行等于(==)比较的结果矩阵。 2.4.5 矩阵的秩与迹1．矩阵的秩矩阵线性无关的行数与列数称为矩阵的秩。在MATLAB中，求矩阵秩的函数是rank(A)。2．矩阵的迹矩阵的迹等于矩阵的对角线元素之和，也等于矩阵的特征值之和。在MATLAB中，求矩阵的迹的函数是trace(A)。 2.4.6 向量和矩阵的范数矩阵或向量的范数用来度量矩阵或向量在某种意义下的长度。范数有多种方法定义，其定义不同，范数值也就不同。 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 矩阵的特征值与特征向量在MATLAB中，计算矩阵A的特征值和特征向量的函数是eig(A)，常用的调用格式有3种：(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的特征值和特征向量。 例2-12 用求特征值的方法解方程。3x5-7x4+5x2+2x-18=0p=[3,-7,0,5,2,-18];A=compan(p); %A的伴随矩阵x1=eig(A) %求A的特征值x2=roots(p) %直接求多项式p的零点 2.5 矩阵的超越函数1．矩阵平方根sqrtmsqrtm(A)计算矩阵A的平方根。2．矩阵对数logmlogm(A)计算矩阵A的自然对数。此函数输入参数的条件与输出结果间的关系和函数sqrtm(A)完全一样 3．矩阵指数expm、expm1、expm2、expm3expm(A)、expm1(A)、expm2(A)、expm3(A)的功能都求矩阵指数eA。4．普通矩阵函数funmfunm(A,‘fun’)用来计算直接作用于矩阵A的由‘fun’指定的超越函数值。当fun取sqrt时，funm(A,‘sqrt’)可以计算矩阵A的平方根，与sqrtm(A)的计算结果一样。 2.6 字符串 在MATLAB中，字符串是用单撇号括起来的字符序列。MATLAB将字符串当作一个行向量，每个元素对应一个字符，其标识方法和数值向量相同。也可以建立多行字符串矩阵。 字符串是以ASCII码形式存储的。abs和double函数都可以用来获取字符串矩阵所对应的ASCII码数值矩阵。相反，char函数可以把ASCII码矩阵转换为字符串矩阵。 例2-13 建立一个字符串向量，然后对该向量做如下处理：(1) 取第1～5个字符组成的子字符串。(2) 将字符串倒过来重新排列。(3) 将字