Matlab：晶体热透镜焦距粗略计算

      固体激光器在高功率运转的条件下，由于增益介质内部的温度梯度分布会产生不同的折射率，这可以等效为插入一个凸透镜。热透镜在畸变中占据着主导地位，热透镜效应产生的畸变将导致其输出功率的不稳定性和多种模式之间的耦合，同时还限制了激光功率的提高，严重影响了固体激光器的性能。 

      所以如果要设计一款激光器，一定要考虑晶体热透镜对于谐振腔带来的影响。其中热透镜焦距严重影响谐振腔的稳定性等。

    本文只是针对激光晶体的热透镜效应做一个粗略的计算。

%测晶体热透镜效应的焦距
%采用钕离子掺杂浓度为0.5% 
% 尺寸为5x5x10 mm^2的Nd:YVO4的晶体，折射角设为10度

K=0.00523;                 %K 为热导系数
W=0.4;                     %W 为激光晶体泵浦光的平均光斑半径
A=3.14;                    %A 为晶体吸收系数
B=(4.7+-0.6)*10^-6;        %B 为热色散系数dn/dT
C=0.2;                     %C 为泵浦功率转化为热的效率20%
L=10;                      %L 为晶体长度
P=1:0.5:70;                %P 为泵浦功率取值范围
f=(pi*K*W^2)./(P*C*B).*(1./(1-exp(-A*L)));     
plot(P,f)
hold on
axis([0 70 0 500]);        %限制X、Y轴取值范围，（0，70）（0，500）
title('热透镜焦距随泵浦功率变化曲线')
xlabel('泵浦功率[W]');ylabel('热焦距[mm]');