矩阵理论第一二章 典型例题
矩阵理论习题与答案

矩阵理论习题与答案矩阵理论习题与答案矩阵理论是线性代数中的重要内容之一,它在数学、工程、计算机科学等领域都有广泛的应用。
为了帮助读者更好地理解和掌握矩阵理论,本文将介绍一些常见的矩阵理论习题,并提供详细的答案解析。
一、基础习题1. 已知矩阵A = [[2, 3], [4, 5]],求A的转置矩阵。
答案:矩阵的转置是将其行和列互换得到的新矩阵。
所以A的转置矩阵为A^T = [[2, 4], [3, 5]]。
2. 已知矩阵B = [[1, 2, 3], [4, 5, 6]],求B的逆矩阵。
答案:逆矩阵是指与原矩阵相乘得到单位矩阵的矩阵。
由于B是一个2×3的矩阵,不是方阵,所以不存在逆矩阵。
3. 已知矩阵C = [[1, 2], [3, 4]],求C的特征值和特征向量。
答案:特征值是矩阵C的特征多项式的根,特征向量是对应于每个特征值的线性方程组的解。
计算特征值和特征向量的步骤如下:首先,计算特征多项式:det(C - λI) = 0,其中I是单位矩阵,λ是特征值。
解特征多项式得到特征值λ1 = 5,λ2 = -1。
然后,将特征值代入线性方程组 (C - λI)x = 0,求解得到特征向量:对于λ1 = 5,解得特征向量v1 = [1, -2]。
对于λ2 = -1,解得特征向量v2 = [1, -1]。
所以C的特征值为λ1 = 5,λ2 = -1,对应的特征向量为v1 = [1, -2],v2 = [1, -1]。
二、进阶习题1. 已知矩阵D = [[1, 2], [3, 4]],求D的奇异值分解。
答案:奇异值分解是将矩阵分解为三个矩阵的乘积,其中一个是正交矩阵,一个是对角矩阵。
计算奇异值分解的步骤如下:首先,计算D的转置矩阵D^T。
然后,计算D和D^T的乘积DD^T,得到一个对称矩阵。
接下来,求解对称矩阵的特征值和特征向量。
将特征值构成对角矩阵Σ,特征向量构成正交矩阵U。
最后,计算D^T和U的乘积D^TU,得到正交矩阵V。
矩阵理论习题答案

