相控阵馈源MVDR波束合成器缓解RFI仿真*

2023-02-01 08:39李健仵鹏段雪峰裴鑫张青松陈卯蒸
天文学报 2023年1期
关键词:射电视场波束

李健 仵鹏 段雪峰 裴鑫 张青松 陈卯蒸

(1 中国科学院新疆天文台 乌鲁木齐 830011)

(2 新疆微波技术重点实验室 乌鲁木齐 830011)

(3 新疆射电天体物理实验室 乌鲁木齐 830011)

(4 中国科学院大学 北京 100049)

1 引言

上世纪末射电工程师开始使用相控阵馈源替代传统喇叭式馈源以获得更大的观测视场,提升了射电望远镜的工作效率,缩短了观测时间[1].国际上美国国家射电天文台(National Radio Astronomy Observatory,NRAO)研制出19个阵元双偶极子PAFs (Phased array feeds)接收机,其数字波束合成器可以实时对波束进行合成[2]; 荷兰射电天文研究所(Netherlands Institute for Radio Astronomy,ASTRON)研制出由144个Vivaldi天线作为阵元的PAFs接收机,结合MaxSNR (Maximum Signal-to-Noise Ratio)和LCMV (Linear Constraint Minimal Variance)波束合成器对连续的大视场进行平滑性优化,提升了波束的旋转对称性,增强了望远镜观测时的稳定性[3]; 澳大利亚国家射电望远镜中心(Australia Telescope National Facility,ATNF)研制出由94×2个贴片偶极子单元组成188阵元PAFs接收机Mk.II,该接收机安装在埃费尔斯伯格的100 m射电望远镜上,经过望远镜校准和波束合成优化基本可以满足实际观测的需求[4].国内由上海天文台率先探索了在大口径射电望远镜使用PAFs接收机系统,该系统实现了对波束的控制功能[5]; 国家天文台和中电五十四所联合成立的射电天文技术联合实验室正在开展平方公里阵列(Square kilometre Array,SKA)先导天线及相控阵馈源等关键技术研究,已经设计出中国SKA项目中的PAFs接收设备样机[6].新疆天文台联合中国科学院和新疆维吾尔族自治区,拟在新疆奇台县建设“110 m口径全向可动射电天文望远镜” (QiTaiTelescope,QTT),初步计划采用20 cm波段(0.7-1.8 GHz)相控阵馈源接收机装配在QTT主焦位置以实现连续的大视场覆盖[7-8].目前国内还没有大口径射电望远镜应用PAFs接收机实现波束合成及RFI (Radio Frequency Interference)缓解的公开研究成果.

PAFs接收机作为一种新型的电磁波接收装置,在全世界范围内都展开了广泛的研究.相对于单波束接收机,虽然拥有良好的主瓣增益和更低的旁瓣响应,但是观测视场小、反射镜利用率低是其主要的缺点.相对于多波束接收机,在继承了单波束接收机高增益的基础上,虽然增大了观测视场、提升了反射镜的利用率,但是受限于馈源物理尺寸的影响,无法形成连续的观测视场是其无法避免的缺点.PAFs接收机使用众多的小型天线对反射镜焦面场的电磁波进行空域采样,返回时域阵列信号,采用灵活的数字波束合成器对望远镜的波束形状、指向和数量等进行控制,从而实现优化照明区域、选择性增强目标方向信号响应,同时抑制干扰方向信号响应,即在空域内实现RFI的缓解为单口径射电望远镜RFI缓解带来新的解决办法.相较于单波束接收机和多波束接收机而言,优化望远镜波束仅仅是天线设计的问题,但是对PAFs接收机而言,它也是一个波束合成器设计的问题,波束合成器通过设计可以实现高灵敏度、给定空域响应或其他设计目标.

Warnick等人在Green Bank Telescope 20 m射电望远镜上初步验证了使用波束合成器可以实现RFI缓解[9].MVDR (Minimum Variance Distortionless Response)波束合成原理可自主识别RFI的方向并给予缓解,在阵列信号处理中有着广泛的应用.Jeffs等[10]和Leshem等[11]在对PAFs接收机信号处理建模时就提到MVDR自适应波束合成器具有很好的应用和拓展前景,并做了一些深入的理论探索和仿真.总体看MVDR波束合成器在PAFs接收机中具有很好的应用前景,特别是在PAFs射电望远镜的大视场波束形成以及大视场下有源RFI缓解等问题带来了新的解决方法.

2 阵列信号处理模型

