pingaoqiang.m

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

wl=1.04e-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)')
hold on;

r=15;
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,'m')
hold on;

r=25;
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,'c')
hold on;

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')
hold on;

r=50;
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,'b')
%hold on;

%r=100;
%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,'k')