有一个zemax的zbf的文件，我现在已经可以完成对这个文件的电场数据的读出，代码如下 function [nx, ny, dx, dy, isPol, lambda, EX, EY, PilotRays] = read_ZBF(filename)
% 接收一个文件名，返回采样点数（nx，ny），网格距离（dx，dy），偏振状态（ispol），
% 波长（lambda），电场数据（Ex，Ey），引导光线的数据（PilotRays），
% 以及照度和相位的图像（intensityImage, phaseImage）
% 打开文件
fid = fopen(filename, 'r');
% 读取第一部分（文件格式，采样点数以及偏振状态）
vec = fread(fid, 9, 'int32'); % 假设整数是32位的，根据实际情况可能需要调整
fileFormat = vec(1);
nx = vec(2);
ny = vec(3);
isPol = vec(4);
% 读取第二部分（根据文件格式）
if fileFormat == 1
vec = fread(fid, 20, 'double');
lambda = vec(9);
PilotRays.posX = vec(3);
PilotRays.posY = vec(6);
PilotRays.rayleighX = vec(4);
PilotRays.rayleighY = vec(7);
PilotRays.waistX = vec(5); % 注意：原代码中waist可能是waistX的笔误，但这里保留原样
PilotRays.waistY = vec(8); % 同上
else
vec = fread(fid, 11, 'double');
lambda = vec(5);
% 对于旧格式，假设PilotRays的某些字段是相同的（这里只是示例，可能需要根据实际情况调整）
PilotRays.posX = vec(3);
PilotRays.posY = vec(3); % 通常posY应该有不同的值，但这里只是示例
PilotRays.rayleighX = vec(4);
PilotRays.rayleighY = vec(4); % 通常rayleighY应该有不同的值，但这里只是示例
PilotRays.waistX = vec(6); % 同上，可能是waistX的笔误
PilotRays.waistY = vec(6); % 同上
end
dx = vec(1);
dy = vec(2);
% 读取电场数据EX
vecEX = fread(fid, nx * ny * 2, 'double');
realEX = reshape(vecEX(1:2:end), [ny, nx]);
imagEX = reshape(vecEX(2:2:end), [ny, nx]);
EX = realEX + 1i * imagEX;
EX = EX'; % 转置以匹配通常的图像格式（行表示y，列表示x）
% 如果光束是偏振的，读取电场数据EY
if isPol == 1
vecEY = fread(fid, nx * ny * 2, 'double');
realEY = reshape(vecEY(1:2:end), [ny, nx]);
imagEY = reshape(vecEY(2:2:end), [ny, nx]);
EY = realEY + 1i * imagEY;
EY = EY'; % 转置
else
EY = zeros(size(EX));
end我现在需要对这个二维的矩阵进行傅里叶变换，我有这个zbf文件的EX和EY电场分布，我需要对这个矩阵进行复电场的描述，然后傅里叶变换，得到矩阵中每一个波数和能量之间的关系，然后我还需要得到角度和能量之间的关系，注意这个角度的话，应该既不是与x轴的夹角，也不是与y轴的夹角，他的夹角应该是与xy平面垂直的z轴的夹角，写一个matlab的代码
|
