当前位置:文档之家› 矩阵理论课件 第五章 特征值的估计与广义逆矩阵

矩阵理论课件 第五章 特征值的估计与广义逆矩阵

4个盖尔圆中只有 G4 是孤立的, G1,G2 ,G3 是连通
的,故结论成立。
定义1 (严格对角占优矩阵)
设 A (aij ),若C满n足n
n
aii aij , i 1, 2, n j 1 ji
则称 A 为(行)对角占优矩阵,若不等式严格成立, 则称 为A(行)严格对角占优矩阵;若 为A行T (严格)对角占优矩阵,则称 A列(严格)对角占
min x : H x min Ax b
xH
2
xC n
2
定义5 (广义逆的一般定义)
设 A C,m若n存在
设 A (aij )nn Rnn(n阶实矩阵),则
Im i
n(n 1)
2
max
1i , jn
cij
例1 估计下面矩阵的特征值的界:
0 0.2 0.1 1 0
A
0.2
0
0.2
2 0.3i
解:
0.1 0.2 0 3 0.3i
B 1 ( A AT ) 0,C 1 ( A AT ) A
的并集内,其中 Ri aij 。 j 1 ji
上述圆盘称为Gerschgorin圆盘,简称盖尔圆。
证明: 设Ax x, x ( x1, x2 , , xn )T
n
aij x j xi , i 1, 2, , n
j 1
n
aij x j ( aii )xi , i 1, 2, , n
1 2 1 0 0 1 2 1 0 0
0
00
1
0
0
00
1
0
2 4 0 0 1 0 0 2 0 1
1 0 1 0 0
0
00
1
0
0 0 2 0 1
1 0 0
P
0
1
0
2 0 1
1 2
Q
0
1
A
Q
E1 Y
X
Z
P
1 0 0
1
0
2 1
1
y
x1 z1
x2 z2

A
3 5
5 3
3 5 E A
5 3 ( 3)2 25 0
1,2 3 5i
2个盖尔圆 3 5

A
2 1
1
4
2 1 E A
1 4
2 1 ( 3)2 0 1,2 3
2个盖尔圆
4 1
1 0.8

