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),[])
运行结果: