|
马上注册,结识高手,享用更多资源,轻松玩转三维网社区。
您需要 登录 才可以下载或查看,没有帐号?注册
x
一个matlab问题
/ r1 t& r" {5 \' R我得到一个矩阵,但是只能求出他的本征值,无法得到他的本征向量
$ t; z4 |' N2 ?' ?( T知道为什么吗?
7 I! [1 ]6 o& }: z我是用matlab算的
" M: K+ }9 {: C) Rsyms j b q p u;
( e9 }( h" R' bx=[0 1;1 0];
9 U) L" K, w* B& \+ N t' xy=[0 -i;i 0];8 i! p9 G, D. b5 k7 y$ B
z=[1 0;0 -1];
% u+ I2 z1 t" _0 ]' I; ^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])
% J( ?" {1 N7 }2 C& U( ~5 |可以得到h矩阵,1 S2 t" l3 [7 p! P) Y& [, J$ X, q
用eig(h) 很快得到了他的本征值
2 [ T! Z5 Z [但是用[v d]=eig(h)时就会出错,9 N. _8 q% }. Q: u2 Z; E
但是我也想知道他的本征向量v/ u: ^+ o' u2 \" ~
错误原因是:
1 g, N$ {- i- a: z& x; u4 z) S7 Z4 H?? Error using ==> sym.eig- N& @6 f2 n' K, u
Error, (in eigenvectors) eigenvects only works for a matrix of rationals, rational functions, algebraic numbers, or algebraic functions at present
' T# z$ S+ V0 L3 r+ x请问为什么?如果能直接帮我解出来将不甚感激! |
|