




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、10/25/20211线性代数第7讲分块矩阵10/25/20212把一个5阶矩阵2110112230001000001000001a用水平和垂直的虚线分成4块.120010021101, 00, 0101222000001aaoi121aaoi10/25/20213把一个mn矩阵a, 在行的方向上分成s块, 在列的方向分成t块, 称为a的st分块矩阵, 记作a=aklst, 其中akt(k=1,2,.,s,l=1,2,.,t)称为a的子块, 它们是各种类型的小矩阵.常用的分块矩阵, 除了上面的4块矩阵, 还有以下几种形式:10/25/20214按行分块1112112122221212,1,2,
2、.nnmmmnmiiiinaaaaaaaaaaaaaaaaaim10/25/20215按列分块1112121222121212,1,2, .sssnnnsjjjnjbbbbbbbb bbbbbbbbjsb10/25/20216当n阶矩阵c中非零元素都集中在主对角线附近, 有时可以分块成对角块矩阵(准对角矩阵)12,ncooococooc其中ci是ri阶方阵(i=1,2,.,m; r1+r2+.+rm=n)10/25/20217例如 12312301000012000000110000112000022000000311001,112,3 .12022bbbbbbb 10/25/20218下面讨
3、论分块矩阵的运算1. 分块矩阵的加法设分块矩阵a=aklst, b=bklst, 如果a与b对应的子块akl和bkl都是同型矩阵, 则a+b=akl+bklst例如11121112111112122122212221212222aabbababaabbabab其中a11与b11, a12与b12, a21与b21, a22与b22分别都是同型小矩阵(子块).10/25/202192. 分块矩阵的数量乘法设分块矩阵a=aklst, h是一个数, 则ha=haklst.10/25/2021103. 分块矩阵的乘法 设a是mn矩阵, b是np矩阵, 如a分块为rs分块矩阵aklrs, b分块为st分
4、块矩阵bklst, 且a的列的分块法和b的行的分块法完全相同, 则10/25/202111可以证明(但略去), 用分块乘法求得的ab与不分块作乘法求得的ab是相等的.12111211112112122221222212121,(1,2, ;1,2, )ststsssstrrrssklr tsklkiilijjjbbbaaajbbbaaajabbbbaaajccca bkr lt列列列行行行10/25/202112例1 将下列5阶矩阵a,b分成4块阵, 并用分块矩阵的乘法计算ab.10000320100100013001,121001000011010010002000100100ab10/25
5、/202113解 由观察, 可将a分成如下4块阵22 31311000001000,121001101020001121120ioaaia10/25/202114将b分块为1223 213201013001,100000100000100320130bibiob10/25/202115则146044042100010001031023021121,311131121321312ibaaibaiboiibiaoiab其中10/25/202116故3201013001240124401164120ab 10/25/202117例2 设a是mn矩阵, b是ns矩阵, b按列分块成1s分块矩阵, 将a
6、看成11分块矩阵, 则 ab=ab1,b2,.,bs=ab1,ab2,.,abs.若已知ab=0, 则显然有abi=0, i=1,2,.,n. 因此, b的每一列都是线性方程组ax=0的解.10/25/202118例3 若n阶矩阵c,d可以分块成同型对角块矩阵, 即1122,mmcdcdcdcd1122mmc dc dcdc d其中ci和di是同阶方阵(i=1,2,.,m), 则10/25/202119例4 证明:若n阶上三角矩阵a可逆, 则其逆a1也是上三角阵.证 对n作数学归纳法, n=1时, a1=1/a, 结论成立.(一阶矩阵可认为是上三角矩阵).假设命题对n1阶可逆上三角阵成立, 考
7、虑n阶情况, 设11121222111000(0,1,2, ),nnnniiaaaaaaaoaaain10/25/202120其中a1是n1阶可逆的上三角阵. 设a的逆矩阵为111212122211112,nnnnnnbbbbbbbbbbbb111111111 11111111111nnnoababioioaboa baboiaab则10/25/202121于是: a1=o=a1o=oa1b1=in1=b1=a11, 根据归纳假设, b1是n1阶上三角矩阵, 因此1111babob是上三角矩阵(其中b11=a111; a111a11).10/25/2021224. 分块矩阵的转置.,.,.,
8、2 , 1;, 2 , 1,2121231322122111232221131211tmtttmttttttttkllkstlkttsklbbbbbbbbaaaaaaaaaaaaaasktlabbaaa则则例如其中的转置矩阵为分块矩阵按行分块10/25/2021235. 可逆分块矩阵的逆矩阵对角块矩阵(准对角矩阵)12maaaa111121maaaa的行列式为|a|=|a1|a2|.|am|, 因此, 对角块矩阵a可逆的充要条件为|ai|0 i=1,2,.,m, 这时10/25/202124例51,aadbdcoba可逆并求证皆为可逆方阵其中设1,xyazt解 |a|=|b|d|0, a可逆.
9、 设其中x与b, t与d分别是同阶方阵, 于是由12boxybxbycdztcxdzcydtiooi10/25/202125bx=i1, 故 x=b1.by=o, 故 y=b1o=o.cx+dz=o, 故 dzcxcb1,zd1cb1.cy+dt=i2,故 dt=i2, t=d1.所以11111.boad cbd12iobxbycxdzcydtoi10/25/202126例如033021,2112,10003010300012100021000123cidbdcoba可逆其中10/25/202127则1001201021001100003/23/10003/13/2122110,2112311
10、11111311dcbdobacbdidb即得10/25/2021286. 分块矩阵的初等变换与分块初等阵这里仅就22分块矩阵为例来讨论. 对于分块矩阵11122122aaaaa12mncoiooioc或可以同样定义它的三种初等行变换和列变换, 并相应地定义三类分块初等矩阵:(i) 分块倍乘阵(c1,c2是可逆阵)10/25/202129(ii) 分块倍加阵43.mmnnioiccioi或.nmoiio(iii)分块对换阵分块初等矩阵是方阵, 它们左乘(或右乘)分块矩阵a(不一定是方阵), 在保证可乘的情况下, 其作用与前述初等矩阵相同.10/25/202130例6 设n阶矩阵a分块表示为11
11、121112,iopa ai解 先对分块阵a作初等变换, 将其化为上三角块矩阵, 为此左乘分块倍加阵.,.,122111212211221122211211aaaaaaaaaaaaaa并求可逆证明可逆和且为方阵其中其中i1,i2为单位阵, 其阶数分别与a11,a22相同.10/25/202131于是1211121221211121221112111212211112111212211112111212212111., 0|0| , 0| , 1|. |,aaaaqaaaaaaaaaaaapaaaaaapbaaaaoaaap记可逆故所以由于记作10/25/202132为求a1, 将b化为对角块矩阵, 为此取211121121121111112111112112112211212)(,)(,iaaoiioqaiqooappcacppacappcqooabpioqaip因此两边取逆得即于是记作10/25/202133记.,12111212211121121112111121221112111211111111222112111aaaaqaaqdqaadqdaaqaaaddddda其中
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
评论
0/150
提交评论