A
1 0.5
0.8
0
E A 0.5
2 0.4 0
( A A)H ,A只要A方程组
有解A,x其则b
x A就是b 它的最小范数解(唯一)。
证明:
x2 2
Ab (En
A A)z
2 2
Ab (En A A)z H Ab (En A A)z
Ab H Ab (En A A)z H (En A A)z
Ab H (En A A)z (En A A)z H Ab
对于方程组 Ax ,b若 rank( A, b) rank( A)
则方程组有解(称为相容方程组,否则称为不相容), 所有解里面2-范数最小的解,称为其最小范数解,即
x* : min x
Ax b
2
定理4 (最小范数解的性质)
通常记为 A{1,3}
若 A是 矩阵 A C的m一个n 广义{1逆} ,且满足
定理5 (最小二乘解的性质) 通常记为 A{1, 4}
若 A是 矩阵 A C的m一个n 广义{1逆} ,且满足
( AA )H ,则AA对
,b C m x Ab
是不相容方程组 Ax 的b最小二乘解(不唯一)。
定义4 (极小最小二乘解)
对于不相容方程组 Ax ,b所有最小二乘解里面
2-范数最小的解,称为极小最小二乘解,即
Ab H (En A A)z 0 (En A A)z H Ab 0
Ab H Ab (En A A)z H (En A A)z
Ab 2 2
(En
A A)z
2 2
Ab 2 2
1 2 1
例7
已知
A
0
0
,
b
0
,求
Ax b
2 4 2
的最小范数解。
0 2
1
0
0 1
1 0 0
1 2y
y
x1 2z1 z1
x2
2z2 z2
0 2
1 0
0
1
A
1
2
y
y
2( x2 2z2
2z2 )
x1 2z1 z1
x2 2z2
z2
A
A
1 2
y
y
2(1 2 y)
2y
( A A)H
1 2y
2(1
2
y)
y
2
y
2
2(1 2 y) y y
上述形式),则 的A所有 广{1}义逆的集合为:
A{1}
Q
Er Y
其中X 分别是 X ,Y , 相应Z阶数P的 任意矩阵。
Z
证明: AGA A PAGAQ PAQ
PAQQ1GP 1PAQ PAQ
Er O
OW
O
Y
X Er
Z
O
O O
Er O
O
O
W
O
O O
1 0 0 0 0 1 0
0 0
1 0
10 01
1 0
0
0
2 0
0 0 0 0 1 1 1
0 0 1 0
P
0
1
0
0
1 0 2 0
0
1
1
1
0 0 1
Q
0
1
1
1 0 0
1 0 0 0 0 1 1 0 0 0 0 1
0
1
00
1
0
0
1
00
1
1
0 0 1 1 0 0 0 0 1 1 0 0
2
2
i 3 0.2 0.6 Rei 3 0 0
Im i
3 2 0.2 2
Imi 3 0.2 0.6
0.3464
§2、圆盘定理 定理1 (圆盘定理)
设 A (aij )nn Cnn ,则 A 的特征值都在复数平面
上的n个圆盘 z aii Ri , i 1, 2, , n
n
0 0 2 1 0 0 0
1 0
1 0
00 10
1 0
0
0
1 0
0 0 0 0 1 1 1
0 0 0 1 0 2 0
1 0
1 0
00 10
1 0
0
0
1 0
0 0 0 0 1 1 1
0 0 1 0 0 1 0
1 0
1 0
00 01
1 0
0
0
2 0
0 0 0 0 1 1 1
v1 v2 vn ,则有下列结论:

i
n
max
1i , jn
aij
;

Re i
n max 1i , jn
bij
;

Im i
n
max
1i , jn
cij
.
其中
1 bij 2 (aij a ji )
1 cij 2 (aij a ji )
i, j 1, 2, , n
推论2 (Bendixson定理)
j 1
ji
n

xk
max
1 i n
xi
0
( akk )xk
akj x j
j 1
jk
n
n
akk xk akj x j akj x j
j 1
j 1
jk
jk
n
akk akj
j 1
xj xk
n
akj Rk
j 1
jk
jk
例2 估计下面矩阵的特征值范围:
1 0.1 0.2 0.3
1,2
1
i 2
0.6
1 0.8
2个盖尔圆
0.5
注:定理2表明由一个盖尔圆组成的连通部分,有且只 有一个特征值。
推论1 如果n阶矩阵 A C nn 的n个盖尔圆两两不相 交,则 A 必相似于对角矩阵。(仅是充分条件) 推论2 如果n阶矩阵 A Rnn 的n个盖尔圆两两不相
交,则 A 的特征值全为实数。
实际计算得到的特征值:
1 3.0253, 2 1.0757, 3 1.1118, 4 3.9893.
定理2 (圆盘定理的推广)
矩阵 A 的任一由k个盖尔圆组成的连通部分里,有且只 有k个特征值(当 A 的主对角线上有相同元素时,则按
重复次数计算,有相同特征值时亦需按重数计算)
例3 应用举例:
对 z Cn , 令 b Az ,则 Ax b 有解 x Gb
AGb b AGAz Az AGA A
设 y是方程组 Ax的解b
AGA A AGAy Ay b AGb b
即 x Gb是方程组 Ax b的解,从而结论成立。
定 设 其理 中A2P(,C{Q分1,}m别若n是广m义阶逆P和的An阶Q计可算,逆)矩E阵r (O即初 等ra变n换k(可A化)为 r
Er O
O
O
W Er
G
Q
Er Y
X
Z
P
0 0 2
例6
求矩阵
A
1
1
0
的一个{1} 广义逆。
0 0 1
1
1
1
相关主题