| 查看: 437 | 回复: 0 | ||
[求助]
急用matlab求解一个四元非线性方程组！能帮忙编程的悬赏100个金币！
|
|
这是一个关于三元A-B-C相图的计算问题， 关系式g_ex1=29433*x1*x2+29241*(x2-x1)*x1*x2-43.65*T*x1*x2-33.85*T*(x2-x1)*x1*x2;这是A-B二元系的过剩自由能表达式 关系式g_ex2=7777*x1*x3+3447*(x3-x1)*x1*x3-21.46*T*x1*x3;这是A-C二元系的过剩自由能表达式 关系式g_ex3=-1906*x2*x3+32*x2*x3*(x3-x2)-0.874*T*x2*x3;这是B-C二元系的过剩自由能表达式 关系式G_ex=x2/(x2+x3)*g_ex1+x3/(x2+x3)*g_ex2+(1-x1)*(1-x1)*g_ex3;这是A-B-C三元系的过剩自由能表达式 关系式G_1=G_ex-x2*diff(G_ex,x2)-x3*diff(G_ex,x3);这是A的偏摩尔过剩自由能表达式 关系式G_2=G_ex+(1-x2)*diff(G_ex,x2)-x3*diff(G_ex,x3);这是B的偏摩尔过剩自由能表达式 关系式G_3=G_ex-x2*diff(G_ex,x2)+(1-x3)*diff(G_ex,x3);这是C的偏摩尔过剩自由能表达式 这是最后得到的等量方程 F1=3796.7+18.64*T-0.0279*T*T+G_1+8.314*T*ln(x1)=0; F2=25563-48.415*T+G_2+8.314*T*ln(x2)=0; F3=15177-26.033*T+G_3+8.314*T*ln(x3)=0; F4=1-x1-x2-x3=0; 其中，温度T的范围可能是从323~834K，而x1,x2,x3的初值分别为0.3,0.4,0.3，如何采用牛顿-拉普森法求解该非线性方程组，其中F4代表的是x1,x2,x3需要满足的关系。 |
回复此楼» 猜你喜欢
青基代表作，AAAI之类的A会的special track在国内认可度高吗？还是归为workshop之流？
已经有3人回复
-
上海工程技术大学【激光智能制造】课题组招收硕士
已经有6人回复
-
带资进组求博导收留
已经有11人回复
-
自荐读博
已经有5人回复
-
求个博导看看
已经有16人回复
-
上海工程技术大学张培磊教授团队招收博士生
已经有4人回复
-
求助院士们，这个如何合成呀
已经有4人回复
-
临港实验室与上科大联培博士招生1名
已经有9人回复
-
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？
已经有6人回复
-
最近几年招的学生写论文不引自己组发的文章
已经有11人回复
找到一些相关的精华帖子，希望有用哦~
初学新手求教。求用Matlab求解一个非线性回归，或请指点门路。
已经有11人回复
- 用1stOpt解出非线性方程组的解作为初值用Matlab求解，解不出来
已经有4人回复
- MATLAB求解五元非线性方程组，在线等
已经有7人回复
- 非线性二阶微分方程组求解，matlab
已经有12人回复
- matlab 求解非线性方程组未果
已经有10人回复
- 用matlab求解非线性方程组说无解，一定是方程组本身无解，还是有可能程序有问题呢？
已经有11人回复
- matlab如何求解一个非线性微分方程组
已经有8人回复
- MATLAB求解非线性方程组
已经有5人回复
- matlab求解非线性方程组，错误提示怎么解决
已经有5人回复
- 请教用matlab求解一个非线性偏微分方程组的数值解
已经有8人回复
- matlab求解非线性方程组，求助！
已经有6人回复
- 求助！matlab用fsolve函数求解非线性方程组的问题！
已经有19人回复
- matlab能否画出这样一个三维球形图？？求指教～～急用！！
已经有10人回复
- 如何使用matlab求解非线性方程组的所有整数解？
已经有9人回复
- matlab求解非线性方程组
已经有16人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- 【求助】matlab求解非线性方程组，并画图处理。要求y,z是实数解！
已经有18人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人