我需要制作一个imshow()图表,每个元素显示极角phi的值(从x轴开始)。我想到了以下代码:
x = np.linspace(-3*a,3*a,1000, dtype='complex')
y = np.linspace(-3*a,3*a,1000, dtype='complex')
X,Y = np.meshgrid(x,y)
rho = np.sqrt(X**2+Y**2)
phi = np.arcsin(Y/rho)
最后一行的问题是,如果我只是天真地使用arcsin,arccos或arctan,那么在某些象限中符号将是不正确的。我需要phi从0到2 * pi,或从-pi到+ pi,但是使用上面的代码,所有值中有一半的符号不正确。
那我该怎么做呢?