一种GNSS 信号模拟器高精度实时伪距计算优化方法

沙海，周兵，仝海波，张国柱，欧钢

(1. 国防科技大学 电子科学与工程学院, 湖南 长沙，410073；2. 北京卫星导航中心，北京，100094)

GNSS 信号模拟器可以同时模拟美国GPS、俄罗斯CLONASS、欧洲GALILEO 和中国BDS 4 种全球卫星导航系统用户接收的信号，是卫星导航系统地面站和多模接收机等设备研制的不可或缺的关键设备。目前，国内外针对GNSS 信号模拟器的硬件平台和软件设计都进行了论述，并取得了一定研究成果[1-6]。其中高精度伪距计算及实时性是GNSS 信号模拟器的研究热点也是难点。由于GNSS 信号模拟器需要模拟全星座信号，其计算量将会成倍增加。同时，考虑仿真环境的逼真度，卫星轨道、电离层特性等也将使用精确模型计算，计算量将大幅度提升。如何在保证伪距精度的条件下，降低伪距计算过程的运算量是本文需要研究的内容。此外，伪距动态模拟常采用高阶累加器DDS 方法，该方法需要求解伪距动态模型参数[6]，传统方法多采用分别计算卫星、用户等各阶导数的方法，精度有限。为此，文献[7]针对高阶导数不易求解的问题，提出了一种基于三阶样条函数的伪距动态模型参数计算方法。但是使用该插值方法时，针对不同的动态环境所需计算的原始伪距点个数及时间间隔是有所差异的，这为进一步减少伪距运算量提供了优化空间。本文作者在分析伪距计算过程的基础上，提出了一种精度较高、运算量相对较小的实时伪距计算优化方法。

1 基本原理

GNSS 信号模拟器伪距计算过程如图1 所示。首先由用户建立仿真场景，包括卫星轨道模型、电离层模型、用户运动轨迹模型等；然后根据仿真模型参数计算固定时间间隔的可见链路伪距、功率衰减值、多普勒等信息；最后，通过插值算法计算伪距动态模型参数，包括伪距、速度、加速度、加加速度等。

图1 伪距计算过程Fig.1 Pseudorange calculation process

1.1 伪距计算

接收机在仿真时刻t 收到卫星在ts时刻发射的卫星信号时，其测得伪距ρ为

其中：R(t,ts)为卫星到接收机的几何距离；Δts为卫星钟差；Δtr为相对论效应延迟；Δτion和Δτtrop分别为电离层和对流层延迟；Δτother为信号模拟器通道零值等其他误差；c 为光速。

由式(1)可以看出：伪距计算本质上为求解卫星发射时刻ts的一元方程求根问题。目前多数文献采用定点迭代法进行求解[1-4]，其计算过程不再累述。使用定点迭代法时，存在算法易发散的问题[8]，但鲜有文献对伪距计算方法的收敛性进行证明，下面将根据卫星导航系统的特征，证明伪距计算方法的收敛性。

将式(1)改写为如下一元方程形式：

对g(ts)求导得：

1.2 伪距动态模型参数计算

为了满足高动态条件下的伪距模拟要求，通常使用伪距的3 阶泰勒展开多项式作为伪距动态模型[9-11]，如式(4)所示。该模型在信号处理中采用多阶DDS 累加器结构实现，由于忽略高阶项，需要周期计算伪距动态模型的各阶项系数。

式中：a(tN)为加速度；j(tN)为加加速度。

计算伪距动态模型参数时，根据已知固定时间间隔的伪距及其变化率，常采用Hermite 插值、样条插值等方法。其中分段三阶样条插值方法无需求解伪距变化率，且精度较高，是一种较理想的伪距动态模型参数计算方法[7]。同时，在工程实现中，GNU 科学库已提供样条插值算法的快速计算函数[12]，使用方便。

2 伪距计算优化方法

高精度GNSS 信号模拟器的伪距计算方法包含从仿真数学模型计算伪距至插值算法得出动态模型参数的整个过程。为了减小运算量，并保证伪距计算的高精度，可以对伪距计算过程进行优化，提出的伪距优化方法主要包含以下3 个方面的改进。

1) 使用插值算法计算卫星位置。高精度GNSS 信号模拟器要求使用轨道动力学积分方法计算卫星轨道，并且伪距计算时需要多次求解任意时刻的卫星位置。若均采用动力学积分方法计算，则运算量较大。一种合理的优化方法是采用动力学积分方法计算一定时间间隔的轨道数据，然后通过插值算法计算任意时刻的卫星位置。

2) 信道延迟量不参与迭代计算。电离层延迟、对流层延迟等信道延迟量与接收机同卫星之间的夹角有关。通常高精度伪距计算时，会迭代重新计算这些延迟量。但是，伪距迭代计算前后，夹角通常变化很小。所以，可以考虑在损失一定伪距计算精度的条件下，在迭代计算过程中不计算信道延迟，以减少运算量的优化方法。

