如何绘制庞加莱截面图

luck0110

照着别人的思路，绘制了达芬方程的仿真图，但是感觉不对；
程序如下：
首先建立了存放函数的M函数；

function dx=mm_4(t,x)
dx=[x(2);cos(t)-x(1)^3+x(1)-4*x(2)];

然后又编写了画制庞加莱截面图的M函数。
function f
tt=2*pi/1.2;
[t,x]=ode45(@mm_4,[0:tt/100:80*tt],[1,1]);
figure
subplot(2,1,1)
plot(t,x(:,1));
title('位移曲线');
xlabel('t');ylabel('x');

subplot(2,2,3)
plot(x(2000:end,1),x(2000:end,2),'-')
xlabel('x');ylabel('dx/dt');
title('相图,r=4,w=1,f=1');

subplot(2,2,4)
i=2000:100:3000 ;
plot(x(i,1),x(i,2),'*')
title('庞加莱截面');


通过修改r的值，得到不同的图形；


通过图片可以看出，r=4，相图为一个封闭图形，庞加莱界面应该是一个点才对；不知道错在哪个地方？
拜托高手帮忙了~谢谢。

luck0110

thdl520

论坛里面的大牛都不知哪里去了  我发了个帖子也是无人问津啊   

haomeng100

貌似画庞加莱截面时，舍去的点太少。

Rotor2014

解决了吗？

dj15175771160

每个周期截取一个点吧