【调剂】北京石油化工学院2024年16个专业接受调剂
查看: 5838  |  回复: 11
本帖产生 1 个 计算强帖 ,点击这里进行查看

yz457694

铁杆木虫 (正式写手)

[交流] 【好贴推荐】Origin共享参数拟合——Wilson汽液相平衡关联参数的拟合已有8人参与

这个方法只适用于自变量和参数相同的情况,可以参见之前解决过的一个问题以及后面(17楼)我给出的例子点击跳转。对于更一般的过程,基本上都可以用matlab编程,通过最小二乘的原理实现,可以参见另外一个例子点击跳转。下面开始此例。
实例说明:(见图1)
实验数据:(见图2)
自定义函数:在Origin中打开Tools--Fitting Function Organizer对话框,添加自定义函数(见图3)
自定义函数中function代码
CODE:
        double gamma1,gamma2,P,ps1,ps2;
        gamma1 = 1/(x+A12*(1-x))*exp((1-x)*(A12/(x+A12*(1-x))-A21/((1-x)+A21*x)));
        gamma2 = 1/((1-x)+A21*x)*exp(x*(A21/((1-x)+A21*x)-A12/(x+A12*(1-x))));
//补充说明一下Antoine方程形式根据手册说明而定,我这里用的就是lg
        ps1=10^((7.33827-(1650.5/(231.48+te))));
        ps2=10^((5.96407-(1200.31/(222.504+te))));
        P = 101.325;
        y1=ps1*x*gamma1/P;
        y2=ps2*(1-x)*gamma2/P;

也就是拟合两个函数y1、y2(对于多函数也可以类似试试),共享参数A12、A21。
关于函数中的常数,可以在Parameter Name里面事先定义,然后固定其值即可,或者直接在function代码中声明,然后使用,比如此处的P。过程中出现的变量也要声明,比如此处的ps1,ps2。
定义好函数之后,就可以在拟合对话框中选择使用自己定义好的函数了
拟合结果: (见图4)可见效果还不错啦
也可以参考一下origin的help文档,同样可以实现这个过程(见图5)
此外,对比用matlab实现,两者结果相同(见图6,自己弄了个gui,哈哈,origin用的LM算法,matlab非线性最小二乘默认的也是这个),所得相图见图7

写的匆忙,仅供交流
【好贴推荐】Origin共享参数拟合——Wilson汽液相平衡关联参数的拟合
图1-说明.jpg


【好贴推荐】Origin共享参数拟合——Wilson汽液相平衡关联参数的拟合-1
图2-data.jpg


【好贴推荐】Origin共享参数拟合——Wilson汽液相平衡关联参数的拟合-2
图3-自定义函数说明.jpg


【好贴推荐】Origin共享参数拟合——Wilson汽液相平衡关联参数的拟合-3
图4-拟合结果.jpg


【好贴推荐】Origin共享参数拟合——Wilson汽液相平衡关联参数的拟合-4
图6-matlab拟合结果.jpg


【好贴推荐】Origin共享参数拟合——Wilson汽液相平衡关联参数的拟合-5
图7-相图.jpg

[ Last edited by 月只蓝 on 2016-12-8 at 15:21 ]
回复此楼

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

yz457694

铁杆木虫 (正式写手)

图5看不清楚,传附件啦,相关文件也传上
[函数文件]Wilson_sharedmodel.FDF
[数据文件]wilson_VLE.opj

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:libolin3@tal.com
  • 附件 1 : origin共享参数拟合.png
  • 2016-10-25 00:16:07, 658.51 K
  • 附件 2 : Wilson_sharedmodel.FDF
  • 2016-10-25 00:17:32, 3.13 K
  • 附件 3 : wilson_VLE.opj
  • 2016-10-25 00:18:34, 198.81 K
2楼2016-10-25 00:18:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yz457694

铁杆木虫 (正式写手)

对应的matlab(不包含GUI)整理见附件
两种数据导入的方式,默认文本导入,把Excel中按顺序输好数据粘贴到文本文件即可。
main.m是主程序,Wilson.m是函数文件,importfile.m是文本数据导入函数
【好贴推荐】Origin共享参数拟合——Wilson汽液相平衡关联参数的拟合-6
Excel_data.jpg


【好贴推荐】Origin共享参数拟合——Wilson汽液相平衡关联参数的拟合-7
txt_data.jpg

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:libolin3@tal.com
  • 附件 1 : data_ethanol&cyclohexane.txt
  • 2017-06-12 08:58:26, 237 bytes
  • 附件 2 : data_ethanol&cyclohexane.xls
  • 2017-06-12 08:58:26, 28 K
  • 附件 3 : importfile.m
  • 2017-06-12 08:58:26, 425 bytes
  • 附件 4 : main.m
  • 2017-06-12 08:58:26, 1.62 K
  • 附件 5 : wilson.m
  • 2017-06-12 08:58:27, 607 bytes
