CN109001770A - 多模卫星导航系统用户层几何精度因子最小值计算方法 - Google Patents
多模卫星导航系统用户层几何精度因子最小值计算方法 Download PDFInfo
- Publication number
- CN109001770A CN109001770A CN201810545454.5A CN201810545454A CN109001770A CN 109001770 A CN109001770 A CN 109001770A CN 201810545454 A CN201810545454 A CN 201810545454A CN 109001770 A CN109001770 A CN 109001770A
- Authority
- CN
- China
- Prior art keywords
- matrix
- navigation system
- satellite navigation
- layer geometric
- precision
- 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
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
- 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/23—Testing, 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)
- Complex Calculations (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Navigation (AREA)
Abstract
本发明公开了一种多模卫星导航系统用户层几何精度因子最小值计算方法,其包括构建多模卫星导航系统中用户层几何矩阵,采用分块矩阵求逆方法对相关矩阵进行处理,构建用户层几何精度因子计算模型,分别计算相关部分的对角元素,构建多模卫星导航系统用户层几何精度因子最小值计算模型。本发明将分块矩阵求逆方法应用于用户层几何精度因子最小值计算,根据几何矩阵对测量矩阵进行分块,再通过分块矩阵求逆方法和伴随矩阵计算相关部分的对角元素,最终获得多模卫星导航系统用户层几何精度因子最小值。
Description
技术领域
本发明属于卫星导航技术领域,具体涉及一种多模卫星导航系统用户层几何精度因子最小值计算方法。
背景技术
全球导航卫星系统(Global Navigation Satellite System,GNSS)广泛应用于载体导航、定位以及授时等诸多领域。随着卫星导航系统的不断完善与发展,将多个单系统组合起来构成多模卫星导航系统,可以有效提高系统定位精度、连续性、可用性以及完好性,从而改善导航定位性能。
在多模卫星导航系统中,几何精度因子(Geometric Dilution of Precision,GDOP)对于选星、定位精度评估以及系统病态性诊断等方面均具有重要意义。通常情况下,GDOP值越小则定位精度越高。因此,计算多模卫星导航系统GDOP最小值具有重要意义。
在多模卫星导航系统中,针对系统之间的时间偏差有不同的处理方法,从而导致几何矩阵也有不同的形式,进而造成GDOP最小值计算方法也有明显区别。在多模卫星导航系统中,由于用户层几何矩阵更为复杂,因此传统方法则不再适用于计算用户层GDOP最小值。
发明内容
本发明的发明目的是:为了解决利用传统方法无法正确计算多模卫星导航系统用户层GDOP最小值的问题,本发明提出了一种多模卫星导航系统用户层几何精度因子最小值计算方法。
本发明的技术方案是:一种多模卫星导航系统用户层几何精度因子最小值计算方法,包括以下步骤:
A、在多模卫星导航系统中构建用户层几何矩阵HU;
B、计算测量矩阵并对其进行分块处理,采用分块矩阵求逆方法对进行求逆,得到表示为:
C、构建用户层几何精度因子计算模型,表示为:
其中,GDOPU为用户层几何精度因子,tr(·)表示矩阵求迹。
D、分别计算tr(φ-1)和tr(ψ-1)的最小值;
E、根据步骤D中tr(φ-1)和tr(ψ-1)的最小值和步骤C中用户层几何精度因子计算模型,得到多模卫星导航系统用户层几何精度因子最小值计算模型,表示为:
其中,ni表示第i个子系统的卫星数目,1≤i≤3,n=n1+n2+n3。
进一步地,所述步骤A中用户层几何矩阵HU表示为:
其中,i=1,2,3表示第i个子系统,1i与0i分别表示单位向量和零向量;Hi为第i个子系统由ni颗卫星的方向余弦向量组成的矩阵。
进一步地,所述用户层几何矩阵中
其中,hj表示第j颗卫星的方向余弦向量,1≤j≤n1+n2+n3。
进一步地,所述步骤B中测量矩阵表示为:
进一步地,所述步骤B中对测量矩阵进行分块处理,得到:
其中,
进一步地,所述步骤B中对分块矩阵进行求逆,得到表示为:
进一步地,所述步骤D中tr(φ-1)最小值计算模型表示为:
进一步地,所述步骤D中tr(ψ-1)最小值计算模型表示为:
本发明的有益效果是:本发明将分块矩阵求逆方法应用于用户层几何精度因子最小值计算,根据几何矩阵对测量矩阵进行分块,再通过分块矩阵求逆方法和伴随矩阵计算相关部分的对角元素,最终获得多模卫星导航系统用户层几何精度因子最小值。
附图说明
图1是本发明的多模卫星导航系统用户层几何精度因子最小值计算方法的流程示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
如图1所示,为本发明的多模卫星导航系统用户层几何精度因子最小值计算方法的流程示意图。一种多模卫星导航系统用户层几何精度因子最小值计算方法,包括以下步骤:
A、在多模卫星导航系统中构建用户层几何矩阵HU;
B、计算测量矩阵并对其进行分块处理,采用分块矩阵求逆方法对进行求逆,得到表示为:
C、构建用户层几何精度因子计算模型,表示为:
其中,GDOPU为用户层几何精度因子,tr(·)表示矩阵求迹。
D、分别计算tr(φ-1)和tr(ψ-1)的最小值;
E、根据步骤D中tr(φ-1)和tr(ψ-1)的最小值和步骤C中用户层几何精度因子计算模型,得到多模卫星导航系统用户层几何精度因子最小值计算模型,表示为:
其中,ni表示第i个子系统的卫星数目,1≤i≤3,n=n1+n2+n3。
在本发明的一个可选实施例中,上述步骤A中用户层几何矩阵HU表示为:
其中,i=1,2,3表示第i个子系统,1i与0i分别表示单位向量和零向量;Hi为第i个子系统由ni颗卫星的方向余弦向量组成的矩阵,具体表示为
其中,ni表示第i个子系统的卫星数目,1≤i≤3,n=n1+n2+n3。
在本发明的一个可选实施例中,上述步骤B中计算测量矩阵表示为:
对测量矩阵进行分块处理,得到:
其中,
再根据分块矩阵求逆方法对分块矩阵进行求逆,得到表示为:
其中,非对角元素对用户层几何精度因子的计算没有影响。
在本发明的一个可选实施例中,上述步骤C构建多模卫星导航系统中用户层几何精度因子的计算模型,表示为:
再根据步骤B中得到的将用户层几何精度因子的计算模型表示为:
从而通过分别计算tr(φ-1)和tr(ψ-1)的最小值,可以获得多模卫星导航系统用户层几何精度因子最小值。
在本发明的一个可选实施例中,上述步骤D中tr(φ-1)最小值计算模型为
tr(φ-1)=tr[(A-UD-1UT)-1]
其中,矩阵A和φ都为对称的非负定矩阵,且A≥φ,则得到
tr[(A-UD-1UT)-1]≥tr(A-1)
当且仅当U=0时,等号成立。
设定λi(A)与aii分别表示矩阵A的特征值与对角元素,则根据矩阵特征值与对角元素之间的关系,可以得到
根据调和不等式,可以得到
当且仅当λ1(A)=λ2(A)=λ3(A)=n/3时,等号成立。
从而将tr(φ-1)最小值计算模型表示为:
得到tr(ψ-1)最小值计算模型为
tr(ψ-1)=tr[(D-UTA-1U)-1]≥tr(D-1)
当且仅当UT=0时,等号成立。
再通过计算矩阵D-1的对角元素计算tr(D-1),令D*与|D|分别矩阵D的伴随矩阵与行列式,得到
其中,行列式伴随矩阵D*的三个对角元素分别为:
得到tr(D-1)为
从而将tr(ψ-1)最小值计算模型表示为:
当且仅当UT=0时,tr(D-1)取得最小值,即tr(ψ-1)也取得最小值。
在本发明的一个可选实施例中,上述步骤E根据步骤D中tr(φ-1)和tr(ψ-1)的最小值和步骤C中用户层几何精度因子计算模型,将用户层几何精度因子计算模型表示为
当且仅当U=UT=0时,等号成立。
从而得到多模卫星导航系统用户层几何精度因子最小值计算模型,表示为:
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。
Claims (8)
1.一种多模卫星导航系统用户层几何精度因子最小值计算方法,其特征在于,包括以下步骤:
A、在多模卫星导航系统中构建用户层几何矩阵HU;
B、计算测量矩阵并对其进行分块处理,采用分块矩阵求逆方法对进行求逆,得到表示为:
C、构建用户层几何精度因子计算模型,表示为:
其中,GDOPU为用户层几何精度因子,tr(·)表示矩阵求迹。
D、分别计算tr(φ-1)和tr(ψ-1)的最小值;
E、根据步骤D中tr(φ-1)和tr(ψ-1)的最小值和步骤C中用户层几何精度因子计算模型,得到多模卫星导航系统用户层几何精度因子最小值计算模型,表示为:
其中,ni表示第i个子系统的卫星数目,1≤i≤3,n=n1+n2+n3。
2.如权利要求1所述的多模卫星导航系统用户层几何精度因子最小值计算方法,其特征在于,所述步骤A中用户层几何矩阵HU表示为:
其中,i=1,2,3表示第i个子系统,1i与0i分别表示单位向量和零向量;Hi为第i个子系统由ni颗卫星的方向余弦向量组成的矩阵。
3.如权利要求2所述的多模卫星导航系统用户层几何精度因子最小值计算方法,其特征在于,所述用户层几何矩阵中
其中,hj表示第j颗卫星的方向余弦向量,1≤j≤n1+n2+n3。
4.如权利要求3所述的多模卫星导航系统用户层几何精度因子最小值计算方法,其特征在于,所述步骤B中测量矩阵表示为:
5.如权利要求4所述的多模卫星导航系统用户层几何精度因子最小值计算方法,其特征在于,所述步骤B中对测量矩阵进行分块处理,得到:
其中,
6.如权利要求5所述的多模卫星导航系统用户层几何精度因子最小值计算方法,其特征在于,所述步骤B中对分块矩阵进行求逆,得到表示为:
7.如权利要求6所述的多模卫星导航系统用户层几何精度因子最小值计算方法,其特征在于,所述步骤D中tr(φ-1)最小值计算模型表示为:
8.如权利要求7所述的多模卫星导航系统用户层几何精度因子最小值计算方法,其特征在于,所述步骤D中tr(ψ-1)最小值计算模型表示为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810545454.5A CN109001770B (zh) | 2018-05-25 | 2018-05-25 | 多模卫星导航系统用户层几何精度因子最小值计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810545454.5A CN109001770B (zh) | 2018-05-25 | 2018-05-25 | 多模卫星导航系统用户层几何精度因子最小值计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109001770A true CN109001770A (zh) | 2018-12-14 |
CN109001770B CN109001770B (zh) | 2021-03-30 |
Family
ID=64573327
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810545454.5A Active CN109001770B (zh) | 2018-05-25 | 2018-05-25 | 多模卫星导航系统用户层几何精度因子最小值计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109001770B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110045396A (zh) * | 2019-05-14 | 2019-07-23 | 成都海亨利科技有限公司 | 一种bd/gps双模接收机gdop最小值计算方法 |
CN111174786A (zh) * | 2020-02-25 | 2020-05-19 | 中南大学 | 一种ins/sar组合导航的几何精度因子计算方法 |
CN113589346A (zh) * | 2021-08-04 | 2021-11-02 | 电子科技大学 | 一种约束条件下系统层与用户层gdop最小值计算方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101149428A (zh) * | 2007-11-15 | 2008-03-26 | 北京航空航天大学 | 一种组合卫星导航系统的快速选星方法 |
CN101571581A (zh) * | 2009-06-15 | 2009-11-04 | 北京航空航天大学 | 全球导航卫星系统快速选星方法 |
CN103592658A (zh) * | 2013-09-30 | 2014-02-19 | 北京大学 | 多模卫星导航系统中基于选星算法的raim新方法 |
CN105589078A (zh) * | 2015-12-11 | 2016-05-18 | 电子科技大学 | 一种双星座组合导航系统几何精度因子计算方法 |
US20170090016A1 (en) * | 2015-09-25 | 2017-03-30 | Texas Instruments Incorporated | Method for Joint Antenna-Array Calibration and Direction of Arrival Estimation for Automotive Applications |
KR101737802B1 (ko) * | 2016-04-12 | 2017-05-19 | 국방과학연구소 | Gps 항재밍을 위한 배열안테나의 시공간 가중치 결정 장치 및 이를 이용한 방법 |
-
2018
- 2018-05-25 CN CN201810545454.5A patent/CN109001770B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101149428A (zh) * | 2007-11-15 | 2008-03-26 | 北京航空航天大学 | 一种组合卫星导航系统的快速选星方法 |
CN101571581A (zh) * | 2009-06-15 | 2009-11-04 | 北京航空航天大学 | 全球导航卫星系统快速选星方法 |
CN103592658A (zh) * | 2013-09-30 | 2014-02-19 | 北京大学 | 多模卫星导航系统中基于选星算法的raim新方法 |
US20170090016A1 (en) * | 2015-09-25 | 2017-03-30 | Texas Instruments Incorporated | Method for Joint Antenna-Array Calibration and Direction of Arrival Estimation for Automotive Applications |
CN105589078A (zh) * | 2015-12-11 | 2016-05-18 | 电子科技大学 | 一种双星座组合导航系统几何精度因子计算方法 |
KR101737802B1 (ko) * | 2016-04-12 | 2017-05-19 | 국방과학연구소 | Gps 항재밍을 위한 배열안테나의 시공간 가중치 결정 장치 및 이를 이용한 방법 |
Non-Patent Citations (4)
Title |
---|
刘洋等: "基于GDOP分析模型的双星无源定位系统精度研究", 《南昌航空大学学报:自然科学版》 * |
应士君等: "BD-2/GPS组合导航系统选星算法及定位分析", 《上海海事大学学报》 * |
杨薛涛等: "北斗卫星导航系统在高纬度地区GDOP最小值分析", 《第八届中国卫星导航学术年会论文集—S02导航与位置服务》 * |
温忠麟等: "《调节效应和中介效应分析》", 30 March 2012 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110045396A (zh) * | 2019-05-14 | 2019-07-23 | 成都海亨利科技有限公司 | 一种bd/gps双模接收机gdop最小值计算方法 |
CN111174786A (zh) * | 2020-02-25 | 2020-05-19 | 中南大学 | 一种ins/sar组合导航的几何精度因子计算方法 |
CN113589346A (zh) * | 2021-08-04 | 2021-11-02 | 电子科技大学 | 一种约束条件下系统层与用户层gdop最小值计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109001770B (zh) | 2021-03-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109001770A (zh) | 多模卫星导航系统用户层几何精度因子最小值计算方法 | |
CN103925925B (zh) | 一种用于多点定位系统的实时高精度位置解算方法 | |
CN111707275B (zh) | 一种定位方法、装置、电子设备及计算机可读存储介质 | |
CN103369466A (zh) | 一种地图匹配辅助室内定位方法 | |
US10060743B2 (en) | Method and apparatus of correcting output value of geomagnetic sensor | |
CN111366156A (zh) | 基于神经网络辅助的变电站巡检机器人导航方法及系统 | |
Xue et al. | Understanding GDOP minimization in GNSS positioning: Infinite solutions, finite solutions and no solution | |
CN109444930A (zh) | 一种基于分步加权最小二乘估计的单点定位的方法及装置 | |
CN105895089A (zh) | 一种语音识别方法及装置 | |
CN111103572B (zh) | 基于分布目标的星载sar极化定标方法及设备 | |
CN105589078A (zh) | 一种双星座组合导航系统几何精度因子计算方法 | |
US6571173B1 (en) | Three-dimensional space curve comparison using spatial angle variance metric and applications thereof | |
CN112822633B (zh) | 一种基于误差补偿的定位方法,系统,终端设备及可读存储介质 | |
CN113466909A (zh) | 一种gnss多频系统部分整周模糊度子集选取方法 | |
EP2879011B1 (en) | On-board estimation of the nadir attitude of an Earth orbiting spacecraft | |
CN108919315B (zh) | 多模接收机系统层几何精度因子最小值计算方法 | |
CN116413758A (zh) | 一种无线电信号辅助城市复杂环境卫星定位的方法 | |
CN112083463A (zh) | 模糊度是否固定正确的检测方法及装置、定位终端 | |
CN106990423B (zh) | 一种卫星导航定位整周模糊度有效性检验方法 | |
CN113870600B (zh) | 车道线信息显示方法、装置、电子设备和计算机可读介质 | |
CN104680564A (zh) | 一种灰度增强层析piv重构方法、装置和设备 | |
US10109060B1 (en) | Uncertainty prediction in passive target motion analysis | |
US11543257B2 (en) | Navigation apparatus and operation method of navigation apparatus | |
CN110045396B (zh) | 一种bd/gps双模接收机gdop最小值计算方法 | |
CN115326081A (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 |