第2章 矩阵运算基础.ppt

虽然逻辑数组的元素为0或1，但它不同与元素为0或1的“0-1数值矩阵” , “0-1数值矩阵”不能起到对矩阵元素标识的作用 5、数值矩阵标识矩阵元素 　　数值矩阵D为m行n列，数值矩阵元素数值k是“1至m×n”之间的正整数，k值可以重复,同时母矩阵A与数据矩阵D结构相同 A(D)的结果是矩阵： “结果矩阵”对应位置元素是其在母矩阵中按单下标标识规定的第k个元素，k为数值矩阵对应位置元素的值 例：2-6_1 例：2-6_1 A=zeros(4,4) A(:)=1:16 %全元素法赋值 L=abs(A)8 % L是“logical array” X=A(L) %将A中被矩阵L标识的所有元素赋给矢量X Dl=[1 9 1 1 1 9 1 1 1 1 5 14 1 1 5 14] Y=A(Dl) %数值矩阵定位 四、矩阵函数 m表示 matrix 矩阵a的对数log(a) logm(a) 矩阵a的逆 inv(a) 矩阵a的伪逆 pinv(a) 矩阵a的的秩 rank(a) 矩阵a的范数 矩阵a的指数ea 方阵a的行列式值 矩阵a的条件数 功 能 norm(a) expm(a) det(a) cond(a) 函数名称 【例2-19】矩阵a为[1 2 3;4 5 6;7 8 9]，计算a的逆矩阵和伪逆矩阵。 a=[1 2 3;4 5 6;7 8 9]; c1=inv(a) %逆 c2=pinv(a) %伪逆 d=det(a) %行列式的值 r=rank(a) %秩 n=norm(a) %范数 e=expm(a) %指数 l=logm(a) %对数 cd=cond(a) %条件数 通用函数(对矩阵元素的函数) 三角函数（数据对象是弧度） sin正弦 asin反正弦　 cos余弦　acos反余胘 tan正切　atan反正切　sec正割　asec反正割 如果x0是角度，则转为弧度：x*pi/180 加例2＿7：x=60 fd=x*pi/180 s=sin(fd) as=asin(s) c=cos(fd) ac=acos(c) t=tan(fd) at=atan(t) 指数函数 exp:以e为底的指数 pow2:2的幂次 log:自然对数　 　　log10:以10为底的对数 log2:以2为底的对数 sqrt:开平方 加例2＿8：t=3 et=exp(t) pt=pow2(t) lt=log(et) l10=log10(1000) l2=log2(pt) sq=sqrt(9) 复数函数 abs：绝对值，复数的模 angle：相位角 conj：复数的共轭 complex(a,b)=a+j*b=a+i*b 构造复数 real：复数实部 imag：复数虚部 isreal：判断对象是否为实数 ,结果为逻辑值 加例2＿9： x=3; y=6 z=complex(x,y) zab=abs(z) zan=angle(z) zconj=conj(z) zi=imag(z) zr=real(z) 取整函数 fix：朝零方向取整 floor：朝负无穷大方向取整，意：地板 ceil：朝正无穷大方向取整，意：天花板 round：四舍五入到最近的整数 sign：符号函数，在信号处理中常用 mod：求余数, mod(x,y)=x -floor(x./y).*y rem：除后求余数，与平时结果相同，是整数也可能是小数 加例2＿10 : p=3.7; n=-3.7 fp=fix(p); fn=fix(n) flp=floor(p); fln=floor(n) cp=ceil(p); cn=ceil(n) rp=round(p); rn=round(n) m1=mod(9,2); m11=9 -floor(9./2).*2 m2=mod(-9,2); m22=-9 -floor(-9./2).*2 r1=rem(9,2); r2=rem(-9,2);r3=rem(9.3,2.1) s1=sign(p); s2=sign(0); s3=sign(n) 五、矩阵分解 特征值分解 　矩阵A的特征矢量ν与特征值λ满足Aν=λν；将特征矢量ν合成矩阵V，特征值λ矩阵Δ对象角上，则AV＝VΔ。 [v,d]=eig(a),v是a的特征矢量矩阵，d是a的特征值矩阵，a×v=v×d d=eig(a),d是a的特征值矢量 [v,d]=eig(a,b),a和b的广义特征值