24小时热门版块排行榜    

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

249038396

新虫 (初入文坛)

[求助] 用1stOpt进行多模型多变量非线性拟合得不到结果 已有1人参与

本人用网上下载的1.5版本的1stOpt进行多模型多变量非线性拟合,但程序部运行。论坛里有人说可能是版本不支持share model的运算。先附上自己写得代码,希望有高版本的虫友帮我看下是代码的问题还是软件版本的问题,谢谢!
Parameters  k1,k2,kco2,ko2,e1,e2,eco2,eo2,a;
Variable    T,P,Y1,YCO2,Y3,YO2,YH2O,Rk1,Rk2;
ConstStr         P3=P*Y3;
ConstStr         P1=P*Y1;
ConstStr         PCO2=P*YCO2 ;
ConstStr         PO2=P*YO2 ;
ConstStr         PH2O=P*YH2O ;
ConstStr         Rg=8.315;
SharedModel;
Function  Rk1=(k1*EXP(-e1/Rg/T)* P3*PO2)/(1+kco2*EXP(-eco2/Rg/T)*PCO2+ko2*EXP(-eo2/Rg/T)*(PO2^0.5)* PH2O);
          Rk2=(k2*EXP(-e2/Rg/T)* P3*(PO2^a))/(1+kco2*EXP(-eco2/Rg/T)*PCO2+ko2*EXP(-eo2/Rg/T)*(PO2^0.5)* PH2O);
Data;
484.9417         2.1163         0.0051         0.0269         0.1650         0.0753         0.0020         2.1457         0.0049
485.0713         2.0952         0.0031         0.0508         0.2872         0.0470         0.0012         1.2896         0.0029
485.0828         2.0876         0.0062         0.0076         0.3143         0.0694         0.0024         2.5848         0.0058
484.9200         2.0992         0.0033         0.0063         0.3170         0.0371         0.0012         1.3700         0.0028
485.0941         2.0947         0.0067         0.0527         0.3440         0.0788         0.0028         2.7487         0.0067
488.1700         2.0973         0.0029         0.0559         0.2270         0.0371         0.0013         1.1490         0.0029
487.9231         2.0977         0.0051         0.0221         0.3454         0.0297         0.0022         1.9418         0.0049
488.1433         2.0967         0.0046         0.0218         0.2561         0.0504         0.0019         1.7768         0.0043
487.9061         2.0960         0.0056         0.0575         0.2841         0.0638         0.0026         2.1888         0.0058
488.5491         2.0973         0.0083         0.0395         0.3130         0.0560         0.0041         3.1633         0.0091
回复此楼

» 猜你喜欢

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

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

zyugiec

金虫 (正式写手)

请问楼主解决了此问题吗?好像是盗版的1stopt无法调用sharemodel功能。
5楼2016-08-17 10:27:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
249038396: 金币+5, ★有帮助, 非常感谢! 2014-01-02 21:34:22
你的Rk2数据与Rk1数据差了将近3个数量级,拟合中RK1几乎要被忽略,建议先调整下。
2楼2013-12-31 15:14:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

249038396

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by dingd at 2013-12-31 15:14:50
你的Rk2数据与Rk1数据差了将近3个数量级,拟合中RK1几乎要被忽略,建议先调整下。

我想知道在不考虑数据合理性的前提下,我的这段程序有问题么?
3楼2014-01-02 09:39:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

引用回帖:
3楼: Originally posted by 249038396 at 2014-01-02 09:39:38
我想知道在不考虑数据合理性的前提下,我的这段程序有问题么?...

没问题吧。
4楼2014-01-02 11:07:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-29 19:27 by cqwQDCxMcL3I
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:54 by yCO1Ll7aHtsw
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:47 by yCO1Ll7aHtsw
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:30 by etmYJ6d2rquH
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:03 by etmYJ6d2rquH
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:04 by etmYJ6d2rquH
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 14:58 by etmYJ6d2rquH
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-29 14:24 by etmYJ6d2rquH
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 IDs3scOF0tjC 2026-09-28 7/350 2026-09-29 14:07 by etmYJ6d2rquH
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 IDs3scOF0tjC 2026-09-28 3/150 2026-09-29 13:56 by etmYJ6d2rquH
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 09:59 by JzYBbHIXWSrW
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 GDBe8tDZqE8z 2026-09-28 4/200 2026-09-29 09:49 by JzYBbHIXWSrW
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:42 by ZvyPK8n6Nfki
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:28 by ZvyPK8n6Nfki
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-28 23:46 by b8fCuTHqckEt
[有机交流] 同一个分子,一条来自文献,一条来自AI——不告诉你答案,你会选哪条? +3 tianxiaoxian 2026-09-28 6/300 2026-09-28 22:44 by tianxiaoxian
[找工作] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 5/250 2026-09-28 22:32 by ez6fGg9abYaj
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 6/300 2026-09-28 22:15 by ez6fGg9abYaj
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-28 19:08 by 9lS3ad5oOymn
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 18:51 by 9lS3ad5oOymn
信息提示
请填处理意见