习 题 一1. 设λ为的任一特征值,则因 λλ22- 为A =-A 22O 的特征值, 故022=-λλ. 即 λ=0或2.2. A ~B , C ~D 时, 分别存在可逆矩阵P 和Q , 使得 P 1-AP =B , Q 1-CQ =D .令T =⎪⎪⎭⎫⎝⎛Q O O P 则 T 是可逆矩阵,且T 1-⎪⎪⎭⎫⎝⎛C O O A T =⎪⎪⎭⎫ ⎝⎛⎪⎪⎭⎫ ⎝⎛⎪⎪⎭⎫ ⎝⎛--Q O O P C O O A Q O O P 11=⎪⎪⎭⎫ ⎝⎛D O O B 3. 设i x 是对应于特征值i λ的特征向量, 则 A i x =i λi x , 用1-A 左乘得i i i x A x 1-λ=.即i i i x x A 11--λ= 故 1-i λ是A 的特征值, i =1,2,, n .4. (1) 可以. A E -λ=)2)(1)(1(-+-λλλ,=P ⎪⎪⎪⎭⎫ ⎝⎛--104003214, ⎪⎪⎪⎭⎫ ⎝⎛-=-2111AP P .(2) 不可以.(3) ⎪⎪⎪⎭⎫ ⎝⎛=110101010P , ⎪⎪⎪⎭⎫⎝⎛=-1221AP P .5. (1) A 的特征值是0, 1, 2. 故A =-(b -a )2=0. 从而 b =a .又11111-λ----λ----λ=-λaa aa A I =)223(22+---a λλλ将λ=1, 2 代入上式求得 a=0.(2) P =⎪⎪⎪⎭⎫ ⎝⎛-101010101.6. A I -λ=)1()2(2+-λλ, A 有特征值 2, 2, -1.λ=2所对应的方程组 (2I -A )x =0 有解向量p 1=⎪⎪⎪⎭⎫ ⎝⎛041, p 2=⎪⎪⎪⎭⎫ ⎝⎛401λ=-1所对应的方程组 (I +A )x =0 有解向量p 3=⎪⎪⎪⎭⎫⎝⎛101令 P =(p ,1p ,2p 3)=⎪⎪⎪⎭⎫ ⎝⎛140004111, 则 P 1-=⎪⎪⎪⎭⎫ ⎝⎛---4416414030121. 于是有A 100=P ⎪⎪⎪⎭⎫ ⎝⎛122100100P 1-=⎪⎪⎪⎭⎫⎝⎛-⋅-⋅-⋅---12412244023012122431100100100100100100100. 7. (1)A I -λ=)1(2+λλ=D 3(λ), λI -A 有2阶子式172111----λ=λ-4λ-4不是D 3(λ)的因子, 所以D 2(λ)=D 1(λ)=1, A 的初等因子为λ-1, 2λ. A 的Jordan 标准形为J =⎪⎪⎪⎭⎫ ⎝⎛-000100001设A 的相似变换矩阵为P =(p 1,p 2,p 3), 则由AP =PJ 得 ⎪⎩⎪⎨⎧==-=23211pAp Ap p Ap 0 解出P =⎪⎪⎪⎭⎫ ⎝⎛-----241231111; (2) 因为),2()1()(23--=λλλD 1)()(12==λλD D ,故A ~J =⎪⎪⎪⎭⎫ ⎝⎛200010011设变换矩阵为 P =(321,,p p p ), 则⎪⎩⎪⎨⎧=+==33212112p Ap p p Ap p Ap ⇒P =⎪⎪⎪⎭⎫ ⎝⎛---502513803 (3) ),2()1()(23-+=-=λλλλA I D ,1)(2+=λλD 1)(1=λD .A 的不变因子是,11=d ,12+=λd )2)(1(3-+=λλdA ~J =⎪⎪⎪⎭⎫ ⎝⎛--211 因为A 可对角化,可分别求出特征值-1,2所对应的三个线性无关的特征向量:当λ=-1时,解方程组 ,0)(=+x A I 求得两个线性无关的特征向量,1011⎪⎪⎪⎭⎫ ⎝⎛-=p ⎪⎪⎪⎭⎫ ⎝⎛-=0122p当λ=2时,解方程组 ,0)2(=-x A I 得⎪⎪⎪⎭⎫ ⎝⎛-=1123p , P =⎪⎪⎪⎭⎫ ⎝⎛---101110221(4) 因⎪⎪⎪⎭⎫ ⎝⎛---+=-41131621λλλλA I ~⎪⎪⎪⎭⎫ ⎝⎛--2)1(11λλ, 故A ~J =⎪⎪⎪⎭⎫ ⎝⎛10111设变换矩阵为P =),,(321p p p , 则⎪⎩⎪⎨⎧+===3232211pp Ap p Ap p Ap 21,p p 是线性方程组 0=-x A I )(的解向量,此方程仴的一般解形为p =⎪⎪⎪⎭⎫ ⎝⎛+-t s t s 3 取⎪⎪⎪⎭⎫ ⎝⎛-=0111p , ⎪⎪⎪⎭⎫ ⎝⎛=1032p为求滿足方程 23)(p p A I -=-的解向量3p , 再取 ,2p p = 根据 ⎪⎪⎪⎭⎫ ⎝⎛------t s t s 3113113622~⎪⎪⎪⎭⎫⎝⎛----t s t s s 00033000311 由此可得 s =t , 从而向量 T 3213),,(x x x =p 的坐标应満足方程s x x x -=-+3213取 T 3)0,0,1(-=p , 最后得P =⎪⎪⎪⎭⎫ ⎝⎛--010001131 8. 设 f (λ)=4322458-++-λλλλ. A 的最小多项式为 12)(3+-=λλλA m ,作带余除法得 f (λ)=(149542235-+-+λλλλ))(λA m +1037242+-λλ, 于是f (A )=I A A 1037242+-=⎪⎪⎪⎭⎫ ⎝⎛----346106195026483.9. A 的最小多项式为 76)(2+-=λλλA m , 设 f(λ)=372919122234+-+-λλλλ,则f (λ)=)()52(2λλA m ++2+λ. 于是 [f (A )]1-=1)2(-+I A .由此求出[f (A )]1-=⎪⎪⎭⎫ ⎝⎛-3217231 10. (1) λI -A =⎪⎪⎪⎭⎫ ⎝⎛---+41131621λλλ标准形⎪⎪⎪⎭⎫ ⎝⎛--2)1(00010001λλ, A 的最小多项式为 2)1(-λ;2) )1)(1(+-λλ; (3) 2λ.11. 将方程组写成矩阵形式:⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎭⎫ ⎝⎛----=⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛321321188034011d d d d d d x x x t x t x t x , ⎪⎪⎪⎭⎫ ⎝⎛=321x x x x , ⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=t x t x t x t d d d d d d d d 321x , A =⎪⎪⎪⎭⎫ ⎝⎛----188034011则有J =PAP 1-=⎪⎪⎪⎭⎫ ⎝⎛-100010011, .其中 P =⎪⎪⎪⎭⎫⎝⎛124012001.令 x =Py , 将原方程组改写成 : ,d d Jy y=t 则⎪⎪⎪⎩⎪⎪⎪⎨⎧-=+==3321211d d d d d d yty y y ty y t y 解此方程组得: y 1=C 1e t +C 2T e t , y 2=C 2e t , y 3=C 3e t -. 于是x =Py =⎪⎪⎪⎭⎫ ⎝⎛++++++-t t t tt t t c )t (c c )t (c c t c c e e 24e 4e 12e 2e e 3212121.12. (1) A 是实对称矩阵. A I -λ=2)1)(10(--λλ,A 有特征值 10, 2, 2.当λ=10时. 对应的齐次线性方程组 (10I -A )x =0的系数矩阵⎪⎪⎪⎭⎫ ⎝⎛--542452228~⎪⎪⎪⎭⎫ ⎝⎛000110102由此求出特征向量p 1=(-1, -2, 2)T , 单位化后得 e 1= (32,32,31--)T . 当λ=1时, 对应的齐次线性方程组 (I -A )x =0的系数矩阵⎪⎪⎪⎭⎫ ⎝⎛-----442442221~⎪⎪⎪⎭⎫ ⎝⎛-000000221 由此求出特征向量 p 2=(-2, 1, 0)T , p 3=(2, 0, 1)T . 单位化后得e 2=(0,51,52-)T , e 3=(535,534,532)T. 令 U =⎪⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛---53503253451325325231, 则 U 1-AU =⎪⎪⎪⎭⎫⎝⎛1110.(2) A 是Hermit 矩阵. 同理可求出相似变换矩阵U =⎪⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛---2121212i 2i 2i 21210, U 1-AU =⎪⎪⎪⎭⎫⎝⎛-22. 13. 若A 是Hermit 正定矩阵,则由定理1.24可知存在n 阶酉矩阵U , 使得U H AU =⎪⎪⎪⎪⎪⎭⎫ ⎝⎛n λλλ21, i λ﹥0, I =1, 2, , n . 于是A =U ⎪⎪⎪⎪⎪⎭⎫ ⎝⎛n λλλ21U H = U ⎪⎪⎪⎪⎪⎭⎫ ⎝⎛n λλλ 21U H U ⎪⎪⎪⎪⎪⎭⎫⎝⎛n λλλ21U H 令B =U ⎪⎪⎪⎪⎪⎭⎫⎝⎛n λλλ21U H 则 A =B 2.反之,当 A =B 2且B 是Hermit 正定矩阵时,则因Hermit 正定矩阵的乘积仍为Hermit 正定矩阵,故A 是Hermit 正定的. 14. (1)⇒(2). 因A 是Hermit 矩阵,则存在酉矩阵U,使得U H AU =diag(n λλλ,,,21 )令x =Uy , 其中 y =e k . 则 x ≠0. 于是x H Ax =y H (U H AU )y =k λ≧0 (k =1, 2, , n ).(2)⇒(3).A =U diag(n λλλ,,,21 )U H =U diag(n λλλ,,,21 )diag(n λλλ,,,21 )U H令 P =diag(n λλλ,,,21 )U H , 则 A =P H P . (3)⇒(1). 任取x ≠0, 有x H Ax =x H P H Px =22Px ≧0.习 题 二1. 1x =01i 42i 1+++-++=7+2,2x =1i)4i(4)2(i)1i)(1(2+-+-+-+=23, ∞x =max {}1i 42i 1,,,-+=4.2. 当 x ≠0时, 有 x ﹥0; 当 x ﹦0时, 显然有 x =0. 对任意∈λC , 有x λ=x nk kk nk kk λξωλλξω==∑∑==1212.为证明三角不等式成立,先证明Minkowski 不等式: 设 1≦p ﹤∞, 则对任意实数 x k ,y k (k =1, 2, , n )有pnk pk k y x 11)(∑=+≦∑∑==+nk ppk nk ppk y x 1111)()(证 当 p =1时,此不等式显然成立. 下设 p ﹥1, 则有∑=+nk pkk y x 1≦∑∑=-=-+++nk p kk k nk p kk k y x y y x x 1111对上式右边的每一个加式分别使用H ölder 不等式, 并由 (p -1)q =p , 得∑=+nk pkky x1≦qnk q p kk pnk pk qnk q p kk pnk pk y x y y x x 11)1(1111)1(11)()()()(∑∑∑∑=-==-=+++=qnk p k k pnk pk pnk p k y x y x 111111)]()()[(∑∑∑===++再用 qnk p k k y x 11)(∑=+ 除上式两边,即得 Minkowski 不等式.现设任意 y =(n ηηη,,,21 )T ∈C n , 则有∑=+=+nk kk k y x 12ηξω=∑=+nk k k k 12)(ηξω≦∑=+nk k k k k 12)(ηωξω≦∑∑==+nk j k nk k k 1212()(ηωξω=y x +.3. (1) 函数的非负性与齐次性是显然的,我们只证三角不等式.利用最大函数的等价定义:max(A , B )=)(21b a b a -++max(),b a y x y x ++≦max(b b a a y x y x ++,)=)(21b b a a b a b a y x y x y y x x --+++++≦)(21b a b a b a b a y y x x y y x x -+-++++ =)(21)(21b a b a b a b a y y y y x x x x -+++-++ =max( b a x x ,)+max( b a y y ,)(2) 只证三角不等式.k 1a y x ++k 2b y x +≦k 1a x +k 1a y +k 2b x +k 2b y =( k 1a x +k 2b x )+( k 1a y +k 2b y ) .4. 218132i 453i 11m +=+++++++=A ;66132i 453i 1222222F =+++++++=A ; 15m =∞A ;=1A 列和范数(最大列模和)=27+;∞A =行和范数(最大行模和)=9 ;5. 非负性: A ≠O 时S 1-AS ≠O , 于是 m 1AS S A -=>0. A =O 时, 显然A =0;齐次性: 设λ∈C , 则 λλλ==-m1)(S A S A m1ASS -=λA ;三角不等式: m11m1)(BSS AS S S B A S B A ---+=+=+≦B A BSS AS S +=+--m 1m 1;相容性: m11m1)(BS ASS S SAB S AB ---==≦m1m1BSS AS S --=A B .6. 因为I n ≠O , 所以n I >0.从而利用矩阵范数的相容性得:n n n I I I =≦n I n I ,即n I ≧1.7. 设 A =(A ij )∈C n n ⨯, x =∈ξξξT 21),,,(n C n , 且 A =ij ji a ,max , 则∑∑=ikk ik Ax ξa 1≦∑∑ikk ik a ξ=∑∑kiik k a ][ξ≦n A ∑kk ξ=∞m A 1x ;∑∑=ikk ikAx 22ξa≦∑∑ikk ika2][ξ=∑∑ikka 22][ξ=n A 2x ≦n A =∞m A 2x .8. 非负性与齐次性是显然的, 我们先证三角不等式和相容性成立. A =(a ij ), B =(b ij )∈C n m ⨯, C =(c st )∈C l n ⨯且 A =ij ji a ,max , B =ij ji a ,max , C =st ts c ,max . 则MBA +=max{m ,n }ij ij ji b a +,max ≦max{m ,n })(m ax ,ij ij ji b a +≦max{m ,n }(A +B )=max{m ,n }A +max{m ,n }B =M M B A +;MAC=max{m ,l }∑kkt ik ti c a ,max ≦max{m ,n }}{max ,∑kkt ik ti c a ≦max{m ,n }}{max 22,∑∑⋅kkt kikti c a (Minkowski 不等式)=max{m ,n }n AC ≦max{m ,n }max{n ,l }AC =M M C A .下证与相应的向量范数的相容性.设 x =∈ξξξT 21),,,(n C n , d =kmax {k ξ}, 则有∑∑=ikk ik a Ax ξ1≦∑∑ikk ik a ξ=∑∑ki ikka)(ξ≦∑kk na ξ=n A ∑kk ξ≦max{m ,n }A ∑kk ξ=1M x A ;2Ax =∑∑ikkik a2ξ≦∑∑ik k ik a 2)(ξ≦∑∑∑ikkkika )(22ξ(H ölder 不等式)=∑∑∑⋅kk ikik a 22ξ≦mn A 2x≦max{m ,n }A 2x =2M x A ;}{max 1∑=∞=n k k ik iAxξa ≦∑=nk k ik ia 1}{max ξ≦}{max 22∑∑⋅kk kik ia ξ≦}max{22nd na i⋅=n AD ≦max{m ,n }AD =∞x A M .9. 只证范数的相容性公理及与向量2–范数的相容性. 设 A =(a ij )∈C n m ⨯, B =(b st )∈C l n ⨯,x =∈ξξξT 21),,,(n C n 且 A =ij ji a ,max , B =st ts b ,max , 则∑=≤≤≤≤=nk ktik lt m i AB11,1Gmaxb aml ≦}{max ,kt kik t i b a ml ∑≦}{max 22,∑∑⋅kkt kikti b a ml (Minkowski 不等式)≦ml n ab =))((b nl a mn =G G B A .∑∑===m i nk k ikAx1212ξa≦∑∑ik k ika2)(ξ≦∑∑∑⋅ikkk ik a )(22ξ (H ölder 不等式)≦∑∑⋅ikkna )(22ξ=mn A 2x=2G x A .10. 利用定理2.12得122H 2===nI UU U.11.A 1-=⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛---0110211214321cond 1(A )=225255111=⋅=-A A ; cond ∞(A )=10251=⋅=∞-∞A A .12.设x 是对应于λ的特征向量, 则A x x m m λ=.又设 v ⋅是C n 上与矩阵范数⋅相容的向量范数,那么vm vm v mx A x x ==λλ≦v m x A因 v x >0, 故由上式可得 mλ≦m A ⇒λ≦m m A .习 题 三1. 2c λc λλ))(2(+-=-A I , 当c λρ=)(﹤1时, 根据定理3.3, A 为收敛矩阵.2. 令S )N (=∑=N0)(k k A , )(lim N N S +∞→=S , 则 0)(lim lim )()()(=-=+∞→+∞→k k k k k S S A .反例: 设 A )(k =k⎪⎪⎭⎫ ⎝⎛0001k, 则因 ∑+∞=01k k发散, 故 ∑+∞=0)(k k A发散, 但)(lim k k A +∞→=O .3. 设 A =⎪⎪⎭⎫⎝⎛6.03.07.01.0, 则 )(A ρ≦=∞A 行和范数=0.9<1, 根据定理3.7,∑∞+=⎪⎪⎭⎫ ⎝⎛06.03.07.01.0k k=(I -A )1-=⎪⎪⎭⎫ ⎝⎛937432.4. 我们用用两种方法求矩阵函数e A : 相似对角化法. 22a λλ+=-A I , a -a i ,i =λ当 =λi a 时, 解方程组 (i a -A )x =0, 得解向量 p 1=(i, 1)T .当 λ=-i a 时, 解方程组 (i a +A )x =0, 得解向量 p 2=(-i, 1)T .令 P =⎪⎪⎭⎫⎝⎛-11i i , 则P 1-=⎪⎪⎭⎫ ⎝⎛-i 1i 1i 21, 于是 e A =P ⎪⎪⎭⎫⎝⎛-a ai 00i P 1-=⎪⎪⎭⎫ ⎝⎛a a a -a cos sin sin cos . 利用待定系数法. 设e λ=(2λ+a 2)q (λ)+r (λ), 且 r (λ)=b 0+b 1λ, 则由⎩⎨⎧=-=+-aaa b b a b b i 10i 10ei e i ⇒b 0=cos a , b 1=a1sin a .于是e A =b 0I +b 1A =cos a ⎪⎪⎭⎫ ⎝⎛11+a 1sin a ⎪⎪⎭⎫ ⎝⎛-a a =⎪⎪⎭⎫ ⎝⎛-a a a a cos sin sin cos . 后一求法显然比前一种方法更简便, 以后我们多用待定系数法. 设f (λ)=cos λ, 或 sin λ则有⎩⎨⎧=-=+a-a b b aa b b sini i sini i 1010 与 ⎩⎨⎧=-=+aa b b aa b b i cos i i cos i 1010 由此可得⎪⎩⎪⎨⎧-==a a b b sini i 010 与 ⎩⎨⎧==0i cos 10b ab 故 (a 2isini a )A =⎪⎪⎭⎫ ⎝⎛-0isini isini 0a a =sin A 与(cosi a )I =⎪⎪⎭⎫⎝⎛a acosi 00cosi =cos A .5. 对A 求得P = ⎪⎪⎪⎭⎫ ⎝⎛--013013111, P 1-=⎪⎪⎪⎭⎫ ⎝⎛-24633011061, P 1-AP =⎪⎪⎪⎭⎫ ⎝⎛-211根据p69方法二,e At =P diag(e t -,e t ,e t 2)P 1-=⎪⎪⎪⎭⎫⎝⎛+--++---------t t t t tt tt t t t t t t e 3e 3e 3e 30e 3e 3e 3e 30e e 3e 2e e 3e 4e 661222tsin A =P diag(sin(-1),sin1,sin2)P 1-=⎪⎪⎪⎭⎫⎝⎛--01sin 601sin 6001sin 42sin 21sin 22sin 42sin 616. D 3(λ)=101011----λλλ=2)1(-λλ, D 2(λ)=D 1(λ)=1, A ~J =⎪⎪⎪⎭⎫⎝⎛000010011.现设r (λ,t )=b 0+b 1λ+b 2λ2, 则有⎪⎩⎪⎨⎧==+=++1e 2e 021210b t b b b b b t t ⇒b 0=1, b 1=2e t -t e t -2, b 2=t e t -e t +1. 于是e t A =r (A , t )=b 0I +b 1A +b 2A 2=I +(2e `t -t e t -2)⎪⎪⎪⎭⎫⎝⎛100100011+(t e t -e t +1)⎪⎪⎪⎭⎫ ⎝⎛100100111=⎪⎪⎪⎭⎫ ⎝⎛-+--tt e 001e 101e e 1e e tt t t t同理,由⎪⎩⎪⎨⎧=-=+=++1sin 2cos 021210b t t b b t b b b ⇒b 0=1, b 1=t sin t +2cos t -2, b 2=1-t sin t -cos t . 将其代入cos A t =b 0I +b 1A +b 2A 2, 求出cos A t =⎪⎪⎪⎭⎫ ⎝⎛----t t t t t t t cos 001cos 10cos sin 11cos cos7. 设 f (A )=∑+∞=0k k A ka ,S N=∑=Nk k A 0k a .则 f (A )=N N S +∞→lim 并且由于(S N)T=T)(∑=N k k k A a =∑=Nk k k A 0T )(a所以, f (A T )=T )(lim N N S +∞→=f (A )T .8, (1) 对A 求得P =⎪⎪⎪⎪⎪⎭⎫ ⎝⎛1111, P 1-=P , J =⎪⎪⎪⎪⎪⎭⎫⎝⎛1111111 则有e t A =P ⎪⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛t t tt t tt ttt t t t t t t e e e e 2e e e 6e 2e 232eP 1-=⎪⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛t ttt t t tt t e e e 2e 60e e e 200e e 000e 232t t t t t t tsin A t =P ⎪⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛---t t t t t t t t t t t t t t t t sin cos sin sin 2cos sin cos 6sin 2cos sin 232P 1-=⎪⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛---t t t t t t t t t t t t t t t t sin cos sin 2cos 6sin cos sin 2sin cos sin 232cos A t =P ⎪⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛-----t t t tt t t t tt t t t t t t cos sin cos cos 2sin cos sin 6cos 2sin cos 232P=⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛-----t t t t t t t t t t t t t t t t cos sin cos 2sin 60cos sin cos 200cos sin 000cos 232(2) 对A 求出P =P 1-=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛0100100000100001, J =⎪⎪⎪⎪⎪⎭⎫⎝⎛--010212 则有e At =P ⎪⎪⎪⎪⎪⎭⎫⎝⎛---11e e e 222t t tt t P 1-=⎪⎪⎪⎪⎪⎭⎫⎝⎛---100010000e 000e e 222t t tt tsin A t =P ⎪⎪⎪⎪⎪⎭⎫⎝⎛--002sin 2cos 2sin t t tt tP 1-=⎪⎪⎪⎪⎪⎭⎫⎝⎛--0000000002sin 0002cos 2sin t t tt tcos A t =P ⎪⎪⎪⎪⎪⎭⎫ ⎝⎛1012cos 2sin 2cos t t t t P 1-=⎪⎪⎪⎪⎪⎭⎫⎝⎛10000100002cos 0002sin 2cos t t t t 9. (1) sin 2A +cos `2A =[)e (e i 21i i A A --]2=[)(e 21i i A A e -+]2=)e e e (e 41)e e e (e 41i 2i 2i 2i 2O O A A O O A A ++++--+---=e O =I(2) sin(A +2πI )=sin A cos(2πI )+cos A sin(2πI )=sin A [I -!21(2πI )2+!41(2πI )4-…]+cos A [2πI -!31(2πI )3+!51(2πI )5-…]= sin A [1-!21(2π)2+!41(2π)4-…]I +cos A [2π-!31(2π)3+!51(2π)5-…]I=sin A cos2π+cos A sin2π (3)的证明同上.(4) 因为 A (2πi I )=(2πi I )A ,所以根据定理3.10可得 e I A i π2+=e A e I πi 2=e A [I +(2πI )+!21(2πi I )2+!31(2πi I )3+…]=e A {[1-!21(2π)2+!41(2π)4-…]+i[2π-!31(2π)3+!51(2π)5-…]}I=e A {cos2π+isin2π}I =e A此题还可用下列方法证明:e I A πi 2+=e ⋅A e I i π2=e ⋅A P ⎪⎪⎪⎪⎪⎭⎫⎝⎛i π2iπ2πi 2e e e P 1-=e ⋅A PIP 1-=e A用同样的方法可证: e I A πi 2-=e A e I πi 2-.10. A T =-A , 根据第7题的结果得 (e A )T =e TA =e A -, 于是有e A (e A )T =e A e TA =e A A -=e O =I11. 因A 是Herm(i A )H =-i A H =-i A , 于是有e A i (e A i )H =e A i e A i -=e O =I12. 根据定理3.13, A 1-tt A e d d =e At , 利用定理3.14得 ⎰tA 0d e ττ=⎰-t A A 01d e d d τττ=A 1-τττd e d d 0A t ⎰=A 1-(e -At I ). 13. t d d A (t )=⎪⎪⎭⎫ ⎝⎛---t t t t sin cos cos sin , t d d (det A (t ))=t d d (1)=0, det(t d dA (t ))=1, A 1-(t )=⎪⎪⎭⎫ ⎝⎛-t t t t cos sin sin cos , t d d A 1-(t )=⎪⎪⎭⎫⎝⎛---t t t t sin cos cos sin14. ⎰t A 0d )(ττ=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛⎰⎰⎰⎰⎰⎰-00d 30d e 2d e d d e d e 002002002t t t t t t τττττττττττττ=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛--+---002301e e1311e e )1(e 212232t t t t t t t t 15. 取 m =2, A (t )=⎪⎪⎭⎫⎝⎛t t t 02, 则 A 2(t )=⎪⎪⎭⎫ ⎝⎛+22340t t t t , t d d (A (t ))2=⎪⎪⎭⎫ ⎝⎛+t t t t 2023423≠2A (t )t d dA (t )=⎪⎪⎭⎫⎝⎛+t t t t 2022423. 困为++==--21)]()[(d d)()]()[(d d )]()()([d d )]([d d m m A A A A A A A A A t t tt t t t t t t t t t m +)(d d)]([1t tt A A m -所以当(t d d A (t ))A (t )=A (t )t d dA (t )时, 有)(d d)]([)(d d )]([)(d d )]([)]([d d 111t tt t t t t t t t t A A A A A A A m m m m ---++= =m [A (t )])(d d1t tA m -16. (1) 设 B =(ij b )n m ⨯, X =(ij ξ)m n ⨯, 则 BX =(∑=nk kj ik 1ξb )m m ⨯,于是有tr(BX )=∑∑∑===++++nk km mk n k kj jk n k k k 11111ξξξb b bijBX ξ∂∂)tr(=ji b (i =1,2,…,n ;j =1,2,…,m ) ⎪⎪⎪⎭⎫ ⎝⎛=mn n m BX X b b b b 1111)(tr(d d=T B 由于 BX 与 T T T )(B X BX =的迹相同,所以T T T ))(tr(d d ))(tr(d d B BX XB X X == (2) 设A =(ij a )n n ⨯,f=tr(AX X T ), 则有⎪⎪⎪⎭⎫ ⎝⎛=nm mn X ξξξξ1111T ,AX =⎪⎪⎪⎪⎪⎭⎫⎝⎛∑∑∑∑k km nk k k nk km k k k k ξξξξa a a a 1111f =∑∑∑∑∑∑++++l kkm lk lm l k kj lk lj l kk lk l ξξξξξξa a a 11)]()([][∑∑∑∑∑∂∂⋅+⋅∂∂=∂∂=∂∂k kj lk l k ijlj kj lk ij lj l k kj lk lj ij ij ξξξξξξξξξξa a a f =∑∑+klj li kkj ik ξξa amn ij X ⨯⎪⎪⎭⎫⎝⎛∂∂=ξff d d =X A A X A AX )(T T +=+ 17. 设A =(ij a )m n ⨯, 则 F (x )=(∑∑∑===nk kn k nk k nk k k 1211,,,a a a 1k ξξξ ),且A d F F F x F nn n n n n n =⎪⎪⎪⎪⎪⎭⎫⎝⎛=⎪⎪⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=a a aa a a a a a 21222211121121d d d d d d d ξξξ 18. ()⎪⎪⎪⎭⎫⎝⎛---------=='t t tt t t tt t t t t t t t t tt AtAt A 222222222e 4e 3e 3e 6e 3e 6e 2e e e 4e e 2e 2e e e 2e e 4e e在上式中令t =0, 则有A =⎪⎪⎪⎭⎫ ⎝⎛---=133131113e OA19. A =⎪⎪⎪⎭⎫ ⎝⎛---502613803, x (0)=⎪⎪⎪⎭⎫⎝⎛111, A 的最小多项式为 2)1()(+=λλϕ. 记f (λ)=t λe ,并设f (λ)=g(λ))(λϕ+)(10λb b +, 则⎩⎨⎧==---tte e 110t b b b ⇒ tt --=+=e ,)1(10t b e t b 于是⎪⎪⎪⎭⎫ ⎝⎛--+=++=---t t t t t t t t 41026138041e e e )1(e t t t At A I , x (t )=Ate x (0)=⎪⎪⎪⎭⎫ ⎝⎛-++-t t t 6191121e t20. A =⎪⎪⎪⎭⎫ ⎝⎛--101024012, f (t )=⎪⎪⎪⎭⎫ ⎝⎛-1e 21t , x (0)=⎪⎪⎪⎭⎫ ⎝⎛-111, =)(λϕdet(λI -A)=23λλ-. 根据O A =)(ϕ,可得; 252423,,A A A A A A ===,….于是23232)!31!21()(!31)(!21)(e A A I A A A I At ++++=++++=t t t t t t=2)1(e A A I t t t --++=⎪⎪⎪⎭⎫⎝⎛---++--t t t e 1e e 210124021t t t t ttx (t )=⎪⎪⎪⎭⎫⎝⎛-=⎪⎪⎪⎭⎫ ⎝⎛+⎪⎪⎪⎭⎫ ⎝⎛-=⎪⎪⎪⎭⎫ ⎝⎛+=+⎰⎰-t t t t f e )1(11]02111[e ]d 021)0([]d )(e )0([e 00At t At tA At x e x ττττ习 题 四1. Doolite 分解的说明,以3阶矩阵为例: 11r 12r 13r 第1框 21l 22r 23r 第2框 31l 32l 33r 第3框 计算方法如下: (ⅰ) 先i 框,后i +1框,先r 后l .第1框中行元素为A 的第1行元素; (ⅱ)第2框中的j r 2为A 中的对应元素j a 2减去第1框中同行的21l 与同列的j r 1之积.第3框中的33r 为A 中的对应元素33a 先减去第1框中同行的31l 与同列的13r 之积,再减去第2框中同行的32l 与同列的23r 之积; (ⅲ)第2框中的32l 为A 中的对应元素32a 先减去第1框中同行的31l 与同列的12r 之积,再除以22r . 计算如下:1 3 02 -3 0 2 2 -6A =⎪⎪⎪⎭⎫ ⎝⎛--⎪⎪⎪⎭⎫ ⎝⎛600030031122012001 2.Crout 分解的说明,以3阶矩阵为例:11l 12u 13u 第1框 21l 22l 23u 第2框 31l 32l 33l 第3框(ⅰ) 先i 框,后i +1框.每框中先l 后r .第1框中的列元素为A 的第1列的对应元素;(ⅱ)第2框中的2i l 为A 中对应元素2i a 减去第1框中同行的1i l 与同列的12u 之积;(ⅲ)第2框中的23u 为A 中的对应元素23a 减去第1框中同行的21l 与同列的13u 之积,再除以22l .第3框中的33l 为A 中的对应元素33a 先减去第1框中同行的31l 与同列的13u 之积,再减去第2框中同行的32l 与同列的23u 之积.计算如下:1 3 02 -3 02 -6 -6A =⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎭⎫ ⎝⎛---100010031662032001 2. 先看下三角矩阵的一种写法:⎪⎪⎪⎭⎫⎝⎛333231222111000a a a a a a =⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛332211223211311121000000101001a a a a a a a a a , ii a ≠0 对本题中的矩阵A 求得Crout 分解为A =⎪⎪⎪⎪⎪⎭⎫⎝⎛--⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛--1002105452115240512005 利用下三角矩阵的写法对上面的分解变形可得A =⎪⎪⎪⎪⎪⎭⎫⎝⎛--⎪⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛--10021054521100051000512540152001 =⎪⎪⎪⎪⎪⎭⎫⎝⎛--⎪⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛--10021054521100510005100051000512540152001=⎪⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛--⎪⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛--10052510545251525405152005 3.对A 的第1列向量)1(β, 构造Householder 矩阵1H 使得 =)1(1βH 12)1(e β, 31C e ∈⎪⎪⎪⎭⎫ ⎝⎛=010)1(β, ⎪⎪⎪⎭⎫ ⎝⎛-=-01112)1()1(e ββ, u =⎪⎪⎪⎭⎫ ⎝⎛-=--01121212)1()1(12)1()1(e e ββββ⎪⎪⎪⎭⎫ ⎝⎛=-=1000010102T 1uu I H , ⎪⎪⎪⎭⎫⎝⎛=2301401111A H , ⎪⎪⎭⎫⎝⎛=23141A对1A 的第1列向量⎪⎪⎭⎫ ⎝⎛=34)2(β, 类似构造Householder 矩阵2H :⎪⎪⎭⎫ ⎝⎛-=--=3110122)2)2(12)2()2ββββe u , 21C e ∈, ⎪⎪⎭⎫ ⎝⎛-=-=4334512T 22uu I H ⎪⎪⎭⎫⎝⎛-=102512A H令12001H H H ⎪⎪⎭⎫⎝⎛=, 则有 ⎪⎪⎪⎭⎫ ⎝⎛-=100250111HA =R 并且⎪⎪⎪⎭⎫ ⎝⎛-⎪⎪⎪⎪⎪⎭⎫ ⎝⎛-=⎪⎪⎭⎫ ⎝⎛=⎪⎪⎭⎫ ⎝⎛==---1002501115453000153540001001T2T 112111R H H R H H R H A =QR4. 对A 的第1列向量⎪⎪⎪⎭⎫⎝⎛=202)1(β, 构造Givens 矩阵⎪⎪⎪⎪⎪⎭⎫⎝⎛-=210210102102113T , ⎪⎪⎪⎪⎭⎫⎝⎛=0022)1(13βT , ⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛-=1132221210220232322A O A T 对1A 的第1列向量⎪⎪⎪⎭⎫⎝⎛-=212)2(β, 构造 ⎪⎪⎪⎪⎭⎫ ⎝⎛-=3223131322~12T , ⎪⎪⎪⎭⎫ ⎝⎛=023~)2(12βT , ⎪⎪⎪⎪⎭⎫⎝⎛=34023723~112A T 令 ⎪⎪⎭⎫ ⎝⎛=12T12~1T O O T , 则有 ⎪⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛==34002372302323221312R A T T . 于是 QR R T T A =⎪⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛--==340023723023232232231213123403223121H13H 125. 设A =),,(i i 0i 0i 0i 1321ααα=⎪⎪⎪⎭⎫ ⎝⎛----, 对向量组321,,ααα施行正交化, 令⎪⎪⎪⎭⎫ ⎝⎛--==0i 111αβ, ⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛--+⎪⎪⎪⎭⎫ ⎝⎛=-=i 212i 0i 12i i 0i ],[],[1111222ββββααβ,⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛--=⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛-⎪⎪⎪⎭⎫ ⎝⎛---⎪⎪⎪⎭⎫ ⎝⎛--=--=323i232i 212i 3i 0i 1211i 0],[],[],[],[222231111333ββββαββββααβ于是⎪⎪⎪⎩⎪⎪⎪⎨⎧++=+-==3213212113i 212iβββαββαβα 写成矩阵行式K ),,(1003i 10212i 1),,(),,(321321321ββββββααα=⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡-= ⎪⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛⎪⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛----=32632316i 203i 612i 316i 21),,(321βββ 最后得A =K ⎪⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛----32632316i 203i 612i 316i 21=⎪⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛-⎪⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛----32006i 630212i 2316i 203i 612i 316i 21=QR 6. 令⎪⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛-==10005152********T T 则⎪⎪⎪⎪⎪⎭⎫⎝⎛=⎪⎪⎪⎭⎫⎝⎛⎪⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛-=011000520550114022011000515*******A T 再令⎪⎪⎪⎪⎪⎭⎫ ⎝⎛-==305061010610305132T T , ⎪⎪⎪⎪⎪⎭⎫⎝⎛-=3010305000061061612A T T 最后令⎪⎪⎪⎭⎫⎝⎛=0101000013T , R A T T T =⎪⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛-=00030103050610616123 A =⎪⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛-⎪⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛---=0003010305061061603056151302625230161H 3H 2H 3R T T T =QR 7. =)1(β(0, 1)T , 12)1(=β, u =2121)1(1)1(=--e e ββ(-1, 1)T ,H 1=⎪⎪⎭⎫⎝⎛=-01102T2uu I , H =⎪⎪⎭⎫⎝⎛1001H 则有HAH T =⎪⎪⎪⎭⎫⎝⎛⎪⎪⎪⎭⎫ ⎝⎛--⎪⎪⎪⎭⎫ ⎝⎛010100001111210121010100001=⎪⎪⎪⎭⎫ ⎝⎛--120111211, H 是Householder 矩阵.同理, 对)1(β, 取 c =0, s =1, T 12=⎪⎪⎭⎫⎝⎛-0110, T =⎪⎪⎭⎫ ⎝⎛12001T , 则 ⎪⎪⎪⎭⎫⎝⎛-⎪⎪⎪⎭⎫ ⎝⎛--⎪⎪⎪⎭⎫ ⎝⎛-=='-0101000011112101210101000011TAT T TA=⎪⎪⎪⎭⎫ ⎝⎛---120111211, T 是Givens 矩阵.8. 对 ⎪⎪⎭⎫⎝⎛=1612)1(β, 计算u =⎪⎪⎭⎫ ⎝⎛-=--2151202021)1(1)1(e e ββ, H =I -2uu T=⎪⎪⎭⎫ ⎝⎛-344351 令 Q =⎪⎪⎭⎫⎝⎛H 001, 则⎪⎪⎪⎭⎫⎝⎛=075075600200200TQAQ同理,对)1(β,为构造Givens 矩阵,令c =53, s =54, ⎪⎪⎪⎪⎭⎫ ⎝⎛-=5354545312T ,则当⎪⎪⎭⎫⎝⎛=12001T T 时,='T TA ⎪⎪⎪⎭⎫ ⎝⎛--075075600200200.1. (1) 对A 施行初等行变换⎪⎪⎪⎭⎫ ⎝⎛----100424201011200010321~⎪⎪⎪⎪⎭⎫ ⎝⎛---142000002102121100111201 S=,1420210011⎪⎪⎪⎪⎭⎫ ⎝⎛-- A =⎪⎪⎭⎫ ⎝⎛-⎪⎪⎪⎭⎫ ⎝⎛-2121101201422021(2) ⎪⎪⎪⎪⎪⎭⎫⎝⎛--------10001111010011110010111100011111~⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛-----11000000001100000210211110021021001 S=⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛---11000011021021021021, A =⎪⎪⎭⎫ ⎝⎛--⎪⎪⎪⎪⎪⎭⎫⎝⎛----1110000111111111(3)⎪⎪⎪⎪⎪⎭⎫ ⎝⎛1000126420100632100101264200016321~⎪⎪⎪⎪⎪⎭⎫⎝⎛---10100000010100000011000000016321 ⎪⎪⎪⎪⎪⎭⎫⎝⎛---=1010010100110001S, ()63212121⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=A 10. (1) ⎪⎪⎪⎭⎫⎝⎛=000000005T A A 的特征值是5,0,0. 分别对应特征向量321,,e e e ,从而V=I,),(11p V =∑=(5), 11AV U =∑1-=⎪⎪⎭⎫ ⎝⎛2151. 令,12512⎪⎪⎭⎫⎝⎛-=U ()21U U U =, 则I U A ⎪⎪⎭⎫⎝⎛=000005(2)⎪⎪⎭⎫⎝⎛=2112T A A 的特征值是,,1321==λλ对应的特征向量分别为TT11,11⎪⎪⎭⎫ ⎝⎛-⎪⎪⎭⎫⎝⎛.于是 ∑=⎪⎪⎭⎫⎝⎛1003, ⎪⎪⎪⎪⎭⎫⎝⎛-=21212121V =1V , 11AV U =∑1-=⎪⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛-06221612161取 ⎪⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛--=3131312U , 构造正交矩阵()21U U U ==⎪⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛---31062312161312161‘所以,A 的奇异值分解为T 001003V U A ⎪⎪⎪⎪⎭⎫ ⎝⎛=11. 根据第一章定理1.5, A A H 的特征值之和为其迹,而由第二章2.7 F-范数的定义A A A A A HH2F )tr(==的特征值之和=∑=ri i 12σ习 题 五1.设x =T 21),,,(n ηηη 为对应于特征值λ的单位特征向量,即(QD )x =λx两边取转置共轭:H H H H x Q D x λ=与上式左乘得2H H λ=Dx D x 即 22222221212n n ηηηd d d λ+++= ,由此立即有 2min iid ≤2λ≤2max i id从而i d imin ≤λ≤i d imax .后一不等式的另一证明:根据定理2.13,λ≤)(QD ρ≤2QD i d imax 最大特征值的H 22.11定理==D D D2. A 的四个盖尔园是 1G : 9-z ≤6, 2G : 8-z ≤2, 3G : 4-z ≤1, 4G : 1-z ≤1.由于4G 是一个单独的连通区域,故其中必有一个实特征值.321G G G ⋃⋃是连通区域,其中恰有三个特征值,因而含有一个实特征值 .3. A 的四个盖尔园:1G 1-z ≤2713, :2G 2-z ≤2713, :3G 3-z ≤2713, :4G 4-z ≤2713 是互相隔离的,并且都在右半平面,从而每个盖尔园中恰有一个特征值且为正实数.4.设 =λβαi +为A 的待征值,则有盖尔园k G ,使得k G ∈λ.若α≤0, 则kk a -α≤βαi )(+-kk a ≤k R 故 kk a +-)(α≤k R ,即 kk a ≤α+kk R ≤kk R , 这与A 是严格对角占优的条件矛盾.5. (1)当两个盖尔园的交集中含有两个特征值时; (2) 当两个盖尔园相切且切点是A 的单特征值时.6. A 的盖尔园 2:1-z G ≤3, 10:2-z G ≤2, 20:3-z G ≤10. 因1G 是与32G G ⋃分离的,故1G 中恰有一个实特征值∈1λ[-1, 5].A 的列盖尔园 :'1G 2-z ≤9, 10:'2-z G ≤4, 20:'3-z G ≤2. 因'3G 是与'2'1G G ⋃分离的,故 '3G 中恰有一个实特征值 ∈3λ[18, 22]. 选取 D =diag(1, 1,21), 则 1-DAD 的盖尔园 ''G 1 : 2-z ≤4, :''2G 10-z ≤3, :''3G20-z ≤5. 这三个盖尔园是相互独立的,故必然有∈1λ[-2, 6], ∈2λ[7, 13], ∈3λ[15, 25]与上面所得的结果对照可知利用Gerschgorin 定理,特征值的最隹估计区间为∈1λ[-1, 5], ∈2λ[7, 13], ∈3λ[18, 22]7. 因为det(λB -A )=)23)(2(422+-=----λλλλλλ所以广义特征值为1λ=2, 2λ=-32.分别求解齐次线性方程组0=-x A B )(1λ , 0=-x A B )(2λ可得对应于1λ与2λ的特征向量分别为⎪⎪⎭⎫⎝⎛121k (01≠k ), ⎪⎪⎭⎫ ⎝⎛-122k (02≠k ) 8. 先证明一个结果:若A 是Hermit 矩阵,n λλ,1分别是A 的最大、最小特征值,则)(m ax )(m ax 112x R x R x x =≠==λ, )(m ax )(m ax 12x R x R =≠==x x n λ事实上,Ax x x x x Axx x x x Axx x x x x H 1H 22H 220H H 002max 11max max )(max =≠≠≠===x R下证1λ>1μ, n λ>n μ. 令 Q =A -B , 则)(m ax m ax H H 1H 1122Qx x Bx x Ax x x x +====λ>Bx x x H 12max ==1μ( Q 正定,Qx x H >0 )同理可证 n λ>n μ.现在设 1<s <n , 则根据定理5.10及上面的结果,有)m ax (m in m ax m in H H H 1021Qx x Bx x Ax x x x P s +====λ>s x x P Bx x μ===H 1021max min 9. 显然,A B 1-的特征值就是A 相对于B 的广义特征值. 设为n λλλ,,,21 且j j j Bq Aq λ=, 0≠j q , j =1, 2, …,n 其中 n q q q ,,,21 是按B 标准正交的广义特征向量. 当 )(1A B -ρ<1时,对任意 x =0≠+++n n q q q c c c 2211)()(2211HH 22H 11H n n n n q q q A q q q Ax x c c c c c c ++++++==))((222111HH 22H 11n n n n n Bq Bq Bq q q q λλλc c c c c c ++++++ =2222211n n c c c λλλ+++ ≤i iλmax )(22221n c c c ++⋅=Bx x A B H 1)(-ρ<Bx x H反之,若对任意 x ≠0, Ax x H <Bx x H 成立,并且 )(1A B -=ρλ,Bq Aq λ=,0≠q ,则取 x=q , 于是有λ=Aq q H <1H =Bq q10. 若λ是BA 的特征值,q 是对应于λ的特征向量,即(BA )q =λq =λIq由此可知,λ是BA 的相对于单位矩阵I 的广义特征值 ,因此BAx x Ix x BAxx x R BA x x I x H 1H H 111222max max )(max )(======λ=)(maxH H 12Ax Bxx x x =≤)(max )(max H 1H 122Ax x Bx x x x == =)()(11A B λλ同理)(m in )(m in )(H H 1122Ax Bxx x x R BA x I x n ====λ≥)(m in )(m in H 1H 122Ax x Bx x x x == =)()(A B n n λλ11. 由于x ≠0时,12)()(==x x R x R ,从而5.24式等价于}0,1)(m in{m ax H 22)(2===-⨯∈x P x x R r n n P r C λ我们约定,下面的最小值都是对12=x 来取的. 令x =Qy , 则y y Ax x x R Qy P x P x P ΛH H H 2H 2H 2m in m in )(m in 0=====由于 n r n Q P ⨯-∈)(H 2C , 则在齐次线性方程组 0=Qy P H 2中,方程的个数小于未知量的个数,根据 Cramer 法则,它必有非零解. 设),,,,0,,0(~1n r r y ηηη +=,(1~2=y )为满足方程的解(容易证明这种形式的解必存在),则)(min ~min 22112~H ~H 2H 2n n r r r r y Q P y Q P y y ηληληλ ++=++==0Λ≤r λ 注意到 ⊆==}1~,~~{2H 2y y Q P y 0}1,{2H 2==y Qy P y 0,从而)(min H 2x R x P 0==)(min H 2y R Qy P 0=≤y y y R y Q P y Q P ΛH ~~~m in )~(m in H 2H 20===≤r λ 特别地,取),,(12n r q q P +=时,根据定理5.9)(min H 2x R x P r 0==λ故(5.24)式成立. 12. 我们约定:以下的最小值是对单位向量来取的,即证},1)(min{max H 22)(20C ===-⨯∈Bx P x x R r n n P r λ成立. 令 x =Qy , 则有y y x R BQy P B Bx P ΛHH2H 2m in )(m in === 设齐次线性方程组 0=BQy P H 2有形如 1~),,,,,0,,0(~21==+y y n r r ηηη 的解(不难证明这样的解一定存在),则因})({}~)(~{H 2H 200=⊆=y BQ P y y BQ P y所以)(min H 2x R B BxP ≤22112H ~~~min H 2n n r r r r y BQ P y y ηληληλ+++=++= Λ0≤r λ 特别地,取 ),,,(21H 2n r r q q q P ++=时,根据定理5.12可得r B Bx P x R λ==)(min H 20由此即知(5.44)成立.习 题 六求广义逆矩阵{1}的一般方法: 1)行变换、列置换法利用行变换矩阵S 和列置换矩阵P , 将矩阵A 化成SAP =⎪⎪⎭⎫⎝⎛O O K I r则。
矩阵理论第一章课后习题答案

1.按通常矩阵的加法及数与矩阵的乘法,下列数域F 上方阵集合是否构成F 上的线性空间:(1)全体形如⎪⎪⎭⎫ ⎝⎛b a-a 0的二阶方阵的集合; (2)全体n 阶对称(或反对称、上三角)矩阵的集合; (3){|0,}n n V X AX X F ⨯==∈(A 为给定的n 阶方阵).解:(1)设⎪⎪⎭⎫ ⎝⎛=111b a-a 0α⎪⎪⎭⎫ ⎝⎛-=222a 0b a β⎪⎪⎭⎫⎝⎛-=3330b a a γ ①αββα+=⎪⎪⎭⎫⎝⎛+⎪⎪⎭⎫ ⎝⎛-=⎪⎪⎭⎫ ⎝⎛+--+=⎪⎪⎭⎫⎝⎛-+⎪⎪⎭⎫ ⎝⎛=+111222212121222111b a -a 0a 00a 0b a -a 0b a b b a a a a b a ②)(0b a -a 0000a 0b a -a 0)(323232111321321321333212121333222111γβαγβα++=⎪⎪⎭⎫⎝⎛+--++⎪⎪⎭⎫ ⎝⎛=⎪⎪⎭⎫ ⎝⎛++---++=⎪⎪⎭⎫⎝⎛-+⎪⎪⎭⎫ ⎝⎛+--+=⎪⎪⎭⎫ ⎝⎛-+⎥⎦⎤⎢⎣⎡⎪⎪⎭⎫⎝⎛-+⎪⎪⎭⎫ ⎝⎛=++b b a a a a b b b a a a a a a b a a b b a a a a b a a b a③存在零向量V ∈0,使得对每个V a ∈,a a =⎪⎪⎭⎫⎝⎛=⎪⎪⎭⎫ ⎝⎛+⎪⎪⎭⎫ ⎝⎛=+111111b a -a 00000b a -a 00④对每个V a ∈,存在负向量a -,使得0b -a a -0b a -a 0)(111111=⎪⎪⎭⎫⎝⎛+⎪⎪⎭⎫ ⎝⎛=-+a a再令F y x ∈,⑤αα)(b a -a 0xyb xya -xya 0yb ya -ya 0b a -a 0)(111111111111xy xy x y x y x =⎪⎪⎭⎫⎝⎛=⎪⎪⎭⎫ ⎝⎛=⎪⎪⎭⎫ ⎝⎛=⎥⎦⎤⎢⎣⎡⎪⎪⎭⎫⎝⎛= ⑥αα=⎪⎪⎭⎫⎝⎛=111b a -a 011⑦βαβαx x b a xb xb xa xa xa xa b b a a a a x b a x x +=⎪⎪⎭⎫⎝⎛-+⎪⎪⎭⎫ ⎝⎛=⎪⎪⎭⎫⎝⎛+--+=⎪⎪⎭⎫ ⎝⎛+--+=⎥⎦⎤⎢⎣⎡⎪⎪⎭⎫⎝⎛-+⎪⎪⎭⎫ ⎝⎛=+222111212121212121222111a 0b a -a 000a 0b a -a 0)(⑧ya xa yb xb yaxa ya xa y x y x +=⎪⎪⎭⎫ ⎝⎛+⎪⎪⎭⎫ ⎝⎛=⎪⎪⎭⎫⎝⎛+--+=⎪⎪⎭⎫ ⎝⎛+=+111111*********yb ya -ya 0xb xa -xa 00b a -a 0)()(α所以全体形如⎪⎪⎭⎫⎝⎛b a -a 0的二阶方阵的集合构成F 上的线性空间。
西北工业大学矩阵论课件PPT第一章例题矩阵的相似变换

2100 3100 2100 3100
2100
例 求解一阶线性常系数微分方程组
ddt x1 2x1 x2 x3
ddt x2 x1 2x2 x3
d dt
x3
x1
x2
2 x3
解令
x
x1 x2 x3
,
dx dt
d dt
d dt
d dt
x1 x2 x3
, A
2 1 1
一次因式方幂的乘积, 并分别写出这些方幂
(相同的按出现的次数计数),称之为A的初等因子,
本题中A的初等因子为
2 和 ( 2)2 第三步:对每个初等因子( i )ri 作出 ri 阶
Jordan块
i
1
i
1
i
ri
ri
所有初等因子对应的Jordan块构成的Jordan矩阵 J
即是A的Jordan标准形。本题中A的Jordan标准形为
1 1
10
1 0 0,
1 0
3 0 ( 3)( 2), 1 2
3
1
1 2,
1
1 1
0 ( 2), 2
1 1 ( 2), 1 0 2,
11
1 2
1 0 ( 1)( 2)
1 2
所以
D2() 2
又 det(I A) ( 2)3,故
D3() ( 2)3
;
1 1 2
解 第一步:对 I A 用初等变换化为Smith
标准形:
3
I A 1
1
3
1
1
1
0
c2 ( 1) c1
1 0
1 2 2 4 4
0
r1( 3) r2
矩阵论的习题集

其中 aij = a ji = 1, (1 ≤ i ≤ n,1 ≤ j ≤ i ) ,其它元素为 0。 ′, ε 2 ′ ,ε3 ′,ε 4 ′ ] = [ε 1 , ε 2 , ε 3 , ε 4 ] A ,可得 6、[解]由 [ε 1 1 0 −1 ′, ε 2 ′ ,ε 3 ′ ,ε 4 ′] = A = [ε 1 , ε 2 , ε 3 , ε 4 ] [ε 1 0 0 1 0 0 0 1 0 2 = 0 0 1 0 0 0 0 3 0 4 0 1 1 2 2 1 0 3 2 1 2 0 4 1 3 2 2 = 3 1 4 1 1 0 3 1 0 2 0 0 1 1 0 2 0 0 1 0 0 0 0 2
2
3、对于 ∀B, C ∈ V 和 ∀λ ∈ F ,满足 BA = AB , CA = AC ,并且 A( B + C ) = AB + AC = BA + CA = ( B + C ) A , A( µB) = µAB = µBA = (µB ) A , 即 B + C ∈ V , µB ∈ V ,从而由第 1.2 节定理 1 可知,V 是 F n×n 的子空间。 满足 trB = 0 , 并且 tr ( B + C ) = trB + trC = 0 , 4、 对于 ∀B, C ∈ V 和 ∀λ ∈ R , trC = 0 , tr (λB) = λ tr ( B) = 0 ,从而由第 1.2 节定理 1 可知,V 是 R 2×2 的子空间。 1 0 0 1 0 0 dim V = 3 ,并且 V 的一组基为 , 0 − 1 和 。 0 0 1 0 5 、 对 于 ∀B, C ∈ V 和 ∀λ ∈ R , 满 足 B = B T , C = C T , 并 且 ( B + C ) T = B T + C T = B + C , (λB) T = λB T = λB ,从而由第 1.2 节定理 1 可知, V 是 R n×n 的子空间。 dim V = n(n + 1) ,并且 V 的一组基为 Vij = (a ij ) n×n , 2
第二章-矩阵(历年真题+答案)

D.
【解析】答案:B 【选择】 【201610】 【2 分】3.设 n 阶矩阵 A,B,C 满足 ABC=E,则 -1 -1 -1 A.A =B C B.A-1= C-1B-1 C.B-1=CA D.B-1=AC 【解析】答案:C 【填空】 【201610】 【2 分】
【解析】答案:-1 【填空】 【201610】 【2 分】
【解析】解:C2=ATBATB= AT(BAT)B=ATEB= ATB
2 1 1 C =A B= 0 0 0 2 1 1
2 T
【填空】 【201610】 【2 分】
0 1 0 1 4 3 -1 -1 【解析】解析:P= ,求得:P = ,则 P AP= 。 1 0 1 0 2 1
1 0 0 1 0 0 P 0 1 0 , Q= 0 0 1 2 0 1 0 1 0
使得 PA=B,BQ=E,即 PAQ=E
1 0 0 1 0 0 1 0 0 所以 A=P-1EQ-1=P-1Q-1= 0 1 0 0 0 1 = 0 0 1 -2 0 1 0 1 0 -2 1 0
足等式 AX B X ,求 X .
【解析】
(提示:另一种简便解法为:此题(A-E)X=B 的形式;构造分块矩阵( A-E,B ) 分 块 矩 阵 进 行 初 等 行 变 换 后 变 为 ( E , X ) 即 可 求 出
3 2 1 1 0 0 1 0 0 3 2 1 3 2 1 , A E,B 2 1 0 2 1 0 0 1 0 4 3 2 故 X= -4 -3 -2 ) 1 0 0 3 2 1 0 0 1 0 0 1 0 0 1
矩阵理论习题解答等材料

西南科技大学研究生试题单(B 卷)(2014级高等工程数学A)第一部分 矩阵理论(共32分)1、(8分)填空题(1)每个n 阶矩阵都相似于一个 矩阵。
(2)n nA C⨯∈,A 为正规矩阵的充要条件是A 对角形矩阵。
(3)正交变换在规范正交基下的矩阵是 矩阵。
(4)A 的最小多项式 A 的零化多项式。
2、(6分) 求4R 的子空间1234123412341234{(,,,)|0},{(,,,)|0}V a a a a a a a a W a a a a a a a a =-+-==+++=的交V W I 的一组基。
3、(8分) 已知111111,012A -⎛⎫ ⎪= ⎪ ⎪-⎝⎭计算5432()2822g A A A A A E =-++-。
4、(10分)求矩阵213121242A -⎛⎫⎪= ⎪ ⎪⎝⎭的Doolittle 分解和LDU 分解。
第二部分 数值分析(共36分)5、 (4分)解答下列各题 设函数2015201420131()5.2015!f x x x x =++,求差商0120142015[2,2,2,2]?f =L 6、(8分)设函数4()f x x =,不直接用拉格朗日插值公式,而用拉格朗日余项公式求出以1,0,1,2x =-为插值节点的三次插值多项式3().L x7、(8分)设有求积公式2120()(0)(1)(2)f x dx af a f a f ≈++⎰试确定系数012,,a a a 使上述公式的代数精度尽量高,且指出其代数精度。
8、(8分)已知方程组123123123102212100.51.931x x x x x x x x x --=⎧⎪-+-=⎨⎪--+=⎩ (1) 构造Jacobi 迭代法的迭代格式,迭代格式是否收敛?说明理由; (2) 取(0)(0,0,0)T x=,用上述迭代法来计算一步迭代值(保留小数点后4位)。
9、(8分)若求解初值问题为24,015(0)1x y y x y y ⎧'=-≤≤⎪⎨⎪=⎩, 试写出Euler 方法求解的迭代格式(0.2)h =,并计算(0.2),(0.4)y y 的值(保留小数点后至少8位)。
线性代数 第一章矩阵 参考答案

0 A2
0 A1
0 I A11r1 , A21r2 I 0 I 0 0 I
0 A11
A2 1 0
P31 习题 1.4 1.按上课要求做,则此题中行阶梯形答案不唯一,行最简形和标准形答案唯一
1 1 1 (1) 0 2 1 0 0 0
法一
2 1 1 B ( A 2 I ) A ,求出 ( A 2 I ) 1 1 1 4 3 4 2 3 3 8 B 1 5 3 1 1 0 2 9 1 6 4 1 2 3 2 12
4.解: 4 X
4 0 0 4 8 (3) 2 14 2 (4) 3 11 5 11 5 4 10 1 1 0 1 7.解: AB ; BA 1 2 0 0 1 2
1 0 (2) 0 0 1 0 (3) 0 0 1 0 (4) 0 0
1 1 0 0 3 2 0 0 1 1 0 0
1 1 0 0 3 1 1 0 2 1 0 0
1 0 1 ,0 1 2 1 0 0 1 1 0 0 0 1 1 1 , 0 0 0 1 0 0 0 0 7 0 1 5 1 , 0 1 1 0 0 0 0 2 1 1 0 5 1 ,0 1 1 0 0 0 0 0 0 0
(法二)
A1 X1 X 2 的逆阵为 B ,则有 0 X 3 X4 A1 X 1 X 2 I 0 0 X X I 0 4 3
A21 。 0
I 0 r1 r2 A2 0 I 0 0 A21 1 所以 A 1 0 A1 A1 0
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
《矩阵理论》第一二章 典型例题
一、 判断题
1.A n 为阶实对称矩阵,n
R x 对中的列向量, ||x |A
x =定义, ||x ||x 则为向量 的范数. ( )
2.设A n 为阶Hermite 矩阵,12,,,n λλλ 是矩阵A 的特征值,则22
2
1
||||n
m i i A λ==∑
.
( )
3. 如果m n A C ⨯∈,且0A ≠,()H
AA AA --=, 则2||||A A n -
=. ( )
4. 若设n x R ∈,则212||||||||||||x x x ≤≤. ( )
5. 设m n
A R
⨯∈的奇异值为12n σσσ≥≥≥ ,则222
1
||||n
i i A σ==∑. ( )
6. 设n n A C ⨯∈,且有某种算子范数||||⋅,使得||||1A <,则11||()||1||||
E A A -->
-,
其中E 为n 阶单位矩阵. ( )
7. 设2H A E uu =-(其中,E 为n 阶单位矩阵,2||||1n u C u ∈=且),则2
||||m A =
( )
8. 设n n A C ⨯∈为正规矩阵,则矩阵的谱半径2()||||r A A =. ( )
9.设n
n C
A ⨯∈可逆,n
n C
B ⨯∈,若对算子范数有1
||||||||1A
B -⋅<,则B A +可逆.
( )
10. 设A 为m n ⨯矩阵,P 为m 阶酉矩阵, 则PA 与A 有相同的奇异值. ( ) 11. 设n n
A C
⨯∈,且A 的所有列和都相等,则()r A A
∞
=. ( )
12. 如果12(,,,) T n
n x x x x C =∈,则1||||m in i i n
x x ≤≤=是向量范数. ( )
13. 设,n n A C ⨯∈则矩阵范数
m A
∞
与向量的1-范数相容. ( )
14、设n n
A C
⨯∈是不可逆矩阵,则对任一自相容矩阵范数 有1I A -≥, 其中I 为单位矩
阵. ( )
二、 设m n
A C
⨯∈,,||||ax ||ij i j
A a =
,证明:
(1)||||A 为矩阵范数; (2)||||A 为与向量2-范数相容.
三、 试证:如果A 为n 阶正规矩阵,且A x x λ=和Ay y μ=,其中,λμ≠,那么x 与y
正交.
四、 (1) 设(1)n n
A C
n ⨯∈>为严格对角占优矩阵,1122
(,,,)nn D diag a a
a = ,其中
(1,2,,)ii a i n = 为A 的对角元,E 为n 阶单位矩阵,则存在一个矩阵范数||||⋅使得
1
()1r E D
A --<.
(2) 设n n
A C
⨯∈, ε为任意给定的正数,()r A 为矩阵的谱半径。
证明:至少存在一个矩
阵范数||||A 使得||||().A r A ε≤+
五.设矩阵U 是酉矩阵, 12diag (,,,)n A a a a = , 证明: U A 的所有特征值λ满足不等式
{||}||{||}max min i i i
i
a a λ≤≤.
六. 设||||a ⋅是n n C ⨯上的相容的矩阵范数, 矩阵,B C 都是n 阶可逆矩阵, 且1
||||a B -及
1
||||a C
-都小于或等于1, 证明: 对任意矩阵n n
A C
⨯∈
||||||||b a A BAC =
定义了n n
C
⨯上的一个相容的矩阵范数.
七.设A 是可逆矩阵, λ是A 的一个特征值, 对于任意的算子范数||||⋅, 证明1
1||||||
A
λ-≥
.
八. 设A 是Hermite 矩阵()H
A
A =,且A
的特征值12n λλλ=== ,证明矩阵A 的
Rayleigh 商恒等于1λ.
九.已知n n C ⨯中的两种矩阵算子范数|| ||a 与|| ||b , 对于任意矩阵n n A C ⨯∈, 验证 ||||||||||||a b A A A =+是n n C ⨯中的相容矩阵范数.
十.设矩阵m n r A C ⨯∈的非零奇异值为12,,,r σσσ (0r >), 求证
1
2
21
||||().r
F i i A σ==∑
十一. 设矩阵n n A C ⨯∈可逆, 矩阵范数||||⋅是n C 上的向量范数||||v ⋅诱导出的算子范数, 令
()L x Ax =, 证明:
||||11
||||1
m ax ||()||||||||||m in ||()||v v v
x v
y L x A A
L y =-==⋅.
证明: 根据算子范数的定义, 有||||1
max ||()||||||x L x A ==,
1
1
1
||||1
||||1
||||||||111||||max
max
||||||||
||||
min ||||
min ||()||
min
||||
y A
x
x y y y y A x y A
Ay x Ay Ay L y y --=-≠≠==≠====
=
,
结论成立.
十二. 设矩阵n n A C ⨯∈为单纯矩阵, 证明: A 的特征值都是实数的充分必要条件是存在正定矩阵n n H C ⨯∈, 使得H A 为Hermite 矩阵.
十三. (1) 设矩阵()ij n n A a ⨯=, 则
,||||m ax ||a ij i j
A n a =⋅
是矩阵范数.
(2) 设,,,n x y p q C ∈, 矩阵H H A xp yq ,x y ,p q =+⊥⊥其中,求2
||||m A .。