CN115598686A - 一种基于伪卫星信号的测向方法 - Google Patents
一种基于伪卫星信号的测向方法 Download PDFInfo
- Publication number
- CN115598686A CN115598686A CN202211278195.7A CN202211278195A CN115598686A CN 115598686 A CN115598686 A CN 115598686A CN 202211278195 A CN202211278195 A CN 202211278195A CN 115598686 A CN115598686 A CN 115598686A
- Authority
- CN
- China
- Prior art keywords
- pseudolite
- double
- pseudo
- difference
- matrix
- 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.)
- Pending
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
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/53—Determining attitude
- G01S19/54—Determining attitude using carrier phase measurements; using long or short baseline interferometry
- G01S19/55—Carrier phase ambiguity resolution; Floating ambiguity; LAMBDA [Least-squares AMBiguity Decorrelation Adjustment] method
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)
Abstract
本发明提供一种基于伪卫星信号的测向方法,包括如下步骤,利用多颗伪卫星发送模拟信号;利用接收机接收伪卫星观测数据,包括伪距、载波相位;对伪卫星观测数据进行周跳检测、载波相位历元差计算以及根据对应的数据统计特征进行可用星判别,去除不可用伪卫星;所述接收机再获取伪卫星的载波相位、伪距;根据接收机获取的伪卫星载波相位、伪距,采用基于MC‑LAMBDA算法的测向方法,得到多基线解向量。本发明增加了可用星的数目,即增加了有效观测量,从而有利于测向求解;伪卫星一个显著的特点就是其高度角很低,且信号不需要通过电离层,将这种低高度角伪卫星与GPS组合定位能够有效地改善几何图形结构,提高测向精度。
Description
技术领域
本发明属于伪卫星测向领域,尤其是涉及一种基于伪卫星信号的测向方法。
背景技术
对于GPS全球卫星定位系统,其定位精度、可用性和可靠性都依赖于所观测到的卫星数目和卫星星座的几何图形结构。在观测条件不理想的情况下,所能观测到的卫星数目和卫星的几何图形结构通常都不理想,也难以满足精密定位的需要。在某些极端条件下,如在室内或地下、隧洞中,则完全接收不到卫星信号,GPS也就无法工作;虽然GPS测量在水平方向上能达到较高的精度,但在垂直方向上定位精度较差,其误差通常为水平定位误差的2~3倍,难以满足一些精密应用的要求。
发明内容
有鉴于此,本发明提供了一种基于伪卫星信号的测向方法,旨在解决在GPS观测中低高度角卫星不足的问题,提高在室内、地下等卫星不足的情况下的测向精度。
为达到上述目的,本发明的技术方案是这样实现的:
一种基于伪卫星信号的测向方法,包括如下步骤
步骤1:利用多颗伪卫星发送模拟信号;
步骤2:利用接收机接收伪卫星观测数据,包括伪距、载波相位;
步骤3:对伪卫星观测数据进行周跳检测、载波相位历元差计算以及根据对应的数据统计特征进行可用星判别,去除不可用伪卫星;
步骤4:所述接收机再获取伪卫星的载波相位、伪距;
步骤5:根据接收机获取的伪卫星载波相位、伪距,采用基于MC-LAMBDA算法的测向方法,得到多基线解向量。
进一步的,所述多颗伪卫星中的每一颗伪卫星根据卫星参数、伪卫星位置和参考时延,模拟卫星的卫星信号并发送伪卫星信号。
进一步的,所述伪卫星观测数据还包括载噪比测量值、多普勒频移测量值。
进一步的,所述步骤5具体包括如下步骤:
步骤51:设某载体上一短基线同时观测两颗伪卫星s1、s2,则双差载波相位和伪距观测方程为:
其中,为双差载波相位观测量;为双差伪距观测值;λ为载波波长;z为双差整周模糊度;载体上天线之间的距离很短,假设同一历元载体上所有天线接收中心到同一颗伪卫星的视向量相同,记为s12,即为s12的转置;bE为基线向量;ε、e分别为双差载波相位和伪距观测噪声;
步骤52:载体上安装m+1个天线,共视伪卫星数为n+1,选择任一天线作为主天线,仰角最高的伪卫星作为参考星(以符号r表示),则基线i的观测方程表示为:
其中,S表示伪卫星视向量;其中,j=1,…,n;A为包含载波波长的设计矩阵;为某一伪卫星相对于仰角最高的伪卫星的双差伪距观测值;为某一伪卫星相对于仰角最高的伪卫星的双差载波相位观测值;v为量测噪声矩阵;D(yi)为yi的方差;为yi的协方差矩阵;zi为双差模糊度向量其中,i=1,…,m;
步骤53:将载体上m条基线观测量组合,得到多基线观测方程:
其中BB表示为BE在载体坐标系中相应的向量,RBB=BE,R为姿态旋转矩阵,满足RTR=I,I表示单位阵;QY为Y的协方差矩阵;Z为双差整周模糊度矩阵;
步骤54:引入向量算子vec,对式(4)进行向量化得到:
vec(Y)=vec(AZ)+vec(SRBB)+vec(V) (5)
根据矩阵乘积向量化公式,由式(5)可以得到:
步骤55:MC-LAMBDA方法解算整周模糊度及姿态矩阵即为解算下面最小化问题:
其中,为双差整周模糊度浮点解的协方差矩阵,vecZ为双差整周模糊度矩阵固定解;分别为双差整周模糊度浮点解和姿态旋转矩阵浮点解及其协方差矩阵,其中,为双差整周模糊度浮点解的协方差矩阵,为双差整周模糊度浮点解和姿态矩阵浮点解的协方差矩阵,为姿态矩阵浮点解和双差整周模糊度浮点解的协方差矩阵,为姿态矩阵浮点解的协方差矩阵;分别为姿态矩阵固定解及其协方差矩阵;
步骤56:采用LAMBDA方法搜索满足下式的所有整数值z:
从而解算整周模糊度;
步骤57:得到多基线解向量,继而得到航向角和俯仰角。
相对于现有技术,本发明所述的一种基于伪卫星信号的测向方法具有以下优势:
1、本发明增加了可用星的数目,即增加了有效观测量,从而有利于提高在室内或地下、隧洞中等卫星信号较弱或没有卫星信号情况下测向的精确度;
2、伪卫星一个显著的特点就是其高度角很低,且信号不需要通过电离层,将这种低高度角伪卫星与GPS组合定位能够有效地改善几何图形结构,提高测向精度;
3、伪卫星可以与GPS组合进行测向,增加了观测量,有效地改善几何图形结构,附加的观测量还有利于增强GPS模糊度的解算,提高测向的精度以及测向的可靠性。
附图说明
构成本发明的一部分的附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1为本发明的基于伪卫星信号的测向方法流程图;
图2为本发明的基于MC-LAMBDA算法的测向流程图。
具体实施方式
需要说明的是,在不冲突的情况下,本发明中的实施例及实施例中的特征可以相互组合。
在本发明的描述中,需要理解的是,术语“中心”、“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,术语“第一”、“第二”等仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”等的特征可以明示或者隐含地包括一个或者更多个该特征。在本发明的描述中,除非另有说明,“多个”的含义是两个或两个以上。
在本发明的描述中,需要说明的是,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以通过具体情况理解上述术语在本发明中的具体含义。
下面将参考附图并结合实施例来详细说明本发明。
如图1所示,本发明提供的一种基于伪卫星信号的测向方法包括按顺序进行的下列步骤:
1)发送单元的多颗伪卫星中的每一颗根据卫星参数、伪卫星位置和参考时延,模拟卫星的卫星信号并发送伪卫星信号,伪卫星可以模拟GNSS卫星信号,例如GPS信号或者北斗信号,并将模拟的卫星信号通过射频天线播发到终端。
2)接收处理单元的接收机接收伪卫星信号,获取所有的伪卫星观测数据,例如伪距、载波相位、载噪比测量值、多普勒频移测量值。
3)接收处理单元的处理器,对伪卫星观测数据进行周跳检测、载波相位历元差计算以及根据对应的数据统计特征进行可用星判别,去除不可用伪卫星。
4)在去除所述不可用卫星后,接收机再获取伪卫星的载波相位、伪距。
5)进行整周模糊度的解算,根据接收机获取的伪卫星的载波相位、伪距,采用基于MC-LAMBDA算法的测向方法。
如图2所示,所述基于MC-LAMBDA算法的测向方法具体步骤如下:
步骤5.1、假设某载体上一短基线同时观测两颗伪卫星s1、s2,则双差载波相位和伪距观测方程为:
其中,为双差载波相位观测量;为双差伪距观测值;λ为载波波长;z为双差整周模糊度;载体上天线之间的距离很短,可假设同一历元载体上所有天线接收中心到同一颗伪卫星的视向量相同,记为s12,即为s12的转置;bE为基线向量;ε、e分别为双差载波相位和伪距观测噪声。
步骤5.2、载体上安装m+1个天线,共视伪卫星数为n+1,选择任一天线作为主天线,仰角最高的伪卫星作为参考星(以符号r表示),则基线i的观测方程表示为:
其中,S表示伪卫星视向量;其中,j=1,…,n;A为包含载波波长的设计矩阵;为某一伪卫星相对于仰角最高的伪卫星的双差伪距观测值;为某一伪卫星相对于仰角最高的伪卫星的双差载波相位观测值;v为量测噪声矩阵;D(yi)为yi的方差;为yi的协方差矩阵;zi为双差模糊度向量其中,i=1,…,m。
步骤5.3、将载体上m条基线观测量组合,得到多基线观测方程:
Y=AZ+SRBB+V
D(vec(Y))=QY,Z∈Zn×m,RTR=I
(4)
其中BB表示为BE在载体坐标系中相应的向量,RBB=BE,R为姿态旋转矩阵,满足RTR=I,I表示单位阵;QY为Y的协方差矩阵;Z为双差整周模糊度矩阵。
步骤5.4、引入向量算子vec,对式(4)进行向量化得到:
vec(Y)=vec(AZ)+vec(SRBB)+vec(V)
(5)根据矩阵乘积向量化公式,由式(5)可以得到:
步骤5.5、MC-LAMBDA方法解算整周模糊度及姿态矩阵即为解算下面最小化问题:
其中,为双差整周模糊度浮点解的协方差矩阵,vecZ为双差整周模糊度矩阵固定解;分别为双差整周模糊度浮点解和姿态旋转矩阵浮点解及其协方差矩阵,其中,为双差整周模糊度浮点解的协方差矩阵,为双差整周模糊度浮点解和姿态矩阵浮点解的协方差矩阵,为姿态矩阵浮点解和双差整周模糊度浮点解的协方差矩阵,为姿态矩阵浮点解的协方差矩阵;分别为姿态矩阵固定解及其协方差矩阵。
步骤5.6、采用LAMBDA方法搜索满足下式的所有整数值z:
从而解算整周模糊度。
步骤5.7、得到高精度多基线解向量,继而计算得到航向角和俯仰角。从而实现了基于伪卫星信号获取高精度测向结果的目的。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (4)
1.一种基于伪卫星信号的测向方法,其特征在于:包括如下步骤
步骤1:利用多颗伪卫星发送模拟信号;
步骤2:利用接收机接收伪卫星观测数据,包括伪距、载波相位;
步骤3:对伪卫星观测数据进行周跳检测、载波相位历元差计算以及根据对应的数据统计特征进行可用星判别,去除不可用伪卫星;
步骤4:所述接收机再获取伪卫星的载波相位、伪距;
步骤5:根据接收机获取的伪卫星载波相位、伪距,采用基于MC-LAMBDA算法的测向方法,得到多基线解向量。
2.根据权利要求1所述的一种基于伪卫星信号的测向方法,其特征在于:所述多颗伪卫星中的每一颗伪卫星根据卫星参数、伪卫星位置和参考时延,模拟卫星的卫星信号并发送伪卫星信号。
3.根据权利要求1所述的一种基于伪卫星信号的测向方法,其特征在于:所述伪卫星观测数据还包括载噪比测量值、多普勒频移测量值。
4.根据权利要求1所述的一种基于伪卫星信号的测向方法,其特征在于:所述步骤5具体包括如下步骤:
步骤51:设某载体上一短基线同时观测两颗伪卫星s1、s2,则双差载波相位和伪距观测方程为:
其中,为双差载波相位观测量;为双差伪距观测值;λ为载波波长;z为双差整周模糊度;载体上天线之间的距离很短,假设同一历元载体上所有天线接收中心到同一颗伪卫星的视向量相同,记为s12,即为s12的转置;bE为基线向量;ε、e分别为双差载波相位和伪距观测噪声;
步骤52:载体上安装m+1个天线,共视伪卫星数为n+1,选择任一天线作为主天线,仰角最高的伪卫星作为参考星(以符号r表示),则基线i的观测方程表示为:
其中,S表示伪卫星视向量;其中,j=1,…,n;A为包含载波波长的设计矩阵;为某一伪卫星相对于仰角最高的伪卫星的双差伪距观测值;为某一伪卫星相对于仰角最高的伪卫星的双差载波相位观测值;v为量测噪声矩阵;D(yi)为yi的方差;Qyi为yi的协方差矩阵;zi为双差模糊度向量其中,i=1,…,m;
步骤53:将载体上m条基线观测量组合,得到多基线观测方程:
其中BB表示为BE在载体坐标系中相应的向量,RBB=BE,R为姿态旋转矩阵,满足RTR=I,I表示单位阵;QY为Y的协方差矩阵;Z为双差整周模糊度矩阵;
步骤54:引入向量算子vec,对式(4)进行向量化得到:
vec(Y)=vec(AZ)+vec(SRBB)+vec(V) (5)
根据矩阵乘积向量化公式,由式(5)可以得到:
步骤55:MC-LAMBDA方法解算整周模糊度及姿态矩阵作为解算下面最小化问题:
其中,为双差整周模糊度浮点解的协方差矩阵,vecZ为双差整周模糊度矩阵固定解;分别为双差整周模糊度浮点解和姿态旋转矩阵浮点解及其协方差矩阵,其中,为双差整周模糊度浮点解的协方差矩阵,为双差整周模糊度浮点解和姿态矩阵浮点解的协方差矩阵,为姿态矩阵浮点解和双差整周模糊度浮点解的协方差矩阵,为姿态矩阵浮点解的协方差矩阵;分别为姿态矩阵固定解及其协方差矩阵;
步骤56:采用LAMBDA方法搜索满足下式的所有整数值z:
从而解算整周模糊度;
步骤57:得到多基线解向量,继而得到航向角和俯仰角。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211278195.7A CN115598686A (zh) | 2022-10-19 | 2022-10-19 | 一种基于伪卫星信号的测向方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211278195.7A CN115598686A (zh) | 2022-10-19 | 2022-10-19 | 一种基于伪卫星信号的测向方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115598686A true CN115598686A (zh) | 2023-01-13 |
Family
ID=84849688
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211278195.7A Pending CN115598686A (zh) | 2022-10-19 | 2022-10-19 | 一种基于伪卫星信号的测向方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115598686A (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101158719A (zh) * | 2007-11-16 | 2008-04-09 | 上海伽利略导航有限公司 | 伪卫星亚分米级室内定位系统及其载波相位定位方法 |
CN102565834A (zh) * | 2011-11-30 | 2012-07-11 | 重庆九洲星熠导航设备有限公司 | 一种单频gps测向系统及其测向定位方法 |
CN105785416A (zh) * | 2016-03-06 | 2016-07-20 | 中国人民解放军国防科学技术大学 | 基线约束下单频单历元gnss快速定向方法 |
CN110058281A (zh) * | 2019-04-29 | 2019-07-26 | 湖南国科微电子股份有限公司 | 动态定位方法及装置 |
CN113267796A (zh) * | 2021-05-13 | 2021-08-17 | 中国人民解放军92859部队 | 一种双天线gnss、rtk定位及测向方法 |
-
2022
- 2022-10-19 CN CN202211278195.7A patent/CN115598686A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101158719A (zh) * | 2007-11-16 | 2008-04-09 | 上海伽利略导航有限公司 | 伪卫星亚分米级室内定位系统及其载波相位定位方法 |
CN102565834A (zh) * | 2011-11-30 | 2012-07-11 | 重庆九洲星熠导航设备有限公司 | 一种单频gps测向系统及其测向定位方法 |
CN105785416A (zh) * | 2016-03-06 | 2016-07-20 | 中国人民解放军国防科学技术大学 | 基线约束下单频单历元gnss快速定向方法 |
CN110058281A (zh) * | 2019-04-29 | 2019-07-26 | 湖南国科微电子股份有限公司 | 动态定位方法及装置 |
CN113267796A (zh) * | 2021-05-13 | 2021-08-17 | 中国人民解放军92859部队 | 一种双天线gnss、rtk定位及测向方法 |
Non-Patent Citations (1)
Title |
---|
王冰 等: ""利用GPS解算载体测姿的算法研究"" * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7292185B2 (en) | Attitude determination exploiting geometry constraints | |
JP7153427B2 (ja) | 衛星測位システムを用いた測位方法および測位装置 | |
US10908293B2 (en) | System and method for detecting false global navigation satellite system satellite signals | |
US4990922A (en) | System and method for measuring ocean surface currents at locations remote from land masses using synthetic aperture radar | |
US5359332A (en) | Determination of phase ambiguities in satellite ranges | |
CN108802782B (zh) | 一种惯导辅助的北斗三频载波相位整周模糊度求解方法 | |
US20030231132A1 (en) | Method and apparatus for selecting optimal satellites in global positioning system | |
CN105334510A (zh) | 一种gnss-r地表探测装置和方法 | |
CN110320540A (zh) | 高精度集中式差分定位方法 | |
CN107703526B (zh) | 基线测向方法、装置和系统 | |
JP2010500578A (ja) | マルチプルgps測定タイプを加重最小二乗解へと融合するための方法 | |
Morales et al. | GNSS vertical dilution of precision reduction using terrestrial signals of opportunity | |
CN103364803A (zh) | 选星方法及应用该选星方法的卫星导航定位方法 | |
CN116879927B (zh) | 基于三天线共线共钟架构的船用卫星罗经艏向测定方法 | |
US10935669B2 (en) | Using SDP relaxation for optimization of the satellites set chosen for positioning | |
CN114895330A (zh) | 一种基于广播星历的单站位移监测方法、设备、存储介质 | |
CN113126022B (zh) | 一种双天线定位测向方法 | |
Lee et al. | Analysis of pseudolite augmentation for GPS airborne applications | |
US6211821B1 (en) | Apparatus and method for determining pitch and azimuth from satellite signals | |
CN115598686A (zh) | 一种基于伪卫星信号的测向方法 | |
RU2422844C2 (ru) | Способ определения угловой ориентации объекта | |
EP3841403A1 (en) | System and method for detecting false global navigation satellite system satellite signals | |
Chen | A study on real-time medium-range carrier phase-based GPS multiple-reference stations. | |
US20100034240A1 (en) | Transmitter-Receiver System | |
CN116413757B (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 |