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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求助！！！1stopt进行常微分方程组参数模拟
汕头大学海洋科学接受调剂
61/1返回列表
查看: 924  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

南方的落落

新虫 (小有名气)

[求助] 求助！！！1stopt进行常微分方程组参数模拟 已有2人参与

求助大神!!!
我在做动力学模型，列出了微分方程组，用1stopt编了如下的程序，但是运行不出来
Parameter k1,k2,k3,k4,k5,k6,q1,q2,q3,q4,q5,q6;
Variable t,y1,y2,y3,y4,y5,y6;
ODEFunction
y1'=-k1*q1*y1;
y2'=k1*q1*y1-(k2*q2+k3*q3+k4*q4+k5*q5)*y2;
y3'=k2*q2*y2;
y4'=k3*q3*y2;
y5'=k4*q4*y2-k6*q6*y5;
y6'=k5*q5*y2+k6*q6*y5;
Data;
0.083333333        93.97083255        7.840492639        9.572094691        3.518630812        9.196398765        9.58599691
0.166666667        90.09133592        18.73857849        9.029740675        1.280160211        6.49678316        7.242624229
0.333333333        86.8961721        29.5306655        8.07535529        0.441393587        3.535430781        5.146967947
0.5        84.31720661        38.82233506        5.671183361        0.279535385        2.541740915        3.738556994
0.75        81.76266667        39.56        3.2        0.25        1.7        2.714213576
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2016-01-20 15:39:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
模型公式不大合理吧，k1和q1，k2和q2，。。。k6和q6，完全可以各合并成一个参数即可，即减少了参数数，也防止参数过拟合而导致结果不唯一
赞一下(1人) 回复此楼
2楼2016-01-20 22:27:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
2楼 dingd 大神说得很对，k1和q1，k2和q2，…，k6和q6，都是可以合并成单个参数的。如果不合并，拟合结果会有无穷多种。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
3楼2016-01-21 09:29:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

南方的落落

新虫 (小有名气)
引用回帖:
2楼: Originally posted by dingd at 2016-01-20 22:27:12
模型公式不大合理吧，k1和q1，k2和q2，。。。k6和q6，完全可以各合并成一个参数即可，即减少了参数数，也防止参数过拟合而导致结果不唯一

谢谢，我把k1和q1合并了，我说的运行不出来是软件没进行运算，我这个是低版本1.5的，不知道是不是软件问题
赞一下 回复此楼
4楼2016-01-21 20:42:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

南方的落落

新虫 (小有名气)
引用回帖:
3楼: Originally posted by 月只蓝 at 2016-01-21 09:29:24
2楼 dingd 大神说得很对，k1和q1，k2和q2，…，k6和q6，都是可以合并成单个参数的。如果不合并，拟合结果会有无穷多种。

谢谢，合并了之后还是一样，软件没进行计算
赞一下 回复此楼
5楼2016-01-21 20:43:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
引用回帖:
5楼: Originally posted by 南方的落落 at 2016-01-21 20:43:52
谢谢，合并了之后还是一样，软件没进行计算...

http://muchong.com/bbs/viewthread.php?tid=9176191

参见该链接中第4条解释。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
6楼2016-01-21 22:03:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 南方的落落 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料工程281还有调剂机会吗 +41 xaw. 2026-04-11 42/2100 2026-04-14 20:05 by honglizhao
[考研] 一志愿211，化学310分，本科重点双非，求调剂 +24 努力奋斗112 2026-04-08 24/1200 2026-04-14 18:02 by lhj2009
[考研] 302分求调剂 +10 凡语祈愿 2026-04-08 11/550 2026-04-14 16:50 by jiangguiquan11
[考研] 材料专业344求调剂 +17 hualkop 2026-04-10 22/1100 2026-04-14 16:21 by sxdj2
[考研] 一志愿西南大学生物学学硕344 求生物学相关调剂/生物与医药 +8 超人不会飞@ 2026-04-08 8/400 2026-04-14 16:06 by zs92450
[考研] 290求调剂 +21 luoziheng 2026-04-10 23/1150 2026-04-14 15:49 by zs92450
[考研] 279求调剂 +12 张番茄不炒蛋 2026-04-11 12/600 2026-04-14 15:38 by zs92450
[考研] 085408光电信息工程专硕355一志愿长春光机所调剂 +6 王ymaa 2026-04-13 13/650 2026-04-14 11:33 by 王ymaa
[考研] 一志愿华中农微生物，288分，三年实验经历 +11 代fish 2026-04-09 11/550 2026-04-12 10:21 by Hayaay
[考研] 343求调剂 +9 王国帅 2026-04-10 9/450 2026-04-11 20:31 by dongdian1
[考研] 085410-273求调剂 +6 X1999 2026-04-10 6/300 2026-04-11 10:32 by Delta2012
[考研] 283求调剂，工科！ +12 苏打水7777 2026-04-08 12/600 2026-04-11 10:28 by 逆水乘风
[考研] 机械专硕270求调剂，接受跨专业 +12 老师看看我吧aba 2026-04-09 14/700 2026-04-11 10:21 by laoshidan
[考研] 284求调剂 +9 让我上岸吧阿西 2026-04-09 11/550 2026-04-10 19:18 by 靖jing
[考研] 初试261 +3 Asht少 2026-04-10 6/300 2026-04-10 16:38 by Asht少
[考研] 085800 能源动力求调剂 +6 阿biu啊啊啊啊啊 2026-04-10 6/300 2026-04-10 15:03 by hemengdong
[考研] 材料专硕(0856) 339分求调剂 +9 哈哈哈鹅哈哈哈 2026-04-09 10/500 2026-04-09 20:01 by Orcid
[论文投稿] 求助文献原文 10+3 18500821399 2026-04-08 3/150 2026-04-09 16:56 by 北京莱茵润色
[考研] 085501机械英二77总分294求调剂，接受跨专业学习 +6 守法公民亓纪 2026-04-08 6/300 2026-04-09 15:55 by wp06
[考研] 调剂 +12 月@163.com 2026-04-08 12/600 2026-04-09 14:27 by rl1980
信息提示
关闭
请填处理意见
关闭