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

频闪法庞加莱截面

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


程序代码：function fga
y0=[-0.05;0;0];
T=4*pi;
ts=;
=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=;
else
   dx=;
end上图分别为庞加莱截面图和相图，始终对不上，请高手指点


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

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

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

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

lihaitao123 发表于 2013-7-17 12:03 static/image/common/back.gif
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

就是不知道那个截面怎么取？
页: [1]
查看完整版本: 频闪法庞加莱截面