小球间完全弹性斜碰（赋Python代码）

完全弹性斜碰：

本文讨论两个考虑质量的小球间完全弹性斜碰。

完全弹性正碰可直接通过动能定理和动量守恒即解决。

完全弹性斜碰稍微复杂一点，具体计算过程如下：

首先拿两小球记为A球和B球，A球B球属性有质量，速度大小，速度方向。

建立坐标系（后面会用到）：

        1）以A球圆心为原点，以水平方向向右为X轴，逆时针偏转PI/2为Y轴，建立A坐标系。

        2）以B球圆心为原点，以水平方向向右为X轴，逆时针偏装PI/2为Y轴，建立B坐标系。

        3）以A球B球圆心连线为X轴，方向由A到B，逆时针偏装PI/2为Y轴，建立AB坐标系。

      

碰撞过程中涉及到的变量：（下文会逐一介绍）

    angleA        VA          VAX         VAY        mA      angleAB

    angleNA      VNA       VNAX      VNAY

    anglePNA    VPNA    VPNAX    VPNAY

    angleZA       VZA       VZAX      VZAY

    angleB        VB          VBX         VBY        mB      angleAB

    angleNB      VNB       VNBX      VNBY

    anglePNB    VPNB    VPNBX    VPNBY

    angleZB       VZB       VZBX      VZBY

碰撞的计算过程：

已知：A球和B球的初速度和方向。

          angleA，VA，angleB，VB。

第零步：给变量改名。

第一步：求XAB的角度，即B相对于A的角度。

第二步：记录XAB坐标系中VB的角度，并求出XAB坐标系中X轴Y轴的分速度。

第三步：XAB坐标系中，Y轴速度不变，X轴速度满足完全弹性正碰。

第四步：合成新速度，并计算XAB中的角度。

第五步：转化为水平坐标系中的速度和角度。

第六步：赋值操作。

第零步：给变量改名（变量多改名便于记忆）

    A球的运动方向记为angleA，速度大小记为VA。

    A球在坐标系A中X轴的角度记为VAX = VA * cos（angleA）

    A球在坐标系A中Y轴的角度记为VAY = VA * sin（angleA）

    同理可得B球的angleB，VB，VBX，VBY。

第一步：求XAB坐标系X轴与水平方向的角度angleAB，并求出A球B球在XAB坐标系中的速度大小和方向

    获得A球B球的圆心实际坐标p