数学建模向量矩阵运算.ppt

,单击此处编辑母版标题样式,单击此处编辑母版文本样式,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,*,数学实验,向量与矩阵运算,向量与矩阵的生成,向量与矩阵运算,向量的生成,直接输入:,a=1,2,3,4,冒号,运,算符,a=1:4,=,a=1,2,3,4,b=0:pi/3:pi,=,b=0,1.0472,2.0944,3.1416,c=6:-2:0,=,c=6,4,2,0,例,：,从矩阵中抽取行或列,向量与矩阵的生成（续）,向量与矩阵运算,矩阵的生成,直接输入:,A,=1,2,3;4,5,6;7,8,9,由向量生成,由函数生成,通过编写,m,文件生成,例,：,x=1,2,3;y=2,3,4;,A=x,y,B=x;y,例,：,C=magic(3),常见矩阵生成函数,zeros,(m,n),生成一个,m,行,n,列的零矩阵，,m=n,时可简写为,zeros(,n,),ones,(m,n),生成一个,m,行,n,列的元素全为 1 的矩阵,m,=,n,时可写为,ones(,n,),eye,(m,n),生成一个主对角线全为 1 的,m,行,n,列矩阵,m=n,时可简写为,eye(,n,)，,即为,n,维单位矩阵,diag,(X),若,X,是矩阵，则,diag(X),为,X,的主对角线向量,若,X,是向量，,diag(X),产生以,X,为主对角线的对角矩阵,tril,(A),提取一个矩阵的下三角部分,triu,(A),提取一个矩阵的上三角部分,rand,(m,n),产生 01 间均匀分布的随机矩阵,m,=,n,时简写为,rand(,n,),randn,(m,n),产生均值为0，方差为1的标准正态分布随机矩阵,m,=,n,时简写为,randn(,n,),其它特殊矩阵生成函数：,magic,、,hilb,、,pascal,矩阵操作,提取矩阵的部分元素：冒号运算符,A(:),A,的所有元素,A(:,:),二维矩阵,A,的所有元素,A(:,k),A,的第,k,列，,A(k,:),A,的第,k,行,A(k:m),A,的第,k,到第,m,个元素,A(:,k:m),A,的第,k,到第,m,列组成的子矩阵,A(:),与,A(:,:),的区别,?,如何获得由,A,的第一、三行和第一、二列组成的子矩阵？,自己动手,矩阵操作,改变矩阵的形状：,reshape,reshape(A,m,n),:,将矩阵元素按,列方向,进行重组,重组后得到的新矩阵的元素个数,必须与原矩阵元素个数相等,！,矩阵操作,查看矩阵的大小：,size,size(A),列出矩阵,A,的,行数和列数,size(A,1),返回矩阵 A 的,行数,size(A,2,),返回矩阵 A 的,列,数,例,：,A=1 2 3;4 5 6,size(A),size(A,1),size(A,2),length(x),返回,向量,X,的,长度,length(A),等价于,max(size(A),矩阵基本运算,矩阵的加减,：,对应分量进行运算,要求参与加减运算的矩阵具有 相同的维数,例,：,A=1 2 3;4 5 6;B=3 2 1;6 5 4,C=A+B;D=A-B;,矩阵的普通乘法,要求参与运算的矩阵满足线性代数中矩阵相乘,的,原则,例,：,A=1 2 3;4 5 6;B=2 1;3 4;,C=A*B,矩阵基本运算,矩阵的,除法,：,/,、,右除和左除,若,A,可逆方阵，则,AB,A,的逆左乘,B,inv(A)*B,B/A,A,的逆右乘,B,B*inv(A),X=AB,A*X=B,X=B/A,X*A=B,通常，矩阵除法可以理解为,当,A,和,B,行数相等,时即可进行,左除,当,A,和,B,列数相等,时即可进行,右除,矩阵的乘方,A,是方阵，,p,是正整数,Ap,表示,A,的,p,次幂，即,p,个,A,相乘。,若,A,是方阵，,p,不是正整数,A,p,的计算涉及到,A,的特征值分解，即若,A=V*D*V,-1,则,A,p=V*(D.p)/V,矩阵的乘方,若,a,是标量，,A,是方阵，且,V,D=eig(A)，,则,aA,V*(aD)/V,若,A,P,均是矩阵，则,AP,无定义,若,a,是标量，,则,矩阵的,Kronecker 乘,积,矩阵,Kronecker,乘积,的定义,设,A,是,n,m,矩阵，,B,是,p,q,矩阵，则,A,与,B,的,kronecker,乘积为：,Kronecker,乘积,的性质,是,npmq,矩阵；,通常,任何两个矩阵都有,Kronecker,乘积,Matlab,中实现两个矩阵,Kronecker,相乘的函数为,kron(A,B),Kronecker,乘积有时也称张量积,矩阵的数组运算,数组运算：,对应元素进行运算,点与算术运算符之间不能有空格！,数组运算包括：,点乘,、,点除,、,点幂,相应的数组运算符为：,“.*”,，,“./”,，,“.”,和,“,.”,参与运算的对象必须具有相同的形状！,例,：,A=1 2 3;4 5 6;B=3 2 1;6 5 4;,C=A.*B;D=A./B;E=A.B;F=A.B;,函数取值,设,x,是变量，,f,是一个函数,当,x=a,是标量时，,f,(x)=f(a),也是一个标量,当,x=a,b,c,是向量时，,f,(x)=f(a),f(b),f(c),函数作用在矩阵上的取值,若,A,是矩阵，则,f,(A),是一个与,A,同形状的矩阵,f,作用在,x,的,每个分量上,函数取值,怎样计算,e,A,?,例,：,x=0:pi/4:pi;A=1 2 3;4 5 6;,y1=sin(x);y2=exp(A);y3=sqrt(A);,例,：,矩阵的超越函数,Matlab,提供了三种矩阵函数,：,expm、sqrtm、logm,详情参见联机帮助（,help expm/sqrtm/logm）,更一般的矩阵函数,：,funm,funm(,A,fun,),参数,fun,的可以是,exp,，log，cos，sin，cosh，sinh,矩阵与数的运算,加减,：,矩阵的每个元素都与数作加减运算,数乘,：,矩阵的每个元素都与数作乘法运算,矩阵除以一个数,：,每个元素都除以这个数,点幂,：,底为矩阵，指数为标量,底为标量，指数为矩阵,数与数组的点幂,x.y,=14,25,36=1,32,729,x.2,=12,22,32=1,4,9,2.x,=,?,例,：,x=1 2 3;y=4 5 6;,2.x;y,=,?,Matlab中的所有,标点符号必须在,英文状态下输入,Matlab,中常见数学函数,sin、cos、tan、cot、sec、csc、,asin、acos、atan、acot、asec、acsc,、,exp、,log,、log2、log10、sqrt,abs、conj、real、imag、sign,fix、floor、ceil、round、mod、rem,max,、,min,、,sum,、,mean,、,sort,、,det、inv、eig、rank、,log,是自然对数，即以,e,为底数,mod(x,y),结果与,y,同号，,rem(x,y),则与,x,同号,max,等函数的参数是矩阵时，是作用在矩阵各列上,上机作业,试分别生成,5,阶的单位阵、,8,阶均匀分布的随机矩阵及其下三角矩阵,生产列向量,x=1,3,5,7,9,29,生成以,x,的前,5,个元素为对角线的矩阵,A,生成一个与,A,同阶的正态分布的随机矩阵,B,输出,A,与,B,的,kronecker,乘积矩阵,C,生成由,B,与其同阶魔方矩阵点乘得到的矩阵,D,生成由,D,的第,2,、,4,、,5,行和第,4,、,1,列组成的子矩阵,E,求出矩阵,E,的最大元素,教材第,53,页，,1,、,2,、,3,、,4,、,5,