|
|
马上注册，结识高手，享用更多资源，轻松玩转三维网社区。
您需要 登录 才可以下载或查看，没有帐号？注册
x
一个matlab问题 ~7 G4 v- y( n5 Z4 `4 t) l3 w1 G
我得到一个矩阵,但是只能求出他的本征值,无法得到他的本征向量# s" e ]. L; t2 T* R- B S
知道为什么吗？8 e2 x5 D; h, t2 P/ q
我是用matlab算的
0 ]- R% ]( S, e0 J2 Vsyms j b q p u;! ?4 W I" w2 t- }# }
x=[0 1;1 0];
% t+ J% C9 X- J5 D4 `1 G4 I+ {y=[0 -i;i 0];5 }# E) T q$ v& y
z=[1 0;0 -1];% D1 Y( O& ^8 r; z: i( C
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])
& U* n8 @/ h; @& G( @4 R! a: k可以得到h矩阵,, P- Q; H2 t' ~6 L& e: ~$ w+ ~( } m
用eig(h) 很快得到了他的本征值
/ }0 D( [% {6 O9 b但是用[v d]=eig(h)时就会出错，
/ I7 |# G2 g. v但是我也想知道他的本征向量v
8 V, X! j* _9 B错误原因是:
S6 P" q2 D9 t4 d6 f" x4 K% w0 H& o？? Error using ==> sym.eig
- q! M% U0 ^/ B S2 H7 AError, (in eigenvectors) eigenvects only works for a matrix of rationals, rational functions, algebraic numbers, or algebraic functions at present0 @ l. K0 s' _" P+ m( M# w
请问为什么？如果能直接帮我解出来将不甚感激！ |
|
[复制链接]
发表于 2007-4-30 01:07:27