%具体的讨论可以见我的个人空间！
%这是duffing方程的Poincare截面的程序，经过试验和调试是没有错误的！
%这是定义的duffing方程的式子！
%要调用poincare.m!
function dy=duffing(t,y);
global det
global w
global r
dy=[y(2);y(1)-y(1)^3-det*y(2)-r*cos(w*t)];




%Duffing 方程的Poincare截面程序！
%经过试验和调试是没有错误的！
%要调用duffing.m!
function  poincare
global det
global r
global w
det=0.25;w=1;
r=0.3;
x0=[0;0];
tspan=[0:0.1*2*pi:50000];
[t,y]=ode45('duffing',tspan,x0);
m=zeros(7000,1);
n=zeros(7000,1);
for i=1:7000
    m(i,1)=y(10*i,1);
    n(i,1)=y(10*i,2);
end
plot(m,n,'k.','markersize',1);

[ 本帖最后由 liliangbiao 于 2008-4-24 15:46 编辑 ]

untitled.PNG