《matlab矩阵运算和数组运算》由会员分享，可在线阅读，更多相关《matlab矩阵运算和数组运算(7页珍藏版)》请在人人文库网上搜索。
1、更新时间：-5-3matlab矩阵运算和数组运算作者：佚名教程来源：网络点击数：1368矩阵运算和数组运算是Matlab的数值运算中的两大类运算。矩阵运算是按矩阵运算法则进行的运算；数 组运算无论是何种运算操作都是对元素逐个进行。矩阵运算和数组运算指令对照汇总矩阵运算指令指令含义数组运算指令指令含义A*矩阵转置A.+B对应兀素相加A+B矩阵相加A.-B对应兀素相减A-B 矩阵相减A.*B冋维数组对应兀素相乘s+B标量加矩阵s.*A A的每个元素乘ss-B,B-s标量矩阵相减A./B A的元素被B的对应元素除A*B矩阵相乘BAA同上A/BA 右除 B s./B, BAss分别被B的兀素除。
2、BAA左除B A.An A的每个兀素自乘n次in v(A)矩阵求逆log(A)对A的每个元素求对数AAn矩阵的n次幕sqrt(A)对A的每个兀素求平方根f(A)求A的各个元素的函数值例：a=1 2 3; 4 5 6; 7 8 9;b=1 2 3;32 1;1 45;c=1 1 1;23 1;1 02;d=a*cA2+bd =3231368279821281291343.4矩阵函数和数组函数3.4.1基本数组函数数组函数是对各个元素的函数设计的。f(.)基本函数表函数名称功能函数名称功能sin 正弦acosh反双曲余弦COS余弦atanh反双曲正切tan正切acoth反双曲余切cot余切asec。
3、h反双曲正割sec正害acsch反双曲余割CSC余割fix朝零方向取整asin反正弦ceil朝正无穷大方向取整acos反余弦floor朝负无穷大方向取整atan反正切round四舍五入到整数atan2四象反正切rem除后取余数acot反余切sign符号函数asec反正割abs绝对值acsc反余割angle复数相角sinh双曲正弦imag复数虚部cosh双曲余弦real复数实部tanh双曲正切conj复数共馳coth双曲余切log 10常用对数sech双曲正割log自然对数csch双曲余割exp指数asinh反双曲正弦aqrt平方根f(.)特殊函数表函数名称功能函数名称功能bessel第一、第二。
4、类Bessel函数erfbeta Beta函数eerfinv逆误差函数误差函数gamma Gamma数ellipk第一、第二类全椭圆积分rat有理近似ellipj Jacobi椭圆函数3.4.2基本矩阵函数基本矩阵函数指令函数指令指令含义函数指令指令含义cond(A)矩阵的条件数(最大奇异值除以最小奇异值)svd(A)矩阵的奇异值分解det(A)方阵的行列式trace(A)矩阵的迹dot(A,B)矩阵的点积expm(A)矩阵指数eig(A)矩阵的特征值expm1(A)用Pade近似求norm(A,1)矩阵1 范数expm2(A)用Taylor级数近似求，精度稍差，但对任何方阵适用norm(A)。
5、矩阵的2-范数expm3(A)用矩阵分解求，仅当独立调整向量数目等于秩时适用norm(A,inf)矩阵的无穷范数logm(A)矩阵对数ln(A)norm(A；fro,)矩阵的f-范数(全部奇异值平方和的正平方根)sqrtm(A)平方根矩阵rank(A)矩阵的秩(非零奇异值的个数)rcond(A)矩阵的倒条件数funm(A,fn) A阵的一般矩阵函数例：注意观察奇异值与矩阵各性质的矢系a=magic(5);s=svd(a) d=det(a),t=trace(a),rk=rank(a),c=cond(a)n1=norm(a,1), n2=norm(a), ninf=n orm(a,inf),nf=。
6、norm(a,fro)7/765.0000 22.5471 21.6874 13.4036 11.9008 d = 5070000t =65rk =5c =5.4618n1 =65n2 =65.0000ninf =65nf =74.33033.5线性方程组的直接解法线性方程组Ax=b, A是的系数矩阵1)n=m且非奇异时，此方程称为“恰定”2)nm时，此方程称为“超定”方程(3)*时，此方程称为“欠定”方程(3.5.1矩阵逆和除法解恰定方程组1)采用求逆运算：x=inv(A)*b采用左除运算：x=Ab说明：方程(Properly Determined Equation)Overdetermin。
7、ed Equation)Underdetermined Equation1、由于MATLAB遵循IEEE算法，所以即使方面给出藩告:“ Warning:Matrix isA阵奇异，该运算也照样进行。但在运算结束时，singular to working precision ” ；另一方面，所得逆阵的元素都是“ Inf ”(无穷大)。1、当A为“病态”时，也给出警告信息。2、在MATLABK inv指令不很有用。MATLA睢荐：尽量使用除运算，少用逆运算。例1 : “求逆”法和“左除”法解恰定方程组的性能对比。为对比两种方法的性能，先用以下指令构造一个条件数很大的高阶恰定方程组。randCsee。
8、d*,);%选定随机种子，目的是可重复产生随机矩阵AA=rand(500)+1 .e8;%rand(500)生成均匀分布的随机矩阵%每个随机矩阵元素加一个数的目的是使A的条件数变大x=ones(500,1); %令解向量x为全1的500元列向量b=A*x;%为使Ax=b方程一致，用A和x生成b向量cond(A)%计算矩阵A的条件数ans =1.7608e+013过程1 : “求逆”法解，恰定方程组的误差、残差和所用计算时间tic%启动记时器(Stopwatch Timer)xi=inv(A)*b; %xi是用“求逆”法解恰定方程组所得的解toe%矢闭计时器，并显示解方程所用的时间eri=nor。
9、m(x-xi) %解向量xi与真解向量的2-范误差rei=norm(A*xi-b)/norm(b) %方程的2-范相对残差 elapsed_time =7.2500eri =0.0066rei =1.5488e-006过程2 : “左除”法解，恰定方程组的误差、残差和所用计算时间tic xd=Ab;% xd是用“左除”法解恰定方程组所得的解toe erd=norm(x-xd)red=no rm(A*xd-b)/no rm(b)elapsed_time =3.3500erd =0.00211.2695e-015说明：1) 计算结果表明：除法求解比求逆求解速度明显快，精度相当；但“除法”的相对残差。
10、几乎是“机器零”，而“逆 阵”法的相对残差高得多。2) MATLAB在设计求逆函数inv时，采用的是Gauss消去法。3) MATLAB在设计“左除”运算解恰定方程时，并不求逆，而是直接采用Gauss消去法求解，有效地减少了残差，所以即便在高条件数下也能得到较好的结果。3.5.2矩阵除法解超定方程组(1) 求正则方程(Normal equations)的解。(2) 用Householder变换(Householder transformation )直接求原超定方程的最小二 乘解。由于第二种方程法采用的是正交变换，据最小二乘理论可知，第二种方法所得的解的准确性、可靠性都比第一种方法好得多。MA。
11、TLAB解超定方程组用的就是第二种方法。例：除运算解超定方程的简单算例。a=1 2 3;4 5 6;7 8 9;10 11 12;b=1 ：4：x=abx = -0.33330.66670.00003.5.3矩阵除法解欠定方程组欠定方程的解是不唯一的。用除法运算所得的解有两个重要特征：(1 )在解中至多有Rank(A)个非零元素；(2)它是这类型解中范数最小的一个。例：除运算解欠定方程的简单算例。a=1 2 3;4 5 6;7 8 9;10 11 12;b=a：C=1 3 3：x=bcx =2.00000.16670-0.1667Zisserman Multiple View Geometry。
12、 in Computer Vision:44 单欣.基于RANSAC法的基本矩阵估计的匹配方法J.上海电 机学院学报，，9(4) : 666945 Xiaoli Liu ，Xiang Peng，Yongkai Yin. A coarse registration method of range image based on SIFTC . Conferenee on Electronic Imaging and Mult imedia Tech no logy，bei jing，46 Yanl i Wan，Zhenj iang Miao. Automatic panoramaimage mosaic and ghost el iminatingC.IEEE InternationalConferenee on Multimedia and Expo，Hannover，Germany，I CME，，v01. 3。
