JP4729197B2 - 物体の姿勢検出装置および整数バイアス再決定方法 - Google Patents
物体の姿勢検出装置および整数バイアス再決定方法 Download PDFInfo
- Publication number
- JP4729197B2 JP4729197B2 JP2001167407A JP2001167407A JP4729197B2 JP 4729197 B2 JP4729197 B2 JP 4729197B2 JP 2001167407 A JP2001167407 A JP 2001167407A JP 2001167407 A JP2001167407 A JP 2001167407A JP 4729197 B2 JP4729197 B2 JP 4729197B2
- Authority
- JP
- Japan
- Prior art keywords
- antenna
- phase difference
- sensor unit
- inertial sensor
- antennas
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Navigation (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Description
【発明の属する技術分野】
この発明は航空機や船舶等の移動体や、海上浮揚構造体、超高層建築物のような姿勢が変化する物体の姿勢を検出する装置に関するものである。
【0002】
【従来の技術】
従来、複数の測位用衛星からの電波を、物体の所定位置に配置された複数のアンテナで受信し、キャリア信号の位相をそれぞれ観測し、複数のアンテナのうち1つのアンテナを基準として他のアンテナの相対位置を求めることによって、受信アンテナの相対位置関係を求め、これによって移動体の3次元姿勢を計測するようにした姿勢検出装置(Attitude Determining System:ADS)が研究開発されている。
【0003】
一方、レートジャイロ等の角速度を検出する慣性センサユニット(以下、IMU:Inertia Measurement Unitという。)を用いて移動体の姿勢を測定するようにした装置も用いられている。
【0004】
【発明が解決しようとする課題】
従来のADSでは、移動体の所定位置に配置した複数のアンテナのうち1つを基準として他のアンテナの相対位置を求めるために、測位用衛星から送信された電波のキャリア位相について2重位相差または1重位相差を観測して2重位相差または1重位相差に含まれる整数バイアスを求めるが、従来の整数バイアス解法では統計確率的に誤りが生じ、正しい整数バイアスを見つける確率は例えば約85〜95%である。また、整数バイアスの決定時間は観測雑音の大きさ、衛星数、およびアンテナ基線長(複数のアンテナのうち1つの基準アンテナと他の相対位置を求めるべきアンテナとを結ぶ直線の長さ)に大きく依存する。通常の雑音レベルにおいて、基線長が1m程度の場合には、平均30秒から数分の時間を要する。また、雑音レベルが高い場合には、たとえば10分程度の長時間をかけても整数バイアスを決定できないという場合もあった。
【0005】
実用の環境下では、測位に用いる衛星の組合せを変更する場合や、雑音によってサイクルスリップ(キャリア再生回路の位相同期ループが追尾不能になることにより、キャリア位相をカウントするカウンタの整数分の値が狂う現象)が生じた場合等、高頻度で整数バイアスを再決定する必要がある。そのため、整数バイアスの決定時間が長いと、ユーザーへの姿勢出力が高頻度に、しかも許容し難い期間、中断してしまう。また、求められた整数バイアスの正解率が低いと、信頼できない姿勢情報を提供してしまう。これらがADSの実用化のネックとなっていた。
【0006】
一方、IMUのみを用いた姿勢検出装置では、x,y,zの3軸成分の角速度を積分することによって三次元の姿勢を検出するものであるため、誤差の累積が問題となっていた。
【0007】
そこで、米国特許5,349,531 では、従来のADSとIMUによる姿勢検出システムとの統合化を図り、ADSによる姿勢出力が中断したときに、IMUの出力に基づく角で補間するようにしている。尚、以降では特に断らない限り、IMUの出力に基づく姿勢はローカル座標系基準の姿勢角を意味するものとし、IMU姿勢角と呼ぶことにする。
【0008】
また、米国特許4,754,280 の姿勢検出装置では、IMUとGPSを統合することによりIMUの誤差を推定し、高レート出力、かつ高精度な姿勢を提供することのみにとどまっている。
【0009】
IMU利用により移動体の姿勢検出を行う装置で、安価な普及タイプのIMUを使用した場合には、バイアス、ドリフト、スケール・ファクタ等のセンサ誤差が大きく、しかも正確な誤差モデリングができないために、これらの誤差を完全に推定し、補正することができない。そのため、例えば姿勢角として約1°を許容誤差とする場合には、IMU姿勢角によるデドレコ(推測航法)可能なバックアップ時間は約30秒〜60秒程度が限界となる。しかし、キャリア位相の整数バイアスのカウント中にサイクルスリップが生じた場合、このような短時間のうちに整数バイアスを再決定することは、従来の方法では困難であり、実用的ではなかった。
【0010】
上記いずれの米国特許も、ADSにおけるキャリア位相差の整数バイアスを求めることに伴う上述の問題を解消するものではない。
【0011】
この発明の目的は、IMU姿勢角を利用して、より短時間に整数バイアスを確実に再決定するようにした物体の姿勢検出装置および整数バイアス再決定方法を提供することにある。
【0012】
【課題を解決するための手段】
この発明は、キャリア位相差の整数バイアスを決定して、複数のアンテナの相対位置を測位している途中でサイクルスリップが生じたとき、IMUの姿勢角を利用して整数バイアスを直接決定することにより、サイクルスリップが頻繁に生じても、ユーザーに姿勢角を中断することなく出力できるようにする。
【0013】
また、この発明は、IMU姿勢角から整数バイアスを決定するために、IMU姿勢角と、アンテナ座標系における基線ベクトルとからローカル座標系における基線ベクトルを求め、アンテナおよび測位用衛星の位置で決まる、アンテナから測位用衛星への視野方向ベクトル、観測による1重位相差または2重位相差、および前記ローカル座標系における基線ベクトルに基づいて決定する。
【0014】
また、この発明は、物体上の、一直線上に並ばないそれぞれ異なった位置に固定配置された複数のアンテナで複数の測位用衛星からの電波をそれぞれ受信し、キャリア位相の1重位相差または2重位相差を観測し、該1重位相差または2重位相差の整数バイアスを初期化時において決定するステップと、物体に配置されたIMU姿勢角と、アンテナ座標系における基線ベクトルとからローカル座標系における基線ベクトルを求めるステップと、前記アンテナおよび測位用衛星の位置で決まるアンテナから測位用衛星への視野方向ベクトル、前記観測による1重位相差または2重位相差、および前記ローカル座標系における基線ベクトルに基づいて、前記1重位相差または2重位相差の整数バイアスを直接再決定するステップとを備える。
【0015】
【発明の実施の形態】
まず、アンテナ座標系とローカル座標系について、および両者の座標変換について、図1を参照して説明する。
【0016】
図1の(A)はローカル座標系を表し、基準アンテナの位置を原点として例えばx軸を北、y軸を東、z軸を鉛直方向(下方)に採る。(x 1 n ,y 1 n ,z 1 n )は第2アンテナの位置、(x 2 n ,y 2 n ,z 2 n )は第3アンテナの位置、L 1 n は第1アンテナから第2アンテナへの基線ベクトル(複数のアンテナのうち1つの基準アンテナと他の相対位置を求めるべきアンテナとを結ぶ直線の長さと方向を表すベクトル)、L 2 n は第1アンテナから第3アンテナへの基線ベクトルである。尚、座標値および基線ベクトルの右肩の添字nはローカル座標系での座標値および基線ベクトルであることを表している。
【0017】
図1の(B)はアンテナ座標系を表し、基準アンテナの位置を原点として例えばx軸を方位(機首または船首の方向)に採り、y軸をx軸方向に直交する水平方向に採り、z軸をx,y軸にそれぞれ直交する方向に採る。(x 1 a ,y 1 a ,z 1 a )は第2アンテナの位置、(x 2 a ,y 2 a ,z 2 a )は第3アンテナの位置、L 1 a は第1アンテナから第2アンテナへの基線ベクトル、L 2 a は第1アンテナから第3アンテナへの基線ベクトルである。尚、座標値および基線ベクトルの右肩の添字aはアンテナ座標系での座標値および基線ベクトルであることを表している。
【0018】
上記アンテナ座標系からローカル座標系への座標変換の変換式は〔数1〕で示される。
【0019】
【数1】
【0020】
ここでL 3 n はL 1 n とL 2 n の外積である。また、L 3 a はL 1 a とL 2 a の外積である。
GPSで計測されるローカル座標系に対するアンテナ座標系の姿勢角を、或る軸回転順序で定義されたオイラー角(θr,θp,θy)で表現すると、〔数1〕における座標変換オペレータC a n は次式で表される。
【0021】
【数2】
【0022】
ここで、θrはロール角、θpはピッチ角、θyはヨー角と言うことができる。尚、GPSによる基線ベクトルは、ローカル座標系で得られるので、姿勢角θyは、真北基準の方位として得られる。また、θrとθpは水平面を基準とする姿勢角として得られる。尚、以降、(θr,θp,θy)をGPS姿勢角と呼ぶ。
【0023】
上記3つの姿勢角はそれぞれ次の関係から求められる。
【0024】
【数3】
【0025】
【数4】
【0026】
【数5】
【0027】
ここで、〔数1〕において、(L 1 a ,L 2 a ,L 3 a )は既知であるから、(L 1 n ,L 2 n ,L 3 n )をキャリア位相差観測量により求めると、C a n が求まる。C a n が求まると、〔数3〕〜〔数5〕で姿勢角が求まる。従って、姿勢角の計測精度は、基線ベクトルの算出精度、すなわち(x 1 n ,y 1 n ,z 1 n ),(x 2 n ,y 2 n ,z 2 n )の算出精度によって決定される。そのため、通常は擬似雑音コードによる測位ではなく、キャリア位相差を観測して相対位置を高精度に計測する。その際、1重位相差方式、2重位相差方式のいずれでも可能であるが、以下2重位相差方式について述べる。
2重位相差方式による観測方程式は次式で与えられる。
【0028】
【数6】
【0029】
ここで、
m:衛星番号
k:衛星数
j:基線ベクトルの番号
λ:受信信号キャリア波長
δρ j m :基準衛星S1 と衛星Sm+1 (m=1,2,3)に対する基線ベクトルL j n の2重距離差
δφ j m :同上に対する2重位相差観測量
δN j m :同上に対する2重整数バイアス差
δe j m :同上に対する2重観測雑音差
この2重位相差を図式化したものを図2に示す。図中の記号の意味は次のとおりである。
【0030】
Δρ j m+1 :衛星sm+1 に対する基線ベクトルL j n の1重距離差
Δφ j m+1 :衛星sm+1 に対する基線ベクトルL j n の1重位相観測量
ΔN j m+1 :衛星sm+1 に対する基線ベクトルL j n の1重整数バイアス差
Δe j m+1 は衛星sm+1 に対する基線ベクトルL j n の1重観測雑音差
また、図中の各記号と前記各変数との関係は次のとおりである。
【0031】
δρ j m =Δρ j m+1 −Δρ j 1
δφ j m =Δφ j m+1 −Δφ j 1
δN j m =ΔN j m+1 −ΔN j 1
δe j m =Δe j m+1 −Δe j 1
〔数6〕において2重距離差δρ j m は、δρ j m =δA j m L j n (j=1,2)で表される。ここで、δA j m は2重位相差における衛星への方向余弦マトリクス差であり、アンテナ位置と衛星位置により決定できる。
【0032】
また、L j n は(x 1 n ,y 1 n ,z 1 n )の関数である。したがって、〔数6〕右辺の未知数はδN j m とL j n すなわち(x j n ,y j n ,z j n )となる。いま、前もって何らかの方法で整数バイアスを求めれば、未知数は(x j n ,y j n ,z j n )の3個であるので、少なくとも3個の観測量δφ j m が得られると、すなわち衛星数が4以上であると、(x j n ,y j n ,z j n )を求めることができる。
【0033】
さて、以下に、IMU姿勢角から整数バイアスを再決定できる根拠を示す。
いま、図1のように定義されたローカル座標系とアンテナ座標系における任意の基線ベクトルLの座標ベクトルの座標位置を、それぞれ、マトリックスXn ,Xa で表記すれば、〔数1〕は次の式で表される。
【0034】
【数7】
Xn =C a n Xa …(7)
したがって、Xn の誤差共分散σP 2 =E[ δXn ( δXn ) T ] は次式で示される。
【0035】
【数8】
【0036】
ここで、σ p 2,σθ2 はそれぞれアンテナ基線ベクトル位置と姿勢角の共分散誤差マトリックスであり、次式で表される。また、HはXn と姿勢角との感度マトリックスである。
【0037】
【数9】
【0038】
基線ベクトル位置の誤差の標準偏差σP と衛星視野方向に射影した値δρとの関係は次式で表される。
【0039】
【数10】
δρ=δAσp …(10)
ここでδAは、〔数6〕の説明中で述べたδA j m (2重位相差における衛星への方向余弦マトリクス差)のサフィックスを省略したものである。
いま、δρの標準偏差をσN とすると、上記〔数8〕と〔数10〕とによって、姿勢角誤差と整数バイアスとの誤差との関係は次式で表される。
【0040】
【数11】
【0041】
今、仮に、各辺が約50cmの三角形の各頂点に配置されたアンテナの架台を、ロール角±20°、ロール角変動周期60秒、ピッチ角±15°、ピッチ角変動周期70秒、ヨー角±30°、ヨー角変動周期180秒で振揺させたとき、そして衛星の方位を360°、仰角を0〜80°の範囲でそれぞれ変化させたとき、IMU姿勢角の誤差が標準偏差で数度以内であれば、衛星視野方向換算の距離誤差は、衛星から送信される電波のキャリア波長の1/4(5cm程度)以下となる。
【0042】
このことは、姿勢角誤差は通常数度以内の精度であるから、上記の条件でIMUの出力に基づいて求めた姿勢角からキャリア位相差のアンビギュイティを丸め込みによって整数化しても、キャリア位相差の整数バイアスを誤りなく決定できることを意味する。
【0043】
次に、IMU姿勢角より整数バイアスを再決定する方法を具体的に記述する。
アンテナ座標系における基線ベクトルLj a(j=1,2)の位置Xa (xj a,yj a,zj a)は既知であるから、IMU姿勢角を用いて〔数2〕、〔数7〕より、ローカル座標系における基線ベクトルLj n (j=1,2)の位置Xn(xj n,yj n,zj n)を求める。
【0044】
この(x j n ,y j n ,z j n )を用いて、次式より各衛星における1重距離差Δρ j m を逆算する。
【0045】
【数12】
Δρ j m =A j m X j n (m=1,2,3,4,j=1,2) …(12)
ここで、A j m は基線ベクトルL j n のm衛星への視野方向ベクトルであるから、ローカル座標系における各衛星位置とアンテナの絶対位置より決定できる。尚、mは使用衛星番号、jは基線ベクトルの番号である。
【0046】
1重位相差の観測方程式は次式で表されるから、1重位相差Δφ j m を観測量とすると、m衛星に対する1重位相差の実数アンビギュイティΔN j m が計算できる。
【0047】
【数13】
λΔN j m =Δρ j m −λΔφ j m +Δe j m …(13)
(m=1,2,3,4,j=1,2)
前述の例示のようにIMU姿勢角誤差の衛星視野方向への距離換算誤差は、通常、キャリア波長の1/4(5cm程度)以下となるので、このΔN j m を丸め込みによって、1重位相差の整数バイアス(1重整数バイアス差)が決定できる。
【0048】
尚、2重位相差方式の場合には、前述の〔数6〕より2重位相差実数アンビギュイティδN j m で計算できる。したがって、1重位相差方式、2重位相差方式のいずれにも適用できる。
【0049】
以上のようにして、IMU姿勢角を用いて整数バイアスを決定することができる。
【0050】
尚、以上に述べた例では、3つのアンテナを用い、そのうちの1つを基準アンテナとし、他の2つのアンテナとの間の組で1重位相差または2重位相差を求めるようにしたが、2つのアンテナのみを設け、その2つのアンテナ間で1重位相差または2重位相差を求めるようにしてもよい。例えば、2つのアンテナ間の基線を船舶の船首尾線上に配置すれば、その2つのアンテナの相対位置から、船首方位(θy)および、ピッチ角(θp)を検出することができる。また、2つのアンテナ間の基線を右舷−左舷方向に配置すれば、その2つのアンテナの相対位置から、船首方位(θy)とロール角(θr)を検出することができる。
【0051】
但し、前者ではロール角(θr)を検出できず、後者ではピッチ角(θp)を検出できない。3次元の姿勢が求まらなければ、〔数1〕および〔数2〕より、この場合は、IMU出力から得られる相対姿勢角を完全にGPS姿勢角で初期化できないために、IMU姿勢角から整数バイアスを求めることができない。そこで、本発明の方法を適用するためには、後述するように、2つのアンテナだけでは検出できない姿勢角を検出するための不足分の姿勢を計測可能な加速度センサ、または傾斜計を併用すればよい。
【0052】
また、4つ以上のアンテナを用いて、そのうちの1つを基準アンテナとし、他の3つのアンテナとの間の組で1重位相差または2重位相差をそれぞれ求めることによって、各アンテナ間の相対位置を求めるようにしてもよい。
【0053】
次に、物体の姿勢検出装置の全体の構成を図3および図4を参照して説明する。
図3は装置のブロック図である。ここでs1,s2・・・skはGPS測位衛星であり、GPSアンテナはこれらの複数の測位用衛星からの電波を受信し、中間周波信号に変換するとともに増幅する。GPS受信機はそれぞれのアンテナで受信した信号から、それぞれのアンテナ位置およびキャリア位相差を演算し、各衛星位置情報と共に姿勢演算処理部へ所定(数十ミリ秒〜数秒)の間隔で送出する。測位により求められるアンテナ位置や衛星のエフェメリス情報等はGPSの座標系に従っているので、GPS受信機部または姿勢演算処理部のいずれかでローカル座標系に変換する。以下の説明ではGPS受信機部で変換するものとする。
【0054】
姿勢演算処理部における衛星計画部は、衛星のエフェメリス情報とアンテナ位置とから、使用する衛星を選定し、その衛星位置や選定された衛星の衛星情報を位相差演算部へ与える。この位相差演算部は、衛星計画部から受け取った情報とGPS受信機より受け取ったキャリア位相差信号から1重(または2重)位相差の観測量を算出し、整数バイアス決定処理部へ送出する。整数バイアス決定部では、該当衛星に対する1重(または2重)位相差の整数バイアスを決定し、その妥当性を確認した後に整数バイアス情報をGPS姿勢演算部へ与える。GPS姿勢演算部では、整数バイアスと位相差演算部より送出されたキャリア位相差観測量とを用いて、公知の方法によって、ローカル座標系における基線ベクトルを算出する。これと既知のアンテナ座標系における基線ベクトルより、〔数1〕〜〔数5〕によって姿勢角(GPS姿勢角)を算出する。
【0055】
図3において、IMUは直交3軸回りの角速度を検出するレートジャイロと、それらの出力を増幅するX軸増幅部、Y軸増幅部、Z軸増幅部を備えていて、姿勢演算処理部のIMU姿勢演算部に対して各軸回りの角速度を与える。具体的には、地球の赤道平面と春分点子午面との交線をx軸、これに直交する東向きをy軸とする右手系の慣性座標軸基準に対するIMU座標軸の角速度ω s isを与える。いま、必要なのは、ローカル座標系に対するIMU座標軸の角速度ω s nsである。ω s isとω s nsとの関係は、ω s ns=ω s is−ω s inで表される。ここで、ω s inは慣性座標軸系基準に対するローカル座標軸の角速度をIMU座標軸系で観測した角速度である。船舶のような低速で移動する物体の場合ではω s is>>ω s inが成立するから、ω s inは無視できる。航空機のような高速移動する物体では、移動体の位置および周知の方法によりω s nsを計算する。
【0056】
IMU姿勢演算部では、周知のオイラー方程式またはGilmore のアルゴリズムによって姿勢角を算出する。尚、IMU姿勢角の初期値はGPS姿勢角を用いて初期化するものとする。
【0057】
IMU/GPS統合化処理演算部では、IMU姿勢角とGPS姿勢角とを統合化して、最終的にロール,ピッチ,ヨーの各姿勢角を出力する。
【0058】
図4は、図3に示した姿勢演算処理部における再初期化判定部およびIMU/GPS統合化処理演算部等での処理手順を示すフローチャートである。まず測位に用いる衛星の選択、および選択した衛星からの電波を受信して求めた観測データを取り込む(n1)。初期状態で未だ整数バイアスを決定していなければ、先ず各アンテナについて基線ベクトルを推定する。(実数解を求める。)この推定には例えばカルマン・フィルタまたは最小2乗法を用いる。続いて、公知の方法により実数解を整数化することによって整数バイアスを決定する(n2→n3→n5→n6→n5・・・)。
【0059】
その後、χ2 検定および基線長が実際の値に一致するか否か等によって、整数バイアス推定の妥当性の検定を行う(n7)。もし、検定結果が妥当でなければ、再び整数バイアスの決定処理を行う(n8→n5・・・)。検定結果が妥当とみなされれば、求められた整数バイアスを既知として基線ベクトルを求め、各アンテナの相対位置からGPS姿勢角を算出する(n9)。続いて、このGPS姿勢角を算定した同時刻におけるIMUの3軸回りの角速度からIMU姿勢角を演算する(n10)。
【0060】
尚、前述のようにIMU姿勢角の初期値はGPS姿勢角を用いて初期化されるので、IMU姿勢角はGPS姿勢角と同じローカル座標系の姿勢角に一致したものとなる。
【0061】
その後、GPS姿勢角とIMU姿勢角との統合化処理を行う。ここでは、GPSダウン時の処理や、GPS姿勢角とIMU姿勢角との差からIMUの誤差補正を行う(n10→n11)。
【0062】
続いて、姿勢演算処理部における再初期化判定部で、サイクルスリップが生じたか否かを判定し、サイクルスリップが生じていなければ、求められた姿勢角を出力する(n12→n13→n15)。
【0063】
サイクルスリップが生じたことを検出すれば、上述した方法により、IMU姿勢角データより整数バイアスを直接決定し、その整数バイアスで基線ベクトルを定める(n13→n14→n4)。このステップn4の基線ベクトル推定とステップn3の基線ベクトル推定との違いは、例えばカルマン・フィルタを用いる場合、ステップn3の基線ベクトル推定でのステート数はx,y,zの3個とm個のアンビギュイティであるが、ステップn4の基線ベクトル推定でのステート数はx,y,zの3個のみである。従って、このステップn4の基線ベクトル推定の処理負荷は小さく、極短時間に基線ベクトルを推定できる。
以降、既に述べた整数バイアスの検定以降の処理を行う(n7→n8→・・・)。
【0064】
尚、通常はステップn14での整数バイアスの決定に誤りがないので、ステップn7,n8の処理を省略して、ステップn4からステップn9へ進むようにしてもよい。
【0065】
以上のようにして、サイクルスリップ、衛星選択の組合せの変更等によって再初期化が必要となった場合に、IMU姿勢角データを基にして整数バイアスを直接決定するようにしたため、整数バイアスの決定に長時間を必要としない。そのため、物体の姿勢角を実質上中断することなく出力でき、高精度な姿勢角情報を連続的にユーザに提供できるようになる。
【0066】
次に、2つのアンテナを用いた物体の姿勢検出装置の構成を図5を参照して説明する。
図3に示した例では、3つのアンテナを用いてGPS姿勢角を求めたが、この図5に示す例では、3つの姿勢角θr,θp,θyのうち、2つ(θy、θrまたはθpのいずれか)の姿勢角を2つのアンテナを用いて求め、2アンテナ方式では算出できないθp、またはθrのいずれかの姿勢角を加速度センサを用いて求める。
【0067】
すなわち、図5において、GPS受信機部には2つのGPSアンテナでGPS衛星からの電波を受信し、キャリア位相を観測する。このGPS受信機部または姿勢演算処理部でθy、θrまたはθpのいずれかの姿勢角を求める。IMUには、レートジャイロの他に1つ以上の加速度センサを備えていて、この加速度センサの検出結果から上記2アンテナ方式では算出できないθpまたはθrのいずれかの姿勢角を求める。
【0068】
初期化時における整数バイアスの決定後にサイクルスリップが生じれば、図3および図4に示した場合と同様の処理によって整数バイアスの再決定を行う。但し、この実施形態ではアンテナが2本しかなく、整数バイアスの再決定の際に、上記2アンテナ方式では算出できないθpまたはθrのいずれかの姿勢角を求めるために、IMU/GPS統合化処理演算部は、加速度センサの出力により求めた姿勢角を利用する。図中のスイッチで示す記号は、この整数バイアスの再決定時にのみオンする、という概念を示している。
【0069】
尚、図5に示した例では、加速度センサの出力に基づいて姿勢角の1つを求めるようにしたが、その他に、傾斜計を用いてもよい。
【0070】
【発明の効果】
この発明によれば、IMUを用いて求めた物体の姿勢角を利用して、複数のアンテナの相対測位に用いるキャリア位相差の整数バイアスを直接的に決定するようにしたため、使用衛星の組み合わせが変更される場合やサイクルスリップの発生時にも、短時間のうちに整数バイアスを再決定することができる。その結果、ユーザーに対してほぼ連続して姿勢情報を出力できるようになる。また、整数バイアスの再決定時の正解率が高まるので、信頼性の高い姿勢情報を出力できるようになる。
【図面の簡単な説明】
【図1】ローカル座標系とアンテナ座標系におけるアンテナ配置の例を示す図
【図2】基線ベクトルと1重位相差等との関係を示す図
【図3】物体の姿勢角検出装置の構成を示すブロック図
【図4】図3におけるIMU/GPS統合化処理演算部等での処理手順を示すフローチャート
【図5】物体の姿勢角検出装置の他の構成を示すブロック図
Claims (5)
- 物体上の、一直線上に並ばないそれぞれ異なった位置に固定配置された少なくとも3つのアンテナと、これらのアンテナで複数の測位用衛星からの電波をそれぞれ受信し、これらのアンテナのうち1つの基準アンテナに対する他のアンテナの相対位置を、キャリア位相の1重位相差または2重位相差を観測するとともに、これらの位相差の整数バイアスを決定することによって求め、該アンテナの相対位置から前記物体の姿勢角を求める手段を備えた装置において、
物体に配置された慣性センサユニットと、該慣性センサユニットの出力により物体の姿勢角を求める手段と、
前記アンテナの相対位置から求められた前記物体の姿勢角で、前記慣性センサユニットの出力により求められる姿勢角を初期化する手段と、
前記慣性センサユニットの出力により求められた姿勢角から前記整数バイアスを直接再決定する手段とを設けた物体の姿勢検出装置。 - 物体上の、離れた位置に固定配置された2つのアンテナと、これらのアンテナで複数の測位用衛星からの電波をそれぞれ受信し、前記2つのアンテナのうち一方の基準アンテナに対する他方のアンテナの相対位置を、キャリア位相の1重位相差または2重位相差を観測するとともに、これらの位相差の整数バイアスを決定することによって求め、該アンテナの相対位置から前記物体の姿勢角を検出する装置において、
物体に配置された慣性センサユニットと、該慣性センサユニットの出力により物体の姿勢角を求める手段と、
前記アンテナの相対位置から求められた前記物体の姿勢角で、前記慣性センサユニットの出力により求められる姿勢角を初期化する手段と、
前記慣性センサユニットの出力により求められた姿勢角から前記整数バイアスを直接再決定する手段とを設けた物体の姿勢検出装置。 - 前記整数バイアスを再決定する手段は、
前記慣性センサユニットの出力により求められた姿勢角と、アンテナ座標系における基線ベクトルとからローカル座標系における基線ベクトルを求め、前記アンテナおよび測位用衛星の位置で決まるアンテナから測位用衛星への視野方向ベクトル、前記観測による1重位相差または2重位相差、および前記ローカル座標系における基線ベクトルに基づいて再決定する請求項1または2に記載の物体の姿勢検出装置。 - 物体上のそれぞれ異なった位置に固定配置された複数のアンテナと、これらのアンテナで複数の測位用衛星からの電波をそれぞれ受信し、前記複数のアンテナの相対位置を、前記アンテナを用いた前記測位用衛星の観測によるキャリア位相差の整数バイアスを決定することによって求め、該アンテナの相対位置から前記物体の姿勢角を検出する装置において、
物体に配置された慣性センサユニットと、該慣性センサユニットの出力により物体の姿勢角を求める手段と、
前記アンテナの相対位置から求められた前記物体の姿勢角で、前記慣性センサユニットの出力により求められる姿勢角を初期化する手段と、
前記慣性センサユニットの出力により求められた姿勢角から前記整数バイアスを直接再決定する手段とを設けた物体の姿勢検出装置。 - 物体上の、一直線上に並ばないそれぞれ異なった位置に固定配置された複数のアンテナで複数の測位用衛星からの電波をそれぞれ受信し、キャリア位相の1重位相差または2重位相差を観測し、該1重位相差または2重位相差の整数バイアスを初期化時において決定するステップと、
前記アンテナの相対位置から前記物体の姿勢角を求めるステップと、
前記アンテナの相対位置から求められた前記物体の姿勢角で、物体に配置された慣性センサユニットの出力により求められる姿勢角を初期化するステップと、
前記慣性センサユニットの出力により求められた姿勢角と、アンテナ座標系における基線ベクトルとからローカル座標系における基線ベクトルを求めるステップと、
前記アンテナおよび測位用衛星の位置で決まるアンテナから測位用衛星への視野方向ベクトル、前記観測による1重位相差または2重位相差、および前記ローカル座標系における基線ベクトルに基づいて、前記1重位相差または2重位相差の整数バイアスを直接再決定するステップとを備えた、物体の姿勢検出装置における整数バイアス再決定方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2001167407A JP4729197B2 (ja) | 2000-06-01 | 2001-06-01 | 物体の姿勢検出装置および整数バイアス再決定方法 |
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2000-164684 | 2000-06-01 | ||
JP2000164684 | 2000-06-01 | ||
JP2000164684 | 2000-06-01 | ||
JP2001167407A JP4729197B2 (ja) | 2000-06-01 | 2001-06-01 | 物体の姿勢検出装置および整数バイアス再決定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2002054946A JP2002054946A (ja) | 2002-02-20 |
JP4729197B2 true JP4729197B2 (ja) | 2011-07-20 |
Family
ID=26593168
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2001167407A Expired - Fee Related JP4729197B2 (ja) | 2000-06-01 | 2001-06-01 | 物体の姿勢検出装置および整数バイアス再決定方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4729197B2 (ja) |
Families Citing this family (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4445313B2 (ja) * | 2004-04-01 | 2010-04-07 | 古野電気株式会社 | 船舶航行支援装置 |
JP4461223B2 (ja) * | 2004-07-07 | 2010-05-12 | 独立行政法人情報通信研究機構 | 3次元配置型gps受信機、それを用いた姿勢測定方法及びシステム |
JP2006126181A (ja) * | 2004-10-01 | 2006-05-18 | Mitsubishi Electric Corp | 移動体姿勢検出装置 |
JP5180447B2 (ja) * | 2005-08-08 | 2013-04-10 | 古野電気株式会社 | キャリア位相相対測位装置及び方法 |
JP5301762B2 (ja) * | 2005-10-07 | 2013-09-25 | 古野電気株式会社 | キャリア位相相対測位装置 |
JP5083749B2 (ja) * | 2006-06-23 | 2012-11-28 | 独立行政法人港湾空港技術研究所 | Gps測位データの誤認識判別方法および誤認識判別装置 |
JP5052845B2 (ja) * | 2006-09-06 | 2012-10-17 | 日本無線株式会社 | 移動体姿勢計測装置 |
JP5084303B2 (ja) * | 2007-03-05 | 2012-11-28 | 日本無線株式会社 | 移動体姿勢計測装置 |
JP5430172B2 (ja) * | 2009-02-20 | 2014-02-26 | 三菱電機株式会社 | 方位算出装置、方位算出装置の方位算出方法および方位算出プログラム |
JP4854778B2 (ja) * | 2009-11-20 | 2012-01-18 | 住友ゴム工業株式会社 | 車両用推定航法装置、車両用推定航法及び車両用推定航法のプログラム |
KR101257935B1 (ko) | 2010-12-06 | 2013-04-23 | 국방과학연구소 | 관성 항법 시스템의 바이어스 추정치를 이용한 정렬 장치 및 그 항법 시스템 |
JP2012202749A (ja) * | 2011-03-24 | 2012-10-22 | Yokogawa Denshikiki Co Ltd | 方位測定装置 |
CN107110978B (zh) * | 2014-12-26 | 2021-02-02 | 古野电气株式会社 | 状态计算装置、状态计算方法及存储介质 |
US10514469B2 (en) | 2014-12-26 | 2019-12-24 | Furuno Electric Co., Ltd. | Attitude angle calculating device, method of calculating attitude angle, and attitude angle calculating program |
RU2608338C1 (ru) * | 2015-09-29 | 2017-01-17 | Акционерное общество "Федеральный научно-производственный центр "Нижегородский научно-исследовательский институт радиотехники" | Устройство обработки сигналов в наземно-космической просветной радиолокационной системе |
JP6886800B2 (ja) * | 2016-11-01 | 2021-06-16 | 三菱重工業株式会社 | 移動体観測システム、移動体観測方法及びプログラム |
US11187812B2 (en) * | 2017-06-29 | 2021-11-30 | Panasonic Intellectual Property Management Co., Ltd. | Positioning method and positioning terminal |
JP7037529B2 (ja) * | 2019-09-25 | 2022-03-16 | 日立建機株式会社 | 作業機械 |
CN115516338A (zh) * | 2020-05-15 | 2022-12-23 | 索尼集团公司 | 信息处理装置及信息处理方法 |
CN117405108B (zh) * | 2023-11-08 | 2024-05-07 | 中国人民解放军63620部队 | 目标姿态测量方法、系统、电子设备和介质 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0688865A (ja) * | 1990-08-20 | 1994-03-29 | Texas Instr Inc <Ti> | 単一受信機使用gps指向又は姿勢システム |
JPH08327374A (ja) * | 1995-05-30 | 1996-12-13 | Koden Electron Co Ltd | Gps受信機 |
JPH09119972A (ja) * | 1995-10-26 | 1997-05-06 | Furuno Electric Co Ltd | 相対測位装置および相対測位方法 |
JPH09154308A (ja) * | 1995-12-12 | 1997-06-17 | Kubota Corp | 作業機械の姿勢制御装置 |
JPH10170629A (ja) * | 1996-12-13 | 1998-06-26 | Taiyo Musen Kk | Gps信号を用いた方位・姿勢測定法 |
JPH10221109A (ja) * | 1997-02-01 | 1998-08-21 | Litton Syst Inc | プラットホームのロール、ピッチおよび機首方位を判断するカルマンフィルタプロセスへの入力のための可観測値を得るための方法および装置 |
JP2001059863A (ja) * | 1999-08-25 | 2001-03-06 | Furuno Electric Co Ltd | 方位測定装置 |
JP2001099910A (ja) * | 1999-10-04 | 2001-04-13 | Japan Aviation Electronics Industry Ltd | Rtk・gps−imu併用測位方法 |
-
2001
- 2001-06-01 JP JP2001167407A patent/JP4729197B2/ja not_active Expired - Fee Related
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0688865A (ja) * | 1990-08-20 | 1994-03-29 | Texas Instr Inc <Ti> | 単一受信機使用gps指向又は姿勢システム |
JPH08327374A (ja) * | 1995-05-30 | 1996-12-13 | Koden Electron Co Ltd | Gps受信機 |
JPH09119972A (ja) * | 1995-10-26 | 1997-05-06 | Furuno Electric Co Ltd | 相対測位装置および相対測位方法 |
JPH09154308A (ja) * | 1995-12-12 | 1997-06-17 | Kubota Corp | 作業機械の姿勢制御装置 |
JPH10170629A (ja) * | 1996-12-13 | 1998-06-26 | Taiyo Musen Kk | Gps信号を用いた方位・姿勢測定法 |
JPH10221109A (ja) * | 1997-02-01 | 1998-08-21 | Litton Syst Inc | プラットホームのロール、ピッチおよび機首方位を判断するカルマンフィルタプロセスへの入力のための可観測値を得るための方法および装置 |
JP2001059863A (ja) * | 1999-08-25 | 2001-03-06 | Furuno Electric Co Ltd | 方位測定装置 |
JP2001099910A (ja) * | 1999-10-04 | 2001-04-13 | Japan Aviation Electronics Industry Ltd | Rtk・gps−imu併用測位方法 |
Also Published As
Publication number | Publication date |
---|---|
JP2002054946A (ja) | 2002-02-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4729197B2 (ja) | 物体の姿勢検出装置および整数バイアス再決定方法 | |
US6424915B1 (en) | System for determining the heading and/or attitude of a body | |
US6792380B2 (en) | Attitude angle detecting apparatus | |
US6611228B2 (en) | Carrier phase-based relative positioning apparatus | |
JP5301762B2 (ja) | キャリア位相相対測位装置 | |
US6496778B1 (en) | Real-time integrated vehicle positioning method and system with differential GPS | |
US7355549B2 (en) | Apparatus and method for carrier phase-based relative positioning | |
US20120086598A1 (en) | Apparatus and methods for driftless attitude determination and reliable localization of vehicles | |
EP0856747A1 (en) | Method and apparatus for attitude determination utilizing an inertial measurement unit and a plurality of satellite transmitters | |
EP2488892B1 (en) | Ultra-short baseline gnss receiver | |
KR100443550B1 (ko) | 오차보정시스템을 구비하는 관성측정유닛-지피에스통합시스템과 미지정수 검색범위 축소방법 및 사이클 슬립검출방법, 및 그를 이용한 항체 위치, 속도,자세측정방법 | |
JP2011185899A (ja) | 位置標定装置、位置標定装置の位置標定方法および位置標定プログラム | |
JP4563157B2 (ja) | 物体の方位および姿勢検出装置 | |
JP2008039691A (ja) | 搬送波位相式測位装置 | |
JP2012098185A (ja) | 方位角推定装置及びプログラム | |
CA2716597A1 (en) | Apparatus and method for driftless attitude determination and reliable localization of vehicles | |
JP5180447B2 (ja) | キャリア位相相対測位装置及び方法 | |
JP2012202749A (ja) | 方位測定装置 | |
JPH07209406A (ja) | 地球上の特定の位置を決定するための方法および装置 | |
CN108351420B (zh) | 用于在惯性测量单元的静态对准期间检测寄生移动的方法、以及相关的检测设备 | |
Kis et al. | Development of state estimation system with INS, magnetometer and carrier phase GPS for vehicle navigation | |
JP7140443B2 (ja) | アンテナ間相対位置推定方法及びアンテナ間相対位置推定プログラム | |
KR100341801B1 (ko) | 다중안테나를이용한도시형차량항법시스템 | |
US9976857B2 (en) | Apparatus for providing an azimuth angle | |
Wang et al. | Development of a low-cost solution for GPS/gyro attitude determination |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20080521 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20100831 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20100921 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20101117 |
|
RD02 | Notification of acceptance of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7422 Effective date: 20101117 |
|
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: 20110412 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20110418 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 4729197 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140422 Year of fee payment: 3 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R255 | Notification that request for automated payment was rejected |
Free format text: JAPANESE INTERMEDIATE CODE: R2525 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
LAPS | Cancellation because of no payment of annual fees |