|
马上注册,结识高手,享用更多资源,轻松玩转三维网社区。
您需要 登录 才可以下载或查看,没有帐号?注册
x
一个matlab问题' |4 g A& f$ l. K; ~) X
我得到一个矩阵,但是只能求出他的本征值,无法得到他的本征向量
8 B3 x) ^3 ^8 t9 d# c) t知道为什么吗?1 b! \7 ?* [9 o, y$ P3 ^, p
我是用matlab算的
1 u m, _: ^+ u, q! |# n0 ysyms j b q p u;. f) N u" s4 e5 c
x=[0 1;1 0];2 S3 ^6 Y& n2 m. b
y=[0 -i;i 0];
; q* s3 v; a p, G8 _: sz=[1 0;0 -1];1 P/ Y }5 d7 M8 w' h" [
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])
( m+ v d5 z' E; l& ^. a7 z可以得到h矩阵,
- }4 V* Q' e( x$ X9 O: U用eig(h) 很快得到了他的本征值1 C. R. T3 O4 k# e y, n' a0 P# B5 o
但是用[v d]=eig(h)时就会出错,
+ {9 m/ O, Z2 D! a% l# o但是我也想知道他的本征向量v
6 K0 ~: o- m! @3 P' Q! ~8 _错误原因是:: a* I4 X9 R0 z; l$ x1 p
?? Error using ==> sym.eig
% @) K' w% g q" w7 K% WError, (in eigenvectors) eigenvects only works for a matrix of rationals, rational functions, algebraic numbers, or algebraic functions at present
$ _& D9 S) Z. p% R2 B/ X9 D$ d" j请问为什么?如果能直接帮我解出来将不甚感激! |
|