|
马上注册,结识高手,享用更多资源,轻松玩转三维网社区。
您需要 登录 才可以下载或查看,没有帐号?注册
x
一个matlab问题
* q/ T8 C; ~' w1 V' |/ A: @我得到一个矩阵,但是只能求出他的本征值,无法得到他的本征向量4 C, C: h; y' ]8 ^
知道为什么吗?- {" P6 X& N* o
我是用matlab算的: |0 G( ] F; l, C
syms j b q p u;
+ M; @/ x1 x! px=[0 1;1 0];
9 K) h- T7 d/ g2 F8 Hy=[0 -i;i 0];% v9 \: m( {% i" r9 x% Q) O
z=[1 0;0 -1];
4 X0 D8 t, M. t, W$ x; Bh=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 R" `6 Z: V5 G, s- b- e可以得到h矩阵,4 [$ G2 |. k- J: [& F8 ~/ a! ?
用eig(h) 很快得到了他的本征值
0 }( `) {' K$ H) \' m但是用[v d]=eig(h)时就会出错,
g0 x: H$ p8 B! j但是我也想知道他的本征向量v. s8 j7 C! B1 L0 D% X4 k
错误原因是:
5 ~' o7 C: `5 [3 Y?? Error using ==> sym.eig
) Q3 p3 P B: [% @- nError, (in eigenvectors) eigenvects only works for a matrix of rationals, rational functions, algebraic numbers, or algebraic functions at present/ O- I# @& r4 K' J
请问为什么?如果能直接帮我解出来将不甚感激! |
|