Matlab：实现菲涅尔圆孔衍射仿真

最新推荐文章于 2024-09-13 11:31:42 发布

Dream_言十

最新推荐文章于 2024-09-13 11:31:42 发布

阅读量1.3w

点赞数 14

分类专栏： Matlab_光学文章标签： matlab 物理光学

本文链接：https://blog.csdn.net/weixin_43487914/article/details/119042145

版权

Matlab_光学专栏收录该内容

23 篇文章 242 订阅

订阅专栏

本文介绍了如何使用Matlab基于光学原理，通过两种方法（菲涅尔衍射积分的离散求和及FFT法）来仿真菲涅尔圆孔衍射现象。代码包括高斯光束产生、圆孔掩膜创建以及菲涅尔积分算法的实现，展示了衍射前后的图像，并提供了相关函数的参考代码。

摘要由CSDN通过智能技术生成

本文基于光学原理，用Matlab实现菲涅尔圆孔衍射现象仿真。

以下使用了两种菲涅尔衍射的数值计算方法：

菲涅尔衍射积分的离散求和
菲涅尔衍射FFT法

Project Code

% 作者：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;