第2章 矩阵和数组及其运算

第二讲数组及矩阵运算数值数组（矩阵）及其运算是MATLAB的核心。

它可以使计算程序简单、易读，使程序更接近于教科书上的数学计算公式，这也是理工科学生最需要的。

另外，它提高了程序的向量化程度，提高计算效率，节省计算机的开销。

注意：矩阵与二维数组在形式上相同，但不同于数组，矩阵具有特定的含义和运算方式。

>>x1＝[0,0.1,0.2,0.3,0.4,0.5]%简单的一维数组>>y1=x1.*exp(-x1)%数组的乘积运算>>x2=rand(1,4) %产生4个随机数，构成行向量>>x＝[1 2 3; 4 5 6; 7 8 9]%二维数组，也可认为是3行2列的矩阵x＝1 2 34 5 67 8 9一、数值数组的创建1. 一维数组的创建逐个元素输入法这是最简单又最通用的方法，用中括号[ ]表示，输入元素以空格或逗号为间隔。

函数生成法利用MATLAB提供的特殊形式的数组/矩阵函数来生成。

特点：数组元素的大小按递增或递减的次序排列，元素之间的“差”是“等步长”的。

作用：用作函数的自变量，for循环中的循环自变量等。

1）“冒号：”生成法x=a ：inc ：b其中：inc 是采样点间隔，即步长，默认值为1。

递增/递减型一维数组的创建2）线性或对数定点法¾x=linspace(a,b,n) 线性等间隔为1×n行数组其中：a，b为左右端点，n为总采样点数，即一维数组的长度，它等价于x=a:(b-a)/(n-1):b¾x=logspace(a,b,n) 对数等间隔为1×n行数组例1：产生一个从0到10的数组，间隔为2 >>a=0:2:10a=0 2 4 6 8 10>>b=linspace(0,10,6)b=0 2 4 6 8 10例2：以对数等分隔>>ak=logspace(0,1,5)ak=0 10 10^2 10^3 10^42. 二维数组的创建小规模数组的直接输入法•整个输入组必须以方括号[ ]为其首尾；•行与行之间必须用分号；或[Enter]键隔离；•数组元素间必须用逗号或空格分隔；•若元素之间有复数，则元素间不能用空格；中规模数组的数组编辑器创建法数组创建之后，将其保存为Matrix.mat文件。

