兰州大学信息科学与工程学院 矩阵理论第 9讲 -1
矩阵理论 -第九讲
兰州大学信息科学与工程学院
2004年
兰州大学信息科学与工程学院 矩阵理论第 9讲 -2
上节内容回顾
? 矩阵的条件数
– 定义矩阵条件数的工程背景
– 矩阵的奇异值
? 矩阵序列
– 矩阵序列收敛的充分必要条件
– 收敛矩阵
? 矩阵级数
– 矩阵级数的绝对收敛的充要条件
绝对收敛 收敛
(){,,0,1,,}k m nA A C k??? ()l im 0k
k AA?? ??
()lim k
k AA?? ?
,mnCR??? ? ?()0 kk A???
lim kk A?? ? 0nnAC??
()0 kk A???
兰州大学信息科学与工程学院 矩阵理论第 9讲 -3
矩阵的幂级数
– 矩阵幂级数
设,,称矩阵级数
为矩阵 A的幂级数
– 方阵幂级数收敛的判别定理
若复变数幂级数 的收敛半径为 r,而矩阵 的谱半径
为,则
1,当 时,方阵幂级数 绝对收敛
2,当 时,方阵幂级数 发散
证明,
1.,取
,使得
( 0,1,,)ka C k??nnAC??
0 kkk aA
?
??
0 kkk aA
?
??
0 kkk aA
?
??
0 kkk az
?
??
nnAC??
()A?
()Ar? ?
()Ar? ?
()Ar? ? 0 ( )rA?? ?? ? ? 0 ?????
()Ar????, nnm CR????
()mA A r??? ? ?
兰州大学信息科学与工程学院 矩阵理论第 9讲 -4
矩阵的幂级数
由于幂级数
收敛,根据正项级数的比较审敛法知矩阵幂级数
绝对收敛
2,由于,设,则
当 时,
由 Jordan定理,,使得
( ( ) )kk k kk k k kmmma A a A a A a A??? ? ? ? ? ? ?
0 ( ( ) ) kkk aA??
?
? ??
0 kkk aA
?
??
( ) m a x jjA??? m a xljj??? () lA???
()Ar? ? l r? ?
11
21
1
( 1 0)i
n
n
P AP J or
??
?
?
?
?
?
?
??
??
? ? ?
??
??
nnnPC???
兰州大学信息科学与工程学院 矩阵理论第 9讲 -5
矩阵的幂级数
矩阵幂级数
的对角线元素为
由于 发散,从而矩阵幂级数 发散
由于矩阵幂级数

具有相同的敛散性,可知
也发散。
– 推论
设幂级数 的收敛半径为 r,。 若
使得,则矩阵幂级数 绝对收敛
0 kkk aJ
?
??
0 ( 1,,)kkjk a j n?
?
? ??
0 kklk a ?
?
?? 0 kkk aJ
?
??
()0 kk A??? ()0 kk P A Q???
0 kkk aA
?
??
0 kkk az
?
??, nnCR????
nnAC??
Ar? 0 kkk aA???
()AA? ?
,nnCR????
兰州大学信息科学与工程学院 矩阵理论第 9讲 -6
矩阵的幂级数
– 举例
判断矩阵幂级数
的敛散性
解:令
>> eig(A)
ans =
0.8333
-0.5000
由于幂级数 的收敛半径为 r = 1
绝对收敛
0
18
216
k
kk
k?
?
???
??????
181
216A
???? ??
???
( ) 0, 8 3 3 3 1A? ??
0 kk kz
?
??
0
18
216
k
kk
k?
?
???
??????
兰州大学信息科学与工程学院 矩阵理论第 9讲 -7
矩阵的幂级数
– Neumann级数收敛充要条件
设,称矩阵幂级数 为 Neumann级数
收敛
并且在此级数收敛时,其和为
证明,
充分性,
幂级数 的收敛半径为 1
必要性:若矩阵幂级数
收敛,记
,,则
nnAC??
0 kk A
?
??
0 kk A
?
?? ( ) 1A? ?
( ) 1A? ?
0 kk kz
?
??
0 kk A
?
??
收敛
0 kkSA
?
?? ? () 0
nnk
kSA?? ?
()lim n
n SS?? ?
0 kk A
?
??
( ) ( 1 ) ( ) ( 1 )l im l im ( ) l im l imn n n n n
n n n nA S S S S
??
? ? ? ? ? ? ? ?? ? ? ? ? 0
( ) 1A? ?
1()IA??
收敛矩阵的充要条件
兰州大学信息科学与工程学院 矩阵理论第 9讲 -8
矩阵的幂级数
当 收敛时,

