CN103545612A - 一种用于动中通系统的卫星波束中心估计方法 - Google Patents

一种用于动中通系统的卫星波束中心估计方法 Download PDF

Info

Publication number
CN103545612A
CN103545612A CN201310517900.9A CN201310517900A CN103545612A CN 103545612 A CN103545612 A CN 103545612A CN 201310517900 A CN201310517900 A CN 201310517900A CN 103545612 A CN103545612 A CN 103545612A
Authority
CN
China
Prior art keywords
scanning
beam center
angle
signal
center position
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.)
Granted
Application number
CN201310517900.9A
Other languages
English (en)
Other versions
CN103545612B (zh
Inventor
周志强
石永生
汪渤
董明杰
高志峰
沈军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN201310517900.9A priority Critical patent/CN103545612B/zh
Publication of CN103545612A publication Critical patent/CN103545612A/zh
Application granted granted Critical
Publication of CN103545612B publication Critical patent/CN103545612B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Radio Relay Systems (AREA)

Abstract

本发明提供一种用于动中通系统的卫星波束中心估计方法,该方法首先对采样得到的信标信号进行适当的滤波处理,通过简单的扫描运动使天线通过卫星波束中心位置,并在扫描过程中采集各点的扫描角度、角速率和信号强度数据,最后利用一阶统计矩准确地计算出扫描到的波束中心位置。本发明解决了现有技术中的动中通系统卫星波束中心扫描计算方法容易受噪声干扰、精度不高、控制过程复杂等问题。本发明的用于动中通系统的卫星波束中心估计方法在噪声干扰下对不同波束宽度的卫星信号都能准确地估计出波束中心位置,具有扫描运动控制简单,计算量小,精度高等优点。

Description

一种用于动中通系统的卫星波束中心估计方法
技术领域
本发明属于动中通系统的天线控制与信号处理领域,尤其涉及一种用于动中通系统的卫星波束中心估计方法。
背景技术
动中通是“移动中的卫星地面站通信系统”的简称。车辆、轮船和飞机等载体通过动中通系统在运动过程中保持天线指向实时跟踪卫星信号方向,不间断地传递语音、数据、图像等多媒体信息,可满足各种军民用应急通信和移动条件下的多媒体通信需要。动中通系统通常包括惯性稳定平台和信号扫描机构两大部分。惯性稳定平台通过陀螺、加速度计等惯性传感器构成的控制系统隔离载体运动,保持天线姿态稳定。信号扫描机构在天线指向偏离卫星波束中心导致卫星信号强度减弱时,控制天线对卫星信号进行扫描,重新确定并控制天线指向信号最大方向。
天线信标接收机输出的信标信号强度反映了当前接收到的卫星信号强度。当天线指向对准卫星波束中心时,信标机输出信号最强。在天线扫描过程中,以信标信号作为参考信号可以判断天线指向相对于卫星波束中心的变化。动中通系统在动态跟踪过程中,通过惯性稳定平台隔离载体姿态的变化,使天线指向始终保持在对准卫星的方向上。但是,由于惯性系统精度和陀螺漂移的影响,仪靠平台的稳定难以满足卫星通信对跟踪精度的要求,在此基础上还需监控卫星信号的变化,通过信号扫描获取跟踪误差,对天线指向进行相应的调整。
常用的信号扫描方式有矩形扫描、三角形扫描、圆锥形扫描等。矩形或三角形轨迹的扫描控制较为简单,圆锥形轨迹的扫描控制则比较复杂,它要求天线方位和俯仰角机械运动相配合,从而使天线极化轴顶点形成圆形轨迹。天线扫描过程中的进动一般分为搜索和调整两种,搜索阶段的每一步对信标电平信号进行取样和保存,经过若干次搜索和对扫描信号值的计算可确定卫早信号方向。根据搜索结果对天线角度进行相应调整,最终使其指向卫星信号最大方向。
在搜索过程中,现有的信号扫描估计方法通常是对各采样点的信号强度进行比较,根据比较结果逐步使天线指向信号最强点。由于信号噪声和信号本身某些特性的影响,这种方式容易出现误动作、不能准确寻找到卫星波束中心等问题。圆锥形扫描方法在计算过程中虽然也进行过一些滤波处理,能够从理论上解调出信号强度最大方向,但在实际中由于受信标信号质量、扫描频率和扫描控制准确度等因素影响,对卫星波束中心的计算精度也并不理想。
发明内容
为解决上述问题,本发明提供一种用于动中通系统的卫星波束中心估计方法,能使动中通系统在噪声干扰下对不同波束宽度的卫星信号通过扫描计算都能准确地估计出波束中心位置,从而提高动中通系统的跟踪精度。
本发明的用于动中通系统的卫星波束中心估计方法,包括:
步骤1,当天线指向偏离卫星波束中心时,控制所述天线朝任一方向进行方位角扫描,并以固定的采样频率获取每一采样时刻ti的信标信号ui,然后对所述信标信号ui根据公式
Figure BDA0000403587660000021
进行滤波处理,其中,扫描初始时刻t0的信标信号记为u0,且
Figure BDA0000403587660000022
α为滤波参数,
Figure BDA0000403587660000023
表示滤波后的信号;
步骤2,保持步骤1中的所述任一方向继续扫描直至
Figure BDA0000403587660000024
并记录扫描至
Figure BDA0000403587660000025
这一采样时刻的信标信号为U0,且U0=ui,然后控制所述天线进行反方向扪描,其中所述ε为设定的信号阈值;
步骤3,保持步骤2中的所述反方向扫描直至
Figure BDA0000403587660000026
则天线已通过方位角波束中心位置,扫描结束;记录扫描至
Figure BDA0000403587660000027
这一采样时刻的信标信号ui、方位角θi和扫描角速率ωi,将满足条件ui>U0的信标信号记为u′i,利用公式
Figure BDA0000403587660000028
计算信标信号关于扫描角度和角速率的一阶统计矩M,利用公式
Figure BDA0000403587660000029
计算相应的信号累加量N,利用公式θs=M/N计算方位角波束中心位置θs
步骤4,根据步骤3中计算获得的方位角波束中心位置θs控制所述天线做姿态调整,使方位角对准方位角波束中心位置θs
步骤5,根据步骤1至步骤4的方法控制所述天线进行俯仰角扫描,计算获得俯仰角波束中心位置θs’,使俯仰角对准俯仰角波束中心位置θs’:
步骤6,根据扫描估计获得的所述方位角波束中心位置θs和所述俯仰角波束中心位置θs’得到最终的卫星波束中心位置。
进一步的,所述步骤1中的滤波参数α的取值范围为0.02到0.05:
所述步骤3中的扫描角速率ωi由方位陀螺或俯仰轴陀螺获得。
本发明的有益效果在于:
本发明所给出的用于动中通系统的卫星波束中心估计方法具有扫描运动控制简单,计算量小,精度高等优点。
本发明利用统计方法估计卫星波束中心位置,比基于信号强度比较和曲线拟合的方法具有更好的鲁棒性和适用性,在噪声干扰下具有更高的估计精度。本发明在统计过程中考虑了扫描角速率信息,避免了扫描角速率的变化对估计精度产生影响。
相比于传统的步进式扫描,本发明方法不容易出现误动作,在噪声干扰下具有良好的鲁棒性和更高的估计精度;与圆锥扫描相比,控制更为简单,估计精度也要更高。
附图说明
图1(a)为本发明的实施例一中扫描过程记录下的信标信号示意图;
图1(b)为本发明的实施例一中扫描估计得到的卫星方位角波束中心位置示意图;
图2(a)为本发明的实施例二中扫描过程记录下的信标信号示意图:
图2(b)为本发明的实施例二中扫描估计得到的卫星方位角波束中心位置示意图。
具体实施方式
现有动中通系统卫星波束中心扫描估计的精度容易受噪声影响,或对扫描控制的精度和速度要求较高。当卫星信号波束较宽,波束中心的信号强度没有明显强于周边信号时,现有方法也很难通过信号强度的比较准确地计算出波束中心位置。本发明首先对采样得到的信标信号进行适当的滤波处理,通过简单的扫描运动使天线通过卫星波束中心位置。扫描过程中采集各点的扫描角度、角速率和信号强度数据,最后利用一阶统计矩准确地估计出扫描到的波束中心位置。
本发明在天线方位和俯仰角上采用相同的方法进行扫描计算,可分别得出方位和俯仰角波束中心位置,利用动中通系统接收亚太五号Ku波段卫星信号。下而以对方位角波束中心扫描估计说明本实施方式的具体步骤:
步骤一、当天线指向偏离卫星波束中心时,控制天线朝某一方向开始进行方位角扫描运动。以固定的采样频率获取每一采样时刻ti的信标信号ui,扫描初始时刻t0的信标信号记为u0。对信标信号ui进行如下的滤波处理:
u ^ i + 1 = ( 1 - α ) u ^ i + α u i + 1 .
式中,
Figure BDA0000403587660000042
α为滤波参数,
Figure BDA0000403587660000043
表示滤波后的信号。
步骤二、保持当前方向继续扫描,直至
Figure BDA0000403587660000044
(ε为设定的信号阈值),表明卫星信号强度有明显下降,控制天线进行反方向扫描。记录当前采样时刻的信标信号为U0(U0=ui)。
步骤三、保持步骤二中的所述反方向扫描直至
Figure BDA0000403587660000045
表明天线已通过波束中心位置,扫描结束。在此过程中记录所有采样时刻的信标信号ui、方位角θi和扫描角速率ωi,将满足条件ui>U0的信标信号记为u′i,按如下公式计算信标信号关于扫描角度和角速率的一阶统计矩M及对应的信号累加量N:
M = Σ u i ′ > U 0 θ i ω i u i ′ ,
N = Σ u i ′ > U 0 ω i u i ′ ,
由公式θs=M/N估计得到方位角上的波束中心位置θs
步骤四、控制天线做姿态调整,使方位角对准波束叫心位置θs
下面对该实施例中的其他技术细节问题作详细说明:
(1)步骤二所述的滤波处理的目的是为了得到更为平滑的滤波信号从而能够在扫描过程中对信标信号强度是否上升或下降做出更准确的判断。滤波参数α的合适取值范围为0.02到0.05。本实施例中α取为0.04。该滤波方法在平滑信标信号的同时,也会造成信标信号滞后。所以,在估计波束中心位置θs时不采用滤波信号
Figure BDA0000403587660000049
而是采用原信号ui
(2)计算一阶统计矩M和信号累加量N时需要考虑扫描角速率ωi。由于采样时间间隔固定,计算公式中加入ωi可以避免扫描角速率变化对波束中心估计精度产生影响,同时也使天线的扫描运动不需要保持恒定角速率,控制更为简单。本实施例中,扫描角速率由方位陀螺给出。
利用本发明方法按上述实施方式对卫星方位角波束中心进行扫描估计实验。图1(a)是.上述步骤三中的扫描过程记录下的信标信号曲线,其中0°角位置是利用频谱仪标定出的理想卫星方位角波束中心位置。从图中可以看出,上述扫描过程经过了卫星方位角波束中心位置,但是信标信号具有较强的噪声干扰,且在波束中心附近未出现明显的波峰,给波束中心位置的精确估计带来了很大的困难。图1(b)中粗曲线为滤波后的信号曲线。该信号曲线较为平滑,虽然比原信号略有滞后,但是利用它可以判断出当前信号强度是否具有明显的上升或下降趋势。图中虚线以上的信号为强度大于U0的有效信号,利用这些信号结合扫描角度和角速率信息,最后计算得出方位角波束中心位置为θs=0.0047°(如图中竖直箭头所示),可以看出本发明方法具有较高的扫描估计精度。
图2(a)为另一实验过程记录下的信标信号曲线。该信标信号在方位角波束中心附近出现了不稳定的跳变,但是并不影响本发明方法对方位角波束中心的扫描。如图2(b)所示,该扫描过程经过了方位角波束中心位置,并在信号强度出现明显下降后停止。图中竖直箭头标示出了扫描估计得到的方位角波束中心位置θs=-0.0190°。该实验表明本发明方法在强干扰作用下仍具有较高的精度。
上述实施例给出的是方位角波束中心的扫描估计方法,而本发明方法是需要对方位和俯仰角波束中心位置采用所述方法分别进行扫描估计的,这样才能得出完整的波束中心位置。采用所述方法对俯仰角波束中心位置扫描估计的具体步骤:
步骤一、当天线方位角对准波束中心位置θs时,控制天线朝某一方向开始进行俯仰角扫描运动。以固定的采样频率获取每一采样时刻ti的信标信号ui,扫描初始时刻t0的信标信号记为u0。对信标信号ui进行如下的滤波处理:
u ^ i + 1 = ( 1 - α ) u ^ i + α u i + 1 .
式中,
Figure BDA0000403587660000052
α为滤波参数,
Figure BDA0000403587660000053
表示滤波后的信号。
步骤二、保持当前方向继续扫描,直至
Figure BDA0000403587660000054
(ε为设定的信号阈值),表明卫星信号强度有明显下降,控制天线进行反方向扫描。记录当前采样时刻的信标信号为U0(U0=ui)。
步骤三、保持步骤二中的所述反方向扫描直至
Figure BDA0000403587660000061
表明天线已通过俯仰角波束中心位置,扫描结束。在此过程中记录所有采样时刻的信标信号ui、俯仰角θ′i和扫描角速率ω′i,将满足条件ui>U0的信标信号记为u′i,按如下公式计算信标信号关于扫描角度和角速率的一阶统计矩M及对应的信号累加量N:
M = Σ u i ′ > U 0 θ i ω i u i ′ ,
N = Σ u i ′ > U 0 ω i u i ′ ,
由公式θ′s=M/N估讣得到俯仰角上的波束中心位置θ′s
步骤四、控制天线做姿态调整,使俯仰角对准波束中心位置θ′s
根据上述对方位角和俯仰角波束中心位置扫描的方法,进行扫描估计获得的所述方位角波束中心位置θs和所述俯仰角波束中心位置θs’来获得最终的卫星波束中心位置。
当然,本发明还可有其他多种实施例,在不背离本发明精神及其实质的情况下,熟悉本领域的技术人员当可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。

Claims (2)

1.一种用于动中通系统的卫星波束中心估计方法,其特征在于,包括:
步骤1,当天线指向偏离卫星波束中心时,控制所述天线朝任一方向进行方位角扫描,并以固定的采样频率获取每一采样时刻ti的信标信号ui,然后对所述信标信号ui根据公式
Figure FDA0000403587650000011
进行滤波处理,其中,扫描初始时刻t0的信标信号记为u0,且
Figure FDA0000403587650000012
α为滤波参数,
Figure FDA0000403587650000013
表示滤波后的信号;
步骤2,保持步骤1中的所述任一方向继续扫描直至
Figure FDA0000403587650000014
并记录扫描至这一采样时刻的信标信号为U0,且U0=ui,然后控制所述天线进行反方向扫描,其中所述ε为设定的信号阈值;
步骤3,保持步骤2中的所述反方向扫描直至
Figure FDA0000403587650000016
则天线已通过方位角波束中心位置,扫描结束;记录扫描至
Figure FDA0000403587650000017
这一采样时刻的信标信号ui、方位角θi和扫描角速率ωi,将满足条件ui>U0的信标信号记为u′i,利用公式
Figure FDA0000403587650000018
计算信标信号关于扫描角度和角速率的一阶统计矩M,利用公式
Figure FDA0000403587650000019
计算相应的信号累加量N,利用公式θs=M/N计算方位角波束中心位置θs
步骤4,根据步骤3中计算获得的方位角波束中心位置θs控制所述天线做姿态调整,使方位角对准方位角波束中心位置θs
步骤5,根据步骤1至步骤4的方法控制所述天线进行俯仰角扫描,计算获得俯仰角波束中心位置θs’,使俯仰角对准俯仰角波束中心位置θs’;
步骤6,根据扫描估计获得的所述方位角波束中心位置θs和所述俯仰角波束中心位置θs’得到最终的卫星波束中心位置。
2.如权利要求1所述的用于动中通系统的卫星波束中心估计方法,其特征在于,所述步骤1中的滤波参数α的取值范围为0.02到0.05;
所述步骤3中的扫描角速率ωi由方位陀螺或俯仰轴陀螺获得。
CN201310517900.9A 2013-10-28 2013-10-28 一种用于动中通系统的卫星波束中心估计方法 Active CN103545612B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310517900.9A CN103545612B (zh) 2013-10-28 2013-10-28 一种用于动中通系统的卫星波束中心估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310517900.9A CN103545612B (zh) 2013-10-28 2013-10-28 一种用于动中通系统的卫星波束中心估计方法

Publications (2)

Publication Number Publication Date
CN103545612A true CN103545612A (zh) 2014-01-29
CN103545612B CN103545612B (zh) 2015-09-09

Family

ID=49968829

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310517900.9A Active CN103545612B (zh) 2013-10-28 2013-10-28 一种用于动中通系统的卫星波束中心估计方法

Country Status (1)

Country Link
CN (1) CN103545612B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111812683A (zh) * 2020-09-11 2020-10-23 四川九洲电器集团有限责任公司 一种低轨卫星实时跟踪方法、装置、设备及介质
CN115986407A (zh) * 2023-03-20 2023-04-18 四川省华盾防务科技股份有限公司 相控阵系统的控制方法、装置及存储介质
CN106954223B (zh) * 2017-04-17 2023-11-14 南京京迪通信设备有限公司 一种动中通端站系统及动中通端站系统的通信方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0755578B1 (en) * 1994-04-08 2002-09-04 Ericsson Inc. Multiple beamwidth phased array
CN1849529A (zh) * 2003-09-10 2006-10-18 皇家飞利浦电子股份有限公司 利用多个同步波束发送的超声空间组合
US20090199642A1 (en) * 2006-07-11 2009-08-13 Hiroyuki Fukutomi Ultrasonic flaw detection apparatus and ultrasonic flaw detection method
CN102859901A (zh) * 2010-04-08 2013-01-02 波音公司 地理位置调节点波束交叠

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0755578B1 (en) * 1994-04-08 2002-09-04 Ericsson Inc. Multiple beamwidth phased array
CN1849529A (zh) * 2003-09-10 2006-10-18 皇家飞利浦电子股份有限公司 利用多个同步波束发送的超声空间组合
US20090199642A1 (en) * 2006-07-11 2009-08-13 Hiroyuki Fukutomi Ultrasonic flaw detection apparatus and ultrasonic flaw detection method
CN102859901A (zh) * 2010-04-08 2013-01-02 波音公司 地理位置调节点波束交叠

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106954223B (zh) * 2017-04-17 2023-11-14 南京京迪通信设备有限公司 一种动中通端站系统及动中通端站系统的通信方法
CN111812683A (zh) * 2020-09-11 2020-10-23 四川九洲电器集团有限责任公司 一种低轨卫星实时跟踪方法、装置、设备及介质
CN111812683B (zh) * 2020-09-11 2020-12-15 四川九洲电器集团有限责任公司 一种低轨卫星实时跟踪方法、装置、设备及介质
CN115986407A (zh) * 2023-03-20 2023-04-18 四川省华盾防务科技股份有限公司 相控阵系统的控制方法、装置及存储介质

Also Published As

Publication number Publication date
CN103545612B (zh) 2015-09-09

Similar Documents

Publication Publication Date Title
CN109120561B (zh) 低轨卫星移动通信系统中频偏校正方法
CN104681986B (zh) 卫星天线在倾斜状态下的寻星方法
CN108061888A (zh) 球面相控阵天线对航天器动态目标的自跟踪方法
CN110441764A (zh) 一种针对单通道接收机在多径信道下稳定跟踪的方法和装置
CN106872967B (zh) 一种基于双基地雷达的动目标检测系统及方法
CN110045361A (zh) 球面相控阵单脉冲数字跟踪系统
CN107271997B (zh) 机载多通道cssar地面运动目标运动参数估计方法
CN110045328B (zh) 基于小型无人机平台的反黑飞无人机探测定位方法
CN112013832B (zh) Victs天线自适应对星跟踪方法
CN103545612B (zh) 一种用于动中通系统的卫星波束中心估计方法
WO2019071917A1 (zh) 卫星跟踪方法
Knoedler et al. Detecting and tracking a small uav in gsm passive radar using track-before-detect
CN109613487A (zh) 一种利用卫星侦测信息提升雷达快速截获重点目标的方法
CN102142853B (zh) 单脉冲跟踪接收机系统误差矩阵补偿方法
CN102253382B (zh) 一种数字中频单脉冲测向方法
CN115336430B (zh) Ka波段底视合成孔径cm级测高装置及测高方法
CN110988858B (zh) 一种双波束微波着陆雷达高精度测距方法及系统
CN106532266A (zh) 一种能够防止vsat天线误指向的卫星捕获方法
US9720077B1 (en) Radio altimeter for detecting accurate height
CN116794616A (zh) 一种基于脱靶量融合计算的无人机校相方法、设备、装置及存储介质
CN115882881A (zh) 一种强干扰下阵列天线时分系统初始角度搜索方法和终端
CN111190178A (zh) 一种相频扫三坐标雷达的多波束顺序比幅测高方法
CN107894589B (zh) 基于双频连续波应答机天线的运载火箭二维姿态测量方法
CN114910944A (zh) 一种相控阵天线跟踪精度测量方法
CN115166715A (zh) 连续波相控阵系统的信号检测和跟踪方法及装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant