24小时热门版块排行榜

返回列表

dana1997$

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 6
帖子: 1
在线: 3.7小时
虫号: 8233083
注册: 2018-03-11
专业: 生态经济学

[求助] 求1stopt高版本代跑微分方程参数估计

最近在研究lotka—volterra方程，找了很多办法解微分方程，发现还是1stopt最好。
代码的格式不一定正确，求帮忙改下代跑程序，感激不尽！！
Parameter a(1:7),b(1:8),c(1:8),d(1:8),e(1:8),f(1:8),g(1:8),h(1:8);
Variable A,B,C,D,E,F,G;
ODEFunction A' = A*(a1+b1*A+b2*B+b3*C+b4*D+b5*E+b6*F+b7*G)+b8;
B'= B*(a2+c1*A+c2*B+c3*C+c4*D+c5*E+c6*F+c7*G)+c8;
C'= C*(a3+d1*A+d2*B+d3*C+d4*D+d5*E+d6*F+d7*G)+d8;
D'= D*(a4+e1*A+e2*B+e3*C+e4*D+e5*E+e6*F+e7*G)+e8;
E'= E*(a5+f1*A+f2*B+f3*C+f4*D+f5*E+f6*F+f7*G)+f8;
F'= F*(a6+g1*A+g2*B+g3*C+g4*D+g5*E+g6*F+g7*G)+g8;
G'= G*(a7+h1*A+h2*B+h3*C+h4*D+h5*E+h6*F+h7*G)+h8;
Data;
//A B C D E F G
0.99062    3.325031321 1.20064    1.505528    0.465544    1.0456    0.465544
1.14162    3.34957033 0.72806    1.02674    0.5293          0.58934    0.5293
1.35414    3.371823121 1.08226    0.924514    0.629788    0.979    0.629788
1.54394    3.358702072 1.05092    1.221278    0.756664    1.5998    0.756664
1.75564    3.386495364 1.12072    1.26153    0.911154    1.7214    0.911154
2.14156    3.41386638 1.0099    1.033504    1.149582    0.9758    1.40958
2.51394    3.441276984 0.97434    1.08504    1.46752    2.0614    1.468836
2.75924    3.469425525 3.23728    0.951092    1.826014    0.87182 1.499694
3.30608    3.532804179 3.35756    1.53954    2.24755    2.1928    1.478168
3.94312    3.554224579 3.482    0.839842    2.933274    1.307    1.514688
4.34662    3.57388128 3.47258    1.519284    3.390792    1.3512    1.531788
4.78092    3.595427135 3.56552    1.11271    3.804548    2.4814    1.517656
5.18538    3.614506304 3.55198    1.285118    4.134148    0.9752    1.542748
5.54818    3.63956211 3.60134    1.42194    4.916166    1.2714    1.543972
6.09814    3.669662833 3.56924    1.69098    5.241606    2.1204    1.532076

回复此楼