基于粒子群算法PSO的宽带消色差超透镜的Matlab核心程序与FDTD仿真

基于粒子群算法PSO的宽带消色差超透镜 matlab核心程序 FDTD仿真

超透镜这玩意儿最近几年火得不行，尤其是宽带消色差的版本。传统透镜堆叠厚度感人，现在用亚波长结构直接在平面上搞定相位调控，想想就刺激。不过要同时搞定多个波长可不是简单活，这时候就得请出优化算法里的老江湖——粒子群算法（PSO）。

先看核心问题：如何在可见光范围内让不同波长的光都能准确聚焦。手工调参基本等于大海捞针，用PSO的关键在于定义好目标函数。咱们的目标是让三个特征波长（比如450nm、550nm、650nm）的相位分布尽可能接近理想情况。

基于粒子群算法PSO的宽带消色差超透镜 matlab核心程序 FDTD仿真

Matlab里搞PSO的套路一般是这样的：先初始化一群"粒子"，每个粒子代表一组超结构参数（比如纳米柱的直径、高度、周期）。这里有个坑，超结构参数的变化范围得设合理，不然算法跑偏到姥姥家去了。比如纳米柱直径一般控制在50nm到250nm之间，周期别超过工作波长。

% 粒子初始化 particle.position = lb + (ub - lb).*rand(1, D); particle.velocity = zeros(1, D); particle.best.position = particle.position; particle.best.cost = inf;

适应度函数才是重头戏，得把FDTD仿真结果揉进去。这里建议把全波仿真结果提前做成查找表，实时跑FDTD会等到地老天荒。举个栗子，可以预存不同结构参数对应的相位延迟，算适应度的时候直接插值调用。

function cost = fitness(x) % x是当前粒子位置（结构参数） phase = query_phase_lookup(x); % 查预存的相位数据 target = calculate_target_phase(); % 理想相位分布 cost = norm(phase - target) + 0.2*max(abs(phase - target)); % 后面这个0.2的惩罚项专治个别波长摆烂 end

粒子更新这块有个小技巧：惯性权重别用固定值。开场给个0.9让大伙多探索，慢慢降到0.4精细调整。这样前期不容易陷在局部最优，后期又能稳定收敛。

w = 0.9 - (0.5 * iter/maxIter); % 线性衰减惯性权重 particle.velocity = w * particle.velocity ... + c1*rand().*(particle.best.position - current_pos) ... + c2*rand().*(global_best.position - current_pos);

跑完优化别急着开香槟，得上FDTD做最终验证。这时候要注意mesh设置，特别是纳米结构边缘的网格得足够细。建议用共形网格技术，既能保证精度又不至于算到明年。

% FDTD仿真设置示例 setup = struct(); setup.mesh_accuracy = 3; setup.mesh_type = 'conformal'; setup.sources = {'type','total-field', 'wavelength', [450e-9, 550e-9, 650e-9]}; setup.monitors = {'xz_plane', 'frequency', 'all'};

最后看结果时重点盯着两个指标：聚焦效率（别低于60%）和焦斑尺寸（要比衍射极限小）。实测中发现某些波长的性能会突然扑街，这时候得回看适应度函数里的权重分配是不是偏心眼了。

搞这玩意最大的感悟就是：算法调参和电磁仿真得玩二人转。有时候看着算法收敛得挺美，一上仿真就翻车，这时候就得在目标函数里加物理约束，比如避开高次模激励的区域。这行当没有银弹，多试错才是王道。