对于行数和列数较高的矩阵A , 运算时常采用分块法
使大矩阵的运算化成小矩阵的运算。
将矩阵A用若干条纵线和横线分成许多个小矩阵 ,
每一个小矩阵称为A的子块
以子块为元素的形式上的矩阵称为分块矩阵
例如将3×4矩阵A=分成子块的分法很多
下面举出三种分块形式:
(i) ,
(ii) ,
(iii)
分法(i)可记为A=
其中A11= ,
A12= ,
A21=(a31 , a32) ,
A22=(a33 , a34)
即A11 , A12 , A21 ,
A22为A的子块 , 而A形式上成为以这些子块为元的分块矩阵
分法(ii)及(iii)的分块矩阵也可类似写出
本章第2节证明公式|AB|=|A||B|时出现的矩阵及
正是分块矩阵
在那里是把四个矩阵拼成一个大矩阵
这与把大矩阵分成多个小矩阵是同一个概念的两个方面
分块矩阵的运算规则与普通矩阵的运算规则相类似 , 分别说明如下:
规则(i) 设矩阵A与B的行数相同、列数相同 , 采用相同的分块法
有A= ,
B=
其中A与B的行数相同、列数相同
那么A+B=
规则(ii)
设A= , λ为数 ,
那么λA=
规则(iii) 设A为m×l矩阵 , B为l×n矩阵
分块成A= ,
B=
其中Ai1 , Ai2 , ⋯ ,
Ait的列数分别等于B1j , B2j , ⋯ ,
Btj的行数
那么AB=
例17设A= ,
B= , 求AB
解: 把A , B分块成
A== ,
B=
则AB==
而A1B11+B21=+=+=
A1+B22=+=
于是AB==
规则(iv)
设A= ,
则AT=
分块矩阵转置 = 分块行列互换 + 每个子块自身转置
例子
设:
我们按虚线分块,上方为 2 行 4 列,分为左右两个子块;
下方为 1 行 4 列,也分为左右两个子块。
更规范地写成
分块矩阵(按列分两块,每块 2
列):
其中,
转置过程:
1. 交换分块位置:
位置放
位置放
位置放
位置放
2. 各子块分别转置:
,
3.
写出转置的分块矩阵:
验证: 直接对原始矩阵
整体转置(不看作分块):
,转置后:
与上面分块方法得到的结果相同。
规则(v) 设A为n阶方阵 , 若A的分块矩阵只有在对角线上有非零子块 ,
其余子块都为零矩阵 , 且在对角线上的子块都是方阵
即A= ,
其中Ai (i=1 , 2 , ⋯ , s)都是方阵
那么称A为分块对角矩阵
分块对角矩阵的行列式的性质:
|A|=|A1||A2|⋯|As|
由此性质可知 , 若|Ai|≠0(i=1 , 2 , ⋯ , s) , 则|A|≠0
并有A-1=
例18设A= , 求A-1
解:
因为A==
A1=(5) ,
=()
,
A2= ,
=
所以A-1=
对矩阵分块时 , 有两种特殊分块法 , 这就是按列分块和按行分块
m×n矩阵A有n列 , 称为矩阵A的n个列向量
若第j列记作aj= ,
则A可按列分块为A=(a1 , a2 , ⋯ , an)
m×n矩阵A有m行 , 称为矩阵A的m个行向量
若第i行记作=(ai1
, ai2 , ⋯ , ain) ,
则A可按行分块为A=
对于矩阵A=(aij)m×s与矩阵B=(bij)s×n的乘积矩阵AB=C=(cij)m×n
若把A按行分成m块 , 把B按列分成n块
便有AB=(b1
, b2 , ⋯ ,
bn)==(cij)m×n
例19试证明矩阵A=O的充分必要条件是方阵ATA=O
证:
必要性:因为A=O,得AT=O ,所以ATA=O⋅O=O
充分性:
设A=(aij)m×n , 把A按列分块为A=(a1 ,
a2 , ⋯ , an),
则ATA=(a1
, a2 , ⋯ ,
an)=,即
ATA 的(i ,
j)元为
因ATA=O ,
有=0(i
, j=1 , 2 , ⋯ , n),
也有==,(j=1
, 2 , ⋯ , n)
因aij为实数,得a1j=a2j=⋯=amj=0
(1 , 2 , ⋯ , n)
所以A=O
本例阐明了矩阵A与方阵ATA之间的一种关系
特别地 , 当A=a为列向量时 , 由于aTa为1×1矩阵 ,
即aTa是一个数
这时 , 本例的结论可叙述为:
aTa=0是列向量a=0的充分必要条件
利用矩阵的按列(按行)分块 , 还可以给出线性方程组的另一矩阵表示形式
重新回到线性方程组 (1)
它的矩阵乘积形式为Am×nxn×1=bm×1
(1')
上式中 , 把A按列分块 , 把x按行分块
由分块矩阵的乘法有=b
即x1a1+x2a2+⋯+xnan=b
(10)
其实把方程组(1)表成
也即是(10)式.
(1)、(1')和(10)是线性方程组(1)的各种变形.
今后 , 它们与(1)将混同使用而不加区分 ,
并都称为线性方程组或线性方程
解与解向量亦不加区别