MATLAB仿真：牛顿环干涉实验的光强分布可视化

1. 牛顿环干涉实验的基本原理牛顿环是光学中一个经典的干涉现象最早由艾萨克·牛顿在1665年发现并描述。这个实验不仅展示了光的波动性也是理解薄膜干涉的绝佳案例。想象一下当你把一块平凸透镜放在平面玻璃上时两者之间会形成一个微小的空气间隙。这个间隙的厚度从中心向外逐渐增加就像是一个微型的碗。当单色光垂直照射这个装置时光线会在空气间隙的上下两个表面分别反射。这两束反射光由于来自同一光源具有相干性会在空间中产生干涉。具体来说在空气间隙上表面玻璃-空气界面反射的光会发生半波损失相位突变π在空气间隙下表面空气-玻璃界面反射的光则不会发生相位突变这两束光相遇时根据它们的光程差会产生相长干涉亮环或相消干涉暗环。光程差δ与空气间隙厚度e的关系为δ2eλ/2其中λ是入射光波长。当δ等于整数倍波长时产生相消干涉暗环半整数倍波长时产生相长干涉亮环。2. MATLAB仿真的核心思路用MATLAB模拟牛顿环实验本质上是要建立一个数学模型来描述上述物理过程。我们需要重点关注以下几个关键点2.1 几何关系建模平凸透镜的曲率半径R与空气间隙厚度e之间存在明确的几何关系。对于距离中心r处的点根据圆的几何性质可以得到 e ≈ r²/(2R) 当r远小于R时2.2 光强分布计算干涉光强I可以用以下公式表示 I I₁ I₂ 2√(I₁I₂)cos(Δφ) 其中Δφ是两束光的相位差与光程差的关系为Δφ2πδ/λ对于牛顿环实验通常假设I₁I₂I₀因此公式简化为 I 4I₀cos²(Δφ/2)2.3 可视化策略我们将采用两种可视化方式二维灰度图像模拟实际观察到的牛顿环图案一维光强曲线展示沿径向的光强分布变化3. MATLAB代码实现详解下面我们分段解析完整的MATLAB实现代码帮助理解每个步骤的物理意义和编程逻辑。clc; clear; close all; % 基本参数设置 R 1e1; % 透镜曲率半径(m) Lambda 0.55e-6; % 光波长(m)对应绿光首先清除工作区并设置基本参数。这里选择R10m的大曲率半径模拟实际实验中常用的平凸透镜。波长选择550nm的绿光这是人眼最敏感的可见光波段。% 创建坐标网格 x linspace(-0.01, 0.01, 1001); % 创建-1cm到1cm的采样区间 y x; % y轴与x轴对称 [X,Y] meshgrid(x); % 生成二维网格坐标 r sqrt(X.^2 Y.^2); % 计算各点到中心的距离建立坐标系是仿真的基础。我们创建一个1cm×1cm的方形区域采样1001个点保证分辨率。meshgrid函数生成二维网格坐标然后计算每个点到中心的径向距离r。% 计算干涉光强分布 phase pi*(r.^2/R Lambda/2)/Lambda; I cos(phase).^2; % 光强分布这里直接计算了干涉光强。相位项中的r²/R对应空气间隙厚度带来的相位差Lambda/2项考虑了半波损失。使用cos²函数得到最终的干涉光强分布。% 可视化牛顿环图案 NCLevels 255; % 灰度等级 colormap(gray(NCLevels)); % 设置灰度色图 Ir I*NCLevels; % 将光强映射到灰度值 image(x, y, Ir); % 显示干涉图样 title(牛顿环单位m); axis square; colorbar;将计算得到的光强映射到0-255的灰度值使用image函数显示。设置方形坐标轴和色标使图像更接近实际观察效果。% 绘制径向光强分布曲线 figure; plot(r(501,:), I(501,:), k, LineWidth, 1.5); title(干涉光强分布图); xlabel(r/m环半径); ylabel(相对光强); grid on;通过提取中间行(r0处)的数据绘制径向光强分布曲线。这条曲线清晰地展示了光强随半径的周期性变化。4. 仿真结果分析与优化运行上述代码后我们会得到两个关键结果4.1 牛顿环干涉图样生成的灰度图像显示出一系列明暗相间的同心圆环中心为暗斑这与理论预测完全一致。环间距随着半径增大而逐渐减小这是牛顿环的典型特征。4.2 光强分布曲线曲线展示了光强随半径的振荡衰减。我们可以观察到中心点(r0)光强为0对应暗斑第一个亮环出现在r≈2.3mm处环间距逐渐减小4.3 参数影响分析通过修改代码中的参数可以研究不同条件对牛顿环的影响% 研究波长影响 Lambda [0.45e-6, 0.55e-6, 0.65e-6]; % 蓝、绿、红光 for i 1:3 phase pi*(r.^2/R Lambda(i)/2)/Lambda(i); I cos(phase).^2; subplot(1,3,i); image(x, y, I*255); title([λ,num2str(Lambda(i)*1e9),nm]); end这段代码比较了不同波长下的牛顿环图案。可以看到波长越长环间距越大这与理论公式Δr≈√(Rλ)一致。5. 常见问题与调试技巧在实际仿真过程中可能会遇到一些典型问题5.1 环纹不清晰如果生成的牛顿环边缘模糊可能原因是采样点数不足增加linspace的第三个参数曲率半径太小尝试增大R值显示范围过大减小x,y的范围5.2 出现锯齿状边缘这是由离散化采样导致的解决方法包括增加采样点数使用抗锯齿算法应用图像平滑滤波器% 抗锯齿处理示例 I_smooth imgaussfilt(I, 1); % 高斯平滑 image(x, y, I_smooth*255);5.3 光强曲线异常如果光强曲线不符合预期检查相位计算是否正确是否使用了正确的三角函数cos²而非sin²坐标范围是否合理6. 扩展应用与进阶方向掌握了基础仿真后可以尝试以下扩展6.1 彩色牛顿环模拟用RGB三色光叠加模拟白光干涉% 三色光叠加 Lambda_R 0.65e-6; Lambda_G 0.55e-6; Lambda_B 0.45e-6; I_R cos(pi*(r.^2/R Lambda_R/2)/Lambda_R).^2; I_G cos(pi*(r.^2/R Lambda_G/2)/Lambda_G).^2; I_B cos(pi*(r.^2/R Lambda_B/2)/Lambda_B).^2; Irgb cat(3, I_R, I_G, I_B); image(x, y, Irgb); axis square;6.2 透镜缺陷分析通过引入非理想因素研究其对干涉图样的影响% 添加透镜表面不平整度 noise 0.1e-6*randn(size(r)); % 纳米级随机起伏 phase pi*((r.^2/R noise) Lambda/2)/Lambda;6.3 动态模拟创建干涉图样随参数变化的动画for R linspace(5,20,50) phase pi*(r.^2/R Lambda/2)/Lambda; I cos(phase).^2; image(x, y, I*255); title([R,num2str(R),m]); pause(0.1); end通过这些扩展可以更深入地理解干涉现象也为实际光学测量应用打下基础。在实验室条件受限时这类仿真尤其有价值。