JP2000275320A - Gps receiver and gps reception method applied to the same - Google Patents

Gps receiver and gps reception method applied to the same

Info

Publication number
JP2000275320A
JP2000275320A JP7848599A JP7848599A JP2000275320A JP 2000275320 A JP2000275320 A JP 2000275320A JP 7848599 A JP7848599 A JP 7848599A JP 7848599 A JP7848599 A JP 7848599A JP 2000275320 A JP2000275320 A JP 2000275320A
Authority
JP
Japan
Prior art keywords
gps
data
error
visible
value
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
Application number
JP7848599A
Other languages
Japanese (ja)
Other versions
JP3544887B2 (en
Inventor
Yuichi Sato
友一 佐藤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Toshiba Corp
Original Assignee
Toshiba Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Toshiba Corp filed Critical Toshiba Corp
Priority to JP7848599A priority Critical patent/JP3544887B2/en
Publication of JP2000275320A publication Critical patent/JP2000275320A/en
Application granted granted Critical
Publication of JP3544887B2 publication Critical patent/JP3544887B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Position Fixing By Use Of Radio Waves (AREA)
  • Navigation (AREA)

Abstract

PROBLEM TO BE SOLVED: To improve the reliability of the navigating/positioning accuracy and navigation/positioning of a GPS receiver. SOLUTION: When the position of a user having a GPS(global positioning system) receiver and the initial estimated value of a clock error are set (step S11), the GPS receiver sets the parameter value indicating the state of its own device (step S12), acquires the observation data of all visible GPS satellites (step S13), and evaluates a Lagrangian factor term (step S14). Then the receiver successively solves navigation equations based on the evaluated results (step S15) and sets the status parameter of its own device by performing the processing in the step S12 based on the calculated results of the equations. Namely, the receiver sets a threshold for evaluating the quality of observation data from the visible GPS satellites.

Description

【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION

【0001】[0001]

【発明の属する技術分野】この発明は、GPS(Global
Positioning System)衛星を利用して、自装置の位置を
計測するGPS受信装置及びこの装置に適用されるGP
S受信方法において、特に、航法もしくは測位精度、及
び航法もしくは測位における信頼性の向上を図るように
したものである。
TECHNICAL FIELD The present invention relates to a GPS (Global
Positioning System) A GPS receiver that measures the position of its own device using a satellite, and a GP applied to this device.
In the S receiving method, particularly, navigation or positioning accuracy and reliability in navigation or positioning are improved.

【0002】[0002]

【従来の技術】周知のように、例えば自動車等の移動体
に搭載されるGPS受信装置は、例えばGPS衛星から
の衛星信号を受信して、この受信信号から移動体の移動
方向及び移動位置を測定している。この種のGPS受信
装置は、4個の可視GPS衛星からの衛星信号が受信可
能である。すなわち、このGPS受信装置では、4個以
上の可視GPS衛星からの衛星信号を受信して連続的に
測位を行なっている時において、測位誤差が最小となる
GPS衛星を予測し、移動に伴って、次々と受信すべき
GPS衛星を切り替えて、常に測位精度を最適にしてお
くようにしている。
2. Description of the Related Art As is well known, a GPS receiver mounted on a mobile body such as an automobile receives a satellite signal from a GPS satellite, for example, and determines a moving direction and a moving position of the mobile body from the received signal. Measuring. This type of GPS receiver can receive satellite signals from four visible GPS satellites. That is, in this GPS receiving apparatus, when satellite signals from four or more visible GPS satellites are received and positioning is continuously performed, a GPS satellite with a minimum positioning error is predicted, and the GPS satellite is moved with movement. The GPS satellites to be received are switched one after another so that the positioning accuracy is always optimized.

【0003】しかし、上記GPS受信装置では、4個の
GPS衛星からの衛星信号により測位を行なっている
が、1つの衛星信号に大きな誤差があると3つの衛星信
号により測位するため測位精度に大きく影響することに
なり、航法もしくは測位における信頼性を低下させてし
まうことになる。また、GPS衛星の切り替わりによ
り、航法もしくは測位データが不連続でもある。なお、
GPS衛星からの衛星信号には、予め誤差が与えられて
おり、また、差分航法を行なうシステムにおいても、例
えば地上局からの軌道データ等に誤差が与えられてい
る。
However, in the above-mentioned GPS receiving apparatus, positioning is performed using satellite signals from four GPS satellites. However, if there is a large error in one satellite signal, positioning is performed using three satellite signals, so that positioning accuracy is large. This has the effect of reducing the reliability of navigation or positioning. In addition, navigation or positioning data is discontinuous due to switching of GPS satellites. In addition,
An error is given in advance to a satellite signal from a GPS satellite, and also in a system that performs differential navigation, an error is given to, for example, orbit data from a ground station.

【0004】従って、近時においては、GPS衛星の切
り替わりに起因する航法もしくは測位データの不連続性
を抑え、精度の向上を図り、さらに差分航法への拡張性
を確保できるようなGPS受信装置の改良が強く望まれ
ている。
[0004] Therefore, recently, a GPS receiving apparatus which can suppress discontinuity of navigation or positioning data due to switching of GPS satellites, improve accuracy, and further ensure scalability to differential navigation. Improvement is strongly desired.

【0005】[0005]

【発明が解決しようとする課題】以上のように、従来の
GPS受信装置では、4個のGPS衛星からの衛星信号
により測位を行なっているが、1つの衛星信号に大きな
誤差があると3つの衛星信号により測位するため測位精
度に大きく影響することになり、航法もしくは測位にお
ける信頼性を低下させるという問題を有している。ま
た、GPS衛星の切り替わりにより、航法もしくは測位
データが不連続になるという問題も有している。
As described above, in the conventional GPS receiver, positioning is performed using satellite signals from four GPS satellites. However, if one satellite signal has a large error, three satellite signals are used. Since positioning is performed by satellite signals, positioning accuracy is greatly affected, and there is a problem that reliability in navigation or positioning is reduced. There is also a problem that the navigation or the positioning data becomes discontinuous due to the switching of the GPS satellites.

【0006】この発明の目的は、選択衛星の切り替わり
に起因する航法/測位のデータの不連続性を最小限に抑
制でき、かつ航法/測位精度及び航法/測位における信
頼性の向上を図り得るGPS受信装置及びこの装置に適
用されるGPS受信方法を提供することにある。
An object of the present invention is to provide a GPS capable of minimizing the discontinuity of navigation / positioning data due to switching of a selected satellite and improving the navigation / positioning accuracy and the reliability in navigation / positioning. An object of the present invention is to provide a receiving device and a GPS receiving method applied to the receiving device.

【0007】[0007]

【課題を解決するための手段】この発明は、GPS衛星
を利用して、自装置の位置を計測するGPS受信装置に
おいて、可視GPS衛星全てから衛星信号を受信する受
信手段と、この受信手段で得られた可視GPS衛星全て
の衛星信号から可視GPS衛星に関する擬似距離データ
もしくは搬送波位相データと、可視GPS衛星の放送暦
データとを含む観測データを抽出する観測データ抽出手
段と、この観測データ抽出手段で得られた可視GPS衛
星全ての観測データを入力して、衛星信号に与えられる
誤差値もしくは自装置に発生する誤差値の推定確率に最
大を与える誤差の推定値を算出する誤差推定値算出手段
と、この誤差推定値算出手段により算出された誤差の推
定値に基づいて、受信手段で得られた可視GPS衛星全
ての観測データの中から必要なGPS衛星からの観測デ
ータを決定する観測データ決定手段とを備えるようにし
たものである。
SUMMARY OF THE INVENTION According to the present invention, there is provided a GPS receiving apparatus for measuring the position of an own apparatus using a GPS satellite, receiving means for receiving satellite signals from all visible GPS satellites, and the receiving means. Observation data extraction means for extracting observation data including pseudorange data or carrier phase data relating to the visible GPS satellites and broadcast calendar data of the visible GPS satellites from the obtained satellite signals of all the visible GPS satellites, and this observation data extraction means Error observation value calculation means for inputting the observation data of all the visible GPS satellites obtained in the above and calculating the error value given to the satellite signal or the error estimation value which gives the maximum to the estimation probability of the error value occurring in the own device. And, based on the error estimation value calculated by the error estimation value calculation means, the observation data of all the visible GPS satellites obtained by the reception means. Observation data from the necessary GPS satellites is obtained as provided with the observed data determining means for determining.

【0008】この構成によれば、可視GPS衛星全てか
らの観測データが取り込まれ、これら観測データを用い
て、衛星信号に与えられる誤差値もしくは自装置に発生
する誤差値の推定確率に最大を与える誤差の推定値が算
出され、この誤差推定値に基づいて、可視GPS衛星全
ての観測データの中から必要なGPS衛星からの観測デ
ータが決定される。すなわち、観測データを入力して、
推定確率が最大となるような誤差の推定値を算出し、こ
の推定値に基づいて推定される観測データと可視GPS
衛星全ての観測データとを比較照合することにより、観
測データの品質を評価でき、この結果を基に所定値以下
のずれの観測データを信頼性の高い観測データとして選
択でき、所定値より大きいずれの観測データに対しては
信頼性の低い観測データと判断し棄却できるようにな
る。
According to this configuration, observation data from all the visible GPS satellites is fetched, and the maximum probability is given to the error value given to the satellite signal or the estimation probability of the error value occurring in the own device using these observation data. An estimated value of the error is calculated, and the required observation data from the GPS satellites is determined from the observation data of all the visible GPS satellites based on the estimated error value. That is, input observation data,
An estimated value of the error that maximizes the estimated probability is calculated, and the observed data and the visible GPS estimated based on the estimated value are calculated.
By comparing and checking the observation data of all satellites, the quality of the observation data can be evaluated, and based on this result, observation data with a deviation of less than a predetermined value can be selected as highly reliable observation data, and deviations larger than the predetermined value can be selected. Observation data can be discarded because it is judged to be unreliable observation data.

【0009】このため、通常の航法においてみられる、
選択衛星の切り替わりに起因する航法/測位のデータの
不連続性が最小限に抑えられ、使用可能な観測データを
最適利用しているので、航法/測位精度及び航法/測位
における信頼性の向上を図ることが可能となる。さらに
GPS衛星等の故障による一部の観測データの品質劣化
に対する耐性が得られる。
For this reason, it can be seen in normal navigation,
Since the discontinuity of navigation / positioning data due to switching of the selected satellite is minimized and the available observation data is optimally used, the navigation / positioning accuracy and the reliability in navigation / positioning are improved. It becomes possible to plan. Furthermore, resistance to quality deterioration of some observation data due to a failure of a GPS satellite or the like can be obtained.

【0010】また、上記構成において、受信手段は、予
め決められた位置に設置された地上局から送出され、地
上局による擬似距離データ及び観測時刻データを含む地
上局信号を受信する手段を有し、観測データ抽出手段
は、受信手段で受信された地上局信号から擬似距離デー
タ及び観測時刻データを抽出する手段を有し、誤差推定
値算出手段は、地上局による擬似距離データ及び観測時
刻データを入力して、地上局による観測衛星位置の誤差
値の推定確率に最大を与える推定値を算出する手段を有
することを特徴とする。このようにすることで、差分航
法を行なう場合に、高精度で、かつ信頼性の高い航法を
実現できる。
In the above configuration, the receiving means has means for receiving a ground station signal transmitted from a ground station installed at a predetermined position and containing pseudorange data and observation time data by the ground station. The observation data extraction means has means for extracting pseudorange data and observation time data from the ground station signal received by the reception means, and the error estimation value calculation means calculates pseudorange data and observation time data by the ground station. The apparatus further comprises means for calculating an estimated value for inputting the maximum value to the estimation probability of the error value of the observation satellite position by the ground station. By doing so, highly accurate and highly reliable navigation can be realized when performing differential navigation.

【0011】[0011]

【発明の実施の形態】以下、この発明の実施形態につい
て図面を参照して詳細に説明する。図1は、この発明に
係るGPS受信装置の一実施形態を示す回路ブロック図
である。図2は、このGPS受信装置と複数のGPS衛
星との位置関係及び測位概念を示す図である。
Embodiments of the present invention will be described below in detail with reference to the drawings. FIG. 1 is a circuit block diagram showing an embodiment of a GPS receiver according to the present invention. FIG. 2 is a diagram showing a positional relationship between the GPS receiver and a plurality of GPS satellites and a positioning concept.

【0012】図1において、例えばGPS衛星から送信
されたGPS衛星信号は、アンテナ11にて受信され、
その直後に低雑音増幅器12にて低雑音増幅されること
によりRF(Radio Frequency) 信号に変換され、乗算器
13の一方の入力端に供給される。この乗算器13の他
方の入力端には、局部発振器14から発生されたGPS
中心周波数を有する局発信号をてい倍器15でてい倍さ
れた第1のローカル信号が供給される。そして、乗算器
13は、RF信号に第1のローカル信号を乗算すること
により、IF(中間周波数)帯まで周波数変換する。
In FIG. 1, for example, a GPS satellite signal transmitted from a GPS satellite is received by an antenna 11,
Immediately thereafter, the signal is low-noise amplified by the low-noise amplifier 12, converted into an RF (Radio Frequency) signal, and supplied to one input terminal of the multiplier 13. The other input terminal of the multiplier 13 has a GPS signal generated from the local oscillator 14 connected thereto.
A first local signal is provided which is obtained by multiplying a local signal having a center frequency by a multiplier 15. Then, the multiplier 13 performs frequency conversion to an IF (intermediate frequency) band by multiplying the RF signal by the first local signal.

【0013】このIF信号は、その時受信することので
きる可視GPS衛星信号で、各GPS衛星を識別するた
めのPN(擬似雑音)コードに応じて予め決められた複
数の振幅が拡散されている。このIF信号は、ベースバ
ンド部101〜10Nにそれぞれ供給される。
The IF signal is a visible GPS satellite signal that can be received at that time, and a plurality of predetermined amplitudes are spread in accordance with a PN (pseudo noise) code for identifying each GPS satellite. This IF signal is supplied to the baseband units 101 to 10N, respectively.

【0014】このうちのベースバンド処理部101にお
いて、入力されたIF信号は、PNコード解除器16の
一方の入力端に供給される。PNコード解除器16の他
方の入力端には、PNコード発生器17により発生され
る自装置に割り当てられたPNコード発振信号が供給さ
れる。そして、PNコード解除器16は、IF信号にP
Nコード発振信号を乗算することにより、逆拡散して、
自装置に割り当てられた各GPS衛星信号を抽出するこ
とができる。よって、ベースバンド処理部101〜10
Nの数がそのGPS受信装置のチャネル数となり最大受
信可能な衛星数である。
In the baseband processing section 101, the input IF signal is supplied to one input terminal of the PN code canceler 16. The other input terminal of the PN code canceler 16 is supplied with a PN code oscillation signal generated by the PN code generator 17 and assigned to the own device. Then, the PN code release unit 16 outputs the P signal to the IF signal.
Despread by multiplying the N code oscillation signal,
Each GPS satellite signal assigned to the own device can be extracted. Therefore, the baseband processing units 101 to 10
The number N is the number of channels of the GPS receiver and is the maximum number of receivable satellites.

【0015】また、PNコード解除器16により抽出さ
れた衛星識別信号は、位相比較器18に供給され、NC
O(Newmedical Control Ocial)19により発生されるG
PS中心周波数を有する第2のローカル信号と位相比較
されることにより、その位相差成分が受信信号として抽
出される。この受信信号は、受信制御部200に供給さ
れる。
The satellite identification signal extracted by the PN code canceler 16 is supplied to a phase comparator 18 where it is transmitted to an NC.
G generated by O (Newmedical Control Ocial) 19
By comparing the phase with the second local signal having the PS center frequency, the phase difference component is extracted as a received signal. This reception signal is supplied to the reception control unit 200.

【0016】受信制御部200は、入力された受信信号
に基づいて、PNコード発生器17に設定するPNコー
ド位相設定値、NCO19に設定する搬送波位相設定値
をリアルタイムで演算し、この演算結果であるPNコー
ド位相設定信号及び搬送波位相設定信号をPNコード発
生器17及びNCO19にそれぞれ与えて、PNコード
解除器16及び位相比較器18にそれぞれ入力される信
号に対する位相追尾を行なうように制御する。また、受
信制御部200は、入力された受信信号に基づいて演算
された擬似距離データもしくは搬送波位相データと、受
信信号のデータ復調により得られる可視GPS衛星の放
送暦データとを含む観測データを測位演算部300に出
力する。この処理は、各ベースバンド部101〜10N
で各GPS衛星信号毎に独立して行なわれる。
The reception control unit 200 calculates a PN code phase set value to be set in the PN code generator 17 and a carrier phase set value to be set in the NCO 19 in real time based on the input received signal. A certain PN code phase setting signal and a carrier phase setting signal are supplied to a PN code generator 17 and an NCO 19, respectively, so as to perform a phase tracking for signals input to a PN code canceler 16 and a phase comparator 18, respectively. In addition, the reception control unit 200 performs positioning of observation data including pseudorange data or carrier phase data calculated based on the input received signal and broadcast GPS data of a visible GPS satellite obtained by demodulating the received signal. Output to the arithmetic unit 300. This processing is performed for each of the baseband units 101 to 10N.
Is performed independently for each GPS satellite signal.

【0017】測位演算部300は、各ベースバンド処理
部101〜10Nに対する観測データに基づいて測位演
算処理を実行することにより、自装置の位置を示す測位
データを求める。ここで得られる測位データは、図2に
示す如く例えば4個のGPS衛星401〜404から見
たGPS受信装置の位置を表す(X,Y,Z)の3次元
情報と、その変化分(速度)を表す(dX/dt,dY
/dt,dZ/dt)である。
The positioning operation section 300 executes positioning operation processing based on observation data for each of the baseband processing sections 101 to 10N, thereby obtaining positioning data indicating the position of the own apparatus. The positioning data obtained here is, as shown in FIG. 2, three-dimensional information (X, Y, Z) indicating the position of the GPS receiver viewed from, for example, four GPS satellites 401 to 404, and the change (speed). (DX / dt, dY)
/ Dt, dZ / dt).

