|
|
马上注册,结识高手,享用更多资源,轻松玩转三维网社区。
您需要 登录 才可以下载或查看,没有帐号?注册
x
一个matlab问题1 S8 U# A& H# Q" R8 A
我得到一个矩阵,但是只能求出他的本征值,无法得到他的本征向量
; _8 t3 X& }/ m3 ^知道为什么吗?8 u0 K. B9 @* z; Y- Q& b, Z
我是用matlab算的, Q: M' B- P5 C4 O3 G2 k( z6 _4 [
syms j b q p u;
3 @ w' L) z! q# N. px=[0 1;1 0];2 w4 p9 s: V( p( J! c- Y
y=[0 -i;i 0];8 z- D+ x6 h) M" S8 M
z=[1 0;0 -1];
5 Z# y: A3 g+ r) l) v1 Ch=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]); h, c, H- F7 R# d3 d
可以得到h矩阵,+ o$ m8 T$ S9 a% k: _" p/ H" ?% O
用eig(h) 很快得到了他的本征值
1 t7 P$ K, W- X% a1 P( l但是用[v d]=eig(h)时就会出错,
0 t8 S+ f/ s3 f0 p7 c但是我也想知道他的本征向量v
6 L- b* ?& _4 O+ g6 \. q错误原因是:4 [2 f6 Z S3 U( ?- I
?? Error using ==> sym.eig
" P# T {/ r' k- D0 vError, (in eigenvectors) eigenvects only works for a matrix of rationals, rational functions, algebraic numbers, or algebraic functions at present: J& ]( ]% ^* {8 a k* U
请问为什么?如果能直接帮我解出来将不甚感激! |
|