第2章矩阵运算及其应用优秀课件.ppt

《第2章矩阵运算及其应用优秀课件.ppt》由会员分享，可在线阅读，更多相关《第2章矩阵运算及其应用优秀课件.ppt（68页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第2章矩阵运算及其应用第1页，本讲稿共68页2.1 矩阵的加减乘法矩阵的加减乘法2.1.1 矩阵的加法矩阵的加法定义定义2.1 设有两个同型的 矩阵，矩阵矩阵A与矩阵与矩阵B的和的和记作 ，规定为：第2页，本讲稿共68页若 ，把 记作 ，称为A的负矩阵。显然有：由此可定义矩阵的减法矩阵的减法为：第3页，本讲稿共68页2.1.2 矩阵的数乘矩阵的数乘定义定义2.2 数 与矩阵 的乘积，简称数乘数乘，记作 或 ，规定为第4页，本讲稿共68页矩阵的加法和数乘统称为矩阵的线性运算，矩阵的线性运算满足下列运算规律（A、B、C是同型矩阵，、是数）（1）加法交换律（2）加法结合律 （3）（4）第5页，本讲稿

2、共68页（5）（6）（7）（8）数乘分配律 第6页，本讲稿共68页2.1.3 矩阵的乘法矩阵的乘法定义定义2.3 设A是矩阵，B是矩阵，那么矩阵矩阵A 和矩阵和矩阵B的乘积的乘积是一个矩阵C，其中记作 C=AB第7页，本讲稿共68页由定义知，只有当第一个矩阵的列数和第二个矩阵的行数相等，即它们的内阶数相等时，两个矩阵才能相乘。乘积矩阵的第 元素等于前一个矩阵的第 行各元素与后一个矩阵的第 列相应元素乘积之和，即：第8页，本讲稿共68页定义定义2.4 对于变量 ，若它们都能由 变量线性表示，即有：（2-1）则称此关系式为变量 到变量 的线性变换线性变换。第9页，本讲稿共68页可以写成输出向量Y等

3、于系数矩阵A左乘输入向量X：第10页，本讲稿共68页例例2.4 式（2-2）给出变量 到变量 的线性变换；式（2-3）给出变量 到变量 的线性变换。请写出变量 到变量 的线性变换。（2-2）（2-3）第11页，本讲稿共68页解：方法一方法一，代换法。将式（2-3）代入式（2-2），得：（2-4）方法二方法二，矩阵运算法。根据矩阵乘法的定义，可以把式（2-2）和式（2-3）分别写为式（2-5）和式（2-6）的矩阵等式：第12页，本讲稿共68页 （2-5）（2-6）把式（2-6）代入式（2-5）中，得：第13页，本讲稿共68页 （2-7）式（2-7）和式（2-4）等价。通过这个例子，可以看出矩阵乘

4、法在线性变换中的运用。第14页，本讲稿共68页有了矩阵乘法的定义后，可以把一般的线性方程组（1-3）写为矩阵形式：（2-8）若用A表示系数矩阵，X表示未知量构成的向量，b表示常数项所构成的向量，则式（2-8）可以化简为：AX=b第15页，本讲稿共68页例例2.5 已知 ，求 AB，BA解：根据矩阵乘法定义，有：第16页，本讲稿共68页由于矩阵有2列，矩阵有3行，所以B不能左乘A。由矩阵乘法定义和前面的例题可以看出：（1）矩阵乘法不满足交换律，即在一般情况 下 （2）不能由 ，推出 或第17页，本讲稿共68页（3）不能由 ，推出 在一般情况下有：矩阵乘法满足下列运算规律：（1）（2）第18页，本

5、讲稿共68页（3），为数（4）（5），其中 为正 整数，必须为方阵。第19页，本讲稿共68页2.1.4 矩阵的转置矩阵的转置定义定义2.5 设是 一个 矩阵，将矩阵 中的行换成同序数的列得到的一个 矩阵，称为矩阵 的转置矩阵转置矩阵，记作 ，或 。例如，则第20页，本讲稿共68页矩阵转置满足以下运算规律（1）（2）（3）（4）第21页，本讲稿共68页在此只证明（4）证明：设 ,记 ，据矩阵乘法定义及矩阵转置定义知：而 的第 行就是 的第 列，为：，的第 列就是 的第 行，为：，因而有第22页，本讲稿共68页即得 ，亦即 。定义定义2.6 如果n阶方阵 满足 ，则称为对称矩阵对称矩阵。如果n阶方

6、阵 满足 ，则称为反对称矩阵反对称矩阵。显然反对称阵的主对角线上元素都是零。第23页，本讲稿共68页2.2 矩阵的逆矩阵的逆2.2.1 逆矩阵的定义逆矩阵的定义定义定义2.7 设 为n阶方阵，若存在n阶方阵 ，使得 ，其中 为n阶单位矩阵，则称 为可逆矩阵可逆矩阵或 是可逆的可逆的，并称 为 的逆矩阵。如果 的逆矩阵为 ，记 ，显然，则 的逆矩阵为 ，记 ，我们也称矩阵 和矩阵 互逆互逆。第24页，本讲稿共68页例例2.7 设 ,分析矩阵 和矩阵 、矩阵 和矩阵 的关系。解：第25页，本讲稿共68页所以，矩阵 和矩阵 互为逆矩阵。矩阵 和矩阵 也互为逆矩阵。第26页，本讲稿共68页2.2.2

