vxingqiang.m

激光器典型腔计算

% 在本程序中，腔由入射腔镜（平面镜），晶体（不考虑热透镜效应），输出镜（凹面，r＝50mm）和SESAM构成
% 晶体的厚度为2mm，对于11mm的短腔来说长度不可忽略，若晶体的热透镜焦距为f，
% 在我的实验的腔中存在两个个束腰，即平面镜处和靠近SESAM处
% 起点从晶体中间开始和SESAM算起
r=50;            % r为输出镜的曲率半径
l1=1.5;          % l1为入射全反镜与晶体表面之间的距离，在我的实验中为1.5mm
l2=19.5;           % l2为晶体表面与凹面镜之间的距离，在我的实验中为24mm
l3=23.2;           % l3为SESAM离输出镜的距离
lm=2;            % lm为晶体的长度，在我的实验中为2mm
wl=1.033e-3;     % 激光波长为1.033um
n=1.74;          % Yb:GdCOB的折射率

T1=[1,l3
    0,1];
T2=[1,0
    -2/r,1];
T3=[1,l2
    0,1];
T4=[1,lm/n
    0,1];
T5=[1,l1
    0,1];
T6=[1,lm/2/n
    0,1];           % T6是晶体一半厚度的传输矩阵

M=T1*T2*T3*T4*T5*T5*T4*T3*T2*T1;
ADhalf=(M(1,1)+M(2,2))/2
wsesam=sqrt(wl/pi*abs(M(1,2))/sqrt(-1*M(1,2)*M(2,1)-((M(1,1)-M(2,2))/2)^2))

Mc=T6*T5*T5*T4*T3*T2*T1*T1*T2*T3*T6;
ADhalfc=(Mc(1,1)+Mc(2,2))/2
wcrytal=sqrt(wl/pi*abs(Mc(1,2))/sqrt(-1*Mc(1,2)*Mc(2,1)-((Mc(1,1)-Mc(2,2))/2)^2))