作业帮 > 数学 > 作业

计算含有符号的矩阵方程组的解,使用mathematica或者matlab都可以.

来源:学生作业帮 编辑:搜搜做题作业网作业帮 分类:数学作业 时间:2024/04/30 02:54:58
计算含有符号的矩阵方程组的解,使用mathematica或者matlab都可以.
我所需要计算的方程组有很多方程,方程中的未知数有很多,而且部分是向量和矩阵形式,好在基本都是线性的,最后会得到一组唯一的解.希望有高手不吝赐教,传授方法,非常感激!这个问题可能能够解决的人会比较少,悬赏分给出去了很可能就流失了,分数什么的就是一个数字而已,我就不到300分,愿意拿去就全部奉上!
方程和变量都在这里了
syms E e q r r1 r2 r3 oX oY oZ pX1 pY1 pZ1 pX2 pY2 pZ2 pX3 pY3 pZ3 xl yl zl xm ym zm xn yn znE=0.06;e=0.06;q=sqrt(3);r=0.0625;r1=0.0000036317;r2=0.000003753;r3=0.0000051272;
M1=[q*E/2;-E/2;0];
M2=[0;E;0];
M3=[-q*E/2;-E/2;0];
m1m=[q*E/2;-E/2;0];
m2m=[0;E;0];
m3m=[-q*E/2;-E/2;0];
o=[oX;oY;oZ];
p1=[pX1;pY1;pZ1];
p2=[pX2;pY2;pZ2];
p3=[pX3;pY3;pZ3];
R=[xl yl zl;xm ym zm;xn yn zn];
m1=R*m1m+o;m2=R*m2m+o;
m3=R*m4m+o;
(r+r1)^2=dot((m1-M1),(m1-M1));
(r+r2)^2=dot((m2-M2),(m2-M2));
(r+r3)^2=dot((m3-M3),(m3-M3));
1=dot(p1,p1);
1=dot(p2,p2);
1=dot(p3,p3);
(q*e)^2=dot((m1-m2),(m1-m2));
(q*e)^2=dot((m3-m2),(m3-m2));
(q*e)^2=dot((m1-m3),(m1-m3));
o=(m1+m2+m3)/3;
m1=M1+(r+r1)*p1;
m2=M2+(r+r2)*p2;
m3=M3+(r+r3)*p3;
dot(p1,(m2-m3))=0;
dot(p2,(m1-m3))=0;dot(p3,(m2-m1))=0;
计算含有符号的矩阵方程组的解,使用mathematica或者matlab都可以.
这种东西你不把实际方程给出来别人没法解答……如果是线性方程组的话,只要你语法不要错,Mathematica的Solve指令一般都可以解的——不过……看你的说法,你在软件方面是0基础?
再问: 好吧,您是老师,0基础就0基础吧,马上奉上题目。
再答: 有几个不明确的地方。 1 你的q是已知数吗? 2 你的小m向量右上角的角标m是什么意思? 3 小o向量右上角的角标M又是什么意思? 4 左侧写着M和m的那个R,这种表示法又有什么特殊含义?右边的九个量和其他未知量的依赖关系是啥? 5 你的式子里应该没有用到张量的求和约定吧? 6 字太小了我也许没太看清楚,不过你是不是把Δr写了两遍?顺便这个和r是什么关系?总不会是散度吧?
再问: 本来想加你百度hi的,可是您没回应,那就把方程贴上来了。第一排最后的znE中间有一个转行,zn是syms的最后一个,E=0.06。我的网络很不给力,所以上起网来有点困难,不好意思啊!您只需要帮我把格式改一下,成为正确的程序即可,谢谢啊!
再答: Malab……算了,这种程度我也懂。你的m1,m2,m3的定义式写错了吧,应该是内积吧?然后,这个方程无解。你仔细检查检查吧。
再问: 可能是我的方程没整理完全,没有解我回头再整理一下,方程挺多的,你能把正确的程序告诉我么?mathematica也可以,不要求解出答案,只想要一个完整的程序。
再答: 发给你也可以,但是你连基础语法都不明白的话……也罢,你多看看帮助吧。直接贴怕被吞,我用私信给你吧,唉,其实很简单的。