QQ登录

只需一步,快速开始

登录 | 注册 | 找回密码

三维网

 找回密码
 注册

QQ登录

只需一步,快速开始

展开

通知     

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

[讨论结束] 一个matlab问题

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

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

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

x
一个matlab问题
- ]; q' a9 a) }我得到一个矩阵,但是只能求出他的本征值,无法得到他的本征向量" O  |0 {6 G# |" z- e3 ~' W
知道为什么吗?: t  n5 a( _8 r7 B( D: j
我是用matlab算的, N  k! b$ \. l( Z/ f  M
syms j b q p u;
5 M+ X+ H1 a, A  Jx=[0 1;1 0];
+ @5 t" f0 V" s- j$ ky=[0 -i;i 0];
2 h+ p9 ]7 P- h$ ez=[1 0;0 -1];8 l# T& w: K/ _2 p0 Y
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])- F. t* t! h  a" k. F" y
可以得到h矩阵,
4 C# l) k) x! X  r用eig(h) 很快得到了他的本征值
: Z0 B' v# n( b0 l但是用[v d]=eig(h)时就会出错,
* x  ^$ v) j' N  n8 q9 F但是我也想知道他的本征向量v' x3 O# Q" M+ Q- G: L$ O7 |) I/ X
错误原因是:0 z9 K; V5 T5 I3 P+ o) |% g) d
?? Error using ==> sym.eig9 c3 d- A! {% X: w& n: W
Error, (in eigenvectors) eigenvects only works for a matrix of rationals, rational functions, algebraic numbers, or algebraic functions at present
# C; L' J, F! s5 t+ z请问为什么?如果能直接帮我解出来将不甚感激!
发表回复
您需要登录后才可以回帖 登录 | 注册

本版积分规则


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

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

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