方慧雯，楊錦宏，賀勝男，衛(wèi)玉嬌，張吳記，汪衛(wèi)華

（安徽大學(xué) 物質(zhì)科學(xué)與信息技術(shù)研究院，安徽 合肥 230601）

隨著技術(shù)的革新，光的衍射現(xiàn)象在生產(chǎn)實(shí)踐中的應(yīng)用越來(lái)越廣泛.多邊形的衍射模擬在光學(xué)領(lǐng)域具有重大意義，可以解釋燈光星芒產(chǎn)生的原因、相機(jī)如何記錄圖樣［1］以及如何在夜晚拍出明顯星芒［2］.光的衍射現(xiàn)象可以用惠更斯-菲涅耳原理和瑞利-索末菲公式進(jìn)行計(jì)算，但數(shù)學(xué)公式計(jì)算繁瑣，不夠直觀(guān).計(jì)算機(jī)模擬則能通過(guò)顯示圖像，直觀(guān)地觀(guān)察衍射效果.且計(jì)算機(jī)模擬不受實(shí)驗(yàn)場(chǎng)地和實(shí)驗(yàn)器材的限制.

目前，有很多學(xué)者已經(jīng)對(duì)不同孔徑的衍射現(xiàn)象進(jìn)行了研究:比如劉有菊等研究了“E”形孔的夫瑯禾費(fèi)衍射現(xiàn)象［3］，宋易知模擬了任意正多邊形小孔的夫瑯禾費(fèi)衍射現(xiàn)象［4］，王海濤等模擬了任意多邊形小孔的菲涅耳衍射［5］，Jeffrey等利用快速菲涅耳衍射算法模擬了孔形為“2”的衍射［6］，Qian等研究了菲涅耳和夫瑯禾費(fèi)衍射的數(shù)值模擬［7］，汪峰等研究了任意形狀孔夫瑯禾費(fèi)計(jì)算機(jī)模擬［8］，戴又善研究了二維小孔夫瑯禾費(fèi)衍射的對(duì)稱(chēng)變換特性［9］.但已有的研究主要針對(duì)的是正多邊形，不能適用于任意亮度孔形，并且現(xiàn)有的數(shù)值模擬不能適用任意距離.由于菲涅耳衍射和夫瑯禾費(fèi)衍射特性差別顯著，需根據(jù)菲涅耳數(shù)進(jìn)行劃分［10］.在上述研究的啟發(fā)下，本文的創(chuàng)新點(diǎn)在于用圖片像素灰度值來(lái)代表小孔內(nèi)不同的亮度值；研究了任意亮度的小孔衍射；并且該程序適用于菲涅耳和夫瑯禾費(fèi)衍射兩種情況.即采用Matlab軟件，利用圖片像素讀取任意亮度孔形，然后對(duì)近場(chǎng)/遠(yuǎn)場(chǎng)衍射進(jìn)行數(shù)值模擬，給出任意形狀、任意亮度孔的衍射圖案來(lái)探尋其中的衍射規(guī)律.

1 衍射的理論與數(shù)值模型

1.1 衍射理論

光的衍射現(xiàn)象是指光在傳播過(guò)程中，遇到障礙物或小孔時(shí)，光將偏離直線(xiàn)傳播的路徑而繞到障礙物后面?zhèn)鞑サ默F(xiàn)象［11］.把波面 S上每個(gè)面積元 ds看作新的波源，如圖1所示，它們均發(fā)出次波，波面前方空間某一點(diǎn)P的振動(dòng)可以由S面上所有面積元所發(fā)出的次波在該點(diǎn)疊加后的合振幅來(lái)表示［12］.在本文中，將小孔均分成一定數(shù)量的方格（即圖像像素），像素灰度值表示孔的亮度.光屏也均分為一定數(shù)量的方格（近似將其看作一點(diǎn)）.光屏上某一點(diǎn)的振動(dòng)可看成小孔內(nèi)所有方格點(diǎn)在該點(diǎn)的合振幅.將光屏上所有方格點(diǎn)的振動(dòng)計(jì)算出來(lái)，就可以顯示該孔的衍射圖像.

圖1 衍射光路示意圖

設(shè)有一個(gè)波長(zhǎng)為λ的平面波入射，衍射平面與光屏之間的距離（焦距）為 f，衍射平面上點(diǎn) P（x1，y1）與光屏上點(diǎn) Q（x2，y2）之間的距離為 r，根據(jù)惠更斯-菲涅耳原理復(fù)數(shù)形式:

當(dāng)忽略時(shí)間t時(shí)則（2）可以簡(jiǎn)化為

1.2 衍射數(shù)值模型

