MATLAB高斯光束透镜焦面衍射数值计算

最新推荐文章于 2024-03-02 17:28:02 发布

「已注销」

最新推荐文章于 2024-03-02 17:28:02 发布

阅读量278

点赞数

分类专栏： MATLAB 数学文章标签： matlab 开发语言

本文链接：https://blog.csdn.net/m0_38127487/article/details/130279807

版权

MATLAB 同时被 2 个专栏收录

57 篇文章 0 订阅

订阅专栏

数学

37 篇文章 0 订阅

订阅专栏

%高斯光束透镜焦面衍射数值计算
clear;
clc;
tic;
n=1.5062;%1064nm 波长折射率，k9 玻璃
d=3; % 透镜中心厚度
RL=0.025e3;% 透镜凸面曲率半径
f=RL/(n-1);%透镜的焦距
R0=1;%入射光束半径
lambda=1.064e-3;
k=2*pi/lambda;
phy=lambda/pi/R0;
z=f;
mr2=51;
ne2=61;
mr0=81;
while sqrt(R0^2+z^2)-sqrt(R0^2*(1-1/mr0)^2+z^2)>lambda/20
mr0=mr0+1;
end
ne0=mr0;
rmax=5*f*phy;
r=linspace(0,rmax,mr2);
eta=linspace(0,2*pi,ne2);
[rho,theta]=meshgrid(r,eta);
[x,y]=pol2cart(theta,rho);
r0=linspace(0,3*R0,mr0);
eta0=linspace(0,2*pi*(ne0-1),ne0-1);
[rho0,theta0]=meshgrid(r0,eta0);
[x0,y0]=pol2cart(theta0,rho0);
deta=3*R0/(mr0-1)*2*pi/(ne0-1);
E2=zeros(size(x));
E1=exp(-(x0.^2+y0.^2)/R0^2);
for gk=1:ne2
for df=1:mr2
Rrho=sqrt((x(gk,df)-x0).^2+(y(gk,df)-y0).^2+z^2);
Rtheta=z./Rrho;
opd=exp(j*k*((n-1)*(sqrt(RL^2-rho0.^2)-(RL-d))+d));
Ep=-j/lambda/2*exp(Rrho*j*k).*(1+Rtheta)./Rrho*deta.*rho0.*opd.*E1;
E2(gk,df)=sum(Ep(:));
end
end
Ie=conj(E2).*E2;
figure;
surf(x,y,Ie);
shading interp;
grid off;
box on;
toc;