| 查看: 1381 | 回复: 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 |
回复此楼» 猜你喜欢
计算开关的响应时间
已经有0人回复
-
今年面上预计又打酱油了
已经有18人回复
-
物理学I论文润色/翻译怎么收费?
已经有196人回复
-
北伊利诺伊大学(Northern Illinois University)机械系招收全奖博士生
已经有0人回复
-
北伊利诺伊大学(Northern Illinois University)机械系招收全奖博士生
已经有0人回复
-
第五届粉末衍射晶体结构解析及Rietveld精修高级研讨会
已经有44人回复
-
FDTD软件里的时域图像画图
已经有2人回复
-
复杂量子系统的新物态和新机理讲习班
已经有0人回复
-
导电原子力显微镜cAFM如何测试mcAFM(即磁导AFM)性能？
已经有0人回复
-
感谢小木虫，已完成目标
已经有0人回复
-
零点能
已经有2人回复
hzlhm
至尊木虫 (著名写手)
- 应助: 373 (硕士)
- 金币: 16237.3
- 红花: 52
- 帖子: 2577
- 在线: 572.8小时
- 虫号: 1517335
- 注册: 2011-11-30
- 性别: GG
- 专业: 常微分方程与动力系统
5楼2021-05-26 21:39:25
独孤神宇
版主 (知名作家)
- 应助: 489 (硕士)
- 贵宾: 0.008
- 金币: 29608.8
- 散金: 802
- 红花: 121
- 沙发: 1
- 帖子: 5354
- 在线: 850.4小时
- 虫号: 3522474
- 注册: 2014-11-06
- 性别: GG
- 专业: 机械动力学
- 管辖: 计算模拟
2楼2021-05-26 11:07:51
 本帖仅楼主可见 |
3楼2021-05-26 11:46:39
|
上边写错了。。是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
50