| 查看: 409 | 回复: 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 |
回复此楼» 猜你喜欢
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有126人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
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