luck0110 发表于 2014-11-28 17:07

如何绘制庞加莱截面图

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

function dx=mm_4(t,x)
dx=;

然后又编写了画制庞加莱截面图的M函数。
function f
tt=2*pi/1.2;
=ode45(@mm_4,,);
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 发表于 2014-12-1 20:30

{:{10}:}

thdl520 发表于 2014-12-3 10:06

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

haomeng100 发表于 2014-12-6 16:38

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

Rotor2014 发表于 2014-12-6 19:56

解决了吗？

dj15175771160 发表于 2014-12-9 14:13

每个周期截取一个点吧
页: [1]
查看完整版本: 如何绘制庞加莱截面图