一个matlab问题

pkuwuli

一个matlab问题
( V# J7 m: E3 m/ }8 M9 O% S我得到一个矩阵,但是只能求出他的本征值,无法得到他的本征向量
知道为什么吗？
; D( [0 I0 x, a) H* D, `我是用matlab算的
5 g) l/ ~; e" F, X9 Gsyms j b q p u;
x=[0 1;1 0];
, J0 A+ z3 v; c& Sy=[0 -i;i 0];
z=[1 0;0 -1];
h=4*j*(kron(x,x)+kron(y,y)+kron(z,z))+u*b*kron(((sin(q)*cos(p))*[0 1;1 0]+(sin(q)*sin(p))*[0 -i;i 0]+cos(q)*[1 0;0 -1]),[1 0;0 1])
可以得到h矩阵,
用eig(h) 很快得到了他的本征值
但是用[v d]=eig(h)时就会出错，
  U# i9 J8 b4 a' m! a: \" E" t但是我也想知道他的本征向量v
; d5 J- L6 S$ }7 [/ x1 y错误原因是:
？? Error using ==> sym.eig
/ ]+ p- K* ]; b) }Error, (in eigenvectors) eigenvects only works for a matrix of rationals, rational functions, algebraic numbers, or algebraic functions at present
请问为什么？如果能直接帮我解出来将不甚感激！