用狭缝法测光学调制传递函数matlab代码

以下是用狭缝法测光学调制传递函数的MATLAB代码：

% 读取原始图像
I = imread('lena.png');

% 转换为灰度图像
I = rgb2gray(I);

% 定义狭缝尺寸及间距
slit_width = 10; % 狭缝宽度
slit_spacing = 20; % 狭缝间距

% 定义狭缝位置
slit_positions = 1:slit_spacing:size(I, 1);

% 计算垂直狭缝的传递函数
M = zeros(size(I));
for i = 1:length(slit_positions)
    slit_position = slit_positions(i);
    slit = zeros(size(I));
    slit(slit_position:slit_position+slit_width-1, :) = 1;
    slit_fft = fftshift(fft2(slit));
    I_fft = fftshift(fft2(I));
    M_fft = I_fft ./ slit_fft;
    M_slit = abs(ifft2(ifftshift(M_fft)));
    M(slit_position:slit_position+slit_width-1, :) = M_slit(slit_position:slit_position+slit_width-1, :);
end

% 显示结果
figure;
subplot(1, 3, 1);
imshow(I);
title('原始图像');
subplot(1, 3, 2);
imshow(M);
title('狭缝法测得的OTF');
subplot(1, 3, 3);
imshow(log(1+abs(M)), []);
title('log(1+abs(OTF))');

需要注意的是，这里计算的是垂直狭缝的传递函数，如果需要计算水平狭缝的传递函数，可以将代码中的行和列互换。另外，这里计算的是光学传递函数的幅值，如果需要计算相位信息，可以将代码中的`abs`函数去掉。