| 查看: 722 | 回复: 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 |
回复此楼» 本帖@通知
» 猜你喜欢
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
物理学I论文润色/翻译怎么收费?
已经有205人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有0人回复
hitcs
铜虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 45.9
- 帖子: 47
- 在线: 18.9小时
- 虫号: 3591713
- 注册: 2014-12-12
- 性别: GG
- 专业: 金属材料的合金相、相变及
2楼2018-04-28 20:20:57
50