1stopt微分方程组参数拟合(各位大哥求帮忙,没有高版本的软件)
初始金币 5
Parameter a11=0.01[0,10000];
Parameter a12=0.01[0,10000];
Parameter a13=0.01[0,10000];
Parameter b1=0.01[0,10000];
Parameter a21=0.01[0,10000];
Parameter a22=0.01[0,10000];
Parameter a23=0.01[0,10000];
Parameter b2=0.01[0,10000];
Parameter a31=0.01[0,10000];
Parameter a32=0.01[0,10000];
Parameter a33=0.01[0,10000];
Parameter b3=0.01[0,10000];
Variable t,x,y;
ODEFunction x'=a11* x+a12*y+a13*z-b1*x*x;
y'= a21* x+a22*y+a23*z-b2*y*y;
z'= a31* x+a32*y+a33*z-b2*y*y;
Data;
1 0 0 0
33 15 12 6
43 39 30 30
72 84 60 63
92 174 102 119
102 334 194 188
109 414 258 284
123 608 452 404 @beefly 返回小木虫查看更多
今日热帖
参考下:
均方差(RMSE): 22.4162459568179
残差平方和(SSE): 7034.83315915174
相关系数(R): 0.99124454515698
相关系数之平方(R^2): 0.982565748303469
修正R平方(Adj. R^2): 0.965161343617583
确定系数(DC): 0.981591902258889
F统计(F-Statistic): -48.1592299042318
参数 最佳估算
-------------------- -------------
a11 0
a12 0
a13 2149.17512435827
b1 182.108564892615
a21 8194.43647500787
a22 708.554507187899
a23 8991.36975189023
b2 1719.88087625676
a31 4602.03187540038
a32 4117.88086212192
a33 9440.06926253831
b3 2765.01241145004
z 初值 0.000129132966324601,
非常感谢!!参数a11为0有点问题啊。
能否再算一下Parameter??a11=0.01
非常感谢楼主,不过参数a11理论上不能是0,能否再重新算一下,谢谢
Parameter a11=0.01[0,10000];
Parameter a12=0.01[0,10000];
Parameter b1=0.01[0,10000];
Parameter a21=0.01[0,10000];
Parameter a22=0.01[0,10000];
Parameter b2=0.01[0,10000];
Variable t,x,y;
ODEFunction x'=a11* x+a12*y-b1*x*x);
y'= a21* x+a22*y-b2*y*y);
Data;
1 0 0
33 12 6
43 30 30
72 60 63
92 102 119
102 194 188
109 258 284
123 452 404 @beefly