|
|
马上注册，结识高手，享用更多资源，轻松玩转三维网社区。
您需要 登录 才可以下载或查看，没有帐号？注册
x
一个matlab问题1 ~, g7 B( x9 J' l% Y/ U3 Z' m
我得到一个矩阵,但是只能求出他的本征值,无法得到他的本征向量+ }' y* l7 `+ w# R$ R
知道为什么吗？# s# K' `; O" Y" q1 }" j. c
我是用matlab算的% g/ l7 l6 i1 B8 }' ^' K
syms j b q p u;) k$ s% F6 r% Y( E& v, g$ T
x=[0 1;1 0];& |6 x9 t4 L# E H8 Q, l
y=[0 -i;i 0];. I* Z( `1 P: c g H$ [1 ~& L: ^' O+ n
z=[1 0;0 -1];
! t0 m4 g/ N4 ^; S' S0 o t. i9 v# _! kh=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])+ b7 C# G/ Z. F) J8 q
可以得到h矩阵,
& p% y e! G3 O- [# Q& `用eig(h) 很快得到了他的本征值2 U5 e& J& n& V, K
但是用[v d]=eig(h)时就会出错，7 T: c3 P8 ?3 r& R: F* B1 h
但是我也想知道他的本征向量v) n/ o3 I0 V, T
错误原因是:
. j V# R: y( P) G' Z- i4 F？? Error using ==> sym.eig0 N1 ?8 c* ]/ ?: f. A" T
Error, (in eigenvectors) eigenvects only works for a matrix of rationals, rational functions, algebraic numbers, or algebraic functions at present# y$ Y# a% G0 o- j4 ]& _
请问为什么？如果能直接帮我解出来将不甚感激！ |
|
[复制链接]
发表于 2007-4-30 01:07:27