闪耀型径向相位光栅检测OAM模式值

最新推荐文章于 2024-07-30 11:28:17 发布

鸡蛋牛肉面

最新推荐文章于 2024-07-30 11:28:17 发布

阅读量874

点赞数 25

分类专栏： matlab 文章标签： matlab

本文链接：https://blog.csdn.net/SmithCGauss/article/details/139219570

版权

matlab 专栏收录该内容

5 篇文章 1 订阅

订阅专栏

1. OAM涡旋光束生成

1.1 拉盖尔-高斯(LG)光束的复振幅表达式:

$U_{pl}^{LG}=A_{LG}\times \sqrt{\frac{2p!}{\pi (p+|l|)!}}\frac{1}{w(z)}\left[\frac{r\sqrt{2}}{w(z)}\right]^{|l|}L_{p}^{|l|}\left[\frac{2r^2}{w^2(z)}\right]\times exp\left[\frac{-r^2}{w^2(z)} \right]\times exp\left[\frac{-ikr^2z}{2 (z^2+z_R^2)} \right]exp\left[i(2p+|l|+1)tan^{-1}\frac{z}{z_R} \right]exp(-il\phi)$
其中 $(r,\phi,z)$ 为柱坐标参数， $A_{LG}$ 为振幅， $l$ 为OAM模式值，表示空间相位的变化， $p$ 为径向指数， $w(z)=w_0\sqrt{1+(z/z_R)^2}$ 是LG光束在 $z$ 处的束腰半径， $w_0$ 是基模高斯光束在 $z = 0$ 处的束腰半径， $w_0=\pi w_0^2/\lambda$ 是瑞丽距离， $\lambda$ 为光束波长， $k=2\pi /\lambda$ 是光束的波数， $(2p+|l|+1)tan^{-1}\frac{z}{z_R}$ 是Gouy相位， $L_{p}^{|l|}$ 为归一化的拉盖尔多项式，其表达式为：
$L_{p}^{|l|}=(-1)^{|l|}\frac{d^{|l|}}{dx^{|l|}}L_{p+|l|}(x)$

1.2 代码

1.2.1 文件"generate_OAM.m"

function result = generate_OAM(mode, w0, z, p, lanbda,X,Y)
    k = 2*pi/lanbda;                     %波数
    Z_R = pi*w0^2/lanbda;                %瑞利长度
    w_z = w0*sqrt(1+(z/Z_R)^2);          %光束在z位置的半径
    [theta,r] = cart2pol(X,Y);
    if mode == 0
        result = zeros(1000);
    else
        result = sqrt(2*factorial(p)/pi/(p+factorial(abs(mode))))*...
            (1/w_z)*(sqrt(2)*r/w_z).^abs(mode).*exp(-r.^2/w_z^2).*...
            laguerre(p,abs(mode),2*r.^2/w_z^2).*exp(-1i*mode*theta).*...
            exp(-1i*k*z).*exp(-1i*k*r.^2*z/2/(z^2+Z_R^2))*...
            exp(-1i*(2*p+abs(mode)+1)*atan(z/Z_R));
    end
end

1.2.2 文件"laguerre.m"

%% 拉盖尔多项式
function result = laguerre(p,l,x)
    result = 0;
    if p == 0
        result = 1;
    elseif p == 1
        result = 1+abs(l)-x;
    else
        result = (1/p)*((2*p+l-1-x).*laguerre(p-1,abs(l),x)-(p+l-1)*laguerre(p-2,abs(l),x));
    end
end

1.2.3 文件"beam_parameter.m"

function [lambda,k,z,w0,x,y] = beam_parameter(N)
    lambda = 0.532e-6;
    k = 2 * pi / lambda; %波数
    z = 0; 
    w0 = 0.005;
    T0 = 4;
    R1 = T0 * w0; %范围
    x = linspace(-R1, R1, N);
    y = linspace(-R1, R1, N);
end

1.3 示例

clc; clear;
N = 1000;
[lambda,k,z,w0,x,y] = beam_parameter(N);

%产生OAM光束
beam = generate_OAM(1,w0,0,0,lambda,x,y) + generate_OAM(4,w0,0,0,lambda,x,y);

%显示光强
figure();
imagesc(x,y,abs(beam).^2)
colormap gray;
axis off;

2. 夫琅禾费远场衍射

$U(x,y)=\frac{exp(ikz)}{i\lambda z}\cdot exp\left[i\frac{k}{2z}(x^2+y^2) \right]\cdot F\left[u(\xi,\eta)\times t(\xi,\eta) \right]$
其中, $(x, y)$ 为坐标平面, $F(\cdot)$ 为傅立叶变换, $u(\xi,\eta)$ 为光束光场, $t(\xi,\eta)$ 为光栅透过率函数.

3. 闪耀型径向相位光栅

3.1 光栅透过率函数

$t(\theta)=exp(i2\pi \theta /P)$
其中, $P$ 是光栅角向周期, $\theta$ 为方位角.

3.2 代码

%%光栅检测法 —— 闪耀型径向相位光栅 (单位：mm)
clc; clear;
N = 2000;

%定义光栅及光束平面
dxs = 10e-3; dys = dxs;              %光源平面采样间隔
R1 = dxs*N/2;                        %范围
x = linspace(-R1/2, R1/2, N);
y = linspace(-R1/2, R1/2, N);
[X, Y] = meshgrid(x, y);

L = R1/20*26;                        %光束偏离中心的距离
theta =atan2(Y, X+L);                %角度
z0 = 5e3;                            %观察屏与光栅的距离

%涡旋光束参数
lambda = 532.8e-6;                   %波长
k=2*pi/lambda;                       %波数
z = 0; w0 = 0.232; p = 0;

%定义观察屏平面
R2 = lambda*z0/dxs;                  %范围
a = linspace(-R2/2, R2/2, N);
b = linspace(-R2/2, R2/2, N);
[A, B] = meshgrid(a, b);

%产生闪耀型径向相位光栅
P = 0.001;
grating = exp(1i*2*pi*theta/P);      %由透过率函数得出

%产生涡旋光束
mode = -5;                           %光束模态, 拓扑荷数
beam1 = generate_OAM(mode, w0, z, p, lambda, X, Y);

%夫琅禾费积分衍射
U1 = fftshift(fft2(grating.*beam1)); %傅里叶变换法
U1 = U1.*dxs*dys;                    %补偿连续傅里叶变换和FFT差值
U1 = U1.*exp(1i*k./(2*z0).*(A.^2+B.^2)).*exp(1i*k*z0)./(1i*z0*lambda);

figure('Name','衍射光斑');
imagesc(abs(U1).^2);
title(['l=' num2str(mode)]);
colormap hot;
axis off;