矩阵的特征值和特征向量二次型优秀PPT.ppt

《矩阵的特征值和特征向量二次型优秀PPT.ppt》由会员分享，可在线阅读，更多相关《矩阵的特征值和特征向量二次型优秀PPT.ppt（58页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、矩阵的特征值和特征向量二次型现在学习的是第1页，共58页实验目的实验目的1、学会用MATLAB软件求矩阵的特征值和特征向量2、学会用MATLAB软件将二次型化为标准型3、通过用MATLAB软件编程来判断二次型的正定性现在学习的是第2页，共58页一、特征值与特征向量一、特征值与特征向量 现在学习的是第3页，共58页现在学习的是第4页，共58页 其中：D为由特征值构成的对角阵，V为由特征向量作为列向量构成的矩阵。且使 AV=VD 成立用Matlab计算特征值和特征向量的命令如下：d=eig(A)仅计算A的特征值(以向量形式d存放)V,D=eig(A)trace(A)计算矩阵A的迹现在学习的是第5页

2、，共58页例例1 1：求方阵的特征值、特征向量和迹解：解：A=2 2-2;2 5-4;-2-4 5;V D=eig(A)trace(A)现在学习的是第6页，共58页V=-0.2981 0.8944 0.3333 -0.5963 -0.4472 0.6667 -0.7454 0 -0.6667D=1.0000 0 0 0 1.0000 0 0 0 10.0000 trace(A)ans=12现在学习的是第7页，共58页答：答：特征值为：现在学习的是第8页，共58页现在学习的是第9页，共58页例例2 2：求方阵的特征值、特征向量和迹解：解：A=4 6 0;-3-5 0;-3-6 1;V D=eig

3、(A)trace(A)现在学习的是第10页，共58页二、矩阵的相似对角化二、矩阵的相似对角化现在学习的是第11页，共58页例例3 3：判断下列方阵是否可对角化。若可对角 化，求出可逆阵P，使P-1AP为对角阵。现在学习的是第12页，共58页解解(1)(1)：A=4 6 0;-3-5 0;-3-6 1;V D=eig(A)rank(V)ans=3答：答：A可对角化，且可对角化，且V=0 0.5774 -0.8944 0 -0.5774 0.4472 1.0000 -0.5774 0D=1 0 0 0 -2 0 0 0 1现在学习的是第13页，共58页 A=0 1 0;-1 2 0;-1 1 1;

4、V D=eig(A)rank(V)ans=2答：答：A不可对角化。不可对角化。解解(2)(2)：V=0 0.6325 0.4511 0 0.6325 0.4511 1.0000 0.4472 0.7701D=1 0 0 0 1 0 0 0 1现在学习的是第14页，共58页 下述函数可用来判断矩阵是否可对角化，若可对角化返回1，否则返回0。现在学习的是第15页，共58页function y=trigle(A)%可对角化返回1 1，否则返回0 0。y=1;c=size(A);if c(1)=c(2)y=0;return;ende=eig(A);n=length(A);while 1 if isem

