24小时热门版块排行榜    

查看: 455  |  回复: 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 的主题更新
信息提示
请填处理意见