帮我编写matlab计算图像mtf曲线函数

好的，下面是使用MATLAB计算图像MTF曲线的函数：

function [freq, mtf] = calcMTF(image)
% image：输入的灰度图像

% 提取灰度图像中的条纹数据
[pixel_num, ~] = size(image);
stripes = image(pixel_num/2-4:pixel_num/2+4, :);
stripe_mean = mean(stripes);

% 对条纹数据进行傅里叶变换，得到频域信息
fft_result = abs(fft(stripe_mean));
fft_result = fft_result(1:floor(length(fft_result)/2));

% 计算频率轴
freq = (0:length(fft_result)-1) / length(fft_result);

% 计算MTF曲线
mtf = fft_result ./ fft_result(1);

% 平滑处理MTF曲线
mtf = smooth(mtf, 5);

% 归一化处理MTF曲线
mtf = mtf / max(mtf);
end

使用方法：

1. 将以上函数保存为.m文件。

2. 准备一张包含灰度条纹的图像。

3. 载入图像，并将其传递给函数进行计算。

   image = imread('test.jpg');
   image = rgb2gray(image);
   [freq, mtf] = calcMTF(image);

4. 可以将计算得到的频率轴和MTF曲线进行可视化。

   plot(freq, mtf);
   xlabel('frequency (cycles/pixel)');
   ylabel('MTF');

希望以上代码可以帮助到你。