QQ登录

只需一步,快速开始

登录 | 注册 | 找回密码

三维网

用户名  找回密码
 注册

QQ登录

只需一步,快速开始

展开

通知     

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

[讨论结束] 一个matlab问题

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

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

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

x
一个matlab问题- W8 H: C' W+ Z6 w6 j% c  U- ~
我得到一个矩阵,但是只能求出他的本征值,无法得到他的本征向量
) b. _) W8 ?. g. w) t3 k知道为什么吗?+ f: F( u3 O7 P3 ~, T& B* `
我是用matlab算的; [% F, A+ k) M+ g7 _" _4 e. f
syms j b q p u;, ~; ^- X8 U/ L: c+ h
x=[0 1;1 0];
- @- ?" S, O* P/ r/ Y2 c" E# j: wy=[0 -i;i 0];
5 R( ^9 r* J  Q7 o; L( k& X. ez=[1 0;0 -1];
0 |% X$ x8 w3 ~. q, U- `5 B% jh=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])+ K7 F$ R: Q* w- ^- N% J
可以得到h矩阵,
9 a2 j' a0 p8 n# p8 G7 l用eig(h) 很快得到了他的本征值
& d5 b( t, |0 |' M3 Z9 O但是用[v d]=eig(h)时就会出错,
- @5 d. l' i& P但是我也想知道他的本征向量v7 B6 q+ d/ K' x' q$ }4 j1 ?
错误原因是:
9 c3 |9 V3 L( w) z, D1 e' e! i?? Error using ==> sym.eig
+ r& I9 F4 w. I5 bError, (in eigenvectors) eigenvects only works for a matrix of rationals, rational functions, algebraic numbers, or algebraic functions at present
+ H" c' v9 L5 m. G) S) c% D请问为什么?如果能直接帮我解出来将不甚感激!
发表回复
您需要登录后才可以回帖 登录 | 注册

本版积分规则

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

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

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