在实验室里搭光路的同行应该都有过这个经历激光器出来的光斑天生是个高斯分布中心亮、边缘暗可好多应用偏偏需要一块照度均匀的平顶光斑。做激光微加工、做光镊阵列、做荧光显微照明谁没为这个发过愁。传统办法是拿透镜组加光阑硬切能量损失大不说边缘还带着一圈圈衍射条纹根本谈不上“平顶”。换成振幅掩膜倒是能整形但光一过掩膜就损耗大半效率实在难看。后来空间光调制器SLM普及了问题才算有了优雅的解法。SLM只调相位不调振幅理论上能量可以几乎无损地通过所有整形任务都落在“如何算出一幅正确的相位图”上。这篇文章就聊聊我用MATLAB实现高斯光束到平顶光束转变的完整过程核心两条路都走一遍一条是大家熟知的GS算法Gerchberg-Saxton迭代逼近最优相位另一条是直接从能量守恒出发解析计算SLM相位分布。两条路没有绝对好坏关键看你的目标光斑是什么形状、对实时性要求有多高。下面我用实际跑过的代码和踩过的坑把整个流程掰开讲清楚。1. 项目全貌从需求到方案的取舍1.1 先搞清楚物理问题再动手先说物理模型。入射的高斯光束在SLM面上的复振幅可以写成E_in(r) E0 * exp(-r^2 / w0^2) * exp(i * φ(r))其中w0是束腰半径φ(r)就是要加载到SLM上的相位分布。光束经过SLM反射后再经过一个透镜在透镜焦面上形成输出光场。透镜的作用就是一个光学傅里叶变换焦面上的复振幅正比于SLM面上复振幅的二维傅里叶变换。我们的目标是让焦面上的光强分布变成一个“帽子”形状中心区域平坦均匀边缘快速下降。常见的两种描述方式一是理想平顶就是一个圆形区域内的均匀光强二是超高斯分布用 I(ρ) I0 * exp(-2(ρ/R)^n) 来近似n越大边缘越陡。为什么非要SLM呢因为单纯用衰减掩膜整形原理上就是把高斯光斑中间强的地方挡掉一部分让整体强度变均匀。但光被挡掉就意味着能量浪费而且掩膜上的硬边会产生严重衍射。相位整形则是完全不同的思路我不减少任何光能量只是改变不同半径处光线的传播方向把中心多余的能量“搬运”到边缘去。这就像把一堆沙子从中间推到四周总量一分没少只是摊平了。1.2 GS算法和直接计算法怎么选GS算法的本质是一个迭代优化过程。它不知道解析解但可以反复在“输入面”和“输出面”之间来回做傅里叶变换每一次都把振幅替换成我们想要的目标值同时保留相位信息。迭代几十次之后相位分布就会收敛到一个让输出光强逼近目标的结果。直接计算法则完全不用迭代。它从能量守恒出发推导出入射高斯光束中每一个半径位置的能量应该被映射到输出面上的哪个半径位置再由几何光学的偏折关系反推SLM上每一点的相位梯度最后积分得到完整的相位分布。这俩怎么选我做了个对比表实际项目里我基本按这个逻辑判断对比维度GS算法直接计算法目标光斑形状任意形状圆形、方形、图案都行基本只适合圆对称目标计算速度迭代几十上百次一般1秒内也能完成一次积分搞定毫秒级实现难度代码逻辑简单但调参有讲究数学推导绕人代码反而简单抗噪声能力迭代过程自带一定的鲁棒性对光束质量和模型偏差比较敏感可扩展性可以加约束、加权重玩法多改个目标形状基本就要重新推公式一句话总结如果你只是想要一个圆形平顶光斑直接计算法又快又稳如果你想要的平顶是方形、环形或者要在平顶基础上叠加复杂图案那就老老实实上GS。我在实际项目里经常两种都跑一遍GS的结果作为精度基准直接计算法用来做实时调参两者配合效率最高。2. 先啃硬骨头GS算法的MATLAB实现全解2.1 把GS算法的“对答案”逻辑吃透GS算法的直觉其实特别朴素。想象你和一个朋友分别在两个房间里中间隔着一堵墙墙上只留一个小孔让你们传递纸条。你手里有一张目标照片输出面的振幅分布但缺一个相位信息你朋友手里有一张已知的照片输入面的振幅分布也缺相位。你们的任务是隔着一堵墙共同猜出一对相位让最终合成的照片符合目标。GS算法的做法就是反复“对答案”从你朋友那边拿到一版带猜测相位的照片通过“傅里叶变换”这个传纸条通道送到你手里你把你手里的目标振幅换上保留传过来的相位再传回去你朋友收到之后把他那边的已知振幅换上保留你传过去的相位再传回来。就这么来回倒腾每一次对答案都在修正相位最后两边同时满意收敛完成。对应到光学系统里“输入面”就是SLM所在的平面“输出面”就是透镜焦面。迭代过程中的两个约束条件非常明确输入面振幅已知就是高斯光束的振幅分布或者说SLM反射后的振幅分布输出面振幅已知就是我们要的平顶光斑的振幅分布。唯一自由的就是两个面上的相位GS算法就是把这个自由度榨干。2.2 MATLAB代码骨架与参数初始化我习惯用512×512的网格来仿真物理尺寸设置成和实际SLM一致。这里给一个我常用的参数配置直接抄作业完全没问题% 基本物理参数 lambda 532e-9; % 波长532nm绿光 N 512; % 网格点数 pixel_size 8e-6; % SLM像素尺寸我这里模拟8um像素的SLM L N * pixel_size; % SLM面物理尺寸 f 300e-3; % 透镜焦距300mm % 构建坐标网格 x linspace(-L/2, L/2, N); [X, Y] meshgrid(x, x); R sqrt(X.^2 Y.^2); % 极坐标半径 % 高斯光束参数 w0 2e-3; % 束腰半径2mm E_in_amp exp(-(R/w0).^2); % 高斯振幅分布 % 目标平顶参数 R_flat 1.5e-3; % 目标平顶半径1.5mm target_amp zeros(N); target_amp(R R_flat) 1; % 理想平顶振幅 % 初始相位 phase zeros(N); % 可以试试随机相位效果不同 U0_in E_in_amp .* exp(1i * phase);GS迭代的核心循环我用的是这种写法iter 100; for k 1:iter % 正向传播从SLM面到焦面傅里叶变换 U_focal fftshift(fft2(ifftshift(U0_in))); % 焦面替换振幅保留相位振幅换成目标 U_focal target_amp .* exp(1i * angle(U_focal)); % 反向传播从焦面回到SLM面逆傅里叶变换 U_slm ifftshift(ifft2(fftshift(U_focal))); % SLM面替换振幅保留相位振幅回到高斯 U0_in E_in_amp .* exp(1i * angle(U_slm)); % 记录收敛曲线方便最后分析 err(k) 1 - abs(sum(sum(conj(U_focal_target) .* U_focal)))^2 ... / (sum(sum(abs(U_focal_target).^2)) * sum(sum(abs(U_focal).^2))); end phase_gs angle(U0_in); % 最终相位这里有几个细节必须强调fftshift和ifftshift的位置不能搞反。如果你用meshgrid生成的坐标是中心为零的对称网格那么输入矩阵的正中心位置对应的就是零频分量傅里叶变换之后要再用fftshift把零频挪回中心。我的习惯是先建好对称坐标正向变换前做ifftshift把零频挪到矩阵左上角fft2之后用fftshift挪回来反向传播对称操作。另一个容易忽略的点是MATLAB的fft2不带物理单位变换结果和网格点数、物理尺寸都存在比例关系。在实际仿真中我们通常只关心相对强度分布所以这个比例系数可以直接省掉。但如果你要把仿真结果和实验绝对强度对比就需要按傅里叶光学的关系补上系数这一步是很多人对不上实验结果的原因。2.3 迭代参数调节心得别让算法白跑几十轮先说迭代次数。我实测下来50次迭代就能看到比较像样的平顶轮廓100次左右基本就稳定了200次以上几乎没有肉眼可见的变化。GS算法的收敛曲线不是直线下降前20次下降最快后面就是在反复修正细节所以没必要盲目堆迭代次数。我习惯设置一个收敛判据比如相邻两次迭代的误差变化小于千分之一就提前退出省时间。再说初始相位。很多人第一次跑GS直接用零相位初始化最后结果也凑合。但实际上初始相位的选择会影响收敛速度和最终质量。我试过三种初始化方式差异不小零相位简单但收敛速度偏慢而且容易陷入局部最优解平顶边缘会出现一些不规则的毛刺随机相位收敛快平顶质量好但每次跑出来的相位都不一样实验复现性差二次相位离焦项相位初始化成圆锥形相当于提前给光束一个发散角收敛速度和随机相位接近而且结果稳定。实际项目中我推荐第三种初始化加一个轻微的二次相位比如phase_init 5e6 * R.^2系数的量级需要根据光斑尺寸试兼顾了收敛速度和复现性。还有一个大坑是目标平顶半径的选取。GS算法不是万能的它只能在傅里叶平面上有效范围内“变魔术”。如果目标平顶半径太大接近甚至超过焦面上有效衍射极限范围算法会拼命振荡牺牲平顶质量换取能量守恒最后出来的光斑边缘一大堆花纹。我的经验值是目标平顶半径不超过SLM有效区域对应衍射极限半径的60%~70%。举个例子对于8um像素、512×512的SLM532nm波长、300mm焦距焦面上可用的最大范围大约在3mm量级那目标平顶半径取1~2mm是合理的超过2.5mm就会出现明显的质量劣化。最后聊一下目标图形的边缘。直接用圆形硬边矩形窗函数构造目标振幅平顶边缘会有明显的振铃现象就是一圈圈明暗交替的纹路。这不是GS算法的bug而是矩形窗函数在频域本身就带有旁瓣。解决办法是把目标边缘换成平滑过渡比如用余弦渐变窗或者超高斯分布让目标强度在边缘连续下降而不是突变。我一般这样处理% 超高斯目标n越大边缘越陡 n 8; target_amp exp(-(R/R_flat*2).^n); target_amp(target_amp 0.01) 0; % 截断太小的值避免无效迭代用了超高斯目标之后仿真的平顶边缘虽然有一点点斜但振铃完全消失而且真实实验里边缘稍微倾斜一点毫无影响。这个技巧直接改善最终光斑质量。3. 换个思路直接算相位分布的推导与实现3.1 从能量守恒到坐标映射先把数学搞顺GS算法虽然好用但终究是迭代的而且每次迭代结果都有微小的随机性。有些场景需要快速得到相位图或者需要精确可重复的结果那就得靠直接计算法。这个方法的起点是能量守恒。设在SLM面上径向位置r处的一个圆环内包含的高斯光束功率为P_in(r) P0 * (1 - exp(-2 * r^2 / w0^2))这其实就是高斯分布从中心到r的累积功率。这个功率经过整形之后应该全部落在输出面焦面上对应半径ρ的圆内。如果目标平顶光束的强度是均匀的I_flat在半径ρ的圆内功率为P_out(ρ) I_flat * π * ρ^2令两者相等取边界条件当r趋于无穷或者说r r_maxSLM有效半径时ρ R_flat也就是目标平顶半径。就能得到输入半径r和输出半径ρ的一一映射ρ(r) R_flat * sqrt(1 - exp(-2 * r^2 / w0^2))这个公式的含义很直观高斯光束中心区域能量密度高对应到输出面上应该被“摊”到比较大的半径范围而高斯光束边缘能量密度低对应输出面上的半径增量就小。你可以把它理解成用一张弹性膜去铺一床被子中心部分拉伸得厉害一些边缘部分几乎不动。如果你用这个映射去统计输出强度会发现每一段输入半径对应的能量经过拉伸或压缩之后强度恰好均匀——这正是能量守恒保证的。这个关系式是整个直接计算法的灵魂后面所有步骤都是它的延续。3.2 从相位梯度到完整相位图有了坐标映射剩下的问题就是怎么让光线按照这个映射从半径r的位置偏转到半径ρ的位置。这要用到几何光学里相位梯度和光线偏折角的关系。在傍轴近似下SLM上某一点处相位梯度 dφ/dr 决定了光线出射后横向偏折的角度θ近似关系是θ ≈ (1/k) * dφ/dr其中k 2π/λ是波数。光束经过焦距为f的透镜后在焦面上横向偏移量ρ和偏折角的关系为ρ ≈ f * θ把两个式子合起来就得到SLM相位梯度和目标映射的关系dφ/dr k * ρ(r) / f (2π/(λf)) * R_flat * sqrt(1 - exp(-2 * r^2 / w0^2))最后一步就是对这个梯度做径向积分得到完整的相位分布φ(r) ∫0^r (2π/(λf)) * R_flat * sqrt(1 - exp(-2 * r^2 / w0^2)) dr这个积分没有初等函数的解析表达但不妨碍我们直接用数值积分算。整个过程在MATLAB里只需要几行% 径向坐标和相位梯度的数值计算 r linspace(0, L/2, N/21); % 只算径向利用圆对称性 rho R_flat * sqrt(1 - exp(-2 * r.^2 / w0^2)); % 相位梯度 dphi_dr (2*pi / (lambda * f)) * rho; % 数值积分得到相位分布 phi_radial cumtrapz(r, dphi_dr); % 扩展到二维 [XX, YY] meshgrid(x, x); RR sqrt(XX.^2 YY.^2); phi2d interp1(r, phi_radial, RR, linear, 0);一个很容易踩的坑是积分区间的处理。r的取值必须覆盖SLM有效面上所有像素对应的径向位置而且要保证最远处的相位梯度不为零。如果你在高斯光束的尾部直接截断截断位置对应的ρ值小于R_flat那平顶的实际半径就会比设计值小。我建议r取到SLM半宽末尾超出w0两三倍以上让映射关系尽量覆盖完整的高斯分布。3.3 直接计算法的适用边界哪里能用哪里别硬用直接计算法推导过程漂亮、计算速度快但它有三个前提偏离了结果就崩第一个是圆对称性。整套推导都建立在高斯光束是圆对称、目标平顶也是圆对称的基础上。如果你要的是方形平顶、六角形平顶或者中心带孔的环形平顶这套公式直接失效只能回到GS算法。第二个是傍轴近似。dφ/dr和偏折角的线性关系只在光线偏折角很小的时候成立。如果设计的目标平顶半径很大导致SLM边缘的相位梯度非常大光线偏折角超出傍轴范围实际光斑就会和仿真结果明显偏移。我一般控制最大相位梯度对应的偏折角不超过10度左右。第三个是模型的理想性。推导假设高斯光束是完美的TEM00模、SLM相位调制是线性的、透镜没有像差。真实实验里这些都不可能完美所以直接计算法的结果只能作为“初始相位”实际光路里往往还要微调。即便如此直接计算法的价值依然很大。它零迭代、毫秒级出结果非常适合做批量扫描参数比如你想看平顶半径从1mm变化到2mm时输出光斑长什么样GS算法要跑一百次迭代直接计算法瞬间就能出全部结果。这种速度优势在实验现场调试时非常宝贵。4. 仿真和实验中间隔了好几个坑4.1 平顶光束质量怎么量化评价很多人仿真跑完感觉光斑“看着挺平”但到底平不平、好到什么程度心里没数。我给一个我自用的评价流程全是数值指标不用拍脑袋均匀度flatness用平顶区域内光强的相对标准差越小越好% 平顶区域掩膜 mask zeros(N); mask(R R_eval) 1; % R_eval是平顶评价区域比设计半径取小一点 % 区域内光强和均匀度 I_focal abs(U_focal).^2; I_region I_focal(mask 1); uniformity 1 - std(I_region) / mean(I_region);均匀度能做到95%以上肉眼看起来已经非常平了90%~95%之间细看能看出一些纹路90%以下基本不合格。边缘陡峭度edge steepness看的是平顶边缘过渡带的宽度。我习惯定义一个“10%~90%边缘宽度”就是从光强最大值的10%上升到90%所经过的横向距离。这个值越小越好一般做到1~2个衍射极限以内就算不错。转换效率efficiency衡量能量利用率公式是平顶区域内功率除以入射光束总功率eta sum(I_focal(mask 1)) / sum(sum(I_focal));注意如果仿真里没有考虑系统损耗GS算法理论上能把效率推到90%以上。如果效率偏低先检查是不是目标平顶半径取得太大或者目标振幅的截断太激进。还有一个指标是相关系数或误差可以直接把仿真光强和目标光强做归一化相关数值越高说明整形精度越好。我常用的是RMS误差的百分比形式diff I_focal / max(I_focal(:)) - target_amp.^2 / max(target_amp(:)).^2; rms_error sqrt(mean(diff(:).^2)) * 100;一般RMS误差控制在5%以内就算可用。把这些指标固化成一个评价函数每次跑完仿真自动输出调参的效率能提升一个量级。4.2 从MATLAB的相位图到SLM硬件的关键环节仿真的相位图是连续的浮点数范围在[-π, π]之间波动或者0到某个大于2π的值。但SLM实际上只能加载有限灰度级的相位处理流程大概是第一步把相位归一化到0~2π之间用mod操作完成这对应SLM的2π相位调制范围phase_map mod(phase_map, 2*pi);第二步把0~2π量化到SLM的灰度级。大多数SLM是8bit也就是256个灰度级gray round(phase_map / (2*pi) * 255); gray uint8(gray);这里有个细节量化误差不可避免理论最大相位误差是半格灰度对应2π/512 ≈ 0.35%的相位误差。对大多数整形任务完全够用但如果要求极高的平顶均匀度就要考虑使用12bit的SLM或者做误差补偿。第三步加载到SLM之后光路里一定要加一个离焦透镜相位或者闪耀光栅相位来分离零级光。因为SLM的像素结构本质上是一个二维光栅不可避免地会产生零级衍射也就是未经调制的直接反射光。如果不把零级光和整形光束在空间上分开你会发现平顶光斑中心总是叠着一个又亮又小的光点。我习惯给相位图叠加一个一维闪耀光栅% 闪耀光栅沿x方向线性相位 grating_phase 2*pi * x / grating_period; % grating_period可按像素数设 phase_map mod(phase_map grating_phase, 2*pi);光栅周期越小零级光和整形光束分离的角度越大但衍射效率会下降且目标光斑可能进入二阶衍射。一般取周期8~16个像素比较平衡。4.3 实际光路调试里的几条经验仿真永远是干干净净的光路里全是“脏活”。我整理几个几乎每次都会踩的实际问题偏振方向必须调对。反射式LCOS SLM对入射光的偏振方向敏感如果偏振方向不对相位调制深度会严重缩水甚至失效。正规做法是在SLM前面加一个半波片旋转光路中激光的偏振方向同时观察SLM的衍射效率调到最大再固定。光束要先扩束。SLM的相位调制是像素级的意味着高斯光束要尽量覆盖SLM的有效区域才能充分利用像素数。如果光束只有1mm粗SLM上只有一小片区域参与调制其他区域都在“空转”整出来的平顶光斑分辨率会严重不足。我的经验是把光束扩束到1/e^2半径占SLM有效宽度一半以上再用仿真里对应的w0值。CCD看焦斑时注意避开饱和。焦面上的激光功率密度很高CCD稍不注意就饱和了一饱和光强的顶点被削平看起来“非常平”实际上完全不能反映真实平顶质量。我会先衰减光功率到CCD响应线性区再开始测量宁可暗一点也不能过曝。透镜焦面位置的确认也要仔细。仿真里假设透镜是理想傅里叶变换实际透镜有像差焦面位置可能有毫米级的偏移。用CCD沿光轴方向扫描找光斑最锐利的那个位置再开始做定量测量。5. 常见问题与排查技巧实录5.1 问题速查表我把这几年来碰到的典型问题整理成一张速查表碰到类似情况直接照着排查现象可能原因解决思路平顶中心凹陷相位初始化不当SLM中心区域死像素换二次相位初始化检查SLM面板平顶边缘振铃目标振幅用硬边窗函数换超高斯目标加边缘过渡带光斑质量远差于仿真光束未扩束有效调制面积不足扩束后重新标定w0平顶光斑叠加亮斑零级衍射未分离叠加闪耀光栅相位平顶半径小于设计值实测量到的w0比仿真值大按实际光斑重新算映射关系灰阶加载后效果不对SLM灰度-相位响应非线性做gamma校准建立查表平顶出现方向性条纹SLM像素栅格衍射干涉换用更大周期的闪耀光栅或离焦同一相位图两次加载效果不同SLM温度漂移或液晶响应不稳定增加预热时间环境恒温平顶边缘粗糙、能量分散目标半径超出衍射极限有效范围减小目标半径或增大透镜焦距1kHz以上刷新频率下效果差SLM响应时间跟不上降低刷新频率或换高速SLM5.2 走了弯路才明白的三件事第一件事w0的取值是整套仿真的地基。我早期用SLM的时候SLM厂家给的像素尺寸、有效面积都清清楚楚但我一直用激光器说明书里的束腰半径做仿真结果实验和仿真怎么都对不上。后来拿刀片法实际测了SLM面上的光斑尺寸才发现经过扩束系统之后光束已经变成了带像差的“准高斯”和理想高斯差不少。从那以后我每次做这类仿真第一件事永远是实测SLM面上的光斑分布拟合出一个真实w0喂给程序。仿真里的美好结果90%取决于输入参数和实验对得上而不是算法本身。第二件事GS算法不是越迭代越好要盯着收敛曲线看。有一次我把迭代次数从100改到500想着反正时间不长结果发现200次之后误差曲线反而缓慢上升。后来查了资料才明白GS算法在某些条件下不单调收敛可能从一个较优解振荡到另一个稍差的解。我现在的习惯是迭代过程中实时记录误差曲线如果发现先降后升就在误差最低点把相位取出来后面再继续迭代的意义不大。第三件事实验上的微小误差仿真里根本看不见。比如SLM面板的相位响应不是完美线性的灰度0和灰度128对应的相位增量不是严格的2π/256液晶的响应和温度有关开机半小时和开机两小时的状态完全不同。所以不要指望实验和仿真完全一致更合理的预期是“趋势一致、数值偏差在可接受范围”。我的做法是在仿真里故意注入一些噪声和误差比如给相位图加1%的随机扰动、把灰度量化位数强制降到6bit看看输出光斑性能退化多少。如果这种“抗噪测试”通过了再上光路心里就比较有底。5.3 一个能让平顶更平的加分技巧最后分享一个我自己调了很长时间才悟到的小技巧属于那种程序没写错、但就是效果差一口气时可以考虑的思路目标平顶半径微补偿。直接计算法和GS算法的推导都默认输出面的坐标系是严格线性缩放的但真实光路里CCD看到的焦面坐标和仿真里假设的傅里叶变换坐标会因为透镜像差产生微小非线性畸变。这种畸变在高斯光斑中心区域几乎看不出来但在边缘位置会让平顶的边缘变得“发毛”。我的解决办法是在仿真阶段就预补偿把目标平顶的半径按一个轻微的非线性函数修正比如让目标半径从边缘向中心略有渐变相当于光路里实际畸变的逆向操作。具体补偿系数要用实验标定一般也就是1%~2%的调整量但就是这点调整能让实验实测的均匀度从93%拉高到97%以上。这种优化思路比换算法本身更实用因为它解决的问题是算法模型之外的现实偏差。做光束整形这几年我的总体体会是GS算法和直接计算法更像是两个工具而不是两个答案。真要做事先拿直接计算法快速摸清楚参数范围再用GS算法精修最后上光路实测迭代。这套组合拳我用了很久稳得很。如果你正准备做类似的项目我的建议是最开始别贪心先把圆对称的平顶光斑做出来均匀度达到95%以上再考虑方形平顶、环形平顶或者动态切换。光学这东西一步一个脚印地验证比一口气追求好看仿真数据要靠谱得多。