| 查看: 623 | 回复: 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 |
回复此楼» 本帖@通知
» 猜你喜欢
2026年循环经济功能材料国际会议（ICFMCE 2026）
已经有0人回复
-
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有230人回复
-
哈尔滨理工大学物理系招收物理学考研调剂
已经有5人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
-
光学工程学硕调剂信息
已经有25人回复
-
欢迎加入课题组
已经有0人回复
-
散金币，求好运，祝面上顺利！
已经有29人回复
hitcs
铜虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 45.9
- 帖子: 47
- 在线: 18.9小时
- 虫号: 3591713
- 注册: 2014-12-12
- 性别: GG
- 专业: 金属材料的合金相、相变及
2楼2018-04-28 20:20:57
