1. 引言
射电望远镜阵列与光学望远镜成像原理并不相同,光学望远镜可直接成像,而前者是通过获取辐射源在空间域的可见度数据再傅里叶反变换转为图像亮度数据,因此最终成图质量几乎取决于空间域的覆盖程度。uv域全覆盖时,射电成像系统才有可能接收到完整的图像,然而这个只有天线无限多时理论上才可实现,但在实际中是不可能实现的;天线阵列分布与uv覆盖息息相关,uv覆盖 [1] 与脏束或psf是傅里叶变换对,不完整的uv覆盖可能导致接收的脏图无法通过clean算法以及其变体的算法来恢复出源图 [2] [3] [4] [5] 。因此诸多天文领域学者通过阵列布局来优化uv覆盖或脏束 [6] 。
模拟退火算法、爬山算法是在元启发算法出现之前的主要解决优化问题的传统优化算法 [7] 。因为这些算法需要从单点出发检索并且需要梯度信息,导致其优化需要大量时间。所以在解决诸多优化问题上都并非良选,如定位优化 [8] 。而元启发算法出现后,许多优化问题得以解决。它们与传统优化算法之间最大的差别是这类算法几乎不需要过多推导 [9] ,并且这类算法能有效避免局部最优。由于元启发算法的简单性、灵活性、低复杂度等特性,使其能应用于各个领域 [10] 。迄今为止开发的突出的自然启发的元启发式算法有萤火虫算法 [11] ,粒子群算法 [9] ,蝙蝠算法 [12] ,以及前沿的、应用领域广泛的、优异的灰狼优化算法 [13] [14] [15] [16] 等。
在很多应用问题上,元启发算法都表现良好,如线性天线阵优化 [17] [18] 、二值化技术 [19] 、路径规划 [20] 、同心圆形天线阵列优化 [21] 等。伴随着此类算法的快速发展,综合孔径天线阵列优化问题也逐渐开始用元启发算法来解决问题。Rahmat-Samii等人在2008年将粒子群优化应用于射电天文学应用的干涉阵列 [6] 。该实验表明在三个臂上不均匀地部署天线相比于均匀地布置天线能提供更好地的uv覆盖。该算法的弊端在于要获取更多的uv覆盖需要对每个观测单独运行。Kiehbadroudinezhad等人在2014年研究了更复杂的问题,即在角度域(l-m域)中的抑制SLL和空间频率域(uv域)中均匀覆盖 [22] 。其使用的是遗传算法进行优化并取得不错的效果,但遗传算法存在稳定性差、复杂度高等问题,而灰狼算法在诸多领域都表现良好 [13] ,因此本文基于灰狼算法进行天线布阵优化研究。
为了解决以上问题,并且考虑现有元启发算法应用到综合孔径二维阵列优化领域的很少,本文首先对现有的VLA_A阵列(由27根天线构成)进行分析,结果发现VLA_A阵列形成的UV Density数值较大。因此本文通过在优化综合孔径二维阵列领域较突出的遗传算法开展27根天线布阵优化,结果表明有所提升。但由于遗传算法存在计算成本高、稳定性差等问题,本文引入了元启发算法中前沿的灰狼算法进行了再次优化,与遗传算法对比,结果和花费时间均比其表现优异。并且考虑到标准元启发算法基本都是采用线性收敛因子,它的缺点是虽然前期探索能力强,但后期局部开发能力弱 [23] [24] ,为了得到前期探索能力强,后期局部开发能力也强的收敛因子,提出了分段收敛因子,并应用到标准灰狼算法上进行改进,最终结果表明,分段收敛因子的灰狼算法优化结果优于标准灰狼算法所优化的结果。
2. 综合孔径成像原理
综合孔径天线阵列的阵型的选择在天线数目不变的情况下,很大程度上取决于采用函数的优化 [1] 。n个干涉仪对应n个空间频率分量,干涉仪的天线之间的基线不同时,其对应的u, v空间频率分量也不同。因此通过增加天线数量能获得更多的空间频率分量。但实际上设计射电望远镜时,由于成本、技术等因素,天线数量都是有限的,要在天数数目有限的情况下获得更多不同的空间频率,天线的布阵尤为重要。
2.1. u, v覆盖
干涉仪的基线不同,对应的空间频率u, v也不同,基线与u, v之间的关系的数学模型如下:
(1)
其中,
、
、
分别为基线的三个分量,相位中心的指向为(h,
),h表示当地的时角,
表示赤纬。空间频率分量w的大小决定着成像视场的扭曲程度 [7] ,因此当存在大视场问题时,必然需要w项修正,大多数方法是将三维数据投影或近似到二维平面上,在二维平面上可以使用二维傅里叶变换,从而返回到现有算法可以解决的情况 [23] 。
综合孔径天线数量为n时,可构成
架干涉仪,在快照模式下这
条基线能获取
个
点。当地球自转时,u和v的值会发生变换,因此利用地球自转获取更多频率信息的工作模式逐渐得以发展,这种工作模式称为“非快照”工作模式。
2.2. 采样函数和脏束
干涉仪的基线矢量对应UV平面上的一个采样点,天线阵列相对应的UV覆盖可定义为该阵列配置参数的空间复相关函数,公式如下:
(2)
其中,f指取遍所有
点,
为有限长的冲击响应。采样函数对应的傅里叶逆变换称为脏束或者点扩展函数(PSF)。PSF代表望远镜系统性能,观测图的质量好坏很大程度上取决于它。在射电天文学中,用于消除望远镜点扩散函数(PSF)影响的最流行算法是CLEAN算法的变体 [2] [3] [4] ,比如Cotton-Schwab CLEAN (CS-Clean),基于峰值信噪比(psnr)来优化最陡下降法在恢复过程中的步长的CLEAN改进方法等。但这些方法都是PSF旁瓣效应在可接受程度上有效,因此脏束最终的形成十分重要,而通过天线阵列排布是改善脏束最直接也是最有效的方法。
3. 本文方法
3.1. 灰狼优化算法
灰狼优化算法(Grey Wolf Optimizer, GWO)是2014年澳大利亚Mirjalili等学者受到灰狼捕食猎物活动的启发而开发的一种元启发优化算法 [13] 。在自然中灰狼群常分为4个等级,第一阶级的灰狼用α表示,第二阶级的灰狼用β表示,第三阶级的灰狼用δ表示,第四阶级的灰狼用ω表示。ω灰狼数量在狼群中占比最多,并且其捕猎行为基本取决于前三阶级灰狼的影响。GWO由于较强的收敛性,被用到诸多优化领域,其捕猎数学模型如下。
1) 包围过程。灰狼发现目标猎物后,带领狼群向猎物方向进行包围,过程数学如下:
(3)
(4)
其中,t为迭代次数;
为目标解在t次迭代后的位置向量;
为灰狼t次迭代后的位置向量;D为灰狼与猎物的距离;A与C表示向量系数,两者定义如下:
(5)
(6)
(7)
其中,
和
是在[0, 1]之间均匀分布的随机数,a取值范围是[0, 2],
是最大迭代次数,本文取值
。
2) 捕猎过程。
头狼α带领β、δ对猎物进行捕猎,捕猎过程中,α、β、δ的行为用数学模型表示如下:
(8)
(9)
(10)
(11)
(12)
(13)
(14)
其中
、
、
分别表示为头狼α、β、δ的位置,
和
分别为α、β、δ的系数。由于
和
是线性收敛因子,这样虽然前期探索能力强,可后期局部开发能力却很弱,因此本文对传统灰狼算法又加以改进。
3.2. 改进的灰狼优化算法
标准灰狼优化算法的线性收敛因子
,其虽然前期探索能力强,但后期局部开发能力弱
[23] [24] ,为了得到前期探索能力强,后期局部开发能力也强的收敛因子,提出了分段收敛因子,其数学模型如公式15所示:
(15)
Figure 1. Comparison data of convergence factor
图1. 收敛因子对比数据
由图1可清晰看到,分段收敛因子和线性因子前期变化速率一样,这保证了两者前期探索能力一样强,而后期分段收敛因子的变化速率和非线性因子变化速率一样比较平缓,使得其后期局部开发能力和非线性因子一样强。这有效的提升了灰狼算法的优化能力。
3.3. 利用改进的鲸鱼算法进行天线阵列排布
VLA_A是一个由27根蝶形天线组成的Y型天线阵列。本文只研究其阵型,因此工作频率、天线孔径直径并未与VLA_A保持一致。
我们的目标是在找到能采样到更完整的频率信息的天线排布。由于天线的排布丰富,遗传算法能有效应用于此,Kiehbadroudinezhad等人利用遗传算法进行30根天线优化排布,最终优化后得到的脏束的SLL和uv覆盖率都还不错。因为遗传算法在元启发算法中其复杂度较高、收敛性一般,而灰狼算法复杂度较低、收敛性较强 [13] ,因此本文采用了灰狼算法应用到该排布问题上,该算法进行优化能很好的解决此类问题。
3.3.1. 初始化
天线所在位置本文用经纬度表示,为了统一单位,将其转为单位米。由经纬度坐标最终转为矩形坐标如下所示:
(16)
其中d是基线距离,az是基线的方位角,lat是天文台维度,el是基线仰角。
灰狼算法和粒子群以及遗传算法获取初始化点的方法都是一样的,如果不给予确定的初始点,都是通过随机初始化获取。
3.3.2. 评估
在空间频率域,我们采用了Babak E. Cohanim [24] 定义的参数M,也称UV Density,作为评估指标。其模型如下所示:
(17)
其中
表示已填充的所有标称网格点,
表示uv点总数,
,M的值越小,uv覆盖率更完整,即M的值越小,设计在优化中表现越好。
4. 结果
该算法在Matlab R2020b中实现,使用的计算机Intel(R) Core(TM) i7-6700 CPU @ 3.40GHz 3.41 GHz,RAM为8GB。
如下图2是VLA_A相关配置:
(a) (b)
Figure 2. (a) VLA_A position distribution; (b) KAT-7 uv coverage
图2. (a) VLA_A位置分布;(b) KAT-7 uv覆盖
如图2所示,VLA_A的uv覆盖并不理想,其UV Density是0.547,这表示其一半以上都是重采样点,这有很大提升空间。由于遗传算法在优化天线阵列方面表现良好,基于此利用遗传算法进行优化结果如图3所示。
(a) (b) (c)
Figure 3. (a) Genetic algorithm optimised curve with M as fitness; (b) Antenna configuration obtained from genetic algorithm optimization; (c) Genetic algorithm optimised antenna corresponding to uv coverage
图3. (a) 遗传算法以M作为适应度优化曲线;(b) 遗传算法优化所得天线配置;(c) 遗传算法优化天线对应uv覆盖
由图3(a)结果可知,使用遗传算法优化27根天线在uv域的UV Density,最终优化结果为0.4017,该结果小于VLA_A在uv域的UV Density,这表明了遗传算法优化所得的天线阵列相比于VLA_A能获取更多的频率信息。但这个值还有降低的空间,并且考虑到灰狼算法的收敛性效果不错,因此利用灰狼算法来进行优化求解,结果如图4所示。
(a) (b) (c)
Figure 4. (a) Optimisation curve of the grey wolf algorithm with M as the fitness; (b) Antenna configurations optimised by the grey wolf algorithm; (c) Optimised antenna corresponding to uv coverage by the grey wolf algorithm
图4. (a) 灰狼算法以M作为适应度优化曲线;(b) 灰狼算法优化所得天线配置;(c) 灰狼算法优化天线对应uv覆盖
由图3(a)和图4(a)可知,使用灰狼算法优化27根天线在uv域的UV Density的值小于使用遗传算法优化27根天线在uv域的UV Density。并且从图3(c)和图4(c)的uv覆盖图也能明显看到图4(c)的uv覆盖重叠采样点少于图3(c)的uv覆盖的重叠采样点。这表明了图4(c)的uv覆盖能获取更多的频率信息。
由于传统灰狼算法采用的是线性收敛因子,局部开发能力较弱,将线性因子优为分段收敛因子进行优化结果如图5所示。
(a) (b) (c)
Figure 5. (a) Improved grey wolf algorithm optimising the curve with M as the fitness; (b) Optimised antenna configuration from the improved grey wolf algorithm; (c) Optimised antenna corresponding to uv coverage by the improved grey wolf algorithm
图5. (a)改进灰狼算法以M作为适应度优化曲线;(b) 改进灰狼算法优化所得天线配置;(c) 改进灰狼算法优化天线对应uv覆盖
对比图5(a)和图4(a),可知,改进后的灰狼算法优化27根天线阵列在uv域的UV Density最终为0.3048,该值低于传统灰狼算法优化所得的UV Density。
将遗传算法和灰狼算法以及改进后灰狼算法优化27跟天线阵列排布最终所花费时间和优化所得指标UV Density记录于表1。
Table 1. UV Density obtained by the three algorithms to optimise the antenna array with their corresponding time spent
表1. 三个算法优化天线阵列所得的UV Density与其对应所花费时间
从表1可清楚看到,遗传算法不仅优化性能表现差于灰狼算法,所花费的时间也高于灰狼算法。而改进后的灰狼算法优化所花费时间虽然增加了一些,但性能相比于传统灰狼算法性能却有所提升。牺牲一些时间来提高优化性能,这是可接受的。
5. 总结
本文提出了一种改进灰狼算法来优化综合孔径天线阵列排布,结果表明改进后的灰狼算法不仅能有效解决天线阵列优化问题,而且优化效果均比标准灰狼算法和遗传算法表现优异。
对于天线的排布优化,元启发算法应用到该领域的很少,其中比较突出的是Kiehbadroudinezhad等人利用遗传算法进行布阵优化研究,其优化的是30根天线组成的阵列,与GMRT阵列作对比得出利用遗传算法优化结果良好。但由于天线布阵研究不同的数量会存在巨大的差异性,并且由于遗传算法存在易陷入局部最优、复杂度较高等问题,因此本文研究27根天线布阵研究时,采用了元启发算法中前沿的、不易陷入局部最优的灰狼算法进行布阵优化,与遗传算法进行对比的同时也与现有的VLA_A阵列进行了对比研究,结果表明标准的灰狼算法能有效解决27根天线布阵优化的问题,并且表现优于遗传算法。之后考虑到标准的元启发算法都是采用线性收敛因子,存在后期局部开发能力弱的缺点,会导致最终优化结果不理想。因此本文提出了分段收敛因子对标准灰狼算法进行改进,结果表明提出的分段收敛灰狼算法优化27根天线布阵的表现优于标准的灰狼算法,其优化的DV Density减少为0.3048。通过改进后的灰狼算法在天线布阵优化的表现,我们相信在其他领域,其也能表现优异。
基金项目
国家自然科学基金(12273007, 11963003, 12242303);贵州省优秀青年科技人才计划项目(黔科合平台人才-YQK[2023]006);贵州省基础研究(自然科学)项目(黔科合基础-ZK[2022]一般143);贵州大学培育项目(贵大培育[2020]76号)。
NOTES
*通讯作者。