MATHCAD解2次方程组的问题
方程在图片里请问像这样的方程MATHCAD是否能解出来？如果可以 请问格式该怎么写？
如果不能 请问该用什么样的软件步骤该如何谢谢 题目十分繁复，是否可以发一个mathcad的xmcd文件的压缩件上来，以节省移植时间。 2# zpz77777
这个就是刚刚那个方程组的谢谢版主帮忙 本帖最后由 zsq-w 于 2009-6-15 19:55 编辑
我初学，随意说几句。
运行的lz文件，显示 solution not found。
lz的方程非常复杂，建议先在草稿纸上化简再在mcad中运算，这样比较好把握。
v1 v2 v3以比值形式出现在方程中，不妨就用比值，另赋变量，这样或许看起来方便点。 4# zsq-w
额 我再去弄弄基本什么方法都用过了 就是说无解。。。。所以才来问问 另外，手工推导一下，说不定你会发现意想不到的一些结论（help if u are doing some research in this field）。我隐隐感觉到这几个未知数之间有一定的对称性和对应关系。 6# zsq-w
.....手工………… 至少是5元2次方程额 这个公式是文献上报道过的只是没有说明详细的计算方法。。。。实在是不知道怎么弄出来的
现在就差解这一步了就是解不出来。。干着急 by the way, where is the value of g12... 可能你误会我的意思了，不是说手工完全解出，而是你的方程组可以简化。
可以肯定的说，你那6个方程组中的delta/RT（就是方程左边那些分式），那6个中间未知数是完全不必要的，最好合并同类项，再来求解。而且合并同类项的时候，我相信你能对你的问题产生更深刻的认识。 用1stOpt试了下，不知对否？delta和theat分别用d和c代替，代码如下，5个方程，6个参数：
Constant V1=18E+6, V2=1.385E+8, V3=5.609E+10;
Constant c3b=0.01, R=8.314, T=298.15, g23=-1.108E+6, g13=0.61;
ConstStr Du1a=(ln(c1a)+1-c1a-V1/V2*c2a-V1/V3*c3a+(g12*c2a+g13*c3a)*(c2a+c3a)-V1/V2*g23*c2a*c3a)*R*T,
Du2a=(ln(c2a)+1-c2a-V2/V1*c1a-V2/V3*c3a+(g12*V2/V1*c1a+g23*c3a)*(c1a+c3a)-V2/V1*g13*c1a*c3a)*R*T,
Du3a=(ln(c3a)+1-c3a-V3/V1*c1a-V3/V2*c2a+(g13*V3/V1*c1a+V3/V2*g23*c2a)*(c2a+c1a)-V3/V1*g12*c1a*c2a)*R*T,
Du1b=(ln(c1b)+1-c1b-V1/V2*c2b-V1/V3*c3b+(g12*c2b+g13*c3b)*(c2b+c3b)-V1/V2*g23*c2b*c3b)*R*T,
Du2b=(ln(c2b)+1-c2b-V2/V1*c1b-V2/V3*c3b+(g12*V2/V1*c1b+g23*c3b)*(c1b+c3b)-V2/V1*g13*c1b*c3b)*R*T,
Du3b=(ln(c3b)+1-c3b-V3/V1*c1b-V3/V2*c2b+(g13*V3/V1*c1b+V3/V2*g23*c2b)*(c2b+c1b)-V3/V1*g12*c1b*c2b)*R*T;
Function du1a=du1b;
du2a=du2b;
du3a=du3b;
c1a+c2a+c3a=1;
c1b+c2b+c3b=1;
可得（解不是唯一的）：
c1a: 0.706028200611506
c2a: 0.283968904297105
c3a: 0.0100000000000069
g12: 2.61001641220209
c1b: 0.706028200611907
c2b: 0.283968904297016
由上面解可看出，解的确是对称的，即c1a=c1b, c2a=c2b, c3a=c3b。如果将此条件带入再求解，就非常容易了。 额。。。。。其实我也发现了。。。g12忘记赋值了解出来确实是两两相等的
我是给5个未知数都赋了初始值 然后求解方程的数值解 用＝号的 确实可以求出来
不过结果跟我想要的答案不太一样罢了
谢谢大家了 本帖最后由 xhxforever 于 2009-6-18 09:42 编辑
10# shamohu
额～～这个算出来的g12 跟我想要赋的值很接近。。。。。
如果直接给g12赋值后能否解出其他6个未知数就是说不给c3b复制的话能解出来不
我不太清楚非线性方程怎么解
如果能解出多个解我是希望能用这些解画出一条线来说明问题的谢谢版主了 如果设g12=0.01，c3b为待求参数，用1stOpt也可很容易得到6个参数的解，并且解也是对称的，但不是唯一的。下面是一组：
c1a: 0.745560900726516
c2a: 0.254434054272094
c3a: 2.75696800642512E-6
c1b: 0.745560899096907
c2b: 0.254434054602625
c3b: 2.7569680300853E-6 13# shamohu
版主能不能传一个 MFF的格式 上来 这个方程组的 谢谢
我想试试你说的那个软件看看怎么求解
页:
[1]