看图说话——布鲁塞尔问题的分叉图

胡晓宇 · 发表于 2008-9-5 21:02

多谢liliangbiao

cam_1980 · 发表于 2008-9-5 21:36

刚看到你的回复，那我也在这里多谢liliangbiao

胡晓宇 · 发表于 2008-9-9 09:08

liliangbiao，我仔细看了一下你的程序，先算了60个周期，并舍弃，将最终值作为初值计算后40个周期，那么后四十个周期每个周期做一次计算是否是必要的？tspan=[i*period:step:(i+1)*period];
[t,Y]=ode45(@duffing,tspan,y0);
能不能改成
tspan=[60*period：step：100*period]
[t,Y]=ode45(@duffing,tspan,y0);
再隔100点取一个点
因为多加一个循环好像速度慢很多，不知道是否可以？

octopussheng · 发表于 2008-9-9 09:56

他的这个程序思想其实是很好的，可以把前面的瞬态解去除，

你的意思是不是前面的60个周期直接不算了，时间t0直接从60周期开始？如果是这样的话就不对了哦！

胡晓宇 · 发表于 2008-9-9 10:06

前四十个周期一样算，并将瞬态去除
只是在进行后60个周期的时候，不一个周期一个周期的循环计算，而是t_span将后60个周期取整体，就是将tspan=[i*period:step:(i+1)*period];改成
tspan=[60*period：step：100*period]

octopussheng · 发表于 2008-9-9 11:06

呵呵，这样应该是可行的！

胡晓宇 · 发表于 2008-9-9 16:45

如果可行，速度又快了一点，我先试算一下看看吧

无水1324 · 发表于 2008-9-12 10:45

应该就是计算的策略不同，所以呢两种应该是一样的

cam_1980 · 发表于 2008-9-16 10:53

确实是一样的结果，我就是这样算出来的，只是在最后绘图的时候再取存储点的后面数据来绘制，这样就能得到轮廓清晰的图。

账号		自动登录	找回密码
密码			我要加入

[稳定性与分岔] 看图说话——布鲁塞尔问题的分叉图

回复 31楼胡晓宇的帖子

回复 30楼 liliangbiao 的帖子

回复 34楼 octopussheng 的帖子

回复 36楼 octopussheng 的帖子

回复 35楼胡晓宇的帖子

回复 37楼胡晓宇的帖子