(57)【要約】【課題】遅延時間の推定精度の劣化を招くことなく、相関行列の次元数Ｍdを小さくして、演算量の低減を図ることができるようにする。【解決手段】ディジタルの受信信号ｓ（１），ｓ（２），・・・，ｓ（Ｎd）を予め設定された複数の遅延時間推定区間＃ｌ（１≦ｌ≦Ｌ）の部分信号に分割し、複数の遅延時間推定区間＃ｌの部分信号にデシメーション処理を施すデシメーション処理部７を設け、遅延時間推定部８−１〜８−Ｌが複数のデシメーション処理信号ｘl（１），ｘl（２），・・・，ｘl（Ｍ）に超分解能遅延時間推定処理を施して、遅延時間推定区間＃ｌに存在している目標信号の遅延時間を推定する。【選択図】図１
2008107242.tif 000002
【特許請求の範囲】
【請求項１】
目標に向けて電波を送信する電波送信手段と、上記電波送信手段から送信されて、上記目標に反射された電波を受信して受信信号を生成する電波受信手段と、上記電波受信手段により生成された受信信号を予め設定された複数の遅延時間推定区間の部分信号に分割し、複数の遅延時間推定区間の部分信号にデシメーション処理を施すデシメーション手段と、上記デシメーション手段による複数のデシメーション処理信号に超分解能遅延時間推定処理を施して、上記遅延時間推定区間に存在している目標信号の遅延時間を推定する遅延時間推定手段とを備えた遅延時間推定装置。
【請求項２】
デシメーション手段は、電波受信手段により生成された受信信号を予め設定された複数の遅延時間推定区間の部分信号に分割する受信信号分割部と、上記受信信号分割部により分割された部分信号にフーリエ変換処理を施して受信信号スペクトルを生成する複数の受信信号スペクトル生成部と、上記受信信号スペクトル生成部により生成された受信信号スペクトルを送信波スペクトルで除算して除算信号を生成する複数の除算信号生成部と、上記除算信号生成部により生成された除算信号にデシメーション処理を施す複数のデシメーションフィルタとから構成されていることを特徴とする請求項１記載の遅延時間推定装置。
【請求項３】
デシメーション手段は、電波受信手段により生成された受信信号を予め設定された複数の遅延時間推定区間の部分信号に分割する受信信号分割部と、上記受信信号分割部により分割された部分信号にフーリエ変換処理を施して受信信号スペクトルを生成する複数の受信信号スペクトル生成部と、上記受信信号スペクトル生成部により生成された受信信号スペクトルを送信波スペクトルで除算して除算信号を生成する複数の除算信号生成部と、上記除算信号生成部により生成された除算信号に逆フーリエ変換処理を施して復調信号を生成する複数の逆フーリエ変換処理部と、上記逆フーリエ変換処理部により生成された復調信号から遅延時間推定区間の信号を切り出す複数の信号切出処理部と、上記信号切出処理部により切り出された遅延時間推定区間の信号にフーリエ変換処理を施す複数のフーリエ変換処理部とから構成されていることを特徴とする請求項１記載の遅延時間推定装置。
【請求項４】
除算信号生成部は、送信波スペクトルの振幅値が予め設定された最小値より小さい場合、上記送信波スペクトルを上記最小値に修正し、修正後の送信波スペクトルで受信信号スペクトルを除算することを特徴とする請求項２または請求項３記載の遅延時間推定装置。
【請求項５】
遅延時間推定手段は、デシメーション手段による複数のデシメーション処理信号間の相関を表す相関行列を生成する相関行列生成部と、上記相関行列生成部により生成された相関行列の固有値と固有ベクトルを算出し、所定の閾値より大きい固有値があれば、目標信号が存在するものと認定して上記固有ベクトルを出力する一方、所定の閾値より大きい固有値がなければ、目標信号が存在していないと認定する目標判別部と、上記目標判別部から出力された固有ベクトルを雑音空間とするＭＵＳＩＣ処理を実施して、目標信号の遅延時間を推定するＭＵＳＩＣ処理部とから構成されていることを特徴とする請求項１記載の遅延時間推定装置。
【請求項６】
目標に向けて電波を送信する電波送信手段と、上記電波送信手段から送信されて、上記目標に反射された電波を受信して受信信号を生成する電波受信手段と、上記電波受信手段により生成された受信信号を予め設定された複数の遅延時間推定区間の部分信号に分割して、複数の遅延時間推定区間の部分信号に目標信号が含まれているか否かを判定し、目標信号が含まれている部分信号にデシメーション処理を施すデシメーション手段と、上記デシメーション手段によるデシメーション処理信号に超分解能遅延時間推定処理を施して、目標信号の遅延時間を推定する遅延時間推定手段とを備えた遅延時間推定装置。
【請求項７】
デシメーション手段は、電波受信手段により生成された受信信号を予め設定された複数の遅延時間推定区間の部分信号に分割する受信信号分割部と、上記受信信号分割部により分割された部分信号にフーリエ変換処理を施して受信信号スペクトルを生成する複数の受信信号スペクトル生成部と、上記複数の受信信号スペクトル生成部により生成された受信信号スペクトルのうち、平均電力が基準電力より大きい受信信号スペクトルには目標信号が含まれていると判定して、その受信信号スペクトルを出力する目標信号検出部と、上記目標信号検出部から出力された受信信号スペクトルを送信波スペクトルで除算して除算信号を生成する除算信号生成部と、上記除算信号生成部により生成された除算信号にデシメーション処理を施すデシメーションフィルタとから構成されていることを特徴とする請求項６記載の遅延時間推定装置。
【請求項８】
電波受信手段を２個実装するとともに、上記２個の電波受信手段により生成された受信信号を最大比合成して、最大比合成後の受信信号をデシメーション手段に出力する最大比合成処理手段を設けたことを特徴とする請求項１から請求項７のうちのいずれか１項記載の遅延時間推定装置。
【請求項９】
電波受信手段は、電波送信手段から送信される電波がパルス単位で符号変調が施されている場合、パルスのチップ幅よりも短い時間間隔で電波をサンプリングして受信信号を生成することを特徴とする請求項１から請求項８のうちのいずれか１項記載の遅延時間推定  装置。
【請求項１０】
目標に向けて電波を送信する電波送信手段と、上記電波送信手段から送信されて、上記目標に反射された電波を受信して受信信号を生成する電波受信手段と、上記電波受信手段が受信信号を生成する毎に、その受信信号をパルス方向に並べて蓄積するとともに、その受信信号を時間方向に並べて蓄積する受信信号蓄積手段と、上記受信信号蓄積手段にパルス方向に並べて蓄積されている受信信号にフーリエ変換処理を施す複数のフーリエ変換処理手段と、上記複数のフーリエ変換処理手段によるフーリエ変換処理信号を解析して、目標信号が存在している遅延時間推定区間を検出する遅延時間推定区間検出手段と、上記受信信号蓄積手段に時間方向に並べて蓄積されている受信信号から上記遅延時間推定区間検出手段により検出された遅延時間推定区間の受信信号を抽出し、上記受信信号にデシメーション処理を施すデシメーション手段と、上記デシメーション手段によるデシメーション処理信号に超分解能遅延時間推定処理を施して、上記遅延時間推定区間に存在している目標信号の遅延時間を推定する遅延時間推定手段とを備えた遅延時間推定装置。
【請求項１１】
目標に向けて電波を送信する電波送信手段と、上記電波送信手段から送信されて、上記目標に反射された電波を受信して受信信号を生成する電波受信手段と、上記電波受信手段が受信信号を生成する毎に、その受信信号をパルス方向に並べて蓄積する受信信号蓄積手段と、上記受信信号蓄積手段にパルス方向に並べて蓄積されている受信信号にフーリエ変換処理を施す複数のフーリエ変換処理手段と、上記複数のフーリエ変換処理手段によるフーリエ変換処理信号を解析して、目標信号が存在している遅延時間推定区間を検出する遅延時間推定区間検出手段と、上記遅延時間推定区間検出手段により検出された遅延時間推定区間のフーリエ変換処理信号にデシメーション処理を施すデシメーション手段と、上記デシメーション手段によるデシメーション処理信号に超分解能遅延時間推定処理を施して、上記遅延時間推定区間に存在している目標信号の遅延時間を推定する遅延時間推定手段とを備えた遅延時間推定装置。
【請求項１２】
遅延時間推定手段は、デシメーション手段による複数のデシメーション処理信号間の相関を表す相関行列を生成する相関行列生成部と、上記相関行列生成部により生成された相関行列の固有値を算出し、その固有値に対応する固有ベクトルを算出する固有ベクトル算出部と、上記固有ベクトル算出部により算出された固有ベクトルを用いてＥＳＰＲＩＴ処理を実施して、目標信号の遅延時間を推定するＥＳＰＲＩＴ処理部とから構成されていることを特徴とする請求項１記載の遅延時間推定装置。
【請求項１３】
遅延時間推定手段は、デシメーション手段による複数のデシメーション処理信号間の相関を表す相関行列を生成する相関行列生成部と、上記相関行列生成部により生成された相関行列の固有値を算出し、その固有値に対応する固有ベクトルを算出する固有ベクトル算出部と、上記固有ベクトル算出部により算出された固有ベクトルを用いて最尤推定処理を実施して、目標信号の遅延時間を推定する最尤推定処理部とから構成されていることを特徴とする請求項１記載の遅延時間推定装置。
【請求項１４】
遅延時間推定手段は、デシメーション手段による複数のデシメーション処理信号を蓄積する信号蓄積部と、上記信号蓄積部に蓄積されているデシメーション処理信号のうち、隣接しているデシメーション処理信号間の相関値を基準にして、目標信号の有無を判定する有無判定部と、上記有無判定部の判定結果が目標信号が有る旨を示していれば、上記信号蓄積部に蓄積されている複数のデシメーション処理信号間の相関を表す相関行列を生成する相関行列生成部と、上記相関行列生成部により生成された相関行列の固有値を算出し、その固有値に対応する固有ベクトルを算出する固有ベクトル算出部と、上記固有ベクトル算出部により算出された固有ベクトルを雑音空間とするＭＵＳＩＣ処理を実施して、目標信号の遅延時間を推定するＭＵＳＩＣ処理部とから構成されていることを特徴とする請求項１記載の遅延時間推定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、例えば、目標を検出するレーダ装置に搭載される遅延時間推定装置に関するものである。
【背景技術】
【０００２】
従来の遅延時間推定装置は、例えば、以下の非特許文献１に開示されており、従来の遅延時間推定装置の概略構成は以下の通りである。
（１）電波を生成する送信機Ａ
（２）送信機Ａにより生成された電波を目標に向けて送信する送信アンテナＢ
（３）目標に反射された電波を受信する受信アンテナＣ
（４）受信アンテナＣの受信波に帯域制限や位相検波を施して受信信号を生成する受信機Ｄ
（５）受信機Ｄにより生成された受信信号をアナログ／ディジタル変換して、ディジタル信号を出力するＡ／Ｄ変換器Ｅ
（６）Ａ／Ｄ変換器Ｅから出力されたディジタル信号に高速フーリエ変換処理であるＦＦＴ（Ｆａｓｔ  Ｆｏｕｒｉｅｒ  Ｔｒａｎｓｆｏｒｍ）処理を施して、受信信号の周波数スペクトルを求めるＦＦＴ処理部Ｆ
（７）ＦＦＴ処理部Ｆにより求められた受信信号の周波数スペクトルを送信波の周波数スペクトルで除算して除算信号を生成する除算信号生成部Ｇ
（８）除算信号生成部Ｇにより生成された除算信号に超分解能遅延時間推定処理を施して、目標信号の遅延時間を推定する超分解能推定処理部Ｈ
【０００３】
次に動作について説明する。
送信機Ａが電波を生成すると、送信アンテナＢが送信機Ａにより生成された電波を目標に向けて送信する。
受信機Ｄは、送信アンテナＢから送信されて、目標に反射された電波を受信アンテナＣが受信すると、受信アンテナＣの受信波に帯域制限や位相検波を施して受信信号を生成する。
ここで、送信アンテナＢから電波が送信されて、受信アンテナＣにより電波が受信されるまでの時間を遅延時間τdとする。
【０００４】
Ａ／Ｄ変換器Ｅは、受信機Ｄが受信信号を生成すると、その受信信号をアナログ／ディジタル変換して、ディジタルの受信信号ｓ（１），ｓ（２），・・・，ｓ（Ｎd）をＦＦＴ処理部Ｆに出力する。
ここで、受信信号ｓ（１），ｓ（２），・・・，ｓ（Ｎd）は、目標の信号成分と受信機Ｄの雑音成分からなる信号である。
【０００５】
ＦＦＴ処理部Ｆは、Ａ／Ｄ変換器Ｅから受信信号ｓ（１），ｓ（２），・・・，ｓ（Ｎd）を受けると、その受信信号ｓ（１），ｓ（２），・・・，ｓ（Ｎd）にＦＦＴ処理を施して、受信信号の周波数スペクトルｙ（１），ｙ（２），・・・，ｙ（Ｎd）を算出する。
除算信号生成部Ｇは、ＦＦＴ処理部Ｆが受信信号の周波数スペクトルｙ（１），ｙ（２），・・・，ｙ（Ｎd）を算出すると、下記の式（１）に示すように、その受信信号の周波数スペクトルｙ（１），ｙ（２），・・・，ｙ（Ｎd）を送信波の周波数スペクトルΓ（１），Γ（２），・・・，Γ（Ｎd）で除算して除算信号ｘ（１），ｘ（２），・・・，ｘ（Ｎd）を生成する。
ｘ（ｉ）＝ｙ（ｉ）／Γ（ｉ）                                    （１）
ただし、１≦ｉ≦Ｎd
【０００６】
超分解能推定処理部Ｈは、除算信号生成部Ｇが除算信号ｘ（１），ｘ（２），・・・，ｘ（Ｎd）を生成すると、その除算信号ｘ（１），ｘ（２），・・・，ｘ（Ｎd）に超分解能遅延時間推定処理を施して、目標信号の遅延時間を推定する。
以下、超分解能推定処理部Ｈの処理内容を具体的に説明する。
【０００７】
超分解能推定処理部Ｈの相関行列生成部Ｈ１は、除算信号生成部Ｇが除算信号ｘ（１），ｘ（２），・・・，ｘ（Ｎd）を生成すると、下記の式（２）を演算して、除算信号ｘ（１），ｘ（２），・・・，ｘ（Ｎd）間の相関を表す相関行列Ｒを生成する。
【数１】
2008107242.tif 000003
式（２）において、Ｍdは相関行列の次元数、ｘ*mはベクトルｘmの共役転置を表している。
【０００８】
超分解能推定処理部Ｈの固有ベクトル算出部Ｈ２は、相関行列生成部Ｈが相関行列Ｒを生成すると、その相関行列Ｒの固有値ｖe（１），ｖe（２），・・・，ｖe（Ｍd）（ｖe（１）＞ｖe（２）＞・・・＞ｖe（Ｍd））と、その固有値ｖe（ｉ）（１≦ｉ≦Ｍd）に対応する固有ベクトルｅ（ｉ）を算出する。
超分解能推定処理部Ｈの固有ベクトル算出部Ｈ２は、上記のようにして、固有値ｖe（ｉ）（１≦ｉ≦Ｍd）に対応する固有ベクトルｅ（ｉ）を算出すると、目標数をＫとして、固有ベクトルｅ（Ｋ＋１），ｅ（Ｋ＋２），・・・，ｅ（Ｍd）をＭＵＳＩＣ処理部Ｈ３に出力する。
【０００９】
超分解能推定処理部ＨのＭＵＳＩＣ処理部Ｈ３は、固有ベクトル算出部Ｈ２から固有ベクトルｅ（Ｋ＋１），ｅ（Ｋ＋２），・・・，ｅ（Ｍd）を受けると、その固有ベクトルｅ（Ｋ＋１），ｅ（Ｋ＋２），・・・，ｅ（Ｍd）を雑音空間とするＭＵＳＩＣ（ＭＵｌｔｉｐｌｅ  ＳＩｇｎａｌ  Ｃｌａｓｓｉｆｉｃａｔｉｏｎ）処理を実施する。
即ち、ＭＵＳＩＣ処理部Ｈ３は、Ｔをサンプリング間隔、τkをｋ番目の目標の遅延時間として、下記の式（３）からステアリングベクトルａ（τ）を算出する。
【数２】
2008107242.tif 000004
そして、ＭＵＳＩＣ処理部Ｈ３は、固有ベクトルｅ（Ｋ＋１），ｅ（Ｋ＋２），・・・，ｅ（Ｍd）の全てに直交するｋ種類のステアリングベクトルａ（τ）に対応する遅延時間τkを求める。
【００１０】
なお、従来の遅延時間推定装置では、十分な遅延時間の推定精度を得るためには、送受信帯域１／Ｔと比較して、相関値の算出幅Ｍd／（ＮdＴ）を大きくとる必要がある。
そのため、相関行列の次元数Ｍdが大きくなり、ＭＵＳＩＣ処理の演算量が大きくなることがある。
例えば、相関行列の次元数Ｍdを入力信号の点数Ｎdの半分とすると、Ｍd＝Ｎd／２となり、ＭＵＳＩＣ処理の演算量は、Ｍd3＝Ｎd3／８のオーダーになる。
【００１１】
【非特許文献１】「アレーアンテナによる適応信号処理」  菊間著  １９９８年  科学技術出版
【発明の開示】
【発明が解決しようとする課題】
【００１２】
従来の遅延時間推定装置は以上のように構成されているので、送受信帯域１／Ｔと比較して、相関値の算出幅Ｍd／（ＮdＴ）を大きくとれば、十分な遅延時間の推定精度を得ることができる。しかし、相関値の算出幅Ｍd／（ＮdＴ）を大きくとると、相関行列の次元数Ｍdが大きくなるため、ＭＵＳＩＣ処理の演算量が膨大なものになるなどの課題があった。
【００１３】
この発明は上記のような課題を解決するためになされたもので、遅延時間の推定精度の劣化を招くことなく、相関行列の次元数Ｍdを小さくして、演算量の低減を図ることができる遅延時間推定装置を得ることを目的とする。
【課題を解決するための手段】
【００１４】
この発明に係る遅延時間推定装置は、電波受信手段により生成された受信信号を予め設定された複数の遅延時間推定区間の部分信号に分割し、複数の遅延時間推定区間の部分信号にデシメーション処理を施すデシメーション手段を設け、遅延時間推定手段がデシメーション手段による複数のデシメーション処理信号に超分解能遅延時間推定処理を施して、遅延時間推定区間に存在している目標信号の遅延時間を推定するようにしたものである。
【発明の効果】
【００１５】
この発明によれば、電波受信手段により生成された受信信号を予め設定された複数の遅延時間推定区間の部分信号に分割し、複数の遅延時間推定区間の部分信号にデシメーション処理を施すデシメーション手段を設け、遅延時間推定手段がデシメーション手段による複数のデシメーション処理信号に超分解能遅延時間推定処理を施して、遅延時間推定区間に存在している目標信号の遅延時間を推定するように構成したので、遅延時間の推定精度の劣化を招くことなく、相関行列の次元数Ｍdを小さくして、演算量の低減を図ることができる効果がある。
【発明を実施するための最良の形態】
【００１６】
実施の形態１．
図１はこの発明の実施の形態１による遅延時間推定装置を示す構成図である。
図において、送信機１は目標３に向けて送信する電波を生成する。
送信アンテナ２は送信機１により生成された電波を目標３に向けて送信する。なお、送信機１及び送信アンテナ２から電波送信手段が構成されている。
受信アンテナ４は送信アンテナ２から送信されて目標３に反射された電波を受信する。
受信機５は受信アンテナ４により受信された電波に帯域制限や位相検波を施して受信信号を生成する。
Ａ／Ｄ変換器６は受信機５により生成された受信信号をアナログ／ディジタル変換して、ディジタルの受信信号ｓ（１），ｓ（２），・・・，ｓ（Ｎd）をデシメーション処理部７に出力する。
なお、受信アンテナ４、受信機５及びＡ／Ｄ変換器６から電波受信手段が構成されている。
【００１７】
デシメーション処理部７はＡ／Ｄ変換器６から出力されたディジタルの受信信号ｓ（１），ｓ（２），・・・，ｓ（Ｎd）を予め設定された複数の遅延時間推定区間＃ｌ（１≦ｌ≦Ｌ）の部分信号に分割し、複数の遅延時間推定区間＃ｌの部分信号にデシメーション処理を施して、複数のデシメーション処理信号ｘl（１），ｘl（２），・・・，ｘl（Ｍ）（１≦ｌ≦Ｌ）を出力する。なお、デシメーション処理部７はデシメーション手段を構成している。
【００１８】
遅延時間推定部８−１〜８−Ｌはデシメーション処理部７による複数のデシメーション処理信号ｘl（１），ｘl（２），・・・，ｘl（Ｍ）（１≦ｌ≦Ｌ）に超分解能遅延時間推定処理を施して、遅延時間推定区間＃ｌに存在している目標信号の遅延時間τを推定する。なお、遅延時間推定部８−１〜８−Ｌは遅延時間推定手段を構成している。
【００１９】
図２はこの発明の実施の形態１による遅延時間推定装置のデシメーション処理部７を示す構成図である。
図において、受信信号分割部１１はＡ／Ｄ変換器６から出力されたディジタルの受信信号ｓ（１），ｓ（２），・・・，ｓ（Ｎd）を予め設定された複数の遅延時間推定区間＃ｌ（１≦ｌ≦Ｌ）の部分信号に分割し、複数の遅延時間推定区間＃ｌの部分信号を受信信号スペクトル生成部１２−１〜１２−Ｌに出力する。
即ち、受信信号分割部１１は所要の最大遅延時間ＮdＴをＬ等分し、そのサンプリング分に受信信号長を付加することにより、Ｎサンプリングの部分信号を生成する。例えば、受信信号スペクトル生成部１２−ｌには、Ｍ＝Ｎd／Ｌとして、部分信号ｓ（１＋Ｍ（ｌ−１）），ｓ（２＋Ｍ（ｌ−１）），・・・，ｓ（Ｎ＋Ｍ（ｌ−１））が出力される。
【００２０】
受信信号スペクトル生成部１２−１〜１２−Ｌは受信信号分割部１１により分割された部分信号ｓ（１＋Ｍ（ｌ−１）），ｓ（２＋Ｍ（ｌ−１）），・・・，ｓ（Ｎ＋Ｍ（ｌ−１））にＦＦＴ処理を施して受信信号スペクトルｙl（１），ｙl（２），・・・，ｙl（Ｎ）を生成する。ただし、１≦ｌ≦Ｌである。
除算信号生成部１３−１〜１３−Ｌは受信信号スペクトル生成部１２−１〜１２−Ｌにより生成された受信信号スペクトルｙl（１），ｙl（２），・・・，ｙl（Ｎ）を送信波スペクトルΓ（１），Γ（２），・・・，Γ（Ｎ）で除算して除算信号ｘ'l（１），ｘ'l（２），・・・，ｘ'l（Ｎ）を生成する。ただし、１≦ｌ≦Ｌである。
デシメーションフィルタ１４−１〜１４−Ｌは除算信号生成部１３−１〜１３−Ｌにより生成された除算信号ｘ'l（１），ｘ'l（２），・・・，ｘ'l（Ｎ）にデシメーション処理を施して、デシメーション処理信号ｘl（１），ｘl（２），・・・，ｘl（Ｍ）を出力する。ただし、１≦ｌ≦Ｌである。
【００２１】
図３はこの発明の実施の形態１による遅延時間推定装置の遅延時間推定部８−ｌ（１≦ｌ≦Ｌ）を示す構成図である。
図において、相関行列生成部２１はデシメーション処理部７によるデシメーション処理信号ｘl（１），ｘl（２），・・・，ｘl（Ｍ）（１≦ｌ≦Ｌ）間の相関を表す相関行列Ｒlを生成する。
固有ベクトル算出部２２は相関行列生成部２１により生成された相関行列Ｒlの固有値ｖe,l（１），ｖe,l（２），・・・，ｖe,l（Ｍd）（ｖe,l（１）＞ｖe,l（２）＞・・・＞ｖe,l（Ｍd））を算出して、その固有値ｖe,l（ｉ）（１≦ｉ≦Ｍd）に対応する固有ベクトルｅl（ｉ）を算出し、目標数をＫlとして、固有ベクトルｅl（Ｋl＋１），ｅl（Ｋl＋２），・・・，ｅl（Ｍd）をＭＵＳＩＣ処理部２３に出力する。
ＭＵＳＩＣ処理部２３は固有ベクトル算出部２２から出力された固有ベクトルｅl（Ｋl＋１），ｅl（Ｋl＋２），・・・，ｅl（Ｍd）を雑音空間とするＭＵＳＩＣ処理を実施し、遅延時間推定区間＃ｌに含まれている目標信号の遅延時間τを推定する。
【００２２】
次に動作について説明する。
送信機１が電波を生成すると、送信アンテナ２が送信機１により生成された電波を目標３に向けて送信する。
受信機５は、送信アンテナ２から送信されて、目標３に反射された電波を受信アンテナ４が受信すると、受信アンテナ４の受信波に帯域制限や位相検波を施して受信信号を生成する。
ここで、送信アンテナ２から電波が送信されて、受信アンテナ４により電波が受信されるまでの時間を遅延時間τdとする。
【００２３】
Ａ／Ｄ変換器６は、受信機５が受信信号を生成すると、その受信信号をアナログ／ディジタル変換して、ディジタルの受信信号ｓ（１），ｓ（２），・・・，ｓ（Ｎd）をデシメーション処理部７に出力する。
ここで、受信信号ｓ（１），ｓ（２），・・・，ｓ（Ｎd）は、目標の信号成分と受信機５の雑音成分からなる信号である。
【００２４】
デシメーション処理部７は、Ａ／Ｄ変換器６からディジタルの受信信号ｓ（１），ｓ（２），・・・，ｓ（Ｎd）を受けると、ディジタルの受信信号ｓ（１），ｓ（２），・・・，ｓ（Ｎd）を予め設定された複数の遅延時間推定区間＃ｌ（１≦ｌ≦Ｌ）の部分信号に分割し、複数の遅延時間推定区間＃ｌの部分信号にデシメーション処理を施して、複数のデシメーション処理信号ｘl（１），ｘl（２），・・・，ｘl（Ｍ）（１≦ｌ≦Ｌ）を出力する。
デシメーション処理部７の具体的な処理内容は次の通りである。
【００２５】
デシメーション処理部７の受信信号分割部１１は、Ａ／Ｄ変換器６からディジタルの受信信号ｓ（１），ｓ（２），・・・，ｓ（Ｎd）を受けると、図４に示すように、ディジタルの受信信号ｓ（１），ｓ（２），・・・，ｓ（Ｎd）を予め設定された複数の遅延時間推定区間＃ｌ（１≦ｌ≦Ｌ）の部分信号に分割する。
即ち、受信信号分割部１１は、所要の最大遅延時間ＮdＴをＬ等分（受信信号スペクトル生成部１２−１〜１２−Ｌの個数分に分割）し、そのサンプリング分に受信信号の信号長を付加することにより、Ｎサンプリングの部分信号を生成する。
これにより、例えば、受信信号スペクトル生成部１２−ｌには、Ｍ＝Ｎd／Ｌとして、部分信号ｓ（１＋Ｍ（ｌ−１）），ｓ（２＋Ｍ（ｌ−１）），・・・，ｓ（Ｎ＋Ｍ（ｌ−１））が出力される。
【００２６】
デシメーション処理部７の受信信号スペクトル生成部１２−１〜１２−Ｌは、受信信号分割部１１から部分信号ｓ（１＋Ｍ（ｌ−１）），ｓ（２＋Ｍ（ｌ−１）），・・・，ｓ（Ｎ＋Ｍ（ｌ−１））を受けると、その部分信号ｓ（１＋Ｍ（ｌ−１）），ｓ（２＋Ｍ（ｌ−１）），・・・，ｓ（Ｎ＋Ｍ（ｌ−１））にＦＦＴ処理を施して受信信号スペクトルｙl（１），ｙl（２），・・・，ｙl（Ｎ）を生成する。ただし、１≦ｌ≦Ｌである。
【００２７】
デシメーション処理部７の除算信号生成部１３−１〜１３−Ｌは、受信信号スペクトル生成部１２−１〜１２−Ｌが受信信号スペクトルｙl（１），ｙl（２），・・・，ｙl（Ｎ）を生成すると、下記の式（４）に示すように、その受信信号スペクトルｙl（１），ｙl（２），・・・，ｙl（Ｎ）を送信波スペクトルΓ（１），Γ（２），・・・，Γ（Ｎ）で除算して除算信号ｘ'l（１），ｘ'l（２），・・・，ｘ'l（Ｎ）を生成する。
ｘ'l（ｉ）＝ｙl（ｉ）／Γ（ｉ）                                （４）
ただし、１≦ｉ≦Ｎ、１≦ｌ≦Ｌである。
【００２８】
デシメーション処理部７のデシメーションフィルタ１４−１〜１４−Ｌは、除算信号生成部１３−１〜１３−Ｌが除算信号ｘ'l（１），ｘ'l（２），・・・，ｘ'l（Ｎ）を生成すると、その除算信号ｘ'l（１），ｘ'l（２），・・・，ｘ'l（Ｎ）にデシメーション処理を施して、デシメーション処理信号ｘl（１），ｘl（２），・・・，ｘl（Ｍ）を出力する。ただし、１≦ｌ≦Ｌである。
【００２９】
ここで、図５はデシメーションフィルタ１４−ｌ（１≦ｌ≦Ｌ）の内部を示す構成図である。
デシメーションフィルタ１４−ｌでは、Ｍ’＝Ｌ／Ｍとして、重みｗ（ｍ）（１≦ｍ≦Ｍ’）が下記の式（５）によって設定されており、下記の式（６）を演算することによりデシメーション処理を実施して、デシメーション処理信号ｘl（ｉ）（１≦ｉ≦Ｍ、１≦ｌ≦Ｌ）を生成する。
【数３】
2008107242.tif 000005
【００３０】
遅延時間推定部８−１〜８−Ｌは、デシメーション処理部７がデシメーション処理信号ｘl（ｉ）（１≦ｉ≦Ｍ、１≦ｌ≦Ｌ）を生成すると、そのデシメーション処理信号ｘl（ｉ）に超分解能遅延時間推定処理を施し、遅延時間推定区間＃ｌに存在している目標数をＫlとして、目標信号の遅延時間τdを推定する。
遅延時間推定部８−１〜８−Ｌの具体的な処理内容は次の通りである。
【００３１】
遅延時間推定部８−ｌ（１≦ｌ≦Ｌ）の相関行列生成部２１は、デシメーション処理部７がデシメーション処理信号ｘl（１），ｘl（２），・・・，ｘl（Ｍ）を生成すると、下記の式（７）を演算して、デシメーション処理信号ｘl（１），ｘl（２），・・・，ｘl（Ｍ）間の相関を表す相関行列Ｒlを生成する。
【数４】
2008107242.tif 000006
式（７）において、Ｍdは相関行列の次元数、ｘ*mはベクトルｘmの共役転置を表している。
【００３２】
遅延時間推定部８−ｌの固有ベクトル算出部２２は、相関行列生成部２１が相関行列Ｒlを生成すると、その相関行列Ｒlの固有値ｖe,l（１），ｖe,l（２），・・・，ｖe,l（Ｍd）（ｖe,l（１）＞ｖe,l（２）＞・・・＞ｖe,l（Ｍd））と、その固有値ｖe,l（ｉ）（１≦ｉ≦Ｍd）に対応する固有ベクトルｅl（ｉ）を算出する。
固有ベクトル算出部２２は、上記のようにして、固有値ｖe,l（ｉ）（１≦ｉ≦Ｍd）に対応する固有ベクトルｅl（ｉ）を算出すると、目標数をＫlとして、固有ベクトルｅl（Ｋ＋１），ｅl（Ｋ＋２），・・・，ｅl（Ｍd）をＭＵＳＩＣ処理部２３に出力する。
【００３３】
遅延時間推定部８−ｌのＭＵＳＩＣ処理部２３は、固有ベクトル算出部２２から固有ベクトルｅl（Ｋl＋１），ｅl（Ｋl＋２），・・・，ｅl（Ｍd）を受けると、その固有ベクトルｅl（Ｋl＋１），ｅl（Ｋl＋２），・・・，ｅl（Ｍd）を雑音空間とするＭＵＳＩＣ処理を実施する。
即ち、ＭＵＳＩＣ処理部２３は、Ｔをサンプリング間隔として下記の式（８）からステアリングベクトルａ（τ）を算出する。
【数５】
2008107242.tif 000007
そして、ＭＵＳＩＣ処理部２３は、固有ベクトルｅl（Ｋl＋１），ｅl（Ｋl＋２），・・・，ｅl（Ｍd）の全てに直交するKl種類のステアリングベクトルａ（τ1）, ａ（τ2）,…, ａ（τKl）を求める。これより、ｋ（1≦k≦Kl）番目の目標の遅延時間τkが求まる。
【００３４】
この実施の形態１では、上記の従来例と比較して、相関行列の次元数Ｍdを小さくすることができるため、処理の演算量を低減することができる。
例えば、相関行列の次元数Ｍdを入力信号の点数Ｍの半分とすると、Ｍd＝Ｍ／２＝Ｎd／２Ｌとなり、ＭＵＳＩＣ処理の演算量は、ＬＭd3＝Ｎd3／８Ｌ2のオーダーになる。
従来例におけるＭＵＳＩＣ処理の演算量は、Ｎd3／８であるため、ＭＵＳＩＣ処理の演算量が１／Ｌ2に低減されている。
【００３５】
以上で明らかなように、この実施の形態１によれば、Ａ／Ｄ変換器６から出力されたディジタルの受信信号ｓ（１），ｓ（２），・・・，ｓ（Ｎd）を予め設定された複数の遅延時間推定区間＃ｌ（１≦ｌ≦Ｌ）の部分信号に分割し、複数の遅延時間推定区間＃ｌの部分信号にデシメーション処理を施して、複数のデシメーション処理信号ｘl（１），ｘl（２），・・・，ｘl（Ｍ）（１≦ｌ≦Ｌ）を出力するデシメーション処理部７を設け、遅延時間推定部８−１〜８−Ｌがデシメーション処理部７による複数のデシメーション処理信号ｘl（１），ｘl（２），・・・，ｘl（Ｍ）（１≦ｌ≦Ｌ）に超分解能遅延時間推定処理を施して、遅延時間推定区間＃ｌに存在している目標信号の遅延時間を推定するように構成したので、遅延時間の推定精度の劣化を招くことなく、相関行列の次元数Ｍdを小さくして、演算量の低減を図ることができる効果を奏する。
【００３６】
実施の形態２．
図６はこの発明の実施の形態２による遅延時間推定装置のデシメーション処理部７を示す構成図であり、図７はデシメーション処理部７における推定区間切り出し型デシメーション部１５−１〜１５−Ｌの内部を示す構成図である。
この発明の実施の形態２による遅延時間推定装置の全体構成は、図１の遅延時間推定装置と同一である。
図６において、推定区間切り出し型デシメーション部１５−１〜１５−Ｌ以外は、図２のデシメーション処理部７と同じである。
【００３７】
図７において、推定区間切り出し型デシメーション部１５−ｌ（１≦ｌ≦Ｌ）のＩＦＦＴ処理部（逆フーリエ変換処理部）３１は除算信号生成部１３−ｌにより生成された除算信号ｘ'l（１），ｘ'l（２），・・・，ｘ'l（Ｎ）にＩＦＦＴ（逆フーリエ変換）処理を施して復調信号ｕl（１），ｕl（２），・・・，ｕl（Ｎ）を生成する。
信号切出処理部である推定区間切り出し処理部３２はＩＦＦＴ処理部３１により生成された復調信号ｕl（１），ｕl（２），・・・，ｕl（Ｎ）から遅延時間推定区間＃ｌの信号ｕl（１），ｕl（２），・・・，ｕl（Ｍ）の切り出しを行う。
フーリエ変換処理部であるＦＦＴ処理部３３は推定区間切り出し処理部３２により切り出された遅延時間推定区間＃ｌの信号ｕl（１），ｕl（２），・・・，ｕl（Ｍ）にＦＦＴ処理を施して、デシメーション処理信号ｘl（１），ｘl（２），・・・，ｘl（Ｍ）を出力する。
【００３８】
上記実施の形態１では、デシメーション処理部７がデシメーションフィルタ１４−１〜１４−Ｌを実装しているものについて示したが、図６に示すように、デシメーションフィルタ１４−１〜１４−Ｌの代わりに、推定区間切り出し型デシメーション部１５−１〜１５−Ｌを実装するようにしても、上記実施の形態１と同様に、相関行列の次元数Ｍdを小さくすることができるため、処理の演算量を低減することができる。
【００３９】
以下、推定区間切り出し型デシメーション部１５−１〜１５−Ｌの処理内容を具体的に説明する。
推定区間切り出し型デシメーション部１５−ｌ（１≦ｌ≦Ｌ）のＩＦＦＴ処理部３１は、上記実施の形態１と同様にして、除算信号生成部１３−ｌが除算信号ｘ'l（１），ｘ'l（２），・・・，ｘ'l（Ｎ）を生成すると、その除算信号ｘ'l（１），ｘ'l（２），・・・，ｘ'l（Ｎ）にＩＦＦＴ処理を施して、復調信号ｕl（１），ｕl（２），・・・，ｕl（Ｎ）を生成する。
【００４０】
推定区間切り出し処理部３２は、ＩＦＦＴ処理部３１が復調信号ｕl（１），ｕl（２），・・・，ｕl（Ｎ）を生成すると、その復調信号ｕl（１），ｕl（２），・・・，ｕl（Ｎ）から遅延時間推定区間＃ｌの信号ｕl（１），ｕl（２），・・・，ｕl（Ｍ）の切り出しを行う（図４を参照）。
ＦＦＴ処理部３３は、推定区間切り出し処理部３２が遅延時間推定区間＃ｌの信号ｕl（１），ｕl（２），・・・，ｕl（Ｍ）の切り出しを行うと、遅延時間推定区間＃ｌの信号ｕl（１），ｕl（２），・・・，ｕl（Ｍ）にＦＦＴ処理を施して、デシメーション処理信号ｘl（１），ｘl（２），・・・，ｘl（Ｍ）を出力する。
以下、上記実施の形態１と同様であるため説明を省略する。
【００４１】
実施の形態３．
図８はこの発明の実施の形態３による遅延時間推定装置のデシメーション処理部７を示す構成図であり、この発明の実施の形態３による遅延時間推定装置の全体構成は、図１の遅延時間推定装置と同一である。
図８において、図２と同一符号は同一または相当部分を示すので説明を省略する。
零割防止型除算信号生成部１６−１〜１６−Ｌは、図２の除算信号生成部１３−１〜１３−Ｌと同様に、受信信号スペクトル生成部１２−１〜１２−Ｌにより生成された受信信号スペクトルｙl（１），ｙl（２），・・・，ｙl（Ｎ）を送信波スペクトルΓ（１），Γ（２），・・・，Γ（Ｎ）で除算して除算信号ｘ'l（１），ｘ'l（２），・・・，ｘ'l（Ｎ）を生成するが、送信波スペクトルΓ（１），Γ（２），・・・，Γ（Ｎ）の振幅値が予め設定された最小値Γminより小さい場合、その送信波スペクトルΓ（１），Γ（２），・・・，Γ（Ｎ）を最小値Γminに修正し、修正後の送信波スペクトルΓ'（１），Γ'（２），・・・，Γ'（Ｎ）で受信信号スペクトルｙl（１），ｙl（２），・・・，ｙl（Ｎ）を除算する。ただし、１≦ｌ≦Ｌである。
【００４２】
上記実施の形態１，２では、除算信号生成部１３−１〜１３−Ｌが受信信号スペクトル生成部１２−１〜１２−Ｌにより生成された受信信号スペクトルｙl（１），ｙl（２），・・・，ｙl（Ｎ）を送信波スペクトルΓ（１），Γ（２），・・・，Γ（Ｎ）で除算して除算信号ｘ'l（１），ｘ'l（２），・・・，ｘ'l（Ｎ）を生成するが、送信波スペクトルΓ（１），Γ（２），・・・，Γ（Ｎ）の振幅値が予め想定している値より小さいと、除算信号ｘ'l（１），ｘ'l（２），・・・，ｘ'l（Ｎ）が突発的に大きな誤差を持つことがある。
【００４３】
そこで、この実施の形態３では、送信波スペクトルΓ（１），Γ（２），・・・，Γ（Ｎ）の振幅値が小さいことに起因する突発的な大きな誤差の発生を防止するようにしている。
即ち、零割防止型除算信号生成部１６−１〜１６−Ｌは、予め、送信波スペクトルの最小値Γminを設定し、受信信号スペクトルｙl（１），ｙl（２），・・・，ｙl（Ｎ）を送信波スペクトルΓ（１），Γ（２），・・・，Γ（Ｎ）で除算するに先立って、受信信号スペクトルｙl（１），ｙl（２），・・・，ｙl（Ｎ）の振幅値を最小値Γminと比較する。
【００４４】
そして、零割防止型除算信号生成部１６−１〜１６−Ｌは、受信信号スペクトルｙl（１），ｙl（２），・・・，ｙl（Ｎ）の振幅値が最小値Γminより小さい場合、下記の式（９）に示すように、その送信波スペクトルΓ（１），Γ（２），・・・，Γ（Ｎ）を最小値Γminに修正する。
【数６】
2008107242.tif 000008
【００４５】
零割防止型除算信号生成部１６−１〜１６−Ｌは、受信信号スペクトルｙl（１），ｙl（２），・・・，ｙl（Ｎ）の振幅値が最小値Γminより小さくなければ、受信信号スペクトルｙl（１），ｙl（２），・・・，ｙl（Ｎ）を送信波スペクトルΓ（１），Γ（２），・・・，Γ（Ｎ）で除算して除算信号ｘ'l（１），ｘ'l（２），・・・，ｘ'l（Ｎ）を生成するが、受信信号スペクトルｙl（１），ｙl（２），・・・，ｙl（Ｎ）の振幅値が最小値Γminより小さければ、修正後の送信波スペクトルΓ'（１），Γ'（２），・・・，Γ'（Ｎ）で受信信号スペクトルｙl（１），ｙl（２），・・・，ｙl（Ｎ）を除算して除算信号ｘ'l（１），ｘ'l（２），・・・，ｘ'l（Ｎ）を生成する。
【００４６】
これにより、この実施の形態３によれば、送信波スペクトルΓ（１），Γ（２），・・・，Γ（Ｎ）の振幅値が予め想定している値より小さい場合でも、除算信号ｘ'l（１），ｘ'l（２），・・・，ｘ'l（Ｎ）が突発的に大きな誤差を持つことを防止することができる効果を奏する。
【００４７】
実施の形態４．
図９はこの発明の実施の形態４による遅延時間推定装置の遅延時間推定部８−ｌ（１≦ｌ≦Ｌ）を示す構成図であり、この発明の実施の形態４による遅延時間推定装置の全体構成は、図１の遅延時間推定装置と同一である。
図９において、図３と同一符号は同一または相当部分を示すので説明を省略する。
目標判別部である目標信号判定型固有ベクトル算出部２４は相関行列生成部２１により生成された相関行列Ｒlの固有値ｖe,l（１），ｖe,l（２），・・・，ｖe,l（Ｍd）（ｖe,l（１）＞ｖe,l（２）＞・・・＞ｖe,l（Ｍd））を算出して、その固有値ｖe,l（ｉ）（１≦ｉ≦Ｍd）に対応する固有ベクトルｅl（ｉ）を算出し、所定の閾値Ｔｈより大きい固有値があれば、目標信号が存在するものと認定して固有ベクトルｅl（ｉ）をＭＵＳＩＣ処理部２３に出力する一方、所定の閾値Ｔｈより大きい固有値がなければ、目標信号が存在していないと認定する。
【００４８】
遅延時間推定装置の遅延時間推定部８−ｌ（１≦ｌ≦Ｌ）以外は、上記実施の形態１〜３と同様であるため、以下、遅延時間推定部８−ｌ（１≦ｌ≦Ｌ）の処理内容を具体的に説明する。
遅延時間推定部８−ｌの相関行列生成部２１は、デシメーション処理部７がデシメーション処理信号ｘl（１），ｘl（２），・・・，ｘl（Ｍ）を生成すると、上記実施の形態１と同様に、デシメーション処理信号ｘl（１），ｘl（２），・・・，ｘl（Ｍ）間の相関を表す相関行列Ｒlを生成する。
【００４９】
遅延時間推定部８−ｌの目標信号判定型固有ベクトル算出部２４は、相関行列生成部２１が相関行列Ｒlを生成すると、その相関行列Ｒlの固有値ｖe,l（１），ｖe,l（２），・・・，ｖe,l（Ｍd）（ｖe,l（１）＞ｖe,l（２）＞・・・＞ｖe,l（Ｍd））と、その固有値ｖe,l（ｉ）（１≦ｉ≦Ｍd）に対応する固有ベクトルｅl（ｉ）を算出する。
目標信号判定型固有ベクトル算出部２４は、上記のようにして、固有値ｖe,l（ｉ）（１≦ｉ≦Ｍd）に対応する固有ベクトルｅl（ｉ）を算出すると、予め設定された閾値Ｔｈと固有値ｖe,l（１），ｖe,l（２），・・・，ｖe,l（Ｍd）を比較し、閾値Ｔｈ以上の固有値を求める。
【００５０】
目標信号判定型固有ベクトル算出部２４は、全ての固有値ｖe,l（１），ｖe,l（２），・・・，ｖe,l（Ｍd）が閾値Ｔｈより小さい場合、目標信号が存在していないと認定して、一連の処理を終了する。この場合、固有ベクトルｅl（ｉ）がＭＵＳＩＣ処理部２３には出力されず、目標信号の遅延時間の推定は行われない。
目標信号判定型固有ベクトル算出部２４は、固有値ｖe,l（１），ｖe,l（２），・・・，ｖe,l（Ｍd）の中で、閾値Ｔｈより大きい固有値がある場合、目標信号が存在するものと認定する。
例えば、ｖe,l（１）＞ｖe,l（２）＞・・・＞ｖe,l（Ｋ）＞Ｔｈであれば、目標数がＫであると認定して、固有ベクトルｅl（Ｋ＋１），ｅl（Ｋ＋２），・・・，ｅl（Ｍd）をＭＵＳＩＣ処理部２３に出力する。
【００５１】
遅延時間推定部８−ｌのＭＵＳＩＣ処理部２３は、目標信号判定型固有ベクトル算出部２４から固有ベクトルｅl（Ｋ＋１），ｅl（Ｋ＋２），・・・，ｅl（Ｍd）を受けると、上記実施の形態１と同様に、その固有ベクトルｅl（Ｋ＋１），ｅl（Ｋ＋２），・・・，ｅ（Ｍd）を雑音空間とするＭＵＳＩＣ処理を実施して、目標信号の遅延時間を推定する。
【００５２】
以上で明らかなように、この実施の形態４によれば、相関行列生成部２１により生成された相関行列Ｒlの固有値ｖe,l（１），ｖe,l（２），・・・，ｖe,l（Ｍd）を算出して、その固有値ｖe,l（ｉ）（１≦ｉ≦Ｍd）に対応する固有ベクトルｅl（ｉ）を算出し、所定の閾値Ｔｈより大きい固有値があれば、目標信号が存在するものと認定して固有ベクトルｅl（ｉ）をＭＵＳＩＣ処理部２３に出力する一方、所定の閾値Ｔｈより大きい固有値がなければ、目標信号が存在していないと認定するように構成したので、目標信号数が未知の場合でも、超分解能遅延時間推定処理を実施することができるとともに、目標信号が存在しないときは無駄な推定処理を省略することができる効果を奏する。
【００５３】
実施の形態５．
図１０はこの発明の実施の形態５による遅延時間推定装置を示す構成図であり、図において、図１と同一符号は同一または相当部分を示すので説明を省略する。
デシメーション処理部９はＡ／Ｄ変換器６から出力されたディジタルの受信信号ｓ（１），ｓ（２），・・・，ｓ（Ｎd）を予め設定された複数の遅延時間推定区間＃ｌ（１≦ｌ≦Ｌ）の部分信号に分割し、複数の遅延時間推定区間＃ｌの部分信号に目標信号が含まれているか否かを判定し、目標信号が含まれている部分信号にデシメーション処理を施して、デシメーション処理信号ｘh（１），ｘh（２），・・・，ｘh（Ｍ）を出力する。なお、デシメーション処理部９はデシメーション手段を構成している。
遅延時間推定部１０はデシメーション処理部９によるデシメーション処理信号ｘh（１），ｘh（２），・・・，ｘh（Ｍ）に超分解能遅延時間推定処理を施して、目標信号の遅延時間τを推定する。なお、遅延時間推定部１０は遅延時間推定手段を構成している。
【００５４】
図１１はこの発明の実施の形態５による遅延時間推定装置のデシメーション処理部９を示す構成図であり、図において、図２と同一符号は同一または相当部分を示すので説明を省略する。
目標信号検出部１７は受信信号スペクトル生成部１２−１〜１２−Ｌにより生成された受信信号スペクトルｙl（１），ｙl（２），・・・，ｙl（Ｎ）（１≦ｌ≦Ｌ）の平均電力Ｐlを算出し、その平均電力Ｐlが基準電力Ｔｈrより大きければ、その受信信号スペクトルｙl（１），ｙl（２），・・・，ｙl（Ｎ）には目標信号が含まれていると判定し、目標信号が含まれている受信信号スペクトルｙh（１），ｙh（２），・・・，ｙh（Ｎ）を除算信号生成部１８に出力する。
【００５５】
除算信号生成部１８は目標信号検出部１７から出力された受信信号スペクトルｙh（１），ｙh（２），・・・，ｙh（Ｎ）を送信波スペクトルΓ（１），Γ（２），・・・，Γ（Ｎ）で除算して除算信号ｘ'h（１），ｘ'h（２），・・・，ｘ'h（Ｎ）を生成する。
デシメーションフィルタ１９は除算信号生成部１８により生成された除算信号ｘ'h（１），ｘ'h（２），・・・，ｘ'h（Ｎ）にデシメーション処理を施して、デシメーション処理信号ｘh（１），ｘh（２），・・・，ｘh（Ｍ）を出力する。
【００５６】
次に動作について説明する。
デシメーション処理部９は、上記実施の形態１と同様にして、Ａ／Ｄ変換器６がディジタルの受信信号ｓ（１），ｓ（２），・・・，ｓ（Ｎd）を出力すると、その受信信号ｓ（１），ｓ（２），・・・，ｓ（Ｎd）を予め設定された複数の遅延時間推定区間＃ｌ（１≦ｌ≦Ｌ）の部分信号に分割し、複数の遅延時間推定区間＃ｌの部分信号に目標信号が含まれているか否かを判定し、目標信号が含まれている部分信号にデシメーション処理を施して、デシメーション処理信号ｘh（１），ｘh（２），・・・，ｘh（Ｍ）を遅延時間推定部１０に出力する。
【００５７】
即ち、デシメーション処理部９の受信信号分割部１１は、Ａ／Ｄ変換器６からディジタルの受信信号ｓ（１），ｓ（２），・・・，ｓ（Ｎd）を受けると、上記実施の形態１と同様に、ディジタルの受信信号ｓ（１），ｓ（２），・・・，ｓ（Ｎd）を予め設定された遅延時間推定区間＃ｌの部分信号ｓ（１＋Ｍ（ｌ−１）），ｓ（２＋Ｍ（ｌ−１）），・・・，ｓ（Ｎ＋Ｍ（ｌ−１））に分割し、その部分信号ｓ（１＋Ｍ（ｌ−１）），ｓ（２＋Ｍ（ｌ−１）），・・・，ｓ（Ｎ＋Ｍ（ｌ−１））を受信信号スペクトル生成部１２−１〜１２−Ｌに出力する。ただし、１≦ｌ≦Ｌである。
【００５８】
デシメーション処理部９の受信信号スペクトル生成部１２−１〜１２−Ｌは、受信信号分割部１１から部分信号ｓ（１＋Ｍ（ｌ−１）），ｓ（２＋Ｍ（ｌ−１）），・・・，ｓ（Ｎ＋Ｍ（ｌ−１））を受けると、上記実施の形態１と同様に、その部分信号ｓ（１＋Ｍ（ｌ−１）），ｓ（２＋Ｍ（ｌ−１）），・・・，ｓ（Ｎ＋Ｍ（ｌ−１））にＦＦＴ処理を施して受信信号スペクトルｙl（１），ｙl（２），・・・，ｙl（Ｎ）を生成する。ただし、１≦ｌ≦Ｌである。
【００５９】
デシメーション処理部９の目標信号検出部１７は、受信信号スペクトル生成部１２−１〜１２−Ｌが受信信号スペクトルｙl（１），ｙl（２），・・・，ｙl（Ｎ）を生成すると、下記の式（１０）に示すように、その受信信号スペクトルｙl（１），ｙl（２），・・・，ｙl（Ｎ）の平均電力Ｐlを算出する。
【数７】
2008107242.tif 000009
ただし、１≦ｌ≦Ｌである。
【００６０】
目標信号検出部１７は、受信信号スペクトルｙl（１），ｙl（２），・・・，ｙl（Ｎ）の平均電力Ｐlと予め設定された基準電力Ｔｈrを比較する。
目標信号検出部１７は、受信信号スペクトルｙl（１），ｙl（２），・・・，ｙl（Ｎ）の平均電力Ｐlが基準電力Ｔｈrより小さければ、その受信信号スペクトルｙl（１），ｙl（２），・・・，ｙl（Ｎ）には目標信号が含まれていないと判定し、その受信信号スペクトルｙl（１），ｙl（２），・・・，ｙl（Ｎ）を除算信号生成部１８には出力しない。
一方、受信信号スペクトルｙl（１），ｙl（２），・・・，ｙl（Ｎ）の平均電力Ｐlが基準電力Ｔｈrより大きければ、その受信信号スペクトルｙl（１），ｙl（２），・・・，ｙl（Ｎ）には目標信号が含まれていると判定し、その受信信号スペクトルｙl（１），ｙl（２），・・・，ｙl（Ｎ）を目標信号が含まれている受信信号スペクトルｙh（１），ｙh（２），・・・，ｙh（Ｎ）として除算信号生成部１８に出力する。
【００６１】
除算信号生成部１８は、目標信号検出部１７から目標信号が含まれている受信信号スペクトルｙh（１），ｙh（２），・・・，ｙh（Ｎ）を受けると、その受信信号スペクトルｙh（１），ｙh（２），・・・，ｙh（Ｎ）を送信波スペクトルΓ（１），Γ（２），・・・，Γ（Ｎ）で除算して除算信号ｘ'h（１），ｘ'h（２），・・・，ｘ'h（Ｎ）を生成する。
デシメーションフィルタ１９は、除算信号生成部１８が除算信号ｘ'h（１），ｘ'h（２），・・・，ｘ'h（Ｎ）を生成すると、その除算信号ｘ'h（１），ｘ'h（２），・・・，ｘ'h（Ｎ）にデシメーション処理を施して、デシメーション処理信号ｘh（１），ｘh（２），・・・，ｘh（Ｍ）を出力する。
【００６２】
遅延時間推定部１０は、デシメーション処理部９から遅延時間推定区間＃ｈのデシメーション処理信号ｘh（１），ｘh（２），・・・，ｘh（Ｍ）を受けると、上記実施の形態１の遅延時間推定部８−１〜８−Ｌと同様に、そのデシメーション処理信号ｘh（１），ｘh（２），・・・，ｘh（Ｍ）に超分解能遅延時間推定処理を施し、遅延時間推定区間＃ｈに存在している目標信号の遅延時間τを推定する。
【００６３】
以上で明らかなように、この実施の形態５によれば、Ａ／Ｄ変換器６から出力されたディジタルの受信信号ｓ（１），ｓ（２），・・・，ｓ（Ｎd）を予め設定された複数の遅延時間推定区間＃ｌ（１≦ｌ≦Ｌ）の部分信号に分割し、複数の遅延時間推定区間＃ｌの部分信号に目標信号が含まれているか否かを判定し、目標信号が含まれている部分信号にデシメーション処理を施して、デシメーション処理信号ｘh（１），ｘh（２），・・・，ｘh（Ｍ）を遅延時間推定部１０に出力するように構成したので、上記実施の形態１と同様の効果を奏する他に、目標信号が存在している遅延時間推定区間が未知の場合でも、超分解能遅延時間推定処理を実施することができる効果を奏する。
【００６４】
実施の形態６．
図１２はこの発明の実施の形態６による遅延時間推定装置を示す構成図であり、図において、図１と同一符号は同一または相当部分を示すので説明を省略する。
受信アンテナ４−１，４−２は送信アンテナ２から送信されて目標３に反射された電波を受信する。
受信機５−１，５−２は受信アンテナ４−１，４−２により受信された電波に帯域制限や位相検波を施して受信信号を生成する。
【００６５】
Ａ／Ｄ変換器６−１，６−２は受信機５−１，５−２により生成された受信信号をアナログ／ディジタル変換して、ディジタルの受信信号ｓ（１），ｓ（２），・・・，ｓ（Ｎd），ｓ'（１），ｓ'（２），・・・，ｓ'（Ｎd）を最大比合成処理部４１に出力する。
なお、受信アンテナ４−１，４−２、受信機５−１，５−２及びＡ／Ｄ変換器６−１，６−２から電波受信手段が構成されている。
最大比合成処理部４１はＡ／Ｄ変換器６−１，６−２から出力されたディジタルの受信信号ｓ（１），ｓ（２），・・・，ｓ（Ｎd），ｓ'（１），ｓ'（２），・・・，ｓ'（Ｎd）を最大比合成して、最大比合成後の受信信号ｓ''（１），ｓ''（２），・・・，ｓ''（Ｎd）をデシメーション処理部７に出力する。なお、最大比合成処理部４１は最大比合成処理手段を構成している。
【００６６】
次に動作について説明する。
送信機１が電波を生成すると、送信アンテナ２が送信機１により生成された電波を目標３に向けて送信する。
受信機５−１，５−２は、送信アンテナ２から送信されて、目標３に反射された電波を受信アンテナ４−１，４−２が受信すると、受信アンテナ４−１，４−２の受信波に帯域制限や位相検波を施して受信信号を生成する。
Ａ／Ｄ変換器６−１，６−２は、受信機５−１，５−２が受信信号を生成すると、その受信信号をアナログ／ディジタル変換して、ディジタルの受信信号ｓ（１），ｓ（２），・・・，ｓ（Ｎd），ｓ'（１），ｓ'（２），・・・，ｓ'（Ｎd）を最大比合成処理部４１に出力する。
【００６７】
最大比合成処理部４１は、Ａ／Ｄ変換器６−１，６−２からディジタルの受信信号ｓ（１），ｓ（２），・・・，ｓ（Ｎd），ｓ'（１），ｓ'（２），・・・，ｓ'（Ｎd）を受けると、その受信信号ｓ（１），ｓ（２），・・・，ｓ（Ｎd）と、受信信号ｓ'（１），ｓ'（２），・・・，ｓ'（Ｎd）とを最大比合成して、最大比合成後の受信信号ｓ''（１），ｓ''（２），・・・，ｓ''（Ｎd）をデシメーション処理部７に出力する。
最大比合成処理部４１の具体的な処理内容は以下の通りである。
【００６８】
最大比合成処理部４１は、Ａ／Ｄ変換器６−１，６−２からディジタルの受信信号ｓ（１），ｓ（２），・・・，ｓ（Ｎd），ｓ'（１），ｓ'（２），・・・，ｓ'（Ｎd）を受けると、下記の式（１１）に示すように、受信信号ｓ（１），ｓ（２），・・・，ｓ（Ｎd）と、受信信号ｓ'（１），ｓ'（２），・・・，ｓ'（Ｎd）との相関行列Ｒrを算出する。ただし、式（１１）において、行列［ｓr（１）・・・ｓr（Ｎd）］Hは、行列［ｓr（１）・・・ｓr（Ｎd）］のエルミート共役を表している。
【数８】
2008107242.tif 000010
【００６９】
最大比合成処理部４１は、２つの受信信号の相関行列Ｒrを算出すると、その相関行列Ｒrの最大固有値に対応する固有ベクトルｅrを算出し、下記の式（１２）に示すように、その固有ベクトルｅr用いて、最大比合成処理信号ｓ''（１），ｓ''（２），・・・，ｓ''（Ｎd）を求め、最大比合成処理信号ｓ''（１），ｓ''（２），・・・，ｓ''（Ｎd）をデシメーション処理部７に出力する。
ｓ''（ｉ）＝ｅrHｓr（ｉ）                                    （１２）
デシメーション処理部７及び遅延時間推定部８−１〜８−Ｌの処理内容は、上記実施の形態１と同様であるため説明を省略する。
【００７０】
以上で明らかなように、この実施の形態６によれば、最大比合成処理部４１がＡ／Ｄ変換器６−１，６−２から出力されたディジタルの受信信号ｓ（１），ｓ（２），・・・，ｓ（Ｎd），ｓ'（１），ｓ'（２），・・・，ｓ'（Ｎd）を最大比合成して、最大比合成後の受信信号ｓ''（１），ｓ''（２），・・・，ｓ''（Ｎd）をデシメーション処理部７に出力するように構成したので、上記実施の形態１と同様の効果を奏する他に、マルチパスフェージング等が発生していても受信信号の電力を改善して、遅延時間の推定精度の劣化を防止することができる効果を奏する。
【００７１】
なお、この実施の形態６では、受信ンテナ４−１，４−２、受信機５−１，５−２及びＡ／Ｄ変換器６−１，６−２を２組用意して、２つの受信信号の最大比合成を行うものについて示したが、２つのＡ／Ｄ変換器６−１，６−２から出力される受信信号のうち、電力が大きい方の受信信号を選択して、その受信信号をデシメーション処理部７に出力する選択性ダイバーシティ方式を採用するようにしてもよい。
【００７２】
実施の形態７．
図１３はこの発明の実施の形態７による遅延時間推定装置を示す構成図であり、図において、図１と同一符号は同一または相当部分を示すので説明を省略する。
オーバーサンプリング型Ａ／Ｄ変換器４２は受信機５により生成された受信信号をアナログ／ディジタル変換する際、送信アンテナ２から送信される電波がパルス単位で符号変調が施されている場合、パルスのチップ幅よりも短い時間間隔で電波をサンプリングしてディジタルの受信信号ｓ（１），ｓ（２），・・・，ｓ（Ｎ'd）を生成する。なお、オーバーサンプリング型Ａ／Ｄ変換器４２は電波受信手段を構成している。
【００７３】
次に動作について説明する。
送信機１は、電波をパルス単位で符号変調し、符号変調後の電波を送信アンテナ２に出力する。
送信アンテナ２は、送信アンテナ２から符号変調後の電波を受けると、その電波を目標３に向けて送信する。
受信機５は、送信アンテナ２から送信されて、目標３に反射された電波を受信アンテナ４が受信すると、受信アンテナ４の受信波に帯域制限や位相検波を施して受信信号を生成する。
【００７４】
オーバーサンプリング型Ａ／Ｄ変換器４２は、受信機５が受信信号を生成すると、図１のＡ／Ｄ変換器６と同様に、アナログの受信信号をアナログ／ディジタル変換するが、その際、パルスのチップ幅よりも短い時間間隔で、アナログの受信信号をサンプリングして、ディジタルの受信信号ｓ（１），ｓ（２），・・・，ｓ（Ｎ'd）を生成する。
ここで、図１４はマルチパスフェージング環境下において、Ａ／Ｄ変換器６における通常のサンプリングと、オーバーサンプリング型Ａ／Ｄ変換器４２におけるオーバーサンプリングとの相違を示す説明図である。
【００７５】
マルチパスフェージング環境下では、Ａ／Ｄ変換器６のように通常のサンプリングを実施すると、図１４に示すように、フェージングしている箇所をサンプリングすることになり、十分な受信信号の電力が得られないことがある。
一方、オーバーサンプリング型Ａ／Ｄ変換器４２のようにオーバーサンプリングを実施すると、図１４に示すように、フェージングしている箇所だけでなく、フェージングしていない箇所もサンプリングすることができるため、十分な受信信号の電力が得られる。
デシメーション処理部７及び遅延時間推定部８−１〜８−Ｌの処理内容は、上記実施の形態１と同様であるため説明を省略する。
【００７６】
以上で明らかなように、この実施の形態７によれば、送信アンテナ２から送信される電波がパルス単位で符号変調が施されている場合、パルスのチップ幅よりも短い時間間隔で電波をサンプリングしてディジタルの受信信号ｓ（１），ｓ（２），・・・，ｓ（Ｎ'd）を生成するように構成したので、上記実施の形態１と同様の効果を奏する他に、マルチパスフェージング等が発生していても、受信信号の電力を改善して、遅延時間の推定精度の劣化を防止することができる効果を奏する。
【００７７】
実施の形態８．
図１５はこの発明の実施の形態８による遅延時間推定装置を示す構成図であり、図において、図１と同一符号は同一または相当部分を示すので説明を省略する。
メモリ回路４３はＡ／Ｄ変換器６からディジタルの受信信号ｓnp（１），ｓnp（２），・・・，ｓnp（Ｎd）（１≦ｎp≦Ｎp）を受ける毎に、その受信信号をパルス方向に並べて蓄積するとともに、その受信信号を時間方向に並べて蓄積する。なお、メモリ回路４３は受信信号蓄積手段を構成している。
【００７８】
パルス方向ＦＦＴ部４４−１〜４４−Ｎdはメモリ回路４３にパルス方向に並べて蓄積されている受信信号ｓ1（ｎd），ｓ2（ｎd），・・・，ｓNp（ｎd）にＦＦＴ処理を施し、パルス方向ＦＦＴ処理信号ｚ1（ｎd），ｚ2（ｎd），・・・，ｚNp（ｎd）を出力する。ただし、１≦ｎd≦Ｎd、１≦ｎp≦Ｎpである。
なお、パルス方向ＦＦＴ部４４−１〜４４−Ｎdはフーリエ変換処理手段を構成している。
遅延時間推定区間検出部４５はパルス方向ＦＦＴ部４４−１〜４４−Ｎdから出力されたパルス方向ＦＦＴ処理信号ｚ1（ｎd），ｚ2（ｎd），・・・，ｚNp（ｎd）を解析して、目標信号が存在している遅延時間推定区間＃ｈを検出する。なお、遅延時間推定区間検出部４５は遅延時間推定区間検出手段を構成している。
【００７９】
デシメーション処理部４６−１〜４６−Ｎpはメモリ回路４３に時間方向に並べて蓄積されている受信信号ｓnp（１），ｓnp（２），・・・，ｓnp（Ｎd）から遅延時間推定区間検出部４５により検出された遅延時間推定区間＃ｈの部分信号ｓnp,h（１），ｓnp,h（２），・・・，ｓnp,h（Ｎ）を抽出し、その遅延時間推定区間＃ｈの部分信号ｓnp,h（１），ｓnp,h（２），・・・，ｓnp,h（Ｎ）にデシメーション処理を施して、デシメーション処理信号ｘh,np（１），ｘh,np（２），・・・，ｘh,np（Ｍ）を出力する。ただし、１≦ｎp≦Ｎpである。なお、デシメーション処理部４６−１〜４６−Ｎpはデシメーション手段を構成している。
遅延時間推定部４７はデシメーション処理部４６−１〜４６−Ｎpによるデシメーション処理信号ｘh,np（１），ｘh,np（２），・・・，ｘh,np（Ｍ）に超分解能遅延時間推定処理を施して、遅延時間推定区間＃ｈに存在している目標信号の遅延時間を推定する。なお、遅延時間推定部４７は遅延時間推定手段を構成している。
【００８０】
図１６はこの発明の実施の形態８による遅延時間推定装置のデシメーション処理部４６−ｎp（１≦ｎp≦Ｎp）を示す構成図である。
図において、信号抽出部５１はメモリ回路４３に時間方向に並べて蓄積されている受信信号ｓnp（１），ｓnp（２），・・・，ｓnp（Ｎd）から遅延時間推定区間検出部４５により検出された遅延時間推定区間＃ｈの部分信号ｓnp,h（１），ｓnp,h（２），・・・，ｓnp,h（Ｎ）を抽出する。ただし、１≦ｎp≦Ｎpである。
受信信号スペクトル生成部５２は信号抽出部５１により抽出された遅延時間推定区間＃ｈの部分信号ｓnp,h（１），ｓnp,h（２），・・・，ｓnp,h（Ｎ）にＦＦＴ処理を施して受信信号スペクトルｙnp,h（１），ｙnp,h（２），・・・，ｙnp,h（Ｎ）を生成する。ただし、１≦ｎp≦Ｎpである。
【００８１】
除算信号生成部５３は受信信号スペクトル生成部５２により生成された受信信号スペクトルｙnp,h（１），ｙnp,h（２），・・・，ｙnp,h（Ｎ）を送信波スペクトルΓ（１），Γ（２），・・・，Γ（Ｎ）で除算して除算信号ｘ'np,h（１），ｘ'np,h（２），・・・，ｘ'np,h（Ｎ）を生成する。ただし、１≦ｎp≦Ｎpである。
デシメーションフィルタ５４は除算信号生成部５３により生成された除算信号ｘ'np,h（１），ｘ'np,h（２），・・・，ｘ'np,h（Ｎ）にデシメーション処理を施して、デシメーション処理信号ｘnp,h（１），ｘnp,h（２），・・・，ｘnp,h（Ｍ）を出力する。ただし、１≦ｎp≦Ｎpである。
【００８２】
図１７はこの発明の実施の形態８による遅延時間推定装置の遅延時間推定部４７を示す構成図であり、図において、図３と同一符号は同一または相当部分を示すので説明を省略する。
パルス方向平均型相関行列生成部２５はデシメーション処理部４６−ｎp（１≦ｎp≦Ｎp）によるデシメーション処理信号ｘnp,h（１），ｘnp,h（２），・・・，ｘnp,h（Ｍ）（１≦ｎp≦Ｎp）間の相関を表す相関行列Ｒを生成する。
【００８３】
次に動作について説明する。
送信機１が電波をＮp回繰り返し生成し、送信アンテナ２が送信機１により生成された電波を目標３に向けてＮp回繰り返し送信する。
受信機５は、送信アンテナ２から送信されて、目標３に反射された電波を受信アンテナ４が受信すると、受信アンテナ４の受信波に帯域制限や位相検波を施して受信信号を生成する。
Ａ／Ｄ変換器６は、受信機５が受信信号を生成すると、その受信信号をアナログ／ディジタル変換して、ディジタルの受信信号ｓnp（１），ｓnp（２），・・・，ｓnp（Ｎd）（１≦ｎp≦Ｎp）をメモリ回路４３に出力する。
【００８４】
メモリ回路４３は、Ａ／Ｄ変換器６からディジタルの受信信号ｓnp（１），ｓnp（２），・・・，ｓnp（Ｎd）（１≦ｎp≦Ｎp）を受ける毎に、その受信信号をパルス方向に並べて蓄積するとともに、その受信信号を時間方向に並べて蓄積する。
なお、メモリ回路４３は、Ｎp回の送受信が終了した時点で、パルス方向に並べて蓄積している受信信号ｓ1（ｎd），ｓ2（ｎd），・・・，ｓNp（ｎd）をパルス方向ＦＦＴ部４４−１〜４４−Ｎdに出力する。ただし、１≦ｎd≦Ｎd、１≦ｎp≦Ｎpである。
【００８５】
パルス方向ＦＦＴ部４４−１〜４４−Ｎdは、メモリ回路４３からパルス方向に並べられている受信信号ｓ1（ｎd），ｓ2（ｎd），・・・，ｓNp（ｎd）を受けると、その受信信号ｓ1（ｎd），ｓ2（ｎd），・・・，ｓNp（ｎd）にＦＦＴ処理を施して、パルス方向ＦＦＴ処理信号ｚ1（ｎd），ｚ2（ｎd），・・・，ｚNp（ｎd）を出力する。
遅延時間推定区間検出部４５は、パルス方向ＦＦＴ部４４−１〜４４−Ｎdからパルス方向ＦＦＴ処理信号ｚ1（ｎd），ｚ2（ｎd），・・・，ｚNp（ｎd）を受けると、そのパルス方向ＦＦＴ処理信号ｚ1（ｎd），ｚ2（ｎd），・・・，ｚNp（ｎd）を解析して、目標信号が存在している遅延時間推定区間＃ｈを検出する。
遅延時間推定区間検出部４５の具体的な処理内容は以下の通りである。
【００８６】
遅延時間推定区間検出部４５は、パルス方向ＦＦＴ部４４−１〜４４−Ｎdからパルス方向ＦＦＴ処理信号ｚ1（ｎd），ｚ2（ｎd），・・・，ｚNp（ｎd）を受けると、下記の式（１３）に示すように、そのパルス方向ＦＦＴ処理信号ｚ1（ｎd），ｚ2（ｎd），・・・，ｚNp（ｎd）の電力Ｐnp（ｌ）を算出する。
【数９】
2008107242.tif 000011
【００８７】
遅延時間推定区間検出部４５は、パルス方向ＦＦＴ処理信号の電力Ｐnp（ｌ）を算出すると、そのパルス方向ＦＦＴ処理信号の電力Ｐnp（ｌ）と予め設定された基準電力Ｔｈrdを比較する。
遅延時間推定区間検出部４５は、パルス方向ＦＦＴ処理信号の電力Ｐnp（ｌ）が基準電力Ｔｈrdより小さければ、当該遅延時間推定区間に目標信号が存在していないと判断するが、パルス方向ＦＦＴ処理信号の電力Ｐnp（ｌ）が基準電力Ｔｈrdより大きければ、当該遅延時間推定区間に目標信号が存在していると判断する。
【００８８】
遅延時間推定区間検出部４５は、例えば、ｎp＝ｎp0、ｌ＝ｈのときのパルス方向ＦＦＴ処理信号の電力Ｐnp0（ｈ）が基準電力Ｔｈrdより大きい場合、パルス方向ＦＦＴ処理信号ｚnp0（１＋Ｍ（ｈ−１）），ｚnp0（２＋Ｍ（ｈ−１）），・・・，ｚnp0（Ｎ＋Ｍ（ｈ−１））に目標信号の成分が含まれているとして、下記の式（１４）によりドップラー周波数ｆdを求め、そのドップラー周波数ｆdを出力する。
【数１０】
2008107242.tif 000012
また、遅延時間推定区間検出部４５は、遅延時間推定区間＃ｈに目標信号が存在している旨をデシメーション処理部４６−１〜４６−Ｎpに通知する。
【００８９】
デシメーション処理部４６−１〜４６−Ｎpは、遅延時間推定区間検出部４５から遅延時間推定区間＃ｈに目標信号が存在している旨の通知を受けると、メモリ回路４３から時間方向に並べられている受信信号ｓnp（１），ｓnp（２），・・・，ｓnp（Ｎd）を入力して、その受信信号ｓnp（１），ｓnp（２），・・・，ｓnp（Ｎd）から遅延時間推定区間＃ｈの部分信号ｓnp,h（１），ｓnp,h（２），・・・，ｓnp,h（Ｎ）を抽出する。ただし、１≦ｎp≦Ｎpである。
即ち、デシメーション処理部４６−１〜４６−Ｎpの信号抽出部５１は、下記の式（１５）に示すように、時間方向に並べられている受信信号ｓnp（１），ｓnp（２），・・・，ｓnp（Ｎd）から遅延時間推定区間＃ｈの部分信号ｓnp,h（１），ｓnp,h（２），・・・，ｓnp,h（Ｎ）を抽出する（図１８を参照）。
【数１１】
2008107242.tif 000013
【００９０】
デシメーション処理部４６−１〜４６−Ｎpの受信信号スペクトル生成部５２は、信号抽出部５１が遅延時間推定区間＃ｈの部分信号ｓnp,h（１），ｓnp,h（２），・・・，ｓnp,h（Ｎ）を抽出すると、その遅延時間推定区間＃ｈの部分信号ｓnp,h（１），ｓnp,h（２），・・・，ｓnp,h（Ｎ）にＦＦＴ処理を施して受信信号スペクトルｙnp,h（１），ｙnp,h（２），・・・，ｙnp,h（Ｎ）を生成する。ただし、１≦ｎp≦Ｎpである。
デシメーション処理部４６−１〜４６−Ｎpの除算信号生成部５３は、受信信号スペクトル生成部５２が受信信号スペクトルｙnp,h（１），ｙnp,h（２），・・・，ｙnp,h（Ｎ）を生成すると、その受信信号スペクトルｙnp,h（１），ｙnp,h（２），・・・，ｙnp,h（Ｎ）を送信波スペクトルΓ（１），Γ（２），・・・，Γ（Ｎ）で除算して除算信号ｘ'np,h（１），ｘ'np,h（２），・・・，ｘ'np,h（Ｎ）を生成する。ただし、１≦ｎp≦Ｎpである。
デシメーション処理部４６−１〜４６−Ｎpのデシメーションフィルタ５４は、除算信号生成部５３が除算信号ｘ'np,h（１），ｘ'np,h（２），・・・，ｘ'np,h（Ｎ）を生成すると、その除算信号ｘ'np,h（１），ｘ'np,h（２），・・・，ｘ'np,h（Ｎ）にデシメーション処理を施して、デシメーション処理信号ｘnp,h（１），ｘnp,h（２），・・・，ｘnp,h（Ｍ）を出力する。ただし、１≦ｎp≦Ｎpである。
【００９１】
遅延時間推定部４７は、デシメーション処理部４６−１〜４６−Ｎpからデシメーション処理信号ｘnp,h（１），ｘnp,h（２），・・・，ｘnp,h（Ｍ）を受けると、そのデシメーション処理信号ｘh,np（１），ｘh,np（２），・・・，ｘh,np（Ｍ）に超分解能遅延時間推定処理を施して、遅延時間推定区間＃ｈに存在している目標信号の遅延時間を推定する。
遅延時間推定部４７の具体的な処理内容は以下の通りである。
【００９２】
遅延時間推定部４７のパルス方向平均型相関行列生成部２５は、デシメーション処理部４６−１〜４６−Ｎpからデシメーション処理信号ｘnp,h（１），ｘnp,h（２），・・・，ｘnp,h（Ｍ）（１≦ｎp≦Ｎp）を受けると、下記の式（１６）に示すように、デシメーション処理信号ｘnp,h（１），ｘnp,h（２），・・・，ｘnp,h（Ｍ）（１≦ｎp≦Ｎp）間の相関を表す相関行列Ｒを生成する。
【数１２】
2008107242.tif 000014
【００９３】
遅延時間推定部４７の固有ベクトル算出部２２は、パルス方向平均型相関行列生成部２５が相関行列Ｒを生成すると、その相関行列Ｒの固有値ｖe（１），ｖe（２），・・・，ｖe（Ｍ）（ｖe（１）＞ｖe（２）＞・・・＞ｖe（Ｍ））と、その固有値ｖe（ｉ）（１≦ｉ≦Ｍ）に対応する固有ベクトルｅ（ｉ）を算出する。
固有ベクトル算出部２２は、上記のようにして、固有値ｖe（ｉ）（１≦ｉ≦Ｍ）に対応する固有ベクトルｅ（ｉ）を算出すると、目標数をＫとして、固有ベクトルｅ（Ｋ＋１），ｅ（Ｋ＋２），・・・，ｅ（Ｍ）をＭＵＳＩＣ処理部２３に出力する。
【００９４】
遅延時間推定部４７のＭＵＳＩＣ処理部２３は、固有ベクトル算出部２２から固有ベクトルｅ（Ｋ＋１），ｅ（Ｋ＋２），・・・，ｅ（Ｍ）を受けると、その固有ベクトルｅ（Ｋ＋１），ｅ（Ｋ＋２），・・・，ｅ（Ｍ）を雑音空間とするＭＵＳＩＣ処理を実施して、遅延時間推定区間＃ｈに存在している目標信号の遅延時間を推定する。
【００９５】
以上で明らかなように、この実施の形態８によれば、目標信号が存在している遅延時間推定区間＃ｈに限定して、超分解能遅延時間推定処理を実施するように構成しているので、上記実施の形態１と同様の効果を奏する他に、低負荷で超分解能遅延時間推定処理を実施することができる効果を奏する。
【００９６】
実施の形態９．
図１９はこの発明の実施の形態９による遅延時間推定装置を示す構成図であり、図において、図１及び図１５と同一符号は同一または相当部分を示すので説明を省略する。
メモリ回路６１はＡ／Ｄ変換器６からディジタルの受信信号ｓnp（１），ｓnp（２），・・・，ｓnp（Ｎd）（１≦ｎp≦Ｎp）を受ける毎に、その受信信号をパルス方向に並べて蓄積する。なお、メモリ回路６１は受信信号蓄積手段を構成している。
遅延時間推定区間検出部６２は図１５の遅延時間推定区間検出部４５と同様の検出処理を実施して、ドップラー周波数ｆdと目標信号が存在している遅延時間推定区間＃ｈを出力する他に、その遅延時間推定区間＃ｈのパルス方向ＦＦＴ処理信号ｚnp0（１＋Ｍ（ｈ−１）），ｚnp0（２＋Ｍ（ｈ−１）），・・・，ｚnp0（Ｎ＋Ｍ（ｈ−１））を出力する。なお、遅延時間推定区間検出部６２は遅延時間推定区間検出手段を構成している。
【００９７】
受信信号スペクトル生成部６３は遅延時間推定区間検出部６２から出力された遅延時間推定区間＃ｈのパルス方向ＦＦＴ処理信号ｚnp0（１＋Ｍ（ｈ−１）），ｚnp0（２＋Ｍ（ｈ−１）），・・・，ｚnp0（Ｎ＋Ｍ（ｈ−１））にＦＦＴ処理を施して受信信号スペクトルｙh（１），ｙh（２），・・・，ｙh（Ｎ）を生成する。
除算信号生成部６４は受信信号スペクトル生成部６３により生成された受信信号スペクトルｙh（１），ｙh（２），・・・，ｙh（Ｎ）を送信波スペクトルΓ（１），Γ（２），・・・，Γ（Ｎ）で除算して除算信号ｘ'h（１），ｘ'h（２），・・・，ｘ'h（Ｎ）を生成する。
【００９８】
デシメーションフィルタ６５は除算信号生成部６４により生成された除算信号ｘ'h（１），ｘ'h（２），・・・，ｘ'h（Ｎ）にデシメーション処理を施して、デシメーション処理信号ｘh（１），ｘh（２），・・・，ｘh（Ｍ）を出力する。
なお、受信信号スペクトル生成部６３、除算信号生成部６４及びデシメーションフィルタ６５からデシメーション手段が構成されている。
遅延時間推定部６６はデシメーションフィルタ６５から出力されたデシメーション処理信号ｘh（１），ｘh（２），・・・，ｘh（Ｍ）に超分解能遅延時間推定処理を施して、遅延時間推定区間＃ｈに存在している目標信号の遅延時間を推定する。なお、遅延時間推定部６６は遅延時間推定手段を構成している。
【００９９】
次に動作について説明する。
送信機１が電波をＮp回繰り返し生成し、送信アンテナ２が送信機１により生成された電波を目標３に向けてＮp回繰り返し送信する。
受信機５は、送信アンテナ２から送信されて、目標３に反射された電波を受信アンテナ４が受信すると、受信アンテナ４の受信波に帯域制限や位相検波を施して受信信号を生成する。
Ａ／Ｄ変換器６は、受信機５が受信信号を生成すると、その受信信号をアナログ／ディジタル変換して、ディジタルの受信信号ｓnp（１），ｓnp（２），・・・，ｓnp（Ｎd）（１≦ｎp≦Ｎp）をメモリ回路４３に出力する。
【０１００】
メモリ回路６１は、Ａ／Ｄ変換器６からディジタルの受信信号ｓnp（１），ｓnp（２），・・・，ｓnp（Ｎd）（１≦ｎp≦Ｎp）を受ける毎に、その受信信号をパルス方向に並べて蓄積する。
なお、メモリ回路６１は、Ｎp回の送受信が終了した時点で、パルス方向に並べて蓄積している受信信号ｓ1（ｎd），ｓ2（ｎd），・・・，ｓNp（ｎd）をパルス方向ＦＦＴ部４４−１〜４４−Ｎdに出力する。ただし、１≦ｎd≦Ｎd、１≦ｎp≦Ｎpである。
【０１０１】
パルス方向ＦＦＴ部４４−１〜４４−Ｎdは、メモリ回路６１からパルス方向に並べられている受信信号ｓ1（ｎd），ｓ2（ｎd），・・・，ｓNp（ｎd）を受けると、上記実施の形態８と同様に、その受信信号ｓ1（ｎd），ｓ2（ｎd），・・・，ｓNp（ｎd）にＦＦＴ処理を施して、パルス方向ＦＦＴ処理信号ｚ1（ｎd），ｚ2（ｎd），・・・，ｚNp（ｎd）を出力する。
遅延時間推定区間検出部６２は、パルス方向ＦＦＴ部４４−１〜４４−Ｎdからパルス方向ＦＦＴ処理信号ｚ1（ｎd），ｚ2（ｎd），・・・，ｚNp（ｎd）を受けると、図１５の遅延時間推定区間検出部４５と同様の検出処理を実施して、目標信号が存在している遅延時間推定区間＃ｈを検出する。
遅延時間推定区間検出部６２は、図１５の遅延時間推定区間検出部４５と同様に、ドップラー周波数ｆdと目標信号が存在している遅延時間推定区間＃ｈを出力する他に、その遅延時間推定区間＃ｈのパルス方向ＦＦＴ処理信号ｚnp0（１＋Ｍ（ｈ−１）），ｚnp0（２＋Ｍ（ｈ−１）），・・・，ｚnp0（Ｎ＋Ｍ（ｈ−１））を受信信号スペクトル生成部６３に出力する。
【０１０２】
受信信号スペクトル生成部６３は、遅延時間推定区間検出部６２から遅延時間推定区間＃ｈのパルス方向ＦＦＴ処理信号ｚnp0（１＋Ｍ（ｈ−１）），ｚnp0（２＋Ｍ（ｈ−１）），・・・，ｚnp0（Ｎ＋Ｍ（ｈ−１））を受けると、そのパルス方向ＦＦＴ処理信号ｚnp0（１＋Ｍ（ｈ−１）），ｚnp0（２＋Ｍ（ｈ−１）），・・・，ｚnp0（Ｎ＋Ｍ（ｈ−１））にＦＦＴ処理を施して受信信号スペクトルｙh（１），ｙh（２），・・・，ｙh（Ｎ）を生成する。
除算信号生成部６４は、受信信号スペクトル生成部６３が受信信号スペクトルｙh（１），ｙh（２），・・・，ｙh（Ｎ）を生成すると、その受信信号スペクトルｙh（１），ｙh（２），・・・，ｙh（Ｎ）を送信波スペクトルΓ（１），Γ（２），・・・，Γ（Ｎ）で除算して除算信号ｘ'h（１），ｘ'h（２），・・・，ｘ'h（Ｎ）を生成する。
【０１０３】
デシメーションフィルタ６５は、除算信号生成部６４が除算信号ｘ'h（１），ｘ'h（２），・・・，ｘ'h（Ｎ）を生成すると、その除算信号ｘ'h（１），ｘ'h（２），・・・，ｘ'h（Ｎ）にデシメーション処理を施して、デシメーション処理信号ｘh（１），ｘh（２），・・・，ｘh（Ｍ）を出力する。
遅延時間推定部６６は、デシメーションフィルタ６５からデシメーション処理信号ｘh（１），ｘh（２），・・・，ｘh（Ｍ）を受けると、そのデシメーション処理信号ｘh（１），ｘh（２），・・・，ｘh（Ｍ）に超分解能遅延時間推定処理を施して、遅延時間推定区間＃ｈに存在している目標信号の遅延時間を推定する。
【０１０４】
以上で明らかなように、この実施の形態９によれば、目標信号が存在している遅延時間推定区間＃ｈに限定して、超分解能遅延時間推定処理を実施するように構成しているので、上記実施の形態１と同様の効果を奏する他に、低負荷で超分解能遅延時間推定処理を実施することができる効果を奏する。
【０１０５】
実施の形態１０．
図２０はこの発明の実施の形態１０による遅延時間推定装置の遅延時間推定部８−ｌ（１≦ｌ≦Ｌ）を示す構成図であり、この発明の実施の形態１０による遅延時間推定装置の全体構成は、図１の遅延時間推定装置と同一である。
図２０において、図３と同一符号は同一または相当部分を示すので説明を省略する。
ＥＳＰＲＩＴ処理部２６は固有ベクトル算出部２２により算出された固有ベクトルｅl（Ｋ＋１），ｅl（Ｋ＋２），・・・，ｅl（Ｍd）を用いて、低負荷の超分解能推定処理であるＥＳＰＲＩＴ（Ｅｓｔｉｍａｔｉｏｎ  ｏｆ  Ｓｉｇｎａｌ  Ｐａｒａｍｅｔｅｒｓ  ｖｉａ  Ｒｏｔａｔｉｏｎａｌ  ＩｎｖａｒｉａｎｃｅＴｅｃｈｎｉｑｕｅｓ）処理を実施して、目標信号の遅延時間を推定する。
【０１０６】
上記実施の形態１では、遅延時間推定部８−ｌ（１≦ｌ≦Ｌ）のＭＵＳＩＣ処理部２３が固有ベクトルｅl（Ｋ＋１），ｅl（Ｋ＋２），・・・，ｅl（Ｍd）を雑音空間とするＭＵＳＩＣ処理を実施して、遅延時間推定区間＃ｌに含まれている目標信号の遅延時間τを推定するものについて示したが、図２０に示すように、ＥＳＰＲＩＴ処理部２６が固有ベクトルｅl（Ｋ＋１），ｅl（Ｋ＋２），・・・，ｅl（Ｍd）を用いて、ＥＳＰＲＩＴ処理を実施することにより、目標信号の遅延時間τを推定するようにしてもよい。
【０１０７】
即ち、遅延時間推定部８−ｌ（１≦ｌ≦Ｌ）のＥＳＰＲＩＴ処理部２６は、固有ベクトル算出部２２が固有ベクトルｅl（Ｋ＋１），ｅl（Ｋ＋２），・・・，ｅl（Ｍd）を算出すると、下記の式（１７）に示すように、相関行列Ｒlの固有値の中で、値が大きい上位Ｋ種類の固有値に対応する固有ベクトルｅl（１），ｅl（２），・・・，ｅl（Ｋl）から行列Ｅlを生成する。
Ｅl＝［ｅl（１）・・・ｅl（Ｋl）］                          （１７）
【０１０８】
ＥＳＰＲＩＴ処理部２６は、固有ベクトルｅl（１），ｅl（２），・・・，ｅl（Ｋl）から行列Ｅlを生成すると、下記の式（１８）に示すように、その行列Ｅlから行列Ψlを算出する。
【数１３】
2008107242.tif 000015
式（１８）において、行列Ｅ'lHは行列Ｅ'lのエルミート共役を表している。
また、行列Ｊ1は（Ｍd−１）行Ｍd列の行列であり、Ｊ1（ｊ，ｋ）は行列Ｊ1のｉ行ｋ列成分を表している。
また、行列Ｊ2は（Ｍd−１）行Ｍd列の行列であり、Ｊ2（ｊ，ｋ）は行列Ｊ2のｉ行ｋ列成分を表している。
【０１０９】
ＥＳＰＲＩＴ処理部２６は、行列Ｅlから行列Ψlを算出すると、下記の式（１９）によって、目標信号の遅延時間τk（ｌ）ハットを求める。
【数１４】
2008107242.tif 000016
式（１９）において、ｖΨ（ｋl）は行列Ψlのｋl番目の固有値、ａｒｇ［ｖΨ（ｋl）］は固有値ｖΨ（ｋl）の偏角を表している。
【０１１０】
以上で明らかなように、この実施の形態１０によれば、ＥＳＰＲＩＴ処理を実施して、目標信号の遅延時間を推定するように構成したので、遅延時間推定部８−ｌ（１≦ｌ≦Ｌ）における超分解能推定処理の負荷を軽減することができる効果を奏する。
【０１１１】
実施の形態１１．
図２１はこの発明の実施の形態１１による遅延時間推定装置の遅延時間推定部８−ｌ（１≦ｌ≦Ｌ）を示す構成図であり、この発明の実施の形態１１による遅延時間推定装置の全体構成は、図１の遅延時間推定装置と同一である。
図２１において、図３と同一符号は同一または相当部分を示すので説明を省略する。
最尤推定処理部２７は固有ベクトル算出部２２により算出された固有ベクトルｅl（Ｋ＋１），ｅl（Ｋ＋２），・・・，ｅl（Ｍd）を用いて最尤推定処理を実施して、目標信号の遅延時間を推定する。
【０１１２】
上記実施の形態１では、遅延時間推定部８−ｌ（１≦ｌ≦Ｌ）のＭＵＳＩＣ処理部２３が固有ベクトルｅl（Ｋ＋１），ｅl（Ｋ＋２），・・・，ｅl（Ｍd）を雑音空間とするＭＵＳＩＣ処理を実施して、遅延時間推定区間＃ｌに含まれている目標信号の遅延時間τを推定するものについて示したが、図２１に示すように、最尤推定処理部２７が固有ベクトルｅl（Ｋ＋１），ｅl（Ｋ＋２），・・・，ｅl（Ｍd）を用いて、最尤推定処理を実施することにより、目標信号の遅延時間τを推定するようにしてもよい。
【０１１３】
即ち、遅延時間推定部８−ｌ（１≦ｌ≦Ｌ）の最尤推定処理部２７は、固有ベクトル算出部２２が固有ベクトルｅl（Ｋ＋１），ｅl（Ｋ＋２），・・・，ｅl（Ｍd）を算出すると、下記の式（２０）に示すように、目標信号の遅延時間τ1，τ2，・・・，τKlを変数とする評価関数Θ（τ1，τ2，・・・，τKl）を設定する。
【数１５】
2008107242.tif 000017
式（２０）において、ａ（τ）は、上記の式（３）で定義しているステアリングベクトルを表している。
【０１１４】
最尤推定処理部２７は、評価関数Θ（τ1，τ2，・・・，τKl）を設定すると、その評価関数Θ（τ1，τ2，・・・，τKl）を最大とするτ1＝τ1（ｌ）ハット，・・・，τKl＝τK（ｌ）ハットを求める。
なお、推定区間ｌに含まれる目標ｋの遅延時間の推定値はτk（ｌ）ハットである。
【０１１５】
以上で明らかなように、この実施の形態１１によれば、最尤推定処理を実施して、目標信号の遅延時間を推定するように構成したので、遅延時間推定部８−ｌ（１≦ｌ≦Ｌ）における超分解能推定処理の負荷を軽減することができる効果を奏する。
【０１１６】
実施の形態１２．
図２２はこの発明の実施の形態１２による遅延時間推定装置の遅延時間推定部８−ｌ（１≦ｌ≦Ｌ）を示す構成図であり、この発明の実施の形態１２による遅延時間推定装置の全体構成は、図１の遅延時間推定装置と同一である。
図２２において、図３と同一符号は同一または相当部分を示すので説明を省略する。
信号蓄積部であるメモリ回路２８はデシメーション処理部７から出力されたデシメーション処理信号ｘl（１），ｘl（２），・・・，ｘl（Ｎd）（１≦ｌ≦Ｌ）を蓄積するとともに、隣接信号相関型目標信号判定部２９から目標信号が有る旨の判定結果を受けると、そのデシメーション処理信号ｘl（１），ｘl（２），・・・，ｘl（Ｎd）（１≦ｌ≦Ｌ）を相関行列生成部２１に出力する。
有無判定部である隣接信号相関型目標信号判定部２９はメモリ回路２８に蓄積されているデシメーション処理信号のうち、隣接しているデシメーション処理信号間の相関値を基準にして、目標信号の有無を判定する。
【０１１７】
次に動作について説明する。
遅延時間推定部８−ｌ（１≦ｌ≦Ｌ）のメモリ回路２８は、上記実施の形態１と同様にして、デシメーション処理部７からデシメーション処理信号ｘl（１），ｘl（２），・・・，ｘl（Ｎd）（１≦ｌ≦Ｌ）を受けると、そのデシメーション処理信号ｘl（１），ｘl（２），・・・，ｘl（Ｎd）を蓄積する。
【０１１８】
遅延時間推定部８−ｌ（１≦ｌ≦Ｌ）の隣接信号相関型目標信号判定部２９は、メモリ回路２８に蓄積されているデシメーション処理信号ｘl（１），ｘl（２），・・・，ｘl（Ｎd）のうち、隣接しているデシメーション処理信号間の相関値を基準にして、目標信号の有無を判定する。
即ち、隣接信号相関型目標信号判定部２９は、下記の式（２１）に示すように、隣接しているデシメーション処理信号間の相関値ｃｏｒ（ｌ）を算出する。
【数１６】
2008107242.tif 000018
式（２１）において、ｘ（ｉ）*はｘ（ｉ）の複素共役を表している。
【０１１９】
隣接信号相関型目標信号判定部２９は、隣接しているデシメーション処理信号間の相関値ｃｏｒ（ｌ）を算出すると、その相関値ｃｏｒ（ｌ）と予め設定されているスレッショルドを比較する。
隣接信号相関型目標信号判定部２９は、相関値ｃｏｒ（ｌ）がスレッショルドより小さければ、目標信号が無いものと判定して、その判定結果をメモリ回路２８に出力する。
一方、相関値ｃｏｒ（ｌ）がスレッショルドより大きければ、目標信号が有るものと判定して、その判定結果をメモリ回路２８に出力する。
【０１２０】
遅延時間推定部８−ｌ（１≦ｌ≦Ｌ）のメモリ回路２８は、隣接信号相関型目標信号判定部２９から目標信号が有る旨の判定結果を受けると、蓄積しているデシメーション処理信号ｘl（１），ｘl（２），・・・，ｘl（Ｍ）（１≦ｌ≦Ｌ）を相関行列生成部２１に出力する。
相関行列生成部２１以降の処理内容は、上記実施の形態１と同様であるため説明を省略する。
【０１２１】
以上で明らかなように、この実施の形態１２によれば、目標信号の有無を判定して、目標信号が有る場合に限り、ＭＵＳＩＣ処理を実施して、目標信号の遅延時間を推定するように構成したので、遅延時間推定部８−ｌ（１≦ｌ≦Ｌ）における超分解能推定処理の負荷を軽減することができる効果を奏する。
【図面の簡単な説明】
【０１２２】
【図１】この発明の実施の形態１による遅延時間推定装置を示す構成図である。
【図２】この発明の実施の形態１による遅延時間推定装置のデシメーション処理部７を示す構成図である。
【図３】この発明の実施の形態１による遅延時間推定装置の遅延時間推定部８−ｌ（１≦ｌ≦Ｌ）を示す構成図である。
【図４】デシメーション処理部７の処理内容を示す説明図である。
【図５】デシメーションフィルタ１４−ｌ（１≦ｌ≦Ｌ）の内部を示す構成図である。
【図６】この発明の実施の形態２による遅延時間推定装置のデシメーション処理部７を示す構成図である。
【図７】デシメーション処理部７における推定区間切り出し型デシメーション部１５−ｌ（１≦ｌ≦Ｌ）の内部を示す構成図である。
【図８】この発明の実施の形態３による遅延時間推定装置のデシメーション処理部７を示す構成図である。
【図９】この発明の実施の形態４による遅延時間推定装置の遅延時間推定部８−ｌ（１≦ｌ≦Ｌ）を示す構成図である。
【図１０】この発明の実施の形態５による遅延時間推定装置を示す構成図である。
【図１１】この発明の実施の形態５による遅延時間推定装置のデシメーション処理部９を示す構成図である。
【図１２】この発明の実施の形態６による遅延時間推定装置を示す構成図である。
【図１３】この発明の実施の形態７による遅延時間推定装置を示す構成図である。
【図１４】マルチパスフェージング環境下において、Ａ／Ｄ変換器６における通常のサンプリングと、オーバーサンプリング型Ａ／Ｄ変換器４２におけるオーバーサンプリングとの相違を示す説明図である。
【図１５】この発明の実施の形態８による遅延時間推定装置を示す構成図である。
【図１６】この発明の実施の形態８による遅延時間推定装置のデシメーション処理部４６−ｎp（１≦ｎp≦Ｎp）を示す構成図である。
【図１７】この発明の実施の形態８による遅延時間推定装置の遅延時間推定部４７を示す構成図である。
【図１８】デシメーション処理部４６−１〜４６−Ｎpの処理内容を示す説明図である。
【図１９】この発明の実施の形態９による遅延時間推定装置を示す構成図である。
【図２０】この発明の実施の形態１０による遅延時間推定装置の遅延時間推定部８−ｌ（１≦ｌ≦Ｌ）を示す構成図である。
【図２１】この発明の実施の形態１１による遅延時間推定装置の遅延時間推定部８−ｌ（１≦ｌ≦Ｌ）を示す構成図である。
【図２２】この発明の実施の形態１２による遅延時間推定装置の遅延時間推定部８−ｌ（１≦ｌ≦Ｌ）を示す構成図である。
【符号の説明】
【０１２３】
１  送信機（電波送信手段）、２  送信アンテナ（電波送信手段）、３  目標、４，４−１，４−２  受信アンテナ（電波受信手段）、５，５−１，５−２  受信機（電波受信手段）、６，６−１，６−２  Ａ／Ｄ変換器（電波受信手段）、７  デシメーション処理部（デシメーション手段）、８−１〜８−Ｌ  遅延時間推定部（遅延時間推定手段）、９  デシメーション処理部（デシメーション手段）、１０  遅延時間推定部（遅延時間推定手段）、１１  受信信号分割部、１２−１〜１２−Ｌ  受信信号スペクトル生成部、１３−１〜１３−Ｌ  除算信号生成部、１４−１〜１４−Ｌ  デシメーションフィルタ、１５−１〜１５−Ｌ  推定区間切り出し型デシメーション部、１６−１〜１６−Ｌ  零割防止型除算信号生成部、１７  目標信号検出部、１８  除算信号生成部、１９  デシメーションフィルタ、２１  相関行列生成部、２２  固有ベクトル算出部、２３  ＭＵＳＩＣ処理部、２４  目標信号判定型固有ベクトル算出部（目標判別部）、２５  パルス方向平均型相関行列生成部、２６  ＥＳＰＲＩＴ処理部、２７  最尤推定処理部、２８  メモリ回路（信号蓄積部）、２９  隣接信号相関型目標信号判定部（有無判定部）、３１  ＩＦＦＴ処理部（逆フーリエ変換処理部）、３２  推定区間切り出し処理部（信号切出処理部）、３３  ＦＦＴ処理部（フーリエ変換処理部）、４１  最大比合成処理部（最大比合成処理手段）、４２  オーバーサンプリング型Ａ／Ｄ変換器（電波受信手段）、４３  メモリ回路（受信信号蓄積手段）、４４−１〜４４−Ｎd  パルス方向ＦＦＴ部（フーリエ変換処理手段）、４５  遅延時間推定区間検出部（遅延時間推定区間検出手段）、４６−１〜４６−Ｎp  デシメーション処理部（デシメーション手段）、４７  遅延時間推定部（遅延時間推定手段）、５１  信号抽出部、５２  受信信号スペクトル生成部、５３  除算信号生成部、５４  デシメーションフィルタ、６１  メモリ回路（受信信号蓄積手段）、６２  遅延時間推定区間検出部（遅延時間推定区間検出手段）、６３  受信信号スペクトル生成部（デシメーション手段）、６４  除算信号生成部（デシメーション手段）、６５  デシメーションフィルタ（デシメーション手段）。
【図１】
2008107242.tif 000019
【図２】
2008107242.tif 000020
【図３】
2008107242.tif 000021
【図４】
2008107242.tif 000022
【図５】
2008107242.tif 000023
【図６】
2008107242.tif 000024
【図７】
2008107242.tif 000025
【図８】
2008107242.tif 000026
【図９】
2008107242.tif 000027
【図１０】
2008107242.tif 000028
【図１１】
2008107242.tif 000029
【図１２】
2008107242.tif 000030
【図１３】
2008107242.tif 000031
【図１４】
2008107242.tif 000032
【図１５】
2008107242.tif 000033
【図１６】
2008107242.tif 000034
【図１７】
2008107242.tif 000035
【図１８】
2008107242.tif 000036
【図１９】
2008107242.tif 000037
【図２０】
2008107242.tif 000038
【図２１】
2008107242.tif 000039
【図２２】
2008107242.tif 000040