|
|
马上注册，结识高手，享用更多资源，轻松玩转三维网社区。
您需要 登录 才可以下载或查看，没有帐号？注册
x
一个matlab问题
9 Z0 w4 W8 s- i5 b4 n9 U V我得到一个矩阵,但是只能求出他的本征值,无法得到他的本征向量: P3 O a% p) \/ Y0 k" X' H
知道为什么吗？
# }; }3 t; R" h a6 e- I- @0 c我是用matlab算的
, p! J' }& O0 esyms j b q p u;2 C5 a, y1 u+ _6 T
x=[0 1;1 0];
$ \8 D* z2 G& A5 o' ty=[0 -i;i 0];6 T% M9 d! T% L
z=[1 0;0 -1];
9 R5 \3 A$ b# [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]). \7 Y. ^# ]5 X! G5 L% y/ ]5 T- y+ p
可以得到h矩阵,
% I. T& U+ P0 t; L# P$ I0 s用eig(h) 很快得到了他的本征值% C3 R/ r- I3 U6 M
但是用[v d]=eig(h)时就会出错，, T F% A8 W1 @! j2 m' l; }8 T$ W5 c
但是我也想知道他的本征向量v& S) [$ K5 M) q8 T. y
错误原因是:
7 Q; x9 t7 d0 L- { A. k？? Error using ==> sym.eig% ^! u5 v4 L9 d$ `0 L7 c
Error, (in eigenvectors) eigenvects only works for a matrix of rationals, rational functions, algebraic numbers, or algebraic functions at present. `/ s: \1 V* i( ] t% N0 k' N) U3 n7 z
请问为什么？如果能直接帮我解出来将不甚感激！ |
|
[复制链接]
发表于 2007-4-30 01:07:27