向各位大侠求助庞加莱投影图的画法

煜宸0922

回复 14 # meiyongyuandeze 的帖子

我知道是选取碰撞面和碰撞后的速度，因为这个在罗老师和一些论文上有说，但具体在我那个程序上实现是我做的总有问题 ，所以才求助，对了，老师，怎样保证是正向穿越呢？

煜宸0922

回复 14 # meiyongyuandeze 的帖子

大侠，我用庞加莱截面法做的图，感觉程序已经对了，而且满足了速度位移的要求，可就是图不对，我到底什么地方出错了，我实在想不通，贴个程序，请大侠指导下
global w
w=3.172;
[t,x,V,P]=changerk144(80000);
ind=find((x(1,:)-0.1)<=1e-4&x(1,:)>0);
V1_penghou=P;
X2 =[];  V2=[];
    for j=1:length(ind)
    if  abs(x(3,ind(j))-V1_penghou)<=1e-4&x(3,ind(j))<0
        X2=[X2,x(2,j)];  
        V2=[V2,x(4,j)];
    end
    end
hold on
plot(X2,V2,'.k')