matlab求解符号变量方程组

clear;
clc;
womg=[0.5 -0.5 0.4]';
j=[1 2.8 2]';
J=diag(j);
syms a11 a12 a13 a21 a22 a23 a31 a32 a33
A=[0 -womg(3) womg(2);womg(3) 0 -womg(1);-womg(2) womg(1) 0];
E=J*womg;
B=[0 -E(3) E(2);E(3) 0 -E(1);-E(2) E(1) 0];
A1=[a11 a12 a13;a21 a22 a23;a31 a32 a33];
E1=[a11 a21 a31;a12 a22 a32;a13 a23 a33];
H=A1*A-A*E1;
% B=sym(D);
% F=H-D;
% F1=eval(F);
eq1=H(1)-B(1);
eq2=H(2)-B(2);
eq3=H(3)-B(3);
eq4=H(4)-B(4);
eq5=H(5)-B(5);
eq6=H(6)-B(6);
eq7=H(7)-B(7);
eq8=H(8)-B(8);
eq9=H(9)-B(9);
S=solve(eq1,eq2,eq3,eq4,eq5,eq6,eq7,eq8,eq9,a11,a12,a13,a21,a22,a23,a31,a32,a33);
B=[S.a11,S.a12,S.a13;S.a21,S.a22,S.a23;S.a31,S.a32,S.a33];
运行的时候提示错误Warning: Explicit solution could not be found.
> In solve at 140
In sym.solve at 49
In fuhao at 25
??? Access to an object's fields is only permitted within its methods.

Error in ==> fuhao at 26
B=[S.a11,S.a12,S.a13;S.a21,S.a22,S.a23;S.a31,S.a32,S.a33];
不知道该怎么改,请高手帮忙

第1个回答  2011-04-03
你的方程组没有解析解
所以出错
你看看前面程序是不是构成方程组有问题追问

这是个九元一次方程组,而且是线性的,应该有解的吧

追答

Explicit solution could not be found
没有解析解

本回答被提问者采纳
相似回答