可逆
由于
0 kk A
?
?? ( ) 1A? ? 0 1 ( )A?? ?? ? ?
0 ????? ( ) 1A????, nnCR????
( ) 1AA??? ? ? IA?
( ) 2()nnS I A I A A A? ? ? ? ? ?
21nnA A A A ?? ? ? ? ?
1nIA???
( ) 1 1 1( ) ( )nnS I A A I A? ? ?? ? ? ?
( ) 1l im ( )n
nS S I A
?
??? ? ?
A是收敛矩阵 lim n
n A?? ? 0
兰州大学信息科学与工程学院 矩阵理论第 9讲 -9
矩阵的幂级数
– 举例

判断矩阵幂级数 的敛散性,若收敛,求其和
解,norm(A,1)
ans = 0.9000
即,所以 绝对收敛
inv(eye(size(A))-A)
ans =
2.0000 1.0000 1.0000
3.1429 4.4286 3.0000
1.4286 1.7857 2.5000
0,2 0,1 0,2
0,5 0,5 0,4
0,1 0,3 0,2
A
??
???
????
0 kk A
?
??
1 0,9 1A ?? 0 kk A
?
??
0 kk A
?
?? ?
兰州大学信息科学与工程学院 矩阵理论第 9讲 -10
矩阵函数
– 定义,
矩阵函数的定义基于收敛的矩阵幂级数 。
收敛于一个唯一的矩阵,即此矩阵幂级数的和 S。 这样,矩
阵幂级数在矩阵 与 之间建立了一个映射,
称此映射为矩阵函数,它是以矩阵为变量(更为确切地,以方阵为变
量)且取值为矩阵(方阵)的一类函数。
称 S为 A在映射 f下的象,记作,
,n n n nf C C???
0 kkk aA
?
??
0 kkk aA
?
??
()S f A?
nnC? nnC?
0 ()!
k z
k
z er
k
?
? ? ? ? ??
21
0
( 1 ) s i n ( )
( 2 1 ) !
k
k
k z z rk
? ?
?
? ? ? ? ?
??
2
0
( 1 ) c o s ( )
( 2 ) !
k
k
k z z rk
?
?
? ? ? ? ??
10 ( 1 ) ( 1 )kk z z r? ?? ? ? ??
1
0
( 1 ) l n ( 1 ) ( 1 )
1
k k
k z z rk
? ?
?
? ? ? ?
??
兰州大学信息科学与工程学院 矩阵理论第 9讲 -11
矩阵函数
相应地,根据矩阵幂级数的收敛准则,将矩阵幂级数
的和分别记为下列矩阵函数
0 ()!
k nn
k
A AC
k
? ?
? ???
21
0
( 1 ) ()
( 2 1 ) !
k
k n n
k A A Ck
? ??
?
? ??
??
2
0
( 1 ) ()
( 2 ) !
k
k n n
k A A Ck
? ?
?
? ???
0 ( ( ) 1 )kk AA ?
?
? ??
1
0
( 1 ) ( ( ) 1 )
1
k k
k AAk ?
? ?
?
? ?
??
1,s i n,c o s,( ),l n ( )Ae A A I A I A???
兰州大学信息科学与工程学院 矩阵理论第 9讲 -12
矩阵函数在矩阵分析中的应用
许多工程问题,常常化为求解一阶常系数微分方程组的问题
– 由线性元件构成的网络状态方程组及输出方程组
– 其它动态系统或受控系统
L
C
R2
R1
u(t) iL iC
x Ax Bu??
y C x D u??
m
y(t)
F(t)
x Ax Bu??
兰州大学信息科学与工程学院 矩阵理论第 9讲 -13
矩阵函数在矩阵分析中的应用
– 离散时间系统
假设上述方程组的初始条件为 或
先考察 u(t) = 0时,的解,这时状态方程组简化为
这相当于求系统的零输入响应
当矩阵 A为数 a时,其解为
可以设想,当 而 时,的解含有
,可以证明 都是收敛的,因而其和是有意义的
离散时间系统 x(n) y(n)
( ) ( 1 ) ( )m n G m n H x n? ? ?
( ) ( ) ( )y n C m n D x n??
x Ax Bu??
x Ax?
00()x t x? 0(0)xx?
0() 0() a t tx t e x??
AtennAC??nxC? x Ax Bu??
tR???
0
1 ()
!
k
k Atk
?
??
0 ()!
k A n n
k
A e A C
k
? ?
? ? ? ??
兰州大学信息科学与工程学院 矩阵理论第 9讲 -14
矩阵函数在矩阵分析中的应用
– 矩阵函数, 及 满足代数三角函数的性质 Euler公式
sinAcosA Ae
s i n ( ) s i nAA? ? ? ?
c o s ( ) c o sAA? ? ?
c o s s i nAe A i A? ? ?
c o s 2i A i AeeA ????
s i n 2i A i AeeA i ????
0 !
ki A k
k
ieA
k
?
?? ?
2 2 1
00
( 1 ) ( 1 )
( 2 ) ! ( 2 1 ) !
c o s s in
kk
kk
kk A i Akk
A i A
?? ?
??
????
?
??
??
c o s ( ) s i n ( )iAe A i A? ? ? ? ?
c os sinA i A??
c o s 2iA iAeeA ??? s in 2
i A i AeeA
i
???
Ae
兰州大学信息科学与工程学院 矩阵理论第 9讲 -15
矩阵函数在矩阵分析中的应用
– 常用矩阵函数的性质
设,且
()A B A B B Ae e e e e?? ? ?
s i n ( ) s i n c o s c o s s i nA B A B A B? ? ? ?
c o s ( ) c o s c o s s i n s i nA B A B A B? ? ? ?
,nnA B C ?? AB BA?
sin 2 2 sin c o sA A A?
22c o s 2 c o s s i nA A A??
00
11( ) ( )
!!
A B k k
kke e A Bkk
??
??? ??
0
1 ()
!
A B k
ke A Bk
??
????
221( ) ( )
2!I A B A A B B A B? ? ? ? ? ? ?
3 2 2 2 2 31 ()
3! A A B A B A B A A B A B B A B B? ? ? ? ? ? ? ?
221( ) ( 2 )
2!I A B A A B B? ? ? ? ? ?
3 2 2 31 ( 3 3 )
3! A A B AB B? ? ? ? ?
充要条件
兰州大学信息科学与工程学院 矩阵理论第 9讲 -16
矩阵函数在矩阵分析中的应用
– 常用矩阵函数的性质

