第二讲 向量与矩阵运算精选文档.ppt

《第二讲 向量与矩阵运算精选文档.ppt》由会员分享，可在线阅读，更多相关《第二讲 向量与矩阵运算精选文档.ppt（18页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第二讲第二讲向量与矩阵运算向量与矩阵运算本讲稿第一页，共十八页a=1:1:4 =a=1,2,3,4b=0:pi/3:pi =b=0,1.0472,2.0944,3.1416c=6:-2:0 =c=6,4,2,0例例：（2）冒号表达式生成）冒号表达式生成利用冒号表达式自动建立一个向量利用冒号表达式自动建立一个向量:冒号表达式可以自动产生一个行向量，一般格式是：冒号表达式可以自动产生一个行向量，一般格式是：X=e1:e2:e3 其中其中e1为初始值，为初始值，e2为步长，为步长，e3为终止值。为终止值。本讲稿第二页，共十八页 （3）用）用linspace函数产生行向量。函数产生行向量。其调用格式为

2、：其调用格式为：linspace(a,b,n)其中其中a和和b是生成向量的第一个和最后一个元素，是生成向量的第一个和最后一个元素，n是元素总数是元素总数(采用等分）采用等分）。a=linspace(1,3,3)a=1,2,3本讲稿第三页，共十八页q 矩阵（矩阵（m行行n列的数表）的生成列的数表）的生成向量与矩阵运算向量与矩阵运算(1)直接输入直接输入:两行之间用;或回车隔开 元素之间用空格或,隔开.A=1,2,3;4,5,6;7,8,9(2)由向量生成由向量生成,每一行为一个向量每一行为一个向量例例：x=1,2,3;y=2,3,4;A=x;yA=X,Y 等于什么等于什么?本讲稿第四页，共十八页

3、(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 维单位矩阵常见矩阵生成函数常见矩阵生成函数本讲稿第五页，共十八页矩阵元素的操作矩阵元素的操作 q 提取矩阵的部分元素：提取矩阵的部分元素：冒号运算符冒号运算符u A(:,k)A的第的第 k 列，列，A(k,:)A的第的第 k 行行 u A(:,k:m)A的第的第

4、k 到第到第 m 列组成的子矩阵列组成的子矩阵u A(k:m,:)A的第的第 k 到第到第 m 行组成的子矩阵行组成的子矩阵u双下标双下标:A(i,j)位于位于A的第的第i行第行第j列的元素列的元素u 单下表单下表:A(i)矩阵矩阵A的第的第i个元素（按列从左到右排列个元素（按列从左到右排列)q元素的标识元素的标识u A(:,k，m)A的第的第 k、m列，列，本讲稿第六页，共十八页向量与矩阵运算向量与矩阵运算na=123n456n789na=n123n456n789na(1,3,:)%第一、三行第一、三行nans=n123n789本讲稿第七页，共十八页矩阵元素的操作矩阵元素的操作 u A(k:

5、m)A的第的第 k 到第到第 m 个元素个元素u删除删除A的第的第km行，构成新矩阵行，构成新矩阵:A(k:m，：，：)=删除删除A的第的第km列，构成新矩阵列，构成新矩阵:A(:,k:m)=将矩阵将矩阵A和和B拼接成新矩阵：拼接成新矩阵：AB；A；B本讲稿第八页，共十八页二、矩阵基本运算二、矩阵基本运算q 矩阵的加减矩阵的加减：对应分量进行运算对应分量进行运算要求参与加减运算的矩阵具有 相同的维数例例：A=1 2 3;4 5 6;B=3 2 1;6 5 4 C=A+B;D=A-B;q 矩阵的普通乘法矩阵的普通乘法要求参与运算的矩阵满足线性代数中矩阵相乘的原则例例：A=1 2 3;4 5 6;

6、B=2 1;3 4;5 6;C=A*B本讲稿第九页，共十八页矩阵基本运算矩阵基本运算q 矩阵的除法矩阵的除法：/、右除和左除右除和左除 若 A 可逆方阵，则AB A 的逆左乘的逆左乘 B =inv(A)*BB/A A 的逆右乘的逆右乘 B B*inv(A)X=AB A*X=BX=B/A X*A=B 通常，矩阵除法可以理解为本讲稿第十页，共十八页矩阵的乘方矩阵的乘方u A 是方阵，p 是正整数Ap 表示 A 的 p 次幂，即 p 个 A 相乘。本讲稿第十一页，共十八页三、矩阵的数组运算三、矩阵的数组运算q 数组运算：数组运算：对应元素进行运算点与算术运算符之间不能有空格！u 数组运算包括：数组运

7、算包括：点乘点乘、点除点除、点幂点幂u 相应的数组运算符为：相应的数组运算符为：“.*”，“./”，“.”和和“.”参与运算的对象必须具有相同的形状！对应元素进行相应的运算例例：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 是一个函数是一个函数u 当当 x=a 是标量时，是标量时，f(x)=f(a)也是一个标量也是一个标量u 当当 x=a,b,c 是向量时，是向量时，f(x)=f(a),f(b),f(c)q 函数作用在矩阵上的取值函数作用在矩阵上的取值u 若若 A

8、 是矩阵，则是矩阵，则 f(A)是一个与是一个与 A 同形状的矩阵同形状的矩阵 f 作用在作用在 x 的的每个分量上每个分量上本讲稿第十三页，共十八页函数取值函数取值例例：x=0:pi/4:pi;A=1 2 3;4 5 6;y1=sin(x);y2=exp(A);y3=sqrt(A);例例：本讲稿第十四页，共十八页1.矩阵的逆矩阵的逆 对于一个方阵对于一个方阵A，如果存在一个与其同阶的方阵，如果存在一个与其同阶的方阵B，使得：使得：AB=BA=I(I为单位矩阵为单位矩阵),则称则称B为为A的逆矩阵，的逆矩阵，当然，当然，A也是也是B的逆矩阵。的逆矩阵。求一个矩阵的逆是一件非常烦琐的工作，容易出

9、错，求一个矩阵的逆是一件非常烦琐的工作，容易出错，但在但在MATLAB中，求一个矩阵的逆非常容易。中，求一个矩阵的逆非常容易。求方阵求方阵A的逆矩阵可调用函数的逆矩阵可调用函数inv(A)。例例 用求逆矩阵的方法解线性方程组。用求逆矩阵的方法解线性方程组。Ax=b 其解为：其解为：x=A-1b本讲稿第十五页，共十八页a=1 2 3;4 5 6;2 3 5;b=inv(a)b=-2.3333 0.3333 1.0000 2.6667 0.3333 -2.0000 -0.6667 -0.3333 1.0000本讲稿第十六页，共十八页2 方阵的行列式方阵的行列式 把一个方阵看作一个行列式，并对其按行

10、列式的规把一个方阵看作一个行列式，并对其按行列式的规则求值，这个值就称为矩阵所对应的行列式的值。则求值，这个值就称为矩阵所对应的行列式的值。在在MATLAB中，求方阵中，求方阵A所对应的行列式的值的函数是所对应的行列式的值的函数是det(A)。3 矩阵的秩矩阵的秩 矩阵线性无关的行数与列数称为矩阵的秩。在矩阵线性无关的行数与列数称为矩阵的秩。在MATLAB中，求矩阵秩的函数是中，求矩阵秩的函数是rank(A)。本讲稿第十七页，共十八页4.矩阵的特征值与特征向量矩阵的特征值与特征向量 在在MATLAB中，计算矩阵中，计算矩阵A的特征值和特征向量的函数是的特征值和特征向量的函数是eig(A)，常用的调用格式有：，常用的调用格式有：(1)E=eig(A)：求矩阵：求矩阵A的全部特征值，构成向量的全部特征值，构成向量E。(2)V,D=eig(A)：求矩阵：求矩阵A的全部特征值，构成对角阵的全部特征值，构成对角阵D，并求，并求A的特征向量构成的特征向量构成V的列向量。的列向量。对于奇异矩阵或长方阵对于奇异矩阵或长方阵b，把同时满足，把同时满足xbx=b和和bxb=x（penrose第一，第二方程）的矩阵第一，第二方程）的矩阵x称为称为b的伪矩阵，的伪矩阵，matlab中求伪矩阵中求伪矩阵x的命令为的命令为pinv（b）；）；是方阵逆阵的推广。是方阵逆阵的推广。本讲稿第十八页，共十八页