7、逆矩阵的性质逆矩阵的性质性质性质1 如果矩阵 可逆，则 的逆矩阵唯一性质性质2 若 和 为同阶方阵，且满足 则 ，即矩阵 和矩阵 互逆。性质性质3 若 可逆，则 也可逆，且 性质性质4 若 可逆，数 ，则 可逆，且 性质性质5 若 、均为 阶可逆方阵，则 也可 逆，且 第27页，本讲稿共68页性质性质6 若 可逆，则 也可逆，且例例2.8 设方阵 满足 ，试证 可逆，并求 。解：根据已知条件，可以得到：则有：所以矩阵 可逆，且 。第28页，本讲稿共68页2.3 矩阵的分块矩阵的分块在矩阵运算中，特别是针对高阶矩阵，常常采用矩阵分块的方法将其简化为较低阶的矩阵运算。用若干条纵线和横线将矩阵分为若

8、干个小矩阵，每一个小矩阵称为的子块子块，以子块为元素的矩阵，称为分块矩阵分块矩阵。第29页，本讲稿共68页比如将43矩阵 分为 ,它们可分别表示为：第30页，本讲稿共68页分块矩阵的运算与普通矩阵类似，1加法运算加法运算设 ，都是 矩阵，且将 ，按完全相同的方法分块：第31页，本讲稿共68页2数乘运算数乘运算设 ,有：3乘法运算乘法运算设 为 矩阵，为 矩阵，将它们分别分块成第32页，本讲稿共68页其中 的列数分别等于 的行数 ，即 可以左乘 。则有：其中 第33页，本讲稿共68页4转置运算转置运算 设 有：注意分块矩阵的转置，不仅要把每个子块内的元素位置转置，而且要要把子块本身的位置转置。第

9、34页，本讲稿共68页5分块对角矩阵分块对角矩阵如果将方阵 分块后，有以下形式：其中主对角线上的子块 均是方阵，而其余子块全是零矩阵，则称 为分块分块对角矩阵对角矩阵，记为 。第35页，本讲稿共68页设有两个同型且分块方法相同的对角矩阵则有 第36页，本讲稿共68页对于上面的分块矩阵 ，若对角线上的所有子块 都可逆，则有：例例2.9 利用分块矩阵的概念，把下列线性方程组写成向量等式。第37页，本讲稿共68页解：线性方程组的矩阵表示为：把系数矩阵按列分成4块：与常数矩阵 分别用向量 和向量 来表示，则有：第38页，本讲稿共68页进而得到向量等式：第39页，本讲稿共68页2.4 初等矩阵初等矩阵定

10、义定义2.8 单位矩阵 经过一次初等变换所得到的矩阵称为初等矩阵初等矩阵或初等方阵初等方阵。前面介绍了三种初等变换，每一种初等变换，都有一个相对应的初等矩阵（1）交换单位矩阵 的 ,两行（或 ,两列），得到的初等矩阵记为 ，即：第40页，本讲稿共68页 （2-12）第41页，本讲稿共68页（2）用一个非零数 乘单位矩阵 的第 行（或第 列），得到的初等矩阵记为 ，即：（2-13）第42页，本讲稿共68页（3）将单位矩阵 第 行的 倍加到第 行上（或将单位矩阵第 列的 倍加到第 列上）得到的初等矩阵记为 ，即：（2-14）第43页，本讲稿共68页例例2.10 设求：E1*A，E2*A，E3*A。

11、第44页，本讲稿共68页解：第45页，本讲稿共68页定理定理2.1 设设 是一个是一个 矩阵，对矩阵，对 施行一施行一次初等行变换，其结果等于在次初等行变换，其结果等于在 的左边乘以的左边乘以相应的相应的 阶初等矩阵；对阶初等矩阵；对 施行一次初等列变施行一次初等列变换，其结果等于在换，其结果等于在 的右边乘以相应的的右边乘以相应的 阶初阶初等矩阵。等矩阵。第46页，本讲稿共68页定理定理2.2 设设 为为 阶方阵，那么下面各命题阶方阵，那么下面各命题等价：等价：（1）是可逆矩阵；是可逆矩阵；（2）线性方程组）线性方程组 只有零解；只有零解；（3）可以经过有限次初等行变换化为单位可以经过有限次

12、初等行变换化为单位 矩阵矩阵 ；（4）可以表示为有限个初等矩阵的乘积。可以表示为有限个初等矩阵的乘积。第47页，本讲稿共68页例例2.11 设判断 、是否可逆，如果可逆，请求之。解：第48页，本讲稿共68页则矩阵 可逆，且其逆为：第49页，本讲稿共68页显然矩阵 通过初等行变换不能化为单位矩阵，则矩阵 不可逆。是降秩的。它通过初等行变换，可以化出一个零行，则其秩为2。故当A不可逆时，(2-15)式应改为：其中是 秩为r的nn方阵，rM=0.1,0.3,0.15;0.3,0.4,0.25;0.1,0.2,0.15;P=4000,4500,4500,4000;2000,2800,2400,2200

