类视网膜阵列结构的光纤相控阵方法研究

刘炜剑，曹 杰，郝 群，张佳利，杨九章，唐鸣元，程 阳

（北京理工大学 光电学院 机器人与系统教育部重点实验室 精密光电测试仪器及技术北京市重点实验室，北京 100081）

引言

光纤相控阵（optical fiber phased array，OFPA）技术具有扫描速度快、指向精度高、空间分辨率高、大功率激光合成以及多目标探测跟踪等突出优点，特别是利用其电控扫描技术，能够实现探测、搜索、跟踪、通信等多功能一体化系统[1-5]，在大功率激光相干合成等领域有着重要的应用前景[6]。

在OFPA 中，出射光束阵列布局对高功率激光相干合成性能有着重要的影响，近年来引起广泛的研究。传统的OFPA 阵列布局主要包括方形、等间距多环形以及不等间距多环形等。研究发现，二维方形阵列布局[7]在远场会存在许多栅瓣和旁瓣，从而影响相干合成性能[8]。由于其自身的结构特性，均匀多环形阵列布局在一定程度上提高了能量集中度[9]，压缩了峰值旁瓣电平[10]（peak side lobe level，PSLL），但是仍然会在远场相干合成过程中带来较多的栅瓣以及旁瓣，从而影响能量合成效率。根据OFPA 理论，为了获得较好的远场相干合成性能，环间阵元间距应小于光束波长的一半[11]。然而，由于光束波长较短，当前的加工制造技术难以使得环间阵元间距满足上述条件。除此之外，另一种可行的方法是通过形成非均匀阵列布局，以减少栅瓣以及旁瓣对于远场相干合成性能的影响[12]。一般通过改变阵列的填充因子来形成非均匀布局。填充因子定义为阵元束腰半径与阵元间距的比值，填充因子越大，远场相干合成性能相对越好。根据定义，改变填充因子形成非均匀阵列布局，既可以通过优化阵元间距实现，也可以通过改变阵元的束腰半径实现[13]。

非均匀多环形OFPA 阵列布局可以通过优化算法实现，目前主要包括基因算法[14]（genetic algorithm，GA）和粒子群优化（particle swarm optimization，PSO）算法[15]。粒子群算法具有较快的搜索速度，搜索性能在一定条件下优于基因算法，因此在优化OFPA 环间阵元间距方面具有广泛应用。PSO 算法优化多环形阵列[16]能够有效地降低栅瓣以及旁瓣电平，充分体现了PSO 算法在曲线阵综合中的优势，相比较于其他传统阵列结构，远场相干合成性能可以得到一定程度的改善。除PSO 算法优化环间阵元间距形成非均匀多环形阵列布局之外，人眼视网膜结构本身也是一种非均匀多环形阵列布局，中央凹区域感光细胞密集，外围区域感光细胞稀疏，并且环间细胞尺寸按一定指数由内向外增长，该结构布局具有一定的冗余压缩作用[17]。受这一特性的启发，为了减少OFPA 远场相干合成过程中存在的栅瓣以及旁瓣，本文在保持中心阵元以及最外环阵元位置不变的情况下，将类视网膜结构以及PSO 算法相结合，设计了一种新型的非均匀多环形OFPA 阵列布局。与OFPA的其他传统阵列布局相比，该方法显著降低远场相干合成的栅瓣以及旁瓣，与此同时，在一定程度上提高了能量集中度，压缩了PSLL。

1 类视网膜多环形结构建模

本文提出的类视网膜阵列结构多环形OFPA由M环，每环N个阵元组成，其中，阵元在每个环上等间距排布，并且有一个阵元位于环中心位置。类视网膜多环形[18]OFPA 阵列布局如图1所示。假设类视网膜多环形OFPA 阵列出射光束束腰半径为wj（j=0，1，2，…，M），环间阵元间距为da(a=1,2,···,M)。

图1 类视网膜多环形OFPA 结构Fig.1 Diagram of retina-like multi-rings OFPA structure

对于只有一个阵元的阵列，m=0。对于m≥1 的OFPA 阵列，其中1≤n≤mn。第m个环上第n个阵元的坐标（xmn，ymn）可以表示为

类视网膜多环形OFPA 阵列的数量N可以表示为

类视网膜多环形OFPA 阵列布局中的每个阵元的出射光束一般服从高斯分布，输出平面远场分布可表示为

式中：Amn为出射光束阵元幅度；ϕmn为出射光束阵元初始相位；*为卷积运算符。

根据弗劳恩霍夫衍射理论，远场分布等同于近场分布的傅里叶变换。因此，远场分布E（x，y）可以表示为

将（3）式代入（4）式可得：

式中：（x，y）是输出平面的坐标；（x′，y′）是远场观察平面的坐标；fx′≈θx/λ以及fy′≈θy/λ是空间频率。

省略常数项，类视网膜多环形OFPA 阵列布局的远场光强度分布可表示为

式中k是波矢。

2 类视网膜多环形结构优化

为了有效抑制远场合成过程中存在的栅瓣以及旁瓣，应用优化算法搜索理想的阵元分布。其中最为常用的是基因算法和粒子群算法。它们都随机初始化相控阵阵列，使用评估函数来衡量阵元个体的优缺点，并基于从评估函数获得的自适应值执行一定的随机搜索。但是，基因算法有交叉、变异等操作，计算复杂度相对较高，并且很难利用阵元个体在解空间中的随机速度来改变搜索进程，因此，不具有较强的随机性。然而，粒子群算法具有较快的搜索速度，并且在搜索性能以及计算复杂度上优于基因算法[19]。因此，本文采用PSO 算法优化环间阵元间距，优化多环形OFPA 环间阵元间距的PSO 算法流程图如图2所示。

