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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求助！！！1stopt进行常微分方程组参数模拟
汕头大学海洋科学接受调剂
61/1返回列表
查看: 925  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料类284调剂 +42 想换手机不想解� 2026-04-08 50/2500 2026-04-15 06:09 by 逍遥三郎
[考研] 调剂求收留 +33 果然有我 2026-04-10 34/1700 2026-04-15 00:37 by zzzggc
[考研] 一志愿鲁东大学071000生物学学硕初试分数276求调剂 +26 慕绝cc 2026-04-09 30/1500 2026-04-14 18:50 by 蔡苏阳
[考研] 071000生物学调剂求助 +18 zzzzwww 2026-04-09 21/1050 2026-04-14 15:39 by zs92450
[考研] 245求调剂 +6 冰糖橘?汽水 2026-04-13 10/500 2026-04-14 10:49 by jyl0317
[考研] 材料085601调剂 +32 何润采123 2026-04-10 34/1700 2026-04-14 08:47 by 木木mumu～
[考研] 考研求调剂 +12 子木呐 2026-04-12 13/650 2026-04-14 01:19 by 王珺璞
[考研] 339求调剂 +4 hanwudada 2026-04-12 4/200 2026-04-13 12:03 by 蓝云思雨
[考研] 322求调剂，08工科 +4 今天是个小号 2026-04-08 4/200 2026-04-13 00:20 by baobaoye
[考研] 求调剂 +16 张番茄不炒蛋 2026-04-10 17/850 2026-04-12 13:58 by 熬夜成！
[考研] 296求调剂 +14 汪！？！ 2026-04-08 15/750 2026-04-11 20:28 by dongdian1
[考研] 农学0904 312求调剂 +6 Say Never 2026-04-10 6/300 2026-04-11 10:33 by wwj2530616
[考研] 283求调剂，工科！ +12 苏打水7777 2026-04-08 12/600 2026-04-11 10:28 by 逆水乘风
[考研] 0854调剂 +8 950824he@ 2026-04-09 8/400 2026-04-11 10:11 by zhq0425
[考研] 263能源动力专硕求调剂 +3 加大号饭盒袋 2026-04-10 3/150 2026-04-10 22:23 by 286640313
[考研] 初试261 +3 Asht少 2026-04-10 6/300 2026-04-10 16:38 by Asht少
[考研] 085601初试330分找调剂 +10 流心奶黄包l 2026-04-09 10/500 2026-04-10 08:14 by Sammy2
[考研] 332，085601求调剂 +12 ydfyh 2026-04-09 14/700 2026-04-09 17:28 by wp06
[考研] 考研求调剂 +4 雯??? 2026-04-08 4/200 2026-04-08 21:44 by 土木硕士招生
[考研] 生物学328分求调剂 +9 闪电kkl 2026-04-08 10/500 2026-04-08 21:42 by liuhuiying09
信息提示
关闭
请填处理意见
关闭