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

24小时热门版块排行榜    

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

phoenix708

铁杆木虫 (著名写手)
引用回帖:
楼: Originally posted by hls85915 at 2012-01-05 21:16:31:
你黏在MATLAB里就OK了

我没有matlab的嘛，我没用过，呵呵！麻烦帮忙运算下，把数据也发上来下，我要自己用origin画图，谢谢
赞一下 回复此楼
21楼2012-01-05 21:19:05
已阅   回复此楼   关注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的回帖

zl030091126
23楼2012-01-06 03:22:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

phoenix708

铁杆木虫 (著名写手)
送鲜花一朵
引用回帖:
楼: Originally posted by dbb627 at 2012-01-05 22:56:01:
[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)+ ...

Thank you!
回复此楼
24楼2012-01-06 09:07:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

笨笨小熊熊
25楼2012-03-06 11:26:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jackzhang19

银虫 (小有名气)
★ ★
cenwanglai(金币+2): ~~ 2012-03-10 08:43:18
function dy=jack(t,y)
dy=zeros(6,1);
dy(1)=1.98*y(5)*(1-y(1))-0.73*y(1)*y(3);
dy(3)=y(2)-0.00207*y(3)/(1+y(3)+2*y(5))-0.59*(y(4)*(1+y(3)+2*y(5))-y(3)*(y(4)+2*y(6)))/(1+y(3)+2*y(5))^2;
dy(5)=1/2*(-y(2)-0.00056*y(5)/(1+y(3)+2*y(5))-0.086*(y(6)*(1+y(3)+2*y(5))-y(5)*(y(4)+2*y(6)))/(1+y(3)+2*y(5))^2);
[t,y] = ode23s(@jack,tspan,[0 ,1.98,0 ,1.0067796593,1,-0.9805981300402]);
plot(t,y(:,1),'r',t,y(:,3),'b',t,y(:,5),'y')
赞一下 回复此楼
路漫漫其修远兮吾将上下而求索
26楼2012-03-07 17:26:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jackzhang19

银虫 (小有名气)
初值中的对x,y,s在零点的导数值设置的比较粗糙你可以再修改一下结果会更好
赞一下 回复此楼
路漫漫其修远兮吾将上下而求索
27楼2012-03-07 17:32:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

中和达人

金虫 (小有名气)
都是牛人啊
回复此楼
务实求新
28楼2012-03-09 11:54:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 phoenix708 的主题更新
283/3返回列表上一页123
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] T^T +6 m584r2z4j3 2026-06-11 7/350 2026-06-13 23:10 by 3eaav2pea1
[找工作] 售SCI文章，我:8.O.5.51O.5.4,科目齐全,可十急 +3 q1zg6x4lf3 2026-06-12 3/150 2026-06-13 19:07 by 1js5n6upxm
[公派出国] 囍 +4 q1zg6x4lf3 2026-06-12 4/200 2026-06-13 19:05 by 1js5n6upxm
[公派出国] π_π +4 m584r2z4j3 2026-06-11 4/200 2026-06-13 18:33 by lkmi35jdzt
[考研] T^T +6 m584r2z4j3 2026-06-11 6/300 2026-06-13 18:29 by lkmi35jdzt
[文学芳草园] 咸菜 +3 myrtle 2026-06-10 3/150 2026-06-13 18:29 by ahsoarli
[考博] π_π +4 m584r2z4j3 2026-06-11 4/200 2026-06-13 18:25 by lkmi35jdzt
[考研] o>_<o +5 m584r2z4j3 2026-06-11 5/250 2026-06-13 18:22 by lkmi35jdzt
[考研] ╭(╯ε╰)╮ +4 m584r2z4j3 2026-06-11 4/200 2026-06-13 14:22 by 80j16ujcu8
[考博] T^T +4 jk31rb6vyz 2026-06-12 4/200 2026-06-13 12:09 by 80j16ujcu8
[考博] T^T +3 m584r2z4j3 2026-06-11 3/150 2026-06-13 12:04 by 80j16ujcu8
[找工作] 售SCI文章，我:8.O.5.51O.5.4,科目齐全,可十急 +4 q1zg6x4lf3 2026-06-12 4/200 2026-06-13 08:53 by t4hz8376wy
[考博] →_→ +3 q1zg6x4lf3 2026-06-12 3/150 2026-06-13 08:49 by t4hz8376wy
[找工作] *囍 +3 m584r2z4j3 2026-06-12 3/150 2026-06-13 08:29 by t4hz8376wy
[基金申请] 面上项目没有好文章就没希望了吗？ +19 1234567wang 2026-06-07 21/1050 2026-06-12 21:24 by littlezl
[考博] [急招]吉林大学招收26年9月入学考核制博士生1名，截止今晚24点前 +3 hehejin 2026-06-10 4/200 2026-06-12 18:24 by 冻雨天
[基金申请] 前几天时间戳更新了 +8 vito刘 2026-06-11 13/650 2026-06-12 14:17 by BlakeReary
[教师之家] AASSDDFF +3 m584r2z4j3 2026-06-11 3/150 2026-06-12 11:28 by zs48m4lz8g
[考博] 还有课题组有博士名额吗 +5 冻雨天 2026-06-07 6/300 2026-06-08 20:59 by RCS_tyro
[考博] 博士申请 +3 张城斐 2026-06-07 3/150 2026-06-08 09:15 by A大朋同学
信息提示
关闭
请填处理意见
关闭