实验四矩阵的特征值与特征向量.docVIP

实验四 矩阵的特征值与特征向量 1、方阵的特征多项式与迹 在MATLAB里，命令poly用于求方阵的特征多项式，它既有数值功能，又有符号功能；trace用于求方阵的迹，它只有数值功能。 2、方阵的特征值与特征向量 MATLAB中roots用于求多项式的根，eig(A)可同时求矩阵A的特征值与特征向量。 例1 求矩阵 的特征多项式，迹，特征值与特征向量。 A=[-1,1,0;-4,3,0;1,0,2]; Q=poly(sym(A)) Q = x^3-4*x^2+5*x-2 T=trace(A) T = 4 p=poly(A) p = 1 -4 5 -2 r=roots(p) r = 2.0000 1.0000 + 0.0000i 1.0000 - 0.0000i [V,D]=eig(A,nobalance) V = 0.5000 -0.5000 -0.0000 1.0000 -1.0000 -0.0000 -0.5000 0.5000 1.0000 D = 1.0000 0 0 0 1.0000 0 0 0 2.0000 3、相似变换 设是可逆阶方阵，是任意阶方阵，则变换称为相似变换。在MATLAB里，写成。 例2 已知4阶魔方矩阵和可逆矩阵,,分别求方阵的特征多项式。 A=magic(4) A = 16 2 3 13 5 11 10 8 9 7 6 12 4 14 15 1 P=[2,0,0,1;3,1,5,3;1,2,9,9;3,3,0,1]; P=sym(P)；A=sym(A) ； B=P\A*P B = [ 5172/175, 3198/175, 678/25, 228/7] [ -2424/175, -1066/175, 1174/25, 204/7] [ -2117/175, -878/175, 342/25, 26/7] [ 3656/175, 1829/175, -431/25, -22/7] PA=poly(A) PA = x^4-34*x^3-80*x^2+2720*x PB=poly(B) PB = x^4-34*x^3-80*x^2+2720*x 实际上此例验证了相似的矩阵具有相同的特征多项式。 4、方阵的相似对角化 在MATLAB里，可用命令eig的符号功能来判断方阵是否可以相似对角化。 例3 已知方阵,请问该矩阵是否可以相似对角化。 A=[1,2,0,0,0;0,1,0,0,0;0,0,3,1,0;0,0,0,3,0;0,0,0,0,2]; n=size(A) n = 5 5 A=sym(A)； [V,D,P]=eig(A) V = [ 1, 0, 0] [ 0, 0, 0] [ 0, 1, 0] [ 0, 0, 0] [ 0, 0, 1] D = [ 1, 0, 0, 0, 0] [ 0, 1, 0, 0, 0] [ 0, 0, 3, 0, 0] [ 0, 0, 0, 3, 0] [ 0, 0, 0, 0, 2] P = 1 3 5 5、 实对称矩阵的正交相似对角化 用命令eig就可以将实对称矩阵正交相似对角化。 例4 已知方阵，求正交矩阵,使为对角矩阵。 A=[1,1,1,0,0;1,1,0,0,0;1,0,3,1,0;0,0,1,3,2;0,0,0,2,2]; [Q,D]=eig(A,nobalance) Q = -0.7116 -0.1834 0.5184 -0.4210 0.1186 0.5839 0.2763 0.7390 -0.1888 0.0300 0.2832 -0.1546 -0.3754 -0.7501 0.4386 -0.2001 0.5952 -0.0310 0.2483 0.7370 0.1803 -0.7154 0.2078 0.4036 0.4995 D = -0.2186 0 0 0 0 0 0.3362