版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

>论坛更新日志 (4428)
>文献求助 (299)
>导师招生 (251)
>虫友互识 (212)
>论文投稿 (166)
>学术会议 (155)
>考博 (151)
>硕博家园 (78)
>教师之家 (71)
>找工作 (70)
>基金申请 (60)
>休闲灌水 (53)
>论文道贺祈福 (47)
>博后之家 (43)
>考研 (40)
>招聘信息布告栏 (36)

返回列表

brqhl_ing

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 262.9
散金: 11
红花: 1
帖子: 115
在线: 19.5小时
虫号: 854693
注册: 2009-09-23
性别: GG
专业: 物理电子学

[交流] 【讨论】大家看看这是怎么了？？？？？？？？？已有1人参与

function dzdt=jisuan(t,x)
dzdt=zeros(3,1);
dx(1)=-8/3*x(1)+x(2)*x(3)+0.5*(x(1)-x(2));
dx(2)=-10*x(2)+10*x(3);
dx(3)=-x(1)*x(2)+28*x(2)-x(3);

[t,y]=ode45('jisuan',[0,10],[0.3,0.8,1.3])

t =

0
0.2500
0.5000
0.7500
1.0000
1.2500
1.5000
1.7500
2.0000
2.2500
2.5000
2.7500
3.0000
3.2500
3.5000
3.7500
4.0000
4.2500
4.5000
4.7500
5.0000
5.2500
5.5000
5.7500
6.0000
6.2500
6.5000
6.7500
7.0000
7.2500
7.5000
7.7500
8.0000
8.2500
8.5000
8.7500
9.0000
9.2500
9.5000
9.7500
10.0000

y =

0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
0.3000 0.8000 1.3000
错误出在哪？？？

回复此楼

» 猜你喜欢

售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急已经有6人回复
售SCI文章，我:8O.5.5.1O.54,科目全,可十急已经有8人回复
售SCI文章，我:8O.5.5.1O.54,科目全,可十急已经有5人回复
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急已经有10人回复
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急已经有6人回复
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急已经有9人回复
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急已经有8人回复
国社科又开始会评了，不知道这次命运如何已经有14人回复
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急已经有8人回复
麻烦专家们看看评委们的意见（F口面上）已经有13人回复

1楼 2010-06-04 22:51:54

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师

应助: 44 (小学生)
金币: 17724.5
红花: 17
帖子: 2413
在线: 546.7小时
虫号: 496517
注册: 2008-01-19
专业: 化学反应工程

★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
robert2020(金币+1):多谢应助！ 2010-06-04 23:20:47
zzuwangshilei(金币+1):多谢应助，还望详细说明下 2010-06-05 10:58:58

function test
clear
clc
x0=eps*ones(1,3);
[t,x]=ode45(@OdeFun,[0:0.02:8],x0);
plot(t,x(:,1),'r.-',t,x(:,2),'b.-',t,x(:,3),'k.-')

function dx = OdeFun(t,x)
dx1=-8/3*x(1)+x(2)*x(3)+0.5*(x(1)-x(2));
dx2= -10*x(2)+10*x(3);
dx3=-x(1)*x(2)+28*x(2)-x(3);
dx=[dx1;dx2;dx3];

赞一下(3人)