基于Matlab的激光脈寬測量實驗

張　剛，焦志勇,黃小燦，劉　鵬

?

基于Matlab的激光脈寬測量實驗

張剛，焦志勇,黃小燦，劉鵬

(中國石油大學(華東) 理學院 物理與光電工程系，山東 青島 266580)

采用樣條插值法，利用Matlab軟件實現(xiàn)了調Q鎖模脈沖包絡寬度自動測量. 使用該程序對GaAs被動調Q鎖模激光器輸出的脈沖數據進行處理，可以直接獲得脈沖包絡寬度，簡化了峰值功率計算的步驟.

調Q鎖模；調Q包絡寬度；Matlab；樣條插值

自1960年梅曼發(fā)明第一臺紅寶石激光器以來，固體激光技術與應用迅速發(fā)展[1-5]. 為了激發(fā)學生對激光應用領域的興趣，大學物理教材中都會對激光知識做詳細介紹，物理實驗中也有針對性地添加激光輸出特性觀測等方面的內容[6-8]. 學生通過實驗獲得激光脈沖模型，測量脈沖寬度，進而獲得脈沖峰值功率，更加直觀地認識脈沖激光. 其中，調Q鎖模技術以其結構簡單、峰值功率高等特點獲得關注. 調Q鎖模脈沖是鎖模激光脈沖的幅度被調Q包絡周期性調制[9-10]，此種狀態(tài)下，激光工作物質能積累更多的反轉粒子數，從而獲得比連續(xù)鎖模激光器更高的峰值功率. 峰值功率的計算依賴于脈沖寬度以及輸出功率的精確測量. 但是，調Q鎖模脈沖具有調Q包絡和鎖模脈沖2個重復頻率，這種不穩(wěn)定性增加了學生實驗過程中精確測量脈沖寬度的難度. 本論文將利用樣條差值描點法，通過Matlab程序實現(xiàn)調Q鎖模激光器中調Q包絡寬度的測量.

1　實驗裝置和輸出脈沖

圖1為實驗中使用的調Q鎖模激光器的結構圖. 尺寸為3 mm×3 mm×8 mm的摻釹釩酸镥(Nd∶LuVO4)晶體充當激光工作物質，其中釹離子(Nd3+)摻雜原子百分數為0.5 %，該晶體左右兩端鍍有808 nm 和1 342 nm的高透膜，晶體側面用錫箔包裹，置于銅塊中心位置，通過水冷方式將溫度控制在20 ℃，以抑制熱透鏡效應.

圖1　GaAs被動調Q鎖模激光器裝置圖

激光腔采用“Z”形腔結構. 輸入鏡M1的曲率半徑為200 mm，泵浦端鍍有808 nm增透膜，另一端鍍有1 342nm高反膜. 折疊鏡M2和M3分別是曲率半徑為500 mm和200 mm的凹面鏡，鍍有1 342 nm高反膜層，充當諧振腔反射鏡. 輸出鏡M4是平鏡，在1 342 nm波段的透過率為8%. M1和M2,M2和M3,M3和M4之間的距離分別為250 mm，650 mm和100 mm. 實驗中使用的被動調制元件是小信號透過率為93.9%的GaAs飽和吸收體，靠近輸出鏡放置.

實驗中選用帶寬為500 MHz的數字示波器(Tektronix Inc., USA)測量和記錄波形，用MAX 500AD激光功率計(Coherent, USA)測量平均輸出功率.

2　實驗數據

實驗中測得的泵浦功率和平均輸出功率的關系如圖2所示. 由實驗結果可知，當泵浦功率為0.9 W時，開始有激光脈沖輸出，當泵浦功率為8 W時，獲得最高的輸出功率為452 mW. 調Q包絡的重復率隨泵浦功率的變化如圖3所示. 由圖3可知，脈沖重復頻率隨泵浦功率增加而增加. 泵浦功率為8 W時，重復頻率為140 kHz. 圖4是在泵浦功率為8 W時，單個調Q包絡的波形圖，鎖模脈沖的調制深度為80%. 通過示波器記錄的數據，逐點讀數，可以獲得此泵浦功率下調Q包絡寬度為75 ns.

圖2　平均輸出功率隨泵浦功率變化

圖3　調Q包絡的重復頻率隨泵浦功率變化

圖4　單個鎖模脈沖波形圖

3　脈寬測量程序

前述脈沖寬度讀取過程中，學生需要逐點讀取數據，工作量較大，數據處理時間長. 為了更加簡潔和精確測量包絡寬度，需要程序自動運算示波器采集的數據，最終得到需要的包絡寬度. Matlab以其功能強大、界面友好的特點，將使用者從繁瑣的底層編程中解放出來，使他們有更多的時間和精力去探究科學問題本身[11-12]. 因此引導學生使用Matlab實現(xiàn)調Q包絡寬度的測量.

