分块矩阵word版本
3.3 分块矩阵
![3.3 分块矩阵](https://img.taocdn.com/s3/m/7e1992387f1922791688e8e1.png)
s1
Csr
其 中Cij Aik Bkj i 1,, s; j 1,, r .
k 1
9
例4 理解常用分块
2 3 4 5
1
AB
4
2 5
3 6
3 4
4 5
5 6
6
7
20
47
26 32 62 77 24
38
92
(1)B 按列分块
AB A B1 B2 B3 B4 AB1 AB2 AB3 AB4
作业: P132 22 24 27 P133 32 33
19
a10
A
0 1
a 0
0 b
0 1
A1
A2
A3
A4 ,
其中A2431
a0 01b
0 1 1 b
1b0
a
A
0 1
0
1 a 0 1
0 0 b 1
0 0 1 b
B1
B2 B3
.
B4
3
二、分块矩阵的运算
1. 加(减)法
设矩阵A和B的行数和列数相同, 采用相同的分
An
Bn
AnBn
16
矩阵分块运算后以块作为元素进行运算的基本原则是: 1° 分块计算与不分块计算的结果相同; 2° 两种计算方法的法则一致; 3 °块也是矩阵,块之间的运算仍然是矩阵的运算.
17
三、 应用
命题1 AB=O充要条件是 B的列向量是齐次线性方程组AX=O的解 .
利用分块矩阵分析式 AmnBns Oms的涵义.
A1t kA11
Ast kAs1
kA12 kAs 2
kA1t
kAst
第一章 矩阵分块矩阵及其运算
![第一章 矩阵分块矩阵及其运算](https://img.taocdn.com/s3/m/b91cb3ef87c24028905fc36a.png)
所以
1 0 1 0
A BA 2B 1B 111 1B21 A 21 E B 22
1 2 1
2 4 1
0 3 3
1
3 1
第一章 矩阵
4. 分块转置
设矩阵A =
A11 A12 … A1r A21 A22 … A2r ………… ,
As1 As2 … Asr
§3 分块矩阵
A11T A21T … As1T
A=
A21 A22 … A2t …………
,
B=
B21 B22 … B2r …………
,
As1 As2 … Ast
Bt1 Bt2 … Btr
其中Ai1, Ai2, …, Ait的列数分别与B1j, B2j, …, Btj的 行数相等.
C11 C12 … C1r
则AB =
C21 C22 … C2r … … ……
从而 X 有 A 2 1 A 2 A 1 -21 11 A 1 11 A O 2 12 ,使AX=E 所以,A可逆,且A-1 =X。
例 利用上例结论求方阵
0 0 1 0 0
0 2 0 0 0
A1Biblioteka 3 110 2 0
0 0 0 的逆阵
1 1
1 0
11
解 记AAA1211
O
A2
2
0 2 0 1
1 0 4 2
0
1
1 0
B11 B21
E B22
所 以ABAE21 O EB B1211 BE22
A2
B11 1B11B21
E
A21B2
2
因为 A21B11B2111 1211 0211 01 03 4211 0112 14
(完整版)分块矩阵及其应用汇总,推荐文档
![(完整版)分块矩阵及其应用汇总,推荐文档](https://img.taocdn.com/s3/m/8225b1cdc5da50e2534d7f42.png)
分块矩阵及其应用徐健,数学计算机科学学院摘要:在高等代数中,分块矩阵是矩阵内容的推广. 一般矩阵元素是数量,而分块矩阵则是将大矩阵分割成小矩形矩阵,它的元素是每个矩阵块.分块矩阵的引进使得矩阵工具的利用更加便利,解决相关问题更加强有力,所以其应用也更广泛. 本文主要研究分块矩阵及其应用,主要应用于计算行列式、解决线性方程组、求矩阵的逆、证明与矩阵秩有关的定理.关键词:分块矩阵;行列式;方程组;矩阵的秩On Block Matrixes and its ApplicationsXu Jian, School of Mathematics and Computer ScienceAbstract In the higher algebra, block matrix is a generalization of matrix content.In general, matrix elements are numbers. However, the block matrix is a large matrix which is divided into some small rectangular matricies, whose elements are matrix blocks. The introduction of the block matrix makes it more convenient to use matrix, and more powerful to solve relevant problems. So the application of the block matrix is much wider. This paper mainly studies the block matrix and its application in the calculation of determinant, such as solving linear equations, calculating inverse matrix, proving theorem related to the rank of matrix , etc.Keywords Block matrix; Determinant; System of equations; Rank of a matrix11 ⎪1 引 言我们在高等代数中接触到矩阵后,学习了矩阵的相关性质,但是对于一些复杂高阶矩阵,我们希望能将问题简化. 考虑将矩阵分割为若干块,并将矩阵的部分性质平移至分块矩阵中,这样的处理往往会使问题简化.定义 1.1 [1] 分块矩阵是把一个大矩阵分割成若干“矩阵的矩阵”,如把 m ⨯ n 矩阵分割为如下形式的矩阵:⎛A 11A ⎫ 1n ⎪A m ⨯n = ⎪A m 1 A m n特别地,对于单位矩阵分块:⎝ ⎭ ⎛E 0 0 ⎫ ⎪ E n ⨯n = 0 0 0 ⎪ 0 E ⎝n n ⎭ 显然,这里我们认识的矩阵元素不再局限于数字,而是一个整体,这里的A 所代表的是大矩阵囊括的小矩阵,而小矩阵一般是我们熟知的常见矩阵.ij依照以上设想,有关矩阵性质的一些问题,我们可以考虑用分块矩阵的思路来解决.2.1 矩阵的相关概念2 分块矩阵在矩阵的学习中,我们学过一些最基本的概念,比如矩阵的行列式、矩阵 的秩、矩阵的逆、初等变换、初等矩阵等等.事实上,我们发现:分块后的矩阵同样用到这些概念.a 11 定义 2.1.1[2]n 级行列式a 21a 12 a 22 a 1n a 2n等于所有取自不同行不同列的a n 1 a n 2a nn 个元素的乘积a 1j a 2ja n j的代数和,这一定义又可写成:12na 11 a 21 a 12a 22a 1na 2n =(-1) (j 1j 2 j n )a aa .a n 1 a n 2a n∑j 1j 2 j n1j 1 2j 2n j n[2]定义 2.1.2向量组的极大无关组所含向量的个数称为这个向量组的的秩.所O I ⎪ ⎪ ⎪1谓矩阵的行秩就是指矩阵的行向量组的秩;矩阵的列秩就是矩阵列向量组的秩. 定义 2.1.3 [2] n 级方阵称为可逆的,如果有n 级方阵 B ,使得A B = A -1 .BA = E (这里 E 是n 级单位矩阵),那么B 就称为 A 的逆矩阵,记为定义 2.1.4 [3] 对分块矩阵施行下列三种初等变换: (1) 互换分块矩阵的某两行(列);(2) 用一个非奇异阵左(右)乘分块矩阵的某一行(列);(3) 用一个非零阵左(右)乘分块矩阵的某一行(列)加至另一行(列)上, 分别称上述三种初等行(列)变换为分块矩阵的初等行(列)变换. 定义 2.1.5 [3] m + n 2 ⨯ 2 ⎛I m O ⎫对 阶单位矩阵作 分块,即I m +n = O I ⎪ ,然后⎝ n ⎭对其作相应的初等变换所得到的矩阵称为分块初等矩阵. 分块矩阵具有以下形式:(1) 分块初等对换阵⎛I n O ;⎫ ⎝ m ⎭⎛P O ⎫ ⎛I m O ⎫(2) 分块初等倍乘阵 0 I ⎪ , ⎪ ;⎝ n ⎭ (3) 分块初等倍加阵⎛I m R 1 ⎫ O I ⎝ 0 Q ⎭ ,⎛I m O ⎫ ; S I ⎝ n ⎭ ⎝ n ⎭其中 P , Q 分别是m 阶和n 阶可逆方阵,且R ∈ R m ⨯n ,S ∈ R n ⨯m为非零阵.2.2 矩阵的运算性质矩阵的运算包括加法、乘法、数乘,这里主要讨论矩阵的运算性质: 定义 2.2.1 [4] 矩阵加法:设A = (a ) , B = (b ) 是两个同型矩阵,ij snij sn则矩阵C = (c i j )= (a i j+ b i j )称为 A 和 B 的和,记为C = A + B .元素全为零的矩阵称为零矩阵,记为O s n ,可简单记为O,对于矩阵 A 、 B ,有:(1) A + O = A(2) A + ( -A ) = 0(3) A - B = A + ( -B )(4) ( A + B ) + C = A + ( B + C )snsnn11 (5)A + B = 定义 2.2.2 [4] B + A矩阵乘法:设A = (a ) ,B = (b ) 是两个不同型矩阵,i k s nk j n m那么矩阵C = A B =(c i j ),称为矩阵 A 与 B 的乘积,其中:smc i j = a i 1b 1j + a i 2b 2j+ a i n b n j= ∑a i k b k jk =1在乘积的定义中,我们要求第二个矩阵行数和第一个矩阵列数相等.特别地,矩阵的乘法和加法满足以下性质:(1) A ( B + C ) = A B + A C(2) ( B + C )A = B A + C A(3) (A B )D =A (B D )⎛k a 11 k a 1k a 1 ⎫定义 2.2.3 [4] 矩阵数乘: k a 21k ak a 2n ⎪ ⎪A = (a ) 与 数 22 ⎪称为矩阵 ⎪⎪ ij sn k a k a k a ⎝ s 1 s 2 s n ⎭k 的数量乘积,记为kA ,有以下性质:(1) 1 * A = A ;(2) k(l A ) = (k l )A ;(3) k ( A + B )= kA + kB ;(4) (k + l )A = kA +lA ; (5) k (A + B ) = kA +kB .2.3 分块矩阵的初等变换性质我们对于分块矩阵,也有其运算性质:设 A 、 B 是m ⨯ n 矩阵,若对它们有相同的划分,也就有:⎛A 11 + B A 1t + B 1t ⎫ ⎪ 加法:A + B = ⎪ . ⎪ A + B A + B ⎪ ⎝ s 1 s 1 st st ⎭乘法:C = A B , 其中:∑ ⎪ 1 C i j = A i 1B 1j + A i 2B 2j+ + A i n B n j⎛k A 11k A 1 ⎫⎪ n= A i k B k j .k =1数乘:k A =⎪ .⎪ k Ak A⎝s 1 s t ⎭总结了矩阵的运算性质,我们主要看看分块矩阵初等变换性质:定义 2.3.1 [2] 由单位矩阵 E 经过一次初等变换得到的矩阵称为初等矩阵. 初等矩阵都是方阵,包括以下三种变换:(1) 互换矩阵 E 的i 行与 j 行的位置; (2) 用数域 P 中的非零数c 乘 E 的i 行; (3) 把矩阵 E 的 j 行的k 倍加到i 行.定义 2.3.2 [5] 将单位矩阵分块,并施行如下三种变换中的一种变换而得到的方阵称为分块初等矩阵:(1) 对调两块同阶的块所在的行或列; (2) 某一块乘以同阶的满秩方阵;(3) 某一块乘以一个矩阵后加到另一行上(假定这种运算可以进行).如:我们对分块矩阵⎛ A B ⎫进行相应变换,只要应用矩阵的计算性质,左乘对⎝C D ⎭ 应分块矩阵: ⎛ O E m ⎫ ⎛ A B ⎫ ⎪⎪⎛C D ⎫ ⎪ ⎝E n O ⎭ ⎝C D ⎭⎝ A B ⎭ ⎛P O ⎫ ⎛ A B ⎫ ⎛P A = P B ⎫ O E ⎪C D ⎪ ⎪⎝ n ⎭ ⎝⎭ ⎝ C D ⎭ ⎛E m O ⎫ ⎛ A B ⎫ ⎛ = A B⎫P E ⎪C D ⎪ ⎪C + P AD + P B⎝ n ⎭ ⎝⎭ ⎝ ⎭2.4 矩阵的分块技巧对矩阵的分块不是唯一的,我们往往根据问题的不同进行不同的分块,分块的合适与否,都对问题的解决至关重要,最常见的有四种分块方法[6] :(1) 列向量分法,即A =(1,⎛ ⎫ ⎪, n ),其中j 为 A 的列向量.(2) 行向量分法,即A = ⎪ ,其中j 为 A 的行向量.⎪ ⎝ m ⎭=1⎪ (3)分两块,即A = (A 1, A 2 ),其中A 1 ,A 2 分别为A 的各若干列作成.或 A = ⎛B ⎫ ,其中B ,B 分别为 A 的若干行作成. B ⎪1 2 ⎝ 2 ⎭⎛C 1 C 2 ⎫(4) 分四块,即A =C C ⎪ .⎝ 3 4 ⎭我们在进行分块时,希望分割的矩阵块尽可能是我们所熟悉的简单矩阵,于是,我们有必要熟悉一些常见的矩阵.2.5 常见的矩阵块我们把高等代数中学习过的一些常见矩阵总结如下: (1) 单位矩阵:对角线元素都为1,其余元素为0 的n 阶方阵. (2) 对角矩阵:对角线之外的元素都为0 的n 阶方阵. (3) 三角矩阵:对角线以上(或以下)元素全为0 的n 阶方阵. (4) 对称矩阵:满足矩阵 A 的转置和 A 相等. (5) 若尔丹(Jordan )块:形如⎛ 0 1 0 0 ⎫ 0 ⎪J ( ,t ) ⎪= ⎪0 0 ⎪ 0 0 0 1 ⎝ ⎭(6) 若尔丹形矩阵:由若干个若尔丹块组成的准对角矩阵, 其一般形状形如:⎛A 1 ⎫⎪ A 2⎪ ⎪ ⎪A ⎪ ⎝n ⎭在复杂矩阵中,找到这些矩阵块,会使计算简化.3.1 行列式计算的应用3 分块矩阵及其应用定理 3.1.1 [2] 拉普拉斯(Laplace )定理:设在行列式 D 中任意取定了k 个 行.由这k 行元素所组成的一切k 级子式与它们的代数余子式的乘积的和等于行列式 D .事实上,行列式计算中的拉普拉斯定理就包括了矩阵分块的思想,它通过取k 级子式的方法,提取出矩阵内的矩阵块. 然而,在行列式计算中,行列式a ⎪ a 按行或列的展开更为常用. 这里,我们最常用到的是取列向量分块和行向量分块.例 3.1.1 [7] :(爪形行列式)计算行列式:a 01 1 1 1 a 10 0 1 0 a 2 0 ,其中a i ≠ 0(i = 1, 2, , n ) .1 0 0 a n解:设Q =A D ,其中A = (a )C B a 1 B =,C = ( 1, 1, , 1)T ,D = ( 1, 1, , 1) .a n因为a i ≠ 0(i = 1, 2, , n ) ,所以 B 是可逆矩阵.-1⎛n 1 ⎫又易知: A - D B C = a 0 - ∑ ⎪ . ⎝ i =1 i ⎭根据分块矩阵乘法: ⎛ E0 ⎫ ⎛ A D ⎫ --1 ⎪ ⎪= ⎛A D ⎫-1 ⎝ C A E ⎭ ⎝C B ⎭ ⎝ 0 B - C A D ⎭A D -1 -1 ⎛ n 1 ⎫则:= AB - C A D =B A - D BC = a a a a-∑ a ⎪C B⎛n 1 ⎫ 12n 0⎝i =1 i ⎭故:原行列式=a 1a 2 a n a 0 - ∑ ⎪ . ⎝ i =1 i ⎭例 3.1.2 [7] :(对角行列式)计算行列式:adH 2n= a d.c bcb解:令⎪ a x A =⎛a ⎫⎪ ,B = ⎛b ⎫⎪ ,C = ⎛ c ⎫ ⎛ ,D = d ⎫⎪ ⎪ ⎪ ⎪ ⎪ a ⎪ b ⎪ c ⎪ d ⎪ ⎝ ⎭ 为n 阶方阵. 由于a ≠ ⎝ ⎭ ⎝ ⎭ ⎝ ⎭ 0,故 A 为可逆方阵.⎛ b - c a -1d⎫⎪ 又易知:B - C A -1D =⎝ b - c a -1d ⎪ b - -1 ⎪ ca d ⎭故 H 2n= A D = C BAB - C A -1D = a n (b - c a -1d )n= (a b - c d )n .例 3.1.3 [8] :设 A 、 B 、C 、 D 都是n 阶矩阵,证明当 AC = CA 时, A 可逆时,有A D= A B - C DC B⎛ A D ⎫ ⎛E -A 1D-⎛ A 0 ⎪ ⎫,证明:若 A 可逆,⎪ ⎪ =-1 ⎝C B ⎭ ⎝OE ⎭ ⎝C B - C A D ⎭A D故:=C BAB - C A -1D = A B - A C A-1D = A B - C D .注意到,这里计算分块矩阵行列式和计算一般数字矩阵行列式有所区别,不是简单的a d c b= a b - c d ,其矩阵块限制条件有所加强. 所以本例告诉我们,在矩阵分块以后,并非所有一般矩阵性质都可以应用到分块矩阵中.3.2 线性方程组的应用对于线性方程组,我们有以下四种表述: (1) 标准型:⎧a 11x 1 + a 12x 2+ + ax = b ⎪ 1nn 1⎨ax + ax + + a x = b ; ⎪a 21 x 1+ 22 2 + + 2n n a x = b ⎩ m1 1 m2 2 m n n m (2) 矩阵型:令A = ⎣a i j ⎦m ⨯n,x = (x 1, x 2, , x n )' ,B = (b 1, b 2, b m )' 方程组可以表述为: Ax = B ;(3) 列向量型:令2⎢a ⎥ ⎝O O⎪ ⎪ ⎪ ⎡a 11 ⎤ ⎢21 ⎥⎡a 12 ⎤⎥ 22 ⎡a 1n ⎤ ⎢ ⎥ = , 1 ⎢ ⎥ 2 = , , ⎢ ⎥= ⎢a 2n ⎥ n ⎢ ⎥ ⎢ ⎥ ⎣a m 1 ⎦ ⎢ ⎥ ⎣a m 2 ⎦ ⎢ ⎥ ⎣a m n ⎦则方程组又可以表述为:x 11 + x22+ + x nn = B ;(4)行向量型: x ' + x ' + + x' = B ' .1 12 2n n可见,矩阵分块为我们解方程组提供了新的思路.事实上,在求齐次线性方程组系数矩阵的秩时,在判断非齐次线性方程组是否有解时,行列向量组的合理应用,使得问题解决更加便捷、明了.例 3.2.1:(齐次线性方程组)求解方程组:⎧ x 1 + 2x 2 2x ⎪ + x + 2x 3 - 2x + x 4 = 0 - 2x = 0 ⎨ 1 x -2x - 4x 3 - 3x 4=0 ⎩ 1 2 3 4 解:对系数矩阵施行行变换,并将结果用分块矩阵表示:⎛1 0 -25 ⎫ - 3⎪ ⎛ 1 2 2 1 ⎪⎫ ⎛ 1 2 2 1 ⎪⎫4 ⎪ ⎛E C ⎫ A = 2 1 -2 -2 0 -3 -6 -4 0 1 2 ⎪ = 2 ⎪ ⎪1 -1 -4 -3⎪ 0 -3 -6 -4⎪ 3 ⎪ 12 ⎭ ⎝ ⎭ ⎝ ⎭ 0 0 0 0 ⎪⎪ ⎝ ⎭R ( A ) = 2,基础解系含4 - 2 = 2 个.而方程又满足:相应的可以取:⎛E 2 C ⎫ ⎛1 ⎫ = ⎛ 0⎫⎪ ,⎝O 1 O 2 ⎭ ⎝2 ⎭⎝ 0⎭⎛ 5 ⎫ 2 3 ⎪ ⎛ -C ⎫⎪⎝ E 2 ⎭⎪ = -2 4 ⎪3 ⎪1 0 ⎪ ⎝ 0 1 ⎭-⎪ 0 3 ⎪⎭⎛ 2 ⎫ ⎛ 5 ⎫3 ⎪有通解: = k + k,其中= -2⎪1, =- ⎪ 4 ⎪ . 1 12 21 ⎪2 ⎪ ⎪ ⎝ 0 ⎭⎪ 1 ⎪ ⎝ ⎭例 3.2.2 [9] :(非齐次线性方程组)求解方程组:⎧⎪ x 1 + 2x 2- 3x 4 + 2x 5 = 1 x - x - 3x + x - 3x = 2 ⎪ ⎨ 1 2 3 4 52x - 3x + 4x - 5x + 2x = 7 ⎪ 9x ⎩ 1= 25 解:我们分别对于方程组的系数矩阵和增广矩阵求秩:r ( A ) = 3,而r ( A ) = 4 , 故r ( A ) ≠ r ( A) . 从而方程组无解. ⎛ Λ45 -b ⎫事实上,我们可以利用分块矩阵叙述:经对分块矩阵 ⎝ E变换,都不能把最后一列变成0 ,所以该方程组无解.例 3.2.3:证明: n 阶方阵 A 的秩为n- 1,则r a n k ( A* )=1首先证明此例需要利用的一个引理: 4进行行列0 引理:A = (a i j )n ⨯n ,B = (b i j )n ⨯n ,r( A ) = r ,A B =0 ,则r ( B ) ≤ n - r证明:对矩阵 B 进行列向量的分块,B = (B 1, B 2, B n ) ,A B = 0 则有:A B i= 0 ,B i 是AX = 0 的解. 而A X =0 基础解系有n - r 个解.故:r ( B ) ≤ n - r 再证明本例: 因为r ( A )= n - 1,则 A = 0 ,A 至少有一个n -1级子式不为零,r a n k ( A* ) ≥ 1.而:A * =AE = 0 .利用引理得:r a n k ( A * ) ≤ 1,故r a n k ( A )=*.51 - 9 x +2 6x - 163 x4 + 2x 52 3 4 5⎝⎪ 1 2= ⎪ ⎪ 得证.3.3 求矩阵逆的应用我们在求矩阵逆的时候包括很多方法:利用定义求逆、利用伴随矩阵求逆、 利用初等变换求逆、混合采用初等行列变换求逆等等.这里我们统一用矩阵分块的思路来求矩阵的逆.例 3.3.1 [6] :设 A 、 B 是n 阶方阵,若 A + B 与 A - B 可逆,试证明: ⎛ A B ⎫可逆,并求其逆矩阵. B A ⎭ ⎪ 解:令D = ⎛ A B ⎫,由假设知 A + B ≠ 0 , A - B ≠ 0B A ⎪ .那么:D =A B⎝ ⎭A +B B =A + BB= A + B A - B ≠ 0 .B AB + A AA - B即 D 可逆. 再令D -1 ⎛D 1= D 2⎫ , 由D -1 = E ,即:可得:D D ⎝ 3 4 ⎭⎛ A B ⎫ ⎛D D ⎫ ⎛E 0 ⎫ ⎪ ⎪⎪ ⎝B A ⎭ ⎝D 3D 4 ⎭ ⎝ 0E ⎭⎪⎧A D 1 + B D 3 = E B D + A D = 0⎪12⎨A D +B D = 0 B D 2 + A D 4 = E ⎩ 2 4将第一行和第二行相加、相减,得:⎪D + D = ( A + B )-1 ⎨1 3⎩D 1 - D 3= ( A - B )-1 解之得:D = 1 ⎡( A + B )-1 + ( A - B )-1 ,D = 1⎡( A + B )-1 - ( A - B )-11 2 ⎣⎦ 2 2 ⎣⎦类似地:D 2所以: = D 3 ,D 4= D 1 .⎛ A B ⎫-11 ⎛( A + B )-1 + ( A - B )-1 ( A + B )-1 - ( A - B )-1 ⎫⎪ = 2 -1 -1 -1-1 ⎪ . ⎝B A ⎭ ⎝( A + B ) - ( A - B )( A + B ) + ( A - B ) ⎭ =⎝⎭ ⎝ - ⎪⎪ ⎪0 例 3.3.2 [6] :已知分块形矩阵M = ⎛ A B ⎫可逆,其中 B 为p ⨯ p 块, C 为C 0 ⎪ ⎝ ⎭q ⨯ q 块,求证: B 与C 都可逆,并求M-1 . 解:由0 ≠M = (-1)p qBC ,则: B ≠0 , C ≠ 0 ,即证 B 、C 都可逆.这里用分块矩阵的广义初等变换来求逆: ⎛ A B E p0 ⎫ → ⎛ A B E 0 ⎫ → ⎛ 0B E -AC -1 ⎫⎪ ⎪ -1 ⎪ -1⎝C 0 0 Eq ⎭ ⎝E 0 0 C ⎭ ⎝E 0 0 E ⎭→ ⎛ 0 E B -1-B -1A C -1 ⎫ → ⎛E 0 0 C-1 ⎫E 0 0 C-1⎪ 0 E B -1-B -1A C -1 ⎪ ⎭-1⎛C -1 ⎫故 :M = B -1-B -1A C-1 ⎪ . ⎝⎭备注:本例和上例属于同一个类型的问题,但我们利用分块矩阵,可以有两种不同的方法来解决,待定系数法和广义初等变换都是求逆的有效方法.值得注意的是,在题目没有直接给出分块矩阵的情况时,我们要学会自己构造:⎛ 1 0 1 ⎫ 例 3.3.3 [10] :求矩阵A = 2 1 0 ⎪的逆矩阵.⎝ ⎭ 解:构造矩阵:⎛ 10 1 1 00⎫⎪⎛ 1 0 1 1 0 0⎫⎪2 0 0 1 -2 -2 1 0 D = ⎛ A E ⎫= -3 1 0 0 1 2 -5 0 0 1⎪ → 0 2 -2 3 0 1⎪ ⎪⎪ ⎪ ⎝E O ⎭6⨯6 1 0 0 0 00 1 0 0 0 0⎪ 1 0 0 0 0 0⎪ 0⎪ 0 1 0 0 0 0⎪0 0 1 0 0 0 0 1 0 0 0 ⎝ ⎭ ⎝ ⎭⎛ 1 0⎫⎪ 00 1⎪ →1 0⎪ ⎛ 1 0 1 1 0 0⎫ 0 1 -2 -2 1 0 0 1⎪ → 1 0⎪⎪ ⎪ 0 0⎪ 0 0⎪ 00⎪ 0 0⎪ ⎝⎭ ⎝ ⎭ 0 1 1 0 1 -2 -2 1 0 2 7 -2 0 0 0 0 1 0 0 0 0 1 0 00 2 7 -2 0 -1 0 0 1 0 0 0 0 1 0 0- - ⎪ ⎝ ⎭ ⎝ ⎭1 ⎛ 1 0 0 1 0 0⎫⎪0 1 0 2 1 0 ⎛ 10 0 1 0 0⎪⎫ 0 1 0 2 1 0 0 0 17 -2 1⎪0 0 2 7 -2 1⎪1 ⎪→ ⎪ → 10 - 0 0 0⎪ .1 0 -1 0 0 0⎪2⎪ 0 1 2 0 0 0⎪ 00 10 01 0 0 0⎪0 0 1 0 0 0⎪⎝所以;⎭⎪⎝2⎭⎛1 0 1 ⎫ ⎛ 5 1 ⎫- 2 ⎪⎛ 1 0 0⎫ - 2 -1 - 2 ⎪ A -1 = 0 1 1 ⎪ -2 1 0⎪ = 5 -1 1 ⎪ . ⎪ ⎪ ⎪ 1 ⎪ 7 -2 17 1 ⎪ 0 0 2 ⎪ ⎝ ⎭ 2 -1 2 ⎪ 此方法在计算上并不简单,但是它把平常的单纯的一种变换变成了两种变换同时应用,把已知的可逆矩阵置于含单位矩阵的分块矩阵中,以此求逆矩阵, 有时比较简单.3.4 矩阵秩基本不等式矩阵理论中, 矩阵的秩是一个重要的概念,而矩阵经过运算后所得新矩阵 的秩往往与原矩阵的秩有一定关系. 现把高等代数书中有关矩阵秩最基本的不等式总结如下:(1)矩阵和的秩不超过两矩阵秩的和.即:设 A 、 B 均为m ⨯ n 矩阵,则:r ( A + B ) ≤ r(A ) + r ( B ) .(2)矩阵乘积的秩不超过各因子的秩.即:设 A 是m ⨯ n 矩阵 , B 是n ⨯ s 矩(3)r ⎛A B ⎫阵,则:r ( A B ) ≤≥ r ( A ) + r ( B ) . m i n {r ( A ) , r ( B )}.(4)r ⎝ 0 C ⎭ ⎪ ⎛A ⎫ ⎪⎪ ≥ A i j .A ⎪ ⎝ m ⎭再来介绍由分块矩阵证明导出的两个基本不等式例 3.4.1[11] :(薛尔弗斯特不等式)设A = (a ) ,B = (b ) ,证明:ij s ⨯nij n ⨯mr a n k ( A B ) ≥ r a n k ( A ) + r a n k ( B ) - n⎪ 证明:由分块矩阵的乘积⎛ E n 0⎪ ⎫ ⎛E B ⎫ ⎪⎛E n -B ⎫⎛E n 0 ⎫ -A E A n0 0 E ⎪ = ⎪0 - ⎝ s ⎭ ⎝ ⎭ ⎝ 知:m ⎭⎝ A B ⎭ r a n k⎛E n B⎫ = r a n k (E ) + r a n k ( -A B ) = n + r a n k ( A B )A 0 ⎪n.⎝ ⎭但,r a n k⎛E nB ⎫ A 0⎪= r a n k⎛B E n ⎫ ≥ r a n k ( A ) + r a n k ( B ) ⎪故:得证.⎝⎭ ⎝ 0 A ⎭.n + r a n k ( A B )≥ r a n k ( A ) + r a n k ( B )备注:在矩阵秩不等式的证明过程中,我们往往会构造如下的分块矩阵: (1) 矩阵不等式中含两个不同矩阵:构造 ⎛A 0 ⎫⎪;⎝ 0 B ⎭(2) 矩阵不等式中含有两个不同矩阵及阶数:构造⎛ A E ⎫ ⎪ 或者 ⎛ A 0 ⎫ ⎪.⎝ 0 B ⎭ ⎝E B ⎭具体分块矩阵的元素则要看题目所给的条件.例 3.4.2 [6] :(Frobenius 不等式)设 A 、 B 、C 是任意3 个矩阵,乘积ABC 有意义,证明:r ( A B C ) ≥ r ( A B ) + r ( B C ) - r ( B )证明:设 B 是n ⨯ m 矩阵,r ( B ) = r那么存在n 阶可逆阵 P , m 阶可逆阵Q ,使B = ⎛Er0⎫ P ⎪ Q .⎝ 0 0⎭把 P 、Q 适当分块:P = (M , S ),Q =⎛N ⎫, 由上式有: T ⎝ ⎭故:r ( A B C )= r ( A M N C ) B = (M , S )⎛E r0⎫ ⎛N ⎫ = M N .⎪ ⎪ ⎝ 0 0⎭ ⎝T ⎭≥ r ( A M ) + r ( N C ) - r0 ≥ r ( A M N ) + r ( M N C ) - r ( B )得证.= r ( A B ) + r ( B C ) - r ( B ) .3.5 矩阵秩不等式证明的应用矩阵基本不等式的证明思路,在一般不等式中也常常用到, 以下例题是对矩阵秩不等式的推广及其应用:例 3.5.1[11] :设 A 为m ⨯ k 矩阵, B 为k ⨯ n 矩阵,则证明:r a n k ( A )+r ank( B ) - k≤ r ank( AB) ≤ m i n {r a n k ( A ) , r a n k ( B )}证明:先证明右边的不等式,由:(A 0)(E k0 B ) = ( A A B ) ;E n可得:⎛E k A E 0⎪ ⎫ ⎛B ⎪⎫ = ⎛ B A B ⎫⎪ ,⎝m ⎭ ⎝ ⎭⎝ ⎭r a n k ( A ) =r ank( A 0) = r a n k ( A A B ) ≥ r a n k ( A B ) ;r a n k ( B ) = r a n k ⎛ B ⎫ = r a n k ⎛ B ⎫≥ r a n k ( A B ) .⎪ ⎪⎝ 0 ⎭ ⎝AB ⎭ 再证左边的不等式.注意到下列事实:⎛E m -A ⎫ ⎛ A 0 ⎫ ⎛E ⎪k -B ⎫ = ⎛ 0 -A B ⎫⎪ 0 E ⎪E B 0E⎪ E 0 ⎝k ⎭ ⎝ k 则:⎭ ⎝ n ⎭⎝ k ⎭0 ⎫⎛ 0r a n k ⎛ A ⎪ = r a n k-A B ⎫ ⎪于是:⎝E kB ⎭ ⎝E k0 ⎭⎛ A 0 ⎫r a n k ( A ) + r ank ( B ) ≤r ank ⎪ = r a n k ( -A B ) + r a n k (E k )= r a n k ( A B ) + k⎝E kB ⎭ 从而: r a n k ( A ) + r a n k ( B ) - k ≤ r a n k ( A B ) .这里也是用到构造矩阵的方法.例 3.5.2 [6] :设n 阶矩阵 A 、 B 可交换,证明:r a n k ( A + B ) ≤ r a n k ( A ) + r a n k ( B ) - r a n k ( A B )→ → , ⎝ ⎭ 解:利用分块初等变换,有:⎛A O ⎫ ⎛A B ⎫ ⎛A + B B ⎫⎪ ⎪⎪ ⎝O B ⎭ ⎝O B ⎭ ⎝ B B ⎭ 因为 AB = BA ,所以:⎛ E O ⎫ ⎛A + B B ⎫ = ⎛A + B B ⎫ .B -A - ⎪ B ⎪ O- ⎪B B A B ⎝ 于是,有:⎭ ⎝ ⎭ ⎝ ⎭r a n k ( A ) + r a n k ( B )= r a n k⎛A + B B ⎫≥ r a n k ⎛A + B B ⎫B ⎪⎝ B ⎭ ⎝ ⎪O-A B ⎭即:r a n k ( A + B )得证.≥ r a n k ( A + B ) + r a n k ( A B ) .≤ r a n k ( A ) + r a n k ( B ) - r a n k ( A B ) .例 3.5.3:设 A 是n 阶方阵,且r ( A ) = r ( A 2 ,证明:对任意自然数k ,有r ( A k ) = r ( A )⎛A 2O ⎫证:构造分块矩阵 O A 2 ⎪,由 Frobenius 不等式: 2 2 2 ⎛A O ⎫ ⎛A 2 -A 3 ⎫ ⎛O -A 3 ⎫ 3 r ( A )+r( A ) ≤ r ⎪ = r A A 2 A O ⎪ = r A O ⎪ = r ( A ) + r ( A ) . 由:r ( A ) = r ( A 2 ) ⎝ ⎭ ⎝ ⎭ ⎝ ⎭所以,r ( A3 ) = r ( A 2 * A )≤ r ( A2 ) .故: r(A 2 ) = r ( A 3 .由此可推得:r ( A3) = r ( A 4) , r ( A4) = r ( A5 ) , .故:对任意自然数k , 有:r ( A k ) = r ( A ) .3.6 综合应用在掌握了分块矩阵的技巧之后,可以由其导出的一个重要的定理:特征多项式的降阶定理,以下主要讨论该定理及其结论的应用.例 3.6.1 [6] :(特征多项式的降阶定理)设 A 是m ⨯n 矩阵, B 是n ⨯ m 矩阵. 证明: AB 的特征多项式f A B ( ) 与 BA 的特征多项式f B A( ) 有如下的关系:nm1 2 s证:先要把上式改写为:n f () =m f () .A BB AnE -m A B =mEE 1 Bn - B A .用构造法,设 ≠ 0 ,令: H =n.A E m⎛1 ⎫ ⎛E 1 B ⎫对 ⎛E n 0⎪ ⎫ E n B ⎪= n ⎪ ⎝ -A E⎪⎪ 1 ⎪ 两边取行列式得: n ⎭ A E⎝ m ⎭ 0 E - ⎝A B ⎪⎭ H = E -1 A B = 1 m E - A B .⎛E 1 B ⎫ ⎛E nm 0 ⎫⎛ 1( ) m1 B ⎫ 再对 n ⎪ -A E ⎪ E - B A ⎪ 两边取行列式得: ⎪ ⎪ = n⎪⎝ A E m ⎭⎝ n ⎭ ⎝ H = E -0 1B A = E m ⎭ 1 n E - B A .故: 1nE n- B A =1Em mn- A B() nmE n - B A = nE m - A B .上述等式是假设了 ≠ 0 ,但是两边均为的n + m 次多项式,有无穷多个值使它们成立(0)≠ ,从而一定是恒等式,即证.这个等式也称为薛尔弗斯特(Sylvester )公式. 以下例题是定理的应用. 例 3.6.2 [6] :设 A 为m ⨯ n 矩阵, B 为n ⨯ m 矩阵,证明: AB 与 BA 有相同的非零特征值.证:由定理:m E - B A = n E - A B .设 E m- A B = m -s (- ) ( - ) ( - ) ,其中12 m ≠么有:0 ,即 AB 有s 个非零特征值:1, 2, , s , 由上面两式,那nE - B A = ( - 1) ( - ) 2 (- )n- s s即证 BA 也只有s 个非零特征值:1, 2, , s .m∑ 例 3.6.3 [6] :设 A 、 B 分别是m ⨯n 和n ⨯ m 矩阵,证明:t r A B = t r B A .解:由上例知,若E - A B = m -s ( - a ) ( - a )m1s其中a 1a 2 a s ≠ 0. 则 AB 的全部特征值为1 = a 1, , s= a s , s +1= = m = 0 ,且:E - B A = n -s ( - a ) ( - a ) .n1s即 BA 的全部特征值为:1 = a 1,2 = a2, ,s +1= = n = 0 .从而 t r A B =sa ii=1=t r B A .可见,在一些问题中,直接利用特征多项式的降阶定理会更加方便处理,这里则要求我们对分块矩阵的了解更加深刻.结论本文主要通过“分块矩阵、分块矩阵及其应用”两个部分,分别简单介绍了分块矩阵的性质概念、导出的定理结论和相关应用.主要是将分块矩阵的技巧和推广做了一个内容的总结.本文简单的将矩阵工具应用于计算行列式、解决线性方程组、求矩阵的逆、证明矩阵秩的相关定理等,对应不同问题也举了几个重要的应用以及它们的综合应用.将以前出现的矩阵思想整体化,并对相关知识也做了一个系统的复习.最后,本文还有一些不足之处,有待于进一步的改善和提高.参考文献[1] 上海交通大学线性代数编写组. 线性代数[M]. 高等教育出版社, 1982. [2] 北京大学. 高等代数{M}. 高等教育出版社, 1998.[3] 高百俊. 分块矩阵的初等变换及其应用[J]. 伊犁师范学院学报, 2007(4):14-18.[4]张红玉, 魏慧敏. 矩阵的研究[M]. ft 西人民出版社, 2010.[5]雷英果. 分块初等方阵及其应用[J].工科数学, 1998, 14(4):150-154. [6]钱吉林. 高等代数题解精粹(第二版)[M]. 中央民族大学出版社, 2010.[7] 王莲花, 李念伟, 梁志新. 分块矩阵在行列式计算中的应用[J]. 河南教育学院学报(自然科学版), 2005, 14(3):12-15.[8] 张贤科, 许甫华. 高等代数学[M]. 清华大学出版社, 1998:91-96.[9]杨子胥. 高等代数习题集[M]. ft东科学技术出版社, 1981.[10]鲁翠仙. 分块矩阵在求矩阵逆的应用[D]. 云南:云南大学数学系数学研究所,2009:14-15.[11]刘丁酉. 高等代数习题精解[M].中国科学技术大学出版社, 1999.“”“”At the end, Xiao Bian gives you a passage. Minand once said, "people who learn to learn are very happy people.". In every wonderful life, learning is an eternal theme. As a professional clerical and teaching position, I understand the importance of continuous learning, "life is diligent, nothing can be gained", only continuous learning can achieve better self. Only by constantly learning and mastering the latest relevant knowledge, can employees from all walks of life keep up with the pace of enterprise development and innovate to meet the needs of the market. This document is also edited by my studio professionals, there may be errors in the document, if there are errors, please correct, thank you!。
word2007如何制作矩阵
![word2007如何制作矩阵](https://img.taocdn.com/s3/m/4a04e7e0710abb68a98271fe910ef12d2af9a972.png)
word2007如何制作矩阵
如今很多人都会使用word2007制作一些数学题,线性代数、数学分析里,经常要用到矩阵,下面是店铺整理的word2007制作矩阵的方法,供您参考。
word2007制作矩阵的方法
打开word文件。点击插入→对象
选择公式编辑3.0,点确定,进入到公式编辑器界面。
点击矩阵图形那里
选择需要的图框
然后选择矩阵的结构,假设我们做一个3*3的矩阵
即出现这Байду номын сангаас效果。
在矩阵里面输入相应的数字。
点击空白处,即出现在word文档里了。
word2007制作矩阵的
分块矩阵
![分块矩阵](https://img.taocdn.com/s3/m/46a4c6758e9951e79b8927df.png)
⎡A A 11 12 ⎢A A ⎢ 21 22 A= ⎢ ⎢ ⎢Ap1 Ap2 ⎣
At ⎤ ⎡B B 1 11 12 ⎥ ⎢B B At ⎥ 2 ⎢ 21 22 , B= ⎥ ⎢ ⎥ ⎢ Apt ⎥ ⎢Bt1 Bt2 ⎣ ⎦
Br ⎤ ⎥ B2r ⎥ ⎥ ⎥ Btr ⎥ ⎦
Aij 为 mi×sj 子 块 ( i=1 , 2 , … , p ; j=1 , 2,…,t),且m1 + m2 + … +mp = m,s1 + s2 + … + st= s; Bjk 为 sj×nk 子 块 ( j=1 , 2 , … , t ; k=1 , 2,…,r),且k1 + k2 + … + kr = n,s1 + s2 + … + st= s。 令AB = C,则
由Laplace定理知,分块对角矩阵的行列式 具有下述性质:|A|=|A1|·|A2|·…·|As| 由此性质可知,若|Ai| ≠ 0(i=1,2,…, s),则|A| ≠ 0,并有
⎡A ⎢ −1 ⎢ 0 A = ⎢ ⎢ ⎢ 0 ⎣
−1 1
0 A
−1 2
0
0 ⎤ ⎥ 0 ⎥ ⎥ −1 ⎥ As ⎥ ⎦
⎡ C11 C12 ⎢C C22 21 AB = C = ⎢ ⎢ ⎢ ⎢C p1 C p 2 ⎣
Cij = Ai1B1j + Ai2B2j + ? + Ait Bt j = ∑AikBkj
k=1 t
C1r ⎤ ⎥ C2 r ⎥ , ⎥ ⎥ C pr ⎥ ⎦
( i =1,2,
, p; j =1,2, , r )
0 0 1 0 0 1 0 0 0 0
第二章分块矩阵(优选)word资料
![第二章分块矩阵(优选)word资料](https://img.taocdn.com/s3/m/11a6bf67a0116c175e0e4863.png)
第二章分块矩阵(优选)word资料分块矩阵1 分块矩阵及其运算对矩阵进行"分块"是处理较高阶矩阵的一种常用技巧,分块矩阵的运算能使矩阵间的一些关系更清楚地反映出来.将n m ⨯矩阵A ()n m kl a ⨯=作如下分块s tm m m n st t t n s n s A A A A A A A A A A }}}2121212221212111⎪⎪⎪⎪⎪⎭⎫⎝⎛==()t s ij A ⨯ 其中ij A 是j i n m ⨯矩阵(s i ,,2,1 =,t j ,,2,1 =),m msi i=∑=1,n n tj j =∑=1.这里A 既是以数kl a 为元素的n m ⨯数字矩阵,又是以矩阵ij A 为元素的t s ⨯分块矩阵.注: 在分块矩阵中,每一行小矩阵有相同的行数,每一列小矩阵有相同的列数,在对一个矩阵进行分块时,一定要注意这一点.不能把分块矩阵简单地理解为“以矩阵为元素的矩阵”.比如: ⎪⎪⎭⎫ ⎝⎛22211211A AA A 其中()()98,7,63,542122211211==⎪⎪⎭⎫ ⎝⎛=⎪⎪⎭⎫ ⎝⎛=A A A A ,不是⎪⎪⎪⎪⎭⎫⎝⎛=987654321A 的分块.分块矩阵的运算数字矩阵的运算可以用分块矩阵进行操作,换句话说,数字矩阵的运算可以转化为分块矩阵的运算.所以分块矩阵的运算不是矩阵的新运算,只不过是数字矩阵运算的一种新的运算方法,其运算的结果与数字矩阵的运算结果是一致的.加法 设()ts ijn m A A ⨯⨯=,()ts ijn m B B ⨯⨯=,其中ij ij B A ,是j i n m ⨯矩阵(s i .,2,1 =;t j ,,2,1 =),m m si i =∑=1,n n tj j =∑=1,则()ts ijij B A B A ⨯+=+.乘法 设()ts ijn m A A ⨯⨯=,()rt jkl n B B ⨯⨯=,其中ij A 是j i n m ⨯矩阵,jk B 是k j l n ⨯矩阵(s i .,2,1 =;t j ,,2,1 =;r k ,,2,1 =),m msi i=∑=1,n n t j j =∑=1,l l rk k =∑=1,则()r s ik C AB ⨯=C =,其中tk it k i k i ik B A B A B A C +++= 2211 (s i .,2,1 =;r k ,,2,1 =).数乘 设()ts ijn m A A ⨯⨯=,则()ts ijkA kA ⨯=.分块矩阵的运算规则可以用一句话概括成“只要运算有意义,分块矩阵的运算可以按照数字矩阵的运算规则进行运算”.这里的“运算有意义”是指运算对矩阵行数、列数的相应要求,它包含三层意思:第一 分块之前的两数字矩阵运算要有意义(相加时两矩阵的行、列数分别相同;相乘时前一矩阵的列数要等于后一矩阵的行数);第二 分块以后的两个分块矩阵的运算要有意义(相加时两分块矩阵的行、列数分别相同;相乘时前一分块矩阵的列数要等于后一分块矩阵的行数);第三 各对子矩阵的运算要有意义(各对相加的子矩阵的行、列数分别相同;各对相乘的子矩阵中前一子矩阵的列数要等于后一子矩阵的行数).要保证运算有意义,只需把握住分块环节即可,具体地说:在用分块方法作矩阵加法时,两数字矩阵行和列的分法应当分别一致;在用分块方法作矩阵乘法时,前一数字矩阵列的分法与后一数字矩阵行的分法应一致.注意:在分块矩阵的乘法中,各子矩阵的先后顺序不能随意颠倒,应与其母矩阵的先后顺序一致.下边证明分块矩阵的上述乘法规则.矩阵()nm ija A ⨯=的位于r i i i ,,,21 行、s j j j ,,,21 列的子块(也称为A 的子矩阵,A叫该子矩阵的母矩阵)记为⎭⎬⎫⎩⎨⎧s r j j j i i i A 2121,以后总假定s r j j j i i i <<<<<< 2121,. 例如:由元素ij a 构成的子块为()⎭⎬⎫⎩⎨⎧=j i A a ij ; 由A 的第i 行元素构成的子块为()⎭⎬⎫⎩⎨⎧=n i A a a a in i i 1221(简记为⎭⎬⎫⎩⎨⎧----i A ); 由A 的第j 列元素构成的子块为⎭⎬⎫⎩⎨⎧=⎪⎪⎪⎪⎪⎭⎫⎝⎛j m A a a a mjjj 1221(简记为⎭⎬⎫⎩⎨⎧----j A ); 由A 的第r i i i ,,,21 行元素构成的子块为⎭⎬⎫⎩⎨⎧n i i i A r 1221(简记为⎭⎬⎫⎩⎨⎧----r i i i A 21);由A 的第s j j j ,,,21 列元素构成的子块为⎭⎬⎫⎩⎨⎧s j j j m A 2112(简记为⎭⎬⎫⎩⎨⎧----s j j j A 21). 由矩阵的乘法可知,AB 的),(j i 元等于A 的第i 行与B 第j 列元素的对应乘积之和,即:))(()(⎭⎬⎫⎩⎨⎧----⎭⎬⎫⎩⎨⎧----=⎭⎬⎫⎩⎨⎧j B i A j i AB .由此易得 ))(()(21212121⎭⎬⎫⎩⎨⎧----⎭⎬⎫⎩⎨⎧----=⎭⎬⎫⎩⎨⎧s r s r j j j B i i i A j j j i i i AB ,即:AB 的位于r i i i ,,,21 行、s j j j ,,,21 列的子块等于A 的第r i i i ,,,21 行元素构成的子块与B 第s j j j ,,,21 列元素构成的子块的乘积. 特别地 B i A B i A i AB )())(()(⎭⎬⎫⎩⎨⎧----=⎭⎬⎫⎩⎨⎧--------⎭⎬⎫⎩⎨⎧----=⎭⎬⎫⎩⎨⎧----,即:AB 的第i 行元素构成的子块等于A 的第i 行元素构成的子块与B 的乘积.)())(()(⎭⎬⎫⎩⎨⎧----=⎭⎬⎫⎩⎨⎧----⎭⎬⎫⎩⎨⎧--------=⎭⎬⎫⎩⎨⎧----j B A j B A j AB ,即:AB 的第j 列元素构成的子块等于A 与B 的第j 列元素构成的子块的乘积.进一步有 ]))[(()(21212121⎭⎬⎫⎩⎨⎧----⎭⎬⎫⎩⎨⎧----=⎭⎬⎫⎩⎨⎧s r s r j j j HK B i i i A j j j i i i HK AB ))()((2121⎭⎬⎫⎩⎨⎧----⎭⎬⎫⎩⎨⎧----=s r j j j K H B i i i A .特别地 ))()(()(⎭⎬⎫⎩⎨⎧----⎭⎬⎫⎩⎨⎧----=⎭⎬⎫⎩⎨⎧j K H B i A j i HK AB ,))(()(HK B i A i HK AB ⎭⎬⎫⎩⎨⎧----=⎭⎬⎫⎩⎨⎧----,))(()(⎭⎬⎫⎩⎨⎧----=⎭⎬⎫⎩⎨⎧----j K H AB j HK AB .下边证明分块矩阵乘法的计算规则:=ABt rtt n n n l tr r r l t l t n st t t n s n s m m m B B B B B B B B B A A A A A A A A A }}}{{{21212121212221212111212221212111⎪⎪⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎪⎪⎭⎫⎝⎛ srm m m l sr r r l s l s C C C C C C C C C }}}2121212221212111⎪⎪⎪⎪⎪⎭⎫ ⎝⎛= C = 其中tk it k i k i ik B A B A B A C +++= 2211 (s i .,2,1 =;r k ,,2,1 =),m m si i =∑=1,n ntj j=∑=1,l l rk k =∑=1.只需证明C 与AB 对应位置的元素相同.设C 的),(j i 元为ij c ,则ij c 必在C 的某个子块th kt h k h k kh B A B A B A C +++= 2211之中,设ij c 是kh C 的),(j i ''元,则 i m m i k '+++=-11 ,j l l j h '+++=-11 ,这里k m i ≤'≤1,h l j ≤'≤1.由于ij c 是th kt h k h k B A B A B A ,,,2211 的),(j i ''元之和.而ph kp B A 的),(j i ''元为kp A 的i '行与ph B 的j '列元素的对应乘积之和:))((⎭⎬⎫⎩⎨⎧'----⎭⎬⎫⎩⎨⎧----'j B i A ph kp .但是,kp A 的第i '行恰为A 的第i 行,kp A 所处的列恰为A 的第p p p p n n n n n n n +++++++++---111111,,2,1 列;ph B 的第j '列恰为B 的第j 列,ph B 所处的行恰为B 的第p p p p n n n n n n n +++++++++---111111,,2,1 行.(参看下图)i i n m m kt p k n st p k kp p k p n p k n k k k pp A A A A A A A A A ''-++----⎪⎪⎪⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛}})1()1()1(1)1(11111}},p p j h l j h n n n pr h p th h p ph h p h l h p l p B B B B B B B B B }}}1111)1()1()1(1)1(1-'-'-⎪⎪⎪⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛++--所以))((⎭⎬⎫⎩⎨⎧'----⎭⎬⎫⎩⎨⎧----'j B i A ph kp ⎪⎪⎪⎪⎪⎭⎫⎝⎛=++++++++++++++++++------j n n n j n n j n n n n n i n n i n n i p p p p p p p p b b b a a a )()2()1()()2()1(111111111111),,,(jn n n n i j n n n n i j n n n n i p p p p p p b a b a b a )()()2()2()1()1(1111111111+++++++++++++++++++=---- , 于是⎥⎦⎤⎢⎣⎡⎭⎬⎫⎩⎨⎧'---⎭⎬⎫⎩⎨⎧---'⎥⎦⎤⎢⎣⎡⎭⎬⎫⎩⎨⎧'---⎭⎬⎫⎩⎨⎧---'⎥⎦⎤⎢⎣⎡⎭⎬⎫⎩⎨⎧'---⎭⎬⎫⎩⎨⎧---'=))(())(())((2211j B i A j B i A j B i A c th kt h k h k ij )()()()()()1()1()()()1()1(1111111121211111j n n n n i j n n n n i j n n n n i j n n i j n in i i t t t t b a b a b a b a b a b a +++++++++++++++++++++++=--nj in j i j i b a b a b a +++= 2211,恰为AB 的),(j i 元,故C AB =.转置 设()ts ijA A ⨯=,则⎪⎪⎪⎪⎪⎭⎫⎝⎛'''''''''='st t ts s A A A A A A A A A A 212221212111对分块对角阵⎪⎪⎪⎪⎪⎭⎫⎝⎛=s A A A A21,其中i A 是方阵,s i ,,2,1 =.有:1、s A A A A 21=;2、A 可逆当且仅当s A A A ,,,21 均可逆,且⎪⎪⎪⎪⎪⎭⎫⎝⎛=----112111s A A A A ; 3、A 的秩等于s A A A ,,,21 的秩之和.(此结论当i A 不是方阵时亦然)利用矩阵分块可以将矩阵乘积中的一些关系反映得更清楚,比如: 设()()l n jk nm ijb B a A ⨯⨯==,.对,A B 作如下分块,A 的每个元素为一块,B 的每一行为一块,可得⎪⎪⎪⎪⎪⎭⎫⎝⎛+++++++++=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎪⎪⎭⎫⎝⎛=n mn m m n n n n n mn m m n n B a B a B a B a B a B a B a B a B a B B B a a aa a aa a a AB22112222121121211121212222111211, 由此可见,AB 的每个行向量是B 的行向量的线性组合.对,A B 作如下分块,A 的每一列为一块,B 的每个元素为一块,可得()⎪⎪⎪⎪⎪⎭⎫⎝⎛=nl n n l l n b b b b b b b b b A A A AB 21222211121121()n nl l l n n n n A b A b A b A b A b A b A b A b A b +++++++++= 221122221121221111由此可见,AB 的每个列向量是A 的列向量的线性组合.当0=AB 时,对,A B 作如下分块,A 整个作为一块,B 的每个列为一块,可得 ()()()()0002121===l l AB AB AB B B B A AB即有0=j AB (l j ,,2,1 =).由此可见,B 的每个列向量是以A 为系数矩阵的齐次线性方程组0=Ax 的解.若对,A B 作下述分块,A 的每一行作为一块,B 整个为一块,可得()⎪⎪⎪⎪⎪⎭⎫⎝⎛=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=0001121 B A B A B A B A A A AB m m即有0=B A i ,亦即0=''i A B (m i ,,2,1 =).由此可见,A 的每个行向量是以B '为系数矩阵的齐次线性方程组0='x B 的解.例1 A 是n m ⨯矩阵,证明:(1) 存在矩阵0≠B 使⇔=0AB 秩n A <)(; (2) 存在矩阵0≠C 使⇔=0CA 秩m A <)(.证明:只证明(1),(2)的证法类似.(1)⇒:0=AB ,B ∴的每个列向量是齐次线性方程组0=Ax 的解.而0≠B ,可知0=Ax 有非零解,所以系数矩阵A 的秩<未知量的个数n ,即:秩n A <)(.⇐:因为秩n A <)(,所以齐次线性方程组0=Ax 有非零解.以0=Ax 的若干非零解为列向量构造矩阵B ,则0≠B 且有分块矩阵的乘法可知0=AB .2 分块矩阵的初等变换类似于数字矩阵,分块矩阵也有其初等变换和初等矩阵,而且它们之间的联系也与数字矩阵中两者之间的联系类似.广义初等变换下述三种变换称为广义初等行(列)变换: 1、对换分块矩阵中两行(列)的位置;2、用一个非退化的矩阵D 左(右)乘分块矩阵某一行(列)中的所有元素; 3、用一个非零矩阵C 左(右)乘分块矩阵的某一行(列)后加于另一行(列).注: ① 行变换左乘,列变换右乘.② 根据所作变换中矩阵运算的需要,对用来左乘或右乘的矩阵应有行、列上的自然要求.比如: 要将分块矩阵j i jt j j it i i m m A A A A A A j i A }}2121⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=的第j 行左乘非零矩阵C 后加于第i 行,为了保证乘法可行,C 的列数就必须等于A 的第j 行子矩阵的行数j m ;而为了保证之后的加法可行,C 的行数就必须等于A 的第i 行子矩阵的行数i m ,所以C 应该是j i m m ⨯矩阵.广义初等矩阵对分块单位矩阵⎪⎪⎪⎪⎪⎭⎫⎝⎛=r n n n n E E E E21 (n nri i=∑=1)施行一次广义初等变换所得到的分块矩阵叫做广义初等矩阵.广义初等矩阵共有六类,它们是:⎪⎪⎭⎫ ⎝⎛j i P :交换分块单位矩阵j i ,两行的位置所得到的广义初等矩阵; ()j iP :交换分块单位矩阵j i ,两列的位置所得到的广义初等矩阵;()i D P )(:用非奇异矩阵D 左乘分块单位矩阵的第i 行所得到的广义初等矩阵; ())(D i P :用非奇异矩阵D 右乘分块单位矩阵的第i 列所得到的广义初等矩阵;⎪⎪⎭⎫ ⎝⎛j C i P )(:将分块单位矩阵的第j 行左乘C 后加于第i 行所得到的广义初等矩阵; ())(C j iP :将分块单位矩阵的第j 列右乘C 后加于第i 列所得到的广义初等矩阵.⎪⎪⎭⎫ ⎝⎛j i P j i E E E E E E E E ji n n n n n n n n r j ij i ji ⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛=+-+-111110()j iP j i E E E E E E E E ji n n n n n n n n r j jj i ii ⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛=+-+-11111()i D P )(i E E DE E in n n n r i i ⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛=+-111())(D i P i E E DE E in n n n r i i ⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛=+-111⎪⎪⎭⎫ ⎝⎛j C i P )(j i E E CE E ji n n n n r j i⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛=1 ())(C j iP j i E E CE E ji n n n n r ji⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛=1 由此可见这六类广义初等矩阵之间的关系:⎪⎪⎭⎫⎝⎛j i P ()j i P '=,()i D P )(=())(D i P ,⎪⎪⎭⎫⎝⎛j C i P )(=())(C i j P .广义初等矩阵的转置:⎪⎪⎭⎫ ⎝⎛'j i P =()j i P ,()j iP '=⎪⎪⎭⎫ ⎝⎛j i P ;()()i D P i D P )()('=',()())()(D i P D i P '=';⎪⎪⎭⎫⎝⎛'=⎪⎪⎭⎫ ⎝⎛'i C j P j C i P )()(,()())()(C i j P C j i P '='.广义初等矩阵的逆:广义初等矩阵的行列式不为零,故可逆,且()j i P j i P =⎪⎪⎭⎫ ⎝⎛-1,()⎪⎪⎭⎫ ⎝⎛=-j i P j iP 1;()()i D P i D P )()(11--=,()())()(11--=D i P D i P ; ⎪⎪⎭⎫ ⎝⎛-=⎪⎪⎭⎫ ⎝⎛-j C i P j C i P )()(1,()())()(1C j i P C j i P -=-.广义初等变换与广义初等矩阵间的关系:对分块矩阵施行一次广义初等行变换相当于在其左边乘上一个相应的广义初等矩阵;对分块矩阵施行一次广义初等列变换相当于在其右边乘上一个相应的广义初等矩阵.注: ① 这里相应的广义初等矩阵是指将所作初等变换作用到分块单位矩阵上所得到的广义初等矩阵.此结果的证明与数字矩阵相应结果的证明相仿.② 对应于n m ⨯数字矩阵A 的分块矩阵s tm m m n st t t n s n s n m A A A A A A A A A A }}}2121212221212111 ⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=⨯,n n m m t j j s i i ==∑∑==11,,相应的分块单位矩阵有两个,它们是⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=s m m m m E E E E 21和⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=t n n n n E E E E 21.对A 施行一次广义初等行变换,等于在A 的左边乘上一个由m E 经过相应的广义初等行变换得到的广义初等矩阵;对A 施行一次广义初等列变换相当于在其右边乘上一个由n E 经过相应的广义初等列变换得到的广义初等矩阵.③ 数字矩阵的初等变换和初等矩阵是广义初等变换和广义初等矩阵的特例,广义初等变换可以用若干次数字矩阵的初等变换来实现,所以在数字矩阵的初等变换下矩阵的不变性质在广义初等变换下也不变.例2 A 是n s ⨯实矩阵,证明:秩-'-)(A A E n 秩s n A A E s -='-)(.证明: 构造矩阵⎪⎪⎭⎫ ⎝⎛'=n s E A A E B ,下边用两种方法化为对角分块矩阵求秩. 用s E 将A '和A 化为0,得⎪⎪⎭⎫ ⎝⎛'-=⎪⎪⎭⎫ ⎝⎛-⎪⎪⎭⎫ ⎝⎛'⎪⎪⎭⎫⎝⎛'-A A E E E A E E A A E E A E n s n s n s n s 0000,所以 秩=)(B 秩=⎪⎪⎭⎫ ⎝⎛'-A A E E n s 00秩+)(s E 秩+='-s A A E n )(秩)(A A E n '-. 用n E 将A '和A 化为0,得⎪⎪⎭⎫ ⎝⎛'-=⎪⎪⎭⎫ ⎝⎛'-⎪⎪⎭⎫ ⎝⎛'⎪⎪⎭⎫ ⎝⎛-n s n s n s n s E A A E E A E E A A E E A E 0000,所以 秩=)(B +n 秩)(A A E s '-.比较上两式即得: 秩-'-)(A A E n 秩s n A A E s -='-)(.例3 C B A ,,是同阶方阵,证明Frobenius 不等式:秩+)(B 秩≥)(ABC 秩+)(AB 秩)(BC .证明:⎪⎪⎭⎫ ⎝⎛-→⎪⎪⎭⎫ ⎝⎛→⎪⎪⎭⎫ ⎝⎛0000AB BC B ABC AB B ABC B, ∴秩=⎪⎪⎭⎫ ⎝⎛ABC B 00秩⎪⎪⎭⎫ ⎝⎛-0AB BC B . 而 秩=⎪⎪⎭⎫ ⎝⎛ABC B 00秩+)(B 秩)(ABC , 秩⎪⎪⎭⎫ ⎝⎛-0AB BC B≥秩+)(AB 秩)(BC -=秩+)(AB 秩)(BC ①. 所以: 秩+)(B 秩≥)(ABC 秩+)(AB 秩)(BC .① 秩≥⎪⎪⎭⎫ ⎝⎛0C B A 秩+)(B 秩)(C 的证明:设 秩1)(r B =,秩2)(r C =,则B 中有一个1r 级子式1B 不为0,C 中有一个2r 级子式1C 不为0.于是⎪⎪⎭⎫⎝⎛0C B A 中由11,C B 所在的行和列确定的21r r +级子式不为0,所以秩=+≥⎪⎪⎭⎫ ⎝⎛210r r C B A 秩+)(B 秩)(C . 对于形如⎪⎪⎭⎫ ⎝⎛⎪⎪⎭⎫ ⎝⎛⎪⎪⎭⎫ ⎝⎛C A B C A B A C B 0,0,0的分块矩阵也有类似的结果.一般地,对三角形分块矩阵此结论亦然.比如:⎪⎪⎪⎪⎪⎭⎫ ⎝⎛rr r r A A A A A A21222111的秩∑=≥r i ii A 1)秩((三角形分块矩阵的秩大于等于其主对角线上块的秩之和).3 标准单位向量n 维列向量j e j ⎪⎪⎪⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=00100 行 n j ,,2,1 =称为(第j 个)n 维标准单位向量.基本性质:设()n m ij a A ⨯=的行向量组为m ααα,,,21 ,列向量组为n βββ,,,21 .1、⎩⎨⎧≠==='j i j i e e ij j i 01δ; 2、i i A e α=',j j Ae β=,ij j i a Ae e =';由此,再利用分块矩阵的乘法规则可得:⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛'''r r i i i i i i A e e e ααα 2121,()()t t j j j j j j e e e A βββ 2121= (1) ()⎭⎬⎫⎩⎨⎧=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛'''t r j j j i i i j j j i i i A e e e A e e e t r 21212121 (2)上述基本性质表明,可以利用标准单位向量将矩阵A 的任意一行、任意一列、任意一个元素、任意一个子矩阵用A 与标准单位向量的乘积表达出来,这对很多问题的讨论将带来帮助.通过以下几个例子能体会到这一点.例4 设A 是n m ⨯矩阵,证明:若对任意n 维列向量x ,恒有0=Ax ,则0=A . 证明 由已知,对标准单位向量n e e e ,,,21 ,有0=j Ae n j ,,2,1 =,而j j Ae β=是A 的第j 个列向量,即A 的每个列向量均为0,故0=A .例5 若对任意n 维列向量x ,恒有0='Ax x ,则A 是反对称矩阵.证明 取i e x =,则ii i i a Ae e ='=0 n i ,,2,1 =.再取j i e e x +=(j i ≠),则有 ji ij i j j i j i j i a a Ae e Ae e e e A e e +='+'=+'+=)()(0,所以ji ij a a -=.总之A 是反对称矩阵.例6 证明:与任意n 阶可逆矩阵可交换的矩阵必是数量矩阵.证明 设()n n ija A ⨯=与任意可逆矩阵可交换,其行向量组和列向量组分别是nαα,,1 和n ββ,,1 .取()⎪⎪⎪⎪⎪⎭⎫ ⎝⎛'''==⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=n n e e e ne e e n B 2121221 则BA AB =,即 ()⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=n n n n αααβββ 212122 比较等号两端矩阵的),(j i 元可得ij ij ja ia =,所以当j i ≠时,0=ij a ,即A 的主对角线以外的元素均为0,故A 为对角矩阵⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=nn a a a A 2211()⎪⎪⎪⎪⎪⎭⎫ ⎝⎛'''==n nn n nn e a e a e a e a e a e a 222111222111. 取 ()n j i j i j i e e e e e e e e P 11111+-+-=⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛''''''''=+-+-n j i j i j i e e e e e e e e 11111),(j i P = 则PA AP =,由(1)式即得()n nn j j j i ii j j j i i i j jj i i i e a e a e a e a e a e a e a e a 111111*********+++---+++---⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛''''''''=+++---+++---n nn j j j i ii j j j i i i j jj i i i e a e a e a e a e a e a e a e a 111111********* 比较等号两端矩阵的),(j i 元可得jj ii a a =,即nn a a a === 2211,所以E a A 11=是数量矩阵.例7 设⎪⎪⎭⎫ ⎝⎛=-⨯0001n n n E A ,证明:⎪⎪⎭⎫ ⎝⎛=-000k n k E A 1,,2,1-=n k ,0=n A . 证明 ()032n e e e A =⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛'''=-1210n e e e ,对k (11-≤≤n k )作归纳: 1=k 时,结论自然成立.假设m k =时结论成立(11-<≤n m ),即⎪⎪⎭⎫ ⎝⎛=-000mn m E A ()0021 n m m e e e ++=. 考虑1+m 的情况,有=+1m A()()003232n m m m n m e A e A e A e e e A =.由性质2,j m e A 是m A 的第j 个列向量,所以()⎪⎪⎭⎫ ⎝⎛==+-+++00000)1(321m n n m m m E e e e A . 又:⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛'=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛'=⎪⎪⎭⎫ ⎝⎛==-00000000001111 A e A e A E A A A n n (A e 1'是A 的第一个行向量故为0).例8 ⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎭⎫ ⎝⎛=-00011000010000100101 n E A 叫n 阶基础循环矩阵.证明:⎪⎪⎭⎫ ⎝⎛=-00k k n k E E A ,1,,2,1-=n k ;n n E A =.证明 ()121-=n n e e e e A ,由本节性质2,有n e Ae =1,1-=j j e Ae n j ,,3,2 =. (*) 下边对k (11-≤≤n k )作归纳:1=k 时,结论自然成立.假设结论对1-k 成立(111-<-≤n k ),即()121321)1(100+-+-+-----=⎪⎪⎭⎫ ⎝⎛=k n n k n k n k k n k e e e e e e E E A . 由分块矩阵的乘法及(*),()121321+-+-+--==k n n k n k n k k e e e e e e A AA A()⎪⎪⎭⎫ ⎝⎛==---+-+-0021121k k n k n n n k n k n E E e e e e e e e . 由(*)()()n n n n n n n E e e e e e e e e A E A AA A ===⎪⎪⎭⎫ ⎝⎛==---12113211010 . 上边例7、例8的结果是应该熟悉的结果.例9 形如⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=---14322154312211321c c c c c c c c c c c c c c c c c c c c C n n n nn n 的矩阵叫循环矩阵.证明:两个循环矩阵的乘积仍是循环矩阵.证明 C 是循环矩阵当且仅当⎪⎪⎭⎫ ⎝⎛+⎪⎪⎭⎫ ⎝⎛++⎪⎪⎭⎫ ⎝⎛+⎪⎪⎭⎫ ⎝⎛+=-----01000000101221223121n n n n n n n E c E E c E E c E c E c C . 利用上题结果,此式可写为11211213121010010010010-------⎪⎪⎭⎫ ⎝⎛+⎪⎪⎭⎫ ⎝⎛++⎪⎪⎭⎫ ⎝⎛+⎪⎪⎭⎫ ⎝⎛+=n n n n n n n n n E c E c E c E c E c C .如果D 是另一个循环矩阵,由上可知D 亦可表为 11211213121010010010010-------⎪⎪⎭⎫ ⎝⎛+⎪⎪⎭⎫ ⎝⎛++⎪⎪⎭⎫ ⎝⎛+⎪⎪⎭⎫ ⎝⎛+=n n n n n n n n n E d E d E d E d E d D .由上题结果,CD 形如 11211213121010010010001-------⎪⎪⎭⎫ ⎝⎛+⎪⎪⎭⎫ ⎝⎛++⎪⎪⎭⎫ ⎝⎛⎪⎪⎭⎫ ⎝⎛+=n n n n n n n n n E b E b E b E b E b CD ,可见CD 是循环矩阵.例10 形如 =F ⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛------121100010001000a a a a n n n 的矩阵叫Frobenius 矩阵.(1) 求k F ,这里n k ≤≤2;(2) 证明:对满足n s <≤1的任意正整数s ,以及任意1+s 个不全为0的数s b b b ,,,10 ,均有00111≠++++--E b F b F b F b s s s s .解 (1) ()αn e e e F 32= ,其中()'---=-11a a a n n α.由基本性质21-=i i Fe e n i ,,3,2 =,n Fe =α (1) 又: ,)(,)(,4312133211221e Fe e F F e F e Fe Fe F e F e Fe =======,如此可得11e F e i i -= n i ,,3,2 =. (2) 设()n k F βββ 21=,则 111)1()()(---====j j k j k j k j Fe e F F Fe F e F βn j ,,3,2 =. (3) 又:当n k =时,αβ)1()2(1111)(====-n n n Fe e FF e F , (4)当1-≤n k 时,1)2(11+==k k e e F β. (5) 于是,当k 给定后,先由(4)或(5)求得1β,再由(3)逐个求得n βββ,,,32 ,从而得到()n k F βββ 21=.(2) 令E b F b F b A s s 01+++= ,下边证明A 有一个列不为0. 101111111Ee b Fe b e F b e F b Ae s s s s ++++=-- ()00010102111)2(≠'=++++=-+ s s s s s b b b e b e b e b e b即A 的第一列01≠Ae ,故0≠A .标准单位向量的定义及结论虽然简单,但在矩阵问题的处理中,却是不可忽视的基本功.(a电枢总成;(b换向器结构电枢轴;2—电枢绕组;3—铁心;4—换向器; 5—换向片;6—轴套;7—压环;8—焊线凸缘S5系列可编序控制器采用模板化结构,由适于各个特殊问题的功能单元组成,图2-1给出S5-115U系统配置。
分 块 矩 阵
![分 块 矩 阵](https://img.taocdn.com/s3/m/62a0e831854769eae009581b6bd97f192279bfe0.png)
,
B
3 0
4 1
1 3
B11 B21
1 0 1
B12 B22
.
所以,
AB
A11 A21
A12 B11
A22
B21
B12 B22
A11B11 A21 B11
A12 B21 A22 B21
A11B12 A21 B12
A12 B22 A22 B22
.
1.1 分块矩阵的运算
1 1
2 3
1,| A3 | 5 ,都不为零,均可逆,故 A 可逆。
又因为
A11
1 3
,
A21
3
1
2 1
,
A31
1 5
,则
.
1.1 分块矩阵的运算
5.例题
由于
A1
A2
A2
A1
A2
A12
A22
,
A3
A3
A32
且
A12
9 ,A22
1 1
2
2
3
3
4
8 11
,A32
C
O
时,有
A O
O 1 A1
B
O
O B 1
.
线性代数
A21
A22
As1 As2
A1r
A2r
Asr
B11 B12
B
B21
B22
Bs1 Bs2
B1r
B2r
Bsr
则 ,
A11 B11
A
B
A21
B21
As1 Bs1
A12 B12 A22 B22
As2 Bs2
A1r B1r
第四节分块矩阵
![第四节分块矩阵](https://img.taocdn.com/s3/m/1e861d1cc5da50e2524d7fd7.png)
0 0 A = 4 1 b
二、分块矩阵的运算规则 加法: 加法:利用分块法对两个同型矩阵进行加法运 算时,两个矩阵必须采用相同的分块法。 算时,两个矩阵必须采用相同的分块法。
a 0 例1 设 A = 0 0 1 a 0 0 0 0 b 1 0 0 , 1 b
数乘: 数乘:数字 k 乘以分块矩阵 A等于 k 乘以 A 的 每一子块。 每一子块。
A11 L A1r kA11 L kA1r M , 则 kA = M 设 A= M M . A L A kA L kA sr sr s1 s1
1 a 0 1
0 0 b 1
0 0 A1 = I 1 b
O A2
单位块
0 0 O= 0 0
a 1 其 中 A = 1 0 a
1 0 I = 0 1
b 1 A = 2 1 b
Ar 1
M , Asr
T A11 L 则 AT = M T A L 1r
M , T Asr
T A1 s1 s
方阵, 分块对角矩阵: 分块对角矩阵: 设A为n阶方阵,若A的分块 矩阵只有主对角线上有非零子块( 矩阵只有主对角线上有非零子块(这些非零子 块必须为方阵),其余子块全为零 方阵),其余子块全为零, 块必须为方阵),其余子块全为零,那么方阵 A就称为分块对角阵. 就称为分块对角阵 分块对角阵.
a 1 B = 0 0 0 a 0 0 0 0 b 1 0 0 0 b
求 A+ B .
解
将 A, B分块
1 a 0 0
0 a 0 0
a 0 A= 0 0
3_4分块矩阵
![3_4分块矩阵](https://img.taocdn.com/s3/m/2d92bc106edb6f1aff001fa7.png)
T
AT A D = T C B D
T
C T B
T
5. 分块对角矩阵
设 A为 n阶矩阵 , 若 A的分块矩阵只有在主对角线上有非零 子块 , 其余子块都为零矩阵 , 且主对角线上子块都是方阵 .即
a 0 A= 0 0
a 0 A1 = 0 0 A1 0 = 0 A , 其中 1 2 b A2 = b 1
1 , a 1 ; b 0 , a 0 ; b
a 1 B= 0 0
0 a 0 0
分块矩阵的结构显得简单且清晰。 分块矩阵的结构显得简单且清晰。
1 0 例1 对于矩阵 A = −1 1
0 0 0 1 0 0 2 1 0 1 0 1
1 0 E2 = 0 1 0 0 0 = 0 0 −1 2 A1 = 1 1
0 0 0 0 1 0 0 1
解:
分块如下: 把A,B分块如下 分块如下
1 0 A= 0 0
E2 0
3 1 2 4 E2 = 0 −1 0 0 0 0 −1 0 1
A12 B11 +B − E2 21 0 E2
2 2 2 1 = 6 3 0 −2
E2 − B11 A− B = 0 − B21
A12 − 0 −2 E 2
3 0 −2 1 −2 1 2 4 = − 6 −3 −2 0 2 0 −2 0
2.数乘分块矩阵 2.数乘分块矩阵
A11 A= M A r1
L L
A1 s B1 1 M , B = M B Ars r1
第3节 分块矩阵(全)
![第3节 分块矩阵(全)](https://img.taocdn.com/s3/m/101511c50c22590102029d19.png)
A12 A22
这是2阶 方阵吗?
a11 a12 A11 a21 a22
为A的子块
A21 a31 a32
A22 a33 a34
矩阵形式上成为以这些子块为元素的分块矩阵。
将一个矩阵分成分块矩阵的方法很多,分块时要注意矩阵的 特点。
例
1 0 A 0 0 0 0 0 1 2 1 0 2 3 E3 0 1 3 4 O 23 0 0 2 0 0 0 0 2
对于分块对角矩阵,可求得
A A1 A2 As
由此可知 A 0 的充分必要条件是 Ai 0 ( i=1,2,…, s)。从而可知分块对角矩阵A可逆的充分必要条件是 Ai(i=1, 2,…,s)均可逆。并且,当A可逆时,有 A11 1 A2 1 A 1 As
T A11 T T A12 A T A 1t
A12 A22 As 2
A1t A2 t Ast
T A21 T A22
则
T A2 t
AsT1 AsT2 T Ast
分块矩阵A的转置,不仅要把分块矩阵A的每一行变为同序 号的列,还要把A的每一个子块 Aij 取转置。
1 0 4 2
0 0 E 0 A1 1
0 1 B11 1 B21 0
O E
1 0 1 2 B 1 0 1 1
E B22
则
E AB A1 O B11 E B21 E B11 B22 A1 B11 B21
例2
设5阶方阵
高等代数第七节分块矩阵
![高等代数第七节分块矩阵](https://img.taocdn.com/s3/m/4855ef26a36925c52cc58bd63186bceb18e8ed14.png)
0 1 b
,
a 0 0 0
B
1 0
a 0
0 b
0 0
0 0 1 b
求 A B, ABA.
解 将 A, B分块
a 1 0 0
A
0 0
0
a 0 0
0 b 1
0 1 b
A1 0
0 ,
A2
a
A1
0
其中
b
A2
1
1 , a 1 ; b
a
B
1 0
0
0 a 0 0
0 0 b 1
0
0 0
*
(
)
A A* 0
(A )
;
0 B B *
B B* 0
(B )
0 A A*
A B* 0
(C )
0 B A*
B A* 0
(D )
0 A B *
分析:根据伴随矩阵公式CC* C E;C* C C 1,由已知分别求C
与C 1即可
2.应用于矩阵的一些运算
解 :C * C C 1, C A 0
Q1B
B1 B2
,B1是u n矩阵,B2是k
u n矩阵,
则
R AB R PAQ Q1B
R
Iu O
O O
B1 B2
R
B1 O
R
B1
3.矩阵秩的不等式证明
另一方面,
RB R
Q1B
R
B1 B2
由秩的不等式性质:
R AB R B1 ≥R B R B2
A2
B2
b 1
1 b b 1
0 2b b 2
1 , 2b
1.4 分块矩阵
![1.4 分块矩阵](https://img.taocdn.com/s3/m/a786a943fe4733687e21aaa0.png)
T A b A
A ( b T A1 )
而 | P || A | 0
( 2) | PQ | A 2 (b T A1 )
Q A (b T A1 )
Q可逆
| Q | 0
T A1 b.
例1 设 A, B 为三阶方阵,且 A (1 , 2 , 3 ), A 1, B (1 , 2 31 ,2 3 1 ), 求 B , 2 A B。
解
B 1 2 31
2 3 1
1 2 2 3 21 2 3 2 2 A B 21 2 2 2 3 1 2 31
例2
求AB. 解
1 0 1 3 0 1 2 1 A 0 0 2 0 0 0 0 2
E2 A O A1 2 E2
1 1 B 0 0
B1 B O
3 0 4 1 0 0 1 2
B2 E2
1.加法与数乘
A11 B11 A1t B1t A B As1 Bs1 Ast Bst
A与B必须用完全相同的分法 。 kA11 kA1t kA kAs1 kAst
3.转置
A11 A1t A As1 Ast
T T A11 As1 AT T T A1t Ast
4 分块对角矩阵 设A为n阶方阵,若A的分块矩阵只有在主对角 线上有非零子块(这些非零子块必须为方阵),其余 子块全为零,那么方阵A就称为分块对角阵. A1 A2 , Ai i 1,2, s 都是方阵. 即如 A A s 1 0 0 0 0 0 1 2 0 0 2 0 0 0 0 3 0 0 都是分块对角阵. 0 1 3 0 0 0 0 0 1 0 0 0 2 1 0 0 2 0 0 0 0 1 5
4分块矩阵
![4分块矩阵](https://img.taocdn.com/s3/m/8e4ee62158fb770bf78a5500.png)
其中
βj ( j = 1, 2, … ,n ).
是 A 的第 j 列.
对于线性方程组
a11x1 a12 x2 a1n xn b1 a x a x a x b 21 1 22 2 2n n 2 am1 x1 am 2 x2 amn xn bm
设 A 和 B 均为 m×n矩阵,分法下:
A11 A1r B11 B1r A B A A B B sr sr s1 s1
其中Aij与Bij的行数相同, 列数 相同, 那末
A11 B11 A B A B s1 s1
4.分ห้องสมุดไป่ตู้矩阵的转置
设分块矩阵
则
A11 A1r A A A sr s1
AT11 AT s1 AT . T T A 1r A sr
5.分块对角矩阵(准 对角矩阵).
设
A1 A
反证法
作业
1.利用逆矩阵解线性方程组:
x1 2 x 2 3 x3 1 2 x1 2 x 2 5 x3 2 3x 5x x 3 1 2 3
2.设
3 4 4 3 A 0 0 0 0
0 0 0 0 8 4 求 A 及A . 2 0 2 2
其运算律与数乘矩阵相同.
3.分块矩阵的乘法. 设A为 m×l 矩阵,B为 l×n矩阵,分块成
A11 A A i1 A s1
A12 Ai2 As2
A1t B11 Ait B Bi1 Ast Bt1
第11节分块矩阵
![第11节分块矩阵](https://img.taocdn.com/s3/m/d8da6f04b52acfc789ebc927.png)
a 0 A = 例1 1 0
1 0 0 B1 a 0 0 = B2 , 0 b 1 B 1 1 b 3
1 0 0
即
a 0 A= 0 0
a 0 0 1 1 b 1 1 b
B1 = B2 B3
1
结论: 结论: 1) 若 A是反对称矩阵 2) 若 A是反对称矩阵 3) 若 A是反对称矩阵 例如
0 1 是反对称矩阵, A= 是反对称矩阵, −1 0 0 3 3A = −3 0
AT = − A
4) 若 A, B都是反对称矩阵 例如
0 1 A= , −1 0
a 0 A = 例4 1 0
1 0 0 a 0 0 = ( A1 , A2 , A3 , A4 ) 0 b 1 1 1 b
说明: 说明: 给了一个矩阵, 给了一个矩阵,可以根据需要把它写成 不同的分块矩阵. 不同的分块矩阵.
其中
a 1 0 0 0 a 0 0 A1 = , A2 = , A3 = , A4 = 1 0 b 1 0 1 1 b
例1 对任一方阵 A 都有 A + AT是对称矩阵, 是对称矩阵,A − AT 是反对称矩阵. 因为 ( A + A ) = A + ( A ) = AT + A = A + AT 证:
T T T T T
例2 证明奇数阶反对称矩阵的行列式等于0. 证明奇数阶反对称矩阵的行列式等于0. 设 A 是 n 阶反对称矩阵, 证: 阶反对称矩阵,则 AT = − A 所以 A = AT = − A = ( −1)n A 当 n奇数时, 奇数时, A = − A 即 A =0 证毕. 证毕.
线性代数分块矩阵
![线性代数分块矩阵](https://img.taocdn.com/s3/m/c7ae152a854769eae009581b6bd97f192279bff5.png)
A1
O
A2
A
,
O
As
其中 Ai i 1, 2,
s 都是方阵, 那末称 A为分块对角矩阵.
6
1
0
⋯
0
0
2
⋯
0
⋯
⋯
⋯
⋯
0
0
⋯
0
A1 B1
0 A2 B2
0
0
1
0
⋯
0
0
2
⋯
0
⋯
⋯
⋯
⋯
0
0
.
1 2 4 1 3 3
A1 B22
,
1 1 2 0 3 1
于是
B11
AB
A1 B11 B21
1
E 1
2
A1 B22
1
0 1 0
2 0 1
.
4 3 3
1 3 1
转
设
1
0
A
1
1
解
0 0 0
1 0 0
,
2 1 0
1 0 1
0
1
1 2
B
1
0
1 1
1 0
0 1
,
4 1
2 0
求 AB .
把A, B分块成
分 块 矩 阵
![分 块 矩 阵](https://img.taocdn.com/s3/m/f4626390dc3383c4bb4cf7ec4afe04a1b071b06e.png)
Ar1
Ar2
A1s
A2
s
Ars
2. 分块矩阵的加法
将m×n 矩阵A 与B 按相同的分块法分别分成r×s的分块矩阵
A11 A12
A
A21
A22
Ar1
Ar 2
A1s
B11 B12
A2 s
,
B
B21B22
Ars
Br1
Br 2
B1s
B2s
Brs
则
A11 B11 A12 B12
3 1
4
0
0
1
在利用分块矩阵的乘法讨论AB 时,下面的特殊情形值得注意。 设A 为m ×l 矩阵,B 为l×n 矩阵,将右矩阵B 按列分块:
B= b11 b12 bn
则
AB= Ab11 Ab12 Abn
若AB=O,则 Ab11 Ab12 Abn O (OO O) ,从而
线性代数
分块矩阵
1
2
3
分块矩阵 的概念
分块矩阵 的运算
分块对角矩阵
1.1 分块矩阵的概念
定义1
用若干条横线与若干条纵线将矩阵分成若干小块,每个小块 称为矩阵的子块;以子块为元素的形式上的矩阵称为分块矩阵。
a b 0 0
例如
A
c
d
0
0
0 0 p q
0
0
r
s
按下述分法分块
a b 0 0
A
Abj O( j 1,2, n)
即 bj ( j 1, 2, n) 是矩阵方程 Aml Xl1 Om1 的解,也就是说 B 的列是 Aml Xl1 Om1 的解。
4. 分块矩阵的转置
将m×n 矩阵A 分成r×s的分块矩阵
2-4-1 分块矩阵
![2-4-1 分块矩阵](https://img.taocdn.com/s3/m/cd768c3351e79b89680226af.png)
Bs1 Bs2 Bsr
其中各对应的子块Aij 与 Bij 有相同的行数和列
数,则 A11 B11
A
B
A21 B21
A12 B12
A22 B22
A1r B1r A2r B2r
As1 Bs1 As2 Bs2 Asr Bsr
(2) 数乘运算
设 为任意的数, 则
所以
1 2 1 0 4 0 0 1
AB
4
1
0
1
5 2 1 2
4
5
2
0
例2 设 求 A1
3 0 0 0 0
0 0 1 0 0
A
0
2
5
0
0
0 0 0 1 0
0
0
0
0
1
解 : 将A分块如下
0 2 0
1 5 0
0 0 1
0 0 0
A1
A2
E2
运算来确定。
二、分块矩阵的运算
矩阵分块后, 把小矩阵当作元素,按普通的 矩阵运算法则进行运算。
(1) 加法
设 A, B 是两个m n矩阵,且用相同的分块 法,得分块矩阵为
A11
A
A21
A12
A22
A1r A2r
B11
B
B21
B12
B22
B1r B2r
As1 As2 Asr
k 1
(4) 转置
设A得分块为 A11
A
A21
A12
A22
A1r A2r
则
As1 As2 Asr
A1T1 A2T1 AsT1
AT
A1T2
第2.5节分块矩阵
![第2.5节分块矩阵](https://img.taocdn.com/s3/m/59dc2d262f60ddccda38a096.png)
都是分块对角阵.
注
(1) A A1 A2 As ; P.21 例 1.20
(3) (3)若 A
A11 1 O A2 1 ( 2 ) 若 每 个 Ai 0, 则 A可 逆 , 且 A . O 1 A1 B1 As
0 0 0 7 5 0 6 , 5E 2 0 5 , 8
0 0 0 , 0 5
3 0 B 0 0 0 3 3E 3 0 0
0 3 0 0
0 0 3 0
0 0 0 1
0 0 9 0 0 1 3 3 0 , B1 9 4 0 3
例 按行、列分块时的乘法运算 :
O 3 A1 5 E2 B1 O O 5 B1
A
按行
分块之后,可以利 用零矩阵和单位 矩阵使运算简化.
a11 a 21 am 1 b11 b21 bs 1
a12 a 22 am 2 b12 b22 bs 2
第2.5节
1. 矩阵的分块 定义
分块矩阵
将 A用若干条横线和纵线分成许多个小矩阵( A的子 块),以子块为元素的形式上的矩阵称为分块矩阵 . 与利将 讨用矩 a11 a12 a13 a14 A1 A2 论已阵 如 A a 21 a 22 a 23 a 24 知适 A A 3 4 a 的当 31 a 32 a33 a 34 性地 质分 a11 a1 2 a13 a14 A1 简块 A2 a 21 a 22 a 23 a 24 化有 运可 A3 a 31 a 32 , A4 a 3 3 a 3 4 算能
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
A12 A22
X X
1 3
X X
2 4
A11
X1 A12 A22 X 3
X
3
A11
X 2 A12 A22 X 4
X
4
E 0
0 E
A11X1 A12 X 3 E
得到4个矩阵方程组
A11
X
2
A12 X 4
0
A22
X
3
0
A22 X 4 E
求解该方程组,得
X 4 A221 X3 0 X1 A1T1 X 2 A111 A12 A221
(2) (解略,请仿(1)方法自行求解)
2. 分块对角矩阵
设A1, A2, … , As均为方阵(不一定同阶),则称下 面的A为分块对角矩阵
A1
A
A2
As
如果矩阵A1, A2, … , As均可逆,则分块对角矩阵A 可逆,且其逆矩阵为
A11
A1
A21
As1
说明:分块对角阵的逆矩阵,与对角矩阵的逆矩
0
0
1 b
A1
A2
A3
a10
A4
,其中A2431
a0 01b
1b0
说明 (1). 矩阵分块时,同一个矩阵可以有不同的 分块方法,应根据需要进行选择。
2、矩阵分块一般形式
矩阵A = ( aij )m×n,在行方向分s块,列方向分t块, 称A为s×t分块矩阵,第k行l列子块Akl是mk×nl阶矩阵。
AtT
二、一些特殊的分块矩阵
1. 2阶分块上(下)三角形矩阵求逆
例2. 求下列2阶分块逆矩阵
(1) A A11
A12 A22
其中A11, A22可逆矩阵
(2) B B21
B12 B22
其中B12, B21可逆矩阵
解(1) :设A的分块逆矩阵为
A1
X1 X3
X X
2 4
AA1 E
AA1 A11
其中B = ( b1, b2, …, bt ), X = ( x1, x2, …, xt ) 例3. 求解下列矩阵方程
1 2 3 1 1 1 X 2 2
0 1
1
B
1
1 1
0 2 0 1
1 0 4 2
0
1
1 0
解:把A, B分块成
11 00
AA
00 1 1
11 2 1
00 00
00 1 0
00 1010
EA1
O E
,
1 0 1 0
B
1 1
1
2 0 1
0 4 2
1 1 0
B11 B21
E
B22
则
AB
E A1
O E
B11 B21
A、B分块如下:
A11 A12 A1t
A
A21 As1
A22 As2
A2t Ast
B11 B12
B
B21 Bs1
B22 Bs 2
B1t B2t Bst
则定义A B Akl Bkl st
注. 分块矩阵运算中,每个子块具有二重性:一 是分块矩阵的元素;二是本身是矩阵。
A11
A12
A
A21 As1
A22 As 2
n1 n2
A1t
m1
各子块行数
A2t m2
Ast
ms
s
mk m
k 1
t
nt 各子块列数
nl n
l 1
说明 (2). 矩阵分块三原则:体现原矩阵特点,依
据问题需要,子块可以作元素运算。
一、分块矩阵的运算规则
1、分块加法
设A、B是m×n阶矩阵,采用相同的分块法分块将
A11 A
As1
A1r
A1T1
, 则 AT ( AT )T
Asr
ij
A1Tr
AsT1
.
AsTr
分外层、内层双重转置
说明:分块转置中,每个子块一方面作为分块阵 元素要转置;另一方面作为矩阵本身也要转置。
特别地,对于列分块矩阵: A1T
A ( A1, A2 , , At ) AT
例如
A
0 1
a 0
0 b
0 1
0 1 1 b
B1 B2 ,
B3
a 1 0 0
即
A
0 00
a
1 1
0
1 1
0 bb
B1 BB32
a 1 0 0
0
1 0
a 0 1
0 b 1
0 b1
A E
O B
,
其中OBEA
ab01 10
01 0ba1
a
0
1 0
1 a 0 1
0 0 b 1
E B22
B11 A1B11
B21
A1
E B22
又
A1
B11
B21
1 1
12
1 1
02
1 1
01
3 0
24
1 1
01
2 1
14
A1
B22
1 1
12
4 2
10 33
13
于是
AB
B11 A1B11
B21
A1
E B22
1 0 1 0
1 2 1
4 4 1
0 3 3
1
3 1
说明 (3). 矩阵分块的目的,是让矩阵的计算过程
分块矩阵
矩阵分块,是矩阵运算的一个重要方法,可将大 规模矩阵的运算化为若干小矩阵进行计算。
一. 分块矩阵的运算规则 二. 分块矩阵的一些例子
1、矩阵分块的方法
在矩阵某些行之间插入横线,某些列之间插入纵 线,将矩阵分割成若干个小矩阵,每个小矩阵称为 矩阵的子块;以子块为元素的矩阵,称为分块矩阵。
a 1 0 0
2、分块数乘 设A是m×n阶矩阵,任意分块,k是常数,则定义
kA kAkl st
3、分块乘法 设A是m×l阶矩阵,
B是l×n阶矩阵, 即A的列数 = B 的行数 分块A = ( Auv )s×r
B = ( Bvw )r×t 即A的列分块法 = B 的行分块法 则A与B的乘积C = ( Cuw ) 是s×t阶分块矩阵,满足
r
Cuw Auv Bvw v1
(u 1, , s; w 1, ,t)
注. 分块矩阵乘积AB中,每个子块:
r
(1)作为分块阵元素参与运算Cuw AuvBvw v 1
(2)作为矩阵也要满足乘法条件 AuvBvw
例1. 用分块矩阵法求AB
1
A
0 1 1
0 1 2 1
0 0 1 0
0
0
阵形式类似。
3. 矩阵乘积AB,A不分块,B按列分块
设矩阵A、B分别是s×n 和n×t 阶矩阵,A不分块,
B按列分块,即
B (1, 2 , , t )
则 AB A(1, 2 , , t )
( A1, A2 , , At )
说明:矩阵方程AX = B 可看成 t 个线性方程组 Ax1 = b1, Ax2 = b2, …, Axt = bt
更简单,计算量更少。
例1的计算量比较: 直接进行矩阵乘积需要的四则运算次数
4 4 (4 3) 112 用分块矩阵进行矩阵乘积需要的四则运算次数
块运算:2 2 (2 1) 12 子块运算:2 2 (2 1) 2 2 2 20
合计32次
4、分块转置 设矩阵A = ( Aij ) 是s×r 阶分块矩阵