【0018】また、図2において、GPS受信装置は、
GPS衛星401を利用して自装置におけるX成分の位
置情報が得られ、GPS衛星402を利用して自装置に
おけるY成分の位置情報が得られ、GPS衛星403を
利用して自装置におけるZ成分の位置情報が得られる。
さらに、GPS衛星404を利用して自装置における動
作時刻情報が得られる。
In FIG. 2, the GPS receiver is
Position information of the X component in the own device is obtained using the GPS satellite 401, position information of the Y component in the own device is obtained using the GPS satellite 402, and Z component in the own device is used using the GPS satellite 403. Is obtained.
Further, the operation time information of the own device can be obtained using the GPS satellite 404.

【0019】以上が一般的なGPS受信装置の構成であ
る。
The above is the configuration of a general GPS receiver.

【0020】次に、上記GPS受信装置における可視G
PS衛星全体からの観測データの処理概念について説明
する。まず、従来のGPS受信装置では、3次元位置情
報(X,Y,Z)と、時間情報(T)とを用いて自装置
の位置を計測していることにより、どれか1つの情報が
間違っている可能性も高く、またデータが不足して大き
い誤差が生じることになる。
Next, the visible G in the GPS receiver is
The concept of processing observation data from the entire PS satellite will be described. First, in the conventional GPS receiving apparatus, one of the information is incorrect because the position of the own apparatus is measured using the three-dimensional position information (X, Y, Z) and the time information (T). It is highly possible that the data is insufficient and a large error occurs due to lack of data.

