ch3Matlab矩阵的运算.ppt

Matlab矩阵的运算,安徽工业大学数理学院侯为根wghou,1.数学运算,2.逻辑运算,三、Matlab基本运算,3.比较运算,4.位运算符,5.集合运算符,一、数学运算,MATLAB的基本算术运算有：(加)、(减)、*(乘)、/(右除)、(左除)、(乘方)、共轭转置（）正号（+）,负号（-）,【注意】运算是在矩阵意义下进行的，单个数据的算术运算只是一种特例。,2.点运算,点运算符有.*、./、.和.。两矩阵进行点运算是指它们的对应元素进行相关运算，要求两矩阵的维参数相同。,1.基本算术运算,1.转置（.）,幂次（.）.共轭转置（）,矩阵幂次（）2.正号（+）,负号（-）3.乘法（.*）,元素右除（./）,元素左除（.）,矩阵乘法（*）,矩阵右除（/）,矩阵左除（）4.加法（+）,减法（-）5.冒号（:）（例如：x=1:2:5等）,数学运算符之优先级:,同一类的运算符均具有相同的优先度(Priority),因此在计算上，是由左至右依次完成,矩阵的数学运算,矩阵的加减与一般标量（Scalar）的加减类似相加或相减的矩阵必需具有相同的维度,A=12345620;B=1324;C=A+BC1=A-B,例1:加减法运算,矩阵与纯量可以直接进行加减，MATLAB会直接将加减应用到每一个元素,A=12321+5,矩阵的加减法运算,矩阵的乘法与除法,标量对矩阵的乘或除，可比照一般写法,A=123;442;C=A/3B=2*A,欲进行矩阵相乘，必需确认第一个矩阵的列数（ColumnDimension）必需等于第二个矩阵的行数（RowDimension）,例2:矩阵的乘法,A=1;2;B=3,4,5;C=A*B,矩阵除法,若A为非奇异方阵，则X=A-1B,若A为非奇异方阵，则X=BA-1,矩阵左除：AX=B，求XMATLAB求解：X=AB,最小二乘解,矩阵右除：XA=B，求XMATLAB求解：X=B/A最小二乘解,注：若A为非奇异方阵，则A-1在Matlab中为inv(A),矩阵的次方运算,矩阵的次方运算，可由“”来达成，但矩阵必需是方阵，其次方运算才有意义,A=magic(3);B=A2,例3：,复数矩阵z，其“共轭转置”矩阵(ConjugateTranspose)可表示成矩阵z,转置和共轭转置矩阵,i=sqrt(-1);%单位虚数z=1+i,2;3,1+2i;w=z%共轭转置(注意z后面的单引号),例4:矩阵共轭转置,i=sqrt(-1);%单位虚数z=1+i,2;3,1+2i;w=z.%单纯转置(注意z后面的句点及单引号),想得到任何矩阵z的转置（Transpose），则可表示成矩阵矩阵z.,例5:矩阵的转置,若z为实数，则z和z.的结果是一样的,4、MATLAB矩阵运算函数,size(A)%矩阵的大小,inv(A)%矩阵的逆,det(A)%矩阵的行列式,rank(A)%矩阵的秩,rref(A)%矩阵列的极大无关组,orth(A)%矩阵化为正交阵,trace(A)%矩阵的迹,chol(A)%三角分解,A=L*LT,lu(A)%A对称，A=L*U,U,S,V=svd(A)%A=U*S*V,qr(A)%A=q*r(q正交阵,r三角阵）,norm(A,1|2|p|inf|fro)%矩阵的模,cond(A)%矩阵的条件数,4、MATLAB矩阵运算函数,eig(A)%A*x=*x,V,D=eig(A)%A*V=V*D,eig(A,B)%A*x=*B*x,V,D=eig(A,B)%A*V=B*V*D,例：设A=pascal(4),B=magic(4)，并使用下列函数进行计算,向量的p-norm,一个向量a的p-norm可以定义为,p=2时，此即为向量a的长度，或称欧氏长度（EuclideanLength）,欲求一向量的p-norm，可使用norm指令norm(x,p),a=34;x=norm(a,1)%x=7y=norm(a,2)%y=5z=norm(a,inf),注意：,矩阵A的p-norm可以定义如下：,例9:normMatrix,norm指令亦可用于计算矩阵的p-norm,A=123;456;789;norm(A,2),Cholesky分解是将一个对称的正定矩阵分解为一个上三角矩阵和其转置的乘积。,X=pascal(5)R=chol(X)C=R*R,LU分解是将原正方(square)矩阵分解成一个上三角形矩阵或是排列(permuted)的上三角形矩阵和一个下三角形矩阵，这样的分解法又称为三角分解。,命令格式为：,L,U=lu(A),L,U,P=lu(A),例如：,B=132;-2-61;257L,U=lu(B)L*U,L,U,P=lu(B)P*L*U,QR分解法是将矩阵分解成一个正规正交矩阵与上三角形矩阵，矩阵A不必为正方矩阵。,A=1325;-2-616;2572Q,R=qr(A)Q*QQ*R,奇异值分解(sigularvaluedecomposition,SVD)是另一种正交矩阵分解法；SVD是最可靠的分解法，但是它比QR分解法要花上近十倍的计算时间。,U,S,V=svd(A)，其中U和V代表二个相互正交矩阵，而S代表一对角矩阵。和QR分解法相同者，原矩阵A不必为正方矩阵。,A=1,2;3,1S,V,D=svd(A),Matlab关于向量运算的函数,x=randi(0,20,1,7),y1=sum(x)y2=sort(x),2、点运算,矩阵对应元素的直接运算,例如,A=123;456;780,B=A.A,a=1:5d1=1./a%求a中元素的倒数d2=a.2%求a中元素的平方d3=2.a%求2的a中元素的幂,例:设,x=1:100;s1=sum(-1).(x+1).*(1./x),x=1:100;s1=(-1).(x+1)*(1./x),方法1,方法2,3MATLAB数组运算函数,1.规则：设X=(xij)mn，则f(X)=(f(xij)mn,sincostancotseccscasinacosatanacotasecacsc,2.常用数学函数,三角与反三角函数,exploglog10sqrtpow2,ceilfixfloorroundremsign,absangleconjimagreal,指数与对数,复变函数,取整函数,cart2sphcart2polpol2cartsph2pol,坐标转换,注：若要对函数进行矩阵运算，方法为：funm(X,f),例：求sin(x)在x=0,/4,/2,3/4,时的值,x=0:pi/4:piy=x.*sin(x.2),例：求y=xsin(x2)在x=0,/4,/2,3/4,时的值,x=0:pi/4:piy=sin(x),矩阵的比较运算,=,ba=b,关系运算符的用法:,关系运算符的用法,A=123;456;78,0b=A5A(b),3、同维矩阵间比较；返回0-1矩阵,2、矩阵与一数比较返回0-1矩阵,例2:给出矩阵A中大与5的元素值,例3:接上例,设B=magic(3),B=magic(3);BA,find指令:可传回非零元素索引,显示一维索引或下标,元素索引,index=find(A5),可以将关系运算符返回的0-1矩阵,来进行矩阵的索（Indexing）求出在A矩阵中，满足AijBij的元素,x=012;003index=find(x)x(index),x(index)则是非零元素所形成的向量,要找出满足某条件的二维索引或下标,2维索引,idx1,idx2=find(x5),idx1和idx2分别是符合条件之矩阵元素的列索引(RowIndex)及行索引(ColumnIndex),例:如何找出介于7和11的元素?,x=magic(5)x(find(7=0).*y;%消去负半波p=sin(pi/3);y2=(y=p)*p+(yp).*y1;%按要求处理第(2)步,MATLAB的位运算符（BitwiseOperators）可对非负整数进行位对位的逻辑运算,位运算符,数字10的二进制表示法是1010，数字4的二进制表示法是0100，因此数字10和数字4的bitor应传回1110的十进制数值，即数字14,数字18的二进制表示法向左平移一个单位，其效果相当乘2,bitor(10,4),bitshift(18,1),MATLAB的集合运算符