版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

南方科技大学公共卫生及应急管理学院2025级博士研究生招生报考通知

返回列表

julescesar

铁杆木虫 (知名作家)

傻子中的牛虻

应助: 85 (初中生)
金币: 7405.9
散金: 89
红花: 5
帖子: 5225
在线: 552.7小时
虫号: 548890
注册: 2008-04-19
性别: GG
专业: 环境地球化学

[求助] 求高版本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

回复此楼

» 猜你喜欢

第一性原理计算方向2026级博士申请 PRB*1，四级484 已经有1人回复
求助VISSIM破解版软件已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有232人回复
求2026年在台湾举行的物理和材料领域国际学术会议信息已经有0人回复
求国际会议网站已经有1人回复
求取一些关于纳米材料和纳米技术相关的英文PPT。已经有0人回复
【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有14人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有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人回复

一个罪犯也有权检举另一个罪犯，以事实为采信依据；反之同理。神州法治move-forward！

1楼2014-04-08 09:03:48

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

julescesar

铁杆木虫 (知名作家)

傻子中的牛虻

应助: 85 (初中生)
金币: 7405.9
散金: 89
红花: 5
帖子: 5225
在线: 552.7小时
虫号: 548890
注册: 2008-04-19
性别: GG
专业: 环境地球化学

请有资源的朋友帮个忙，谢谢

赞一下

回复此楼

一个罪犯也有权检举另一个罪犯，以事实为采信依据；反之同理。神州法治move-forward！

2楼2014-04-08 11:21:39

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

julescesar

铁杆木虫 (知名作家)

傻子中的牛虻

应助: 85 (初中生)
金币: 7405.9
散金: 89
红花: 5
帖子: 5225
在线: 552.7小时
虫号: 548890
注册: 2008-04-19
性别: GG
专业: 环境地球化学

求脚印。。。。

回复此楼

一个罪犯也有权检举另一个罪犯，以事实为采信依据；反之同理。神州法治move-forward！

3楼2014-04-08 12:24:12

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

julescesar

铁杆木虫 (知名作家)

4楼2014-04-08 15:09:13

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

计算强帖: 4
应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.5小时
虫号: 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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

julescesar

铁杆木虫 (知名作家)

傻子中的牛虻

应助: 85 (初中生)
金币: 7405.9
散金: 89
红花: 5
帖子: 5225
在线: 552.7小时
虫号: 548890
注册: 2008-04-19
性别: GG
专业: 环境地球化学

引用回帖:

5楼: Originally posted by dingd at 2014-04-08 15:30:28
将就参考下：

参数最佳估算
-------------------- -------------
k1 2.0920685752888
k2 2.03145877517177E-9
k3 0.248632668721487
k4 1.35390776847532E-17
k5 6.03374135447111E-9
k6 3.61 ...

谢谢版主。。。
能不能把你修正过的程序代码和计算后的实测值与模拟值比较部分的数据贴出来？谢谢啦

赞一下

回复此楼

一个罪犯也有权检举另一个罪犯，以事实为采信依据；反之同理。神州法治move-forward！

6楼2014-04-08 15:52:30

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