Matlab菲涅尔基尔霍夫衍射公式

weixin_50538161

已于 2024-06-29 22:26:57 修改

阅读量354

点赞数 7

分类专栏： Matlab 光束文章标签： matlab 开发语言

于 2024-06-29 22:25:54 首次发布

本文链接：https://blog.csdn.net/weixin_50538161/article/details/140071737

版权

光束同时被 2 个专栏收录

5 篇文章 0 订阅

订阅专栏

Matlab

4 篇文章 0 订阅

订阅专栏

1. 代码

clear all;
close all;
lambda = 632.8e-9; % 光波的波长
k = 2 * pi / lambda; % 波数
L = 10e-3; % 采样区域大小单位m
N = 1024;% 采样点数
a = 1e-3; % 光阑半径单位m
z = 1; % 衍射距离

% 定义计算[函数]
fresnel_diffraction = @(x, y) (1 / (1i * lambda * z)) * exp(1i * (2 * pi / lambda * z + pi / lambda * (x^2 + y^2) / z)) * exp(1i * pi / lambda * (x^2 + y^2) / z);

% 定义计算区域
% 生成采样点
x = linspace(-L/2, L/2, N);
y = linspace(-L/2, L/2, N);
[X, Y] = meshgrid(x, y);

% 计算衍射场
diffraction_field = fresnel_diffraction(X, Y);

% 绘制衍射图像
m = 4;
for i = 1:m
    figure(i);
    I = abs(diffraction_field).^2 .* (2 * besselj(i, k * a * sqrt(X.^2 + Y.^2) / z) ./ (k * a * sqrt(X.^2 + Y.^2) / z)).^2;

    imagesc(x, y, I);
    axis equal tight;
    colormap('jet');
    xlabel('x (m)');
    ylabel('y (m)');
    title([num2str(i),'阶besselj光束'],'菲涅尔基尔霍夫衍射');
end