【0021】そこで、この発明は、可視GPS衛星全て
から衛星信号を取り込み、多数決原理を用いて品質の良
い観測データを決定するものである。すなわち、受信す
べき可視GPS衛星全てからの衛星信号には、各GPS
衛星で予め与えられる誤差やGPS受信装置内で発生す
る誤差が含まれている。
Therefore, the present invention fetches satellite signals from all visible GPS satellites and determines high quality observation data using the majority rule. That is, satellite signals from all visible GPS satellites to be received include each GPS
The error includes an error given in advance by a satellite and an error generated in the GPS receiver.

【0022】測位演算部300は、受信制御部200か
ら与えられる可視GPS衛星全ての観測データを入力し
て、衛星信号に与えられる誤差値もしくは自装置に発生
する誤差値の推定確率に最大を与える誤差の推定値を算
出する。すなわち、誤差推定値を求める式として、以下
の式(1)に示すようなTOTAL Lagrangian(トータルラ
グランジュ)を導入し、 TOTAL Lagrangian:L(Δz”G,Δz’)=L1+L2+L3+L4 …(1) 以下の式(2)に示すこのラグランジュに対する変分方
程式を解くことにより、推定確率に最大を与えるユーザ
位置/クロック誤差の推定値を算出する。
The positioning operation unit 300 receives the observation data of all the visible GPS satellites supplied from the reception control unit 200 and gives the maximum to the error value given to the satellite signal or the estimation probability of the error value generated in the own device. Calculate the estimated value of the error. That is, as an expression for obtaining the error estimates, introduced TOTAL Lagrangian (total Lagrangian) as shown in the following equation (1), TOTAL Lagrangian: L (Δz "G, Δz ') = L1 + L2 + L3 + L4 ... (1) below By solving the variational equation for this Lagrange shown in Expression (2), an estimated value of the user position / clock error that gives the maximum estimated probability is calculated.

【0023】 航法方程式:∂L/∂Δz”G =0 ∂L/∂Δz’=0…(2) なお、Δz”Gは、可視GPS衛星♯Gの位置/クロッ
ク誤差に関する放送暦データからの偏差にSAGNAC
効果、つまり地球の自転による補正係数を加えた値を示
す。また、Δz’は、ユーザの位置の更新前位置からの
偏差及び自装置のクロック誤差を表す4次元変数で、こ
れが推定対象である。
Navigation equation: ∂L / ∂Δz ″ G = 0 ∂L / ∂Δz ′ = 0 (2) where Δz ″ G is a value obtained from broadcast calendar data on the position / clock error of the visible GPS satellite ♯G. SAGNAC for deviation
It shows the effect, that is, the value to which the correction coefficient by the rotation of the earth is added. Δz ′ is a four-dimensional variable representing the deviation of the user's position from the pre-update position and the clock error of the own device, which is the estimation target.

【0024】次に、トータルラグランジュに使用される
上記L1,L2,L3,L4それぞれの項の意味及び作
用について説明する。まず、L1はGPS受信装置位置
/クロック偏差のアプリオリ推定のラグランジュの式で
以下の式(3)となる。
Next, the meaning and operation of each of the above L1, L2, L3 and L4 terms used in the total Lagrange will be described. First, L1 is a Lagrangian equation of a priori estimation of the GPS receiver position / clock deviation, and is given by the following equation (3).

【数1】 (Equation 1)

【0025】ここで、Sは観測データ取得/更新前の共
分散行列である。つまり、L1は観測が行なわれる直前
の時点でのGPS受信装置の状態に関する推定を表して
いる。
Here, S is a covariance matrix before acquisition / update of observation data. That is, L1 represents the estimation of the state of the GPS receiver immediately before the observation is performed.

【0026】L2は、GPS受信装置の位置推定のラグ
ランジュの式で以下の式(4)となる。
L2 is the Lagrange's equation for estimating the position of the GPS receiver, and is expressed by the following equation (4).

【数2】 (Equation 2)

【0027】ここで、σGRはGPS受信装置によるGP
S衛星♯Gの受信装置誤差レベルである。εGはデータ
採用のフラグで、GPS衛星♯Gのデータを採用するか
否かに従って1または0をとる。VG は、時間成分1
と、空間成分として推定更新前のGPS受信装置の位置
から見たGPS衛星♯Gの方向の単位ベクトルからなる
4次元ベクトルである。ΔyG は、観測データ(擬似距
離)に電離層伝搬遅延補正と対流圏伝搬遅延補正を考慮
して計算される誘導観測量である。ΣGは、可視GPS
衛星全ての和を表している。すなわち、L2は、可視G
PS衛星に対する観測データに起因する状態推定を表し
ていることになる。
Here, σ GR is the GP obtained by the GPS receiver.
This is the receiver error level of the S satellite $ G. ε G is a data adoption flag, which takes 1 or 0 according to whether or not to adopt data of the GPS satellite ΔG. V G is the time component 1
And a four-dimensional vector composed of unit vectors in the direction of the GPS satellite ♯G as viewed from the position of the GPS receiver before the estimation update as a spatial component. Δy G is a guided observation amount calculated in consideration of ionospheric propagation delay correction and tropospheric propagation delay correction to observation data (pseudorange). Σ G is visible GPS
Represents the sum of all satellites. That is, L2 is the visible G
This indicates the state estimation resulting from the observation data for the PS satellite.

【0028】次に、L3は、衛星位置/クロック誤差に
関するアプリオリ推定のラグランジュの式で以下の式
(5)となる。
Next, L3 is the Lagrange's equation for a priori estimation relating to the satellite position / clock error, and is given by the following equation (5).

【数3】 (Equation 3)

【0029】ここで、σGEはGPS衛星♯Gのエフェメ
リス精度を表している。すなわち、GPS衛星の放送暦
が提供する衛星位置情報は、完全に正確なわけではな
く、これには精度情報も付加されている。これより、衛
星放送暦から計算した衛星位置/クロック誤差について
も、精度情報に従った重みをかけてこれを推定変数とし
ている。
Here, σ GE represents the ephemeris accuracy of the GPS satellite ΔG. That is, the satellite position information provided by the GPS satellite broadcast calendar is not completely accurate, and accuracy information is also added thereto. Thus, the satellite position / clock error calculated from the satellite broadcast calendar is also weighted according to the accuracy information and used as an estimated variable.

【0030】次に、L4は、予め決められた位置に設置
された地上局による観測の衛星位置推定のラグランジュ
の式で以下の式(6)となる。
Next, L4 is a Lagrange's equation for estimating the satellite position observed by the ground station installed at a predetermined position, and is given by the following equation (6).

【数4】 (Equation 4)

【0031】ここで、σGBは、例えば地上局BによるG
PS衛星♯Gの受信機誤差レベルを表している。すなわ
ち、L4は差分航法を行なう場合に導入される。差分航
法を行なわない場合には、L4は0となる。なお、ΣGB
は、地上局Bにおける可視GPS衛星全ての和である。
Here, σ GB is, for example, G
It represents the receiver error level of PS satellite $ G. That is, L4 is introduced when performing differential navigation. When the differential navigation is not performed, L4 becomes 0. Σ GB
Is the sum of all visible GPS satellites at ground station B.

【0032】次に、上記ラグランジュの式を利用した測
位演算部300の処理について図3〜図5のフローチャ
ートを参照して説明する。図3において、測位演算部3
00は、まず、GPS受信装置を所有したユーザの位置
及びクロック誤差の初期推定値が設定されると(ステッ
プS11)、自装置の状態を示すパラメータ値を設定し
(ステップS12)、可視GPS衛星全ての観測データ
を取得し(ステップS13)、ラグランジュ要素項の評
価を行なう(ステップS14)。すなわち、上記L1、
L2、L3、L4の項内のパラメータ値の評価を行な
う。そして、この評価結果に基づいて、航法方程式の逐
次解法を行ない(ステップS15)、この算出結果に基
づいて、ステップS12の処理で自装置の状態パラメー
タ設定を行なう。つまり、可視GPS衛星からの観測デ
ータの品質評価を行なうためのしきい値を設定する。
Next, the processing of the positioning operation unit 300 using the Lagrangian equation will be described with reference to the flowcharts of FIGS. In FIG. 3, the positioning calculation unit 3
First, when an initial estimated value of the position and clock error of the user who owns the GPS receiving device is set (step S11), a parameter value indicating the state of the own device is set (step S12), and the visible GPS satellite is set. All observation data is acquired (step S13), and the Lagrange element term is evaluated (step S14). That is, the above L1,
The parameter values in the terms L2, L3, L4 are evaluated. Then, based on the evaluation result, a sequential solution of the navigation equation is performed (step S15), and based on the calculation result, the state parameter of the own device is set in the process of step S12. That is, a threshold for evaluating the quality of observation data from the visible GPS satellites is set.

【0033】図4は、上記ステップS14の処理の詳細
を示している。すなわち、測位演算部300は、上記ス
テップS13で可視GPS衛星全ての観測データを入力
した時点で、ラグランジュの要素項L1、L2、L3、
L4が所定値以上であるか否かを判定し(ステップS1
41)、所定値を満たしていない場合に(NG)、2回
連続NGであるか否かを判定し(ステップS142)、
2回連続NGである場合に(YES)、可視GPS衛星
数が4より大きいか否かを判定する。ここで、4より大
きい数である場合に(YES)、ラグランジュ要素項L
2、L4のパラメータ値εGを0に設定して(ステップ
S144)、観測データの棄却を図るようにする。ま
た、上記ステップS141で評価値が所定値以上である
場合(GOOD)や、ステップS142でNGが1回だ
けである場合(NO)、ステップS143で可視GPS
衛星数が4以下である場合(NO)には、εGを1に設
定する(ステップS145)。
FIG. 4 shows the details of the processing in step S14. That is, the positioning calculation unit 300, at the time when the observation data of all the visible GPS satellites is input in step S13, the Lagrange element terms L1, L2, L3,
It is determined whether or not L4 is equal to or greater than a predetermined value (step S1).
41), when the predetermined value is not satisfied (NG), it is determined whether or not two consecutive NGs (step S142).
In the case of two consecutive NGs (YES), it is determined whether or not the number of visible GPS satellites is greater than four. If the number is greater than 4 (YES), the Lagrange element term L
2. The parameter value ε G of L4 is set to 0 (step S144), and the observation data is rejected. If the evaluation value is equal to or more than the predetermined value in step S141 (GOOD) or if the NG is only once in step S142 (NO), the visible GPS is determined in step S143.
If the number of satellites is 4 or less (NO), ε G is set to 1 (step S145).

