《第2章矩阵运算及其应用精选文档.ppt》由会员分享，可在线阅读，更多相关《第2章矩阵运算及其应用精选文档.ppt（68页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、第2章矩阵运算及其应用本讲稿第一页，共六十八页2.1 矩阵的加减乘法矩阵的加减乘法2.1.1 矩阵的加法矩阵的加法定义定义2.1 设有两个同型的 矩阵，矩阵矩阵A与矩阵与矩阵B的和的和记作 ，规定为：本讲稿第二页，共六十八页若 ，把 记作 ，称为A的负矩阵。显然有：由此可定义矩阵的减法矩阵的减法为：本讲稿第三页，共六十八页2.1.2 矩阵的数乘矩阵的数乘定义定义2.2 数 与矩阵 的乘积，简称数乘数乘，记作 或 ，规定为本讲稿第四页，共六十八页矩阵的加法和数乘统称为矩阵的线性运算，矩阵的线性运算满足下列运算规律（A、B、C是同型矩阵，、是数）（1）加法交换律（2）加法结合律 （3）（4）本讲稿
2、第五页，共六十八页（5）（6）（7）（8）数乘分配律 本讲稿第六页，共六十八页2.1.3 矩阵的乘法矩阵的乘法定义定义2.3 设A是矩阵，B是矩阵，那么矩阵矩阵A 和矩阵和矩阵B的乘积的乘积是一个矩阵C，其中记作 C=AB本讲稿第七页，共六十八页由定义知，只有当第一个矩阵的列数和第二个矩阵的行数相等，即它们的内阶数相等时，两个矩阵才能相乘。乘积矩阵的第 元素等于前一个矩阵的第 行各元素与后一个矩阵的第 列相应元素乘积之和，即：本讲稿第八页，共六十八页定义定义2.4 对于变量 ，若它们都能由 变量线性表示，即有：（2-1）则称此关系式为变量 到变量 的线性变换线性变换。本讲稿第九页，共六十八页可
3、以写成输出向量Y等于系数矩阵A左乘输入向量X：本讲稿第十页，共六十八页例例2.4 式（2-2）给出变量 到变量 的线性变换；式（2-3）给出变量 到变量 的线性变换。请写出变量 到变量 的线性变换。（2-2）（2-3）本讲稿第十一页，共六十八页解：方法一方法一，代换法。将式（2-3）代入式（2-2），得：（2-4）方法二方法二，矩阵运算法。根据矩阵乘法的定义，可以把式（2-2）和式（2-3）分别写为式（2-5）和式（2-6）的矩阵等式：本讲稿第十二页，共六十八页 （2-5）（2-6）把式（2-6）代入式（2-5）中，得：本讲稿第十三页，共六十八页 （2-7）式（2-7）和式（2-4）等价。通过
4、这个例子，可以看出矩阵乘法在线性变换中的运用。本讲稿第十四页，共六十八页有了矩阵乘法的定义后，可以把一般的线性方程组（1-3）写为矩阵形式：（2-8）若用A表示系数矩阵，X表示未知量构成的向量，b表示常数项所构成的向量，则式（2-8）可以化简为：AX=b本讲稿第十五页，共六十八页例例2.5 已知 ，求 AB，BA解：根据矩阵乘法定义，有：本讲稿第十六页，共六十八页由于矩阵有2列，矩阵有3行，所以B不能左乘A。由矩阵乘法定义和前面的例题可以看出：（1）矩阵乘法不满足交换律，即在一般情况 下 （2）不能由 ，推出 或本讲稿第十七页，共六十八页（3）不能由 ，推出 在一般情况下有：矩阵乘法满足下列运
5、算规律：（1）（2）本讲稿第十八页，共六十八页（3），为数（4）（5），其中 为正 整数，必须为方阵。本讲稿第十九页，共六十八页2.1.4 矩阵的转置矩阵的转置定义定义2.5 设是 一个 矩阵，将矩阵 中的行换成同序数的列得到的一个 矩阵，称为矩阵 的转置矩阵转置矩阵，记作 ，或 。例如，则本讲稿第二十页，共六十八页矩阵转置满足以下运算规律（1）（2）（3）（4）本讲稿第二十一页，共六十八页在此只证明（4）证明：设 ,记 ，据矩阵乘法定义及矩阵转置定义知：而 的第 行就是 的第 列，为：，的第 列就是 的第 行，为：，因而有本讲稿第二十二页，共六十八页即得 ，亦即 。定义定义2.6 如果n阶方
6、阵 满足 ，则称为对称矩阵对称矩阵。如果n阶方阵 满足 ，则称为反对称矩阵反对称矩阵。显然反对称阵的主对角线上元素都是零。本讲稿第二十三页，共六十八页2.2 矩阵的逆矩阵的逆2.2.1 逆矩阵的定义逆矩阵的定义定义定义2.7 设 为n阶方阵，若存在n阶方阵 ，使得 ，其中 为n阶单位矩阵，则称 为可逆矩阵可逆矩阵或 是可逆的可逆的，并称 为 的逆矩阵。如果 的逆矩阵为 ，记 ，显然，则 的逆矩阵为 ，记 ，我们也称矩阵 和矩阵 互逆互逆。本讲稿第二十四页，共六十八页例例2.7 设 ,分析矩阵 和矩阵 、矩阵 和矩阵 的关系。解：本讲稿第二十五页，共六十八页所以，矩阵 和矩阵 互为逆矩阵。矩阵
7、和矩阵 也互为逆矩阵。本讲稿第二十六页，共六十八页2.2.2 逆矩阵的性质逆矩阵的性质性质性质1 如果矩阵 可逆，则 的逆矩阵唯一性质性质2 若 和 为同阶方阵，且满足 则 ，即矩阵 和矩阵 互逆。性质性质3 若 可逆，则 也可逆，且 性质性质4 若 可逆，数 ，则 可逆，且 性质性质5 若 、均为 阶可逆方阵，则 也可 逆，且 本讲稿第二十七页，共六十八页性质性质6 若 可逆，则 也可逆，且例例2.8 设方阵 满足 ，试证 可逆，并求 。解：根据已知条件，可以得到：则有：所以矩阵 可逆，且 。本讲稿第二十八页，共六十八页2.3 矩阵的分块矩阵的分块在矩阵运算中，特别是针对高阶矩阵，常常采用矩
8、阵分块的方法将其简化为较低阶的矩阵运算。用若干条纵线和横线将矩阵分为若干个小矩阵，每一个小矩阵称为的子块子块，以子块为元素的矩阵，称为分块矩阵分块矩阵。本讲稿第二十九页，共六十八页比如将43矩阵 分为 ,它们可分别表示为：本讲稿第三十页，共六十八页分块矩阵的运算与普通矩阵类似，1加法运算加法运算设 ，都是 矩阵，且将 ，按完全相同的方法分块：本讲稿第三十一页，共六十八页2数乘运算数乘运算设 ,有：3乘法运算乘法运算设 为 矩阵，为 矩阵，将它们分别分块成本讲稿第三十二页，共六十八页其中 的列数分别等于 的行数 ，即 可以左乘 。则有：其中 本讲稿第三十三页，共六十八页4转置运算转置运算 设 有
9、：注意分块矩阵的转置，不仅要把每个子块内的元素位置转置，而且要要把子块本身的位置转置。本讲稿第三十四页，共六十八页5分块对角矩阵分块对角矩阵如果将方阵 分块后，有以下形式：其中主对角线上的子块 均是方阵，而其余子块全是零矩阵，则称 为分块分块对角矩阵对角矩阵，记为 。本讲稿第三十五页，共六十八页设有两个同型且分块方法相同的对角矩阵则有 本讲稿第三十六页，共六十八页对于上面的分块矩阵 ，若对角线上的所有子块 都可逆，则有：例例2.9 利用分块矩阵的概念，把下列线性方程组写成向量等式。本讲稿第三十七页，共六十八页解：线性方程组的矩阵表示为：把系数矩阵按列分成4块：与常数矩阵 分别用向量 和向量 来
10、表示，则有：本讲稿第三十八页，共六十八页进而得到向量等式：本讲稿第三十九页，共六十八页2.4 初等矩阵初等矩阵定义定义2.8 单位矩阵 经过一次初等变换所得到的矩阵称为初等矩阵初等矩阵或初等方阵初等方阵。前面介绍了三种初等变换，每一种初等变换，都有一个相对应的初等矩阵（1）交换单位矩阵 的 ,两行（或 ,两列），得到的初等矩阵记为 ，即：本讲稿第四十页，共六十八页 （2-12）本讲稿第四十一页，共六十八页（2）用一个非零数 乘单位矩阵 的第 行（或第 列），得到的初等矩阵记为 ，即：（2-13）本讲稿第四十二页，共六十八页（3）将单位矩阵 第 行的 倍加到第 行上（或将单位矩阵第 列的 倍加到
11、第 列上）得到的初等矩阵记为 ，即：（2-14）本讲稿第四十三页，共六十八页例例2.10 设求：E1*A，E2*A，E3*A。本讲稿第四十四页，共六十八页解：本讲稿第四十五页，共六十八页定理定理2.1 设设 是一个是一个 矩阵，对矩阵，对 施行一施行一次初等行变换，其结果等于在次初等行变换，其结果等于在 的左边乘以的左边乘以相应的相应的 阶初等矩阵；对阶初等矩阵；对 施行一次初等列变施行一次初等列变换，其结果等于在换，其结果等于在 的右边乘以相应的的右边乘以相应的 阶初阶初等矩阵。等矩阵。本讲稿第四十六页，共六十八页定理定理2.2 设设 为为 阶方阵，那么下面各命题阶方阵，那么下面各命题等价：
12、等价：（1）是可逆矩阵；是可逆矩阵；（2）线性方程组）线性方程组 只有零解；只有零解；（3）可以经过有限次初等行变换化为单位可以经过有限次初等行变换化为单位 矩阵矩阵 ；（4）可以表示为有限个初等矩阵的乘积。可以表示为有限个初等矩阵的乘积。本讲稿第四十七页，共六十八页例例2.11 设判断 、是否可逆，如果可逆，请求之。解：本讲稿第四十八页，共六十八页则矩阵 可逆，且其逆为：本讲稿第四十九页，共六十八页显然矩阵 通过初等行变换不能化为单位矩阵，则矩阵 不可逆。是降秩的。它通过初等行变换，可以化出一个零行，则其秩为2。故当A不可逆时，(2-15)式应改为：其中是 秩为r的nn方阵，rM=0.1,0
13、.3,0.15;0.3,0.4,0.25;0.1,0.2,0.15;P=4000,4500,4500,4000;2000,2800,2400,2200;5800,6200,6000,6000;Q=M*P Q=1870 2220 2070 1960 3450 4020 3810 3580 1670 1940 1830 1740本讲稿第五十三页，共六十八页为了进一步计算矩阵Q的每一行和每一列的和，可以继续键入：Q*ones(4,1)ans=8120 14860 7180ones(1,3)*Qans=6990 8180 7710 7280并可以继续算出全年的总成本：ans*ones(4,1)ans=
14、30160 本讲稿第五十四页，共六十八页根据以上计算结果，可以完成每季度总成本分类表，如表2.8所示。表2.8 每季度总成本分类表成本（元）夏 秋 冬 春 全年 原材料 18702220 2070 19608120劳动 34504020 3810 3580 14860企业管理费 16701940 1830 17407180总成本（元）69908180 7710 7280 30160本讲稿第五十五页，共六十八页2.5.2 特殊矩阵的生成特殊矩阵的生成例例2.13 在Matlab环境下生成矩阵X：矩阵X有相同的10行，每一行都是公差为1的等差数列。解：令则 ，就实现了矩阵赋值。本讲稿第五十六页，共
15、六十八页键入MATLAB语句：v1=-10:10;v2=ones(1,10)X=v2*v1例例2.14 在Matlab环境下生成范德蒙矩阵。解：这里用了Matlab的符号运算功能。键入：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%两个方阵的元素群求幂 本讲稿第五十七页，共六十八页程序的运行结果为：Matlab内置的范德蒙矩阵生成函数vander.m是不能用符号表示的，只能产生数值矩阵。本讲稿第五十八页，共六十八页2.5.3 逆矩阵的求解逆矩阵的求解例例
16、2.15 设 试求其逆阵解：当矩阵的阶数较高时，利用Matlab辅助计算就尤显重要。用Matlab来求矩阵的逆，其方法很多。首先在Matlab环境下键入：A=3,0,3,-6;5,-1,1,-5;-3,1,4,-9;1,-3,4,-4;本讲稿第五十九页，共六十八页方法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
17、 -0.2424 -0.1515 0.0303本讲稿第六十页，共六十八页例例2.16 求矩阵 的逆。解：矩阵求逆命令inv也可以用符号变量。在Matlab环境下，键入：syms a 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)本讲稿第六十一页，共六十八页2.5.4 图及其矩阵表述图及其矩阵表述例如，图2.1为1，2，3，4四个城市之间的空运航线的有向图。图2.1 航线图它可以用下列航路矩阵表示：本讲稿第六十二页，共六十八页2.5.5 网络的矩阵分割和连接网络的矩阵分割和连接在电路设
18、计中，经常要把复杂的电路分割为局部电路，每一个电路都用一个网络“黑盒子”来表示。“黑盒子”的输入为u1，i1，输出为u2，i2，其输入输出关系用矩阵A来表示（如图2.2所示）：图2.2 单个子网络模型 本讲稿第六十三页，共六十八页以图2.3为例，把两个电阻组成的分压电路分成两个串接的子网络。第一个子网络包含电阻R1，第二个子网络包含电阻R2，列出第一个子网络的电路方程为：图2.3 两个子网络串联模型本讲稿第六十四页，共六十八页写成矩阵方程为：同样可列出第二个子网络的电路方程，写成矩阵方程为：本讲稿第六十五页，共六十八页从上分别得到两个子网络的传输矩阵整个电路的传输矩阵为两者的乘积本讲稿第六十六页，共六十八页2.5.6 弹性梁的柔度矩阵弹性梁的柔度矩阵设简支梁如图2.4所示，在梁的三个位置分别施加力f1，f2和f3后，在该处产生的综合变形（通常称为挠度）为图示的y1，y2和y3。图2.4 简支梁在三个点的力和变形（挠度）本讲稿第六十七页，共六十八页根据虎克定律，在材料未失去弹性的范围内，三个力与它引起的三个变形都呈线性关系，可以写成矩阵形式：用向量和矩阵符号表示为：yDf D中的各元素为挠度元素，这些元素的值越大表明这个梁愈柔软，所以矩阵D被称作柔度矩阵。柔度矩阵D的逆就是刚度矩阵K。本讲稿第六十八页，共六十八页