基于星载可移点波束天线的地面辐射源定位方法 |
|||||||
申请号 | CN201610034916.8 | 申请日 | 2016-01-19 | 公开(公告)号 | CN105676166A | 公开(公告)日 | 2016-06-15 |
申请人 | 中国人民解放军理工大学; | 发明人 | 张更新; 钟锡健; 谢智东; 边东明; 王闯; 孔博; 张威; 葛欣星; | ||||
摘要 | 本 发明 公开一种基于星载可移点波束天线的地面 辐射 源 定位 方法,包括:(10)根据天线方向图函数及天线增益与 俯仰 角 和方位角之间的联系,构建天线方向图指向模型;(20)通过全球波束对该地面辐射源 信号 参数认知,通过可移点波束扫描得到该地面辐射源预估 位置 ;(30)根据预估计的辐射源位置,多次测量辐射信号强度,得到辐射信号强度的空间分布;(40)根据天线方向图指向模型和辐射信号强度空间分布,建立定向方程组,获得辐射源位置指向;(50)根据辐射源位置指向,结合地理信息与星历信息,得到地面辐射源的经纬度。本发明射源定位方法,简单易行,无需额外星载设备。 | ||||||
权利要求 | 1.一种基于星载可移点波束天线的地面辐射源定位方法,其特征在于,包括以下步骤: |
||||||
说明书全文 | 基于星载可移点波束天线的地面辐射源定位方法技术领域背景技术[0002] 由于卫星系统具有覆盖能力强、传输距离长、不受地理条件限制、机动灵活等独特的优点,卫星系统开始逐渐进入民用和国防领域,并在几十年中得到了迅速的发展,具有非常好的发展前景。然而,随着卫星通信的发展,卫星天线口径大、接收灵敏度高,卫星又处于开放的空间,GEO卫星通信系统更容易受到各种辐射源有意或无意的射频干扰。随着各国卫星业务的迅速增长,轨道与频率资源日趋紧张,卫星所处电磁环境日益恶化,卫星受到不明辐射源干扰事件的数量与日俱增,造成了诸多恶劣影响以及经济上的损失。为排除或降低不明辐射源干扰事件对卫星通信的影响,必须及时快速定位辐射源。 [0003] 常用的卫星辐射源定位方法有测时差定位和测向定位。测时差定位是随着多平台通信技术的发展和到达时间差技术的进步而发展起来的,并成为现代高精度辐射源定位的主要方法之一,主要有三星时差定位法、双星时差/频差联合定位法。一般需要由两颗或两颗以上的卫星组网同步相互配合实现。由于卫星轨道资源的有限性,对于GEO卫星而言,系统中仅包括有限的几颗卫星,很多情况下只有一颗卫星,很难找到协助定位的邻星,难以基于双星或三星干扰源定位体制实现对辐射源的定位。 [0004] 目前的单星侧向定位方法,主要包括:相位干涉仪测向、比相测向、比幅测向、时差测向、空间谱估计测向等。但这些单星侧向定位方法都必须额外增加星载定向设备,并且对卫星的姿态测量精度有较高的要求。 [0005] 总之,现有卫星辐射源定位技术存在的问题是:单星辐射源定位方法需要额外增加星载设备,导致发射成本高、星上资源消耗多,不易实现。。 发明内容[0006] 本发明的目的在于提供一种基于星载可移点波束天线的地面辐射源单星定位方法,利用星载可移点波束获取辐射源的信号强度空间分布,以解方程组的形式获取辐射源位置,方法简单,容易实现,无需额外增加星载设备。 [0007] 实现本发明目的的技术解决方案为: [0008] 一种基于星载可移点波束天线的地面辐射源定位方法,包括以下步骤: [0010] (20)预估辐射源位置:将星载全球波束天线接收到的地面辐射源信号进行参数认知,利用可移点波束在全球波束覆盖范围内交叉扫描,根据可移点波束接收到的信号强弱变化预估该地面辐射位置; [0011] (30)获取辐射源信号强度:根据预估计的辐射源位置,将可移点波束指向移至辐射源位置附近,并多次移动可移点波束指向,多次测量辐射信号强度,选取接收信号强度较强且中心指向位置成等边三角形的波束,结合指向信息,得到辐射信号强度的空间分布; [0012] (40)求解辐射源位置指向:根据天线方向图指向模型和辐射信号强度空间分布,建立定向方程组,求解定向方程组,获得辐射源位置指向; [0013] (50)获取辐射源经纬度:根据辐射源位置指向,结合地理信息与星历信息,得到地面辐射源的经纬度。 [0014] 本发明与现有技术相比,其显著优点: [0015] 1、定位方法易于实现:由于定位方程组形式简单,所需的数据容易获得,使得该定位方法容易实现; [0016] 2、发射成本低、星上资源消耗少:由于无需额外的星载设备用于辐射源定位,节约了卫星发射成本,降低了星上资源消耗。 [0017] 3、利用可移点波束实现辐射源定位,灵活性强。 附图说明[0019] 图1是本发明基于星载可移点波束天线的地面辐射源单星定位方法所涉及的流程图。 [0020] 图2是本发明基于星载可移点波束天线的地面辐射源单星定位方法的原理图。 [0021] 图3是定义α角和β角及与俯仰角 和方位角θ的转换示意图。 [0022] 图4是本发明方法中天线方向图指向模型所得指向与波束增益的映射图。 [0024] 图6是本发明方法在天线方向图为理想函数时求辐射源位置指向仿真结果图。 [0025] 图7是本发明方法根据数据库拟合天线方向图函数的仿真结果图。 [0026] 图8是本发明方法拟合天线方向图函数相对于理想天线方向图函数绝对误差图。 [0027] 图9是表1不同拟合粒度的拟合方向图函数误差情况表。 [0028] 图10是表2本发明方法根据拟合函数求解位置指向结果表。 [0029] 图11是表3本发明方法根据任意位置指向求经纬度结果表。 具体实施方式[0030] 如图1、2所示,本发明基于星载可移点波束天线的地面辐射源定位方法,其包括以下步骤: [0031] (10)构建方向图指向模型:根据天线方向图函数及天线增益与俯仰角和方位角之间的联系,构建天线方向图指向模型; [0032] 所述(10)构建方向图指向模型步骤包括: [0033] (11)获取天线方向图函数:若天线增益与指向关系为数据库形式储存,则根据数据库对天线方向图进行拟合,得到星载可移点波束天线方向图函数;否则,根据可移点波束天线的设计特性及参数,获取星载可移点波束天线方向图函数; [0034] 所述(11)获取天线方向图函数步骤包括: [0035] (111)天线方向图拟合:若天线增益与指向关系为数据库形式储存,则,[0036] 拟合函数 [0037] [0038] 其中, 表示一个拟合区域上的近似函数,b为基函数, [0039] bT=(1,x,x2), [0040] A(x)=BTW(x)B, [0041] D(x)=BTW(x), [0042] fT=(f1,f2,…,fi), [0043] 其中, [0044] [0045] [0047] (112)方向图函数获取:针对可移点波束,假设波束中心增益为G0,那么波束增益G近似计算为, [0048] [0049] 其中,u=2.07123sin(φ)/sin(φ3dB),J1和J3分别是1阶和3阶第一类贝塞尔函数,[0050] [0051] 其中,φ3dB是天线增益相对波束中心3dB衰减处所对应的角度,φ表示辐射源信号入射方向与可移点波束某一指向中心之间的夹角;D为天线口径,λ为辐射信号的波长,η为天线效率。 [0052] (12)构建方向图指向模型:根据可移点波束天线方向图函数,结合可移点波束中心指向,将天线增益与俯仰角和方位角相关联,得到反映指向与天线增益之间函数关系的方向图指向模型。 [0053] 所述(12)构建方向图指向模型步骤具体为: [0054] [0055] 其中, [0056] [0057] [0058] [0059] 为俯仰角,θ为方位角,α为任意指向位置在平行于卫星运行方向上与星下点的夹角,β为任意指向位置在垂直于卫星运行方向上与星下点的夹角,α0和β0为可移点波束某一中心指向。 [0060] (20)预估辐射源位置:将星载全球波束天线接收到的地面辐射源信号进行参数认知,利用可移点波束在全球波束覆盖范围内交叉扫描,根据可移点波束接收到的信号强弱变化预估该地面辐射位置; [0061] (30)获取辐射源信号强度:根据预估计的辐射源位置,将可移点波束指向移至辐射源位置附近,并多次移动可移点波束指向,多次测量辐射信号强度,选取接收信号强度较强且中心指向位置成等边三角形的波束,结合指向信息,得到辐射信号强度的空间分布; [0062] 所述(30)获取辐射源信号强度步骤包括: [0063] (31)测量辐射信号强度:将可移点波束移至地面辐射源预估位置附近,结合主波束认知的该地面辐射源信号参数信息,多次测量该地面辐射源的辐射信号强度; [0064] (32)获得辐射信号强度空间分布:结合接收信号强度、且中心指向位置成等边三角形的波束的指向信息,得到辐射信号强度在可移点波束天线环境下的空间分布。 [0065] (40)求解辐射源位置指向:根据天线方向图指向模型和辐射信号强度空间分布,建立定向方程组,求解定向方程组,获得辐射源位置指向; [0066] 所述(40)求解辐射源位置指向步骤包括: [0067] (41)建立定向方程组:利用天线方向图指向模型,结合辐射信号强度空间分布,建立由辐射源位置指向、可移点波束各指向接收辐射信号强度和天线方向图函数组成的定向方程组; [0068] 所述(41)建立定向方程组步骤具体为: [0069] 假设地面辐射源辐射强度为t,Lf为自由空间传输损耗,则可移点波束指向A接收的辐射信号强度uA表示为: [0070] t-Lf+G(φA)=uA, [0071] 由于φA为关于α和β的函数,将G(φA)记为GA(α,β),得: [0072] t-Lf+GA(α,β)=uA, [0073] 同理可得关于波束指向B和C时的链路计算方程: [0074] [0075] 其中,Gi(α,β)(i=A,B,C)为辐射源对应可移点波束三个指向获得的增益,uA、uB和uC为可移点波束接收辐射信号的强度值; [0076] 联立上述式子,化简得: [0077] [0078] 记: [0079] [0080] 求解辐射源位置指向α,β,相当于求解f=|f1|+|f2|取最小值时的α,β。 [0081] (42)获取辐射源位置指向:运用粒子群搜索算法,对定向方程组搜索求解,得到辐射源位置指向。 [0082] 如图5所示,所述(42)获取辐射源位置指向步骤包括: [0083] (421)初始化:在问题空间中随机产生粒子的位置与速度,位置即(α,β); [0084] (422)评价种群:计算每个粒子的适应度函数值,即粒子(α,β)对应的f的值; [0085] (423)更新最优:1)比较粒子适应度函数值与它的个体最优值Pbest,如果优于Pbest,则将其Pbest位置设置为当前粒子位置;2)比较粒子适应度函数值与群体最优值Gbest,如果优于Gbest,则设置Gbest位置为当前粒子位置; [0086] (424)更新粒子:按照式 和式 改变粒子的速度和位置; [0087] (425)停止条件:循环回到步骤2,直到满足最大的迭代次数,得到最优适应度函 数值fbest和最优粒子位置(α,β)best, [0088] 其中,(α,β)为辐射源位置指向,个体极值Pbest是指个体粒子搜索到的适应度值最优位置,群体极值Gbest是指种群中的所有粒子搜索到的适应度最优位置,假设在一个D维的搜索空间中,由n个粒子组成的种群X=(x1,x2…xn),其中第i个粒子表示为一个D维量Xi=(xi1,xi2…xiD)T,代表第i个粒子在D维搜索空间中的位置,亦代表问题的一个潜在解,根据目标函数即可计算出每个粒子位置Xi对应的适应度值,第i个粒子的速度为Vi=(Vi1,Vi2,…,ViD)T,其个体极值为Pi=(Pi1,Pi2,…,PiD)T,种群的群体极值为Pg=(Pg1,Pg2,…,PgD)T, [0089] 在每次迭代过程中,粒子通过个体极值和群体极值更新自身的速度和位置,即[0090] [0091] [0092] 其中,ω是惯性权重;d=1,2,…,D;i=1,2,…,n;k为当前迭代次数;Vid为粒子的速度;c1和c2是非负的常数,称为加速因子;r1和r2是分布在[0,1]区间的随机数。 [0093] 通过以上搜索算法步骤,对f=|f1|+|f2|进行求解,得到辐射源位置指向(α,β)。 [0094] (50)获取辐射源经纬度:根据辐射源位置指向,结合地理信息与星历信息,得到地面辐射源的经纬度。 [0095] 所述(50)获取辐射源经纬度步骤包括: [0096] (51)α角和β角与俯仰角 和方位角θ的转换: [0097] 根据几何模型, [0098] 俯仰角 [0099] 方位角 [0100] (52)求任意指向的经纬度:根据卫星的轨道运行参数,将卫星与地心地固坐标相关联;根据求得的辐射源位置指向的俯仰角和方位角,结合地面信息和卫星信息,得到辐射源指向点在地面的坐标,进而求出其经纬度。 [0101] 图3是定义α角和β角及与俯仰角 和方位角θ的转换示意图。由图中的集合模型,根据三角函数关系,可得: [0102] 俯仰角 [0103] 方位角 [0104] 其中,S为卫星所在点,O为星下点,A为指向位置,B为A在X轴上的投影。 [0105] 图4为本发明方法中天线方向图指向模型所得指向与波束增益的映射图,其中可移点波束中心指向的α角和β角假设为(5,0)。 [0106] 图6为用本发明方法对天线方向图为确定的理想函数,可移点波束中心指向的α角和β角分别为(5.0000,0.0000)、(5.3631,1.8871)、(6.8159,0.6289),辐射源真实位置指向的α角和β角为(6.0000,1.0000)的条件下进行定向仿真,运用粒子群搜索算法进行搜索求解,仿真次数100次取平均值的最优个体fbest与进化代数关系变化曲线及求解结果。数据表明,在理想情况下,本发明方法的定向精度非常高。 [0107] 图7为用本发明运用移动最小二乘法对数据库形式的天线方向图拟合结果,其中可移点波束中心指向的α角和β角设为(5.0000,0.0000),拟合数据库来源于理想方向图函数,数据库的粒度为指向的0.1°。 [0108] 图8为拟合函数相对于理想函数的绝对误差(dB),结果表明在天线方向图函数的主瓣边缘滚降较快的地方,误差会相对较大,但在与波束中心夹角小于2°的拟合函数误差非常小,在0.2dB以内。 [0109] 图9为检验拟合函数误差情况。将拟合的两种误差量定义为: [0110] [0111] [0112] 其中,f(xi)为原函数, 为拟合函数。由仿真结果表明:拟合函数的误差取决于数据库的粒度,数据库的粒度越小,拟合函数的误差越小,精度越高;数据库的粒度越大,拟合函数的误差越大,精度越低。 [0113] 图10为天线方向图为数据库形式时的定向结果,其中拟合方向图的拟合粒度为0.1°。数据表明误差基本在0.015°以内,转换为地面误差,距离约在10km以内,定向误差主要来源于拟合的天线方向图,若数据库的粒度更小,拟合精度会更高,定向精度也会更高。 [0114] 图11为本发明方法根据位置指向求经纬度结果,仿真条件为:以地球同步轨道卫星为例,将星下点的经纬度设为99.85°E、0°N,卫星轨道半长轴42164.137km、偏心率0°、轨道倾角0°。 [0115] 由于定位方程组形式简单,所需的数据容易获得,使得本发明的单星定位方法容易实现;而且,由于无需额外的星载设备用于辐射源定位,节约了卫星发射成本,降低了星上资源消耗;利用可移点波束实现辐射源定位,弥补固定点波束的不足,灵活性强。 |