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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 其他 » 微分方程组求解，求助！！！
51/1返回列表
查看: 919  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zhl717

铜虫 (小有名气)

[求助] 微分方程组求解，求助！！！ 已有1人参与

哪位大神能用四阶龙格库特法帮忙解下微分方程组，感激不尽！！！


dx/dt=a*x(1-x/b)*(1+s/c)^(-1)
dp/dt=d*(dx/dt)+e*[s/(s+f)]*[(1+s/g)^(-1)]*x
-ds/dt=h*dx/dt+i*dp/dt+jx

t是自变量，x,p,s是因变量，a、b、c、d、e、f、g、h、i、j为参数
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
加油！！！
1楼 2015-07-25 10:24:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★
zhl717: 金币+1, ★有帮助 2015-08-20 09:51:33
引用回帖:
3楼: Originally posted by zhl717 at 2015-07-25 13:47:01
初值有，但是a-j参数没有数值，我最终的目就是拟合曲线求解参数，是不是没办法处理？...

可以用高版本的stopt软件来做，不过你的参数个数比较多，得用高版本的1stopt软件。
MATLAB可以参见：http://muchong.com/bbs/viewthread.php?tid=6425538
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
4楼2015-07-25 14:25:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

月只蓝

主管区长 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
龙格库塔法本身是用来求初值问题数值解的，楼主没给出初值，a~j参数的具体数值，以及要求解的t的范围，目前无法求解。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
2楼2015-07-25 11:14:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhl717

铜虫 (小有名气)
引用回帖:
2楼: Originally posted by 月只蓝 at 2015-07-25 11:14:32
龙格库塔法本身是用来求初值问题数值解的，楼主没给出初值，a~j参数的具体数值，以及要求解的t的范围，目前无法求解。

初值有，但是a-j参数没有数值，我最终的目就是拟合曲线求解参数，是不是没办法处理？
赞一下 回复此楼
加油！！！
3楼2015-07-25 13:47:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhl717

铜虫 (小有名气)
引用回帖:
4楼: Originally posted by 月只蓝 at 2015-07-25 14:25:06
可以用高版本的stopt软件来做，不过你的参数个数比较多，得用高版本的1stopt软件。
MATLAB可以参见：http://muchong.com/bbs/viewthread.php?tid=6425538...

大神，你能帮我写一下吗，我对matlab不了解，感激不尽
赞一下 回复此楼
加油！！！
5楼2015-07-25 14:52:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +4 hKe7U1kad5Sg 2026-10-31 9/450 2026-11-02 03:24 by xzcv0Tehd8lQ
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 xvJS1BYDxvBS 2026-10-31 7/350 2026-11-02 03:18 by xzcv0Tehd8lQ
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 xvJS1BYDxvBS 2026-10-31 9/450 2026-11-02 03:11 by xzcv0Tehd8lQ
[考研] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-02 02:58 by xzcv0Tehd8lQ
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-02 02:47 by xzcv0Tehd8lQ
[硕博家园] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-02 00:47 by xzcv0Tehd8lQ
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 hKe7U1kad5Sg 2026-10-31 8/400 2026-11-02 00:37 by xzcv0Tehd8lQ
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 22:24 by xzcv0Tehd8lQ
[基金申请] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 22:11 by xzcv0Tehd8lQ
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 20:03 by xzcv0Tehd8lQ
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 16:48 by BCdvBAEpt3zm
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +5 X5bHxOULxFqb 2026-10-30 22/1100 2026-11-01 16:35 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 xvJS1BYDxvBS 2026-10-31 5/250 2026-11-01 12:51 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 12:38 by BCdvBAEpt3zm
[考博] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 10/500 2026-11-01 12:27 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 1uaZt2kAZj2S 2026-10-30 17/850 2026-11-01 12:24 by BCdvBAEpt3zm
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-01 12:14 by BCdvBAEpt3zm
[硕博家园] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 12:14 by BCdvBAEpt3zm
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 4/200 2026-11-01 11:01 by BCdvBAEpt3zm
[硕博家园] 申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请 +3 小小阿光 2026-10-27 3/150 2026-10-29 23:56 by 起风了的快乐
信息提示
关闭
请填处理意见
关闭