第四章 矩阵的标准型PPT幻灯片
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
最后,根据 J j ( i ) 的结构,设
p ij (p i(1 j),p i(2 j),L ,p i(n jij))
由 A pij pijJj( i),可知
(
A
i
I
)
p
( i
1 j
)
(
A
iI
)
p
( i
2 j
)
p
( i
1 j
)
L
(
A
iI
)
p ( ni j ) ij
p ( ni j 1) ij
称 J 为 A 的Jordan标准型。并称方阵
i
Ji (i )
1
i O
O
1
, i 1,2,L,s
i mimi
为 m i 阶Jordan 块。
数学系 李继根(jgli@ecust.edu.cn)
定理 2 设 ACnn。如果 A 的特征多项式可
分解因式为 () ( 1 ) m 1 L (s ) m s
JA A1(2)
A2(1)
因为特征值 `1 2 为单根,所以 A1(2) 2
并从 (A2I)x解得对应的特征向量为
1 (1,2,4)T
数学系 李继根(jgli@ecust.edu.cn)
对于二重特征值 `2,3 1 ,由 (AI)x
只解得唯一的特征向量为
2 (1,1,1)T
因此 A 2 ( 1 ) 中只有一个Jordan块,即
i
阶的Jordan子矩阵,有 k
个
i
阶数为 n i j ( n i 1 n i 2 L n i k i n i ) 的
Jordan块,即 A i ( i ) d i a g ( J 1 ( i ) , J 2 ( i ) , L , J k i ( i ) )
数学系 李继根(jgli@ecust.edu.cn)
其中, pi(1 j)(j1 ,2,L,ki)是矩阵 A 关于特征
值 i 的一个特征向量,
的广义特征向量,称
p
( i
n j
i
j
)
pi(2j),L,
为 i 的
p(ni j ) 则称为 ij
n i j 级根向量。
iห้องสมุดไป่ตู้
数学系 李继根(jgli@ecust.edu.cn)
当所有的 n i j 1 时,可知 k i n i ,此时矩阵没
则 V 可分解成不变子空间的直和
V N 1 排 N 2L ? N s
这里 N iK er((T iE )m i)
数学系 李继根(jgli@ecust.edu.cn)
适当选取每个子空间 N i 的基(称为Jordan基), 则每个子空间的Jordan基合并起来即为 V 的Jordan
基,并且 V 在该Jordan基下的矩阵为块对角阵 J d i a g ( J 1 ( 1 ) , J 2 ( 2 ) , L , J s ( s ) )
x&P1x&P1APP1xP1Bu
AxBu
yCPxDuCxDu
显然,最简单的 A 就是 A 的Jordan标准型。此时
虽然没有实现状态变量间的完全解耦,但也达到了可
能达到的最简耦合形式。因此线性变换就是状态空间
的基底变换,其目的在于寻找描述同一系统的运动行
为的尽可能简单的状态空间描述。
数学系 李继根(jgli@ecust.edu.cn)
( m 1 m 2 L m s n )
则 A 可经过相似变换化成唯一的 Jordan标准型 J (不计Jordan块的排列次序),即存在可逆矩阵(称为 Jordan变换矩阵) PCnn 使
P1APJ
APJP 或者 A 有Jordan分解
1
数学系 李继根(jgli@ecust.edu.cn)
2 (1,2,1)T
对重根有几个特 征向量,就有几 个约旦块
因此 A 2 ( 1 ) 中只有一个Jordan块,即
1 1
A2
(1)
0
1
求解 (AI)2,可得所需的广义特征向量
(0,1,1)T
数学系 李继根(jgli@ecust.edu.cn)
综合上述,可得
0 1 0
2 0 0
P0 2 1, JA0 1 1
1 1 1
0 0 1
数学系 李继根(jgli@ecust.edu.cn)
所以原方程组变为
d dy tP 1d dx tP 1A xP 1A PyJAy
即 d dy t12y1, d dy t2y2y3, d dy t3y3
解得
y 1 c 1 e 2 t,y 2 c 2 e t c 3 te 2 t,y 3 c 3 e t,
数学系 李继根(jgli@ecust.edu.cn)
最后,由可逆线性变换 x P y 得原方程组的解
xx12
c2et c3t et 2c2et c3(2t
1)et
x3 c1e2t c2et c3(t 1)et
数学系 李继根(jgli@ecust.edu.cn)
例 5 现代控制理论中,线性定常系统(Linear time
2 1 1 1
A
2
1
3
2
1 1 0 1
1
1
2
2
数学系 李继根(jgli@ecust.edu.cn)
解: A 的特征值为 `1 0 ,`2`3`4 1 ,则
JA
A1(2)
A2(1)
因为特征值 `1 0 为单根,所以 A1(0) 0
并从 (A0I)x解得对应的特征向量为
1(1,3,1,2)T
有广义特征向量, p i 的各列是 i 的线性无关的特
征向量,因此Jordan块 J j(i)(j 1 ,L ,k i;
i1,L, t)都是一阶的,此时Jordan标准型为
J d ia g (1 4 1 ,L 2,4 3 1 ,1 4 2 ,L 2,4 3 2,L ,1 4 t,L 2,4 3 t)
求下列状态方程的约当标准型:
0 1 0 0 x&AxBu0 0 1x0u
2 3 0 1
这里矩阵 A 是特征多项式 | I A|的友矩阵。
数学系 李继根(jgli@ecust.edu.cn)
解:|I A | 3 3 2 ( 2 ) ( 1 ) 2 0
A 的特征值为 `12, `2`31,故设
二、 Jordan标准型的一种简易求法
把 A 的同一个特征值的若干个Jordan块排列在一起,
就得到Jordan标准型 J A d i a g ( A 1 ( 1 ) , A 2 ( 2 ) , L , A t ( t ) )
(n 1 n 2 L n t n )
其中 A i ( i ) 是 n
解: A
特征值为 `12,`2`31,所以设
JA
A1(2)
A2(1)
因为特征值 `1 2 为单根,所以 A1(2) 2
并从 (A2I)x解得对应的特征向量为
1 (0,0,1)T
数学系 李继根(jgli@ecust.edu.cn)
对于二重特征值 `2 `3 1 ,由 (AI)x
只解得唯一的特征向量为
数学系 李继根(jgli@ecust.edu.cn)
一、 Jordan标准型的概念
定理 1 设 T 是复数域 C 上的线性空间 V 上的 线性变换 。令 T 在V 的一组基下的矩阵表示 为 A ,如果 A 的特征多项式可分解因式为
() ( 1 ) m 1 L (s ) m s
( m 1 m 2 L m s n )
因此经过可逆线性变换 x& P x 后,系统矩阵 A 和
控制矩阵 B 分别为
2 0 0
A P1AP0 1
1
J
0 0 1
2
B
P 1B
1 9
1
.
1
数学系 李继根(jgli@ecust.edu.cn)
例 6 求矩阵 A 的 Jordan标准型 J A 和相应的
Jordan变换矩阵 P ,其中
这里 x(x 1,x 2,x 3)T ,d d x t(d d x t1,d d x t2,d d x t3)T ,
1 1 0
A
4
3
0
1 0 2
数学系 李继根(jgli@ecust.edu.cn)
由上例,存在可逆线性变换 x P y 使得
P1APJA
其中
0 1 0
2 0 0
P0 2 1, JA0 1 1
n 1
n 2
n t
即矩阵 A 是可对角化矩阵。显然正规矩阵是一类最
特殊的可对角化矩阵。
数学系 李继根(jgli@ecust.edu.cn)
例 3 求矩阵 A 的 Jordan标准型 J A 和相应的
Jordan变换矩阵 P ,其中
1 1 0
A
4
3
0
1 0 2
数学系 李继根(jgli@ecust.edu.cn)
数学系 李继根(jgli@ecust.edu.cn)
解这个方程组,可得到Jordan链
{pi(1 j),pi(2 j),L,pi(n jij)}
这个名称也可以这样理解:
p i ( n j i j ) A i I p i ( n j i j 1 ) A i I L A i I p i ( 1 j ) A i I
p i (p i1 ,p i2 ,L ,p ik i)
其中
p ij(j 1 ,2 ,L ,k i)是 n ni
阶矩阵。
j
由 Api piA i(i) ,可知
A p i j p i j J j ( i )( j 1 , 2 , L , k i )
数学系 李继根(jgli@ecust.edu.cn)
1 1
A2(1)
0
1
求解 (AI)2,可得所需的广义特征向量
(1,0,1)T
数学系 李继根(jgli@ecust.edu.cn)
综合上述,可得
1 1 1
2 0 0
P2 1 0, JA0 1 1
4 1 1
0 0 1
1 2 1
P
1
1 9
2
6
5 3
2
3
数学系 李继根(jgli@ecust.edu.cn)
数学系 李继根(jgli@ecust.edu.cn)
第四章 矩阵的标准型
数学系 李继根(jgli@ecust.edu.cn)
标准型的理论源自矩阵的相似性,因为相 似矩阵有许多相似不变量:特征多项式、特征 值(包括代n 数重数和几何重数)、行列式、迹 及秩等,并且特征向量也可以借助于可逆的相 似变换矩阵互相求出。这自然导出了寻找相似 矩阵集合中的“代表矩阵”的问题。“代表矩 阵”当然越简单越好。对于可对角化矩阵, “代表矩阵”就是特征值组成的对角矩阵。特 别地,对于正规矩阵,可逆的相似变换矩阵特 殊化为酉矩阵或正交矩阵。但是令人非常遗憾 的是:一般矩阵未必与对角矩阵相似!!!
数学系 李继根(jgli@ecust.edu.cn)
对于三重特征值 `2`3`41 ,由
(A(1)I)x解得两个特征向量为
2 ( 1 ,0 ,0 , 1 ) T ,3 ( 0 ,1 ,0 , 1 ) T
因此A 2 ( 1 ) 中有两个Jordan块,即
1 1 1
A2(1)
1 或 1
invariant , LTI )的状态空间描述为
x& Ax Bu
y
Cx
Du
这里矩阵 A 表示了系统内部状态变量之间的联系,
称为系统矩阵;矩阵 B 称为输入矩阵或控制矩阵;
矩阵 C 称为输出矩阵或观测矩阵;矩阵 D 称为直
接观测矩阵。
数学系 李继根(jgli@ecust.edu.cn)
做可逆线性变换 x& P x ,则
1 1 1
0 0 1
数学系 李继根(jgli@ecust.edu.cn)
例 4 用 Jordan标准型理论求解线性微分方程组
d x1 dt
x1
x2
d d
x2 t
4 x1
3 x2
d x3 dt
x1
2 x3
数学系 李继根(jgli@ecust.edu.cn)
解: 方程组的矩阵形式为
dx Ax dt
根据 J A 的结构,将Jordan变换矩阵 P 列分块为
P (p 1 ,p 2 ,L ,p t) 其中 p i 是 n n i 阶的矩阵。 由 APPJA ,可知
A p i p iA i(i)( i 1 ,2 ,L ,t)
数学系 李继根(jgli@ecust.edu.cn)
进一步,根据 A i ( i ) 的结构,将 p i 列分块为
1 1 1
数学系 李继根(jgli@ecust.edu.cn)
求解 (A+I)2,无解!! 求解 (A+I)3,可得所需的广义特征向量
(1, 0,1,0)T
综合上述,可得
1 1 0 1
P
3
0 1 0,
1 0 0 1
2 1 1 0
0
JA
1
数学系 李继根(jgli@ecust.edu.cn)
§1、矩阵的Jordan标准型
由于一般矩阵与对角矩阵不相似,因此我 们“退而求其次”,寻找“几乎对角的” 矩阵。这就引出了矩阵在相似下的各种标 准型问题,其中Jordan标准型是最接近 对角的矩阵,只在第1条对角线上取1或0。 弄清楚了矩阵相似的本质,理论上、计算 上以及应用上的许多问题就容易处理了, 当然花费也大了。