矩阵分析所有习题及标准答案ppt课件

1、习题习题3-13-1已知已知A A CnCn n n是正定是正定HermiteHermite矩阵矩阵, , , ,Cn.Cn.定义内积定义内积 ( ( , , )=)= A A * *. .试证它试证它是内积是内积; ;写出相应的写出相应的C-SC-S不等式不等式: :Cauchy-Schwarz:Cauchy-Schwarz不等式：不等式： *,()(),;TAAAA *(, )( , );kk Ak *(, )()( , ) ( , );AAA *( , ) 0; ( , )0,0 ( A).A 因 正定|( , )| 11nnnnnni ijji ijji ijjijijijxa yxa

2、 xya y习题习题3-3(1)3-3(1)#3-3(1):#3-3(1):已知已知A= ,A= ,试求试求U U UnUn n n使使U U* *AU=RAU=R为为上三角矩阵上三角矩阵. .解解:det(:det( E-A)=(E-A)=( +1)3+1)3给出给出 =-1=-1是是A A的的3 3重特征值重特征值. .显然显然, , 1=(0,1,0)T1=(0,1,0)T是是A A的一个特征向量的一个特征向量. .作酉矩阵作酉矩阵V=(V=( 1,1, 2,2, 3),3), 2=(1,0,0)T,2=(1,0,0)T, 3=(0,0,1)T,3=(0,0,1)T,那么那么 V V*

3、*AV= AV= 子矩阵子矩阵A1A1的特征值仍是的特征值仍是-1,-1,对应的单位特征向量对应的单位特征向量是是 1=(-2/1=(-2/ 5,1/5,1/ 5)T,5)T,作作2 2阶酉矩阵阶酉矩阵W1=(W1=( 1,1, 2),2), 2=(1/2=(1/ 5,2/5,2/ 5)T,5)T,则则W1W1* *A1W1=A1W1=作作3 3阶酉矩阵阶酉矩阵W=diag(1,W1),U=VW,W=diag(1,W1),U=VW,那么那么 U U* *AU=AU=为上三角矩阵为上三角矩阵. .5026138035283,0063152083063111AA101011001010530110

4、06311*111WAW21010105521U=VW= 10001055001120125555 习题习题3-93-9#3-9:#3-9:若若S,TS,T分别为实对称分别为实对称, ,反实对称矩阵反实对称矩阵, ,则则A=(E+T+iS)(E-T-iS)-1A=(E+T+iS)(E-T-iS)-1为酉矩阵为酉矩阵. . 证证: :A A* *A=(E-T-iS)A=(E-T-iS)* *)-1(E+T+iS)-1(E+T+iS)* *(E+T+iS)(E-T-(E+T+iS)(E-T-iS)-1iS)-1 =(E+T+iS)-1(E-(T+iS)(E+(T+iS)(E-T-iS)- =(E+