图2 多环形OFPA 环间阵元间距的PSO 算法流程图Fig.2 Flow chart of PSO algorithm for multi-rings OFPA inter-ring array spacing

3 仿真结果与分析

3.1 仿真设置

为了研究本文提出的PSO 算法优化类视网膜多环形OFPA 阵列布局的远场相干合成性能，进行相应的数值模拟仿真。设置关键参数如下：阵元总体数量N= 79，初始惯性权重系数ω=0.9，然后逐渐收敛至0.4[20]，阵元加速常数c1=c2=2，阵元速度变化区间vmax=0.01，vmin=-0.01，发射光波长（λ）设置为当前OFPA 技术最常用的1.064 μm波长。在均匀间距多环形阵列以及PSO 算法优化多环形阵列中，将束腰半径w设置为1.20 μm，因此，79 阵元总的束腰半径为94.8 μm。为了确保总的束腰半径不变，在类视网膜多环形阵列布局中，将中间环（第3 环）的束腰半径设置为1.20 μm，环间增长系数q为1.02。因此，设置各环束腰半径数值如表1所示。类视网膜多环形OFPA 阵列的总束腰半径为95.64 μm，与均匀间距多环形阵列以及PSO 算法优化多环形阵列相比，总束腰半径几乎没有变化。此外，将初始环间阵元间距设置为d= 4.8 μm。然后，在中心阵元以及最外环阵元位置不变的情况下，使用PSO 算法优化其他环间阵元间距。为确保优化的OFPA 是圆形对称的，仅优化相邻环之间的间隔，每环中的阵元仍然均匀分布。

表1 各环阵元束腰半径数值Table 1 Value of waist radius of each ring array element

3.2 结果对比分析

根据上述仿真设置，分别对5 种不同的OFPA阵列布局进行远场相干合成性能仿真，仿真结果如图3所示。评价参数对比如表2所示。可以看出，本文提出的方法在性能上相对传统方法有明显提升，其中，能量集中度具有一定程度的提高，峰值旁瓣电平也有较好的压缩，具体表现在以下2 个方面。1）能量集中度方面，传统方形阵列结构为0.562，如图3（a）所示；普通的多环形阵列结构相对有所提高，能量集中度为0.697，如图3（b）所示；PSO 算法优化多环形阵列结构，能量集中度为0.803，如图3（c）所示；单一的类视网膜多环形阵列结构，能量集中度为0.809，如图3（d）所示；本文提出的阵列结构如图3（e）所示，相比前者，在能量集中度方面分别提升了0.359、0.224、0.118、0.112。可以看出，相比较于方形阵列，均匀多环形阵列由于自身的环形对称结构，远场相干合成所含的栅瓣以及旁瓣相对较少，能量集中度已有较大的提升。一方面，在此基础上，利用粒子群算法优化，改变环间阵元间距形成非均匀阵列布局，能量集中度有进一步的提升；另一方面，改变环间阵元束腰半径形成的单一类视网膜多环形阵列结构，相比较传统均匀多环形阵列布局，在能量集中度方面也有一定程度的提高；最后，本文提出的将粒子群优化算法与类视网膜阵列结构两者进行结合形成的新型非均匀阵列结构，相比其他几种传统的阵列结构，能够显著提升远场相干合成的能量集中度。从能量集中度方面的提升趋势可以看出，PSO 算法优化类视网膜多环形阵列相比较于方形阵列结构，能量集中度提升最大，均匀多环形阵列结构次之，PSO 算法优化多环形阵列结构和单一类视网膜多环形阵列结构远场相干合成能量集中度接近，对比前两种传统阵列结构提升相对较小，从而验证了本文所提方法的可行性。2）PSLL方面，前4 种阵列布局分别为0.212，0.153，0.111，0.097，本文提出的PSO 算法优化类视网膜多环形阵列为0.043，相比前4 种阵列布局，在PSLL 方面分别压缩了0.169、0.11、0.075、0.069。同样，与能量集中度类似，PSO 算法优化类视网膜多环形阵列相对于其他几种传统阵列布局，在PSLL 方面也有较大幅度的压缩，从而在仿真层面也验证了本文方法的优越性。

表2 5 种不同OFPA 阵列布局评价参数对比Table 2 Comparison of evaluation parameters for five different OFPA array layouts

图3 5 种不同OFPA 阵列布局及其归一化强度的仿真Fig.3 Simulation of five different OFPA array layouts and normalized intensity

通过PSO 算法优化后，类视网膜多环形OFPA环间阵元间距分别为5.58 μm、4.17 μm、5.69 μm、4.35 μm、4.33 μm 以及4.69 μm，分布如图4所示。

图4 PSO 算法优化后，类视网膜多环形OFPA 环间阵元分布Fig.4 Distribution of retina-like multi-rings OFPA interring array element after optimization by PSO algorithm

4 结论

本文提出了一种类视网膜结构阵列布局的OFPA，在保持中心阵元以及最外环阵元位置不变的情况下，利用粒子群优化算法优化其他环间阵元间距，使得这种新型的非均匀间距阵元布局打破了光束波长的限制。与传统的OFPA 阵列布局相比，PSO 算法优化后的类视网膜多环形OFPA 阵列布局在一定程度上减少了远场合成过程中产生的栅瓣以及旁瓣，与此同时，提高了能量集中度，压缩了PSLL，该方法可为提升光束相干合成性能提供指导，有利于其在诸如激光通信、激光雷达等领域的应用。