运用对比度法测光学调制传毒函数matlab仿真代码

以下是运用对比度法测光学调制传递函数的 Matlab 仿真代码：

% 定义物理常数
lambda = 633e-9; % 波长
k = 2*pi/lambda; % 波数
NA = 0.2; % 数值孔径
a = 5.5e-6; % 光栅常数
f = 1e-3; % 透镜焦距
L = 2*a*f/lambda/NA; % 半长轴

% 定义光栅函数
N = 1024; % 光栅大小
x = linspace(-a/2, a/2, N);
[X, Y] = meshgrid(x);
grating = 0.5*(1+sign(X));

% 定义透镜函数
R = 0.5*f/NA; % 透镜半径
r = linspace(-R, R, N);
[Rho, Phi] = meshgrid(r, r);
lens = double(Rho <= R);

% 定义传递函数
P = fft2(ifftshift(grating.*lens));
H = exp(1j*k*L*sqrt(1-(lambda*Rho/f).^2));
Q = P.*H;

% 计算幅度传递函数
Q1 = abs(Q);
Q2 = abs(fft2(lens));
C = Q1./Q2;

% 显示结果
figure;
subplot(221); imshow(grating, []); title('光栅函数');
subplot(222); imshow(lens, []); title('透镜函数');
subplot(223); imshow(log10(abs(H)), []); title('相位传递函数');
subplot(224); imshow(C, []); title('幅度传递函数');

代码首先定义了一些物理常数，包括波长、波数、数值孔径、光栅常数、透镜焦距等。然后定义了一个光栅函数，使用 sign 函数实现了一个周期为 $a$ 的矩形波形。接着定义了一个透镜函数，使用一个半径为 $R$ 的圆形窗口来限制透镜的作用范围。通过傅里叶变换计算出传递函数 $Q$，并计算出幅度传递函数 $C=|Q|/|H|$，最后在一个图像窗口中显示结果。