CN110658536B - 一种卫星接收机导航输出信号时序测试方法 - Google Patents

一种卫星接收机导航输出信号时序测试方法 Download PDF

Info

Publication number
CN110658536B
CN110658536B CN201810700852.XA CN201810700852A CN110658536B CN 110658536 B CN110658536 B CN 110658536B CN 201810700852 A CN201810700852 A CN 201810700852A CN 110658536 B CN110658536 B CN 110658536B
Authority
CN
China
Prior art keywords
signal
sampling
test
navigation
time sequence
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
Application number
CN201810700852.XA
Other languages
English (en)
Other versions
CN110658536A (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 Automation Control Equipment Institute BACEI
Original Assignee
Beijing Automation Control Equipment Institute BACEI
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 Automation Control Equipment Institute BACEI filed Critical Beijing Automation Control Equipment Institute BACEI
Priority to CN201810700852.XA priority Critical patent/CN110658536B/zh
Publication of CN110658536A publication Critical patent/CN110658536A/zh
Application granted granted Critical
Publication of CN110658536B publication Critical patent/CN110658536B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/23Testing, monitoring, correcting or calibrating of receiver elements

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)
  • Radio Relay Systems (AREA)

Abstract

本发明公布了一种卫星接收机导航输出信号的时序测试方法,属于室温工作环境下卫星接收机信号时序测试领域。针对卫星接收机信号时序测试工作,本发明提高了科研生产的自动化水平,本发明摆脱了对示波器的依赖,开发了一种新的测试计算方法:首先,通过时域采样得到数字信号形式的数据;然后,利用数字信号的特征和待测项目设计一种计算方法;最后,在测试计算机上实现了一键自动化计算、显示测试结果。本发明既可以用于检验卫星导航接收机输出时序是否满足设计要求,也能作为完备的测试方法应用到科研生产测试工作中,与现有测试计算方法相比具有测试时间花销小、测试流程自动化程度高和测试结果精度高等优点。

Description

