(完整版)逆矩阵的几种求法与解析(很全很经典)

整理文档很辛苦,赏杯茶钱您下走!

免费阅读已结束,点击下载阅读编辑剩下 ...

阅读已结束,您可以下载文档离线阅读编辑

资源描述

逆矩阵的几种求法与解析矩阵是线性代数的主要内容,很多实际问题用矩阵的思想去解既简单又快捷.逆矩阵又是矩阵理论的很重要的内容,逆矩阵的求法自然也就成为线性代数研究的主要内容之一.本文将给出几种求逆矩阵的方法.1.利用定义求逆矩阵定义:设A、B都是n阶方阵,如果存在n阶方阵B使得AB=BA=E,则称A为可逆矩阵,而称B为A的逆矩阵.下面举例说明这种方法的应用.例1求证:如果方阵A满足Ak=0,那么EA是可逆矩阵,且(E-A)1=E+A+A2+…+A1K证明因为E与A可以交换,所以(E-A)(E+A+A2+…+A1K)=E-AK,因AK=0,于是得(E-A)(E+A+A2+…+A1K)=E,同理可得(E+A+A2+…+A1K)(E-A)=E,因此E-A是可逆矩阵,且(E-A)1=E+A+A2+…+A1K.同理可以证明(E+A)也可逆,且(E+A)1=E-A+A2+…+(-1)1KA1K.由此可知,只要满足AK=0,就可以利用此题求出一类矩阵EA的逆矩阵.例2设A=0000300000200010,求E-A的逆矩阵.分析由于A中有许多元素为零,考虑AK是否为零矩阵,若为零矩阵,则可以采用例2的方法求E-A的逆矩阵.解容易验证1A2=0000000060000200,A3=0000000000006000,A4=0而(E-A)(E+A+A2+A3)=E,所以(E-A)1=E+A+A2+A3=1000310062106211.2.初等变换法求元素为具体数字的矩阵的逆矩阵,常用初等变换法.如果A可逆,则A可通过初等变换,化为单位矩阵I,即存在初等矩阵SPPP,,21使(1)sppp21A=I,用A1右乘上式两端,得:(2)sppp21I=A1比较(1)(2)两式,可以看到当A通过初等变换化为单位矩阵的同时,对单位矩阵I作同样的初等变换,就化为A的逆矩阵A1.用矩阵表示(AI)初等行变换为(IA1),就是求逆矩阵的初等行变换法,它是实际应用中比较简单的一种方法.需要注意的是,在作初等变换时只允许作行初等变换.同样,只用列初等变换也可以求逆矩阵.例1求矩阵A的逆矩阵.已知A=521310132.解[AI]1005210103100011320011320103101005213/16/16/11000103101005213/16/16/110012/32/10103/46/136/10012故A1=3/16/16/112/32/13/46/136/1.在事先不知道n阶矩阵是否可逆的情况下,也可以直接用此方法.如果在初等变换过程中发现左边的矩阵有一行元素全为0,则意味着A不可逆,因为此时表明A=0,则A1不存在.例2求A=987654321.解[AE]=1009870106540013211071260014630001321121000014630001321.由于左端矩阵中有一行元素全为0,于是它不可逆,因此A不可逆.3.伴随阵法定理n阶矩阵A=[aij]为可逆的充分必要条件是A非奇异.且A1=A1nnnnnnAAAAAAAAA.....................212221212111其中Aij是A中元素aij的代数余子式.矩阵nnnnnnAAAAAAAAA.....................212221212111称为矩阵A的伴随矩阵,记作A3,于是有A1=A1A3.证明必要性:设A可逆,由AA1=I,有1AA=I,则A1A=I,所以A0,即A为非奇异.充分性:设A为非奇异,存在矩阵3B=A1nnnnnnAAAAAAAAA.....................212221212111,其中AB=nnnnnnaaaaaaaaa.....................212222111211A1nnnnnnAAAAAAAAA.....................212221212111=A1AAAA...00.........0...00...0=1...00...1......0...100...01=I同理可证BA=I.由此可知,若A可逆,则A1=A1A3.用此方法求逆矩阵,对于小型矩阵,特别是二阶方阵求逆既方便、快阵,又有规律可循.因为二阶可逆矩阵的伴随矩阵,只需要将主对角线元素的位置互换,次对角线的元素变号即可.若可逆矩阵是三阶或三阶以上矩阵,在求逆矩阵的过程中,需要求9个或9个以上代数余子式,还要计算一个三阶或三阶以上行列式,工作量大且中途难免出现符号及计算的差错.对于求出的逆矩阵是否正确,一般要通过AA1=I来检验.一旦发现错误,必须对每一计算逐一排查.4.分块矩阵求逆法4.1.准对角形矩阵的求逆命题设A11、A22都是非奇异矩阵,且A11为n阶方阵,A22为m阶方阵221100AA12211100AA证明因为A=221100AA=11A22A0,所以A可逆.设A1=WZYX,于是有WZYX221100AA=mnII00,4其中XA11=In,YA22=0,ZA11=0,WA22=Im.又因为A11、A22都可逆,用A111、A122分别右乘上面左右两组等式得:X=A111,Y=0,Z=0,W=A122故A21=12211100AA把上述结论推广到每一个子块都是非奇异矩阵的准对角形状矩阵中去,即:121...kAAA=11211...kAAA4.2.准三角形矩阵求逆命题设A11、A22都是非奇异矩阵,则有12212110AAA=122122121111110AAAAA证明因为2212110AAAIAAI012111=221100AA两边求逆得1121110IAAI12212110AAA=12211100AA所以12212110AAA=IAAI01211112211100AA=122122121111110AAAAA同理可证12221110AAA=122122211111110AAAAA此方法适用于大型且能化成对角子块阵或三角块阵的矩阵.是特殊方阵求逆的一种方法,并且在求逆矩阵之前,首先要将已给定矩阵进行合理分块后方能使用.5.恒等变形法5恒等变形法求逆矩阵的理论依据为逆矩阵的定义,此方法也常用与矩阵的理论推导上.就是通过恒等变形把要求的值化简出来,题目中的逆矩阵可以不求,利用AA1=E,把题目中的逆矩阵化简掉。例1计算(A+4E)T(4E-A)1(16E-A2)的行列式,其中A=141021001解令)16()4()4(21AEAEEAT=DD=)16()4()4(21AEAEEAT=)4)(4()4()4(1AEAEAEAET=TAEAE)4)(4(=2)4(AE.虽然题目中出现了(4E-A)1.但是经过化简之后不再出现此式,因此得D=24AE=22500.例2已知n阶矩阵A满足A2+2A-3E=0.求证:A+4E可逆并求出A+4E的逆.证明把A2+2A-3E=0变形为A2+2A-8E=5E,即(A+4E)(A-2E)=-5E,可得(A+4E)(-A/5+2E/5)=E,所以存在一个矩阵B=-A/5+2E/5,使(A+4E)B=E,由定义得A+4E可逆,且(A+4E)1=B=-A/5+2E/5.另外,有些计算命题中虽出现逆矩阵,但通过适当的矩阵运算可消去,因而不必急于求出逆矩阵.6.利用线性方程组求逆矩阵若n阶矩阵A可逆,则AA1=E,于是A1的第i列是线性方程组AX=E的解,i=1,2,…,n,E是第i个分量是I的单位向量.因此,我们可以去解线性方程组AX=B,其中B=(b1,b2,…,bn)T,然后把所求的解的公式中的b1,b2,…,bn分别用E1=(1,0,0,…,0),E2=(0,1,0,…,0),……,6En=(0,0,0,…,1)代替,便可以求得A1的第1,2,…n列,这种方法在某些时候可能比初等变换法求逆矩阵稍微简一点.下面例子说明该方法的应用.例求矩阵A=3000013000013000013000013的逆矩阵.解设X=(x1,x2,x3,x4,x5)T,B=(b1,b2,b3,b4,b5)T解方程组AX=B,即:5545434323212133333bxbxxbxxbxxbxx解得:51554245432335432234254322314513)3(3)33(3)333(3)3333(3bxbbxbbbxbbbbxbbbbbx然后把B=(b1,b2,…,bn)列,分别用E1=(1,0,0,…,0),E2=(0,1,0,…,0),……,En=(0,0,0,…,1)代入,得到矩阵A1的第1,2,3,4,5列,分别用x1=(13,0,0,0,0)T,x2=(32,13,0,0,0)T,x3=(33,32,13,0,0)T,7x4=(34,33,32,13,0)T,x5=(35,34,33,32,13)TA1=1213214321543213000033000333003333033333.这种方法特别适用于线性方程组AX=B比较容易求解的情形,也是很多工程类问题的解决方法.以上各种求逆方法只是我的一些粗浅的认识,也许有很多的不当之处,我希望我的这篇文章能给大家带来帮助,能帮助我们更快更准地解决好繁琐的求逆矩阵问题.同时,它还是我们更好的学习线性代数的必备基础知识,认真掌握它,可供我们以后继续在数学方面深造打下坚实的基础.但我很希望各位老师和同学给于指导.能使我的这篇文章更加完善和实用.参考文献[1]北京大学数学系几何与代数教研室代数小组.高等代数[M].北京:高等教育出版社,2001.[2]杨明顺.三角矩阵求逆的一种方法[J].渭南师范学院学报,2003.[3]丘维声.高等代数[M].北京:高等教育出版社,2001.[4]杨子胥.高等代数习题集[M].济南:山东科学技术出版社,1984.[5]赵树原.线性代数[M].北京:中国人民大学出版社,1997.[6]李宗铎.求逆矩阵的一个方法[J].数学通报,1983.[7]贺福利等.关于矩阵对角化的几个条件[J].高等函授学报(自然科学版),2004,(1)[8]张禾瑞.郝炳新.高等代数[M].北京:高等教育出版社.1999.[9]王永葆.线性代数[M].长春:东北大学出版社.2001.[10]同济大学遍.线性代数(第二版).北京:高等教育出版社,1982.[11]王萼芳,丘维声编,高等代数讲义.北京大学出版社,1983.[13]华东师范大学数学系编.数学分析.人民教育出版社,1980[14]杜汉玲求逆矩阵的方法与解析高等函授学报(自然科学版)第17卷第4期2004年8月8[15

1 / 9
下载文档,编辑使用

©2015-2020 m.777doc.com 三七文档.

备案号:鲁ICP备2024069028号-1 客服联系 QQ:2149211541

×
保存成功