| 查看: 977 | 回复: 1 | |||
wws1203金虫 (小有名气)
|
[求助]
多元微分方程组参数求解 已有1人参与
|
|
MATLAB程序求解 微分方程组如下： dx(1)/dt=x(8)*x(4)/x(11)-x(7)*x(1)/x(11)-x(1)*x(2)*k(1)*exp(0.098*(x(10)-15))/(k(2)*(x(1)+exp(0.051*(x(10)-1.58))))-k(3)*x(1)-0.18*mpower(k(4),x(10)-20)*x(1)*x(2)/(k(5)+x(1))-1/(1+mpower(10,10.068-0.033*x(10)-x(9)))+k(6)+x(2)*mpower(k(7),x(10)-25); dx(2)/dt=x(8)*x(5)/x(11)-x(7)*x(2)/x(11)+0.18*mpower(k(4),x(10)-20)*x(1)*x(2)/(x(1)+k(5))+0.18*mpower(k(4),x(10)-20)*x(3)*x(2)/(k(8)+x(3))+k(9)-x(2)*mpower(k(7),x(10)-25)-k(10)*x(2); dx(3)=x(8)*x(6)/x(11)-x(7)*x(3)/x(11)+x(1)*x(2)*k(1)*exp(0.098*(x(10)-15))/(k(2)*(x(1)+exp(0.051*(x(10)-1.58))))-x(3)*k(3)-0.18*mpower(k(4),x(10)-20)*x(3)*x(2)/(k(8)+x(3))-k(10)*x(3)*mpower(k(12),x(10)-20); 现有想x(1)---x(11)随时间变化的值，求方程组中的参数k(1)---k(10);求出参数出给定x(4)--x(11),反演出x(1)--x(3); t x(1) x(2) x(3) x(4) x(5) x(6) x(7) x(8) x(9) x(10) x(11) 0 3.289 7.852 0.535 0.904 14.926 0.282 0.019 0.050 7.860 31.600 0.690 2 2.909 9.055 0.518 0.460 61.763 0.315 0.029 0.072 7.710 32.000 0.720 4 3.441 4.177 0.518 0.807 10.849 0.352 0.035 0.094 7.690 33.200 0.700 7 2.563 7.692 0.502 0.330 15.841 0.315 0.035 0.026 7.780 30.900 0.690 9 2.324 5.497 0.435 3.712 4.145 0.348 0.026 0.037 8.040 29.300 0.670 11 2.248 7.585 0.390 1.956 5.412 0.431 0.020 0.020 8.020 29.300 0.668 13 2.216 8.663 0.452 0.644 8.672 0.410 0.012 0.017 8.150 29.100 0.674 |
回复此楼» 收录本帖的淘帖专辑推荐
 常微分方程 |
» 猜你喜欢
各位大佬，求一份最新Wien2k版本作为学习用途
已经有0人回复
-
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后
已经有553人回复
-
物理学I论文润色/翻译怎么收费?
已经有273人回复
-
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有130人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
2楼2018-03-25 21:10:27
