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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 会用matlab的同志帮帮忙哈！
53/1返回列表
查看: 1384  |  回复: 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的回帖

dbb627

荣誉版主 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
phoenix708(金币+20): ★★★很有帮助 答案很完善！ 2012-01-06 09:07:06
cenwanglai(金币+10): 真是不知道如何感谢你~ 2012-01-11 20:47:04
CODE:
odefun=@(t,x,dx)[dx(1)-1.98*x(3)*(1-x(1))+0.73*x(1)*x(2)
    dx(2)-(dx(1)-4.6*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)+1.25*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')
display('t  X   Y  S')
[t‘ x’]

t      X          Y            S
       0         0            0             1.0000
    0.0000    0.0000    0.0000    1.0000
    0.0000    0.0001    0.0001    1.0000
    0.0001    0.0002    0.0002    0.9999
    0.0002    0.0004    0.0004    0.9998
    0.0004    0.0008    0.0008    0.9995
    0.0006    0.0012    0.0012    0.9993
    0.0008    0.0016    0.0016    0.9991
    0.0012    0.0024    0.0023    0.9986
    0.0016    0.0032    0.0031    0.9981
    0.0024    0.0048    0.0047    0.9971
    0.0041    0.0080    0.0078    0.9952
    0.0073    0.0144    0.0140    0.9914
    0.0139    0.0269    0.0261    0.9839
    0.0204    0.0392    0.0378    0.9766
    0.0270    0.0512    0.0492    0.9694
    0.0335    0.0630    0.0602    0.9623
    0.0466    0.0858    0.0812    0.9484
    0.0597    0.1078    0.1009    0.9351
    0.0858    0.1491    0.1367    0.9098
    0.1120    0.1873    0.1681    0.8861
    0.1382    0.2226    0.1954    0.8640
    0.1643    0.2553    0.2190    0.8432
    0.1905    0.2856    0.2394    0.8237
    0.2429    0.3399    0.2715    0.7881
    0.2952    0.3868    0.2940    0.7564
    0.3475    0.4277    0.3085    0.7278
    0.3999    0.4634    0.3167    0.7020
    0.4522    0.4948    0.3198    0.6785
    0.5045    0.5226    0.3187    0.6570
    0.5569    0.5473    0.3144    0.6370
    0.6615    0.5892    0.2988    0.6011
    0.7662    0.6235    0.2773    0.5693
    0.8709    0.6521    0.2530    0.5406
    0.9755    0.6766    0.2278    0.5142
    1.0802    0.6978    0.2031    0.4896
    1.1849    0.7164    0.1797    0.4665
    1.2895    0.7330    0.1581    0.4446
    1.3942    0.7480    0.1384    0.4237
    1.4989    0.7616    0.1207    0.4038
    1.6035    0.7739    0.1051    0.3848
    1.7082    0.7852    0.0914    0.3665
    1.8129    0.7956    0.0795    0.3490
    1.9175    0.8052    0.0691    0.3321
    2.0222    0.8140    0.0601    0.3159
    2.1269    0.8221    0.0524    0.3004
    2.3362    0.8365    0.0400    0.2712
    2.5246    0.8476    0.0316    0.2470
    2.7130    0.8572    0.0252    0.2245
    2.8826    0.8648    0.0207    0.2058
    3.0000    0.8695    0.0182    0.1936
赞一下(1人) 回复此楼

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

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.
22楼2012-01-05 22:56:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 28 个回答

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的回帖
查看全部 28 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿天津大学339材料与化工求调剂 +3 江往卖鱼 2026-03-26 3/150 2026-03-26 09:42 by 王小欠i
[考研] 303求调剂 +6 蓝山月 2026-03-25 6/300 2026-03-25 22:47 by 418490947
[考研] 材料求调剂 +4 .m.. 2026-03-25 4/200 2026-03-25 21:30 by peike
[考研] 考研一志愿苏州大学初始315（英一）求调剂 +3 sbdksD 2026-03-24 4/200 2026-03-25 18:16 by xcjcqu
[考研] 招08考数学 +8 laoshidan 2026-03-20 17/850 2026-03-25 17:52 by 一个红太阳
[考研] 材料277求调剂 +4 min3 2026-03-24 4/200 2026-03-25 15:29 by fch1983
[考研] 0854电子信息求调剂 +7 α____ 2026-03-22 9/450 2026-03-25 13:37 by α____
[考研] 281求调剂 +4 Koxui 2026-03-24 5/250 2026-03-25 11:38 by userper
[考研] 一志愿吉林大学材料与化工303分求调剂 +4 为学666 2026-03-24 4/200 2026-03-25 11:27 by BruceLiu320
[考研] 一志愿北化315 求调剂 +3 akrrain 2026-03-24 3/150 2026-03-24 19:35 by 了了了了。。
[考研] 080500求调剂 +3 zzzzfan 2026-03-24 3/150 2026-03-24 16:38 by barlinike
[考研] 【双一流院校新能源、环境材料，材料加工与模拟招收大量调剂】 +4 Higraduate 2026-03-22 7/350 2026-03-24 11:23 by 种大树
[考研] 344求调剂 +3 desto 2026-03-24 3/150 2026-03-24 10:09 by 搏击518
[考研] 298求调剂 +8 上岸6666@ 2026-03-20 8/400 2026-03-23 11:02 by laoshidan
[考研] 石河子大学（211、双一流）硕博研究生长期招生公告 +3 李子目 2026-03-22 3/150 2026-03-22 21:01 by 怎么释怀
[考研] 材料与化工085600，总分304，本科有两篇sci参与，求调剂 +4 幸运的酱酱 2026-03-22 5/250 2026-03-22 20:15 by edmund7
[考研] 求调剂院校信息 +6 CX 330 2026-03-21 6/300 2026-03-22 15:25 by 无懈可击111
[考研] 求调剂 +5 Zhangbod 2026-03-21 7/350 2026-03-22 13:13 by Zhangbod
[考研] 材料与化工（0856）304求 B区 调剂 +3 邱gl 2026-03-21 3/150 2026-03-21 13:47 by lature00
[考研] 261求B区调剂，科研经历丰富 +3 牛奶很忙 2026-03-20 4/200 2026-03-20 19:34 by JourneyLucky
信息提示
关闭
请填处理意见
关闭