PAFs接收机可以理解为采用众多小型天线按照一定规则排布的阵列电磁波接收系统.因此PAFs接收机生成信号的处理方法满足阵列信号处理的基本模型.

2.1 平面阵列时域响应模型

在理想情况下,假设射电望远镜位于自由空间内,接收到的电磁波信号未受到遮挡与散射等影响,并且假设PAFs各阵元在工作频带内具有相同、平坦且各向同性的接收灵敏度.通常射电望远镜的主焦方向会对准目标源,而干扰源的电磁辐射往往从射电望远镜的旁瓣进入,射电望远镜观测示意图如图1所示.

图1 射电望远镜观测示意图Fig.1 Schematic diagram of radio telescope observations

假设有D+1个平面电磁波分别从球面角Ωd(d=0,1,··· ,D)入射到馈源,且存在加性背景噪声,则相对于参考阵元的第m个阵元观测到的信号波形可以描述为[12]:

式中,t为时间信息,M表示阵元数量,τm(Ωd)是信号沿Ωd方向到达第m个阵元相对于参考阵元的传播时间延迟,sd(t) (d=0,1,··· ,D)为在空间参考阵元测量的第D+1个信号波形,其中d=0对应于期望信号,d=1,2,··· ,D对应于D个干扰,nm(t)表示第m个阵元接收到的满足广义平稳特性的综合背景噪声.

2.2 MVDR波束合成器

Capon[13]于1969年提出了MVDR波束合成器,它可以使目标方向信号无失真通过的同时,最大限度地抑制波束输出中的干扰和噪声功率.1972年,Frost[14]提出了LCMV波束合成器,将期望信号不受影响的这一约束扩展为一组约束,由于RFI的方向无法提前预知,不能很好地实现提前约束.Max-SNR波束合成器可以实现输出信噪比最大化,但对干扰的方位和缓解没有识别和处理能力,更适合在没有RFI的环境下使用.

对PAFs各阵元接收到的信号进行相位补偿和线性加权后,可以实现射电望远镜的空域滤波能力,即形成一定方向性的接收,让目标方向入射的信号通过,同时抑制其他方向入射的干扰与噪声,这种信号处理技术称为波束合成.在射电天文中,通常希望最大限度地提高目标方向的电磁波响应,同时尽可能地减小其他方向的电磁波响应.MVDR波束合成器可以使目标方向信号在无失真输出的状态下,使望远镜输出的噪声与干扰功率达到统计意义下的最小化,从而实现对干扰和噪声的抑制,减小系统输出中噪声与干扰的响应.PAFs接收机中MVDR波束合成器原理框架如图2所示,其中LNA_1,LNA_2,···,LNAM表示各阵元对应的低噪声放大器,x1(ω),x2(ω),··· ,xM(ω)为各阵元窄带频域信号,ω为频域信息.阵元信号经过一定量的样本统计后传入MVDR运算器,再经过运算器处理后返回相应的匹配系数矩阵wH=[w1,w2,··· ,wM]H,其中wM表示第M个阵元的权重,H表示共轭转置,X(ω)=[x1(ω),x2(ω),··· ,xM(ω)]为PAFs阵元输出矩阵,最后这些系数与对应的统计样本加权后实现波束合成,即y(ω),波束合成器完整表达式为[10]:

图2 MVDR波束合成器框架图Fig.2 MVDR beamformer frame diagram

3 仿真

新疆天文台拟在新疆奇台县建设“110 m口径全向可动射电天文望远镜”,依据其科学目标,计划在QTT主焦点配备一台20 cm波段(0.7-1.8 GHz)PAFs接收机,以实现连续的高灵敏度大视场观测.根据项目进展规划,该台PAFs接收机采用96个阵元组成四边形平面阵列,阵元排布方式如图3所示[15].

3.1 仿真环境

使用Python仿真阵元输出信号和波束合成过程,并对结果进行可视化处理,仿真过程主要分为3个部分.

(1)信号发生器: 生成含有信号、干扰和综合噪声的PAFs阵列复指数信号流.对PAFs建立坐标系如图4所示,将阵元平面置于yoz平面,θ与φ分别为信号相对于PAFs的水平方位角与垂直俯仰角,dym为第m个阵元相对于参考阵元沿着y轴的距离,dzm为第m个阵元相对于参考阵元沿z轴方向上的距离.设目标所在的方向余弦表示为(cosθcosφ,cosθsinφ,sinθ),则平面电磁波到达各阵元之间的波程差决定了信号接收过程中的时间延迟,因此,第m个阵元相对于参考阵元沿y轴和z轴的时间延迟可分别表示为:

