我正在尝试绘制具有可变入射角(0 到 180 度)的垂直波和平行波的相移。波从折射率为 1.33 传播到折射率为 1.5 的介质。
我使用了以下方程:方程理论<-- 第 18 页
我使用了以下代码:
def Phase(theta):
n=1.5/1.33
Shift=np.sqrt(np.sin(theta*np.pi/180)**2-n**2)
Shift=Shift/np.cos(theta*np.pi/180)
Shift=2*np.degrees(np.arctan(Shift))
return Shift
print(Phase(x))
x=np.linspace(0,180,30)
问题是我得到[ nan nan nan nan nan nan nan nan nan nan]
了回报。