余弦光栅matlab

余弦光栅是一种一维光栅结构，在MATLAB中可以通过以下代码实现：

% 定义参数
a = 200; % 光栅周期
d = 5; % 光栅宽度
c = 1000; % 图像尺寸
r = 1; % 图像高度

% 创建一个全零的图像
U = zeros(c, r);

% 绘制余弦光栅
for i = 1:a
    for j = 0:(c/d - 1)
        U(i+j*d, :) = cos(2*pi*i/a);
    end
end

% 显示图像
imshow(U);

这段代码中，我们首先定义了光栅周期`a`、光栅宽度`d`、图像尺寸`c`和图像高度`r`。然后创建一个全零的图像`U`。接着使用嵌套的循环来绘制余弦光栅，通过在每个光栅周期内将对应位置的像素值设为余弦函数值。最后使用`imshow`函数显示图像。

相关问题

余弦光栅衍射仿真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');

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

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

好看一点的三维余弦光栅衍射的MATLAB仿真代码

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

clear all;
close all;
clc;

%% 设置参数
wavelength = 0.5e-6; % 波长
k = 2*pi/wavelength; % 波数
L = 10e-3; % 衍射屏尺寸
N = 128; % 格点数
dx = L/N; % 格点间距
z = 1; % 衍射距离
R = 1.5e-3; % 光阑半径
M = 512; % 接收屏尺寸
dy = L/M; % 接收屏格点间距

%% 创建衍射屏
[X,Y] = meshgrid((-N/2:N/2-1)*dx);
R = sqrt(X.^2+Y.^2);
phase = exp(1i*k*z).*exp(-1i*pi*R.^2/z/wavelength);
phase(R>R(end)/2) = 0;

%% 创建接收屏
[Xr,Yr] = meshgrid((-M/2:M/2-1)*dy);
Rr = sqrt(Xr.^2+Yr.^2);
intensity = abs(fftshift(fft2(phase))).^2;

%% 绘图
figure;
subplot(1,2,1);
imagesc(abs(phase));
colormap gray;
axis equal tight;
title('衍射屏');

subplot(1,2,2);
imagesc(Rr,Rr,intensity);
colormap jet;
axis equal tight;
title('接收屏');

该代码首先设置了一些参数，如波长、衍射屏尺寸等。然后创建了衍射屏和接收屏，并计算了光强。最后使用subplot函数绘制了衍射屏和接收屏的图像。