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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 会用matlab的同志帮帮忙哈！
282/3返回列表上一页123下一页
查看: 1357  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 317一志愿华南理工电气工程求调剂 +3 Soliloquy_Q 2026-02-28 3/150 2026-02-28 21:43 by xyx2012xyx
[考研] 304求调剂 +3 52hz~~ 2026-02-28 4/200 2026-02-28 21:41 by gaoxiaoniuma
[考研] 295求调剂 +5 19171856320 2026-02-28 5/250 2026-02-28 21:39 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
[考研] 085600材料工程一志愿中科大总分312求调剂 +8 吃宵夜1 2026-02-28 10/500 2026-02-28 20:27 by L135790
[考研] 298求调剂 +8 人间唯你是清欢 2026-02-28 11/550 2026-02-28 20:26 by L135790
[基金申请] 面上模板改不了页边距吧？ +5 ieewxg 2026-02-25 5/250 2026-02-28 20:11 by iwuli
[考研] 276求调剂 +3 路lyh123 2026-02-28 4/200 2026-02-28 19:45 by 路lyh123
[考博] 博士推荐 +5 花儿笑? 2026-02-21 6/300 2026-02-28 18:53 by nxgogo
[考研] 化工专硕348，一志愿985求调剂 +3 弗格个 2026-02-28 5/250 2026-02-28 17:04 by sandychj
[考博] 博士自荐 +3 kkluvs 2026-02-28 3/150 2026-02-28 16:59 by StarAura
[高分子] 求环氧树脂研发1名 +3 孙xc 2026-02-25 11/550 2026-02-28 16:57 by ichall
[考研] 265分求调剂不调专业和学校有行学上就 +4 礼堂丁真258 2026-02-28 6/300 2026-02-28 16:18 by 求调剂zz
[考研] 寻找调剂 +3 LYidhsjabdj 2026-02-28 3/150 2026-02-28 12:59 by miniwendy
[考研] 304求调剂 +5 曼殊2266 2026-02-28 6/300 2026-02-28 12:44 by 迷糊CCPs
[考研] 272求调剂 +3 田智友 2026-02-28 3/150 2026-02-28 12:31 by 王加浩to
[硕博家园] 【博士招生】太原理工大学2026化工博士 +4 N1ce_try 2026-02-24 8/400 2026-02-26 08:40 by N1ce_try
信息提示
关闭
请填处理意见
关闭