Multi-target Interference Localization Using Single Satellite Multi-beam Antenna Based on Compressive Sensing
-
摘要: 针对卫星干扰处理中的多目标定位问题,该文提出基于压缩感知的定位方法。该方法利用目标的空间稀疏性,以及多波束天线在不同信号源方向上的增益不同,仅需要测量接收信号强度便可实现多个干扰的位置识别。研究结果表明,定位性能与节点分布、目标个数、波束覆盖半径、判决门限有关。在给定参数及原对偶内点算法下,该方法可实现1~4个干扰源的空域定位,在信噪比为20 dB时定位精度达到7.7 km,优于经典的旋转干涉仪和空间谱估计测向方法。Abstract: To cope with the issue of locating multi-target in mitigating satellite interference, a localization method is proposed based on Compressive Sensing (CS). The sources of satellite interference can be identified by using Received Signal Strength (RSS) measurement only, relying on the spatial sparsity of the target source and the fact that multi-beam antenna has different gain at the position of interference. The conclusions show that positioning performance is related to node distribution, target number, coverage radius and decision threshold. Furthermore, over the Primal-Dual Interior Point (PDIP) algorithm, the simulation result represents that the target number is four under certain conditions, and the position accuracy is closed to 7.7 km with SNR of 20 dB. In addition, the study result also confirms that the proposed algorithm is better than the classic methods of Rotating Interferometer (RI) and Direction Of Arrival (DOA) estimation
-
1. 引言
随着在轨卫星数量的不断增多及卫星频谱资源越来越紧张,干扰问题日益突出,严重影响卫星通信网络正常稳定运行。为了减小和消除干扰造成的恶劣影响,通常需要采用卫星干扰源定位技术[1,2]。单星定位系统仅利用1颗观测卫星实施干扰定位,可摆脱双星定位中对邻星和卫星星历等遴选条件的限制[3,4]。
常用的单星定位体制有:(1)单星仅测向:利用相位干涉仪或空间谱估计算法对来波方向进行估计。其中,干涉仪法的测向范围受基线长度限制,容易产生相位模糊;同时由于其设备复杂、造价高,目前常应用于侦查卫星中[5]。空间谱估计方法可实现高测向灵敏度和准确度,但对星载天线阵列有很高的要求[6]。(2)利用运动学原理:星上接收机在不同位置通过对同一干扰源测量得到的多个多普勒频率或脉冲到达时间进行定位[7]。该技术难以有效应对多目标问题,且对卫星姿态及同步要求很高,因此只适用于低轨卫星或倾斜轨道卫星。(3)幅度比较式测向:利用天线的方向特性,按照接收信号幅度的不同,测定不同来波方向。该体制实现相对简单,但定位精度过度依赖天线方向图的测量精度,且无法同时处理多个同频干扰[8]。
近年来,压缩感知(Compressive Sensing, CS)理论[9]的兴起为我们提供了一种全新的视角去看待多目标问题。研究表明,基于CS的多目标定位方法已在无线传感器网络(Wire Sensor Networks, WSN)中得到了广泛应用,并凸显出比传统方法更好的定位性能[10-12]。但是该方法需建立大量的无线链路以实现对定位区域的充分覆盖,而卫星多波束天线的覆盖特性恰好为二者的结合提供了基础。
鉴于此,本文针对多目标定位问题,提出一种将压缩感知理论应用于多波束卫星的干扰定位方法,给出了明确的系统模型和算法流程;推导了定位误差的理论值并分析了影响误差的主要因素;通过软件仿真,验证了所提方法具备较好的定位性能。
2. 多目标定位模型
2.1 节点分布特性
空间网格化是构建多目标定位模型的基础。不失一般性,设已知存在K个目标干扰源的2维感知区域被离散成N个网格节点,M个卫星波束各自对应1个传感器,用以测量其覆盖范围内各节点处参考信号的接收信号强度(Received Signal Strength, RSS)。以圆形波束为例,其覆盖单元可用正六边形的蜂窝结构(等效边长为l),且各图形之间不存在重叠或者空隙,图1所示为19个蜂窝网波束下的四色频分复用方式。
将整个卫星感知区域近似为圆形(易知其半径
R≥√19⋅l ),由于地面参考站数量有限,大量参考信号由己方移动用户提供,故将任一参考信号坐标(X,Y) 看作随机过程,从而有X(t)=I⋅cos(2πt+ω)Y(t)=I⋅sin(2πt+ω)} (1) 其中,
−∞<t<∞ , I和ω 是相互独立的随机变量,且分别服从区间[0,R] 和[0,2π] 上的均匀分布。在空间网格点数目足够大的条件下,可以用节点坐标近似估计目标的位置来对K个干扰源进行定位。2.2 压缩感知定位原理
由于
K≪N ,所以干扰源的位置在空域上具有稀疏性。根据CS理论,多目标问题可以抽象成通过M维接收信号强度的有噪测量(卫星信道中近似为高斯白噪声),重构出N维空间中K稀疏位置矢量的近似问题。基本模型为y=FYx+e=Ax+e (2) 其中,e为噪声矢量。其他参数意义如下:
(1) 干扰源位置矢量
x∈RN×1+ :x=(x1,x2,···, xN)T 为K稀疏矢量,编码了干扰源在动态节点中的近似位置,其非零元素表示与参考信号发射功率有关的正数k。(2) 基于RSS的稀疏基矩阵
Y∈RN×N :{{Y}} \!= \left( {\varphi _1}, {\varphi _2}, ··· ,{\varphi _N} \right) 的每一列φj 代表所有网格点对节点j处信号的RSS测量。(3) 测量矩阵
F∈RM×N :该矩阵表示节点对应传感器的部署方案。与地面WSN不同的是,多波束卫星传感器的位置并不是随机均匀散布在感知区域内,而是固定在转发器上,成为卫星通信有效载荷的一部分,宏观上相当于将WSN由2维变换为3维,利用多波束天线的方向特性起到分集接收的效果。(4) 感知矩阵
A∈RM×N :该矩阵由基矩阵和测量矩阵相乘得到,其任一元素ai,j 表示波束i处传感器对节点j处信号的RSS测量,可表示为RSSi,j=P⋅GL⋅At (3) 其中,P为信号发射功率,
At 为雨衰、大气吸收、电离层闪烁等引起的附加衰减,与实时链路状态相关。G表示波束增益,可近似计算为[13]G=G0⋅[J1(t)2t+36J3(t)t3]2 (4) 其中,
t=2.07123sinθ/sinθ3dB ,θ 为发射信号入射方向与波束中心的夹角,θ3dB 为波束中心增益3 dB衰减处对应的角度,G0为波束中心增益,J1, J3分别表示1阶、3阶第1类贝塞尔函数。L表示自由空间传播损耗,可近似计算为L=(4πdfc)2 (5) 其中,c代表光速,d为链路长度,f为载波频率。
(5) 观测矢量
y∈RM×1 :观测矢量为M个传感器的实际RSS测量结果。由于多波束天线具有高灵敏度,点波束旁瓣值较高,导致邻近波束之间有较大的重叠区域。当系统受到干扰时,除主要受干扰波束接收到干扰信号外,其邻近波束也会接收到该信号(即带外信号)。因此,y=(y1,y2,···,yM)T 中的任一元素yi代表第i个传感器感知到的所有目标的功率叠加结果。2.3 算法流程
图2展示了基于CS的多目标定位方法总体框架结构,分为离线和在线两个阶段:
(1) 离线阶段
步骤 1 建立坐标系。以1号蜂窝单元中心为原点建立平面坐标系,并以此推算各波束中心指向坐标;
步骤 2 设计节点分布方案。按照2.1节方法在圆形覆盖区域内随机产生N个动态节点,并记录所有节点坐标;
步骤 3 构建感知矩阵。将参考信号置于各节点处,根据多波束天线方向图的分布特性,各点波束在不同信号源方向上的天线增益不同,使得所接收到的信号强度有较大的差异,因此按照式(3)计算各传感器处的RSS,具体实施为:归一化参考功率,则干扰功率等于功率系数k;入射角度、链路长度可结合三角函数关系,由各节点坐标、星地距离、地球曲率等联合计算得到;考虑到目标频率的随机性会形成感知矩阵的不唯一性,故将载波频率按照复用方式分别代入,对测量结果取平均值。
(2) 在线阶段
步骤 4 测量RSS。当系统检测到干扰时,对各传感器接收到的所有干扰信号进行观测,并统计功率叠加结果;
步骤 5 重构稀疏矢量。利用稀疏恢复算法(如正交匹配追踪(Orthogonal Matching Pursuit, OMP)、基追踪(Basis Pursuit, BP)、原对偶内点(Primal-Dual Interior Point, PDIP)算法等[14,15])求解目标干扰的位置矢量信息;
步骤 6 确定干扰位置。以通过门限判决的有效值对应的节点作为目标位置的近似估计,完成定位。
3. 定位性能分析
3.1 随机分布误差
由于采用动态网格化的定位方法,故当定位成功时,估计误差
ε 可等效为任意两个相邻参考信号之间距离的期望。不妨设任意相邻节点的坐标为X1=I⋅cos(ω),X2=(I+ΔI)⋅cos(ω+Δω)Y1=I⋅sin(ω),Y2=(I+ΔI)⋅sin(ω+Δω)} (6) 其中,
ΔI 和Δω 分别服从区间[−2R/N,2R/N] 和[−4π/N,4π/N] 上的均匀分布。根据概率分布理论,有ε2=E[(X2−X1)2+(Y2−Y1)2]−D(ε)=E[(I+ΔI)2+I2−2(I+ΔI)I⋅cos(Δω)]−D(ε) (7) 其中,
D(ε) 为ε 的方差。已知下列算式成立E(I2)=R∫0x2Rdx=R23 E(ΔI2)=2RN∫−2RNx24RNdx=4R23N2 E(ΔI)=2RN∫−2RNx4RNdx=0 E[cos(Δω)]=4πN∫−4πNcos(x)8πNdx=sin(4π/N)4π/N 将式(8)代入式(7)可得
ε<R√43N2+23⋅[1−sin(4π/N)4π/N] (9) 由于
D(ε) 随N增大而减小,因此N越大,ε 越趋近于右边表达式。3.2 稀疏重构算法的影响
考虑到重构算法对定位性能的影响,还需从以下因素进行分析:
(1) 节点数N:在多目标定位问题中,基于CS的方法将参数空间离散化以获得目标所在节点的近似估计,显然N越大,节点与目标之间的误差就越小,式(9)也佐证了这一点;但同时N过大会导致稀疏基矩阵中各元素间的相关性增强,影响算法性能。
(2) 干扰数K:从压缩重构的角度看,给定节点分布方案确定了A的构成,从而进一步确定了能够重构出的最大目标数K,换言之,传感器个数M需满足[16]
M≥O[K⋅lg(N/K)] (10) 因此对于生成波束不变的卫星系统,K存在理论上限。
(3) 覆盖半径r:点波束覆盖半径越大,其与相邻波束之间的重叠区域就越大,不同传感器对同一干扰提供的信息就越多,但由于目标在感知区域内是稀疏分布的,故这个信息量是有限的;若波束间重叠区域很小,则各传感器近似独立工作,使信号重构效果变差。
(4) 判决门限T:通过对
ˆx 中的非零元素设定阈值,以减小因噪声和算法误差引起的“扰动项”。需要注意的是,由于目标数量和发射功率是未知的,利用重构算法进行信息恢复时,可能将一处干扰的RSS识别为多个目标共同作用的结果,因此随着T增大,ˆx 中可用信息的数量会减少,影响对定位结果的判断。综上,给出下述评价指标作为衡量定位性能的依据:
定位成功率
ρ :以估计位置为圆心,当实际目标处于任一搜索半径范围内视为定位成功,即能够被正确识别的干扰数所占的比例称为定位成功率。定位误差
δ :在定位成功前提下,满足搜索范围重叠区最小处的质心与实际目标的距离视为定位误差。4. 仿真实验
仿真参数配置如表1所示,其中频率设置的标准是3 dB波束宽度满足至少能覆盖1个蜂窝单元;循环迭代次数为500,每次执行时随机生成N个节点坐标和K个干扰坐标,同时随机分配频率和功率;不同仿真条件下只对相应参数作动态调整,其余不变;不考虑附加衰减影响,最后对统计结果进行仿真分析。
表 1 仿真参数配置固定参数 取值 调整参数 取值 波束数M 19 载波频率f(GHz) 3.60~3.75 天线口径D(m) 14 功率系数k 1~5 天线效率η 0.6 节点数N 400 蜂窝边长l(km) 240 干扰数K 4 区域半径R(km) 1050 覆盖半径r(km) 4l 搜索半径s(km) 30 判决门限T 0.2 参考功率P0 1 信噪比(dB) 20 4.1 计算干扰识别概率
图3为不同重构算法下干扰数识别概率的统计结果。通过计算干扰识别数的期望,可以得到在4目标随机分布条件下,OMP, BP和PDIP算法执行时成功定位的概率分别达到0.66, 0.70和0.73,充分验证了本文提出方法的可靠性,其中PDIP算法的性能相对最优;但针对上述重构算法的比较,并非本文研究重点,这里不再赘述。实际应用时可在条件允许情况下对RSS进行连续采样,从而利用多次观测降低随机误差对定位结果的影响。
4.2 不同参数下的性能验证
采用PDIP算法,图4为不同节点数下定位成功率、定位误差随干扰数变化的仿真结果。从图中可以看到,定位性能随干扰个数的增大而显著下降,这是因为位置矢量的稀疏度随之增大,算法的信号重构能力下降。由式(10)可知,对于M = 19和N = 100~400的部署方案,K的理论上限分别为7, 5, 4和4。图4(a)中,相同干扰数下N越大定位成功率越低,但当K分别取理论上限时各情形的
ρ 值相近;当K超过理论上限时,定位成功率均低于70%。图4(b)中,当K小于理论上限时,各情形的定位误差均满足式(9),且相同干扰数下N越大定位误差越小;但是随着K增大,受限于式(10)所示约束条件,N为400时δ 值的增长趋势甚至超过其他情形。上述仿真结论均与第3节相关理论一致,充分验证了分析的合理性。综上所述,当波束数(传感器数)固定时,节点数越大可获得更精确的估计结果,越小则能进行更多目标的定位,因此对N的选取应根据实际情况进行折衷考量。图5为不同覆盖半径下定位成功率、定位误差随判决门限变化的仿真结果。从图中可以看到,相同判决门限下,定位性能随覆盖半径增大而显著上升,这是因为各传感器之间协作程度随之增大,提高了信号重构的准确度;但当
r≥3l 时存在一个相对稳定的上界,这是由目标的K-稀疏分布决定的。图5(a)中,定位成功率随判决门限增大而缓慢下降,但在r≥2l 的情形中均保持在60%以上。图5(b)中,定位误差不随判决门限的变化而显著改变;当r = l时,δ 值在15 km上下波动,已经大于式(9)给出的理论上限,说明此时定位模型几近失效,根据误差定义,仿真曲线反映的即是搜索半径值的平均统计特性。上述仿真结论均与第3节相关理论一致,充分验证了分析的合理性。综上所述,判决门限对定位性能的影响有限,而波束覆盖半径的大小则是实现压缩感知多目标定位的关键因素。根据仿真结论将各评价指标与关键参数的大致对应关系汇总如下(见表2)。其中,“+”表示成正比,“–”表示成反比,“/”表示无直接关系。
表 2 评价指标与关键参数对应关系定位成功率ρ 定位误差δ 节点数N – – 干扰数K – + 覆盖半径r + – 判决门限T – / 4.3 与其他卫星定位方法的比较
将本文方法与相位干涉仪、空间谱估计(又称波达方向(Direction Of Arrival, DOA)估计)测向方法进行性能比较。针对传统干涉仪只适用于单目标情形,考虑利用旋转干涉仪(Rotating Interferometer, RI)处理多目标问题,即通过对干涉仪在不同旋转角度下的采样数据进行预处理,构造虚拟阵列数据对多目标参数进行估计;DOA估计中采用经典的MUSIC(MUltiple SIgnal Classification)算法,接收端考虑19馈源的3层环阵排布方式。由于RI和DOA均为针对来波角度参量(俯仰角、方位角)的估计方法,为便于比较,需要对测向结果进行转换。
假设确定干扰坐标的随机变量为
I 和ω ,测向时俯仰角、方位角误差分别为α 和β ,则估计坐标的参量可表示为ˆI=H⋅tan(ξ+α)ˆω=ω+β} (11) 其中,
tan(ξ)=I/H , H表示星地距离,这里取值为36000 km。利用式(11)得到的坐标可计算求解定位误差。此外,为更直观地描述低信噪比下的定位性能,本例中将搜索半径设置为500 km。图6为不同定位方法下定位误差随信噪比(Signal-to-Noise Ratio, SNR)变化的仿真结果。可以看到,随着信噪比的增大,目标参数的定位误差逐渐减小,3种方法在SNR = 20 dB时的定位精度分别达到30.2 km, 18.9 km和7.7 km;其中本文方法的定位性能相对最优,且在低信噪比段(0~5 dB)仍具备较强的竞争力。具体原因是RI和DOA法利用阵元测量相位差时,受通信卫星轨道高度的影响,其基线长度难以保证有效的测量精度;DOA法测向时还利用了信号的幅度参数,因此其定位性能优于RI法;而所提方法基于“同一地理区域信号在同一时间帧中具有功率变化的正相关”这一特征,在构造感知矩阵时能够最大化减少低信噪比的影响,只是在信号重构时会引入与噪声相关的“扰动项”,但当SNR > 5 dB时误差曲线趋于平缓,定位精度较高。
5. 结束语
本文提出将压缩感知相关理论应用于多波束卫星的干扰定位场景中,旨在解决空分频率复用下多个干扰目标位置的识别问题。该方法利用星载多波束天线的方向特性,通过动态节点化参考信号,只需测量接收信号强度及求解一个稀疏恢复问题,便可重构出包含多个目标位置信息的稀疏矢量。研究结果表明,定位性能与节点分布、干扰个数、波束覆盖半径、判决门限等诸多因素有关;在一定参数条件下,该方法在定位成功率及定位误差方面均能取得较好的定位效果,且其定位精度明显优于相位干涉仪和空间谱估计方法,因此对卫星干扰定位技术应用具有一定的参考价值。在此基础上,之后的研究工作可结合坐标系转换、链路衰减适配等方面深入展开,以进一步简化模型参数、提高定位性能。
-
表 1 仿真参数配置
固定参数 取值 调整参数 取值 波束数M 19 载波频率f(GHz) 3.60~3.75 天线口径D(m) 14 功率系数k 1~5 天线效率η 0.6 节点数N 400 蜂窝边长l(km) 240 干扰数K 4 区域半径R(km) 1050 覆盖半径r(km) 4l 搜索半径s(km) 30 判决门限T 0.2 参考功率P0 1 信噪比(dB) 20 表 2 评价指标与关键参数对应关系
定位成功率ρ 定位误差δ 节点数N – – 干扰数K – + 覆盖半径r + – 判决门限T – / -
[1] 易克初, 李怡, 孙晨华, 等. 卫星通信的近期发展与前景展望[J]. 通信学报, 2015, 36(6): 157–172. doi: 10.11959/j.issn.1000-436x.2015223YI Kechu, LI Yi, SUN Chenhua, et al. Recent development and its prospect of satellite communications[J]. Journal on Communications, 2015, 36(6): 157–172. doi: 10.11959/j.issn.1000-436x.2015223 [2] 郝才勇. 卫星干扰处理技术综述[J]. 电信科学, 2017, 8(1): 106–113. doi: 10.11959/j.issn.1000-0801.2017018HAO Caiyong. A survey of mitigating satellite interference technology[J]. Telecommunications Science, 2017, 8(1): 106–113. doi: 10.11959/j.issn.1000-0801.2017018 [3] 孙霆, 董春曦, 毛昱. 一种基于半定松弛技术的TDOA-FDOA无源定位算法[J]. 电子与信息学报, 2020, 42(7): 1599–1605. doi: 10.11999/JEIT190435SUN Ting, DONG Chunxi, and MAO Yu. A TDOA-FDOA passive positioning algorithm based on the semi-definite relaxation technique[J]. Journal of Electronics &Information Technology, 2020, 42(7): 1599–1605. doi: 10.11999/JEIT190435 [4] 孙光才, 王裕旗, 高昭昭, 等. 一种基于短合成孔径的双星干涉精确定位方法[J]. 电子与信息学报, 2020, 42(2): 472–479. doi: 10.11999/JEIT180940SUN Guangcai, WANG Yuqi, GAO Zhaozhao, et al. A dual satellite interferometric precise localization method based on short synthetic aperture[J]. Journal of Electronics &Information Technology, 2020, 42(2): 472–479. doi: 10.11999/JEIT180940 [5] 张金秀, 陶海红, 王渊. 一种基于双基线旋转的改进干涉仪定位算法[J]. 北京理工大学学报, 2018, 38(3): 320–324. doi: 10.15918/j.tbit1001-0645.2018.03.016ZHANG Jinxiu, TAO Haihong, and WANG Yuan. An improved locating algorithm based on double baselines rotating interferometer[J]. Transactions of Beijing Institute of Technology, 2018, 38(3): 320–324. doi: 10.15918/j.tbit1001-0645.2018.03.016 [6] ADEOGUN R O. A robust music based scheme for interference location in satellite systems with multibeam antennas[J]. International Journal of Computer Applications, 2013, 82(12): 1–6. doi: 10.5120/14165-2322 [7] 徐义, 郭福成, 冯道旺. 一种单星仅测TOA无源定位方法[J]. 宇航学报, 2010, 31(2): 502–508. doi: 10.3873/j.issn.1000-1328.2010.02.031XU Yi, GUO Fucheng, and FENG Daowang. A new satellite passive localization method using TOA measurement only[J]. Journal of Astronautics, 2010, 31(2): 502–508. doi: 10.3873/j.issn.1000-1328.2010.02.031 [8] HE Chao, XIE Zhidong, BIAN Dongming, et al. Study of interference localization using single satellite based on signal strength distribution in multi-beam antenna for satellite communications system[J]. International Journal of Distributed Sensor Networks, 2018, 14(5): 16–27. doi: 10.1177/1550147718774015 [9] CANDES E J and WAKIN M B. An introduction to compressive sampling[J]. IEEE Signal Processing Magazine, 2008, 25(2): 21–30. doi: 10.1109/MSP.2007.914731 [10] CEVHER V, DUARTE M F, and BARANIUK R G. Distributed target localization via spatial sparsity[C]. The 2008 16th European Signal Processing Conference, Lausanne, Switzerland: IEEE, 2008: 1–5. doi: 10.5281/zenodo.41257. [11] 余东平, 郭艳, 李宁, 等. 压缩感知多目标无源定位中的字典适配方法[J]. 电子与信息学报, 2019, 41(4): 865–871. doi: 10.11999/JEIT180531YU Dongping, GUO Yan, LI Ning, et al. Dictionary refinement method for compressive sensing based multi-target device-free localization[J]. Journal of Electronics &Information Technology, 2019, 41(4): 865–871. doi: 10.11999/JEIT180531 [12] 游康勇, 杨立山, 郭文彬. 无线传感器网络下基于压缩感知的多目标分层贪婪匹配定位[J]. 自动化学报, 2019, 45(3): 480–489. doi: 10.16383/j.aas.2018.c170237YOU Kangyong, YANG Lishan, and GUO Wenbin. Hierarchical greedy matching pursuit for multi-target localization in wireless sensor networks using compressive sensing[J]. Acta Automatica Sinica, 2019, 45(3): 480–489. doi: 10.16383/j.aas.2018.c170237 [13] CAINI C, CORAZZA G E, FALCIASECCA G, et al. A spectrum- and power-efficient EHF mobile satellite system to be integrated with terrestrial cellular systems[J]. IEEE Journal on Selected Areas in Communications, 1992, 10(8): 1315–1325. doi: 10.1109/49.166759 [14] CANDES E J and TAO T. Decoding by linear programming[J]. IEEE Transactions on Information Theory, 2005, 51(12): 4203–4215. doi: 10.1109/TIT.2005.858979 [15] BOYD S and VANDENBERGHE L. Convex Optimization[M]. Cambridge: Cambridge University Press, 2004. doi: 10.1017/CBO9780511804441. [16] CANDÈS E J, ROMBERG J K, and TAO T. Stable signal recovery from incomplete and inaccurate measurements[J]. Communications on Pure and Applied Mathematics, 2006, 59(8): 1207–1223. doi: 10.1002/cpa.20124 期刊类型引用(2)
1. 朱江,王彦敏. 无人机数据采集系统中基于动态干扰的轨迹优化方案. 数据采集与处理. 2024(05): 1271-1286 . 百度学术
2. 吴涛,胡艳霞,田甜,张健铤. GNSS干扰定位技术分析. 全球定位系统. 2023(05): 103-111 . 百度学术
其他类型引用(2)
-