任志浩 馬琨

摘 要： 根據(jù)等厚干涉的基本理論，利用Matlab平臺(tái)編程進(jìn)行仿真，模擬劈尖干涉實(shí)驗(yàn)在不同參數(shù)條件下的光強(qiáng)分布情況。并用Matlab自帶的Peaks函數(shù)模擬利用劈尖干涉對平面玻璃表面進(jìn)行缺陷檢測的干涉條紋。仿真所得到的干涉圖像細(xì)致逼真，并且能夠動(dòng)態(tài)地反映出干涉條紋隨參數(shù)的改變而發(fā)生變化的過程，使得物理規(guī)律形象直觀。Matlab模擬出的干涉圖像與實(shí)驗(yàn)結(jié)果非常吻合，還可以靈活地調(diào)節(jié)參數(shù)，大大擴(kuò)展了劈尖干涉問題的研究途徑與方法，有助于學(xué)生更加深刻地理解物理光學(xué)的現(xiàn)象和規(guī)律，也為光學(xué)的理論教學(xué)和實(shí)驗(yàn)方式提供了便利。

關(guān)鍵詞： Matlab； 劈尖干涉； 數(shù)值模擬； 表面檢測

中圖分類號(hào)： TN012?34； O436.1 文獻(xiàn)標(biāo)識(shí)碼： A 文章編號(hào)： 1004?373X（2016）06?0108?02

Simulation of wedge interference based on Matlab

REN Zhihao， MA Kun

（Faculty of Science， Kunming University of Science and Technology， Kunming 650500， China）

Abstract： Based on the fundamental theory of equal thickness interference， the light intensity distribution at different parameters in the wedge interference experiments is simulated by means of Matlab platform programming. The interference fringes conducting the defect detection of plain glass surface by wedge interference are simulated with Peaks function of Matlab. The interference image obtained from simulation is meticulous and realistic， and can dynamically reflect the changing process of interference fringes with the change of parameters， which makes the physics rules visual and intuitive. The interference image simulated by Matlab is identical with experimental result. The parameters can be set flexibly， which greatly extends the research approach and method of wedge interference subject， is helpful for students to profoundly understand the phenomenon and rules of physical optics， and provides the convenience for optical theory teaching and experimental mode selection.

Keywords： Matlab； wedge interference； numerical simulation； surface detection

0 引 言

波動(dòng)光學(xué)中的等厚干涉理論在大學(xué)物理課程的教學(xué)中占有很重要的地位。目前，國內(nèi)大部分的物理教材都把這部分內(nèi)容作為思考題或課后習(xí)題供給學(xué)生參考。但由于等厚干涉的理論比較復(fù)雜和抽象，一直以來都是教學(xué)過程中的重點(diǎn)和難點(diǎn)。此外，由于實(shí)驗(yàn)條件的限制，學(xué)生很難在課堂上通過實(shí)驗(yàn)來觀察干涉現(xiàn)象，也就無法獲得直觀的認(rèn)識(shí)。隨著計(jì)算機(jī)技術(shù)的飛速發(fā)展，計(jì)算機(jī)輔助教學(xué)（CAI）越來越受到人們的青睞，這為解決上述難題提供了途徑。

Math Works公司開發(fā)的Matlab（Matrix Laboratory）軟件是當(dāng)今國際上公認(rèn)的在數(shù)值計(jì)算和可視化圖形圖像處理領(lǐng)域最為優(yōu)秀的工程軟件之一，廣泛地應(yīng)用于數(shù)學(xué)分析、系統(tǒng)仿真、信息光學(xué)、數(shù)字圖像處理等領(lǐng)域[1]。Matlab強(qiáng)大的繪圖和動(dòng)畫功能可以生動(dòng)形象且準(zhǔn)確地反應(yīng)出連續(xù)變化的物理光學(xué)現(xiàn)象。因此，將Matlab與光學(xué)教學(xué)有機(jī)地結(jié)合起來[2?4]，有助于加深學(xué)生對基本概念和現(xiàn)象的理解以及認(rèn)識(shí)，激發(fā)學(xué)生的學(xué)習(xí)熱情，同時(shí)也能夠使教學(xué)內(nèi)容更加豐富，從而提高教學(xué)質(zhì)量。