將程序劃分為4個功能模塊：數據選擇模塊、數據揀峰模塊、數據模擬模塊和脈寬輸出模塊.

3.1數據選擇模塊

在利用Matlab程序計算脈沖寬度時，希望獲得具有通用性的程序，而不是僅僅針對某個特定數據文件. 不論示波器何時獲取的脈沖數據，只需要在數據提取時選擇相應文件即可，而不必更改源代碼. 要完成這樣的要求，數據選擇模塊必不可少. 因此在該部分，使用“uigetfile”函數來完成文件選擇任務.

相應程序源代碼如下：

1- clc

2- clear all

3-

4- [f,p]=uigetfile(′*.csv′,′選擇文件′);

5- if f

6- y=xlsread(strcat(p,f));

7- end

8- t=1∶1∶500;

9- z=y(∶,2);

10- ys=z;

3.2數據揀峰模塊

類似于圖4所示調Q鎖模脈沖，調Q包絡內含有許多鎖模脈沖，想要獲取脈沖寬度，需要知道調Q包絡的寬度. 需要對原始數據處理，找出每個鎖模脈沖的峰值，然后選擇合適差值方式處理數據，為繪制調Q包絡打下基礎.

具體說來，第一部分首先用嵌套的“diff”函數對脈沖數據二階求導，再根據導數的正負值用于尋找峰谷點. 第二部分對數據進行插值處理. 為了獲得更多的數據點，也為了獲得更好的包絡線，使用3次樣條插值[13-15]. 第三部分用循環(huán)嵌套約束插值結果帶來的誤差. 這主要是因為在峰值部分有1個拐點，插值過程中就可能插入與原始數據最大值偏差較大的數據點，這樣就可能導致最終繪制出的包絡線有較大的誤差. 為此需要對這兩點之間的插值進行約束限制. 代碼如下：

11- %% 循環(huán)判斷

12- for i=1

13- while 1

14- max_location=find(diff(sign(diff(z)))==-2)+1;

15- z_max=z(max_location);

16- T_max=t(max_location);

17- min_location=find(diff(sign(diff(z)))==2)+1;

18- z_min=z(min_location);

19- T_min=t(min_location);

20- zmax=spline(T_max,z_max,t);

21- zmin=spline(T_min,z_min,t);

22- figure(2)

23- baoluo2=(zmax+zmin)/2;

24- z=baoluo2;

25- if max(zmax)>0.9*max(ys)

26- if max(zmax)<1.2*max(ys)

27- break;

28- end

29- end

30- end

31- end

3.3數據模擬模塊

在Matlab中，繪制二維圖像最基本的是plot函數. 因此在本程序中，使用該函數完成包絡線的繪制.

32- hold on

33- plot(zmax,′k′)

34- % plot(zmin,′g′)

35- % plot(baoluo2,′r′)

36- plot(ys,′*′)

37- maxm=max(zmax);

38- minn=min(zmax);

39- mid=(maxm-minn)/2+minn;

40- eps=0.007;

41- j=1;

42- for i=1∶length(zmax)

43- if abs(zmax(i)-mid)

44- midy(j)=zmax(i);

45-j=j+1;

46- end

47- end

48- y1=midy′;

49- m=length(midy);

50- for j=1∶m

51-ind(j)=find(midy(j)==zmax);

52- end

53- x1=ind′;

54- m=length(x1);

55- for i=1∶m-1

56-if abs(x1(i)-x1(i+1))>10

57- x(1)=x1(i);

58-x(2)=x1(i+1);

59-y3(1)=y1(i);

60-y3(2)=y1(i+1);

61- end

62- end

63-

其中，默認狀態(tài)下Matlab軟件只能顯示1條命令繪制的圖像，后續(xù)繪制的圖像會替換掉前圖. 為避免該現(xiàn)象，本段程序使用了“hold on”命令，該命令可以完成圖形保持動作，命令執(zhí)行過程中，新舊兩圖同時顯示. 這樣，最終完成圖像繪制.

3.4脈寬輸出模塊

前3個模塊運行結束后，包絡繪制完成. 數據中點的位置坐標也得到. 最后模塊的功能是求出脈寬T，并將結果輸出. 該部分源代碼如下：

