JP2000111630A - 電波到来方向推定方法及びその装置 - Google Patents
電波到来方向推定方法及びその装置Info
- Publication number
- JP2000111630A JP2000111630A JP10284433A JP28443398A JP2000111630A JP 2000111630 A JP2000111630 A JP 2000111630A JP 10284433 A JP10284433 A JP 10284433A JP 28443398 A JP28443398 A JP 28443398A JP 2000111630 A JP2000111630 A JP 2000111630A
- Authority
- JP
- Japan
- Prior art keywords
- signal
- arrival
- signals
- radio wave
- normalized
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 239000011159 matrix material Substances 0.000 claims abstract description 73
- 238000012545 processing Methods 0.000 claims abstract description 52
- 238000000926 separation method Methods 0.000 claims abstract description 6
- 238000012937 correction Methods 0.000 claims description 25
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 claims description 23
- 238000012935 Averaging Methods 0.000 claims description 7
- 238000011156 evaluation Methods 0.000 claims description 6
- 230000001934 delay Effects 0.000 abstract description 6
- 238000010586 diagram Methods 0.000 description 36
- 230000001186 cumulative effect Effects 0.000 description 13
- 230000005540 biological transmission Effects 0.000 description 6
- 238000010606 normalization Methods 0.000 description 5
- 239000000284 extract Substances 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 239000010432 diamond Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000004913 activation Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
(57)【要約】
【課題】 到来信号遅延時間のオフセットをなくし、到
来角推定部８で常時各到来信号の到来方向を正確に推定
できる電波到来方向推定装置を提供する。 【解決手段】 近接配置のＮ個のアンテナ１1 〜１4 、
受信信号からＮ個の正規化信号を得る正規化信号発生部
２1 〜２4 、各正規化信号から自己相関行列を得る自己
相関行列発生部３1 〜３4 、平均化自己相関行列を得る
平均化自己相関行列発生部４、平均化自己相関行列を処
理し、到来信号数（Ｋ個）と各到来信号遅延時間を発生
する高分解能信号処理部５、各正規化信号からＭ個の相
互相関行列を発生する相互相関行列発生部６、各相互相
関行列と到来信号数と各到来信号遅延時間からＫ個の到
来信号を各Ｍ通り発生する到来信号発生部７、各到来信
号遅延時間を、各Ｍ通り発生した各到来信号振幅に基づ
き誤差評価を行い、誤差が最小な各補正到来信号遅延時
間を到来信号発生部７に加える遅延時間補正部２６、各
Ｍ通り発生したＫ個の到来信号位相から各到来信号の到
来方向を推定する到来角推定部８を備える。
来角推定部８で常時各到来信号の到来方向を正確に推定
できる電波到来方向推定装置を提供する。 【解決手段】 近接配置のＮ個のアンテナ１1 〜１4 、
受信信号からＮ個の正規化信号を得る正規化信号発生部
２1 〜２4 、各正規化信号から自己相関行列を得る自己
相関行列発生部３1 〜３4 、平均化自己相関行列を得る
平均化自己相関行列発生部４、平均化自己相関行列を処
理し、到来信号数（Ｋ個）と各到来信号遅延時間を発生
する高分解能信号処理部５、各正規化信号からＭ個の相
互相関行列を発生する相互相関行列発生部６、各相互相
関行列と到来信号数と各到来信号遅延時間からＫ個の到
来信号を各Ｍ通り発生する到来信号発生部７、各到来信
号遅延時間を、各Ｍ通り発生した各到来信号振幅に基づ
き誤差評価を行い、誤差が最小な各補正到来信号遅延時
間を到来信号発生部７に加える遅延時間補正部２６、各
Ｍ通り発生したＫ個の到来信号位相から各到来信号の到
来方向を推定する到来角推定部８を備える。
Description
【０００１】
【発明の属する技術分野】本発明は、電波到来方向推定
方法及びその装置に係わり、特に、Ｎ（ここで、Ｎは３
以上の整数）個のアンテナで捉えた信号電波の中で時間
的に近接して到来する複数の信号成分のそれぞれの到来
方向を安定に、かつ、正確に推定することが可能な電波
到来方向推定方法及びその装置に関する。
方法及びその装置に係わり、特に、Ｎ（ここで、Ｎは３
以上の整数）個のアンテナで捉えた信号電波の中で時間
的に近接して到来する複数の信号成分のそれぞれの到来
方向を安定に、かつ、正確に推定することが可能な電波
到来方向推定方法及びその装置に関する。
【０００２】
【従来の技術】従来、所定領域内を適宜移動する移動体
の現在位置を追尾するために、電波を利用した移動体追
尾方式が知られている。この移動体追尾方式は、移動体
が信号電波を送信する送信機を携帯し、送信機から放射
された信号電波を受ける複数のアンテナを備えた基地局
を所要領域の近傍に配置しているものである。基地局
は、時間的に近接して到来する複数の信号電波の中で最
も早く到来する信号電波成分、即ち、到来信号電波成分
の中で最も遅延時間の短い信号成分を主信号成分として
抽出し、この主信号成分の到来方向を求めることによ
り、信号電波の発信源の位置を特定することができる。
そして、信号電波の発信源が移動体の携帯する送信機で
あれば、移動体の現在位置を追尾することができる。
の現在位置を追尾するために、電波を利用した移動体追
尾方式が知られている。この移動体追尾方式は、移動体
が信号電波を送信する送信機を携帯し、送信機から放射
された信号電波を受ける複数のアンテナを備えた基地局
を所要領域の近傍に配置しているものである。基地局
は、時間的に近接して到来する複数の信号電波の中で最
も早く到来する信号電波成分、即ち、到来信号電波成分
の中で最も遅延時間の短い信号成分を主信号成分として
抽出し、この主信号成分の到来方向を求めることによ
り、信号電波の発信源の位置を特定することができる。
そして、信号電波の発信源が移動体の携帯する送信機で
あれば、移動体の現在位置を追尾することができる。
【０００３】この場合、このような移動体追尾方式にお
いては、移動体が携帯する送信機と基地局との間で送受
信される信号変調方式に、種々の信号変調方式を用いる
ことが可能であるが、その信号変調方式の中の１つにＰ
Ｎ（疑似ランダムノイズ）符号を用いた拡散変調方式が
ある。
いては、移動体が携帯する送信機と基地局との間で送受
信される信号変調方式に、種々の信号変調方式を用いる
ことが可能であるが、その信号変調方式の中の１つにＰ
Ｎ（疑似ランダムノイズ）符号を用いた拡散変調方式が
ある。
【０００４】ここで、図９は、前記移動体追尾方式に用
いられる既知の電波到来方向推定装置の構成の一例を示
すブロック図であり、４個のアンテナを用いる場合の構
成例を示すものである。
いられる既知の電波到来方向推定装置の構成の一例を示
すブロック図であり、４個のアンテナを用いる場合の構
成例を示すものである。
【０００５】図９に示されるように、既知の電波到来方
向推定装置は、４個のアンテナ１1乃至１4 と、４個の
正規化信号発生部２1 乃至２4 と、４個の自己相関行列
発生部３1 乃至３4 と、平均化自己相関行列発生部４
と、高分解能信号処理部５と、相互相関行列発生部６
と、到来信号発生部７と、到来角推定部８と、信号出力
端子９とからなっている。
向推定装置は、４個のアンテナ１1乃至１4 と、４個の
正規化信号発生部２1 乃至２4 と、４個の自己相関行列
発生部３1 乃至３4 と、平均化自己相関行列発生部４
と、高分解能信号処理部５と、相互相関行列発生部６
と、到来信号発生部７と、到来角推定部８と、信号出力
端子９とからなっている。
【０００６】そして、正規化信号発生部２1 は、入力が
アンテナ１1 に接続され、出力が自己相関行列発生部３
1 の入力と相互相関行列発生部６の対応する入力に接続
される。正規化信号発生部２2 は、入力がアンテナ１2
に接続され、出力が自己相関行列発生部３2 の入力と相
互相関行列発生部６の対応する入力に接続される。正規
化信号発生部２3 は、入力がアンテナ１3 に接続され、
出力が自己相関行列発生部３3 の入力と相互相関行列発
生部６の対応する入力に接続される。正規化信号発生部
２4 は、入力がアンテナ１4 に接続され、出力が自己相
関行列発生部３4 の入力と相互相関行列発生部６の対応
する入力に接続される。平均化自己相関行列発生部４
は、４つの入力がそれぞれ自己相関行列発生部３1 乃至
３4 の各出力に接続され、出力が高分解能信号処理部５
の入力に接続される。到来信号発生部７は、６つの入力
がそれぞれ相互相関行列発生部６の対応する出力に接続
されるとともに、２つの入力が高分解能信号処理部５の
対応する２つの出力（到来信号数出力、信号遅延時間出
力）に接続され、６つの出力が到来角推定部８の対応す
る６つの入力にそれぞれ接続される。到来角推定部８
は、出力が信号出力端子９に接続される。この場合、高
分解能信号処理部５には、例えば、マルチプル・シグナ
ル・クラシフィケーション（ＭＵＳＩＣ）法を用いた高
時間分解能信号処理装置である。
アンテナ１1 に接続され、出力が自己相関行列発生部３
1 の入力と相互相関行列発生部６の対応する入力に接続
される。正規化信号発生部２2 は、入力がアンテナ１2
に接続され、出力が自己相関行列発生部３2 の入力と相
互相関行列発生部６の対応する入力に接続される。正規
化信号発生部２3 は、入力がアンテナ１3 に接続され、
出力が自己相関行列発生部３3 の入力と相互相関行列発
生部６の対応する入力に接続される。正規化信号発生部
２4 は、入力がアンテナ１4 に接続され、出力が自己相
関行列発生部３4 の入力と相互相関行列発生部６の対応
する入力に接続される。平均化自己相関行列発生部４
は、４つの入力がそれぞれ自己相関行列発生部３1 乃至
３4 の各出力に接続され、出力が高分解能信号処理部５
の入力に接続される。到来信号発生部７は、６つの入力
がそれぞれ相互相関行列発生部６の対応する出力に接続
されるとともに、２つの入力が高分解能信号処理部５の
対応する２つの出力（到来信号数出力、信号遅延時間出
力）に接続され、６つの出力が到来角推定部８の対応す
る６つの入力にそれぞれ接続される。到来角推定部８
は、出力が信号出力端子９に接続される。この場合、高
分解能信号処理部５には、例えば、マルチプル・シグナ
ル・クラシフィケーション（ＭＵＳＩＣ）法を用いた高
時間分解能信号処理装置である。
【０００７】また、図２は、図９に図示された各正規化
信号発生部２1 乃至２4 の構成の一例を示すブロック図
であり、後述する本発明の電波到来方向推定装置に用い
られる各正規化信号発生部の構成と同じ構成のものであ
る。
信号発生部２1 乃至２4 の構成の一例を示すブロック図
であり、後述する本発明の電波到来方向推定装置に用い
られる各正規化信号発生部の構成と同じ構成のものであ
る。
【０００８】図２に示すように、各正規化信号発生部２
1 乃至２4 は、同じ構成のもので、受信部１０と、第１
フーリエ変換部１１と、第１フィルタ手段１２と、参照
信号発生部１３と、第２フーリエ変換部１４と、第２フ
ィルタ手段１５と、除算部１６と、雑音除去処理部１７
と、受信信号入力端子１８と、正規化信号出力端子１９
とからなっている。この場合、受信部１０は、ベースバ
ンド（以下、ＢＢという）信号発生部２０と、アナログ
−ディジタル（以下、Ａ／Ｄという）変換部２１と、メ
モリ２２とからなり、雑音除去処理部１７は、周波数窓
乗算部２３と、フーリエ変換部２４と、雑音除去部２５
と、逆フーリエ変換部２６と、周波数窓除算部２７とか
らなる。
1 乃至２4 は、同じ構成のもので、受信部１０と、第１
フーリエ変換部１１と、第１フィルタ手段１２と、参照
信号発生部１３と、第２フーリエ変換部１４と、第２フ
ィルタ手段１５と、除算部１６と、雑音除去処理部１７
と、受信信号入力端子１８と、正規化信号出力端子１９
とからなっている。この場合、受信部１０は、ベースバ
ンド（以下、ＢＢという）信号発生部２０と、アナログ
−ディジタル（以下、Ａ／Ｄという）変換部２１と、メ
モリ２２とからなり、雑音除去処理部１７は、周波数窓
乗算部２３と、フーリエ変換部２４と、雑音除去部２５
と、逆フーリエ変換部２６と、周波数窓除算部２７とか
らなる。
【０００９】そして、受信部１０は、入力が受信信号入
力端子１８に接続され、出力が第１フーリエ変換部１１
の入力に接続される。第１フィルタ手段１２は、入力が
第１フーリエ変換部１１の出力に接続され、出力が除算
部１６の第１入力に接続される。参照信号発生部１３
は、出力が第２フーリエ変換部１４の入力に接続され
る。第２フィルタ手段１５は、入力が第２フーリエ変換
部１４の出力に接続され、出力が除算部１６の第２入力
に接続される。雑音除去処理部１７は、入力が除算部１
６の出力に接続され、出力が正規化信号出力端子１９に
接続される。
力端子１８に接続され、出力が第１フーリエ変換部１１
の入力に接続される。第１フィルタ手段１２は、入力が
第１フーリエ変換部１１の出力に接続され、出力が除算
部１６の第１入力に接続される。参照信号発生部１３
は、出力が第２フーリエ変換部１４の入力に接続され
る。第２フィルタ手段１５は、入力が第２フーリエ変換
部１４の出力に接続され、出力が除算部１６の第２入力
に接続される。雑音除去処理部１７は、入力が除算部１
６の出力に接続され、出力が正規化信号出力端子１９に
接続される。
【００１０】また、受信部１０において、ＢＢ信号発生
部２０は、入力が受信部１０の入力に接続され、出力が
Ａ／Ｄ変換部２１の入力に接続される。メモリ２２は、
入力がＡ／Ｄ変換部２１の出力に接続され、出力が受信
部１０の出力に接続される。雑音除去処理部１７におい
て、周波数窓乗算部２３は、入力が雑音除去処理部１７
の入力に接続され、出力がフーリエ変換部２４の入力に
接続される。雑音除去部２５は、入力がフーリエ変換部
２４の出力に接続され、出力が逆フーリエ変換部２６の
入力に接続される。周波数窓除算部２７は、入力が逆フ
ーリエ変換部２６の出力に接続され、出力が雑音除去処
理部１７の出力に接続される。
部２０は、入力が受信部１０の入力に接続され、出力が
Ａ／Ｄ変換部２１の入力に接続される。メモリ２２は、
入力がＡ／Ｄ変換部２１の出力に接続され、出力が受信
部１０の出力に接続される。雑音除去処理部１７におい
て、周波数窓乗算部２３は、入力が雑音除去処理部１７
の入力に接続され、出力がフーリエ変換部２４の入力に
接続される。雑音除去部２５は、入力がフーリエ変換部
２４の出力に接続され、出力が逆フーリエ変換部２６の
入力に接続される。周波数窓除算部２７は、入力が逆フ
ーリエ変換部２６の出力に接続され、出力が雑音除去処
理部１７の出力に接続される。
【００１１】続く、図１０は、図９に図示された４個の
アンテナ１1 乃至１4 の配置状態と信号電波の到来状態
の一例を示す構成図である。
アンテナ１1 乃至１4 の配置状態と信号電波の到来状態
の一例を示す構成図である。
【００１２】図１０に示されるように、４個のアンテナ
１1 乃至１4 は、基準点Ｏを通って延びるＸ軸と、基準
点Ｏを通ってＸ軸と直交する方向に延びるＹ軸からなる
平面上の、基準点Ｏを中心とする半径ｒ’（ここで、
ｒ’は、信号電波の波長をλとした場合、０．３６λに
選ばれる）の円周上に等角度間隔に配置され、４素子円
形アレーアンテナを構成する。
１1 乃至１4 は、基準点Ｏを通って延びるＸ軸と、基準
点Ｏを通ってＸ軸と直交する方向に延びるＹ軸からなる
平面上の、基準点Ｏを中心とする半径ｒ’（ここで、
ｒ’は、信号電波の波長をλとした場合、０．３６λに
選ばれる）の円周上に等角度間隔に配置され、４素子円
形アレーアンテナを構成する。
【００１３】この場合、信号電波はＰＮ（疑似ランダム
ノイズ）符号で拡散変調された信号を含む伝搬波であっ
て、その第１波は、前記平面に平行で、Ｘ軸に対して５
０°の方位角方向から到来し、その第２波は、前記平面
に平行で、Ｘ軸に対して０°の方位角方向から到来する
ものとし、第１波に対して第２波が相対的に０．２５チ
ップ（Ｔｃ）だけ遅延している例を示している。
ノイズ）符号で拡散変調された信号を含む伝搬波であっ
て、その第１波は、前記平面に平行で、Ｘ軸に対して５
０°の方位角方向から到来し、その第２波は、前記平面
に平行で、Ｘ軸に対して０°の方位角方向から到来する
ものとし、第１波に対して第２波が相対的に０．２５チ
ップ（Ｔｃ）だけ遅延している例を示している。
【００１４】また、図１１は、送信側において、送信デ
ータ（周期Ｔ）がＰＮ符号で拡散変調され、さらにその
帯域を制限した帯域制限拡散変調信号の一例を示す信号
波形図である。
ータ（周期Ｔ）がＰＮ符号で拡散変調され、さらにその
帯域を制限した帯域制限拡散変調信号の一例を示す信号
波形図である。
【００１５】図１２は、４個の正規化信号発生部２1 乃
至２4 の中の１つ、例えば、正規化信号発生部２1 のＢ
Ｂ信号発生部２０から出力されるＢＢ拡散変調信号の一
例を示す信号波形図である。図中、実線で示される曲線
ａは同相チャネルの信号、点線で示される曲線ｂは直交
チャネルの信号である。
至２4 の中の１つ、例えば、正規化信号発生部２1 のＢ
Ｂ信号発生部２０から出力されるＢＢ拡散変調信号の一
例を示す信号波形図である。図中、実線で示される曲線
ａは同相チャネルの信号、点線で示される曲線ｂは直交
チャネルの信号である。
【００１６】図１３は、正規化信号発生部２1 の第２フ
ーリエ変換部１４から出力される周波数領域参照信号の
一例を示す信号波形図である。図中、周波数帯域幅Ａの
部分は、第２フィルタ手段１５で抽出される周波数領域
である。
ーリエ変換部１４から出力される周波数領域参照信号の
一例を示す信号波形図である。図中、周波数帯域幅Ａの
部分は、第２フィルタ手段１５で抽出される周波数領域
である。
【００１７】図１４は、正規化信号発生部２1 の第１フ
ーリエ変換部１１から出力される周波数領域受信信号の
一例を示す信号波形図である。図中、周波数帯域幅Ｂの
部分は、第１フィルタ手段１２で抽出される周波数領域
である。
ーリエ変換部１１から出力される周波数領域受信信号の
一例を示す信号波形図である。図中、周波数帯域幅Ｂの
部分は、第１フィルタ手段１２で抽出される周波数領域
である。
【００１８】図１５は、雑音除去処理部１７に入力され
る雑音成分を含む正規化信号、及び、雑音除去処理部１
７から出力される雑音成分が除かれた正規化信号の各一
例を示す信号波形図である。図中、点線で示される曲線
ａは雑音成分を含んだ正規化信号、実線で示される曲線
ｂは雑音成分を除去した正規化信号である。
る雑音成分を含む正規化信号、及び、雑音除去処理部１
７から出力される雑音成分が除かれた正規化信号の各一
例を示す信号波形図である。図中、点線で示される曲線
ａは雑音成分を含んだ正規化信号、実線で示される曲線
ｂは雑音成分を除去した正規化信号である。
【００１９】図１６は、高分解能信号処理部５から出力
される相関計量信号の一例を示す信号波形図である。
される相関計量信号の一例を示す信号波形図である。
【００２０】前記構成による既知の電波到来方向推定装
置の動作について、図２、図９乃至図１６を参照して説
明する。
置の動作について、図２、図９乃至図１６を参照して説
明する。
【００２１】この場合、送信側（図示なし）から送信さ
れる信号電波は、送信データをＰＳＫ（位相シフトキー
イング）等の変調を行って形成した１次変調信号にＰＮ
符号を乗算して拡散変調信号を形成し、さらにこの拡散
変調信号の帯域を制限して図１１に示されるような帯域
制限拡散変調信号を形成し、この帯域制限拡散変調信号
が送信信号に周波数変換されて送信されるものである。
れる信号電波は、送信データをＰＳＫ（位相シフトキー
イング）等の変調を行って形成した１次変調信号にＰＮ
符号を乗算して拡散変調信号を形成し、さらにこの拡散
変調信号の帯域を制限して図１１に示されるような帯域
制限拡散変調信号を形成し、この帯域制限拡散変調信号
が送信信号に周波数変換されて送信されるものである。
【００２２】また、この電波到来方向推定装置の動作説
明の中の４個の正規化信号発生部２1 乃至２4 の動作に
ついては、それぞれ、４個のアンテナ１1 乃至１4 が受
信する信号電波に応じて取扱う信号波形が異なるだけで
同じ動作が行なわれる。このため、ここでは正規化信号
発生部２1 の動作だけを説明し、他の正規化信号発生部
２2 乃至２4 の動作の説明は省略する。
明の中の４個の正規化信号発生部２1 乃至２4 の動作に
ついては、それぞれ、４個のアンテナ１1 乃至１4 が受
信する信号電波に応じて取扱う信号波形が異なるだけで
同じ動作が行なわれる。このため、ここでは正規化信号
発生部２1 の動作だけを説明し、他の正規化信号発生部
２2 乃至２4 の動作の説明は省略する。
【００２３】送信側から図１１に示されるような拡散変
調された信号が周波数変換されて送信されると、図１０
に示されるように、４個のアンテナ１1 乃至１4 は、そ
れぞれ、第１波及び第２波を含む信号電波を捉える。
調された信号が周波数変換されて送信されると、図１０
に示されるように、４個のアンテナ１1 乃至１4 は、そ
れぞれ、第１波及び第２波を含む信号電波を捉える。
【００２４】４個のアンテナ１1 乃至１4 で捉えられた
信号電波は、それぞれ、受信信号として４個の正規化信
号発生部２1 乃至２4 の対応する受信部１０に供給され
る。このとき、正規化信号発生部２1 においては、以下
に述べるような動作が行なわれる。
信号電波は、それぞれ、受信信号として４個の正規化信
号発生部２1 乃至２4 の対応する受信部１０に供給され
る。このとき、正規化信号発生部２1 においては、以下
に述べるような動作が行なわれる。
【００２５】まず、受信部１０のＢＢ信号発生部２０
は、図１２に示されるようなＢＢアナログ拡散変調信号
を再生し、Ａ／Ｄ変換部２１は、ＢＢアナログ拡散変調
信号をＢＢディジタル拡散変調信号に変換してメモリ２
２に出力する。
は、図１２に示されるようなＢＢアナログ拡散変調信号
を再生し、Ａ／Ｄ変換部２１は、ＢＢアナログ拡散変調
信号をＢＢディジタル拡散変調信号に変換してメモリ２
２に出力する。
【００２６】次に、第１フーリエ変換部１１は、ＢＢデ
ィジタル拡散変調信号をフーリエ変換して図１４に示さ
れるような信号波形の周波数領域受信信号に変換する。
ィジタル拡散変調信号をフーリエ変換して図１４に示さ
れるような信号波形の周波数領域受信信号に変換する。
【００２７】次いで、第１フィルタ手段１２は、周波数
領域受信信号の周波数帯域を図１４に図示された帯域幅
Ｂの範囲になるように制限し、除算部１６に供給する。
領域受信信号の周波数帯域を図１４に図示された帯域幅
Ｂの範囲になるように制限し、除算部１６に供給する。
【００２８】一方、参照信号発生部１３は、参照信号と
して図１１に図示されたような帯域制限ＰＮ符号を発生
する。この帯域制限ＰＮ符号は、送信側において送信デ
ータを拡散変調する際に用いたＰＮ符号と同じ符号を帯
域制限して得たものであり、１送信データ期間（１Ｔ）
における帯域制限拡散変調信号と同じである。
して図１１に図示されたような帯域制限ＰＮ符号を発生
する。この帯域制限ＰＮ符号は、送信側において送信デ
ータを拡散変調する際に用いたＰＮ符号と同じ符号を帯
域制限して得たものであり、１送信データ期間（１Ｔ）
における帯域制限拡散変調信号と同じである。
【００２９】次に、第２フーリエ変換部１４は、参照信
号をフーリエ変換して図１３に示されるような信号波形
の周波数領域参照信号に変換する。
号をフーリエ変換して図１３に示されるような信号波形
の周波数領域参照信号に変換する。
【００３０】次いで、第２フィルタ手段１５は、周波数
領域参照信号の周波数帯域を図１３に図示された帯域幅
Ａの範囲になるように制限し、除算部１６に供給する。
この場合、、帯域幅Ａと帯域幅Ｂとは同じ帯域幅で、と
もに（１／Ｔｃ）になるように選ばれる。
領域参照信号の周波数帯域を図１３に図示された帯域幅
Ａの範囲になるように制限し、除算部１６に供給する。
この場合、、帯域幅Ａと帯域幅Ｂとは同じ帯域幅で、と
もに（１／Ｔｃ）になるように選ばれる。
【００３１】除算部１６は、周波数領域受信信号に対し
て周波数領域参照信号を基準とした比を各周波数毎に計
算し、図１５の点線で示される曲線ａのような雑音成分
を含んだ正規化信号を発生する。
て周波数領域参照信号を基準とした比を各周波数毎に計
算し、図１５の点線で示される曲線ａのような雑音成分
を含んだ正規化信号を発生する。
【００３２】次に、雑音処理除去部１７は、周波数窓乗
算部２３で雑音成分を含んだ正規化信号に周波数窓信号
を乗算し、フーリエ変換部２４で周波数窓信号を乗算し
た正規化信号をフーリエ変換して時間領域信号に変換
し、雑音除去部２５で時間領域信号中の雑音成分を除去
し、逆フーリエ変換部２６で時間領域信号をフーリエ変
換して周波数領域の正規化信号に再変換し、周波数窓除
算部２７で正規化信号を周波数窓信号で除算して、図１
５の実線で示される曲線ｂのような雑音成分を含まない
正規化信号Ｘ1 を形成し、正規化信号出力端子１９から
出力される。
算部２３で雑音成分を含んだ正規化信号に周波数窓信号
を乗算し、フーリエ変換部２４で周波数窓信号を乗算し
た正規化信号をフーリエ変換して時間領域信号に変換
し、雑音除去部２５で時間領域信号中の雑音成分を除去
し、逆フーリエ変換部２６で時間領域信号をフーリエ変
換して周波数領域の正規化信号に再変換し、周波数窓除
算部２７で正規化信号を周波数窓信号で除算して、図１
５の実線で示される曲線ｂのような雑音成分を含まない
正規化信号Ｘ1 を形成し、正規化信号出力端子１９から
出力される。
【００３３】以上の説明は、正規化信号発生部２1 の動
作に関するものであるが、他の正規化信号発生部２2 乃
至２4 においても同様の動作が行なわれ、それぞれ雑音
成分を含まない正規化信号Ｘ2 乃至Ｘ4 が形成される。
作に関するものであるが、他の正規化信号発生部２2 乃
至２4 においても同様の動作が行なわれ、それぞれ雑音
成分を含まない正規化信号Ｘ2 乃至Ｘ4 が形成される。
【００３４】続いて、図９において、４個の自己相関行
列発生部３1 乃至３4 は、対応する４つの正規化信号発
生部２1 乃至２4 が発生した正規化信号Ｘ1 乃至Ｘ4 を
各別に受け、対応する自己相関行列を発生する。一例と
して、正規化信号発生部２1は、正規化信号Ｘ1 と正規
化信号Ｘ1 の複素共役信号との積を周波数成分毎に計算
する等の信号処理を行って自己相関行列を発生する。他
の正規化信号発生部２2 乃至２4 も同様にして自己相関
行列を発生する。
列発生部３1 乃至３4 は、対応する４つの正規化信号発
生部２1 乃至２4 が発生した正規化信号Ｘ1 乃至Ｘ4 を
各別に受け、対応する自己相関行列を発生する。一例と
して、正規化信号発生部２1は、正規化信号Ｘ1 と正規
化信号Ｘ1 の複素共役信号との積を周波数成分毎に計算
する等の信号処理を行って自己相関行列を発生する。他
の正規化信号発生部２2 乃至２4 も同様にして自己相関
行列を発生する。
【００３５】平均化自己相関行列発生部４は、４個の自
己相関行列発生部３1 乃至３4 がそれぞれ発生する４個
の自己相関行列を、自己相関行列の要素毎に平均をとっ
て求めた平均化自己相関行列を発生し、高分解能信号処
理部５に供給する。
己相関行列発生部３1 乃至３4 がそれぞれ発生する４個
の自己相関行列を、自己相関行列の要素毎に平均をとっ
て求めた平均化自己相関行列を発生し、高分解能信号処
理部５に供給する。
【００３６】高分解能信号処理部５は、平均化自己相関
行列をＭＵＳＩＣ法を用いて高い時間分解能で解析し、
図１６に図示されるような相関計量信号を形成し、得ら
れた相関計量信号に基づいて相関計量値が所定値以上あ
るピーク計量値の数を計数して到来信号数を定める。ま
た、そのピーク計量値をとる際の遅延時間を計測し、各
到来信号遅延時間を定めて、到来信号数及び各到来信号
遅延時間を到来信号発生部７に供給する。図１６に図示
された例の場合は、到来信号数（ピーク計量値の数）は
２で、各到来信号遅延時間は第１波が１６．０５Ｔｃ、
第２波が１６．３６Ｔｃであり、第１波に対する第２波
の相対信号遅延時間は０．３１Ｔｃである。
行列をＭＵＳＩＣ法を用いて高い時間分解能で解析し、
図１６に図示されるような相関計量信号を形成し、得ら
れた相関計量信号に基づいて相関計量値が所定値以上あ
るピーク計量値の数を計数して到来信号数を定める。ま
た、そのピーク計量値をとる際の遅延時間を計測し、各
到来信号遅延時間を定めて、到来信号数及び各到来信号
遅延時間を到来信号発生部７に供給する。図１６に図示
された例の場合は、到来信号数（ピーク計量値の数）は
２で、各到来信号遅延時間は第１波が１６．０５Ｔｃ、
第２波が１６．３６Ｔｃであり、第１波に対する第２波
の相対信号遅延時間は０．３１Ｔｃである。
【００３７】また、相互相関行列発生部６は、４個の正
規化信号発生部２1 乃至２4 が発生した正規化信号Ｘ1
乃至Ｘ4 を受け、４個の正規化信号Ｘ1 乃至Ｘ4 の中の
それぞれ異なる第１及び第２の正規化信号の６つの組み
合わせに対して、第１の正規化信号と第２の正規化信号
の複素共役信号との積を周波数成分毎に計算して得た６
つの相互相関行列を発生する。ここで、４個の正規化信
号Ｘ1 乃至Ｘ4 の任意の２つの組合せの数は６になるの
で、６個の相互相関行列が発生される。
規化信号発生部２1 乃至２4 が発生した正規化信号Ｘ1
乃至Ｘ4 を受け、４個の正規化信号Ｘ1 乃至Ｘ4 の中の
それぞれ異なる第１及び第２の正規化信号の６つの組み
合わせに対して、第１の正規化信号と第２の正規化信号
の複素共役信号との積を周波数成分毎に計算して得た６
つの相互相関行列を発生する。ここで、４個の正規化信
号Ｘ1 乃至Ｘ4 の任意の２つの組合せの数は６になるの
で、６個の相互相関行列が発生される。
【００３８】到来信号発生部７は、相互相関行列発生部
６が発生した６つの相互相関行列及び高分解能信号処理
部５が発生する到来信号数（Ｋ個）及びＫ個の到来信号
遅延時間を受け、相互相関行列単位に到来信号数Ｋ個分
の到来信号を計算し、Ｋ個の到来信号のそれぞれについ
て６通りの到来信号を発生し、到来角推定部８に供給す
る。この場合、６通りの到来信号には、到来信号の振幅
の情報と、到来信号の到来天頂角及び到来方位角の情報
とがそれぞれ到来信号数Ｋ個分だけ含まれている。な
お、図１６に図示されている例の場合には、Ｋ＝２とな
る。
６が発生した６つの相互相関行列及び高分解能信号処理
部５が発生する到来信号数（Ｋ個）及びＫ個の到来信号
遅延時間を受け、相互相関行列単位に到来信号数Ｋ個分
の到来信号を計算し、Ｋ個の到来信号のそれぞれについ
て６通りの到来信号を発生し、到来角推定部８に供給す
る。この場合、６通りの到来信号には、到来信号の振幅
の情報と、到来信号の到来天頂角及び到来方位角の情報
とがそれぞれ到来信号数Ｋ個分だけ含まれている。な
お、図１６に図示されている例の場合には、Ｋ＝２とな
る。
【００３９】到来角推定部８は、到来信号発生部７がＫ
個の到来信号のそれぞれについて発生した６通りの到来
信号を受け、６通りの到来信号の位相を相互に比較し、
到来信号毎にその到来方向（到来天頂角及び到来方位
角）を計算し、計算結果を信号出力端子９から出力す
る。この場合、図１６に図示された例のように、到来信
号数が２個であると推定されたとき（Ｋ＝２）に、この
２個の到来信号の到来角を推定すれば、第１波の到来方
位角は５０．７°になり、第２波の到来方位角は６．２
°になる。
個の到来信号のそれぞれについて発生した６通りの到来
信号を受け、６通りの到来信号の位相を相互に比較し、
到来信号毎にその到来方向（到来天頂角及び到来方位
角）を計算し、計算結果を信号出力端子９から出力す
る。この場合、図１６に図示された例のように、到来信
号数が２個であると推定されたとき（Ｋ＝２）に、この
２個の到来信号の到来角を推定すれば、第１波の到来方
位角は５０．７°になり、第２波の到来方位角は６．２
°になる。
【００４０】このように、既知の電波到来方向推定装置
は、第１波、第２波等の到来信号を抽出するとともに、
それら到来信号の到来方向を求めることができる。
は、第１波、第２波等の到来信号を抽出するとともに、
それら到来信号の到来方向を求めることができる。
【００４１】
【発明が解決しようとする課題】前記既知の電波到来方
向推定装置は、時間的に近接して到来する信号成分を抽
出できるだけでなく、その到来方向を比較的正確に知る
ことができるものであるが、正規化信号発生部２1 乃至
２4 が発生する正規化信号Ｘ1 乃至Ｘ4 に雑音成分が含
まれているため、このような正規化信号Ｘ1 乃至Ｘ4 に
基づいて高分解能信号処理部５で到来信号遅延時間を求
めた場合、正規化信号Ｘ1 乃至Ｘ4 中の雑音成分の影響
を受けて、到来信号遅延時間にオフセット成分が加わ
り、到来信号発生部７において常時正確な到来信号を発
生させることができなくなる。
向推定装置は、時間的に近接して到来する信号成分を抽
出できるだけでなく、その到来方向を比較的正確に知る
ことができるものであるが、正規化信号発生部２1 乃至
２4 が発生する正規化信号Ｘ1 乃至Ｘ4 に雑音成分が含
まれているため、このような正規化信号Ｘ1 乃至Ｘ4 に
基づいて高分解能信号処理部５で到来信号遅延時間を求
めた場合、正規化信号Ｘ1 乃至Ｘ4 中の雑音成分の影響
を受けて、到来信号遅延時間にオフセット成分が加わ
り、到来信号発生部７において常時正確な到来信号を発
生させることができなくなる。
【００４２】その結果、前記既知の電波到来方向推定装
置は、主信号成分の到来方向を正しく求めることができ
ない場合が生じたり、主信号成分以外の他の信号成分の
到来方向を正しく求めることができない場合が生じると
いう問題を有している。
置は、主信号成分の到来方向を正しく求めることができ
ない場合が生じたり、主信号成分以外の他の信号成分の
到来方向を正しく求めることができない場合が生じると
いう問題を有している。
【００４３】最も、前記既知の電波到来方向推定装置
は、正規化信号発生部２1 乃至２4 内に雑音除去処理部
１７を配置し、正規化信号Ｘ1 乃至Ｘ4 中の雑音成分を
雑音除去処理部１７によって除去するようにしている
が、このような雑音除去処理部１７を配置しても、ＭＵ
ＳＩＣ法のような手法を用いる高分解能信号処理によっ
て得られる信号遅延時間には、オフセットが含まれると
いう本質的な問題がある。
は、正規化信号発生部２1 乃至２4 内に雑音除去処理部
１７を配置し、正規化信号Ｘ1 乃至Ｘ4 中の雑音成分を
雑音除去処理部１７によって除去するようにしている
が、このような雑音除去処理部１７を配置しても、ＭＵ
ＳＩＣ法のような手法を用いる高分解能信号処理によっ
て得られる信号遅延時間には、オフセットが含まれると
いう本質的な問題がある。
【００４４】従って、到来信号発生部７において常時正
確な到来信号を発生させることができない。
確な到来信号を発生させることができない。
【００４５】その結果、前記既知の電波到来方向推定装
置は、依然として、雑音除去処理部１７を配置しても、
主信号成分の到来方向を正しく求めることができない場
合が生じたり、主信号成分以外の他の信号成分の到来方
向を正しく求めることができない場合が生じるという問
題を有している。
置は、依然として、雑音除去処理部１７を配置しても、
主信号成分の到来方向を正しく求めることができない場
合が生じたり、主信号成分以外の他の信号成分の到来方
向を正しく求めることができない場合が生じるという問
題を有している。
【００４６】例えば、先に示した具体例の場合におい
て、第２波の到来方位角は本来０°であるのに対して、
前記既知の電波到来方向推定装置で推定した第２波の到
来方位角が６．２°になり、大きな方位角誤差が導入さ
れることになる。
て、第２波の到来方位角は本来０°であるのに対して、
前記既知の電波到来方向推定装置で推定した第２波の到
来方位角が６．２°になり、大きな方位角誤差が導入さ
れることになる。
【００４７】本発明は、このような問題点を解決するも
ので、その目的は、到来信号遅延時間に加わるオフセッ
ト成分をなくし、到来角推定部において常時正確な到来
信号を発生させ、その到来方向を正確に推定することを
可能にした電波到来方向推定方法及びその装置を提供す
ることにある。
ので、その目的は、到来信号遅延時間に加わるオフセッ
ト成分をなくし、到来角推定部において常時正確な到来
信号を発生させ、その到来方向を正確に推定することを
可能にした電波到来方向推定方法及びその装置を提供す
ることにある。
【００４８】
【課題を解決するための手段】前記目的を達成するため
に、本発明による電波到来方向推定方法は、Ｎ（３以上
の整数）個の正規化信号からＮ個の自己相関行列及びＭ
（２以上の整数）個の相互相関行列を発生し、Ｎ個の自
己相関行列を平均化した平均化自己相関行列を発生し、
平均化自己相関行列を高い時間分解能で信号分離処理し
て到来信号数（Ｋ個）及びＫ個の到来信号遅延時間を発
生し、Ｍ個の相互相関行列と到来信号数とＫ個の補正到
来信号遅延時間とを用いて、Ｋ個の到来信号を各Ｍ通り
発生し、各Ｍ通り発生したＫ個の到来信号の位相に基づ
いてＫ個の到来信号の到来方向を推定するもので、Ｋ個
の補正到来信号遅延時間は、各Ｍ通り発生したＫ個の到
来信号の振幅に基づいた誤差評価を行い、誤差が最小に
なるようにＫ個の到来信号遅延時間を補正したものであ
る第１の手段を具備する。
に、本発明による電波到来方向推定方法は、Ｎ（３以上
の整数）個の正規化信号からＮ個の自己相関行列及びＭ
（２以上の整数）個の相互相関行列を発生し、Ｎ個の自
己相関行列を平均化した平均化自己相関行列を発生し、
平均化自己相関行列を高い時間分解能で信号分離処理し
て到来信号数（Ｋ個）及びＫ個の到来信号遅延時間を発
生し、Ｍ個の相互相関行列と到来信号数とＫ個の補正到
来信号遅延時間とを用いて、Ｋ個の到来信号を各Ｍ通り
発生し、各Ｍ通り発生したＫ個の到来信号の位相に基づ
いてＫ個の到来信号の到来方向を推定するもので、Ｋ個
の補正到来信号遅延時間は、各Ｍ通り発生したＫ個の到
来信号の振幅に基づいた誤差評価を行い、誤差が最小に
なるようにＫ個の到来信号遅延時間を補正したものであ
る第１の手段を具備する。
【００４９】また、前記目的を達成するために、本発明
による電波到来方向推定装置は、Ｎ個の正規化信号を発
生するＮ個の正規化信号発生部と、Ｎ個の正規化信号か
らそれぞれ自己相関行列を発生するＮ個の自己相関行列
発生部と、Ｎ個の自己相関行列の平均化自己相関行列を
発生する平均化自己相関行列発生部と、平均化自己相関
行列を高い時間分解能で信号分離処理して到来信号数
（Ｋ個）及びＫ個の到来信号遅延時間を発生する高分解
能信号処理部と、Ｎ個の正規化信号からＭ個の相互相関
行列を発生する相互相関行列発生部と、Ｍ個の相互相関
行列と到来信号数とＫ個の補正到来信号遅延時間とを用
いてＫ個の到来信号を各Ｍ通り発生する到来信号発生部
と、高分解能信号処理部から出力されるＫ個の到来信号
遅延時間を、各Ｍ通り発生したＫ個の到来信号の振幅に
基づいた誤差評価を行い、誤差が最小になるようにＫ個
の到来信号遅延時間を補正したＫ個の補正到来信号遅延
時間を到来信号発生部に供給する遅延時間補正部と、各
Ｍ通り発生したＫ個の到来信号の位相に基づいてＫ個の
到来信号の到来方向を推定する到来角推定部とを備えた
第２の手段を具備する。
による電波到来方向推定装置は、Ｎ個の正規化信号を発
生するＮ個の正規化信号発生部と、Ｎ個の正規化信号か
らそれぞれ自己相関行列を発生するＮ個の自己相関行列
発生部と、Ｎ個の自己相関行列の平均化自己相関行列を
発生する平均化自己相関行列発生部と、平均化自己相関
行列を高い時間分解能で信号分離処理して到来信号数
（Ｋ個）及びＫ個の到来信号遅延時間を発生する高分解
能信号処理部と、Ｎ個の正規化信号からＭ個の相互相関
行列を発生する相互相関行列発生部と、Ｍ個の相互相関
行列と到来信号数とＫ個の補正到来信号遅延時間とを用
いてＫ個の到来信号を各Ｍ通り発生する到来信号発生部
と、高分解能信号処理部から出力されるＫ個の到来信号
遅延時間を、各Ｍ通り発生したＫ個の到来信号の振幅に
基づいた誤差評価を行い、誤差が最小になるようにＫ個
の到来信号遅延時間を補正したＫ個の補正到来信号遅延
時間を到来信号発生部に供給する遅延時間補正部と、各
Ｍ通り発生したＫ個の到来信号の位相に基づいてＫ個の
到来信号の到来方向を推定する到来角推定部とを備えた
第２の手段を具備する。
【００５０】前記第１の手段及び第２の手段によれば、
高分解能信号処理部から出力されたＫ個の到来信号数及
びＫ個の到来信号遅延時間を到来信号発生部に供給し、
到来信号発生部においてＫ個の到来信号をＭ通り発生す
る際に、各Ｍ通りのＫ個の到来信号の振幅に基づいて誤
差評価を行い、誤差が最小になるようにＫ個の到来信号
遅延時間を補正し、各補正到来信号遅延時間に変換した
後で到来信号発生部に供給し、到来信号発生部において
は、この補正到来信号遅延時間を用いて相互相関行列単
位毎にそれぞれＫ個の到来信号をＭ通り発生させてい
る。このため、正規化信号発生部から出力されるＮ個の
正規化信号中の雑音成分に起因した到来信号遅延時間の
オフセットや高分解能信号処理部から出力されるＫ個の
到来信号遅延時間のオフセットは、Ｋ個の到来信号遅延
時間をＫ個の補正到来信号遅延時間に補正することによ
り実質的になくなり、その結果、到来信号発生部からＭ
通り発生されるＫ個の到来信号中に誤差成分が含まれる
ことがなく、到来角推定部において、常時、Ｋ個の到来
信号の到来方向を正確に推定することができる。
高分解能信号処理部から出力されたＫ個の到来信号数及
びＫ個の到来信号遅延時間を到来信号発生部に供給し、
到来信号発生部においてＫ個の到来信号をＭ通り発生す
る際に、各Ｍ通りのＫ個の到来信号の振幅に基づいて誤
差評価を行い、誤差が最小になるようにＫ個の到来信号
遅延時間を補正し、各補正到来信号遅延時間に変換した
後で到来信号発生部に供給し、到来信号発生部において
は、この補正到来信号遅延時間を用いて相互相関行列単
位毎にそれぞれＫ個の到来信号をＭ通り発生させてい
る。このため、正規化信号発生部から出力されるＮ個の
正規化信号中の雑音成分に起因した到来信号遅延時間の
オフセットや高分解能信号処理部から出力されるＫ個の
到来信号遅延時間のオフセットは、Ｋ個の到来信号遅延
時間をＫ個の補正到来信号遅延時間に補正することによ
り実質的になくなり、その結果、到来信号発生部からＭ
通り発生されるＫ個の到来信号中に誤差成分が含まれる
ことがなく、到来角推定部において、常時、Ｋ個の到来
信号の到来方向を正確に推定することができる。
【００５１】
【発明の実施の形態】本発明の第１の実施の形態におい
て、電波到来方向推定方法は、近接配置した３以上の整
数（Ｎ）個のアンテナで捉えた信号電波をＮ個の正規化
信号に変換し、Ｎ個の正規化信号からＮ個の自己相関行
列及び２以上の整数（Ｍ）個の相互相関行列を発生し、
Ｎ個の自己相関行列を平均化した平均化自己相関行列を
発生し、平均化自己相関行列を高い時間分解能で信号分
離処理して到来信号数（Ｋ個）及びＫ個の到来信号遅延
時間を発生し、Ｍ個の相互相関行列と到来信号数とＫ個
の補正到来信号遅延時間とを用いて、Ｋ個の到来信号を
各Ｍ通り発生し、各Ｍ通り発生したＫ個の到来信号の位
相に基づいてＫ個の到来信号の到来方向を推定するもの
で、Ｋ個の補正到来信号遅延時間は、各Ｍ通り発生した
のＫ個の到来信号の振幅に基づいた誤差評価を行い、誤
差が最小になるようにＫ個の到来信号遅延時間を補正し
たものである。
て、電波到来方向推定方法は、近接配置した３以上の整
数（Ｎ）個のアンテナで捉えた信号電波をＮ個の正規化
信号に変換し、Ｎ個の正規化信号からＮ個の自己相関行
列及び２以上の整数（Ｍ）個の相互相関行列を発生し、
Ｎ個の自己相関行列を平均化した平均化自己相関行列を
発生し、平均化自己相関行列を高い時間分解能で信号分
離処理して到来信号数（Ｋ個）及びＫ個の到来信号遅延
時間を発生し、Ｍ個の相互相関行列と到来信号数とＫ個
の補正到来信号遅延時間とを用いて、Ｋ個の到来信号を
各Ｍ通り発生し、各Ｍ通り発生したＫ個の到来信号の位
相に基づいてＫ個の到来信号の到来方向を推定するもの
で、Ｋ個の補正到来信号遅延時間は、各Ｍ通り発生した
のＫ個の到来信号の振幅に基づいた誤差評価を行い、誤
差が最小になるようにＫ個の到来信号遅延時間を補正し
たものである。
【００５２】本発明の第１の実施の形態の１つの例にお
いて、電波到来方向推定方法は、正規化信号が各信号電
波を周波数変換して発生したベースバンド信号と基準と
なる参照信号との比を求めて得たものである。
いて、電波到来方向推定方法は、正規化信号が各信号電
波を周波数変換して発生したベースバンド信号と基準と
なる参照信号との比を求めて得たものである。
【００５３】本発明の第１の実施の形態の他の１つの例
において、電波到来方向推定方法は、正規化信号がベー
スバンド信号と参照信号との比を周波数領域で求めて得
たものである。
において、電波到来方向推定方法は、正規化信号がベー
スバンド信号と参照信号との比を周波数領域で求めて得
たものである。
【００５４】本発明の第１の実施の形態のさらに他の１
つの例において、電波到来方向推定方法は、ベースバン
ド信号が所定の符号で拡散変調された信号を含み、参照
信号が所定の符号と高い相関を持つ信号である。
つの例において、電波到来方向推定方法は、ベースバン
ド信号が所定の符号で拡散変調された信号を含み、参照
信号が所定の符号と高い相関を持つ信号である。
【００５５】本発明の第１の実施の形態のさらに他の１
つの例において、電波到来方向推定方法は、
つの例において、電波到来方向推定方法は、
【００５６】
【数１】
【００５７】により求めているものである。
【００５８】本発明の第１の実施の形態のさらに他の１
つの例において、電波到来方向推定方法は、Ｍ個の相互
相関行列を、Ｎ個の正規化信号の組み合わせによって発
生させることが可能な最大｛Ｎ（Ｎ−１）／２｝個の相
互相関行列からＭ個を選択して得たもので、選択された
Ｍ個の相互相関行列が、Ｎ個のアンテナの中で比較的離
れて配置されている２つのアンテナの正規化信号に対応
した相互相関行列からなるものである。
つの例において、電波到来方向推定方法は、Ｍ個の相互
相関行列を、Ｎ個の正規化信号の組み合わせによって発
生させることが可能な最大｛Ｎ（Ｎ−１）／２｝個の相
互相関行列からＭ個を選択して得たもので、選択された
Ｍ個の相互相関行列が、Ｎ個のアンテナの中で比較的離
れて配置されている２つのアンテナの正規化信号に対応
した相互相関行列からなるものである。
【００５９】本発明の第２の実施の形態において、電波
到来方向推定装置は、近接配置した３以上の整数（Ｎ）
個のアンテナと、Ｎ個のアンテナで捉えた信号電波の受
信によりＮ個の正規化信号を発生するＮ個の正規化信号
発生部と、Ｎ個の正規化信号からそれぞれ自己相関行列
を発生するＮ個の自己相関行列発生部と、Ｎ個の自己相
関行列の平均化自己相関行列を発生する平均化自己相関
行列発生部と、平均化自己相関行列を高い時間分解能で
信号分離処理して到来信号数（Ｋ個）及びＫ個の到来信
号遅延時間を発生する高分解能信号処理部と、Ｎ個の正
規化信号から２以上の整数（Ｍ）個の相互相関行列を発
生する相互相関行列発生部と、Ｍ個の相互相関行列と到
来信号数とＫ個の補正到来信号遅延時間とを用いて、Ｋ
個の到来信号を各Ｍ通り発生する到来信号発生部と、高
分解能信号処理部から出力されるＫ個の到来信号遅延時
間を、各Ｍ通り発生したＫ個の到来信号の振幅に基づい
た誤差評価を行い、誤差が最小になるようにＫ個の到来
信号遅延時間を補正したＫ個の補正到来信号遅延時間を
到来信号発生部に供給する遅延時間補正部と、各Ｍ通り
発生したＫ個の到来信号の位相に基づいてＫ個の到来信
号の到来方向を推定する到来角推定部とを備えているも
のである。
到来方向推定装置は、近接配置した３以上の整数（Ｎ）
個のアンテナと、Ｎ個のアンテナで捉えた信号電波の受
信によりＮ個の正規化信号を発生するＮ個の正規化信号
発生部と、Ｎ個の正規化信号からそれぞれ自己相関行列
を発生するＮ個の自己相関行列発生部と、Ｎ個の自己相
関行列の平均化自己相関行列を発生する平均化自己相関
行列発生部と、平均化自己相関行列を高い時間分解能で
信号分離処理して到来信号数（Ｋ個）及びＫ個の到来信
号遅延時間を発生する高分解能信号処理部と、Ｎ個の正
規化信号から２以上の整数（Ｍ）個の相互相関行列を発
生する相互相関行列発生部と、Ｍ個の相互相関行列と到
来信号数とＫ個の補正到来信号遅延時間とを用いて、Ｋ
個の到来信号を各Ｍ通り発生する到来信号発生部と、高
分解能信号処理部から出力されるＫ個の到来信号遅延時
間を、各Ｍ通り発生したＫ個の到来信号の振幅に基づい
た誤差評価を行い、誤差が最小になるようにＫ個の到来
信号遅延時間を補正したＫ個の補正到来信号遅延時間を
到来信号発生部に供給する遅延時間補正部と、各Ｍ通り
発生したＫ個の到来信号の位相に基づいてＫ個の到来信
号の到来方向を推定する到来角推定部とを備えているも
のである。
【００６０】本発明の第２の実施の形態の１つの例にお
いて、電波到来方向推定装置は、Ｎ個の正規化信号発生
部が、それぞれ、アンテナで捉えた各信号電波の受信信
号を周波数変換してベースバンド信号を発生するベース
バンド信号発生手段と、基準となる参照信号を発生する
参照信号発生手段と、ベースバンド信号と参照信号との
比をとって正規化信号を発生する除算手段とを備えてい
るものである。
いて、電波到来方向推定装置は、Ｎ個の正規化信号発生
部が、それぞれ、アンテナで捉えた各信号電波の受信信
号を周波数変換してベースバンド信号を発生するベース
バンド信号発生手段と、基準となる参照信号を発生する
参照信号発生手段と、ベースバンド信号と参照信号との
比をとって正規化信号を発生する除算手段とを備えてい
るものである。
【００６１】本発明の第２の実施の形態の他の１つの例
において、電波到来方向推定装置は、Ｎ個の正規化信号
発生部が、それぞれ、除算手段でベースバンド信号と参
照信号との比を周波数領域でとって正規化信号を発生す
るものである。
において、電波到来方向推定装置は、Ｎ個の正規化信号
発生部が、それぞれ、除算手段でベースバンド信号と参
照信号との比を周波数領域でとって正規化信号を発生す
るものである。
【００６２】本発明の第２の実施の形態のさらに他の１
つの例において、電波到来方向推定装置は、ベースバン
ド信号発生手段が所定の符号で拡散変調された信号を含
むベースバンド信号を発生し、参照信号発生手段が所定
の符号と高い相関を持つ参照信号を発生するものであ
る。
つの例において、電波到来方向推定装置は、ベースバン
ド信号発生手段が所定の符号で拡散変調された信号を含
むベースバンド信号を発生し、参照信号発生手段が所定
の符号と高い相関を持つ参照信号を発生するものであ
る。
【００６３】本発明の第２の実施の形態のさらに他の１
つの例において、電波到来方向推定装置は、遅延時間補
正部が、
つの例において、電波到来方向推定装置は、遅延時間補
正部が、
【００６４】
【数１】
【００６５】により求めることにより、到来信号遅延時
間の補正を行い、各到来信号遅延時間の補正を行ってい
るものである。
間の補正を行い、各到来信号遅延時間の補正を行ってい
るものである。
【００６６】本発明の第２の実施の形態のさらに他の１
つの例において、電波到来方向推定装置は、相互相関行
列発生部を、Ｎ個の正規化信号の組み合わせによって発
生させることが可能な最大｛Ｎ（Ｎ−１）／２｝個の相
互相関行列からＭ個の相互相関行列を選択して得たもの
で、選択されたＭ個の相互相関行列が、Ｎ個のアンテナ
の中で比較的離れて配置されている２つのアンテナの正
規化信号に対応したものである。
つの例において、電波到来方向推定装置は、相互相関行
列発生部を、Ｎ個の正規化信号の組み合わせによって発
生させることが可能な最大｛Ｎ（Ｎ−１）／２｝個の相
互相関行列からＭ個の相互相関行列を選択して得たもの
で、選択されたＭ個の相互相関行列が、Ｎ個のアンテナ
の中で比較的離れて配置されている２つのアンテナの正
規化信号に対応したものである。
【００６７】これらの本発明の実施の形態によれば、高
分解能信号処理部から出力されたＫ個の到来信号数及び
Ｋ個の到来信号遅延時間を到来信号発生部に供給し、到
来信号発生部において相互相関行列単位毎にそれぞれＫ
個の到来信号をＭ通り発生する際に、Ｋ個の到来信号遅
延時間を、Ｍ通り発生するＫ個の到来信号の振幅に基づ
き、所定ステップに従って、それらの誤差が最小になる
ように順次補正してそれぞれ対応する補正到来信号遅延
時間に変換した後で到来信号発生部に供給し、到来信号
発生部においては、この補正到来信号遅延時間を用いて
Ｋ個の到来信号をＭ通り発生させるようにしているの
で、Ｎ個の正規化信号発生部から出力されるＮ個の正規
化信号中の雑音成分に起因した到来信号遅延時間のオフ
セットや高分解能信号処理部から出力される各到来信号
遅延時間のオフセットが、Ｋ個の到来信号遅延時間をＫ
個の補正到来信号遅延時間に補正することにより実質的
になくすことができる。その結果、到来信号発生部から
Ｍ通り発生されるＫ個の到来信号成分中に含まれる誤差
成分がなくなり、到来角推定部において、常時、Ｋ個の
到来信号の到来方向を正確に推定することができる。
分解能信号処理部から出力されたＫ個の到来信号数及び
Ｋ個の到来信号遅延時間を到来信号発生部に供給し、到
来信号発生部において相互相関行列単位毎にそれぞれＫ
個の到来信号をＭ通り発生する際に、Ｋ個の到来信号遅
延時間を、Ｍ通り発生するＫ個の到来信号の振幅に基づ
き、所定ステップに従って、それらの誤差が最小になる
ように順次補正してそれぞれ対応する補正到来信号遅延
時間に変換した後で到来信号発生部に供給し、到来信号
発生部においては、この補正到来信号遅延時間を用いて
Ｋ個の到来信号をＭ通り発生させるようにしているの
で、Ｎ個の正規化信号発生部から出力されるＮ個の正規
化信号中の雑音成分に起因した到来信号遅延時間のオフ
セットや高分解能信号処理部から出力される各到来信号
遅延時間のオフセットが、Ｋ個の到来信号遅延時間をＫ
個の補正到来信号遅延時間に補正することにより実質的
になくすことができる。その結果、到来信号発生部から
Ｍ通り発生されるＫ個の到来信号成分中に含まれる誤差
成分がなくなり、到来角推定部において、常時、Ｋ個の
到来信号の到来方向を正確に推定することができる。
【００６８】
【実施例】以下、本発明の実施例を図面を参照して説明
する。
する。
【００６９】図１は、本発明に係わる電波到来方向推定
装置の第１実施例の構成を示すブロック図である。
装置の第１実施例の構成を示すブロック図である。
【００７０】図１に示されるように、第１実施例の電波
到来方向推定装置と、図９に示されている既知の電波到
来方向推定装置とを比べた場合、既知の電波到来方向推
定装置に用いられていない遅延時間補正部２６を用いて
いる点において構成が異なっているだけであって、その
他に、第１実施例の電波到来方向推定装置と既知の電波
到来方向推定装置との間に構成上の違いはない。
到来方向推定装置と、図９に示されている既知の電波到
来方向推定装置とを比べた場合、既知の電波到来方向推
定装置に用いられていない遅延時間補正部２６を用いて
いる点において構成が異なっているだけであって、その
他に、第１実施例の電波到来方向推定装置と既知の電波
到来方向推定装置との間に構成上の違いはない。
【００７１】この場合、本実施例の電波到来方向推定装
置に用いられている遅延時間補正部２６は、２個の入力
及び６個の入力と、１個の出力とを有し、２個の入力が
高分解能信号処理部５の到来信号数出力と各信号遅延時
間出力にそれぞれ接続され、６個の入力が到来信号発生
部７の６個の出力にそれぞれ接続され、１個の出力が到
来信号発生部７の遅延時間入力に接続される。
置に用いられている遅延時間補正部２６は、２個の入力
及び６個の入力と、１個の出力とを有し、２個の入力が
高分解能信号処理部５の到来信号数出力と各信号遅延時
間出力にそれぞれ接続され、６個の入力が到来信号発生
部７の６個の出力にそれぞれ接続され、１個の出力が到
来信号発生部７の遅延時間入力に接続される。
【００７２】なお、図１においては、図９に示された構
成要素と同じ構成要素については同じ符号を付け、図１
に示された本実施例の電波到来方向推定装置の構成につ
いては、前の説明と重複するので、これ以上の詳細な説
明を省略する。
成要素と同じ構成要素については同じ符号を付け、図１
に示された本実施例の電波到来方向推定装置の構成につ
いては、前の説明と重複するので、これ以上の詳細な説
明を省略する。
【００７３】また、図２は、図１に図示された各正規化
信号発生部２1 乃至２4 の構成の一例を示すブロック図
であって、図９に図示された既知の電波到来方向推定装
置における各正規化信号発生部２1 乃至２4 と同じ構成
のものである。
信号発生部２1 乃至２4 の構成の一例を示すブロック図
であって、図９に図示された既知の電波到来方向推定装
置における各正規化信号発生部２1 乃至２4 と同じ構成
のものである。
【００７４】図２に示された本実施例の各正規化信号発
生部２1 乃至２4 の構成についても、前の説明と重複す
るので、これ以上の詳細な説明を省略する。
生部２1 乃至２4 の構成についても、前の説明と重複す
るので、これ以上の詳細な説明を省略する。
【００７５】次に、図３は、本実施例の電波到来方向推
定装置の遅延時間補正部２６で実行される、到来電波の
第１波及び第２波の到来信号遅延時間の各変化に対応し
た累積正規化２乗誤差の変化の一例を示す等高線図であ
り、平面方向の２軸がＰＮ符号の１チップ（Ｔｃ）を時
間単位とした第１波及び第２波の到来信号遅延時間を示
し、高さ方向の軸が累積正規化２乗誤差を示すものであ
る。
定装置の遅延時間補正部２６で実行される、到来電波の
第１波及び第２波の到来信号遅延時間の各変化に対応し
た累積正規化２乗誤差の変化の一例を示す等高線図であ
り、平面方向の２軸がＰＮ符号の１チップ（Ｔｃ）を時
間単位とした第１波及び第２波の到来信号遅延時間を示
し、高さ方向の軸が累積正規化２乗誤差を示すものであ
る。
【００７６】前記構成による第１実施例の電波到来方向
推定装置の動作を説明する。
推定装置の動作を説明する。
【００７７】ただし、この動作説明において、既知の電
波到来方向推定装置の動作と同じ動作が行なわれる部分
については、説明が重複するので、その説明を省略す
る。
波到来方向推定装置の動作と同じ動作が行なわれる部分
については、説明が重複するので、その説明を省略す
る。
【００７８】本実施例の電波到来方向推定装置におい
て、送信側から送信される信号電波は、既知の電波到来
方向推定装置で用いられる信号電波と同じである。
て、送信側から送信される信号電波は、既知の電波到来
方向推定装置で用いられる信号電波と同じである。
【００７９】また、本実施例の電波到来方向推定装置に
おける４個のアンテナ１1 乃至１4の配置は、図１０に
示されるように、既知の電波到来方向推定装置における
４個のアンテナ１1 乃至１4 の配置と同じである。そし
て、本実施例の電波到来方向推定装置においても、第１
波は、４個のアンテナ１1 乃至１4 の配置面に平行に、
かつ、Ｘ軸に対して５０°の方位角方向から到来し、第
２波は、４個のアンテナ１1 乃至１4 の配置面に平行
に、かつ、Ｘ軸に対して０°の方位角方向から到来し、
第１波と第２波の電力が等しく、第１波に対して第２波
の相対的に０．２５チップ（Ｔｃ）だけ遅延しているも
のとする。
おける４個のアンテナ１1 乃至１4の配置は、図１０に
示されるように、既知の電波到来方向推定装置における
４個のアンテナ１1 乃至１4 の配置と同じである。そし
て、本実施例の電波到来方向推定装置においても、第１
波は、４個のアンテナ１1 乃至１4 の配置面に平行に、
かつ、Ｘ軸に対して５０°の方位角方向から到来し、第
２波は、４個のアンテナ１1 乃至１4 の配置面に平行
に、かつ、Ｘ軸に対して０°の方位角方向から到来し、
第１波と第２波の電力が等しく、第１波に対して第２波
の相対的に０．２５チップ（Ｔｃ）だけ遅延しているも
のとする。
【００８０】さらに、４個の正規化信号発生部２1 乃至
２4 の動作、４個の自己相関行列発生部３1 乃至３4 の
動作、平均化自己相関行列発生部４の動作、高分解能信
号処理部５の動作、相互相関行列発生部６の動作、到来
信号発生部７の動作、それに到来角推定部８の動作は、
いずれも、既知の電波到来方向推定装置における４個の
正規化信号発生部２1 乃至２4 の動作、４個の自己相関
行列発生部３1 乃至３4 の動作、平均化自己相関行列発
生部４の動作、高分解能信号処理部５の動作、相互相関
行列発生部６の動作、到来信号発生部７の動作、それに
到来角推定部８の動作と同じである。
２4 の動作、４個の自己相関行列発生部３1 乃至３4 の
動作、平均化自己相関行列発生部４の動作、高分解能信
号処理部５の動作、相互相関行列発生部６の動作、到来
信号発生部７の動作、それに到来角推定部８の動作は、
いずれも、既知の電波到来方向推定装置における４個の
正規化信号発生部２1 乃至２4 の動作、４個の自己相関
行列発生部３1 乃至３4 の動作、平均化自己相関行列発
生部４の動作、高分解能信号処理部５の動作、相互相関
行列発生部６の動作、到来信号発生部７の動作、それに
到来角推定部８の動作と同じである。
【００８１】高分解能信号処理部５は、前述のようにＭ
ＵＳＩＣ法を用いて信号処理を行うもので、到来信号数
と各信号遅延時間を出力し、相互相関行列発生部６は、
６個の相互相関行列を出力する。また、到来信号発生部
７は、６通りの到来信号を発生し、到来角推定部８は、
６通りの到来信号の位相を相互に比較し、最も遅延時間
の短い第１波及びそれに続く第２波等を抽出し、抽出し
た第１波及び第２波等の到来方向（到来天頂角及び到来
方位角）を計算し、計算結果を信号出力端子９から出力
する。
ＵＳＩＣ法を用いて信号処理を行うもので、到来信号数
と各信号遅延時間を出力し、相互相関行列発生部６は、
６個の相互相関行列を出力する。また、到来信号発生部
７は、６通りの到来信号を発生し、到来角推定部８は、
６通りの到来信号の位相を相互に比較し、最も遅延時間
の短い第１波及びそれに続く第２波等を抽出し、抽出し
た第１波及び第２波等の到来方向（到来天頂角及び到来
方位角）を計算し、計算結果を信号出力端子９から出力
する。
【００８２】第１実施例の電波到来方向推定装置は、使
用しているアンテナ１1 乃至１4 が４個であるので、相
互相関行列発生部６から出力される相互相関行列の組み
合わせは全部で６通りである。このとき、遅延時間補正
部２６は、各到来信号遅延時間を、次式（１）に従って
その累積正規化２乗誤差が最小になるように修正し、補
正到来信号遅延時間を発生し、到来信号発生部７に供給
する。
用しているアンテナ１1 乃至１4 が４個であるので、相
互相関行列発生部６から出力される相互相関行列の組み
合わせは全部で６通りである。このとき、遅延時間補正
部２６は、各到来信号遅延時間を、次式（１）に従って
その累積正規化２乗誤差が最小になるように修正し、補
正到来信号遅延時間を発生し、到来信号発生部７に供給
する。
【００８３】
【数２】
【００８４】また、アンテナ１1 乃至１4 が４個（Ｎ＝
４）である場合、（ｉ、ｊ）の組み合わせは（２、
１）、（３、１）、（４、１）、（３、２）、（４、
２）、（４、３）の６通り（Ｍ＝６）である。
４）である場合、（ｉ、ｊ）の組み合わせは（２、
１）、（３、１）、（４、１）、（３、２）、（４、
２）、（４、３）の６通り（Ｍ＝６）である。
【００８５】第１実施例の電波到来方向推定装置は、図
９に図示された既知の電波到来方向推定装置に遅延時間
補正部２６を付加した構成のものであり、遅延時間補正
部２６において各到来信号遅延時間を各補正到来信号遅
延時間に修正しているので、既知の電波到来方向推定装
置に比べて各到来信号の到来方向の推定精度を高めるこ
とができる。
９に図示された既知の電波到来方向推定装置に遅延時間
補正部２６を付加した構成のものであり、遅延時間補正
部２６において各到来信号遅延時間を各補正到来信号遅
延時間に修正しているので、既知の電波到来方向推定装
置に比べて各到来信号の到来方向の推定精度を高めるこ
とができる。
【００８６】ここで、到来信号遅延時間を修正して補正
到来信号遅延時間を得る状態を、図３に図示された等高
線図を用いて説明する。
到来信号遅延時間を得る状態を、図３に図示された等高
線図を用いて説明する。
【００８７】なお、図３において、各格子の一辺は、遅
延時間を調整する際の時間ステップであり、本実施例に
おいては時間ステップを０．０１２５チップ（Ｔｃ）に
選んでいる。
延時間を調整する際の時間ステップであり、本実施例に
おいては時間ステップを０．０１２５チップ（Ｔｃ）に
選んでいる。
【００８８】最初に、高分解能信号処理部５から出力さ
れた第１波の到来信号遅延時間は１６．０５チップ（Ｔ
ｃ）、第２波の到来信号遅延時間は１６．３６チップ
（Ｔｃ）であるので、累積正規化２乗誤差は図３の黒丸
の位置にある。このとき、その位置の累積正規化２乗誤
差は最小になっていない。
れた第１波の到来信号遅延時間は１６．０５チップ（Ｔ
ｃ）、第２波の到来信号遅延時間は１６．３６チップ
（Ｔｃ）であるので、累積正規化２乗誤差は図３の黒丸
の位置にある。このとき、その位置の累積正規化２乗誤
差は最小になっていない。
【００８９】そこで、遅延時間補正部２６は、第１波及
び第２波の到来信号遅延時間を前記時間ステップ単位で
変化させ、変化させる毎に得られた到来信号遅延時間を
補正到来信号遅延時間として到来信号発生部７に供給す
る。このとき、遅延時間補正部２６は、到来信号発生部
７から１つの到来信号当たり６通り出力された到来信号
を用い、再び、前記式（１）に従った誤差評価を行う。
即ち、第１波及び第２波の信号遅延時間を、累積正規化
２乗誤差が最小になる方向に前記時間ステップ単位で変
化させる。
び第２波の到来信号遅延時間を前記時間ステップ単位で
変化させ、変化させる毎に得られた到来信号遅延時間を
補正到来信号遅延時間として到来信号発生部７に供給す
る。このとき、遅延時間補正部２６は、到来信号発生部
７から１つの到来信号当たり６通り出力された到来信号
を用い、再び、前記式（１）に従った誤差評価を行う。
即ち、第１波及び第２波の信号遅延時間を、累積正規化
２乗誤差が最小になる方向に前記時間ステップ単位で変
化させる。
【００９０】このような一連の動作は、到来信号の振幅
に関する累積正規化２乗誤差が最小になるまで続行さ
れ、最終的に累積正規化２乗誤差が図３の白丸の位置に
きたときに、到来信号の振幅に関する累積正規化２乗誤
差が最小になったものと判断する。このとき、第１波の
信号遅延時間は１６．１２５チップ（Ｔｃ）、第２波の
信号遅延時間は１６．３７５チップ（Ｔｃ）であると推
定され、この白丸の位置における第１波及び第２波の各
信号遅延時間の推定誤差はともにゼロである。また、第
１波に対する第２波の相対的な信号遅延時間は０．２５
チップ（Ｔｃ）であって、初期設定した時間に一致して
いる。さらに、各到来信号遅延時間を遅延時間補正部２
６で補正した補正到来信号遅延時間を用いた場合、信号
振幅誤差は、第１波において＋０．１５ｄＢであり、第
２波において−０．０５ｄＢであるので、殆んど誤差を
有していない。
に関する累積正規化２乗誤差が最小になるまで続行さ
れ、最終的に累積正規化２乗誤差が図３の白丸の位置に
きたときに、到来信号の振幅に関する累積正規化２乗誤
差が最小になったものと判断する。このとき、第１波の
信号遅延時間は１６．１２５チップ（Ｔｃ）、第２波の
信号遅延時間は１６．３７５チップ（Ｔｃ）であると推
定され、この白丸の位置における第１波及び第２波の各
信号遅延時間の推定誤差はともにゼロである。また、第
１波に対する第２波の相対的な信号遅延時間は０．２５
チップ（Ｔｃ）であって、初期設定した時間に一致して
いる。さらに、各到来信号遅延時間を遅延時間補正部２
６で補正した補正到来信号遅延時間を用いた場合、信号
振幅誤差は、第１波において＋０．１５ｄＢであり、第
２波において−０．０５ｄＢであるので、殆んど誤差を
有していない。
【００９１】このように、第１実施例の電波到来方向推
定装置によれば、高分解能信号処理部５から出力された
到来信号遅延時間を遅延時間補正部２６において修正
し、到来信号発生部７においては補正到来信号遅延時間
を用いて到来信号毎に６通りの到来信号を形成し、到来
角推定部８においては１個の到来信号あたり６通りの到
来信号を用いて遅延時間が最も短い到来信号、即ち、第
１波の到来角の推定を行ったところ、図１０に示される
ように、第１波の本来の到来角５０°に対し、推定到来
角として４９．５°が得られ、その推定方位誤差が０．
５°になり、また、第２波の本来の到来角０°に対し、
推定到来角として０．３°が得られ、その推定方位誤差
が０．３°になり、既知の電波到来方向推定装置と比べ
て推定方位誤差を極めて少なくすることができるもので
ある。
定装置によれば、高分解能信号処理部５から出力された
到来信号遅延時間を遅延時間補正部２６において修正
し、到来信号発生部７においては補正到来信号遅延時間
を用いて到来信号毎に６通りの到来信号を形成し、到来
角推定部８においては１個の到来信号あたり６通りの到
来信号を用いて遅延時間が最も短い到来信号、即ち、第
１波の到来角の推定を行ったところ、図１０に示される
ように、第１波の本来の到来角５０°に対し、推定到来
角として４９．５°が得られ、その推定方位誤差が０．
５°になり、また、第２波の本来の到来角０°に対し、
推定到来角として０．３°が得られ、その推定方位誤差
が０．３°になり、既知の電波到来方向推定装置と比べ
て推定方位誤差を極めて少なくすることができるもので
ある。
【００９２】次に、図４は、本発明に係わる電波到来方
向推定装置の第２実施例の構成を示すブロック図であ
る。
向推定装置の第２実施例の構成を示すブロック図であ
る。
【００９３】図４に図示の第２実施例の電波到来方向推
定装置は、最大のアンテナ開口位置にある２つのアンテ
ナ、例えば、図１０に図示されたアンテナ１1 、１3 及
びアンテナ１2 、１4 の２通りの組み合わせに対応した
２つの相互相関行列を選択した例を示すものである。
定装置は、最大のアンテナ開口位置にある２つのアンテ
ナ、例えば、図１０に図示されたアンテナ１1 、１3 及
びアンテナ１2 、１4 の２通りの組み合わせに対応した
２つの相互相関行列を選択した例を示すものである。
【００９４】第２実施例の電波到来方向推定装置は、第
１実施例の電波到来方向推定装置において用いられてい
る相関相関行列発生部６、到来信号発生部７、到来角推
定部８、遅延時間補正部２６の代わりに、それぞれ相関
相関行列発生部６’、到来信号発生部７’、到来角推定
部８’、遅延時間補正部２６’を用いているものであ
る。この場合、相関相関行列発生部６と相関相関行列発
生部６’、到来信号発生部７と到来信号発生部７’、到
来角推定部８と到来角推定部８’、遅延時間補正部２６
と遅延時間補正部２６’は、それぞれ、同じ位置に同じ
ように接続され、ともに同じ働きをするものである。
１実施例の電波到来方向推定装置において用いられてい
る相関相関行列発生部６、到来信号発生部７、到来角推
定部８、遅延時間補正部２６の代わりに、それぞれ相関
相関行列発生部６’、到来信号発生部７’、到来角推定
部８’、遅延時間補正部２６’を用いているものであ
る。この場合、相関相関行列発生部６と相関相関行列発
生部６’、到来信号発生部７と到来信号発生部７’、到
来角推定部８と到来角推定部８’、遅延時間補正部２６
と遅延時間補正部２６’は、それぞれ、同じ位置に同じ
ように接続され、ともに同じ働きをするものである。
【００９５】なお、図４においては、図１に示された構
成要素と同じ構成要素については同じ符号を付け、図４
に示された第２実施例の電波到来方向推定装置の構成に
ついては、前述の第１実施例の電波到来方向推定装置の
構成と同じであってその説明と重複するので、これ以上
の詳細な説明を省略する。
成要素と同じ構成要素については同じ符号を付け、図４
に示された第２実施例の電波到来方向推定装置の構成に
ついては、前述の第１実施例の電波到来方向推定装置の
構成と同じであってその説明と重複するので、これ以上
の詳細な説明を省略する。
【００９６】第２実施例の電波到来方向推定装置は、最
大の開口位置にある２つのアンテナの組み合わせは、ア
ンテナ１1 とアンテナ１3 及びアンテナ１2 とアンテナ
１4の２通り（Ｍ＝２）である。このとき、遅延時間補
正部２６’は、各到来信号遅延時間を、次式（３）に従
ってその累積正規化２乗誤差が最小になるように修正
し、補正到来信号遅延時間を発生し、到来信号発生部
７’に供給する。
大の開口位置にある２つのアンテナの組み合わせは、ア
ンテナ１1 とアンテナ１3 及びアンテナ１2 とアンテナ
１4の２通り（Ｍ＝２）である。このとき、遅延時間補
正部２６’は、各到来信号遅延時間を、次式（３）に従
ってその累積正規化２乗誤差が最小になるように修正
し、補正到来信号遅延時間を発生し、到来信号発生部
７’に供給する。
【００９７】
【数３】
【００９８】なお、式（３）において、前記式（１）に
おける相互相関行列の数を６から２に変更した場合、式
（３）は式（１）と同じ形になる。
おける相互相関行列の数を６から２に変更した場合、式
（３）は式（１）と同じ形になる。
【００９９】図５は、第２実施例において、到来信号遅
延時間を修正して補正到来信号遅延時間を得る状態を表
す等高線図である。
延時間を修正して補正到来信号遅延時間を得る状態を表
す等高線図である。
【０１００】図５に図示の等高線図を用い、到来信号遅
延時間を修正して補正到来信号遅延時間を得る状態を説
明する。
延時間を修正して補正到来信号遅延時間を得る状態を説
明する。
【０１０１】なお、図５においても、各格子の一辺は、
遅延時間を調整する際の時間ステップであり、本実施例
においては時間ステップを０．０１２５チップ（Ｔｃ）
に選んでいる。
遅延時間を調整する際の時間ステップであり、本実施例
においては時間ステップを０．０１２５チップ（Ｔｃ）
に選んでいる。
【０１０２】始めに、高分解能信号処理部５から出力さ
れた第１波の到来信号遅延時間は１６．０５チップ（Ｔ
ｃ）、第２波の到来信号遅延時間は１６．３６チップ
（Ｔｃ）であるので、累積正規化２乗誤差は図５の黒丸
の位置にある。このときに、その位置の累積正規化２乗
誤差は最小になっていない。
れた第１波の到来信号遅延時間は１６．０５チップ（Ｔ
ｃ）、第２波の到来信号遅延時間は１６．３６チップ
（Ｔｃ）であるので、累積正規化２乗誤差は図５の黒丸
の位置にある。このときに、その位置の累積正規化２乗
誤差は最小になっていない。
【０１０３】遅延時間補正部２６’は、第１波及び第２
波の到来信号遅延時間を前記時間ステップ単位で変化さ
せ、変化させる毎に得られた到来信号遅延時間を補正到
来信号遅延時間として到来信号発生部７’に供給する。
このとき、遅延時間補正部２６’は、到来信号発生部
７’から１つの到来信号当たり２通り出力された到来信
号を用い、再び、前記式（３）に従った誤差評価を行
う。即ち、第１波及び第２波の信号遅延時間を、累積正
規化２乗誤差が最小になる方向に前記時間ステップ単位
で変化させる。
波の到来信号遅延時間を前記時間ステップ単位で変化さ
せ、変化させる毎に得られた到来信号遅延時間を補正到
来信号遅延時間として到来信号発生部７’に供給する。
このとき、遅延時間補正部２６’は、到来信号発生部
７’から１つの到来信号当たり２通り出力された到来信
号を用い、再び、前記式（３）に従った誤差評価を行
う。即ち、第１波及び第２波の信号遅延時間を、累積正
規化２乗誤差が最小になる方向に前記時間ステップ単位
で変化させる。
【０１０４】このような一連の動作は、到来信号の振幅
に関する累積正規化２乗誤差が最小になるまで続行さ
れ、最終的に累積正規化２乗誤差が図５の白丸の位置に
きたときに、到来信号の振幅に関する累積正規化２乗誤
差が最小になったものと判断する。このとき、第１波の
信号遅延時間は１６．１２５チップ（Ｔｃ）、第２波の
信号遅延時間は１６．３７５チップ（Ｔｃ）であると推
定され、この白丸の位置における第１波及び第２波の各
信号遅延時間の推定誤差はともにゼロである。また、第
１波に対する第２波の相対的な信号遅延時間は０．２５
チップ（Ｔｃ）であって、初期設定した時間に一致して
いる。さらに、各到来信号遅延時間を遅延時間補正部２
６’で補正した補正到来信号遅延時間を用いた場合、信
号振幅誤差は、第１波において＋０．１０ｄＢであり、
第２波において−０．１０ｄＢであるので、殆んど誤差
がない。
に関する累積正規化２乗誤差が最小になるまで続行さ
れ、最終的に累積正規化２乗誤差が図５の白丸の位置に
きたときに、到来信号の振幅に関する累積正規化２乗誤
差が最小になったものと判断する。このとき、第１波の
信号遅延時間は１６．１２５チップ（Ｔｃ）、第２波の
信号遅延時間は１６．３７５チップ（Ｔｃ）であると推
定され、この白丸の位置における第１波及び第２波の各
信号遅延時間の推定誤差はともにゼロである。また、第
１波に対する第２波の相対的な信号遅延時間は０．２５
チップ（Ｔｃ）であって、初期設定した時間に一致して
いる。さらに、各到来信号遅延時間を遅延時間補正部２
６’で補正した補正到来信号遅延時間を用いた場合、信
号振幅誤差は、第１波において＋０．１０ｄＢであり、
第２波において−０．１０ｄＢであるので、殆んど誤差
がない。
【０１０５】ここで、図６は、第２実施例における遅延
時間補正部２６’及び既知の電波到来方向推定装置にお
ける高分解能信号処理部４で得られる信号遅延時間の推
定誤差の変化状態の一例を示す特性図であって、第２波
の到来方位角を０°に固定し、第１波の到来方位角を０
°から１８０°まで変化させた際に得られた特性を示す
ものである。
時間補正部２６’及び既知の電波到来方向推定装置にお
ける高分解能信号処理部４で得られる信号遅延時間の推
定誤差の変化状態の一例を示す特性図であって、第２波
の到来方位角を０°に固定し、第１波の到来方位角を０
°から１８０°まで変化させた際に得られた特性を示す
ものである。
【０１０６】また、図７は、第２実施例及び既知の電波
到来方向推定装置において、到来信号発生部７’、７で
得られる信号振幅誤差の変化状態の一例を示す特性図で
あり、第２波の到来方位角を０°に固定し、第１波の到
来方位角を０°から１８０°まで変化させた際に得られ
た特性を示すものである。
到来方向推定装置において、到来信号発生部７’、７で
得られる信号振幅誤差の変化状態の一例を示す特性図で
あり、第２波の到来方位角を０°に固定し、第１波の到
来方位角を０°から１８０°まで変化させた際に得られ
た特性を示すものである。
【０１０７】さらに、図８は、本実施例及び既知の電波
到来方向推定装置において、到来角推定部８’、８で得
られる到来方位角の推定誤差の変化状態の一例を示す特
性図であリ、第２波の到来方位角を０°に固定し、第１
波の到来方位角を０°から１８０°まで変化させた際に
得られた特性を示すものである。
到来方向推定装置において、到来角推定部８’、８で得
られる到来方位角の推定誤差の変化状態の一例を示す特
性図であリ、第２波の到来方位角を０°に固定し、第１
波の到来方位角を０°から１８０°まで変化させた際に
得られた特性を示すものである。
【０１０８】図６乃至図８において、白菱形点を結んだ
曲線ａは、第２実施例の電波到来方向推定装置における
特性であり、黒菱形点を結んだ曲線ｂは、既知の電波到
来方向推定装置における特性である。
曲線ａは、第２実施例の電波到来方向推定装置における
特性であり、黒菱形点を結んだ曲線ｂは、既知の電波到
来方向推定装置における特性である。
【０１０９】図６乃至図８を用いて、第２実施例の電波
到来方向推定装置と既知の電波到来方向推定装置との各
種特性の比較を行う。
到来方向推定装置と既知の電波到来方向推定装置との各
種特性の比較を行う。
【０１１０】まず、図６に図示されるように、第２波の
到来角を０°に固定し、第１波の到来角を０°から１８
０°まで変化させたときに得られる信号遅延時間誤差
は、既知の電波到来方向推定装置の最大誤差が０．１９
Ｔｃ（曲線ｂ参照）であるのに対し、第２実施例の電波
到来方向推定装置の最大誤差が０．０５Ｔｃ（曲線ａ参
照）になり、本実施例の電波到来方向推定装置の信号遅
延時間誤差は大幅に低減されている。
到来角を０°に固定し、第１波の到来角を０°から１８
０°まで変化させたときに得られる信号遅延時間誤差
は、既知の電波到来方向推定装置の最大誤差が０．１９
Ｔｃ（曲線ｂ参照）であるのに対し、第２実施例の電波
到来方向推定装置の最大誤差が０．０５Ｔｃ（曲線ａ参
照）になり、本実施例の電波到来方向推定装置の信号遅
延時間誤差は大幅に低減されている。
【０１１１】次に、図７に図示されるように、第２波の
到来角を０°に固定し、第１波の到来角を０°から１８
０°まで変化させたときに得られる信号振幅誤差は、既
知の電波到来方向推定装置の最大誤差が１８．２ｄＢ
（曲線ｂ参照）であるのに対し、第２実施例の電波到来
方向推定装置の最大誤差が４．４ｄＢ（曲線ａ参照）に
なり、本実施例の電波到来方向推定装置の信号振幅誤差
は大幅に低減されている。
到来角を０°に固定し、第１波の到来角を０°から１８
０°まで変化させたときに得られる信号振幅誤差は、既
知の電波到来方向推定装置の最大誤差が１８．２ｄＢ
（曲線ｂ参照）であるのに対し、第２実施例の電波到来
方向推定装置の最大誤差が４．４ｄＢ（曲線ａ参照）に
なり、本実施例の電波到来方向推定装置の信号振幅誤差
は大幅に低減されている。
【０１１２】次いで、図８に示されるように、第２波の
到来角を０°に固定し、第１波の到来角を０°から１８
０°まで変化させたときに得られる到来方位角の誤差
は、既知の電波到来方向推定装置の最大誤差が１５３．
０°（曲線ｂ参照）であるのに対し、第２実施例の電波
到来方向推定装置の最大誤差が３．２°（曲線ａ参照）
になり、本実施例の電波到来方向推定装置の到来方位角
の誤差についても、大幅に低減されている。
到来角を０°に固定し、第１波の到来角を０°から１８
０°まで変化させたときに得られる到来方位角の誤差
は、既知の電波到来方向推定装置の最大誤差が１５３．
０°（曲線ｂ参照）であるのに対し、第２実施例の電波
到来方向推定装置の最大誤差が３．２°（曲線ａ参照）
になり、本実施例の電波到来方向推定装置の到来方位角
の誤差についても、大幅に低減されている。
【０１１３】このように、第２実施例の電波到来方向推
定装置によれば、高分解能信号処理部５から出力された
到来信号遅延時間を遅延時間補正部２６’で補正到来信
号遅延時間に修正し、到来信号発生部７’はこの補正到
来信号遅延時間を用いて相互相関行列単位で到来信号毎
に最大の開口位置にある２つのアンテナの組み合わせに
よる２通りの到来信号を形成し、到来角推定部８’にお
いては２通りの到来信号を用いて遅延時間が最も短い到
来信号、即ち、第１波の到来角の推定を行ったところ、
図１０に図示されるように、第１波の本来の到来角５０
°に対し、推定到来角として４９．５°が得られてその
推定方位誤差が０．５°になり、また、第２波の本来の
到来角０°に対し、推定到来角として０．３°が得られ
てその推定方位誤差が０．３°になり、第１実施例と同
様の推定方位誤差を極めて少なくすることができる。
定装置によれば、高分解能信号処理部５から出力された
到来信号遅延時間を遅延時間補正部２６’で補正到来信
号遅延時間に修正し、到来信号発生部７’はこの補正到
来信号遅延時間を用いて相互相関行列単位で到来信号毎
に最大の開口位置にある２つのアンテナの組み合わせに
よる２通りの到来信号を形成し、到来角推定部８’にお
いては２通りの到来信号を用いて遅延時間が最も短い到
来信号、即ち、第１波の到来角の推定を行ったところ、
図１０に図示されるように、第１波の本来の到来角５０
°に対し、推定到来角として４９．５°が得られてその
推定方位誤差が０．５°になり、また、第２波の本来の
到来角０°に対し、推定到来角として０．３°が得られ
てその推定方位誤差が０．３°になり、第１実施例と同
様の推定方位誤差を極めて少なくすることができる。
【０１１４】また、第２実施例の電波到来方向推定装置
によれば、第１実施例の電波到来方向推定装置と同様の
推定方位誤差で各信号電波の到来方向を推定することが
できるので、第１実施例の電波到来方向推定装置に比べ
て相互相関行列の数を低減することができ、その分、各
部の計算量を減らすことができる。
によれば、第１実施例の電波到来方向推定装置と同様の
推定方位誤差で各信号電波の到来方向を推定することが
できるので、第１実施例の電波到来方向推定装置に比べ
て相互相関行列の数を低減することができ、その分、各
部の計算量を減らすことができる。
【０１１５】なお、前記各実施例においては、信号電波
として第１波及び第２波からなる２波が到来した例を挙
げて説明したが、本発明の電波到来方向推定装置におい
て到来角の推定処理可能な信号電波の数は２つである場
合に限られるものでなく、３つまたはそれ以上の信号電
波が到来した場合においても、同様に動作させることが
可能である。
として第１波及び第２波からなる２波が到来した例を挙
げて説明したが、本発明の電波到来方向推定装置におい
て到来角の推定処理可能な信号電波の数は２つである場
合に限られるものでなく、３つまたはそれ以上の信号電
波が到来した場合においても、同様に動作させることが
可能である。
【０１１６】また、前記各実施例においては、４個のア
ンテナ１1 乃至１4 を用い、それに関連して４個の正規
化信号発生部２’1 乃至２’4 及び４個の自己相関行列
発生部３1 乃至３4 を用いた例を挙げて説明したが、本
発明の電波到来方向推定装置において用いられるアンテ
ナ、正規化信号発生部、自己相関行列発生部の数は４個
である場合に限られるものでなく、複数であれば３個以
上の数を用いるようにしてもよい。
ンテナ１1 乃至１4 を用い、それに関連して４個の正規
化信号発生部２’1 乃至２’4 及び４個の自己相関行列
発生部３1 乃至３4 を用いた例を挙げて説明したが、本
発明の電波到来方向推定装置において用いられるアンテ
ナ、正規化信号発生部、自己相関行列発生部の数は４個
である場合に限られるものでなく、複数であれば３個以
上の数を用いるようにしてもよい。
【０１１７】さらに、前記各実施例においては、４個の
アンテナ１1 乃至１4 が円形に配置された例を挙げて説
明したが、本発明の電波到来方向推定装置におけるアン
テナの配置は円形である場合に限られるものでなく、格
子状の配置やそれ以外の形状の配置であってもよい。
アンテナ１1 乃至１4 が円形に配置された例を挙げて説
明したが、本発明の電波到来方向推定装置におけるアン
テナの配置は円形である場合に限られるものでなく、格
子状の配置やそれ以外の形状の配置であってもよい。
【０１１８】この他に、前記各実施例においては、各正
規化信号発生部２1 乃至２4 が図２に図示される構成の
ものを用いた例を挙げて説明したが、本発明の電波到来
方向推定装置に用いられる各正規化信号発生部は、図２
に図示されるような構成のものに限られるものでなく、
その構成に類似の構成を有し、類似の動作を行うもので
あれば、他の構成のものを用いるようにしてもよい。
規化信号発生部２1 乃至２4 が図２に図示される構成の
ものを用いた例を挙げて説明したが、本発明の電波到来
方向推定装置に用いられる各正規化信号発生部は、図２
に図示されるような構成のものに限られるものでなく、
その構成に類似の構成を有し、類似の動作を行うもので
あれば、他の構成のものを用いるようにしてもよい。
【０１１９】また、前記各実施例においては、高分解能
信号処理部５で行われる信号処理にＭＵＳＩＣ法を用い
た例を挙げて説明したが、本発明の電波到来方向推定装
置に用いられる高分解能信号処理部５は、ＭＵＳＩＣ法
を用いたものに限られるものでなく、高い時間的信号分
解能を有するものであれば、他の方法を用いたものであ
ってもよい。
信号処理部５で行われる信号処理にＭＵＳＩＣ法を用い
た例を挙げて説明したが、本発明の電波到来方向推定装
置に用いられる高分解能信号処理部５は、ＭＵＳＩＣ法
を用いたものに限られるものでなく、高い時間的信号分
解能を有するものであれば、他の方法を用いたものであ
ってもよい。
【０１２０】
【発明の効果】以上のように、本発明によれば、高分解
能信号処理部から出力された到来信号数（Ｋ個）及びＫ
個の到来信号遅延時間を到来信号発生部に供給し、到来
信号発生部においてＫ個の到来信号をそれぞれＭ通り発
生する際に、各Ｍ通り発生したＫ個の到来信号の振幅に
基づき、所定ステップに従って、それらの累積正規化２
乗誤差が最小になるようにＫ個の到来信号遅延時間を順
次補正してＫ個の補正到来信号遅延時間に変換した後で
到来信号発生部に供給し、到来信号発生部においては、
この補正到来信号遅延時間を用いてＫ個の到来信号をそ
れぞれＭ通り発生させるようにしているので、Ｎ個の正
規化信号発生部から出力されるＮ個の正規化信号中の雑
音成分に起因した到来信号遅延時間のオフセットや高分
解能信号処理部から出力されるＫ個の到来信号遅延時間
のオフセットを、各到来信号遅延時間を各補正到来信号
遅延時間に補正することにより実質的になくすことがで
き、その結果、到来信号発生部からそれぞれＭ通り発生
されるＫ個の到来信号に含まれる誤差成分が除去され、
到来角推定部において、常時、Ｋ個の到来信号の到来方
向を正確に推定できるという効果がある。
能信号処理部から出力された到来信号数（Ｋ個）及びＫ
個の到来信号遅延時間を到来信号発生部に供給し、到来
信号発生部においてＫ個の到来信号をそれぞれＭ通り発
生する際に、各Ｍ通り発生したＫ個の到来信号の振幅に
基づき、所定ステップに従って、それらの累積正規化２
乗誤差が最小になるようにＫ個の到来信号遅延時間を順
次補正してＫ個の補正到来信号遅延時間に変換した後で
到来信号発生部に供給し、到来信号発生部においては、
この補正到来信号遅延時間を用いてＫ個の到来信号をそ
れぞれＭ通り発生させるようにしているので、Ｎ個の正
規化信号発生部から出力されるＮ個の正規化信号中の雑
音成分に起因した到来信号遅延時間のオフセットや高分
解能信号処理部から出力されるＫ個の到来信号遅延時間
のオフセットを、各到来信号遅延時間を各補正到来信号
遅延時間に補正することにより実質的になくすことがで
き、その結果、到来信号発生部からそれぞれＭ通り発生
されるＫ個の到来信号に含まれる誤差成分が除去され、
到来角推定部において、常時、Ｋ個の到来信号の到来方
向を正確に推定できるという効果がある。
【図１】本発明による電波到来方向推定装置の第１実施
例の構成を示すブロック図である。
例の構成を示すブロック図である。
【図２】第１及び第２実施例における正規化信号発生部
の構成の一例を示すブロック図である。
の構成の一例を示すブロック図である。
【図３】第１実施例の電波到来方向推定装置に到来した
第１波及び第２波の信号電波の信号遅延時間に対する累
積正規化２乗誤差の一例を示す等高線図である。
第１波及び第２波の信号電波の信号遅延時間に対する累
積正規化２乗誤差の一例を示す等高線図である。
【図４】本発明による電波到来方向推定装置の第２実施
例の構成を示すブロック図である。
例の構成を示すブロック図である。
【図５】第２実施例の電波到来方向推定装置に到来した
第１波及び第２波の信号電波の信号遅延時間に対する累
積正規化２乗誤差の一例を示す等高線図である。
第１波及び第２波の信号電波の信号遅延時間に対する累
積正規化２乗誤差の一例を示す等高線図である。
【図６】第２実施例及び既知の電波到来方向推定装置に
おける第１波の到来方位角に対する遅延時間推定誤差の
変化状態の一例を示す特性図である。
おける第１波の到来方位角に対する遅延時間推定誤差の
変化状態の一例を示す特性図である。
【図７】第２実施例及び既知の電波到来方向推定装置に
おける第１波の到来方位角に対する信号振幅誤差の変化
状態の一例を示す特性図である。
おける第１波の到来方位角に対する信号振幅誤差の変化
状態の一例を示す特性図である。
【図８】第２実施例及び既知の電波到来方向推定装置に
おける第１波の到来方位角に対する到来方位角の推定誤
差の変化状態の一例を示す特性図である。
おける第１波の到来方位角に対する到来方位角の推定誤
差の変化状態の一例を示す特性図である。
【図９】既知の電波到来方向推定装置の構成の一例を示
すブロック図である。
すブロック図である。
【図１０】４つのアンテナの配置状態と信号電波の到来
状態の一例を示す構成図である。
状態の一例を示す構成図である。
【図１１】送信側でＰＮ符号により拡散変調されさらに
帯域制限して得た帯域制限拡散変調信号及び参照信号発
生部が発生する参照信号の一例を示す信号波形図であ
る。
帯域制限して得た帯域制限拡散変調信号及び参照信号発
生部が発生する参照信号の一例を示す信号波形図であ
る。
【図１２】ベースバンド（ＢＢ）拡散変調信号の一例を
示す信号波形図である。
示す信号波形図である。
【図１３】周波数領域参照信号の一例を示す信号波形図
である。
である。
【図１４】周波数領域受信信号の一例を示す信号波形図
である。
である。
【図１５】雑音成分を含んだ正規化信号及び雑音成分を
除去した正規化信号の一例を示す信号波形図である。
除去した正規化信号の一例を示す信号波形図である。
【図１６】高分解能信号処理部から出力される遅延時間
を含んだ相関計量信号の一例を示す信号波形図である。
を含んだ相関計量信号の一例を示す信号波形図である。
１1 〜１4 アンテナ ２1 〜２4 正規化信号発生部 ３1 〜３4 自己相関行列発生部 ４ 平均化自己相関行列発生部 ５ 高分解能信号処理部 ６、６’ 相互相関行列発生部 ７、７’ 到来信号発生部 ８、８’ 到来角推定部 ９ 信号出力端子 １０ 受信部 １１ 第１フーリエ変換部 １２ 第１フィルタ手段 １３ 参照信号発生部 １４ 第２フーリエ変換部 １５ 第２フィルタ手段 １６ 除算部 １７ 雑音除去処理部 １８ ベースバンド（ＢＢ）信号発生部 １９ アナログ−ディジタル（Ａ／Ｄ）変換部 ２０ メモリ ２１ フーリエ変換部 ２２ 雑音除去部 ２３ 逆フーリエ変換部 ２４ 受信信号入力端子 ２５ 正規化信号出力端子 ２６、２６’ 遅延時間補正部 ２７ 周波数窓乗算部 ２８ 周波数窓除算部
Claims (12)
- 【請求項１】 近接配置した３以上の整数（Ｎ）個のア
ンテナで捉えた信号電波をＮ個の正規化信号に変換し、
前記Ｎ個の正規化信号からＮ個の自己相関行列及び２以
上の整数（Ｍ）個の相互相関行列を発生し、前記Ｎ個の
自己相関行列を平均化した平均化自己相関行列を発生
し、前記平均化自己相関行列を高い時間分解能で信号分
離処理して到来信号数（Ｋ個）及びＫ個の到来信号遅延
時間を発生し、前記Ｍ個の相互相関行列と前記到来信号
数とＫ個の補正到来信号遅延時間とを用いて、Ｋ個の到
来信号を各Ｍ通り発生し、前記各Ｍ通り発生したＫ個の
到来信号の位相に基づいて前記Ｋ個の到来信号の到来方
向を推定するもので、前記Ｋ個の補正到来信号遅延時間
は、前記各Ｍ通り発生したＫ個の到来信号の振幅に基づ
いた誤差評価を行い、誤差が最小になるように前記Ｋ個
の到来信号遅延時間を補正したものであることを特徴と
する電波到来方向推定方法。 - 【請求項２】 前記Ｎ個の正規化信号は、各信号電波を
周波数変換して発生したベースバンド信号と、基準とな
る参照信号との比を求めて得たものであることを特徴と
する請求項１に記載の電波到来方向推定方法。 - 【請求項３】 前記Ｎ個の正規化信号は、前記ベースバ
ンド信号と前記参照信号との比を周波数領域で求めて得
たものであることを特徴とする請求項２に記載の電波到
来方向推定方法。 - 【請求項４】 前記ベースバンド信号は、所定の符号で
拡散変調された信号を含み、前記参照信号は、前記所定
の符号と高い相関を持つ信号であることを特徴とする請
求項２または３に記載の電波到来方向推定方法。 - 【請求項５】 【数１】 により求めることを特徴とする請求項１に記載の電波到
来方向推定方法。 - 【請求項６】 前記Ｍ個の相互相関行列は、前記Ｎ個の
正規化信号の組み合わせによって発生させることが可能
な最大｛Ｎ（Ｎ−１）／２｝個の相互相関行列からＭ個
が選択されたもので、前記選択されたＭ個の相互相関行
列は、前記Ｎ個のアンテナの中で比較的離れて配置され
ている２つのアンテナの正規化信号に対応した相互相関
行列であることを特徴とする請求項１に記載の電波到来
方向推定方法。 - 【請求項７】 近接配置した３以上の整数（Ｎ）個のア
ンテナと、前記Ｎ個のアンテナで捉えた信号電波の受信
によりＮ個の正規化信号を発生するＮ個の正規化信号発
生部と、前記Ｎ個の正規化信号からそれぞれ自己相関行
列を発生するＮ個の自己相関行列発生部と、前記Ｎ個の
自己相関行列の平均化自己相関行列を発生する平均化自
己相関行列発生部と、前記平均化自己相関行列を高い時
間分解能で信号分離処理して到来信号数（Ｋ個）及びＫ
個の到来信号遅延時間を発生する高分解能信号処理部
と、前記Ｎ個の正規化信号から２以上の整数（Ｍ）個の
相互相関行列を発生する相互相関行列発生部と、前記Ｍ
個の相互相関行列と前記到来信号数とＫ個の補正到来信
号遅延時間とを用いて、Ｋ個の到来信号を各Ｍ通り発生
する到来信号発生部と、前記高分解能信号処理部から出
力される前記Ｋ個の到来信号遅延時間を、前記各Ｍ通り
発生したＫ個の到来信号の振幅に基づいた誤差評価を行
い、誤差が最小になるように前記Ｋ個の到来信号遅延時
間を補正した前記Ｋ個の補正到来信号遅延時間を前記到
来信号発生部に供給する遅延時間補正部と、前記各Ｍ通
り発生したＫ個の到来信号の位相に基づいて前記Ｋ個の
到来信号の到来方向を推定する到来角推定部とを備えて
いることを特徴とする電波到来方向推定装置。 - 【請求項８】 前記Ｎ個の正規化信号発生部は、それぞ
れ、前記アンテナで捉えた各信号電波の受信信号を周波
数変換してベースバンド信号を発生するベースバンド信
号発生手段と、基準となる参照信号を発生する参照信号
発生手段と、前記ベースバンド信号と前記参照信号との
比をとって正規化信号を発生する除算手段とを備えてい
ることを特徴とする請求項７に記載の電波到来方向推定
装置。 - 【請求項９】 前記Ｎ個の正規化信号発生部は、それぞ
れ、前記除算手段で前記ベースバンド信号と前記参照信
号との比を周波数領域でとって正規化信号を発生するも
のであることを特徴とする請求項８に記載の電波到来方
向推定装置。 - 【請求項１０】 前記ベースバンド信号発生手段は、所
定の符号で拡散変調された信号を含むベースバンド信号
を発生し、前記参照信号発生手段は、前記所定の符号と
高い相関を持つ参照信号を発生することを特徴とする請
求項８または９に記載の電波到来方向推定装置。 - 【請求項１１】 前記遅延時間補正部は、 【数１】 により求めることにより、到来信号遅延時間の補正を行
うものであることを特徴とする請求項７に記載の電波到
来方向推定装置。 - 【請求項１２】 前記相互相関行列発生部は、前記Ｎ個
の正規化信号の組み合わせによって発生させることが可
能な最大｛Ｎ（Ｎ−１）／２｝個の相互相関行列からＭ
個の相互相関行列を選択するもので、前記選択されたＭ
個の相互相関行列は、前記Ｎ個のアンテナの中で比較的
離れて配置されている２つのアンテナの正規化信号に対
応したものであることを特徴とする請求項７に記載の電
波到来方向推定装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP28443398A JP3600459B2 (ja) | 1998-10-06 | 1998-10-06 | 電波到来方向推定方法及びその装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP28443398A JP3600459B2 (ja) | 1998-10-06 | 1998-10-06 | 電波到来方向推定方法及びその装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2000111630A true JP2000111630A (ja) | 2000-04-21 |
| JP3600459B2 JP3600459B2 (ja) | 2004-12-15 |
Family
ID=17678497
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP28443398A Expired - Fee Related JP3600459B2 (ja) | 1998-10-06 | 1998-10-06 | 電波到来方向推定方法及びその装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3600459B2 (ja) |
Cited By (23)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3424648B2 (ja) | 2000-04-05 | 2003-07-07 | 日本電気株式会社 | 電波環境分析装置及びコンピュータ読み取り可能な記憶媒体 |
| WO2004104620A1 (ja) * | 2003-05-22 | 2004-12-02 | Fujitsu Limited | 固有値分解を利用しない信号到来方向推定手法および受信ビーム形成装置 |
| JP2006242895A (ja) * | 2005-03-07 | 2006-09-14 | Toshiba Corp | 電波誘導装置 |
| KR100785056B1 (ko) | 2006-11-29 | 2007-12-12 | 국방과학연구소 | 상관 벡터를 이용한 방향 탐지 방법 |
| WO2008143142A1 (ja) * | 2007-05-18 | 2008-11-27 | Tokyo Electron Device Limited | 音源定位装置、音源定位方法、及び、プログラム |
| JP2009523648A (ja) * | 2006-01-19 | 2009-06-25 | テールズ | 編隊を組む宇宙船グループ内の一宇宙船用の、出力測定による相対位置制御装置 |
| JP2016511403A (ja) * | 2013-02-25 | 2016-04-14 | コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． | 微分到来角を用いた自律的方向探知 |
| KR101991844B1 (ko) * | 2018-06-20 | 2019-06-21 | 국방과학연구소 | 지연 시간 추정 장치 및 방법 |
| JP2020176902A (ja) * | 2019-04-17 | 2020-10-29 | 株式会社Ihi | 時間差測定装置及び到来方向推定装置 |
| US11451287B1 (en) * | 2021-03-16 | 2022-09-20 | Pivotal Commware, Inc. | Multipath filtering for wireless RF signals |
| US11497050B2 (en) | 2021-01-26 | 2022-11-08 | Pivotal Commware, Inc. | Smart repeater systems |
| US11563279B2 (en) | 2020-01-03 | 2023-01-24 | Pivotal Commware, Inc. | Dual polarization patch antenna system |
| US11670849B2 (en) | 2020-04-13 | 2023-06-06 | Pivotal Commware, Inc. | Aimable beam antenna system |
| US11706722B2 (en) | 2018-03-19 | 2023-07-18 | Pivotal Commware, Inc. | Communication of wireless signals through physical barriers |
| US11757180B2 (en) | 2019-02-20 | 2023-09-12 | Pivotal Commware, Inc. | Switchable patch antenna |
| US11843955B2 (en) | 2021-01-15 | 2023-12-12 | Pivotal Commware, Inc. | Installation of repeaters for a millimeter wave communications network |
| US11844050B2 (en) | 2020-09-08 | 2023-12-12 | Pivotal Commware, Inc. | Installation and activation of RF communication devices for wireless networks |
| US11848478B2 (en) | 2019-02-05 | 2023-12-19 | Pivotal Commware, Inc. | Thermal compensation for a holographic beam forming antenna |
| US11929822B2 (en) | 2021-07-07 | 2024-03-12 | Pivotal Commware, Inc. | Multipath repeater systems |
| US11937199B2 (en) | 2022-04-18 | 2024-03-19 | Pivotal Commware, Inc. | Time-division-duplex repeaters with global navigation satellite system timing recovery |
| US11968593B2 (en) | 2020-08-03 | 2024-04-23 | Pivotal Commware, Inc. | Wireless communication network management for user devices based on real time mapping |
| US11973568B2 (en) | 2020-05-27 | 2024-04-30 | Pivotal Commware, Inc. | RF signal repeater device management for 5G wireless networks |
| US12185453B2 (en) | 2021-10-26 | 2024-12-31 | Pivotal Commware, Inc. | RF absorbing structures |
-
1998
- 1998-10-06 JP JP28443398A patent/JP3600459B2/ja not_active Expired - Fee Related
Cited By (30)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3424648B2 (ja) | 2000-04-05 | 2003-07-07 | 日本電気株式会社 | 電波環境分析装置及びコンピュータ読み取り可能な記憶媒体 |
| WO2004104620A1 (ja) * | 2003-05-22 | 2004-12-02 | Fujitsu Limited | 固有値分解を利用しない信号到来方向推定手法および受信ビーム形成装置 |
| US7068221B2 (en) | 2003-05-22 | 2006-06-27 | Fujitsu Limited | Technique for direction-of-arrival estimation without eigendecomposition and its application to beamforming at base station |
| JP2006242895A (ja) * | 2005-03-07 | 2006-09-14 | Toshiba Corp | 電波誘導装置 |
| JP2009523648A (ja) * | 2006-01-19 | 2009-06-25 | テールズ | 編隊を組む宇宙船グループ内の一宇宙船用の、出力測定による相対位置制御装置 |
| KR100785056B1 (ko) | 2006-11-29 | 2007-12-12 | 국방과학연구소 | 상관 벡터를 이용한 방향 탐지 방법 |
| WO2008143142A1 (ja) * | 2007-05-18 | 2008-11-27 | Tokyo Electron Device Limited | 音源定位装置、音源定位方法、及び、プログラム |
| JP2016511403A (ja) * | 2013-02-25 | 2016-04-14 | コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． | 微分到来角を用いた自律的方向探知 |
| US12425987B2 (en) | 2018-03-19 | 2025-09-23 | Pivotal Commware, Inc. | Communication of wireless signals through physical barriers |
| US11706722B2 (en) | 2018-03-19 | 2023-07-18 | Pivotal Commware, Inc. | Communication of wireless signals through physical barriers |
| KR101991844B1 (ko) * | 2018-06-20 | 2019-06-21 | 국방과학연구소 | 지연 시간 추정 장치 및 방법 |
| US11848478B2 (en) | 2019-02-05 | 2023-12-19 | Pivotal Commware, Inc. | Thermal compensation for a holographic beam forming antenna |
| US12362472B2 (en) | 2019-02-20 | 2025-07-15 | Pivotal Commware, Inc. | Switchable patch antenna |
| US11757180B2 (en) | 2019-02-20 | 2023-09-12 | Pivotal Commware, Inc. | Switchable patch antenna |
| JP7267825B2 (ja) | 2019-04-17 | 2023-05-02 | 株式会社Ihi | 時間差測定装置及び到来方向推定装置 |
| JP2020176902A (ja) * | 2019-04-17 | 2020-10-29 | 株式会社Ihi | 時間差測定装置及び到来方向推定装置 |
| US11563279B2 (en) | 2020-01-03 | 2023-01-24 | Pivotal Commware, Inc. | Dual polarization patch antenna system |
| US11670849B2 (en) | 2020-04-13 | 2023-06-06 | Pivotal Commware, Inc. | Aimable beam antenna system |
| US11973568B2 (en) | 2020-05-27 | 2024-04-30 | Pivotal Commware, Inc. | RF signal repeater device management for 5G wireless networks |
| US11968593B2 (en) | 2020-08-03 | 2024-04-23 | Pivotal Commware, Inc. | Wireless communication network management for user devices based on real time mapping |
| US11844050B2 (en) | 2020-09-08 | 2023-12-12 | Pivotal Commware, Inc. | Installation and activation of RF communication devices for wireless networks |
| US11843955B2 (en) | 2021-01-15 | 2023-12-12 | Pivotal Commware, Inc. | Installation of repeaters for a millimeter wave communications network |
| US12010703B2 (en) | 2021-01-26 | 2024-06-11 | Pivotal Commware, Inc. | Smart repeater systems |
| US11497050B2 (en) | 2021-01-26 | 2022-11-08 | Pivotal Commware, Inc. | Smart repeater systems |
| US20220302992A1 (en) * | 2021-03-16 | 2022-09-22 | Pivotal Commware, Inc. | Multipath filtering for wireless rf signals |
| US11451287B1 (en) * | 2021-03-16 | 2022-09-20 | Pivotal Commware, Inc. | Multipath filtering for wireless RF signals |
| US11929822B2 (en) | 2021-07-07 | 2024-03-12 | Pivotal Commware, Inc. | Multipath repeater systems |
| US12185453B2 (en) | 2021-10-26 | 2024-12-31 | Pivotal Commware, Inc. | RF absorbing structures |
| US11937199B2 (en) | 2022-04-18 | 2024-03-19 | Pivotal Commware, Inc. | Time-division-duplex repeaters with global navigation satellite system timing recovery |
| US12495377B2 (en) | 2022-04-18 | 2025-12-09 | Pivotal Commware, Inc. | Time-division-duplex repeaters with global navigation satellite system timing recovery |
Also Published As
| Publication number | Publication date |
|---|---|
| JP3600459B2 (ja) | 2004-12-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3600459B2 (ja) | 電波到来方向推定方法及びその装置 | |
| US6529745B1 (en) | Radio wave arrival direction estimating antenna apparatus | |
| US6009334A (en) | Method and system for determining position of mobile radio terminals | |
| JP3673700B2 (ja) | スペクトル拡散信号を用いた測距及び位置測定方法、その方法を行う装置 | |
| US20040102157A1 (en) | Wireless LAN with distributed access points for space management | |
| JP2003004834A (ja) | 信号の到来方向を推定するための方法 | |
| US6931244B2 (en) | Radio equipment capable of real time change of antenna directivity and doppler frequency estimating circuit used for the radio equipment | |
| JPH11215543A (ja) | セルラ基地局およびそれに搭載される位置標定装置 | |
| KR20010062214A (ko) | 개선된 신호 포착과 프로세싱을 지닌 코드분할 다중접속시스템 및 오퍼레이션 방법 | |
| US6208842B1 (en) | Method and apparatus for estimating a channel parameter | |
| CN100444531C (zh) | 接收装置 | |
| JPH11340884A (ja) | ダイバーシチ受信下での信号合成方法及び装置 | |
| JP2973416B1 (ja) | Ｒａｋｅ受信回路 | |
| WO1998008319A1 (en) | Rake receiver for spread spectrum signal demodulation | |
| JP2000147083A (ja) | 到来角測定方法及び到来角測定器 | |
| JP2000286629A (ja) | 無線送信装置及び送信指向性調整方法 | |
| JP3630581B2 (ja) | 拡散変調信号受信装置 | |
| JP3639106B2 (ja) | 拡散変調信号受信装置 | |
| JP4193939B2 (ja) | アレーアンテナ、到来波推定装置、及び平面アレーの合成方法 | |
| JP2002204193A (ja) | 移動通信システム | |
| JP3639120B2 (ja) | 拡散変調信号受信装置 | |
| US6807223B2 (en) | Method of performing code synchronization, and receiver | |
| JP7638460B2 (ja) | 電波到来角推定装置、電波到来角推定方法、制御回路および記憶媒体 | |
| JPH11118898A (ja) | 電波方位計測システム | |
| JP2002359584A (ja) | 無線基地局装置及び指向性送信方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20040810 |
|
| 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: 20040831 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20040916 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20070924 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080924 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080924 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090924 Year of fee payment: 5 |
|
| LAPS | Cancellation because of no payment of annual fees |