作为通过光圈模拟衍射的长程序的一部分,我正在摆弄试图让我的结果plt.imshow()
的原点位于绘图的中心,即我希望改变轴
相关的代码部分是:
n=40
lam=0.0006
k=(2*np.pi)/lam
z=float(input("Type screen distance, z in mm: "))
rho=float(input("Type rho in mm: "))
delta = 2/n
intensity = np.zeros((n+1,n+1))
for i in range(n+1):
x=-1+(i*delta)
for j in range(n+1):
y =-1+(j*delta)
intensity[i,j] = (abs(square_2dsimpson_eval(-rho/2,rho/2,n)))**2
plt.imshow(intensity)
plt.show()