CN111044972A - 基于gnss精密时间同步的飞行器时差定位方法及系统 - Google Patents
基于gnss精密时间同步的飞行器时差定位方法及系统 Download PDFInfo
- Publication number
- CN111044972A CN111044972A CN201911276822.1A CN201911276822A CN111044972A CN 111044972 A CN111044972 A CN 111044972A CN 201911276822 A CN201911276822 A CN 201911276822A CN 111044972 A CN111044972 A CN 111044972A
- Authority
- CN
- China
- Prior art keywords
- base station
- gnss
- difference
- time
- observation
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/02—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
- G01S5/06—Position of source determined by co-ordinating a plurality of position lines defined by path-difference measurements
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/24—Acquisition or tracking or demodulation of signals transmitted by the system
- G01S19/25—Acquisition or tracking or demodulation of signals transmitted by the system involving aiding data received from a cooperating element, e.g. assisted GPS
- G01S19/256—Acquisition or tracking or demodulation of signals transmitted by the system involving aiding data received from a cooperating element, e.g. assisted GPS relating to timing, e.g. time of week, code phase, timing offset
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W4/00—Services specially adapted for wireless communication networks; Facilities therefor
- H04W4/02—Services making use of location information
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明实施例提供一种基于GNSS精密时间同步的飞行器时差定位方法及系统,该方法包括:基于GNSS精密时间,对多个基站进行时间同步处理,以使得每个基站之间时间精密同步;根据时间精密同步后的基站中GNSS接收机中的时间,获取目标飞行器信号源对应各个基站的时差观测值,所述时差观测值为所述目标飞行器信号源发送的同一帧无线电信号到达各个基站的时差;根据所述时差观测值以及基站的已知坐标,获取所述目标飞行器信号源的坐标,以用于对所述目标飞行器进行定位。本发明实施例提高了飞行器的定位精度。
Description
技术领域
本发明涉及无线电测量和定位技术领域,尤其涉及一种基于GNSS精密时间同步的飞行器时差定位方法及系统。
背景技术
到达时差(Time Difference of Arrival,简称TDOA)定位技术在无线电监测和航空领域具有十分重要的作用,具有精度高、覆盖范围广和无源的特点。
现有基于TDOA的时差定位方法,通过从信号源发出的无线电信号到达基站(4处以上)的时间差,再根据基站的位置信息,计算出信号源的位置坐标,而定位精度主要由时间同步精度和无线电信号到达不同基站的时差测量精度决定。然而,在现有的时差定位方法中,对于接收基站的时间同步的问题,暂未提供有效的解决方案,导致基站间的时间同步效果较差,降低了飞行器的定位精度。
因此,现在亟需一种基于GNSS精密时间同步的飞行器时差定位及系统来解决上述问题。
发明内容
针对现有技术存在的问题,本发明实施例提供一种基于GNSS精密时间同步的飞行器时差定位方法及系统。
第一方面,本发明实施例提供了一种基于GNSS精密时间同步的飞行器时差定位方法,包括:
基于GNSS精密时间,对多个基站进行时间同步处理,以使得所有基站之间的时间精密同步;
根据时间精密同步后的基站中GNSS接收机中的时间,获取目标飞行器信号源对应各个基站的时差观测值,所述时差观测值为所述目标飞行器信号源发送的同一帧无线电信号到达各个基站的时差;
根据所述时差观测值,获取所述目标飞行器信号源的坐标,以用于对所述目标飞行器进行定位。
进一步地,所述基于GNSS精密时间,对多个基站进行时间同步处理,以使得每个基站之间的时间精密同步,包括:
获取每个基站的GNSS观测数据,并通过探测粗差,将GNSS观测数据中存在粗差的GNSS卫星剔除,以使得每个基站的GNSS观测数据来源于不存在粗差的GNSS卫星,所述GNSS观测数据包括双频伪距和载波相位观测值;
根据剔除粗差后每个基站中的GNSS观测数据,获取每两个基站之间的GNSS共视卫星,并根据所述GNSS共视卫星的观测数据,获取所述GNSS共视卫星的站间单差无电离层组合伪距观测值和站间单差无电离层组合载波相位观测值,以用于构建单差观测方程;
将单差GNSS接收机钟差、天顶对流层湿延迟和站间单差无电离层组合模糊度项作为参考估计,并根据所述站间单差无电离层组合伪距观测值和所述站间单差无电离层组合载波相位观测值,构建误差方程;
根据法方程和所述误差方程,获取单差GNSS接收机钟差,并将其中任一基站的GNSS接收机时间作为参考时间,以使得各个基站的时间精密同步。
进一步地,在所述根据剔除粗差后每个基站中的GNSS观测数据,获取每两个基站之间的GNSS共视卫星,并根据所述GNSS共视卫星的观测数据,获取所述GNSS共视卫星的站间单差无电离层组合伪距观测值和站间单差无电离层组合载波相位观测值,以用于构建单差观测方程之后,所述方法还包括:
对所述站间单差无电离层组合伪距观测值和所述站间单差无电离层组合载波相位观测值中的误差进行修正,以根据误差修正后的站间单差无电离层组合伪距观测值和站间单差无电离层组合载波相位观测值,构建误差方程。
进一步地,所述对所述站间单差无电离层组合伪距观测值和所述站间单差无电离层组合载波相位观测值中的误差进行修正,包括:
对所述站间单差无电离层组合伪距观测值和所述站间单差无电离层组合载波相位观测值中的地球自转效应误差、卫星天线相位中心偏差、卫星天线相位中心变化、接收机天线相位中心偏差、接收机天线相位中心变化、相对论效应、引力延迟、相位缠绕和对流层干延迟进行修正,得到误差修正后的站间单差无电离层组合伪距观测值和站间单差无电离层组合载波相位观测值。
进一步地,所述单差观测方程为:
其中,表示第i个GNSS卫星下基站a和基站b的站间单差无电离层组合伪距观测值;表示第i个GNSS卫星下基站a和基站b的站间单差卫地距;tab表示基站a和基站b之间的单差接收机钟差;表示第i个GNSS卫星下基站a和基站b的站间单差对流层延迟;表示第i个GNSS卫星下基站a和基站b的站间单差无电离层组合载波相位观测值,单位为米;表示第i个GNSS卫星下基站a和基站b的单差无电离层组合模糊度项,单位为米;a和b表示GNSS接收机标识,i表示GNSS卫星标识;
其中,表示第i个GNSS卫星下基站a的L1载波的伪距观测值;表示第i个GNSS卫星下基站b的L1载波的伪距观测值;表示第i个GNSS卫星下基站a的L2载波的伪距观测值;表示第i个GNSS卫星下基站b的L2载波的伪距观测值;表示第i个GNSS卫星下基站a的L1载波的载波相位观测值;表示第i个GNSS卫星下基站b的L1载波的载波相位观测值;表示第i个GNSS卫星下基站a的L2载波的载波相位观测值;表示第i个GNSS卫星下基站b的L2载波的载波相位观测值;ρ表示卫地距;λ1和λ2分别表示L1载波和L2载波的波长;表示第i个GNSS卫星下基站a的L1载波的模糊度;表示第i个GNSS卫星下基站b的L1载波的模糊度;表示第i个GNSS卫星下基站a的L2载波的模糊度;表示第i个GNSS卫星下基站b的L2载波的模糊度;表示第i个GNSS卫星下基站a的对流层干延迟投影函数;表示第i个GNSS卫星下基站b的对流层干延迟投影函数;表示第i个GNSS卫星下基站a的对流层湿延迟投影函数;表示第i个GNSS卫星下基站b的对流层湿延迟投影函数;TD,a表示基站a的天顶对流层干延迟;TD,b表示基站b的天顶对流层干延迟;TW,a表示基站a的天顶对流层湿延迟;TW,b表示基站b的天顶对流层湿延迟;f1表示L1载波频率,f2表示L2载波频率,(xa,ya,za)表示基站a的GNSS接收机坐标;(xb,yb,zb)表示基站b的GNSS接收机坐标;表示基站a对应的第i个GNSS卫星的坐标;表示基站b对应的第i个GNSS卫星的坐标。
进一步地,所述误差方程为:
V=BX-l,W;
其中,B表示设计矩阵,X表示待估参数向量,l表示观测值向量,V表示观测值残差向量,W表示观测值向量权阵,所述观测值向量权阵根据观测值先验精度获取得到的;
其中,和分别表示误差修正后的站间单差无电离层组合伪距观测值和站间单差无电离层组合载波相位观测值;n表示GNSS卫星总数;表示第n个GNSS卫星对应的基站a的非差无电离层组合伪距观测值先验方差,表示第n个GNSS卫星对应的基站b的非差无电离层组合伪距观测值先验方差,表示第n个GNSS卫星对应的基站a的非差无电离层组合载波相位观测值先验方差,表示第n个GNSS卫星对应的基站b的非差无电离层组合载波相位观测值先验方差。
进一步地,所述根据所述时差观测值,获取所述目标飞行器信号源的信号源坐标,以用于对所述目标飞行器信号源进行定位,包括:
根据接收目标飞行器信号源的基站对应的坐标信息,获取所述目标飞行器信号源的概略坐标;
将任一接收目标飞行器信号源的基站k作为参考基站,根据时差观测值和基站d坐标信息,d,k∈K,d≠k,构建观测方程:
其中,(xd,yd,zd)表示第d个基站的坐标信息,(x,y,z)表示目标飞行器信号源的坐标信息,c表示光速;Δtd表示当基站k为参考基站时,目标飞行器信号源发送的同一帧无线电信号到达第d个基站与参考基站k的时差;
根据接收目标飞行器信号源的基站总数K,构建K-1个观测方程,并根据泰勒级展开法和所述目标飞行器信号源的概略坐标,对K-1个观测方程进行求解,得到所述目标飞行器信号源的信号源坐标,以用于对所述目标飞行器信号源进行定位。
第二方面,本发明实施例提供了一种基于GNSS精密时间同步的飞行器时差定位系统,包括:
精密时间同步模块,基于GNSS精密时间,对多个基站进行时间同步处理,以使得所有基站之间的时间精密同步;
时差获取模块,用于根据时间精密同步后的基站中GNSS接收机中的时间,获取目标飞行器信号源对应各个基站的时差观测值,所述时差观测值为所述目标飞行器信号源发送的同一帧无线电信号到达各个基站的时差;
信号源定位模块,用于根据所述时差观测值,获取所述目标飞行器信号源的坐标,以用于对所述目标飞行器进行定位。
第三方面,本发明实施例提供一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现如第一方面所提供的方法的步骤。
第四方面,本发明实施例提供一种非暂态计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时实现如第一方面所提供的方法的步骤。
本发明实施例提供的基于GNSS精密时间同步的飞行器时差定位方法及系统,该方法不需要原子钟,即可实现基站间精密时间同步,解决了时差定位系统中基站间精密时间同步的问题,进而实现高精度时差定位,提高了飞行器的定位精度,降低了GNSS星历误差的影响。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的基于GNSS精密时间同步的飞行器时差定位方法的流程示意图;
图2为本发明实施例提供的基于GNSS精密时间同步的时差定位示意图;
图3为本发明实施例提供的IGS跟踪站中PTBB-POTS接收机钟时间比对误差的示意图;
图4为本发明实施例提供的飞行器信号源定位结果的仿真示意图;
图5为本发明实施例提供的飞行器信号源的三维定位误差示意图;
图6为本发明实施例提供的基于GNSS精密时间同步的飞行器时差定位系统的结构示意图;
图7为本发明实施例提供的电子设备结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
图1为本发明实施例提供的基于GNSS精密时间同步的飞行器时差定位方法的流程示意图,如图1所示,本发明实施例提供了一种基于GNSS精密时间同步的飞行器时差定位方法,包括:
步骤101,基于GNSS精密时间,对多个基站进行时间同步处理,以使得所有基站之间的时间精密同步;
步骤102,根据时间精密同步后的基站中GNSS接收机中的时间,获取目标飞行器信号源对应各个基站的时差观测值,所述时差观测值为所述目标飞行器信号源发送的同一帧无线电信号到达各个基站的时差;
步骤103,根据所述时差观测值,获取所述目标飞行器信号源的坐标,以用于对所述目标飞行器进行定位。
在本发明实施例中,图2为本发明实施例提供的基于GNSS精密时间同步的时差定位示意图,可参考图2所示,在每个定位基站中均安装有高精度GNSS接收机,基于4个以上基站中GNSS(Global Navigation Satellite System)接收机采集的数据,通过GNSS精密时间比对,实现定位基站间精密时间同步。然后,以完成时间同步后的定位基站GNSS接收机钟的时间为基准,测量目标飞行器信号源发出的同一帧无线电信号到达不同基站的时差观测值;最后,数据中心根据时差观测值,计算目标飞行器信号源的坐标,从而实现目标飞行器信号源的定位。
本发明实施例提供的基于GNSS精密时间同步的飞行器时差定位方法,该方法不需要原子钟,即可实现基站间精密时间同步,解决了时差定位系统中基站间精密时间同步的问题,进而实现高精度时差定位,提高了飞行器的定位精度,降低了GNSS星历误差的影响。
在上述实施例的基础上,所述基于GNSS精密时间,对多个基站进行时间同步处理,以使得每个基站之间的时间精密同步,包括:
获取每个基站的GNSS观测数据,并通过探测粗差,将GNSS观测数据中存在粗差的GNSS卫星剔除,以使得每个基站的GNSS观测数据来源于不存在粗差的GNSS卫星,所述GNSS观测数据包括双频伪距和载波相位观测值;
根据剔除粗差后每个基站中的GNSS观测数据,以两处位置的基站为一组,获取每两个基站之间的GNSS共视卫星,并根据所述GNSS共视卫星的观测数据,获取所述GNSS共视卫星的站间单差无电离层组合伪距观测值和站间单差无电离层组合载波相位观测值,以用于构建单差观测方程;具体地,在本发明实施例中,所述单差观测方程为:
其中,表示第i个GNSS卫星下基站a和基站b的站间单差无电离层组合伪距观测值;表示第i个GNSS卫星下基站a和基站b的站间单差卫地距;tab表示基站a和基站b之间的单差接收机钟差;表示第i个GNSS卫星下基站a和基站b的站间单差对流层延迟;表示第i个GNSS卫星下基站a和基站b的站间单差无电离层组合载波相位观测值,单位为米;表示第i个GNSS卫星下基站a和基站b的单差无电离层组合模糊度项,单位为米;a和b表示GNSS接收机标识,i表示GNSS卫星标识;
其中,表示第i个GNSS卫星下基站a的L1载波的伪距观测值;表示第i个GNSS卫星下基站b的L1载波的伪距观测值;表示第i个GNSS卫星下基站a的L2载波的伪距观测值;表示第i个GNSS卫星下基站b的L2载波的伪距观测值;表示第i个GNSS卫星下基站a的L1载波的载波相位观测值;表示第i个GNSS卫星下基站b的L1载波的载波相位观测值;表示第i个GNSS卫星下基站a的L2载波的载波相位观测值;表示第i个GNSS卫星下基站b的L2载波的载波相位观测值;ρ表示卫地距;λ1和λ2分别表示L1载波和L2载波的波长;表示第i个GNSS卫星下基站a的L1载波的模糊度;表示第i个GNSS卫星下基站b的L1载波的模糊度;表示第i个GNSS卫星下基站a的L2载波的模糊度;表示第i个GNSS卫星下基站b的L2载波的模糊度;表示第i个GNSS卫星下基站a的对流层干延迟投影函数;表示第i个GNSS卫星下基站b的对流层干延迟投影函数;表示第i个GNSS卫星下基站a的对流层湿延迟投影函数;表示第i个GNSS卫星下基站b的对流层湿延迟投影函数;TD,a表示基站a的天顶对流层干延迟;TD,b表示基站b的天顶对流层干延迟;TW,a表示基站a的天顶对流层湿延迟;TW,b表示基站b的天顶对流层湿延迟;f1表示L1载波频率,f2表示L2载波频率,(xa,ya,za)表示基站a的GNSS接收机坐标;(xb,yb,zb)表示基站b的GNSS接收机坐标;表示基站a对应的第i个GNSS卫星的坐标;表示基站b对应的第i个GNSS卫星的坐标。其中,GNSS接收机坐标是根据长时间连续的GNSS观测数据提前计算得到的;GNSS卫星的坐标是通过IGS超快速预报星历计算得到的。
将单差GNSS接收机钟差、天顶对流层湿延迟和站间单差无电离层组合模糊度项作为参考估计,并根据所述站间单差无电离层组合伪距观测值和所述站间单差无电离层组合载波相位观测值,构建误差方程;
优选地,在上述实施例的基础上,在所述根据剔除粗差后每个基站中的GNSS观测数据,获取每两个基站之间的GNSS共视卫星,并根据所述GNSS共视卫星的观测数据,获取所述GNSS共视卫星的站间单差无电离层组合伪距观测值和站间单差无电离层组合载波相位观测值,以用于构建单差观测方程之后,所述方法还包括:
对所述站间单差无电离层组合伪距观测值和所述站间单差无电离层组合载波相位观测值中的误差进行修正,以根据误差修正后的站间单差无电离层组合伪距观测值和站间单差无电离层组合载波相位观测值,构建误差方程。
进一步地,在上述实施例的基础上,所述对所述站间单差无电离层组合伪距观测值和所述站间单差无电离层组合载波相位观测值中的误差进行修正,包括:
对所述站间单差无电离层组合伪距观测值和所述站间单差无电离层组合载波相位观测值中的地球自转效应误差、卫星天线相位中心偏差、卫星天线相位中心变化、接收机天线相位中心偏差、接收机天线相位中心变化、相对论效应、引力延迟、相位缠绕和对流层干延迟进行修正,得到误差修正后的站间单差无电离层组合伪距观测值和站间单差无电离层组合载波相位观测值。
在本发明实施例中,所述误差方程为:
V=BX-l,W;
其中,B表示设计矩阵,X表示待估参数向量,l表示观测值向量,V表示观测值残差向量,W表示观测值向量权阵,所述观测值向量权阵根据观测值先验精度获取得到的;
其中,和分别表示误差修正后的站间单差无电离层组合伪距观测值和站间单差无电离层组合载波相位观测值;n表示GNSS卫星总数;表示第n个GNSS卫星对应的基站a的非差无电离层组合伪距观测值先验方差,表示第n个GNSS卫星对应的基站b的非差无电离层组合伪距观测值先验方差,表示第n个GNSS卫星对应的基站a的非差无电离层组合载波相位观测值先验方差,表示第n个GNSS卫星对应的基站b的非差无电离层组合载波相位观测值先验方差。在本发明实施例中,非差无电离层组合伪距观测值先验方差和非差无电离层组合载波相位观测值先验方差是根据非差伪距观测值和载波相位观测值的先验标准方差和卫星高度角计算得到的。
根据法方程和所述误差方程,获取单差GNSS接收机钟差,并将任一基站的GNSS接收机时间作为参考时间,以使得各个基站的时间精密同步。
进一步地,在本发明实施例中,获取第m个历元的法方程:
Gm=NmX;
然后,探测载波相位观测值周跳,若为首历元,则不需探测周跳;进一步地,通过法方程相加,获取整体法方程,需要说明的是,法方程相加过程中,需要考虑前后两个历元待估参数的变化,若为首历元,则不需进行法方程相加。具体地,第m(m≠1)个历元的整体法方程为:
Gm=Gm-1+Gm;
Nm=Nm-1+Nm;
进一步地,对上述整体法方程进行求解,得到该方程的解为:
方程一协方差阵为:
从而获取到每个基站的GNSS接收机钟差之差。最后,以其中一个基站中的GNSS接收机钟为基准,根据GNSS接收机钟差之差,调整其他基站的GNSS接收机钟的时间,以使得所有基站间的时间同步。
在上述实施例的基础上,所述根据所述时差观测值,获取所述目标飞行器信号源的信号源坐标,以用于对所述目标飞行器信号源进行定位,包括:
步骤S1,根据接收目标飞行器信号源的基站对应的坐标信息,获取所述目标飞行器信号源的概略坐标(X0,Y0,Z0);
步骤S2,将任一接收目标飞行器信号源的基站k作为参考基站,根据时差观测值和基站d坐标信息,d,k∈K,d≠k,对于基站d构建观测方程:
其中,(xd,yd,zd)表示第d个基站的坐标信息,(x,y,z)表示目标飞行器信号源的坐标信息,c表示光速;Δtd表示当基站k为参考基站时,目标飞行器信号源发送的同一帧无线电信号到达第d个基站与参考基站k的时差;
步骤S3,根据接收目标飞行器信号源的基站总数K,构建K-1个上述观测方程,并根据泰勒级展开法和所述目标飞行器信号源的概略坐标(X0,Y0,Z0),对K-1个观测方程构建的方程组进行求解,当时,计算过程结束,得到所述目标飞行器信号源的信号源坐标(x,y,z),以用于对所述目标飞行器信号源进行定位;否则,令(x,y,z)=(X0,Y0,Z0),再次执行步骤S2至S3。
在本发明一实施例中,图3为本发明实施例提供的IGS跟踪站中PTBB-POTS接收机钟时间比对误差的示意图,可参考图3所示,从IGS服务器下载PTBB和POTS跟踪站在2019年第117天的数据,测站间距约为178公里,数据采样间隔是30秒。计算得到接收机钟差之差后,与IGS提供的事后精密钟差进行对比(IGS事后精密钟差的标称精度为0.075纳秒)。如图3所示,本发明实施例提供的方法,时间传递精度优于0.3纳秒,统计得到GNSS接收机钟差之差的均方根误差(Root Mean Square,简称RMS)为0.091纳秒。
进一步地,在本发明另一实施例中,通过以4个基站定位1架飞行器的模拟实验为例进行说明。在本发明实施例中,假设基站的坐标分别为:基站1(0km,-5km,0km),基站2(10km,0.1km,0.01km),基站3(20km,-0.1km,0.005km),基站4(30km,5km,0.002km);飞行器信号源坐标为(25km,0km,4km)。假定基站坐标误差为10米,基站间时间同步误差与到达时差测量误差的综合误差为30纳秒。在该条件下,随机模拟10000次。图4为本发明实施例提供的飞行器信号源定位结果的仿真示意图,可参考图4所示,本发明实施例提供的飞行器信号源的定位方法估计出目标飞行器的位置,和目标飞行器的实际位置基本一致;图5为本发明实施例提供的飞行器信号源的三维定位误差示意图,可参考图5所示,信号源三维定位精度基本优于250米。
图6为本发明实施例提供的基于GNSS精密时间同步的飞行器时差定位系统的结构示意图,如图6所示,本发明实施例提供了一种基于GNSS精密时间同步的飞行器时差定位系统,包括精密时间同步模块601、时差获取模块602和信号源定位模块603,其中,精密时间同步模块601基于GNSS精密时间,对多个基站进行时间同步处理,以使得所有基站之间的时间精密同步;时差获取模块602用于根据时间精密同步后的基站中GNSS接收机中的时间,获取目标飞行器信号源对应各个基站的时差观测值,所述时差观测值为所述目标飞行器信号源发送的同一帧无线电信号到达各个基站的时差;信号源定位模块603用于根据所述时差观测值,获取所述目标飞行器信号源的坐标,以用于对所述目标飞行器进行定位。
本发明实施例提供的基于GNSS精密时间同步的飞行器时差定位系统,该方法不需要原子钟,即可实现基站间精密时间同步,解决了时差定位系统中基站间精密时间同步的问题,进而实现高精度时差定位,提高了飞行器的定位精度,降低了GNSS星历误差的影响。
本发明实施例提供的系统是用于执行上述各方法实施例的,具体流程和详细内容请参照上述实施例,此处不再赘述。
图7为本发明实施例提供的电子设备结构示意图,参照图7,该电子设备可以包括:处理器(processor)701、通信接口(Communications Interface)702、存储器(memory)703和通信总线704,其中,处理器701,通信接口702,存储器703通过通信总线704完成相互间的通信。处理器701可以调用存储器703中的逻辑指令,以执行如下方法:基于GNSS精密时间,对多个基站进行时间同步处理,以使得所有基站之间的时间精密同步;根据时间精密同步后的基站中GNSS接收机中的时间,获取目标飞行器信号源对应各个基站的时差观测值,所述时差观测值为所述目标飞行器信号源发送的同一帧无线电信号到达各个基站的时差;根据所述时差观测值,获取所述目标飞行器信号源的坐标,以用于对所述目标飞行器进行定位。
此外,上述的存储器703中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
另一方面,本发明实施例还提供一种非暂态计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时实现以执行上述各实施例提供的基于GNSS精密时间同步的飞行器时差定位方法,例如包括:基于GNSS精密时间,对多个基站进行时间同步处理,以使得所有基站之间的时间精密同步;根据时间精密同步后的基站中GNSS接收机中的时间,获取目标飞行器信号源对应各个基站的时差观测值,所述时差观测值为所述目标飞行器信号源发送的同一帧无线电信号到达各个基站的时差;根据所述时差观测值,获取所述目标飞行器信号源的坐标,以用于对所述目标飞行器进行定位。
以上所描述的装置实施例仅仅是示意性的,其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下,即可以理解并实施。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件。基于这样的理解,上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在计算机可读存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (10)
1.一种基于GNSS精密时间同步的飞行器时差定位方法,其特征在于,包括:
基于GNSS精密时间,对多个基站进行时间同步处理,以使得所有基站之间的时间精密同步;
根据时间精密同步后的基站中GNSS接收机中的时间,获取目标飞行器信号源对应各个基站的时差观测值,所述时差观测值为所述目标飞行器信号源发送的同一帧无线电信号到达各个基站的时差;
根据所述时差观测值以及基站的已知坐标,获取所述目标飞行器信号源的坐标,以用于对所述目标飞行器进行定位。
2.根据权利要求1所述的基于GNSS精密时间同步的飞行器时差定位方法,其特征在于,所述基于GNSS精密时间,对多个基站进行时间同步处理,以使得每个基站之间的时间精密同步,包括:
获取每个基站的GNSS观测数据,并通过探测粗差,将GNSS观测数据中存在粗差的GNSS卫星剔除,以使得每个基站的GNSS观测数据来源于不存在粗差的GNSS卫星,所述GNSS观测数据包括双频伪距和载波相位观测值;
根据剔除粗差后每个基站中的GNSS观测数据,获取每两个基站之间的GNSS共视卫星,并根据所述GNSS共视卫星的观测数据,获取所述GNSS共视卫星的站间单差无电离层组合伪距观测值和站间单差无电离层组合载波相位观测值,以用于构建单差观测方程;
将单差GNSS接收机钟差、天顶对流层湿延迟和站间单差无电离层组合模糊度项作为参考估计,并根据所述站间单差无电离层组合伪距观测值和所述站间单差无电离层组合载波相位观测值,构建误差方程;
根据法方程和所述误差方程,获取单差GNSS接收机钟差,并将其中任一基站的GNSS接收机时间作为参考时间,以使得各个基站的时间精密同步。
3.根据权利要求2所述的基于GNSS精密时间同步的飞行器时差定位方法,其特征在于,在所述根据剔除粗差后每个基站中的GNSS观测数据,获取每两个基站之间的GNSS共视卫星,并根据所述GNSS共视卫星的观测数据,获取所述GNSS共视卫星的站间单差无电离层组合伪距观测值和站间单差无电离层组合载波相位观测值,以用于构建单差观测方程之后,所述方法还包括:
对所述站间单差无电离层组合伪距观测值和所述站间单差无电离层组合载波相位观测值中的误差进行修正,以根据误差修正后的站间单差无电离层组合伪距观测值和站间单差无电离层组合载波相位观测值,构建误差方程。
4.根据权利要求3所述的基于GNSS精密时间同步的飞行器时差定位方法,其特征在于,所述对所述站间单差无电离层组合伪距观测值和所述站间单差无电离层组合载波相位观测值中的误差进行修正,包括:
对所述站间单差无电离层组合伪距观测值和所述站间单差无电离层组合载波相位观测值中的地球自转效应误差、卫星天线相位中心偏差、卫星天线相位中心变化、接收机天线相位中心偏差、接收机天线相位中心变化、相对论效应、引力延迟、相位缠绕和对流层干延迟进行修正,得到误差修正后的站间单差无电离层组合伪距观测值和站间单差无电离层组合载波相位观测值。
5.根据权利要求2所述的基于GNSS精密时间同步的飞行器时差定位方法,其特征在于,所述单差观测方程为:
其中,表示第i个GNSS卫星下基站a和基站b的站间单差无电离层组合伪距观测值;表示第i个GNSS卫星下基站a和基站b的站间单差卫地距;tab表示基站a和基站b之间的单差接收机钟差;表示第i个GNSS卫星下基站a和基站b的站间单差对流层延迟;表示第i个GNSS卫星下基站a和基站b的站间单差无电离层组合载波相位观测值,单位为米;表示第i个GNSS卫星下基站a和基站b的单差无电离层组合模糊度项,单位为米;a和b表示GNSS接收机标识,i表示GNSS卫星标识;
其中,表示第i个GNSS卫星下基站a的L1载波的伪距观测值;表示第i个GNSS卫星下基站b的L1载波的伪距观测值;表示第i个GNSS卫星下基站a的L2载波的伪距观测值;表示第i个GNSS卫星下基站b的L2载波的伪距观测值;表示第i个GNSS卫星下基站a的L1载波的载波相位观测值;表示第i个GNSS卫星下基站b的L1载波的载波相位观测值;表示第i个GNSS卫星下基站a的L2载波的载波相位观测值;表示第i个GNSS卫星下基站b的L2载波的载波相位观测值;ρ表示卫地距;λ1和λ2分别表示L1载波和L2载波的波长;表示第i个GNSS卫星下基站a的L1载波的模糊度;表示第i个GNSS卫星下基站b的L1载波的模糊度;表示第i个GNSS卫星下基站a的L2载波的模糊度;表示第i个GNSS卫星下基站b的L2载波的模糊度;表示第i个GNSS卫星下基站a的对流层干延迟投影函数;表示第i个GNSS卫星下基站b的对流层干延迟投影函数;表示第i个GNSS卫星下基站a的对流层湿延迟投影函数;表示第i个GNSS卫星下基站b的对流层湿延迟投影函数;TD,a表示基站a的天顶对流层干延迟;TD,b表示基站b的天顶对流层干延迟;TW,a表示基站a的天顶对流层湿延迟;TW,b表示基站b的天顶对流层湿延迟;f1表示L1载波频率,f2表示L2载波频率,(xa,ya,za)表示基站a的GNSS接收机坐标;(xb,yb,zb)表示基站b的GNSS接收机坐标;表示基站a对应的第i个GNSS卫星的坐标;表示基站b对应的第i个GNSS卫星的坐标。
6.根据权利要求5所述的基于GNSS精密时间同步的飞行器时差定位方法,其特征在于,所述误差方程为:
V=BX-l, W;
其中,B表示设计矩阵,X表示待估参数向量,l表示观测值向量,V表示观测值残差向量,W表示观测值向量权阵,所述观测值向量权阵根据观测值先验精度获取得到的;
7.根据权利要求6所述的基于GNSS精密时间同步的飞行器时差定位方法,其特征在于,所述根据所述时差观测值,获取所述目标飞行器信号源的信号源坐标,以用于对所述目标飞行器信号源进行定位,包括:
根据接收目标飞行器信号源的基站对应的坐标信息,获取所述目标飞行器信号源的概略坐标;
将任一接收目标飞行器信号源的基站k作为参考基站,根据时差观测值和基站d坐标信息,d,k∈K,d≠k,构建观测方程:
其中,(xd,yd,zd)表示第d个基站的坐标信息,(x,y,z)表示目标飞行器信号源的坐标信息,c表示光速;Δtd表示当基站k为参考基站时,目标飞行器信号源发送的同一帧无线电信号到达第d个基站与参考基站k的时差;
根据接收目标飞行器信号源的基站总数K,构建K-1个观测方程,并根据泰勒级展开法和所述目标飞行器信号源的概略坐标,对K-1个观测方程进行求解,得到所述目标飞行器信号源的信号源坐标,以用于对所述目标飞行器信号源进行定位。
8.一种基于GNSS精密时间同步的飞行器时差定位系统,其特征在于,包括:
精密时间同步模块,基于GNSS精密时间,对多个基站进行时间同步处理,以使得所有基站之间的时间精密同步;
时差获取模块,用于根据时间精密同步后的基站中GNSS接收机中的时间,获取目标飞行器信号源对应各个基站的时差观测值,所述时差观测值为所述目标飞行器信号源发送的同一帧无线电信号到达各个基站的时差;
信号源定位模块,用于根据所述时差观测值,获取所述目标飞行器信号源的坐标,以用于对所述目标飞行器进行定位。
9.一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时实现如权利要求1至7任一项所述基于GNSS精密时间同步的飞行器时差定位方法的步骤。
10.一种非暂态计算机可读存储介质,其上存储有计算机程序,其特征在于,该计算机程序被处理器执行时实现如权利要求1至7任一项所述基于GNSS精密时间同步的飞行器时差定位方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911276822.1A CN111044972B (zh) | 2019-12-12 | 2019-12-12 | 基于gnss精密时间同步的飞行器时差定位方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911276822.1A CN111044972B (zh) | 2019-12-12 | 2019-12-12 | 基于gnss精密时间同步的飞行器时差定位方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111044972A true CN111044972A (zh) | 2020-04-21 |
CN111044972B CN111044972B (zh) | 2021-01-05 |
Family
ID=70235912
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911276822.1A Active CN111044972B (zh) | 2019-12-12 | 2019-12-12 | 基于gnss精密时间同步的飞行器时差定位方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111044972B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111983650A (zh) * | 2020-07-07 | 2020-11-24 | 北京自动化控制设备研究所 | 一种基于gnss的高精度时间传递方法 |
CN112284195A (zh) * | 2020-09-19 | 2021-01-29 | 中国航空工业集团公司洛阳电光设备研究所 | 一种基于相对卫星导航的对地精确打击设计方法 |
CN113625223A (zh) * | 2020-05-08 | 2021-11-09 | 大唐高鸿数据网络技术股份有限公司 | 一种定位方法及终端设备 |
CN115933356A (zh) * | 2023-01-09 | 2023-04-07 | 北京航空航天大学 | 一种虚拟原子钟的高精度时间同步系统和方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020183070A1 (en) * | 2001-05-30 | 2002-12-05 | L. Scott Bloebaum | Method and system for GPS bit-edge synchronization in the presence of burst-mode interference |
CN102165332A (zh) * | 2008-07-24 | 2011-08-24 | Ses阿斯特拉有限公司 | 空间飞行器位置估计系统和方法 |
CN202421491U (zh) * | 2012-02-03 | 2012-09-05 | 中国电子科技集团公司第三十八研究所 | 基于分布式时钟的多边定位系统 |
CN108445518A (zh) * | 2018-03-16 | 2018-08-24 | 中国科学院数学与系统科学研究院 | 一种基于双差模糊度固定解约束的gnss精密时间传递方法 |
CN108519586A (zh) * | 2018-04-03 | 2018-09-11 | 芜湖泰贺知信息系统有限公司 | 一种分布式无源雷达系统及其目标定位方法 |
CN109655846A (zh) * | 2019-01-30 | 2019-04-19 | 清华大学 | 一种多站差分后处理高精度时间同步方法及系统 |
CN110471087A (zh) * | 2019-08-13 | 2019-11-19 | 北京电子工程总体研究所 | 一种空间飞行器的时间漂移计算方法和系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107153209B (zh) * | 2017-07-06 | 2019-07-30 | 武汉大学 | 一种短弧段低轨导航卫星实时精密定轨方法 |
-
2019
- 2019-12-12 CN CN201911276822.1A patent/CN111044972B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020183070A1 (en) * | 2001-05-30 | 2002-12-05 | L. Scott Bloebaum | Method and system for GPS bit-edge synchronization in the presence of burst-mode interference |
CN102165332A (zh) * | 2008-07-24 | 2011-08-24 | Ses阿斯特拉有限公司 | 空间飞行器位置估计系统和方法 |
CN202421491U (zh) * | 2012-02-03 | 2012-09-05 | 中国电子科技集团公司第三十八研究所 | 基于分布式时钟的多边定位系统 |
CN108445518A (zh) * | 2018-03-16 | 2018-08-24 | 中国科学院数学与系统科学研究院 | 一种基于双差模糊度固定解约束的gnss精密时间传递方法 |
CN108519586A (zh) * | 2018-04-03 | 2018-09-11 | 芜湖泰贺知信息系统有限公司 | 一种分布式无源雷达系统及其目标定位方法 |
CN109655846A (zh) * | 2019-01-30 | 2019-04-19 | 清华大学 | 一种多站差分后处理高精度时间同步方法及系统 |
CN110471087A (zh) * | 2019-08-13 | 2019-11-19 | 北京电子工程总体研究所 | 一种空间飞行器的时间漂移计算方法和系统 |
Non-Patent Citations (2)
Title |
---|
张建立 等: "时差定位系统精度分析及工程应用", 《电子技术与软件工程》 * |
解楠: "多站时差定位精度分析", 《信息与电子工程》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113625223A (zh) * | 2020-05-08 | 2021-11-09 | 大唐高鸿数据网络技术股份有限公司 | 一种定位方法及终端设备 |
CN113625223B (zh) * | 2020-05-08 | 2024-04-30 | 中信科智联科技有限公司 | 一种定位方法及终端设备 |
CN111983650A (zh) * | 2020-07-07 | 2020-11-24 | 北京自动化控制设备研究所 | 一种基于gnss的高精度时间传递方法 |
CN111983650B (zh) * | 2020-07-07 | 2022-07-15 | 北京自动化控制设备研究所 | 一种基于gnss的高精度时间传递方法 |
CN112284195A (zh) * | 2020-09-19 | 2021-01-29 | 中国航空工业集团公司洛阳电光设备研究所 | 一种基于相对卫星导航的对地精确打击设计方法 |
CN112284195B (zh) * | 2020-09-19 | 2022-11-22 | 中国航空工业集团公司洛阳电光设备研究所 | 一种基于相对卫星导航的对地精确打击设计方法 |
CN115933356A (zh) * | 2023-01-09 | 2023-04-07 | 北京航空航天大学 | 一种虚拟原子钟的高精度时间同步系统和方法 |
CN115933356B (zh) * | 2023-01-09 | 2023-08-22 | 北京航空航天大学 | 一种虚拟原子钟的高精度时间同步系统和方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111044972B (zh) | 2021-01-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111044972B (zh) | 基于gnss精密时间同步的飞行器时差定位方法及系统 | |
CN111045034B (zh) | 基于广播星历的gnss多系统实时精密时间传递方法及系统 | |
CN108931915B (zh) | 利用导航卫星的授时方法和装置、计算机可读存储介质 | |
CN109655846B (zh) | 一种多站差分后处理高精度时间同步方法及系统 | |
CN109581452B (zh) | 一种gnss参考站载波相位整周模糊度解算方法 | |
CN108226985B (zh) | 基于精密单点定位的列车组合导航方法 | |
CN108415049A (zh) | 提高网络rtk双差宽巷模糊度固定正确率的方法 | |
US7576690B2 (en) | Position determination with reference data outage | |
CN108196272A (zh) | 一种基于实时精密单点定位的卫星导航定位装置及方法 | |
EP3109672B1 (en) | Gnss receiver with a capability to resolve ambiguities using an uncombined formulation | |
CN108254773A (zh) | 一种多gnss的实时钟差解算方法 | |
CN105589087A (zh) | 一种高精度实时卫星定位装置及其方法 | |
GB2487256A (en) | Location fix from an unknown position | |
CN110412634B (zh) | 基于参考站网的伪距差分定位方法及装置 | |
CN111381264A (zh) | 网络rtk中长基线模糊度固定方法和平台 | |
CN112146557A (zh) | 一种基于gnss的实时桥梁变形监测系统及方法 | |
CN103235321A (zh) | Gps单站伪距定位精密授时方法 | |
CN111290004A (zh) | 伪距差分定位方法、装置、电子设备及存储介质 | |
CN117388881B (zh) | 一种低轨卫星的星载原子钟向UTC(k)的溯源方法及系统 | |
CN112462399A (zh) | 一种基于接收机的伪距偏差计算及改正的方法及装置 | |
CN114994727A (zh) | 一种实现高精度时间校准及卫星定位的设备 | |
CN110068848B (zh) | 一种高性能rtk处理技术方法 | |
CN105549048A (zh) | 基于三角形接收机阵列的高精度卫星定位装置及其方法 | |
CN111983650A (zh) | 一种基于gnss的高精度时间传递方法 | |
CN111913199B (zh) | 基于变分模态分解的手机gnss数据噪声提取方法 |
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 |