MATLAB源程序代码分享：杨氏双孔干涉实验的MATLAB计算模拟

House

MATLAB源程序代码分享：杨氏双孔干涉实验的MATLAB计算模拟
: L9 q6 D3 S7 j0 h0 S( S%% 定义实验的参数
clear;clc;close all
lambda=3000e-9;     % 光的波长
( f: R9 E: X# p# V5 l- X8 D" r7 F5 p. bd=7e-5;             % 双孔的距离
D=1.2;              % 屏幕距双孔的距离

%% 定义观察屏幕的范围, 并计算各个观测点的坐标值
xa=-2;           % 最小的横坐标值
xb=2;            % 最大的横坐标值
% ?* H2 F: p/ q; ?3 x# L5 h' O( eya=-0.5;         % 最小的纵坐标值
6 M8 f2 N2 I& f, h  d: i6 q/ myb=0.5;          % 最大的纵坐标值
n1=1000;         % x 方向等分份数
6 [9 e; U" F- c. \) F+ yn2=1000;         % y 方向等分份数
' B- o9 t/ c% h6 W9 x
x=linspace(xa,xb,n1);     % x 坐标
6 P+ H( y9 o" S3 f& ~y=linspace(ya,yb,n2);     % y 坐标
[xx,yy]=meshgrid(x,y);    % 网格化处理, 得到观测点的坐标值

; `6 q( b% Z% x# B6 {, [7 s9 B%% 计算各个观测点处的光强
k=2*pi/lambda;                     % 计算波数
r1=sqrt((xx-d/2).^2+yy.^2+D^2);    % 观测点到孔 1 的距离
r2=sqrt((xx+d/2).^2+yy.^2+D^2);    % 观测点到孔 2 的距离
+ H) T: n" }; G7 U# UI=(cos(k*r1)./r1+cos(k*r2)./r2).^2+(sin(k*r1)./r1+sin(k*r2)./r2).^2;   % 计算观测点处的光强
I=I/(max(max(I)));     % 归一化处理
5 D4 b7 d* S% E5 V
$ m! N7 F8 _+ I%% 绘制杨氏双孔干涉条纹
figure
set(gcf,'units','normalized','position',[0.2 0.2 0.6 0.6]);  % 设置 figure 窗口的位置和尺寸
subplot(2,1,1)
image([xa xb],[ya yb],I*255)     % 显示条纹在屏幕上的二维图像
colormap(gray(255))              % 图像显示为黑白, 图像的灰度值为 255 级
2 C0 \: ?1 k' v0 h1 X' \; Jaxis equal
8 D/ D/ P; _9 G0 `9 V2 ~3 h# @axis([xa xb ya yb])
. F/ c. J  F* y1 B. G" Fxlabel('x (m)')
3 ]; P; K. B4 B$ eylabel('y (m)')
title('杨氏双孔干涉条纹')
9 c1 y3 o7 i, s% s8 l
%% 绘制光强分布的三维图
: Y6 e& L+ X3 ]: t/ Y; L! qsubplot(2,1,2)
7 l& I1 `: a) Q) K& G* R9 J; Fmesh(xx,yy,I)
7 x" v/ ~4 Q' a2 |5 _7 I1 \& yxlabel('x (m)')
ylabel('y (m)')
zlabel('光强')
title('光强分布')

wu68aq

杨氏双孔干涉实验的MATLAB计算模拟1 N5。