【0034】次に、図5は、上記ステップS15の処理
の詳細を示している。すなわち、測位演算部300は、
上記式(2)示すような航法方程式を解いて、推定確率
が最大となる誤差推定値Δz’,Δz”Gを求める(ス
テップS151)。そして、これら誤差推定値Δz’,
Δz”Gに基づいて推定される観測データに入力された
観測データが収束するか否かを判定し(ステップS15
2)、収束する場合(YES)には、GPS受信装置に
設定された状態パラメータの更新を行なう(ステップS
153)。また、ステップS152において、収束しな
い場合(NO)、航法方程式の逐次解法を3回より多く
実行したか否かを判定し(ステップS154)、3回よ
り多い場合(YES)には、異常処理を行なう(ステッ
プS155)。
FIG. 5 shows the details of the processing in step S15. That is, the positioning calculation unit 300
By solving the navigation equation as shown in the above equation (2), error estimation values Δz ′, Δz ″ G at which the estimation probability becomes the maximum are obtained (step S151).
It is determined whether the observation data input to the observation data estimated based on Δz ″ G converges (step S15).
2) In the case of convergence (YES), the state parameters set in the GPS receiver are updated (step S).
153). In step S152, if the convergence does not occur (NO), it is determined whether the sequential solution of the navigation equation has been performed more than three times (step S154). (Step S155).

