【光学】模拟拉盖尔高斯【含Matlab源码 2167期】

⛄一、获取代码方式

获取代码方式1：
完整代码已上传我的资源：【光学】基于matlab模拟拉盖尔高斯【含Matlab源码 2167期】
点击上面蓝色字体，直接付费下载，即可。

获取代码方式2：
付费专栏Matlab物理应用（初级版）

备注：
点击上面蓝色字体付费专栏Matlab物理应用（初级版），扫描上面二维码，付费29.9元订阅海神之光博客付费专栏Matlab物理应用（初级版），凭支付凭证，私信博主，可免费获得1份本博客上传CSDN资源代码（有效期为订阅日起，三天内有效）；
点击CSDN资源下载链接：1份本博客上传CSDN资源代码

⛄二、拉盖尔高斯光束简介

LGp, l光束是一种发散的高斯光束, 是亥姆霍兹方程在缓变振幅近似下的一个特解, 具体表达式为

式 (1) 中: (r, φ, z) 为场点的柱坐标;k为波数;Ap, l为振幅常数;w0为z=0处的光斑半径;w为拉盖尔-高斯光束的光斑半径;R为拉盖尔-高斯光束

p与亮光环的个数有关, 文中研究p=0, 仅有1个亮光环的情况;l为拉盖尔多项式的阶数, 与光环的半径成正比关系, LGp, l光束具有轨道角动量, 可以旋转被捕获的微粒.随着LGp, l光束阶数l的增加, 角动量也会增大.利用高阶LGp, l光束可获得较大的角动量, 但是阶数l的增加可能会导致捕获效率的下降.因此, 阶数的选取很重要.

⛄三、部分源代码

%%%LG光束的光强
clear
N = 512;
delta =0.003;
% c=-1:0.01:1;

% [x,y]=meshgrid(c,r);
[theta,rho]=cart2pol(x,y);
w0=0.3; %束腰半径
m=4;p=0;z=0; %确定光束的阶数m,模数p，传输距离z
lambda=0.632e-3;k=2pi/lambda;
zr=piw0^2/lambda;
wz=w0sqrt(1+(z/zr)^2);
A=2factorial§/(pifactorial(m+p)); %归一化因子
u=(sqrt(2)rho/w0).^m.*exp(-rho.2/w0^{2).*exp(1i*m*theta)*sqrt(2/factorial(m)/pi).*Laguerre(p,m,2*rho.}2/w0^2);
%aaa= A/wz^2.*(2*rho.2/wz^2).m.(Laguerre(p,m,2rho.^2/wz2)).^{2.*exp(-2*rho.}2/wz^2).exp(1im*theta); %光场表达式
beta=10;
I2=u.*conj(u);
Ie=I2/max(max(I2));
%画光强相位
figure
surf(x,y,Ie);
shading interp %色彩插值处理
zlabel(‘intensity’);

figure
pha=rem(angle(u)+pi2,2pi);
surf(x,y,pha);
shading interp
%colormap gray;

%画光强相位像素图
figure
imshow(Ie)
figure
vortex_phase=rem(angle(u)+2pi,2pi);
imshow(vortex_phase,[])

⛄四、运行结果

⛄五、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1]门云阁.MATLAB物理计算与可视化[M].清华大学出版社，2013.
[2]王娟,任洪亮,周业鹏.拉盖尔-高斯光束光镊捕获性质[J].华侨大学学报(自然科学版). 2014,35(06)

3 备注
简介此部分摘自互联网，仅供参考，若侵权，联系删除