若, 是 A的特征值,则矩阵函数 的特
征值为
由 Jordan定理,,使得
eI??0
nnAC??
trd e t AAee??
1()AAee????
nnAC??
12,,n? ? ? ()fA
12( ),( ),( )nf f f? ? ?
11
21
1
( 1 0)i
n
n
P AP J or
??
?
?
?
?
?
?
??
??
? ? ?
??
??
nnnPC???
兰州大学信息科学与工程学院 矩阵理论第 9讲 -17
矩阵函数在矩阵分析中的应用
再设,可得如下矩阵幂级数
收敛。由于
即矩阵幂级数 收敛,由于 的对角线元素为
所以,这些复数项幂级数收敛,且
0 ( 1,,)kkjk a j n?
?
? ??
1100( ) ( )kkkkkkP a A P a P A P?????? ???
0() kkkf A a A
?
?? ?
10 ()kkk P a A P? ???
1
0
0
() kkk
k
kk
a P A P
aJ
? ?
?
?
?
?
?
?
?
0 kkk aJ
?
??
0 ( ) ( 1,,)kk j ik a f j n??
?
? ???
0 kkk aJ
?
??
10 0 0d e t( ) ( ) ( )k k kk k k nk k kI a J a a? ? ? ? ?? ? ?? ? ?? ? ? ?? ? ?
12( ( ) ) ( ( ) ) ( ( ) )nf f f? ? ? ? ? ?? ? ? ? ? ? ?
110( ( ) ) ( )kkkP a A P P f A P???????
兰州大学信息科学与工程学院 矩阵理论第 9讲 -18
矩阵函数在矩阵分析中的应用
由于相似矩阵具有相同的特征值,所以 的特征值为
由此,的特征值为
0( ) ~ kkkf A a J
?
??
12( ),( ),( )nf f f? ? ?
()fA
Ae
12,,ne e e ???
12d e t nA ? ? ?? ? ? ?
1 1 2 2 1 2tr n n nA a a a ? ? ?? ? ? ? ? ? ? ?
12d e t nAe e e e ???? ? ? ?
12
tr
n
A
e
e
? ? ?? ? ??
? ?0 Ae 可逆
AAe e e I? ??0 1()AAee???
兰州大学信息科学与工程学院 矩阵理论第 9讲 -19
矩阵函数在矩阵分析中的应用
– 需要注意的几点
除非 A为对角矩阵
1 1 1
1
n
n n n
aa
A
aa
??
???
????
111
1
n
n nn
aa
A
aa
ee
e
ee
??
??
?
??
??
11 1
1
s i n s i n
sin
s i n s i n
n
n nn
aa
A
aa
??
???
??
兰州大学信息科学与工程学院 矩阵理论第 9讲 -20
矩阵函数在矩阵分析中的应用
? 的求法举例
– 利用 Hamilton-Cayley定理
已知

