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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 会用matlab的同志帮帮忙哈！
282/3返回列表上一页123下一页
查看: 1371  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学工程085602 305分求调剂 +10 RichLi_ 2026-03-25 10/500 2026-03-26 02:17 by BruceLiu320
[考研] 一志愿上海交大生物与医药专硕324分，求调剂 +6 jiajunX 2026-03-22 6/300 2026-03-25 23:05 by licg0208
[考研] 316求调剂 +9 梁茜雯 2026-03-19 9/450 2026-03-25 23:01 by peike
[考研] 材料与化工328分调剂 +6 。，。，。，。i 2026-03-23 6/300 2026-03-25 22:30 by 418490947
[考研] 299求调剂 +7 某某某某位 2026-03-21 8/400 2026-03-25 20:34 by 热情沙漠
[考研] 【2026考研调剂】制药工程 284分 求相关专业调剂名额 +4 袁奂奂 2026-03-25 8/400 2026-03-25 14:32 by lbsjt
[考研] 求调剂323材料与化工 +4 1124361 2026-03-24 4/200 2026-03-25 11:19 by shulmg
[考研] 299求调剂 +7 shxchem 2026-03-20 9/450 2026-03-25 10:41 by lbsjt
[考研] 一志愿 西北大学 ，070300化学学硕，总分287，双非一本，求调剂。 +4 晨昏线与星海 2026-03-20 4/200 2026-03-25 10:16 by allen-yin
[考研] B区考研调剂 +4 yqdszhdap－ 2026-03-22 5/250 2026-03-25 08:51 by baoball
[考研] 材料292调剂 +8 橘颂思美人 2026-03-23 8/400 2026-03-24 16:33 by laoshidan
[考研] 321求调剂 +4 Ymlll 2026-03-24 4/200 2026-03-24 14:44 by sprinining
[考研] 一志愿北京化工大学材料与化工 264分各科过A区国家线 +3 哈哈157349 2026-03-21 3/150 2026-03-24 14:11 by zhyzzh
[考研] 一志愿中南化学（0703）总分337求调剂 +9 niko- 2026-03-19 10/500 2026-03-22 16:08 by ColorlessPI
[考研] 260求调剂 +3 朱芷琳 2026-03-20 4/200 2026-03-22 15:12 by 朱芷琳
[考研] 求调剂 +7 Auroracx 2026-03-22 7/350 2026-03-22 12:38 by 素颜倾城1988
[考研] 296求调剂 +4 www_q 2026-03-20 4/200 2026-03-21 17:26 by 学员8dgXkO
[基金申请] 学校已经提交到NSFC，还能修改吗？ 40+4 babangida 2026-03-19 9/450 2026-03-21 16:12 by babangida
[考研] 329求调剂 +9 想上学吖吖 2026-03-19 9/450 2026-03-20 22:01 by luoyongfeng
[考研] 招收调剂硕士 +4 lidianxing 2026-03-19 12/600 2026-03-20 12:25 by lidianxing
信息提示
关闭
请填处理意见
关闭