式中,c为电磁波传递速度,则方向Ω(θ,φ)的平面电磁波到达第m个阵元的时间与到达参考阵元的时间延迟可表示为:

代入(1)式可得平面电磁波在各个阵元相对于参考阵元的时域响应信号.

(2)波束合成器: 使用MVDR波束合成器计算统计样本下最优匹配系数组wMVDR,实现波束合成.在信号发生器的参数基础上输入主瓣方向参数: Ω0(θ0,φ0),生成波束导向矢量a(Ω0),波束导向矢量应当为(4)式所对应时间上的补偿,具体频域表达式为:

式中,exp(·)=e(·)表示指数函数,i=是虚数单位.

使用不依赖于样本统计的匹配系数作为固定波束合成器,该波束合成器用作与MVDR波束合成器性能对比,其恒定匹配系数wCONST为:[1,1,··· ,1] (此时波束指向为Ω(0°,0°),即是指向望远镜焦点的方向).

(3)结果可视化: 对波束合成器的输出信号遍历成像区域,生成望远镜二维响应图像.可视化主要参数为: 成像范围,图像分辨率.

参考图3的阵元排布方式确定对应的PAFs阵元坐标,PAFs的窄带工作频率为1.5 GHz,通常设定阵元间距与信号波长比例系数r=0.5,确定样本采样精度SP=4,样本总量N=1024.假定在理想反射条件下PAFs完整覆盖电磁波在抛物面反射形成的平面电磁场,因此使用所有阵元进行波束合成.定义信噪比SNR=其中Ps为样本信号功率,Pn为样本噪声功率.向信号发生器输入SNR=0 dB,θ=0°,φ=0°的信号参数,并对生成的所有阵元信号使用固定波束合成器进行波束合成,归一化后可得望远镜的波束方向图,可视化结果如图5所示.

图3 QTT PAFs 96阵元布局图(0.7-1.8 GHz)Fig.3 The layout of 96 elements of PAFs in QTT (0.7-1.8 GHz)

图4 阵列建模示意图Fig.4 Schematic diagram of element modeling

图5 QTT PAFs 96阵元在1.5 GHz下的方向图Fig.5 The patterning of PAFs with 96 elements of QTT at 1.5 GHz

3.2 有源RFI缓解仿真结果

向波束生成器同时输入两种不同参数的信号,分别作为一个目标源和一个干扰源,两个源的具体参数如表1所示.固定波束合成器的响应结果如图6(a)所示.将生成的N个样本量作为MVDR波束合成器的统计输入[10],MVDR波束合成器的输出即为该样本下满足RFI缓解的合成信号,在该合成信号下的望远镜响应如图6(b)所示.

表1 信号源与干扰源参数表Table 1 Signal source and interference sourceparameters list

图6(a)分析可得: 在固定波束合成器的成像范围内,干扰源的响应位于以φ=-0.25°和θ=-0.25°为中心的区域内,在分辨率内完全被SNR ≥45 dB的响应区间覆盖,外围主要被SNR ≥30 dB的响应区间覆盖; 信号源的响应被设定于φ=0°和θ=0°为中心的区域内,由于该区域的响应与图中大部分SNR ≥30 dB的响应区域处于同一动态范围,无法辨别是信号响应还是噪声和干扰引起的综合响应.干扰位置处有着明亮且最大范围的响应,目标位置处的信号响应类型无法识别.

图6(b)分析可得: 在MVDR波束合成器的成像范围内,干扰源的响应位于φ=-0.25°与θ=-0.25°的中心区域内,在分辨率内以SNR ≥15 dB的成像范围占据主要区域,外围更大的响应区域处于SNR ≥0 dB的响应范围; 反观位于φ=0°与θ=0°为中心的信号区域,在有效的分辨范围内主要以SNR ≥15 dB为主要响应区域,甚至于在最中心可看到SNR ≥30 dB的小范围响应区域.

图6 (a)有源干扰下固定波束合成器响应图,(b)有源干扰下MVDR波束合成器响应图.Fig.6 (a) Response diagram of conventional beamformer under active interference,(b) Response diagram of MVDR beamformer under active interference.

