phase函数计算方向场