我正在模拟一个通过两个针孔的不同电荷状态的离子束。我使用来自外部函数的随机角度来评估沿轴的速度分量。 为了评估粒子运动，我正在求解微分方程：
O=[0 0 0]';
f = @(t,ys) [ys(4:6); O];
options=odeset('RelTol',1e-7,'Events',@(t,ys)Event_Stop_Sorgente(t,ys,Pinhole1));
[t,ys] = ode23t(f,tspan,y_sorgente(i,:),options);
只要光束没有角度扩展且初始条件向量为y_sorgente = [0 0 0 0 0 v_z]，它就可以正常工作。我认为它也适用于角展度，因为轨迹是直线，正如它们应该的那样。
外部函数从一个sferic曲面得到一些点，并使用它们的坐标作为方向余弦来从模数中获得速度分量。
我遇到的问题是，只有少数电荷状态为1+的粒子可以通过针孔，其他所有粒子都不能，并且没有任何意义。即使每个粒子具有相同的速度模量，我也有同样的问题。我做错了什么，但我无法理解什么，有什么建议吗？ 如果你想检查一下，我可以提供完整的代码。