| 查看: 1195 | 回复: 7 | ||
julescesar铁杆木虫 (知名作家)
傻子中的牛虻
|
[求助]
求高版本1stopt跑一个程序 已有1人参与
|
|
程序如下，如果不合语法，请高手帮忙修正后再跑，谢谢啦 Title "lump6"; Parameter k1,k2,k3,k4,k5,k6,k7,k8,k9,k10,k11,k12,k13,k14; ParameterDomain=[0,]; Variable t,y(6),z(6),p(6),q(6),r(6); ODEFunction y1'=-(k1+k2+k3+k4+k5)*y1/(y1/500+y2/300+y3/150+y4/50+y5/20+y6/400)/(1+50*k6)/(1+0.5*k7); y2'=(k1*y1-(k8+k9+k10+k11)*y2)/(y1/500+y2/300+y3/150+y4/50+y5/20+y6/400)/(1+50*k6)/(1+0.5*k7); y3'=(k2*y1+k8*y2-(k12+k13+k14)*y3)/(y1/500+y2/300+y3/150+y4/50+y5/20+y6/400)/(1+50*k6)/(1+0.5*k7); y4'=(k3*y1+k9*y2+k12*y3)/(y1/500+y2/300+y3/150+y4/50+y5/20+y6/400)/(1+50*k6)/(1+0.5*k7); y5'=(k4*y1+k10*y2+k13*y3)/(y1/500+y2/300+y3/150+y4/50+y5/20+y6/400)/(1+50*k6)/(1+0.5*k7); y6=100-y1-y2-y3-y4-y5; z1'=-(k1+k2+k3+k4+k5)*z1/(z1/500+z2/300+z3/150+z4/50+z5/20+z6/400)/(1+50*k6)/(1+0.4*k7); z2'=(k1*z1-(k8+k9+k10+k11)*z2)/(z1/500+z2/300+z3/150+z4/50+z5/20+z6/400)/(1+50*k6)/(1+0.4*k7); z3'=(k2*z1+k8*z2-(k12+k13+k14)*z3)/(z1/500+z2/300+z3/150+z4/50+z5/20+z6/400)/(1+50*k6)/(1+0.4*k7); z4'=(k3*z1+k9*z2+k12*z3)/(z1/500+z2/300+z3/150+z4/50+z5/20+z6/400)/(1+50*k6)/(1+0.4*k7); z5'=(k4*z1+k10*z2+k13*z3)/(z1/500+z2/300+z3/150+z4/50+z5/20+z6/400)/(1+50*k6)/(1+0.4*k7); z6=100-z1-z2-z3-z4-z5; p1'=-(k1+k2+k3+k4+k5)*p1/(p1/500+p2/300+p3/150+p4/50+p5/20+p6/400)/(1+50*k6)/(1+0.3*k7); p2'=(k1*p1-(k8+k9+k10+k11)*p2)/(p1/500+p2/300+p3/150+p4/50+p5/20+p6/400)/(1+50*k6)/(1+0.3*k7); p3'=(k2*p1+k8*p2-(k12+k13+k14)*p3)/(p1/500+p2/300+p3/150+p4/50+p5/20+p6/400)/(1+50*k6)/(1+0.3*k7); p4'=(k3*p1+k9*p2+k12*p3)/(p1/500+p2/300+p3/150+p4/50+p5/20+p6/400)/(1+50*k6)/(1+0.3*k7); p5'=(k4*p1+k10*p2+k13*p3)/(p1/500+p2/300+p3/150+p4/50+p5/20+p6/400)/(1+50*k6)/(1+0.3*k7); p6=100-p1-p2-p3-p4-p5; q1'=-(k1+k2+k3+k4+k5)*q1/(q1/500+q2/300+q3/150+q4/50+q5/20+q6/400)/(1+50*k6)/(1+0.2*k7); q2'=(k1*q1-(k8+k9+k10+k11)*q2)/(q1/500+q2/300+q3/150+q4/50+q5/20+q6/400)/(1+50*k6)/(1+0.2*k7); q3'=(k2*q1+k8*q2-(k12+k13+k14)*q3)/(q1/500+q2/300+q3/150+q4/50+q5/20+q6/400)/(1+50*k6)/(1+0.2*k7); q4'=(k3*q1+k9*q2+k12*q3)/(q1/500+q2/300+q3/150+q4/50+q5/20+q6/400)/(1+50*k6)/(1+0.2*k7); q5'=(k4*q1+k10*q2+k13*q3)/(q1/500+q2/300+q3/150+q4/50+q5/20+q6/400)/(1+50*k6)/(1+0.2*k7); q6=100-q1-q2-q3-q4-q5; r1'=-(k1+k2+k3+k4+k5)*r1/(r1/500+r2/300+r3/150+r4/50+r5/20+r6/400)/(1+50*k6)/(1+0.1*k7); r2'=(k1*r1-(k8+k9+k10+k11)*r2)/(r1/500+r2/300+r3/150+r4/50+r5/20+r6/400)/(1+50*k6)/(1+0.1*k7); r3'=(k2*r1+k8*r2-(k12+k13+k14)*r3)/(r1/500+r2/300+r3/150+r4/50+r5/20+r6/400)/(1+50*k6)/(1+0.1*k7); r4'=(k3*r1+k9*r2+k12*r3)/(r1/500+r2/300+r3/150+r4/50+r5/20+r6/400)/(1+50*k6)/(1+0.1*k7); r5'=(k4*r1+k10*r2+k13*r3)/(r1/500+r2/300+r3/150+r4/50+r5/20+r6/400)/(1+50*k6)/(1+0.1*k7); r6=100-r1-r2-r3-r4-r5; Data; t y1 y2 y3 y4 y5 y6 z1 z2 z3 z4 z5 z6 p1 p2 p3 p4 p5 p6 q1 q2 q3 q4 q5 q6 r1 r2 r3 r4 r5 r6 0 93.7 6.3 0 0 0 0 93.7 6.3 0 0 0 0 93.7 6.3 0 0 0 0 93.7 6.3 0 0 0 0 93.7 6.3 0 0 0 0 1 5.34 16 47 20 2 9.66 4.95 15 48 19 3 10.05 4.42 13 48 20 3 11.58 3.79 13 47 21 3.3 11.91 3.58 12 47 22 3.3 12.12 |
回复此楼» 猜你喜欢
帮我的英语口语老师找学生
已经有4人回复
-
计算三阶力常数用thirdorder扩胞时参数如何选择？临近原子数怎么取值
已经有7人回复
-
物理学I论文润色/翻译怎么收费?
已经有150人回复
-
2026年循环经济功能材料国际会议（ICFMCE 2026）
已经有0人回复
-
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）
已经有0人回复
-
哈尔滨理工大学物理系招收物理学考研调剂
已经有5人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 1stopt解方程组
已经有12人回复
- 用1stopt软件做微分方程组的参数估计
已经有13人回复
- 求用高版本的1stOpt拟合方程
已经有8人回复
- 有1stopt的虫友给跑下程序，求！
已经有5人回复
- 我的MS低版本打不开高版本，怎么办啊，急求助！！！
已经有9人回复
- 谁用1stopt这个软件啊，离散数据拟合方程，方程未知
已经有5人回复
- 求大神帮忙用1stOpt 拟合复数数据
已经有13人回复
- rmse，sse计算求助
已经有3人回复
- 求高版本1stopt拟合，
已经有11人回复
- 1stopt中使用sharedmodel拟合曲线出现问题，求指示！！
已经有6人回复
- 微分方程组参数拟合的问题（多参数动力系统参数估计）
已经有11人回复
- 求大哥大姐帮跑一下1stopt~新手啊，版本低啊
已经有9人回复
- 求大侠帮助跑一下1stopt
已经有12人回复
- 谁有1stopt，帮忙跑一下啊
已经有6人回复
- 哪位能用1stopt的高版本帮我走个模型啊！谢谢了！我的1.5的不能跑两模型共参数啊
已经有3人回复
- 哪位大哥能帮我用1stopt3.0以上版本跑一个拟合程序
已经有11人回复
julescesar
铁杆木虫 (知名作家)
傻子中的牛虻
- 应助: 85 (初中生)
- 金币: 7405.9
- 散金: 89
- 红花: 5
- 帖子: 5225
- 在线: 552.7小时
- 虫号: 548890
- 注册: 2008-04-19
- 性别: GG
- 专业: 环境地球化学
2楼2014-04-08 11:21:39
julescesar
铁杆木虫 (知名作家)
傻子中的牛虻
- 应助: 85 (初中生)
- 金币: 7405.9
- 散金: 89
- 红花: 5
- 帖子: 5225
- 在线: 552.7小时
- 虫号: 548890
- 注册: 2008-04-19
- 性别: GG
- 专业: 环境地球化学
3楼2014-04-08 12:24:12
julescesar
铁杆木虫 (知名作家)
 |
