线性代数课件:2-4分块矩阵
- 格式:ppt
- 大小:1.30 MB
- 文档页数:25
§4 矩阵分块法本节我们将介绍矩阵运算的一种有用的技巧——矩阵的分块,这种技巧在处理某些较高阶的矩阵时常常被用到。
一、分块矩阵的概念设A 是一个矩阵,我们在它的行或列之间加上一些直线,把这个矩阵分成若干个小块,例如,设A 是一个43⨯矩阵111213212223313233414243a a a a a a A a a a a a a ⎛⎫ ⎪ ⎪= ⎪ ⎪⎝⎭ 我们可以把它分成如下的四块111213212223313233414243a a a a a a A a a a a a a ⎛⎫ ⎪ ⎪= ⎪ ⎪⎝⎭用这种方法被分成若干个小块的矩阵称为分块矩阵,每一个小块称为A 的一个子块。
在一个分块矩阵中,每一个小块也可以看成是一个矩阵。
例如,上面的分块矩阵A 是由以下四个矩阵组成的111121a A a ⎛⎫= ⎪⎝⎭ 1213122223a a A a a ⎛⎫= ⎪⎝⎭ 312141a A a ⎛⎫= ⎪⎝⎭ 3233224243a a A a a ⎛⎫= ⎪⎝⎭我们可以把A 简单地写成11122122A A A A A ⎛⎫=⎪⎝⎭对一个矩阵来讲,可以有各种不同的分法。
二、分块矩阵的运算规则分块矩阵的运算规则与普通矩阵的运算规则相类似,分别说明如下:(1)分块矩阵的加法设()ij m n A a ⨯=,()ij m n B b ⨯=,采用同样的分块方法得1111r s sr A A A A A ⎛⎫ ⎪= ⎪ ⎪⎝⎭ , 1111r s sr B B B B B ⎛⎫ ⎪= ⎪ ⎪⎝⎭其中ij A 与ij B 的行数与列数都相同,则11111111r r s s sr sr A B A B A B A B A B ++⎛⎫ ⎪+= ⎪ ⎪++⎝⎭(2)数乘分块矩阵设1111r s sr A A A A A ⎛⎫⎪= ⎪ ⎪⎝⎭ ,λ为实数,则1111r s sr A A A A A λλλλλ⎛⎫⎪= ⎪ ⎪⎝⎭(3)分块矩阵的乘法设()ij m l A a ⨯=,()ij l n B b ⨯=,分别分块成1111t s st A A A A A ⎛⎫ ⎪= ⎪ ⎪⎝⎭ , 1111r t tr B B B B B ⎛⎫⎪= ⎪ ⎪⎝⎭其中12,,i i it A A A (1,2,,i s = )的列数分别等于12,,,j j t j B B B (1,2,,j r = )的行数,则1111r s sr C C AB C A ⎛⎫ ⎪= ⎪ ⎪⎝⎭其中1tij ik kj k C A B ==∑(1,2,,i s = ,1,2,,j r = )例1 设1000010012101101A ⎛⎫ ⎪⎪= ⎪- ⎪⎝⎭, 1010120110411120B ⎛⎫ ⎪-⎪= ⎪ ⎪--⎝⎭求乘积AB解 为了求乘积AB ,我们可以对A 、B 进行如下的分块1000010012101101A ⎛⎫ ⎪ ⎪= ⎪- ⎪⎝⎭1E O A E ⎛⎫= ⎪⎝⎭,1010120110411120B ⎛⎫ ⎪- ⎪= ⎪ ⎪--⎝⎭112122B E B B ⎛⎫= ⎪⎝⎭按分块矩阵的乘法可得11111212211121122E O B E B EAB A E B B A B B A B ⎛⎫⎛⎫⎛⎫== ⎪⎪⎪++⎝⎭⎝⎭⎝⎭而 11121121010111211A B B -⎛⎫⎛⎫⎛⎫+=+⎪ ⎪ ⎪---⎝⎭⎝⎭⎝⎭2411-⎛⎫= ⎪-⎝⎭122124133112031A B -⎛⎫⎛⎫⎛⎫+=+= ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭故 1010120124331131AB ⎛⎫ ⎪- ⎪= ⎪- ⎪-⎝⎭(4)分块矩阵的转置设1111r s sr A A A A A ⎛⎫⎪= ⎪ ⎪⎝⎭ , 则1111T T s T T T r srA A A A A ⎛⎫ ⎪= ⎪ ⎪⎝⎭(5)分块对角阵在n 阶方阵A 的分块矩阵中,如果只有在主对角线上有非零的小方阵,而其余子块均为零矩阵,即12s A A A A ⎛⎫ ⎪⎪= ⎪ ⎪⎝⎭则A 称为分块对角阵。