24小时热门版块排行榜    

查看: 582  |  回复: 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的回帖

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的回帖

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的回帖
相关版块跳转 我要订阅楼主 hyj_swjtu 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 2026国自然函评费到账 +9 羊腰板 2026-08-21 9/450 2026-08-22 00:03 by anji999
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 让我中一个面上吧! +11 大萍1987 2026-08-20 13/650 2026-08-21 14:48 by 20120902066
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 什么时候开奖? +7 CrisMessi 2026-08-18 7/350 2026-08-21 12:15 by weiyin
[基金申请] 时间戳又变了8-15 +15 archvillain 2026-08-15 27/1350 2026-08-21 11:32 by tim76
[基金申请] 今天基金会出结果吗?20260819 +15 kkkl_v 2026-08-19 16/800 2026-08-21 11:22 by 365372687
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
信息提示
请填处理意见