第四章 矩阵的特征值与特征向量问题
计算方法(5)第四章 矩阵特征值和特征向量的计算
2 ,比值越小,收敛越快。 1
几点说明: 1)如果u0的选取恰恰使得1 0, 乘幂法计算仍能进行。 因为计算过程中舍入误差的影响,迭代若干次后,必然会 产生一个向量uk , 它在x1方向上的分量不为零,可以说实际 中出现1 0的可能性几乎为零。
2)因uk 1k1 x1 , 计算过程中可能会出现溢出( 1 1) 或成为0( 1 1)的情形。解决方法:每次迭代所求的向量 都要规范化。因此,乘幂法实际使用的计算公式是
n
uk j i k 由lim( ) 0 1 = lim ,(j 1, 2, k k u 1 k 1 j
n)
uk j 1 uk 1 j
按上面式子计算矩阵A 按模最大的特征值与相应的 特征向量的方法称为乘幂法。 乘幂法的收敛速度依赖 于比值
/1 ,当比值接近于1时,乘幂法收敛很慢。 比值 2
乘幂法加速有多种,重点介绍原点平移法。
矩阵A与A pI 的特征值有以下关系:若i 是A 的特 征值,则i p就是A pI 的特征值,而且相应的特征向 量不变。如果用矩阵A pI 按uk Auk 1计算,则有 uk ( A pI )uk 1
a22 a11 2 1 t t 0 a12
解二次方程可以求出t , 容易得到 c 1 1 t
2Байду номын сангаас
s ct
(1) a 0 c s T T 11 A1 RAR , R (1) s c 0 a 22 R T的两个列向量是相应的特征向量。
再将A 2中(2,3)和(3,2)位置上的元素变成0, 0 0 1 取R 2 0 0.975 0.226 0 0.226 0.975
线性代数第四章矩阵的特征值
令 P ( p1 p2 L pn ), 则P 可逆,且
AP ( Ap1 Ap2 L Apn ) (1 p1 2 p2 L n pn )
1
( p1 p2 L
pn
)
2
O
P,
n
2. 求出矩阵A对应于所有特征值的特征向量
若A有一个t重特征值,对应的特征向量在线性 无关的意义下小于t,则A不与对角矩阵相似。
3.写出对角矩阵和相似变换矩阵。 特征值和特征向量的对应.
1. 求出n阶矩阵A的所有特征值 2. 求出矩阵A对应于所有特征值的特征向量 3.写出对角矩阵和相似变换矩阵。
3 1
的λ都是方阵A的特征值.
定义4.2 设A为n阶矩阵,含有未知量λ的矩阵λI-A
称为A的特征矩阵,其行列式
I A
为λ的n次多项式,称为A的特征多项式, I A 0
称为A的特征方程.
求n阶矩阵的特征值和特征向量的步骤:
1. 由矩阵A的特征方程 I A 0 求出A的特征值 1,2 ,L s (s n 2k )
所以 P 1 AP , 即A与对角矩阵Λ相似.
定理的证明告诉我们,如果n阶矩阵A与对角矩 阵Λ相似,则Λ的主对角线上的元素就是A的全部
特征值.相似矩阵P的列是对应于Λ对角线上 元素的特征向量。
推论 若n阶矩阵A有n个两两不同的特征值,则
A必与对角矩阵Λ相似
推论 若n阶矩阵A有n个特征值,则可相似对 角化<==>A的任ti重特征值有对应ti个线性无
A
4 1
3 0
0 2
4第四章 矩阵的特征值与特征向量
THANKS
演示完毕 感谢观看
定义9 如果 n 阶矩阵 A 满足 , AT A I(即A1 AT) 则称 A 为正交矩阵,简称正交阵.
第四章
4.3
相似矩阵与矩阵的对角化
一、相似矩阵
定理2
n 阶矩阵 A 相似于对角阵(即矩阵 A 可对角化)的充分必要条件是矩阵 A 有 n 个线
性无关的特征向量. 推论 若 n 阶矩阵 A 有 n 个不同的特征值,则矩阵 A 与对角阵相似. 证 设 A 有 n 个不同的特征值 1, 2 , , n ,对应的特征向量 1 , 2 , , n . 根据§4.1定理3,向量组 1 , 2 , , n 线性无关,所以由定理2可知, A 与对角阵 相似,其中
线性代数
主编 金桂堂
普通高等教育”十三五“规划教材
肖马成
矩阵的特征值 与特征向量
第四章
本章主要讨论方阵的特征值与特征向量、方阵的相 似对角化等问题,然后介绍向量空间、基与维数, 以及向量的内积、长度及正交等知识.
第四章
4.2
n维向量空间
一、向量空间与子空间
第四章
ห้องสมุดไป่ตู้
4.2
n维向量空间
四、正交矩阵
第四章 矩阵的特征值和特征向量
第四章 矩阵的特征值和特征向量例1 求下列矩阵的特征值与特征向量⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=163053064A ,并判断它能否相似对角化。
若能,求可逆阵P ,使∧=-AP P 1(对角阵)。
例2 已知三阶方阵A 的三个特征值为4,3,2-,则1-A 的特征值为_______,TA 的特征值为_______,*A 的特征值为_______,E A A 232+-的特征值为_______例3 设矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=0011100y x A 有三个线性无关的特征向量,则y x ,应满足条件_______ 例5 已知矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=x A 10200002与⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=10000002y B 相似,则____________==y x 例6 设n 阶方阵A 满足0232=+-I A A ,求A 的特征值例7 已知向量T k )1,,1(=ξ是矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=211121112A 的逆矩阵1-A 的特征向量,求常数k例8 设A 为非零方阵,且0=mA (m 为某自然数),证明:A 不能与对角阵相似 例9 设n 阶方阵A 满足01072=+-I A A ,求证:A 相似于一个对角矩阵结论 总结1 n 阶方阵A 有n 个特征值,它们的和等于A 的主对角线元素之和(即A 的逆trA ),它们的乘积等于A 的行列式A2 如果mλλ,,1 是方阵A 的特征值,m P P ,,1 是与之对应的特征向量,如mλλ,,1 互不相等时,m P P ,,1 线性无关3 如果n 阶方阵A 与B 相似,则A 与B 有相同的特征多项式,从而有相同的特征值4 如果n 阶方阵A 与对角阵∧相似,则∧的主对角线元素就是A 的n 个特征值5 n 阶方阵A 与对角阵∧相似,即A 可相似对角化的充要条件是A 有n 个线性无关的特征向量6 如果n 阶方阵A 的n 个特征值互不相等,则A 与对角阵相似,即A 可相似对角化7 实对称矩阵的特征值全为实数8 实对称矩阵的不同特征值对应的特征向量相互正交9 对实对称矩阵n n A A ⨯=,必存在正交矩阵P ,使∧=-AP P 1,其中∧是以A 的n 个特征值为主对角线元素的对角阵10 方阵A 可逆的充要条件是A 的特征值全不为零习 题一 填空题1 设A 为3阶矩阵,其特征值为2,1,3-,则A =________ 1-A 的特征值为________,E A A +-322的特征值为________2 如果二阶矩阵⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡=4231,127B x y A 相似,则 __________==y x 3 若n 阶可逆阵A 的每行元素之和是)0(≠a a ,则数________一定是E A +-12的特征值4 设三阶矩阵A 有3个属于特征值λ的线性无关的特征向量,则______=A5 若E A =2,则A 的特征值为________6 设n 阶方阵A 的n 个特征值为n ,,2,1 ,则_______=+I A7 设⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=101120101A 2≥n ,则_______21=--n n A A 8 ⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡-=61000512141A 则 ______lim =∞→n n A二 选择题1 设矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=100321z y x A A 的特征值为3,2,1,则( ) A )8,4,2===z y x B) R z y x ∈==,4,1 C) R z y x ∈=-=,2,2 D) 3,4,1===z y x2 已知矩阵⎪⎪⎭⎫⎝⎛-x 123022有一特征向量⎪⎪⎭⎫⎝⎛-35,则)(=xA) 18- B) 16- C) 14- D) 12-3 设⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=53342111x A A 有特征值2,621==λλ (二重 ),且A 有三个线性无关的特征向量,则______=xA) 2 B) 2- C) 4 D) 4-4 若B A ~(等价 ),则有( )A )B I A I -=-λλ B) B A =C) 对于λ,矩阵A 与B 有相同的特征值与特征向量 D) A 与B 均与一对角矩阵相似 5 已知矩阵A 的各列元素之和为3,则( )A) A 有一个特征值为3,并对应一个特征向量T )1,,1,1( B) A 有一个特征值为3,并不一定对应有特征向量T )1,,1,1( C) 3不一定是A 的特征值 D) A 是否有特征值不能确定 6 设A 是三阶矩阵,有特征值2,1,1-,则下列矩阵中可逆的是( ) A) A I - B) A I + C) A I -2 D) A I +2 三 解答题1. 设三阶矩阵A 的特征值为3,2,1321===λλλ,对应的特征向量依次为:T )1,1,1(1=ξ,T )4,2,1(2=ξ,T )2,3,1(3=ξ,又向量T )3,1,1(=β1) 将β 用321,,ξξξ线性表示 2) 求βnA(n 为自然数)2 已知⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=06303012x A 有3个线性无关的特征向量,求100A3. 设⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=122212221A 求A 的特征值与对应的特征向量,A 是否对角阵相似。
(线性代数)第四章 矩阵的特征值和特征向量
∴η1
=
a2 − a1
1 0 0 0 ,η 2 = 1 ,L ,η n −1 = 0 M M M 0 1 0
对应λ=0的 =0的 特征向量为 k1η1 + L + kn −1η n −1 , k ,L , k 不全 n −1 1
第四章 矩阵的特征值和特征向量
§4.1 相似矩阵
§4.1 相似矩阵 一. 相似矩阵的定义和性质 AP= 都是n阶方阵, 若有可逆矩阵P 设A, B都是n阶方阵, 若有可逆矩阵P, 使P−1AP=B, 则称矩阵A 相似. 记为A 相似变换矩阵. 则称矩阵A与B相似. 记为A~B. P为相似变换矩阵. 相似是相抵的特例 相似必相抵,反之不然. 特例: 注1: 相似是相抵的特例: 相似必相抵,反之不然. 注2: 矩阵间的相似关系是一种等价关系 (1) 反身性: A~A; 反身性: P−1AP =B (2) 对称性: A~B ⇒ B~A; 对称性: PBP−1 =A (3) 传递性: A~B, B~C ⇒ A~C. 传递性: 相抵关系下的不变量: 相抵关系下的不变量:矩阵的秩 相似关系下的不变量: 相似关系下的不变量: 矩阵的秩
第四章 矩阵的特征值和特征向量
§4.2 特征值与特征向量
解: |λE–A| = (λ+1)(λ –2)2. +1)( 所以A 所以A的特征值为λ1= –1, λ2= λ3= 2. (–E–A)x = 0的基础解系: ξ1=(1,0,1)T. 的基础解系: 对应于λ1= –1的特征向量为kξ1 (0≠k∈R). 的特征向量为k (0≠ (2E–A)x = 0的基础解系: (2E 的基础解系: ξ2=(0, 1, –1)T, ξ3=(1, 0, 4)T. =2的特征向量为 的特征向量为k 对应于λ2=λ3 =2的特征向量为k2ξ2 +k3ξ3 (k2, k3不同时为零). 不同时为零).
分析04-矩阵的特征值
矩阵的特征值与特征向量的计算
产生迭代向量序列
(由x 的某一分量的相邻二次结果之比 可得出1),而相应的特征向量为 x ( k 1) 。
实际上, 由式(4-1)可得 :
x ( k 1) Ax( k ) Ak 1 x ( 0) Ak 1 ( i ui )
i 1 k i k 1ui 11 1u1 2k 1u2 n k 1un i 2 n i 1 n n
当n不大时,如n4 解特征方程,可求出全部特征值 (n 3较难)当 n较大(n>5),计算量会增大得惊人, 且不可能求得准确结果,还可能出现不稳定,所以当n稍 大一般不直接求解特征方程,而根据实际问题的需要,介 绍相应的一些行之有效的数值解法
第四章 矩阵的特征值与特征向量的计算
4-4
W
矩阵的特征值与特征向量的计算概述(续1)
x(k+1)为1对应的特征向量收敛到1u1+…+mum
W Y
两 点 注 释(续2)
( k 1) xi ( xi k )
x
( k 1)
( k 1) 1
x
( 0)
可构造向量序列
所以乘幂法实际上是,对于给定的初始向量 x ( 0) ( 零向量)由迭代法:
x
( k 1)
第四章
W Y
u
i 1
n
i i
1u1 2u 2 nU n
x
( k 1)
Ax
(k )
Ax
(k )
(k 0,1,, )
(4 -1)
5. 若 , 输出 , x, 停机; 否则, 转6
6. 若k<N,置k+1k, ,转3;否则, 输出失败信息,停机。
概率论与数理统计4-1矩阵的特征值与特征向量
k 2 p2 k 3 p3
( k 2 , k 3 不同时为 ). 0
求矩阵特征值与特征向量的步骤:
1. 计算A的特征多项式 | A E | 或 | E A | ;
2. 求特征方程 | A E | 0 或 | E A | 0 的全部根
1 , 2 , , n , 就是A的全部特征值 ;
1 0 1 ~ 0 1 0 , 0 0 0
故对应于1 1的全体特征向量为 k p1 ( k 0).
当2 3 2时, 解方程 A 2 E x 0.由
4 1 1 4 1 1 A 2 E 0 0 0 ~ 0 0 0 , 4 1 1 0 0 0 得基础解系为: 0 1 p2 1 , p3 0 , 1 4 所以对应于 2 3 2的全部特征向量为:
推广
. 是A 的特征值
m
m
例3 设λ是方阵A的特征值, 证明
2 是 A 2 的特征值; (1) 1 是A 1的特征值. (2) 当A可逆时,
m m . 是A 的特征值
2 当A可逆时, 0,
1
1
由Ax x可得
1
A Ax A x A x
A x x
解
2 A E 0 4
1 2 1
2
1 0 3
( 1) 2 , 2 令 ( 1) 2 0
得A的特征值为1 1, 2 3 2.
当1 1时, 解方程 A E x 0.由
1 1 1 A E 0 3 0 4 1 4 1 得基础解系 p1 0 , 1
《线性代数》第四章第二节 方阵的特征值与特征向量
若P是与对应的特征向量,则显然k 0时, kP也是与对应的特征向量.
6.属于同一特征值的特征向量的非零线性组合 仍是属于这个特征值的特征向量.
例1
设
A
=
−2 0
1 2
1 0,
求A的特征值与特征向量.
− 4 1 3
分析:
1.特征方程的根就是特征值;
2. (A-E)x=0的通解(去掉零解)就是特征值对应
所以对应于 2 = 3 = 2的全部特征向量为 :
k2 p2 + k3 p3 (k2 , k3不同时为0).
例2 证明:若 是矩阵A的特征值,x 是A的属于 的特征向量,则
(1) m是Am的特征值(m是任意常数).
(2) 当A可逆时,−1是A−1的特征值.
证明 (1) Ax = x A(Ax) = A(x) = (Ax) = (x) A2 x = 2 x
有x.
3. A − E = 0 为A的特征方程。
a11 −
a21
an1
a12
a22 −
an2
a1n
a2n
=0
ann −
记 f ( ) = A − E ,它是的n次多项式, 称其
为方阵A的 特征多项式 .
( ) 4. 设 n阶方阵A = aij 的特征值为1, 2 ,,
n ,则有 (1) 1 + 2 + + n = a11 + a22 + + ann; (2) 12 n = A .
将1 = 2 = 1代入(A − E )x = 0,
解之得基础解系
− 2
1 = 1 ,
0
第四章 特征值与特征向量(研究生矩阵论)
A( A x 3Ax 4x) 0, ( A 4E)(A x Ax) 0, ( A E)(A x 4 Ax) 0.
线性无关,所以
2
由于
x, Ax, A x
2
A x 3 Ax 4 x 0,
A x Ax 0, A x 4 Ax 0. 据此立得: A 0, A 4E 0, A E 0. 故 A 的三个特征值为0,-1,4. 于是
A
有特征值
1
以及相应的特征向量 的一组基,设为
n
1. 将 1 扩充成 C n 1, 2 ,, n . 令
则因 A1 11 , A i b ji j (2 i n),
j 1
1 b12 b1n 0 b22 b2 n A1 , A 2 ,, A n 1 , 2 ,, n 0 b b n2 nn 若令 P 1 , 2 ,, n ,
【命题4.2.3】设
n 阶矩阵
A 相似于对角矩阵
f ( x)
n 个特征值为 1, 2 ,, n, 是一多项式,则 f ( A) 的 n 个特征值为
A
的
f (1 ), f (2 ),, f (n ).
【推论4.2.1】设 则对
f ( x)
是一多项式,若
f ( A) 0,
A
的任一特征值
a11 0 0 a11 0 a13 a11 a12 0 a12 0 a21 a23 a21 a22 0 a13 0 a31 0 a33 a31 a32
a12 a13 a11 a12 a13 0 a22 a23 a21 a22 a23 0 a23 a33 a31 a32 a33
数值计算方法第04章矩阵特征值与特征向量的计算
• 计算出k=2时的x和y。 • (保留四位有效数字)
22
二、幂法的加速
因为幂法的收敛速度是线性的,而且依赖 于比值 2 /1 ,当比值接近于1时,幂法收敛 很慢。幂法加速有多种,介绍两种。
23
幂法的加速—原点移位法 应用幂法计算矩阵A的主特征值的收敛速度主要
26
4 14 0 , 2.9, 用原点移位法求矩 例:A 5 13 0 0 1 0 2.8 -4 阵A的按模最大的特征值,要求误差不超过10 。 解:取x (0) (1,1,1)T , 按x ( k 1) ( A pI )x (k )进行计算 0 6.9 14 A 0 I 5 10.1 0 0 0.1 1 (3.1000568, 2.214326, 0.9687661) 4 3.1000568
在一定条件下, 当k充分大时: 相应的特征向量为:
x 1 x
x
( k 1)
( k 1 ) i (k ) i
10
幂法的理论依据 对任意向量x(0), 有 x ( 0 ) i ui , 设1不为零.
i 1 n
x
( k 1 )
Ax
n i 1
(k )
A
k 1
x
(0) n
1 Ak 1 i ui i k i ui i 1
k 1 1
2 k 1 n k 1 1u1 ( ) 2 u2 ( ) n un 1 1
k 1 1 1u1
故 1 xi( k 1) xi( k ) x(k+1)为1的特征向量的近似向量(除一个因子外).
第四章 方阵的特征值和特征向量
4.1.3 反幂法 由Axi=ixi易推得A-1xi=(1/i)xi ,若有
| 1 | | 2 | | 3 | | n |,
则1/n是A-1的按模最大的特征值,我们只要求出A-1的按模最大的 特征值,也就求出了A的按模最小的特征值.为了避免求逆阵,我们 用解方程组的方法构造如下算法:
5 结束
2. 我们假设在(4.3)中α1≠0,这在选择u0时,也无法判断,但这往往不 影响幂法的成功使用.因为若选u0,使α1=0,由于舍入误差的影响, 在迭代某一步会产生uk,它在x1方向上的分量不为零,这时以后的 迭代仍会收敛. 3. 我们假设了 | 1 | | 2 | | 3 | | n |,
u k 1 1 x1
k
2 1
1
Au k
k 1 1
1 x1 1 1 x1 1 u k ,
k 1
不是零向量,
即uk为1的近似的特征向量. 2 结束
实际计算时,为防止uk的模过大或过小,以致产生计算机运算的 上下溢出,通常每次迭代都对uk进行归一化,使‖ uk ‖∞=1,因此 以上幂法公式改进为:
y k 1 u k 1 u k Ay k 1 u k 1
k 1, 2 ,
( 4 .4 )
此时uk仍收敛于1对应的特征向量。1可用如下公式计算:
1
ak a k 1 (4 .7 )
其中ak 是uk 的绝对值最大的分量,a k 1 是yk-1 的绝对值最大 的分量。
cos 为实对称阵, U sin sin cos
1. 二阶实对称矩阵的对角化
设
为二阶旋转矩阵,容易验证U正交。 12 结束
第四章-矩阵的特征值与特征向量问题讲解
Ax 2 x
1 2 x 0,
则x 0, 与定义矛盾.
12
注记
4. 若λ是矩阵A的r重特征值,对应λ有s个线性 无关的特征向量,则1≤s≤r; 若A为实对称矩阵,则对应特征值λ 恰有r 个线性无 关的特征向量。
5. 实对称矩阵的特征值是实数,属于不同特 征值的特征向量正交。
13
注记
6. 设 n阶方阵 A aij 的特征值为1, 2 ,, n ,记:
定义:设A是n阶方阵, 是一复数,如果方程 Ax x
存在非零解向量,则称 为方阵A的特征值, 相应的非零解向量x 称为与特征值 对应的特征向量, 此特征值与特征向量x称为一特征对, P(A )=det(I A)称为矩阵A的特征多项式。
4
注记
1. 特征向量x 0, 特征值问题是对方阵而 言的. 2. n阶方阵A的特征值,就是使齐次线性方程组
0.2 0.3 0.1 4
G1 = {z:|z – 1| 0.6};G2 = {z:|z – 3| 0.8}; G3 = {z:|z + 1| 1.8};G4 = {z:|z + 4| 0.6}。
G4
G1
G2
G3
注:定理推断A的n个特征值全落在n个盖氏圆
上,但未说明每个圆盘内都有一个特征值。
20
对应的特征值1,2,…,n,满足
|1| > |2| … |n|
(4.1.1)
26
1.基本思想
因为{v1,v2,…,vn}为Cn的一组基,故:
任给x(0) 0,
n
x (0) aivi
所以有:
i 1
n
n
Ak x(0) Ak ( aivi ) ai Akvi
矩阵的特征值与特征向量
1
所以,A 的特征值为 1 2 , 2 3 1,
7
当 1 2 时, 解方程组 ( 2 I A ) x 0 ,
1 x1 2 1 x 2 0, 即 1 2 x 3 1 解之得基础解系为 p 1 1 , 1 2 1 1 1
故 1是 A1的特征值, 且 x 也是 A1对应于1的特征向量.
24
性质2 矩阵 A 和 AT 的特征值相同. 证 因为 IAT = ( I)TAT = ( IA)T 所以 det ( IA) = det ( IAT)
因此, A 和AT 有完全相同的特征值.
补充 性质 设 是方阵 A 的特征值.设
(*)式中不含 的常数项为
a 11 a 21 a n1 a 12 a 22 an2 a1n a2n a nn
21
( 1) A ,
n
即 c n ( 1) A
n
f 所以, ( ) I A ( 1 )( 2 ) ( n )
的全部特征向量.
9
例
2 设矩阵 A 2 0
2 1 2
0 2 , 求 A 的特征值. 0
解 A 的特征多项式为
2 I A
2 0 2 0 2 ( 2 )( 1 )( 4 ),
1
2
所以,A 的特征值为 1 2 , 2 1 , 3 4 . 特征值的计算不容易!!
0 A 1 1 1 0 1 1 1 0
所以 k 1 p 1 是对应于 1 2 的全部特征向量;
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
则显然有A(1) = A,A(0) = D, 易知A()的特征多项式的系数是的多项式, 从而A()的特征值1(),2(),…,n()为的连续函 数。
22
G ( ) { z : z a | a | | a | } G ,( 0 1 ) i ii ij ij i
1 1
1
1
1
的特征向量 .
★ 特征值和特征向量的性质:
A或 AT
特征值
kA k
x
对应特征 向量 x
x
A m x
m
A
1
A
*
f ( A)
f ()
1/
A/
x
x
x
9
,x 依次是与之对应的特 向量。如果 , , , m 1 2 m 各不相等 , 则 x ,x , ,x 线性无关 。 1 2 m
式 ,当各 ,该行列式不等于 0 ,从而该矩阵 i不相等时 k x , k x , , k x 0 , 0 , , 0 , 可逆 .于是有 1 1 2 2 m m
但 x 0 ,故 即 k x 0 j 1 , 2 , , m . k 0 j 1 , 2 , , m . j j j j
第四章 矩阵的特征值与 特征向量问题
1
第三章 矩阵的特征值与特征向量
4.1 幂法与反幂法 4.2 Jacobi方法(重点) 4.3 多项式方法求特征值问题(自学) 4.4 QR算法 (重点)
Givens矩阵; Householder矩阵; Gram-Schmidt正交化方法
2
3
a kx jj
j k
Gk Gi
i 1
19
n
例 估计方阵A特征值的范围
0.1 0.2 1 0.5 3 0.1 A 1 0.3 1 0.2 0.3 0.1 0.3 0.2 0.5 4
解:
G1 = {z:|z – 1| 0.6};G2 = {z:|z – 3| 0.8}; G3 = {z:|z + 1| 1.8};G4 = {z:|z + 4| 0.6}。
为S,它与后n – k个圆盘严格分离,显然,A()的前k
个盖氏圆盘与后n – k个圆盘严格分离。 当 = 0时,A(0) = D的前k个特征值刚好落在前k个圆 盘G1,…,Gk中,而另n – k个特征值则在区域S之外,
G4 G3
G1
G2
注:定理推断A的n个特征值全落在n个盖氏圆 上,但未说明每个圆盘内都有一个特征值。
20
盖氏圆的连通部分
称相交盖氏圆之并构成的连通部分为连通部分。 孤立的盖氏圆本身也为一个连通部分。
定理 若由A的k个盖氏圆组成的连通部分,含且 仅含A的k个特征值。
21
定理 若由A的k个盖氏圆组成的连通部分,含且仅含A 的k个特征值。 证明: 令D = diag(a11,a22,…,ann),M = A – D,记
的特征向量 , 即有 1 2
因为 , 如果设 x 同时是 A 的属于特征 , 的 1 2
Ax x , Ax x 1 2 x x x 0 , 1 2 1 2 由于 0 , 则 x0 , 与定义矛盾 . 1 2
12
概述
定义: 设 A是 n阶方阵 , 是一复数,如果方程 Ax x 存在非零解向量 ,则称 为方阵 A的 特征值 , 相应的非零解向量 x 称为与特征值 对应的 特征向量 , 此特征值 与特征向量 x称为一 特征对 , P ( )= det (I A)称为矩阵 A的 特征多项式 。 是 A 的特征值 .
m m ( 1 ) 是 A 的特征值 m 是任意常 .
证明: 1 Ax x
2 2 A x x A Ax A x Ax x
m m m 2 A x x 再继续施行上述步骤 次,就得
3 . 是矩阵 A 的一个特征值,则一定 是 A λI0, 的根,因此又称 特征根 。若 是 A I 0 的 k 重根 , 则称 为 A 的 k 重特征值 ( 根 ) 。
5
重数:
det( I A ) ( 1 ) ( 2 ) ( P )
k x k x k x 0 .
k k 1 1 1 2 2 2
k x k x k x 0 , 1 1 1 2 2 2 m m m
k m m m
10 k 1 , 2 , , m 1
把上列各式合写成矩阵形式,得 m 1 1 1 1 m 1 1 2 2 k x , k x , , k x 0 , 0 , , 0 1 1 2 2 m m m 1 1 m m 上式等号左端第二个矩 阵的行列式为范德蒙 列
a 0 a a 11 12 1 n a a 0 a 21 22 2 n A ( ) D M ( 0 1 ) a a 0 n n n 1 a n 2
15
Jordan分解定理
设 A C n n , 有 r 个互不相同的特征值 其重数分别为 矩阵 P C
n n 1
( i i 1, , r )
n( i i 1 , , r ),则一定存在非奇异 使得 P AP diag ( J ( 1 ) J ( r )) J 1 i
n1 n2 nP
其中:n1 n 2 n p n ; i j (i j ) 称 n i 为 i的 代数重数 (简称 重数); m i n rank ( i I A )为 i的 几何重数 ( m i n i ); 如果 n i=1,则称 i 为 A 的一个 单特征值 ; 否则称 i 为 A 的一个 重特征值 ; 如果 m i= n i,则称该特征值 i 为 A 的一个 半单特征值 ; 如果 A 的所有特征值都是半单 的,则称 A 是 非亏损 的; A 是非亏损的充要条件是 A 有 n 个线性无关的特征向量 。
m m m m 故 是矩阵 A 的特征值 , 且 x 是 A 对应于 的
征向量 .
8
2 当 A 可逆时 , 0 ,由 Ax x 可得
1 1 1 A Ax A x A x
x A x A x x
1 1
故 是矩阵 A 的特征值 , 且 x 是 A 对应
盖氏圆
设A = [aij]nn,称由不等式
z a ii a ij
j 1 j i
n
所确定的复区域为A的第i个盖氏圆,记为Gi:
i = 1,2,…,n。 G { z : z a a } i ii ij
j 1 j i
18
n
Gerschgorin圆盘定理
定理 若为A的特征值,则 G i
Schur分解定理
n n n n 设 A C ,则 存 在 酉 矩 阵 U C 使 得 H U A U T
其 中是 T 上 三 角 矩 阵 , 且 选 择 适 当 矩 阵 U ,可 使 T 的 主 对 角 线 元 素 按 任 意 指 定 次 序 排 列 。
17
特征值估计
粗略估计有(A) ||A||; 可将复平面上的特征值一个个用圆盘围起。
注记
4. 若λ是矩阵A的r重特征值,对应λ有s个线性 无关的特征向量,则1≤s≤r; 若A为实对称矩阵,则对应特征值λ 恰有r 个线性无 关的特征向量。
5.
实对称矩阵的特征值是实数,属于不同特 征值的特征向量正交。
13
注记
n
6 .设 n 阶方阵 A a 的特征值为 , i j 1, 2, n,记:
I A( ( ( 1 )s 1)( 2) n) k
k k 0 n n 1 n 2 n s s ( 1 ) A 1 2
n k
其中 :s A 中一切 k 阶主子式之和。 k为
( 1 ) a a a = tr ( A ) ; 1 2 n 11 22 nn
4
注记
1 . 特征向量 x 0 ,特征值问题是对方阵而 言的 .
2 .n 阶 方 阵 A 的 特 征 值 ,就 是 使 齐 次 线 性 方 程 组 Ix 0 有 非 零 解 的 值 ,即 满 足 方 程 A A I 0 的 都 是 矩 阵 A 的 特 征 值 .
11
所以向量组 x , x , , x 线性无关 . 1 2 m
注记
1. 2.
属于不同特征值的特征向量是线性无关的.
属于同一特征值的特征向量的非零线性组合 仍是属于这个特征值的特征向量. 3. 矩阵的特征向量总是相对于矩阵的特征值而 言的,一个特征值具有的特征向量不唯一; 一个特征向量不能属于不同的特征值.
i 1 n
证明:设Ax = x (x 0),若k使得 x max x k i x
1 i n
因为
a
j 1
n
kj
x j xk
n
( a )x a xj kk k kj
n |x | j a | | | a | |a | kk kj kj |x j k x j 1 j 1 k k| n n j k j k
, , , 是方阵 A 的 m 个特征值 ,x ,x , 定理 : 设 1 2 m 1 2
证明: 设有常数 k , k , , k 使 1 2 m