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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 会用matlab的同志帮帮忙哈！
281/3返回列表123下一页
查看: 1361  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上模板改不了页边距吧？ +5 ieewxg 2026-02-25 6/300 2026-03-01 00:10 by addressing
[考研] 307求调剂 +4 73372112 2026-02-28 6/300 2026-03-01 00:04 by ll247
[考研] 317一志愿华南理工电气工程求调剂 +5 Soliloquy_Q 2026-02-28 8/400 2026-02-28 23:36 by xyx2012xyx
[考研] 292求调剂 +3 yhk_819 2026-02-28 3/150 2026-02-28 21:57 by gaoxiaoniuma
[考研] 290求调剂 +5 材料专硕调剂； 2026-02-28 6/300 2026-02-28 21:40 by gaoxiaoniuma
[考博] 26申博 +4 想申博！ 2026-02-26 4/200 2026-02-28 21:37 by limorning
[考研] 材料学调剂 +5 提神豆沙包 2026-02-28 5/250 2026-02-28 21:34 by gaoxiaoniuma
[考研] 311求调剂 +8 南迦720 2026-02-28 8/400 2026-02-28 21:30 by gaoxiaoniuma
[考研] 材料类求调剂 +6 wana_kiko 2026-02-28 6/300 2026-02-28 21:20 by gaoxiaoniuma
[考研] 求调剂 +4 repeatt?t 2026-02-28 4/200 2026-02-28 21:16 by gaoxiaoniuma
[考研] 284求调剂 +4 天下熯 2026-02-28 4/200 2026-02-28 21:13 by gaoxiaoniuma
[考研] 高分子化学与物理调剂 +4 好好好1233 2026-02-28 7/350 2026-02-28 20:42 by 好好好1233
[考研] 298求调剂 +8 人间唯你是清欢 2026-02-28 11/550 2026-02-28 20:26 by L135790
[考研] 285求调剂 +5 满头大汗的学生 2026-02-28 5/250 2026-02-28 18:10 by 材料专硕调剂；
[高分子] 求环氧树脂研发1名 +3 孙xc 2026-02-25 11/550 2026-02-28 16:57 by ichall
[考研] 0856调剂 +3 刘梦微 2026-02-28 3/150 2026-02-28 13:22 by houyaoxu
[考研] 寻找调剂 +3 LYidhsjabdj 2026-02-28 3/150 2026-02-28 12:59 by miniwendy
[硕博家园] 博士自荐 +6 科研狗111 2026-02-26 9/450 2026-02-28 12:32 by seaskyy
[考研] 272求调剂 +3 田智友 2026-02-28 3/150 2026-02-28 12:31 by 王加浩to
[基金申请] 面上可以超过30页吧？ +12 阿拉贡aragon 2026-02-22 13/650 2026-02-26 22:09 by Hahaxia
信息提示
关闭
请填处理意见
关闭