单缝弗朗禾费衍射光强公式代码MATLAB

以下是单缝弗朗禾费衍射光强公式的MATLAB代码实现：

% 定义常量
lambda = 632.8e-9;    % 波长
b = 1e-3;             % 单缝宽度
L = 1;                % 距离
k = 2*pi/lambda;      % 波数

% 定义计算函数
f = @(theta) (sin(k*b/2*sin(theta))./(k*b/2*sin(theta))).^2;
I = @(theta) f(theta).*cos(theta);

% 定义计算范围
theta = linspace(-pi/2, pi/2, 10000);

% 计算光强分布
I_theta = I(theta);

% 绘制图像
plot(theta, I_theta);
xlabel('θ (rad)');
ylabel('I(\theta)');
title('Single Slit Fraunhofer Diffraction Pattern');

在运行此代码时，请注意将波长、单缝宽度和距离调整为您需要计算的值。