快捷導(dǎo)航

Python光學(xué)仿真學(xué)習(xí)衍射算法初步理解

更新時間：2021年10月20日 11:21:19 作者：微小冷

這篇文章主要為大家介紹了Python光學(xué)仿真學(xué)習(xí)中對衍射算法的初步理解，有需要的朋友可以借鑒參考下，希望能夠有所幫助，祝大家多多進(jìn)步

對衍射最經(jīng)典的解釋是Huygens-Fresnel原理，Huygens認(rèn)為波陣面上每一點都會成為新的波源，這些子波源的相互干涉就形成了衍射。這顯然是一種離散的觀點，仿佛是專門為程序員準(zhǔn)備的一樣。

假設(shè)一束光打在一個方形孔上，這個方形孔被細(xì)分成 n×n個網(wǎng)格，那么每個網(wǎng)格都相當(dāng)于是一個小孔，而這些小孔的互相干涉，即為衍射。隨著網(wǎng)格不斷被細(xì)分，最終可以逼近真實的衍射情形。那么，假設(shè)矩孔處為等相位面，其網(wǎng)格坐標(biāo)為 (i,j)，到衍射屏距離為 d d d,那么對于衍射屏上任意一點P(x,y)，其光強為

在這里插入圖片描述

同理，我們剛剛寫下的平面矩孔光強疊加也出現(xiàn)了問題，如果我們默認(rèn)矩孔上每個格點都是一個點光源，那么打在衍射屏上之后，應(yīng)該遵從球面波的衰減原則。又因為這種假設(shè)其實忽略了從光源射到矩孔過程中的光線的傳播方向，所以應(yīng)該有一個傾斜因子，即 ( i , j ) 點打在 ( x , y )點的光強為

在這里插入圖片描述

#基爾霍夫衍射,衍射屏坐標(biāo)范圍-dGrid:dGrid,光源坐標(biāo)(0,0)
#簡單的矩孔衍射,dSource為光源到小孔的距離；dScreen為衍射屏到小孔距離
#dHole為矩孔網(wǎng)格尺寸；dGrid為衍射屏網(wǎng)格尺寸；nGrid為網(wǎng)格數(shù)目
def squareDiff(dSource=1,dScreen=1.4,dWave=1e-6,
            dHole=3e-5,nGrid=100,dGrid=1e-5):
    nX,nY = nGrid*np.array([1,1])
    axisX = np.arange(-nX,nX+1)*dGrid       
    axisY = np.arange(-nY,nY+1)*dGrid
    xAxis,yAxis = np.meshgrid(axisX,axisY)  #此為衍射屏的x坐標(biāo)
    axisX = np.arange(-nX,nX+1)*dHole       
    axisY = np.arange(-nY,nY+1)*dHole
    xHole,yHole = np.meshgrid(axisX,axisY)  #此為矩孔的x坐標(biāo)
    dArrS = np.sqrt(xHole**2+yHole**2+dSource**2)   #孔平面到光源的距離
    nSide = int(nGrid*2+1)              #格點個數(shù)
    pane = np.zeros([nSide,nSide])      #衍射屏強度
    for m in range(nSide):
        for n in range(nSide):
            dArr = np.sqrt((xHole-xAxis[m,n])**2+(yHole-yAxis[m,n])**2+dScreen**2)
            pane[m,n] = np.sum(
                np.cos(np.pi*2*(dArr+dArrS)/dWave)/dArr/dArrS)
    pane = np.abs(pane)
    pane = pane/np.max(np.array(pane))
    fig = plt.figure()
    ax = axd(fig)
    ax.plot_surface(xAxis,yAxis,pane)
    plt.show()
    return pane, xAxis, yAxis

其結(jié)果為

在這里插入圖片描述

故可定義矩陣索引

#輸入對于M*M矩陣的第一個值到N*N矩陣的距離，返回M(m,n)的距離矩陣
def getDisMat(dMat,N,m,n):
    dMat = np.mat(dMat)
    A = dMat[1:m,1:n]
    B = dMat[1:m,0:N-n+1]
    C = dMat[0:N-m+1,1:n]
    D = dMat[0:N-m+1,0:N-n+1]
    return np.vstack((np.hstack((np.flip(A),np.flipud(B))),np.hstack((np.fliplr(C),D))))    #stack矩陣拼接，flip翻轉(zhuǎn)矩陣

相應(yīng)地算法改為(其他位置不變)

    dArrS = np.sqrt(xHole**2+yHole**2+dSource**2) #孔平面到光源的距離
    dScreen = np.sqrt(                             #衍射平上第(0,0)個點的距離矩陣
        (xHole-xAxis[0,0])**2+(yHole-yAxis[0,0])**2+dScreen**2)
    nSide = int(nGrid*2+1)              #格點個數(shù)
    pane = np.zeros([nSide,nSide])      #衍射屏強度
    for m in range(nSide):
        for n in range(nSide):
            dArr = getDisMat(dScreen,nSide,m,n)

但這種優(yōu)化是極其有限的，這是直觀無腦的思維方式所帶來的麻煩。在接下來的傅里葉光學(xué)中，我們將繼續(xù)處理衍射的計算問題。

以上就是Python光學(xué)仿真學(xué)習(xí)衍射算法初步理解的詳細(xì)內(nèi)容，更多關(guān)于Python光學(xué)仿真衍射算法的資料請關(guān)注腳本之家其它相關(guān)文章！

您可能感興趣的文章: