CN105372692B - 一种北斗测姿接收机的快速整周模糊度方法 - Google Patents
一种北斗测姿接收机的快速整周模糊度方法 Download PDFInfo
- Publication number
- CN105372692B CN105372692B CN201410395770.0A CN201410395770A CN105372692B CN 105372692 B CN105372692 B CN 105372692B CN 201410395770 A CN201410395770 A CN 201410395770A CN 105372692 B CN105372692 B CN 105372692B
- Authority
- CN
- China
- Prior art keywords
- satellite
- carrier phase
- msub
- fuzziness
- pseudorange
- 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
Abstract
本发明属于无线通信技术领域,具体涉及一种北斗测姿接收机的快速整周模糊度方法。它包括:步骤一:观测数据的读取;包括卫星的三维坐标位置、卫星的伪距、载波相位、载波相位测量噪声值,步骤二:基准星选取,步骤三:求差,用其余卫星的载波相位与基准星的载波相位求差,步骤四:求模糊度备选集合,步骤五:模糊度确定,计算,并取最小值对应的模糊度初始值,作为本方法的输出值。本发明的效果是:本发明提供一种基于北斗导航系统的测姿接收机中使用的快速求解整周模糊度的方法。该方法不仅提高了计算效率和缩短了整周模糊度的时间,而且提高了整周模糊度的成功率。
Description
技术领域
本发明属于无线通信技术领域,具体涉及一种北斗测姿接收机的快速整周模糊度方法。
背景技术
随着卫星导航及相关技术的进步,导航卫星系统的应用日益广泛,其经济效益越来越明显,军事和民用上的作用也更为突出。导航卫星系统是以空间卫星为基础的无线电导航与定位系统,能为全球任意地点、任意多个用户(包括空中、陆上、海上以及外层空间的用户)全天候、实时、连续地提供精确的三维位置、三维速度和时间信息。目前,卫星定位系统主要包括美国的GPS(Global Positioning System),俄罗斯的GLONASS卫星定位系统,欧洲的伽利略(GALILEO)和我国的北斗(COMPASS)导航卫星系统。
随着北斗系统三步走的发展策略,目前北斗已经发布了ICD,并且自2012年末开始提供了区域性的导航服务。北斗(COMPASS)导航系统有B1、B2、B3三个频段,免费为用户提供可用性、可持续性和健全的导航服务,尤其在定位精度方面,可以达到厘米级。这些性能应用广泛,可以为科学家检测地壳的微小移动;测绘人员确定边界;为道路、桥梁建造测定距离等。
利用北斗(COMPASS)精确定位的同时,也能通过载波相位测量载体姿态。测量载体姿态的功能不仅可用在军事上,如进行武器(例如野战炮,陆用战车,战略导弹)的瞄准、炮兵阵地联测中的赋予基准射向、舰船载武器的初始对准等,也可以用在民用上,如飞机的航向和俯仰的,汽车的行驶方向等。
利用北斗(COMPASS)进行载体姿态测量是利用测得的载波相位信息,通过干涉测量技术实现,因此存在相位模糊问题,如何正确而又快速地解算相位模糊是北斗姿态测量的一项关键技术。因而利用载波相位进行实时动态测量的核心是运动中解算整周模糊度。
在动态实时条件下,需要快速可靠地在线模糊度求解。基于整周模糊度求解,己经提出了许多与此有关的算法,例如在观测域里搜索的模糊求解、在位置域里搜索的模糊度求解等,但是这些模糊度方法在求解速度上无法满足要求。具体的说,上述模糊度求解过程需要复杂的计算过程,对运算量的要求高,基本无法达到实时解算,实时使用的情况。
发明内容
本发明的目的是针对现有技术的缺陷,提供一种北斗测姿接收机的快速整周模糊度方法。
本发明是这样实现的:一种北斗测姿接收机的快速整周模糊度方法,包括下述步骤:
步骤一:观测数据的读取
需要读取的数据包括:卫星的三维坐标位置、卫星的伪距、载波相位、载波相位测量噪声值,
所述的卫星的三维坐标位置、卫星的伪距、载波相位、载波相位测量噪声值均由卫星导航解算给出,
步骤二:基准星选取
任意时刻,在天空中选取一颗卫星作为基准星,所述的基准星是步骤一中的卫星之一,
步骤三:求差
用其余卫星的载波相位与基准星的载波相位求差,求得的载波相位差用表示,其中i表示卫星序号,
用其余卫星的位置与基准星的位置求差,求得的位置差用hi表示,该参数表征了两颗卫星到观测点之间的距离差,
用其余卫星的伪距与基准星的伪距求差,求得的伪距差用Si表示,其中i表示卫星序号,
步骤四:求模糊度备选集合
用下述公式计算模糊度初始值Ni
其中,Si是步骤三球的伪距差,λ是载波波长,
本步骤得到集合N={Ni},其中i是卫星序号,
步骤五:模糊度确定
用下述公式计算
其中,hi是步骤三求得的位置差,λ是载波波长,Ni是步骤四求得的模糊度初始值,ni是步骤一中输入的载波相位测量噪声值,
计算所有卫星的σi之后,取σi最小值对应的模糊度初始值Ni作为本方法的输出值。
如上所述的一种北斗测姿接收机的快速整周模糊度方法,其中,步骤一至少读取5颗卫星的数据。
如上所述的一种北斗测姿接收机的快速整周模糊度方法,其中,步骤二中的基准星是相对地面观测点仰角最高的卫星。
本发明的效果是:本发明提供一种基于北斗导航系统的测姿接收机中使用的快速求解整周模糊度的方法。该方法不仅提高了计算效率和缩短了整周模糊度的时间,而且提高了整周模糊度的成功率。
具体实施方式
一种北斗测姿接收机的快速整周模糊度方法,包括下述步骤:
步骤一:观测数据的读取
需要读取的数据包括:卫星的三维坐标位置、卫星的伪距、载波相位、载波相位测量噪声值。
上述数据均是通过卫星导航解算得到的。
所述的载波相位测量噪声值用ni表示,i表示卫星序号,该载波相位测量噪声值对于不同卫星的数值不同,由卫星自身的导航解算给出。
本步骤至少读取5颗卫星的数据。
步骤二:基准星选取
任意时刻,在天空中选取一颗卫星作为基准星,该基准星相对地面观测点仰角最高。所述的基准星是步骤一中的卫星之一。
对于任意地面观测点,仰角最高的卫星是可以唯一确定的。
步骤三:求差
用其余卫星的载波相位与基准星的载波相位求差,求得的载波相位差用表示,其中i表示卫星序号。
用其余卫星的位置与基准星的位置求差,求得的位置差用hi表示,该参数表征了两颗卫星到观测点之间的距离差。
用其余卫星的伪距与基准星的伪距求差,求得的伪距差用Si表示,其中i表示卫星序号。
步骤四:求模糊度备选集合
用下述公式计算模糊度初始值Ni
其中,Si是步骤三球的伪距差,λ是载波波长。
本步骤得到集合N={Ni},其中i是卫星序号。
步骤五:模糊度确定
用下述公式计算
其中,hi是步骤三求得的位置差,λ是载波波长,Ni是步骤四求得的模糊度初始值,ni是步骤一中输入的载波相位测量噪声值。
计算所有卫星的σi之后,取σi最小值对应的模糊度初始值Ni作为本方法的输出值。
Claims (3)
1.一种北斗测姿接收机的快速整周模糊度方法,其特征在于,包括下述步骤:
步骤一:观测数据的读取
需要读取的数据包括:卫星的三维坐标位置、卫星的伪距、载波相位、载波相位测量噪声值,
所述的卫星的三维坐标位置、卫星的伪距、载波相位、载波相位测量噪声值均由卫星导航解算给出,
步骤二:基准星选取
任意时刻,在天空中选取一颗卫星作为基准星,所述的基准星是步骤一中的卫星之一,
步骤三:求差
用其余卫星的载波相位与基准星的载波相位求差,求得的载波相位差用表示,其中i表示卫星序号,
用其余卫星的位置与基准星的位置求差,求得的位置差用hi表示,
用其余卫星的伪距与基准星的伪距求差,求得的伪距差用Si表示,其中i表示卫星序号,
步骤四:求模糊度备选集合
用下述公式计算模糊度初始值Ni
<mrow>
<msub>
<mi>N</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>S</mi>
<mi>i</mi>
</msub>
<mi>&lambda;</mi>
</mfrac>
</mrow>
其中,Si是步骤三的伪距差,λ是载波波长,
本步骤得到集合N={Ni},其中i是卫星序号,
步骤五:模糊度确定
用下述公式计算
<mrow>
<msub>
<mi>&sigma;</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>h</mi>
<mi>i</mi>
</msub>
<mi>&lambda;</mi>
</mfrac>
<mo>-</mo>
<msub>
<mi>N</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>n</mi>
<mi>i</mi>
</msub>
</mrow>
其中,hi是步骤三求得的位置差,λ是载波波长,Ni是步骤四求得的模糊度初始值,ni是步骤一中输入的载波相位测量噪声值,
计算所有卫星的σi之后,取σi最小值对应的模糊度初始值Ni作为本方法的输出值。
2.如权利要求1所述的一种北斗测姿接收机的快速整周模糊度方法,其特征在于:步骤一至少读取5颗卫星的数据。
3.如权利要求1所述的一种北斗测姿接收机的快速整周模糊度方法,其特征在于:步骤二中的基准星是相对地面观测点仰角最高的卫星。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410395770.0A CN105372692B (zh) | 2014-08-13 | 2014-08-13 | 一种北斗测姿接收机的快速整周模糊度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410395770.0A CN105372692B (zh) | 2014-08-13 | 2014-08-13 | 一种北斗测姿接收机的快速整周模糊度方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105372692A CN105372692A (zh) | 2016-03-02 |
CN105372692B true CN105372692B (zh) | 2018-05-25 |
Family
ID=55375040
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410395770.0A Active CN105372692B (zh) | 2014-08-13 | 2014-08-13 | 一种北斗测姿接收机的快速整周模糊度方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105372692B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108463741A (zh) * | 2015-08-17 | 2018-08-28 | 加利福尼亚大学董事会 | 多历元gnss载波相位整数解析 |
CN106526629B (zh) * | 2016-12-23 | 2019-01-22 | 广州海格通信集团股份有限公司 | 卫星导航设备及其定向方法及装置 |
CN108490474A (zh) * | 2018-01-31 | 2018-09-04 | 中国航天电子技术研究院 | 基于阵列天线求解整周模糊度以实现单频测姿的方法 |
CN108919819B (zh) * | 2018-05-22 | 2023-10-27 | 北京和协导航科技有限公司 | 无人机导航通讯一体化系统及方法 |
CN111505694A (zh) * | 2020-05-07 | 2020-08-07 | 中航机载系统共性技术有限公司 | 一种面向机载的bds-3三天线阵多频点测姿方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101526608A (zh) * | 2009-04-14 | 2009-09-09 | 中国科学院国家天文台 | 卫星导航定位中三频码波伪距和载波相位的组合方法 |
CN101770033A (zh) * | 2010-02-08 | 2010-07-07 | 东南大学 | 连续运行参考站系统站间整周模糊度网络固定方法 |
CN101825717A (zh) * | 2010-04-16 | 2010-09-08 | 北京航空航天大学 | 一种基于载波平滑码伪距技术的动态定姿方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7102563B2 (en) * | 2004-02-26 | 2006-09-05 | Topcon Gps Llc | Methods and apparatuses of estimating the position of a mobile user in a system of satellite differential navigation |
JP5760001B2 (ja) * | 2009-11-17 | 2015-08-05 | トプコン ポジショニング システムズ, インク. | 全地球的航法衛星システム受信機における異常測定値の検出および補正ならびにアンビギュイティの決定 |
NL2009695C2 (en) * | 2012-10-25 | 2014-05-06 | Fugro N V | Ppp-rtk method and system for gnss signal based position determination. |
-
2014
- 2014-08-13 CN CN201410395770.0A patent/CN105372692B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101526608A (zh) * | 2009-04-14 | 2009-09-09 | 中国科学院国家天文台 | 卫星导航定位中三频码波伪距和载波相位的组合方法 |
CN101770033A (zh) * | 2010-02-08 | 2010-07-07 | 东南大学 | 连续运行参考站系统站间整周模糊度网络固定方法 |
CN101825717A (zh) * | 2010-04-16 | 2010-09-08 | 北京航空航天大学 | 一种基于载波平滑码伪距技术的动态定姿方法 |
Non-Patent Citations (2)
Title |
---|
GNSS整周模糊度实时解算方法研究;李娜;《中国优秀硕士学位论文全文数据库 基础科学辑》;20140715(第07期);全文 * |
一种实时 GPS 姿态测量中的整周模糊度的解算方法;彭晓刚等;《测绘通报》;20111231(第5期);第16-18页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105372692A (zh) | 2016-03-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106990424B (zh) | 一种双天线gps测姿方法 | |
Noureldin et al. | Fundamentals of inertial navigation, satellite-based positioning and their integration | |
CN107003386B (zh) | 一种卫星导航测姿方法和装置及无人机 | |
CN107272039B (zh) | 一种基于双天线gps的定位测姿方法 | |
CN105372692B (zh) | 一种北斗测姿接收机的快速整周模糊度方法 | |
CN108871336B (zh) | 一种车辆位置估算系统及方法 | |
CN106646564B (zh) | 一种基于低轨卫星增强导航方法 | |
US6281841B1 (en) | Direction determining apparatus | |
CN108226980A (zh) | 基于惯性测量单元的差分gnss与ins自适应紧耦合导航方法 | |
CN104807479A (zh) | 一种基于主惯导姿态变化量辅助的惯导对准性能评估方法 | |
CN103760585B (zh) | 一种适用林区的星‑地结合定位方法 | |
CN103033822B (zh) | 移动信息确定装置、方法以及接收机 | |
Eliasson | A Kalman filter approach to reduce position error for pedestrian applications in areas of bad GPS reception | |
Bitner et al. | Multipath and spoofing detection using angle of arrival in a multi-antenna system | |
Angrisano et al. | Using local redundancy to improve GNSS absolute positioning in harsh scenario | |
CN108151765A (zh) | 一种在线实时估计补偿磁强计误差的定位测姿方法 | |
Deng et al. | Underwater map-matching aided inertial navigation system based on multi-geophysical information | |
CN108106597A (zh) | 全捷联激光导引头在目标出线性视场情况下角度测量方法 | |
Groves et al. | Enhancing micro air vehicle navigation in dense urban areas using 3D mapping aided GNSS | |
Huang et al. | An online calibration method for the DVL scale factor error based on a circular track | |
CN106123894B (zh) | 基于干涉条纹匹配的InSAR/INS组合导航方法 | |
RU2414685C1 (ru) | Способ определения угловой ориентации объекта | |
JP2017009561A (ja) | 測量装置、測量方法、及び測量プログラム | |
Gordeev et al. | The influence of a geometric factor on a satellite positioning accuracy | |
Projović | Application of relative positioning in topographic survey preparations on a full basis in artillery |
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 |