一种卫星接收机导航输出信号时序测试方法
技术领域
本发明属于信号时序测试技术领域,具体涉及一种卫星接收机导航输出信号时序测试方法。
背景技术
随着卫星导航技术和卫星导航接收机的普遍应用,对卫星导航信号定位精度的要求越来越高,而卫星接收机导航输出信号中同步脉冲信号和导航数据信号的时序关系对卫星导航信号定位精度至关重要。
传统的卫星接收机导航输出信号的时序测试工作主要通过示波器完成,由于输出数据信号的最小宽度在微秒级,而且同步脉冲信号的脉冲宽度在毫秒级,精度在微秒级,这对测试仪器精度和测试人员操作熟练度都有着很高的要求。测试人员需要用示波器通过测试接口反复调整信号波形显示,使得测试结果的可靠度降低且时间花销较大,在一定程度上延缓了科研生产任务。为了提高科研生产的自动化水平,针对卫星接收机信号时序测试工作,本发明摆脱了对示波器的依赖,开发了一种新的测试、计算方法:首先针对卫星接收机导航输出信号设计一种时域采样方法;然后设计一种量化编码方法;最后利用数字信号的特征和待测项目设计一种计算方法,并在测试计算机上实现了自动化计算、显示测试结果。
发明内容
本发明的目的在于:构建一种新的测试、计算方法,在测试时间花销小、测试流程自动化程度高的基础上,对卫星接收机导航输出信号中同步脉冲信号和导航数据信号的时序进行高精度的测试检验。
为了实现这一目的,本发明采取的技术方案是:
一种卫星接收机导航输出信号时序测试方法,本方法包括四步;
第一步、时域采样;
按照待测卫星接收机导航输出信号时序的指标值,设置信号采样设备的采样周期T。
第二步、量化编码;
根据卫星接收机导航输出信号的标称电平值,设置量化编码设备的逻辑判别门限Vth,量化编码设备把高于Vth的离散值判别为逻辑值1,低于Vth的离散值判别为逻辑值0,按照下式量化方法:
Figure GDA0003852209410000021
Figure GDA0003852209410000022
分别得到对应的两路数字信号
Figure GDA0003852209410000023
Figure GDA0003852209410000024
Figure GDA0003852209410000025
并保存为文本形式的数据;
第三步、计算方法
通过计算机读取第二步中保存文本形式的数据,统计同步脉冲采样点的逻辑值连续为1的点数M,结合第一步中时域采样的采样周期T,可以得到此路同步脉冲信号的脉冲宽度τh=MT;
统计同步脉冲采样点的逻辑值连续为1的片段之间连续的逻辑值0的点数N,可以得到数据更新率fs=NT;
在同步脉冲采样点的逻辑值由连续的1变为0时做一个标记k1,从标记k1开始统计逻辑值的点数,直到导航数据信号采样点的逻辑值由连续的0变为1时做另一标记k2,得到同步脉冲信号中脉冲与导航数据信号中数据包之间的间隔为τhf=(k2-k1)T;
在导航数据信号采样点的逻辑值由0变成1后统计逻辑值连续为1的点数N1,紧接着统计逻辑值连续为0的点数N2,得到此路导航数据信号的高电平宽度τf1=N1T,低电平宽度τf0=N2T,延续上述过程n次,然后取平均值得到此路导航数据信号的波特率
Figure GDA0003852209410000031
一种卫星接收机导航输出信号时序测试方法,所述第一步时域采样,为了有效测量产品的数据更新率fs、脉冲与数据间隔、波特率等参数,必须要确保截取的信号片段包含一个完整的脉冲和数据包,即截取的信号片段的宽度τcut需满足:
Figure GDA0003852209410000032
现实中信号采样设备由于其自身存储有限的原因,其采样点数最大值为Smax,结合采样周期可得到信号采样设备所能截取到的信号片段长度τcut满足下式:
τcut=TSmax
结合上述内容,可知信号采样设备的采样周期T的设置需要满足:
Figure GDA0003852209410000033
然后,用信号采样设备对同步脉冲信号h(t)和导航数据信号f(t)进行时域采样,分别得到两路时域离散信号
Figure GDA0003852209410000034
Figure GDA0003852209410000035
表示为:
Figure GDA0003852209410000036
Figure GDA0003852209410000037
本发明的有益效果为:本发明既可以用于检验卫星导航接收机输出时序是否满足设计要求等设计工作,也能作为完备的测试方法应用到科研生产测试工作中。本发明构建了一种新的测试、计算方法,与现有测试、计算方法相比具有测试时间花销小、测试流程自动化程度高和计算结果精度高等优点。
附图说明
图1为本发明的卫星接收机导航输出信号时序测试方法流程图。
具体实施方式
下面结合附图及具体实施例对本发明详细说明。
一种卫星接收机导航输出信号时序测试方法包括四步,第一步,对同步脉冲信号h(t)和导航数据信号f(t)进行时域采样,分别得到两路时域离散信号
Figure GDA0003852209410000041
Figure GDA0003852209410000042
表示为:
Figure GDA0003852209410000043
Figure GDA0003852209410000044
其中,T代表时域采样的采样周期,按照待测卫星接收机导航输出信号时序的具体指标值来设置的T取值,δ(·)为单位冲激函数且表达式为:
Figure GDA0003852209410000045
第二步,对两路时域离散信号
Figure GDA0003852209410000046
Figure GDA0003852209410000047
通过逻辑门限判别的方式进行量化编码,得到数字信号表示式为:
Figure GDA0003852209410000048
Figure GDA0003852209410000049
其中
Figure GDA00038522094100000410
Figure GDA00038522094100000411
表示整数域,u(·)为单位阶跃函数且表达式为:
Figure GDA00038522094100000412
第三步,设计一种算法,计算第二步中数字信号特征,得到的同步脉冲信号的脉冲宽度τh、导航数据信号的波特率B、同步脉冲信号中脉冲与导航数据信号中数据包之间的间隔τhf,数学表达式为:
τh=MT
Figure GDA0003852209410000051
τhf=(k2-k1)T
其中,M代表同步脉冲信号的逻辑值中连续为1的点数;逻辑值由0变成1后统计逻辑值连续为1的点数N1,紧接着统计逻辑值连续为0的点数N2,延续上述过程n次,Ni表示其中的第i(i=1,2,...n)次;k1代表在同步脉冲信号的逻辑值由连续的1变为0时的标记量,k2代表从标记k1到导航输出数据信号的逻辑值由连续的0变为1时的标记量。
具体实施例中,本方法分为四个步骤,如下:
第一步、时域采样
按照待测卫星接收机导航输出信号时序的指标值,设置信号采样设备的采样周期T。为了有效测量产品的数据更新率fs(同步脉冲信号频率)、脉冲与数据间隔、波特率等参数,必须要确保截取的信号片段包含一个完整的脉冲和数据包,即截取的信号片段的宽度τcut需满足:
Figure GDA0003852209410000052
现实中信号采样设备由于其自身存储有限等原因,其采样点数最大值为Smax,结合采样周期可得到信号采样设备所能截取到的信号片段长度τcut满足下式:
τcut=TSmax
结合上述内容,可知信号采样设备的采样周期T的设置需要满足:
Figure GDA0003852209410000053
然后,用信号采样设备对同步脉冲信号h(t)和导航数据信号f(t)进行时域采样,分别得到两路时域离散信号
Figure GDA0003852209410000061
Figure GDA0003852209410000062
表示为:
Figure GDA0003852209410000063
Figure GDA0003852209410000064
第二步、量化编码;
根据卫星接收机导航输出信号的标称电平值,设置量化编码设备的逻辑判别门限Vth,量化编码设备把高于Vth的离散值判别为逻辑值1,低于Vth的离散值判别为逻辑值0,按照下式量化方法:
Figure GDA0003852209410000065
Figure GDA0003852209410000066
分别得到对应的两路数字信号
Figure GDA0003852209410000067
Figure GDA0003852209410000068
Figure GDA0003852209410000069
并保存为文本形式的数据;
第三步、计算方法
通过计算机读取第二步中保存文本形式的数据,其数据结构如下表所示:
采样点 同步脉冲采样点逻辑值 导航数据信号采样点逻辑值
0 0 0
1 1 0
2 1 0
k 0 1
k+1 0 1
S<sub>max</sub> 1 1
由上表可知,统计同步脉冲采样点的逻辑值连续为1的点数M,结合第一步中时域采样的采样周期T,可以得到此路同步脉冲信号的脉冲宽度τh=MT;
由上表可知,统计同步脉冲采样点的逻辑值连续为1的片段之间连续的逻辑值0的点数N,可以得到数据更新率fs=NT;
由上表可知,在同步脉冲采样点的逻辑值由连续的1变为0时做一个标记k1,从标记k1开始统计逻辑值的点数,直到导航数据信号采样点的逻辑值由连续的0变为1时做另一标记k2,得到同步脉冲信号中脉冲与导航数据信号中数据包之间的间隔为τhf=(k2-k1)T;
由上表可知,在导航数据信号采样点的逻辑值由0变成1后统计逻辑值连续为1的点数N1,紧接着统计逻辑值连续为0的点数N2,得到此路导航数据信号的高电平宽度τf1=N1T,低电平宽度τf0=N2T,延续上述过程n次,然后取平均值得到此路导航数据信号的波特率
Figure GDA0003852209410000071