5、pty(e)return;endd=e(1);f=sum(abs(e-d)10*eps);g=n-rank(A-d*eye(n);if f=g y=0;return;end e(find(abs(e-d)10*eps)=;end现在学习的是第16页，共58页function y=trigle(A)%可对角化返回1 1，否则返回0 0。y=1;c=size(A);if c(1)=c(2)y=0;returnende=eig(A);n=length(A);while 1 if isempty(e)%若为空阵则为真若为空阵则为真现在学习的是第17页，共58页return;endd=e(1);f=su

6、m(abs(e-d)10*eps);%特征值d的代数重数 g=n-rank(A-d*eye(n);%特征值d的几何重数 if f=g y=0;return;end e(find(abs(e-d)A=4-3 1 2;5-8 5 4;6-12 8 5;1-3 2 2 trigle(A)ans=0 A=1 1 1 1;1 1 1 1;1 1 1 1;1 1 1 1;trigle(A)ans=1答：答：A不不可对角化。可对角化。P D=eig(A)解解(2)(2)：现在学习的是第20页，共58页答：答：A可对角化，且可对角化，且P=-0.5000 0.2113 0.2887 0.7887 0.5000

7、 0.7887 -0.2887 0.2113 0.5000 -0.5774 -0.2887 0.5774 0.5000 0 0.8660 0D=-2.0000 0 0 0 0 2.0000 0 0 0 0 2.0000 0 0 0 0 2.0000现在学习的是第21页，共58页二、二、二次型化标准型二次型化标准型现在学习的是第22页，共58页现在学习的是第23页，共58页例例5 5：判断下列矩阵是否对称A=1 3 4 6;3 7 9 5;4 9 4 1;6 5 1 0;B=A;if(A=B)fprintf(A是对称矩阵是对称矩阵)else if(A=-B)fprintf(A是反对称矩阵是反对称

8、矩阵)else fprintf(A既不是对称矩阵，也不是反对称矩阵既不是对称矩阵，也不是反对称矩阵)endendA A是对称矩阵是对称矩阵解：解：现在学习的是第24页，共58页现在学习的是第25页，共58页Matlab中二次型化成标准形的命令为：P,T =schur (A)其中：A 二次型矩阵(即实对称矩阵)；T 为 A 的特征值所构成的对角形矩阵；P 为 T 对应的正交变换的正交矩阵，P 的列向量为 A的特征值所对应的特征向量 P,T =eig (A)现在学习的是第26页，共58页例例6 6：求一个正交变换，将二次型解：解：该二次型所对应的矩阵为化成标准形现在学习的是第27页，共58页 A=

9、1 1 0 1;1 1 1 0;0 1 1 1;-1 0 1 1;P,T =schur(A)P=-0.5000 0.7071 0.0000 0.5000 0.5000 -0.0000 0.7071 0.5000 0.5000 0.7071 0.0000 -0.5000 -0.5000 0 0.7071 -0.5000 P,T =eig(A)现在学习的是第28页，共58页T=-1.0000 0 0 0 0 1.0000 0 0 0 0 1.0000 0 0 0 0 3.0000答：所作的正交变换为：二次型的标准型为：现在学习的是第29页，共58页例例7 7：求一个正交变换，将二次型解：解：该二次

10、型所对应的矩阵为化成标准形现在学习的是第30页，共58页 A=4 2 2;-2 1 1/2;2 1/2 1;P,T =schur(A)P=0.5458 -0.0000 0.8379 0.5925 0.7071 -0.3859 -0.5925 0.7071 0.3859现在学习的是第31页，共58页T=-0.3423 0 0 0 0.5000 0 0 0 5.8423答：所作的正交变换为：二次型的标准型为：现在学习的是第32页，共58页三、三、正定二次型的判定正定二次型的判定现在学习的是第33页，共58页现在学习的是第34页，共58页1.1.顺序主子式判断法顺序主子式判断法 求二次型求二次型 F

11、=XAX 的矩阵的矩阵 A 的各阶顺序的各阶顺序 主子式主子式 Di (i=1,2,3.);判断判断 Di 是否大于是否大于0.程序：建立函数文件程序：建立函数文件 shxu.m现在学习的是第35页，共58页function C,M=shxu(A)%C为为A的各阶顺序主子式组成的向量的各阶顺序主子式组成的向量%M为判定向量为判定向量:if C(i)0,then M(i)=1;%others M(i)=0 n=size(A);C=;M=;for i=1:n(1)A1=A(1:i,1:i);D=det(A1);C=C D;if D0 m=1;else m=0;end M=M,m;end 现在学习的

12、是第36页，共58页2 2、特征值判别法、特征值判别法 求二次型求二次型 f=XAX 的矩阵的矩阵 A 的全部特征的全部特征 值值 （i=1,2,i=1,2,）；）；判断判断 是否大于是否大于 0.程序：建立函数文件程序：建立函数文件 tezh.m现在学习的是第37页，共58页 function T,M =tezh(A)n=size(A);T=(eig(A);M=;for i=1:n(1)if T(i)0 m=1;else m=0;end M=M,m;end现在学习的是第38页，共58页例例8 8 判定下列二次型是否正定 解解 二次型矩阵现在学习的是第39页，共58页方法一方法一 顺序主子式顺

13、序主子式 A=1 1 2 1;-1 3 0 3;2 0 9 6;1 3 6 19;C,M=shxu(A)答：此二次型是正定的。C=1 2 6 24 M=1 1 1 1现在学习的是第40页，共58页方法二方法二 特征值法特征值法 T=0.0643 2.2421 7.4945 22.1991 M=1 1 1 1 A=1 1 2 1;-1 3 0 3;2 0 9 6;1 3 6 19 T,M =tezh(A)答：此二次型是正定的。现在学习的是第41页，共58页例例9 9 判定下列二次型是否正定 解解 二次型矩阵现在学习的是第42页，共58页方法一方法一 顺序主子式顺序主子式 A=9 6 24;-6

14、130 30;24 30 71;C,M=shxu(A)答：此二次型是正定的。C=9 1134 6174 M=1 1 1 现在学习的是第43页，共58页方法二方法二 特征值法特征值法 T=0.6576 65.0894 144.2530M=1 1 1 A=9 6 24;-6 130 30;24 30 71;T,M =tezh(A);答：此二次型是正定的。现在学习的是第44页，共58页例例1010 判定下列二次型是否正定 解解 二次型矩阵现在学习的是第45页，共58页方法一方法一 顺序主子式顺序主子式 A=10 4 12;4 2 14;12 14 1;C,M=shxu(A)答：此二次型不是正定的。C

15、=10 4 -3588M=1 1 0现在学习的是第46页，共58页方法二方法二 特征值法特征值法T=-17.4209 10.1708 20.2501M=0 1 1 A=10 4 12;4 2 14;12 14 1;T,M =tezh(A)答：此二次型不是正定的。现在学习的是第47页，共58页现在学习的是第48页，共58页function C,M=shxuf(A)%C为为A的各阶顺序主子式组成的向量的各阶顺序主子式组成的向量%M为判定向量为判定向量:if C(i)0,then M(i)=1;if C(i)0 m=1;elseif D0 m=-1；else m=0;end M=M,m;end 现在

16、学习的是第49页，共58页 function T,M =tezhf(A)n=size(A);T=(eig(A);M=;for i=1:n(1)if T(i)0 m=1;elseif T(i)A=-1 1 2 1;1 3 0 3;-2 0 9 6;-1 3 6-19;C,M=shxuf(A)答：此二次型是负定的。C=-1 2 -6 24M=-1 1 -1 1现在学习的是第53页，共58页方法二方法二 特征值法特征值法T=-22.1991 -7.4945 -2.2421 -0.0643M=-1 -1 -1 -1 A=-1 1 2 1;1 3 0 3;-2 0 9 6;-1 3 6-19;T,M =

17、tezhf(A)答：此二次型是负定的。现在学习的是第54页，共58页1、已知矩阵(1)求矩阵A的特征值;(2)求矩阵A的特征值对应的全部特征向量.习题习题现在学习的是第55页，共58页2 2 判断下列方阵是否可对角化,若可对角 化，求出可逆阵P，使P-1AP为对角阵。现在学习的是第56页，共58页3、已知二次型 f=x1x2+x2x3+x3x4+x4x1 (1)写出二次型矩阵A;(2)用正交变换将二次型化为标准形,并 写出所作的正交变换;现在学习的是第57页，共58页5、判别下列二次型是否为负定负定二次型 （用两种方法求，写出程序）4、判别下列二次型是否为正定正定二次型 （用两种方法求，写出程序）现在学习的是第58页，共58页