Matlab：用空间光调制器产生不同拓扑荷的涡旋光束的相息图

clc
clear all
close all
%%
c=-516:516;
r=-516:516;
l=1;
[x,y]=meshgrid(c,r);
[theta,r]=cart2pol(x,y);
figure;
for l=1:8
    subplot(2,4,l)
    g=mod(l*theta,2*pi);
    imshow(g,[])
    l=l+1;
end
figure;
for l=1:8
    subplot(2,4,l)
    g=mod(-1*l*theta,2*pi);
    imshow(g,[])
    l=l+1;
end
%imshow(mod(-2*theta,2*pi),[])

运行结果：