帮我生成根据边缘扩展函数得到MTF的matlab程序

边缘扩展函数在图像处理中常用于计算像片的调制传递函数（Modulation Transfer Function, MTF），它描述了图像系统的分辨率性能。以下是一个简单的MATLAB程序，用于从边缘扩展函数计算MTF：

% 假设边缘扩展函数为EdgeResponse，单位通常为弧度
EdgeResponse = ... % 请替换为实际的边缘响应数据

% 定义像素大小和最大角度范围
pixelSize = ...; % 像素尺寸，单位通常是弧度/meter
maxAngle = pi / 2; % 角度范围上限

% 确定频率间隔 (假设是等间隔)
numSamples = 1000; % 分辨率采样点数
freqInterval = maxAngle / numSamples;

% 计算空间频率范围
spaceFreqRange = (-maxAngle : freqInterval : maxAngle) * pixelSize;

% 使用傅里叶变换计算MTF
MTF = abs(fftshift(fft(EdgeResponse))) ./ (length(EdgeResponse) / 2);

% 窗口函数（如Hanning窗，有助于减少高频噪声）
MTF = hanning(length(MTF)) .* MTF;

% 对结果取对数并归一化到0-1范围
MTF = log10(max(MTF(:))) - log10(min(MTF(:)));

% 打印或保存MTF结果
disp(['MTF: ' mat2str(MTF)]);
save('MTF_data.mat', 'MTF', 'spaceFreqRange');

%