64- disp('中點坐標[x,y3]為：′)

65- disp(unique([x,y3],′rows′))

66- disp(′該脈沖的脈寬為：′)

67- T=(x(2)-x(1))

68- plot(x,y3,′r*′)

69- hold off

4　程序結果處理

圖5為程序運行結果，可得脈沖包絡的寬度為76 nm，與前述結果(75 nm)相比，偏差僅為1.3%，符合實驗對數據精度的要求. 根據前面所得數據，通過公式Pp=Pa/(FT)，可計算出泵浦功率為8 W時脈沖峰值功率為42.5 W. 其中，Pa是平均輸出功率，F(xiàn)是調Q包絡的重復率，T是調Q包絡的寬度.

圖5　插值描點法程序運行結果

5　結束語

調Q鎖模脈沖包絡寬度讀取過程繁瑣，易錯. Matlab強大的數值處理功能,為調Q鎖模激光器的脈寬測量提供了便利. 通過編寫Matlab程序, 可以一鍵獲取調Q包絡的寬度. 實驗中，GaAs被動調Q鎖模脈沖寬度為76 ns，峰值功率為42.5 W. 實驗改進并使用Matlab程序后，計算過程更簡便，提高了學生實驗積極性.

[1]Maiman T H. Stimulated optical radiation in ruby [J]. Nature, 1960,187(4736):493-494.

[2]Newman R. Excitation of the Nd3+fluorescence in CaWO4by recombination radiation in GaAs [J]. J. Appl. Phys., 1963,34(2):437-438.

[3]Fields R A, Brinbaum M, Fincher C L. Highly efficient Nd:YVO4diode-laser end-pumped laser [J]. Appl. Phys. Lett., 1987,51(23):1885-1886.

[4]Barnes N P．Solid-state lasers from all efficiency persptive [J]. IEEE Journal of Selected Topics in Quantum Electronics, 2007,13(3):435-446.

[5]Han Shuo, Li Xianlei, Xu Honghao, et al. Graphene Q-switched 0.9μm Nd∶La0.11Y0.89VO4laser [J]. Chinese Optics Letters, 2014,12(1):33-35.

[6]吳弘,張禮,陳杰. 基于FPGA的邁克耳孫干涉儀測量激光波長[J]. 物理實驗,2015,35(10):32-36.

[7]楊燕婷,王敏,周譽昌. 用超聲光柵測量氦氖激光的波長[J]. 大學物理,2009,28(7):43-44.

[8]李莉,齊曉慧,劉秉琦,等. LD泵浦Nd∶YAG連續(xù)激光器轉換性能實驗研究[J]. 物理實驗,2009,29(11):20-23.

[9]Datta P K, Mukhopadhyay S, Das S K. Enhancement of stability and efficiency of a nonlinear mirror mode-locked Nd∶YVO4oscillator by an active Q-switch [J]. Opt. Express, 2004,12(17):4041-4046.

[10]Theobald C, Weitz M, Knappe R, et al. Stable Q-switch mode-locking of Nd∶YVO4lasers with a semiconductor saturable absorber [J]. Appl. Phys. B, 2008,92(1):1-3.

[11]莊建，青莉. 基于Matlab的光學實驗仿真平臺 [J]. 成都大學學報(自然科學版),2006,25(1):70-72.

[12]陳守滿. MATLAB在光學諧振腔設計中的應用[J]. 安康師專學報, 2004,16(6):71-73.

[13]Mckinley S, Levine M. Cubic spline interpolation [J]. Numer. Math. J. Chinese Univ., 1999,64(1):44-56.

[14]Cuche E, Marquet P, Depeursinge C. Aperture apodization using cubic spline interpolation: application in digital holographic microscopy [J]. Optics Communications, 2000,182(1-3):59-69.

[15]胡俊華，田錦會. 基于樣條插值的第二代小波設計與應用[J]. 計算機與現(xiàn)代化,2013(7):76-79.

[責任編輯：任德香]

Measuring the envelope width of Q-switched laser pulse based on Matlab

ZHANG Gang, JIAO Zhi-yong, HUANG Xiao-can, LIU Peng

(Department of Physics and Optoelectronic Engineering, College of Science,China University of Petroleum (East China)， Qingdao 266580, China)

Based on spline interpolation and Matlab, the automatic measurement of Q-switched and mode-locked laser pulse width was realized. The pulse width could be obtained by processing the pulse data stored in oscilloscope directly. The method simplified the calculation steps of the peak power of the pulse.

Q-switched and mode-locked; Q-switched envelope width; Matlab; spline interpolation

2016-05-20

山東省自然科學基金(No.ZR2014FP008)；山東省研究生教育創(chuàng)新計劃(No.SDYY14141)

張剛(1984-)，男，山東泰安人，中國石油大學(華東)理學院物理與光電工程系講師,博士，主要研究方向為激光物理與技術.

TN248.1

A

1005-4642(2016)10-0001-04