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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求助！！！1stopt进行常微分方程组参数模拟
61/1返回列表
查看: 902  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085602 289分求调剂 +3 WWW西西弗斯 2026-03-24 3/150 2026-03-24 14:26 by JourneyLucky
[考研] 材料292调剂 +6 橘颂思美人 2026-03-23 6/300 2026-03-24 13:32 by 连洋洋洋
[考研] 求调剂 +5 研研，接电话 2026-03-24 6/300 2026-03-24 11:39 by 星空星月
[考研] 一志愿南航材料专317分求调剂 +4 炸呀炸呀炸薯条 2026-03-23 4/200 2026-03-24 07:32 by wangy0907
[考研] 求调剂 +7 十三加油 2026-03-21 7/350 2026-03-23 23:48 by 热情沙漠
[考研] 一志愿211 初试270分 求调剂 +4 谷雨上岸 2026-03-23 5/250 2026-03-23 21:18 by 不惑可乐
[考研] 336化工调剂 +4 王大坦1 2026-03-23 5/250 2026-03-23 18:32 by allen-yin
[考研] 招08考数学 +6 laoshidan 2026-03-20 14/700 2026-03-23 14:37 by 15614359529
[考研] 315分，诚求调剂，材料与化工085600 +3 13756423260 2026-03-22 3/150 2026-03-22 20:11 by edmund7
[考研] 一志愿华中农业071010，总分320求调剂 +5 困困困困坤坤 2026-03-20 6/300 2026-03-22 17:41 by hxsm
[考研] 298求调剂一志愿211 +3 上岸6666@ 2026-03-20 3/150 2026-03-22 15:50 by ColorlessPI
[考研] 311求调剂 +3 勇敢的小吴 2026-03-20 3/150 2026-03-21 17:40 by ColorlessPI
[考研] 296求调剂 +4 www_q 2026-03-20 4/200 2026-03-21 17:26 by 学员8dgXkO
[考研] 299求调剂 +6 △小透明* 2026-03-17 6/300 2026-03-21 02:42 by JourneyLucky
[考研] 330求调剂 +4 小材化本科 2026-03-18 4/200 2026-03-20 23:13 by JourneyLucky
[考研] 一志愿武汉理工材料工程专硕调剂 +9 Doleres 2026-03-19 9/450 2026-03-20 22:36 by JourneyLucky
[考研] 295复试调剂 +8 简木ChuFront 2026-03-19 8/400 2026-03-20 20:44 by zhukairuo
[考研] 材料学求调剂 +4 Stella_Yao 2026-03-20 4/200 2026-03-20 20:28 by ms629
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
[考研] 288求调剂，一志愿华南理工大学071005 +5 ioodiiij 2026-03-17 5/250 2026-03-19 18:22 by zcl123
信息提示
关闭
请填处理意见
关闭