一个matlab问题

pkuwuli

一个matlab问题
% p$ g  d4 i  `1 n; r我得到一个矩阵,但是只能求出他的本征值,无法得到他的本征向量
知道为什么吗？
* E' Z1 K4 g; [$ ^! Q我是用matlab算的
& |3 {! t( J. C8 Fsyms j b q p u;
2 s9 U9 n- n- t# D7 `% |x=[0 1;1 0];
y=[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])
% \2 l' [; X3 G* i8 X可以得到h矩阵,
! F% b( O. r2 O用eig(h) 很快得到了他的本征值
7 o0 m( ~) }3 H8 p! |; c2 u, h但是用[v d]=eig(h)时就会出错，
" u3 ^2 O$ E) _; P8 u6 C0 E7 I) g' g但是我也想知道他的本征向量v
错误原因是:
？? Error using ==> sym.eig
0 b% U/ [$ \( j8 E. o# e0 `3 aError, (in eigenvectors) eigenvects only works for a matrix of rationals, rational functions, algebraic numbers, or algebraic functions at present
' K' @6 i* t# |1 G9 }3 X请问为什么？如果能直接帮我解出来将不甚感激！