基模高斯光束的MATLAB仿真模拟

光电GGBond

已于 2024-01-12 20:39:32 修改

阅读量1.8k

点赞数 21

文章标签： matlab 算法

于 2023-12-23 22:14:07 首次发布

本文链接：https://blog.csdn.net/jinverynb/article/details/135173896

版权

基模高斯光束是各种高斯光束的基础，其复振幅分布呈高斯分布。

下面我将对运用matlab模拟基模高斯光束进行讲解：

公式部分

首先需要了解一些高斯光束的基本参量：

光束的输入波长 --------------------------- $\lambda$ （lambda）
光束的传播距离---------------------------- $z$
传播常数，又称波数---------------------- $k=2\pi/\lambda$
束腰半径------------------------------------- $w_0$
常用：r(86.5%)=1；r(99%)=1.5；r(50%)=0.6，最常用1/e^2半径，r(x%)表示包含光束总功率x%的半径相比于1/e^2半径的比值。
r瑞利长度或共焦参数--------------------- $z_R=kw_0^2/2$
光斑范围半径，传播到z处的束宽---- $w(z)=w_0\cdot \sqrt{1+(z/z_R)^2}$
高斯球面波的波前曲率半径------------ $R(z)=z+z_R^2/z$
古伊相位------------------------------------ $\vartheta (z)=arctan(z/z_R)$

确定了基本参量后便可以带入基模高斯光束的的完整电场表达式进行求解。

公式如下：

光强为电场强度与其复共轭的乘积：

代码部分

编程时需要先实现基本参量的输入及简单计算：

%% 基本参量
lambda=input('入射光波波长(mm)=');  %光束的输入波长mm
z=input('传播距离=');               %光束传播距离
k=2*pi/lambda;                     %传播常数，又称波数
w0=2;                              %高斯光束的中心束腰半径，这里设为2mm
ZR=k*w0^2/2;                       %瑞利长度或共焦参数
wz=w0*sqrt(1+(z/ZR)^2);            %光斑范围半径，传播到z处的束宽
Rz=z+ZR*ZR/z;                      %高斯球面波的波前曲率半径
Phi=atan(z/ZR);                    %古伊相位

然后建立坐标系（极坐标系）：

%% 建立极坐标y
L=10;                              %x，y方向范围，[-L,L]
N=1001                             %取样数（必须为奇数）
x=linspace(-L,L,N);                %x方向分割
y=linspace(-L,L,N);                %y方向分割
[x1,y1]=meshgrid(x,y);             %生成二维网格坐标
[phi,r]=cart2pol(x1,y1);           %转化成极坐标

根据公式表达基模高斯光束：

%% 基模高斯光束的表达式
E_Gauss=(w0./wz).*exp(-(r.^2./wz^2))...
.*exp(1i.*(k.*z+k.*r.^2./(2*Rz)-Phi));%基模高斯光束的完整的电场表达式
I_Gauss=E_Gauss.*conj(E_Gauss);       %基模高斯光束的光强
normI_Gauss=I_Gauss/max(max(I_Gauss));%光强的归一化处理（normalization process）

最后绘制高斯光束的图像：

figure

subplot(2,2,1)
plot(x,normI_Gauss((N+1)/2,:),'linewidth',1.5); %归一化强度一维图像
set(gca,'fontname','times new roman','fontsize',15); %坐标轴设置
xlabel('\itx','fontname','times new roman','fontsize',20);
ylabel('归一化强度','fontname','宋体','fontsize',20); 

subplot(2,2,2)                                      
pcolor(x1,y1,normI_Gauss);                           %归一化强度二维图像           
set(gca,'fontname','times new roman','fontsize',15);     
xlabel('\itx','fontname','times new roman','fontsize',20);
ylabel('\ity','fontname','times new roman','fontsize',20); 
clim([0,1]);                                         %色彩映射范围
shading interp;
colormap hot;
colorbar; 

subplot(2,2,3)
surf(x1,y1,normI_Gauss);                             %归一化强度三维图像
set(gca,'fontname','times new roman','fontsize',15);     
xlabel('\itx','fontname','times new roman','fontsize',20); 
ylabel('\ity','fontname','times new roman','fontsize',20);
clim([0,1]);
shading interp;
colormap hot;
colorbar; 

theAxes=axis;