首先用繪圖軟件繪制尺寸合適的任意孔形，然后用imread函數(shù)讀出所繪制的孔形，在程序中將小孔的部分設(shè)為像素灰度值，其余部分設(shè)為0，進(jìn)行讀取，編寫(xiě)程序.此程序適用于任意孔徑、任意亮度、任意距離的衍射，對(duì)于不同的孔形或不同的亮度只需調(diào)換相應(yīng)的孔形圖片；對(duì)于不同距離的衍射，只需更改焦距f.在此基礎(chǔ)上，調(diào)換合適的光屏大小，孔徑大小，得到衍射明顯的衍射圖像.計(jì)算源程序附文末，程序處理小孔圖像衍射的大致步驟是:1）讀取任意小孔的圖片；2）對(duì)衍射平面和光屏建立合理的坐標(biāo)系；3）計(jì)算衍射平面和光屏上任意兩點(diǎn)的距離；4）計(jì)算其復(fù)振幅和光強(qiáng)；5）畫(huà)出光強(qiáng)圖像和衍射圖.

2 衍射模擬

2.1 衍射模擬與理論計(jì)算對(duì)比驗(yàn)證

以圓孔為例，按照半波帶法，當(dāng)入射光波長(zhǎng) λ和圓孔半徑 R一定時(shí)，焦距 f變化，屏中心點(diǎn)O（0，0）對(duì)應(yīng)的半波帶數(shù)k會(huì)變化，并且該點(diǎn)的明暗隨著k值的變化而變化.k為偶數(shù)時(shí)，為暗條紋；k為奇數(shù)時(shí)，為亮條紋.k＝1是菲涅耳衍射和夫瑯禾費(fèi)衍射的分界點(diǎn)，當(dāng)k＞1時(shí)是菲涅耳衍射，k＜1是夫瑯禾費(fèi)衍射.故當(dāng)半波帶數(shù) k≤1時(shí)，衍射圖樣中心均為亮斑［11］.并且從中心向外是明暗相間的圓環(huán)，且 k越大（即z越?。靼淡h(huán)越多越密.k值由下式?jīng)Q定:

當(dāng) f＝0.2 m，λ＝500 nm，衍射圓孔半徑 R＝1 mm，光屏寬度h＝0.4 mm時(shí)，由圖2可見(jiàn)，整個(gè)光強(qiáng)分布以圓心為對(duì)稱(chēng)點(diǎn)，并且圓心點(diǎn)光強(qiáng)是減弱的（即產(chǎn)生暗斑）.將數(shù)據(jù)帶入式（4）計(jì)算得出 k＝10，即中心光強(qiáng)減弱，與程序模擬相符.

圖2 f＝0.2 m，λ＝500 nm，R＝1 mm圓孔衍射光強(qiáng)圖樣

當(dāng) f＝2 m，λ＝500 nm，衍射圓孔半徑 R＝1 mm，光屏寬度h＝4 mm時(shí)，由圖3可見(jiàn)，整個(gè)光強(qiáng)分布以圓心為對(duì)稱(chēng)點(diǎn)，并且圓心點(diǎn)光強(qiáng)是增強(qiáng)的（即產(chǎn)生亮斑）.將數(shù)據(jù)帶入式（4）計(jì)算得出 k＝1，即中心光強(qiáng)增強(qiáng)，與程序模擬相符.

當(dāng) f＝10 m，λ＝500 nm，衍射圓孔半徑 R＝1 mm，光屏寬度h＝10 mm時(shí)，由圖4可見(jiàn)，整個(gè)光強(qiáng)分布以圓心為對(duì)稱(chēng)點(diǎn)，并且圓心點(diǎn)光強(qiáng)是增強(qiáng)的（即產(chǎn)生亮斑）.將數(shù)據(jù)帶入式（4）計(jì)算得出 k＝0.2＜1，即中心光強(qiáng)增強(qiáng)，與程序模擬相符.

綜上所述，理論計(jì)算結(jié)果與程序模擬結(jié)果基本一致，證明模擬程序正確.由計(jì)算可知圖2是近場(chǎng)衍射圖像，圖3是近場(chǎng)衍射和遠(yuǎn)場(chǎng)衍射的分界點(diǎn)圖像，圖4為遠(yuǎn)場(chǎng)衍射圖像.從上述圖中可以得出，遠(yuǎn)場(chǎng)衍射圖像的中央亮斑比近場(chǎng)衍射的亮斑大，并且隨著k值逐漸減小，暗紋的光強(qiáng)趨于零，且在菲涅耳衍射時(shí)，中心的光強(qiáng)可能減弱.

圖3 f＝2 m，λ＝500 nm，R＝1 mm圓孔衍射光強(qiáng)圖樣

圖 4 f＝10 m，λ＝500 nm，R＝1 mm圓孔衍射光強(qiáng)圖樣

2.2 任意孔形及任意亮度的衍射模擬結(jié)果

運(yùn)用程序數(shù)值模擬正三角形和正六邊形的夫瑯禾費(fèi)衍射，驗(yàn)證正多邊形孔衍射特性，并對(duì)亮度不同的不規(guī)則小孔衍射特性進(jìn)行分析.

