24小时热门版块排行榜    

查看: 1563  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

sh073

木虫 (正式写手)

[求助] 求高版本1stOpT跑代码 已有2人参与

Title "Phase equilibrium stage2";
Parameters K1[0,],K2[0,],n[1,10,0],m[0,10,0],y[0,1];
Variable caq,corg,cn235,chol,D;
ShareModel;
Function D=(n-1)*K1*(cn235-y*(corg-cn235)/(n-1))*caq^(n-2)+cn235/caq+m*K2*(chol-(1-y)*(corg-cn235)/m)*caq^(m-1);
         D=(n-1)*K1*(cn235-y*(corg-cn235)/(n-1))*caq^(n-2)/y+cn235/caq;

Data;
0.2598        0.4132        0.4107        3.654770278        1.590454196
0.5372        0.4606        0.4107        3.654770278        0.857408786
0.7053        0.4742        0.4107        3.654770278        0.672338012
0.9979        0.515        0.4107        3.654770278        0.516083776
1.585        0.5818        0.4107        3.654770278        0.367066246
2.083        0.618        0.4107        3.654770278        0.29668747
2.628        0.666        0.4107        3.654770278        0.253424658
3.036        0.7122        0.4107        3.654770278        0.23458498
3.4        0.7447        0.4107        3.654770278        0.219029412
1.81        0.4605        0.308025        3.654770278        0.25441989
1.761        0.5421        0.3593625        3.654770278        0.307836457
1.704        0.6005        0.4107        3.654770278        0.352406103
1.62        0.6531        0.4620375        3.654770278        0.403148148
1.574        0.7144        0.513375        3.654770278        0.453875476
1.707        0.5721        0.4107        2.741077708        0.335149385
1.693        0.5836        0.4107        3.197923993        0.344713526
1.691        0.5889        0.4107        3.654770278        0.34825547
1.702        0.6025        0.4107        4.111616563        0.3539953
1.673        0.6143        0.4107        4.568462847        0.367184698
3.991        0.8286        0.4107        3.654770278        0.207617139
4.675        0.929        0.4107        3.654770278        0.198716578
5.226        1.081        0.4107        3.654770278        0.206850364
6.1732        1.324        0.4107        3.654770278        0.214475475
6.614        1.437        0.4107        3.654770278        0.217266405
7.184        1.618        0.4107        3.654770278        0.225222717
8.707        2.103        0.4107        3.654770278        0.241529804
6.683        1.227        0.308025        3.654770278        0.18360018
6.667        1.328        0.3593625        3.654770278        0.19919004
6.683        1.43        0.4107        3.654770278        0.213975759
6.532        1.49        0.4620375        3.654770278        0.228107777
6.586        1.568        0.513375        3.654770278        0.238080777
6.653        1.303        0.4107        2.741077708        0.195851496
6.622        1.35        0.4107        3.197923993        0.203865902
6.563        1.426        0.4107        3.654770278        0.217278684
6.511        1.522        0.4107        4.111616563        0.233758255
6.545        1.62        0.4107        4.568462847        0.247517189
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
sh073: 金币+30, ★★★★★最佳答案, 谢谢了 2014-11-04 21:04:19
修改后代码如下:
CODE:
Title "Phase equilibrium stage2";
Parameters K1=[0,],K2=[0,],n=[1,10,0],m=[0,10,0],y=[0,1];
Variable caq,corg,cn235,chol,D,DD;
SharedModel;
Function D=(n-1)*K1*(cn235-y*(corg-cn235)/(n-1))*caq^(n-2)+cn235/caq+m*K2*(chol-(1-y)*(corg-cn235)/m)*caq^(m-1);
         DD=(n-1)*K1*(cn235-y*(corg-cn235)/(n-1))*caq^(n-2)/y+cn235/caq;

Data;
0.2598        0.4132        0.4107        3.654770278        1.590454196        1.590454196
0.5372        0.4606        0.4107        3.654770278        0.857408786        0.857408786
0.7053        0.4742        0.4107        3.654770278        0.672338012        0.672338012
0.9979        0.515        0.4107        3.654770278        0.516083776        0.516083776
1.585        0.5818        0.4107        3.654770278        0.367066246        0.367066246
2.083        0.618        0.4107        3.654770278        0.29668747        0.29668747
2.628        0.666        0.4107        3.654770278        0.253424658        0.253424658
3.036        0.7122        0.4107        3.654770278        0.23458498        0.23458498
3.4        0.7447        0.4107        3.654770278        0.219029412        0.219029412
1.81        0.4605        0.308025        3.654770278        0.25441989        0.25441989
1.761        0.5421        0.3593625        3.654770278        0.307836457        0.307836457
1.704        0.6005        0.4107        3.654770278        0.352406103        0.352406103
1.62        0.6531        0.4620375        3.654770278        0.403148148        0.403148148
1.574        0.7144        0.513375        3.654770278        0.453875476        0.453875476
1.707        0.5721        0.4107        2.741077708        0.335149385        0.335149385
1.693        0.5836        0.4107        3.197923993        0.344713526        0.344713526
1.691        0.5889        0.4107        3.654770278        0.34825547        0.34825547
1.702        0.6025        0.4107        4.111616563        0.3539953        0.3539953
1.673        0.6143        0.4107        4.568462847        0.367184698        0.367184698
3.991        0.8286        0.4107        3.654770278        0.207617139        0.207617139
4.675        0.929        0.4107        3.654770278        0.198716578        0.198716578
5.226        1.081        0.4107        3.654770278        0.206850364        0.206850364
6.1732        1.324        0.4107        3.654770278        0.214475475        0.214475475
6.614        1.437        0.4107        3.654770278        0.217266405        0.217266405
7.184        1.618        0.4107        3.654770278        0.225222717        0.225222717
8.707        2.103        0.4107        3.654770278        0.241529804        0.241529804
6.683        1.227        0.308025        3.654770278        0.18360018        0.18360018
6.667        1.328        0.3593625        3.654770278        0.19919004        0.19919004
6.683        1.43        0.4107        3.654770278        0.213975759        0.213975759
6.532        1.49        0.4620375        3.654770278        0.228107777        0.228107777
6.586        1.568        0.513375        3.654770278        0.238080777        0.238080777
6.653        1.303        0.4107        2.741077708        0.195851496        0.195851496
6.622        1.35        0.4107        3.197923993        0.203865902        0.203865902
6.563        1.426        0.4107        3.654770278        0.217278684        0.217278684
6.511        1.522        0.4107        4.111616563        0.233758255        0.233758255
6.545        1.62        0.4107        4.568462847        0.247517189        0.247517189

均方差(RMSE):0.0350538182185307
残差平方和(SSE):0.0884714523622415
相关系数(R): 0.992402867842004
相关系数之平方(R^2): 0.984863452101034
确定系数(DC): 0.980778115771102
F统计(F-Statistic): 589.665203190385

参数                  最佳估算
--------------------        -------------
k1        0.0162982778480216
k2        0.0260910853201875
n        3
m        1
y        0.319189098340218
9楼2014-11-04 20:13:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

月只蓝

主管区长 (职业作家)

2楼2014-11-03 14:15:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

水晶杯喝茶

新虫 (正式写手)

高版本好像很贵,这个低版本能跑出来,想不通为啥要用高版本,能解释一下吗?
3楼2014-11-03 15:34:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

引用回帖:
3楼: Originally posted by 水晶杯喝茶 at 2014-11-03 15:34:08
高版本好像很贵,这个低版本能跑出来,想不通为啥要用高版本,能解释一下吗?

其一,低版本ShareModel 用不了。
其二,高版本拟合结果更好、更稳定。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
4楼2014-11-03 15:47:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-30 03:26 by ZPa0EcMwuECS
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-30 02:29 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-30 02:18 by ZPa0EcMwuECS
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 3/150 2026-08-30 01:58 by ZPa0EcMwuECS
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:58 by ZPa0EcMwuECS
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:36 by ZPa0EcMwuECS
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 01:04 by ZPa0EcMwuECS
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 ASdOkHsho7FD 2026-08-28 6/300 2026-08-30 01:03 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 12/600 2026-08-30 00:31 by ZPa0EcMwuECS
[基金申请] 有没有仍没收到信息的 +5 德尚中行 2026-08-27 6/300 2026-08-29 18:10 by manplx
[基金申请] 为什么到现在没收到通知? +4 tannykie 2026-08-29 4/200 2026-08-29 17:14 by lmz0216
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
信息提示
请填处理意见