clc;clear all; close all

y=0:.25:45;
x=1:.25:60;

[x2,y2] = meshgrid(x,y);

y_dir = sin(2*pi/10*y)';
%  x_dir = ones(1,length(x));
x_dir = sin(2*pi/10*x);

Z = y_dir*x_dir;

figure
subplot(2,1,1)
imagesc(Z)

subplot(2,1,2)
FFT2 = fft2(Z);
imagesc(fftshift(abs(FFT2)))









