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

House

0 r& A. p9 z9 S* {! X; u- b  m: U% |MATLAB源程序代码分享：杨氏双孔干涉实验的MATLAB计算模拟
%% 定义实验的参数
clear;clc;close all
lambda=3000e-9;     % 光的波长
d=7e-5;             % 双孔的距离
. o* |. M4 o1 @- r5 b- G3 YD=1.2;              % 屏幕距双孔的距离

%% 定义观察屏幕的范围, 并计算各个观测点的坐标值
, M* j& p( T( }4 G8 zxa=-2;           % 最小的横坐标值
xb=2;            % 最大的横坐标值
! g. ]1 k2 E: _% [3 aya=-0.5;         % 最小的纵坐标值
% i5 i8 W: t; c: m0 Q7 syb=0.5;          % 最大的纵坐标值
n1=1000;         % x 方向等分份数
n2=1000;         % y 方向等分份数

x=linspace(xa,xb,n1);     % x 坐标
y=linspace(ya,yb,n2);     % y 坐标
, N7 _5 C6 {  {1 L0 ^[xx,yy]=meshgrid(x,y);    % 网格化处理, 得到观测点的坐标值

%% 计算各个观测点处的光强
k=2*pi/lambda;                     % 计算波数
9 h& J& Y1 J2 [% o: T! ~$ ]$ Yr1=sqrt((xx-d/2).^2+yy.^2+D^2);    % 观测点到孔 1 的距离
! J! }5 C0 {7 z$ `2 }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;   % 计算观测点处的光强
  s. e9 U1 ^9 u4 d0 h2 d* GI=I/(max(max(I)));     % 归一化处理
; x; z, p5 R; j" K* d
7 b( [0 y& z- u' g; r& p%% 绘制杨氏双孔干涉条纹
+ l9 Z" S! D( Z  }figure
set(gcf,'units','normalized','position',[0.2 0.2 0.6 0.6]);  % 设置 figure 窗口的位置和尺寸
$ @  S# t; M, Z$ R& zsubplot(2,1,1)
0 G& h- j, _9 c3 Bimage([xa xb],[ya yb],I*255)     % 显示条纹在屏幕上的二维图像
colormap(gray(255))              % 图像显示为黑白, 图像的灰度值为 255 级
. Z2 H6 S/ r0 paxis equal
axis([xa xb ya yb])
xlabel('x (m)')
ylabel('y (m)')
title('杨氏双孔干涉条纹')

%% 绘制光强分布的三维图
subplot(2,1,2)
mesh(xx,yy,I)
xlabel('x (m)')
ylabel('y (m)')
7 P  `  r9 W" w) Zzlabel('光强')
title('光强分布')

wu68aq

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