24小时热门版块排行榜

>论坛更新日志 (3206)
>虫友互识 (604)
>导师招生 (221)
>文献求助 (168)
>考博 (165)
>休闲灌水 (160)
>论文投稿 (112)
>博后之家 (79)
>硕博家园 (76)
>考研 (71)
>教师之家 (43)
>招聘信息布告栏 (42)
>基金申请 (34)
>精细化工 (34)
>公派出国 (27)
>找工作 (25)

返回列表

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

回复此楼