CN105676166B - 基于星载可移点波束天线的地面辐射源定位方法 - Google Patents
基于星载可移点波束天线的地面辐射源定位方法 Download PDFInfo
- Publication number
- CN105676166B CN105676166B CN201610034916.8A CN201610034916A CN105676166B CN 105676166 B CN105676166 B CN 105676166B CN 201610034916 A CN201610034916 A CN 201610034916A CN 105676166 B CN105676166 B CN 105676166B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- radiation source
- antenna
- pointing
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 40
- 230000004807 localization Effects 0.000 title claims abstract 4
- 230000005855 radiation Effects 0.000 claims abstract description 139
- 230000006870 function Effects 0.000 claims description 67
- 239000002245 particle Substances 0.000 claims description 46
- 238000010586 diagram Methods 0.000 claims description 21
- 238000010845 search algorithm Methods 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 4
- 230000005540 biological transmission Effects 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000019771 cognition Effects 0.000 claims description 3
- 230000001133 acceleration Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 claims description 2
- 238000013461 design Methods 0.000 claims description 2
- 230000001149 cognitive effect Effects 0.000 abstract 1
- 235000019580 granularity Nutrition 0.000 description 7
- 238000004088 simulation Methods 0.000 description 6
- 238000004891 communication Methods 0.000 description 5
- 238000005259 measurement Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000013507 mapping Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S1/00—Beacons or beacon systems transmitting signals having a characteristic or characteristics capable of being detected by non-directional receivers and defining directions, positions, or position lines fixed relatively to the beacon transmitters; Receivers co-operating therewith
- G01S1/02—Beacons or beacon systems transmitting signals having a characteristic or characteristics capable of being detected by non-directional receivers and defining directions, positions, or position lines fixed relatively to the beacon transmitters; Receivers co-operating therewith using radio waves
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Geophysics And Detection Of Objects (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明公开一种基于星载可移点波束天线的地面辐射源定位方法，包括：(10)根据天线方向图函数及天线增益与俯仰角和方位角之间的联系，构建天线方向图指向模型；(20)通过全球波束对该地面辐射源信号参数认知，通过可移点波束扫描得到该地面辐射源预估位置；(30)根据预估计的辐射源位置，多次测量辐射信号强度，得到辐射信号强度的空间分布；(40)根据天线方向图指向模型和辐射信号强度空间分布，建立定向方程组，获得辐射源位置指向；(50)根据辐射源位置指向，结合地理信息与星历信息，得到地面辐射源的经纬度。本发明射源定位方法，简单易行，无需额外星载设备。
Description
技术领域
本发明属于无线通信技术领域，特别是一种发射成本低、星上资源消耗少、定位精度较高的基于星载可移点波束天线的地面辐射源定位方法。
背景技术
由于卫星系统具有覆盖能力强、传输距离长、不受地理条件限制、机动灵活等独特的优点，卫星系统开始逐渐进入民用和国防领域，并在几十年中得到了迅速的发展，具有非常好的发展前景。然而，随着卫星通信的发展，卫星天线口径大、接收灵敏度高，卫星又处于开放的空间，GEO卫星通信系统更容易受到各种辐射源有意或无意的射频干扰。随着各国卫星业务的迅速增长，轨道与频率资源日趋紧张，卫星所处电磁环境日益恶化，卫星受到不明辐射源干扰事件的数量与日俱增，造成了诸多恶劣影响以及经济上的损失。为排除或降低不明辐射源干扰事件对卫星通信的影响，必须及时快速定位辐射源。
常用的卫星辐射源定位方法有测时差定位和测向定位。测时差定位是随着多平台通信技术的发展和到达时间差技术的进步而发展起来的，并成为现代高精度辐射源定位的主要方法之一，主要有三星时差定位法、双星时差/频差联合定位法。一般需要由两颗或两颗以上的卫星组网同步相互配合实现。由于卫星轨道资源的有限性，对于GEO卫星而言，系统中仅包括有限的几颗卫星，很多情况下只有一颗卫星，很难找到协助定位的邻星，难以基于双星或三星干扰源定位体制实现对辐射源的定位。
目前的单星侧向定位方法，主要包括：相位干涉仪测向、比相测向、比幅测向、时差测向、空间谱估计测向等。但这些单星侧向定位方法都必须额外增加星载定向设备，并且对卫星的姿态测量精度有较高的要求。
总之，现有卫星辐射源定位技术存在的问题是：单星辐射源定位方法需要额外增加星载设备，导致发射成本高、星上资源消耗多，不易实现。。
发明内容
本发明的目的在于提供一种基于星载可移点波束天线的地面辐射源定位方法，利用星载可移点波束获取辐射源的信号强度空间分布，以解方程组的形式获取辐射源位置，方法简单，容易实现，无需额外增加星载设备。
实现本发明目的的技术解决方案为：
一种基于星载可移点波束天线的地面辐射源定位方法，包括以下步骤：
(10)构建方向图指向模型：根据天线方向图函数及天线增益与俯仰角和方位角之间的联系，构建天线方向图指向模型；
(20)预估辐射源位置：将星载全球波束天线接收到的地面辐射源信号进行参数认知，利用可移点波束在全球波束覆盖范围内交叉扫描，根据可移点波束接收到的信号强弱变化预估该地面辐射位置；
(30)获取辐射源信号强度：根据预估计的辐射源位置，将可移点波束指向移至辐射源位置附近，并多次移动可移点波束指向，多次测量辐射信号强度，选取接收信号强度较强且中心指向位置成等边三角形的波束，结合指向信息，得到辐射信号强度的空间分布；
(40)求解辐射源位置指向：根据天线方向图指向模型和辐射信号强度空间分布，建立定向方程组，求解定向方程组，获得辐射源位置指向；
(50)获取辐射源经纬度：根据辐射源位置指向，结合地理信息与星历信息，得到地面辐射源的经纬度。
本发明与现有技术相比，其显著优点：
1、定位方法易于实现：由于定位方程组形式简单，所需的数据容易获得，使得该定位方法容易实现；
2、发射成本低、星上资源消耗少：由于无需额外的星载设备用于辐射源定位，节约了卫星发射成本，降低了星上资源消耗。
3、利用可移点波束实现辐射源定位，灵活性强。
下面结合附图和具体实施方式对本发明作进一步的详细描述。
附图说明
图1是本发明基于星载可移点波束天线的地面辐射源定位方法所涉及的流程图。
图2是本发明基于星载可移点波束天线的地面辐射源定位方法的原理图。
图3是定义α角和β角及与俯仰角和方位角θ的转换示意图。
图4是本发明方法中天线方向图指向模型所得指向与波束增益的映射图。
图5是求解辐射源位置指向的粒子群搜索算法程序流程图。
图6是本发明方法在天线方向图为理想函数时求辐射源位置指向仿真结果图。
图7是本发明方法根据数据库拟合天线方向图函数的仿真结果图。
图8是本发明方法拟合天线方向图函数相对于理想天线方向图函数绝对误差图。
图9是表1不同拟合粒度的拟合方向图函数误差情况表。
图10是表2本发明方法根据拟合函数求解位置指向结果表。
图11是表3本发明方法根据任意位置指向求经纬度结果表。
具体实施方式
如图1、2所示，本发明基于星载可移点波束天线的地面辐射源定位方法，其包括以下步骤：
(10)构建方向图指向模型：根据天线方向图函数及天线增益与俯仰角和方位角之间的联系，构建天线方向图指向模型；
所述(10)构建方向图指向模型步骤包括：
(11)获取天线方向图函数：若天线增益与指向关系为数据库形式储存，则根据数据库对天线方向图进行拟合，得到星载可移点波束天线方向图函数；否则，根据可移点波束天线的设计特性及参数，获取星载可移点波束天线方向图函数；
所述(11)获取天线方向图函数步骤包括：
(111)天线方向图拟合：若天线增益与指向关系为数据库形式储存，则，
拟合函数
其中，表示一个拟合区域上的近似函数，b为基函数，
bT＝(1,x,x2)，
A(x)＝BTW(x)B，
D(x)＝BTW(x)，
fT＝(f1，f2，…，fi)，
其中，
n是x近影响节点数；m是基函数项数；w(x-xi)是节点xi具有紧支撑性质的光滑连续权函数，在xi紧支撑域内部和边界的wi＝w(x-xi)＞0,在其外部wi＝0；fi是xi节点值；
(112)方向图函数获取：针对可移点波束，假设波束中心增益为G0，那么波束增益G近似计算为，
其中，u＝2.07123sin(φ)/sin(φ3dB)，J1(u)和J3(u)分别是1阶和3阶第一类贝塞尔函数，
其中，φ3dB是天线增益相对波束中心3dB衰减处所对应的角度，φ表示辐射源信号入射方向与可移点波束某一指向中心之间的夹角；D为天线口径，λ为辐射信号的波长，η为天线效率。
(12)构建方向图指向模型：根据可移点波束天线方向图函数，结合可移点波束中心指向，将天线增益与俯仰角和方位角相关联，得到反映指向与天线增益之间函数关系的方向图指向模型。
所述(12)构建方向图指向模型步骤具体为：
其中，
为俯仰角，θ为方位角，α为任意指向位置在平行于卫星运行方向上与星下点的夹角，β为任意指向位置在垂直于卫星运行方向上与星下点的夹角，(α0，β0)为可移点波束某一中心指向。
(20)预估辐射源位置：将星载全球波束天线接收到的地面辐射源信号进行参数认知，利用可移点波束在全球波束覆盖范围内交叉扫描，根据可移点波束接收到的信号强弱变化预估该地面辐射位置；
(30)获取辐射源信号强度：根据预估计的辐射源位置，将可移点波束指向移至辐射源位置附近，并多次移动可移点波束指向，多次测量辐射信号强度，选取接收信号强度较强且中心指向位置成等边三角形的波束，结合指向信息，得到辐射信号强度的空间分布；
所述(30)获取辐射源信号强度步骤包括：
(31)测量辐射信号强度：将可移点波束移至地面辐射源预估位置附近，结合主波束认知的该地面辐射源信号参数信息，多次测量该地面辐射源的辐射信号强度；
(32)获得辐射信号强度空间分布：结合接收信号强度、且中心指向位置成等边三角形的波束的指向信息，得到辐射信号强度在可移点波束天线环境下的空间分布。
(40)求解辐射源位置指向：根据天线方向图指向模型和辐射信号强度空间分布，建立定向方程组，求解定向方程组，获得辐射源位置指向；
所述(40)求解辐射源位置指向步骤包括：
(41)建立定向方程组：利用天线方向图指向模型，结合辐射信号强度空间分布，建立由辐射源位置指向、可移点波束各指向接收辐射信号强度和天线方向图函数组成的定向方程组；
所述(41)建立定向方程组步骤具体为：
假设地面辐射源辐射强度为t，Lf为自由空间传输损耗，则可移点波束指向A接收的辐射信号强度uA表示为：
t-Lf+G(φA)＝uA，
由于φA为关于α和β的函数，将G(φA)记为GA(α,β)，得：
t-Lf+GA(α,β)＝uA，
同理可得关于波束指向B和C时的链路计算方程：
其中，Gi(α,β)(i＝A,B,C)为辐射源对应可移点波束三个指向获得的增益，uA、uB和uC为可移点波束接收辐射信号的强度值；
联立上述式子，化简得：
记：
求解辐射源位置指向(α,β)，相当于求解f＝|f1|+|f2|取最小值时的(α,β)。
(42)获取辐射源位置指向：运用粒子群搜索算法，对定向方程组搜索求解，得到辐射源位置指向。
如图5所示，所述(42)获取辐射源位置指向步骤包括：
(421)初始化：在问题空间中随机产生粒子的位置与速度，位置即(α,β)；
(422)评价种群：计算每个粒子的适应度函数值，即粒子(α,β)对应的f的值；
(423)更新最优：1)比较粒子适应度函数值与它的个体最优值Pbest，如果优于Pbest，则将其Pbest位置设置为当前粒子位置；2)比较粒子适应度函数值与群体最优值Gbest，如果优于Gbest，则设置Gbest位置为当前粒子位置；
(424)更新粒子：按照式和式改变粒子的速度和位置；
(425)停止条件：循环回到(422)评价种群步骤，直到满足最大的迭代次数，得到最优适应度函数值fbest和最优粒子位置(α,β)best，
其中，(α,β)为辐射源位置指向，个体极值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，
在每次迭代过程中，粒子通过个体极值和群体极值更新自身的速度和位置，即
其中，ω是惯性权重；d＝1,2，…,D；i＝1,2,…,n；k为当前迭代次数；Vid为粒子的速度；c1和c2是非负的常数，称为加速因子；r1和r2是分布在[0,1]区间的随机数。
通过以上搜索算法步骤，对f＝|f1|+|f2|进行求解，得到辐射源位置指向(α,β)。
(50)获取辐射源经纬度：根据辐射源位置指向，结合地理信息与星历信息，得到地面辐射源的经纬度。
所述(50)获取辐射源经纬度步骤包括：
(51)α角和β角与俯仰角和方位角θ的转换：
根据几何模型，
俯仰角
方位角
(52)求任意指向的经纬度：根据卫星的轨道运行参数，将卫星与地心地固坐标相关联；根据求得的辐射源位置指向的俯仰角和方位角，结合地面信息和卫星信息，得到辐射源指向点在地面的坐标，进而求出其经纬度。
图3是定义α角和β角及与俯仰角和方位角θ的转换示意图。由图中的集合模型，根据三角函数关系，可得：
俯仰角
方位角
其中，S为卫星所在点，O为星下点，A为指向位置，B为A在X轴上的投影。
图4为本发明方法中天线方向图指向模型所得指向与波束增益的映射图，其中可移点波束中心指向的α角和β角假设为(5,0)。
图6为用本发明方法对天线方向图为确定的理想函数，可移点波束中心指向的α角和β角分别为(5.0000,0.0000)、(5.3631,1.8871)、(6.8159,0.6289)，辐射源真实位置指向的α角和β角为(6.0000,1.0000)的条件下进行定向仿真，运用粒子群搜索算法进行搜索求解，仿真次数100次取平均值的最优个体fbest与进化代数关系变化曲线及求解结果。数据表明，在理想情况下，本发明方法的定向精度非常高。
图7为用本发明运用移动最小二乘法对数据库形式的天线方向图拟合结果，其中可移点波束中心指向的α角和β角设为(5.0000,0.0000)，拟合数据库来源于理想方向图函数，数据库的粒度为指向的0.1°。
图8为拟合函数相对于理想函数的绝对误差(dB)，结果表明在天线方向图函数的主瓣边缘滚降较快的地方，误差会相对较大，但在与波束中心夹角小于2°的拟合函数误差非常小，在0.2dB以内。
图9为检验拟合函数误差情况。将拟合的两种误差量定义为：
其中，f(xi)为原函数，为拟合函数。由仿真结果表明：拟合函数的误差取决于数据库的粒度，数据库的粒度越小，拟合函数的误差越小，精度越高；数据库的粒度越大，拟合函数的误差越大，精度越低。
图10为天线方向图为数据库形式时的定向结果，其中拟合方向图的拟合粒度为0.1°。数据表明误差基本在0.015°以内，转换为地面误差，距离约在10km以内，定向误差主要来源于拟合的天线方向图，若数据库的粒度更小，拟合精度会更高，定向精度也会更高。
图11为本发明方法根据位置指向求经纬度结果，仿真条件为：以地球同步轨道卫星为例，将星下点的经纬度设为99.85°E、0°N，卫星轨道半长轴42164.137km、偏心率0°、轨道倾角0°。
由于定位方程组形式简单，所需的数据容易获得，使得本发明的单星定位方法容易实现；而且，由于无需额外的星载设备用于辐射源定位，节约了卫星发射成本，降低了星上资源消耗；利用可移点波束实现辐射源定位，弥补固定点波束的不足，灵活性强。
Claims (9)
1.一种基于星载可移点波束天线的地面辐射源定位方法，其特征在于，包括以下步骤：
(10)构建方向图指向模型：根据天线方向图函数及天线增益与俯仰角和方位角之间的联系，构建天线方向图指向模型；
(20)预估辐射源位置：将星载全球波束天线接收到的地面辐射源信号进行参数认知，利用可移点波束在全球波束覆盖范围内交叉扫描，根据可移点波束接收到的信号强弱变化预估该地面辐射位置；
(30)获取辐射源信号强度：根据预估计的辐射源位置，将可移点波束指向移至辐射源位置附近，并多次移动可移点波束指向，多次测量辐射信号强度，选取接收信号强度较强且中心指向位置成等边三角形的波束，结合指向信息，得到辐射信号强度的空间分布；
(40)求解辐射源位置指向：根据天线方向图指向模型和辐射信号强度空间分布，建立定向方程组，求解定向方程组，获得辐射源位置指向；
(50)获取辐射源经纬度：根据辐射源位置指向，结合地理信息与星历信息，得到地面辐射源的经纬度。
2.权利要求1所述的地面辐射源定位方法，其特征在于，所述(10)构建方向图指向模型步骤包括：
(11)获取天线方向图函数：若天线增益与指向关系为数据库形式储存，则根据数据库对天线方向图进行拟合，得到星载可移点波束天线方向图函数；否则，根据可移点波束天线的设计特性及参数，获取星载可移点波束天线各波束天线方向图函数；
(12)构建方向图指向模型：根据可移点波束天线方向图函数，结合可移点波束中心指向，将天线增益与俯仰角和方位角相关联，得到反映指向与天线增益之间函数关系的方向图指向模型。
3.权利要求2所述的地面辐射源定位方法，其特征在于，所述(11)获取天线方向图函数步骤包括：
(111)天线方向图拟合：若天线增益与指向关系为数据库形式储存，则，
拟合函数
<mrow>
<mover>
<mi>F</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mi>b</mi>
<mi>T</mi>
</msup>
<msup>
<mi>A</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mi>D</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mi>f</mi>
<mo>,</mo>
</mrow>
其中，表示一个拟合区域上的近似函数，b为基函数，
bT＝(1,x,x2)，
A(x)＝BTW(x)B，
D(x)＝BTW(x)，
fT＝(f1，f2，…，fi)，
其中，
n是x近影响节点数；m是基函数项数；w(x-xi)是节点xi具有紧支撑性质的光滑连续权函数，在xi紧支撑域内部和边界的wi＝w(x-xi)＞0,在其外部wi＝0；fi是xi节点值；
(112)方向图函数获取：针对可移点波束，假设波束中心增益为G0，那么波束增益G近似计算为，
<mrow>
<mi>G</mi>
<mo>=</mo>
<msub>
<mi>G</mi>
<mn>0</mn>
</msub>
<mo>*</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mfrac>
<mrow>
<msub>
<mi>J</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>2</mn>
<mi>u</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mn>36</mn>
<mfrac>
<mrow>
<msub>
<mi>J</mi>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
</mrow>
<msup>
<mi>u</mi>
<mn>3</mn>
</msup>
</mfrac>
<mo>&rsqb;</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>,</mo>
</mrow>
其中，u＝2.07123sin(φ)/sin(φ3dB)，J1(u)和J3(u)分别是1阶和3阶第一类贝塞尔函数，
<mrow>
<msub>
<mi>G</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<mi>&pi;</mi>
<mn>2</mn>
</msup>
<msup>
<mi>D</mi>
<mn>2</mn>
</msup>
<mi>&eta;</mi>
</mrow>
<msup>
<mi>&lambda;</mi>
<mn>2</mn>
</msup>
</mfrac>
<mo>,</mo>
</mrow>
其中，φ3dB是天线增益相对波束中心3dB衰减处所对应的角度，φ表示辐射源信号入射方向与可移点波束某一指向中心之间的夹角；D为天线口径，λ为辐射信号的波长，η为天线效率，波束中心增益为G0、波束增益G的单位为dBi。
4.权利要求3所述的地面辐射源定位方法，其特征在于，所述(12)构建方向图指向模型步骤具体为：
其中，
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&theta;</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mi>t</mi>
<mi>a</mi>
<mi>n</mi>
<mi>&alpha;</mi>
</mrow>
<msqrt>
<mrow>
<msup>
<mi>tan</mi>
<mn>2</mn>
</msup>
<mi>&alpha;</mi>
<mo>+</mo>
<msup>
<mi>tan</mi>
<mn>2</mn>
</msup>
<mi>&beta;</mi>
</mrow>
</msqrt>
</mfrac>
<mo>,</mo>
</mrow>
为俯仰角，θ为方位角，α为任意指向位置在平行于卫星运行方向上与星下点的夹角，β为任意指向位置在垂直于卫星运行方向上与星下点的夹角，(α0，β0)为可移点波束某一中心指向。
5.权利要求1所述的地面辐射源定位方法，其特征在于，所述(30)获取辐射源信号强度步骤包括：
(31)测量辐射信号强度：将可移点波束移至地面辐射源预估位置附近，结合主波束认知的该地面辐射源信号参数信息，多次测量该地面辐射源的辐射信号强度；
(32)获得辐射信号强度空间分布：结合接收信号强度强、且中心指向位置成等边三角形的波束的指向信息，得到辐射信号强度在可移点波束天线环境下的空间分布。
6.权利要求4所述的地面辐射源定位方法，其特征在于，所述(40)求解辐射源位置指向步骤包括：
(41)建立定向方程组：利用天线方向图指向模型，结合辐射信号强度空间分布，建立由辐射源位置指向、可移点波束各指向接收辐射信号强度和天线方向图函数组成的定向方程组；
(42)获取辐射源位置指向：运用粒子群搜索算法，对定向方程组搜索求解，得到辐射源位置指向。
7.权利要求6所述的地面辐射源定位方法，其特征在于，所述(41)建立定向方程组步骤具体为：
假设地面辐射源辐射强度为t，Lf为自由空间传输损耗，则可移点波束指向A接收的辐射信号强度uA表示为：
t-Lf+G(φA)＝uA，
由于φA为关于α和β的函数，将G(φA)记为GA(α,β)，得：
t-Lf+GA(α,β)＝uA，
同理可得关于波束指向B和C时的链路计算方程：
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>t</mi>
<mo>-</mo>
<msub>
<mi>L</mi>
<mi>f</mi>
</msub>
<mo>+</mo>
<msub>
<mi>G</mi>
<mi>B</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>,</mo>
<mi>&beta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>u</mi>
<mi>B</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>t</mi>
<mo>-</mo>
<msub>
<mi>L</mi>
<mi>f</mi>
</msub>
<mo>+</mo>
<msub>
<mi>G</mi>
<mi>C</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>,</mo>
<mi>&beta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>u</mi>
<mi>C</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
其中，Gi(α,β)(i＝A,B,C)为辐射源对应可移点波束三个指向获得的增益，uA、uB和uC为可移点波束接收辐射信号的强度值；
联立上述式子，化简得：
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>G</mi>
<mi>A</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>,</mo>
<mi>&beta;</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>G</mi>
<mi>B</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>,</mo>
<mi>&beta;</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>u</mi>
<mi>A</mi>
</msub>
<mo>-</mo>
<msub>
<mi>u</mi>
<mi>B</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>G</mi>
<mi>A</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>,</mo>
<mi>&beta;</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>G</mi>
<mi>C</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>,</mo>
<mi>&beta;</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>u</mi>
<mi>A</mi>
</msub>
<mo>-</mo>
<msub>
<mi>u</mi>
<mi>C</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
记：
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<msub>
<mi>G</mi>
<mi>A</mi>
</msub>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>,</mo>
<mi>&beta;</mi>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>G</mi>
<mi>B</mi>
</msub>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>,</mo>
<mi>&beta;</mi>
<mo>)</mo>
<mo>-</mo>
<mo>(</mo>
<msub>
<mi>u</mi>
<mi>A</mi>
</msub>
<mo>-</mo>
<msub>
<mi>u</mi>
<mi>B</mi>
</msub>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<msub>
<mi>G</mi>
<mi>A</mi>
</msub>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>,</mo>
<mi>&beta;</mi>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>G</mi>
<mi>C</mi>
</msub>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>,</mo>
<mi>&beta;</mi>
<mo>)</mo>
<mo>-</mo>
<mo>(</mo>
<msub>
<mi>u</mi>
<mi>A</mi>
</msub>
<mo>-</mo>
<msub>
<mi>u</mi>
<mi>C</mi>
</msub>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
求解辐射源位置指向(α,β)，相当于求解f＝|f1|+|f2|取最小值时的(α,β)。
8.根据权利要求7所述的地面辐射源定位方法，其特征在于，所述(42)获取辐射源位置指向步骤包括：
(421)初始化：在问题空间中随机产生粒子的位置与速度，位置即(α,β)；
(422)评价种群：计算每个粒子的适应度函数值，即粒子(α,β)对应的f的值；
(423)更新最优：1)比较粒子适应度函数值与它的个体最优值Pbest，如果优于Pbest，则将其Pbest位置设置为当前粒子位置；2)比较粒子适应度函数值与群体最优值Gbest，如果优于Gbest，则设置Gbest位置为当前粒子位置；
(424)更新粒子：按照式和式改变粒子的速度和位置；
(425)停止条件：循环回到(422)评价种群步骤，直到满足最大的迭代次数，得到最优适应度函数值fbest和最优粒子位置(α,β)best，
其中，(α,β)为辐射源位置指向，个体极值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，
在每次迭代过程中，粒子通过个体极值和群体极值更新自身的速度和位置，即
<mrow>
<msubsup>
<mi>V</mi>
<mrow>
<mi>i</mi>
<mi>d</mi>
</mrow>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>&omega;V</mi>
<mrow>
<mi>i</mi>
<mi>d</mi>
</mrow>
<mi>k</mi>
</msubsup>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
<msub>
<mi>r</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mi>d</mi>
</mrow>
<mi>k</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>d</mi>
</mrow>
<mi>k</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>2</mn>
</msub>
<msub>
<mi>r</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>g</mi>
<mi>d</mi>
</mrow>
<mi>k</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>d</mi>
</mrow>
<mi>k</mi>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>d</mi>
</mrow>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>d</mi>
</mrow>
<mi>k</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>V</mi>
<mrow>
<mi>i</mi>
<mi>d</mi>
</mrow>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
</mrow>
其中，ω是惯性权重；d＝1,2，…,D；i＝1,2,…,n；k为当前迭代次数；Vid为粒子的速度；c1和c2是非负的常数，称为加速因子；r1和r2是分布在[0,1]区间的随机数。
9.根据权利要求8所述的地面辐射源定位方法，其特征在于，所述(50)获取辐射源经纬度步骤包括：
(51)α角和β角与俯仰角和方位角θ的转换：
根据几何模型，
俯仰角
方位角
(52)求任意指向的经纬度：根据卫星的轨道运行参数，将卫星与地心地固坐标相关联；根据求得的辐射源位置指向的俯仰角和方位角，结合地面信息和卫星信息，得到辐射源指向点在地面的坐标，进而求出其经纬度。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610034916.8A CN105676166B (zh) | 2016-01-19 | 2016-01-19 | 基于星载可移点波束天线的地面辐射源定位方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610034916.8A CN105676166B (zh) | 2016-01-19 | 2016-01-19 | 基于星载可移点波束天线的地面辐射源定位方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105676166A CN105676166A (zh) | 2016-06-15 |
| CN105676166B true CN105676166B (zh) | 2017-09-15 |
Family
ID=56301590
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610034916.8A Active CN105676166B (zh) | 2016-01-19 | 2016-01-19 | 基于星载可移点波束天线的地面辐射源定位方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105676166B (zh) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108549050B (zh) * | 2018-03-27 | 2022-03-04 | 天津工业大学 | 一种基于定向充电时间的可充电节点定位方法 |
| CN110531315B (zh) * | 2019-07-15 | 2020-08-21 | 南京航空航天大学 | 一种基于信号强度变化率的卫星干扰源直接定位方法及定位装置 |
| CN112346003B (zh) * | 2020-10-20 | 2023-12-29 | 西安空间无线电技术研究所 | 一种基于等势寻优的单波束测向系统 |
| CN112671453A (zh) * | 2020-12-18 | 2021-04-16 | 航天恒星科技有限公司 | 低轨星座的电磁态势建模方法和装置 |
| CN113189625B (zh) * | 2021-04-29 | 2024-05-14 | 中国电子科技集团公司第三十六研究所 | 基于单星干涉仪测向系统的误差校正方法、装置及卫星 |
| CN114035150B (zh) * | 2021-07-28 | 2023-07-18 | 中国电子科技集团公司第二十九研究所 | 一种基于无人机升空平台的无线电射频源测向装置及定位方法 |
| CN113884978A (zh) * | 2021-09-24 | 2022-01-04 | 正成卫星网络集团有限公司 | 一种victs天线波束指向测试方法 |
| CN114390432B (zh) * | 2022-01-14 | 2024-06-25 | 中国人民解放军61096部队 | 干扰源定位方法及装置 |
| CN115826004B (zh) * | 2022-10-24 | 2023-12-22 | 湖南六九零六信息科技股份有限公司 | 一种基于二维角度及时差联合的三星协同直接定位方法 |
| CN115986407B (zh) * | 2023-03-20 | 2023-05-26 | 四川省华盾防务科技股份有限公司 | 相控阵系统的控制方法、装置及存储介质 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5917446A (en) * | 1995-11-08 | 1999-06-29 | The Charles Stark Draper Laboratory, Inc. | Radio-wave reception system using inertial data in the receiver beamforming operation |
| CN101765785A (zh) * | 2007-05-21 | 2010-06-30 | 空间数码系统公司 | 用于多波束全球定位系统(gps)接收器的数字波束成形装置和技术 |
| CN103746757A (zh) * | 2014-01-08 | 2014-04-23 | 中国人民解放军理工大学 | 一种基于星载多波束天线的单星干扰源定位方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR100951618B1 (ko) * | 2008-02-19 | 2010-04-09 | 한국과학기술원 | 광주파수 발생기를 이용한 절대거리 측정방법 및 시스템 |
-
2016
- 2016-01-19 CN CN201610034916.8A patent/CN105676166B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5917446A (en) * | 1995-11-08 | 1999-06-29 | The Charles Stark Draper Laboratory, Inc. | Radio-wave reception system using inertial data in the receiver beamforming operation |
| CN101765785A (zh) * | 2007-05-21 | 2010-06-30 | 空间数码系统公司 | 用于多波束全球定位系统(gps)接收器的数字波束成形装置和技术 |
| CN103746757A (zh) * | 2014-01-08 | 2014-04-23 | 中国人民解放军理工大学 | 一种基于星载多波束天线的单星干扰源定位方法 |
Non-Patent Citations (2)
| Title |
|---|
| Particle swarm optimization with varying inertia weight for solving nonlinear optimization;PANDEY B B.el;《The 2015 International Conference on EESCO》;20150125;第1-5页 * |
| 一种基于星载多波束天线的单星干扰源定位方法;余强等;《通信技术》;20150131;第48卷(第1期);第37-40、45页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105676166A (zh) | 2016-06-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105676166B (zh) | 基于星载可移点波束天线的地面辐射源定位方法 | |
| CN105589086B (zh) | 基于信号强度多波束空间分布的地面辐射源单星定位方法 | |
| CN103746757B (zh) | 一种基于星载多波束天线的单星干扰源定位方法 | |
| Lazarev et al. | Positioning for location-aware beamforming in 5G ultra-dense networks | |
| GB2534894A (en) | Direction finding using signal power | |
| CN106840212A (zh) | 基于地面激光光斑质心位置的星载激光在轨几何检校方法 | |
| CN109975755B (zh) | 一种校正源存在条件下的短波多站直接定位方法 | |
| CN110531315B (zh) | 一种基于信号强度变化率的卫星干扰源直接定位方法及定位装置 | |
| Klus et al. | Neural network fingerprinting and GNSS data fusion for improved localization in 5G | |
| CN103323856A (zh) | 基于高轨三星时差体制的非合作无线电信号源定位方法 | |
| Aitmagambetov et al. | Energy budget and methods for determining coordinates for a radiomonitoring system based on a small spacecraft | |
| RU2613369C1 (ru) | Способ навигации летательного аппарата с использованием высокоточного одноэтапного пеленгатора и адресно-ответной пакетной цифровой радиолинии в дкмв диапазоне | |
| CN110300901A (zh) | 用于确定地球大气层中的电子密度分布的方法 | |
| CN102023301A (zh) | 一种中轨道卫星搜救系统的选星方法 | |
| CN117031519A (zh) | 一种多目标三维协同定位方法及系统 | |
| CN108828513B (zh) | 基于多监测点电波传播衰减等差线交叉的信号源定位方法 | |
| CN114155347B (zh) | 基于菲波那契网格的星载无源定位误差分布图投影方法 | |
| CN114993299A (zh) | 一种基于超宽带的无人机辅助定位方法及系统 | |
| CN111025358B (zh) | 基于导航卫星信号单短基线的定向方法 | |
| Li et al. | An effective integrated communication and localization method based on digital phased array antenna | |
| Dong et al. | Passive localization for frequency hopping signal emitter based on synthetic aperture principle | |
| Hao et al. | A new passive localization method of the interference source for satellite communications | |
| CN107272024B (zh) | 一种飞行器在高动态环境下的基准站布设方法 | |
| Kendra et al. | Motion-Extended Array Synthesis—Part II: Experimental Validation | |
| RU2606958C1 (ru) | Способ определения местоположения станции сети связи vsat |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |