pingao_with_thermaleffect.m

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

wl=1.054e-3;         % wl为波长，单位为mm
r=10;                % r为输出镜的曲率半径
z=2;                 % z为晶体离全反镜的距离
l=0.001:0.1:r-0.001;  % l为腔长，单位为mm。    
%这里l从0.001开始到r-0.001结束，是为了不让l和r-l等于0，否则第9句分母imag(1./qz)会出现有0的情况
w0=sqrt(wl.*sqrt(l.*(r-l))./pi);
q0=i*pi.*w0.^2./wl;
qz=q0+z;
wz=sqrt(-1*wl/pi./imag(1./qz));
plot(l,wz,'r')
xlabel('腔长l(mm)')
ylabel('晶体上的光斑w(z) (mm)')
%ylabel('腰斑w0 (mm)')
r=38;
l=0.001:0.1:r-0.001;
w0=sqrt(wl.*sqrt(l.*(r-l))./pi);
q0=i*pi.*w0.^2./wl;
qz=q0+z;
wz=sqrt(-1*wl/pi./imag(1./qz));
plot(l,wz,'g')