Performing the transform on the below logic

I have the below function for which I need to perform spatial Fourier transform and plot its spatial frequency spectrum for various values of y.

How to do this using matlab?

Y1 = acos(x/2+y) + asin(g/3+x);

%g is constant

[AndresVar]

for various values of y & x?

you can create a meshgrid of y's and x's

evaluate the function f(x,y)

then use fft2: 2-D fast Fourier transform - MATLAB fft2 (mathworks.com)

Answer 1

[AndresVar]

Make a grid, evaluate the function and use fft2 to get the spatial spectrum.

Here is an example, note function might be different fom what you posted.

Note fftshift was used to shift the spectrum to the center, you can try without.

g=10; % for example
x=linspace(1,4,10); % for example
y=linspace(1,4,10); % for example

[X,Y]=meshgrid(x,y); % a grid

F = (atan(g/2+X)./Y+atan(g/2-X)./Y)/pi; % evaluate the grid

F_PS = abs(fftshift(fft2(F))); % power spectrum, shifted 

imagesc(F_PS);

[Amy Topaz]

Thank you