clear;
wl=1.05e-3;         % wl為Yb:Na:CaF2波長(zhǎng)，單位為mm

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

r=25;                % r為輸出鏡的曲率半徑
l1=1;                % l1為晶體離全反鏡的距離
lc=9;                %晶體長(zhǎng)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;        %此時(shí)算的是晶體離泵浦光最近的右端面
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;        %此時(shí)算的是晶體離泵浦光最遠(yuǎn)的左端面
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))