24小时热门版块排行榜    

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

hyj_swjtu

新虫 (初入文坛)

[求助] 求高版本1stopt代跑程序 已有1人参与

求各位大佬能帮我用高版本1stopt拟合一下4参数微分方程
Parameter r1,r2,c1,c2;
Variable x,p,y;
ODEFunction y''=(r2*p+298.15-y-(r1*c1+r2*c1+r2*c2)*y')/(r1*r2*c1*c2);
Data;
0                                  111.381838113994 298.150000000000
0.0100000000000000 111.503353400036 298.150000000000
0.0200000000000000 111.624868686079 298.151827753426
0.0300000000000000 111.746383972121 298.155412534919
0.0400000000000000 111.867899258163 298.160686241319
0.0500000000000000 111.989414544206 298.167583294755
0.0600000000000000 112.110929830248 298.176040549312
0.0700000000000000 112.232445116290 298.185997201158
0.0800000000000000 112.353960402332 298.197394701979
0.0900000000000000 112.475475688375 298.210176675616
0.100000000000000 112.596990974417 298.224288837791
0.110000000000000 112.679178726144 298.239678918785
0.120000000000000 112.761366477871 298.256296588981
0.130000000000000 112.843554229598 298.274092741802
0.140000000000000 112.925741981325 298.293020095686
0.150000000000000 113.007929733053 298.313033126625
0.160000000000000 113.090117484780 298.334088003196
0.170000000000000 113.172305236507 298.356142523997
0.180000000000000 113.254492988234 298.379156057391
0.190000000000000 113.336680739961 298.403089483476
0.200000000000000 113.418868491688 298.427905138199
0.210000000000000 113.267199148875 298.453566759541
0.220000000000000 113.115529806062 298.480039435681
0.230000000000000 112.963860463250 298.507285717534
0.240000000000000 112.812191120437 298.535269560909
0.250000000000000 112.660521777624 298.563956274539
0.260000000000000 112.508852434811 298.593312470037
0.270000000000000 112.357183091998 298.623306013696
0.280000000000000 112.205513749186 298.653905980069
0.290000000000000 112.053844406373 298.685082607273
0.300000000000000 111.902175063560 298.716807253934
0.310000000000000 111.738586966646 298.749052357731
0.320000000000000 111.574998869732 298.781791395464
0.330000000000000 111.411410772817 298.814998649022
0.340000000000000 111.247822675903 298.848649371716
0.350000000000000 111.084234578989 298.882719752328
0.360000000000000 110.920646482075 298.917186880485
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hyj_swjtu

新虫 (初入文坛)

我添加添加了约束条件,希望有好心人帮忙代跑,谢谢
Parameter r1,r2,c1,c2;
Variable x,p,y;
ODEFunction y\'\'=(r2*p+298.15-y-(r1*c1+r2*c1+r2*c2)*y\')/(r1*r2*c1*c2);
r1>0;
r2>0;
c1>0;
c2>0;
r1<1;
r2<1;
c1<1;
Data;
0                                  111.381838113994 298.150000000000
0.0100000000000000 111.503353400036 298.150000000000
0.0200000000000000 111.624868686079 298.151827753426
0.0300000000000000 111.746383972121 298.155412534919
0.0400000000000000 111.867899258163 298.160686241319
0.0500000000000000 111.989414544206 298.167583294755
0.0600000000000000 112.110929830248 298.176040549312
0.0700000000000000 112.232445116290 298.185997201158
0.0800000000000000 112.353960402332 298.197394701979
0.0900000000000000 112.475475688375 298.210176675616
0.100000000000000 112.596990974417 298.224288837791
0.110000000000000 112.679178726144 298.239678918785
0.120000000000000 112.761366477871 298.256296588981
0.130000000000000 112.843554229598 298.274092741802
0.140000000000000 112.925741981325 298.293020095686
0.150000000000000 113.007929733053 298.313033126625
0.160000000000000 113.090117484780 298.334088003196
0.170000000000000 113.172305236507 298.356142523997
0.180000000000000 113.254492988234 298.379156057391
0.190000000000000 113.336680739961 298.403089483476
0.200000000000000 113.418868491688 298.427905138199
0.210000000000000 113.267199148875 298.453566759541
0.220000000000000 113.115529806062 298.480039435681
0.230000000000000 112.963860463250 298.507285717534
0.240000000000000 112.812191120437 298.535269560909
0.250000000000000 112.660521777624 298.563956274539
0.260000000000000 112.508852434811 298.593312470037
0.270000000000000 112.357183091998 298.623306013696
0.280000000000000 112.205513749186 298.653905980069
0.290000000000000 112.053844406373 298.685082607273
0.300000000000000 111.902175063560 298.716807253934
0.310000000000000 111.738586966646 298.749052357731
0.320000000000000 111.574998869732 298.781791395464
0.330000000000000 111.411410772817 298.814998649022
0.340000000000000 111.247822675903 298.848649371716
0.350000000000000 111.084234578989 298.882719752328
0.360000000000000 110.920646482075 298.917186880485
3楼2020-07-31 15:56:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
独孤神宇: 金币+2, 鼓励交流 2020-07-28 12:33:03
hyj_swjtu: 金币+5, ★有帮助, 计算结果和理论值差距比较大。。。 2020-07-29 16:56:41
多解吧:
1:
均方差(RMSE): 6.80722680798251E-6
残差平方和(SSR): 1.66818012535137E-9
相关系数(R): 0.999999999598921
相关系数之平方(R^2): 0.999999999197842
修正R平方(Adj. R^2): 0.999999999149226
确定系数(DC): 0.999999999193842
F统计(F-Statistic): 9613467604.93016

参数                  最佳估算
--------------------        -------------
r1        -0.033947317590326
r2        0.182708353590288
c1        -6.044507004901
c2        29.1575066853031
y' 初值         -0.093603801762556

2:
均方差(RMSE): 6.80722680798251E-6
残差平方和(SSR): 1.66818012535137E-9
相关系数(R): 0.999999999598921
相关系数之平方(R^2): 0.999999999197842
修正R平方(Adj. R^2): 0.999999999149226
确定系数(DC): 0.999999999193842
F统计(F-Statistic): 9613467604.93016

参数                  最佳估算
--------------------        -------------
r1        -0.033947317590326
r2        0.182708353590288
c1        -6.044507004901
c2        29.1575066853031
y' 初值         -0.093603801762556

。。。。
2楼2020-07-28 10:02:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 6/300 2026-10-31 23:23 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 xvJS1BYDxvBS 2026-10-31 7/350 2026-10-31 23:18 by YXM1FDCSE8Bn
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 5/250 2026-10-31 23:10 by YXM1FDCSE8Bn
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 12/600 2026-10-31 22:52 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 X5bHxOULxFqb 2026-10-30 14/700 2026-10-31 22:46 by YXM1FDCSE8Bn
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 X5bHxOULxFqb 2026-10-30 19/950 2026-10-31 22:44 by YXM1FDCSE8Bn
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 X5bHxOULxFqb 2026-10-30 24/1200 2026-10-31 22:39 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 X5bHxOULxFqb 2026-10-30 11/550 2026-10-31 22:39 by YXM1FDCSE8Bn
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 5/250 2026-10-31 21:59 by YXM1FDCSE8Bn
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 4/200 2026-10-31 21:51 by YXM1FDCSE8Bn
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 4/200 2026-10-31 21:51 by YXM1FDCSE8Bn
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 10/500 2026-10-31 21:41 by YXM1FDCSE8Bn
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 1uaZt2kAZj2S 2026-10-30 11/550 2026-10-31 21:25 by YXM1FDCSE8Bn
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 X5bHxOULxFqb 2026-10-30 18/900 2026-10-31 17:05 by hKe7U1kad5Sg
[考博] 2027博士申请-临床检验诊断学、病原生物学、微生物方向 50+3 yyh@ 2026-10-28 4/200 2026-10-31 12:47 by 可淡不可忘
[硕博家园] 越来越担心硕士毕业后没有竞争力了 +8 zuesli 2026-10-27 14/700 2026-10-30 22:44 by 好运连连公主
[基金申请] 论文的挂名怎么评判贡献?基金申请中,排名第一的挂名文章 +3 biily 2026-10-25 3/150 2026-10-30 18:33 by 旁观2020
[硕博家园] 申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请 +3 小小阿光 2026-10-27 3/150 2026-10-29 23:56 by 起风了的快乐
[文学芳草园] 闲聊 +3 myrtle 2026-10-25 6/300 2026-10-29 15:28 by myrtle
[考博] 申博求助 20+3 Katrinae 2026-10-26 4/200 2026-10-27 09:56 by 北京莱茵编辑
信息提示
请填处理意见