CN106772466B - 一种基于形状特征搜索的近地卫星目标自动捕获算法 - Google Patents
一种基于形状特征搜索的近地卫星目标自动捕获算法 Download PDFInfo
- Publication number
- CN106772466B CN106772466B CN201611009936.6A CN201611009936A CN106772466B CN 106772466 B CN106772466 B CN 106772466B CN 201611009936 A CN201611009936 A CN 201611009936A CN 106772466 B CN106772466 B CN 106772466B
- Authority
- CN
- China
- Prior art keywords
- antenna
- point
- search
- curve
- target
- 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.)
- Expired - Fee Related
Links
- 238000010586 diagram Methods 0.000 claims description 18
- 230000005855 radiation Effects 0.000 claims description 13
- 230000003068 static effect Effects 0.000 claims description 4
- 238000000034 method Methods 0.000 abstract description 9
- 230000004069 differentiation Effects 0.000 abstract description 2
- 238000012795 verification Methods 0.000 abstract description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000010845 search algorithm 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
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/24—Acquisition or tracking or demodulation of signals transmitted by the system
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明提供了一种基于形状特征搜索的近地卫星目标自动捕获算法,利用建立天线的方向图模型,然后将目标与天线相对运动形成的电压的变化率曲线作为待检测曲线,在天线的方向图模型中搜索最符合该曲线特征的起始点和终止点,最后利用终止点确定目标与天线电轴的相对位置,从而控制天线完成目标的捕获跟踪。实际验证结果表明,当目标在天线副瓣范围内运动时,该方法能自动完成目标的捕获跟踪,而无需人工的扫描干预,同时该方法也可用于主副瓣判别。该方法能有效提高天线的搜索范围,缩短捕获时间,提高自动化成功率。
Description
技术领域
本发明属于在轨航天器测控管理领域,适用于动态目标的智能搜索和捕获。
背景技术
大部分近地卫星具有机动性强、轨道高度低、衰减快的特点,地面测控站为了获得更精确的轨道测量数据,常使用大口径的天线,但天线口径的增大也就意味着主瓣宽度的变窄,特别是当大动态目标或当航天器失去控制轨道变得不可准确预测时,目标的搜索和捕获成为难题。此时完成目标的捕获跟踪通常的做法是依靠操作人员对目标进行人工扫描搜索,这样更多依赖于操作员的经验积累,并且操作时长不确定,成功与否不确定,极大影响了正常测控流程的进行。
发明内容
为了克服现有技术的不足,本发明提供一种基于形状特征的捕获算法SASG(Search Algorithm on Sharp Graph),建立天线TE21跟踪模的方向图模型,然后根据目标运动轨迹和接收电压的变化曲线,获得目标相对于测控站的真实位置,利用形状特征匹配模型对其定位,最后驱动天线向目标运动完成搜索捕获,本发明能有效提高天线的搜索范围,缩短捕获时间,提高捕获自动化成功率。
本发明解决其技术问题所采用的技术方案包括以下步骤:
步骤一,收集测控天线A的特征参数T={f,r},f是该天线的接收工作频率,r是该天线的半径;
步骤二,根据特征参数T生成天线A的方向图P={(a1,e1|g11),…,(ai,ej|gij),…,(am,en|gmn)},其中,m是天线方向图方位方向的点数量,n是天线方向图俯仰方向的点数量,ai是天线方向图在方位方向的坐标,ej是天线方向图在俯仰方向的坐标,gij是该点对应的增益值,一个方位、俯仰坐标对(ai,ej)对应该点的天线增益值gij;从方向图P中,获取(ai,ej)点对应的天线增益gij=FP(ai,ej);
步骤三,根据方向图P生成搜索结构S={(a1,e1|g11 0,g11 1,g11 2,g11 3,g11 4),(ai,ej|gij 0,gij 1,gij 2,gij 3,gij 4),…,(am,en|gmn 0,gmn 1,gmn 2,gmn 3,gmn 4)},对于S中每一个点(ai,ej),其对应的gij 0=FP(ai,ej),gij 1=FP(ai+1,ej)-FP(ai,ej),gij 2=FP(ai-1,ej)-FP(ai,ej),gij 3=FP(ai,ej+1)-FP(ai,ej),gij 4=FP(ai,ej-1)-FP(ai,ej),即gij 0是(ai,ej)点的增益值,gij 1、gij 2、gij 3、gij 4是该点与其四个相邻点之间的增益差值,表示(ai,ej)点向四个不同方向继续搜索得到的增益差值,
步骤四,保持天线A静止,目标与天线相对运动形成的接收信噪比曲线点集为Vcurve={v1,…,vp},p是接收信噪比曲线中点的数量;
步骤五,将Vcurve信噪比序列转换为增益差值序列Dcurve={d1,…,dp-1},其中di=vi+1-vi,0≤i≤p-1;
步骤六,在天线A的方向图搜索结构S中查找一段曲线Psearch={(ai,ej),…,(ai+p-1,ej+p-1)},(ai,ej)为搜索路径的起点,(ai+p-1,ej+p-1)为搜索路径的终点,使得点集Psearch中所有点形成的搜索结果集合Dsearch={dij,…,di+p-1,j+p-1}与Dcurve的欧拉距离最小,即Dis(Dsearch,Dcurve)取得最小值;
步骤七,则控制天线向点(ai+p-1,ej+p-1)运动搜索目标。
本发明的有益效果是:利用建立天线的方向图模型,然后将目标与天线相对运动形成的电压的变化率曲线作为待检测曲线,在天线的方向图模型中搜索最符合该曲线特征的起始点和终止点,最后利用终止点确定目标与天线电轴的相对位置,从而控制天线完成目标的捕获跟踪。实际验证结果表明,当目标在天线副瓣范围内运动时,该方法能自动完成目标的捕获跟踪,而无需人工的扫描干预,同时该方法也可用于主副瓣判别。该方法能有效提高天线的搜索范围,缩短捕获时间,提高自动化成功率。
附图说明
图1是本发明的方法流程图;
图2是球坐标系下目标的位置表示。
图3是天线和路方向图;
图4是跟踪曲线在搜素结构中的搜索结果;
图5是天线跟踪目标曲线图。
具体实施方式
下面结合附图和实施例对本发明进一步说明,本发明包括但不仅限于下述实施例。
本发明的技术方案包括以下步骤:
步骤一:收集天线参数
选取一台测控天线A,其特征参数为T={f,r},f是该天线的接收工作频率,r是该天线的半径;
步骤二:生成天线方向图
根据特征参数T生成天线A的方向图P={(a1,e1|g11),…,(ai,ej|gij),…,(am,en|gmn)},其中,m是天线方向图方位方向的点数量,n是天线方向图俯仰方向的点数量,ai是天线方向图在方位方向的坐标,ej是天线方向图在俯仰方向的坐标,gij是该点对应的增益值,一个方位、俯仰坐标对(ai,ej)对应该点的天线增益值gij。从方向图P中,获取(ai,ej)点对应的天线增益gij值使用Fp表示,即gij=FP(ai,ej);
步骤三:生成搜索结构
根据方向图P生成搜索结构S={(a1,e1|g11 0,g11 1,g11 2,g11 3,g11 4),(ai,ej|gij 0,gij 1,gij 2,gij 3,gij 4),…,(am,en|gmn 0,gmn 1,gmn 2,gmn 3,gmn 4)},对于S中每一个点(ai,ej),其对应的gij 0=FP(ai,ej),gij 1=FP(ai+1,ej)-FP(ai,ej),gij 2=FP(ai-1,ej)-FP(ai,ej),gij 3=FP(ai,ej+1)-FP(ai,ej),gij 4=FP(ai,ej-1)-FP(ai,ej),即gij 0是(ai,ej)点的增益值,gij 1,gij 2,gij 3,gij 4是该点与其四个相邻点之间的增益差值。gij 1,gij 2,gij 3,gij 4表示(ai,ej)点向四个不同方向继续搜索得到的增益差值,例如当(ai,ej)点确定向方位正方向继续搜索时,得到该点的搜索结果dij=gij 1。
步骤四:收集跟踪曲线
天线A静止,目标与天线相对运动形成的接收信噪比曲线点集为Vcurve={v1,…,vp},p是接收信噪比曲线中点的数量;
步骤五:曲线形状转换
将Vcurve信噪比序列转换为增益差值序列Dcurve={d1,…,dp-1},其中di=vi+1-vi,0≤i≤p-1。
步骤六:相似序列搜索
在天线A的方向图搜索结构S中查找一段曲线Psearch={(ai,ej),…,(ai+p-1,ej+p-1)},(ai,ej)为搜索路径的起点,(ai+p-1,ej+p-1)为搜索路径的终点,使得点集Psearch中所有点形成的搜索结果集合Dsearch={dij,…,di+p-1,j+p-1}与Dcurve的欧拉距离最小,即Dis(Dsearch,Dcurve)取得最小值;
步骤七:驱动天线跟踪
在得到了搜索点集Psearch后,就能判断目标与天线的相对位置最可能在(ai+p-1,ej+p-1)点,则控制天线向该点运动搜索目标。
本发明的实施例包括以下步骤,算法流程如图1所示:
步骤一:天线参数收集
为检验算法的可行性,选取参数半径r=7.5m,接收频率f=2200MHZ的天线对卫星目标进行跟踪实验分析。
步骤二:天线方向图生成
TE11模和TE21模在球坐标下的远场辐射场强可用(1)-(4)式表示,方向图可用(5)式表示。
其中Jn(u)是n阶第一类Bessel函数,u=kasinθ,并且c为光速,f为接收频率,a为天线口径,θ,Φ为球坐标下目标相对天线的夹角,如图2所示。
同时由于(1)-(4)式是在球坐标系下的计算结果,需要将其转换为天线基座坐标系下的表达方式,如(6)和(7)式所示。
根据天线参数生成和路TE11模的方向图如图3所示,x轴为天线指向方位角度值,y轴为天线指向俯仰角度值,z轴为该点对应的天线增益值。在图3中示例的天线电轴指向方位角为22°,俯仰角为15°,在算法搜索过程中该角度值为天线静止等待目标运动时的实际指向角度值,可取为天线运动范围内的任意值。
步骤三:生成搜索结构
在生成天线的方向图G后,即可根据每一个点的增益值计算得到搜索结构S。对方向图中的每一个点分别计算与其相邻的方位正方向、方位负方向、俯仰正方向、俯仰负方向四个点之间的增益差值。
步骤四:收集跟踪曲线
首先天线转为指向状态,收集目标相对于天线运动的信噪比曲线V。
步骤五:曲线形状转换
将接收到的信噪比曲线中的每一个点依次与后一个点作差,得到增益差值序列Dcurve。
步骤六:相似序列搜索
在搜索结构S中查找与Dcurve最匹配的路径,算法步骤如下:
(1)将S中每个节点Pij=(ai,ei)作为搜索起点,i=1:m,j=1:n;
(2)搜索路径集Psearch和累计差值Disij清空初始化,将pij加入Psearch;
(3)对Dcurve中每个节点dk进行以下匹配,k=1:q;
(4)选取点pij四向增益差值gij 1,gij 2,gij 3,gij 4中与dk差值最小的值dij作为搜索结果,加入Dsearch,pij更新为搜索前进方向的下一个节点,将新的pij加入Psearch;
(5)Disij累加运算,Disij=Disij+dij;
(6)如果Dcurve中仍有节点未进行匹配,跳转到(3);
(7)Disij计算完成,如果S中仍有节点未作为搜索起始点,跳转到(1);
(8)查找具有最小值的Dismin=min{Disij,i=1:m,j=1:n}。
Dismin就是搜索结构S中与Dcurve最小的搜索路径累计差值,其对应的搜索路径点集Psearch,就是天线指向静止期间,目标相对于天线的运动轨迹,如图4中所示。
步骤七:驱动天线跟踪
查找到配置的搜索路径路径后,就得到目标相对于天线电轴的方位、俯仰位置差值,最后控制天线向目标运动进行捕获跟踪,如图5所示,A点为目标相对于天线的位置,B点为天线完成跟踪目标位置。整个算法流程结束。
Claims (1)
1.一种基于形状特征搜索的近地卫星目标自动捕获算法,其特征在于包括下述步骤:
步骤一,收集测控天线A的特征参数T={f,r},f是该天线的接收工作频率,r是该天线的半径;
步骤二,根据特征参数T生成天线A的方向图P={(a1,e1|g11),…,(ai,ej|gij),…,(am,en|gmn)},其中,m是天线方向图方位方向的点数量,n是天线方向图俯仰方向的点数量,ai是天线方向图在方位方向的坐标,ej是天线方向图在俯仰方向的坐标,gij是该点对应的增益值,一个方位、俯仰坐标对(ai,ej)对应该点的天线增益值gij;从方向图P中,获取(ai,ej)点对应的天线增益gij=FP(ai,ej);
步骤三,根据方向图P生成搜索结构S={(a1,e1|g11 0,g11 1,g11 2,g11 3,g11 4),(ai,ej|gij 0,gij 1,gij 2,gij 3,gij 4),…,(am,en|gmn 0,gmn 1,gmn 2,gmn 3,gmn 4)},对于S中每一个点(ai,ej),其对应的gij 0=FP(ai,ej),gij 1=FP(ai+1,ej)-FP(ai,ej),gij 2=FP(ai-1,ej)-FP(ai,ej),gij 3=FP(ai,ej+1)-FP(ai,ej),gij 4=FP(ai,ej-1)-FP(ai,ej),即gij 0是(ai,ej)点的增益值,gij 1、gij 2、gij 3、gij 4是该点与其四个相邻点之间的增益差值,表示(ai,ej)点向四个不同方向继续搜索得到的增益差值,
步骤四,保持天线A静止,目标与天线相对运动形成的接收信噪比曲线点集为Vcurve={v1,…,vp},p是接收信噪比曲线中点的数量;
步骤五,将Vcurve信噪比序列转换为增益差值序列Dcurve={d1,…,dp-1},其中di=vi+1-vi,0≤i≤p-1;
步骤六,在天线A的方向图搜索结构S中查找一段曲线Psearch={(ai,ej),…,(ai+p-1,ej+p-1)},(ai,ej)为搜索路径的起点,(ai+p-1,ej+p-1)为搜索路径的终点,使得点集Psearch中所有点形成的搜索结果集合Dsearch={dij,…,di+p-1,j+p-1}与Dcurve的欧拉距离最小,即Dis(Dsearch,Dcurve)取得最小值;
步骤七,则控制天线向点(ai+p-1,ej+p-1)运动搜索目标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611009936.6A CN106772466B (zh) | 2016-11-17 | 2016-11-17 | 一种基于形状特征搜索的近地卫星目标自动捕获算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611009936.6A CN106772466B (zh) | 2016-11-17 | 2016-11-17 | 一种基于形状特征搜索的近地卫星目标自动捕获算法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106772466A CN106772466A (zh) | 2017-05-31 |
CN106772466B true CN106772466B (zh) | 2019-01-11 |
Family
ID=58968654
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611009936.6A Expired - Fee Related CN106772466B (zh) | 2016-11-17 | 2016-11-17 | 一种基于形状特征搜索的近地卫星目标自动捕获算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106772466B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108828511B (zh) * | 2018-06-20 | 2022-04-01 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 全空域多波束协同的目标搜索方法 |
CN110380770B (zh) | 2019-06-10 | 2021-07-06 | 浙江大学 | 一种低轨移动卫星通信网络的自适应对星方法 |
CN110703283B (zh) * | 2019-08-15 | 2021-08-31 | 北京航天驭星科技有限公司 | 一箭多星任务中的卫星识别方法 |
CN111313959B (zh) * | 2020-02-20 | 2021-11-09 | 中国西安卫星测控中心 | 一种北斗三号相控阵捷变星间链路遥测下行路径搜索方法 |
CN113391138B (zh) * | 2020-03-13 | 2022-08-30 | 中国人民解放军63756部队 | 基于跟踪轨迹拟合的天线副瓣识别及其自动转主瓣方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101339435A (zh) * | 2008-08-08 | 2009-01-07 | 上海杰盛无线通讯设备有限公司 | 一种搜寻、追踪卫星的方法 |
CN103337704A (zh) * | 2013-02-19 | 2013-10-02 | 中国人民解放军66061部队 | Ku波段卫星通信便携站天线自动对星装置和方法 |
CN104659488A (zh) * | 2015-01-29 | 2015-05-27 | 四川安迪科技实业有限公司 | 一种卫星天线自动对星方法 |
-
2016
- 2016-11-17 CN CN201611009936.6A patent/CN106772466B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101339435A (zh) * | 2008-08-08 | 2009-01-07 | 上海杰盛无线通讯设备有限公司 | 一种搜寻、追踪卫星的方法 |
CN103337704A (zh) * | 2013-02-19 | 2013-10-02 | 中国人民解放军66061部队 | Ku波段卫星通信便携站天线自动对星装置和方法 |
CN104659488A (zh) * | 2015-01-29 | 2015-05-27 | 四川安迪科技实业有限公司 | 一种卫星天线自动对星方法 |
Non-Patent Citations (2)
Title |
---|
中继卫星天线指向控制策略研究;孙小松 等;《航空学报》;20040731;第25卷(第4期);第376-380页 |
基于梯度法的对空中目标搜索与跟踪;李正周 等;《系统仿真学报》;20070331;第19卷(第5期);第1078-1081页 |
Also Published As
Publication number | Publication date |
---|---|
CN106772466A (zh) | 2017-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106772466B (zh) | 一种基于形状特征搜索的近地卫星目标自动捕获算法 | |
CN108493610B (zh) | 一种相控阵天线自动对星方法及装置 | |
CN105180942B (zh) | 一种无人船自主导航方法及装置 | |
CN109813319B (zh) | 一种基于slam建图的开环优化方法及系统 | |
CN104777478B (zh) | 一种相控阵雷达搜索捕获目标的方法 | |
CN104869639B (zh) | 一种室内定位方法和装置 | |
US20170329351A1 (en) | Apparatus-assisted sensor data collection | |
CN1230692A (zh) | 气象观测雷达系统 | |
CN110958575A (zh) | 一种基于WiFi融合预测的定位方法及系统 | |
CN106028446B (zh) | 室内停车场定位方法 | |
WO2020248443A1 (zh) | 一种基于二维相控天线阵列的快速精确波束跟踪方法 | |
CN112432649A (zh) | 一种引入威胁因子的启发式无人机蜂群航迹规划方法 | |
CN107181521B (zh) | 星间链路天线扫描捕获装置及方法 | |
CN113890598A (zh) | 基于相控阵的混合扫描天线跟踪方法、系统、终端及介质 | |
CN110488809A (zh) | 一种室内移动机器人的自主建图方法和装置 | |
CN108267736A (zh) | 一种geo sar凝视成像模式方位模糊度确定方法 | |
CN105223591B (zh) | 一种用于减弱远近效应的伪卫星阵列天线方向图综合计算方法 | |
Nguyen et al. | Low speed vehicle localization using wifi fingerprinting | |
JP2009019984A (ja) | 目標観測レーダ装置及び目標追尾方法 | |
Zhou et al. | RSSI-based mobile robot navigation in grid-pattern wireless sensor network | |
CN110309627A (zh) | 小天体表面弹跳移动的可达区获取方法 | |
KR101598023B1 (ko) | 지향 벡터를 이용한 위성 추적 장치 및 그 방법 | |
RU2680969C1 (ru) | Способ навигации летательных аппаратов | |
CN114302359A (zh) | 一种基于WiFi-PDR融合的高精度室内定位方法 | |
Gawronski et al. | Antenna scanning techniques for estimation of spacecraft position |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190111 |