【正文】 分布 39。)。 xlabel(39。x軸 /mm39。) ylabel(39。y軸 /mm39。) zlabel(39。溫度分布 39。) KTP倍頻晶體溫度分布的平面圖： c=。%KTP晶體對基頻光的吸收率 w=。 %基頻光的束腰半徑 p=5。 %基頻光強度 b=。 %KTP晶體的徑向熱導率 a=15。%KTP冷卻溫度 q=p*c/(pi*w^2)。 %熱功率密度 x=0::3。 %倍頻晶體截面在 x軸上的長度 y=0::3。 %倍頻晶體截面在 y軸上的長度 [x,y]=meshgrid(x,y)。 %在 [0,]*[0,]區(qū)域生成網(wǎng)格坐標 r=(()./10).^2+(()./10).^2。 nx=length(x)。 ny=length(y)。 for i=1:nx。 for j=1:ny。 if sqrt(r(i,j))wamp。sqrt(r(i,j))w z(i,j)=(w^(i,j)).*q/(4*b)q*w^2/(2*b)*log(w/)+a。 else z(i,j)=q*w^2/(2*b)*log(sqrt(r(i,j))/)+a。 end end end surface(x,y,z)。 shading interp colorbar title(39。KTP倍頻晶體的溫度分布 39。)。 xlabel(39。x軸 /mm39。) ylabel(39。y軸 /mm39。) KTP倍頻晶體的二維溫度分布： c=。%KTP晶體對基頻光的吸收率 w=。 %基頻光的束腰半徑 p=5。 %基頻光強度 b=。 %KTP晶體的徑向熱導率 a=15。%KTP冷卻溫度 q=p*c/(pi*w^2)。 %熱功率密度 y=0::3。 %倍頻晶體截面在 y軸上的長度 r=(()./10).^2。 西安工業(yè)大學畢業(yè)設(shè)計（論文） 29 for i=1:length(r) if sqrt(r(i))wamp。sqrt(r(i))w z(i)=(w^(i)).*q/(4*b)q*w^2/(2*b)*log(w/)+a。 else z(i)=q*w^2/(2*b)*log(sqrt(r(i))/)+a。 end end plot(y,z,39。b39。) title(39。x=a/2處 KTP倍頻晶體的溫度分布 39。)。 xlabel(39。y軸 /mm39。) ylabel(39。溫度 39。) gtext(39。5w39。)。 legend(39。功率 5w39。) KTP 倍頻晶體的功率變化時的溫度分布： c=。%KTP晶體對基頻光的吸收率 w=。 %基頻光的束腰半徑 a=15。%KTP冷卻溫度 p=5。 %基頻光強度 p1=10。 %基頻光強度 p2=15。 %基頻光強度 p3=20。 %基頻光強度 b=。 %KTP晶體的徑向熱導率 q=p*c/(pi*w^2)。 %基頻光強為 5w時的熱功率密度 q1=p1*c/(pi*w^2)。 %基頻光強為 10w時的熱功率密度 q2=p2*c/(pi*w^2)。 %基頻光強為 15w時的熱功率密度 q3=p3*c/(pi*w^2)。 %基頻光強為 20w時的熱功率密度 y=0::3。 r=(()./10).^2。 for i=1:length(r) if sqrt(r(i))wamp。sqrt(r(i))w z(i)=(w^(i)).*q/(4*b)q*w^2/(2*b)*log(w/)+a。 else z(i)=q*w^2/(2*b)*log(sqrt(r(i))/)+a。 end end for j=1:length(r) if sqrt(r(j))wamp。sqrt(r(j))w z1(j)=(w^(j)).*q1/(4*b)q1*w^2/(2*b)*log(w/)+a。 else z1(j)=q1*w^2/(2*b)*log(sqrt(r(j))/)+a。 end end for k=1:length(r) 西安工業(yè)大學畢業(yè)設(shè)計（論文） 30 if sqrt(r(k))wamp。sqrt(r(k))w z2(k)=(w^(k)).*q2/(4*b)q2*w^2/(2*b)*log(w/)+a。 else z2(k)=q2*w^2/(2*b)*log(sqrt(r(k))/)+a。 end end for n=1:length(r) if sqrt(r(n))wamp。sqrt(r(n))w z3(n)=(w^(n)).*q3/(4*b)q3*w^2/(2*b)*log(w/)+a。 else z3(n)=q3*w^2/(2*b)*log(sqrt(r(n))/)+a。 end end plot(y,z,39。b39。,y,z1,39。r39。,y,z2,39。g:39。,y,z3,39。k.39。) title(39。不同功率下 x=a/2處 KTP倍頻晶體的溫度分布 39。)。 xlabel(39。y軸 /mm39。) ylabel(39。溫度 39。) gtext(39。5w39。)。 gtext(39。10w39。)。 gtext(39。15w39。)。 gtext(39。20w39。)。 legend(39。功率 5w39。,39。功率 10w39。,39。功率 15w39。,39。功率 20w39。) KTP 倍頻晶體的基頻光半 徑變化時的溫度分布： c=。 %KTP晶體對基頻光的吸收率 a=15。%KTP冷卻溫度 w=。 %基頻光的束腰半徑 w1=。 %基頻光的束腰半徑 w2=。 %基頻光的束腰半徑 w3=。 %基頻光的束腰半徑 p=5。 %基頻光強度 b=。 %KTP晶體的徑向熱導率 q=p*c/(pi*w^2)。 %基頻光強為 5w時的熱功率密度 q1=p*c/(pi*w1^2)。 %基頻光強為 5w時的熱功率密度 q2=p*c/(pi*w2^2)。 %基頻光強為 5w時的熱功率密度 q3=p*c/(pi*w3^2)。 %基頻光強為 5w時的熱功率密度 y=0::3。 r=(()./10).^2。 for i=1:length(r) if sqrt(r(i))wamp。sqrt(r(i))w z(i)=(w^(i)).*q/(4*b)q*w^2/(2*b)*log(w/)+a。 else z(i)=q*w^2/(2*b)*log(sqrt(r(i))/)+a。 end 西安工業(yè)大學畢業(yè)設(shè)計（論文） 31 end for j=1:length(r) if sqrt(r(j))w1amp。sqrt(r(j))w1 z1(j)=(w1^(j)).*q1/(4*b)q1*w1^2/(2*b)*log(w1/)+a。 else z1(j)=q1*w1^2/(2*b)*log(sqrt(r(j))/)+a。 end end for k=1:length(r) if sqrt(r(k))w2amp。sqrt(r(k))w2 z2(k)=(w2^(k)).*q2/(4*b)q2*w2^2/(2*b)*log(w2/)+a。 else z2(k)=q2*w2^2/(2*b)*log(sqrt(r(k))/)+a。 end end for n=1:length(r) if sqrt(r(n))w3amp。sqrt(r(n))w3 z3(n)=(w3^(n)).*q3/(4*b)q3*w3^2/(2*b)*log(w3/)+a。 else z3(n)=q3*w3^2/(2*b)*log(sqrt(r(n))/)+a。 end end plot(y,z,39。b39。,y,z1,39。r39。,y,z2,39。g:39。,y,z3,39。k.39。) title(39。不同基頻光半徑下 x=a/2處 KTP倍頻晶體的溫度分布 39。)。 xlabel(39。y軸 /mm39。) ylabel(39。溫度 39。) gtext(39。39。)。 gtext(39。39。)。 gtext(39。39。)。 gtext(39。39。)。 legend(39。激光半徑為 39。,39。激光半徑為 39。,39。激光半徑為 39。, 39。激光半徑為 39。) KTP 倍頻晶體冷卻溫度變化時的溫度 分布： c=。%KTP晶體對基頻光的吸收率 w=。 %基頻光的束腰半徑 p=5。 %基頻光強度 b=。 %KTP晶體的徑向熱導率 q=p*c/(pi*w^2)。 %熱功率密度 y=0::3。 r=(()./10).^2。 for i=1:length(r) if sqrt(r(i))wamp。sqrt(r(i))w z(i)=(w^(i)).*q/(4*b)q*w^2/(2*b)*log(w/)+15。 else 西安工業(yè)大學畢業(yè)設(shè)計（論文） 32 z(i)=q*w^2/(2*b)*log(sqrt(r(i))/)+15。 end end for j=1:length(r) if sqrt(r(j))wamp。sqrt(r(j))w z1(j)=(w^(j)).*q/(4*b)q*w^2/(2*b)*log(w/)+17。 else z1(j)=q*w^2/(2*b)*log(sqrt(r(j))/)+17。 end end for k=1:length(r) if sqrt(r(k))wamp。sqrt(r(k))w z2(k)=(w^(k)).*q/(4*b)q*w^2/(2*b)*log(w/)+20。 else z2(k)=q*w^2/(2*b)*log(sqrt(r(k))/)+20。 end end for n=1:length(r) if sqrt(r(n))wamp。sqrt(r(n))w z3(n)=(w^(n)).*q/(4*b)q*w^2/(2*b)*log(w/)+23。 else z3(n)=q*w^2/(2*b)*log(sqrt(r(n))/)+23。 end end plot(y,z,39。b39。,y,z1,39。g39。,y,z2,39。r39。,y,z3,39。c39。) axis([0 0 28]) title(39。不同冷卻溫度下 x=a/2處 KTP倍頻晶體的溫度分布 39。)。 xlabel(39。y軸 /mm39。) ylabel(39。溫度 39。) gtext(39。15度 39。)。 gtext(39。17度 39。)。 gtext(39。20度 39。)。 gtext(39。23度 39。)。 BBO 晶體的熱效應(yīng) BBO倍頻晶體的三維溫度分布 ： c=。%BBO晶體對基頻光的吸收率 w=。 %基頻光的束腰半徑 p=5。 %基頻光強度 a=15。%BBO冷卻溫度 b=。 %BBO晶體的徑向熱導率 q=p*c/(pi*w^2)。 %熱功率密度 x=0::3。 %倍頻晶體截面在 x軸上的長度 西安工業(yè)大學畢業(yè)設(shè)計（論文） 33 y=0::3。 %倍頻晶體截面在 y軸上的長度 [x,y]=meshgrid(x,y)。 %在 [0,]*[0,]區(qū)域生成網(wǎng)格坐標 r=(()./10).^2+(()./10).^2。 nx=length(x)。 ny=length(y)。 for i=1:nx。 for j=1:ny。 if sqrt(r(i,j))wamp。sqrt(r(i,j))w z(i,j)=(w^(i,j)).*q/(4*b)q*w^2/(2*b)*log(w/)+a。 else z(i,j)=q*w^2/(2*b)*log(sqrt(r(i,j))/)+a。 end end end mesh(x,y,z) colorbar title(39。 BBO倍頻晶體的溫度分布 39。)。 xlabel(39。x軸 /mm39。) ylabel(39。y軸 /mm39。) zlabel(39。溫度分布 39。) BBO倍頻晶體的平面溫度分布： c=。%BBO晶體對基頻光的吸收率 w=。 %基頻光的束腰半徑 p=5。%基頻光強度 b=。 %BBO晶體的徑向熱導率 a=15。%BBO冷卻溫度 q=p*c/(pi*w^2)。 %熱功率密度 x=0::3。 %倍頻晶體截面在 x軸上的長度 y=0::3。 %倍頻晶體截面在 y軸上的長度 [x,y]=meshgrid(x,y)。 %在 [0,]*[0,]區(qū)域生成網(wǎng)格坐標 r=(()./10).^2+(()./