N体仿真缩放与误差容限适配方法咨询

一、误差容限的缩放逻辑 #

首先明确Matlab ODE45的误差控制规则：求解器会确保每一步的局部误差不超过 max(relTol * |y|, absTol)，其中y是包含粒子位置、速度的状态向量。这里的threshold = absTol / relTol是关键阈值——当|y| > threshold时，相对误差容限relTol主导误差控制；当|y| < threshold时，绝对误差容限absTol主导。

针对你的缩放场景（比如将原系统尺寸从10缩放到1，对应状态向量y的缩放因子k=10，即y' = y/k），要保证缩放前后误差控制等效，需遵循以下原则：

1. 相对误差容限relTol无需调整：相对误差是相对于状态值的比例，缩放后|y'| = |y|/k，relTol * |y'| = relTol * |y|/k，而缩放后的局部误差e' = e/k（与状态向量同比例缩放），因此原系统的相对误差控制逻辑在缩放后依然成立，无需修改relTol。
2. 绝对误差容限absTol需按缩放因子反比例调整：原系统中absTol控制绝对误差上限e ≤ absTol，缩放后绝对误差变为e' = e/k，因此缩放后的absTol'应设为absTol / k，这样e' ≤ absTol'等价于e ≤ k*absTol' = absTol，和原系统的绝对误差控制一致。

举个例子，原系统relTol=1e-3、absTol=1e-6（threshold=1e-3），缩放因子k=10时，应设置relTol=1e-3、absTol=1e-7。这样既保持了原系统的误差控制强度，又不会因缩放导致absTol相对过大/过小，避免求解器陷入极小步长。

你提到的“同比例调整容限保持threshold不变”（比如relTol=1e-4、absTol=1e-7），会让误差控制更严格——relTol*|y'| = 1e-4*(|y|/10) = 1e-5|y|，远小于原系统的1e-3|y|，反而迫使求解器用更小的步长，加重你遇到的问题，并不推荐。

二、仅使用单一误差容限的合理性 #

在N体仿真场景下，仅依赖relTol或absTol都存在明显缺陷：

• 仅用relTol（设absTol为极大值）：当粒子靠近时，位置或相对速度的量级会变得极小（比如从10降到1e-2），此时relTol*|y|会变得非常小（比如1e-3*1e-2=1e-5），求解器会为满足这个极小的误差上限而选择过小的步长，导致仿真效率骤降。
• 仅用absTol（设relTol为极大值）：当粒子处于大尺度运动状态时，|y|很大（比如10），absTol相对于|y|可以忽略，此时误差控制完全由absTol主导，相当于强制固定绝对误差上限，会导致大尺度运动的相对误差过大，仿真结果失真。

因此，必须同时使用relTol和absTol，两者配合才能覆盖N体系统中状态向量量级变化的所有场景——relTol保证大状态值下的相对精度，absTol保证小状态值下的绝对精度。

三、关于“将relTol和absTol均设为1”的问题 #

这个设置完全不合理，它相当于允许每一步的局部误差达到状态值本身或1，会导致仿真结果严重偏离真实动力学行为，只是暂时规避了极小步长的问题，本质上是牺牲所有精度换取速度，完全无法保证缩放前后系统行为的一致性。

内容的提问来源于stack exchange，提问作者Paul Aner