⛄一、获取代码方式
获取代码方式1:
完整代码已上传我的资源:【光学】基于matlab拉盖尔高斯光束光场仿真【含Matlab源码 2471期】
点击上面蓝色字体,直接付费下载,即可。
获取代码方式2:
付费专栏Matlab物理应用(初级版)
备注:
点击上面蓝色字体付费专栏Matlab物理应用(初级版),扫描上面二维码,付费29.9元订阅海神之光博客付费专栏Matlab物理应用(初级版),凭支付凭证,私信博主,可免费获得1份本博客上传CSDN资源代码(有效期为订阅日起,三天内有效);
点击CSDN资源下载链接:1份本博客上传CSDN资源代码
⛄二、部分源代码
%%%LG光束的光强
clear
N = 512;
delta =0.003;
% c=-1:0.01:1;
% r=-1:0.01:1;
x=(-N/2:N/2-1)delta;
y=(-N/2:N/2-1)delta;
[x,y]=meshgrid(x,y);
% [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/w02).*exp(1i*m*theta)*sqrt(2/factorial(m)/pi).*Laguerre(p,m,2*rho.2/w0^2);
%aaa= A/wz2.*(2*rho.2/wz2).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.
3 备注
简介此部分摘自互联网,仅供参考,若侵权,联系删除