3) 不同动态条件下的伪距时间间隔。根据接收机的运动速度，可以将信号模拟器的动态特性分为静态、中低动态和高动态3 种环境。通常伪距动态模型参数的更新频率固定不变，因此，为了满足高动态特性要求，需要计算时间间隔较小的伪距，以提高插值算法精度。然而，如果统一采用高动态条件下的伪距计算时间间隔，对于其他动态环境，将造成浪费。为此，针对不同动态特性选择合适的时间间隔计算伪距，是一种合理的优化方法。

3 仿真论证及性能分析

根据伪距计算优化方法提出的3 个改进方式，下面通过仿真实验，论证优化算法的具体步骤及其性能。

3.1 卫星轨道插值计算

采用卫星动力学模型计算卫星轨道时，可以准确计算出卫星的位置和速度，因此选择精度较高的分段三阶Hermite 插值算法作为卫星轨道的插值算法[13]。文献[13]指出当卫星轨道误差小于1 cm 时，对定位结果的影响可以忽略，因此，插值算法的精度要求满足小于1 cm。

仿真时，卫星动力学积分法使用的摄动力模型如表1 所示。实验方法为使用卫星动力学模型计算1 h，时间间隔5 ms 的原始轨道数据。然后分别选择不同时间间隔下的轨道数据作为样本数据，并使用Hermite插值算法得到5 ms 的插值数据。通过与原始轨道数据比较，得出最优时间间隔。仿真结果如图2 所示。

由图2 可以看出：采用分段3 阶Hermite 插值算法计算卫星轨道，当时间间隔取120 s 时，最大位置误差为4.46 mm，最大速度误差为0.44 mm/s，伪距计算误差小于6×10-7m。同时与动力学积分法计算相同时间长度的卫星轨道耗时比较，计算效率提高16 倍。

表1 摄动力模型Table 1 Model of forces

图2 卫星轨道插值算法仿真结果Fig.2 Simulation results of satellite orbits interpolation algorithm

3.2 信道延迟量影响

信道延迟量主要包括电离层延迟和对流层延迟。为了分析迭代计算伪距时，忽略信道延迟量对计算精度的影响，采用如下仿真方法。分别计算5 h 时间内的传统高精度伪距计算方法和优化方法的伪距值，并作差得到优化算法的误差。其中电离层延迟使用Klobuchar 模型，对流层延迟使用Hopfield 模型，仿真结果如图3 所示。

图3 误差项影响仿真结果Fig.3 Simulation results of error affecting

由图3 可以看出：伪距误差随着仰角的减少逐渐增大，最大误差在3×10-5m 以内。同时，通过统计1 000 次重复计算相同时刻伪距的计算时间，得出优化算法计算效率提高了约1 倍。

3.3 不同动态条件下的伪距时间间隔

根据已知时间间隔的伪距，采用分段三阶样条插值算法计算伪距动态模型参数时，针对不同的动态条件，可以选择不同的伪距时间间隔。实验仿真时，分别选择静态、中低动态和高动态3 种场景。其中，在中低动态条件下，接收机的运动速度为1 000 m/s。在高动态条件下，选择JPL 实验室推荐的一种具有典型意义的高动态运动模型[14]，该高动态环境的接收机速度的变化曲线如图4 所示。

图4 高动态环境速度变化曲线Fig.4 Curve of high dynamic environment velocity

表2 不同动态条件下的插值算法精度Table 2 Pseudorange accuracy under different dynamic environments

仿真时，接收机的初始位置设置为仿真开始时刻卫星星下点位置，并且接收机始终沿该时刻垂直高度方向运动。分别选择不同时间间隔条件的原始伪距数据，然后使用样条函数计算1 ms 间隔的伪距值。通过与原始数据比较，得出不同时间间隔条件下的伪距动态模型参数计算精度，仿真结果如表2 所示。

由表2 可见：对于静态条件，选择时间间隔20 s时，伪距计算精度就可以达到2.43×10-5m；对于中低动态时，选择时间间隔2 s 时，伪距计算精度可以达到2.96×10-4m；对于高动态条件时，时间间隔选择0.1 s 时，伪距计算精度可以达到3.71×10-4m。

4 结论

1) 根据卫星导航系统的特征及定点迭代法收敛定理，证明了伪距计算时的算法收敛性问题。

2) 为了减少伪距计算过程的运算量，提出了一种优化算法，其具体步骤为：①使用分段三阶Herimite插值算法计算任意时刻的卫星轨道，已知卫星轨道数据的时间间隔为120 s；②迭代计算卫星发射时刻过程时，忽略信道延迟量的影响；③使用三阶样条函数计算伪距动态模型参数时，对于静态条件，已知伪距时间间隔为20 s；对于中低动态条件，时间间隔为2 s；对于高动态条件，时间间隔为0.1 s。

