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

House

: Z$ F! I* {5 W2 y. JMATLAB源程序代码分享：杨氏双孔干涉实验的MATLAB计算模拟
' J' B* U# g, [+ w! H. v4 a0 V$ P%% 定义实验的参数
  ]( d, f* O7 w8 i+ A9 a. G7 @clear;clc;close all
8 E$ f! ]' s2 i3 Ulambda=3000e-9;     % 光的波长
d=7e-5;             % 双孔的距离
4 K6 ~6 c. E4 QD=1.2;              % 屏幕距双孔的距离
' i8 n8 w0 g1 G; o' H
%% 定义观察屏幕的范围, 并计算各个观测点的坐标值
$ k+ J. O/ p2 Oxa=-2;           % 最小的横坐标值
xb=2;            % 最大的横坐标值
ya=-0.5;         % 最小的纵坐标值
1 x7 Z  h/ V+ C& iyb=0.5;          % 最大的纵坐标值
n1=1000;         % x 方向等分份数
n2=1000;         % y 方向等分份数

  x9 A: L) U8 g2 t4 i6 ~x=linspace(xa,xb,n1);     % x 坐标
y=linspace(ya,yb,n2);     % y 坐标
) P' I/ x" N# h; ]: C8 l) z[xx,yy]=meshgrid(x,y);    % 网格化处理, 得到观测点的坐标值
2 v+ d& D1 H- C# I1 E
9 |1 [  i/ j+ ?# E% b%% 计算各个观测点处的光强
; C- ]% I, P. |k=2*pi/lambda;                     % 计算波数
r1=sqrt((xx-d/2).^2+yy.^2+D^2);    % 观测点到孔 1 的距离
7 X7 k4 a$ d$ F  |r2=sqrt((xx+d/2).^2+yy.^2+D^2);    % 观测点到孔 2 的距离
I=(cos(k*r1)./r1+cos(k*r2)./r2).^2+(sin(k*r1)./r1+sin(k*r2)./r2).^2;   % 计算观测点处的光强
2 Z' i2 U! h  U% [' D% WI=I/(max(max(I)));     % 归一化处理
2 J+ `' L4 W8 Z8 D, m+ ~% Z6 G9 S% h
. q" T. S" v& x: j%% 绘制杨氏双孔干涉条纹
, U8 z* K0 q+ [, w7 {figure
) _4 [6 m1 @3 O3 ~set(gcf,'units','normalized','position',[0.2 0.2 0.6 0.6]);  % 设置 figure 窗口的位置和尺寸
subplot(2,1,1)
9 I+ }; @, x0 O1 C# N0 I/ Oimage([xa xb],[ya yb],I*255)     % 显示条纹在屏幕上的二维图像
: w; a- V; i2 o; R9 M2 rcolormap(gray(255))              % 图像显示为黑白, 图像的灰度值为 255 级
axis equal
axis([xa xb ya yb])
0 f3 B5 ^% [' F3 dxlabel('x (m)')
ylabel('y (m)')
title('杨氏双孔干涉条纹')
% p2 {* V, Q; [) W2 W
! X5 F1 D" l1 A1 i% C%% 绘制光强分布的三维图
  c# J0 \9 y% Osubplot(2,1,2)
mesh(xx,yy,I)
6 ]8 M( s: S+ N$ Nxlabel('x (m)')
ylabel('y (m)')
zlabel('光强')
& c. |, K9 S  Ptitle('光强分布')

wu68aq

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