MATLAB Fraunhofer Diffraction Example（Analytic Results）

clear all;clc;

cm = 1e-2;
mm = 1e-3;
um = 1e-6;
nm = 1e-9;
%% 
L = 1*cm;                % 边长 
M = 512;                 % 边长方向上的采样点数 
dx = L/M;                % 采样间隔
x = -L/2:dx:L/2-dx; y=x; % 采样点的坐标
[X,Y] = meshgrid(x,y); 


disp(strcat('x方向的采样间隔(像素尺寸)是：',num2str(dx/um),' um'));
%%
w = 25*um;          % x half-width  这个是圆孔的半径 r  
lambda = 633*nm;    % wavelength    波长
z = 5*cm;           % prop distance 传播距离 m
k =2*pi/lambda;     % wavenumber    
lz = lambda*z; 


NF = w^2/(lambda*z);
disp(strcat('菲涅尔数 NF：',num2str(NF)));
%%
%irradiance 
I2=(w^2/lz)^2.*(jinc(w/lz*sqrt(X.^2+Y.^2))).^2; 
%%
figure(1)  %irradiance image 
imagesc(x./cm,y./cm,nthroot(I2,3)); 
xlabel('x (cm)'); ylabel('y (cm)'); 
colormap('gray'); 
axis square; 
axis xy; 


figure(2) %x-axis profile 
plot(x./cm,I2(M/2+1,:)); 
xlabel('x(cm)'); ylab