| 查看: 240 | 回复: 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 ] |
回复此楼» 猜你喜欢
FileCode能看出啥？
已经有19人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有32人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有5人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有10人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有13人回复
-
是这周出结果还是下周出结果？
已经有3人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有9人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有4人回复
-
帮忙看看fileCode
已经有13人回复
★ ★ ★ ★ ★ ★
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
