24小时热门版块排行榜    

查看: 1312  |  回复: 27

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


» 本帖已获得的红花(最新10朵)

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的回帖
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 的主题更新
信息提示
请填处理意见