【0035】以上のように、上記実施形態によれば、測
位演算部300にて可視GPS衛星全てからの観測デー
タが取り込まれ、これら観測データを用いて、衛星信号
に与えられる誤差値もしくは自装置に発生する誤差値の
推定確率に最大を与える誤差推定値が算出され、この誤
差推定値に基づいて、可視GPS衛星全ての観測データ
から必要なGPS衛星からの観測データが決定される。
すなわち、測位演算部300は、可視GPS衛星全ての
観測データを入力して、誤差値の推定確率に最大を与え
る誤差推定値を算出し、この誤差推定値に基づいて推定
される観測データと実際の可視GPS衛星全ての観測デ
ータとを比較照合することにより、観測データの品質を
評価でき、この結果を基に所定値以下のずれの観測デー
タを信頼性の高い観測データとして選択でき、所定値よ
り大きいずれの観測データに対しては信頼性の低い観測
データと判断し棄却できるようになる。
As described above, according to the above-described embodiment, the positioning calculation unit 300 fetches observation data from all visible GPS satellites, and uses these observation data to obtain an error value given to a satellite signal or its own device. An error estimation value that gives the maximum to the estimation probability of the error value that occurs in is calculated, and based on the error estimation value, required observation data from GPS satellites is determined from observation data of all visible GPS satellites.
That is, the positioning calculation unit 300 receives the observation data of all the visible GPS satellites, calculates an error estimation value that gives the maximum to the estimation probability of the error value, and compares the observation data estimated based on the error estimation value with the actual observation data. The quality of the observation data can be evaluated by comparing and checking the observation data of all visible GPS satellites, and based on the result, the observation data with a deviation of a predetermined value or less can be selected as highly reliable observation data, and the predetermined value can be selected. Observation data with a larger deviation can be judged as unreliable observation data and rejected.

