写一个磁异常离散傅里叶正逆变换matlab程序，2M+1、dx和2N+1、dy分别为在空间域沿x和y方向的上的采样点数、采样点距，j、k分别为空间域沿x和y方向上的采样点的序号，m、n分别为波数域中沿u、v方向上的采样点的序号，1/Lx和1/Ly分别为波数域u、v方向上的采样间隔

以下是磁异常离散傅里叶正逆变换 Matlab 程序：

% 磁异常离散傅里叶正逆变换程序
% 2M+1、dx和2N+1、dy分别为在空间域沿x和y方向的上的采样点数、采样点距
% j、k分别为空间域沿x和y方向上的采样点的序号
% m、n分别为波数域中沿u、v方向上的采样点的序号
% 1/Lx和1/Ly分别为波数域u、v方向上的采样间隔

% 设置采样点数和采样间隔
M = 100; % x方向采样点数
N = 100; % y方向采样点数
dx = 1; % x方向采样点距
dy = 1; % y方向采样点距
Lx = M*dx; % x方向采样间隔
Ly = N*dy; % y方向采样间隔

% 生成磁异常数据
[X,Y] = meshgrid(-M*dx:dx:M*dx, -N*dy:dy:N*dy);
Z = exp(-(X.^2+Y.^2)/(2*10^2)); % 高斯形状
Z = Z.*cos(2*pi*(X+Y)/100); % 旋转

% 磁异常离散傅里叶正变换
u = (-M:M)/Lx; % u方向采样间隔
v = (-N:N)/Ly; % v方向采样间隔
[U,V] = meshgrid(u,v);
F = fft2(Z)/(2*M+1)/(2*N+1); % 离散傅里叶正变换结果

% 磁异常离散傅里叶逆变换
Z_recon = real(ifft2(F)*(2*M+1)*(2*N+1)); % 离散傅里叶逆变换结果

% 绘制磁异常原始数据和重构数据的图像
figure;
subplot(1,2,1);
imagesc(X(1,:), Y(:,1), Z);
axis equal tight;
colormap(jet);
colorbar;
title('磁异常原始数据');
subplot(1,2,2);
imagesc(X(1,:), Y(:,1), Z_recon);
axis equal tight;
colormap(jet);
colorbar;
title('磁异常重构数据');