| 查看: 458 | 回复: 2 | |||
[求助]
求高版本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 |
回复此楼» 猜你喜欢
帮我的英语口语老师找学生
已经有4人回复
-
计算三阶力常数用thirdorder扩胞时参数如何选择？临近原子数怎么取值
已经有7人回复
-
物理学I论文润色/翻译怎么收费?
已经有200人回复
-
2026年循环经济功能材料国际会议（ICFMCE 2026）
已经有0人回复
-
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）
已经有0人回复
-
哈尔滨理工大学物理系招收物理学考研调剂
已经有5人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
独孤神宇: 金币+2, 鼓励交流 2020-07-28 12:33:03
hyj_swjtu: 金币+5, ★有帮助, 计算结果和理论值差距比较大。。。 2020-07-29 16:56:41
感谢参与，应助指数 +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
|
我添加添加了约束条件，希望有好心人帮忙代跑，谢谢 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
