一个matlab问题

pkuwuli

一个matlab问题
) m7 W0 g: V1 `- a  _1 }( y3 M6 V2 y. H我得到一个矩阵,但是只能求出他的本征值,无法得到他的本征向量
4 O  x: m% E% b" [1 e知道为什么吗？
我是用matlab算的
' C! J  G& K3 {) csyms j b q p u;
x=[0 1;1 0];
/ y* U9 d2 R! O7 b$ N% ?# Ly=[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矩阵,
$ Z; J: k% }$ |* z& {" o3 }用eig(h) 很快得到了他的本征值
$ ^. ^1 f% `) x1 c2 W7 i! q& a3 r# h但是用[v d]=eig(h)时就会出错，
- T3 M( A2 P, @+ }但是我也想知道他的本征向量v
错误原因是:
5 O, l5 |+ j, ~9 x% E？? Error using ==> sym.eig
Error, (in eigenvectors) eigenvects only works for a matrix of rationals, rational functions, algebraic numbers, or algebraic functions at present
6 i) i' d# i/ c6 Z) U请问为什么？如果能直接帮我解出来将不甚感激！