奇怪復相位顯示

Question

使用以下代碼：

from numpy import * 
from matplotlib.pyplot import * 

Radius=10 
N=1024 
dx=2*Radius/N 
dy=dx 
x=r_[-Radius:Radius:dx] 
y=r_[-Radius:Radius:dy] 
X, Y = meshgrid(x,y) 
R = sqrt(X**2+Y**2) 
PHI = arctan2(Y,X) 

ringthing = R < Radius 
ring = zeros((2,N,N),dtype=complex) 
ring[0] = ringthing 
ring[1] = ringthing*exp(1j*PHI) 

f=fig() 
p1=f.add_subplot(121) 
p1.imshow(angle(ring[0])) 
p2=f.add_subplot(122) 
p2.imshow(angle(ring[1])) 

f.show() 

第二圖像的左下方被標記爲紅色（相等於PI）沒有明顯的理由。爲什麼是這樣？

Answer 1

問題是圓外的值是零，復角的零不明確（它是一個奇點）。在浮點算術鍛鍊，使得一些地區，他們評估爲0和他人-0這是通過運行

from __future__ import division 

from numpy import * 

Radius=10 
N=1024 
dx=2*Radius/N 
dy=dx 
x=r_[-Radius:Radius:dx] 
y=r_[-Radius:Radius:dy] 
X, Y = meshgrid(x,y) 
R = sqrt(X**2+Y**2) 
PHI = arctan2(Y,X) 

ringthing = R < Radius 
ring = zeros((2,N,N),dtype=complex) 
ring[0] = ringthing 
ring[1] = ringthing*exp(1j*PHI) 

print ring[1][-1, 0], angle(ring[1][-1, 0]) 
print ring[1][0, -1], angle(ring[1][0, -1]) 

與輸出

(-0+0j) 3.14159265359 
0j 0.0 

一個解決方案解決此是設置可見圓圈外的所有值都明確爲零。