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

多谢liliangbiao

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

回复 31楼 胡晓宇 的帖子

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

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

回复 30楼 liliangbiao 的帖子

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

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

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

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

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

回复 34楼 octopussheng 的帖子

前四十个周期一样算，并将瞬态去除
只是在进行后60个周期的时候，不一个周期一个周期的循环计算，而是t_span将后60个周期取整体，就是将tspan=;改成
tspan=

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

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

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

回复 36楼 octopussheng 的帖子

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

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

回复 35楼 胡晓宇 的帖子

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

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

回复 37楼 胡晓宇 的帖子

确实是一样的结果，我就是这样算出来的，只是在最后绘图的时候再取存储点的后面数据来绘制，这样就能得到轮廓清晰的图。
页: 1 2 [3]
查看完整版本: 看图说话——布鲁塞尔问题的分叉图