JP3711538B2 - Inertial altitude measurement method - Google Patents

Inertial altitude measurement method Download PDF

Info

Publication number
JP3711538B2
JP3711538B2 JP2002188623A JP2002188623A JP3711538B2 JP 3711538 B2 JP3711538 B2 JP 3711538B2 JP 2002188623 A JP2002188623 A JP 2002188623A JP 2002188623 A JP2002188623 A JP 2002188623A JP 3711538 B2 JP3711538 B2 JP 3711538B2
Authority
JP
Japan
Prior art keywords
altitude
gps
inertial
barometric
output
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
Application number
JP2002188623A
Other languages
Japanese (ja)
Other versions
JP2004028911A (en
Inventor
雅喜 山田
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.)
Japan Aviation Electronics Industry Ltd
Original Assignee
Japan Aviation Electronics Industry Ltd
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 Japan Aviation Electronics Industry Ltd filed Critical Japan Aviation Electronics Industry Ltd
Priority to JP2002188623A priority Critical patent/JP3711538B2/en
Publication of JP2004028911A publication Critical patent/JP2004028911A/en
Application granted granted Critical
Publication of JP3711538B2 publication Critical patent/JP3711538B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Landscapes

  • Traffic Control Systems (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Description

【0001】
【発明の属する技術分野】
本発明は、鉛直方向加速度の積分処理による高度を気圧高度で補正する高度測定方式に関し、特に気圧高度の取得ができなくなった時における測定値の階段的な変動を抑制し、連続的な高度測定を可能にする慣性高度測定方式に関する。
【0002】
【従来の技術】
航空機等の飛翔体の慣性高度を測定する手段、またはその測定高度を補正する従来の手段としては、例えば特開平11ー211472号公報に開示された上下動変位測定装置がある。この上下動変位測定装置は、上下加速度測定手段、上下動変位演算手段および高度測定手段を備えており、上下動変位演算手段において、高度測定手段による測定高度を基準高度として、上下加速度測定手段の加速度出力を負帰還ループを介して積分し、その積分出力を上下動変位測定データとしている。
【0003】
図2は、その上下動変位測定装置における上下動変位演算手段の構成を示すブロック図である。図2の上下動変位演算手段は、減算回路23、25及び31と、積分回路24、26及び29と、乗算回路27および28と、加算回路30とで構成されている。鉛直方向加速度117は前記上下加速度測定手段の出力であり、気圧高度120は前記高度測定手段の出力である。
【0004】
その鉛直方向加速度117は、減算回路23に入力され、気圧高度120が減算回路31に入力される。減算回路23の出力は積分回路24で積分される。積分回路24の出力は、減算回路25を経由して、鉛直方向速度118として積分回路26に入力され、積分される。積分回路26の積分出力は、上下動変位119として出力されるとともに、減算回路31に入力される。減算回路31は、上下動変位119と気圧高度120との差分121を生成する。差分121は、乗算回路28、積分回路29および加算回路30を介して、帰還信号124として減算回路23に対し負帰還されるとともに、乗算回路27を介して、帰還信号125として減算回路25に負帰還される。
【0005】
乗算回路28、積分回路29および加算回路30並びに乗算回路27を含む帰還回路は、上下動変位119を気圧高度120に移行させるように作用している。この帰還回路では、負帰還ループが構成されており、発振を生じさせない値に帰還定数が設定してある。この負帰還ループの安定時においては、減算回路31における差分121はゼロに移行する。差分121がゼロに移行することにより、上下動変位119では、鉛直方向加速度117における加速度計のバイアス誤差の影響が除かれる。かくして、図2の上下動変位演算手段において、気圧高度120を基準高度とする上下動変位119が高い精度で測定される。
【0006】
【発明が解決しようとする課題】
特開平11ー211472号公報の上下動変位測定装置においては、基準高度として気圧高度を用いる場合に、気圧高度計が無効になると、前述の負帰還ループを構成することが不可能となり、鉛直方向加速度117の積分演算の過程において、鉛直方向加速度計のバイアス誤差が時間の経過とともに累積され、上下動変位119の値が発散する状態となり、該上下動変位の測定を的確に行うことが不可能になる。また、気圧高度に代えて、基準高度としてGPS高度を用いる場合には、該GPS高度の精度は気圧高度に比較して誤差が大きく、安定した精度の高い高度測定を行うことができないという欠点とともに、上記のように気圧高度を基準高度として用いる場合と同様に、GPS受信機が無効となった場合には、測定高度が発散して測定不可能になるという欠点がある。
【0007】
本発明の目的は、上記の欠点を排除し、気圧高度の有効度とGPS高度の有効度とを比較照合して、基準高度として気圧高度を優先し、より有効な高度の方を選択して、常に安定に慣性高度を測定することを可能にするとともに、基準高度として気圧高度が選択されている動作時に、気圧高度の有効度が低下して基準高度をGPS高度に切替えるような事態においても、急激な基準高度の変動を抑制して緩やかにGPS高度に移行させることにより、高度測定を安定且つ円滑に行うことのできる慣性高度測定方式を提供することにある。
【0008】
【課題を解決するための手段】
前述の課題を解決するために本発明は次の手段を提供する。
【0009】
(1)加速度計出力の垂直方向加速度および基準高度を処理する演算処理手段および該基準高度として気圧高度を出力する気圧高度計とを備え、該垂直方向加速度を該基準高度にダンピングさせることにより慣性高度を生成するフィードバックループを該演算処理手段に含み、該気圧高度が有効なときにだけに該基準高度として該気圧高度を選択する慣性高度測定方式において、
GPS受信機からGPS高度を受け、
前記気圧高度が有効な状態では、該気圧高度を前記基準高度として前記演算処理手段に供給して得た前記慣性高度を基準として前記GPS高度のバイアス誤差を推定し、
前記気圧高度が無効な状態になったとき、該気圧高度が無効な状態になる直前の前記バイアス誤差により前記GPS高度を補正したGPS補正高度を生成し、該GPS補正高度を前記基準高度として選択する
ことを特徴とする慣性高度測定方式。
【0010】
(2)加速度計から出力される垂直方向加速度と基準高度とを処理することにより慣性高度を測定するために、第1、第2及び第3の減算手段、第1及び第2の積分手段並びに第1及び第2の負帰還手段を備え、該第1の負帰還手段は該第3の減算手段の出力に基づき第1の帰還値を生成し、該第2の負帰還手段は該第3の減算手段の出力に基づき第2の帰還値を生成し、該第1の減算手段は該垂直方向加速度と該第1の帰還値との減算をし、該第1の積分手段は該第1の減算手段の出力に関する時間積分をし、該第2の減算手段は該第1の積分手段の出力と該第2の帰還値との減算をし、該第2の積分手段は該第2の減算手段の出力に関する時間積分値を該慣性高度として出力し、該第3の減算手段は該慣性高度と該基準高度との減算をする慣性高度測定方式において、
気圧高度を出力する気圧高度計と、GPS高度を出力するGPS受信機と、GPS高度バイアス推定手段と、基準高度選択手段とを備え、
前記GPS高度バイアス推定手段は、前記慣性高度および前記GPS高度を入力し、該慣性高度に対する該GPS高度のバイアス誤差を推定し、該GPS高度を該バイアス誤差分だけ補正したGPS補正高度を生成し、
前記基準高度選択手段は、前記気圧高度が有効であるときは前記基準高度として該気圧高度を選択し、該気圧高度が無効であるときは該基準高度として前記GPS補正高度を選択する
ことを特徴とする慣性高度測定方式。
【0011】
(3)前記GPS高度バイアス推定手段は、第4の減算手段と、第2の切替手段と、第3の積分手段と、第5の減算手段とを有してなり、
前記第4の減算手段は前記慣性高度と前記GPS補正高度との減算をし、
前記第2の切替手段は、前記基準高度選択手段が前記気圧高度を選択するときは前記第4の減算手段の出力をそのまま出力端に接続し、前記基準高度選択手段が前記GPS補正高度を選択するときは出力端を開放端とし、
前記第3の積分手段は、前記第2の切替手段の出力端の信号に関する時間積分をし、積分値を前記GPS高度のバイアス誤差として出力し、
前記第5の減算手段は、前記GPS高度と該バイアス誤差との減算をし、減算値を前記GPS補正高度として出力する
ことを特徴とする前記(2)に記載の慣性高度測定方式。
【0012】
【発明の実施の形態】
次に、本発明の慣性高度測定方式について実施の形態を挙げ、一層具体的に説明する。本実施形態は、基準高度として、気圧高度とGPS高度の有効度を比較照合して、有効度の高い方の測定高度を選択し、負帰還ループを含む2重積分機構による鉛直方向加速度の積分処理を介して、基準高度の値に追随する状態で慣性高度を測定することを特徴としており、基準高度としての気圧高度とGPS高度との相互補完作用により、高度の測定が安定に行われる。また基準高度を気圧高度からGPS高度に切替える際においては、基準高度の急激な変動が回避されて円滑に高度測定が行われる。
【0013】
図1は、本発明の一実施形態の構成を示すブロック図である。図1に示されるように、本実施の形態は、主要構成要素として、慣性高度補正部1と、GPS高度バイアス推定部2と、基準高度選択部3と、慣性加速度計4と、気圧高度計5と、GPS受信機6とを備えて構成される。慣性高度補正部1は、減算回路8,10,12と、積分回路9,11と、乗算回路13,14と、切替回路15と、積分回路16と加算回路17とを備えてなる。積分回路9及び11は2重積分機構をなしている。乗算回路14、切替回路15、積分回路16および加算回路17は第1の帰還回路をなし、乗算回路13は第2の帰還回路をなし、第1の帰還回路および第2の帰還回路でもって負帰還ループを形成している。
【0014】
GPS高度バイアス推定部2は、減算回路18,21と、切替回路19と、積分回路20とで構成される。このGPS高度バイアス推定部2は、負帰還ループとして形成されており、基準高度114として気圧高度104が選択されている間にGPS高度誤差を推定し、積分回路20に保持し、GPS高度誤差111として積分回路20から出力する。基準高度選択部3は、気圧高度計5とGPS受信機6による各測定高度の有効度を判定する高度計判定回路7と、基準高度を選択して出力する切替回路22とを備えて構成される。
【0015】
次に、基準高度114として気圧高度104が選択されているときに、基準高度を気圧高度104からGPS高度112へ切替える際の動作について説明する。
【0016】
図1において、気圧高度計5から出力される気圧高度104と、GPS受信機6から出力されるGPS高度112は、高度計判定回路7に入力される。高度計判定回路7は、気圧高度104とGPS高度112とのうちのどちらを基準高度114として選択するべきかを指定する基準高度選択信号116を生成する。基準高度選択信号116は切替回路15、19及び22に入力される。基準高度選択信号116が論理値“1”とき、切替回路15、19及び22の各可動切片は接点aに接続され、基準高度選択信号116が論理値“0”とき、切替回路15、19及び22の各可動切片は接点bに接続される。
【0017】
本実施の形態に電源が投入された当初(起動時)においては、高度計判定回路7は、気圧高度104の精度がGPS高度112の精度より高いものと想定し、気圧高度104の有効度の方がGPS高度112の有効度より高く、気圧高度104を基準高度114として選択されるべきであると判定する。そこで、起動時には、高度計判定回路7は、気圧高度104を基準高度114として選択させる論理値“1”に基準高度選択信号116を設定する。起動時の基準高度選択信号116を受けた切替回路15、19及び22においては、図1に示されるように、可動切片が全て接点aの側に接続される。このとき、基準高度選択部3においては切替回路22が基準高度114として気圧高度104を選択し、GPS高度バイアス推定部2においては、GPS高度112の高度誤差を推定して保持する負帰還ループが形成される。また慣性高度補正部1においては、切替回路15が減算回路12の出力の差分105を積分回路16へ導く。このとき、前記第1の帰還回路は、乗算回路14、積分回路16および加算回路17により構成される。
【0018】
慣性加速度計4により検出された鉛直方向加速度101は、慣性高度補正部1の減算回路8に入力される。慣性高度補正部1は、鉛直方向加速度101に対応する積分手段として積分回路9及び11を含む。終段の積分回路11から出力される慣性高度103は、減算回路12に入力されるとともに、GPS高度バイアス推定部2の減算回路18に入力される。減算回路12は、慣性高度103と基準高度104との間の差分出力105を生成する。差分出力105は、第1および第2の帰還回路でなる前記負帰還ループに入力される。乗算回路14、積分回路16および加算回路17でなる第1の帰還回路による帰還信号108は減算回路8に負帰還される。乗算回路13でなる第2の帰還回路による帰還信号109は減算回路10に負帰還入力される。これらの第1および第2の帰還回路では、当該負帰還ループが発振しないように帰還利得が設定してある。この負帰還ループにおける減算回路10からは慣性速度102が出力される。また、その負帰還ループにおける積分回路11からは、本実施の形態の測定対象である慣性高度103が、基準高度の気圧高度104に追随する状態で常時安定に出力される。
【0019】
GPS高度バイアス推定部2においては、慣性高度103が減算回路18に入力されるとともに、GPS高度112が減算回路21に入力される。このとき、減算回路18、切替回路19、積分回路20および減算回路21が負帰還ループを構成しているので、積分回路20からは、慣性高度103に対するGPS高度112の推定誤差としてGPS高度誤差111が出力されるとともに、減算回路21からは、慣性高度103に対して推定されるGPS補正高度113が出力される。GPS高度誤差111は、慣性高度103を基準にして推定したGPS高度112のバイアス誤差であり、時間の経過やGPS受信機6の緯度経度の変動があってもさして変動しない。GPS補正高度113は、GPS高度112をGPS高度誤差111で補正した高度である。
【0020】
GPS高度バイアス推定部2において負帰還ループが構成され、GPS高度誤差111が安定に生成されている状態で、もし基準高度選択信号116が論理値“1”から論理値“0”に変化し、切替回路19の可動切片が接点aから接点bへ切り換わったならば、積分回路20の入力110は開放となるので、積分回路20はその時のデータ(GPS高度誤差111)を保持する。
【0021】
このように、基準高度として気圧高度104が選択されている状態においては、慣性高度補正部1からは、慣性高度103が、気圧高度104に追随する状態で安定に出力され、またGPS高度バイアス推定部2においては、推定誤差としてのGPS高度誤差111と、推定されるGPS補正高度113が生成される。
【0022】
ここまで、気圧高度104がGPS高度112より有効であると高度計判定回路7が判定した状態の動作について説明した。この状態において、例えば、気圧高度104がステップ的に変動し、気圧高度104の信頼性が疑われる等の事情が発生したり、或いは気圧高度測定限界にまで高度が上昇し、気圧高度104よりGPS高度112の方が有効であると高度計判定回路7が判定したとする。このとき、基準高度選択信号116の論理値は“0”となり、切替回路15、19及び22における可動切片は全て接点bの側に切り換わる。この切り換わりにより、GPS高度バイアス推定部2においては前記負帰還ループが解消され、切り換わりの時点におけるGPS高度誤差111は保持される。減算回路21は、GPS高度112をGPS高度誤差111で補正したGPS補正高度113を生成し続ける。
【0023】
このときのGPS高度誤差111は、気圧高度104を基準高度114として選択したときにおける慣性高度103に基づき生成したGPS高度112のバイアス誤差であるから、そのGPS高度誤差111によりGPS高度112を補正して得たGPS補正高度113は切替時の気圧高度104に連続しており、切替後のGPS補正高度113はGPS高度112に追随する。したがって、基準高度114の値は、気圧高度104からGPS補正高度113に切り換わるとき、連続し、階段的に変動することはない。以後、慣性高度補正部1は、GPS補正高度113を基準高度104として補正した慣性高度103を出力する。勿論、切替において基準高度104が連続するから、慣性高度103も切替において連続し、階段状に不連続に変動することはない。
GPS補正高度113は、切替回路22を経由して、基準高度114として慣性高度補正部1の減算回路12に入力される。慣性高度補正部1においては、切替回路15の可動切片が接点bに切り換わるので、積分回路16が第1の帰還回路から除去される。そこで、積分回路11の出力とGPS補正高度113の差分105は、第1の帰還回路では、乗算回路14で所定の利得を乗じられ、乗算値107となる。乗算値107は、加算回路17を経由して、そのまま帰還信号108として減算回路8に対し負帰還入力される。乗算回路13の帰還信号109は、基準高度114の切替前と同じく、減算回路10に負帰還入力される。積分回路11からは、GPS補正高度113を基準高度114とし、積分回路16を除いた新たな負帰還ループにより生成した慣性高度103が出力される。
【0024】
【発明の効果】
以上に実施の形態を挙げ具体的に説明したように、本発明によれば、気圧高度又はGPS補正高度の一方を基準高度として選択できるようにし、基準高度として気圧高度を選択して慣性高度を測定している間に、その慣性高度に対して推定されるGPS高度のバイアス誤差(GPS高度誤差111)を継続的に生成し、保持し、かつGPS高度を該バイアス誤差により補正したGPS補正高度(113)を常時生成することにより、基準高度を気圧高度からGPS補正高度に切替える事態においても、該基準高度の値が連続し、階段的に急激な変動することなく、ひいては慣性高度の測定を安定且つ円滑に行うことができる慣性高度測定方式を得ることができる。
【図面の簡単な説明】
【図1】本発明の一実施の形態を示すブロック図である。
【図2】従来例を示すブロック図である。
【符号の説明】
1 慣性高度補正部
2 GPS高度バイアス推定部 3 基準高度選択部
4 慣性加速度計
5 気圧高度計
6 GPS受信機
7 高度計判定回路
8,10,12,18,21,23,25,31 減算回路
9,11,16,20,24,26,29 積分回路
13,14,27,28 乗算回路
15,19,22 切替回路
17,30 加算回路
101,117 鉛直方向加速度
102,118 慣性速度
103,119 慣性高度
104 気圧高度
105 慣性高度103と基準高度104との間の差分出力
106 積分回路16出力
107 乗算回路14の乗算値
108 帰還信号(加算回路17の加算値)
109 乗算回路13の乗算値
110 切替回路19出力
111 GPS高度誤差
112 GPS高度
113 GPS補正高度
114 基準高度
113 GPS補正高度
116 基準高度選択信号
119 上下動変位
120 基準高度
121 慣性高度119と基準高度120との間の差分出力
124 帰還信号
125 帰還信号(乗算回路27出力)
[0001]
BACKGROUND OF THE INVENTION
The present invention relates to an altitude measurement method that corrects altitude by integration of vertical acceleration with barometric altitude, and in particular, suppresses stepwise fluctuations in measured values when the barometric altitude cannot be acquired, thereby continuously measuring altitude. The present invention relates to an inertial altitude measurement method that enables
[0002]
[Prior art]
As a means for measuring the inertial altitude of a flying object such as an aircraft or a conventional means for correcting the measured altitude, there is a vertical movement displacement measuring apparatus disclosed in, for example, Japanese Patent Application Laid-Open No. 11-21472. This vertical movement displacement measuring device includes vertical acceleration measuring means, vertical movement displacement calculating means, and altitude measuring means. In the vertical movement displacement calculating means, the vertical height measuring means uses the height measured by the height measuring means as a reference height. The acceleration output is integrated via a negative feedback loop, and the integrated output is used as vertical movement displacement measurement data.
[0003]
FIG. 2 is a block diagram showing the configuration of the vertical movement displacement calculating means in the vertical movement displacement measuring apparatus. 2 includes subtracting circuits 23, 25, and 31, integrating circuits 24, 26, and 29, multiplying circuits 27 and 28, and an adding circuit 30. A vertical acceleration 117 is an output of the vertical acceleration measuring means, and a barometric altitude 120 is an output of the altitude measuring means.
[0004]
The vertical acceleration 117 is input to the subtraction circuit 23, and the atmospheric pressure altitude 120 is input to the subtraction circuit 31. The output of the subtracting circuit 23 is integrated by the integrating circuit 24. The output of the integrating circuit 24 is input to the integrating circuit 26 via the subtracting circuit 25 as the vertical speed 118 and integrated. The integration output of the integration circuit 26 is output as the vertical movement displacement 119 and also input to the subtraction circuit 31. The subtraction circuit 31 generates a difference 121 between the vertical movement displacement 119 and the atmospheric pressure altitude 120. The difference 121 is negatively fed back as a feedback signal 124 to the subtraction circuit 23 via the multiplication circuit 28, the integration circuit 29 and the addition circuit 30, and is negatively fed to the subtraction circuit 25 as a feedback signal 125 via the multiplication circuit 27. Returned.
[0005]
The feedback circuit including the multiplier circuit 28, the integrator circuit 29, the adder circuit 30, and the multiplier circuit 27 operates to shift the vertical movement displacement 119 to the atmospheric pressure altitude 120. In this feedback circuit, a negative feedback loop is formed, and the feedback constant is set to a value that does not cause oscillation. When the negative feedback loop is stable, the difference 121 in the subtraction circuit 31 shifts to zero. By shifting the difference 121 to zero, the vertical movement displacement 119 eliminates the influence of the accelerometer bias error in the vertical acceleration 117. Thus, in the vertical movement displacement calculation means of FIG. 2, the vertical movement displacement 119 with the atmospheric pressure altitude 120 as the reference altitude is measured with high accuracy.
[0006]
[Problems to be solved by the invention]
In the vertical displacement measuring apparatus disclosed in Japanese Patent Application Laid-Open No. 11-211472, when the barometric altitude is used as the reference altitude, if the barometric altimeter becomes invalid, the negative feedback loop cannot be formed, and the vertical acceleration In the process of integral calculation 117, the bias error of the vertical accelerometer is accumulated over time, and the value of vertical displacement 119 diverges, making it impossible to accurately measure the vertical displacement. Become. In addition, when the GPS altitude is used as the reference altitude instead of the barometric altitude, the accuracy of the GPS altitude is large compared to the barometric altitude, with the disadvantage that stable and highly accurate altitude measurement cannot be performed. Similarly to the case where the barometric altitude is used as the reference altitude as described above, there is a drawback that when the GPS receiver becomes invalid, the measurement altitude diverges and measurement is impossible.
[0007]
The object of the present invention is to eliminate the above-mentioned drawbacks, compare and collate the effectiveness of the barometric altitude with the effectiveness of the GPS altitude, give priority to the barometric altitude as the reference altitude, and select the more effective altitude. In addition to making it possible to measure the inertial altitude stably at all times, even when the atmospheric altitude is selected as the reference altitude, the effectiveness of the atmospheric altitude decreases and the reference altitude is switched to the GPS altitude. Another object of the present invention is to provide an inertial altitude measurement method capable of performing altitude measurement stably and smoothly by suppressing a sudden change in reference altitude and gradually shifting to a GPS altitude.
[0008]
[Means for Solving the Problems]
In order to solve the above-mentioned problems, the present invention provides the following means.
[0009]
(1) Comprising arithmetic processing means for processing the vertical acceleration and reference altitude of the accelerometer output, and a barometric altimeter that outputs the barometric altitude as the reference altitude, and by damping the vertical acceleration to the reference altitude, In an inertial altitude measurement method that includes a feedback loop that generates a pressure altitude and that selects the barometric altitude as the reference altitude only when the barometric altitude is valid,
Receive the GPS altitude from the GPS receiver,
In a state where the atmospheric pressure altitude is valid, a bias error of the GPS altitude is estimated based on the inertial altitude obtained by supplying the atmospheric altitude to the arithmetic processing unit as the reference altitude,
When the barometric altitude becomes invalid, a GPS corrected altitude is generated by correcting the GPS altitude by the bias error immediately before the barometric altitude becomes invalid, and the GPS corrected altitude is selected as the reference altitude. An inertial altitude measurement system characterized by
[0010]
(2) first, second and third subtracting means, first and second integrating means for measuring the inertial height by processing the vertical acceleration output from the accelerometer and the reference height; and First and second negative feedback means are provided, the first negative feedback means generates a first feedback value based on the output of the third subtraction means, and the second negative feedback means is the third negative feedback means. A second feedback value is generated based on the output of the subtracting means, the first subtracting means subtracts the vertical acceleration and the first feedback value, and the first integrating means is the first integrating means. The second subtracting means subtracts the output of the first integrating means and the second feedback value, and the second integrating means A time integral value relating to the output of the subtracting means is output as the inertia height, and the third subtracting means subtracts the inertia height from the reference height. In the inertial altitude measurement method,
A barometric altimeter that outputs a barometric altitude; a GPS receiver that outputs a GPS altitude; a GPS altitude bias estimating means; and a reference altitude selecting means.
The GPS altitude bias estimation means inputs the inertial altitude and the GPS altitude, estimates a bias error of the GPS altitude relative to the inertial altitude, and generates a GPS corrected altitude in which the GPS altitude is corrected by the bias error. ,
The reference altitude selecting means selects the barometric altitude as the reference altitude when the barometric altitude is valid, and selects the GPS corrected altitude as the reference altitude when the barometric altitude is invalid. Inertial altitude measurement method.
[0011]
(3) The GPS altitude bias estimation means includes a fourth subtraction means, a second switching means, a third integration means, and a fifth subtraction means.
The fourth subtracting means subtracts the inertial altitude and the GPS corrected altitude,
The second switching means connects the output of the fourth subtracting means to the output terminal as it is when the reference altitude selecting means selects the barometric altitude, and the reference altitude selecting means selects the GPS corrected altitude. The output end is the open end,
The third integrating means performs time integration with respect to the signal at the output terminal of the second switching means, and outputs an integrated value as a bias error of the GPS altitude.
The inertial altitude measurement method according to (2), wherein the fifth subtracting unit subtracts the GPS altitude from the bias error and outputs a subtraction value as the GPS corrected altitude.
[0012]
DETAILED DESCRIPTION OF THE INVENTION
Next, the inertial altitude measurement method of the present invention will be described more specifically with reference to an embodiment. In the present embodiment, as the reference altitude, the effectiveness of the barometric altitude and the GPS altitude are compared and selected, the measurement altitude having the higher effectiveness is selected, and the vertical acceleration is integrated by the double integration mechanism including the negative feedback loop. Through the processing, the inertial altitude is measured in a state of following the value of the reference altitude, and the altitude measurement is stably performed by the mutual complementary action of the barometric altitude as the reference altitude and the GPS altitude. Further, when the reference altitude is switched from the barometric altitude to the GPS altitude, the altitude measurement is smoothly performed while avoiding a sudden change in the reference altitude.
[0013]
FIG. 1 is a block diagram showing a configuration of an embodiment of the present invention. As shown in FIG. 1, the present embodiment includes, as main components, an inertial altitude correction unit 1, a GPS altitude bias estimation unit 2, a reference altitude selection unit 3, an inertial accelerometer 4, and a barometric altimeter 5. And a GPS receiver 6. The inertial altitude correction unit 1 includes subtracting circuits 8, 10 and 12, integrating circuits 9 and 11, multiplying circuits 13 and 14, a switching circuit 15, an integrating circuit 16 and an adding circuit 17. The integration circuits 9 and 11 form a double integration mechanism. The multiplication circuit 14, the switching circuit 15, the integration circuit 16 and the addition circuit 17 constitute a first feedback circuit, the multiplication circuit 13 constitutes a second feedback circuit, and the first feedback circuit and the second feedback circuit are negative. A feedback loop is formed.
[0014]
The GPS altitude bias estimation unit 2 includes subtraction circuits 18 and 21, a switching circuit 19, and an integration circuit 20. The GPS altitude bias estimation unit 2 is formed as a negative feedback loop, estimates a GPS altitude error while the atmospheric pressure altitude 104 is selected as the reference altitude 114, holds the GPS altitude error in the integrating circuit 20, and stores the GPS altitude error 111. Is output from the integrating circuit 20 as follows. The reference altitude selection unit 3 includes an altimeter determination circuit 7 that determines the effectiveness of each altitude measured by the barometric altimeter 5 and the GPS receiver 6, and a switching circuit 22 that selects and outputs the reference altitude.
[0015]
Next, the operation when the reference altitude 104 is switched from the atmospheric pressure altitude 104 to the GPS altitude 112 when the atmospheric altitude 104 is selected as the reference altitude 114 will be described.
[0016]
In FIG. 1, the barometric altitude 104 output from the barometric altimeter 5 and the GPS altitude 112 output from the GPS receiver 6 are input to the altimeter determination circuit 7. The altimeter determination circuit 7 generates a reference altitude selection signal 116 that specifies which one of the barometric altitude 104 and the GPS altitude 112 should be selected as the reference altitude 114. The reference altitude selection signal 116 is input to the switching circuits 15, 19 and 22. When the reference altitude selection signal 116 is a logical value “1”, each movable segment of the switching circuits 15, 19 and 22 is connected to the contact a, and when the reference altitude selection signal 116 is a logical value “0”, the switching circuits 15, 19 and Each movable piece 22 is connected to a contact b.
[0017]
At the beginning (at start-up) when the power is turned on in this embodiment, the altimeter determination circuit 7 assumes that the accuracy of the barometric altitude 104 is higher than the accuracy of the GPS altitude 112, and the effectiveness of the barometric altitude 104 Is higher than the effectiveness of the GPS altitude 112, and the barometric altitude 104 should be selected as the reference altitude 114. Therefore, at the time of activation, the altimeter determination circuit 7 sets the reference altitude selection signal 116 to a logical value “1” that causes the barometric altitude 104 to be selected as the reference altitude 114. In the switching circuits 15, 19, and 22 that have received the reference altitude selection signal 116 at the time of activation, as shown in FIG. 1, all the movable segments are connected to the contact a side. At this time, in the reference altitude selection unit 3, the switching circuit 22 selects the barometric altitude 104 as the reference altitude 114, and in the GPS altitude bias estimation unit 2, there is a negative feedback loop that estimates and holds the altitude error of the GPS altitude 112. It is formed. In the inertial altitude correction unit 1, the switching circuit 15 guides the difference 105 of the output of the subtraction circuit 12 to the integration circuit 16. At this time, the first feedback circuit includes a multiplication circuit 14, an integration circuit 16, and an addition circuit 17.
[0018]
The vertical acceleration 101 detected by the inertial accelerometer 4 is input to the subtraction circuit 8 of the inertial altitude correction unit 1. The inertial altitude correction unit 1 includes integration circuits 9 and 11 as integration means corresponding to the vertical acceleration 101. The inertia altitude 103 output from the final stage integration circuit 11 is input to the subtraction circuit 12 and also input to the subtraction circuit 18 of the GPS altitude bias estimation unit 2. The subtraction circuit 12 generates a difference output 105 between the inertial altitude 103 and the reference altitude 104. The differential output 105 is input to the negative feedback loop composed of the first and second feedback circuits. The feedback signal 108 from the first feedback circuit including the multiplication circuit 14, the integration circuit 16 and the addition circuit 17 is negatively fed back to the subtraction circuit 8. The feedback signal 109 from the second feedback circuit composed of the multiplication circuit 13 is input to the subtraction circuit 10 as negative feedback. In these first and second feedback circuits, the feedback gain is set so that the negative feedback loop does not oscillate. The inertia speed 102 is output from the subtraction circuit 10 in the negative feedback loop. In addition, from the integration circuit 11 in the negative feedback loop, the inertial altitude 103 which is the measurement target of the present embodiment is always stably output in a state of following the atmospheric altitude 104 of the reference altitude.
[0019]
In the GPS altitude bias estimation unit 2, the inertia altitude 103 is input to the subtraction circuit 18 and the GPS altitude 112 is input to the subtraction circuit 21. At this time, since the subtracting circuit 18, the switching circuit 19, the integrating circuit 20, and the subtracting circuit 21 form a negative feedback loop, the integrating circuit 20 outputs a GPS altitude error 111 as an estimated error of the GPS altitude 112 with respect to the inertial altitude 103. Is output, and the GPS correction altitude 113 estimated with respect to the inertial altitude 103 is output from the subtraction circuit 21. The GPS altitude error 111 is a bias error of the GPS altitude 112 estimated on the basis of the inertial altitude 103, and does not vary even if time passes or the latitude and longitude of the GPS receiver 6 vary. The GPS corrected altitude 113 is an altitude obtained by correcting the GPS altitude 112 with the GPS altitude error 111.
[0020]
In a state where a negative feedback loop is configured in the GPS altitude bias estimation unit 2 and the GPS altitude error 111 is stably generated, the reference altitude selection signal 116 changes from the logical value “1” to the logical value “0”. If the movable piece of the switching circuit 19 is switched from the contact point a to the contact point b, the input 110 of the integration circuit 20 is opened, so that the integration circuit 20 holds the data at that time (GPS altitude error 111).
[0021]
As described above, in a state in which the atmospheric pressure altitude 104 is selected as the reference altitude, the inertial altitude correcting unit 1 stably outputs the inertial altitude 103 in a state of following the atmospheric altitude 104, and GPS altitude bias estimation. In the unit 2, a GPS altitude error 111 as an estimation error and an estimated GPS correction altitude 113 are generated.
[0022]
So far, the operation in a state where the altimeter determination circuit 7 determines that the atmospheric pressure altitude 104 is more effective than the GPS altitude 112 has been described. In this state, for example, the atmospheric pressure altitude 104 fluctuates in a stepwise manner and the reliability of the atmospheric pressure altitude 104 is suspected, or the altitude rises to the atmospheric pressure altitude measurement limit. Assume that the altimeter determination circuit 7 determines that the altitude 112 is more effective. At this time, the logical value of the reference altitude selection signal 116 becomes “0”, and all the movable segments in the switching circuits 15, 19 and 22 are switched to the contact b side. By this switching, the GPS altitude bias estimation unit 2 eliminates the negative feedback loop, and the GPS altitude error 111 at the time of switching is retained. The subtraction circuit 21 continues to generate a GPS corrected altitude 113 obtained by correcting the GPS altitude 112 with the GPS altitude error 111.
[0023]
Since the GPS altitude error 111 at this time is a bias error of the GPS altitude 112 generated based on the inertial altitude 103 when the barometric altitude 104 is selected as the reference altitude 114, the GPS altitude error 111 is used to correct the GPS altitude 112. The GPS corrected altitude 113 obtained in this way is continuous with the atmospheric pressure altitude 104 at the time of switching, and the GPS corrected altitude 113 after the switching follows the GPS altitude 112. Therefore, the value of the reference altitude 114 is continuous and does not change stepwise when the barometric altitude 104 is switched to the GPS corrected altitude 113. Thereafter, the inertial altitude correction unit 1 outputs the inertial altitude 103 corrected using the GPS corrected altitude 113 as the reference altitude 104. Of course, since the reference altitude 104 continues in switching, the inertial altitude 103 also continues in switching and does not fluctuate discontinuously in a staircase pattern.
The GPS corrected altitude 113 is input to the subtraction circuit 12 of the inertial altitude correction unit 1 as the reference altitude 114 via the switching circuit 22. In the inertial altitude correction unit 1, since the movable piece of the switching circuit 15 is switched to the contact point b, the integrating circuit 16 is removed from the first feedback circuit. Therefore, the difference 105 between the output of the integrating circuit 11 and the GPS correction altitude 113 is multiplied by a predetermined gain in the multiplication circuit 14 in the first feedback circuit, and becomes a multiplication value 107. The multiplication value 107 is input as negative feedback to the subtraction circuit 8 as it is through the adder circuit 17 as a feedback signal 108. The feedback signal 109 of the multiplication circuit 13 is negatively fed back to the subtraction circuit 10 as before the reference altitude 114 is switched. The integrating circuit 11 outputs the inertial altitude 103 generated by a new negative feedback loop excluding the integrating circuit 16 with the GPS corrected altitude 113 as the reference altitude 114.
[0024]
【The invention's effect】
As described above in detail with reference to the embodiment, according to the present invention, either the atmospheric pressure altitude or the GPS corrected altitude can be selected as the reference altitude, and the atmospheric altitude is selected by selecting the atmospheric altitude as the reference altitude. While measuring, a GPS corrected altitude is obtained by continuously generating and maintaining a GPS altitude bias error (GPS altitude error 111) estimated for the inertial altitude and correcting the GPS altitude by the bias error. By constantly generating (113), even when the reference altitude is switched from the barometric altitude to the GPS corrected altitude, the value of the reference altitude is continuous, and the inertial altitude can be measured without abrupt changes in steps. An inertial altitude measurement method that can be performed stably and smoothly can be obtained.
[Brief description of the drawings]
FIG. 1 is a block diagram showing an embodiment of the present invention.
FIG. 2 is a block diagram showing a conventional example.
[Explanation of symbols]
DESCRIPTION OF SYMBOLS 1 Inertial altitude correction | amendment part 2 GPS altitude bias estimation part 3 Reference | standard altitude selection part 4 Inertial accelerometer 5 Barometric altimeter 6 GPS receiver 7 Altimeter determination circuit 8, 10, 12, 18, 21, 23, 25, 31 Subtraction circuit 9, 11, 16, 20, 24, 26, 29 Integration circuit 13, 14, 27, 28 Multiplication circuit 15, 19, 22 Switching circuit 17, 30 Addition circuit 101, 117 Vertical acceleration 102, 118 Inertial speed 103, 119 Inertial altitude 104 Barometric altitude 105 Difference output 106 between inertial altitude 103 and reference altitude 104 Integration circuit 16 output 107 Multiplication value 108 of multiplication circuit 14 Feedback signal (addition value of addition circuit 17)
109 Multiplying value 110 of multiplication circuit 13 Switching circuit 19 Output 111 GPS altitude error 112 GPS altitude 113 GPS correction altitude 114 Reference altitude 113 GPS correction altitude 116 Reference altitude selection signal 119 Vertical movement displacement 120 Reference altitude 121 Inertial altitude 119 and reference altitude 120 Difference output 124 between and feedback signal 125 feedback signal (multiplier circuit 27 output)

Claims (3)

加速度計出力の垂直方向加速度および基準高度を処理する演算処理手段および該基準高度として気圧高度を出力する気圧高度計とを備え、該垂直方向加速度を該基準高度にダンピングさせることにより慣性高度を生成するフィードバックループを該演算処理手段に含み、該気圧高度が有効なときにだけに該基準高度として該気圧高度を選択する慣性高度測定方式において、
GPS受信機からGPS高度を受け、
前記気圧高度が有効な状態では、該気圧高度を前記基準高度として前記演算処理手段に供給して得た前記慣性高度を基準として前記GPS高度のバイアス誤差を推定し、
前記気圧高度が無効な状態になったとき、該気圧高度が無効な状態になる直前の前記バイアス誤差により前記GPS高度を補正したGPS補正高度を生成し、該GPS補正高度を前記基準高度として選択する
ことを特徴とする慣性高度測定方式。
Computational processing means for processing vertical acceleration and reference height of accelerometer output and a barometric altimeter that outputs barometric height as the reference height, and generating inertial height by damping the vertical acceleration to the reference height In an inertial altitude measurement method that includes a feedback loop in the processing means and selects the barometric altitude as the reference altitude only when the barometric altitude is valid,
Receive the GPS altitude from the GPS receiver,
In a state where the atmospheric pressure altitude is valid, a bias error of the GPS altitude is estimated based on the inertial altitude obtained by supplying the atmospheric altitude to the arithmetic processing unit as the reference altitude,
When the barometric altitude becomes invalid, a GPS corrected altitude is generated by correcting the GPS altitude by the bias error immediately before the barometric altitude becomes invalid, and the GPS corrected altitude is selected as the reference altitude. An inertial altitude measurement system characterized by
加速度計から出力される垂直方向加速度と基準高度とを処理することにより慣性高度を測定するために、第1、第2及び第3の減算手段、第1及び第2の積分手段並びに第1及び第2の負帰還手段を備え、該第1の負帰還手段は該第3の減算手段の出力に基づき第1の帰還値を生成し、該第2の負帰還手段は該第3の減算手段の出力に基づき第2の帰還値を生成し、該第1の減算手段は該垂直方向加速度と該第1の帰還値との減算をし、該第1の積分手段は該第1の減算手段の出力に関する時間積分をし、該第2の減算手段は該第1の積分手段の出力と該第2の帰還値との減算をし、該第2の積分手段は該第2の減算手段の出力に関する時間積分値を該慣性高度として出力し、該第3の減算手段は該慣性高度と該基準高度との減算をする慣性高度測定方式において、
気圧高度を出力する気圧高度計と、GPS高度を出力するGPS受信機と、GPS高度バイアス推定手段と、基準高度選択手段とを備え、
前記GPS高度バイアス推定手段は、前記慣性高度および前記GPS高度を入力し、該慣性高度に対する該GPS高度のバイアス誤差を推定し、該GPS高度を該バイアス誤差分だけ補正したGPS補正高度を生成し、
前記基準高度選択手段は、前記気圧高度が有効であるときは前記基準高度として該気圧高度を選択し、該気圧高度が無効であるときは該基準高度として前記GPS補正高度を選択する
ことを特徴とする慣性高度測定方式。
In order to measure the inertial height by processing the vertical acceleration output from the accelerometer and the reference height, first, second and third subtracting means, first and second integrating means, and first and Second negative feedback means, wherein the first negative feedback means generates a first feedback value based on the output of the third subtraction means, and the second negative feedback means is the third subtraction means. Based on the output of the second feedback value, the first subtracting means subtracts the vertical acceleration and the first feedback value, and the first integrating means is the first subtracting means. The second subtracting means subtracts the output of the first integrating means and the second feedback value, and the second integrating means is the second subtracting means of the second subtracting means. The time integral value relating to the output is output as the inertia height, and the third subtracting means is an inertia for subtracting the inertia height from the reference height. In time measurement method,
A barometric altimeter that outputs a barometric altitude; a GPS receiver that outputs a GPS altitude; a GPS altitude bias estimating means; and a reference altitude selecting means.
The GPS altitude bias estimation means inputs the inertial altitude and the GPS altitude, estimates a bias error of the GPS altitude relative to the inertial altitude, and generates a GPS corrected altitude in which the GPS altitude is corrected by the bias error. ,
The reference altitude selecting means selects the barometric altitude as the reference altitude when the barometric altitude is valid, and selects the GPS corrected altitude as the reference altitude when the barometric altitude is invalid. Inertial altitude measurement method.
前記GPS高度バイアス推定手段は、第4の減算手段と、第2の切替手段と、第3の積分手段と、第5の減算手段とを有してなり、
前記第4の減算手段は前記慣性高度と前記GPS補正高度との減算をし、
前記第2の切替手段は、前記基準高度選択手段が前記気圧高度を選択するときは前記第4の減算手段の出力をそのまま出力端に接続し、前記基準高度選択手段が前記GPS補正高度を選択するときは出力端を開放端とし、
前記第3の積分手段は、前記第2の切替手段の出力端の信号に関する時間積分をし、積分値を前記GPS高度のバイアス誤差として出力し、
前記第5の減算手段は、前記GPS高度と該バイアス誤差との減算をし、減算値を前記GPS補正高度として出力する
ことを特徴とする請求項2に記載の慣性高度測定方式。
The GPS altitude bias estimation means includes a fourth subtraction means, a second switching means, a third integration means, and a fifth subtraction means.
The fourth subtracting means subtracts the inertial altitude and the GPS corrected altitude,
The second switching means connects the output of the fourth subtracting means to the output terminal as it is when the reference altitude selecting means selects the barometric altitude, and the reference altitude selecting means selects the GPS corrected altitude. The output end is the open end,
The third integrating means performs time integration with respect to the signal at the output terminal of the second switching means, and outputs an integrated value as a bias error of the GPS altitude.
The inertial altitude measurement method according to claim 2, wherein the fifth subtracting unit subtracts the GPS altitude from the bias error and outputs a subtraction value as the GPS corrected altitude.
JP2002188623A 2002-06-27 2002-06-27 Inertial altitude measurement method Expired - Fee Related JP3711538B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2002188623A JP3711538B2 (en) 2002-06-27 2002-06-27 Inertial altitude measurement method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2002188623A JP3711538B2 (en) 2002-06-27 2002-06-27 Inertial altitude measurement method

Publications (2)

Publication Number Publication Date
JP2004028911A JP2004028911A (en) 2004-01-29
JP3711538B2 true JP3711538B2 (en) 2005-11-02

Family

ID=31183320

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2002188623A Expired - Fee Related JP3711538B2 (en) 2002-06-27 2002-06-27 Inertial altitude measurement method

Country Status (1)

Country Link
JP (1) JP3711538B2 (en)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009103488A (en) * 2007-10-22 2009-05-14 Ihi Aerospace Co Ltd Positioning device of missile
WO2009104254A1 (en) * 2008-02-20 2009-08-27 パイオニア株式会社 Position recognition device, method for correction, and program for correction
JP5094589B2 (en) * 2008-06-25 2012-12-12 ヤフー株式会社 Current position estimation apparatus, method and system
WO2013122716A2 (en) * 2012-01-20 2013-08-22 Sandel Avionics, Inc. Pressure altitude stabilization
KR101366015B1 (en) 2012-07-06 2014-02-24 한국전자통신연구원 Apparatus and method for measuring altitude
JP6318499B2 (en) * 2013-08-27 2018-05-09 株式会社ニコン Imaging device
JP6477950B2 (en) * 2018-04-03 2019-03-06 株式会社ニコン Imaging device
JP7145655B2 (en) * 2018-06-20 2022-10-03 三菱電機株式会社 Altitude measurement device and altitude measurement program
CN112882072B (en) * 2021-01-15 2022-03-08 中国人民解放军海军航空大学 Aircraft inertia and GPS hybrid altitude measurement method based on error feedback

Also Published As

Publication number Publication date
JP2004028911A (en) 2004-01-29

Similar Documents

Publication Publication Date Title
CN108241161B (en) Distributed Kalman filter architecture for carrier range ambiguity estimation
JP5046687B2 (en) Method and apparatus for implementing an iterative extended Kalman filter in a navigation system
US9915550B2 (en) Method and apparatus for data fusion of a three-axis magnetometer and three axis accelerometer
JP3711538B2 (en) Inertial altitude measurement method
US8922429B2 (en) Satellite navigation device
US7103477B1 (en) Self-calibration for an inertial instrument based on real time bias estimator
JPH0291513A (en) Method and device for correcting zero point of gyro
JPH01503170A (en) Vertical position stabilization control circuit in inertial navigation system
WO2009140184A1 (en) System and method of navigation based on state estimation using a stepped filter
CA2329262A1 (en) System and method for improving the accuracy of pressure altitude determinations in an inertial navigation system
JP2001174275A (en) Hybrid navigation and its apparatus
CN115047505A (en) GNSS positioning method and navigation method based on carrier phase differential assistance
JPH06317428A (en) Inertial navigation method
JPH07123763A (en) Speed estimating observer
JP3169213B2 (en) Moving speed detecting method and device, vehicle slip angle detecting device
JP6178776B2 (en) Speed measuring device
GB2534417B (en) Multiple sensor integration
JP2001343444A (en) Positioning device and method, and medium in which positioning program is recorded
JP6512155B2 (en) Output device
JP2008216050A (en) System and method for sensor drift compensation
US8200725B2 (en) Arithmetic processing system and method thereof
JPS61155813A (en) Estimating device for position of moving body
JPH0746874A (en) Rotational speed controller for motor
JP2006275860A (en) Automatic compensation system and method for gps cycle slip, and gps receiver and its program and recording medium
JP5914316B2 (en) Direction measuring device

Legal Events

Date Code Title Description
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20050210

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: 20050726

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20050805

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

Ref document number: 3711538

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080826

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090826

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090826

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100826

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100826

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100826

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110826

Year of fee payment: 6

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110826

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120826

Year of fee payment: 7

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120826

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120826

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130826

Year of fee payment: 8

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130826

Year of fee payment: 8

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

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

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