Question

所以我正在一个项目工作，我需要从服务器检索粒子数据（可能会有几千个对象）（为了同步客户端），我正在努力寻找如何计算多个粒子的轨道倾角。

到目前为止，我有轨道总半径，宽度（粒子的起点和终点）

 var asteroidsOrbit = Common.random( 90, 100),
     rotation = Math.PI / 2, // Not sure where to use this yet
     radians = Common.random(0, 360) * Math.PI / 180;

我也通过设置粒子在所需位置的位置来渲染轨道：

     position_x: Math.sin(radians) * asteroidsOrbit,
     position_y: Common.random(-10, 10),
     position_z: Math.cos(radians) * asteroidsOrbit

我注意到默认情况下，轨道垂直于X平面，表现得如此。但是我应该如何应用rotation，因为它可能在x平面上旋转20度？

Answer 1

我在Python中制作了一些代码，将旋转应用于粒子在3D空间中的速度，可能会有所帮助：

def make_rotation_matrix(theta=0.,phi=0.,cos_theta=False,cos_phi=False):
    """
    Returns a rotation matrix for given angles.
    Makes the rotation matrix for a given theta angle around the y axis (polar
    angle) followed by another rotation around the z axis by a given phi angle
    (azimuthal angle).

    Parameters
    ----------
    theta : float
        Polar angle, between 0 and pi (or its cosine between -1 and 1).
    phi : float
        Azimuthal angle, between 0 and 2pi.
    cos_theta : boolean
        Wether to consider theta as an angle (False) or its cosine (True).
    cos_phi : boolean
        Wether to consider phi as an angle (False) or its cosine (True).

    Returns
    ----------
    out : np.array
        Rotation matrix (shape == (3,3)).

    Examples
    --------
    >>> make_rotation_matrix(0,10)
    array([[-0.83907153,  0.54402111, -0.        ],
           [-0.54402111, -0.83907153, -0.        ],
           [-0.        ,  0.        ,  1.        ]])
    >>> make_rotation_matrix(10,0)
    array([[-0.83907153, -0.        , -0.54402111],
           [-0.        ,  1.        , -0.        ],
           [ 0.54402111,  0.        , -0.83907153]])
    >>> make_rotation_matrix(0.1,0,True)
    array([[ 0.1       , -0.        ,  0.99498744],
           [ 0.        ,  1.        ,  0.        ],
           [-0.99498744,  0.        ,  0.1       ]])
    """
    if cos_theta:
        cos_theta = theta
        sin_theta = np.sqrt(1. - cos_theta**2)
    else:
        cos_theta = np.cos(theta)
        sin_theta = np.sin(theta)
    if cos_phi:
        cos_phi = phi
        sin_phi = np.sqrt(1. - cos_phi**2)
    else:
        cos_phi = np.cos(phi)
        sin_phi = np.sin(phi)
    return np.array(
    [[cos_theta*cos_phi ,  -sin_phi  ,  sin_theta*cos_phi],
     [cos_theta*sin_phi ,   cos_phi  ,  sin_theta*sin_phi],
     [     -sin_theta   ,      0     ,  cos_theta        ]])


def skew_symmetric_cross_product(vector):
    """
    Returns the skew symmetric cross product of given vector.

    Parameters
    ----------
    vector : iterable
        Iterable object, must have at least 3 elements, corresponding to
        the vectors x, y and z components.

    Returns
    ----------
    out : np.array
        Skew symmetric cross product of the given vector (shape == (3,3)).

    Examples
    --------
    >>> skew_symmetric_cross_product([1,0,0])
    array([[ 0.,  0.,  0.],
           [ 0.,  0., -1.],
           [ 0.,  1.,  0.]])
    >>> skew_symmetric_cross_product([0,1,0])
    array([[ 0.,  0.,  1.],
           [ 0.,  0.,  0.],
           [-1.,  0.,  0.]])
    >>> skew_symmetric_cross_product([0,0,1])
    array([[ 0., -1.,  0.],
           [ 1.,  0.,  0.],
           [ 0.,  0.,  0.]])
    """
    return np.array([[ 0.         , -vector[2] , vector[1]  ],
                     [ vector[2]  , 0.         , -vector[0] ],
                     [ -vector[1] , vector[0]  , 0.         ]])


def rotation_a2b(a,b):
    """
    Returns the rotation matrix that transforms vector a in b.
    From
    http://math.stackexchange.com/questions/180418/

    Parameters
    ----------
    a : iterable
        Iterable object, must have at least 3 elements, corresponding to
        the vectors x, y and z components.

    b : iterable
        Similar to a.

    Returns
    ----------
    out : np.array
        Rotation matrix that transforms a in b (shape == (3,3)).

    Examples
    --------
    >>> rotation_a2b([1,0,0],[1,0,0])
    array([[ 1.,  0.,  0.],
           [ 0.,  1.,  0.],
           [ 0.,  0.,  1.]])
    >>> rotation_a2b([1,0,0],[0,1,0])
    array([[ 0., -1.,  0.],
           [ 1.,  0.,  0.],
           [ 0.,  0.,  1.]])
    >>> rotation_a2b([1,0,0],[0,0,1])
    array([[ 0.,  0., -1.],
           [ 0.,  1.,  0.],
           [ 1.,  0.,  0.]])
    """
    if np.all(a==b):
        return np.identity(3)
    v = np.cross(a,b)
    v_sscp = skew_symmetric_cross_product(v)
    s = np.sqrt(np.dot(v,v))
    c = np.dot(a,b)
    return np.identity(3) + v_sscp +\
           (((1-c)/(s**2))*np.dot(v_sscp,v_sscp))

摘自[编辑]来自我的报告：

由于旋转位于粒子的速度（vp）附近，因此R生成的矩阵make_rotation_marix无法直接应用于vp - 我们必须先更改使速度与（0,0,1）向量共线的基础，应用旋转，然后将基础更改回原始。由于第一次基础变化的结果是（0,0，sp），sp是粒子的速度（velocitie的幅度），所以跳过第一个基础变化并且R适用于（0,0，sp）。

为了改变基础，我们使用矩阵，它将（0,0,1）变回初始粒子速度，由函数Rotation_a2b提供。矩阵Ra2b将向量a转换为b，因此调用a =（0,0,1）和b等于标准化粒子初始速度。我们对速度进行标准化的事实意味着基础的变化保留了向量的大小 - 或者可以使用a =（0,0，sp）。

三维空间中的轨道倾角

1 个答案: