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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求助！！！1stopt进行常微分方程组参数模拟
61/1返回列表
查看: 1023  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +14 爱看书的可乐 2026-08-26 16/800 2026-08-26 17:22 by skyjudy
[基金申请] 基金不中，共勉 +3 eulota 2026-08-26 3/150 2026-08-26 16:17 by 路yyyy
[基金申请] 我不理解！ +8 Edward_pc 2026-08-26 13/650 2026-08-26 15:37 by Tide man
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +6 process2012 2026-08-23 9/450 2026-08-26 15:33 by process2012
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 只有每年这种时候来逛逛小木虫 +26 yaoyewhu2008 2026-08-20 28/1400 2026-08-26 13:49 by zzuzxg
[基金申请] 科研孤儿太难了 +21 我4大白菜 2026-08-20 22/1100 2026-08-26 11:29 by nlgza
[基金申请] 怎么看青基中了没有啊 +4 叶九微 2026-08-26 4/200 2026-08-26 10:52 by xiacongxu
[基金申请] 国合现在查不到了吗？ +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 牛来！米来！面来！ +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 让我中一个面上吧！ +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 时间戳今天，20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
关闭
请填处理意见
关闭