§1 矩阵及其运算一、矩阵的基本概念（必考）矩阵，是由m*n个数组成的一个m行n列的矩形表格，通常用大写字母表示，组成矩阵的每一个数，均称为矩阵的元素，通常用小写字母其元素表示，其中下标都是正整数，他们表示该元素在矩阵中的位置.比如,或表示一个m*n 矩阵，下标ij 表示元素位于该矩阵的第行、第列.元素全为零的矩阵称为零矩阵. 特别地，一个m*1矩阵，也称为一个 m维列向量；而一个 1*n矩阵B=(b1,b2,…，bn)，也称为一个 n维行向量.当一个矩阵的行数m与烈数n 相等时，该矩阵称为一个 n阶方阵.若一个n阶方阵的主对角线上的元素都是，而其余元素都是零，则称为单位矩阵，记为，即： .单位矩阵与实数中的‘1’的运算相近.如一个阶方阵的主对角线上（下）方的元素都是零，则称为下（上）三角矩阵是一个阶下三角矩阵.例题：1.A既是上三角矩阵，又是下三角矩阵，则A必是对角矩阵2.两矩阵既可相加又可相乘的充要条件是两矩阵为同阶方阵.3.A=(l≠n)，则A的主对角线上个元素的和为 (设矩阵为2行3列的矩阵，找规律)二、矩阵的运算1、矩阵的加法：如果是两个同型矩阵（即它们具有相同的行数和列数，比如说），则定义它们的和仍为与它们同型的矩阵（即），的元素为和对应元素的和，即：.给定矩阵，我们定义其负矩阵为： .这样我们可以定义同型矩阵的减法为： .由于矩阵的加法运算归结为其元素的加法运算，容易验证，矩阵的加法满足下列运算律：(1)交换律：； (2)结合律：；(3)存在零元：;(4)存在负元：.2 、数与矩阵的乘法的运算律：（1)；（2)；（3)；（4) .3 、矩阵的乘法（必考）设为距阵，为距阵，则矩阵可以左乘矩阵（注意：距阵的列数等与矩阵的行数），所得的积为一个距阵，即,其中,并且(即左行乘右列)矩阵的乘法满足下列运算律（假定下面的运算均有意义）：(1）结合律：； (2）左分配律：；(3）右分配律：；(4）数与矩阵乘法的结合律：；(5）单位矩阵的存在性：.若为阶方阵，则对任意正整数，我们定义：，并规定：由于矩阵乘法满足结合律，我们有：， .注意：矩阵的乘法与通常数的乘法有很大区别，特别应该注意的是：（必考重要）（1)矩阵乘法不满足交换律：一般来讲即便有意义，也未必有意义；倘使都有意义，二者也未必相等.正是由于这个原因，一般来讲,在实数中的某些运算不再适应，如,，反过来，这些公式成立的条件又恰是A、B 可逆.例：A,B,C 是同阶矩阵，A ≠0,若AB=BC,必有B=C,则A满足可逆（2)两个非零矩阵的乘积可能是零矩阵，即未必能推出或者. 同理，A ≠0，B ≠0，而AB却肯能等于0.例题：（选择题5、6）（3)矩阵的乘法不满足消去律：如果并且，未必有 .4 、矩阵的转置：定义：设为矩阵，我们定义的转置为一个矩阵，并用表示的转置，即：.矩阵的转置运算满足下列运算律：（1)；（2)；（3)；（4) (重要).5、对称矩阵：n 阶方阵若满足条件：，则称为对称矩阵；若满足条件：，则称为反对称矩阵.若设，则为对称矩阵，当且仅当对任意的成立；为反对称矩阵，当且仅当对任意的成立.从而反对称矩阵对角线上的元素必为零.对称矩阵具有如下性质：（1)对于任意矩阵，为阶对称矩阵；而为阶对称矩阵；（2)两个同阶（反）对称矩阵的和，仍为（反）对称矩阵；（3)如果两个同阶（反）对称矩阵可交换，即，则它们的乘积必为对称矩阵，即.运算性质：1） （2） （3）（4） （5）三、逆矩阵1.定义 对于n 阶矩阵A ，如果存在n 阶矩阵B ，使得E BA AB ==．则A 称为可逆矩阵或非奇异矩阵．B 称为A 的逆矩阵，．由定义可得，A 与B 一定是同阶的，而且A 如果可逆，则A 的逆矩阵是唯一的．这是因为（反证法），如果1B 、2B 都是A 的逆矩阵，则有E A B AB ==11，E A B AB ==22，那么22212111)()(B EB B A B AB B E B B =====所以逆矩阵是唯一的．我们把矩阵A 的逆矩阵记作1-A ．逆矩阵有下列性质： （1）如果A 可逆，则1-A 也可逆，且A A =--11)(．由可逆的定义，显然有A 与1-A 是互逆的． （2）如果A 、B 是两个同阶可逆矩阵，则)(AB 也可逆，且111)(---=A B AB ．（必考重点） 这是因为 E A A AEA ABB A A B AB =⋅===------111111)())((E B B EB B B A A B AB A B ====------111111)())((，所以111)(---=A B AB ．（必考重点）这个结论也可以推广到有限个可逆矩阵想乘的情形. （3）可逆矩阵A 的转置矩阵T A 也是可逆矩阵，且T T A A )()(11--=．这是因为E E A A A A T T TT===--)()(11，E E AA A A T T T T ===--)()(11所以 T TA A )()(11--=．（4）如果A 是可逆矩阵，则有11--=A A ．这是因为E AA=-1，两边取行列式有 11=⋅-A A ，所以111--==A AA ． 矩阵可逆的条件（1）n 阶方阵A 可逆的充分必要条件是| A | ≠ 0（也即r （A ）= n ）；（2）n 阶方阵A 可逆的充分必要条件是A 可以通过初等变换（特别是只通过初等行（列）变换）化为n 阶单位矩阵；（3）n 阶方阵A 可逆的充分必要条件是A 可以写成一些初等矩阵的乘积；（4）n 阶方阵A 可逆的充分必要条件是A 的n 个特征值不为零；（5）对于n 阶方阵A ，若存在n 阶方阵B 使得AB = E （或BA = E ），则A 可逆，且A -1= B. 逆矩阵的有关结论及运算必考 ——求法方法1 定义法：设A 是数域P 上的一个n 阶方阵，如果存在P 上的n 阶方阵B ，使得AB = BA= E ，则称A 是可逆的，又称B 为A 的逆矩阵.当矩阵A 可逆时，逆矩阵由A 惟一确定，记为A -1.例1：设A 为n 阶矩阵，且满足22A - 3A + 5E = 0，求A -1.【解】22 2 -12A - 3A + 5E = 02A - 3A = - 5E23-A - A =E 552323A (- A - E) = - A - E = E555523A A = - A - E55∴∴∴∴可逆且方法 2 伴随矩阵法：A -1= 1|A|A*.定理n 阶矩阵A = a ij 为可逆的充分必要条件是A 非奇异.且11211122221121n n nnnn A A A A A A A A A A A -⎛⎫ ⎪ ⎪=⎪ ⎪⎝⎭其中A ij 是|A|中元素a ij 的代数余子式.矩阵112111222212n n nnnn A A A A A A A A A ⎛⎫ ⎪ ⎪⎪ ⎪⎝⎭称为矩阵A 的伴随矩阵，记作A*，于是有A -1=1|A|A*. 注 ①对于阶数较低（一般不超过3阶）或元素的代数余子式易于计算的矩阵可用此法求其逆矩阵.注意A* = （A ji ）n ×n 元素的位置及符号.特别对于2阶方阵11122122a a A a a ⎛⎫= ⎪⎝⎭，其伴随矩阵22122111*a a A a a -⎛⎫=⎪-⎝⎭,即伴随矩阵具有“主对角元素互换，次对角元素变号”的规律.②对于分块矩阵A B C D ⎛⎫⎪⎝⎭不能按上述规律求伴随矩阵.例2：已知101A=210325⎛⎫ ⎪ ⎪ ⎪--⎝⎭，求A -1.【解】 ∵| A | = 2 ≠ 0 ∴A 可逆.由已知得111213212223313233A = - 5, A = 10, A = 7A = 2, A = - 2, A = - 2A = - 1, A = 2, A = 1 ， A -1= 1|A| A* = 5115212211022511272171122⎛⎫-- ⎪--⎛⎫ ⎪⎪-=- ⎪ ⎪ ⎪ ⎪-⎝⎭- ⎪⎝⎭方法3 初等变换法：注 ①对于阶数较高（n ≥3）的矩阵，采用初等行变换法求逆矩阵一般比用伴随矩阵法简便.在用上述方法求逆矩阵时，只允许施行初等行变换.②也可以利用1E A E A -⎛⎫⎛⎫−−−−→⎪ ⎪⎝⎭⎝⎭初等列变换求得A 的逆矩阵. ③当矩阵A 可逆时，可利用求解求得A -1B 和CA -1.这一方法的优点是不需求出A 的逆矩阵和进行矩阵乘法，仅通过初等变换即求出了A -1B 或CA -1.例3：:用初等行变换求矩阵231A 013125⎛⎫⎪= ⎪ ⎪⎝⎭的逆矩阵.【解】()231100125001125001A E 01301001301001301012500123110000611212500112500101301001301001910211100166311341006631310122111001663⎛⎫⎛⎫⎛⎫ ⎪ ⎪ ⎪=→→ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪--⎝⎭⎝⎭⎝⎭⎛⎫ ⎪⎛⎫⎪ ⎪→→ ⎪ ⎪ ⎪ ⎪---⎝⎭-- ⎪⎝⎭⎛--→---⎝⎫⎪⎪⎪⎪ ⎪⎪ ⎪⎭1113410066313A 010********1663-⎛⎫--⎪ ⎪ ⎪=- ⎪ ⎪ ⎪-- ⎪⎝⎭故 方法4 用分块矩阵求逆矩阵：设A 、B 分别为P 、Q 阶可逆矩阵，则：1111111111111111A A 000B 0C O A A A CB A O A O BD B O B B DA B B O A O B B O AO ----------------⎛⎫⎛⎫⎛⎫-⎛⎫⎛⎫⎛⎫=== ⎪ ⎪ ⎪ ⎪ ⎪ ⎪-⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭⎛⎫⎛⎫= ⎪ ⎪⎝⎭⎝⎭例4：已知0052002112001100A ⎛⎫⎪ ⎪=⎪-⎪⎝⎭，求A -1.【解】 将A 分块如下：12005200211200110O A A A O ⎛⎫ ⎪ ⎪⎛⎫⎪== ⎪⎪⎝⎭- ⎪ ⎪⎝⎭其中 125212,2111A A -⎛⎫⎛⎫== ⎪ ⎪⎝⎭⎝⎭可求得 1*1*1122121212111,2511||||3A A A A A A ---⎛⎫⎛⎫==== ⎪ ⎪--⎝⎭⎝⎭ 从而11211120033110331200250O A A A O ---⎛⎫ ⎪ ⎪ ⎪-⎛⎫ ⎪== ⎪⎪⎝⎭ ⎪ ⎪- ⎪-⎝⎭方法5 恒等变形法求逆矩阵：有些计算命题表面上与求逆矩阵无关,但实质上只有求出矩 阵的逆矩阵才能算出来,而求逆矩阵须对所给的矩阵等式恒等变 形,且常变形为两矩阵的乘积等于单位矩阵的等式.例8 已知，且，试求．解 由题设条件得3.伴随矩阵 如果n 阶矩阵A 的行列式0≠A ,则称A 是非奇异的(或非退化的).否则,称A 是奇异的(或退化的).（n 阶矩阵A 可逆的充要条件是:|A|≠0）设n n ij a A ⨯=)(,ij A 是A 中元素)21(n j i a ij ，，，， =的代数余子式.矩阵 ⎪⎪⎪⎪⎪⎭⎫⎝⎛=nn n n n n A A A A A A A A A A 212221212111*（顺序变化，重点）称为A 的伴随矩阵. 矩阵n n ij a A ⨯=)(为可逆矩阵的充分必要条件是A 为非奇异矩阵,并且当A 可逆时,有*11A AA =-，伴随矩阵 例1. 已知矩阵⎪⎪⎪⎭⎫ ⎝⎛---=313132121A 判断A 是否可逆,如果可逆,求1-A .解: 因为01313132121≠=---=A ,所以A 可逆.又.13221)1(11211)1(;11312)1(71321)1(;63311)1(53112)1(;11332)1(93312)1(;83113)1(333323321331322322221221311321121111=---==-==---=-=--=-=--=-=---==--==--==---=+++++++++A A A A A A A A A所以 ⎪⎪⎪⎭⎫⎝⎛---==-1711691581*1A A A 四、分块矩阵一、分块矩阵的概念对于行数和列数较高的矩阵, 为了简化运算，经常采用分块法，使大矩阵的运算化成若干小矩阵间的运算，同时也使原矩阵的结构显得简单而清晰. 具体做法是：将大矩阵用若干条纵线和横线分成多个小矩阵. 每个小矩阵称为A 的子块, 以子块为元素的形式上的矩阵称为分块矩阵.矩阵的分块有多种方式，可根据具体需要而定注：一个矩阵也可看作以n m ⨯个元素为1阶子块的分块矩阵. 二、分块矩阵的运算分块矩阵的运算与普通矩阵的运算规则相似. 分块时要注意，运算的两矩阵按块能运算，并且参与运算的子块也能运算，即，内外都能运算.1. 设矩阵A 与B 的行数相同、列数相同，采用相同的分块法, 若,,11111111⎪⎪⎪⎭⎫⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛=st s t st s t B B B B B A A A A A其中ij A 与ij B 的行数相同、列数相同, 则.11111111⎪⎪⎪⎭⎫ ⎝⎛++++=+st st s s t t B A B A B A B A B A2.设,1111⎪⎪⎪⎭⎫ ⎝⎛=st s t A A A A Ak 为数, 则.1111⎪⎪⎪⎭⎫ ⎝⎛=st s t kA kA kA kA kA 3．设A 为l m ⨯矩阵, B 为n l ⨯矩阵, 分块成,,11111111⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛=tr t r st s t B B B B B A A A A A其中pt p p A A A ,,,21 的列数分别等于tq q q B B B ,,,21 的行数, 则,1111⎪⎪⎪⎭⎫ ⎝⎛=sr s r C C C C AB 其中).,,2,1;,,2,1(1r q s p B A C t k kqpk pq ===∑=4. 分块矩阵的转置设,1111⎪⎪⎪⎭⎫ ⎝⎛=st s t A A A A A则.1111⎪⎪⎪⎪⎭⎫ ⎝⎛=T st T tT s T TA A A A A 5. 设A 为n 阶矩阵, 若A 的分块矩阵只有在对角线上有非零子块, 其余子块都为零矩阵, 且在对角线上的子块都是方阵, 即⎪⎪⎪⎪⎪⎭⎫⎝⎛=s A O A O A A21, 其中),,2,1(s i A i =都是方阵, 则称A 为分块对角矩阵.分块对角矩阵具有以下性质：(1) 若 ),,2,1(0||s i A i =≠，则0||≠A ，且|;|||||||21s A A A A =(2) .112111⎪⎪⎪⎪⎪⎭⎫⎝⎛=----s A O A O A A(3) 同结构的对角分块矩阵的和、差、积、商仍是对角分块矩阵. 且运算表现为对应子块的运算。

12m m mna a a 矩阵。

为了表示它是一个整体，总是加一个括号将它界起来，并通常用大写字母表示它。

记做12m m mn a a a ⎥⎦12m m mn a a a a ⎛⎪⎭。

切记不允许使用111212122212n n m m mna a a a a a a a a =A 。

矩阵的横向称行，纵向称列。

矩阵中的每个数称为元素，所有元素都是实数的矩阵称为实矩阵，所有元素都是复数的矩阵称为复矩阵。

本课中的矩阵除特殊说明外，都指12n n nn a a a ⎥⎦不是方阵没有主对角线。

在方阵中，00nn a ⎥⎦11212212000n n nn a a a a a a ⎤⎥⎥⎥⎥⎦（主对角线以上均为零）1122000000nn a aa ⎡⎤⎢⎥⎢⎥⎥⎥⎦（既}nn a .对角元素为1的对角矩阵，记作Ｅ 或001⎡⎢⎥⎦()11a ，此时矩阵退化为一个数矩阵的引进为许多实际的问题研究提供方便。

a x +)1(+⨯n 矩阵：12m m mnm a b a a a b ⎥⎦任何一个方程组都可以用这样一个矩阵来描述；反之，一个矩阵也完全刻划了一个方122m m m mn mn b a b a b ⎥+++⎦⎥⎦⎤⎢⎣⎡-=4012B ，计算 B A +。

122m m m mn mn b a b a b ⎥---⎦与矩阵n m ij a A ⨯=}{的乘积（称之为数乘），12m m mn a a a λλ⎥⎦以上运算称为矩阵的线性运算，它满足下列运算法则：n b ⎪⎭上述几个例子显示，当有意义时，不一定有意义（例6），即便有相同的阶数，也不一定相等（例A = O 或Ba x +12m m mn a a a ⎥⎦为系数矩阵； m b ⎥⎦，称b 为常数项矩阵；12n x x x ⎡⎢⎢=⎥⎦X = b 。

四、矩阵的转置 5 (转置矩阵12m m mn a a a ⎥⎦12nnmn a a a ⎢⎥⎣⎦矩阵，称它为A 的转置矩阵，记作TA 。