QQ登录

只需一步,快速开始

登录 | 注册 | 找回密码

三维网

 找回密码
 注册

QQ登录

只需一步,快速开始

展开

通知     

全站
goto3d 说: 此次SW竞赛获奖名单公布如下,抱歉晚了,版主最近太忙:一等奖:塔山817;二等奖:a9041、飞鱼;三等奖:wx_dfA5IKla、xwj960414、bzlgl、hklecon;请以上各位和版主联系,领取奖金!!!
2022-03-11
系统
[系统通知] 平台第一个项目外包——项目拼多;正式上线,欢迎各单位个人有外包、设计、采购、加工需求的,在此寻找更牛的解决方案
2021-07-01
查看: 1266|回复: 0
收起左侧

[讨论结束] 一个matlab问题

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

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

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

x
一个matlab问题
* q/ T8 C; ~' w1 V' |/ A: @我得到一个矩阵,但是只能求出他的本征值,无法得到他的本征向量4 C, C: h; y' ]8 ^
知道为什么吗?- {" P6 X& N* o
我是用matlab算的: |0 G( ]  F; l, C
syms j b q p u;
+ M; @/ x1 x! px=[0 1;1 0];
9 K) h- T7 d/ g2 F8 Hy=[0 -i;i 0];% v9 \: m( {% i" r9 x% Q) O
z=[1 0;0 -1];
4 X0 D8 t, M. t, W$ x; Bh=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 R" `6 Z: V5 G, s- b- e可以得到h矩阵,4 [$ G2 |. k- J: [& F8 ~/ a! ?
用eig(h) 很快得到了他的本征值
0 }( `) {' K$ H) \' m但是用[v d]=eig(h)时就会出错,
  g0 x: H$ p8 B! j但是我也想知道他的本征向量v. s8 j7 C! B1 L0 D% X4 k
错误原因是:
5 ~' o7 C: `5 [3 Y?? Error using ==> sym.eig
) Q3 p3 P  B: [% @- nError, (in eigenvectors) eigenvects only works for a matrix of rationals, rational functions, algebraic numbers, or algebraic functions at present/ O- I# @& r4 K' J
请问为什么?如果能直接帮我解出来将不甚感激!
发表回复
您需要登录后才可以回帖 登录 | 注册

本版积分规则

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

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

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