5、T+iS)-1(E-(T+iS)(E+(T+iS)(E-T-iS)-1 1 =(E+T+iS)-1(E+T+iS)(E-T-iS)(E-T-iS)-1 =(E+T+iS)-1(E+T+iS)(E-T-iS)(E-T-iS)-1 =E =E注注: :可以不证可以不证 AAAA* *=E;=E;(E-(T+iS)(E+(T+iS)=(E+(T+iS)(E-(T+iS)(E-(T+iS)(E+(T+iS)=(E+(T+iS)(E-(T+iS) =(E+T+iS)(E-T-iS) =(E+T+iS)(E-T-iS)习题习题3-123-12设设A,BA,B均是正规矩阵均是正规矩阵, ,试证试证:A:A与

6、与B B酉酉相似的充要条件是相似的充要条件是A A与与B B的特征值相同的特征值相同 证证: :充分性：因为充分性：因为A,BA,B是正规矩阵是正规矩阵, ,所以存在所以存在U,VU,V UnUn n n 使得使得 A=Udiag(A=Udiag( 1,1, n)Un)U* *, B=Vdiag(, B=Vdiag( 1,1, n)Vn)V* *, , 其中其中 1,1, n n是是A,BA,B的特征值集合的特征值集合. .于是于是B=VUB=VU* *AUVAUV* *=W=W* *AW, W=UVAW, W=UV* * UnUn n n即得证即得证A A与与B B酉相似酉相似. . 必要性

7、必要性: :显然显然, ,因为因为, ,相似矩阵有相同的特征值相似矩阵有相同的特征值. . 习题习题3-133-13#3-13:#3-13:若若A A HnHn n,A2=A,n,A2=A,则存在则存在U U UnUn n n使得使得 U U* *AU=diag(Er,0),r=rank(A).AU=diag(Er,0),r=rank(A).证证: :存在存在U U UnUn n n使得使得 A=Udiag(A=Udiag( 1,1, n)Un)U* *, (, (* *) )其中其中 1,1, n n是是A A的特征值的任意排列的特征值的任意排列. . A2=A A2=A 和和 A2=Udi

8、ag(A2=Udiag( 1,1, n)Un)U* *Udiag(Udiag( 1,1, n)Un)U* * =Udiag( =Udiag( 12,12, n2)Un2)U* * i2=i2= i,i,即即 i i 0,1,i=1,n,.0,1,i=1,n,.取取 1,1, n n的排列使特征值的排列使特征值0 0全排在后面全排在后面, ,那么那么( (* *) )式即给出所需答案式即给出所需答案. .习题习题3-143-14#3-14:#3-14:若若A A HmHm n,A2=E,n,A2=E,则存在则存在U U UnUn n n使得使得 U U* *AU=diag(Er,-En-r).A

9、U=diag(Er,-En-r).证证: :存在存在U U UnUn n n使得使得 A=Udiag(A=Udiag( 1,1, n)Un)U* *, (, (* *) )其中其中 1,1, n n是是A A的特征值的任意排列的特征值的任意排列. . A2=E=Udiag(1,1)U A2=E=Udiag(1,1)U* * 和和 A2=Udiag(A2=Udiag( 1,1, n)Un)U* *Udiag(Udiag( 1,1, n)Un)U* * =Udiag( =Udiag( 12,12, n2)Un2)U* * i2=1,i2=1,即即 i=i= 1,i=1,n,.1,i=1,n,.取取

10、 1,1, n n的排列使特征值的排列使特征值1(1(设共有设共有r r个个) )全排在全排在前面前面, ,那么那么( (* *) )式即给出所需答案式即给出所需答案. .习题习题3-163-16#3-16:#3-16:设若设若A,BA,B HnHn n,n,且且A A为正定为正定HermiteHermite矩阵矩阵, , 试证试证:AB:AB与与BABA的特征值都是实数的特征值都是实数. .证证1:1:由定理由定理3.9.4,A1/23.9.4,A1/2是正定矩阵是正定矩阵, ,于是于是A-1/2(AB)A1/2=A1/2BA1/2=MA-1/2(AB)A1/2=A1/2BA1/2=M Hm

11、Hm n,n,即即ABAB相似于一个相似于一个HermiteHermite矩阵矩阵M.M. (AB)=(AB)= (M)(M) R,R,得证得证ABAB的特征值都是实数的特征值都是实数. .又又 A1/2(BA)A-1/2=A1/2BA1/2=MA1/2(BA)A-1/2=A1/2BA1/2=M HmHm n,n,即即BABA相似于一个相似于一个HermiteHermite矩阵矩阵M.M. (BA)=(BA)= (M)(M) R,R,得证得证BABA的特征值都是实数的特征值都是实数. .#3-16:#3-16:设若设若A,BA,B HmHm n,n,且且A A正定正定, ,试证试证:AB:AB

12、与与BABA的的特征值都是实数特征值都是实数. .证证2:2:由定理由定理3.9.1,PAP3.9.1,PAP* *=E,=E,那么那么PABP-1=PAPPABP-1=PAP* *(P(P* *)-1BP-1=(P)-1BP-1=(P* *)-1BP-1=M)-1BP-1=M HmHm n,n,即即ABAB相似于一个相似于一个HermiteHermite矩阵矩阵M.M. (AB)=(AB)= (M)(M) R,R,得证得证ABAB的特征值都是实数的特征值都是实数. .又又因因BABA的非零特征值与的非零特征值与ABAB的非零特征值完全相的非零特征值完全相同同, ,故故BABA的特征值也都是实

13、数的特征值也都是实数. .证证3:det(3:det( E-AB)=det(A(E-AB)=det(A( A-1-B)A-1-B) =det A det( =det A det( A-1-B)=0.A-1-B)=0.但但det A 0,det A 0,和和det(det( A-1-B)=0A-1-B)=0的根全为实数的根全为实数( (见见例例.1的相关证明的相关证明) )习题习题3-193-19设设A A是正定是正定HermiteHermite矩阵且矩阵且A A UnUn n,n,则则A=E A=E 证证: :存在存在U U UnUn n n使得使得 A=Udiag(A=Udi

14、ag( 1,1, n)Un)U* *, (, (* *) )其中其中 1,1, n n是是A A的特征值的任意排列的特征值的任意排列. . A A 是正定蕴含是正定蕴含 i0,i=1,n i0,i=1,n A A UnUn n n 蕴含蕴含| | i|=1,i=1,ni|=1,i=1,n 因而因而 i=1,i=1,ni=1,i=1,n A=Udiag( A=Udiag( 1,1, n)Un)U* *=UEU=UEU* *=UU=UU* *=E.=E.习题习题3-20 3-20 试证试证: :两个半正定矩阵之和是半两个半正定矩阵之和是半正定正定; ;半正定矩阵与正定矩阵之和是正定矩半正定矩阵与正

15、定矩阵之和是正定矩阵阵解解: : 设设A,BA,B HnHn n n 分别是半正定矩阵分别是半正定矩阵, ,正定矩阵正定矩阵. .那么那么A A* *=A&B=A&B* *=B =B (A+B) (A+B)* *=A+B =A+B HnHn n n x x Cn,xCn,x* *AxAx 0,x0,x* *BxBx 0 0 x x Cn,xCn,x* *(A+B)x(A+B)x 0 0 A+B A+B是半正定是半正定HermiteHermite矩阵矩阵. . 0 0 x x Cn,xCn,x* *AxAx 0,x0,x* *Bx0 Bx0 0 0 x x Cn,xCn,x* *

16、(A+B)x=x(A+B)x=x* *Ax+xAx+x* *Bx0Bx0 A+B A+B是正定是正定HermiteHermite矩阵矩阵. .习题习题3-223-22设设A,BA,B均是正规矩阵均是正规矩阵, ,试证试证:A:A与与B B相似的充要条件是相似的充要条件是A A与与B B酉相似酉相似证证: :因为因为A,BA,B是正规矩阵是正规矩阵, ,所以存在所以存在U,VU,V UnUn n n 使使得得 A=Udiag(A=Udiag( 1,1, n)Un)U* *, B=Vdiag(, B=Vdiag( 1,1, n)Vn)V* *, , 其中其中 1, 1, n,n, 1,1, n n

17、分别是分别是A,BA,B的特征值的特征值集合的任意排列集合的任意排列. .必要性：若必要性：若A A与与B B相似相似, ,那么那么 i=i= i,i=1,n,i,i=1,n,于是于是B=VUB=VU* *AUVAUV* *=W=W* *AW, W=UVAW, W=UV* * UnUn n n即得证即得证A A与与B B酉相似酉相似. . 充分性充分性: :显然显然, ,因为因为, ,酉相似必然相似酉相似必然相似. . 习题习题3-233-23设设A A* *=A.=A.试证试证: :总存在总存在t0,t0,使使得得A+tEA+tE是正定是正定;A-tE;A-tE是负定是负定证证: :因为因为

18、A A是是HermiteHermite矩阵矩阵, ,所以存在所以存在U U UnUn n n 使得使得 A=Udiag(A=Udiag( 1,1, n)Un)U* *, , 其中其中 1, 1, n n是是A A的特征值并且全为实数的特征值并且全为实数. .令令tMax|tMax| 1|,|1|,| n|,n|,于是于是,A+tE,A+tE是是HermiteHermite矩阵矩阵并且特征值全为正数，即得证并且特征值全为正数，即得证A+tEA+tE是正定是正定HermiteHermite矩阵矩阵. A. AtEtE是是HermiteHermite矩阵矩阵并且特征值全为负数，即得证并且特征值全为负

19、数，即得证A AtEtE是负定是负定HermiteHermite矩阵矩阵. .习题习题3-253-25#3-25:A#3-25:A* *=-A(A=-A(A SHnSHn n)n) U=(A+E)(A-E)- U=(A+E)(A-E)-1 1 UnUn n.n.(A(A SHnSHn n nA A E E的特征值全不为的特征值全不为0,0,从而从而A A E E可逆可逆) )解解: U: U* *=U-1=U-1 (A-E) (A-E)* *)-1(A+E)-1(A+E)* *=(A-E)(A+E)-1=(A-E)(A+E)-1 (-A-E)-1(-A+E)=(A-E)(A+E)-1 (-A-

20、E)-1(-A+E)=(A-E)(A+E)-1 (A+E)-1(A-E)=(A-E)(A+E)-1 (A+E)-1(A-E)=(A-E)(A+E)-1 (A-E)(A+E)=(A+E)(A-E) (A-E)(A+E)=(A+E)(A-E) A2-E=A2-E A2-E=A2-E 因最后一式恒成立因最后一式恒成立, ,得证得证U U* *=U-1,=U-1,从而从而 U=(A+E)(A-E)-1U=(A+E)(A-E)-1 UnUn n.n.习题习题3-263-26设设A A为正规矩阵特征值为为正规矩阵特征值为 1, 1, n.n.试证试证:A:A* *A A的特征值为的特征值为| | 1|2,

21、|1|2,| n|2.n|2.证证: :因为因为A A是正规矩阵是正规矩阵, ,所以存在所以存在U U UnUn n n 使得使得 A=Udiag(A=Udiag( 1,1, n)Un)U* *, , 其中其中 1, 1, n n是是A A的特征值的特征值. .于是于是, ,A A* *A=Udiag(|A=Udiag(| 1|2,|1|2,| n|2)Un|2)U* *. .因对角矩阵因对角矩阵diag(|diag(| 1|2,|1|2,| n|2)n|2)酉相似于酉相似于A A* *A,A,故故A A* *A A的特征值为的特征值为 | | 1|2,|1|2,| n|2n|2习题习题3-2

22、73-27#3-27(1):A#3-27(1):A* *A,AAA,AA* *都是半正定都是半正定HermiteHermite矩阵矩阵. . (2): (2):若若A A CmCm n,n,则则A A* *A,AAA,AA* *的非零特征值相的非零特征值相同同( (它们的谱可能不一样它们的谱可能不一样) )证证:(1): (A:(1): (A* *A)A)* *=A=A* *A,(AAA,(AA* *) )* *=AA=AA* *. . x x Cn,xCn,x* *(A(A* *A)x =(Ax)A)x =(Ax)* *Ax=(Ax,Ax)Ax=(Ax,Ax) 0.0. (2): (2):

23、对对AAAA* *的任意非零特征值的任意非零特征值 有有AAAA* *x=x= x,xx,x 0.0. 于是于是 A A* *A(AA(A* *x)=x)= (A(A* *x).x). 因因 x x 0,0,故故A A* *x x 0,0,从而得证从而得证AAAA* *的任意非零特的任意非零特征值征值 也是也是A A* *A A的非零特征值的非零特征值. . 同理可证同理可证:A:A* *A A的任意非零特征值的任意非零特征值 也是也是AAAA* *的非的非零特征值零特征值. .习题习题3-27(2)3-27(2)另一解法另一解法证证: :不难验证下列矩阵等式不难验证下列矩阵等式: : 因因S

24、= S= 可逆可逆, ,故故从而从而det(det( E-AAE-AA* *)=0)=0与与det(det( E-AE-A* *A)=0A)=0有相同非零有相同非零解解, ,得证得证AAAA* *与与A A* *A A有相同的非零特征值有相同的非零特征值. .AAAEAEAAAAAAAAEAEAAAnmnm*0000nmEAEAAASAAASAAA*1*000000习题习题3-283-28设设A A为正规矩阵为正规矩阵. .试证试证: :若若Ar=0,Ar=0,则则A=0.A=0.若若A2=A,A2=A,则则A A* *=A.=A.证证: :因为因为A A是正规矩阵是正规矩阵, ,所以存在所以

25、存在U U UnUn n n 使得使得 A=Udiag(A=Udiag( 1,1, n)Un)U* *, , 其中其中 1, 1, n n是是A A的特征值的特征值. .于是于是, ,Ar=Udiag(Ar=Udiag( 1r,1r, nr)Unr)U* *=0=0蕴涵蕴涵 ir=0,i=1,n.ir=0,i=1,n.后者又蕴涵后者又蕴涵 1=1= n=0. n=0. A=Udiag(0,0)U A=Udiag(0,0)U* *=0. =0. 假设假设 A2=A, A2=A, 那么那么 i2=i2= i,i=1,n. i,i=1,n. 后者又蕴后者又蕴涵涵 i=0i=0或或1, i=1,n,(

26、1, i=1,n,(即正规矩阵即正规矩阵A A的特征的特征值全为实数值全为实数).). A A* *=Udiag(=Udiag( 1,1, n)Un)U* *=A. =A. 习题习题3-303-30#3-30:#3-30:若若A A CnCn n,n,则则A A可唯一地写为可唯一地写为A=B+C,A=B+C,其中其中B B HnHn n,Cn,C SHnSHn n.n.证证: :存在性存在性 取取 B=(1/2)(A+AB=(1/2)(A+A* *),C=(1/2)(A-A),C=(1/2)(A-A* *), ), 则显然则显然B,CB,C分别是分别是HermiteHermite矩阵和反矩阵和

27、反HermiteHermite矩阵矩阵, ,并且满足并且满足A=B+C.A=B+C. 唯一性唯一性 假设假设 A=B+C,A=B+C,其中其中B B HnHn n,Cn,C SHnSHn n,n,那那么么A A* *=(B+C)=(B+C)* *=B=B* *+C+C* *=B-C.=B-C.于是于是 B=(1/2)(A+AB=(1/2)(A+A* *),C=(1/2)(A-A),C=(1/2)(A-A* *). ). 证毕证毕注注: :令令T=-iC,T=-iC,则则T T* *=iC=iC* *=i(-C)=T,=i(-C)=T,即即T T HnHn n.n.由此由此推出推出:A:A可唯一

28、地写为可唯一地写为A=B+iT,A=B+iT,其中其中B,TB,T HnHn n.n.习题习题3 3* *1 1试证：向量长度的齐次性试证：向量长度的齐次性#3#3* *1:1:试证试证证证: :令令 =(a1,an)T ,=(a1,an)T ,那么那么 k k =(a1,an)T =(a1,an)T ,nkkkCC2222111nnniiiiiikkakakak习题习题3 3* *2 2试证：在酉空间试证：在酉空间V V中成立广义中成立广义商高定理商高定理#3#3* *2:2:试证试证 1,1, k k V &(V &( i,i, j)=0,j)=0, i i j j 或等价

29、地或等价地( ( 1+1+ k,k, 1+1+ k)=(k)=( 1,1, 1)+(1)+( k,k, k k) )证证: :对对k k用归纳法证明用归纳法证明.k=2.k=2时时, ,有有 ( ( 1+1+ 2,2, 1+1+ 2)2=(2)2=( 1,1, 1)+(1)+( 1,1, 2)+(2)+( 2,2, 1)+(1)+( 2,2, 2)2) =( =( 1,1, 1)+(1)+( 2,2, 2)2)若若k-1k-1时结论成立时结论成立, ,那么那么 ( ( 1+1+ k-1,k-1, k)=0k)=0( ( 1+1+ k,k, 1+1+ k)=(k)=( 1+1+ k-k-1)+1

30、)+ k,(k,( 1+1+ k-1)+k-1)+ k)k) =( =( 1+1+ k-1,k-1, 1+1+ k-1)+(k-1)+( k,k, k)k) =( =( 1,1, 1)+(1)+( k,k, k)+(k)+( k,k, k)k)22211.kk 习题习题3 3* *3 3令令 1=(1,1,1,1)T,1=(1,1,1,1)T, 2=(3,3,-1,-2=(3,3,-1,-1)T, 1)T, 3=(-2,0,6,8)T,3=(-2,0,6,8)T,求求SpanSpan 1,1, 2,2, 33的标正基的标正基解解: : 1,1, 2,2, 3 3就是所要求的标正基就是所要求的标

31、正基. .11(1,1,1,1) ;T2122111( , )(2,2, 2, 2);( , )T 32313212211( ,)( , )( 1,1, 1,1) .( ,)( , )T 1111 1 1 1( , , , ) ;2 2 2 2T2221 111( , ,) ;2 222TT)21,21,21,21(333习题习题3 3* *5(i)5(i)用归纳法证明用归纳法证明1+3+5+(2n-1)2=n21+3+5+(2n-1)2=n2证证: :对对k k用归纳法证明用归纳法证明.k=1.k=1时结论显然成立时结论显然成立. . 若若n-1n-1时结论成立时结论成立1+3+5+(2n-

32、3)=(n-1)21+3+5+(2n-3)=(n-1)2那么那么 1+3+5+(2n-1)2 1+3+5+(2n-1)2 =1+3+5+(2n- =1+3+5+(2n-3)+(2n-1) 3)+(2n-1) =(n-1)2+(2n-1) =(n-1)2+(2n-1) =n2-2n+1+2n-1 =n2-2n+1+2n-1 =n2 =n2习题习题3 3* *6 6试证试证: : 为正规矩阵为正规矩阵解解所以所以A A为正规矩阵为正规矩阵. .易见易见:A:A不是对角阵且不是对角阵且A A* * A A和和A A* * -A-A因而因而,A,A不是不是HermiteHermite矩阵矩阵, ,也不

33、是反也不是反HermiteHermite矩阵矩阵. .0010 ,10 01iAiii 00000010102 00 0 10 0 10 0 2iiiAAiiiAAii 习题习题3 3* *7 7证明证明: :对任意正定矩阵对任意正定矩阵A,A,任意任意正整数正整数k k 都有正定矩阵都有正定矩阵S S 使使 Sk=ASk=A证证: :因为因为A A是正定矩阵是正定矩阵, ,所以存在所以存在U U UnUn n n 使使得得 A=Udiag(A=Udiag( 1,1, n)Un)U* *, , 其中其中 1, 1, n n是全为正数是全为正数. .令令S=Udiag(S=Udiag( 11/k

34、,11/k, n1/k)Un1/k)U* *, , 其中其中 i1/ki1/k是正数是正数 i i的的k k次算术根次算术根, ,也全为正也全为正数数. .由此推出由此推出: Sk=A,: Sk=A,并且并且S S酉相似于对角酉相似于对角元全为正数的对角矩阵元全为正数的对角矩阵, ,从而得证从而得证S S是正是正定定HermiteHermite矩阵矩阵习题习题4-1(1)4-1(1)4-1:4-1:求求 A= A= 的满秩分解的满秩分解. .解解1: A 1: A = C = C A=BC, B=(A5,A3,A1)= A=BC, B=(A5,A3,A1)=020210114011050121

35、311415213212011210114013212020210114015092111211221020210114011050习题习题4-1(1)4-1(1)4-1:4-1:求求 A= A= 的满秩分解的满秩分解. .解解2: A 2: A = C = C A=BC, B=(A1,A2,A3)= A=BC, B=(A1,A2,A3)=54511001011022201121311415213212011210114013212110501011012131131152212545151515258100010001习题习题4-1(2)4-1(2)4-1(2):4-1(2):求求 A= A

36、= 的满秩分解的满秩分解. .解解: A : A = C = C A=BC, B=(A1,A3)= A=BC, B=(A1,A3)=000001111001011131321111001011111101111001011121001习题习题4-24-2求求 A= A= 的奇异值分解的奇异值分解. .解解: : A A的奇异值是的奇异值是: : 2,1; 2,1; =diag(=diag( 2,1)2,1) AA AA* *的对应于特征值的对应于特征值2,12,1的单位特征向量是的单位特征向量是(1/(1/ 2,1/2,1/ 2,0)T, (1,0,0)T2,0)T, (1,0,0)T0100

37、0;10002121212121211UU100101)2)(1(100011011|;100011011*AEAAVUAV1001100100010001121212111A A的奇异值分解是的奇异值分解是: :10010010020100021212121UDVA100110021000A2121*11VU或习题习题4 4* *1A1A与与B B酉等价酉等价A A与与B B奇异值相同奇异值相同 必要性必要性: A=UBV : A=UBV AA AA* *=UBVV=UBVV* *B B* *U U* *=UBB=UBB* *U U* * BBBB* * AA AA* *与与BBBB* *有

38、相同的特征值集有相同的特征值集, ,得证得证A A与与B B有相有相同的奇异值集同的奇异值集. . 充分性充分性: :作作A,BA,B的奇异值分解的奇异值分解A=UDVA=UDV* *,B=U1DV1,B=U1DV1* *,D=diag(,D=diag( ,0),0),其中其中, , 是由它们的全部正奇异值组成的正对角是由它们的全部正奇异值组成的正对角矩阵矩阵. .于是于是U U* *AV=D=U1AV=D=U1* *BV1 BV1 A=(UU1 A=(UU1* *)B(V1V)B(V1V* *) )因酉矩阵的乘积因酉矩阵的乘积 UU1UU1* *,V1V,V1V* * 仍为酉矩阵仍为酉矩阵,

39、 ,故上故上式表明式表明A A酉等价于酉等价于B.B.习题习题4 4* *2 24 4* *2: 2: 设设A A CrmCrm n,Un,U UmUm m,Vm,V UnUn n n使使B=UB=U* *AV=diag(AV=diag( ,0),0), =diag(b1,br=diag(b1,br), (), (* *) ) 那么那么|b1|,|br|b1|,|br|为为A A的全部正奇异值的全部正奇异值. . 证证: U: U* *AAAA* *U=BBU=BB* *=diag(=diag(* *,0) ,0) 写成写成 2 2不对！不对！ =diag(|b1|2,|br|2,0,0)=d

40、iag(|b1|2,|br|2,0,0) AA AA* * |b1|,|br| |b1|,|br|为为A A的全部正奇异值的全部正奇异值. .奇异值分解定理另一奇异值分解定理另一( (更强更强) )表述表述定理定理: : 令令 1,1, r r为为A A CrmCrm n n的全部正奇异值的全部正奇异值; ; =diag(=diag( 1,1, r),r),则有则有U U UmUm m,Vm,V UnUn n n使使 U U* *AV= =DAV= =D CrmCrm n n ( (* *) ) 反之反之, ,若有若有U U UmUm m,Vm,V UnUn n n使使( (* *) )成立成

41、立, ,其中其中 =diag(d1,dr),=diag(d1,dr), i,di0,i,di0,则则d1,drd1,dr为为A A的全部正奇异值的全部正奇异值.(.(奇异值分解的某种唯一奇异值分解的某种唯一性性) )证证: AA: AA* *=U V=U V* *V UV U* *=U U=U U* * diag(d12,dr2,0,0) diag(d12,dr2,0,0) d1,dr d1,dr为为A A的全部正奇异值的全部正奇异值. .注注: :后半部等价于补充题后半部等价于补充题4 4* *2.2.00000020000004 4* *3 3已知已知A A奇异值求奇异值求AT,AAT,A

42、* *,A-1,A-1的奇的奇异值异值补充题补充题4 4* *3: 3: 令令 1,1, r r为为A A CrmCrm n n的全部正奇的全部正奇异值异值; ; =diag(=diag( 1,1, r),r),则有则有U U UmUm m,Vm,V UnUn n n使使 A=U VA=U V* *=Udiag(=Udiag( ,0)V,0)V* * ( (* *) ) 易见易见 A A* *=Vdiag(=Vdiag( ,0)U,0)U* *AT=(Udiag(AT=(Udiag( ,0)V,0)V* *)T=(V)T=(V* *)Tdiag()Tdiag( ,0)UT,0)UT 1,1,

43、r r为为A A* *,AT, ,AT, 的全部正奇异值的全部正奇异值( (利用利用奇异值分解定理的更强表述奇异值分解定理的更强表述).).A-1=(UA-1=(U V V* *)-1=V)-1=V -1U-1U* *=Vdiag(=Vdiag( 1-1,1-1, n-1)Un-1)U* * 1-1,1-1, n-1n-1为为A-1A-1的全部正奇异值的全部正奇异值. .000*000VUAA习题习题#5-1(2)#5-1(2)试证试证: : x,yx,y V,xV,x yy |x-y|.|x-y|.证：首先证：首先x=(x-y)+yx=(x-y)+y x-y+yx-y+y x-y x-y x

44、-y.x-y.其次其次x-y=-(y-x)=y-xx-y=-(y-x)=y-x y-x= -(x-y)y-x= -(x-y) x-y x-y |x-y|.|x-y|.此外此外 x+y=x-(-y)x+y=x-(-y) |x-y|=|x-y|x-y|=|x-y| x x yy |x-y|.|x-y|.习题习题#5-2#5-2试证试证A= n maxi,j|aij|A= n maxi,j|aij|是矩是矩阵范数阵范数 A=(aij)A=(aij) CnCn n n证证: : 非负性非负性, ,齐次性显然齐次性显然 三角不等式三角不等式: :A+B= n maxi,j|aij+bij|A+B= n

45、maxi,j|aij+bij| n maxi,j|aij|+n n maxi,j|aij|+n maxi,j|bij|=A+Bmaxi,j|bij|=A+B 相容性相容性: :AB= n maxi,j|ai1b1j+ainbnj|AB= n maxi,j|ai1b1j+ainbnj| n2 maxi,t|ait| maxtj|btj| n2 maxi,t|ait| maxtj|btj| =n maxi,j|aij|(n =n maxi,j|aij|(n maxi,j|bij|)=ABmaxi,j|bij|)=AB习题习题#5-3#5-3设设是诱导范数是诱导范数detAdetA 0 0 试证试证

46、: : A A CnCn n,A-1n,A-1 A-1A-1和和 A-1-1= minxA-1-1= minx 0(Ax/x).0(Ax/x).证证: 1=E=AA-1: 1=E=AA-1 AA-1AA-1 detA detA 0 0 A0 A0 A-1 A-1 1/A=A-1.1/A=A-1.A-1= maxxA-1= maxx 0(A-1x/x)0(A-1x/x) = maxy = maxy 0(y/Ay) y=A-1x0(y/Ay) y=A-1x 0 0 x x 0 0 = maxy = maxy 0(1/(Ay/y)0(1/(Ay/y) = 1/miny = 1/miny 0(Ay/y

47、)0(Ay/y) A-1-1= minx A-1-1= minx 0(Ax/x).0(Ax/x).同一向量的三种范数之间的大小关系同一向量的三种范数之间的大小关系习题习题#5-4:#5-4:对对n n维线性空间的任意向量维线性空间的任意向量x x成立成立 xx x2 x2 x1 x1 nxnx nx2 nx2 nx1 nx1 n2xn2x 证证: :xx = max|x1|,|xn|= max|x1|,|xn| ( ( i=1n|xi|2)1/2 = x2i=1n|xi|2)1/2 = x2 (|x1|+|xn|)2)1/2 = (|x1|+|xn|)2)1/2 = x1x1 n max|x1

48、|,|xn| = n max|x1|,|xn| = nxnx 习题习题#5-6A#5-6A CnCn n n是正定矩阵是正定矩阵,x,x CnCn 证明证明:x=(x:x=(x* *Ax)1/2 Ax)1/2 是向量范数是向量范数. .解解1:1:因因A A是正定是正定HermiteHermite矩阵矩阵A,A,故存在可逆故存在可逆矩阵矩阵B B使得使得A=BA=B* *B.B.则则x x的上述表示式可写的上述表示式可写为为: :x=(xx=(x* *Ax)1/2 =(Bx)Ax)1/2 =(Bx)* *(Bx)1/2 (Bx)1/2 =Bx2=Bx2 其中其中2 2 是向量是向量2-2-范数

49、范数. .再注意可逆矩阵再注意可逆矩阵B B的性质的性质:x=0 :x=0 Bx=0, Bx=0,即可直接推出非即可直接推出非负性负性. .kx=B(kx)2=|k|Bx2=|k|xkx=B(kx)2=|k|Bx2=|k|x 推出齐次性推出齐次性; ;三角不等式则由下式推出三角不等式则由下式推出: :x+y=B(x+y)2x+y=B(x+y)2 Bx2+By2Bx2+By2#5-6 A#5-6 A正定正定, ,定义定义x x Cn,x=(xCn,x=(x* *Ax)1/2Ax)1/2试证试证: : 是一个向量范数是一个向量范数. .解解2:2:验证矩阵范数验证矩阵范数3 3条公理成立条公理成立

50、. .前两条显然成立前两条显然成立. .只须证三角不等式只须证三角不等式. . x+y2=(x+y) x+y2=(x+y)* *A(x+y)=(xA(x+y)=(x* *+y+y* *)(Ax+Ay)(Ax+Ay) =x =x* *Ax+yAx+y* *Ay+xAy+x* *Ay+yAy+y* *AxAx =x2+y2+2Re(x =x2+y2+2Re(x* *Ay)Ay)令令B B为为A A的正定的正定HermiteHermite平方根平方根:A=BB,:A=BB,那么那么 x x* *Ay=xAy=x* *BBy=(Bx)BBy=(Bx)* *(By)=(Bx,By) (By)=(Bx,B

51、y) 标准内积标准内积由由Cauchy-SchwarzCauchy-Schwarz不等式不等式 |2Re(x|2Re(x* *Ay)| Ay)| 2|x 2|x* *Ay|Ay| 2(Bx,Bx)1/2(By,By)1/2 = 2(Bx,Bx)1/2(By,By)1/2 = 2xy2xy x+y2 x+y2 (x+y)2, (x+y)2, 得证所需结论得证所需结论. .习题习题#5-7#5-7试找一个收敛的试找一个收敛的2 2阶可逆方阵序列其极限矩阶可逆方阵序列其极限矩阵不可逆阵不可逆 解解: :下列矩阵序列满足所提条件下列矩阵序列满足所提条件: :AkAk的行列式都大于的行列式都大于0,0,

52、故可逆故可逆, ,但极限矩阵是但极限矩阵是行列式不为行列式不为0 0的不可逆矩阵的不可逆矩阵: :,.3 , 2 , 1;11kAkk01)(lim1lim1kkkkA习题习题#5-9 #5-9 计算矩阵幂级数计算矩阵幂级数 试计算幂级数试计算幂级数: : 解解1:1:利用利用JordanJordan标准形标准形B=Pdiag(.5,-.3)P-1,P=B=Pdiag(.5,-.3)P-1,P=解解2:2:利用谱半径小于利用谱半径小于1 1的矩阵性质的矩阵性质, , (B)=0.51.(B)=0.51. E+ E+ k=1k=1 Bk=(E-B)-1=Bk=(E-B)-1= 答案是答案是 k=

53、1k=1 Bk =Bk =解解3: 3: 也可利用也可利用 (B)(B) B1=BB1=B =0.91=0.91=R.:21=R. 所以所以, ,此矩阵幂级数发散此矩阵幂级数发散. .(2):(2):解解: :因因A1=MAX0.9,0.8,0.9=0.91=R,A1=MAX0.9,0.8,0.9=0.91=R=MAX1.1,0.9,0.6=1.11=R2011,0AAkk04.02.01.01.07.08.03.00,0AkAkk补充题补充题5 5* *5 5下列矩阵幂级数是否绝对收敛下列矩阵幂级数是否绝对收敛? ?(3)(3)解解1:1:此矩阵幂级数对应幂级数的收敛半径此矩阵幂级数对应幂级

54、数的收敛半径因因AA =MAX1.7,1.9=1.9R,=MAX1.7,1.9=1.9R.(A1=2.3R发散？）发散？）解解2:2:此矩阵幂级数等价于此矩阵幂级数等价于而的矩阵幂级数绝对收敛而的矩阵幂级数绝对收敛(B(B =0.951).=0.951).8 .01 .15 .02 .1,20AAkkk2lim12121kkkR4 . 055. 05 . 26 . 0,22000BBAAkkkkkkk习题习题#6-5#6-5求已知矩阵求已知矩阵A A的最小多项式的最小多项式 知知 A= A= 解解I:I:解解II: II: A(A( )=dn()=dn( )=Dn()=Dn( )/Dn-1()

55、/Dn-1( ) ) =( =( -1)3/(-1)3/( -1)=(-1)=( -1)2 -1)2 1114320013) 1(111432001det)det()(AED22) 1()(,0, 0321532110321532110)(AEAEA故因11143,1142,1132,*01gcd1nD习题习题#6-5#6-5求已知矩阵求已知矩阵A A的最小多项式的最小多项式 知知 A= A= 解解I:I:因因 A+E A+E 和和 A-2EA-2E都都 0,0,并且并且(A-2E)(A+E)= 0,(A-2E)(A+E)= 0,故故 A(A( )=()=( -2)(-2)( +1)+1)01

56、11011102) 1)(2(1001) 2(1111111) 2(111111|AE0111111111211121112)(2(EAEA习题习题#6-5#6-5求已知矩阵求已知矩阵A A的最小多项式的最小多项式 知知 A= A= 解解II:II: A(A( )=dn()=dn( )=Dn()=Dn( )/Dn-1()/Dn-1( ) ) =( =( -2)(-2)( +1)2/(+1)2/( +1)=(+1)=( -2)(-2)( -1)-1)011101110112001101011001EA 习题习题#6-6#6-6已知矩阵已知矩阵A A求求f(A)f(A)的的JordanJordan

57、表示表示式式 知知 A= A= 解解: :因因 (A-E)(A-2E) (A-E)(A-2E) 0, 0,故故 A(A( )=()=( -1)(-1)( -2)2,-2)2,从而得从而得A A的初等因子为的初等因子为: : -1,(-1,( -2)2.-2)2.设变设变换矩阵为换矩阵为P=(P=( 1,1, 2,2, 3),3),那么那么 A(A( 1,1, 2,2, 3)=(3)=( 1,1, 2,2, 3)3) 给出给出(A-E)(A-E) 1=0,(A-2E)1=0,(A-2E) 2=0,(A-2E)2=0,(A-2E) 3=3= 2 2 解这些方程组求得解这些方程组求得 P=(P=(

58、1,1, 2,2, 3)=3)=2001200012001210002) 2)(1(200121001|AE100011001习题习题#6-6#6-6续续: : :)(005/1)()()(00)(42121214141arctgarctgarctgarctgarctgAarctg11100011001)2()2()2() 1 (100011001)2()2()2() 1 ()(ffffPffffPAf)2(00)2()2()2() 1 (00) 1 (100011001)2(00)2()2() 1 (00) 1 (fffffffffffttttttAtAeeeeeeeeeeeeee22222

59、22200)00,00)00注注: : f(x)=arctg(x/4) f(x)=arctg(x/4) f f (x)= (x)= 51204)2(,1644114122fxx补充题补充题#6#6* *1 1 已知已知A A和和p(p( ),),求求p(A)p(A)知知 A= A= p(p( )=)= 4-24-2 3+3+ -1,f(-1,f( )=)= 12-412-4 11+411+4 10-10- +3+3解解I:I:易见的特征多项式易见的特征多项式D(D( )=()=( -2)3.-2)3.(A-2E)2=0&A-2E(A-2E)2=0&A-2E 0 0 A(A( )

60、=()=( -2)2=-2)2= 2-42-4 +4+4p(p( )=()=( 2+22+2 +4)(+4)( 2-42-4 +4)+9+4)+9 -17-17 p(A)= 0+9A-17E = p(A)= 0+9A-17E =f(f( )=)= 10(10( 2-42-4 +4)-+4)- +3+3 p(A)= 0-A+3E = p(A)= 0-A+3E =1099989001311111002011121001解解II:II:由由D(D( )=()=( -2)3.-2)3.和和 A(A( )=()=( -2)2=-2)2= 2-2-4 4 +4+4A A有有JordanJordan标准形标准形 并有变换矩阵并有变换矩阵P P满足满足