矩阵特征值问题
(完整版)线性代数第五章特征值、特征向量试题及答案
第五章 特征值和特征向量一、特征值与特征向量定义1:设A 是n 阶矩阵,λ为一个数,若存在非零向量α,使λαα=A ,则称数λ为矩阵A 的特征值,非零向量α为矩阵A 的对应于特征值λ的特征向量。
定义2:()E A f λλ-=,称为矩阵A 的特征多项式,)(λf =0E A λ-=,称为矩阵A 的特征方程,特征方程的根称为矩阵A 的特征根 矩阵E A λ-称为矩阵A 的特征矩阵齐次方程组(0)=-X E A λ称为矩阵A 的特征方程组。
性质1:对等式λαα=A 作恒等变形,得(0)=-αλE A ,于是特征向量α是齐次方程组(0)=-X E A λ的非零解向量,由齐次线性方程组有非零解的充要条件知其系数行列式为零,即0=-E A λ,说明A 的特征值λ为0E A λ-=的根。
由此得到对特征向量和特征值的另一种认识:(1)λ是A 的特征值⇔0=-E A λ,即(λE -A )不可逆.(2)α是属于λ的特征向量⇔α是齐次方程组(0)=-X E A λ的非零解.计算特征值和特征向量的具体步骤为: (1)计算A 的特征多项式,()E A f λλ-=(2)求特征方程)(λf =0E A λ-=的全部根,他们就是A 的全部特征值;(3)然后对每个特征值λ,求齐次方程组(0)=-X E A λ的非零解,即属于λ的特征向量.性质2:n 阶矩阵A 的相异特征值m λλλ 21,所对应的特征向量21,ξξ……ξ线性无关性质3:设λ1,λ2,…,λn 是A 的全体特征值,则从特征多项式的结构可得到:(1)λ1+λ2+…+λ n =tr(A )( A 的迹数,即主对角线上元素之和). (2)λ1λ2…λn =|A |.性质4:如果λ是A 的特征值,则(1)f(λ)是A 的多项式f(A )的特征值.(2)如果A 可逆,则1/λ是A -1的特征值; |A |/λ是A *的特征值. 即: 如果A 的特征值是λ1,λ2,…,λn ,则 (1)f(A )的特征值是f(λ1),f(λ2),…,f(λn ).(2)如果A 可逆,则A -1的特征值是1/λ1,1/λ2,…,1/λn ; 因为A AA =*,A *的特征值是|A |/λ1,|A |/λ2,…,|A |/λn .性质5:如果α是A 的特征向量,特征值为λ,即λαα=A 则(1)α也是A 的任何多项式f(A )的特征向量,特征值为f(λ);(2)如果A 可逆,则α也是A -1的特征向量,特征值为1/λ;α也是A *的特征向量,特征值为|A |/λ 。
矩阵特征值问题的数值方法.
矩阵特征值问题的数值方法矩阵特征值设A 是n 阶矩阵,x 是非零列向量. 如果有数λ 存在,满足那么,称x 是矩阵A 关于特征值λ的特征向量. 很显然一般地有主特征值的乘幂迭代法设n 阶矩阵A 的n 个特征值按模从大到小排序为:n 其对应的n 个线性无关的特征向量分别为:设是任意一个非零的n 维向量,则:假设,构造一个向量序列:则:或者:当时:如果是矩阵A 的关于特征值的一个特征向量,特征值个特征那么对于任意一个给定的,也是特征值的特征向量。
所以,是对主特征值对应的特征向量的近似。
如果则会变得很大或者如果,则会变得很大,或者如果,则会变得非常小,在实际计算中,为避免这种情况的出现需对做归一化处理况的出现,需对做归一化处理:由:左乘得:所以主特征值的近似值所以主特征值的近似值:残余误差向量定义为:当迭代次数充分大时,残余误差将充分小。
逆乘幂法:类似地,也可以求模最小特征值和对应的特征向量特征向量。
上述问题的主特征值问题就是矩阵A 的模最小特征值问题。
结果,逆乘幂法的迭代公式为:在实际应用中,无需计算逆矩阵,但需求解线性系统实对称矩阵的基本定理:对实对称矩阵A ,一定存在一个正交相似变换使得为对角矩阵且其对角矩阵P ,使得:为对角矩阵,且其对角的特征值元素为矩阵A 的特征值。
相似变换:相似变换保持矩阵特征值(但不是特征向量)不变不变。
(证明略)正交相似变换:中。
正交相似变换的例子—坐标旋转:叫旋转矩阵。
容易验证:。
适当选择旋转角,可消去xy 项—得到对角阵D 。
矩阵特征值问题的数值方法实对称矩阵的基本定理再看下面的例子:令:O 平面的坐标旋转变换适当同样地有:。
则是在x-O-z 平面的坐标旋转变换。
适当x z —D 。
选择旋转角可消去z 项得到对角阵实对称矩阵的Jacobi 方法:全部特征值和特征向量根据实对称矩阵的基本定理,求得矩阵A 的全部特征值的关键是找到正交相似变换矩阵P 使部特征值的关键,是找到正交相似变换矩阵P ,使得为对角阵。
矩阵的特征值问题求解
矩阵特征值问题求解矩阵在数学和工程领域有着广泛的应用,而研究矩阵的特征值是其中一个重要的问题。
矩阵的特征值对于矩阵的性质和行为具有重要的影响,因此求解矩阵的特征值是一项非常重要的任务。
什么是特征值和特征向量在矩阵理论中,矩阵A的特征值(eigenvalue)是一个数λ,满足方程$A\\mathbf{v} = \\lambda\\mathbf{v}$的向量$\\mathbf{v}$存在且不为零。
其中,$\\mathbf{v}$被称为对应于特征值$\\lambda$的特征向量(eigenvector)。
特征值和特征向量的求解是矩阵理论和线性代数中的重要问题之一。
特征值问题的求解方法1. 特征值分解我们可以通过特征值分解的方法求解矩阵的特征值。
给定一个方阵A,我们可以将其表示为$A=Q\\Lambda Q^{-1}$的形式,其中Q是由A的特征向量所组成的矩阵,Λ是由A的特征值所组成的对角矩阵。
2. 特征多项式特征值问题的另一种求解方法是通过矩阵的特征多项式。
特征多项式是关于矩阵A的一个多项式,它的根就是矩阵A的特征值。
通过求解特征多项式的根,我们可以得到矩阵的特征值。
3. 幂法幂法是一种常用的求解特征值问题的迭代方法。
通过不断的迭代计算$A\\mathbf{v}^{(k)}$,其中$\\mathbf{v}^{(k)}$是第k次迭代得到的特征向量,我们可以逐渐逼近矩阵的特征值和特征向量。
应用和意义矩阵的特征值问题求解在计算机图形学、信号处理、物理学等领域都有着重要的应用和意义。
通过求解矩阵的特征值,我们可以分析矩阵的性质、系统的稳定性以及模式识别等问题,为我们深入理解和应用矩阵提供了重要的工具和方法。
综上所述,矩阵的特征值问题求解是一个具有重要意义和广泛应用的问题,通过不同的方法和技术,我们可以有效地求解矩阵的特征值和特征向量,为我们更好地理解和利用矩阵提供了重要的支持。
数值分析-第7章 矩阵特征值问题的数值解法n
7
9 11 12
6.104716
6.026349 6.006637 6.003327
(-0.450275, -0.322058, 1.0)
(-0.445914, -0.318617, 1.0) (-0.444814, -0.31775, 1.0) (-0.444630, -0.317606, 1.0)
其中i为A的特征值,P的各列为相应于i的特征向量。
P -1 AP D
2
n
2
定理7.1.3 ARnn,1, …, n为A的特征值,则
(1)A的迹数等于特征值之和,即 tr ( A) aii i
i 1 i 1
n
n
(2)A的行列式值等于全体特征值之积,即
1 xi(k +1) / xi(k )
i 1,2,, n
可见,当k充分大时, ( k ) 近似于主特征值, ( k +1) 与x ( k )的对应非零分量的比值 x x 近似于主特征值。
在实际计算中需要对计算结果进行规 , 范化。因为当 1 1时,x (k ) 趋于零, 当1 1时, x ( k )的非零分量趋于无穷。 从而计算时会出现下溢 或上溢。
特征值的范围. 解 我们先分别求出各个圆盘区域。 D1 = {z:|z – 1|£0.6};D2 = {z:|z – 3|£0.8} D3 = {z:|z + 1|£1.8};D4 = {z:|z + 4|£0.6}. 易见D2和D4为 弧立圆盘分别 包含A的两个实 特征值.
第七章—矩阵的特征值问题
1 p仍是B的主特征值,且使 2 p 2 1 p 1
对B应用幂法,使得在计算B的主特征值1 p的过程中 得到加速。这种方法通常称为原点平移法。
若A的特征值满足
1 2
希望p
n,
且2,n能估计时,我们就能确定P的近似值。
2 n
2
使得应用幂法计算加速。
a21 an1 n (a11 a22 为A的特征多项式.
a12 a 22 an 2
a1n a2 n
a nn
ann ) n 1 (次级 n 2的项)
A的特征方程
( ) det( I A) 0 (1.1) 的根称为A的特征值. ( A)表示A的所有特征值的集合. 设为A的特征值, 相应的齐次方程组 ( I A )x 0 (1.2) 的非零解x称为A的对应于的特征向量.
x1 [1 -1 1]T , x2 [1 0 -1]T , x3 [1 2 1]T
定理1 设是矩阵A R nn的特征值, x是对应的非零特征 向量,则 (1) c是cA的特征值(常数c 0); (2) p为A pI的特征值,即 (A pI)x ( p)x; (3) k 是A k的特征值,即 A k x k x; 1 1 1 1 (4) 设A非奇异,则 0且 为A 的特征值,即 A x x.
物理、力学和工程技术中很多问题在数学上都归结为求矩阵 的特征值问题。例如,振动问题(大型桥梁或建筑物的振动、 机械的振动、电磁震荡等),结构屈曲,物理学中的某些临界 值的确定。它们都归结为下述数学问题。
定义1 已知A (aij ) nn , 则称
a 11 ( ) det( I A)
用初等变换求矩阵的特征值问题
参考文献 :
[ 1 ] 王萼芳 , 石生明. 高等代数 (第三版) [ M ] . 北京 : 高等教育出版社 , 2003. 1872193 , 2902299 [ 2 ] 张禾瑞 , 郝炳新. 高等代数 (第四版) [ M ] . 北京 : 高等教育出版社 , 1999. 2782287 [ 3 ] 李星 , 李宏伟. 高等代数学习指导与习题解析 [ M ] . 武昌 : 华中科技大学出版社 , 2005. 3002310 [ 4 ] 罗芳 , 康淑瑰. 求矩阵特征值的一种分治法 [J ] . 雁北师范学院学报 , 2003 , 19 (2) : 9211 [ 5 ] 陈泽安. 求矩阵特征值与特征向量的新方法 [J ] . 长沙通信职业技术学院学报 , 2003 , 2 (1) : 66269
角线上的元素就为所求的矩阵 A 的特征值.
·11 ·
2006 年 12 月 河北北方学院学报 (自然科学版) 第 6 期
3
列的
3 2
倍加到
2
列
1 2
1 2
23 2
4 列的 1 倍加到 2 列
0 1 2 21
0000
0000
所以矩阵
k
引理 1[1] : 对一个 s ×n 的矩阵 A 作一初等行变换就相当于 A 左乘一个相应的 s ×s 初等矩阵 ; 对 A 作
一初等列变换就相当于 A 右乘一个相应的 n ×n 初等矩阵. 且任意一个可逆矩阵都可以表示为一系列初等
矩阵的乘积.
定义 3[1] : 设 A , B 是两个 n 阶矩阵 , 如果存在 P 上的 n 阶可逆矩阵 X , 使得 B = X21 A X , 则称 A 相
A
的全部特征值为
1 2
第五章矩阵特征值问题同步复习
第五章矩阵特征值问题同步复习第五章矩阵特征值问题一、内容提要§5.1 特征值与特征向量1.定义设A 为阶方阵,如果存在数n λ以及一个非零n 维列向量ξ,使得关系式λξξ=A 成立,则称λ为A 的一个特征值,非零向量ξ为A 的属于特征值λ的特征向量。
2.求特征值和特征向量的步骤:(1)计算特征多项式A I -λ;(2)求A 的特征方程A I -λ=0的全部根,它们就是A 的所有特征值;(3)对于A 的每一个特征值λ,求解齐次线性方程组()0=。
设它的一个-X A I λ基础解系为,,,,21r n -ξξξ (其中)(A I r r -=λ),则A 的属于λ的全部特征向量为,2211r n r n k k k --+++ξξξ其中是不全为零的任意数。
r n -21k k k ,,,3.性质● 方阵A 与其转置矩阵T A 有相同的特征多项式,从而有相同的特征值;● )(21A tr n =+++λλλ , A n =λλλ 21;● 可逆矩阵A 与1-A 的特征值互为倒数;● 设λ是矩阵A 的特征值,)(x g 是一个多项式,则)(λg 是)(A g 的特征值;● 如果n 阶矩阵A 有n 个不同的特征值,则A 有n 个线性无关的特征向量;● 设s λλλ,,,21 是矩阵A 的s 个互不相同的特征值,而i in i i ααα,,,21 是A 的分别对应于特征值i λ的线性无关的特征向量组,则向量组111211,,,n ααα ; 222221,,,n ααα ; ...; ssn s s ααα,,,21 线性无关.§5.2 矩阵的相似性1.定义设A ,都是阶方阵,如果阶可逆矩阵B n P ,使B AP P =-1,则称矩阵A 与相似,记为B B A ~。
如果P 为正交矩阵,则称A 与B 正交相似。
2.命题相似矩阵有相同的特征多项式,从而有相同的特征值,相同的行列式和迹。
3.对角化的条件(1)充要条件:n 阶方阵A 可对角化的充分必要条件是A 有n 个线性无关的特征向量。
jacobi方法求特征值和特征向量 例题
一、引言Jacobi方法是一种用于计算矩阵特征值和特征向量的迭代数值方法。
它是数值线性代数中的重要算法之一,广泛应用于科学计算、工程技术和金融领域。
本文将通过一个例题来介绍Jacobi方法的原理和求解过程,并分析其在实际问题中的应用。
二、Jacobi方法的原理Jacobi方法是一种通过迭代对矩阵进行相似变换,使得原矩阵逐步转化为对角矩阵的方法。
通过数值迭代,可以逐步逼近矩阵的特征值和对应的特征向量。
其基本原理如下:1. 对称矩阵特征值问题:对于对称矩阵A,存在一个正交矩阵P,使得P^T * A * P = D,其中D为对角矩阵,其对角线上的元素为A的特征值。
所以我们可以通过迭代找到P,使得P逼近正交矩阵,从而逼近A的特征值和特征向量。
2. Jacobi迭代:Jacobi方法的基本思想是通过正交相似变换,逐步将矩阵对角化。
具体来说,对于矩阵A,找到一个旋转矩阵G,使得A' = G^T * A * G为对角矩阵,然后递归地对A'进行相似变换,直到达到精度要求。
三、Jacobi方法求解特征值和特征向量的例题考虑以下矩阵A:A = [[4, -2, 2],[-2, 5, -1],[2, -1, 3]]我们将通过Jacobi方法来计算矩阵A的特征值和特征向量。
1. 对称化矩阵我们需要对矩阵A进行对称化处理。
对称化的思路是找到正交矩阵P,使得P^T * A * P = D,其中D为对角矩阵。
我们可以通过迭代找到逼近P的矩阵序列,直到达到一定的精度。
2. Jacobi迭代在Jacobi迭代的过程中,我们需要找到一个旋转矩阵G,使得A' =G^T * A * G为对角矩阵。
具体的迭代过程是:找到矩阵A中绝对值最大的非对角元素a[i][j],然后构造一个旋转矩阵G,将a[i][j]置零。
通过迭代地对A'进行相似变换,最终使得A'的非对角元素逼近零,即达到对角化的目的。
3. 计算特征值和特征向量经过一定次数的Jacobi迭代后,得到了对称矩阵A的对角化矩阵D和正交矩阵P。
矩阵特征值求解
矩阵特征值求解的分值算法12组1.1矩阵计算的基本问题(1)求解线性方程组的问题•即给定一个n阶非奇异矩阵A和n维向量b,求一个n维向量X,使得Ax =b (1. 1. 1 )(2)线性最小二乘问题,即给定一个mx n阶矩阵A和m维向量b ,求一个n维向量X,使得|AX -b| =min{ | Ay -比严R n} (1.1.2 )(3)矩阵的特征问题,即给定一个n阶实(复)矩阵A,求它的部分或全部特征值以及对应的特征向量,也就是求解方程(1. 1. 3 )一对解(4 X),其中R(C), x- R n(C n),即A为矩阵A的特征值,X为矩阵Ax = ZxA的属于特征值A的特征向量。
在工程上,矩阵的特征值具有广泛的应用,如大型桥梁或建筑物的振动问题:机械和机件的振动问题;飞机机翼的颤振问题;无线电电子学及光学系统的电磁振动问题;调节系统的自振问题以及声学和超声学系统的振动问题•又如天文、地震、信息系统、经济学中的一些问题都与矩阵的特征值问题密切相关。
在科学上,计算流体力学、统计计算、量子力学、化学工程和网络排队的马尔可夫链模拟等实际问题,最后也都要归结为矩阵的特征值问题.由于特征值问题在许多科学和工程领域中具有广泛的应用,因此对矩阵的特征值问题的求解理论研究算法的开发软件的制作等是当今计算数学和科学与工程计算研究领域的重大课题,国际上这方面的研究工作十分活跃。
1.2矩阵的特征值问题研究现状及算法概述对一个nxn阶实(复)矩阵A,它的特征值问题,即求方程(1.1.3)式的非平凡解,是数值线性代数的一个中心问题•这一问题的内在非线性给计算特征值带来许多计算问题•为了求(1.1.3)式中的A ,—个简单的想法就是显式地求解特征方程det (A 一几I)二0 (121 ) 除非对于个别的特殊矩阵,由于特征方程的系数不能够用稳定的数值方法由行列式的计算来求得,既使能精确计算出特征方程的系数,在有限精度下,其特征多项式f〃)二det(A-ZJ)的根可能对多项式的系数非常敏感能•因此,这个方法只在理论上是有意义的,实际计算中对一般矩阵是不可行的数 _ . _ . 人较大,则行列式det (A -几I)的计算量将非常大;其次,根据•首先,右矩即AfbJ阳数大于四的多项式求根不存在一种通用的方法,基于上述原Galois理论对于次因,人们只能寻求其它途径•因此,如何有效地!精确地求解’矩阵特征值问题,就成为数值线性代数领域的一个中心问题.目前,求解矩阵特征值问题的方法有两大类:一类称为变换方法,另一类称为向量迭代方法•变换方法是直接对原矩阵进行处理,通过一系列相似变换,使之变换成 一个易于求解特征值的形式,如Jacobi 算法,Givens 算法,QR 算法等。
第五章 矩阵的特征值与特征向量 习题
第五章 矩阵的特征值与特征向量 习题1. 试用施密特法把下列向量组正交化:(1)⎪⎪⎪⎭⎫ ⎝⎛=931421111) , ,(321a a a ;(2)⎪⎪⎪⎪⎪⎭⎫ ⎝⎛---=011101110111) , ,(321a a a . 2. 设x 为n 维列向量, x T x =1, 令H =E -2xx T , 证明H 是对称的正交阵. 3. 求下列矩阵的特征值和特征向量:(1)⎪⎪⎪⎭⎫ ⎝⎛----201335212; (2)⎪⎪⎪⎭⎫ ⎝⎛633312321.4. 设A 为n 阶矩阵, 证明A T 与A 的特征值相同.5. 设λ≠0是m 阶矩阵A m ⨯n B n ⨯m 的特征值, 证明λ也是n 阶矩阵BA 的特征值.6. 已知3阶矩阵A 的特征值为1, 2, 3, 求|A 3-5A 2+7A |.7. 已知3阶矩阵A 的特征值为1, 2, -3, 求|A *+3A +2E |.8. 设矩阵⎪⎪⎪⎭⎫ ⎝⎛=50413102x A 可相似对角化, 求x .9. 已知p =(1, 1, -1)T 是矩阵⎪⎪⎪⎭⎫ ⎝⎛---=2135212b a A 的一个特征向量.(1)求参数a , b 及特征向量p 所对应的特征值;(2)问A 能不能相似对角化?并说明理由.10. 试求一个正交的相似变换矩阵, 将对称阵⎪⎪⎪⎭⎫ ⎝⎛----020212022化为对角阵.11. 设矩阵⎪⎪⎪⎭⎫ ⎝⎛------=12422421x A 与⎪⎪⎪⎭⎫ ⎝⎛-=Λy 45相似, 求x , y ; 并求一个正交阵P , 使P -1AP =Λ.12. 设3阶方阵A 的特征值为λ1=2, λ2=-2, λ3=1; 对应的特征向量依次为p 1=(0, 1, 1)T , p 2=(1, 1, 1)T , p 3=(1, 1, 0)T , 求A .13. 设3阶对称矩阵A 的特征值λ1=6, λ2=3, λ3=3, 与特征值λ1=6对应的特征向量为p 1=(1, 1, 1)T , 求A .14. 设⎪⎪⎪⎭⎫ ⎝⎛-=340430241A , 求A 100.。
第五章 特征值和特征向量、矩阵的对角化 扩展例题及求解
[例 1]设 A , B 是 n 阶方阵,且 R(A) R(B) n ,证明 A , B 有公共的特征值,有公共的
特征向量。 [分析]此类题型综合程度较高,主要根据已知条件推导出对证明有用的结论,需要对知识点 熟练掌握。
[证] R(A) R(B) n R(A) n , R(B) n AX 0 , BX 0 都有非零解。
i 1
i 1
a1
a2
n
a3
an
n
从而特征值 ai2 的特征向量为 knn ,( kn 0 )。 i 1
[例]6 证明:设 A 是 n 阶正交矩阵, 是 A 的特征值,则 1 。
[分析]本题要利用正交矩阵的性质来证明。要证 1 ,可证 1。 [证]设 是对应于 的特征向量,则
[例 8]设 A 为 3 阶方阵, 为 3 维向量, , A , A2 线性无关, A3 4A 3A2 ,
计算行列式 2A2 3E 。
[分析]这里根据已知条件,直接计算方阵的行列式比较困难,可以采取特征值和特征向量的
性质来计算,即设 n 阶方阵 A 的特征值为 1, 2, , n ,则 A 12 n 。
不妨设 , 分别为 AX 0 , BX 0 的非零解,则
A 0 , B 0
由特征值、特征向量的定义, 0 是 A , B 的特征值, , 分别为 A , B 对应于特征值 0
的特征向量。
要证有公共的特征向量,只要证明 AX 0 ,BX 0 有公共的非零解即可。事实上,AX 0 ,
乘 ( A 4E)1 得: A(A E) 0 ,即 A2 A 0 ,从而 , A , A2 线性相关,矛
方法求矩阵全部特征值
方法求矩阵全部特征值要求矩阵全部特征值的方法有多种,其中最常用的方法是使用特征值分解或者通过求解矩阵的特征多项式来得到。
特征值分解是一种将矩阵表示为特征向量和特征值的形式的方法。
对于一个nxn的矩阵A,特征向量x满足Ax=λx,其中λ为特征值。
特征向量x可以通过求解方程(A-λI)x=0来获得,其中I为单位矩阵。
步骤如下:1. 对于给定矩阵A,求解特征多项式det(A-λI)=0,可以得到一个关于λ的n次方程。
2.解这个n次方程,求得n个特征值λ1,λ2,...,λn。
这些特征值可能是重复的。
3. 对于每个特征值λi,解方程(A-λiI)x=0,可以得到对应的特征向量xi。
4.可以验证Ax=λx是否成立来验证特征值和特征向量的正确性。
特征值分解的优点是准确性高,能得到精确的特征值和特征向量。
但是它的计算量较大,对于大型矩阵来说可能需要较长的计算时间。
除了特征值分解,还有一些其他的方法可以用来求解矩阵的全部特征值。
以下是一些常用的方法:1. 幂迭代法(Power Iteration Method):该方法通过反复迭代矩阵与一个初始向量的乘积来不断逼近最大的特征值。
通过迭代,该方法可以找到矩阵的一个特征值及其对应的特征向量。
2. 反幂迭代法(Inverse Power Iteration Method):该方法是幂迭代法的变种,用来求解最小特征值及其对应的特征向量。
3. QR迭代法(QR Iteration Method):该方法通过迭代进行QR分解,逐渐将矩阵转化为上三角矩阵,在迭代的过程中得出矩阵的全部特征值。
4. 特征值转换法(Eigenvalue Transformation Method):该方法通过变换矩阵的形式,转化为带有一些特殊特征的矩阵,例如Hessenberg矩阵或Schur矩阵,从而更容易求解特征值。
这些方法各有优缺点,应根据具体情况选择合适的方法。
另外,对于一些特殊类型的矩阵,例如对称矩阵或正定矩阵,还可以使用更为高效的方法来求解特征值。
第五章 矩阵的特征值和特征向量
第五章 矩阵的特征值和特征向量习题一 矩阵的特征值和特征向量一、填空题1.A 为n 阶方阵,Ax =0有非零解,则A 必有一特征值为________. 2.若λ0为A 的特征值,则A k (k 为正整数)有特征值为________.3.若α为A 的特征向量,则________为P -1AP 的特征向量. 4.n 阶矩阵A 与_____________有相同的特征值. 二、计算题1.设A =⎪⎪⎪⎭⎫ ⎝⎛-----122212221 (1) 试求矩阵A 的特征值;(2) 利用(1)的结果,求矩阵E +A -1的特征值,其中E 是三阶单位矩阵.2.求矩阵A =⎪⎪⎪⎭⎫ ⎝⎛----632223221的实特征值及对应的特征向量.三、证明题1.设A 满足A 2-3A +2E =0,证明其特征值只能取值1或2.2.若n 阶矩阵A ,存在自然数m ,使得0=mA ,则A 的特征值是0.3.如果A 可逆,λ是A 的特征值,则1-λ是1-A 的特征值.4.证明:)()(),()()(A kTr kA Tr B Tr A Tr B A Tr =+=+.习题二 相似矩阵和矩阵可对角化一、填空题1.若A ~kE ,则A =________.2.若n 阶方阵A 与B 相似,且A 2=A ,则B 2=________.. 3.已知A =⎪⎪⎪⎭⎫⎝⎛----533242111,B =⎪⎪⎪⎭⎫ ⎝⎛20002000λ 且A ~B ,则λ=________.4.A 可对角化当且仅当 . 5.n 阶矩阵A 有n 个互不相同的特征值是A 可对角化的___________. 6.判别矩阵A 可对角化的方法是 .二、 1.设A =[a ij ]为三角矩阵,且对角线元素互不相等.试指出A 是否有与它相似的对角矩阵,并说明理由.2.矩阵A =⎪⎪⎪⎭⎫ ⎝⎛--314020112能否对角化?若能,求可逆矩阵P ,使P -1AP 为对角矩阵.三、判别下列矩阵是否可对角化⎪⎪⎪⎭⎫ ⎝⎛=001010100A ⎪⎪⎪⎭⎫ ⎝⎛---=031302120B四、矩阵A=⎪⎪⎪⎭⎫⎝⎛-113222x和B=⎪⎪⎪⎭⎫⎝⎛-y21是相似矩阵.求x与y;习题三实对称矩阵的对角化一、求正交矩阵T,使ATT1-为对角矩阵.①⎪⎪⎪⎭⎫⎝⎛----=34243222A②⎪⎪⎪⎭⎫⎝⎛----=1222223B③⎪⎪⎪⎪⎪⎭⎫⎝⎛----=411141141114C二、设实对称矩阵A =⎪⎪⎪⎭⎫ ⎝⎛-124222421,求可逆矩阵Q ,使Q -1AQ 为对角矩阵.三、已知三阶方阵A 的特征值为1,-1,2,设矩阵B =A 3-5A 2试求:(1) 矩阵B 的特征值及与其相似的对角阵;(2) 行列式|B |和|A -5E |.四、设A =⎪⎪⎪⎭⎫ ⎝⎛201021313,求 (1) A 的所有特征值与特征向量;(2) 判断A 能否对角化,若能对角化,则求出可逆矩阵P ,使A 化为对角形矩阵; (3) 计算A m .综合复习题一、空题与选择题1.矩阵________20222002⎪⎪⎭⎫ ⎝⎛⎪⎪⎭⎫ ⎝⎛与.2.设),(,F M T A n ∈其中T 可逆,则k AT T k _(__________)(1=-为非负整数), ][)(_,__________)(1x F x f AT T f ∈=- .(][x F 表示数域F 的全体多项式,)(F M n 表示全体n 阶矩阵) 3.相似矩阵有________秩,有相同秩的矩阵_________相似.4.设⎪⎪⎪⎭⎫ ⎝⎛=411205123A 的三个特征值为321,,λλλ则 .________.____________321321==++λλλλλλ 5.设)(x f 是方阵A 的特征多项式,则_______;)(=A f若B A ~,则)(B f = _________.6.下面四个命题中原命题和逆命题都正确的是( ) (A )相似矩阵有相同的特征多项式;(B )设σ是数域F 上向量空间的一个线性变换.A 是σ关于V 的一个基的矩阵,如果λ是σ的特征根,那么λ是A 的特征根;(C )n 维向量空间的一个线性变换关于V 的两个基的矩阵是相似矩阵; (D )设)(F M A n ∈,若)(x f A 在数域F 内有单根,则A 可对角化. 7.下列三个矩阵中⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛=a a a A a a a A a a a A 001001,001000,000000321 ① 21~A A ; ② 31~A A ; ③ 32~A A ;④ 321,,A A A 中两两都不相似.(A )① 正确; (B) ②正确; (C) ③ 正确 ; (D) ④ 正确. 8.设A 是n 阶矩阵,那么① 在复数域C 上A 一定与某一对角矩阵相似; ② 在C 上A 一定与某一上三角矩阵相似;③ 在C 上A 一定与某一下三角矩阵相似.(A )① 正确; (B) ②,③正确; (C) ①, ② 正确 ; (D) ①,②,③正确. 9.下列矩阵中,不可对角化的仅是(A)⎪⎪⎭⎫ ⎝⎛--0280; (B) ⎪⎪⎭⎫ ⎝⎛1111; (C) ⎪⎪⎭⎫ ⎝⎛---1101; (D)⎪⎪⎭⎫ ⎝⎛-3210. 10.设,0),(,,≠∈T F M T B A n 且B A ,在F 上均可对角化,则① B A + 可对角化; ②AB 可对角化; ③AT T 1-可对角化;④ T B T m 1-可对角化. *N m ∈(*N 表示全体正整数) (A ),②正确; ( B) ③,④正确; (C) ①,②,③,④正确 ; (D) ① 正确. 二、计算与证明题1.求下列矩阵的全部特征值与特征向量(1)⎪⎪⎪⎭⎫ ⎝⎛-=200210311A (2)⎪⎪⎪⎭⎫⎝⎛=624232426A (3)⎪⎪⎪⎭⎫ ⎝⎛=633312321A (4)⎪⎪⎪⎭⎫⎝⎛--=201034011A2.找出1题中可对角化的矩阵A ,并求可逆矩阵X 使AX X1-为对角矩阵.3.求正交矩阵T 使AT T1-为对角矩阵(1) ⎪⎪⎪⎭⎫ ⎝⎛----=542452222A (2)⎪⎪⎪⎭⎫⎝⎛----=342432220A (3)⎪⎪⎪⎪⎪⎭⎫ ⎝⎛----=1333313333133331A (4)⎪⎪⎪⎪⎪⎭⎫ ⎝⎛----=1101111001111011A 4.试证:矩阵A 可逆的充分必要条件是:它的特征值都不等于零.5.设n 阶可逆矩阵A 的特征值是n λλλ,,,21 ,证明:1-A 的特征值为11211,,,---n λλλ . 6.如果任一个n 维非零向量都是n 阶矩阵A 的特征向量,试证明A 是一个数量矩阵. 7.A 是一个n 阶实对称矩阵,试证:如果0λ是A 的k 重特征值,则矩阵A E -0λ的秩等于k n -.自测题一、填空题1.若A 为n 阶矩阵,0=AX 有非零解,则A 必有一特征值为__________. 2.若0λ是A 特征值,则kA (k 为正整数)有特征值为____________.3.若α为A 的特征向量,则AP P 1-的特征向量为_____________.4.若n 阶矩阵A 有n 个属于特征值λ的线性无关的特征向量,则A =__________.5.已知三阶矩阵A 的三个特征值为1,2,3,则1_____;-=A A 的特征值为___________.6.n 阶零矩阵的全部特征向量是___________. 7.若kE A ~,则=A ______________.8.若n 阶矩阵A 与B 相似,且A A =2,则=2B ___________.9.已知⎪⎪⎪⎭⎫⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛----=20002000,533242111λB A 且B A ~,则._______=λ 10.三阶矩阵A 的三个互异特征值为321,,λλλ,它们对应的特征列向量分别为 ,,,321ααα则矩阵(,,,321ααα)的秩为__________.二、选择题1.设λ=2是非奇异矩阵A 的特征值,则矩阵12)31(-A 有一特征值等于( ).(a ) 34 (b ) 43 (c ) 21 (d ) 412.若n 阶矩阵A 的任意一行中n 个元素的和都是a ,则A 的一个特征值为( ).(a ) a (b ) –a (c ) 0 (d ) a -13.设A 是n 阶矩阵,λ1,λ2是A 的特征值,α1,α2是A 的分别对应于λ1,λ2的特征向量,则( ).(a ) λ1=λ2时,α1,α2一定成比例 (b ) λ1=λ2时,α1,α2一定不成比例(c ) λ1≠λ2时,α1,α2一定成比例 (d ) λ1≠λ2时,α1,α2一定不成比例4.设n 阶矩阵A 与B 相似,则( )(a ) λE -A =λE -B (b ) |λE -A |=|λE -B |(c ) |λE -A |~λE -B (d ) A 与B 都相似于一个对角矩阵D5.n 阶方阵A 具有n 个特征值是A 与对角矩阵相似的( )(a ) 充分必要条件 (b ) 充分而非必要条件 (c ) 必要而非充分条件 (d ) 既非充分也非必要条件6.矩阵A =⎪⎪⎪⎭⎫ ⎝⎛300030000与下列哪个矩阵相似( )(a ) ⎪⎪⎪⎭⎫ ⎝⎛000030300 (b ) ⎪⎪⎪⎭⎫ ⎝⎛300130010 (c ) ⎪⎪⎪⎭⎫ ⎝⎛300000003 (d ) ⎪⎪⎪⎭⎫ ⎝⎛030300010 7.n 阶矩阵与对角矩阵相似的充分必要条件是( ).(a ) A 有n 个不全相同的特征值 (b ) A T 有n 个不全相同的特征值 (c ) A 有n 个不相同的特征值 (d ) A 有n 个线性无关的特征向量8.n 阶方阵A 与某对角矩阵相似,则( ).(a ) 方阵A 的秩等于n (b ) 方阵A 有n 个不同的特征值(c ) 方阵A 一定是对称矩阵 (d ) 方阵A 有n 个线性无关的特征向量9.λ1,λ2是n 阶矩阵A 的特征值,X 1,X 2是相应于λ1,λ2的特征向量,对于不全为零的常数c 1,c 2:( )(a ) 当λ1≠λ2时,则c 1X 1+ c 2X 2必为A 特征向量(b ) 当λ1≠λ2时,则X 1,X 2是A 相应于λ1,λ2唯一的两个线性无关的特征向量 (c ) 当λ1=λ2时,则c 1X 1+ c 2X 2必为A 特征向量(d ) 当λ1=λ2时,则X 1,X 2必为A 相应于λ1,λ2的线性无关的特征向量 10.设n 阶矩阵A 为满秩矩阵,则A ( )(a ) 必有n 个线性无关的特征值 (b ) 必有n 个线性无关的特征向量 (c ) 必相似于一满秩的对角矩阵 (d ) 特征值必不为零 三、计算题1.设⎪⎪⎪⎭⎫ ⎝⎛-----=122212221A (1) 试求矩阵A 的特征值;(2)利用(1)的结果,求1-+A E 的特征值.2.求矩阵⎪⎪⎪⎭⎫ ⎝⎛----=632223221A 的特征值及特征向量.3.设实对称矩阵⎪⎪⎪⎭⎫ ⎝⎛-=124222421A ,求可逆矩阵Q 使AQ Q 1-为对角矩阵. 4.设A 为n 阶实矩阵,满足0,<=A E AA T,试求A 的伴随矩阵*A 的一个特征值.5.已知三阶矩阵A 的特征值为1,1-,2,矩阵235A A B -=,试求 (1) 矩阵B 的特征值和与B 相似的对角矩阵;(2) 行列式B 和E A 5-.6.设⎪⎪⎪⎭⎫ ⎝⎛=201021313A ,求 (1)A 的所有特征值与特征向量;(2)判别A 能否对角化,若能对角化,则求出可逆矩阵P ,使AP P 1-为对角矩阵; (3)计算mA .四、证明题1.若n 阶矩阵A 满足A A =2,则A 的特征值仅能是0或1.2.若n 阶矩阵A 满足I A =2,则A 的特征值仅能是1或1-.3.设A 满足0232=+-E A A ,证明:A 的特征值只能是1或2.4.设A 是实数域上奇数阶方阵,且0>A ,证明:A 有正特征值. 5.设][)(),(x F x f F M A n ∈∈,A 在F 上可对角化,证明:)(A f 在F 上可对角化.二次型习题一 二次型及表示方法一、填空题1.二次型f (x 1,x 2,x 3,x 4)=x 12+2x 22+3x 32+4x 1x 2+2x 2x 3________.2. 矩阵A =⎪⎪⎪⎭⎫ ⎝⎛--314122421对应的二次型是________.3.),(21x x q =⎪⎪⎭⎫ ⎝⎛⎪⎪⎭⎫ ⎝⎛-21211222)(x x x x 的矩阵为__________. 4.二次型),,,(21n x x x q 经过__________的线性替换总可以化为标准形2222211nn y c y c y c +++ .5.n 阶对称矩阵同时实行行和列的初等变换总可化为_______矩阵. 二、写出下列各二次型的矩阵1.23322231212138232x x x x x x x x x ++-+-2.243231212x x x x x x x ++-三、写出下列对称矩阵所对应的二次型1.⎪⎪⎪⎪⎪⎭⎫ ⎝⎛--------=012320113113221233121A2.⎪⎪⎪⎭⎫ ⎝⎛=011121110B四、对于对称矩阵A 与B ,求出可逆矩阵C ,使B AC C T=⎪⎪⎪⎭⎫ ⎝⎛=011121110A ⎪⎪⎪⎭⎫ ⎝⎛=011101112B习题二 化二次型为标准型一、用配方法化下列二次型为标准型. 1.31212322214245x x x x x x x -+-+2. 32312164x x x x x x +-二、用初等变化的方法求一奇异矩阵C ,使AC C T为对角矩阵.⎪⎪⎪⎭⎫ ⎝⎛=310102021A三、用初等变换法将二次型f (x 1,x 2,x 3,x 4)=x 12+x 22+x 32+x 42+2x 1x 2+2x 2x 3+2x 3x 4化为规范形,并求所作的非退化变换矩阵,且用矩阵验算结果.四、求一正交矩阵P ,使AP P T为对角矩阵.⎪⎪⎪⎪⎪⎭⎫⎝⎛--------=1132112332112311A四、试用配方法将二次型f (x 1,x 2,x 3)=x 12+x 22+3x 32+4x 1x 2+2x 1x 3+2x 2x 3化为标准形(平方和)和规范形.习题三 正定二次型一、填空题1.实二次型f (x 1,x 2,x 3)=x 12-x 22+3x 32的秩为________,正惯性指数为________,负惯性指数为________.2.设n 阶实对称矩阵A 的特征值分别为1,2,…,n ,则当t ________时,tE -A 为正定矩阵.3. 若n 阶实对称矩阵A 的秩为r (<n )且A 2=A ,则是________矩阵(正定、半正定,…),正惯性指数为________.4.____二次型),,,(21n x x x q 成为正定的,如果对于任意一组),,,(21n c c c ______ 都有),,,(21n c c c q _________.5. 5. 对称矩阵A 正定当且仅当A 与_________矩阵合同.6.实对称矩阵A 正定当且仅当A 的一切顺序主子式__________或者A 的一切主子式 ______________.7. 7. 对称矩阵的特征根都是_____________.二、计算题:1.求α的值,使二次型为正定.(1) 3231212322214225x x x x x ax x x x +-+++(2) 3231212322212245x x x x x x ax x x --+++2.设矩阵A =⎪⎪⎪⎭⎫ ⎝⎛101020101,矩阵B =(kE +A )2,其中k 为实数,E 为单位矩阵.求对角矩阵Λ,使B 与Λ相似,并求k 为何值时,B 为正定矩阵.(1) (1) 3.设A 1~A 1和B 1~B 2.试证⎪⎪⎭⎫ ⎝⎛11B A ~⎪⎪⎭⎫ ⎝⎛22B A判断三元二次型f = x 12+5x 22+x 32+4x 1x 2-4x 2x 3的正定性.三、证明题:1.A 是n 阶实对称矩阵,AB +B T A 是正定矩阵,证明A 可逆.2.设A 是n 阶正定矩阵,证明|A +2E |>2n .3.令A =⎪⎪⎭⎫ ⎝⎛21A O O A , B =⎪⎪⎭⎫ ⎝⎛21B O O B ,如果1A 与1B 合同,2A 与2B 合同,则A 与B 合同.4.证明:实二次型),,,(21n x x x q 负定的充分必要条件是它的矩阵A 的奇数阶顺序主 子式全小于零,偶数阶顺序主子式全大于零.自测题一、填空题1.二次型322123222143212432,,,(x x x x x x x x x x x f ++++=)_______. 2.矩阵⎪⎪⎪⎭⎫ ⎝⎛--=314122421A 对应的二次型是_____________________. 3.二次型),,(321x x x f =31212322212224x x x tx x x x ++++是正定的,那么t 应满足不等式_________.4.二次型),,(321x x x f =2322213x x x +-的秩为__________.正惯性指数为__________,负惯性指数为__________.5.设n 阶实对称矩阵A 的特征值分别为n ,,2,1 ,则当t =______时,A tE -为正定矩阵.6.若n 阶实对称矩阵A 的秩为)(n r <且A A =2,则是_______矩阵,正惯性指数为___________.7.二次型的规范形由_____________唯一确定;复二次型的规范形由____唯一确定.8.实对称矩阵A 正定的充分必要条件是它的特征值___________. 9.若A 是实对称矩阵且可逆,则将Ax x f T =化为y A y f T 1-= 的线性变换为_____________.10.设A 为n 阶实对称矩阵,那么T AA 是_______(对称、非对称、对角).二、选择题i. i. 1.设A ,B 均为n 阶方阵,x =(x 1,x 2,…,x n )T ,且X T AX = X T BX ,当( )时,A =B .(a ) 秩(A )=秩(B ) (b ) A T =A(c ) B T =B (d ) A T =A 且B T =Bii. ii. 2.实二次型f (x 1,x 2,x 3,x 4)= X T AX 为正定的充分必要条件是( ).(a ) |A |>0 (b ) 存在n 阶可逆矩阵C ,使A =C T C(c ) 负惯性指数为零 (d ) 对于某一x =(x 1,x 2,…,x n )T ≠0,有X T AX >0.iii.iii. 3.实二次型f (x 1,x 2,x 3,x 4)= x 12+2x 1x 2+tx 22+3x 32,当t =( )时,其秩为2. (a ) 0 (b ) 1 (c ) 2 (d ) 3 iv. iv. 4.设A ,B 为同阶可逆矩阵,则( )(a ) AB =BA(b ) 存在可逆矩阵P ,使P -1AP =B(c ) 存在可逆矩阵C ,使C T AC =B(d ) 存在可逆矩阵P 和Q ,使P AQ =Bv.v. 5.设A 为正定矩阵,则下列矩阵不一定是正定的是( ) (a ) A T (b )A -1 (c ) A +E (d ) A -E vi.vi. 6.设A 是一个三阶实矩阵,如果对任一三维列向量X ,都有X T AX =0,那么( ). (a ) |A |=0 (b ) |A |>0 (c ) |A |<0 (d ) 以上都不是 vii. vii. 7.n 阶实对称矩阵A 为正定矩阵的充分必要条件是( ).(a ) 所有k 阶子式为正(k =1,2,…,n )(b ) A 的所有特征值非负(c ) A -1为正定矩阵 (d ) 秩(A )=nviii. viii. 8.设A ,B 都是n 阶实对称矩阵,且都正定,那么AB 是( )(a ) 实对称矩阵 (b ) 正定矩阵 (c ) 可逆矩阵 (d ) 正交矩阵ix. ix. 9.下列矩阵为正定的是( ).(a ) ⎪⎪⎪⎭⎫ ⎝⎛200032021 (b ) ⎪⎪⎪⎭⎫ ⎝⎛200042021 (c ) ⎪⎪⎪⎭⎫ ⎝⎛---200052021 (d ) ⎪⎪⎪⎭⎫ ⎝⎛520210002 x.x. 10.设A 、B 是n 阶正定矩阵,则( )是正定矩阵. (a ) A *+B * (b ) A *-B * (c ) A *B * (d ) k 1A *+k 2B *三、对二次型32212221442x x x x x x f --+=分别作下列两个非退化线性替换. (1) ⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎭⎫ ⎝⎛--=⎪⎪⎪⎭⎫ ⎝⎛321321*********y y y x x x (2) ⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎪⎪⎭⎫ ⎝⎛--=⎪⎪⎪⎭⎫ ⎝⎛32132121001101121y y y x x x四、试用配方法将二次型3231212322213212243),,(x x x x x x x x x x x x f +++++=化为标准形(平方和)和规范形.五、用初等变换法将二次型43322142322214321222),,,(x x x x x x x x x x x x x x f -+++++= 化为标准形,求所作的非退化矩阵,并用矩阵验算结果.六、已知二次型)0(2332),,(32232221321>+++==a x ax x x x x x x f ,通过正交变换化成标准形23222152y y y f ++=,求参数a 及所用正交变换矩阵. 七、设矩阵⎪⎪⎪⎭⎫ ⎝⎛=101020101A ,矩阵2)(A kE B +=,其中k 为实数,E 为单位矩阵,求 对角矩阵A ,使B 与A 相似,并求k 为何值时,B 为正定矩阵.八、设1A 与2A 相似,1B 与2B 相似.试证⎪⎪⎭⎫ ⎝⎛11B A 与⎪⎪⎭⎫ ⎝⎛22B A . 九、判断三元二次型3221232221445x x x x x x x f -+++=的正定性. 十、A 是n 阶实对称矩阵,A B AB T +是正定矩阵,证明:A 可逆.十一、设A 是n 阶正定矩阵,证明:n E A 22>+.。
特征值问题与特征方程
特征值问题与特征方程特征值问题与特征方程是线性代数中重要的概念和工具。
在本文中,将详细介绍特征值问题与特征方程的定义、性质和解法。
一、特征值问题的定义在线性代数中,特征值问题是研究线性变换或矩阵对向量空间中的向量进行操作时的一个基本问题。
特征值问题可以描述为:给定一个线性变换或矩阵A,寻找一个非零向量v以及一个数λ,使得下式成立:Av = λv其中,A是一个n阶方阵,v是n维非零向量,λ是常数。
二、特征值问题的解法为了解决特征值问题,我们需要求解特征方程。
特征方程是通过对特征值问题进行变形得到的一个方程。
假设λ是A的一个特征值,v是与λ相应的特征向量。
那么我们有:(A - λI)v = 0其中,I是单位矩阵。
根据线性代数的基本定理,当且仅当(A - λI)的行列式为零时,方程(A - λI)v = 0有非零解v。
因此,我们可以将上述方程转化为一个特征方程。
三、特征方程的定义与求解特征方程是由特征值问题转化得到的一个方程。
具体来说,对于一个n阶方阵A,其特征方程可以表示为:det(A - λI) = 0其中,det表示矩阵的行列式。
解特征方程可以得到A的所有特征值。
解特征方程的方法有很多,常见的有代数学、数值法和矩阵迭代法等。
代数学方法是通过对特定类型的矩阵应用代数定理和求根公式来解特征方程。
数值法是通过数值计算方法来近似求解特征方程。
矩阵迭代法是通过迭代计算来逼近特征值和特征向量。
四、特征值问题的性质特征值问题具有许多重要的性质和应用。
以下是一些常见的性质:1. 特征值可以是实数,也可以是复数。
2. 特征值对应的特征向量是线性无关的。
3. 特征值和特征向量是矩阵的固有性质,不随矩阵变换而变化。
4. 特征值问题在许多领域中都有广泛的应用,例如物理、工程、计算机科学等。
总结:特征值问题与特征方程是线性代数中的重要概念。
通过求解特征方程,我们可以找到矩阵的特征值和特征向量。
特征值问题具有许多重要的性质和应用,对于理解矩阵的特性和解决实际问题具有重要意义。
第六章_特征值问题与矩阵变换
⎛ − 1 1 0⎞ ⎟ ⎜ 例2 求矩阵 A = ⎜ − 4 3 0 ⎟的特征值和特征向量 . ⎜ 1 0 2⎟ ⎠ ⎝
解
A的特征多项式为 −1− λ 1 0
2
3−λ 0 = ( 2 − λ ) (1− λ ) , 1 0 2−λ 所以A的特征值为 λ 1 = 2, λ 2 = λ 3 = 1. −4 当 λ 1 = 2时, 解方程( A − 2 E ) x = 0.由
A − λE =
⎛ − 3 1 0⎞ ⎜ ⎟ A − 2E = ⎜ − 4 1 0⎟ ⎜ 1 0 0⎟ ⎠ ⎝ 得基础解系
⎛ 1 0 0⎞ ⎜ ⎟ ~ ⎜ 0 1 0 ⎟, ⎜ 0 0 0⎟ ⎠ ⎝ ⎛ 0⎞ ⎜ ⎟ p1 = ⎜ 0 ⎟ , ⎜ 1⎟ ⎝ ⎠
所以kp1(k ≠ 0)是对应于 1 = 2的全部特征向量 λ
若A与B相似 , B与C相似 , 则A与C相似 .
结论.n维线性性空间V上的一些线性变换σ在V的 不同基下的矩阵是相似矩阵。
二、相似矩阵与相似变换的性质
定理6.6:
;
求齐次线性方程组( A − λ E ) x = 0 的一个基础 x 解系
η 1 ,η 2 ,
,η t
可得 A 的属于特征值 λ 的全部特征向量 k 1η 1 + k 2η 2 + + k t η t 其中 k 1 , k 2 , , k t 为不全为零的常数 .
注、 n 次多项式的求根 问题一般并不容易, 在实际问题中常常应用 近似计算公式来求 特征值
6.2、矩阵的相似变换
(一)、相似变换与相似矩阵的性质
一、相似变换与相似矩阵概念
定义1 设A, B都是 n阶矩阵, 若有可逆矩阵 P , 使
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
§1、特征值的估计
由于工程计算中求矩阵尤其是高阶矩阵的 精确特征值通常比较困难,而许多情况下我们 只需要知道特征值在什么范围内变化或者落在 什么区域内,例如判断方阵的幂级数是否收敛 只要看方阵的特征值的模或谱半径是否小于1, 因此特征值的估计就显得尤其必要,这方面的 理论在特征值问题中相当经典。
由于
实际上是 的
一个
维子空间,因此我们希望将
搜索极值的空间放大到任意
维子空
间 。而增大后的集合的极大值不会比原集
合的小,极小值也不会比原集合大。
58
设有 则
,并假定
,即
59
并且当
时等号成立。因此
60
一般地,我们有
定理4 (Courant-Fischer)设
是
Hermite矩阵,其特征值为
,则
存在Hermite矩阵特征值的极值原理
48
一、 Rayleigh商
二次型
,如果存在
,那么
所以如果
,我们自然也希望
49
定义1 设
是Hermite矩阵,称
为矩阵 的Rayleigh商。 注意到
因此我们可以把对 在单位球面
的极性的讨论限定 上。
50
单位球面 是闭集,又因为
是 的连续
函数,因此根据多元函数的最值定理,
在 上存在最大值和最小值。由于特征值与
对于广义特征值问题
,可以通过
适当选择位移(shift)或极点(pole) ,再通过 求逆,将之转化为SEP:
这种方法的优点是特征向量不变,矩阵 奇 异时也可以使用,并且在求解邻近 的特征 值或绝对值很小的特征值时效率较高。缺点仍 然是 一般不是特殊矩阵。
35
例 3 (广义特征值问题的Cayley变换法)
3
一、从特征值问题的稳定性说起
工程计算中,求解特征值问题
的特征对
时,由于数据往往带有误差,
因此我们计算出的特征对
,实际上是
扰动后的特征值问题
的解。这里
4
我们希望知道矩阵元素的变化对特征对的影响。
由于我们一般只知道
或 的某个上界,
因此有必要研究如何利用这样的上界,尽可能
准确地估计 与 、 与 之间的差距,从
13
例 6 矩阵
经过对角相似变换
后,得
14
三个行Gerschgorin圆分别收缩为:
15
Gerschgorin定理与对角占优矩阵有密切关系。
定义7 对方阵
,如果
则称矩阵 为按行对角占优矩阵。如果
则称矩阵 为按行严格对角占优矩阵。
16
定理8 (Levy--Desplanques) 严格对角占优矩阵是可逆矩阵。
33
例 1 (广义特征值问题的直接变换法)
对于广义特征值问题
,可以两边
做 的逆变换,将之转化为SEP:
这种方法的优点是特征向量不变,但缺点是矩 阵 奇异时不能使用,并且当矩阵 是正定 Hermite阵时,矩阵 一般不再是对称矩阵, 因此不是保结构的算法,从而使计算复杂。
34
例 2 (广义特征值问题的位移求逆法)
特征值问题
1
特征值问题是线性代数的研究重点,在理论和应 用上都非常重要。 理论上 ,矩阵的特征值就是线性算子的谱。因此 可以从泛函分析里找到理论的支撑和生长点。 应用上,常微分方程(ODE)和偏微分方程(PDE)中 许多问题都可以转化为矩阵特征值问题。 矩阵特征值问题的算法也是高性能计算机的主要 计算任务之一。可大致分为求解稠密、中小型矩 阵全部特征值的变换类方法和求解稀疏、大型矩 阵部分特征值的投影类方法。
与
是等
价的。注意到等价矩阵的特征对之间的关系:
44
虽然矩阵束
也存在与Jordan标准型类似
的Weierstrass标准型及Kronecher标准型,但也
同Jordan标准型一样,在数值计算上存在困难。
从数值观点看,更吸引人的是Moler和Stewart
(1973)提出的广义Schur分解:
对矩阵
,存在酉矩阵
是原广义特征值问
41
证明:由于 是Hermite正定矩阵,所以有
再根据
是Hermite矩阵,所以有酉相似
令
,则有
42
最后根据
因此
这说明 广义特征值
,得
是
的特征值,因此也是
的特征值。
43
例 6 (基于广义Schur分解的QZ算法****)
对于广义特征值问题
,我们的目标
是寻找可逆矩阵
,使得
均为标准型。我们称
证明:令
,则矩阵
的元素
因此
,所以矩阵
可逆,即矩阵 可逆。
17
根据定理8,严格对角占优矩阵 没有零特 征值,而
这说明矩阵 的特征值 可能满足
由此,我们可以将Gerschgorin定理看成定理8 的“推论” 。
18
事实上,设矩阵 的特征值 不属于 的 Gerschgorin区域 ,则有
因此矩阵
严格对角占优,根据定理8
37
设
是广义特征值问题
的
特征对,
是标准特征值问题
的特征对。因为矩阵 是Hermite阵,所以有
完备的标准正交特征向量系,如果 都是单
位特征向量,那么
38
由于
这说明,广义特征向量 对矩阵 是正交的, 或称特征向量矩阵 是 正交(共轭)的 又由于 所以 这说明,广义特征向量 对矩阵 是加权正 交的,或称特征向量矩阵 是加权 正交 (共轭)的。
30
一、从 矩阵的视角看特征值问题
标准特征值问题 (SEP)
,即为
按此视角,广义特征值问题(GEP)
,
即为
类似地,我们有二次特征值问题(QEP)
31
这里系数矩阵为 矩阵,将矩阵元素即 多项式的次数推广到 次,可得多项式特征值 问题(PEP):
更进一步,从线性推广到非线性,我们有非线 性特征值问题(NEP):
的特征值分别为
。令
(1)对 的任意特征值 征值 ,使得
(2)存在
的排列
,存在 的特 ,使得
6
遗憾的是矩阵的特征向量一般不是矩阵元素的 连续函数,因此不一定是稳定的。
例 2 矩阵
的特征值为
,特征向量为
和
。而
的特征值为 ,
特征向量为
和
。矩阵 的
特征向量在
处不连续。
7
二、盖尔(Gerschgorin)定理
对于广义特征值问题
,如果已经
得到特征值的近似值 ,那么通过Cayley变换,
可以将之转化为SEP:
显然 这种方法实质上仍然是位移求逆法(SI)。
36
例 4 (广义特征值问题的Cholesky分解法)
对于广义特征值问题
,当 是正
定Hermite矩阵时,可以通过Cholesky分解
将之转化为SEP:
矩阵 仍然是Hermite矩阵,因此此算法是保 结构算法,并且说明特征值全是实数。
(1)矩阵 的特征值都位于其行盖尔区域内;
(2)若矩阵 有 个盖尔圆构成的并集 是
连通区域,并且与其余
个盖尔圆均不相
交,则 中恰好有 的 个特征值。
10
(1)的证明:
设 有特征对
,这里Leabharlann ,则令,则
,因此
从而
11
例 5 矩阵
的三个行Gerschgorin圆分别是:
12
因为相似变换不改变特征值,为了得到特征值
换先将 化为上Hessenberg矩阵 ,将 化
为上三角阵 ;
(2)QZ迭代,即对
使用一步位移QR迭
代,将 化为上三角阵或准上三角阵 ,同
时将 化成上三角阵 。这是本算法的核心
部分。
(3)计算
的特征对,并据此求出原
GEP的特征对。
47
§3、Rayleigh商和广义Rayleigh商
在物理、信息等学科及理论研究中,经常会遇 到Hermite矩阵的二次型函数的商的最大化或 最小化。这种商包括一个Hermite矩阵的 Rayleigh商和两个Hermite矩阵的广义Rayleigh 商。L.Rayleigh在1930年代研究振动系统的小 振荡时,为了找到合适的广义坐标,提出的一 种特殊形式的商,被后人称为Rayleigh商。
39
由于特征向量 相互正交,因此它们也可以看
成 维内积空间 的一组基,只是内积空间
的内积必须定义为 内积
,即
这样,对任意
,可得
因此得展开式
40
例 5 (广义特征值问题的同时合同对角化)
对于广义特征值问题
,如果
均为Hermite矩阵,并且 还是正定矩阵,那 么存在可逆矩阵 ,使得
这里 题的特征值。
是Hermite ,则
54
其余特征值是否有类似结论呢?
设
,并假定
,即
由于
,因此前面的证明可修改为
55
类似地,假定
,即
由于
,因此前面的证明可修改为
56
综上,一般地,我们有
定理3 设
是Hermite矩阵,其特征
值为
,相应的标准正交特征向
量为
。令
则
57
遗憾的是,定理2中的公式实用价值不大。因为 它们将 的计算限定在求Rayleigh商在 的局 部极值上,因此必须先求出 ,这在数值上是 比较困难的。
NEP的主要算法都是基于求解非线性方程组的 Newton法及其变体,目前尚不成熟。
32
二、广义特征值问题
结构动力分析、信号处理、通信等学科中经常 需要求解广义特征值: