一个matlab问题

pkuwuli

一个matlab问题
4 Q  Z" F6 W9 M+ q我得到一个矩阵,但是只能求出他的本征值,无法得到他的本征向量
0 ~$ K% `7 A: x' k9 h知道为什么吗？
0 _, p/ p% c  c8 \8 ?& w我是用matlab算的
syms j b q p u;
: r0 @5 R- _3 K  I8 f1 px=[0 1;1 0];
4 j* p+ G* O, `+ g5 [# y% ry=[0 -i;i 0];
2 x2 O; K9 t  R2 P& c5 o& f% kz=[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) 很快得到了他的本征值
9 H, p# W. v5 t( |但是用[v d]=eig(h)时就会出错，
. w0 S$ j; s( j7 I  h0 I$ Q1 y# m但是我也想知道他的本征向量v
错误原因是:
) B$ F0 ~0 F/ O/ p: m2 X？? Error using ==> sym.eig
Error, (in eigenvectors) eigenvects only works for a matrix of rationals, rational functions, algebraic numbers, or algebraic functions at present
请问为什么？如果能直接帮我解出来将不甚感激！