Claims (2)

1.一种卫星接收机导航输出信号时序测试方法,其特征在于,本方法包括四步;
第一步、时域采样;
按照待测卫星接收机导航输出信号时序的指标值,设置信号采样设备的采样周期T;
第二步、量化编码;
根据卫星接收机导航输出信号的标称电平值,设置量化编码设备的逻辑判别门限Vth,量化编码设备把高于Vth的离散值判别为逻辑值1,低于Vth的离散值判别为逻辑值0,按照下式量化方法:
Figure FDA0003852209400000011
分别得到对应的两路数字信号
Figure FDA0003852209400000012
Figure FDA0003852209400000013
并保存为文本形式的数据;
Figure FDA0003852209400000014
为同步脉冲的时域离散信号,
Figure FDA0003852209400000015
为导航数据的时域离散信号;
第三步、计算方法
通过计算机读取第二步中保存文本形式的数据,统计同步脉冲采样点的逻辑值连续为1的点数M,结合第一步中时域采样的采样周期T,可以得到此路同步脉冲信号的脉冲宽度τh=MT;
统计同步脉冲采样点的逻辑值连续为1的片段之间连续的逻辑值0的点数N,可以得到数据更新率fs=NT;
在同步脉冲采样点的逻辑值由连续的1变为0时做一个标记k1,从标记k1开始统计逻辑值的点数,直到导航数据信号采样点的逻辑值由连续的0变为1时做另一标记k2,得到同步脉冲信号中脉冲与导航数据信号中数据包之间的间隔为τhf=(k2-k1)T;
在导航数据信号采样点的逻辑值由0变成1后统计逻辑值连续为1的点数N1,紧接着统计逻辑值连续为0的点数N2,得到此路导航数据信号的高电平宽度τf1=N1T,低电平宽度τf0=N2T,延续上述过程n次,然后取平均值得到此路导航数据信号的波特率
Figure FDA0003852209400000021
2.如权利要求1所述的一种卫星接收机导航输出信号时序测试方法,其特征在于:所述第一步时域采样,为了有效测量产品的数据更新率fs、脉冲与数据间隔、波特率参数,必须要确保截取的信号片段包含一个完整的脉冲和数据包,即截取的信号片段的宽度τcut需满足:
Figure FDA0003852209400000022
现实中信号采样设备由于其自身存储有限的原因,其采样点数最大值为Smax,结合采样周期可得到信号采样设备所能截取到的信号片段长度τcut满足下式:
τcut=TSmax
结合上述内容,可知信号采样设备的采样周期T的设置需要满足:
Figure FDA0003852209400000023
然后,用信号采样设备对同步脉冲信号h(t)和导航数据信号f(t)进行时域采样,分别得到两路时域离散信号
Figure FDA0003852209400000024
Figure FDA0003852209400000025
表示为:
Figure FDA0003852209400000026
CN201810700852.XA 2018-06-29 2018-06-29 一种卫星接收机导航输出信号时序测试方法 Active CN110658536B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810700852.XA CN110658536B (zh) 2018-06-29 2018-06-29 一种卫星接收机导航输出信号时序测试方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810700852.XA CN110658536B (zh) 2018-06-29 2018-06-29 一种卫星接收机导航输出信号时序测试方法

