写一个磁异常离散傅里叶正逆变换matlab程序，2M+1、dx和2N+1、dy分别为在空间域沿x和y方向的上的采样点数、采样点距

以下是磁异常离散傅里叶正逆变换的 MATLAB 程序，其中2M+1和2N+1分别为在空间域沿x和y方向的采样点数，dx和dy为采样点的距离。

%% 磁异常离散傅里叶正逆变换程序

% 输入磁异常数据
data = input('请输入磁异常数据：');

% 输入采样点数和采样点距离
dx = input('请输入采样点距离dx：');
dy = input('请输入采样点距离dy：');
M = input('请输入采样点数M：');
N = input('请输入采样点数N：');

% 计算采样区域大小
Lx = 2 * M + 1;
Ly = 2 * N + 1;

% 计算频域采样点数
kx = -M:M;
ky = -N:N;

% 计算频域采样点间隔
dkx = 2 * pi / (Lx * dx);
dky = 2 * pi / (Ly * dy);

% 计算频域网格
[kx_grid, ky_grid] = meshgrid(kx * dkx, ky * dky);

% 计算磁异常离散傅里叶变换
f = fftshift(fft2(data)) * dx * dy;

% 计算磁异常离散傅里叶逆变换
data_recon = ifft2(ifftshift(f)) / (dx * dy);

% 显示结果
figure;
subplot(1, 2, 1);
imagesc(data);
title('原始数据');
xlabel('x');
ylabel('y');
colorbar;
axis equal tight;
subplot(1, 2, 2);
imagesc(real(data_recon));
title('逆变换结果');
xlabel('x');
ylabel('y');
colorbar;
axis equal tight;

% 保存结果
save('磁异常离散傅里叶正逆变换结果.mat', 'data', 'data_recon', 'dx', 'dy', 'M', 'N');