用matlab写一段求矩阵特征值和特征向量的代码

举一个简单的例子供参考一下
>>A=sym('[1 1/5 1/3 1/3;5 1 5 5;3 1/5 1 1/3;3 1/5 3 1 ]')

A =
[ 1, 1/5, 1/3, 1/3]
[ 5, 1, 5, 5]
[ 3, 1/5, 1, 1/3]
[ 3, 1/5, 3, 1]

>> [V,D]=eig(A)

V =

[ 1/3*i*3^(1/2), -1/3*i*3^(1/2), 1/3*3^(1/2), -1/3*3^(1/2)]
[ -5, -5, 5, 5]
[ 1, 1, 1, 1]
[ -i*3^(1/2), i*3^(1/2), 3^(1/2), -3^(1/2)]

D =

[ 2/3*i*3^(1/2), 0, 0, 0]
[ 0, -2/3*i*3^(1/2), 0, 0]
[ 0, 0, 2+4/3*3^(1/2), 0]
[ 0, 0, 0, 2-4/3*3^(1/2)]

其它的道理一样。
[V,D] = EIG(X) produces a diagonal matrix D of eigenvalues and a full matrix V whose columns are the corresponding eigenvectors so that X*V = V*D.

V是特征向量，D是特征值

mark先