10楼2017-06-12 09:02:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

cj129658

新虫 (初入文坛)

3楼2016-10-25 10:09:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
送红花一朵
4楼2016-10-26 09:55:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

旭日永恒

木虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
送红花一朵
感谢楼主,分享是一种美德,我要向你学习
我愿修禅,得证菩提。
5楼2016-12-08 14:02:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

瑞小仙儿

铁虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
送红花一朵
楼主方便也附上Matlab代码吗?
6楼2017-01-07 10:31:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zoe070721

银虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
腻害.......学习了,正在努力中
7楼2017-05-03 09:32:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小灵鱼儿

金虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
7楼: Originally posted by zoe070721 at 2017-05-03 09:32:38
腻害.......学习了,正在努力中

请问你有MATLAB代码吗,我QQ1005661393,急求呀
不打无准备之战
8楼2017-06-11 15:25:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yz457694

铁杆木虫 (正式写手)

引用回帖:
8楼: Originally posted by 小灵鱼儿 at 2017-06-11 15:25:52
请问你有MATLAB代码吗,我QQ1005661393,急求呀...

明天整理一下我再传上来吧

发自小木虫Android客户端
9楼2017-06-11 23:58:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yz457694 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 颜宁:基础研究应顶天立地,立地即从需求端作为问题的出发点,立项背景就是面向真实需求 +3 zju2000 2024-04-18 3/150 2024-04-20 00:48 by cococat951
[教师之家] 发现督导在听课,需要跟TA招呼示意吗?有同学看手机、课件有疏漏算大问题吗? +4 河西夜郎 2024-04-19 4/200 2024-04-20 00:17 by LittleBush
[教师之家] 国内科研环境差导致国内科研水平差的根源 +3 zju2000 2024-04-18 3/150 2024-04-19 21:46 by songwz
[论文投稿] 无under review,直接DIP,咋么办。 +5 lizhengke06 2024-04-19 7/350 2024-04-19 20:14 by otoo
[育儿交流] 小孩7岁,上一年级,不自信,学习较差 10+7 ljaxida 2024-04-16 34/1700 2024-04-19 18:49 by 水可南阳
[教师之家] 某系主任被同校教授举报包养情妇、长期嫖娼 +15 zju2000 2024-04-13 17/850 2024-04-19 18:17 by ZHONGWU_U
[找工作] 事业单位还是大学好? +17 青萍之沫 2024-04-16 18/900 2024-04-19 17:45 by charles-c
[考博] 山东大学化学与化工学院刘磊课题组博士研究生招生启事 +3 sdorganic 2024-04-17 5/250 2024-04-19 17:04 by 368ghnf
[找工作] 杭州国企和浙江高校如何选择? +15 restart2024 2024-04-15 21/1050 2024-04-19 16:05 by l419110028
[基金申请] 下雨了 +13 zju2000 2024-04-16 19/950 2024-04-19 09:24 by duxin_30
[有机交流] 怎么清洗烧瓶 20+5 ww34523 2024-04-16 6/300 2024-04-17 15:20 by 591950582
[基金申请] 请问教育部人文社科难度大吗 +9 锦衣卫寒战 2024-04-15 14/700 2024-04-16 23:10 by 锦衣卫寒战
[基金申请] 迟国泰通过向学生发放劳务费再回收的方式套取科学基金重点项目 +6 babu2015 2024-04-13 7/350 2024-04-16 20:32 by sundiv
[考研] 347求调剂 +3 寒辰ovo 2024-04-15 7/350 2024-04-16 19:05 by 寒辰ovo
[考博] 25年申博求助,射频方向 +4 SherlockAH 2024-04-15 10/500 2024-04-16 18:06 by ljysudaee
[考研] 294求调剂 +3 694062003 2024-04-15 4/200 2024-04-16 15:01 by 邹邹哈哈
[考研] 329求调剂 +18 王郁洁哈哈哈 2024-04-14 26/1300 2024-04-15 19:10 by mumin1990
[考研] 334求调剂 +4 学药救人 2024-04-14 4/200 2024-04-15 15:05 by hunanzang
[考研] 287求调剂 +6 南沨 2024-04-14 6/300 2024-04-14 23:08 by lincunhui
[考研] 309求调剂 +3 留不住放不下 2024-04-14 6/300 2024-04-14 15:28 by 19902169
信息提示
请填处理意见