频闪法庞加莱截面

飞翼 · 发表于 2013-7-16 20:28

您需要登录才可以下载或查看，没有账号？我要加入

x

本帖最后由牛小贱于 2014-7-8 10:01 编辑

程序代码：

function fga
y0=[-0.05;0;0];
T=4*pi;
ts=[0:T/1000:500*T];
[t,x]=ode45(@dsg,ts,y0);
figure(1)
plot(x(40000:end,1),x(40000:end,2))
figure(2)
plot(x(40000:1000:end,1),x(40000:1000:end,2))
function dx=dsg(t,x)
u=0.0213;
a=20;
B=0.5;
w=0.5;
g=0.54;
if x(1)<0
dx=[x(2);-2*u*x(2)-(1+108/35*g^2)*x(1)-108/35*g*x(1)^2-36/35*x(1)^3-g-36/35*g^3+B*cos(x(3));w];
else
dx=[x(2);-2*u*x(2)-(1+108/35*g^2)*x(1)-108/35*g*x(1)^2-36/35*x(1)^3-g-36/35*g^3+B*cos(x(3))-a*x(1);w];
end

复制代码

上图分别为庞加莱截面图和相图，始终对不上，请高手指点

lihaitao123 · 发表于 2013-7-17 12:02

本帖最后由牛小贱于 2014-7-8 10:02 编辑

复制代码

飞翼 · 发表于 2013-7-17 20:02

lihaitao123 发表于 2013-7-17 12:03

登录/注册后可看大图

plot(x(40000:1000:end,1),x(40000:1000:end,2),'k.')

确实行了，但是不知道为啥？

zhimu · 发表于 2014-7-7 23:38

你这好像不是庞加莱截面图啊

Rotor2014 · 发表于 2014-12-22 11:34

zhimu 发表于 2014-7-7 23:38
你这好像不是庞加莱截面图啊

上面修改后的对吗

张某某zyl · 发表于 2015-3-4 20:30

就是不知道那个截面怎么取？

账号		自动登录	找回密码
密码			我要加入

[编程技巧] 频闪法庞加莱截面