特征值和特征向量、若当标准型.docVIP

特征值与特征向量、若当标准型[实验目的] 1.?????? 掌握特征值与特征向量基本概念及其性质 2.?????? 掌握若当标准型的基本概念 3.?????? 学习、掌握MATLAB软件有关的命令。 [实验准备] 1．求一个方阵的特征值与特征向量可以使用函数eig( )。 d=eig(A), 返回A所有特征值组成的向量d. [V,D]= eig(A), 返回A所有特征值组成的矩阵D和特征向量组成的矩阵V.() d=eig(A,B), 返回复数矩阵A+Bi所有特征值组成的向量d. [V,D]= eig(A,B), 返回复数矩阵A+Bi所有特征值组成的矩阵D和特征向量组成的矩阵V. 2．trace( ) 求矩阵的迹 3．sum(A,dim) 矩阵元素求和函数，dim=1则按列求和，dim=2则按行求和。 prod(A,dim) 矩阵元素求积函数，dim=1则按列求积，dim=2则按行求积。 4．jordan( ) 求若当标准型 [实验内容] 一、特征值与特征向量 例1 求矩阵A=的特征值与特征向量。? 相应的MATLAB代码及计算结果如下: Clear A= [1 2 2;2 1 2; 2 2 1]; [V,D]=eig(A) V = 0.6015 0.5522 0.5774 0.1775 -0.7970 0.5774 -0.7789 0.2448 0.5774 D = -1.0000 0 0 0 -1.0000 0 0 0 5.0000 d=eig(A) d = -1.0000 -1.0000 5.0000 inv(V)*A*V ans = -1.0000 0.0000 -0.0000 -0.0000 -1.0000 -0.0000 0.0000 0.0000 5.0000 说明：（1）方阵A的特征值为， （2）特征值-1的特征向量为V中的第1、2列，特征值5的特征向量为V中的第3列。每一个特征向量的范数为1，且彼此之间正交。 二、矩阵的迹 矩阵的迹可用函数trace( )来求。 例2 求例1中矩阵A的迹，并验证 相应的MATLAB代码及计算结果如下: clear A= [1 2 2;2 1 2; 2 2 1]; a=trace(A) a = 3 d=eig(A) d = -1.0000 -1.0000 5.0000 b=sum(d,1) b = 3 e=det(A) e = 5 f=prod(d,1) f = 5 说明：（1）函数sum(A,dim)是矩阵元素求和函数，dim=1则按列求和，dim=2则按行求和。 （2）函数prod(A,dim)是矩阵元素求积函数，dim=1则按列求积，dim=2则按行求积。 三、若当标准形 若当标准型可用函数jordan( ) 来求。 J = jordan(A), 其中J为A的若当标准型. 例3 求矩阵A=的若当标准型。 ? 相应的MATLAB代码及计算结果如下: clear A=[-1 -2 6;-1 0 3;-1 -1 4] A = -1 -2 6 -1 0 3 -1 -1 4 J=jordan(A) J = 1 1 0 0 1 0 0 0 1 说明：Matlab中若当块是按上三角形定义的。 作业: P324 19#