Publications (2)

Publication Number Publication Date
CN110658536A CN110658536A (zh) 2020-01-07
CN110658536B true CN110658536B (zh) 2023-02-10

Family

ID=69026962

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810700852.XA Active CN110658536B (zh) 2018-06-29 2018-06-29 一种卫星接收机导航输出信号时序测试方法

Country Status (1)

Country Link
CN (1) CN110658536B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005003807A1 (en) * 2003-07-03 2005-01-13 Nemerix Sa Gps receiver with fast acquisition time
CN1590961A (zh) * 2003-09-04 2005-03-09 伊莱德公司 组合导航卫星接收机和通信设备
CN101908895A (zh) * 2009-06-08 2010-12-08 泰纳股份公司 接收机和接收信号的方法
CN103308930A (zh) * 2013-05-24 2013-09-18 北京东方计量测试研究所 卫星导航信号模拟器伪距精度测量方法
CN104035104A (zh) * 2014-06-26 2014-09-10 桂林电子科技大学 一种抗欺骗伪卫星系统及信号产生方法
CN106656451A (zh) * 2016-12-29 2017-05-10 中国科学院西安光学精密机械研究所 一种基于卫星授时系统的守时、授时精度测试装置及方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005003807A1 (en) * 2003-07-03 2005-01-13 Nemerix Sa Gps receiver with fast acquisition time
CN1590961A (zh) * 2003-09-04 2005-03-09 伊莱德公司 组合导航卫星接收机和通信设备
CN101908895A (zh) * 2009-06-08 2010-12-08 泰纳股份公司 接收机和接收信号的方法
CN103308930A (zh) * 2013-05-24 2013-09-18 北京东方计量测试研究所 卫星导航信号模拟器伪距精度测量方法
CN104035104A (zh) * 2014-06-26 2014-09-10 桂林电子科技大学 一种抗欺骗伪卫星系统及信号产生方法
CN106656451A (zh) * 2016-12-29 2017-05-10 中国科学院西安光学精密机械研究所 一种基于卫星授时系统的守时、授时精度测试装置及方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
刘进 等.低信噪比下卫星通信信号的非线性增强算法.《舰船电子对抗》.2018,第41卷(第3期),第79-83页. *

Also Published As

Publication number Publication date
CN110658536A (zh) 2020-01-07

Similar Documents

Publication Publication Date Title
CN100575965C (zh) 一种单片机测量频率的方法
CN106680782B (zh) 基于脉冲时延匹配的雷达重频分选方法
CN109407500B (zh) 一种基于fpga的时间间隔测量方法
CN104502684A (zh) 一种全数字化峰值到达时刻鉴别方法
CN104393981B (zh) 一种多路测量数据并行的时间标记方法及系统
CN110554364B (zh) 一种雷达信号脉冲到达时间测量精度的测试系统及方法
WO2022170933A1 (zh) 激光测距误差校正方法、装置、电子设备和存储介质
CN110703583A (zh) 基于soc的多通道高精度大量程时间数字转换器
CN106679694A (zh) 塔康信标模拟器空‑空应答延时测量精度标校装置及方法
CN114355174B (zh) 一种进位链延时测量校准方法及装置
CN104931785B (zh) 一种基于众数法的pcb特性阻抗的测试方法
CN100465944C (zh) 一种补偿示波器时基抖动的方法
CN209486181U (zh) 一种准确测量频率的系统
CN110658536B (zh) 一种卫星接收机导航输出信号时序测试方法
CN106855990B (zh) 核电机组仪表通道测量误差论证方法
CN102193029A (zh) 非常规采样时间短期频率稳定度测量方法
CN106093905B (zh) 基于apFFT相位差的精密测量雷达频偏分类测频方法
CN104111481A (zh) 同步时钟相位差测量系统和方法
CN109814135B (zh) 基于gnss设备的时间对准方法
CN109145258B (zh) 基于非线性拟合的威布尔分布参数置信区间估计方法
CN103529687B (zh) 脉冲时间间隔测量装置
CN107908133A (zh) 一种频率自适应采集分发方法
CN211180593U (zh) 时频检测设备
CN103308083B (zh) 一种基于能量分布的布里渊谱识别方法
CN109257033B (zh) 高精度步进延迟系统

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