QQ登录

只需一步,快速开始

登录 | 注册 | 找回密码

三维网

 找回密码
 注册

QQ登录

只需一步,快速开始

展开

通知     

查看: 1295|回复: 0
收起左侧

[讨论结束] 一个matlab问题

 关闭 [复制链接]
发表于 2007-4-30 01:07:27 | 显示全部楼层 |阅读模式

马上注册,结识高手,享用更多资源,轻松玩转三维网社区。

您需要 登录 才可以下载或查看,没有帐号?注册

x
一个matlab问题' |4 g  A& f$ l. K; ~) X
我得到一个矩阵,但是只能求出他的本征值,无法得到他的本征向量
8 B3 x) ^3 ^8 t9 d# c) t知道为什么吗?1 b! \7 ?* [9 o, y$ P3 ^, p
我是用matlab算的
1 u  m, _: ^+ u, q! |# n0 ysyms j b q p u;. f) N  u" s4 e5 c
x=[0 1;1 0];2 S3 ^6 Y& n2 m. b
y=[0 -i;i 0];
; q* s3 v; a  p, G8 _: sz=[1 0;0 -1];1 P/ Y  }5 d7 M8 w' h" [
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])
( m+ v  d5 z' E; l& ^. a7 z可以得到h矩阵,
- }4 V* Q' e( x$ X9 O: U用eig(h) 很快得到了他的本征值1 C. R. T3 O4 k# e  y, n' a0 P# B5 o
但是用[v d]=eig(h)时就会出错,
+ {9 m/ O, Z2 D! a% l# o但是我也想知道他的本征向量v
6 K0 ~: o- m! @3 P' Q! ~8 _错误原因是:: a* I4 X9 R0 z; l$ x1 p
?? Error using ==> sym.eig
% @) K' w% g  q" w7 K% WError, (in eigenvectors) eigenvects only works for a matrix of rationals, rational functions, algebraic numbers, or algebraic functions at present
$ _& D9 S) Z. p% R2 B/ X9 D$ d" j请问为什么?如果能直接帮我解出来将不甚感激!
发表回复
您需要登录后才可以回帖 登录 | 注册

本版积分规则

Licensed Copyright © 2016-2020 http://www.3dportal.cn/ All Rights Reserved 京 ICP备13008828号

小黑屋|手机版|Archiver|三维网 ( 京ICP备2023026364号-1 )

快速回复 返回顶部 返回列表