|
|
马上注册,结识高手,享用更多资源,轻松玩转三维网社区。
您需要 登录 才可以下载或查看,没有帐号?注册
x
一个matlab问题
$ H2 M, P! L% y; y我得到一个矩阵,但是只能求出他的本征值,无法得到他的本征向量( R# D: o0 k. p1 v% _
知道为什么吗?
6 f4 v' e0 J* {$ ~我是用matlab算的
" n0 Q9 U0 Y3 c% o. E1 j! ksyms j b q p u;/ V9 V3 k0 b- T. ?7 `) e
x=[0 1;1 0]; B) M2 Z* \" x! H
y=[0 -i;i 0];4 e3 C0 n+ ?4 u- }: n5 m( N" J
z=[1 0;0 -1];5 N U2 v* w! C8 {
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])
9 h* K3 ?: ?# ^7 Z可以得到h矩阵,
, ^) ^% Y3 u0 @/ c) {; r8 \( n用eig(h) 很快得到了他的本征值' [' I' E" c* s7 N
但是用[v d]=eig(h)时就会出错,
7 y; g. j- H% [ `但是我也想知道他的本征向量v( |& A' L& L0 g: ~! y! p( {, g; \
错误原因是:3 v$ l5 ?4 y8 e: N, [% R! d) X8 O
?? Error using ==> sym.eig
$ w# W" {( S$ y# YError, (in eigenvectors) eigenvects only works for a matrix of rationals, rational functions, algebraic numbers, or algebraic functions at present: U( N- I" f% [+ ^
请问为什么?如果能直接帮我解出来将不甚感激! |
|