| 查看: 704 | 回复: 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 |
回复此楼» 本帖@通知
» 猜你喜欢
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有130人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
物理学I论文润色/翻译怎么收费?
已经有214人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
hitcs
铜虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 45.9
- 帖子: 47
- 在线: 18.9小时
- 虫号: 3591713
- 注册: 2014-12-12
- 性别: GG
- 专业: 金属材料的合金相、相变及
2楼2018-04-28 20:20:57
