matlab:夫琅禾费衍射简单仿真

clc
clear all
close all
%% 光源部分
N = 200;
Diameter = 10;
x = linspace(-5,5,N); y = linspace(-5,5,N);
[x1,y1] = meshgrid(x,y);
E1 = ones(N);
E1(sqrt(x1.^2+y1.^2)>Diameter/2) = 0;
I1 = E1.*conj(E1);
I1 = I1/(max(max(I1)));
figure;mesh(x1,y1,I1)
set(gca,'fontname','times new roman','fontsize',16);
title('光源光强分布','fontname','华文中宋','fontsize',16);
xlabel('x/mm','fontname','times new roman','fontsize',16);
ylabel('y/mm','fontname','times new roman','fontsize',16);
zlabel('归一化强度','fontname','华文中宋','fontsize',16);
%% 夫琅禾费衍射(单位：mm)
f = 5000;            %%透镜焦距
lambda = 1064e-6;   %%波长
k = 2*pi/lambda;    
z = 1000;
A = 0;  B = f;  C = -1/f;   D = 1-1/f;  %%ABCD定律
xoo = linspace(-1,1,N);  yoo = linspace(-1,1,N);
[x2,y2] = meshgrid(xoo,yoo);
[THETA,R] = cart2pol(x2,y2);
E2 = -1i*k*Diameter^2/8/f*exp(1i*k*z).*exp(1i*k/2/f.*R.^2).*(2*(besselj(1,k*Diameter*R/2*f))./(k*Diameter*R/2*f));
I2 = E2.*conj(E2);  
I2 = I2/max(max(I2));
figure;mesh(x2,y2,I2)
set(gca,'fontname','times new roman','fontsize',16);
title('夫琅禾费衍射后的光强分布','fontname','华文中宋','fontsize',16);
xlabel('x/mm','fontname','times new roman','fontsize',16);
ylabel('y/mm','fontname','times new roman','fontsize',16);
zlabel('归一化强度','fontname','华文中宋','fontsize',16);

运行结果：


参考文献：
[1]吕百达. 激光光学：光束描述、传输变换与光腔技术物理[M]. 高等教育出版社, 2003.