JP3216561B2 - 方向探知装置 - Google Patents

方向探知装置

Info

Publication number
JP3216561B2
JP3216561B2 JP07018297A JP7018297A JP3216561B2 JP 3216561 B2 JP3216561 B2 JP 3216561B2 JP 07018297 A JP07018297 A JP 07018297A JP 7018297 A JP7018297 A JP 7018297A JP 3216561 B2 JP3216561 B2 JP 3216561B2
Authority
JP
Japan
Prior art keywords
calculator
projection length
covariance matrix
target
mode vector
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
JP07018297A
Other languages
English (en)
Other versions
JPH10268016A (ja
Inventor
由昌 大橋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Mitsubishi Electric Corp
Original Assignee
Mitsubishi Electric Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Mitsubishi Electric Corp filed Critical Mitsubishi Electric Corp
Priority to JP07018297A priority Critical patent/JP3216561B2/ja
Publication of JPH10268016A publication Critical patent/JPH10268016A/ja
Application granted granted Critical
Publication of JP3216561B2 publication Critical patent/JP3216561B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Description

【発明の詳細な説明】
【0001】
【発明の属する技術分野】この発明は目標の位置する方
向を高い測角精度で求める方向探知装置に関するもので
ある。
【0002】
【従来の技術】図2は従来の方向探知装置を示すもの
で、この図はRalph O.Schmidt著「Mu
ltiple Emitter Location a
nd Signal Parameter Estim
ation」IEEE Trans.on Anten
na and Propagation,vol. A
P−34,no.3,March 1986 および大
橋、比嘉、佐藤、瀬尾、柳沢、立木著「スーパーレゾリ
ューションESPRITアルゴリズムのレーダへの応用
の検討」電子情報通信学会、宇宙航空エレクトロニクス
研究会、SANE97、2月、1997年に記述されて
いる内容にもとづいた方向探知装置の構成図である。
【0003】図2において、1は送信パルスの発生およ
び受信サンプリングのタイミングを指示するトリガパル
スを発生する制御器、2は制御器1で発生されたトリガ
パルスにもとづいて送信パルス幅τの送信パルスを発生
する送信機、3は送信機2で発生された送信パルスを電
波として放射する送信アンテナ、4−1,4−2,・・
・,4−Mは目標で反射された電波を受信する同一仕様
のM個の受信アンテナ、12−1,12−2,・・・,
12−Mはそれぞれ上記受信アンテナ4−1,4−2,
・・・,4−Mで受信された電波に対し増幅、周波数変
換を行うと共に制御器からのトリガパルスにもとづいて
サンプリングする周波数帯域幅1/τの同一仕様の狭帯
域受信機、6はM個の上記狭帯域受信機12−1,12
−2,・・・,12−Mの出力より共分散行列を算出す
る共分散行列算出器、7は上記共分散行列より固有値、
固有ベクトルを算出する固有値・固有ベクトル算出器、
8は目標を仮想したとき予想されるモードベクトルを算
出するモードベクトル算出器、9は固有値・固有ベクト
ル算出器7で算出された固有値および固有ベクトルにも
とづいて張られる雑音部分空間にモードベクトル算出器
8で算出したモードベクトルを射影して、その射影長を
算出する雑音部分空間射影長算出器、10は仮想目標を
空間上で走査して雑音部分空間射影長算出器9より出力
される射影長が極小になるときの仮想目標の位置する方
向を算出する射影長極小値算出器、11は狭帯域受信機
12、共分散行列算出器6、固有値・固有ベクトル算出
器7、モードベクトル算出器8、雑音部分空間射影長算
出器9、射影長極小値算出器10から構成される信号処
理装置である。
【0004】図3において13は目標の例として航空機
を示したものである。
【0005】次に動作について説明する。制御器1より
出力される繰返し周波数fPRF のトリガパルスにもとづ
いて送信機2よりRF(Radio Frequenc
y)f0 で、かつ送信パルス幅τの送信パルスが送信さ
れ、送信アンテナ3より電波として空中に放射される。
【0006】送信アンテナ3より放射された電波が目標
に照射され、目標で反射された上記電波は受信アンテナ
4−1,4−2,・・・,4−Mで受信される。
【0007】受信アンテナ4−1,4−2,・・・,4
−Mで受信された上記目標からの電波はそれぞれ狭帯域
受信機12−1,12−2,・・・,12−Mで増幅、
周波数変換された後制御器1からのトリガパルスに同期
してサンプリングされる。
【0008】M個の狭帯域受信機12の出力はそれぞれ
共分散行列算出器6に入力される。ここで、M個の狭帯
域受信機12の出力をS1 (k),S2 (k),・・
・,SM (K)とするとき入力ベクトルは数1で表され
る。ただし、kはサンプリング番号を表わす。
【0009】
【数1】
【0010】共分散行列を求めるにあたり、P回のスナ
ップショットをP個のパルスヒットにより行うものとす
れば共分散行列算出器6では数2に示す共分散行列が算
出される。
【0011】
【数2】
【0012】ただし、*は複素共役またはエルミート共
役を表わす。
【0013】固有値・固有ベクトル算出器7では、得ら
れた共分散行列にもとづいてM個の固有値が算出される
と共に、それぞれの固有値に対応して固有ベクトルが算
出され、これらの固有値および固有ベクトルは雑音部分
空間射影長算出器9に入力される。上記の固有値をλ
1 ,λ2 ,・・・,λM とし、対応する固有ベクトルを
1 ,X2 ,・・・,XM とすればこれらの固有値およ
び固有ベクトルは数2を用いて数3により与えられる。
【0014】
【数3】
【0015】得られた共分散行列は正定値行列であるこ
とから固有値はすべて零より大きくなる。狭帯域受信機
12の受信機雑音は互いに無相関で、かつ、その標準偏
差σはすべて等しく、さらに、目標からの信号と受信機
雑音が無相関のとき数4が成り立つ。
【0016】
【数4】
【0017】固有値λ1 に対応する固有ベクトルをX1
とし、λ2 ,λ3 ,・・・,λM 対応する固有ベクトル
をX2 ,X3 ,・・・,XM とすればX1 にもとづく信
号部分空間とX2 ,X3 ,・・・,XM により張られる
雑音部分空間は互いに直交補空間となる。
【0018】さて、モードベクトル算出器8にはM個の
受信アンテナ4が配列されているとき、ある角度θの方
向に目標が存在すると仮定した場合におけるM個の受信
アンテナ4にもとづくデータ、即ち、モードベクトルが
蓄えられており、通常、この角度θは所定の範囲にわた
っている。
【0019】従って、モードベクトル算出器8よりある
角度θにおけるモードベクトルが発生され、雑音部分空
間射影長算出器9に入力される。一方、雑音部分空間射
影長算出器9では、入力された固有値、固有ベクトルに
もとづいて雑音部分空間が張られ、モードベクトル発生
器8から出力されたモードベクトルが雑音部分空間に射
影される。このモードベクトルをa(θ)とおけば数5
で与えられる射影長Ra (θ)が雑音部分空間射影長算
出器9より出力される。
【0020】
【数5】
【0021】射影長極小値算出器10ではモードベクト
ルa(θ)の射影長Ra (θ)に対しその極小値を与え
るθの値θ、が求められる。
【0022】θ1 は目標の存在する方向の角度推定値と
なる。
【0023】
【発明が解決しようとする課題】従来の方向探知装置は
以上のように構成されているのでP回のスナップショッ
ト数を得るにはP個のパルスヒットが必要となり、目標
の角度推定値を得るのに時間が長くなる問題があった。
【0024】この発明は上記のような問題点を解決する
ためになされたもので、P回のスナップショット数を得
るのにPより少ないパルスヒットにより目標の角度推定
値を得ることを目標とする。
【0025】
【課題を解決するための手段】第1の発明による方向探
知装置はトリガパルスを発生する制御器と、上記トリガ
パルスにもとづいて送信パルス幅τの送信パルスを発生
する送信機と、上記送信パルスを電波として放射する送
信アンテナと、目標で反射された電波を受信する同一仕
様のM個の受信アンテナと、上記受信アンテナで受信さ
れた受信波を増幅、周波数変換および複数サンプリング
する周波数帯域2/τ以上の同一仕様のM個の広帯域受
信機と、上記広帯域受信機からの出力にもとづいて共分
散行列を算出する共分散行列算出器と、上記共分散行列
算出器より出力される共分散行列にもとづいて固有値お
よび固有ベクトルを算出する固有値・固有ベクトル算出
器と、上記受信アンテナの配列にもとづき仮想目標に対
応するモードベクトルを算出するモードベクトル算出器
と、上記固有値・固有ベクトル算出器より出力される固
有値および固有ベクトルにもとづいて張られる雑音部分
空間に上記モードベクトル算出器より出力されるモード
ベクトルを射影し、その射影長を算出する雑音部分空間
射影長算出器と、上記雑音部分空間射影長算出器より出
力される上記射影長が極小となるときの上記仮想目標の
位置する方向を求め上記方向を出力する射影長極小値算
出器を備えるものである。
【0026】
【発明の実施の形態】
実施の形態1.図1はこの発明の実施の形態1を示す構
成図であり、図において、1はトリガパルスを発生する
制御器、2は送信パルス幅τを送信する送信機、3は送
信パルスを放射する送信アンテナ、4−1,4−2,・
・・,4−Mは目標で反射された電波を受信する受信ア
ンテナ、5−1,5−2,・・・,5−Mは受信波を増
幅、周波数変換および複数サンプリングする周波数帯域
2/τ以上の同一仕様の広帯域受信機、6は共分散行列
を算出する共分散行列算出器、7は共分散行列より固有
値および固有ベクトルを算出する固有値・固有ベクトル
算出器、8は受信アンテナ4の配列にもとづき仮想目標
に対応するモードベクトルを算出するモードベクトル算
出器、9は固有値・固有ベクトル算出器7より出力され
る固有値および固有ベクトルにもとづいて張られる雑音
部分空間にモードベクトル算出器より出力されるモード
ベクトルを射影し、その射影長を算出する雑音部分空間
射影長算出器、10は雑音部分空間射影長算出器9より
出力される上記射影長が極小となるときの仮想目標の位
置する方向を算出する射影長極小値算出器である。
【0027】次に動作について説明する。制御器1より
出力される繰返し周波数fPRF のトリガパルスにもとづ
いて送信機2よりRFf0 で、かつ送信パルス幅での送
信パルスが送信され、送信アンテナ3より電波として空
中に放射される。
【0028】送信アンテナ3より放射された電波が目標
に照射され、目標で反射された上記電波は受信アンテナ
4−1,4−2,・・・,4−Mで受信される。
【0029】受信アンテナ4−1,4−2,・・・,4
−Mで受信された上記目標からの信号はそれぞれ広帯域
受信機5−1,5−2,・・・,5−Mで増幅、周波数
変換された後制御器1からのトリガパルスに同期して複
数回サンプリングされる。
【0030】例えば、広帯域受信機5の周波数帯域幅を
送信パルス幅τに対しN/τとする。ただし、Nは2以
上の整数である。
【0031】このとき、広帯域受信機5それぞれの受信
機雑音は時間軸上でτ/Nごとに無相関となる。従っ
て、目標からの反射波である1パルス内をτ/Nごとに
N回サンプリングする。
【0032】上記のサンプリングを目標からの反射波の
L個の受信パルスに対し行えばM個の広帯域受信機5の
受信機雑音は互いに無相関であり、かつこれらの受信機
雑音と目標信号は無相関であることからスナップショッ
ト数はNLとなる。
【0033】M個の広帯域受信機5によるスナップショ
ット数NLのデータはそれぞれ共分散行列算出器6に入
力される。M個の広帯域受信機5の出力をS1 (k),
2 (k),・・・,SM (k)とするとき入力ベクト
ルは数1で表され、上記のスナップショット数NLに対
応する共分散行列は数6で表される。
【0034】
【数6】
【0035】共分散行列算出器6では数6に示す共分散
行列が算出され、固有値・固有ベクトル算出器7に出力
される。
【0036】固有値・固有ベクトル算出器7では、得ら
れた共分散行列にもとづいてM個の固有値が算出される
と共に、それぞれの固有値に対応して固有ベクトルが算
出され、これらの固有値および固有ベクトルは雑音部分
空間射影長算出器9に入力される。上記の固有値をλ
1 、λ2 ,・・・,λM とし、対応する固有ベクトルを
1 ,X2 ,・・・,XM とすればこれらの固有値およ
び固有ベクトルは数6を用いて数3により与えられる。
【0037】得られた共分散行列は正定値行列であるこ
とから固有値はすべて零より大きくなる。広帯域受信機
5の受信機雑音は互いに無相関で、かつその標準偏差を
σとすれば数4が成り立つ。ただし、信号と受信機雑音
は無相関と仮定している。
【0038】固有値λ1 に対応する固有ベクトルをX1
とし、λ2 ,λ3 ,・・・,λM に対応する固有ベクト
ルをX2 、X3 ,・・・,XM とすればX1 にもとづく
信号部分空間とX2 ,X3 ,・・・,XM により張られ
る雑音部分空間は互いに直交補空間となる。
【0039】モードベクトル算出器8にはM個の受信ア
ンテナ4が配列されているので、ある角度θの方向に目
標が存在すると仮定した場合におけるM個の受信アンテ
ナ4にもとづくデータ、即ち、モードベクトルが蓄えら
れており、通常、この角度θは所定の範囲にわたってい
る。
【0040】従って、モードベクトル算出器8よりある
角度θにおけるモードベクトルが発生され、雑音部分空
間射影長算出器9に入力される。一方、雑音部分空間射
影長算出器9では、入力された固有値、固有ベクトルに
もとづいて雑音部分空間が張られ、モードベクトル発生
器8から出力されたモードベクトルが雑音部分空間に射
影される。このモードベクトルをa(θ)とおけば数5
で与えられる射影長が雑音部分空間射影長算出器9より
出力される。
【0041】射影長極小値算出器10ではモードベクト
ルa(θ)の射影長に対しその極小値を与えるθの値θ
1 が求められ、このθ1 が目標の存在する方向の角度推
定値となる。
【0042】
【発明の効果】第1の発明によれば、送信パルスの帯域
に比べ広い帯域を持つ広帯域受信機を用いることにより
目標に対するパルスヒットの数より多いスナップショッ
ト数を得ることができ、目標を測角するときに必要な目
標の観測時間を短縮することができる。
【図面の簡単な説明】
【図1】 この発明による方向探知装置の実施の形態1
を示す図である。
【図2】 従来の方向探知装置の構成を示す図である。
【図3】 方向探知装置の一般的な運用例を示す図であ
る。
【符号の説明】 1 制御器、2 送信機、3 送信アンテナ、4 受信
アンテナ、5 広帯域受信機、6 共分散行列算出器、
7 固有値・固有ベクトル算出器、8 モードベクトル
算出器、9 雑音部分空間射影長算出器、10 射影長
極小値算出器、11 信号処理装置、12 狭帯域受信
機、13 目標。
───────────────────────────────────────────────────── フロントページの続き (58)調査した分野(Int.Cl.7,DB名) G01S 3/00 - 3/74 G01S 13/00 - 13/95

Claims (1)

    (57)【特許請求の範囲】
  1. 【請求項1】 送信トリガパルスを発生する制御器と、
    上記制御器からのトリガパルスにもとづいて送信パルス
    幅τを発生する送信機と、上記送信機からの送信パルス
    幅を電波として放射する送信アンテナと、目標で反射さ
    れた上記送信アンテナからの電波を受信する同一仕様の
    複数の受信アンテナと、上記複数の受信アンテナそれぞ
    れに対応して接続されて上記制御器からのトリガパルス
    に同期して、上記目標からの信号を増幅、周波数変換お
    よびN回(N≧2)サンプリングする周波数帯域/τ
    の同一仕様の複数の広帯域受信機と、上記広帯域受信機
    からの出力にもとづいて共分散行列を算出する共分散行
    列算出器と、上記共分散行列算出器より出力される共分
    散行列にもとづいて固有値および固有ベクトルを算出す
    る固有値・固有ベクトル算出器と、上記受信アンテナの
    配列にもとづき、かつ仮想目標の位置する方向の関数で
    表されるモードベクトルを算出するモードベクトル算出
    器と、上記固有値・固有ベクトル算出器より出力される
    固有値および固有ベクトルにもづいて張られる雑音部分
    空間に上記モードベクトル算出器より出力されるモード
    ベクトルを射影し、その射影長を算出する雑音部分空間
    射影長算出器と、上記雑音部分空間射影長算出器より出
    力される上記射影長が極小となるときの上記仮想目標の
    位置する方向を求め上記方向を出力する射影長極小値算
    出器とで構成したことを特徴とする方向探知装置。
JP07018297A 1997-03-24 1997-03-24 方向探知装置 Expired - Fee Related JP3216561B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP07018297A JP3216561B2 (ja) 1997-03-24 1997-03-24 方向探知装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP07018297A JP3216561B2 (ja) 1997-03-24 1997-03-24 方向探知装置

Publications (2)

Publication Number Publication Date
JPH10268016A JPH10268016A (ja) 1998-10-09
JP3216561B2 true JP3216561B2 (ja) 2001-10-09

Family

ID=13424135

Family Applications (1)

Application Number Title Priority Date Filing Date
JP07018297A Expired - Fee Related JP3216561B2 (ja) 1997-03-24 1997-03-24 方向探知装置

Country Status (1)

Country Link
JP (1) JP3216561B2 (ja)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5701106B2 (ja) * 2011-03-04 2015-04-15 富士通テン株式会社 レーダ装置及び該レーダ装置の到来角算出方法

Also Published As

Publication number Publication date
JPH10268016A (ja) 1998-10-09

Similar Documents

Publication Publication Date Title
US10866304B1 (en) Signal detection and denoising systems
JP5617334B2 (ja) レーダ装置及び目標探知方法
JP5631763B2 (ja) Mimoレーダシステム、送信装置、受信装置及びmimoレーダ信号処理方法
JP5600866B2 (ja) 探知測距装置および探知測距方法
US11269070B2 (en) Radar apparatus
EP2200384B1 (en) Ranging diversity-reception method and receiver
US10523288B2 (en) Wireless receiver, wireless receiving method, and wireless system
US11454702B2 (en) Synthetic aperture radar method and synthetic aperture radar device
US20220003835A1 (en) Radar apparatus and method for determining range side lobe
US6446025B1 (en) Multiple propagation wave parameter measuring method and apparatus and machine-readable recording medium recording multiple propagation wave parameter measuring program
JP4494233B2 (ja) 到来時刻推定装置
US6549762B1 (en) Method for estimating arrival direction of desired wave
JP3216561B2 (ja) 方向探知装置
WO2023056991A2 (zh) Mimo传感器及到达角近似度判断方法和目标信息匹配方法
JPH11133130A (ja) 電波の到来方向の検出装置及びfmーcwレーダ
JP6164936B2 (ja) レーダ装置
JP3484995B2 (ja) 瞬時パッシブ距離測定装置
US20220321164A1 (en) Angle information estimation of ultra-wideband wireless signals
JP2937068B2 (ja) 方位検出装置
JP7413850B2 (ja) 物体位置の角度推定装置及び方法、並びにレーダ装置
US20240134006A1 (en) Mimo sensor, method for determining direction-of-arrival approximation degree, and target information matching method
JP6989663B1 (ja) レーダ装置
JP7396128B2 (ja) アンテナ装置の較正装置、較正システム、及び較正方法
WO2023210451A1 (ja) 測位システムおよびそれを備える車両、並びに測位方法
JP3220899B2 (ja) 方向探知装置

Legal Events

Date Code Title Description
LAPS Cancellation because of no payment of annual fees