余弦光栅衍射仿真MATLAB代码

以下是一个简单的 MATLAB 代码，用于余弦光栅衍射的仿真：

% 定义光栅参数
p = 10e-6; % 光栅周期
d = p/2; % 光栅线宽
L = 5*p; % 光栅大小
N = 1000; % 采样点数
x = linspace(-L/2, L/2, N); % 生成采样点

% 定义输入波前
wavelength = 633e-9; % 波长
k = 2*pi/wavelength; % 波数
z = 1; % 衍射距离
u0 = exp(1i*k*z) .* exp(1i*k*x.^2/(2*z)); % 输入波前

% 定义光栅函数
g = zeros(size(x));
g(mod(x,p) < d/2) = 1;

% 计算输出波前
u1 = u0 .* g;
u2 = fftshift(fft(u1));
I = abs(u2).^2;

% 绘制结果
figure;
subplot(2,1,1);
plot(x, abs(u1).^2);
title('Input Intensity');
xlabel('Position (m)');
ylabel('Intensity');
subplot(2,1,2);
plot(x, I);
title('Output Intensity');
xlabel('Position (m)');
ylabel('Intensity');

在这个代码中，我们首先定义了光栅的参数，包括光栅周期、光栅线宽、光栅大小和采样点数。然后我们定义了输入波前，这里使用了一个高斯波前。接着我们定义了光栅函数，这里使用了一个简单的方形函数。最后，我们计算了输出波前，并绘制了输入和输出的强度分布。

需要注意的是，这个代码只是一个简单的示例，实际的光栅衍射仿真需要考虑更多的因素，如光栅的方向、光栅的反射率、衍射距离等等。