|
|
马上注册,结识高手,享用更多资源,轻松玩转三维网社区。
您需要 登录 才可以下载或查看,没有帐号?注册
x
一个matlab问题
! s/ C$ u; v( z/ m3 [! A1 k' F我得到一个矩阵,但是只能求出他的本征值,无法得到他的本征向量+ n% @! d+ \ w$ D/ G; a! m, m
知道为什么吗?
6 @6 l2 }6 ^- V* T0 H我是用matlab算的
+ Z6 E' ~0 k; N! x$ Gsyms j b q p u;/ d/ p# w# K8 s3 G
x=[0 1;1 0];9 p9 n$ t* {4 K- F& ^
y=[0 -i;i 0];
9 _7 X; X- o8 G4 dz=[1 0;0 -1];
7 F+ f0 p- R& u- }( ]& C: v1 N. \. \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])& s* F" q/ `+ ~/ z# l! F4 `
可以得到h矩阵,
9 e* f: ^2 `. e, l- l9 K用eig(h) 很快得到了他的本征值
3 l/ N* K) k1 x- e但是用[v d]=eig(h)时就会出错,
) Q3 S3 ]& i' z% O8 d- \6 x但是我也想知道他的本征向量v
, H: u# y' T0 k- v4 ]$ T% }错误原因是:
4 p5 v3 x# r; ?7 W6 R' q F?? Error using ==> sym.eig( k/ h; g# v+ x' i: S& R; J& Q
Error, (in eigenvectors) eigenvects only works for a matrix of rationals, rational functions, algebraic numbers, or algebraic functions at present
6 K! t% {0 H/ u! Q3 K2 P请问为什么?如果能直接帮我解出来将不甚感激! |
|