【正文】 %倍頻晶體截面在 y軸上的長度 [x,y]=meshgrid(x,y)。 %熱功率密度 x=0::3。 %BBO晶體的徑向熱導(dǎo)率 a=15。 %基頻光的束腰半徑 p=5。) BBO倍頻晶體的平面溫度分布： c=。) zlabel(39。) ylabel(39。 xlabel(39。 BBO倍頻晶體的溫度分布 39。 else z(i,j)=q*w^2/(2*b)*log(sqrt(r(i,j))/)+a。 if sqrt(r(i,j))wamp。 for i=1:nx。 nx=length(x)。 %倍頻晶體截面在 y軸上的長度 [x,y]=meshgrid(x,y)。 %熱功率密度 x=0::3。%BBO冷卻溫度 b=。 %基頻光的束腰半徑 p=5。 BBO 晶體的熱效應(yīng) BBO倍頻晶體的三維溫度分布 ： c=。23度 39。)。 gtext(39。17度 39。)。) gtext(39。) ylabel(39。 xlabel(39。不同冷卻溫度下 x=a/2處 KTP倍頻晶體的溫度分布 39。c39。r39。g39。b39。 else z3(n)=q*w^2/(2*b)*log(sqrt(r(n))/)+23。 end end for n=1:length(r) if sqrt(r(n))wamp。sqrt(r(k))w z2(k)=(w^(k)).*q/(4*b)q*w^2/(2*b)*log(w/)+20。 else z1(j)=q*w^2/(2*b)*log(sqrt(r(j))/)+17。 end end for j=1:length(r) if sqrt(r(j))wamp。sqrt(r(i))w z(i)=(w^(i)).*q/(4*b)q*w^2/(2*b)*log(w/)+15。 r=(()./10).^2。 %KTP晶體的徑向熱導(dǎo)率 q=p*c/(pi*w^2)。 %基頻光的束腰半徑 p=5。) KTP 倍頻晶體冷卻溫度變化時的溫度 分布： c=。, 39。,39。,39。 legend(39。39。)。 gtext(39。39。)。) gtext(39。) ylabel(39。 xlabel(39。不同基頻光半徑下 x=a/2處 KTP倍頻晶體的溫度分布 39。k.39。g:39。r39。b39。 else z3(n)=q3*w3^2/(2*b)*log(sqrt(r(n))/)+a。 end end for n=1:length(r) if sqrt(r(n))w3amp。sqrt(r(k))w2 z2(k)=(w2^(k)).*q2/(4*b)q2*w2^2/(2*b)*log(w2/)+a。 else z1(j)=q1*w1^2/(2*b)*log(sqrt(r(j))/)+a。 end 西安工業(yè)大學(xué)畢業(yè)設(shè)計（論文） 31 end for j=1:length(r) if sqrt(r(j))w1amp。sqrt(r(i))w z(i)=(w^(i)).*q/(4*b)q*w^2/(2*b)*log(w/)+a。 r=(()./10).^2。 %基頻光強為 5w時的熱功率密度 q3=p*c/(pi*w3^2)。 %基頻光強為 5w時的熱功率密度 q1=p*c/(pi*w1^2)。 %基頻光強度 b=。 %基頻光的束腰半徑 w3=。 %基頻光的束腰半徑 w1=。 %KTP晶體對基頻光的吸收率 a=15。功率 20w39。功率 15w39。功率 10w39。功率 5w39。)。 gtext(39。15w39。)。 gtext(39。5w39。溫度 39。y軸 /mm39。)。) title(39。,y,z3,39。,y,z2,39。,y,z1,39。 end end plot(y,z,39。sqrt(r(n))w z3(n)=(w^(n)).*q3/(4*b)q3*w^2/(2*b)*log(w/)+a。 else z2(k)=q2*w^2/(2*b)*log(sqrt(r(k))/)+a。 end end for k=1:length(r) 西安工業(yè)大學(xué)畢業(yè)設(shè)計（論文） 30 if sqrt(r(k))wamp。sqrt(r(j))w z1(j)=(w^(j)).*q1/(4*b)q1*w^2/(2*b)*log(w/)+a。 else z(i)=q*w^2/(2*b)*log(sqrt(r(i))/)+a。 for i=1:length(r) if sqrt(r(i))wamp。 %基頻光強為 20w時的熱功率密度 y=0::3。 %基頻光強為 10w時的熱功率密度 q2=p2*c/(pi*w^2)。 %KTP晶體的徑向熱導(dǎo)率 q=p*c/(pi*w^2)。 %基頻光強度 p3=20。 %基頻光強度 p1=10。 %基頻光的束腰半徑 a=15。) KTP 倍頻晶體的功率變化時的溫度分布： c=。 legend(39。5w39。溫度 39。y軸 /mm39。)。) title(39。 end end plot(y,z,39。sqrt(r(i))w z(i)=(w^(i)).*q/(4*b)q*w^2/(2*b)*log(w/)+a。 %倍頻晶體截面在 y軸上的長度 r=(()./10).^2。%KTP冷卻溫度 q=p*c/(pi*w^2)。 %基頻光強度 b=。%KTP晶體對基頻光的吸收率 w=。y軸 /mm39。x軸 /mm39。)。 shading interp colorbar title(39。 else z(i,j)=q*w^2/(2*b)*log(sqrt(r(i,j))/)+a。 if sqrt(r(i,j))wamp。 for i=1:nx。 nx=length(x)。 %倍頻晶體截面在 y軸上的長度 [x,y]=meshgrid(x,y)。 %熱功率密度 x=0::3。 %KTP晶體的徑向熱導(dǎo)率 a=15。 %基頻光的束腰半徑 p=5。) KTP倍頻晶體溫度分布的平面圖： c=。) zlabel(39。) ylabel(39。 xlabel(39。 KTP倍頻晶體的溫度分布 39。 else z(i,j)=q*w^2/(2*b)*log(sqrt(r(i,j))/)+a。 if sqrt(r(i,j))wamp。 for i=1:nx。 nx=length(x)。 %倍頻晶體截面在 y軸上的長度 [x,y]=meshgrid(x,y)。 %熱功率密度 x=0::3。 %KTP晶體的徑向熱導(dǎo)率 a=15。 %基頻光的束腰半徑 p=5。) KTP 晶體的溫度分布 KTP 倍頻晶體的三維溫度分布： c=。 legend(39。5w39。)。 title(39。光強強度 39。)。 xlabel(39。 z=I*exp((2*((x(a/2)).^2+(y(b/2)).^2))/w^2)。 x=。 %倍頻晶體截面在 x、 y軸上的長度 I0=5。%束腰半徑 a=3。西安工業(yè)大學(xué)畢業(yè)設(shè)計（論文） 27 光強的二維分布圖： c=。KTP倍頻晶體的光強分布 39。)。 ylabel(39。x軸 39。 surface(x,y,z)。 [x,y]=meshgrid(x,y)。 x=0::3。 %倍頻晶體截面在 x、 y軸上的長度 I0=5。%束腰半徑 a=3。 光強分布平面圖： c=。KTP倍頻晶體的光強分布 39。)。 zlabel(39。y軸 /mm39。)。 %繪 出光強分布圖 colorbar xlabel(39。 %在 [0,3]*[0,3]區(qū)域生成網(wǎng)格坐標 z=I*exp((2*((x(a/2)).^2+(y(b/2)).^2))/w^2)。 y=0::3。 %基頻光強 I=I0(I0*c/10)。%倍頻晶體截面在 x軸上的長度 b=3。%KTP晶體對基頻光的吸收率 w=。 畢業(yè)設(shè)計（論文）與資料若有不實之處，本人承擔(dān)一切相關(guān)責(zé)任。盡我所知，除了文中特別加以標注和致謝的地方外，畢業(yè)設(shè)計（論文）中不包含其他人已經(jīng)發(fā)表或撰寫過的成果，不包含他人已申請學(xué)位或其他用途使用過的成果。學(xué)校有權(quán)保留送交的畢業(yè)設(shè)計（論文）的原文或復(fù)印件，允許畢業(yè)設(shè)計（論文）被查閱和借閱；學(xué)?？梢怨籍厴I(yè)設(shè)計（論文）的全部或部分內(nèi)容，可以采用影印、縮印或其他復(fù)制手段保存畢業(yè)設(shè)計（論文）。 再一次向所有關(guān)心、幫助我的人表示最誠摯的謝意 ! 畢業(yè)設(shè)計 (論文 )知識產(chǎn)權(quán)聲明 24 畢業(yè)設(shè)計（論文）知識產(chǎn)權(quán)聲明 本人完全了解西安工業(yè)大學(xué)有關(guān)保護知識產(chǎn)權(quán)的規(guī)定，即：本科學(xué)生在校攻讀學(xué)士學(xué)位期間畢業(yè)設(shè)計（論文）工作的知識產(chǎn)權(quán)屬于西安工業(yè)大學(xué)。 劉蓉 老師淵博的學(xué)識、敏銳的思維、民主而嚴謹?shù)淖黠L(fēng)，使學(xué)生收益匪淺，終生難忘。 她在學(xué)習(xí)以及生活上對我的指導(dǎo)和幫助，教會我的不僅是知識，更讓我明白很多為人處事的道理，她 追求真理、獻身科學(xué)、嚴以律己、寬已待人的崇高品質(zhì)對學(xué)生將是永遠的鞭策。emy, et al., Intracavity frequency tripl ing of activelymodelocked diodepumped Nd:YAG laser[J], Appl. Phys. B., 1999, 69(1):8991 致謝 23 致 謝 時間轉(zhuǎn)眼間進入 2022 年，我的大學(xué)生活即將結(jié)束。e, V. Couderc, A. Barth180。 參考文獻 21 參考文獻 [1] 李隆 ,史彭 ,白晉濤等 ,高功率 Nd:YVO4/KTP 腔內(nèi)倍頻晶體溫度分布的半解析熱分析 [J],光 學(xué)技術(shù) , 2022,165(1):4447. [2] 干福熹主編 ,陳創(chuàng)天 ,第 11 章 —非線性光學(xué)晶體材料 ,信息材料 [M], 天津 :天津大學(xué)出版社 ,2022,537598. [3] Fan T Y, Byer R L. Diode laser pumped solidstate lasers[J].IEEE J Quantum Electron, 1988, 24(6): 895912. [4] 姚建銓著 ,非線性光學(xué)頻率變換及激光調(diào)諧技術(shù) [M],北京 :科學(xué)出版社 ,1995,69?108. [5] 李隆 ,劉小建 ,聶建萍 ,等 . 全固態(tài)多模綠光激光器腔內(nèi)倍頻 KTP 晶體熱效應(yīng) ,紅外與激光工程 ,2022,40( 5) :831. [6] 滑文強 . LD 端面泵浦激光熱效應(yīng)研究 .碩士學(xué)位論文 . 山東 :山東師范大學(xué) ,2022. [7] 高俊超 .固體激光介質(zhì)熱效應(yīng)的理論研究 [D].華中科技大 ,2022: 13. [8] 呂百達 .固體激光器件 [M].北京 :北京郵電大學(xué)出版社 ,2022:25. [9] SRIVASTAVA G P, MOHAN S. GUPTA S Tempperature Profile and Inhomgeneity Effect on Second Harmonic Generation using Focuseed Laser Beam[J]. ,1976,23(9):753759. [10] 克希耐爾 .固體激光工程 [M].北京 :科學(xué)出版社 ,1983: 379. [11] 克希耐爾 .固體激光工程 [M].第五版 .北京 :科學(xué)出版社 ,2022:357364. [12] 張玲 .LD 側(cè)面泵浦 Nd:YAG 激光器的熱效應(yīng)研究 [J]. 激光與紅外 ,2022,33(1):3739. [13] 張玲 .LD 側(cè)面泵浦 Nd:YAG 激光器的熱效應(yīng)研究 [J].北方交通大學(xué)學(xué)報 ,2022,26(6):4547. [14] 張行愚 . 激光二極管抽運的激光器熱透鏡效應(yīng)的研究 [J]. 中國激光 ,2022,27(9):7477. [15] ZHIGANG G. Detailed Investigation of Thermal Effects in Longitudinally DiodePumped Nd:YVO4 Laser[J]. IEEE J. ,2022,39(8):979987. [16] 程兆谷 ,程亞 ,王潤溫 ,等 .高功率橫流 CO2激光器高反膜耦合窗口 ZnSe熱形變理論和實驗研究 [J].光子學(xué)報 ,1995,1(5):7882. [17] 盛朝霞 ,李鳳敏 ,郭力 .高功率激光器輸出窗口熱效應(yīng)的數(shù)值分析 [J].激光雜志 ,2022,26(6):2526. [18] 史彭 ,李隆 ,甘安生 . 端面抽運矩形截面 Nd:GdVO4 晶體熱效應(yīng)研究 [J].中國激光 ,2022,32(7):923928. [19] 史彭 ,李隆 ,甘安生 ,等 . LD 斷面抽運矩形截面 YAGNd:YA