1）f＝2 m，λ＝500 nm，三角孔邊長(zhǎng) a＝1 mm，光屏寬度h＝6 mm的衍射如圖5所示.在垂直于三角形的3條邊的方向和3個(gè)角的方向上衍射現(xiàn)象均明顯，但邊方向衍射更為顯著，距離中心越遠(yuǎn)，衍射現(xiàn)象越弱.三角形的中央光強(qiáng)最強(qiáng)（即圖5的中央亮斑）.

2）f＝2 m，λ＝500 nm，六邊形孔邊長(zhǎng) a＝1 mm，光屏寬度h＝6 mm時(shí)的衍射如圖6所示.在垂直于六邊形各邊的方向衍射現(xiàn)象明顯，產(chǎn)生明暗相間的條紋，距離中心越遠(yuǎn)，衍射現(xiàn)象越弱，中央的衍射光強(qiáng)最強(qiáng).

圖 5 f＝2 m，λ＝500 nm，a＝1 mm三角形孔衍射圖樣

圖6 f＝2 m，λ＝500 nm，a＝1 mm六邊形孔衍射圖樣

3）f＝10 m，λ＝500 nm，圓孔半徑 R＝1 mm，光屏寬度h＝10 mm時(shí)的衍射如圖7所示.圓心光強(qiáng)最強(qiáng)，產(chǎn)生明暗相間的圓環(huán)，距離中心越遠(yuǎn)，衍射現(xiàn)象越弱.在夫瑯禾費(fèi)衍射中，以第一暗環(huán)為范圍的中央亮斑的光強(qiáng)占整個(gè)入射光束光強(qiáng)的84%，這個(gè)中央光斑稱(chēng)為艾里斑，由公式可得艾里斑的線(xiàn)半徑為

由圖7（a）可以看出艾里斑的線(xiàn)半徑大致為3.6 mm，與理論計(jì)算的結(jié)果大致相符.

圖7 f＝10 m，λ＝500 nm，R＝1 mm圓形孔衍射圖樣

4）f＝2 m，λ＝500 nm，面積為1 mm2的衍射平面上不規(guī)則小孔，光屏寬度h＝10 mm時(shí)的衍射現(xiàn)象如圖8所示.在垂直于各邊的方向上衍射現(xiàn)象明顯，產(chǎn)生明暗相間的條紋，且距離中心越遠(yuǎn)，衍射現(xiàn)象越弱.

圖8 f＝2 m，λ＝500 nm，面積為 1 mm2異形孔衍射圖樣

5）任意形狀任意亮度的光源也能通過(guò)該程序得到衍射圖像.如圖 9所示，f＝2 m，λ＝500 nm，正方形邊長(zhǎng)a＝1 mm（孔內(nèi)不同地方亮度值不同），光屏寬度h＝10 mm時(shí)的衍射圖像中可以看出，最大光強(qiáng)在正方孔中央，亮度值不同衍射強(qiáng)弱不同.

圖9 f＝2 m，λ＝500 nm，邊長(zhǎng)為1 mm不同亮度正方形孔衍射圖樣

3 結(jié)論

本文基于惠更斯-菲涅耳原理，用 Matlab進(jìn)行模擬仿真任意形狀任意亮度孔的衍射圖像，得到以下結(jié)論:1）同一衍射孔，在改變衍射面與光屏的距離時(shí)得到不同的衍射圖像.當(dāng)中心半波帶k為偶數(shù)時(shí)，中心光強(qiáng)減弱；k為奇數(shù)時(shí)，中心光強(qiáng)加強(qiáng)；k小于1時(shí)，一直處于加強(qiáng)狀態(tài).隨著焦距慢慢增大，幾何中心光強(qiáng)處于減弱狀態(tài)的情況就越來(lái)越少，最后一直處于加強(qiáng)狀態(tài)［13］.整個(gè)過(guò)程可以看做菲涅耳衍射到夫瑯禾費(fèi)衍射的一個(gè)過(guò)渡［13］.夫瑯禾費(fèi)衍射圖像隨著k值的減小，中央亮斑越來(lái)越大，其暗紋的光強(qiáng)也逐漸趨向于零.2）正多邊形衍射圖像成中心對(duì)稱(chēng)分布，每條邊的垂直方向衍射現(xiàn)象明顯，邊數(shù)越多，光強(qiáng)增強(qiáng)的方向越密集，并且在各角的方位上衍射呈條帶分布（正多邊形的邊數(shù)n為奇數(shù)時(shí)呈亮條帶，n為偶數(shù)時(shí)呈暗條帶）［5］.可以推斷出當(dāng)邊數(shù)趨向于無(wú)窮時(shí)，衍射圖像為圓孔衍射的圖像.且可以看出當(dāng)正多邊形的邊數(shù)為奇數(shù)時(shí)，在垂直邊方向上的衍射比角方向上的衍射更顯著.3）在遠(yuǎn)場(chǎng)衍射時(shí)，任意形狀任意亮度的小孔，中央光強(qiáng)最強(qiáng)，周?chē)难苌涔鈴?qiáng)與小孔各個(gè)位置的亮度有關(guān)，且遵循隨距離增大光強(qiáng)減弱的規(guī)律.