如何编写 poincare 的截面图程序

zhailiangjun

我這里有一個4個變量的微分方程组zhzfun，变量为y（1），y（2），y（3），y（4），我要画出
X=(y（1）,y（2）)在y(3)=0截面处的poincare截面图。程序应该怎么写啊，请各位高手帮帮忙。谢谢了。

[ 本帖最后由 eight 于 2008-1-21 18:57 编辑 ]

无水1324

%    |x-b    when x>b
% y= | x     when -b<x<b
%    |x+b    when x<-b
clc;
b=1;
ff=inline('[y(2);-(y(1)>b)*(y(1)-b)-(y(1)<-b)*(y(1)+b)-(-b<y(1)&y(1)<b)*y(1)]',...
    'x','y','flag','b');
[t,y]=ode45(ff,[0,40],[2,0],[],b);
close all;
subplot(131);plot(t,y(:,1),'k',t,y(:,2),'k:','markersize',2);
subplot(132);plot(y(:,1),y(:,2),'k');title('相图');xx=xlim;
q=find(abs(y(:,1))<1e-1);
subplot(133);plot(y(q,1),y(q,2),'k.','markersize',6);xlim(xx);
title('庞克来相图')
set(gcf,'Position',[16 135 995 327])


此程序转自一个论坛，具体是哪个我不知道了，感谢原创者，发在这里对楼主应该有帮助。，还有你的短信我也看到了

zhailiangjun

谢谢，无水大哥