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

登录注册

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

24小时热门版块排行榜

返回列表

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

回复此楼

» 猜你喜欢

瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向已经有130人回复
PRB投稿，Acknowledgment sent to author状态十天了已经有3人回复
物理学I论文润色/翻译怎么收费? 已经有286人回复
PVA溶解已经有2人回复
基金申请已经有44人回复
CSC与新西兰维多利亚大学PhD奖学金项目已经有0人回复
新西兰Robinson研究所招聘CSC公派访问人员已经有0人回复
帮我的英语口语老师找学生已经有0人回复
什么时候开奖？已经有8人回复
散金币祈福已经有94人回复

» 本主题相关价值贴推荐，对您同样有帮助:

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.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

已阅回复此楼关注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的回帖

julescesar

铁杆木虫 (知名作家)

傻子中的牛虻

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

引用回帖:

k2的数值太小了点儿，感觉不大对

赞一下

回复此楼

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

7楼2014-04-08 16:15:22

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

julescesar

铁杆木虫 (知名作家)

傻子中的牛虻

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

引用回帖:

请版主指点一二，小虫外行

赞一下

回复此楼

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

8楼2014-04-08 21:55:34

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

相关版块跳转我要订阅楼主 julescesar 的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3	2JOx3r2CYEgw 2026-08-22	4/200	2026-08-23 03:55 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3	2JOx3r2CYEgw 2026-08-22	7/350	2026-08-23 03:55 by OEbVnUOu01ol
[论文投稿] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3	2JOx3r2CYEgw 2026-08-22	5/250	2026-08-23 03:52 by OEbVnUOu01ol
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3	2JOx3r2CYEgw 2026-08-21	9/450	2026-08-23 03:31 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3	QTy3jDtz1uLt 2026-08-21	9/450	2026-08-23 03:04 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4	QTy3jDtz1uLt 2026-08-21	11/550	2026-08-23 02:55 by OEbVnUOu01ol
[基金申请] 什么时候开奖？ +7	CrisMessi 2026-08-18	8/400	2026-08-22 18:29 by 淀粉搬运工
[基金申请] 人气不行了 +8	fansofjerry 2026-08-21	8/400	2026-08-22 16:30 by zyqchem
[基金申请] 时间戳今天，20号变了 +5	archvillain 2026-08-20	5/250	2026-08-22 06:12 by hui_daxiao
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +16	医学老男孩 2026-08-20	20/1000	2026-08-21 21:13 by Ldrop2023
[基金申请] 科研孤儿太难了 +17	我4大白菜 2026-08-20	18/900	2026-08-21 20:57 by zhangev
[基金申请] 看来今天不会放榜了？ +8	chengyan1220 2026-08-21	11/550	2026-08-21 17:52 by dcqxinyang
[基金申请] 让我中一个面上吧！ +11	大萍1987 2026-08-20	13/650	2026-08-21 14:48 by 20120902066
[基金申请] 应该是下周三26日公布了吧？ +4	哈哈蛤？ 2026-08-21	4/200	2026-08-21 10:58 by Vivilian
[论文投稿] 投稿咨询 +5	wwm09 2026-08-17	7/350	2026-08-21 10:11 by 期刊论文帮手
[基金申请] 基金啊基金 +4	longfie172 2026-08-20	4/200	2026-08-21 08:58 by mark mao
[基金申请] 估计是周四 +3	archvillain 2026-08-18	3/150	2026-08-21 01:48 by jnhyjjm
[基金申请] 重要消息，中午系统在维护 +11	yuleib84 2026-08-18	12/600	2026-08-20 11:09 by xskun
[基金申请] 今天放榜吗？ +14	布布和一二 2026-08-19	15/750	2026-08-19 18:07 by gltch
[基金申请] 朋友圈看到的 +6	wangzilk 2026-08-18	8/400	2026-08-19 10:55 by Haru815

24小时热门版块排行榜

[求助] 求高版本1stopt跑一个程序 已有1人参与

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

【答案】应助回帖

[求助] 求高版本1stopt跑一个程序已有1人参与