Labfans是一个针对大学生、工程师和科研工作者的技术社区。 论坛首页 | 联系我们(Contact Us)
MATLAB爱好者论坛-LabFans.com
返回   MATLAB爱好者论坛-LabFans.com > 工程数学软件 > MATLAB论坛
MATLAB论坛 一切MATLAB相关问题在此讨论。
回复
 
主题工具 显示模式
旧 2013-04-23, 17:47   #1
ziyoubianjun
初级会员
 
注册日期: 2013-04-23
帖子: 2
声望力: 0
ziyoubianjun 正向着好的方向发展
警告 急求一道求解矩阵特征值的问题!!!!

我想求解A矩阵的特征值,目的是求它的对角矩阵,我编写的代码如下:
clc;
syms J1;
syms J2;
syms J3;
syms x1 x2 x3 x4 x5;
syms m g l;
A=[0,(J2-J3)*x3/(2*J1),(J2-J3)*x2/(2*J1),0,-m*g*l/(J1),0;
(J3-J1)*x3/(2*J2),0,(J3-J1)*x1/(2*J2),m*g*l/J2,0,0;
0,0,0,0,0,0;
0,0,x5/2,0,x3/2,-x2;
0,0,-x4/2,-x3/2,0,x1;
-x5/2,x4/2,0,-x2/2,x1/2,0];
[V,D]=eig(A);
最后总是运行不出结果,总是显示:
??? Error using ==> mupadmex
Error in MuPAD command: cannot compute the explicit representation of the eigenvalues; use 'numeric::eigenvectors'
[linalg::eigenvectors]

Error in ==> sym.sym>sym.mupadmexnout at 2003
out = mupadmex(fcn,args{:});

Error in ==> sym.eig at 68
[V,D,p] = mupadmexnout('mllib::eigenvectors',A);

Error in ==> Untitled at 13
[V,D]=eig(A);


哪位大仙知道这事怎么回事啊!!!急求!!!
ziyoubianjun 当前离线   回复时引用此帖
旧 2013-04-23, 22:21   #2
ziyoubianjun
初级会员
 
注册日期: 2013-04-23
帖子: 2
声望力: 0
ziyoubianjun 正向着好的方向发展
警告 回复: 急求一道求解矩阵特征值的问题!!!!

有没有大师解答啊!急死我了!~~~
ziyoubianjun 当前离线   回复时引用此帖
回复

主题工具
显示模式

发帖规则
不可以发表新主题
不可以发表回复
不可以上传附件
不可以编辑自己的帖子

启用 BB 代码
论坛启用 表情符号
论坛启用 [IMG] 代码
论坛禁用 HTML 代码



所有时间均为北京时间。现在的时间是 18:30


Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.