1. 球面波干涉法检测
clc; clear;
N = 2000;
lambda = 632e-6;
k=2*pi/lambda;
z = 0; w0 = 0.5; p = 0;
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);
theta =atan2(Y,X);
r = sqrt(X.^2 + Y.^2);
beam = exp(1i * theta * m);
beam0 = exp(-1i*k/200*r);
I = (beam + beam0).*conj(beam + beam0);
figure();
imagesc(I);
colormap gray;
axis off;
2. 平面波干涉法检测
clc; clear;
N = 2000;
lambda = 632e-6;
k=2*pi/lambda;
z = 0; w0 = 0.5; p = 0;
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);
theta =atan2(Y,X);
r = sqrt(X.^2 + Y.^2);
beam = exp(1i * theta * m);
beam1 = exp(-1i*k*x);
I = (beam + beam1).*conj(beam + beam1);
figure();
imagesc(I)
colormap gray;
axis off;
3. 共轭干涉法检测
3.1 函数"lg_oam.m"
function [LG_OAM]=lg_oam(L,w0,P, lamda, x,y)
[px,py]=meshgrid(x,y);
r=sqrt(px.^2+py.^2);
phai=atan2(py,px);
z=0;
k0=2*pi./lamda;
zR=pi*w0^2/lamda;
wz=w0.*sqrt((z.^2+zR.^2)./zR.^2);
Clp=sqrt(2.*factorial(P)./(pi.*factorial(P+abs(L))));
X=2*(r.^2)./wz^2;
LG_poly=0;
for n=1:P+1
l=n-1;
LG_poly=LG_poly+(-1)^l*((factorial(P+abs(L)))/((factorial(P-l))*(factorial(abs(L)+l))*(factorial(l)))).*X.^l;
end
E1=wz^(-1).*(sqrt(2)*sqrt(r.^2)/wz).^abs(L).*LG_poly.*exp(-(r.^2)./wz^2);
E2=exp((-1i*k0*(r.^2)*z)/(2*(z^2+zR^2)));
E3=exp(-1i.*(2.*P+abs(L)+1).*atan(z/zR));
LG_OAM=Clp.*E1.*E2.*E3.*exp(-1i.*L.*phai);
end
3.2 函数"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
3.3 示例
clc; clear;
N = 1000;
[lambda,k,z,w0,x,y] = beam_parameter(N);
beam1 = lg_oam(2, w0, 0, lambda, x, y);
beam2 = lg_oam(-2, w0, 0, lambda, x, y);
I = (beam1 + beam2).*conj(beam1 + beam2);
figure();
imagesc(I);
colormap grap;
axis off;
4. Reference
- 陈瑞山. 涡旋光束的高效产生与检测. 2018. 中国科学技术大学, PhD dissertation.