我知道是选取碰撞面和碰撞后的速度，因为这个在罗老师和一些论文上有说，但具体在我那个程序上实现是我做的总有问题 ，所以才求助，对了，老师，怎样保证是正向穿越呢？ 回复 14 # meiyongyuandeze 的帖子
大侠，我用庞加莱截面法做的图，感觉程序已经对了，而且满足了速度位移的要求，可就是图不对，我到底什么地方出错了，我实在想不通，贴个程序，请大侠指导下
global w
w=3.172;
=changerk144(80000);
ind=find((x(1,:)-0.1)<=1e-4&x(1,:)>0);
V1_penghou=P;
X2 =[];V2=[];
for j=1:length(ind)
ifabs(x(3,ind(j))-V1_penghou)<=1e-4&x(3,ind(j))<0
X2=;
V2=;
end
end
hold on
plot(X2,V2,'.k')
页:
1
[2]