pingao_thermalense_jisuanw.m

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

%本程序与caculate_thermalense.m程序配合使用，本程序是第一个程序
clear;
wl=1.033e-3;          % wl为Yb:GdCOB的波长，单位为mm
n=1.74;              % Yb:GdCOB的折射率

r=50;                % r为输出镜的曲率半径
l1=1;                % l1为晶体离入射全反镜的距离
lc=2;                %晶体长2mm
l3=7.5;
syms f;

T1=[1,l1
    0,1];
T2=[1,lc/n/2
    0,1];
T3=[1,0
    -1/f,1];
T4=[1,lc/n/2
    0,1];
T5=[1,l3
    0,1];
T6=[1,0
    -2/r,1];

M=T3*T4*T5*T6*T5*T4*T3*T2*T1*T1*T2;
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))