| 查看: 1368 | 回复: 6 | |||
| 【悬赏金币】回答本帖问题，作者09天岚将赠送您 5 个金币 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
NRTL方程拟合敛到参数全为0已有2人参与
|
|||
|
各位大佬们，我在1stopt上拟合固液相平衡的NRTL方程，为什么收敛到参数全为0啊。。。有没有大佬帮忙看看方程哪里出了问题。。。谢谢了！ 数据是找的文献的数据。 拿这套数据拟合其他模型都没什么问题。。 parameters a,b,c,r; variable y,x; function ln(r)=exp(((1-x)^2)*((b/(8.314*y)*exp(c*a/(8.314*y))/((x+(1-x)*exp(c*a/(8.314*y)))^2))+((a/(8.314*y))*exp(c*b/(8.314*y)))/(((1-x)+(1-x)*exp(c*b/(8.314*y)))^2))); data; //y,x 278.15 0.001541 283.15 0.002004 288.15 0.002704 293.15 0.003392 298.15 0.004415 303.15 0.005561 308.15 0.007136 313.15 0.009151 318.15 0.01113 图为NRTL方程原形式。 急急急 求求大家帮忙分析分析！  PWN3F1)AM5{5W4{)S}0X1OB.png |
回复此楼» 猜你喜欢
SCAPS-1D学习交流
已经有2人回复
-
深圳大学跟哈工大联合博士
已经有13人回复
-
物理学I论文润色/翻译怎么收费?
已经有75人回复
-
Journal of alloys and compounds期刊
已经有2人回复
-
请教一个概率问题
已经有0人回复
-
铁氧体材料在太赫兹波段的响应时间
已经有0人回复
-
optisworks 的软件和教程
已经有0人回复
-
求助FAPbI3 cif文件
已经有0人回复
-
计算开关的响应时间
已经有0人回复
-
今年面上预计又打酱油了
已经有8人回复
|
上边写错了。。是function r=exp(((1-x)^2)*((b/(8.314*y)*exp(c*a/(8.314*y))/((x+(1-x)*exp(c*a/(8.314*y)))^2))+((a/(8.314*y))*exp(c*b/(8.314*y)))/(((1-x)+(1-x)*exp(c*b/(8.314*y)))^2))); 发自小木虫Android客户端 |
4楼2021-05-26 12:31:26
独孤神宇
版主 (知名作家)
- 应助: 488 (硕士)
- 贵宾: 0.008
- 金币: 29561.8
- 散金: 802
- 红花: 120
- 沙发: 1
- 帖子: 5342
- 在线: 850.2小时
- 虫号: 3522474
- 注册: 2014-11-06
- 性别: GG
- 专业: 机械动力学
- 管辖: 计算模拟
2楼2021-05-26 11:07:51
 本帖仅楼主可见 |
3楼2021-05-26 11:46:39
hzlhm
至尊木虫 (著名写手)
- 应助: 373 (硕士)
- 金币: 16197.8
- 红花: 52
- 帖子: 2571
- 在线: 572.3小时
- 虫号: 1517335
- 注册: 2011-11-30
- 性别: GG
- 专业: 常微分方程与动力系统
5楼2021-05-26 21:39:25
5