3) 使用伪距计算优化算法时，计算误差小于1 mm(即伪距损失精度1 mm)，计算效率大幅度提高。

致谢：感谢中国科学院上海天文台黄勇博士在卫星轨道积分算法方面提供的帮助。

[1] 赵军祥. 高动态智能GPS 卫星信号模拟器软件数学模型研究[D]. 北京： 北京航空航天大学电子信息工程学院, 2003：55-67.ZHAO Junxiang. The study on mathematical model of high dynamic and intelligent GPS satellite signal simulator software[D]. Beijing： Beihang University. School of Electronic and Information Engineering,2003：55-67.

[2] 吕志成. 高动态卫星导航信号模拟器软件研究[D]. 长沙： 国防科学技术大学电子科学与工程学院,2006：23-30.LÜ Zhicheng. The research of the high dynamic satellite signal simulator software[D]. Changsha： School of National University of Defense Technology. School of Electronic Science and Engineering,2006：23-30.

[3] 李成军, 陆明泉, 冯振明.GPS 中频信号模拟器的数学模型及实现[J]. 清华大学学报(自然科学版), 2008, 48(10)： 1582-1585.LI Chengjun, LU Mingquan, FENG Zhenming. Mathematical model and realization for GPS IF signal simulator[J]. Journal of Tsinghua University(Science and Technology), 2008, 48(10)：1582-1585.

[4] 罗益鸿. 导航卫星信号模拟器软件设计与实现[D]. 长沙： 国防科学技术大学机电工程与自动化学院,2008：29-43.LUO Yihong. The software design and implementation of navigation satellite signal simulator[D]. Changsha： National University of Defense Technology. College of Mechatronics Engineering and Automation,2008：29-43.

[5] 沙海, 吕志成, 张国柱, 等. 多系统卫星导航信号模拟源的系统结构与软件设计[C]// 第二届卫星导航学术年会. 上海,2011：341-345.SHA Hai, LÜ Zhicheng, ZHANG Guozhu, et al. System architecture and software design of multi-system satellite navigation signal simulator[C]// The Second China Satellite Navigation Conference.Shanghai,2011：341-345.

[6] Song Y,Zhou H,Zeng T,et al.Algorithm and realization of high dynamic satellite signal doppler simulation based on FPGA[C]//Proceedings of the 2010 International Technical Meeting of Institute of Navigation.San Diego,USA,2010：1044-1050.

[7] 刘旻, 吴嗣亮. 实时高动态GNSS 信号模拟器高精度伪距生成方法[J]. 北京理工大学学报,2011,31(9)：1053-1057.LIU Min, WU Siliang. Method of high accuracy pseudorange generation for real-time high dynamic GNSS signal simulator[J].Transactions of Beijing Institute of Technology, 2011, 31(9)：1053-1057.

[8] 丁丽娟. 数值计算方法[M]. 北京： 北京理工大学出版社,1997：184-187.DING Lijuan. Numerical computation method[M]. Beijing：Beijing Institute of Technology Press,1997：184-187.

[9] 齐巍, 常青, 洪晟. 导航信号模拟器可见卫星同步调度算法[J]. 沈阳工业大学学报,2011,33(4)：428-433.QI Wei, CHANG Qing, HONG Sheng. Algorithm for synchronous scheduling of visible satellite based on navigation signal simulator[J]. Journal of Shenyang University of Technology,2011,33(4)：428-433.

[10] ZHAO Yun, KOU Yanhong, HUANG Zhigang. Multi-Constellation GNSS signal simulator and its software implementation[C]// 22nd International Meeting of the Satellite Division of the Institute of Navigation. Savannah GA, 2009：2213-2219.

[11] 宋媛媛, 曾大治, 曾涛. 基于三阶DDS 的卫星信号多普勒模拟方法[J]. 北京理工大学学报,2010,30(10)：1213-1216.SONG Yuanyuan, ZENG Dazhi, ZENG Tao. Satellite signal doppler simulation method based on third-order DDS[J].Transactions of Beijing Institute of Technology, 2010, 30(10)：1213-1216.

[12] Mark Galassi. GNU scientific library[EB/OL]. [2009-08-25].http：//www.gnu.org.

[13] Korvenoja P,Piche R.Efficient satellite orbit approximation[C]//The 13th International Technical Meeting of the Satellite Division of The Institute of Navigation. Salt Lake City, USA,Sept,2000：1930-1937.

[14] 向洋. 高动态GPS 载波跟踪技术研究[D]. 武汉： 华中科技大学电信系,2010：22-52.XIANG Yang. Study on high dynamic GPS carrier tracking techniques[D]. Wuhan： Huazhong University of Science and Technology. Department of Electronic and Information Engineering,2010：22-52.