基于MATLAB的諧振腔穩(wěn)定性分析和高斯光束傳輸特性計(jì)算

1、 中北大學(xué)課程設(shè)計(jì)說(shuō)明書(shū) 2014/2015 學(xué)年第 一 學(xué)期學(xué) 院： 信息與通信工程學(xué)院 專 業(yè)： 光電信息工程 學(xué) 生 姓 名： 趙策 學(xué) 號(hào)： 1105104138 課程設(shè)計(jì)題目： 基于MATLAB的諧振腔穩(wěn)定性分析和 高斯光束傳輸特性計(jì)算 起 迄 日 期： 2015年1月 12日2015年1月30日 指 導(dǎo) 教 師： 王小燕 中北大學(xué)課程設(shè)計(jì)任務(wù)書(shū) 2014/2015 學(xué)年第 一 學(xué)期學(xué) 院： 信息與通信工程學(xué)院 專 業(yè)： 光電信息工程 學(xué) 生 姓 名： 趙策 學(xué) 號(hào)： 1105104138 課程設(shè)計(jì)題目： 基于MATLAB的諧振腔穩(wěn)定性分析和 高斯光束傳輸特性計(jì)算 起 迄 日 期：

2、2015年1月 12日2015年1月30日 指 導(dǎo) 教 師： 王小燕 課 程 設(shè) 計(jì) 任 務(wù) 書(shū)1 設(shè)計(jì)目的：在學(xué)習(xí)專業(yè)基礎(chǔ)課和專業(yè)課的基礎(chǔ)上，主要對(duì)激光原理，激光技術(shù)課程中出現(xiàn)的諸多理論模型進(jìn)行數(shù)值求解，學(xué)會(huì)將MATLAB用于光學(xué)仿真中，鍛煉運(yùn)用數(shù)值分析方法解決專業(yè)問(wèn)題的能力。本設(shè)計(jì)的主要目的如下：1.理解光束在自由空間傳輸?shù)腁BCD傳輸規(guī)律2.在MATLAB中運(yùn)用ABCD傳輸規(guī)律實(shí)現(xiàn)光線在諧振腔內(nèi)傳輸?shù)能壽E，考察諧振腔的穩(wěn)定性。3.理解高斯光束的q參數(shù)傳輸規(guī)律2設(shè)計(jì)內(nèi)容和要求（包括原始數(shù)據(jù)、技術(shù)參數(shù)、條件、設(shè)計(jì)要求等）：設(shè)計(jì)內(nèi)容和技術(shù)參數(shù)：1.查閱資料，理解諧振腔穩(wěn)定振蕩的條件2.如下圖所

3、示的諧振腔，用Matlab程序計(jì)算光線在腔內(nèi)的軌跡，演示腔的穩(wěn)定時(shí)光線在腔內(nèi)往返次數(shù)增加時(shí)光線軌跡。初始光線任意選擇。R1=2m，R2=1m，L=0.8m3.在如圖所示的平凹諧振腔內(nèi)，插入透鏡，分析透鏡放置什么位置時(shí)腔是穩(wěn)定腔。在穩(wěn)定腔的情況下，演示在腔內(nèi)往返100次以上的光線軌跡。中L R1=1000mmF=50mm4,上圖中，計(jì)算自在現(xiàn)高斯光束的q參數(shù)，并演示往返一周腔內(nèi)光斑半徑沿軸線的自在現(xiàn)曲線，取=0.5um3設(shè)計(jì)工作任務(wù)及工作量的要求包括課程設(shè)計(jì)計(jì)算說(shuō)明書(shū)(論文)、圖紙、實(shí)物樣品等：1） 課程設(shè)計(jì)說(shuō)明書(shū)一份；2） 電子文檔（說(shuō)明書(shū)、設(shè)計(jì)程序）一份4主要參考文獻(xiàn)：l 要求按國(guó)標(biāo)GB 7

4、71487文后參考文獻(xiàn)著錄規(guī)則書(shū)寫(xiě)，例：1 傅承義，陳運(yùn)泰，祁貴中.地球物理學(xué)基礎(chǔ).北京：科學(xué)出版社，1985至少五篇5設(shè)計(jì)成果形式及要求：1） 課程說(shuō)明書(shū)打印，并裝訂；2） 必要的光路圖、設(shè)計(jì)程序及仿真結(jié)果圖。6工作計(jì)劃及進(jìn)度：2015年 1月12日1月16日：了解設(shè)計(jì)題目及查找資料；1月17日1月21日：理解原理，確定各題目要求計(jì)算相關(guān)參數(shù)；2015年 1月22日1月23日: 結(jié)合各題目確定具體設(shè)計(jì)方案；2015年 1月24日2015年01月28日: 結(jié)合要求具體設(shè)計(jì)并編程仿真；2015年 01月28日01月30日: 整理課程設(shè)計(jì)說(shuō)明書(shū)；01月30日: 答辯或成績(jī)考核。 學(xué)科部審查意見(jiàn)：

