MATLAB矩阵及其运算.pptVIP

(5) 伴随矩阵MATLAB生成伴随矩阵的函数是compan(p)，其中p是一个多项式的系数向量，高次幂系数排在前，低次幂排在后。例如，为了求多项式的x3-7x+6的伴随矩阵，可使用命令：p=[1,0,-7,6];compan(p) 第三十页，共八十一页。 (6) 帕斯卡矩阵 二次项(x+y)n展开后的系数随n的增大组成一个三角形表，称为杨辉三角形。由杨辉三角形表组成的矩阵称为帕斯卡(Pascal)矩阵。函数pascal(n)生成一个n阶帕斯卡矩阵。 第三十一页，共八十一页。 例2-7 求(x+y)5的展开式。在MATLAB命令窗口，输入命令： pascal(6)矩阵次对角线上的元素1,5,10,10,5,1即为展开式的系数。 第三十二页，共八十一页。 2.3 MATLAB运算 2.3.1算术运算1．基本算术运算 MATLAB的基本算术运算有：＋(加)、－(减)、*(乘)、/(右除)、\(左除)、^(乘方)。注意，运算是在矩阵意义下进行的，单个数据的算术运算只是一种特例。 第三十三页，共八十一页。 (1) 矩阵加减运算 假定有两个矩阵A和B，则可以由A+B和A-B实现矩阵的加减运算。运算规则是：若A和B矩阵的维数相同，则可以执行矩阵的加减运算，A和B矩阵的相应元素相加减。如果A与B的维数不相同，则MATLAB将给出错误信息，提示用户两个矩阵的维数不匹配。 第三十四页，共八十一页。 (2) 矩阵乘法  假定有两个矩阵A和B，若A为m×n矩阵，B为n×p矩阵，则C=A*B为m×p矩阵。 第三十五页，共八十一页。 (3) 矩阵除法 在MATLAB中，有两种矩阵除法运算：\和/，分别表示左除和右除。如果A矩阵是非奇异方阵，则A\B和B/A运算可以实现。A\B等效于A的逆左乘B矩阵，也就是inv(A)*B，而B/A等效于A矩阵的逆右乘B矩阵，也就是B*inv(A)。 对于含有标量的运算，两种除法运算的结果相同，如3/4和4\3有相同的值，都等于0.75。又如，设a=[10.5,25]，则a/5=5\a=[2.1000 5.0000]。对于矩阵来说，左除和右除表示两种不同的除数矩阵和被除数矩阵的关系。对于矩阵运算，一般A\B≠B/A。 第三十六页，共八十一页。 (4) 矩阵的乘方 一个矩阵的乘方运算可以表示成A^x，要求A为方阵，x为标量。 第三十七页，共八十一页。 2．点运算 在MATLAB中，有一种特殊的运算，因为其运算符是在有关算术运算符前面加点，所以叫点运算。点运算符有.*、./、.\和.^。两矩阵进行点运算是指它们的对应元素进行相关运算，要求两矩阵的维参数相同。 第三十八页，共八十一页。 2.3.2 关系运算 MATLAB提供了6种关系运算符：(小于)、=(小于或等于)、(大于)、=(大于或等于)、==(等于)、～=(不等于)。它们的含义不难理解，但要注意其书写方法与数学中的不等式符号不尽相同。 第三十九页，共八十一页。 关系运算符的运算法则为： (1) 当两个比较量是标量时，直接比较两数的大小。若关系成立，关系表达式结果为1，否则为0。 (2) 当参与比较的量是两个维数相同的矩阵时，比较是对两矩阵相同位置的元素按标量关系运算规则逐个进行，并给出元素比较结果。最终的关系运算的结果是一个维数与原矩阵相同的矩阵，它的元素由0或1组成。 第四十页，共八十一页。 (3) 当参与比较的一个是标量，而另一个是矩阵时，则把标量与矩阵的每一个元素按标量关系运算规则逐个比较，并给出元素比较结果。最终的关系运算的结果是一个维数与原矩阵相同的矩阵，它的元素由0或1组成。 第四十一页，共八十一页。 例2-8 产生5阶随机方阵A，其元素为[10,90]区间的随机整数，然后判断A的元素是否能被3整除。 (1) 生成5阶随机方阵A。A=fix((90-10+1)*rand(5)+10) (2) 判断A的元素是否可以被3整除。 P= rem(A,3)= =0其中，rem(A,3)是矩阵A的每个元素除以3的余数矩阵。此时，0被扩展为与A同维数的零矩阵，P是进行等于(= =)比较的结果矩阵。 第四十二页，共八十一页。 2.3.3 逻辑运算 MATLAB提供了3种逻辑运算符：(与)、|(或)和～(非)。 逻辑运算的运算法则为： (1) 在逻辑运算中，确认非零元素为真，用1表示，零元素为假，用0表示。 (2) 设参与逻辑运算