本文针对单棒磁天线在接收罗兰C信号时,会因来波方向的不同而产生相位翻转的特性,设计了基于幅度叠加的信号合成算法,并实现了罗兰C信号的合成处理。该算法在FPGA中实现简单,处理过程快速,具有较强的使用价值,已经在所研制的全向磁天线中得到了使用。 The single stick H-field antenna induce reversed phase as the waves directions are different in receiving Loran C signal. The signal synthesized algorithm based on amplitude superposition is designed, which actualize the synthesization of Loran C signal. This algorithm which process fast, can easily realized in FPGA, and has been applied in Loran C omni-directional H-field antenna.
罗兰C导航系统是一种中远程低频无线电导航系统,属于陆基、脉冲相位调制导航方式[
我国关于罗兰C磁天线的研究尚处于起步阶段,而国外关于磁天线的研究显著早于我国。以早在上世纪七十年代,美国就在其拖曳天线上采用了磁天线技术[
单棒磁天线具有方向性,要想实现全向收信,必须采用多根磁棒同时接收信号并进行信号合成处理。而单棒磁天线在接收信号时,对立波瓣信号具有相位相反的特性,这就给信号合成带来了麻烦。本文主要针对全向磁天线设计环节中的信号合成部分展开研究。
全向磁天线原理结构如图1所示。
两根正交放置的单棒磁天线感应到的信号经前置低噪声放大器放大、模拟带通滤波器滤波、自动增益放大器放大等信号调理过程后,对两个通道的模拟信号进行模数转换,完成双通道信号的数字量化;数字量化后的两路数字信号由高性能可编程逻辑器件FPGA进行处理,分别完成信号的数字滤波、数字陷波、数字移相等过程,完成两路信号的数字调理过程;通过FPGA器件的控制和处理,对两路信号进行搜索和识别,针对信号的相位变化特点,采用文中设计的信号合成方法进行信号合成;最后将合成后的信号进行数模转换,输出的模拟合成信号提供给罗兰C接收机进行定位解算。对于罗兰C数字接收机来说,合成信号可以直接供给接收机进行定位解算,不需进行数模转换。
对罗兰C脉冲组信号中的每个脉冲载频起始相位进行严格规定,使其按照一定的格式排列,称为脉冲组相位编码。通过相位编码,可实现罗兰C系统全自动接收并可有效抑制天波干扰。
罗兰C系统脉冲组采用二相二周期互补编码方式,其中“二相”是指载频起始相位有两种,载波起始相位为0 rad时,为正相位编码,用“+”表示,起始相位为π rad时,为负相位编码,用“–”表示;“二周期”是指编码周期等于两个GRI长度;“互补”是指第一个周期脉冲组中偶数位脉冲的相位与第二个周期脉冲组中相应的偶数位脉冲的相位相差π rad,而奇数位脉冲相位相同;主、副台信号的相位编码是不同的,各主台、各副台信号的相位编码是相同的。罗兰C脉冲群中各脉冲的相位编码如表1所列[7,8]。
图1. 全向磁天线结构设计图
表1. 罗兰C信号相位编码
罗兰C信号在到达单棒磁天线对立模型波瓣时有一个180˚的相位改变,这种相位特性就要求在设计全向磁天线时,必须考虑每根单棒磁天线接收信号的相位特性,即导航台发射信号的来波方向处于单棒磁天线不同的波瓣区域时,必须针对每根单棒磁天线接收到的信号,选择合适的相位进行信号合成处理。
对于十字交叉环路布置的全向磁天线来说,其接收信号的全向区域可分为四个象限,如图2所示,处于不同象限的信号在进行信号合成时,必须进行适当的相位变换[
Ⅰ象限:A + B Ⅱ象限:A – B
Ⅲ象限:–A – B Ⅳ象限:–A + B
全向磁天线两路信号合成不仅要完成信号幅度的合成,还要完成信号的相位合成。由于单棒磁天线的方向性和相位特性,使得两根磁棒接收相同台站发射信号的相位会发生180˚的改变,由此导致罗兰C台链信号的合成变得复杂[
针对这种情况,本文设计了基于幅度叠加的全向磁天线信号合成方法。该方法在信号合成过程中消除了相位变化的影响,只需简单改变接收机内部处理程序,就能实现合成信号的搜索、跟踪、识别和定位解算。
幅度叠加信号合成算法原理如图3所示。
两根正交放置的单棒磁天线感应到的信号经前置低噪声放大器放大、模拟带通滤波器滤波、自动增益放大器放大等信号调理过程后,对两个通道的模拟信号进行模数转换,完成双通道信号的数字量化;数字量化后的两路数字信号由高性能可编程逻辑器件FPGA进行处理,分别完成信号的数字滤波、数字陷波、数字移相等过程,完成两路信号的数字调理过程。在对信号调理的基础上,完成信号合成的第一步,对信号取绝对值,完成|dataA|和|dataB|运算,此举可将信号的相位特性消除;然后完成信号合成的第二步,幅度叠加,完成|dataA| + |dataB|运算,得到幅度叠加合成信号;合成信号输出供罗兰C接收机使用。
假设罗兰C台链分布情况及磁天线的放置情况如图4所示。主台信号、副台2信号到达磁天线时,位于磁天线的第Ⅰ象限;副台1信号到达磁天线时,位于磁天线的第Ⅳ象限。各台站发射信号的波形图和相位
图2. 全向磁天线象限分布
图3. 幅度叠加信号合成流程
图4. 台站及磁天线布置情况
编码情况如图5所示,向上突起表示载波起始相位为0 rad,为正相位编码,用“+”表示;向下突起表示载波起始相位为π rad时,为负相位编码,用“–”表示。
对于A、B磁棒来说,由于单棒磁天线的方向性和相位特性,其接收三个台站的信号的幅度和相位编码都会发生改变,其信号幅度和相位编码情况如图6所示。
信号合成的第一步:对信号取绝对值。
完成|dataA|和|dataB|运算,此举可将信号的相位特性消除,所得的结果如图7中所示。
此时罗兰C信号的波形已经发生了改变,信号的相位特性已经不再包含在信号波形中,包络形状也不再是“鱼”形,其波形如图8所示。对于现有接收机来说,这种包络形状难以完成信号的搜索、跟踪和定位解算,必须根据信号的形状设计一种新的算法,才能使用要求。
图5. 发射信号的波形图和相位编码情况图
图6. 信号dataA,dataB幅度和相位情况
图7. |dataA|和|dataB|波形图
图8. data同|data|波形比较
信号合成的第二步:幅度叠加。
完成|dataA| + |dataB|运算,所得的结果如图9中所示。
利用自制的罗兰C全向磁天线接收东海台链的导航台信号,GRI = 83900 μs,该台链的台站和接收地点的分布情况及磁天线的放置情况如图10所示。
由图中可以看出,接收信号地点距离主台——宣城台最近,距离第一副台——饶平台较远,距离第二副台——荣成台最远;磁天线接收信号时,宣城台和荣成台信号位于全向磁天线的第一象限,饶平台信号位于第二象限。
双通道采集的数据dataA和dataB的原始数据及其相位情况如图11所示,采样率为1 M,数据长度为
图9. 幅度叠加合成信号
图10. 台站、接收地点及磁天线位置分布图
图11. dataA和dataB原始信号
10 M,图中显示信号只为一个GRI的信号。由于各台链之间信号的分布关系,在该段GRI信号中,会有其他台链导航台发射的信号,如宣城台是一个双工台,既发射东海台链的主台信号,又发射北海台链的副1台信号。图中对东海台链的信号进行了标注,信号合成过程中,只是针对东海台链的信号进行处理,其他台链的信号可以不予考虑。
信号合成的第一步:对信号取绝对值。
完成|dataA|和|dataB|运算,此举可将信号的相位特性消除,所得的结果如图12中所示。
信号合成的第二步:幅度叠加。
完成|dataA| + |dataB|运算,所得的结果如图13中所示。
图12. |dataA| + |dataB|信号波形
图13. 幅度叠加合成信号波形
为了衡量合成信号定位结果的精确性,实现对幅度叠加合成信号的时差解算,本文在现有接收机时差测量程序的基础上进行修改,以满足幅度叠加合成信号的时差测量要求,采用该程序算得幅度叠加合成信号的定位时差分别为TD1 = 14988.559097 μs、TD2 = 33426.903251 μs。多种罗兰C接收机的多次试验结果表明,本地8390台链的实际测量时差为TD1 = 14988.5 μs,TD2 = 33427 μs,二者分别相差0.059097 μs与–0.096749 μs,充分证明基于幅度叠加的信号合成方法的正确性与有效性。
单棒磁天线由于其接收信号的方向性,导致其无法完成全方位接收某一台链三个导航台的信号,要实现全方位接收罗兰C信号,必须采用多根不同方向放置的磁棒接收信号,然后进行信号合成。由于单棒磁天线的相位特性,多根不同方向放置的磁棒在接收某一台站的信号时,相位会不同,如果采用模拟合成信号的方法,难以解决同一台站信号相位相反的问题,这就必须采用数字信号合成方法。本章提出了幅度叠加信号合成算法,在信号合成过程中消除了相位变化的影响,只需简单改变接收机内部处理程序,就能实现合成信号的搜索、跟踪、识别和定位解算。通过仿真数据和实测数据的验证,本章提出的信号合成算法都能获得较好的效果。该算法在FPGA中实现简单,处理过程快速,具有较强的使用价值,已经在所研制的全向磁天线中得到了使用。
[
[
[
[
[
[
[
[
[