5、簽字： 年 月 日摘要此次課程設(shè)計(jì)主要針對(duì)激光原理、激光技術(shù)課程中出現(xiàn)的諸多理論模型進(jìn)行數(shù)值求解，通過(guò)MATLAB軟件進(jìn)行仿真驗(yàn)證，從而鍛煉運(yùn)用數(shù)值分析方法解決專業(yè)問(wèn)題的能力,進(jìn)一步學(xué)習(xí)高斯光束的特性和傳播規(guī)律。關(guān)鍵詞：激光原理 高斯光束 MATLAB仿真設(shè)計(jì)AbstractThis course design principle of laser, the laser technology course many theoretical models appeared in numerical simulation by MATLAB software and ability to exer

6、cise using numerical methods to solve professional problems, further study characteristics and propagation of Gaussian beams. Key words: Principle of laser Gaussian beam The MATLAB simulation目 錄一、總體設(shè)計(jì)方案：7二、共軸球面腔穩(wěn)定性分析：7三、任務(wù)一具體設(shè)計(jì)：8四、任務(wù)二具體設(shè)計(jì)：9五、任務(wù)三具體設(shè)計(jì)：11六、心得體會(huì)：12一、總體設(shè)計(jì)方案：此次課程設(shè)計(jì)主要為了完成諧振腔穩(wěn)定性的分析、高斯光束的傳播規(guī)

7、律，理解諧振腔穩(wěn)定振蕩的條件，利用穩(wěn)定震蕩的條件區(qū)別穩(wěn)定腔、臨界腔、非穩(wěn)腔，以及利用Q參數(shù)來(lái)描述高斯光束的基本特征和傳輸規(guī)律。本設(shè)計(jì)借助MATLAB軟件進(jìn)行仿真設(shè)計(jì)和驗(yàn)證，快捷方便，可視性性強(qiáng)。根據(jù)三個(gè)任務(wù)不同的參數(shù)在MATLAB中進(jìn)行編寫(xiě)程序，通過(guò)仿真結(jié)果判斷是否達(dá)到預(yù)定的要求，通過(guò)不斷嘗試，找到最理想?yún)?shù)完成實(shí)驗(yàn)任務(wù)。二、共軸球面腔穩(wěn)定性分析：利用幾何光學(xué)的光線矩陣分析法，分析共軸球面腔的穩(wěn)定性條件非常方便，用幾何光學(xué)的方法分析諧振腔的實(shí)質(zhì)是研究光線在腔內(nèi)往復(fù)反射的過(guò)程。如圖1所示的共軸球面腔，由曲率半徑為R1和R2兩個(gè)球面鏡M1和M2構(gòu)成，腔長(zhǎng)為L(zhǎng)。圖1 光線在共軸球面腔中的往返傳播我們

8、來(lái)分析光線的傳播過(guò)程：光線從M1面出發(fā)，向M2面方向行進(jìn)，我們用矩陣（1.1）來(lái)描述當(dāng)光線在自由空間中行進(jìn)距離為L(zhǎng)時(shí)所引起的坐標(biāo)變換。1L01 （1.1）當(dāng)光線傳播到M2鏡面時(shí)，在球面鏡上發(fā)生反射時(shí)，根據(jù)球面鏡對(duì)旁軸光線的反射規(guī)律，我們用矩陣（1.2）表示，稱為球面鏡反射矩陣。其中R為球面鏡的曲率半徑，而F=R/2為球面鏡對(duì)旁軸光線的焦距。T=10-2R1=10-1F1 （1.2）當(dāng)光線完成腔內(nèi)一次往返，總的變換為公式（1.3）。T=ABCD=10-2R11L0110-2R11L01 （1.3）我們關(guān)注的是在什么情況下旁軸光線能在腔內(nèi)往返任意多次而不致橫向逸出腔外，我們根據(jù)傳播變換公式可得出共

9、軸球面腔的穩(wěn)定性條件，如公式（1.4）。-1<12A+D<1 （1.4）滿足條件，就保證了旁軸光線能在腔中往返多次而不致于從側(cè)面逸出（只要鏡面橫向尺寸足夠大）。反之，如果不滿足上述條件，旁軸光線在腔內(nèi)經(jīng)歷有限次往返后必將橫向逸出腔外。三、任務(wù)一具體設(shè)計(jì)：1、用MATLAB仿真給定光學(xué)諧振腔的穩(wěn)定性和往返光學(xué)軌跡。具體編程如下：L=input('input L:='); r0=3;theta0=0.01; %初始化光線方程相關(guān)參數(shù) R1=200;R2=100; C=r0;theta0; T1=1,L;0,1; T2=1,0;-2/R2,1; T3=T1; T4=1,0