解,
由 Hamilton-Cayley定理
0
0A
?
?
??? ??
???
Ate
Ate
22( ) d e t ( )IA? ? ? ? ?? ? ? ?
()A? ? 0 22AI???0
2 2 2( ) ( 1 )k k k kA I I??? ? ? ?
2 1 2 2 1( ) ( 1 )k k k k AA I A??
?
??? ? ? ?
( 1,2,)k ?
2 4 3 5
0
( ) ( ) ( ) ( ) ( )( 1 ) ( 1 )
! 2 ! 4 ! 3 ! 5 !
kAt
k
A t t t t t AeI
k
? ? ? ?
?
?
?? ? ? ? ? ? ? ? ??
c o s s i nc o s s i n
s i n c o s
ttAt I t
tt
????
???
??? ? ? ? ??
???
兰州大学信息科学与工程学院 矩阵理论第 9讲 -21
矩阵函数在矩阵分析中的应用
? 的求法举例
– 利用相似对角化

同理
Ate
1 12(,,)nP A P d i a g ? ? ?? ? ? ?
100( ) ( )kkkkkkf A a A a P P?? ???? ? ???
1
120 0 0
1
12
d ia g (,,,)
d ia g ( ( ),( ),,( ) )
k k k
k k k nk k k
n
P a a a P
P f f f P
? ? ?
? ? ?
? ? ? ?
? ? ?
?
?
?
? ? ?
112( ) d i a g ( ( ),( ),,( ) )nf A t P f t f t f t P? ? ? ??
4 6 0
3 5 0
3 6 1
A
??
??? ? ?
??????
兰州大学信息科学与工程学院 矩阵理论第 9讲 -22
矩阵函数在矩阵分析中的应用
? 的求法举例
– 利用 Jordan标准形
Jordan块的幂
Ate
1
1
ii
i
i rr
i
i
JC
?
?
?
?
??
??
??
??
??
nnAC?? nn
nPC???
1
21
s
J
J
P A P J
J
?
??
??
??
??
??
( 1 )
( ) ( ) ( )
1 ! 2 ! ( 1 ) !
()
1!
()
2!
()
1!
i
i
rk k k
k
i
k
k
k
k
i
k
k
k
r
J
??
? ? ?
?
?
?
?
?
?
?
?
?
??? ??
??
?
??
?
??
? ??
??
?? ?
??
??
兰州大学信息科学与工程学院 矩阵理论第 9讲 -23
矩阵函数在矩阵分析中的应用
1 ( 1 )2
2
00
( ) ( ) ( )
1 ! 2 ! ( 1 ) !
()
1!
()
()
2!
()
1!
ii
i
rrk k k
k
i
k
k
kk
k
i k i k
kk k
k
k
t
t t t
r
t
f J t a J t a
t
t
??
? ? ?
?
?
?
?
?
?
?
??
??
??
?
?? ? ??
??
?
??
?
??
?? ??
??
?? ?
??
??
??
兰州大学信息科学与工程学院 矩阵理论第 9讲 -24
矩阵函数在矩阵分析中的应用
1 ( 1 )2
0 0 0 0
00
2
00
0
0
( ) ( ) ( )
1 ! 2 ! ( 1 ) !
()
1!
()
()
2!
()
1!
ii
i
rrk k k
k
k k k k
k k k k i
k
k
kk
kk
k
k
i kk
kk
k
k
k
k
k
k
t
t t t
a a a a
r
t
aa
t
f J t aa
t
a
a
??
? ? ?
?
?
?
?
?
?
?
??? ? ? ?
? ? ? ?
??
??
??
??
?
?
?
?
?
?? ? ??
??
?
??
?
??
??
? ??
??
?
??
??
??
? ? ? ?
??
??
?
?
兰州大学信息科学与工程学院 矩阵理论第 9讲 -25
矩阵函数在矩阵分析中的应用
1 ( 1 )2
2
( ) ( ) ( )
()
1 ! 2 ! ( 1 ) !
()
()
1!
()
()
()
2!
()
1!
()
ii
i
rr
i
i
t
tf t f t f
f
r
tf
f
f J t
tf
f
tf
f
??
? ? ?
?
?
?
?
?
?
?
??
?
??? ??
??
?
?? ?
??
? ??
??
???
??
??
兰州大学信息科学与工程学院 矩阵理论第 9讲 -26
矩阵函数在矩阵分析中的应用
100( ) ( )k k k kkkkkf A t a A t a P J P t?? ???????
1
0
1
0
12
0
0
1
1
()
()
()
kk
k
k
kk
k
k
kk
k
k
kk
ks
k
s
P a J t P
a J t
a J t
PP
a J t
f J t
PP
f J t
?
?
?
?
?
?
?
?
?
?
?
?
??
??
???
??
??
??
??
?
??
??
??
?
?
?
?
兰州大学信息科学与工程学院 矩阵理论第 9讲 -27
矩阵函数在矩阵分析中的应用
已知