1 劈尖干涉原理

如圖1所示，兩塊平面玻璃，一塊玻璃水平放置，另一塊玻璃與水平放置的玻璃形成一定的夾角。這時(shí)，在兩塊玻璃之間形成空氣劈尖。兩玻璃片的交線稱為棱邊，在平行于棱邊的線上，劈尖的厚度是相等的。

當(dāng)單色平面光垂直入射劈尖裝置時(shí)，AB面的反射光和AC面的反射光會(huì)發(fā)生干涉。設(shè)空氣的折射率為n0，則在劈尖上下表面反射的兩束光之間的光程差為[5]：

[δ=2n0d+λ2] （1）

圖1 空氣劈尖示意圖

由于從空氣劈尖的上表面和從空氣劈尖的下表面反射的情況不同，所以存在附加的光程差，即半波損失。因此，在厚度為d處出現(xiàn)亮暗紋的條件為[5]：

[δ=2n0d+λ2=kλ， k=1，2，…亮紋δ=2n0d+λ2=2k+1λ2， k=1，2，…暗紋] （2）

假設(shè)反射光的光強(qiáng)相同均為I0，則干涉光的強(qiáng)度可以表示為[6]：

[I=I0cos2π2n0dλ+12] （3）

2 數(shù)值模擬

設(shè)空氣的折射率n0=1，反射光的光強(qiáng)I0設(shè)為1。分別改變劈尖的夾角和入射光的波長，就能得到干涉條紋分布的變化情況。相對于以往的資料[7?10]中采用二維平面圖顯示結(jié)果，本文采用三維圖像，這樣更加直觀。在Matlab中編程，還可以將反映連續(xù)變化過程的關(guān)鍵幀組合起來以GIF的格式輸出。這樣就可以脫離Matlab的平臺(tái)，直接將圖片添加到PPT中以動(dòng)態(tài)圖片的形式放映，更加便捷。

2.1 改變劈尖夾角的大小

如圖2所示，入射光為波長λ=589.3 nm的鈉黃光，劈尖夾角依次增大。從圖2中可以發(fā)現(xiàn)，干涉條紋隨著夾角的增大變得越來越密。由式（2）可知這與厚度d有關(guān)，實(shí)際中常用這種方法測量金屬細(xì)絲的直徑等[11]。

圖2 不同劈尖夾角的干涉條紋圖

2.2 改變?nèi)肷涔獠ㄩL

如圖3所示，劈尖夾角固定不變，入射光波長依次增加。由圖3可知，干涉條紋隨著波長的增加變得越來越疏。由式（2）可知，隨著波長的增加，相同區(qū)域內(nèi)能顯示的條紋級(jí)數(shù)就越少。利用這種方法可以測出入射光的波長。此外，利用光的三原色疊加原理能夠細(xì)致的仿真出白光干涉的圖像。

圖3 不同入射光波長的干涉條紋圖

3 缺陷檢測的仿真

劈尖干涉的一個(gè)典型應(yīng)用就是檢測元件表面的缺陷。大學(xué)物理課中也有類似的習(xí)題，要求學(xué)生根據(jù)條紋的形狀來判斷缺陷的凹凸性。但是，有特定形狀缺陷的元件不易加工，因而很難用實(shí)驗(yàn)去觀察得出結(jié)論。然而用Matlab編程能夠比較容易地模擬出元件表面有凹陷或凸起的缺陷時(shí)產(chǎn)生的干涉圖。如圖4所示，劈尖干涉裝置虛線框中的元件表面存在缺陷，缺陷形狀用Matlab自帶的Peaks函數(shù)來模擬。這樣模擬得到的缺陷既有凹陷部分也有凸起部分，而且有不同程度的凹陷和凸起，有利于對比。由于缺陷的存在，仿真得到的干涉條紋如圖5所示。

圖4 元件表面缺陷示意圖

圖5 元件表面有缺陷時(shí)的干涉條紋圖

從圖5中可以清楚地觀察到缺陷凸起的部分所形成的干涉條紋向右邊彎曲了，條紋也相應(yīng)的變密了。缺陷部分凸起的越厲害條紋彎曲的程度也就越大，凹陷部分的情況也類似。沒有缺陷的部分干涉條紋則和棱邊平行，寬度完全相等。仿真得到的干涉圖像與實(shí)驗(yàn)結(jié)果是一致的，但是清晰度更高，更沒有實(shí)驗(yàn)條件的諸多限制。

4 結(jié) 語

文中僅以劈尖干涉為例進(jìn)行仿真，事實(shí)上利用Matlab能夠?qū)Υ蠖鄶?shù)的光學(xué)問題進(jìn)行數(shù)值模擬。這樣可以很容易地得到不同參數(shù)設(shè)置下的精確結(jié)果，而且調(diào)節(jié)起來十分靈活，克服了實(shí)驗(yàn)中的種種限制。數(shù)值仿真得到的圖像能將抽象的物理規(guī)律很直觀地展現(xiàn)出來，易于學(xué)生理解和認(rèn)識(shí)，這對于提高教學(xué)質(zhì)量起到了很好的輔助作用。要值得注意的是，數(shù)值仿真的前提條件一般都非常理想化，因而并不能夠完全代替實(shí)際的實(shí)驗(yàn)操作。實(shí)驗(yàn)操作是培養(yǎng)動(dòng)手能力的關(guān)鍵環(huán)節(jié)，只有將模擬仿真和實(shí)驗(yàn)操作結(jié)合起來才能相得益彰。

致謝：在此特別感謝理學(xué)院錢曉凡教授。錢教授的大學(xué)物理課幽默風(fēng)趣，深受學(xué)生的喜愛，筆者也受益匪淺。

參考文獻(xiàn)

[1] 張志涌.精通Matlab R2011a[M].北京：北京航空航天大學(xué)出版社，2011.

[2] 周群益，侯兆陽，劉讓蘇.Matlab可視化大學(xué)物理學(xué)[M].北京：清華大學(xué)出版社，2011.

[3] 錢曉凡.信息光學(xué)數(shù)字實(shí)驗(yàn)室[M].Matlab版.北京：科學(xué)出版社，2015.

[4] 歐攀.高等光學(xué)仿真：光波導(dǎo)·激光[M].Matlab版.北京：北京航空航天大學(xué)出版社，2011.

[5] 梁銓廷.物理光學(xué)[M].4版.北京：電子工業(yè)出版社，2012.

[6] BORN M， WOLF E. Principles of optics [M]. 7th ed. Cambridge： Cambridge University Press， 2005.

[7] 李繼軍，程慶華，楊長銘.基于Matlab GUI的薄膜干涉仿真[J].洛陽師范學(xué)院學(xué)報(bào)，2011，30（11）：16?19.

[8] 劉建龍，高揚(yáng).等厚干涉實(shí)驗(yàn)中缺陷檢測的數(shù)值模擬[J].物理實(shí)驗(yàn)，2015，35（1）：23?26.

[9] 張晶，李建映，鄧婷，等.基于Matlab的物理光學(xué)實(shí)驗(yàn)教學(xué)研究[J].廣西師范學(xué)院學(xué)報(bào)（自然科學(xué)版），2014，31（1）：122?127.

[10] 李小燕，馮卓宏，邱俊才.Matlab在大學(xué)物理教學(xué)中的應(yīng)用[J].實(shí)驗(yàn)技術(shù)與管理，2010，27（11）：124?126.

[11] 尹志勇，賈同福，楊洪志，等.用劈尖干涉測金屬線脹系數(shù)[J].大學(xué)物理實(shí)驗(yàn)，2010，23（2）：45?46.