CN104049528B - 一种北斗授时方法及卫星导航接收机 - Google Patents
一种北斗授时方法及卫星导航接收机 Download PDFInfo
- Publication number
- CN104049528B CN104049528B CN201410307499.0A CN201410307499A CN104049528B CN 104049528 B CN104049528 B CN 104049528B CN 201410307499 A CN201410307499 A CN 201410307499A CN 104049528 B CN104049528 B CN 104049528B
- Authority
- CN
- China
- Prior art keywords
- receiver
- delta
- centerdot
- satellite
- rho
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000005070 sampling Methods 0.000 claims abstract description 53
- 238000004364 calculation method Methods 0.000 claims description 20
- 238000009499 grossing Methods 0.000 claims description 14
- 230000009471 action Effects 0.000 claims description 12
- 230000005540 biological transmission Effects 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 6
- 230000008054 signal transmission Effects 0.000 claims description 5
- 230000001360 synchronised effect Effects 0.000 abstract description 2
- 239000000243 solution Substances 0.000 description 14
- 238000004891 communication Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 239000013078 crystal Substances 0.000 description 2
- 238000010276 construction Methods 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000012895 dilution Substances 0.000 description 1
- 238000010790 dilution Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
一种北斗授时方法及卫星导航接收机;所述方法包括:产生采样脉冲;利用采样脉冲采集所接收的卫星信号中的观测信息;所述观测信息包括卫星信号发射时刻、多普勒频移、载波相位;利用所述卫星信号发射时刻和载波相位组成的伪距,构成接收机位置解算方程组;在所述接收机位置解算方程组中添加一个高程方程,解算得到接收机的钟差;利用所述观测信息中的多普勒频移构成接收机速度解算方程组,在所述接收机速度解算方程组中添加一个高程方程,解算得到接收机的钟漂;根据所述钟差和钟漂调节所述采样脉冲的相位和频率,使采样脉冲与标准时间同步。本发明能够提高卫星导航定位系统的授时精度,拓宽授时范围。
Description
技术领域
本发明涉及卫星导航领域,尤其涉及一种北斗授时方法及卫星导航接收机。
背景技术
高精度时间基准已经成为通信、电力、广播电视、安防监控、工业控制等领域的基础保障平台之一。卫星导航定位系统可提供高精度、全天候的定位和授时服务。授时性能优异;高精度、低成本、安全可靠、覆盖范围广。精确授时不仅在相对论研究、地球动力学研究等基础研究领域有重要的作用,而且在航空航天、深空通讯、信息高速公路、导航通信、电力传输等国防和国民经济建设中也有广泛的应用。卫星导航接收机采用授时算法能够使得接收机的时间与标准时间同步。
但是目前有些区域接收卫星信号的情况并不好,会导致定位系统的DOP(Dilutionof Precision,精度因子)值比较差,而影响授时精度。甚至当只接收3颗卫星时,可能会导致不能进行定位而不能完成授时。
发明内容
本发明要解决的技术问题是如何提高卫星导航定位系统的授时精度,拓宽授时范围。
为了解决上述问题,本发明提供了一种北斗授时方法,包括:
产生采样脉冲;
利用采样脉冲采集所接收的卫星信号中的观测信息;所述观测信息包括卫星信号发射时刻、多普勒频移、载波相位;
利用所述卫星信号发射时刻和载波相位组成的伪距,构成接收机位置解算方程组;在所述接收机位置解算方程组中添加一个高程方程,解算得到接收机的钟差;
利用所述观测信息中的多普勒频移构成接收机速度解算方程组,在所述接收机速度解算方程组中添加一个高程方程,解算得到接收机的钟漂;
根据所述钟差和钟漂调节所述采样脉冲的相位和频率,使采样脉冲与标准时间同步。
可选地,所述利用卫星信号发射时刻和载波相位组成的伪距,构成接收机位置解算方程组的步骤包括:
利用卫星信号接收时刻减去卫星信号发射时刻,得到原始伪距;
通过下式计算平滑后的伪距:
ρs,i=α×ρCA,i+(1-α)×(ρs,i-1+Δρφ,i)
其中,ρs,i是平滑后的伪距;α为平滑系数;ρCA,i是原始伪距;ρs,i-1是前一次平滑后的伪距;Δρφ,i是伪距变化量;
所述伪距变化量为:
λ是所接收的卫星信号的波长;φi是本次采样得到的载波相位值;φi-1是前一次采样得到的载波相位值;fIF是所接收的卫星信号的载波中频频率;为接收机钟漂;X是采样时刻锁存的载波相位计数值;fL1是所接收的卫星信号的载波射频频率;
利用平滑后的伪距,构成解算接收机位置的线性方程组如下:
ρi为第i颗卫星的平滑后的伪距,i的取值从1到n,n为接收机所能连接的卫星的颗数;(xu,yu,zu)为接收机的位置坐标;(xi,yi,zi)为第i颗卫星的位置坐标;tu为接收机钟差;ri为接收机到第i颗卫星的距离;c表示光速;带^标示的表示该变量为系统估计值。
可选地,利用观测信息中的多普勒频移所构成的接收机速度的解算方程组为如下的线性方程组:
其中,
fd是经过卫星钟漂修正的多普勒频移观测值;fTj是经过卫星钟漂修正的卫星发射信号的射频频率;()为接收机的速度;(vxj,vyj,vzj)为卫星速度;为接收机的钟漂。
可选地,解算得到接收机的钟差的步骤包括:
采用最小二乘法通过下式得到接收机的钟差:
ΔX=(HTH)-1HTΔρ
接收机钟差为
通过最小二乘法得到接收机的钟漂,公式如下:
g=(HTH)-1Hd
其中,
可选地,在所述接收机位置解算方程组中所添加的高程方程为:
Δρ=axΔxu+ayΔyu+azΔzu
其中,
CGS2000_A=6378137.
e2=CGS2000_F*(2-CGS2000_F)
CGS2000_F=(1./298.257222101)
lat为接收机的纬度,在每次迭代过程中由接收机近似位置经过转换获得;
alt为接收机的高度;
增加的H中的一行为[ax ay az 0];
在所述接收机速度解算方程组中所添加的高程方程为:
dj=0
增加的H中的一行为[ax ay az 0]。
本发明还提供了一种卫星导航接收机,包括:
脉冲生成模块,用于产生采样脉冲;
采集单元,用于利用所述采样脉冲采集所接收的卫星信号中的观测信息;所述观测信息包括卫星信号发射时刻、多普勒频移、载波相位;
钟差解算单元,用于利用所述卫星信号发射时刻和载波相位组成的伪距,构成接收机位置解算方程组;在所述接收机位置解算方程组中添加一个高程方程,解算得到接收机的钟差;
钟漂解算单元,用于利用所述观测信息中的多普勒频移构成接收机速度解算方程组,在所述接收机速度解算方程组中添加一个高程方程,解算得到接收机的钟漂;
调整单元,用于根据所述钟差和钟漂调节所述采样脉冲的相位和频率,使采样脉冲与标准时间同步。
可选地,所述钟差解算单元利用卫星信号发射时刻和载波相位组成的伪距,构成接收机位置解算方程组是指:
所述钟差解算单元利用卫星信号接收时刻减去卫星信号发射时刻,得到原始伪距;通过下式计算平滑后的伪距:
ρs,i=α×ρCA,i+(1-α)×(ρs,i-1+Δρφ,i)
其中,ρs,i是平滑后的伪距;α为平滑系数;ρCA,i是原始伪距;ρs,i-1是前一次平滑后的伪距;Δρφ,i是伪距变化量;
所述伪距变化量为:
λ是所接收的卫星信号的波长;φi是本次采样得到的载波相位值;φi-1是前一次采样得到的载波相位值;fIF是所接收的卫星信号的载波中频频率;为接收机钟漂;X是采样时刻锁存的载波相位计数值;fL1是所接收的卫星信号的载波射频频率;
利用平滑后的伪距,构成解算接收机位置的线性方程组如下:
ρi为第i颗卫星的平滑后的伪距,i的取值从1到n,n为接收机所能连接的卫星的颗数;(xu,yu,zu)为接收机的位置坐标;(xi,yi,zi)为第i颗卫星的位置坐标;tu为接收机钟差;ri为接收机到第i颗卫星的距离;c表示光速;带^标示的表示该变量为系统估计值。
可选地,所述钟漂解算单元利用观测信息中的多普勒频移所构成的接收机速度的解算方程组为如下的线性方程组:
其中,
fd是经过卫星钟漂修正的多普勒频移观测值;fTj是经过卫星钟漂修正的卫星发射信号的射频频率;()为接收机的速度;(vxj,vyj,vzj)为卫星速度;为接收机的钟漂。
可选地,所述钟差解算单元解算得到接收机的钟差是指:
所述钟差解算单元采用最小二乘法通过下式得到接收机的钟差:
ΔX=(HTH)-1HTΔρ
接收机钟差为
所述钟漂解算单元通过最小二乘法得到接收机的钟漂,公式如下:
g=(HTH)-1Hd
其中,
可选地,所述钟差解算单元在所述接收机位置解算方程组中所添加的高程方程为:
Δρ=axΔxu+ayΔyu+azΔzu
其中,
CGS2000_A=6378137.
e2=CGS2000_F*(2-CGS2000_F)
CGS2000_F=(1./298.257222101)
lat为接收机的纬度,在每次迭代过程中由接收机近似位置经过转换获得;
alt为接收机的高度;
增加的H中的一行为[ax ay az 0];
所述钟漂解算单元在所述接收机速度解算方程组中所添加的高程方程为:dj=0
增加的H中的一行为[ax ay az 0]。
本发明采取基于高程进行授时解算,改善了定位系统的DOP值,在已知接收机高度的时候可以提高授时精度,甚至当只接收3颗卫星时,也可以完成授时,拓宽了授时的范围;其中,接收机的高程可以通过很多途径获取。
附图说明
图1为本发明实施例一的一种北斗授时方法的流程示意图。
具体实施方式
下面将结合附图及实施例对本发明的技术方案进行更详细的说明。
需要说明的是,如果不冲突,本发明实施例以及实施例中的各个特征可以相互结合,均在本发明的保护范围之内。另外,虽然在流程图中示出了逻辑顺序,但是在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤。
实施例一,一种北斗授时方法,如图1所示,包括:
产生采样脉冲;
利用采样脉冲采集所接收的卫星信号中的观测信息;所述观测信息包括卫星信号发射时刻、多普勒频移、载波相位;
利用所述卫星信号发射时刻和载波相位组成的伪距,构成接收机位置解算方程组;在所述接收机位置解算方程组中添加一个高程方程,解算得到接收机的钟差;
利用所述观测信息中的多普勒频移构成接收机速度解算方程组,在所述接收机速度解算方程组中添加一个高程方程,解算得到接收机的钟漂;
根据所述钟差和钟漂调节所述采样脉冲的相位和频率,使采样脉冲与标准时间同步。
本实施例的上述步骤是周期性执行的,将不断进行采集、解算和调节;其中解算钟差和钟漂的步骤不分先后,也可以并行执行。
本实施例的一种实施方式中,所述利用卫星信号发射时刻和载波相位组成的伪距,构成接收机位置解算方程组的步骤具体可以包括:
利用卫星信号接收时刻减去卫星信号发射时刻,得到原始伪距;
利用载波相位对原始伪距进行平滑后得到平滑后的伪距;
利用平滑后的伪距,构成解算接收机位置的线性方程组。
本实施方式中,利用载波相位对所述原始伪距进行平滑,得到平滑后的伪距的步骤具体可以包括:
通过下式计算平滑后的伪距:
ρs,i=α×ρCA,i+(1-α)×(ρs,i-1+Δρφ,i)
其中,ρs,i是平滑后的伪距;α为平滑系数;ρCA,i是原始伪距;ρs,i-1是前一次平滑后的伪距;Δρφ,i是伪距变化量;
所述伪距变化量为:
λ是所接收的卫星信号的波长;φi是本次采样得到的载波相位值;φi-1是前一次采样得到的载波相位值;fIF是所接收的卫星信号的载波中频频率;为接收机钟漂;X是采样时刻锁存的载波相位计数值;fL1是所接收的卫星信号的载波射频频率。
本实施方式中,所述α可以但不限于为0.05。
本实施方式中,利用平滑后的伪距所构成的解算接收机位置的解算方程组具体可以为如下的线性方程组:
ρi为第i颗卫星的平滑后的伪距,i的取值从1到n,n为接收机所能连接的卫星的颗数;(xu,yu,zu)为接收机的位置坐标;(xi,yi,zi)为第i颗卫星的位置坐标;tu为接收机钟差;ri为接收机到第i颗卫星的距离;c表示光速;带^标示的表示该变量为系统估计值。
本实施方式中,解算得到接收机的钟差的步骤具体可以包括:
采用最小二乘法通过下式得到接收机的钟差:
ΔX=(HTH)-1HTΔρ
接收机钟差为
本实施方式中,在所述接收机位置解算方程组中所添加的高程方程具体可以为:
Δρ=axΔxu+ayΔyu+azΔzu
其中,
CGS2000_A=6378137.
e2=CGS2000_F*(2-CGS2000_F)
CGS2000_F=(1./298.257222101)
lat为接收机的纬度,在每次迭代过程中由接收机近似位置经过转换获得。
alt为接收机的高度。
增加的H中的一行为[ax ay az 0]。
本实施例的一种实施方式中,利用观测信息中的多普勒频移所构成的接收机速度的线性方程组具体可以如下:
其中,
fd是经过卫星钟漂修正的多普勒频移观测值;fTj是经过卫星钟漂修正的卫星发射信号的射频频率;()为接收机的速度;(vxj,vyj,vzj)为卫星速度;为接收机的钟漂。
该线性方程组可以通过最小二乘法得到解,公式如下:
g=(HTH)-1Hd
其中,
系数矩阵H同前文。
本实施方式中,在所述接收机速度解算方程组中所添加的高程方程为:
dj=0
增加的H中的一行为[ax ay az 0],(ax,ay,az)同前文。
本实施例中,采样脉冲可以通过DDS(Direct Digital Synthesizer,直接数字式频率合成器)技术中的相位累加器的方式产生。相位累加器的位数取64位。相位累加器位数的选取取决于FPGA系统主时钟频率和接收机晶振的频率精度。如果接收机晶振的频率精确度能达到0.1~1ppm。可见累加器位数取32位的话,系统主时钟是40MHz时,相位累加器的频率分辨率为0.00931322574615478515625Hz,不能满足精度要求。
FPGA中的频率控制字寄存器将接收到的频率控制字送入相位累加器,每次FPGA系统主时钟的上升沿到来时,相位累加器就会在相位初值的基础上对频率控制字进行累加。每当相位累加器的最高位溢出时,就会产生一个采样脉冲。
如图1所示,利用平滑后的伪距,构成接收机位置解算方程组,得到接收机钟差。而采取基于高程的授时算法时,要添加一个基于高程的方程。这样可以在已知接收机高程的情况下,完成3颗卫星授时,还可以改善一定条件下的卫星定位系统的DOP值,得到更加准确的接收机钟差。公式如下:
本实施方式中,在得到接收机的钟差和钟漂后,利用钟差去调节采样脉冲发生器的相位,而利用钟漂去调节采样脉冲发生器的频率。将相位值转化为频率值,频率控制字通过下面公式得到:
f为FPGA的系统时钟频率;
t为采样脉冲的周期;
tu为接收机的钟差;
在本实施例中,用钟漂去调节采样脉冲的频率,主要通过调节计算频率控制字的FPGA系统时钟频率来达到调节采样脉冲频率的目的,即将上面公式中的f改为
实施例二、一种卫星导航接收机,包括:
脉冲生成模块,用于产生采样脉冲;
采集单元,用于利用所述采样脉冲采集所接收的卫星信号中的观测信息;所述观测信息包括卫星信号发射时刻、多普勒频移、载波相位;
钟差解算单元,用于利用所述卫星信号发射时刻和载波相位组成的伪距,构成接收机位置解算方程组;在所述接收机位置解算方程组中添加一个高程方程,解算得到接收机的钟差;
钟漂解算单元,用于利用所述观测信息中的多普勒频移构成接收机速度解算方程组,在所述接收机速度解算方程组中添加一个高程方程,解算得到接收机的钟漂;
调整单元,用于根据所述钟差和钟漂调节所述采样脉冲的相位和频率,使采样脉冲与标准时间同步。
本实施例的一种实施方式中,所述钟差解算单元利用卫星信号发射时刻和载波相位组成的伪距,构成接收机位置解算方程组具体可以是指:
所述钟差解算单元利用卫星信号接收时刻减去卫星信号发射时刻,得到原始伪距;通过下式计算平滑后的伪距:
ρs,i=α×ρCA,i+(1-α)×(ρs,i-1+Δρφ,i)
其中,ρs,i是平滑后的伪距;α为平滑系数;ρCA,i是原始伪距;ρs,i-1是前一次平滑后的伪距;Δρφ,i是伪距变化量;
所述伪距变化量为:
λ是所接收的卫星信号的波长;φi是本次采样得到的载波相位值;φi-1是前一次采样得到的载波相位值;fIF是所接收的卫星信号的载波中频频率;为接收机钟漂;X是采样时刻锁存的载波相位计数值;fL1是所接收的卫星信号的载波射频频率;
利用平滑后的伪距,构成解算接收机位置的线性方程组如下:
ρi为第i颗卫星的平滑后的伪距,i的取值从1到n,n为接收机所能连接的卫星的颗数;(xu,yu,zu)为接收机的位置坐标;(xi,yi,zi)为第i颗卫星的位置坐标;tu为接收机钟差;ri为接收机到第i颗卫星的距离;c表示光速;带^标示的表示该变量为系统估计值。
本实施方式中,,所述钟漂解算单元利用观测信息中的多普勒频移所构成的接收机速度的解算方程组具体可以为如下的线性方程组:
其中,
fd是经过卫星钟漂修正的多普勒频移观测值;fTj是经过卫星钟漂修正的卫星发射信号的射频频率;()为接收机的速度;(vxj,vyj,vzj)为卫星速度;为接收机的钟漂。
本实施方式的一种备选方案中,所述钟差解算单元解算得到接收机的钟差具体可以是指:
所述钟差解算单元采用最小二乘法通过下式得到接收机的钟差:
ΔX=(HTH)-1HTΔρ
接收机钟差为
所述钟漂解算单元通过最小二乘法得到接收机的钟漂,公式可以如下:
g=(HTH)-1Hd
其中,
本备选方案中,所述钟差解算单元在所述接收机位置解算方程组中所添加的高程方程具体可以为:
Δρ=axΔxu+ayΔyu+azΔzu
其中,
CGS2000_A=6378137.
e2=CGS2000_F*(2-CGS2000_F)
CGS2000_F=(1./298.257222101)
lat为接收机的纬度,在每次迭代过程中由接收机近似位置经过转换获得;
alt为接收机的高度;
增加的H中的一行为[ax ay az 0];
所述钟漂解算单元在所述接收机速度解算方程组中所添加的高程方程为:dj=0;
增加的H中的一行为[ax ay az 0]。
本领域普通技术人员可以理解上述方法中的全部或部分步骤可通过程序来指令相关硬件完成,所述程序可以存储于计算机可读存储介质中,如只读存储器、磁盘或光盘等。可选地,上述实施例的全部或部分步骤也可以使用一个或多个集成电路来实现。相应地,上述实施例中的各模块/单元可以采用硬件的形式实现,也可以采用软件功能模块的形式实现。本发明不限制于任何特定形式的硬件和软件的结合。
当然,本发明还可有其他多种实施例,在不背离本发明精神及其实质的情况下,熟悉本领域的技术人员当可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于本发明的权利要求的保护范围。
Claims (10)
1.一种北斗授时方法,包括:
产生采样脉冲;
利用采样脉冲采集所接收的卫星信号中的观测信息;所述观测信息包括卫星信号发射时刻、多普勒频移、载波相位;
利用所述卫星信号发射时刻和载波相位组成的伪距,构成接收机位置解算方程组;在所述接收机位置解算方程组中添加一个高程方程,解算得到接收机的钟差;
利用所述观测信息中的多普勒频移构成接收机速度解算方程组,在所述接收机速度解算方程组中添加一个高程方程,解算得到接收机的钟漂;
根据所述钟差和钟漂调节所述采样脉冲的相位和频率,使采样脉冲与标准时间同步。
2.如权利要求1所述的方法,其特征在于,所述利用卫星信号发射时刻和载波相位组成的伪距,构成接收机位置解算方程组的步骤包括:
利用卫星信号接收时刻减去卫星信号发射时刻,得到原始伪距;
通过下式计算平滑后的伪距:
ρs,i=α×ρCA,i+(1-α)×(ρs,i-1+Δρφ,i)
其中,ρs,i是平滑后的伪距;α为平滑系数;ρCA,i是原始伪距;ρs,i-1是前一次平滑后的伪距;Δρφ,i是伪距变化量;
所述伪距变化量为:
λ是所接收的卫星信号的波长;φi是本次采样得到的载波相位值;φi-1是前一次采样得到的载波相位值;fIF是所接收的卫星信号的载波中频频率;为接收机钟漂;X是采样时刻锁存的载波相位计数值;fL1是所接收的卫星信号的载波射频频率;
利用平滑后的伪距,构成解算接收机位置的线性方程组如下:
ρi为第i颗卫星的平滑后的伪距,i的取值从1到n,n为接收机所能连接的卫星的颗数;(xu,yu,zu)为接收机的位置坐标;(xi,yi,zi)为第i颗卫星的位置坐标;tu为接收机钟差;ri为接收机到第i颗卫星的距离;c表示光速;带^标示的表示该变量为系统估计值。
3.如权利要求2所述的方法,其特征在于,利用观测信息中的多普勒频移所构成的接收机速度的解算方程组为如下的线性方程组:
其中,
fd是经过卫星钟漂修正的多普勒频移观测值;fTj是经过卫星钟漂修正的卫星发射信号的射频频率;为接收机的速度;(vxj,vyj,vzj)为卫星速度;为接收机的钟漂。
4.如权利要求3所述的方法,其特征在于,解算得到接收机的钟差的步骤包括:
采用最小二乘法通过下式得到接收机的钟差:
ΔX=(HTH)-1HTΔρ
接收机钟差为
通过最小二乘法得到接收机的钟漂,公式如下:
g=(HTH)-1Hd
其中,
5.如权利要求4所述的方法,其特征在于,在所述接收机位置解算方程组中所添加的高程方程为:
Δρ=axΔxu+ayΔyu+azΔzu
其中,
CGS2000_A=6378137.
e2=CGS2000_F*(2-CGS2000_F)
CGS2000_F=(1./298.257222101)
lat为接收机的纬度,在每次迭代过程中由接收机近似位置经过转换获得;
alt为接收机的高度;
增加的H中的一行为[ax ay az 0];
在所述接收机速度解算方程组中所添加的高程方程为:
dj=0
增加的H中的一行为[ax ay az 0]。
6.一种卫星导航接收机,其特征在于,包括:
脉冲生成模块,用于产生采样脉冲;
采集单元,用于利用所述采样脉冲采集所接收的卫星信号中的观测信息;所述观测信息包括卫星信号发射时刻、多普勒频移、载波相位;
钟差解算单元,用于利用所述卫星信号发射时刻和载波相位组成的伪距,构成接收机位置解算方程组;在所述接收机位置解算方程组中添加一个高程方程,解算得到接收机的钟差;
钟漂解算单元,用于利用所述观测信息中的多普勒频移构成接收机速度解算方程组,在所述接收机速度解算方程组中添加一个高程方程,解算得到接收机的钟漂;
调整单元,用于根据所述钟差和钟漂调节所述采样脉冲的相位和频率,使采样脉冲与标准时间同步。
7.如权利要求6所述的卫星导航接收机,其特征在于,所述钟差解算单元利用卫星信号发射时刻和载波相位组成的伪距,构成接收机位置解算方程组是指:
所述钟差解算单元利用卫星信号接收时刻减去卫星信号发射时刻,得到原始伪距;通过下式计算平滑后的伪距:
ρs,i=α×ρCA,i+(1-α)×(ρs,i-1+Δρφ,i)
其中,ρs,i是平滑后的伪距;α为平滑系数;ρCA,i是原始伪距;ρs,i-1是前一次平滑后的伪距;Δρφ,i是伪距变化量;
所述伪距变化量为:
λ是所接收的卫星信号的波长;φi是本次采样得到的载波相位值;φi-1是前一次采样得到的载波相位值;fIF是所接收的卫星信号的载波中频频率;为接收机钟漂;X是采样时刻锁存的载波相位计数值;fL1是所接收的卫星信号的载波射频频率;
利用平滑后的伪距,构成解算接收机位置的线性方程组如下:
ρi为第i颗卫星的平滑后的伪距,i的取值从1到n,n为接收机所能连接的卫星的颗数;(xu,yu,zu)为接收机的位置坐标;(xi,yi,zi)为第i颗卫星的位置坐标;tu为接收机钟差;ri为接收机到第i颗卫星的距离;c表示光速;带^标示的表示该变量为系统估计值。
8.如权利要求7所述的卫星导航接收机,其特征在于,所述钟漂解算单元利用观测信息中的多普勒频移所构成的接收机速度的解算方程组为如下的线性方程组:
其中,
fd是经过卫星钟漂修正的多普勒频移观测值;fTj是经过卫星钟漂修正的卫星发射信号的射频频率;为接收机的速度;(vxj,vyj,vzj)为卫星速度;为接收机的钟漂。
9.如权利要求8所述的卫星导航接收机,其特征在于,所述钟差解算单元解算得到接收机的钟差是指:
所述钟差解算单元采用最小二乘法通过下式得到接收机的钟差:
ΔX=(HTH)-1HTΔρ
接收机钟差为
所述钟漂解算单元通过最小二乘法得到接收机的钟漂,公式如下:
g=(HTH)-1Hd
其中,
10.如权利要求6所述的卫星导航接收机,其特征在于,所述钟差解算单元在所述接收机位置解算方程组中所添加的高程方程为:
Δρ=axΔxu+ayΔyu+azΔzu
其中,
CGS2000_A=6378137.
e2=CGS2000_F*(2-CGS2000_F)
CGS2000_F=(1./298.257222101)
lat为接收机的纬度,在每次迭代过程中由接收机近似位置经过转换获得;
alt为接收机的高度;
增加的H中的一行为[ax ay az 0];
所述钟漂解算单元在所述接收机速度解算方程组中所添加的高程方程为:dj=0
增加的H中的一行为[ax ay az 0]。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410307499.0A CN104049528B (zh) | 2014-06-30 | 2014-06-30 | 一种北斗授时方法及卫星导航接收机 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410307499.0A CN104049528B (zh) | 2014-06-30 | 2014-06-30 | 一种北斗授时方法及卫星导航接收机 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104049528A CN104049528A (zh) | 2014-09-17 |
CN104049528B true CN104049528B (zh) | 2017-05-17 |
Family
ID=51502559
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410307499.0A Active CN104049528B (zh) | 2014-06-30 | 2014-06-30 | 一种北斗授时方法及卫星导航接收机 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104049528B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104678408B (zh) * | 2015-02-05 | 2017-02-22 | 清华大学 | 星载导航接收机授时方法和授时型星载导航接收机以及星载导航应用系统 |
CN105871495B (zh) * | 2015-12-29 | 2018-04-17 | 中国科学院国家天文台 | 一种时间同步方法、通信地面站和用户终端 |
CN105929192B (zh) * | 2016-04-13 | 2018-10-30 | 中国电子科技集团公司第五十四研究所 | 一种基于gnss自主测速的测风装置及测风方法 |
CN109388591B (zh) * | 2017-08-07 | 2021-07-16 | 北京北斗星通导航技术股份有限公司 | 一种matlab应用系统与stk应用系统接口适配器和适配方法 |
CN108196288A (zh) * | 2018-01-11 | 2018-06-22 | 南京理工大学 | 基于微惯性、芯片原子钟辅助北斗接收机的重定位方法 |
CN108646543B (zh) * | 2018-06-15 | 2020-11-10 | 西安空间无线电技术研究所 | 一种具有高稳定度性能的gnss驯服时钟方法 |
CN109001972B (zh) * | 2018-08-13 | 2020-06-12 | 中国科学院国家授时中心 | 一种北斗广域授时系统与方法 |
CN110837221B (zh) * | 2019-11-27 | 2021-01-19 | 深圳市力合微电子股份有限公司 | 一种有效提高授时可靠性与连续性的方法 |
CN111398994B (zh) * | 2020-04-26 | 2020-10-30 | 清华大学 | 中轨通信卫星定位授时方法和装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103235321A (zh) * | 2013-03-29 | 2013-08-07 | 东南大学 | Gps单站伪距定位精密授时方法 |
CN103698781A (zh) * | 2013-12-27 | 2014-04-02 | 北京北斗星通导航技术股份有限公司 | 一种接收机观测信息的提取方法和装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7564406B2 (en) * | 2006-11-10 | 2009-07-21 | Sirf Technology, Inc. | Method and apparatus in standalone positioning without broadcast ephemeris |
US8497801B2 (en) * | 2007-02-05 | 2013-07-30 | Qualcomm Incorporated | Prediction refresh method for ephemeris extensions |
-
2014
- 2014-06-30 CN CN201410307499.0A patent/CN104049528B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103235321A (zh) * | 2013-03-29 | 2013-08-07 | 东南大学 | Gps单站伪距定位精密授时方法 |
CN103698781A (zh) * | 2013-12-27 | 2014-04-02 | 北京北斗星通导航技术股份有限公司 | 一种接收机观测信息的提取方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN104049528A (zh) | 2014-09-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104049528B (zh) | 一种北斗授时方法及卫星导航接收机 | |
Zhao et al. | An open source GPS/GNSS vector tracking loop-implementation, filter tuning, and results | |
CN103235327B (zh) | 一种gnss/mins超深组合导航方法、系统及装置 | |
CN103105615B (zh) | 卫星导航信号的虚假检测方法及卫星导航定位接收机 | |
US9172524B2 (en) | Method of generating unambiguous correlation function for TMBOC (6,1,4/33)signal based on partial correlation functions, apparatus for tracking TMBOC signal, and satellite navigation signal receiver system using the same | |
CN118409339A (zh) | 近地轨道卫星通信系统中的位置估计 | |
Musumeci et al. | Design of a High Sensitivity GNSS receiver for Lunar missions | |
CN103823206B (zh) | 一种基于导航卫星的星地双基地sar时频同步方法 | |
CN111149018A (zh) | 用于校准系统参数的方法和系统 | |
CN104849731A (zh) | 一种天线阵元通道的校准方法、装置及接收机 | |
EP3141931B1 (en) | Tracking of signals with at least one subcarrier | |
Qin et al. | Performance assessment of a low‐cost inertial measurement unit based ultra‐tight global navigation satellite system/inertial navigation system integration for high dynamic applications | |
JP2006349672A (ja) | Gpsr多周計測型電離層遅延補正装置、gpsr多周計測型電離層遅延補正方法及びプログラム | |
CN103809192A (zh) | 一种gnss接收机的动态校正算法 | |
Jiao et al. | Carrier phase ranging with DTMB signals for urban pedestrian localization and GNSS aiding | |
CN105699992A (zh) | 高动态gnss载波信号跟踪方法及系统 | |
CN103116171A (zh) | 卫星导航信号的虚假检测方法及卫星导航定位接收机 | |
CN104309817B (zh) | 基于多台并址接收机的北斗导航卫星区域定轨方法 | |
KR101838406B1 (ko) | Boc 신호를 추적하는 방법 및 boc 신호 추적 장치 | |
CN106154300B (zh) | 一种惯性/卫星超紧组合实现方法 | |
CN105425257A (zh) | 一种高动态gnss载波信号的跟踪方法及系统 | |
JP6047944B2 (ja) | 受信装置及び相関積算処理方法 | |
CN103116172A (zh) | 卫星导航信号的虚假检测方法及卫星导航定位接收机 | |
CN106125117A (zh) | 一种惯性/卫星超紧组合本地信号控制量生成方法 | |
CN108885266A (zh) | 用于利用蜂窝cdma信号导航的sdr |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |