| 查看: 155 | 回复: 1 | |||
| 当前主题已经存档。 | |||
jily木虫 (正式写手)
|
[交流]
Mathematica求解的问题
|
||
|
一个有三个方程的方程组，其中包含 dw，dT，dN和dq四个变量，需要借助于Mathematica化简该方程组，使之成为分别以dw/dq=..., dT/dq=..., dN/dq=...表示的三个方程，具体如附件的图片所示。 按照作者所提供的Mathematica源程序运行，总是提示所输入的公式存在“is not a well-formed equation.”错误。由于本人对该软件完全不懂，求高人指点。 顺便问一下，Matlab里有能实现上述功能的函数？？若有，还望指点一二。谢谢！ 附上Mathematica源程序： Clear[w,q,y,a,b,u,T,p] -u*dN*a*Cos[w]*Cos[q]+dN*Sin[q]+u*dN*b*Sin[w]*Cos[q]-T*dq*Cos[w]*Sin[q]+dT*Cos[w]*Cos[q]-T*dw*Sin[w]*Cos[q]==0 u*dN*a*Sin[w]+u*dN*b*Cos[w]-T*dw*Cos[w]-dT*Sin[w]==0 u*dN*a*Cos[w]*Sin[q]+dN*Cos[q]+u*dN*b*Sin[w]*Sin[q]-T*dq*Cos[w]*Cos[q]+T*dw*Sin[w]*Sin[q]-dT*Cos[w]*Sin[q]==0 Solve[{%#,%#,%#},{dw,dT,dN}] Simplify[%#]  [ Last edited by jily on 2009-12-21 at 19:18 ] |
回复此楼» 猜你喜欢
假如你的研究生提出不合理要求
已经有8人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
实验室接单子
已经有4人回复
-
全日制（定向）博士
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
★ ★ ★ ★ ★ ★
jily(金币+6,VIP+0):谢谢关注！同时鼓励新虫。:-) 这个问题已通过其他途径解决了。 1-11 13:29
jily(金币+6,VIP+0):谢谢关注！同时鼓励新虫。:-) 这个问题已通过其他途径解决了。 1-11 13:29
|
感觉是solve命令中的问题。试试下面的命令，看能不能运算 先让3个等式都除以dq Clear[w,q,y,a,b,u,T,p] (*Comment: dNq=dN/dq, dwq=dw/dq,dTq=dT/dq*) eq1=-u*dNq*a*Cos[w]*Cos[q]+dNq*Sin[q]+u*dNq*b*Sin[w]*Cos[q]-T*Cos[w]*Sin[q]+dTq*Cos[w]*Cos[q]-T*dwq*Sin[w]*Cos[q]; eq2=u*dNq*a*Sin[w]+u*dNq*b*Cos[w]-T*dwq*Cos[w]-dTq*Sin[w]; eq3=u*dNq*a*Cos[w]*Sin[q]+dNq*Cos[q]+u*dNq*b*Sin[w]*Sin[q]-T*Cos[w]*Cos[q]+T*dwq*Sin[w]*Sin[q]-dTq*Cos[w]*Sin[q]; Solve[{eq1==0,eq2==0,eq3==0},{dNq,dwq,dTq}] PS: symbolic的运算Matlab要用symbol定义变量，不如直接用Mathematica或者Maple，除非没有这些软件。 |
2楼2009-12-31 00:11:09