MATLAB 之特征值与特征向量、jordan标准形

实验六  特征值与特征向量、若当标准形

【实验目的】

1．了解特征值与特征向量基本概念及其性质；

2．了解若当标准型的基本概念；

3．学习、掌握MATLAB软件有关的命令。

【实验准备】

1．特征多项式

设A为n阶方阵, 如果数“ ”和n维列向量x使得关系式 成立, 则称 为方阵A的特征值, 非零向量x称为A对应于特征值“ ”的特征向量。

poly(A)，返回矩阵A的特征多项式的向量表示形式，例如：

>> clear

>> A=[1 0;2 3];

>> p=poly(A)         %矩阵A的特征多项式的向量表示形式

p =

     1    -4     3

>> f=poly2str(p,’x’)     %矩阵A的特征多项式

f =

   x^2 – 4 x + 3

或者由定义出发，计算特征多项式.例如：

>> clear

>> A=[1 0;2 3];

>> E=eye(2);        %2阶单位阵

>> syms x

>> f=det(x*E-A)      %矩阵A的特征多项式

f =

 (x-1)*(x-3)

2．特征值与特征向量eigenvalue

求一个方阵的特征值与特征向量可以使用函数eig( ).

d=eig(A),  返回A所有特征值组成的列向量d.

[V,D]= eig(A),  返回A所有特征值组成的矩阵D和特征向量组成的矩阵V.

[V,D]= eigs(A), 返回A所有特征值(按大小次序)组成的对角矩阵D和特征向量组成的矩阵V，且满足D=V-1AV.

d=eig(A,B), 返回复数矩阵A+Bi所有特征值组成的向量d.

[V,D]= eig(A,B), 返回复数矩阵A+Bi所有特征值组成的矩阵D和特征向量组成的矩阵V.

例如：

>> clear

（>> format）（>> format rat）

>> A=[0 1 0 0;1 0 0 0;0 0 0 1;0 0 1 0];

>> d=eig(A)     %求矩阵A的特征值

d =

     1

    -1

     1

    -1

>> %特征值以列向量的形式输出，例如：

>> [V,D]=eig(A)      %求矩阵A的特征值与特征向量所组成的矩阵

V =

   -0.7071         0         0    0.7071

    0.7071         0         0    0.7071

         0   -0.7071    0.7071         0

         0    0.7071    0.7071         0

D =

    -1     0     0     0

     0    -1     0     0

     0     0     1     0

     0     0     0     1

>>%说明（1）矩阵D的主对角线上的元素为特征值，