24小时热门版块排行榜    

Znn3bq.jpeg
查看: 961  |  回复: 5

南方的落落

新虫 (小有名气)

[求助] 求助!!!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
回复此楼

» 猜你喜欢

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

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

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
模型公式不大合理吧,k1和q1,k2和q2,。。。k6和q6,完全可以各合并成一个参数即可,即减少了参数数,也防止参数过拟合而导致结果不唯一
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的回帖
相关版块跳转 我要订阅楼主 南方的落落 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 上海大学实验技术岗位非升即走 +4 嘻嘻哈哈乐呵呵 2026-05-15 4/200 2026-05-15 22:17 by waitnick
[论文投稿] 有带发论文的吗 +4 山楂之术 2026-05-09 4/200 2026-05-15 15:40 by 妹子不好惹
[文学芳草园] 风把牡丹吹跑了 +5 myrtle 2026-05-12 9/450 2026-05-15 15:27 by myrtle
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 l7k6xnh0yc 2026-05-14 3/150 2026-05-15 12:03 by onwj4wpxp2
[基金申请] 精华III评审感受-评审感受-评审感受 +14 ferrarichen 2026-05-11 18/900 2026-05-15 11:12 by cmhchen
[基金申请] 这年头没有找到涵评专家,还有中面上的可能吗 +9 dd921ww 2026-05-12 10/500 2026-05-15 10:41 by muyiliuhui
[考博] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 l7k6xnh0yc 2026-05-14 3/150 2026-05-15 09:23 by onwj4wpxp2
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 cjf4bx70cj 2026-05-14 4/200 2026-05-15 09:03 by gagyerk94e
[教师之家] 教学课件你会给同学吗 +8 硕士研究生吗 2026-05-13 8/400 2026-05-14 22:23 by 常规沥青
[基金申请] 重磅!青年科学基金项目(C类)资助增幅预计超过50% +5 水和泥不是水泥 2026-05-13 7/350 2026-05-14 20:57 by 水和泥不是水泥
[有机交流] 求助2,4-二氯-5-嘧啶甲醛的合成方法 20+3 光吃不拉 2026-05-14 5/250 2026-05-14 20:15 by 一切都是空工
[高分子] 本人最近太闲了,谁有问题可以提,每天会统一回复 +8 一切都是空工 2026-05-12 19/950 2026-05-14 20:03 by 一切都是空工
[考博] 申博自荐 +4 食品的橙子 2026-05-09 6/300 2026-05-14 16:05 by great1919
[考博] 材料类只有一篇综述能申博么 +4 乐逍遥谷 2026-05-13 4/200 2026-05-14 12:05 by zhyzzh
[基金申请] 请问大佬b0816评完了吗 +3 市民华南虎 2026-05-12 7/350 2026-05-14 07:41 by 市民华南虎
[硕博家园] 导师各种操作恶心咋办 +11 苍白的小青天 2026-05-09 13/650 2026-05-13 17:11 by 六两废铜
[论文投稿] 求助大佬sci投稿哪个好中 +3 江沅188 2026-05-12 4/200 2026-05-13 14:35 by 江沅188
[考博] 西南大学考核制博士 +3 lijunjie84 2026-05-11 6/300 2026-05-12 18:09 by lijunjie84
[文学芳草园] 窗边初夏的小雨 +7 阿美_Lml888 2026-05-09 10/500 2026-05-12 15:27 by 阿美_Lml888
[考博] 现在不知道怎么办,感觉很痛苦 +4 qweww 2026-05-11 5/250 2026-05-11 20:23 by Oversize
信息提示
请填处理意见