chapter3 MATLAB数值运算.pptx

2016-12-15基于MATLAB的控制系统计算机辅助设计与仿真1第3章MATLAB数值运算基于MATLAB的控制系统计算机辅助设计与仿真2016-12-15基于MATLAB的控制系统计算机辅助设计与仿真23.1矩阵和数组MATLAB具有强大的数值运算能力。数值运算主要指数组、矩阵和多项式的运算。MATLAB数值运算的表达式变量中不允许有未定义的自由变量。在MATLAB语言中，二维矩阵和数组在形式上一致，但它们的运算规则不同。矩阵运算是按照线性代数的规则进行运算。数组运算是两个数组按元素位置进行的逐元运算，即元素对元素的运算，数组中的元素都是平等的，要求被操作的两个数组中元素的行数和列数必须相等。数组运算符是在矩阵运算符前面加一点“.”（加减法除外），如“.*”、“./”和“.＾”等。向量是数组的一种特例，一行多列的数组为行向量，一列多行的数组为列向量。3.1.1矩阵的代数运算MATLAB算术运算（Arithmeticoperators）符关于算术运算符的用法可参见“docarithmeticoperators”html帮助。2016-12-15基于MATLAB的控制系统计算机辅助设计与仿真3MatrixOperatorsArrayOperators符号函数功能符号函数功能+Plus加+Unaryplus一元加法-Minus减-Unaryminus一元减法ctranspose矩阵转置.transpose数组转置*mtimes矩阵乘.*times数组乘\mldivide矩阵左除.\ldivide数组左除/mrdivide矩阵右除./rdivide数组右除^mpower矩阵乘方.^power数组乘方2016-12-15基于MATLAB的控制系统计算机辅助设计与仿真4常用的十进制符号如小数点、负号等，在MATLAB中也可以同样使用。表示10的幂次用符号e或E，如：3、-99、0.0001、1.6e-20、6.2e23。2016-12-15基于MATLAB的控制系统计算机辅助设计与仿真5 1.矩阵转置（共轭转置）运算符：n×m矩阵A的转置矩阵B定义为bji=aij，i=1,…,n，j=1,…,m，B为m×n矩阵，记为B=AT。如果A为复数矩阵，则A为复数矩阵的共轭转置（Hermit转置）。例： a=[2+3i4-4i5-2i;7+5i6-9i4+3i] b=a %矩阵转置 c=a. %数组转置 d=[245;764] e=d %矩阵转置 f=d. %数组转置2016-12-15基于MATLAB的控制系统计算机辅助设计与仿真6 2.矩阵的加、减运算符：＋、－只有同阶的两个矩阵或者一个矩阵和一个标量才能进行加减运算。两个同阶矩阵的加减运算是将对应元素相加减，而矩阵和标量的加减运算是将矩阵中每一个元素与标量进行加减运算。如果两个矩阵维数不匹配将提示出错信息例： A=[123;456],B=A+1 A= 123 456 B= 234 5672016-12-15基于MATLAB的控制系统计算机辅助设计与仿真7 3.矩阵的乘法（叉乘）运算符：*两矩阵的乘积为两个矩阵线性代数乘积，表示为C=A*B，其中前一矩阵的列数等于后一矩阵的行数，或其一为标量（称两矩阵的维数相容）。矩阵和标量相乘是将矩阵中的每个元素分别与标量相乘。矩阵乘法不满足交换律，即A*B≠B*A。例： A=[123;456;789];B=[41;75;03];C=3; D=A*B %A*B≠B*A E=A*C2016-12-15基于MATLAB的控制系统计算机辅助设计与仿真8 4.矩阵的除法在线性代数中，只有矩阵的逆运算，而没有矩阵的除法运算。在MATLAB中，矩阵的除法用来求线性方程的解。MATLAB定义了矩阵的左除“\”和右除“/”两种除法运算，表示了两种不同的除数矩阵和被除数矩阵的关系。对于标量，两种除法运算的结果相同，如1/4和4\1有相同的值0.25。2016-12-15基于MATLAB的控制系统计算机辅助设计与仿真9（1）矩阵左除运算符：\ X=A\B是线性方程组A*X=B的解。若A为非奇异方阵，则X=A\B=A-1*B=inv(A)*B；当方程是欠定或超定情况时，A\B所对应的是方程最小二乘解；若X不存在或不唯一，则产生一警告信息。矩阵A可以是矩形矩阵（即非方形矩阵），但此时方程组必须相容。 （2）矩阵右除运算符：/ X=B/A是线性方程组X*A=B的解。A为非奇异方阵时，X=B/A=B*A-1=B*inv(A)。通常A\B≠B/A2016-12-15基于MATLAB的控制系统计算机辅助设计与仿真10inv函数求矩阵的逆附录1：奇异矩阵定义： 行列式等于0的矩阵称为奇异矩阵。奇异矩阵的判断方法：首先，看这个矩阵是不是方阵（方阵是行数和列