简介杨氏双缝干涉MATLAB仿真资源面向光学工程与应用物理专业学生以及MATLAB仿真入门者专注于解决干涉实验的教学演示与可视化问题。压缩包共3个文件包含2个.m脚本和1个.fig图形界面整体仅63KB结构简练。脚本部分负责构建双缝干涉模型基于波动理论逐点计算观察屏上的相位差与光强分布并绘制出明暗相间的干涉条纹GUI部分则提供直观的交互界面支持通过滑动条动态调节波长、缝间距等关键参数实时刷新图案帮助理解条纹间距和亮度随条件的变化。资源在光强计算中引入菲涅尔积分的数值实现思路完整涵盖了从光源与双缝参数设定、干涉图样计算到结果可视化的标准流程可作为工程光学课程设计或光学现象仿真演示的参考模板。目前已有2036人学习浏览对于希望复现经典杨氏双缝实验并快速上手MATLAB编程与界面设计的学习者这份资料提供了低门槛、高实用性的入门路径。1. 杨氏双缝干涉的MATLAB仿真从光程差公式到条纹图的完整落地做工程光学课程设计或者项目复现时“杨氏双缝干涉”听起来是个再经典不过的题目可真正上手用MATLAB把它仿真出来不少人第一步就卡住了——光强公式写在纸上没问题一进代码就不知道坐标网格怎么建、条纹方向为什么转了90度、暗纹位置为什么总差半个周期。这套杨氏双缝干涉MATLAB仿真资源不是只给你一个能出图的脚本而是把干涉条件、光程差计算、条纹可见度、参数扫描这套完整链路都拆开了新手能照着一步步跑出和教材一致的干涉图熟手也能直接拿里面的参数化函数去改波长、缝间距、屏距做定量分析。它解决的就是从“看懂公式”到“画出对的图”这中间最容易翻车的几段路。2. 干涉条件与光程差为什么条纹是等间距的仿真前要先想清楚这几件事2.1 杨氏双缝的物理模型两个相干点源怎么在屏幕上叠加杨氏双缝干涉的核心是两个相干点源在观察屏上叠加。教材里通常把缝宽忽略不计两个缝 S1 和 S2 之间的距离为 d缝到屏幕的距离为 D屏幕上任意一点 P 到两个缝的光程差 Δ d·sinθ ≈ d·x / D其中 x 是 P 点相对屏幕中心的坐标。这个近似成立的条件是 D 远大于 d 和 x也就是夫琅禾费近似。做仿真时如果缝间距和屏幕尺寸不在这个量级用精确的几何光程差公式算才靠谱Δ sqrt(D² (x - d/2)²) - sqrt(D² (x d/2)²)这个精确公式在仿真初期看着麻烦但你会发现当 D 1 m、d 1 mm 时近似公式和精确公式计算出的条纹间距差异在纳米量级肉眼完全看不出来。可是如果你把缝间距加大到厘米级、屏幕距离缩短到 10 cm 量级近似的误差就会直接导致暗纹位置偏掉好几个像素这是第一个容易踩的坑。2.2 光强公式的MATLAB表达从标量叠加到可视化两个相干点源在屏幕上某一点的光强 I I1 I2 2·sqrt(I1·I2)·cos(δ)δ 是相位差等于 2π·Δ / λ。当两个缝的宽度相同、照明光强相同时I1 I2公式简化为 I 2·I0·(1 cos(δ))或者写成 I 4·I0·cos²(δ/2)。这个简化公式在 MATLAB 里三行就能表达但很多人直接照抄公式后发现一个问题MATLAB 默认的三角函数输入是弧度而有些教材给的 δ 是角度制混用之后条纹周期就会对不上。建议代码里统一先算光程差 Δ单位用米再统一转相位差 δ 2·π·Δ / λ最后用 cos 函数这样不管波长单位是纳米还是微米只要统一换算成米就不会错。2.3 仿真参数的量纲和数量级一份能直接抄的初始化清单参数符号典型值单位说明波长λ632.8e-9m氦氖激光也可设 500e-9缝间距d1e-3m决定条纹密度可调屏幕距离D1.0m决定条纹间距尺度屏幕宽度L0.02m约 ±10 mm 范围采样点数N2001个奇数确保中心点被采样缝宽a0m理想化处理后续可扩展采样点数选奇数这一点特别重要如果选偶数屏幕中心 x 0 这个点可能落在两个采样点之间光强最大值的位置会偏移画出来的条纹中心就不是最亮的这在数值仿真里属于那种“看着对但总差一点”的玄学问题。3. 光强分布计算一维快速验证版本与二维条纹图版本3.1 一维光强分布脚本用 20 行代码先验证公式是否正确先写一个一维版本目的是快速验证干涉公式的相位关系是否正确避免一上来就画二维图出了问题还不好定位。% 杨氏双缝干涉一维光强分布 lambda 632.8e-9; % 波长 632.8 nm单位换算成米 d 1e-3; % 缝间距 1 mm D 1.0; % 缝到屏幕距离 1 m L 0.02; % 屏幕观察范围 ±10 mm N 2001; % 采样点数奇数保证中心点被采样 x linspace(-L/2, L/2, N); % 屏幕坐标向量单位 m % 精确光程差计算 delta sqrt(D^2 (x - d/2).^2) - sqrt(D^2 (x d/2).^2); % 相位差与光强 phi 2 * pi * delta / lambda; % 相位差弧度制 I 4 * cos(phi / 2).^2; % 归一化光强最大值为 4 % 绘图 plot(x * 1000, I, b-, LineWidth, 1.2); xlabel(屏幕位置 x (mm)); ylabel(归一化光强 I/I_0); title(杨氏双缝干涉一维光强分布); grid on;这个脚本用了精确光程差公式而不是近似公式 d·x / D好处是普适性更强。有人说近似公式就够了但当你把缝间距调到 5 mm 再来看两种算法的差异就会变得肉眼可见。用 linspace 生成坐标向量时N 选 2001 这个奇数确保 x 0 这个点是采样点之一否则中心亮纹的光强最大值会被低估看起来整个图案的对称性都会被破坏。相位差计算统一用弧度制避免了角度制混用导致条纹周期变化的低级错误。运行后如果看到一条等间距的余弦平方曲线说明公式部分没问题可以进二维版本了。这段代码跑完你还能顺手验证一个结论把缝间距 d 从 1 mm 改成 2 mm条纹间距会缩小一半。这就是教材里说的“缝间距越大条纹越密”现在你有了亲手验证的工具了。3.2 二维干涉条纹图从一维扩展到平面难点在网格构建二维版本的作用是生成一张能直接放进课程设计报告里的干涉条纹图。和一维版本的区别在于屏幕坐标变成二维平面 (x, y)需要用 meshgrid 生成网格矩阵。% 杨氏双缝干涉二维条纹图 lambda 500e-9; % 波长 500 nm绿色光 d 2e-3; % 缝间距 2 mm D 1.0; % 屏幕距离 L 0.01; % 观察范围 ±5 mm更精细的条纹 N 1000; % 每个维度的采样点数 x linspace(-L/2, L/2, N); y linspace(-L/2, L/2, N); [X, Y] meshgrid(x, y); % 生成二维网格坐标 % 光程差只和 x 坐标有关与 y 无关 delta sqrt(D^2 (X - d/2).^2 Y.^2) - sqrt(D^2 (X d/2).^2 Y.^2); phi 2 * pi * delta / lambda; I 4 * cos(phi / 2).^2; % 干涉项这里为了看清条纹没有加衍射包络 % 显示 figure; imagesc(x * 1000, y * 1000, I); axis xy; % 翻转 y 轴方向让图像和实际对应 colormap(gray); axis equal tight; xlabel(x (mm)); ylabel(y (mm)); title(杨氏双缝干涉条纹无衍射包络); colorbar;这里有个关键的几何问题两个缝沿 y 方向排列还是沿 x 方向排列这决定了条纹的方向。如果缝沿着 y 方向排列光程差只取决于 x 坐标那么干涉条纹应该是一组竖直条纹——同一 x 坐标上的所有 y 值光强相同。代码里 delta 同时依赖 X 和 Y是因为光程差实际上是空间距离差两个缝在 x 方向上有偏移 d/2这样处理物理上更严格。运行后你会看到竖直条纹但注意一个问题这套代码生成的是“纯粹的干涉条纹”没有叠加单缝衍射的包络调制也就是说所有亮纹的亮度都是均匀的。这在做原理验证时可以但如果你的课程设计要求模拟真实的双缝干涉——也就是所谓“缺级”现象那就还得把单缝衍射因子乘进去下一节说这个。3.3 加上单缝衍射包络把“理想条纹”变成“真实干涉图”真实的杨氏双缝实验中每个缝本身有一定的宽度光线经过单缝时会先发生衍射之后两束衍射光再干涉。结果是干涉条纹的亮度被单缝衍射的包络调制某些本该是亮纹的位置变成暗纹这就是“缺级”现象。代码里要做的就是在原来的干涉项上乘以一个单缝衍射因子。% 双缝干涉同时考虑单缝衍射包络 lambda 500e-9; d 2e-3; % 缝间距 a 0.1e-3; % 单缝宽度 0.1 mm D 1.0; L 0.02; N 2001; x linspace(-L/2, L/2, N); [X, Y] meshgrid(x, x); % 干涉项 delta_interf sqrt(D^2 (X - d/2).^2 Y.^2) - sqrt(D^2 (X d/2).^2 Y.^2); phi_interf 2 * pi * delta_interf / lambda; I_interf 4 * cos(phi_interf / 2).^2; % 单缝衍射项用 sinc 函数描述注意 MATLAB 的 sinc 定义是 sin(pi x)/(pi x) beta pi * a * X / (lambda * D); I_diff (sin(beta) ./ beta).^2; I_diff(beta 0) 1; % 处理 beta0 处的除零问题 % 总光强 干涉项 × 衍射包络 I_total I_interf .* I_diff; imagesc(x * 1000, x * 1000, I_total); axis xy; colormap(gray); axis equal tight; xlabel(x (mm)); ylabel(y (mm)); title(双缝干涉含单缝衍射包络);衍射因子用 sin(beta)/betaMATLAB 里可以直接用 sinc 函数但要注意 MATLAB 内置的 sinc(x) 的定义是 sin(πx)/(πx)如果你的 β 已经是 π·a·x/(λD)直接用 sin(beta)./beta 更不容易出错。除零问题要提前处理因为 x 0 时 beta 0对应的是衍射主极大的位置光强应该取极限值 1这在代码里有单独一行处理。跑完这段你能明显看到条纹亮度从中心往两边衰减并且在某些高阶位置出现“该亮却暗”的缺级这就是衍射调制干涉的典型特征。4. 参数扫描与三维可视化波长、缝间距对条纹的定量影响4.1 波长与缝间距的参数扫描条纹间距的线性规律把公式 I ∝ cos²(πdx/(λD)) 和条纹间距公式 Δx λD/d 放在一起看一个直接的推论是条纹间距和波长成正比、和缝间距成反比。用 MATLAB 的 for 循环扫描多组参数把不同参数下的光强分布曲线画在同一张图上这个结论就能看得非常直观。% 参数扫描不同缝间距下的光强分布对比 lambda 632.8e-9; D 1.0; L 0.02; N 2001; x linspace(-L/2, L/2, N); d_list [0.5e-3, 1e-3, 2e-3]; % 三组不同的缝间距 colors [r, g, b]; legend_str {}; figure; hold on; for i 1:length(d_list) d d_list(i); delta sqrt(D^2 (x - d/2).^2) - sqrt(D^2 (x d/2).^2); phi 2 * pi * delta / lambda; I 4 * cos(phi / 2).^2; plot(x * 1000, I, colors(i), LineWidth, 1.2); legend_str{i} sprintf(d %.1f mm, d * 1000); end xlabel(x (mm)); ylabel(归一化光强); legend(legend_str); grid on; hold off;这段代码能帮你目测出条纹周期和缝间距的关系d 从 0.5 mm 加大到 2 mm条纹间距依次减半而且中心亮纹的位置不变。如果要做更严谨的验证可以再加一段循环用 findpeaks 函数提取相邻极大值之间的 x 坐标差再和理论值 λD/d 做对比误差通常在 0.1% 以内——这也能顺便检验你光程差公式算得对不对。4.2 用 surf 和 mesh 展示三维光强分布从截面到立体视角很多课程设计的要求里有一条“绘制三维光强分布图”。实现起来就是在已有二维光强矩阵 I 的基础上用 surf 函数把 x、y、I 三个维度画出来z 轴表示光强。这里有一个容易翻车的点直接用 surf(I) 而不设置坐标轴范围出来的图坐标单位是像素而不是毫米报告里没法看。% 三维光强分布可视化 figure; surf(x * 1000, x * 1000, I_total, EdgeColor, none); xlabel(x (mm)); ylabel(y (mm)); zlabel(光强 I); title(双缝干涉三维光强分布); colormap(jet); view(45, 30); % 设置视角方位角 45 度仰角 30 度 colorbar; % 同时绘制中心截面曲线 figure; I_center I_total(round(N/2), :); % 取中间一行 plot(x * 1000, I_center, b-, LineWidth, 1.5); xlabel(x (mm)); ylabel(光强 I); title(y 0 中心截面光强分布); grid on;view(45, 30) 这个视角参数可以自己调数值只影响观看角度不影响结果。surf 的 EdgeColor 设成 none 是为了去掉网格线否则条纹密集时整个图会被黑色网格线覆盖什么都看不清。中心截面的提取用的是矩阵索引 I_total(round(N/2), :)因为 meshgrid 生成的矩阵行对应 y 方向中间行就是 y 0 的截面这个细节能避免你截错方向拿到一条和预期垂直的曲线。4.3 影响条纹清晰度的隐藏因素屏幕范围与采样密度的匹配很多人仿真出来的条纹图“糊成一片”根本原因不是公式错了而是屏幕范围和采样密度不匹配。屏幕取 ±10 mm、采样点取 500 个那么每个像素代表 0.04 mm 的物理宽度如果条纹间距是 0.316 mm632.8 nm 波长、1 mm 缝距、1 m 屏距每个条纹大约只有 8 个像素勉强能分辨。但如果你把屏幕范围扩到 ±20 mm、采样点还是 500 个每个条纹就只剩下 4 个像素条纹边缘的余弦变化被严重欠采样看起来就是模糊的。这类分辨率问题值得记住。一个实用的经验是先估算条纹间距 Δx λD/d再让屏幕范围和采样点数的比值远小于 Δx也就是每个条纹至少 20 个像素这样出来的条纹边缘才够锐利。代码里可以先算好这个值再选 N避免反复试错。5. 常见问题与避坑指南条纹方向、坐标网格与数据精度5.1 条纹方向转了90度meshgrid 使用不当导致的经典翻车现象明明用的是竖直双缝的模型画出来的条纹却是水平方向的或者图案整体转了个 90 度。原因meshgrid 生成的矩阵中X 矩阵的每一行是相同的每一列才是在 x 方向上变化。如果光程差计算的是 dx 而画图时用的是 Y 矩阵或者反过来条纹方向就会跟着错掉。还有一种情况是 imagesc 默认的 y 轴方向是从上往下如果不加 axis xy图像看起来会是上下颠倒的。解决光程差计算前先想清楚你的缝是沿哪个方向排列的——如果缝沿 y 方向排列典型的竖直双缝光程差应当与 X 矩阵相关反过来就是与 Y 矩阵相关。画图时固定加 axis xy再检查一下图的坐标标注确认 x 轴方向与物理模型的水平方向一致。5.2 暗纹位置总是差半个条纹三角函数参数的单位制问题现象条纹间距看起来对极大值的位置也对但暗纹的位置和理论值偏差半个条纹周期。原因相位差计算时混用了角度制和弧度制。cos 函数默认输入是弧度如果你把 δ 算成了角度制那么 cos(45) 会被当成 cos(45 弧度) 而不是 cos(45 度)结果自然错乱。另一种可能是光程差公式用近似时忽略了符号方向暗纹判断公式里用了 d·x/D 但 x 的零点位置取偏了。解决统一用弧度制全代码里只保留一个转换入口phi 2·pi·delta / lambda。千万不要在某个子表达式里用 deg2rad、某个地方又直接用数值。检查暗纹条件的定义暗纹是相位差为 π 的奇数倍的位置对应 cos² 函数取零的点验证方法是提取 I 矩阵里最小值位置和理论计算值做对比看差值是否小于一个采样间隔。5.3 中心点和两侧的光强不对称采样网格没有关于零点对称现象干涉图案中心左边几个条纹的亮度、宽度和右边不完全一致肉眼能看出来偏了。原因linspace 的起点和终点选取导致网格不对称。比如用 linspace(-0.01, 0.01, 2000)采样点是 2000 个但没有一个点精确落在 x 0 上光强最大值和最小值的位置都会偏离对称中心。严格来说这不是物理问题而是数值采样的瑕疵。解决N 取奇数用 linspace(-L/2, L/2, N) 的时候N 为奇数保证中间点就是 0。习惯上我在几个仿真实战里都会检查 length(x) 的奇偶性如果是偶数直接 N1 或者 N-1 调整。另外diff(x) 要均匀如果 meshgrid 之后发现图像有轻微畸变多半是坐标向量的问题。5.4 三维图一片漆黑或者全是白色数据范围与 colormap 的映射问题现象用 surf 画光强分布时图全黑或全白细节完全看不见。原因surf 会把数据的最小值映射到 colormap 的最下端最大值映射到最上端。如果光强矩阵里有一个特别大的异常值比如某个点的 beta0 没有处理导致除零得到 Inf整个色标都会被拉坏其余数据挤在很窄的范围内看不出层次。解决显示前用 min、max 检查数据范围处理 Inf 和 NaN。用 caxis 限制颜色映射范围比如 caxis([0, 1])把注意力集中在光强的有效范围。确认单缝衍射因子的除零操作是否生效可以输出 I_diff(beta 0) 的值正常应当是 1 而非 NaN。5.5 条纹“糊”成一整片屏幕大小与采样数的平衡点现象仿真出来的干涉条纹边缘模糊亮暗过渡区占据半个条纹宽度看起来不锐利。原因欠采样。一个条纹周期内只有 3-5 个采样点余弦变化被粗暴地线性插值图案当然不锐利。有时候是屏幕范围设置过大导致单位像素对应的物理尺寸远大于条纹间距。解决先用理论公式算条纹间距 Δx λD/d再按“每个条纹至少 20 个采样点”的原则确定 N。比如 Δx 0.5 mm屏幕宽度 20 mm至少需要 40 个条纹周期的分辨能力也就是 N 40 × 20 800实际取 2001 留出余量。采样点过多也有问题计算变慢且图片文件变大N 3000 以上对于二维图会有明显的计算开销自己平衡一下。6. 进阶玩法把仿真代码改成参数化函数一键生成课程设计要用的图表到这一步你已经能用脚本跑出正确的干涉条纹了但脚本的复用性差——每改一个波长就得去翻代码改参数效率太低。进阶做法是把它封装成一个函数输入波长、缝间距、屏距、屏幕范围、采样点数直接输出光强矩阵和几张标准图。这样课程设计里要对比不同参数下的条纹只需一行函数调用就行。function I yds_double_slit(lambda, d, a, D, L, N) % 杨氏双缝干涉仿真函数 % 输入: % lambda - 波长 (m) % d - 缝间距 (m) % a - 单缝宽度 (m)设为 0 表示不考虑衍射包络 % D - 缝到屏幕距离 (m) % L - 屏幕观察范围宽度 (m) % N - 采样点数 (建议奇数) % 输出: % I - 光强矩阵 (N x N) x linspace(-L/2, L/2, N); [X, Y] meshgrid(x, x); delta sqrt(D^2 (X - d/2).^2 Y.^2) - sqrt(D^2 (X d/2).^2 Y.^2); phi 2 * pi * delta / lambda; I 4 * cos(phi / 2).^2; if a 0 beta pi * a * X / (lambda * D); I_diff (sin(beta) ./ beta).^2; I_diff(beta 0) 1; I I .* I_diff; end % 归一化到 0~1 范围方便显示和比较 I I / max(I(:)); % 绘图 figure; imagesc(x * 1000, x * 1000, I); axis xy; colormap(gray); axis equal tight; xlabel(x (mm)); ylabel(y (mm)); title(sprintf(双缝干涉: \\lambda %.0f nm, d %.2f mm, lambda*1e9, d*1e3)); end函数里保留 a 作为单缝宽度的入参传 0 时跳过衍射包络计算这样既能拿来做纯干涉演示也能模拟真实实验的缺级效果。输出做了归一化不同参数下光强对比更公平。调用方式很简单I1 yds_double_slit(632.8e-9, 1e-3, 0, 1.0, 0.02, 2001); I2 yds_double_slit(500e-9, 1e-3, 0, 1.0, 0.02, 2001);这样两行就能对比红光和绿光的条纹密度差异课程设计报告里需要的“不同波长干涉图对比”也就有了素材。更进一步你还能在这个函数基础上加一个 for 循环生成一组不同缝间距的图并保存成 PNG报告插图和答辩 PPT 一次搞定。做这组仿真时我有一个一直坚持的习惯每次改动物理参数后先用理论公式手算一个关键值比如中心亮纹两侧第一个暗纹的位置再去看代码输出对不对得上。这个习惯帮我抓出过好几次坐标单位换算的疏漏——有一次把微米直接当成米用画面里的条纹密到肉眼完全看不出间隔检查理论值才发现差了三个数量级。从那以后我每次仿真都会强制走一遍“手算理论值 → 跑代码 → 对比误差”的流程整套做下来既快又稳。希望帮到你。本文还有配套的精品资源点击获取