24小时热门版块排行榜

>论坛更新日志 (2981)
>导师招生 (880)
>虫友互识 (575)
>休闲灌水 (423)
>论文投稿 (153)
>招聘信息布告栏 (149)
>文献求助 (125)
>博后之家 (95)
>学术会议 (62)
>硕博家园 (50)
>考博 (43)
>教师之家 (42)
>考研 (40)
>公派出国 (30)
>基金申请 (24)
>找工作 (19)

返回列表

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

回复此楼