代码拉取完成,页面将自动刷新
clear,clc;
%波长(典型值:600)
Lambda = input('输入光的波长(单位为nm):');
Lambda = Lambda * 1e-9;
%缝间隙(典型值:1)
d = input('输入两个缝的间距(单位为mm):');
d = d * 1e-3;
%缝到屏的距离(典型值:1)
D = input('输入缝到屏的距离(单位为m):');
%设置空间坐标的边界
Xmax = 3 * Lambda * D/d;
Ymax = Xmax;
Tmax = Xmax;
%坐标取样1001份
x = linspace(-Xmax, Xmax, 1001);
y = linspace(-Ymax, Ymax, 1001);
t = linspace(-Tmax, Tmax, 1001);
%生成二维网格
[x, y] = meshgrid(x, y);
%计算复振幅
r1 = sqrt((x+d/2).^2+D^2);
r2 = sqrt((x-d/2).^2+D^2);
phi = 2*pi*(r2-r1)/Lambda;
E1 = 1;
E2 = 1 * exp(-1i * phi);
E = E1 + E2;
%合成波光强
I = abs(E) .^ 2;
I = I ./ 4;
%取y轴的中线处的光强
I0 = I(501, :);
%合成光强灰度分布图
subplot(2,1,1);
imshow(I);
xlabel('空间坐标x');
ylabel('空间坐标y');
title('杨氏双缝干涉合成光强空间分布');
% 合成光强函数分布图
subplot(2,1,2);
%设置坐标轴
axis([-Tmax,Tmax,0,1]);
plot(t, I0);
xlabel('空间坐标x'), ylabel('合成光强度');
title('杨氏双缝干涉合成光强空间分布函数');
此处可能存在不合适展示的内容,页面不予展示。您可通过相关编辑功能自查并修改。
如您确认内容无涉及 不当用语 / 纯广告导流 / 暴力 / 低俗色情 / 侵权 / 盗版 / 虚假 / 无价值内容或违法国家有关法律法规的内容,可点击提交进行申诉,我们将尽快为您处理。