| 查看: 713 | 回复: 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 |
回复此楼» 本帖@通知
» 猜你喜欢
青基已中
已经有2人回复
-
散金币祈福
已经有186人回复
-
物理学I论文润色/翻译怎么收费?
已经有241人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有109人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有64人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有4人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有257人回复
hitcs
铜虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 45.9
- 帖子: 47
- 在线: 18.9小时
- 虫号: 3591713
- 注册: 2014-12-12
- 性别: GG
- 专业: 金属材料的合金相、相变及
2楼2018-04-28 20:20:57
10