[imgid=0]理想的光学涡旋，或者纯粹的光学涡旋：

clc

clear

close all;

M=512;N=512;  [x,y]=meshgrid(-M/2+1:M/2,-N/2+1:N/2);

%============生成光学涡旋相位图==========%

theta=atan2(x,y);

L=30;

r=sqrt(x.^2+y.^2);

a=0.01;

vortex=exp(1i*L*theta).*cyl(x,y,N);

% vortex=(sqrt(2)*r/a)^L*exp(-r.^2/a)*exp(-1i*L*theta);

vortex_I=abs(vortex).^2;

vortex_phase=angle(vortex);

vortex_counter=vortex/max(max(vortex));%归一化

figure(1)

imshow(vortex_phase);

F=fftshift(fft2(vortex_counter));

F1=abs(F)/max(max(abs(F)));

figure(2)

imshow(F1);

figure(3)

plot(r,F1);