Labfans是一个针对大学生、工程师和科研工作者的技术社区。 | 论坛首页 | 联系我们(Contact Us) |
![]() |
|
![]() |
#1 |
初级会员
注册日期: 2013-04-23
帖子: 2
声望力: 0 ![]() |
![]()
我想求解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); 哪位大仙知道这事怎么回事啊!!!急求!!! |
![]() |
![]() |
![]() |
#2 |
初级会员
注册日期: 2013-04-23
帖子: 2
声望力: 0 ![]() |
![]()
有没有大师解答啊!急死我了!~~~
|
![]() |
![]() |