24小时热门版块排行榜

返回列表

luck0110

金虫 (初入文坛)

应助: 0 (幼儿园)
金币: 745.3
帖子: 28
在线: 23.7小时
虫号: 1992625
注册: 2012-09-11
性别: GG
专业: 机械测试理论与技术

[求助] 如何绘制庞加莱截面图已有1人参与

照着别人的思路，绘制了达芬方程的仿真图，但是感觉不对；
程序如下：
首先建立了存放函数的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，相图为一个封闭图形，庞加莱界面应该是一个点才对；不知道错在哪个地方？
拜托高手帮忙了~谢谢。

1.JPG