如何编写 poincare 的截面图程序
我這里有一個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 编辑 ]
回复 楼主 的帖子
% |x-b when x>b% y= | x when -b<x<b
% |x+b when x<-b
clc;
b=1;
ff=inline('',...
'x','y','flag','b');
=ode45(ff,,,[],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',)
此程序转自一个论坛，具体是哪个我不知道了，感谢原创者，发在这里对楼主应该有帮助。，还有你的短信我也看到了
回复 2楼 的帖子
谢谢，无水大哥
页:
[1]