QQ登录

只需一步,快速开始

登录 | 注册 | 找回密码

三维网

 找回密码
 注册

QQ登录

只需一步,快速开始

展开

通知     

全站
10天前
查看: 1298|回复: 0
收起左侧

[讨论结束] 一个matlab问题

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

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

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

x
一个matlab问题/ m2 W7 Z2 g  K* z3 S7 l' `: ^
我得到一个矩阵,但是只能求出他的本征值,无法得到他的本征向量
! g& E3 z! ]! n9 O' i* Q知道为什么吗?: q! F  Y, }, J
我是用matlab算的
+ p, T  h5 M4 {) _& W* G* Lsyms j b q p u;
% V6 U6 B; @0 V& j8 A; P8 f& ox=[0 1;1 0];
! Q1 h1 c* f& f# a6 Y6 D8 d- a4 gy=[0 -i;i 0];  Z6 K( Y( V- O/ a
z=[1 0;0 -1];
: I6 m0 n6 X! p9 ~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])$ x7 I0 m8 `, H3 |# O
可以得到h矩阵,* I% E( C! o9 A% K
用eig(h) 很快得到了他的本征值
9 z& a: {+ V+ G0 U  E( b8 j. i但是用[v d]=eig(h)时就会出错,! E0 h9 f( K, f
但是我也想知道他的本征向量v
% m5 `% z2 |& z$ X错误原因是:+ W+ {0 r* ]. L0 N; ^; M  m- x
?? Error using ==> sym.eig
$ z4 D+ f0 C# DError, (in eigenvectors) eigenvects only works for a matrix of rationals, rational functions, algebraic numbers, or algebraic functions at present5 j$ I! K4 w: a  ^. R" M
请问为什么?如果能直接帮我解出来将不甚感激!
发表回复
您需要登录后才可以回帖 登录 | 注册

本版积分规则

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

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

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