【0036】このため、通常の航法においてみられる、
選択衛星の切り替わりに起因する航法/測位のデータの
不連続性が最小限に抑えられ、使用可能な観測データを
最適利用しているので航法/測位精度及び航法/測位に
おける信頼性の向上を図ることが可能となる。さらにマ
ルチパス等に起因する一部の観測データの品質劣化に対
する耐性が得られる。
For this reason, it can be seen in normal navigation,
Discontinuity of navigation / positioning data due to switching of selected satellites is minimized, and optimal use of available observation data improves navigation / positioning accuracy and reliability in navigation / positioning. It becomes possible. Furthermore, resistance to quality deterioration of some observation data due to multipath or the like can be obtained.

【0037】また、差分航法を行なう場合には、上記式
(1)に示すように、ラグランジュ要素項L4を付加す
ることで確保されるため、差分航法への拡張が簡単とな
る。さらに、上記式(1)を構成する各項が観測残差の
2次形式で与えられるため、可視GPS衛星数が4より
大きい場合に、εG=0とすることで、信頼性の低い観
測データを棄却することが可能となる。
Further, when performing the differential navigation, since it is secured by adding the Lagrange element term L4 as shown in the above equation (1), the extension to the differential navigation is simplified. Further, since each term constituting the above equation (1) is given in a quadratic form of the observation residual, when the number of visible GPS satellites is larger than 4, ε G = 0 is set, so that low reliability observation is performed. Data can be rejected.

【0038】また、上記実施形態によるGPS受信装置
を用いれば、現在使用されている自動車、船舶の航法、
及び土木、建築、測量の測位等に即座に実施可能とな
り、高精度な測位、航法を実現できるようになる。
Further, if the GPS receiver according to the above embodiment is used, the navigation of the currently used automobile and ship,
In addition, the present invention can be immediately implemented for civil engineering, construction, surveying positioning, and the like, and highly accurate positioning and navigation can be realized.

【0039】なお、上記実施形態では、誤差の推定値を
算出する演算アルゴリズムとして、トータルラグランジ
ュの方程式を利用しているが、この他に推定誤差値を最
小にする最小自乗法を用いても同様に実施できる。
In the above embodiment, the total Lagrangian equation is used as an operation algorithm for calculating the estimated value of the error, but the least square method for minimizing the estimated error value is also used in the same manner. Can be implemented.

【0040】その他、GPS受信装置の構成、誤差の推
定値を算出する演算アルゴリズム等についても、この発
明の要旨を逸脱しない範囲で種々変形して実施できる。
In addition, the configuration of the GPS receiver, the operation algorithm for calculating the estimated value of the error, and the like can be variously modified without departing from the gist of the present invention.

【0041】[0041]

【発明の効果】以上詳述したようにこの発明によれば、
選択衛星の切り替わりに起因する航法/測位のデータの
不連続性を最小限に抑制でき、かつ航法/測位精度及び
航法/測位における信頼性の向上を図り得るGPS受信
装置及びこの装置に適用されるGPS受信方法を提供す
ることができる。
As described in detail above, according to the present invention,
The present invention is applied to a GPS receiver and a GPS receiver capable of minimizing discontinuity of navigation / positioning data due to switching of a selected satellite and improving navigation / positioning accuracy and reliability in navigation / positioning. A GPS receiving method can be provided.

【図面の簡単な説明】[Brief description of the drawings]

【図1】この発明に係るGPS受信装置の一実施形態を
示す回路ブロック図。
FIG. 1 is a circuit block diagram showing an embodiment of a GPS receiver according to the present invention.

【図2】同実施形態におけるGPS受信装置と複数のG
PS衛星との位置関係及び測位概念を示す図。
FIG. 2 shows a GPS receiver and a plurality of G's in the embodiment.
The figure which shows the positional relationship with a PS satellite, and a positioning concept.

【図3】上記図1に示した測位演算部の動作を示すフロ
ーチャート。
FIG. 3 is a flowchart showing the operation of the positioning calculation unit shown in FIG. 1;

【図4】上記図3に示したステップS14の詳細な処理
を示すフローチャート。
FIG. 4 is a flowchart showing a detailed process of step S14 shown in FIG. 3;

【図5】上記図3に示したステップS15の詳細な処理
を示すフローチャート。
FIG. 5 is a flowchart showing a detailed process of step S15 shown in FIG. 3;

【符号の説明】[Explanation of symbols]

11…アンテナ、12…低雑音増幅器、13…乗算器、
14…局部発振器、15…てい倍器、16…PNコード
解除器、17…PNコード発生器、18…位相比較器、
19…NCO(Newmedical Control Ocial)、101〜1
0N…ベースバンド処理部、200…受信制御部、30
0…測位演算部、401〜404…GPS衛星。
11 antenna, 12 low-noise amplifier, 13 multiplier
14 local oscillator, 15 multiplier, 16 PN code canceller, 17 PN code generator, 18 phase comparator,
19 ... NCO (Newmedical Control Ocial), 101-1
0N: baseband processing unit, 200: reception control unit, 30
0: Positioning calculation unit, 401 to 404: GPS satellites.

Claims (14)

