JP3501332B2 - 電波方向探知方法 - Google Patents
電波方向探知方法Info
- Publication number
- JP3501332B2 JP3501332B2 JP00889097A JP889097A JP3501332B2 JP 3501332 B2 JP3501332 B2 JP 3501332B2 JP 00889097 A JP00889097 A JP 00889097A JP 889097 A JP889097 A JP 889097A JP 3501332 B2 JP3501332 B2 JP 3501332B2
- Authority
- JP
- Japan
- Prior art keywords
- steering vector
- radio wave
- antenna
- wave
- calculating
- 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
Landscapes
- Radar Systems Or Details Thereof (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Description
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】この発明はミュージック法を
用いた方位測定方法、特に遠距離からの到来電波の方位
測定における評価関数の演算に用いるステアリングベク
トルを求める方法に関する。 【０００２】 【従来の技術】ミュージック法を用いた方向測定方法は
特願平４−２０５５２１号「方向探知機」、IEEE, Tran
s. ANTENA AND PROPAGATION, vol. AP-34, No.3, Marc
h, 1986, P.276 〜280 などで提案されている。この方
位測定方法は次式で示される評価関数ＰMU（φ，θ）に
おけるピークから電波の到来方向φm ，θm （φ：仰
角，θ：水平角）を求める。 ＰMU（φ，θ）≡a(φ，θ）* a(φ，θ)/(a（φ，θ）* ＥN ＥN * a(φ，θ)) ・・・ (1) ここでａ（φ，θ）はステアリングベクトル又はアンテ
ナ電圧ベクトルと呼ばれ、用いられるアンテナ装置に固
有な値であって、周波数が決まると、φ，θの関数とし
て与えられる。ＥN は雑音固有ベクトルからなる行列と
呼ばれ、各アンテナ素子の受信信号から求まる値、つま
り実測値と対応した値である。 *は複素共役転置を表わ
している。このようにミュージック法を用いた方向測定
方法においては、測定値より求めた固有ベクトルＥとス
テアリングベクトルａ（φ，θ）の内積の逆数を、各
（φ，θ）について演算して、ピークが得られた各
（φ，θ）からそれぞれ電波が到来していると判定して
いる。 【０００３】ところでステアリングベクトルａ（φ，
θ）はアンテナ装置に固有の値であってアンテナ素子間
に相互干渉がなければ理論的に求まる、つまり計算によ
り求まる。しかし、実際には、特にＶＨＦ帯やＵＨＦ帯
の円形等間隔アレーアンテナでは色々な制約からアンテ
ナ半径を大きくすることができず、アンテナ素子間の相
互干渉が無視できず、この干渉をも含めてステアリング
ベクトルを正確に計算により求めることは困難であり、
実測により求めていた。 【０００４】ＶＨＦ帯、ＵＨＦ帯においては、電離層反
射が生じないため仰角φの測定にあまり意味がなく、通
常水平方位角θを測定すればよい。従って、対象アンテ
ナ装置と送信機を適当に離して設け、対象アンテナ装置
を単位角度Δθずつ回転させて送信機よりの電波を受信
し、評価関数をＰMU（θ）を計算し、既知方位θからス
テアリングベクトルを求め、これを各周波数ごとに求め
る。このようにして求めたステアリングベクトルａ
（θ）を、実際の到来電波の方位測定に用いる。 【０００５】 【発明が解決しようとする課題】このように実測により
求めたステアリングベクトル（以下実測ステアリングベ
クトルと記す）を用いても、特に遠距離よりの到来電波
の方位測定に比較的大きな誤差が生じることが判明し
た。この原因を追求した結果、図２に示すように、電波
源１１から距離ｒ１離れた円形アレイアンテナ中心点Ｐ
１と、電波源１１及び点Ｐ１を結線と垂直に、点Ｐ１か
らｄ１だけ離れた点Ｐ２のアンテナ素子とにおける電波
源１１よりの受信電波の位相に差が生じる。つまり電波
源１１と点Ｐ２間の距離ｒ２とｒ１との差Δｒ＝√（ｒ
１2 ＋ｄ１2 ）−ｒ１だけ点Ｐ２に到達する電波経路が
長くなり、それだけ位相が遅れたものとなる。 【０００６】一方、電波源１１が点Ｐ１，Ｐ２から十分
遠くに在る場合は点Ｐ１，Ｐ２での各受信電波の位相差
が著しく小さくなる。ちなみにｒ１＝５０ｍ，ｄ1 ＝
０．３６ｍ，波長λ＝０．６ｍとするとλ／Δｒ≒４６
３と大きな値となるが、ｒ１＝５０ｍ，ｄ１＝２．５
ｍ，λ＝０．６ｍとするとλ／Δｒ≒２７となる。実測
ステアリングベクトルの測定は、多重反射波の影響、外
来電波の影響などがないようにする点から、アンテナ装
置と送信機間の距離は例えば５０ｍ程度が限度であり、
十分離すことはできない。つまり実測ステアリングベク
トルは電波の経路差が無視できない近距離における測定
であり、即ち、球面波での測定結果である。 【０００７】一方、電波の経路差が無視できる遠方から
の到来電波方位の測定においては、到来電波は平面波と
なっている。このように球面波で求めた実測ステアリン
グベクトルを、平面波の到来電波の方位測定に用いる点
に問題があると考えるに到った。 【０００８】 【課題を解決するための手段】この発明によれば使用ア
ンテナに対し、近距離での実測ステアリングベクトル、
つまりアンテナ素子間における経路差にもとずく位相差
（λを３６０°とした時のΔｒ／λ）が例えば５°以上
での実測ステアリングベクトルを求め、この実測ステア
リングベクトルからアンテナ相互干渉行列（較正行列）
を求め、このアンテナ相互干渉行列と、理論から求めた
平面波ステアリングベクトルとから、上記近距離よりも
遠距離からの到来電波、つまり経路差にもとづく位相差
が５°以下の到来電波の方位測定の評価関数の演算に用
いるステアリングベクトルを求める。 【０００９】つまり、実測ステアリングベクトルａ
R （θ）は理論上求めたステアリングベクトルａ（θ）
とアンテナ相互干渉行列Ｍとの積で表わせることが知ら
れている。アンテナ素子間の相互干渉は到来波が平面波
でも球面波でも同様と考えられ、つまりアンテナ相互干
渉行列Ｍは到来波が球面波であるか、平面波であるかに
関係しないと考えられる。即ち近距離実測ステアリング
ベクトルａN （θ）は球面波の理論上のステアリングベ
クトルａS （θ）とアンテナ相互干渉行列Ｍとの積と等
しい。 【００１０】 ａN （θ）＝ＭａS （θ） ・・・ (2) 同様に遠距離実測ステアリングベクトルａF （θ）と、
平面波の理論上のステアリングベクトルａP （θ）との
間に次の関係がある。 ａF （θ）＝ＭａP （θ） ・・・ (3) よって、Ｍを求めれば、式（３）によりａF （θ）を求
めることができる。アンテナ相互干渉行列Ｍを求める演
算は例えば論文「Sensor array calibration in the pr
esence of mutual coupling and unknown sensor gains
and phase」（ELECTRONICS LETTERS 3rd March 1994 v
ol.30, No.5, P.373) に記載されている。またＮ本の円
形アレーアンテナの場合、各アンテナ素子の座標を
Ｘm ，Ｙm （原点はアレーアンテナの中心Ｐ１）とする
と、平面波の理論上のステアリングベクトルａP （θ）
は次式で求まる。 【００１１】ａP （θ）＝［exp(ｊ２πｄ１（Ｘ1 sin
θ＋Ｙ1 cos θ）／λ），・・・・，exp(ｊ２πｄ１
（ＸN sin θ＋ＹN cos θ）／λ）］T 【００１２】 【発明の実施の形態】この発明の計算機シミュレーショ
ンによる実験例を示す。アンテナとして、素子が７本の
等間隔、半径ｄ１＝２．５ｍの円形アレーを用い、使用
周波数λ＝５００ＭＨｚ、到来電波数を１波、到来方位
を１００度とした。近距離（５０ｍ）の実測ステアリン
グベクトルａN （θ）を用いた場合、遠距離（無限大）
の実測ステアリングベクトルを用いた場合、この発明を
適用して求めた変換ステアリングベクトルａF （θ）＝
ＭａP （θ）を用いた場合の各評価関数ＰMUの演算結果
を図１にそれぞれ曲線２１，２２，２３として示す。 【００１３】この図１から明らかなように、従来技術で
ある曲線２１ではピークが小さく、正しく方位を測定で
きない。一方、この発明を適用した曲線２３は、遠距離
実測ステアリングベクトルを用いた曲線２２と、特にピ
ーク部分でほぼ一致し、高い鋭いピークが生じており、
この発明により正確に方位を測定できることが理解され
る。この発明の方法は、低い周波数帯ではアームが長く
なり、アームの強度を充分にしなければならないなど、
構造設計が難しい。価格も高くなるなどの点からＶＨＦ
帯、ＵＨＦ帯に適用して好ましいものである。 【００１４】 【発明の効果】以上述べたように、比較的狭い場所でも
近距離実測ステアリングベクトルａN（θ）は正しく求
めることができるため、この発明によれば経路差Δｒが
位相で表わして５度程度以上になるような条件で求めた
実測ステアリングベクトルを遠距離ステアリングベクト
ルａF （θ）に変換し、これを用いて評価関数の演算を
行うことにより、経路差Δｒが５度以下の小さい値とな
る遠方からの到来電波の方位を正確に測定することがで
きる。
用いた方位測定方法、特に遠距離からの到来電波の方位
測定における評価関数の演算に用いるステアリングベク
トルを求める方法に関する。 【０００２】 【従来の技術】ミュージック法を用いた方向測定方法は
特願平４−２０５５２１号「方向探知機」、IEEE, Tran
s. ANTENA AND PROPAGATION, vol. AP-34, No.3, Marc
h, 1986, P.276 〜280 などで提案されている。この方
位測定方法は次式で示される評価関数ＰMU（φ，θ）に
おけるピークから電波の到来方向φm ，θm （φ：仰
角，θ：水平角）を求める。 ＰMU（φ，θ）≡a(φ，θ）* a(φ，θ)/(a（φ，θ）* ＥN ＥN * a(φ，θ)) ・・・ (1) ここでａ（φ，θ）はステアリングベクトル又はアンテ
ナ電圧ベクトルと呼ばれ、用いられるアンテナ装置に固
有な値であって、周波数が決まると、φ，θの関数とし
て与えられる。ＥN は雑音固有ベクトルからなる行列と
呼ばれ、各アンテナ素子の受信信号から求まる値、つま
り実測値と対応した値である。 *は複素共役転置を表わ
している。このようにミュージック法を用いた方向測定
方法においては、測定値より求めた固有ベクトルＥとス
テアリングベクトルａ（φ，θ）の内積の逆数を、各
（φ，θ）について演算して、ピークが得られた各
（φ，θ）からそれぞれ電波が到来していると判定して
いる。 【０００３】ところでステアリングベクトルａ（φ，
θ）はアンテナ装置に固有の値であってアンテナ素子間
に相互干渉がなければ理論的に求まる、つまり計算によ
り求まる。しかし、実際には、特にＶＨＦ帯やＵＨＦ帯
の円形等間隔アレーアンテナでは色々な制約からアンテ
ナ半径を大きくすることができず、アンテナ素子間の相
互干渉が無視できず、この干渉をも含めてステアリング
ベクトルを正確に計算により求めることは困難であり、
実測により求めていた。 【０００４】ＶＨＦ帯、ＵＨＦ帯においては、電離層反
射が生じないため仰角φの測定にあまり意味がなく、通
常水平方位角θを測定すればよい。従って、対象アンテ
ナ装置と送信機を適当に離して設け、対象アンテナ装置
を単位角度Δθずつ回転させて送信機よりの電波を受信
し、評価関数をＰMU（θ）を計算し、既知方位θからス
テアリングベクトルを求め、これを各周波数ごとに求め
る。このようにして求めたステアリングベクトルａ
（θ）を、実際の到来電波の方位測定に用いる。 【０００５】 【発明が解決しようとする課題】このように実測により
求めたステアリングベクトル（以下実測ステアリングベ
クトルと記す）を用いても、特に遠距離よりの到来電波
の方位測定に比較的大きな誤差が生じることが判明し
た。この原因を追求した結果、図２に示すように、電波
源１１から距離ｒ１離れた円形アレイアンテナ中心点Ｐ
１と、電波源１１及び点Ｐ１を結線と垂直に、点Ｐ１か
らｄ１だけ離れた点Ｐ２のアンテナ素子とにおける電波
源１１よりの受信電波の位相に差が生じる。つまり電波
源１１と点Ｐ２間の距離ｒ２とｒ１との差Δｒ＝√（ｒ
１2 ＋ｄ１2 ）−ｒ１だけ点Ｐ２に到達する電波経路が
長くなり、それだけ位相が遅れたものとなる。 【０００６】一方、電波源１１が点Ｐ１，Ｐ２から十分
遠くに在る場合は点Ｐ１，Ｐ２での各受信電波の位相差
が著しく小さくなる。ちなみにｒ１＝５０ｍ，ｄ1 ＝
０．３６ｍ，波長λ＝０．６ｍとするとλ／Δｒ≒４６
３と大きな値となるが、ｒ１＝５０ｍ，ｄ１＝２．５
ｍ，λ＝０．６ｍとするとλ／Δｒ≒２７となる。実測
ステアリングベクトルの測定は、多重反射波の影響、外
来電波の影響などがないようにする点から、アンテナ装
置と送信機間の距離は例えば５０ｍ程度が限度であり、
十分離すことはできない。つまり実測ステアリングベク
トルは電波の経路差が無視できない近距離における測定
であり、即ち、球面波での測定結果である。 【０００７】一方、電波の経路差が無視できる遠方から
の到来電波方位の測定においては、到来電波は平面波と
なっている。このように球面波で求めた実測ステアリン
グベクトルを、平面波の到来電波の方位測定に用いる点
に問題があると考えるに到った。 【０００８】 【課題を解決するための手段】この発明によれば使用ア
ンテナに対し、近距離での実測ステアリングベクトル、
つまりアンテナ素子間における経路差にもとずく位相差
（λを３６０°とした時のΔｒ／λ）が例えば５°以上
での実測ステアリングベクトルを求め、この実測ステア
リングベクトルからアンテナ相互干渉行列（較正行列）
を求め、このアンテナ相互干渉行列と、理論から求めた
平面波ステアリングベクトルとから、上記近距離よりも
遠距離からの到来電波、つまり経路差にもとづく位相差
が５°以下の到来電波の方位測定の評価関数の演算に用
いるステアリングベクトルを求める。 【０００９】つまり、実測ステアリングベクトルａ
R （θ）は理論上求めたステアリングベクトルａ（θ）
とアンテナ相互干渉行列Ｍとの積で表わせることが知ら
れている。アンテナ素子間の相互干渉は到来波が平面波
でも球面波でも同様と考えられ、つまりアンテナ相互干
渉行列Ｍは到来波が球面波であるか、平面波であるかに
関係しないと考えられる。即ち近距離実測ステアリング
ベクトルａN （θ）は球面波の理論上のステアリングベ
クトルａS （θ）とアンテナ相互干渉行列Ｍとの積と等
しい。 【００１０】 ａN （θ）＝ＭａS （θ） ・・・ (2) 同様に遠距離実測ステアリングベクトルａF （θ）と、
平面波の理論上のステアリングベクトルａP （θ）との
間に次の関係がある。 ａF （θ）＝ＭａP （θ） ・・・ (3) よって、Ｍを求めれば、式（３）によりａF （θ）を求
めることができる。アンテナ相互干渉行列Ｍを求める演
算は例えば論文「Sensor array calibration in the pr
esence of mutual coupling and unknown sensor gains
and phase」（ELECTRONICS LETTERS 3rd March 1994 v
ol.30, No.5, P.373) に記載されている。またＮ本の円
形アレーアンテナの場合、各アンテナ素子の座標を
Ｘm ，Ｙm （原点はアレーアンテナの中心Ｐ１）とする
と、平面波の理論上のステアリングベクトルａP （θ）
は次式で求まる。 【００１１】ａP （θ）＝［exp(ｊ２πｄ１（Ｘ1 sin
θ＋Ｙ1 cos θ）／λ），・・・・，exp(ｊ２πｄ１
（ＸN sin θ＋ＹN cos θ）／λ）］T 【００１２】 【発明の実施の形態】この発明の計算機シミュレーショ
ンによる実験例を示す。アンテナとして、素子が７本の
等間隔、半径ｄ１＝２．５ｍの円形アレーを用い、使用
周波数λ＝５００ＭＨｚ、到来電波数を１波、到来方位
を１００度とした。近距離（５０ｍ）の実測ステアリン
グベクトルａN （θ）を用いた場合、遠距離（無限大）
の実測ステアリングベクトルを用いた場合、この発明を
適用して求めた変換ステアリングベクトルａF （θ）＝
ＭａP （θ）を用いた場合の各評価関数ＰMUの演算結果
を図１にそれぞれ曲線２１，２２，２３として示す。 【００１３】この図１から明らかなように、従来技術で
ある曲線２１ではピークが小さく、正しく方位を測定で
きない。一方、この発明を適用した曲線２３は、遠距離
実測ステアリングベクトルを用いた曲線２２と、特にピ
ーク部分でほぼ一致し、高い鋭いピークが生じており、
この発明により正確に方位を測定できることが理解され
る。この発明の方法は、低い周波数帯ではアームが長く
なり、アームの強度を充分にしなければならないなど、
構造設計が難しい。価格も高くなるなどの点からＶＨＦ
帯、ＵＨＦ帯に適用して好ましいものである。 【００１４】 【発明の効果】以上述べたように、比較的狭い場所でも
近距離実測ステアリングベクトルａN（θ）は正しく求
めることができるため、この発明によれば経路差Δｒが
位相で表わして５度程度以上になるような条件で求めた
実測ステアリングベクトルを遠距離ステアリングベクト
ルａF （θ）に変換し、これを用いて評価関数の演算を
行うことにより、経路差Δｒが５度以下の小さい値とな
る遠方からの到来電波の方位を正確に測定することがで
きる。
【図面の簡単な説明】
【図１】この発明の効果を示す計算機シミューレション
結果を示す図。 【図２】近距離到来電波の場合の経路差を示す図。
結果を示す図。 【図２】近距離到来電波の場合の経路差を示す図。
─────────────────────────────────────────────────────
フロントページの続き
(56)参考文献 特開 平６−317529（ＪＰ，Ａ)
特開 平９−219616（ＪＰ，Ａ)
(58)調査した分野(Int.Cl.7，ＤＢ名)
G01S 3/00 - 3/74
Claims (1)
- (57)【特許請求の範囲】 【請求項１】 複数のアンテナからの各受信信号を、そ
れぞれフーリエ変換し、これらフーリエ変換出力をミュ
ージック法で処理して上記アンテナに受信された電波の
到来方向を演算する電波方向探知方法において、 上記ミュージック法における評価関数の演算に用いるス
テアリングベクトルを、近距離で実測して求め、 その実測ステアリングベクトルからアンテナ相互干渉行
列を求め、 上記アンテナ相互干渉行列と、理論上の平面波ステアリ
ングベクトルとから、上記近距離よりも遠距離からの到
来電波の方位測定における上記評価関数の演算に用いる
ステアリングベクトルを求めることを特徴とする電波方
向探知方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP00889097A JP3501332B2 (ja) | 1997-01-21 | 1997-01-21 | 電波方向探知方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP00889097A JP3501332B2 (ja) | 1997-01-21 | 1997-01-21 | 電波方向探知方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH10206517A JPH10206517A (ja) | 1998-08-07 |
| JP3501332B2 true JP3501332B2 (ja) | 2004-03-02 |
Family
ID=11705282
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP00889097A Expired - Fee Related JP3501332B2 (ja) | 1997-01-21 | 1997-01-21 | 電波方向探知方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3501332B2 (ja) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007256171A (ja) * | 2006-03-24 | 2007-10-04 | Nec Corp | ミリ波画像処理装置及びミリ波画像処理方法 |
-
1997
- 1997-01-21 JP JP00889097A patent/JP3501332B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JPH10206517A (ja) | 1998-08-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Friedlander | Localization of signals in the near-field of an antenna array | |
| EP1055268B1 (en) | Determining angles-of-arrival and polarization of rf signals | |
| US7880684B2 (en) | Small aperture broadband localizing system | |
| US20140104111A1 (en) | Method for calculating spacing ratio of interferometer array antenna for direction finder | |
| CN110673086A (zh) | 一种基于数字阵列雷达的二维角度超分辨方法 | |
| US6407702B1 (en) | Method and system for obtaining direction of an electromagnetic wave | |
| WO2000039601A1 (en) | Compensation of direction finding estimates for polarimetric errors | |
| JP2000121716A (ja) | 電波到来方向推定装置 | |
| US12032053B2 (en) | Angle measuring device, angle measuring method, and in-vehicle device | |
| US20190064339A1 (en) | Angle Estimating Method and Radar System | |
| US6195043B1 (en) | AOA estimation and polarization induced phase correction using a conformal array of titled antenna elements | |
| JP6584714B1 (ja) | レーダ装置及び目標測角方法 | |
| KR100527848B1 (ko) | 근역 측정을 위한 안테나 정합 장치 | |
| US20230314127A1 (en) | Angle width estimation device, angle width estimation method, and target angle measurement device | |
| JP3501332B2 (ja) | 電波方向探知方法 | |
| Park et al. | Dynamically reconfigurable direction-finding antenna array for unmanned arial systems | |
| US9939513B2 (en) | Apparatus and method for finding hybrid direction using two baselines | |
| Tan et al. | Ambiguity in MUSIC and ESPRIT for direction of arrival estimation | |
| US8004463B2 (en) | Systems and methods for determining direction-of-arrival | |
| JP2982751B2 (ja) | 位置標定方法及び装置 | |
| JP2002107440A (ja) | 方向検出方法及び方向検出装置 | |
| JP3875198B2 (ja) | 電波到来方向探知方法及び装置 | |
| CN109490875B (zh) | 角度估测方法及雷达系统 | |
| JPS61284680A (ja) | 方位測定装置 | |
| US11733367B1 (en) | Chromatic correlation interferometry direction finding |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| RD01 | Notification of change of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7426 Effective date: 20031125 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20031127 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081212 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081212 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091212 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091212 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101212 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101212 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111212 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111212 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121212 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131212 Year of fee payment: 10 |
|
| LAPS | Cancellation because of no payment of annual fees |