Hi Ben,
This may be a stupid question, but in one over your papers you define the wave vector
as
nx=(-Nx/2:Nx/2-1)';
ny=(-Ny/2:Ny/2-1)';
[x,y]=meshgrid(nx*dx,ny*dy);
[kx,ky]=meshgrid(nx*2*pi/(Nx*dx),ny*2*pi/(Ny*dy));
k=sqrt(kx.^2+ky.^2);
can you explain why you did this? i don't quite understand how it works.
thank you, jeremy
Hi Jeremy,
Please can you be a bit more specific about what it is you don't understand. Is it that there are two slightly different uses of the word 'vector' here?
kx and ky are vectors in the Matlab sense of 'vector' meaning a list of numbers. They contain a discretised version of the wavenumber components in the x and y directions of the wavevector (kx,ky), where wavevector here means 'vector' in the mathematical sense meaning a quantity with direction and magnitude. k is the magnitude of the wavevector (kx,ky).
The relationship between the discretised spatial coordinates (x,y) and the discretised wavenumber components (the coordinates in wavenumber space) (kx,ky) is fixed by the requirements of the Fast Fourier Transform.
Hope that helps,
Ben
You must log in to post.