10、;-2/R1,1; T=T4*T3*T2*T1; S=(T(1,1)+T(2,2)/2 %穩(wěn)定性判據(jù) for k=1:100 %利用for語(yǔ)句完成光線往返100次 x=0:.1:L; y=C(1)+C(2).*x; %.*是矩陣相乘的意思（相同位置的數(shù)相乘） plot(x,y,'k'),hold on %光線從M1到M2的軌跡 C=T1*C; C=T2*C; x=L:-.1:0; y=C(1)-C(2).*(x-L); plot(x,y,'r'),hold on %光線從M1到M2的軌跡 C=T3*C; C=T4*C; 2、仿真結(jié)果：當(dāng)L=80時(shí)，S=-0.76

11、,所以為穩(wěn)定腔。仿真圖如圖2所示。當(dāng)L=100時(shí)，S=1為臨界腔。仿真圖如圖3所示。圖2 穩(wěn)定腔往返光線軌跡（S=-0.76） 圖3 臨界腔往返光學(xué)軌跡(S=-1)四、任務(wù)二具體設(shè)計(jì)：1、在平凹透鏡中間插入透鏡，利用MATLAB確定透鏡的位置使得整個(gè)系統(tǒng)為穩(wěn)定光學(xué)諧振腔。根據(jù)系統(tǒng)參數(shù)，編寫(xiě)代碼如下：L1=input('input L1:='); r0=3;theta0=0.01; %初始化光線方程的相關(guān)參數(shù) R=1000;F=50;L=800; C=r0;theta0; T1=1,L1;0,1; T2=1,0;-1/F,1; T3=1,L;0,1;T4=1,0;-2/R,1;

12、T=T1*T2*T3*T4*T3*T2*T1; S=(T(1,1)+T(2,2)/2 %穩(wěn)定性判據(jù) for k=1:100 x=0:.1:L1; y=C(1)+C(2).*x; plot(x,y,'r'),hold on %從平面鏡到薄透鏡的光線軌跡 C=T1*C; C=T2*C; x=L1:.1:(L1+L); y=C(1)+C(2).*(x-L1); plot(x,y,'r'),hold on %從薄透鏡到球面鏡的光線軌跡 C=T3*C; C=T4*C; x=(L1+L):-.1:L1; y=C(1)-C(2)*(x-L1-L); plot(x,y,'

13、;r'),hold on %從球面鏡返回薄透鏡的光線軌跡 C=T3*C; C=T2*C; x=L1:-.1:0; y=C(1)-C(2)*(x-L1); plot(x,y,'r'),hold on %從薄透鏡返回平面鏡的光線軌跡 C=T1*C; 2、仿真結(jié)果：當(dāng)L1=49時(shí)，S=0.35，為穩(wěn)定腔，光線往返150次軌跡圖如圖4所示。當(dāng)L1=40或52.9時(shí)，S=±1，此時(shí)該諧振腔為臨界腔，光線往返150次軌跡圖如圖5所示。故，L的范圍應(yīng)該是40mm<L1<52.9mm。圖4 穩(wěn)定模式（S=0.35）圖5 臨界模式（S=-1）五、任務(wù)三具體設(shè)計(jì)：1、

14、高斯光束q參數(shù)：是將高斯光束基本特征的兩個(gè)參數(shù)w(z)和R(z)統(tǒng)一在一個(gè)表達(dá)式中，一旦知道某一位置高斯光束q參數(shù)就能知道該位置處的w(z)和R(z)。計(jì)算公式見(jiàn)公式（1.5）。1q(z)=1R(z)-iw2(z) （1.5）2、高斯光束的自再現(xiàn)變換：如果高斯光束通過(guò)透鏡后結(jié)構(gòu)不發(fā)生變化，即參數(shù)w0或f不變，則稱這樣的變換為自再現(xiàn)變換。根據(jù)自再現(xiàn)模定義，穩(wěn)定腔的任意高斯模在腔內(nèi)往返一周后應(yīng)能重現(xiàn)其自身。自在現(xiàn)q參數(shù)的計(jì)算方法參考公式（1.6）。1q=D-A2C±i1-A+D2/4C （1.6）3、要描繪高斯光束在橫截面上的場(chǎng)振幅分布及光斑半徑的變化情況，可根據(jù)公式（1.7）來(lái)求得。W0為腰斑半徑。wz=w01+(zf)2 （1.7）可見(jiàn)，光斑半徑隨坐標(biāo)Z按雙曲線規(guī)律擴(kuò)展。其中w0可用公式（1.8）來(lái)求解。其中f為共焦參數(shù)。w0=f （1.8）4、仿真結(jié)果如圖6所示。圖6 往返一周腔內(nèi)光斑半徑沿軸線的自在現(xiàn)曲線六、心得體會(huì)：此次課程設(shè)計(jì)內(nèi)容相對(duì)簡(jiǎn)單，工程量小，完成效率較高，質(zhì)量較好。通過(guò)