【特許請求の範囲】[Claims] 【請求項1】 GPS(Global Positioning System)衛
星を利用して、自装置の位置を計測するGPS受信装置
において、 可視GPS衛星全てからの衛星信号を受信する受信手段
と、 この受信手段で得られた可視GPS衛星全ての衛星信号
から可視GPS衛星に関する擬似距離データもしくは搬
送波位相データと、可視GPS衛星の放送暦データとを
含む観測データを抽出する観測データ抽出手段と、 この観測データ抽出手段で得られた可視GPS衛星全て
の観測データを入力して、前記衛星信号に与えられる誤
差値もしくは自装置に発生する誤差値の推定確率に最大
を与える誤差の推定値を算出する誤差推定値算出手段
と、 この誤差推定値算出手段により算出された推定値に基づ
いて、前記受信手段で得られた可視GPS衛星全ての観
測データの中から必要なGPS衛星の観測データを決定
する観測データ決定手段とを具備してなることを特徴と
するGPS受信装置。
1. A GPS receiving device that measures the position of its own device using a GPS (Global Positioning System) satellite, receiving means for receiving satellite signals from all visible GPS satellites, Observation data extraction means for extracting observation data including pseudorange data or carrier phase data relating to the visible GPS satellites and broadcast calendar data of the visible GPS satellites from all satellite signals of the visible GPS satellites, Error estimation value calculation means for inputting observation data of all the visible GPS satellites and calculating an error value given to the satellite signal or an error estimation value which gives a maximum to an estimation probability of an error value occurring in the own device; Based on the estimated value calculated by the error estimated value calculating means, the observation data of all visible GPS satellites obtained by the receiving means GPS receiver apparatus characterized by comprising; and a monitoring data determining means for determining the observation data of al necessary GPS satellites.
【請求項2】 前記誤差推定値算出手段は、可視GPS
衛星全ての観測データを入力して、自装置の位置の更新
前位置からの偏差及び自装置のクロック誤差を示す推定
値を算出する手段を有することを特徴とする請求項1記
載のGPS受信装置。
2. The method according to claim 1, wherein the error estimation value calculating means is a visible GPS.
2. The GPS receiving apparatus according to claim 1, further comprising means for inputting observation data of all satellites and calculating an estimated value indicating a deviation of the position of the own apparatus from the pre-update position and a clock error of the own apparatus. .
【請求項3】 前記誤差推定値算出手段は、可視GPS
衛星全ての観測データを入力して、可視GPS衛星の位
置・クロック誤差に関する放送暦データからの偏差に地
球の自転による補正係数を加えた値を示す推定値を算出
する手段を有することを特徴とする請求項1記載のGP
S受信装置。
3. The apparatus according to claim 2, wherein said error estimation value calculating means is a visible GPS.
Means for inputting observation data of all the satellites and calculating an estimated value indicating a value obtained by adding a correction coefficient due to the rotation of the earth to a deviation from the broadcast calendar data regarding the position / clock error of the visible GPS satellite. The GP according to claim 1,
S receiver.
【請求項4】 前記誤差推定値算出手段は、さらに可視
GPS衛星の位置・クロックの誤差値の推定確率に最大
を与える推定値を算出する手段を有することを特徴とす
る請求項1記載のGPS受信装置。
4. The GPS according to claim 1, wherein said error estimated value calculating means further comprises means for calculating an estimated value that gives a maximum to an estimated probability of an error value of a position / clock of a visible GPS satellite. Receiver.
【請求項5】 前記受信手段は、予め決められた位置に
設置された地上局から送出され、地上局による擬似距離
データ及び観測時刻データを含む地上局信号を受信する
手段を有し、 前記観測データ抽出手段は、前記受信手段で受信された
前記地上局信号から前記擬似距離データ及び観測時刻デ
ータを抽出する手段を有し、 前記誤差推定値算出手段は、前記地上局による擬似距離
データ及び観測時刻データを入力して、前記地上局によ
る観測衛星位置の誤差値の推定確率に最大を与える推定
値を算出する手段を有することを特徴とする請求項1記
載のGPS受信装置。
5. The receiving means has means for receiving a ground station signal transmitted from a ground station installed at a predetermined position and including pseudorange data and observation time data by the ground station. The data extracting means has means for extracting the pseudorange data and the observation time data from the ground station signal received by the receiving means, and the error estimation value calculating means has the pseudorange data and the observation data obtained by the ground station. 2. The GPS receiving apparatus according to claim 1, further comprising means for inputting time data and calculating an estimated value which gives a maximum to an estimated probability of an error value of an observation satellite position by the ground station.
【請求項6】 前記誤差推定値算出手段は、可視GPS
衛星全ての観測データを入力して、ラグランジュの方程
式を用いた演算アルゴリズムにより、誤差の推定値を算
出することを特徴とする請求項1記載のGPS受信装
置。
6. The method according to claim 6, wherein the error estimation value calculating means includes a visible GPS
2. The GPS receiving apparatus according to claim 1, wherein the observation data of all the satellites is input, and an estimated value of the error is calculated by an arithmetic algorithm using Lagrange's equation.
【請求項7】 前記誤差推定値算出手段は、可視GPS
衛星全ての観測データを入力して、最小自乗法を用いた
演算アルゴリズムにより、誤差の推定値を算出すること
を特徴とする請求項1記載のGPS受信装置。
7. The apparatus according to claim 1, wherein said error estimation value calculating means includes a visible GPS
2. The GPS receiving apparatus according to claim 1, wherein the observation data of all the satellites is input, and an estimated value of the error is calculated by an arithmetic algorithm using a least square method.
【請求項8】 GPS衛星を利用して、自装置の位置を
計測するGPS受信装置に適用されるGPS受信方法に
おいて、 可視GPS衛星全てから衛星信号を受信し、 得られた可視GPS衛星全ての衛星信号から可視GPS
衛星に関する擬似距離データもしくは搬送波位相データ
と、可視GPS衛星の放送暦データとを含む観測データ
を抽出し、 得られた可視GPS衛星全ての観測データを用いて、前
記衛星信号に与えられる推定誤差値もしくは自装置に発
生する誤差値の推定確率に最大を与える推定値を算出
し、 この算出された推定値に基づいて、前記可視GPS衛星
全ての観測データの中から必要なGPS衛星の観測デー
タを決定することを特徴とするGPS受信方法。
8. A GPS receiving method applied to a GPS receiving device that measures the position of its own device using GPS satellites, wherein satellite signals are received from all visible GPS satellites, and all of the obtained visible GPS satellites are received. Visible GPS from satellite signals
Extraction of observation data including pseudorange data or carrier phase data relating to satellites and broadcast calendar data of visible GPS satellites, and using the obtained observation data of all visible GPS satellites, an estimated error value given to the satellite signal Alternatively, an estimated value that gives the maximum to the estimated probability of the error value generated in the own device is calculated. Based on the calculated estimated value, observation data of the necessary GPS satellites is obtained from observation data of all the visible GPS satellites. A GPS receiving method, characterized by determining.
【請求項9】 前記推定値は、GPS受信装置の位置の
更新前位置からの偏差及びGPS受信装置のクロック誤
差を示すことを特徴とする請求項8記載のGPS受信方
法。
9. The GPS receiving method according to claim 8, wherein the estimated value indicates a deviation of the position of the GPS receiver from the pre-update position and a clock error of the GPS receiver.
【請求項10】 前記推定値は、可視GPS衛星の位置
・クロック誤差に関する放送暦データからの偏差に地球
の自転による補正係数を加えた値を示すことを特徴とす
る請求項8記載のGPS受信方法。
10. The GPS receiving apparatus according to claim 8, wherein the estimated value indicates a value obtained by adding a correction coefficient due to the rotation of the earth to a deviation from broadcast calendar data relating to a position / clock error of a visible GPS satellite. Method.
【請求項11】 前記推定値は、さらに可視GPS衛星
の位置・クロックの誤差値の推定確率が最大となるよう
に算出されることを特徴とする請求項8記載のGPS受
信方法。
11. The GPS receiving method according to claim 8, wherein the estimated value is further calculated such that the estimated probability of the error value of the position / clock of the visible GPS satellite is maximized.
【請求項12】 前記GPS受信装置が予め決められた
位置に設置された地上局から送出され、地上局による擬
似距離データ及び観測時刻データを含む地上局信号を受
信する場合に、 受信された前記地上局信号から前記擬似距離データ及び
観測時刻データを抽出し、 この抽出した地上局による擬似距離データ及び観測時刻
データを入力して、前記地上局による観測衛星位置の誤
差値の推定確率に最大を与える推定値を算出することを
特徴とする請求項8記載のGPS受信方法。
12. When the GPS receiving apparatus receives a ground station signal including pseudo distance data and observation time data transmitted from a ground station installed at a predetermined position and received by the ground station, the GPS receiver receives the received GPS signal. The pseudo-range data and observation time data are extracted from the ground station signal, and the extracted pseudo-range data and observation time data by the ground station are input, and the estimation probability of the error value of the observation satellite position by the ground station is maximized. The GPS receiving method according to claim 8, wherein an estimated value to be given is calculated.
【請求項13】 前記推定値は、可視GPS衛星全ての
観測データを入力して、ラグランジュの方程式を用いた
演算アルゴリズムにより、算出されることを特徴とする
請求項8記載のGPS受信方法。
13. The GPS receiving method according to claim 8, wherein the estimated value is calculated by inputting observation data of all visible GPS satellites and using an arithmetic algorithm using Lagrange's equation.
【請求項14】 前記推定値は、可視GPS衛星全ての
観測データを入力して、最小自乗法を用いた演算アルゴ
リズムにより、算出されることを特徴とする請求項8記
載のGPS受信方法。
14. The GPS receiving method according to claim 8, wherein the estimated value is calculated by inputting observation data of all visible GPS satellites and using an arithmetic algorithm using a least squares method.
JP7848599A 1999-03-23 1999-03-23 GPS receiving device and GPS receiving method applied to this device Expired - Fee Related JP3544887B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP7848599A JP3544887B2 (en) 1999-03-23 1999-03-23 GPS receiving device and GPS receiving method applied to this device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP7848599A JP3544887B2 (en) 1999-03-23 1999-03-23 GPS receiving device and GPS receiving method applied to this device