4楼2014-04-08 15:09:13
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
感谢参与，应助指数 +1
|
将就参考下： 参数 最佳估算 -------------------- ------------- k1 2.0920685752888 k2 2.03145877517177E-9 k3 0.248632668721487 k4 1.35390776847532E-17 k5 6.03374135447111E-9 k6 3.61924495624649E-9 k7 0.689252026420965 k8 2.16983628672973 k9 2.98161186429235E-10 k10 4.39628051135967E-8 k11 1.57735517991293E-9 k12 0.278468868858522 k13 0.0748917320017696 k14 0.283316627567864 |
5楼2014-04-08 15:30:28
julescesar
铁杆木虫 (知名作家)
傻子中的牛虻
- 应助: 85 (初中生)
- 金币: 7405.9
- 散金: 89
- 红花: 5
- 帖子: 5225
- 在线: 552.7小时
- 虫号: 548890
- 注册: 2008-04-19
- 性别: GG
- 专业: 环境地球化学
6楼2014-04-08 15:52:30
julescesar
铁杆木虫 (知名作家)
傻子中的牛虻
- 应助: 85 (初中生)
- 金币: 7405.9
- 散金: 89
- 红花: 5
- 帖子: 5225
- 在线: 552.7小时
- 虫号: 548890
- 注册: 2008-04-19
- 性别: GG
- 专业: 环境地球化学
7楼2014-04-08 16:15:22
julescesar
铁杆木虫 (知名作家)
傻子中的牛虻
- 应助: 85 (初中生)
- 金币: 7405.9
- 散金: 89
- 红花: 5
- 帖子: 5225
- 在线: 552.7小时
- 虫号: 548890
- 注册: 2008-04-19
- 性别: GG
- 专业: 环境地球化学
8楼2014-04-08 21:55:34
20