在matlab中，可以用eig函数计算矩阵的特征值和特征向量。举例如下：
>> a = [1, 2, 3; 4, 5, 6; 7, 8, 9] % 原始数据矩阵
a =
1 2 3
4 5 6
7 8 9
>> [V, D] = eig(a) % 特征值分解，其中V的每一列表示矩阵a的一个特征向量，D是一个对角矩阵，对角线上的元素表示矩阵a的特征值
V =
-0.2320 -0.7858 0.4082
-0.5253 -0.0868 -0.8165
-0.8187 0.6123 0.4082
D =
16.1168 0 0
0 -1.1168 0
0 0 -0.0000