Publications (2)

Publication Number Publication Date
JP2000275320A true JP2000275320A (en) 2000-10-06
JP3544887B2 JP3544887B2 (en) 2004-07-21

Family

ID=13663297

Family Applications (1)

Application Number Title Priority Date Filing Date
JP7848599A Expired - Fee Related JP3544887B2 (en) 1999-03-23 1999-03-23 GPS receiving device and GPS receiving method applied to this device

Country Status (1)

Country Link
JP (1) JP3544887B2 (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007248345A (en) * 2006-03-17 2007-09-27 Toyota Motor Corp Positioning method and positioning device
CN103900581A (en) * 2014-04-04 2014-07-02 哈尔滨工程大学 MIMU and GPS combined pedestrian navigation method based on augmented lagrangian condition
CN108563108A (en) * 2017-12-23 2018-09-21 航天恒星科技有限公司 A kind of GNSS satellite clock health status monitoring method

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007248345A (en) * 2006-03-17 2007-09-27 Toyota Motor Corp Positioning method and positioning device
CN103900581A (en) * 2014-04-04 2014-07-02 哈尔滨工程大学 MIMU and GPS combined pedestrian navigation method based on augmented lagrangian condition
CN103900581B (en) * 2014-04-04 2017-01-18 哈尔滨工程大学 MIMU and GPS combined pedestrian navigation method based on augmented lagrangian condition
CN108563108A (en) * 2017-12-23 2018-09-21 航天恒星科技有限公司 A kind of GNSS satellite clock health status monitoring method

Also Published As

Publication number Publication date
JP3544887B2 (en) 2004-07-21

Similar Documents

Publication Publication Date Title
US6950059B2 (en) Position estimation using a network of a global-positioning receivers
EP1498745B1 (en) Enhanced real time kinematics determination method and apparatus
EP2044457B1 (en) A method for increasing the reliability of position information when transitioning from a regional, wide-area, or global carrier-phase differential navigation (wadgps) to a local real-time kinematic (rtk) navigation system
RU2543618C2 (en) Method to determine position by linkage of measurements
EP2864809B9 (en) Selection of a subset of global navigation satellite system measurements based on prediction of accuracy of target parameters
EP2864810B1 (en) Selection of a subset of global navigation satellite system measurements based on relation between shifts in target parameters and sum of residuals
KR100916875B1 (en) Method and apparatus for determining location using a coarse position estimate
US6587692B1 (en) Location determination using weighted ridge regression
CN107710017A (en) For the satellite navigation receiver and method switched between real time kinematics pattern and relative positioning mode
EP3109672B1 (en) Gnss receiver with a capability to resolve ambiguities using an uncombined formulation
US20240142639A1 (en) Terminal positioning method and apparatus, device, and medium
US20080180315A1 (en) Methods and systems for position estimation using satellite signals over multiple receive signal instances
EP3106897A1 (en) Gnss receiver with an on-board capability to implement an optimal error correction mode
Bourdeau et al. Constructive use of GNSS NLOS-multipath: Augmenting the navigation Kalman filter with a 3D model of the environment
US10194269B2 (en) Systems and methods for using doppler measurements to estimate a position of a receiver
Houzeng et al. GPS/BDS/INS tightly coupled integration accuracy improvement using an improved adaptive interacting multiple model with classified measurement update
US11209552B2 (en) Method and apparatus for improving the quality of position determination using GNSS data
WO2015145718A1 (en) Positioning device
EP1506430A1 (en) Determining the position of a receiver and/ or the system time of a positioning system
JP5636552B2 (en) How to optimize the acquisition of spread spectrum signals from satellites by mobile receivers
CN113109840A (en) Ionosphere TEC real-time measurement method based on GNSS receiver
JP4592526B2 (en) Positioning system
US20030054759A1 (en) Integrity monitoring for geo-location systems
JP2009079975A (en) Positioning system
JP3544887B2 (en) GPS receiving device and GPS receiving method applied to this device

Legal Events

Date Code Title Description
A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20040202

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20040330

A61 First payment of annual fees (during grant procedure)

Effective date: 20040406

Free format text: JAPANESE INTERMEDIATE CODE: A61

LAPS Cancellation because of no payment of annual fees