13、;5800,6200,6000,6000;Q=M*P Q=1870 2220 2070 1960 3450 4020 3810 3580 1670 1940 1830 1740第53页，本讲稿共68页为了进一步计算矩阵Q的每一行和每一列的和，可以继续键入：Q*ones(4,1)ans=8120 14860 7180ones(1,3)*Qans=6990 8180 7710 7280并可以继续算出全年的总成本：ans*ones(4,1)ans=30160 第54页，本讲稿共68页根据以上计算结果，可以完成每季度总成本分类表，如表2.8所示。表2.8 每季度总成本分类表成本（元）夏 秋 冬 春 全

14、年 原材料 18702220 2070 19608120劳动 34504020 3810 3580 14860企业管理费 16701940 1830 17407180总成本（元）69908180 7710 7280 30160第55页，本讲稿共68页2.5.2 特殊矩阵的生成特殊矩阵的生成例例2.13 在Matlab环境下生成矩阵X：矩阵X有相同的10行，每一行都是公差为1的等差数列。解：令则 ，就实现了矩阵赋值。第56页，本讲稿共68页键入MATLAB语句：v1=-10:10;v2=ones(1,10)X=v2*v1例例2.14 在Matlab环境下生成范德蒙矩阵。解：这里用了Matlab的

15、符号运算功能。键入：syms x1 x2 x3 x4 real%令x1 x2 x3 x4为实数符号变量x=x1,x2,x3,x4;y=0:3;A=x*ones(1,4)B=(ones(4,1)*yV=A.B%两个方阵的元素群求幂 第57页，本讲稿共68页程序的运行结果为：Matlab内置的范德蒙矩阵生成函数vander.m是不能用符号表示的，只能产生数值矩阵。第58页，本讲稿共68页2.5.3 逆矩阵的求解逆矩阵的求解例例2.15 设 试求其逆阵解：当矩阵的阶数较高时，利用Matlab辅助计算就尤显重要。用Matlab来求矩阵的逆，其方法很多。首先在Matlab环境下键入：A=3,0,3,-6

16、;5,-1,1,-5;-3,1,4,-9;1,-3,4,-4;第59页，本讲稿共68页方法1,A-1,方法2,inv(A),方法3,Aeye(4),方法4,U=rref(A,eye(4);U(:,5:8)运行结果都为：ans=0.2323 -0.0101 -0.1313 -0.0404 0.5354 -0.3131 -0.0707 -0.2525 0.5859 -0.4747 -0.1717 0.1010 0.2424 -0.2424 -0.1515 0.0303第60页，本讲稿共68页例例2.16 求矩阵 的逆。解：矩阵求逆命令inv也可以用符号变量。在Matlab环境下，键入：syms a

17、 b c d,A=a,b;c,d,V=inv(A)结果为：V =d/(a*d-b*c),-b/(a*d-b*c)-c/(a*d-b*c),a/(a*d-b*c)第61页，本讲稿共68页2.5.4 图及其矩阵表述图及其矩阵表述例如，图2.1为1，2，3，4四个城市之间的空运航线的有向图。图2.1 航线图它可以用下列航路矩阵表示：第62页，本讲稿共68页2.5.5 网络的矩阵分割和连接网络的矩阵分割和连接在电路设计中，经常要把复杂的电路分割为局部电路，每一个电路都用一个网络“黑盒子”来表示。“黑盒子”的输入为u1，i1，输出为u2，i2，其输入输出关系用矩阵A来表示（如图2.2所示）：图2.2 单

18、个子网络模型 第63页，本讲稿共68页以图2.3为例，把两个电阻组成的分压电路分成两个串接的子网络。第一个子网络包含电阻R1，第二个子网络包含电阻R2，列出第一个子网络的电路方程为：图2.3 两个子网络串联模型第64页，本讲稿共68页写成矩阵方程为：同样可列出第二个子网络的电路方程，写成矩阵方程为：第65页，本讲稿共68页从上分别得到两个子网络的传输矩阵整个电路的传输矩阵为两者的乘积第66页，本讲稿共68页2.5.6 弹性梁的柔度矩阵弹性梁的柔度矩阵设简支梁如图2.4所示，在梁的三个位置分别施加力f1，f2和f3后，在该处产生的综合变形（通常称为挠度）为图示的y1，y2和y3。图2.4 简支梁在三个点的力和变形（挠度）第67页，本讲稿共68页根据虎克定律，在材料未失去弹性的范围内，三个力与它引起的三个变形都呈线性关系，可以写成矩阵形式：用向量和矩阵符号表示为：yDf D中的各元素为挠度元素，这些元素的值越大表明这个梁愈柔软，所以矩阵D被称作柔度矩阵。柔度矩阵D的逆就是刚度矩阵K。第68页，本讲稿共68页