一个matlab问题

pkuwuli

一个matlab问题
我得到一个矩阵,但是只能求出他的本征值,无法得到他的本征向量
: S# P. P% A/ j. l知道为什么吗？
我是用matlab算的
syms j b q p u;
5 X3 j0 T: ?+ O3 l' }. }) T! nx=[0 1;1 0];
9 B1 A" B' r$ n% K% }+ v" ^; Ry=[0 -i;i 0];
& e. F  \) T% y# f  Yz=[1 0;0 -1];
+ H; `. b. o8 ]* H, S2 {) \' Fh=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])
9 ?3 F, |' u) i0 V- D& Q可以得到h矩阵,
0 z2 z! H7 `5 ?8 Z0 l: X( i5 j用eig(h) 很快得到了他的本征值
3 k1 z) Z6 a. J: {但是用[v d]=eig(h)时就会出错，
" a5 K% z; S" Z: j' f但是我也想知道他的本征向量v
" T* o% ^0 C4 \6 s; G7 D错误原因是:
？? Error using ==> sym.eig
$ ^6 Q; Z& F% _9 c' U: n( F6 s0 [Error, (in eigenvectors) eigenvects only works for a matrix of rationals, rational functions, algebraic numbers, or algebraic functions at present
- R% v* C4 E- O0 n0 ^4 s请问为什么？如果能直接帮我解出来将不甚感激！