在下面的代码中，我计算了旋转速度作为半径（V_r）的函数，用于xy坐标的笛卡尔网格，将xy坐标转换为极坐标r。
然后我通过取arctan（y / x）来计算极坐标phi，然后使用cos（phi）计算y中的V_r分量* V_r = V_y：
res = 100
x = np.linspace(-1.2,1.2,res)
y = np.linspace(-1.2,1.2,res)
r = np.sqrt(x**2 + y**2)
phi = np.arctan2(y, x)
cosphi = abs(np.cos(phi))
# Velocity Const.
V0 = 200.
# Constants for scaling r
rpe = 0.164
alpha = 0.002
V_r = V0 * (1 - np.exp(-r/rpe)) * (1 + (alpha * r / rpe))
V_los = V_r * abs(np.cos(phi))
plt.figure()
plt.plot(r, V_los)
无需实际创建数组x和y，转换为r和phi。我知道我可以简单地输入极坐标。
然而，我努力达到相同的结果，即在phi = 0（即cos（phi）= 1）时V_r = V_y，或cos（phi）= 0，因此V_y = 0。
如何输入数组r和phi以达到相同的效果？