| 查看: 631 | 回复: 1 | ||
[求助]
求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 |
回复此楼» 本帖@通知
» 猜你喜欢
［调剂信息］211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有139人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有2人回复
-
070200求调剂，一志愿某211，288分
已经有18人回复
-
法国博士后职位
已经有0人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
广州大学光电信息工程专业调剂，招收物理学专业学生
已经有1人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有0人回复
hitcs
铜虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 45.9
- 帖子: 47
- 在线: 18.9小时
- 虫号: 3591713
- 注册: 2014-12-12
- 性别: GG
- 专业: 金属材料的合金相、相变及
2楼2018-04-28 20:20:57