从SNR响应分布变化可得,干扰位置处的SNR得到了明显的抑制与缓解,SNR的抑制最大可达45 dB,最低也有30 dB; 目标位置处的SNR整体得到了相对较小的抑制,最大抑制仅为15 dB,中心点处实现了0 dB抑制,即无畸变响应.两个源相比较,虽然SNR均得到了抑制,但是等价于信号处的源得到了约25 dB的增强.从成像效果变化可得,图像整体能量分布下降,干扰与噪声的功率得到明显的缓解,目标信号得到明显增强,从原来的淹没状态中凸显而出,且全图能量最强点位于目标方向处,成像噪声整体得到改善.

在该仿真中RFI的方向是除焦点方向以外随机设定的,具备方向上的随机性,从仿真结果可得MVDR波束合成器拥有时变RFI缓解能力,这是因为MVDR波束合成器是建立在一定时间样本统计的基础上,所以缓解RFI在时间上的分辨率与样本统计的时间分辨率具有一致性.

3.3 噪声累积引起的RFI缓解仿真结果

当目标源非常微弱,且无有源干扰时,PAFs接收机系统的每个阵元接收到的信号通常被淹没在噪声里,MVDR波束合成器可以在一定程度内提取出每个阵元中的目标信号进行合成,从而增强目标信号响应,同时抑制波束合成时阵元之间加性噪声累积引起的干扰响应.

向波束生成器输入SNR=-30 dB、θ=0°、φ=0°的信号作为阵元输出信号[10],固定波束合成器的响应结果如图7(a)所示,将生成的N个样本量作为MVDR波束合成器的统计输入,MVDR波束合成器的输出即为该样本下满足RFI缓解的合成信号,在该合成信号下的望远镜响应如图7(b)所示.

图7(a)分析可得: 在固定波束合成器的成像中,有许多大小不一,呈现出随机分布的响应区域,这些响应处于SNR ≥20 dB的范围内,并且位于以φ=0°和θ=0°为中心的信号区域也处于该范围内,因此该响应是随机分布所引起的还是目标源的响应无法区分.主要原因是由于各阵元之间加性噪声在波束合成时累积引起的响应,这种响应在微弱的天文信号观测中会表现为一种无规律,相对强的干扰分布.

图7(b)分析可得:在MVDR波束合成器的成像中,目标源响应位于φ=0°和θ=0°为中心的区域内,在分辨率内以SNR ≥0 dB的动态范围为主要响应,在最中心处更有一部分成像区域处于SNR ≥15 dB的动态范围内.从全图来看,除目标区域以外均无其余响应高于SNR ≥15 dB的动态范围,加性噪声引起的干扰得到了明显的抑制.

图7 (a)无源干扰下固定波束合成器响应图,(b)无源干扰下MVDR波束合成器响应图.Fig.7 (a) Response diagram of conventional beamformer under passive interference,(b) Response diagram of MVDR beamformer under passive interference.

经过MVDR波束合成器处理后,从SNR动态变化范围可知,目标源的SNR处于最高动态范围,噪声累计引起的干扰全部得到了抑制; 从成像效果变化可知,图像整体能量分布下降目标信号得以凸显,满足基本处理预期.

4 总结与展望

综上仿真结果,MVDR波束合成器在基于一定量的样本统计后可以自适应地缓解RFI,提升望远镜的灵敏度,在强RFI干扰下,MVDR波束合成器可以缓解RFI的响应,提亮目标源.在进行弱源且无有源RFI的观测条件下,加性噪声在固定波束合成时降低了望远镜的灵敏度,MVDR波束合成器对这些噪声引起的干扰有很好的抑制效果.

本文围绕PAFs接收机系统单波束合成优化展开研究,根据波束合成原理以及阵列信号处理方法,对有源干扰源进行了MVDR波束合成优化,同时也对加性噪声在PAFs波束合成中的累积干扰进行了MVDR波束合成优化,针对波束合成结果均做了2维成像分析.MVDR波束合成器对非主瓣内的干扰有很强的识别和抑制能力,可以提升在RFI环境下射电望远镜的观测能力,同时单波束的自适应优化也为PAFs连续大视场优化及相应的RFI缓解提供了基础.

猜你喜欢
射电视场波束
一种晶圆自动光学检测系统的混合路径规划算法
谁能抓住“神秘天神”——快速射电暴?
射电星系
美国的绿岸射电望远镜
一种基于基准视场扩散拼接的全景图像投影方法
基于共形超表面的波束聚焦研究
超波束技术在岸基光纤阵中的应用
毫米波大规模阵列天线波束扫描研究*
医用内窥镜矩形视场下入瞳视场角的测试方法研究
Optimization of Structure Layout for Three-FOV Star Sensors*