|
|
马上注册，结识高手，享用更多资源，轻松玩转三维网社区。
您需要 登录 才可以下载或查看，没有帐号？注册
x
一个matlab问题
0 [/ p9 a2 ?+ v& d( ]$ ~# z我得到一个矩阵,但是只能求出他的本征值,无法得到他的本征向量3 q5 n9 g p. ~+ F5 ]( e
知道为什么吗？3 @3 m. z; R% L2 v9 `# \, j8 i( {6 ~
我是用matlab算的
# P6 j9 [' K$ m/ Q( a" V3 Jsyms j b q p u;
- y, R G; e% fx=[0 1;1 0];% R% s: _: m/ U U$ h0 N( h# D
y=[0 -i;i 0];- B$ T$ F# u# b
z=[1 0;0 -1];
2 w! T/ O: p* M) A- X9 D7 nh=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 h4 ]/ M4 M% l8 R# a7 k可以得到h矩阵,
1 r! U) j$ B& D. p$ L用eig(h) 很快得到了他的本征值
* U1 S- t6 k. B; h, u+ f9 s但是用[v d]=eig(h)时就会出错，) [0 E3 k- i- n% e
但是我也想知道他的本征向量v" C$ v9 p! ~6 ]7 q- R# D8 _% z' M% q; o7 `
错误原因是:
1 l6 p& n" K- N5 D1 d2 E7 m？? Error using ==> sym.eig
+ G t( w% r% \) v8 F5 r" n: B) fError, (in eigenvectors) eigenvects only works for a matrix of rationals, rational functions, algebraic numbers, or algebraic functions at present
/ Y, P9 z( B; k. B请问为什么？如果能直接帮我解出来将不甚感激！ |
|
[复制链接]
发表于 2007-4-30 01:07:27