Z=[]; [T,Y]=ode45('Chen',[0,5000],[1;1;1;35;3;28]); for k=1:length(Y) if abs(Y(k,1))<1e-2 Z=[Z Y(k,2)+i*Y(k,3)]; end end plot(Z,'.','markersize',2) title('Chen 系统的 Poincare 映像 x=0') xlabel('y'),ylabel('z')
Z=[]; [T,Y]=ode45('Chen',[0,5000],[1;1;1;35;3;28]); for k=1:length(Y) if abs(Y(k,3)-28)<1e-2 Z=[Z Y(k,1)+i*Y(k,2)]; end end plot(Z,'.','markersize',2) title('Chen 系统的 Poincare 映像 z=c') xlabel('x'),ylabel('y')
5 comments:
建议误差设定值要足够小，否则失真：）
不知是否同意！
http://complex.blog.163.com/
你可以试试，我取过小的误差，但得到的点太少以至于看不出什么来，要得到更多的点只有加大迭代的步数，但运行挺费时间的，为快点得结果，感觉只要趋势差不多也算过得去了，所以误差随便取大了些.
:)
以前作过这些东西，看样子你误解了“截面”的面的含义，严格讲不能太“厚”，呵呵。
我的Blog搬家了，不介意我加了你的连接吧：
http://
www.blog.163.com/mail2liujie/
另外，取截面有讲究，比如不稳定的平衡点处，。。。。。。位置不对当然点就反映不出特征来。。。。。。
-----------------------
jerry
后面我想了下，应该修改下程序，不过现在忙于考研只好等到周五周六晚上不熄灯的时候去弄了。
另外当然不介意你加了我的连接，应该表示感谢支持才行。
Post a Comment