本文基于光学原理,用Matlab实现菲涅尔圆孔衍射现象仿真。
以下使用了两种菲涅尔衍射的数值计算方法:
% 作者:ZQJ
% 日期:2021.7.22 星期四
%**********************基于光场传输模拟菲涅尔圆孔衍射*************************
%% 菲涅尔衍射积分*******************************************
clear,clc,close all;
allmasks = f_Masks;
func_translight = f_Lightfield_transmission;
% ********************** 系统参数设定
lamda = 1550e-9; % 波长
N = 200; % 一边像素数目
x_length = 0.01; % 矩阵长度
w0 = 1e-3; % 高斯光束腰半径
Gs_z = 1e-3; % 高斯光传输距离
% ********************** 菲涅尔衍射
mask_hole = allmasks.circle_hole(x_length/N,N,0.08/N);
figure; subplot(1,2,1), imagesc(mask_hole); colormap gray;
Gaussian_I = f_Gaussian_beams(w0,Gs_z,lamda,x_length,x_length,N,N);
subplot(1,2,2),imagesc(Gaussian_I), colormap hot;
E0 = Gaussian_I.*mask_hole;
E1 = func_translight.Fresnel_integrate(lamda,E0,x_length,x_length,x_length,x_length,N,N,0.5);
figure,imagesc(abs(E1)), colormap hot;
%% 菲涅尔衍射FFT法******************************************
clear,clc,close all;
allmasks = f_Masks;
func_translight = f_Lightfield_transmission;
% ********************** 系统参数设定
lamda = 1550e-9; % 波长
N = 200; % 一边像素数目
x_length = 0.01; % 矩阵长度
w0 = 1e-3; % 高斯光束腰半径
Gs_z = 1e-3; % 高斯光传输距离
% ********************** 菲涅尔衍射
mask_hole = allmasks.circle_hole(x_length/N,N,0.08/N);
figure; subplot(1,2,1), imagesc(mask_hole); colormap gray;
Gaussian_I = f_Gaussian_beams(w0,Gs_z,lamda,x_length,x_length,N,N);
subplot(1,2,2),imagesc(Gaussian_I), colormap hot;
E0 = Gaussian_I.*mask_hole;
E1 = func_translight.FFT_(lamda,E0,x_length,x_length,0.5);
figure,imagesc(abs(E1)), colormap hot;
专栏内容供作者本人或大家学习使用,多多指教 ~