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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 会用matlab的同志帮帮忙哈！
281/3返回列表123下一页
查看: 1385  |  回复: 27
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

phoenix708

铁杆木虫 (著名写手)

[求助] 会用matlab的同志帮帮忙哈！

dX/dt=1.98S(1-X)-0.73XY
dY/dt=dX/dt-0.00207Y/(1+Y+2S)-0.059d[Y/(1+Y+2S)]/dt
2dS/dt=-dX/dt-0.00056S/(1+Y+2S)-0.086d[S/(1+Y+2S)]/dt
当t=0时,X=0,Y=0,S=1,
其中 t的最大值为3，
求X，Y, S与t的关系图
谢谢!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-01-05 11:47:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

scumat

金虫 (正式写手)
问题很难，20个金币少了。
赞一下 回复此楼
2楼2012-01-05 16:40:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

phoenix708

铁杆木虫 (著名写手)
引用回帖:
楼: Originally posted by scumat at 2012-01-05 16:40:17:
问题很难，20个金币少了。

你能解吗？能解就可以加金币
赞一下 回复此楼
3楼2012-01-05 16:46:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

phoenix708

铁杆木虫 (著名写手)
我把方程组化简了一下，不知对否？懂得虫友帮忙运算下，应该不难！x'=1.98S*(1-X)-0.73*X*Y;                       
y'=x'-0.00207Y/(1+Y+2*S)-0.059*(y'/(1+y+2*s)+y*(y'+2*s')/(1+y+2*s)^2);         
s'=0.5*(-x'-0.00056S/(1+Y+2*S)-0.086*(s'/(1+y+2*s)+s*(y'+2*s')/(1+y+2*s)^2))
其中y'=dy/dt, s'=ds/dt，x'=dx/dt
赞一下 回复此楼
4楼2012-01-05 17:49:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dbb627

荣誉版主 (著名写手)

【答案】应助回帖

感谢参与，应助指数 +1
>> syms t
>> Y=sym('Y(t)')

Y =

Y(t)

>> S=sym('S(t)')

S =

S(t)

>> F=Y/(1+Y+2*S)

F =

Y(t)/(2*S(t) + Y(t) + 1)

>> dF=diff(F,t)

dF =

diff(Y(t), t)/(2*S(t) + Y(t) + 1) - (Y(t)*(2*diff(S(t), t) + diff(Y(t), t)))/(2*S(t) + Y(t) + 1)^2

1：先将公式中的项“d[Y/(1+Y+2S)]/dt”求导展开变为：
    y'/(1+y+2*s)+y*(y'+2*s')/(1+y+2*s)^2，其中y'=dy/dt, s'=ds/dt
2：同理将“d[S/(1+Y+2S)]/dt”变为：
   s'/(1+y+2*s)+s*(y'+2*s')/(1+y+2*s)^2
用matlab ode15i解隐式微分方程
CODE:
odefun=@(t,x,dx)[dx(1)-1.98*x(3)*(1-x(1))+0.73*x(1)*x(2)
    dx(2)-(dx(1)-0.00207*x(2)/(1+x(2)+2*x(3))-0.059*(dx(2)/(1+x(2)+2*x(3))+x(2)*(dx(2)+2*dx(3))/(1+x(2)+2*x(3))^2))
    2*dx(3)+dx(1)+0.00056*x(3)/(1+x(2)+2*x(3))+0.086*(dx(3)/(1+x(2)+2*x(3))+x(3)*(dx(2)+2*dx(3))/(1+x(2)+2*x(3))^2)];
t0=0;
x0=[0 0 1]';
fix_x0=ones(3,1);
dx0=[0 1 1]';
fix_dx0=zeros(3,1);
[x02,dx02]=decic(odefun,t0,x0,fix_x0,dx0,fix_dx0);
solution=ode15i(odefun,[0 3],x02,dx02);
t=solution.x;
x=solution.y;
plot(t,x)
legend('X','Y','S')

[ Last edited by dbb627 on 2012-1-5 at 18:29 ]
赞一下 回复此楼
The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
5楼2012-01-05 18:23:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dbb627

荣誉版主 (著名写手)

【答案】应助回帖


phoenix708(金币+30): ★★★很有帮助 2012-01-05 20:02:35
做图如下




[ Last edited by dbb627 on 2012-1-5 at 18:30 ]
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

phoenix708 phoenix708
The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
6楼2012-01-05 18:25:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

芝芝花

铜虫 (初入文坛)
引用回帖:
6楼: Originally posted by dbb627 at 2012-01-05 18:25:08:
做图如下

[ Last edited by dbb627 on 2012-1-5 at 18:30 ]

Thank you!
回复此楼
7楼2012-01-05 19:23:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

scumat

金虫 (正式写手)
引用回帖:
6楼: Originally posted by dbb627 at 2012-01-05 18:25:08:
做图如下

[ Last edited by dbb627 on 2012-1-5 at 18:30 ]

佩服！可以得全分。
回复此楼
8楼2012-01-05 19:40:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

phoenix708

铁杆木虫 (著名写手)
送鲜花一朵
引用回帖:
楼: Originally posted by dbb627 at 2012-01-05 18:25:08:
做图如下

[ Last edited by dbb627 on 2012-1-5 at 18:30 ]

很厉害啊，不过这个结果和我想要的结果有差距，我纠正了第二个和第三个方程的系数，肯定是对的，劳驾大侠最后帮我计算一次！多谢多谢！
dX/dt=1.98S(1-X)-0.73XY
dY/dt=dX/dt-4.6Y/(1+Y+2S)-0.059d[Y/(1+Y+2S)]/dt
2dS/dt=-dX/dt-1.25S/(1+Y+2S)-0.086d[S/(1+Y+2S)]/dt
当t=0时,X=0,Y=0,S=1,
其中 t的最大值为3，
求X，Y, S与t的关系图
谢谢!
赞一下 回复此楼
9楼2012-01-05 19:52:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kingslcq

新虫 (正式写手)

爱科学的PH.D & Prof.

【答案】应助回帖

感谢参与，应助指数 +1
phoenix708(金币+1): 2012-01-05 20:00:10
http://muchong.com/bbs/viewthread.php?tid=4017234
赞一下 回复此楼
10楼2012-01-05 19:54:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 phoenix708 的主题更新
281/3返回列表123下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 274求调剂 +11 顾九笙要谦虚 2026-03-24 16/800 2026-03-26 09:55 by stone_128
[考研] 【2026考研调剂】制药工程 284分 求相关专业调剂名额 +4 袁奂奂 2026-03-25 8/400 2026-03-25 14:32 by lbsjt
[考研] 359求调剂 +3 王了个楠 2026-03-25 3/150 2026-03-25 12:50 by Dyhoer
[考研] 考研化学308分求调剂 +10 你好明天你好 2026-03-23 11/550 2026-03-25 10:23 by userper
[考研] 0703化学求调剂 +6 奶油草莓. 2026-03-22 7/350 2026-03-25 10:00 by shangxh
[考研] 07化学280分求调剂 +7 722865 2026-03-23 7/350 2026-03-25 09:29 by aa331100
[考研] 0856材料专硕353求调剂 +6 NIFFFfff 2026-03-20 6/300 2026-03-24 21:02 by hello七七
[考研] 306求0703调剂一志愿华中师范 +10 纸鱼ly 2026-03-21 11/550 2026-03-24 17:22 by qingfeng258
[考研] 求调剂 +6 研研，接电话 2026-03-24 7/350 2026-03-24 17:01 by barlinike
[考研] 291求调剂 +3 HanBeiNingZC 2026-03-24 3/150 2026-03-24 16:34 by barlinike
[考研] 085404电子信息284分求调剂 +4 13659058978 2026-03-24 4/200 2026-03-24 12:15 by syl20081243
[论文投稿] 急发核心期刊论文 +3 贤达问津 2026-03-23 5/250 2026-03-23 17:13 by 妹子不好惹
[考研] 328求调剂 +4 LHHL66 2026-03-23 4/200 2026-03-23 14:55 by lbsjt
[考研] 一志愿070300浙大化学358分，求调剂！ +4 酥酥鱼.. 2026-03-21 4/200 2026-03-23 08:12 by Iveryant
[考研] 石河子大学（211、双一流）硕博研究生长期招生公告 +3 李子目 2026-03-22 3/150 2026-03-22 21:01 by 怎么释怀
[考研] 一志愿华中科技大学071000，求调剂 +4 沿岸有贝壳6 2026-03-21 4/200 2026-03-22 07:21 by ilovexiaobin
[考研] 一志愿东华大学控制学硕320求调剂 +3 Grand777 2026-03-21 3/150 2026-03-21 19:23 by 简之-
[考研] 求调剂 +3 白QF 2026-03-21 3/150 2026-03-21 13:12 by zhukairuo
[考研] 332求调剂 +3 凤凰院丁真 2026-03-20 3/150 2026-03-21 10:27 by luoyongfeng
[考研] 一志愿西安交通大学 学硕 354求调剂211或者双一流 +3 我想要读研究生 2026-03-20 3/150 2026-03-20 20:13 by JourneyLucky
信息提示
关闭
请填处理意见
关闭