1 0 1
1 2 0
4 0 3
A
???
???
?????
Ae
1
2
1
tt
J J t t
t
t
t
ee
e e e
e
?
?
??
??
??
??
??
11
1
2
J
??
???
??
??
兰州大学信息科学与工程学院 矩阵理论第 9讲 -28
矩阵函数在矩阵分析中的应用
? 的求法举例
– 利用矩阵多项式
为计算
利用多项式的带余除法,将 表示为如下形式,
由 Hamilton-Cayley定理

因为
所以
Ate
0() kkkkf A t a A t
?
?? ?
()ft?
( ) (,) ( ) (,) d e g (,)f t g t r t r t n? ? ? ? ? ?? ? ?
()A? ? 0
( ) (,)f A t r A t?
11 1 0(,) ( ) ( ) ( )nnr t b t b t b t? ? ???? ? ? ?
() ( ) 0 ( 0,1,,1 ; 1,2,,)l ii l r i s?? ? ? ? ?
1212( ) ( ) ( ) ( ) srrr s? ? ? ? ? ? ? ?? ? ? ?
( ) (,) ( 0,1,,1 ; 1,2,,)
i i
ll
ill
ddf t r t l r i s
dd ??
??
??
?? ?
?
? ? ? ?
兰州大学信息科学与工程学院 矩阵理论第 9讲 -29
矩阵函数在矩阵分析中的应用
已知

1 0 1
1 2 0
4 0 3
A
???
???
?????
Ate
() tf t e?? ?
2( ) ( 1 ) ( 2 )? ? ? ?? ? ?
22 1 0(,) ( ) ( ) ( )r t b t b t b t? ? ?? ? ?
2 1 0( 1,) ( ) ( ) ( ) tr t b t b t b t e? ? ? ?
21( 1,) 2 ( ) ( ) tr t b t b t t e? ? ? ?
22 1 0( 2,) 4 ( ) 2 ( ) ( ) tr t b t b t b t e? ? ? ?
22 1 0( ) ( ) ( )Ate b t A b t A b t I? ? ?
利用矩阵 A的最小多项式,
还可以进一步降低 的
次数,从而使待定的系数减
少,进一步化简计算。
(,)rt?
兰州大学信息科学与工程学院 矩阵理论第 9讲 -30
矩阵的微分和积分
– 以函数为元素的矩阵 ——函数矩阵
? 函数矩阵的微分和积分
– 泛函
? 数量函数对矩阵变量的导数
– 向量值函数或矩阵值函数对向量变量或矩阵变量的导数
**********************************************************************
– 函数矩阵的微分和积分
定义
以变量 t的函数为元素的矩阵
是定义在 [a,b]上的,若
在 [a,b]上连续、可微、可积,若每个 在 [a,b]上连续、可
微、可积
,mnf C F? ?
( ) ( ( ) )ij m nA t a t ??
[,]t a b?()At
()At ()ijat
兰州大学信息科学与工程学院 矩阵理论第 9讲 -31