简介面向光学、物理及信号处理方向本科与硕士研究生的Matlab仿真资料包聚焦涡旋光波与平面光波干涉的数值模拟适用于毕业设计、课程项目或科研预研内容适合光学、光通信、光镊等领域的初学者快速入门。压缩包共含2个文件包括1个Matlab脚本.m和1张干涉结果示意图.png整体仅61KB轻量易用。脚本基于Matlab 2014/2019a/2021a编写可直接运行并生成干涉图样通过调整轨道角动量模式数等关键参数可对比不同涡旋光束与平面波干涉条纹的差异直观理解涡旋光相位分布、拓扑荷数与干涉图样之间的对应关系。资源内附带运行结果图方便对照校验仿真效果免去调试困扰代码结构简洁便于在此基础上修改参数、继续扩展为更复杂的光场仿真研究也适合作为光学实验教学的辅助素材。已有234人学习下载对需要快速上手光学仿真或完成相关课题设计的同学是一份简洁实用的参考样例同时也有助于深入理解涡旋光的轨道角动量特性及其在光通信、微粒操控等场景中的应用基础。 拿到“基于Matlab模拟涡旋光波和平面光波干涉仿真.zip”这种压缩包大多数人的第一反应是直接解压、运行、看效果。但作为常年和光学仿真打交道的人我建议你多花半小时把里面的物理模型和代码逻辑捋一遍。这个包的核心就是一件事用Matlab分别生成涡旋光波和平面光波的复振幅让它们在空域里叠加再通过强度公式算出干涉图样。项目看起来不大却把光学涡旋、轨道角动量、干涉测量这几个概念很好地串了起来。如果你正在准备光学实验预演、课程设计或者毕业论文里的仿真部分这个项目是很合适的起点。这篇博文就替你把这个zip包拆开讲清楚每一步在算什么以及我实际运行时踩过哪些坑。1. 项目背景与物理原理拆解1.1 涡旋光波和平面光波是什么涡旋光波的核心特征是它的波前呈螺旋状而不是一个平面。描述这个螺旋程度的物理量叫拓扑荷一般用字母 l 表示。可以简单理解成一个波长内相位绕中心轴旋转了多少个2π。l1就是旋转一圈l2就是旋转两圈。旋转导致光束中心出现相位奇异点那里相位没有确定值强度也降为零所以涡旋光中心总有一个暗核。平面光波就好理解得多等相位面是垂直于传播方向的平面相位在横截面上线性变化。两种光叠加之后相位分布产生差异最终在探测器上形成明暗交替的条纹。这些条纹不是随便出现的它正好记录了两束光的相位差在空间上的分布。对于涡旋光来说相位分布里带着 l 的信息所以干涉图样的形状会直接暴露拓扑荷的取值这也是为什么很多实验室用干涉法来检测涡旋光束。1.2 干涉的本质与仿真意义干涉本质上是复振幅相加后取模平方。设涡旋光复振幅为 U_vortex平面光复振幅为 U_plane干涉强度为I(x,y) |U_vortex(x,y) U_plane(x,y)|^2为什么要用复数去表示光因为光不仅自然界有强度还携带着相位。Matlab里的复数变量天然适合做这种运算直接调 abs() 和 .^2 就能得到强度。仿真相对实验最大的优势是参数可调你可以任意修改拓扑荷 l、束腰半径、参考光倾斜角度、波长甚至把 l 设成负数看条纹反向分叉的样子而不需要重新搭光路。这对理解涡旋光、轨道角动量检测、甚至空间光通信里的复用原理都有立竿见影的启发作用。2. Matlab仿真模型搭建与参数设计2.1 网格与光束定义仿真第一步是构建空间坐标系。最常见的做法是用 meshgrid 生成二维网格然后用 cart2pol 把直角坐标换算成极坐标。之所以要极坐标是因为涡旋光的相位分布天然用极角 θ 描述相位因子是 exp(i * l * θ)。注意坐标的单位会影响后续参数设置。如果你严格按物理单位来比如坐标用米波长用532nm那么波数 k 2π/λ 会非常非常大。此时只要参考光稍微偏一个很小的角度条纹密度就可能高到超出采样极限。所以我更推荐在这个仿真里做“归一化坐标”处理把网格范围设定为 [-1,1]然后直接给参考光加一个空间相位倾斜因子比如 exp(i2π(fxx fyy))其中 fx 和 fy 是空间频率控制条纹密度。这样数值范围可控物理含义也清晰。2.2 干涉叠加与强度计算核心计算其实只有三行定义涡旋场、定义平面波参考场、叠加取模平方。但参数设计直接决定最终图画得美不美。涡旋场一般用拉盖尔-高斯模式近似比如U_vortex (r/w)^|l| * exp(-r^2/w^2) * exp(ilθ)这里 w 是光斑束腰半径r 是离中心的距离括号里的 (r/w)^|l| 保证了中心处振幅为零。当 l0 时它自动退化为普通高斯光斑。平面参考场可以加一个线性相位因子来模拟倾斜入射U_plane exp(i2π(fxx fyy))如果你希望参考光和涡旋光完全同轴就把 fx 和 fy 设为 0这样观察到的是螺旋条纹如果稍微偏一点就能看到经典的叉形条纹。两种图样对应着不同的干涉检测模式都值得跑一遍。3. 核心代码解析与实操步骤3.1 参数初始化打开zip包里的主脚本很可能就是一个类似下面的结构。第一步当然是清空环境并设置参数。clear; clc; close all; % 网格参数 N 1024; % 网格数建议至少512 L 2; % 网格总宽度归一化单位 x linspace(-L/2, L/2, N); [X, Y] meshgrid(x, x); [Theta, R] cart2pol(X, Y); % 涡旋光参数 l 1; % 拓扑荷 w 0.8; % 高斯束腰半径 % 参考平面波参数 fx 0.8; % x方向空间频率控制条纹密度 fy 0; % y方向空间频率这里的 N、l、fx、fy 就是你要反复调节的核心变量。N 太低条纹会变形太高计算和绘图变慢l 决定干涉条纹的分叉数fx 决定条纹密度。我一般先把 N 设成512跑通确认没问题再上1024出高清图。3.2 核心仿真代码这里直接给出核心计算片段这是整个项目最精华的部分。% 涡旋光复振幅使用拉盖尔-高斯近似 amplitude (R / w) .^ abs(l) .* exp(-R.^2 / w^2); vortex amplitude .* exp(1i * l * Theta); % 平面参考光复振幅 plane exp(1i * 2 * pi * (fx .* X fy .* Y)); % 干涉强度 I abs(vortex plane).^2;有没有注意到 vortex 中心位置amplitude 计算时 R/w 在 R0 时会变成0的绝对值次方当 l0 时会有 0^0 的情况Matlab会把这个当作1所以 l0 时不会出错。但当 l0 时R/w 在中心为0振幅为0这是涡旋光暗核的正确体现也会规避相位奇点处的数值问题。3.3 绘制干涉图样绘图这一步决定了你最后呈现出来的效果。建议至少画三个子图强度图、相位图、中心一维剖面。figure(Position, [100 100 1400 400]); % 灰度强度图 subplot(1,3,1); imagesc(x, x, I); axis image; colormap(gray); xlabel(x); ylabel(y); title([干涉强度图, l , num2str(l)]); % 相位图 subplot(1,3,2); imagesc(x, x, angle(vortex plane)); axis image; colormap(hsv); xlabel(x); ylabel(y); title(合成波前相位); % 中心横向剖面 subplot(1,3,3); plot(x, I(:, round(N/2)), LineWidth, 1.5); xlabel(x); ylabel(强度); title(y0 强度剖面); grid on;相位图用 hsv 色图可以清楚看到螺旋相位线的位置以及平面波带来的倾斜相位如何与螺旋相位叠加。剖面图则帮你确认暗核是否存在——在中心位置强度应该有个明显的凹陷这是涡旋光干涉图样的标志性特征。4. 仿真结果分析与物理验证4.1 拓扑荷l对干涉条纹的决定性影响把 l 从0改到1、2、3重新运行你会看到干涉条纹分叉数量与 l 严格相等。l0 时涡旋光退化成了普通高斯光干涉图样就是一组竖直的平行条纹和两束平面波干涉没有任何区别。l1 时条纹中央出现一个叉口上下两部分条纹错开一个周期。l2 时出现两个叉口l3 时是三个。这个特征在光学检测里经常被用来反向推算未知涡旋光束的拓扑荷。我把不同 l 下的特征整理成了一张表方便你对照检查自己的仿真结果。拓扑荷 l同轴干涉图样离轴干涉图样0无涡旋均匀光斑/平行条纹等间距平行直条纹1单臂阿基米德螺旋条纹单叉形条纹叉口指向由 l 正负决定2双臂螺旋条纹双叉形条纹像两把叉子并排-1反向螺旋叉口方向与 l1 相反注意拓扑荷的正负在强度图里往往不容易直接分辨因为强度是模平方运算会丢掉绝对相位符号。但如果你把干涉条纹叉口的方向和参考光倾斜方向结合起来看就能区分。仿真时可以改变 l-1 再运行一次对比两张图你会看到叉口的弯曲方向反过来了。4.2 参考光倾斜角度与条纹拓扑结构参考光的 fx、fy 直接影响条纹周期和走向。fx 越大条纹越密fx 为正条纹沿 x 方向倾斜。当你把 fx 设为0、fy 也设为0时干涉图样会从叉形条纹变成连续螺旋条纹。这个转变很有意思因为同轴干涉时某一点的相位差就是涡旋相位 lθ等相位差的轨迹满足 θ 常数 m2π/l也就是辐射状的螺旋线所以在极坐标下看条纹是一圈圈从中心伸出来的螺旋臂。实际操作中你会发现 fx 设置过大时条纹密度可能超过奈奎斯特采样率出现摩尔纹需要加密网格或者减小 fx。反过来fx 太小也不行叉形区域会变成几乎贴着中心的短线条无法清晰显示分叉。我个人的经验是在归一化坐标 [-1,1] 下fx 取 0.5 到 1.5 这个范围图样最舒适。5. 常见问题与调试心得5.1 相位奇点处NaN与暗核噪声仿真中最容易遇到的问题是中心区域出现奇怪的NaN或不连续花纹。原因通常出在极角 θ 的计算和振幅归一化上。e^ (ilθ) 本身在中心是良定义的但如果你把振幅直接写成了 (R.^abs(l))当 R0 且 l0 时振幅为0没有毛病但如果你用 R.^abs(l) 去除其他量R0 就会出现 0/0。所以建议在振幅计算时始终加上一个极小量比如 R/(weps)或者直接用 (R.^abs(l)) 乘以衰减项避免产生NaN。另一个暗核噪声的来源是灰度图显示范围。涡旋光中心强度很低周围强度很高直接 imagesc 会把暗部压没。可以调整显示范围比如 imagesc(x, x, I, [0 max(I(:))*0.8]); 这样暗核的边缘会更清晰。5.2 条纹走向与meshgrid/ndgrid的坑Matlab里 meshgrid 产生矩阵时X 的每一行是相同的Y 的每一列是相同的而 ndgrid 则刚好相反。如果你原本代码用的是 meshgrid后来为了性能改用 ndgrid会导致 X 和 Y 互换干涉条纹的方向会旋转90度拓扑荷的符号也会反转。这个问题很隐蔽排查起来特别费时间。我的习惯是统一用 meshgrid并且开头注释里写清楚避免自己隔几天再回来忘了。另外imagesc 的坐标轴方向可能导致条纹看起来镜像翻转。如果希望 y 轴正方向朝上记得加一句 axis xy; 否则默认的 axis image 会让 y 轴反向螺旋条纹看起来像反了容易被误判为 l 符号出问题。5.3 参数微调与程序包使用建议这个zip包里的脚本通常可以直接运行但我不建议只满足于得到一张默认图。你可以做这样几个实验保持参考光不变把 l 设成 -3 到 3保存所有结果做一个对比图保持 l1把 fx 从0逐渐增大到2观察条纹密度变化把平面参考光也换成不同半径的高斯光看看条纹对比度如何随两束光的强度比变化。改成参数扫描时可以套一个 for 循环用 sgtitle 给每个子图标上 l 的参数值生成一张网格图。这几乎是很多光学组汇报PPT里的标准插图方式操作起来也没难度figure(Position, [50 50 900 600]); for l -3:3 % ...重新计算涡旋场和干涉强度... subplot(2, 4, l 4); % 中间留个空位 imagesc(x, x, I); axis image; colormap(gray); title([l , num2str(l)]); end最后分享一个小技巧仿真完成后用 exportgraphics(gcf, fork_fringe_l1.png, Resolution, 300) 导出图片而不是直接用截图。高分辨率图放到论文或者报告里几乎没有锯齿清晰度完全够用。这个操作在最近的Matlab版本里都很稳定比 print 命令更直观。把这些边角功夫做好了整个仿真项目才算真正落到了可展示、可汇报的层面。本文还有配套的精品资源点击获取
Matlab涡旋光与平面光干涉仿真：原理、代码与调试指南
简介面向光学、物理及信号处理方向本科与硕士研究生的Matlab仿真资料包聚焦涡旋光波与平面光波干涉的数值模拟适用于毕业设计、课程项目或科研预研内容适合光学、光通信、光镊等领域的初学者快速入门。压缩包共含2个文件包括1个Matlab脚本.m和1张干涉结果示意图.png整体仅61KB轻量易用。脚本基于Matlab 2014/2019a/2021a编写可直接运行并生成干涉图样通过调整轨道角动量模式数等关键参数可对比不同涡旋光束与平面波干涉条纹的差异直观理解涡旋光相位分布、拓扑荷数与干涉图样之间的对应关系。资源内附带运行结果图方便对照校验仿真效果免去调试困扰代码结构简洁便于在此基础上修改参数、继续扩展为更复杂的光场仿真研究也适合作为光学实验教学的辅助素材。已有234人学习下载对需要快速上手光学仿真或完成相关课题设计的同学是一份简洁实用的参考样例同时也有助于深入理解涡旋光的轨道角动量特性及其在光通信、微粒操控等场景中的应用基础。 拿到“基于Matlab模拟涡旋光波和平面光波干涉仿真.zip”这种压缩包大多数人的第一反应是直接解压、运行、看效果。但作为常年和光学仿真打交道的人我建议你多花半小时把里面的物理模型和代码逻辑捋一遍。这个包的核心就是一件事用Matlab分别生成涡旋光波和平面光波的复振幅让它们在空域里叠加再通过强度公式算出干涉图样。项目看起来不大却把光学涡旋、轨道角动量、干涉测量这几个概念很好地串了起来。如果你正在准备光学实验预演、课程设计或者毕业论文里的仿真部分这个项目是很合适的起点。这篇博文就替你把这个zip包拆开讲清楚每一步在算什么以及我实际运行时踩过哪些坑。1. 项目背景与物理原理拆解1.1 涡旋光波和平面光波是什么涡旋光波的核心特征是它的波前呈螺旋状而不是一个平面。描述这个螺旋程度的物理量叫拓扑荷一般用字母 l 表示。可以简单理解成一个波长内相位绕中心轴旋转了多少个2π。l1就是旋转一圈l2就是旋转两圈。旋转导致光束中心出现相位奇异点那里相位没有确定值强度也降为零所以涡旋光中心总有一个暗核。平面光波就好理解得多等相位面是垂直于传播方向的平面相位在横截面上线性变化。两种光叠加之后相位分布产生差异最终在探测器上形成明暗交替的条纹。这些条纹不是随便出现的它正好记录了两束光的相位差在空间上的分布。对于涡旋光来说相位分布里带着 l 的信息所以干涉图样的形状会直接暴露拓扑荷的取值这也是为什么很多实验室用干涉法来检测涡旋光束。1.2 干涉的本质与仿真意义干涉本质上是复振幅相加后取模平方。设涡旋光复振幅为 U_vortex平面光复振幅为 U_plane干涉强度为I(x,y) |U_vortex(x,y) U_plane(x,y)|^2为什么要用复数去表示光因为光不仅自然界有强度还携带着相位。Matlab里的复数变量天然适合做这种运算直接调 abs() 和 .^2 就能得到强度。仿真相对实验最大的优势是参数可调你可以任意修改拓扑荷 l、束腰半径、参考光倾斜角度、波长甚至把 l 设成负数看条纹反向分叉的样子而不需要重新搭光路。这对理解涡旋光、轨道角动量检测、甚至空间光通信里的复用原理都有立竿见影的启发作用。2. Matlab仿真模型搭建与参数设计2.1 网格与光束定义仿真第一步是构建空间坐标系。最常见的做法是用 meshgrid 生成二维网格然后用 cart2pol 把直角坐标换算成极坐标。之所以要极坐标是因为涡旋光的相位分布天然用极角 θ 描述相位因子是 exp(i * l * θ)。注意坐标的单位会影响后续参数设置。如果你严格按物理单位来比如坐标用米波长用532nm那么波数 k 2π/λ 会非常非常大。此时只要参考光稍微偏一个很小的角度条纹密度就可能高到超出采样极限。所以我更推荐在这个仿真里做“归一化坐标”处理把网格范围设定为 [-1,1]然后直接给参考光加一个空间相位倾斜因子比如 exp(i2π(fxx fyy))其中 fx 和 fy 是空间频率控制条纹密度。这样数值范围可控物理含义也清晰。2.2 干涉叠加与强度计算核心计算其实只有三行定义涡旋场、定义平面波参考场、叠加取模平方。但参数设计直接决定最终图画得美不美。涡旋场一般用拉盖尔-高斯模式近似比如U_vortex (r/w)^|l| * exp(-r^2/w^2) * exp(ilθ)这里 w 是光斑束腰半径r 是离中心的距离括号里的 (r/w)^|l| 保证了中心处振幅为零。当 l0 时它自动退化为普通高斯光斑。平面参考场可以加一个线性相位因子来模拟倾斜入射U_plane exp(i2π(fxx fyy))如果你希望参考光和涡旋光完全同轴就把 fx 和 fy 设为 0这样观察到的是螺旋条纹如果稍微偏一点就能看到经典的叉形条纹。两种图样对应着不同的干涉检测模式都值得跑一遍。3. 核心代码解析与实操步骤3.1 参数初始化打开zip包里的主脚本很可能就是一个类似下面的结构。第一步当然是清空环境并设置参数。clear; clc; close all; % 网格参数 N 1024; % 网格数建议至少512 L 2; % 网格总宽度归一化单位 x linspace(-L/2, L/2, N); [X, Y] meshgrid(x, x); [Theta, R] cart2pol(X, Y); % 涡旋光参数 l 1; % 拓扑荷 w 0.8; % 高斯束腰半径 % 参考平面波参数 fx 0.8; % x方向空间频率控制条纹密度 fy 0; % y方向空间频率这里的 N、l、fx、fy 就是你要反复调节的核心变量。N 太低条纹会变形太高计算和绘图变慢l 决定干涉条纹的分叉数fx 决定条纹密度。我一般先把 N 设成512跑通确认没问题再上1024出高清图。3.2 核心仿真代码这里直接给出核心计算片段这是整个项目最精华的部分。% 涡旋光复振幅使用拉盖尔-高斯近似 amplitude (R / w) .^ abs(l) .* exp(-R.^2 / w^2); vortex amplitude .* exp(1i * l * Theta); % 平面参考光复振幅 plane exp(1i * 2 * pi * (fx .* X fy .* Y)); % 干涉强度 I abs(vortex plane).^2;有没有注意到 vortex 中心位置amplitude 计算时 R/w 在 R0 时会变成0的绝对值次方当 l0 时会有 0^0 的情况Matlab会把这个当作1所以 l0 时不会出错。但当 l0 时R/w 在中心为0振幅为0这是涡旋光暗核的正确体现也会规避相位奇点处的数值问题。3.3 绘制干涉图样绘图这一步决定了你最后呈现出来的效果。建议至少画三个子图强度图、相位图、中心一维剖面。figure(Position, [100 100 1400 400]); % 灰度强度图 subplot(1,3,1); imagesc(x, x, I); axis image; colormap(gray); xlabel(x); ylabel(y); title([干涉强度图, l , num2str(l)]); % 相位图 subplot(1,3,2); imagesc(x, x, angle(vortex plane)); axis image; colormap(hsv); xlabel(x); ylabel(y); title(合成波前相位); % 中心横向剖面 subplot(1,3,3); plot(x, I(:, round(N/2)), LineWidth, 1.5); xlabel(x); ylabel(强度); title(y0 强度剖面); grid on;相位图用 hsv 色图可以清楚看到螺旋相位线的位置以及平面波带来的倾斜相位如何与螺旋相位叠加。剖面图则帮你确认暗核是否存在——在中心位置强度应该有个明显的凹陷这是涡旋光干涉图样的标志性特征。4. 仿真结果分析与物理验证4.1 拓扑荷l对干涉条纹的决定性影响把 l 从0改到1、2、3重新运行你会看到干涉条纹分叉数量与 l 严格相等。l0 时涡旋光退化成了普通高斯光干涉图样就是一组竖直的平行条纹和两束平面波干涉没有任何区别。l1 时条纹中央出现一个叉口上下两部分条纹错开一个周期。l2 时出现两个叉口l3 时是三个。这个特征在光学检测里经常被用来反向推算未知涡旋光束的拓扑荷。我把不同 l 下的特征整理成了一张表方便你对照检查自己的仿真结果。拓扑荷 l同轴干涉图样离轴干涉图样0无涡旋均匀光斑/平行条纹等间距平行直条纹1单臂阿基米德螺旋条纹单叉形条纹叉口指向由 l 正负决定2双臂螺旋条纹双叉形条纹像两把叉子并排-1反向螺旋叉口方向与 l1 相反注意拓扑荷的正负在强度图里往往不容易直接分辨因为强度是模平方运算会丢掉绝对相位符号。但如果你把干涉条纹叉口的方向和参考光倾斜方向结合起来看就能区分。仿真时可以改变 l-1 再运行一次对比两张图你会看到叉口的弯曲方向反过来了。4.2 参考光倾斜角度与条纹拓扑结构参考光的 fx、fy 直接影响条纹周期和走向。fx 越大条纹越密fx 为正条纹沿 x 方向倾斜。当你把 fx 设为0、fy 也设为0时干涉图样会从叉形条纹变成连续螺旋条纹。这个转变很有意思因为同轴干涉时某一点的相位差就是涡旋相位 lθ等相位差的轨迹满足 θ 常数 m2π/l也就是辐射状的螺旋线所以在极坐标下看条纹是一圈圈从中心伸出来的螺旋臂。实际操作中你会发现 fx 设置过大时条纹密度可能超过奈奎斯特采样率出现摩尔纹需要加密网格或者减小 fx。反过来fx 太小也不行叉形区域会变成几乎贴着中心的短线条无法清晰显示分叉。我个人的经验是在归一化坐标 [-1,1] 下fx 取 0.5 到 1.5 这个范围图样最舒适。5. 常见问题与调试心得5.1 相位奇点处NaN与暗核噪声仿真中最容易遇到的问题是中心区域出现奇怪的NaN或不连续花纹。原因通常出在极角 θ 的计算和振幅归一化上。e^ (ilθ) 本身在中心是良定义的但如果你把振幅直接写成了 (R.^abs(l))当 R0 且 l0 时振幅为0没有毛病但如果你用 R.^abs(l) 去除其他量R0 就会出现 0/0。所以建议在振幅计算时始终加上一个极小量比如 R/(weps)或者直接用 (R.^abs(l)) 乘以衰减项避免产生NaN。另一个暗核噪声的来源是灰度图显示范围。涡旋光中心强度很低周围强度很高直接 imagesc 会把暗部压没。可以调整显示范围比如 imagesc(x, x, I, [0 max(I(:))*0.8]); 这样暗核的边缘会更清晰。5.2 条纹走向与meshgrid/ndgrid的坑Matlab里 meshgrid 产生矩阵时X 的每一行是相同的Y 的每一列是相同的而 ndgrid 则刚好相反。如果你原本代码用的是 meshgrid后来为了性能改用 ndgrid会导致 X 和 Y 互换干涉条纹的方向会旋转90度拓扑荷的符号也会反转。这个问题很隐蔽排查起来特别费时间。我的习惯是统一用 meshgrid并且开头注释里写清楚避免自己隔几天再回来忘了。另外imagesc 的坐标轴方向可能导致条纹看起来镜像翻转。如果希望 y 轴正方向朝上记得加一句 axis xy; 否则默认的 axis image 会让 y 轴反向螺旋条纹看起来像反了容易被误判为 l 符号出问题。5.3 参数微调与程序包使用建议这个zip包里的脚本通常可以直接运行但我不建议只满足于得到一张默认图。你可以做这样几个实验保持参考光不变把 l 设成 -3 到 3保存所有结果做一个对比图保持 l1把 fx 从0逐渐增大到2观察条纹密度变化把平面参考光也换成不同半径的高斯光看看条纹对比度如何随两束光的强度比变化。改成参数扫描时可以套一个 for 循环用 sgtitle 给每个子图标上 l 的参数值生成一张网格图。这几乎是很多光学组汇报PPT里的标准插图方式操作起来也没难度figure(Position, [50 50 900 600]); for l -3:3 % ...重新计算涡旋场和干涉强度... subplot(2, 4, l 4); % 中间留个空位 imagesc(x, x, I); axis image; colormap(gray); title([l , num2str(l)]); end最后分享一个小技巧仿真完成后用 exportgraphics(gcf, fork_fringe_l1.png, Resolution, 300) 导出图片而不是直接用截图。高分辨率图放到论文或者报告里几乎没有锯齿清晰度完全够用。这个操作在最近的Matlab版本里都很稳定比 print 命令更直观。把这些边角功夫做好了整个仿真项目才算真正落到了可展示、可汇报的层面。本文还有配套的精品资源点击获取
