|
|
马上注册，结识高手，享用更多资源，轻松玩转三维网社区。
您需要 登录 才可以下载或查看，没有帐号？注册
x
一个matlab问题
. y/ j: F% K. J1 \& r我得到一个矩阵,但是只能求出他的本征值,无法得到他的本征向量2 J8 d+ _/ I e9 v
知道为什么吗？& K6 S+ H0 m8 }, U" F9 t- C* K3 x
我是用matlab算的
* w h6 B C0 k$ Qsyms j b q p u;5 Z. c9 a2 C( |" c3 z* Q
x=[0 1;1 0];7 e D2 I, Z$ {
y=[0 -i;i 0];: Y1 e# k" T; M! g8 ?6 N) j% ^7 n( Q
z=[1 0;0 -1];
) B( }+ ~* b+ b: m3 jh=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])4 T# c' p8 _" S( D$ Y: _
可以得到h矩阵,/ R8 \& I y/ ~0 q
用eig(h) 很快得到了他的本征值! O0 s+ g! @- g4 c
但是用[v d]=eig(h)时就会出错，. b! J# }5 x) e: p
但是我也想知道他的本征向量v7 w" f: u; X3 v4 n- p7 b9 L; w# H5 Q/ s
错误原因是:
8 p4 f" P0 d) ~* j. a; U4 b( Z3 i？? Error using ==> sym.eig
# e7 s! S! j, S' V4 d( qError, (in eigenvectors) eigenvects only works for a matrix of rationals, rational functions, algebraic numbers, or algebraic functions at present
3 o( w3 w, B. _/ b; n请问为什么？如果能直接帮我解出来将不甚感激！ |
|
[复制链接]
发表于 2007-4-30 01:07:27