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

24小时热门版块排行榜    

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

dingd

铁杆木虫 (职业作家)
该问题与回答过的“芝芝花”的问题几乎一样啊，参见：
1：http://muchong.com/bbs/viewthread.php?tid=3977878&fpage=1
2：http://muchong.com/bbs/viewthread.php?tid=4003917&fpage=1

是同一个人还是在做相同研究的人？
赞一下 回复此楼

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

phoenix708
11楼2012-01-05 19:55:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

phoenix708

铁杆木虫 (著名写手)
引用回帖:
楼: Originally posted by dingd at 2012-01-05 19:55:42:
该问题与回答过的“芝芝花”的问题几乎一样啊，参见：
1：http://muchong.com/bbs/viewthread.php?tid=3977878&fpage=1
2：[url]http://muchong.com/bbs/viewthread.php?tid=4003917&fpage=1[/u ...

同一个课题组的，最近共同完成点数据，不好意思！打扰大家了，这个计算好了就完了！
赞一下 回复此楼
12楼2012-01-05 19:58:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

phoenix708

铁杆木虫 (著名写手)
送鲜花一朵
引用回帖:
楼: Originally posted by dingd at 2012-01-05 19:55:42:
该问题与回答过的“芝芝花”的问题几乎一样啊，参见：
1：http://muchong.com/bbs/viewthread.php?tid=3977878&fpage=1
2：[url]http://muchong.com/bbs/viewthread.php?tid=4003917&fpage=1[/u ...

同学推荐说这个论坛里高手很多，大侠能帮忙算一下吗，最近一直在算这方面的东西，计算工作马上就要完了，辛苦大家了，谢谢！
赞一下 回复此楼
13楼2012-01-05 20:06:35
已阅   回复此楼   关注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 ]

版主好人做到底了，劳驾再帮我计算一下，这是这段时间以来修改的最后结果，打扰大家了，谢谢！把计算后的结果即：图中对应的x,y,s,t的值也发上来下，因为要自己画图！对不起，打扰了！
赞一下 回复此楼
14楼2012-01-05 20:32:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
phoenix708(金币+10): ★有帮助 很热心的虫友，好人啊！ 2012-01-05 20:45:36
6楼给的图应该是对的吧？1stOpt求解代码及结果如下。如果用的频繁，让导师弄个高版本的1stOpt，方便。
CODE:
Variable t=[0:0.1:3],x=0,y=0,s=1;
Plot t[x],x,y,s[y2];
ODEFunction 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))

t        x(t)        y(t)        s(t)
0        0        0        1
0.1        0.172015852079919        0.168692414507187        0.91520287028323
0.2        0.301008914823512        0.295181857886797        0.851612843575276
0.3        0.398410079826341        0.390681232086829        0.803595330591451
0.4        0.472340820022573        0.463157198420494        0.767147038400838
0.5        0.528679098258018        0.518375712926536        0.739370558584698
0.6        0.57174089232378        0.560570533610007        0.718138426693297
0.7        0.604731032316469        0.592885451556486        0.701870893031543
0.8        0.630050124550394        0.61767535899965        0.68938463184859
0.9        0.649508597277285        0.636716213168225        0.679787232494888
1        0.664478996144567        0.651354428633045        0.672402128859789
1.1        0.676006150214996        0.662614944014651        0.666714295511898
1.2        0.684887929226023        0.671280440772958        0.662330440202897
1.3        0.691735117178001        0.677950076423071        0.658949490500964
1.4        0.697016164438162        0.683083384430652        0.656340534220312
1.5        0.70109088349502        0.687033328390857        0.654326209833513
1.6        0.704235956743193        0.690071324190412        0.652770132963963
1.7        0.706664329451613        0.692406263797917        0.651567337044101
1.8        0.708540003734212        0.69419902762291        0.650636980949813
1.9        0.709989354031696        0.695573584610092        0.649916771280749
2        0.711109799846648        0.696625499895661        0.649358687325529
2.1        0.711976463612829        0.697428465951553        0.648925699211464
2.2        0.712647288165669        0.698039322651576        0.648589245360499
2.3        0.713166974013066        0.698501919601665        0.648327291698946
2.4        0.713570010858805        0.698850089962322        0.648122837335628
2.5        0.7138830130987        0.699109941491468        0.647962763331418
2.6        0.714126519906615        0.699301622367755        0.647836945387419
2.7        0.714316383138106        0.699440682676658        0.647737569709322
2.8        0.714464837723565        0.699539124429341        0.647658605381513
2.9        0.714581327361879        0.699606211539531        0.647595397360315
3        0.714673141559223        0.699649094736309        0.647544352460397
赞一下 回复此楼
15楼2012-01-05 20:37:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

phoenix708

铁杆木虫 (著名写手)
引用回帖:
楼: Originally posted by dingd at 2012-01-05 20:37:38:
6楼给的图应该是对的吧？1stOpt求解代码及结果如下。如果用的频繁，让导师弄个高版本的1stOpt，方便。

[code]
Variable t=[0:0.1:3],x=0,y=0,s=1;
Plot t[x],x,y,s[y2];
ODEFunction x'=1.98S*(1-X)-0.73*X ...

版主是回复对的，不过这个结果和我想要的结果有差距，我纠正了第二个和第三个方程的系数，劳驾大侠最后帮我计算一次！多谢多谢！
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的关系图
谢谢!
赞一下 回复此楼
16楼2012-01-05 20:44:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hls85915

荣誉版主 (文坛精英)

科研工具版——沙漠之胡

优秀版主优秀版主

【答案】应助回帖

感谢参与，应助指数 +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
赞一下 回复此楼
哈哈
17楼2012-01-05 20:58:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

phoenix708

铁杆木虫 (著名写手)
引用回帖:
楼: Originally posted by hls85915 at 2012-01-05 20:58:47:
>> 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 ...

还要有图哈，不然看不懂！
t的最大值为3，
求X，Y, S与t的关系图（t为横坐标）
赞一下 回复此楼
18楼2012-01-05 21:15:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hls85915

荣誉版主 (文坛精英)

科研工具版——沙漠之胡

优秀版主优秀版主
引用回帖:
18楼: Originally posted by phoenix708 at 2012-01-05 21:15:10:
还要有图哈，不然看不懂！
t的最大值为3，
求X，Y, S与t的关系图（t为横坐标）

你黏在MATLAB里就OK了
赞一下 回复此楼
哈哈
19楼2012-01-05 21:16:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

phoenix708

铁杆木虫 (著名写手)
引用回帖:
楼: Originally posted by hls85915 at 2012-01-05 20:58:47:
>> 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 ...

注意我的第二个方程和第三个方程的系数是改了的哟！见15楼的方程
赞一下 回复此楼
20楼2012-01-05 21:17:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 phoenix708 的主题更新
282/3返回列表上一页123下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