CN109856646A - 一种高精度卫星导航信号间相对时延控制方法 - Google Patents
一种高精度卫星导航信号间相对时延控制方法 Download PDFInfo
- Publication number
- CN109856646A CN109856646A CN201811586780.7A CN201811586780A CN109856646A CN 109856646 A CN109856646 A CN 109856646A CN 201811586780 A CN201811586780 A CN 201811586780A CN 109856646 A CN109856646 A CN 109856646A
- Authority
- CN
- China
- Prior art keywords
- time delay
- satellite navigation
- pseudo noise
- noise code
- navigation signals
- 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
Links
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明涉及一种高精度卫星导航信号间相对时延控制方法,包括步骤如下:根据两路卫星导航信号生成的时钟频率,确定初步调整阈值;在这两个时钟频率下生成两路伪随机码,根据伪随机码时延差测量结果,通过调整启动信号对生成的伪随机码进行延迟控制;根据调整后时延差测量结果,计算滤波冲激响应以及时延精细调整所需冲激响应参数;将所述滤波冲激响应参数与第一路卫星导航信号卷积,将所述时延精细调整所需冲激响应参数与第二路卫星导航信号卷积;再对所述两路卫星导航信号进行数模转换,即生成所述相对时延精细调整的卫星导航信号。
Description
技术领域
本发明涉及卫星导航信号生成技术,特别涉及一种高精度卫星导航信号间相对时延同步控制方法。
背景技术
导航信号由多个频点组成,不同频点的基带信号元素不同要求必须在不同时钟域下生成信号,从而导致不同时钟域下信号存在时延差。如何通过相对时延控制缩小该时延差是信号生成必须解决的问题。
当前,一般采用的办法是使用第三方启动信号,但是这样存在亚稳态问题,最多只能将相对时延控制在一个时钟周期之内,不满足同步性要求。
发明内容
本发明所要解决的技术问题是:本发明结合卫星导航信号生成的特点,提出一种高精度卫星导航信号相对时延控制方法,结合了卫星导航信号特有的性质,以此提高卫星发射端信号质量,保证高精度定位领域应用性能。
本发明采用的技术方案是:一种高精度卫星导航信号间相对时延控制方法,包括步骤如下:
步骤1、根据两路卫星导航信号生成的时钟频率f1、f2,确定初步调整阈值;
步骤2、在上述两个时钟频率f1、f2下生成两路伪随机码,根据伪随机码时延差测量结果,通过调整启动信号对生成的伪随机码进行延迟控制,获得伪随机码延迟控制量;
步骤3、根据调整后时延差测量结果,计算滤波冲激响应参数以及时延精细调整所需冲激响应参数;
步骤4、将所述滤波冲激响应参数与第一路卫星导航信号卷积,将所述时延精细调整所需冲激响应参数与第二路卫星导航信号卷积;
步骤5、对步骤4中获得的两路卫星导航信号进行数模转换,生成精细同步的卫星导航信号。
所述初步调整阈值是根据两路卫星导航信号生成的时钟频率的最大公约数确定的,方法如下:
两路卫星导航信号生成的两个时钟频率f1、f2的最大公约数频率为fy,其周期为Ty,Ty=M·T1=N·T2;
其中,M和N为正整数;T1为时钟频率f1对应的周期,T2为时钟频率f2对应的周期;f1<f2;
计算
[x(1),ind(1)]=min{abs[1T1-[0 1 2 … N]·T2]},
[x(2),ind(2)]=min{abs[2T1-[0 1 2 … N]·T2]},
[x(M-1),ind(M-1)]=min{abs[(M-1)T1-[0 1 2 … N]·T2]};
其中,min{}表示取集合中的绝对值最小值及其索引,abs[]表示取绝对值;x()表示取得的绝对值最小值,ind()表示该绝对值最小值的索引;
计算[y,indy]=min{x(1),x(2),x(3),...x(M-1)},则初步调整阈值为y/2。
所述伪随机码或卫星导航信号在监测到启动信号第一个上升沿时开始生成,调整启动信号指延迟一个或者多个时钟周期。
所述步骤2中,延迟控制量的计算方法如下:
(2.1)在f1、f2两个时钟域下生成两路伪随机码,以一个第三方钟为测量时刻测量两路伪随机码时延差,通过调整启动信号使f1时钟域下的伪随机码超前f2时钟域下的伪随机码;
(2.2)记录f1时钟域下的伪随机码与f2时钟域下的伪随机码的时延值t1,判断t1,是否小于y/2,小于则调整结束,不小于则进入下一步;
(2.3)延迟f2时钟域下的伪随机码一个时钟周期,记录f1时钟域下的伪随机码与f2时钟域下的伪随机码的时延值t2,判断t2是否小于y/2,小于则调整结束,不小于则进入下一步;
(2.4)重复(2.3),直至f2时钟域下的伪随机码延迟了N次;
(2.5)延迟f1时钟域下伪随机码一个时钟周期,重复步骤(2.2)~步骤(2.4);
(2.6)重复步骤(2.5),直至f1时钟域下的伪随机码延迟了M次;
(2.7)记录最终小于y/2的f1时钟域下的伪随机码与f2时钟域下的伪随机码的时延值τ,得到延迟控制量。
所述步骤3中,滤波冲激响应参数以及时延精细调整所需冲激响应参数的计算方法如下:
f1时钟域下进行时延调整的冲激响应,得到滤波冲激响应参数h1(n):
其中,K≥50,K=kM,k为正整数,n∈{0,1,2,...,2K},fc1≤f1/2;fc1为第一路卫星导航信号的单边带宽;
f2时钟域下进行时延精确调整的冲激响应,得到时延精细调整所需冲激响应参数h2(n-τ):
其中,fc2为第二路卫星导航信号的单边带宽。
所述步骤4中的卷积所用结构为FIR结构,在各自时钟域实现。
本发明与现有技术相比具有如下有益效果:
通过采用本发明的技术方案,通过进行在线精细测量,可实时获取两路导航信号精确时延;通过延迟控制和时延精确调整,结合闭环和开环环节,将两个时钟域调制信号相对时延可以控制在较小范围,可小于0.1ns。这在高精度双频精密单点定位中,可以有效提升定位精度和收敛时间,使信号可直接应用于无人驾驶等实时应用领域。
附图说明
图1是根据本发明的卫星导航信号间相对时延控制方法的示意图。
图2是根据本发明的实施例计算得到的冲激响应以及时延精细调整所需冲激响应参数。
具体实施方式
下面结合附图和具体实施例对本发明作进一步详细的描述:
如图1所示,是根据本发明的实施例的方法实施框图。一种高精度卫星导航信号间相对时延控制方法,包括步骤如下:
步骤1、根据所要生成的两个时钟频率f1、f2确定初步调整阈值;
初步调整阈值是根据两路卫星导航信号生成的时钟频率的最大公约数确定的,其中,
两路卫星导航信号生成的时钟频率为
f1,其周期为T1,
f2,其周期为T2,且f1<f2
两个频率的最大公约数频率为fy,其周期为Ty。因此,Ty=M·T1=N·T2,其中,M和N为正整数。
计算
[x(1),ind(1)]=min{abs[1T1-[0 1 2 … N]·T2]}
[x(2),ind(2)]=min{abs[2T1-[0 1 2 … N]·T2]}
……
[x(M-1),ind(M-1)]=min{abs[(M-1)T1-[0 1 2 … N]·T2]}其中,min{}表示取集合中的绝对值最小值及其索引,abs[]表示取绝对值。x()表示取得的绝对值最小值,ind()表示该绝对值最小值的索引;
计算
[y,indy]=min{x(1),x(2),x(3),...x(M-1)}
y/2即为初步调整阈值阈值。
步骤2、在这两个时钟频率下生成两路伪随机码,使用第三方时钟对两路伪随机码进行相位测量以确定时延差;
步骤3、调整启动信号使f1域下伪码超前f2域下伪码;
所述启动信号就是一个该时钟域下的高电平脉冲,持续一个时钟周期,其中,伪随机码或卫星导航信号在监测到启动信号第一个上升沿时开始生成,调整启动信号即是延迟一个或者多个时钟周期。
步骤4、记录f1域下伪码-f2域下伪随机码时延差值为t1,判断是否小于y/2,小于则调整结束,不小于则进入下一步;
步骤5、延迟f2时钟域下伪随机码一个时钟周期,记录f1域下伪码-f2域下伪码时延值为t2,判断是否小于y/2,小于则调整结束,不小于则进入下一步;
步骤6、返回步骤5,直至f2时钟域下伪随机码延迟了N次;
步骤7、延迟f1时钟域下伪随机码一个时钟周期,重复步骤4~步骤6;
步骤8、返回步骤7,直至f1时钟域下伪随机码延迟了M次;
步骤9、记录最终小于y/2的f1域下伪码-f2域下伪随机码时延值为τ;
步骤10、在启动信号后进行两路卫星导航信号生成;
步骤11、实时计算滤波冲激响应参数以及时延精细调整所需冲激响应参数,按照数字计算有效位数进行量化和截断;
滤波冲激响应参数以及时延精细调整所需冲激响应参数是通过实时计算得到的,其中,
f1时钟域下进行时延调整的冲激响应,得到滤波冲激响应参数h1(n):
其中,K≥50,K=kM,k为正整数,n∈{0,1,2,...,2K},fc1≤f1/2,fc1为第一路卫星导航信号的单边带宽;
f2时钟域下进行时延精确调整的冲激响应,得到时延精细调整所需冲激响应参数h2(n-τ):
其中,fc2≤f2/2,fc2为第二路卫星导航信号的单边带宽。
应用时需要按照数字计算有效位数进行量化和截断使用。
步骤12、将所述滤波冲激响应参数与第一路卫星导航信号卷积,将所述时延精细调整所需冲激响应参数与第二路卫星导航信号卷积;
卷积所用结构为FIR(Finite Impulse Response有限冲激响应)结构,在各自时钟域实现。
步骤13、将卫星导航信号数模转换输出。
实施例:
本发明提供了一种卫星导航信号间相对时延控制方法,包括步骤如下:
设两个信号生成时钟频率为122.76MHz(周期约为T1=8.146ns)和163.68MHz(周期约为T2=6.109ns),则其最大公约数频率为40.92MHz,则Ty=M·T1=N·T2=3·T1=4·T2=24.44ns。
计算
[x(1),ind(1)]=min{abs[1T1-[0 1 2 … N]·T2]}
[x(2),ind(2)]=min{abs[2T1-[0 1 2 … N]·T2]}
……
[x(M-1),ind(M-1)]=min{abs[(M-1)T1-[0 1 2 … N]·T2]}
有表:
其中第一行绝对值最小值为2.037,第二行绝对值最小值为-2.035,因此该步调整的阈值为2.035/2=1.017ns,即最大精度为1.017ns。
在这两个时钟频率下生成两路伪随机码,使用第三方时钟对两路伪随机码进行相位测量以确定时延差;
调整启动信号使两路伪随机码时延差小于所算阈值;
实时计算滤波冲激响应以及时延精细调整所需冲激响应参数,按照数字计算有效位数进行量化和截断;
将所述滤波冲激响应参数与第一路卫星导航信号卷积,将所述时延精细调整所需冲激响应参数与第二路卫星导航信号卷积;
最后将卫星导航信号数模转换输出。
图2是根据本发明的实施例计算得到的冲激响应以及时延精细调整所需冲激响应参数。
以上所述,仅为本发明最佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。
本发明说明书中未作详细描述的内容属于本领域技术人员的公知技术。
Claims (6)
1.一种高精度卫星导航信号间相对时延控制方法,其特征在于,包括步骤如下:
步骤1、根据两路卫星导航信号生成的时钟频率f1、f2,确定初步调整阈值;
步骤2、在上述两个时钟频率f1、f2下生成两路伪随机码,根据伪随机码时延差测量结果,通过调整启动信号对生成的伪随机码进行延迟控制,获得伪随机码延迟控制量;
步骤3、根据调整后时延差测量结果,计算滤波冲激响应参数以及时延精细调整所需冲激响应参数;
步骤4、将所述滤波冲激响应参数与第一路卫星导航信号卷积,将所述时延精细调整所需冲激响应参数与第二路卫星导航信号卷积;
步骤5、对步骤4中获得的两路卫星导航信号进行数模转换,生成精细同步的卫星导航信号。
2.根据权利要求1所述的一种高精度卫星导航信号间相对时延控制方法,其特征在于,所述初步调整阈值是根据两路卫星导航信号生成的时钟频率的最大公约数确定的,方法如下:
两路卫星导航信号生成的两个时钟频率f1、f2的最大公约数频率为fy,其周期为Ty,Ty=M·T1=N·T2;
其中,M和N为正整数;T1为时钟频率f1对应的周期,T2为时钟频率f2对应的周期;f1<f2;
计算
[x(1),ind(1)]=min{abs[1T1-[012…N]·T2]},
[x(2),ind(2)]=min{abs[2T1-[012…N]·T2]},
[x(M-1),ind(M-1)]=min{abs[(M-1)T1-[012…N]·T2]};
其中,min{ }表示取集合中的绝对值最小值及其索引,abs[]表示取绝对值;x()表示取得的绝对值最小值,ind()表示该绝对值最小值的索引;
计算[y,indy]=min{x(1),x(2),x(3),...x(M-1)},则初步调整阈值为y/2。
3.根据权利要求1或2所述的一种高精度卫星导航信号间相对时延控制方法,其特征在于,所述伪随机码或卫星导航信号在监测到启动信号第一个上升沿时开始生成,调整启动信号指延迟一个或者多个时钟周期。
4.根据权利要求3所述的一种高精度卫星导航信号间相对时延控制方法,其特征在于,所述步骤2中,延迟控制量的计算方法如下:
(2.1)在f1、f2两个时钟域下生成两路伪随机码,以一个第三方钟为测量时刻测量两路伪随机码时延差,通过调整启动信号使f1时钟域下的伪随机码超前f2时钟域下的伪随机码;
(2.2)记录f1时钟域下的伪随机码与f2时钟域下的伪随机码的时延值t1,判断t1,是否小于y/2,小于则调整结束,不小于则进入下一步;
(2.3)延迟f2时钟域下的伪随机码一个时钟周期,记录f1时钟域下的伪随机码与f2时钟域下的伪随机码的时延值t2,判断t2是否小于y/2,小于则调整结束,不小于则进入下一步;
(2.4)重复(2.3),直至f2时钟域下的伪随机码延迟了N次;
(2.5)延迟f1时钟域下伪随机码一个时钟周期,重复步骤(2.2)~步骤(2.4);
(2.6)重复步骤(2.5),直至f1时钟域下的伪随机码延迟了M次;
(2.7)记录最终小于y/2的f1时钟域下的伪随机码与f2时钟域下的伪随机码的时延值τ,得到延迟控制量。
5.根据权利要求4所述的一种高精度卫星导航信号间相对时延控制方法,其特征在于:所述步骤3中,滤波冲激响应参数以及时延精细调整所需冲激响应参数的计算方法如下:
f1时钟域下进行时延调整的冲激响应,得到滤波冲激响应参数h1(n):
其中,K≥50,K=kM,k为正整数,n∈{0,1,2,...,2K},fc1≤f1/2;fc1为第一路卫星导航信号的单边带宽;
f2时钟域下进行时延精确调整的冲激响应,得到时延精细调整所需冲激响应参数h2(n-τ):
其中,fc2为第二路卫星导航信号的单边带宽。
6.根据权利要求5所述的一种高精度卫星导航信号间相对时延控制方法,其特征在于,所述步骤4中的卷积所用结构为FIR结构,在各自时钟域实现。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811586780.7A CN109856646B (zh) | 2018-12-24 | 2018-12-24 | 一种高精度卫星导航信号间相对时延控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811586780.7A CN109856646B (zh) | 2018-12-24 | 2018-12-24 | 一种高精度卫星导航信号间相对时延控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109856646A true CN109856646A (zh) | 2019-06-07 |
CN109856646B CN109856646B (zh) | 2021-03-26 |
Family
ID=66892227
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811586780.7A Active CN109856646B (zh) | 2018-12-24 | 2018-12-24 | 一种高精度卫星导航信号间相对时延控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109856646B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111045042A (zh) * | 2019-12-20 | 2020-04-21 | 西安空间无线电技术研究所 | 一种基于“云-端”架构的ppp-rtk增强方法与系统 |
CN113703021A (zh) * | 2021-07-29 | 2021-11-26 | 西安空间无线电技术研究所 | 一种基于码伪距的秒级实时高精度定位方法与系统 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1148344A1 (en) * | 2000-04-20 | 2001-10-24 | Lucent Technologies Inc. | Positioning of a wireless terminal with satellite positioning signals or base station signals |
CN204256171U (zh) * | 2014-12-05 | 2015-04-08 | 中国电子科技集团公司第五十四研究所 | 卫星导航源多频时延测量基带处理电路 |
CN105607096A (zh) * | 2015-08-31 | 2016-05-25 | 中国电子科技集团公司第三十六研究所 | 一种双星时差频差定位方法和定位装置 |
CN106371110A (zh) * | 2016-08-16 | 2017-02-01 | 上海航天测控通信研究所 | 一种gnss‑r双补时延干涉处理系统及方法 |
CN107271955A (zh) * | 2017-07-25 | 2017-10-20 | 西安电子科技大学 | 一种宽带线性调频信号的时差和尺度差估计方法 |
CN107272033A (zh) * | 2017-05-27 | 2017-10-20 | 金华航大北斗应用技术有限公司 | 一种北斗meo卫星机载反射信号的随机时延建模方法 |
CN109039421A (zh) * | 2018-06-28 | 2018-12-18 | 上海卫星工程研究所 | 用于多载荷数据相对传输时延测试的方法 |
-
2018
- 2018-12-24 CN CN201811586780.7A patent/CN109856646B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1148344A1 (en) * | 2000-04-20 | 2001-10-24 | Lucent Technologies Inc. | Positioning of a wireless terminal with satellite positioning signals or base station signals |
CN204256171U (zh) * | 2014-12-05 | 2015-04-08 | 中国电子科技集团公司第五十四研究所 | 卫星导航源多频时延测量基带处理电路 |
CN105607096A (zh) * | 2015-08-31 | 2016-05-25 | 中国电子科技集团公司第三十六研究所 | 一种双星时差频差定位方法和定位装置 |
CN106371110A (zh) * | 2016-08-16 | 2017-02-01 | 上海航天测控通信研究所 | 一种gnss‑r双补时延干涉处理系统及方法 |
CN107272033A (zh) * | 2017-05-27 | 2017-10-20 | 金华航大北斗应用技术有限公司 | 一种北斗meo卫星机载反射信号的随机时延建模方法 |
CN107271955A (zh) * | 2017-07-25 | 2017-10-20 | 西安电子科技大学 | 一种宽带线性调频信号的时差和尺度差估计方法 |
CN109039421A (zh) * | 2018-06-28 | 2018-12-18 | 上海卫星工程研究所 | 用于多载荷数据相对传输时延测试的方法 |
Non-Patent Citations (1)
Title |
---|
刘潇 等: "基于时间序列的导航信号群时延特性检测算法", 《第九届中国卫星导航学术年会》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111045042A (zh) * | 2019-12-20 | 2020-04-21 | 西安空间无线电技术研究所 | 一种基于“云-端”架构的ppp-rtk增强方法与系统 |
CN113703021A (zh) * | 2021-07-29 | 2021-11-26 | 西安空间无线电技术研究所 | 一种基于码伪距的秒级实时高精度定位方法与系统 |
CN113703021B (zh) * | 2021-07-29 | 2023-09-29 | 西安空间无线电技术研究所 | 一种基于码伪距的秒级实时高精度定位方法与系统 |
Also Published As
Publication number | Publication date |
---|---|
CN109856646B (zh) | 2021-03-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7042958B2 (en) | Digital time alignment in a polar modulator | |
CN102223322B (zh) | 一种频率偏差的估计方法及装置 | |
TWI300661B (en) | Apparatus for digital timing adjustment of rake receiver and method thereof | |
CN104300969B (zh) | 一种基于全数字锁相环的高精度同步时钟实现方法 | |
CN109856646A (zh) | 一种高精度卫星导航信号间相对时延控制方法 | |
US10148345B2 (en) | Accurate timing distribution by high-frequency radio | |
CN204631463U (zh) | 一种辐射源监测定位的gps时钟同步系统 | |
CN103702409B (zh) | 一种td‑lte网络优化设备的同步检测方法 | |
CN108768577B (zh) | 一种基于ptp时间同步信号的通信网络授时方法及系统 | |
CN102981551B (zh) | 一种实时时钟温度补偿系统及方法 | |
CN111580380B (zh) | 一种提高gnss同步授时精度方法 | |
CN112087407B (zh) | 基于采样率的动态调整的组合多普勒估计方法 | |
CN107505832A (zh) | 一种高精度授时系统 | |
CN104135326B (zh) | 适用于edge系统8psk信号的射频测试装置与方法 | |
CN107064918B (zh) | 一种高精度雷达测距方法 | |
CN109856615B (zh) | 一种基于css技术的测距方法及系统 | |
CN104052524A (zh) | 一种gps辅助的快速跳频同步方法 | |
CN106169949B (zh) | 一种基带信号位同步时钟宽频自适应提取装置及方法 | |
CN107948111B (zh) | Ofdm系统的采样频偏矫正方法 | |
CN106922015A (zh) | 无线通信设备及其频率同步方法 | |
CN108401582B (zh) | 一种适应tdma/fdma体制的geo卫星移动通信系统初始频率同步方法 | |
CN110515102A (zh) | 一种基于卫星共视观测驾驭晶振的ns级时间复现方法 | |
CN109120260A (zh) | 一种基于asic-tdc的时钟模块高精度鉴相系统及方法 | |
CN110971325B (zh) | 基于位同步的时间传递方法 | |
CN103441978A (zh) | 一种基于ofdm调制的电力线载波通信方法 |
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 |