24小时热门版块排行榜    

查看: 532  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-30 19:07 by xiangy672
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 4/200 2026-08-30 18:51 by hZiFeudZyoGR
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 gy1nBQXYQJqL 2026-08-29 5/250 2026-08-30 18:29 by hZiFeudZyoGR
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 18:17 by hZiFeudZyoGR
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +6 ASdOkHsho7FD 2026-08-28 10/500 2026-08-30 17:36 by hZiFeudZyoGR
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +7 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 17:34 by hZiFeudZyoGR
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +7 ASdOkHsho7FD 2026-08-28 14/700 2026-08-30 17:01 by hZiFeudZyoGR
[基金申请] 面上意见出来了 +9 黄鸟于飞Chao 2026-08-29 18/900 2026-08-30 16:47 by 黄鸟于飞Chao
[基金申请] 能否申诉? +4 echo8914667 2026-08-30 4/200 2026-08-30 16:03 by sundiv
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 gy1nBQXYQJqL 2026-08-29 8/400 2026-08-30 12:14 by l0VvVHGBGRLv
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +8 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 11:40 by l0VvVHGBGRLv
[基金申请] 基金不中,共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见