Pattern8

Shanshan Wang
on 22 Oct 2022
- 9
- 14
- 1
- 0
- 173
n = 2^10;
M = zeros(n);
I = 1:n;
x = I-n/2;
y = n/2-I;
[X,Y] = meshgrid(x,y);
R = 10;
A = (X.^2 +Y.^2 <= R^2);
M(A) = 5;
DP = fft2(fftshift(M));
imagesc(abs(log2(DP)))
colormap jet
axis image off