登 录
註 冊
论坛
微波仿真网
注册
登录论坛可查看更多信息
微波仿真论坛
>
程序
>
一个matlab问题
发帖
回复
2705
阅读
0
回复
[
求助
]
一个matlab问题
离线
pkuwuli
UID :2199
注册:
2007-04-26
登录:
2023-09-01
发帖:
174
等级:
仿真三级
0楼
发表于: 2007-04-30 00:52:48
— 本帖被 casey 从 其它软件 移动到本区(2008-05-28) —
我得到一个矩阵,但是只能求出他的本征值,无法得到他的本征向量
(f#b7O-Wn
知道为什么吗?
f.xA_Y>
我是用matlab算的
VaLs`q&3>
syms j b q p u;
[hh/1[
x=[0 1;1 0];
.*5 Z"Q['G
y=[0 -i;i 0];
dh -,E
z=[1 0;0 -1];
Es4qPB`g.
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])
]?LB?:6
可以得到h矩阵,
vjUp *R>h
用eig(h) 很快得到了他的本征值
\ |!\V
但是用[v d]=eig(h)时就会出错,
eJ@~o{,?>
但是我也想知道他的本征向量v
9pD 7 f`
WAPhv-6
错误原因是:
L%\b' fs
8P: spD0
??? Error using ==> sym.eig
K8[DZ)rO;Z
Error, (in eigenvectors) eigenvects only works for a matrix of rationals, rational functions, algebraic numbers, or algebraic functions at ..
^@6q
{/8Q)2*>0
未注册仅能浏览
部分内容
,查看
全部内容及附件
请先
登录
或
注册
共
条评分
发帖
回复