计算机模与拟第6课 .ppt

例2 已知矩阵A和B，分别计算A-B、 A+B、A+B’、A+6。 A=[1 2 3 4;5 6 7 8;9 10 11 12] B=[2 3 7 8;1 5 3 7;4 1 6 9] 注意：如果A和B中有一个是标量，则该标量与矩阵的每个元素进行加减运算。 3.2.3 特殊矩阵1．通用的特殊矩阵常用的产生通用特殊矩阵的函数有：zeros：产生全0矩阵(零矩阵)。ones：产生全1矩阵(幺矩阵)。eye：产生单位矩阵。rand：产生0～1间均匀分布的随机矩阵。randn：产生均值为0，方差为1的标准正态分布随机矩阵。 2. 用于专门学科的特殊矩阵 这类矩阵的元素值具有一定的规律，在有关学科中 是很有用的。如魔方矩阵、希尔伯特矩阵等。 例1 特殊矩阵示例。 zeros(3） zeros(3,2） ones(4,3） ones(3） rand(3） magic(5) hilb(5) 3.3 矩阵运算  3.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 A=pascal(3) B=magic(3) X=A*B Y=B*A (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。 对于矩阵来说，左除和右除表示两种不同的除数矩阵和被除数矩阵的关系。对于矩阵运算，一般A\B≠B/A。 例3 (4) 矩阵的乘方(幂) 一个矩阵的乘方运算可以表示成A^x，要求A为方阵，x为标量。例4 2．点运算 在MATLAB中，有一种特殊的运算，因为其运算符是在有关算术运算符前面加点，所以叫点运算。 点运算符有.*、./、.\和.^。两矩阵进行点运算是指它们的对应元素进行相关运算，要求两矩阵的维参数相同。 例5 3.3.2 关系运算 MATLAB提供了6种关系运算符：(小于)、=(小于或等于)、(大于)、=(大于或等于)、= =(等于)、～=(不等于)。 它们的含义不难理解，但要注意其书写方法与数学中的不等式符号不尽相同。 关系运算符的运算法则为： (1) 当两个比较量是标量时，直接比较两数的大小。若关系成立，关系表达式结果为1，否则为0。 (2) 当参与比较的量是两个维数相同的矩阵时，比较是对两矩阵相同位置的元素按标量关系运算规则逐个进行，并给出元素比较结果。最终的关系运算的结果是一个维数与原矩阵相同的矩阵，它的元素由0或1组成。 (3) 当参与比较的一个是标量，而另一个是矩阵时，则把标量与矩阵的每一个元素按标量关系运算规则逐个比较，并给出元素比较结果。最终的关系运算的结果是一个维数与原矩阵相同的矩阵，它的元素由0或1组成。 例6 产生5阶随机方阵A，其元素为[10,90]区间的随机整数，然后判断A的元素是否能被3整除。 (1) 生成5阶随机方阵A。A=fix((90-10+1)*rand(5)+10) (2) 判断A的元素是否可以被3整除。 rem(A,3)==0其中，rem(A,3)是矩阵A的每个元素除以3的余数矩阵。 3.3.3 逻辑运算 MATLAB提供了3种逻辑运算符：(与)、|(或)和