24小时热门版块排行榜

>论坛更新日志 (4504)
>虫友互识 (528)
>文献求助 (248)
>导师招生 (196)
>基金申请 (187)
>休闲灌水 (130)
>博后之家 (110)
>考博 (109)
>找工作 (97)
>招聘信息布告栏 (87)
>论文道贺祈福 (86)
>论文投稿 (83)
>硕博家园 (78)
>公派出国 (78)
>考研 (60)
>教师之家 (54)

返回列表

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

回复此楼