假如你是天线测试领域的算法工程师，假如你有了平面近场测试数据，请给出详细matlab程序的近远场变换算法

近场数据经过FFT变换后可以得到频域数据，采用快速傅里叶变换FFT可以快速地将离散时域信号转换成离散频域信号。然后可以通过一个数学变换将频域数据转换为远场数据，这个过程就是近场到远场的变换。下面是一个简单的matlab程序实现：

%% 近场到远场变换
clc;
clear all;
 
%% 参数设置
f = 10e9; % 频率10GHz
c = 3e8; % 光速
lambda = c / f; % 波长
d = 2*lambda; % 接收天线与测量平面的距离
k = 2*pi/lambda; % 波数
L = 2; % 测量平面边长
Nx = 256; % 测量平面x方向上的采样点数
Ny = 256; % 测量平面y方向上的采样点数
x = linspace(-L/2,L/2,Nx); % 测量平面x方向上的采样点坐标
y = linspace(-L/2,L/2,Ny); % 测量平面y方向上的采样点坐标
[X,Y] = meshgrid(x,y); % 测量平面上的网格坐标
dx = x(2)-x(1); % x方向上的采样间距
dy = y(2)-y(1); % y方向上的采样间距
 
%% 读取近场数据
data = load('nearfield_data.txt');
Ex = reshape(data(:,1),Nx,Ny).'; % 近场电场x分量
Ey = reshape(data(:,2),Nx,Ny).'; % 近场电场y分量
Ez = reshape(data(:,3),Nx,Ny).'; % 近场电场z分量
 
%% FFT变换
fEx = fft2(Ex); % 近场电场x分量的频域数据
fEy = fft2(Ey); % 近场电场y分量的频域数据
fEz = fft2(Ez); % 近场电场z分量的频域数据
 
%% 近场到远场变换
kx = 2*pi/L*[-Nx/2:Nx/2-1]; % x方向上的波数
ky = 2*pi/L*[-Ny/2:Ny/2-1]; % y方向上的波数
[Kx,Ky] = meshgrid(kx,ky); % 波数网格坐标
kz = sqrt(k^2 - Kx.^2 - Ky.^2); % z方向上的波数
fHx = 1j*kz./k.*(Kx.*fEx + Ky.*fEy); % 远场磁场x分量的频域数据
fHy = 1j*kz./k.*(Kx.*fEy - Ky.*fEx); % 远场磁场y分量的频域数据
fHz = -kz./k.*fEz; % 远场磁场z分量的频域数据
 
%% 逆FFT变换
Hx = ifft2(fHx); % 远场磁场x分量的时域数据
Hy = ifft2(fHy); % 远场磁场y分量的时域数据
Hz = ifft2(fHz); % 远场磁场z分量的时域数据
 
%% 计算远场强度
E = sqrt(abs(Hx).^2 + abs(Hy).^2 + abs(Hz).^2); % 远场电场强度
E_dB = 10*log10(E/max(max(E))); % 远场电场强度的dB值
 
%% 画图
figure;
imagesc(x,y,E_dB);
axis equal tight;
colorbar;
xlabel('x/m');
ylabel('y/m');
title('远场电场强度分布图');

这个程序实现了平面近场测试数据的近远场变换算法，并可以画出远场电场强度分布图。需要注意的是，这个程序仅适用于平面波测试，对于球面波测试或其他形状的测量平面需要进行适当的修改。