24小时热门版块排行榜    

查看: 551  |  回复: 2

hitsz.civil

金虫 (正式写手)

[求助] 求助下列微分方程如何拟合出系数y"=(v-b*y'-c*y-d*y^3-e*y^5)/a;具体见内容 已有1人参与

各位高手,请问如何拟合出y"=(v-b*y'-c*y-d*y^3-e*y^5)/a;中的a,b,c,d,e五个系数,已知t,v,y,非常感谢!下面是我尝试的计算,但没有结果,请各位高手z指教,谢谢
Parameters a,b,c,d,e;
Variable t,v,y;
ODEFunction y"=(v-b*y'-c*y-d*y^3-e*y^5)/a;
Data;
0.004         62.85         9.89
0.048         59.98         9.66
0.075         55.99         9.50
0.095         52.00         9.12
0.113         47.71         8.89
0.127         43.99         8.51
0.140         40.00         8.13
0.153         35.98         7.74
0.166         31.71         7.13
0.177         28.00         6.99
0.188         24.01         6.52
0.198         19.98         5.99
0.209         15.99         5.53
0.219         12.00         5.07
0.230         7.98         4.62
0.239         4.30         4.08
0.250         0.00         3.17
0.256         -2.28         2.63
0.260         -3.99         1.86
0.265         -6.00         0.88
0.270         -8.01         0.03
0.281         -12.00         -0.04
0.292         -16.30         -0.20
0.296         -18.00         -0.95
0.298         -18.58         -1.80
0.304         -20.87         -1.34
0.308         -22.57         -2.17
0.313         -24.28         -2.56
0.319         -26.29         -3.93
0.323         -27.72         -4.93
0.327         -29.15         -5.92
0.330         -30.28         -6.69
0.335         -31.99         -7.37
0.348         -36.28         -7.75
0.360         -40.00         -8.28
0.373         -43.99         -8.59
0.388         -48.01         -8.89
0.406         -52.28         -9.20
0.425         -55.99         -9.50
0.452         -59.98         -9.74
0.500         -63.43         -9.97
0.548         -59.98         -9.74
0.575         -55.99         -9.50
0.594         -52.28         -9.20
0.612         -48.01         -8.89
0.627         -43.99         -8.59
0.640         -40.00         -8.28
0.652         -36.28         -7.75
0.665         -31.99         -7.37
0.676         -28.27         -6.99
0.687         -24.28         -6.38
0.698         -20.02         -5.84
0.709         -15.99         -5.46
0.719         -12.00         -4.70
0.730         -8.01         -4.16
0.740         -3.99         -3.86
0.750         0.00         -3.09
0.756         2.28         -2.17
0.760         3.99         -1.03
0.767         6.85         0.11
0.775         9.99         0.11
0.786         13.98         0.42
0.791         15.99         0.88
0.796         18.00         1.79
0.799         18.86         2.25
0.802         20.29         1.72
0.807         22.00         2.40
0.812         24.01         2.48
0.818         25.99         3.55
0.824         28.00         4.69
0.828         29.70         6.14
0.834         31.71         7.13
0.847         35.98         7.74
0.860         40.00         8.13
0.873         43.99         8.51
0.887         47.71         8.89
0.905         52.00         9.12
0.925         55.99         9.50
0.996         62.85         9.89
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
参考:

均方差(RMSE): 0.451508938019568
残差平方和(SSE): 15.9011050467015
相关系数(R): 0.997535757851112
相关系数之平方(R^2): 0.995077588191592
修正R平方(Adj. R^2): 0.9949463238767
确定系数(DC): 0.995034595750716
F统计(F-Statistic): 2886.39879977383

参数                  最佳估算
--------------------        -------------
a        0.000647685287294533
b        0.105079723338619
c        2.60999578795763
d        0.0238747667505231
e        0.000159801268513268
y' 初值         -730.05176280801
求助下列微分方程如何拟合出系数y"=(v-b*y'-c*y-d*y^3-e*y^5)/a;具体见内容
c6.jpg

2楼2016-10-25 14:16:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hitsz.civil

金虫 (正式写手)

引用回帖:
2楼: Originally posted by dingd at 2016-10-25 14:16:29
参考:

均方差(RMSE): 0.451508938019568
残差平方和(SSE): 15.9011050467015
相关系数(R): 0.997535757851112
相关系数之平方(R^2): 0.995077588191592
修正R平方(Adj. R^2): 0.9949463238767
确定系数(DC ...

非常感谢dingd!
3楼2016-10-25 15:12:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hitsz.civil 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:42 by ZvyPK8n6Nfki
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:34 by ZvyPK8n6Nfki
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:28 by ZvyPK8n6Nfki
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:18 by ZvyPK8n6Nfki
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-28 23:46 by b8fCuTHqckEt
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 GDBe8tDZqE8z 2026-09-28 4/200 2026-09-28 23:19 by b8fCuTHqckEt
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-28 23:12 by b8fCuTHqckEt
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 5/250 2026-09-28 22:45 by ez6fGg9abYaj
[有机交流] 同一个分子,一条来自文献,一条来自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文章T0P,我:8O.551.O54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 22:28 by ez6fGg9abYaj
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 3/150 2026-09-28 22:21 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
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 17:30 by ZYXdhzDAy9ZX
[教师之家] 某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意! +7 zju2000 2026-09-26 8/400 2026-09-28 16:42 by beefly
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-28 15:32 by ZYXdhzDAy9ZX
信息提示
请填处理意见