pingaoqiang_longcrystal.m

来自「激光器典型腔计算」· M 代码 · 共 33 行

clear;
wl=1.05e-3;         % wl为Yb:Na:CaF2波长，单位为mm

n=1.74;              %注意这是Yb:GdCOB的折射率，我不知道Yb:Na:CaF2的折射率才用别的晶体的折射率代替的

r=25;                % r为输出镜的曲率半径
l1=1;                % l1为晶体离全反镜的距离
lc=9;                %晶体长9mm
syms l3;

T1=[1,l1
    0,1];
T2=[1,lc/n
    0,1];
T3=[1,l3
    0,1];
T4=[1,0
    -2/r,1];
T5=[1,lc/n/4
    0,1];
T6=[1,lc/n/4*3
    0,1];
M=T5*T1*T1*T2*T3*T4*T3*T6;
ADhalf=(M(1,1)+M(2,2))/2
w=sqrt(wl/pi*abs(M(1,2))/sqrt(-1*M(1,2)*M(2,1)-((M(1,1)-M(2,2))/2)^2))

Mr=T1*T1*T2*T3*T4*T3*T2;        %此时算的是晶体离泵浦光最近的右端面
ADhalfr=(Mr(1,1)+Mr(2,2))/2;
wr=sqrt(wl/pi*abs(Mr(1,2))/sqrt(-1*Mr(1,2)*Mr(2,1)-((Mr(1,1)-Mr(2,2))/2)^2))

Ml=T2*T1*T1*T2*T3*T4*T3;        %此时算的是晶体离泵浦光最远的左端面
ADhalfl=(Ml(1,1)+Ml(2,2))/2;
wl=sqrt(wl/pi*abs(Ml(1,2))/sqrt(-1*Ml(1,2)*Ml(2,1)-((Ml(1,1)-Ml(2,2))/2)^2))