应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜     石溪大学接受考研调剂申请>

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 理论与计算方法 » NRTL方程拟合敛到参数全为0
老同志的第一次SCI投稿经历，对新人极具参考价值！
51/1返回列表
查看: 1367  |  回复: 6
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者09天岚将赠送您 5 个金币
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

09天岚

铁虫 (初入文坛)

[求助] 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方程原形式。
急急急   求求大家帮忙分析分析！

NRTL方程拟合敛到参数全为0
PWN3F1)AM5{5W4{)S}0X1OB.png
回复此楼

» 猜你喜欢
1楼2021-05-26 10:59:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hzlhm

至尊木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
独孤神宇: 金币+5, 鼓励交流 2021-05-26 21:52:01
引用回帖:
4楼: Originally posted by 09天岚 at 2021-05-26 12:31:26
上边写错了。。是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)));
...

利用matlab的lsqcurvefit函数，可以得到
a= 0.088241733690097801 ，b= -0.052836750140336099 ，c= -0.024187951865466500 ，r=0
NRTL方程拟合敛到参数全为0-1
123a.jpg
赞一下 回复此楼
QQ：2120156492
5楼2021-05-26 21:39:25
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

独孤神宇

版主 (知名作家)

【答案】应助回帖

感谢参与，应助指数 +1
方程写成y=……
因变量在方程左边，变量，数据也要对应

发自小木虫Android客户端
赞一下 回复此楼
数值计算
2楼2021-05-26 11:07:51
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名
本帖仅楼主可见
3楼2021-05-26 11:46:39
 已阅   申请计算强帖   回复此楼   编辑   查看我的主页

09天岚

铁虫 (初入文坛)
上边写错了。。是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
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫 SCI英文润色翻译-EnPapers美团队-先修改后付款
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国家重点实验室基金，通知去答辩。希望大不？ 20+3 qweasd12345 2024-04-26 7/350 2024-04-30 16:33 by sang222
[硕博家园] 失之东隅 收之桑榆 +14 半世迷离 2024-04-27 25/1250 2024-04-30 16:28 by yinshan1006
[考博] 申请25年博士，求推荐学校及博导 +4 明开夜合20 2024-04-28 7/350 2024-04-30 14:23 by 明开夜合20
[教师之家] 一般发什么文章才能中青基啊 +6 保持肌肤 2024-04-29 6/300 2024-04-30 13:22 by ECOblance
[材料综合] 为什么石墨样氮化碳会拍出来这么诡异的TEM照片 20+3 Hikarico 2024-04-26 8/400 2024-04-30 12:55 by lxs7405
[找工作] 江苏理工学院要慎选 +8 jjchenshui 2024-04-27 10/500 2024-04-30 10:29 by 半生梦君
[分析] C18色谱柱可以分析单糖，二糖，低聚糖（3-7）吗 +5 备安网公 2024-04-26 6/300 2024-04-30 09:25 by shuiliu
[硕博家园] 博士后和科研助理好难找啊 +10 找工作啊 2024-04-28 10/500 2024-04-29 18:13 by 安塔瓦拉多
[考博] 博士招生 +6 zx179 2024-04-24 9/450 2024-04-29 17:53 by 263905569
[论文投稿] optics express投稿求助 +4 带带大周 2024-04-25 6/300 2024-04-29 16:52 by laozhao02
[考博] 25年博士申请 +9 Changzixuan 2024-04-25 17/850 2024-04-29 15:10 by Napoleonsky
[有机交流] 盐酸乐卡地平结晶问题 5+3 沉梦昂志光 2024-04-24 6/300 2024-04-29 14:34 by a229762843
[考研] 312求调剂 +4 平平无奇小Q 2024-04-26 5/250 2024-04-29 07:33 by 晓目崇
[论文投稿] with editor 两个月了，什么原因？ +9 yiersan9 2024-04-24 22/1100 2024-04-28 17:12 by ds521chun
[论文投稿] Jounal of Ethnopharmacology投稿 +3 lay木子禾 2024-04-25 4/200 2024-04-28 16:58 by hanyunch
[有机交流] 以2氟甲苯为原料做2氟4溴甲苯 +3 高树桩 2024-04-24 7/350 2024-04-28 14:37 by 高树桩
[考博] 材料方向24博士申请/一作SCI三篇 +3 白天不碰 2024-04-24 5/250 2024-04-27 17:41 by 安塔瓦拉多
[基金申请] 两类问题算是白选了~ +8 jurkat.1640 2024-04-23 13/650 2024-04-27 12:03 by 淀粉搬运工
[论文投稿] 投稿RSC旗下杂志突然间看不到投稿状态，里面投稿的文章信息不见了，有遇到过的吗？ 50+3 dlying 2024-04-26 3/150 2024-04-27 09:58 by bobvan
[考研] 381求调剂 +4 小刺猬987654321 2024-04-25 6/300 2024-04-26 10:57 by czl12138
信息提示
关闭
请填处理意见
关闭