MATLAB不同光束质量高斯光束的聚焦

%不同光束质量高斯光束的聚焦
clear;clc;
lambda=1.064e-3;
w0=0.2;
R0=1.0e30;
Mf1=[1,0;-1/50,1];
z=linspace(0,300,1000);
Mp2=input('输入光束质量因子Mp2=')%输入光束质量因子
q0=1/(1/R0-j*lambda*Mp2/pi/w0^2);
L1=100; %透镜的位置
wz=zeros(size(z));
for gk=1:1000
if z(gk)<=L1
M=[1,z(gk);0,1];
q=(M(1,1)*q0+M(1,2))/(M(2,1)*q0+M(2,2));
wz(gk)=sqrt(-1/imag(1/q)/pi*lambda*Mp2);
elseif z(gk)>L1
M=[1,z(gk)-L1;0,1]*Mf1*[1,L1;0,1];
q=(M(1,1)*q0+M(1,2))/(M(2,1)*q0+M(2,2));
wz(gk)=sqrt(-1/imag(1/q)/pi*lambda*Mp2);
end
end
plot(z,wz,'b',z,-wz,'b');
title('不同光束质量高斯光束的聚焦 ');
xlabel('z/mm');
ylabel('Wz/mm');
hold on;

输入光束质量因子Mp2=3

Mp2 =

 3

本程序摘自《MATLAB仿真及其在光学课程中的应用》一书第六章