JP4001745B2 - Target detection device - Google Patents
Target detection device Download PDFInfo
- Publication number
- JP4001745B2 JP4001745B2 JP2001397919A JP2001397919A JP4001745B2 JP 4001745 B2 JP4001745 B2 JP 4001745B2 JP 2001397919 A JP2001397919 A JP 2001397919A JP 2001397919 A JP2001397919 A JP 2001397919A JP 4001745 B2 JP4001745 B2 JP 4001745B2
- Authority
- JP
- Japan
- Prior art keywords
- integration
- circuit
- target
- range bin
- path
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Description
【0001】
【発明の属する技術分野】
この発明はレーダ受信信号から目標検出処理を行う目標検出装置に関するものである。
【0002】
【従来の技術】
図17は従来の目標検出装置の構成を示すブロック図である。図において、1は送信信号を変調し変調信号を出力する変調回路、2は変調回路1からの変調信号により送信ビームを形成する送信機、3は送信器2からの送信ビームを電波として発射する送信アンテナ、4は送信された電波が目標に反射して生成される目標信号を受信する受信アンテナ、5は受信された目標信号に対して帯域制限、位相検波、増幅を行いアナログ信号を出力する受信機、6は受信機5からのアナログ信号をディジタル信号に変換するA/D変換器、7は変調回路1からの変調信号とA/D変換器6からのディジタル信号を入力しレンジビン毎の復調信号を出力する復調回路である。
【0003】
また、図17において、8は復調回路7からのレンジビン毎の復調信号を検波して、レンジビン毎の検波信号を出力する検波回路、24は検波回路8からのレンジビン毎の検波信号を入力し、初期レンジビンを基点とする目標が移動すると想定される積分路を設定し、設定された積分路におけるインコヒーレント積分処理を行う移動量補償積分回路、10は移動量補償積分回路24によるインコヒーレント積分処理結果と、雑音を誤って目標と判定する誤警報確率を基準に予め設定されている固定スレッショルドとを比較して目標信号を判定する警報回路、11は警報回路10により判定された目標信号の追尾処理を行い目標の位置を推定する追尾回路である。
【0004】
図18は復調回路7の内部構成を示すブロック図である。図において、25は変調回路1からの変調信号を入力してレンジビン毎に所定時間ずつ遅延させた参照信号を生成し、A/D変換器6からのディジタル信号と生成した参照信号とを乗算して乗算結果を出力する参照信号乗算回路、26は参照信号乗算回路25から出力されたレンジビン毎の信号についてコヒーレント積分を行い、レンジビン毎に復調信号を出力する復調フィルタである。
【0005】
次に動作について説明する。
変調回路1は送信信号を変調して変調信号を出力し、送信機2は変調回路1からの変調信号により送信ビームを形成し、送信アンテナ3は送信機2からの送信ビームを電波として発射する。送信信号をs0 (t)、変調周期をT[s]とすると、送信信号s0 (t)(t:時刻)は次の(1)式で示される。
s0 (t)=s0 (t+T) (1)
【0006】
受信アンテナ4は送信された電波が目標に反射して生成される目標信号を受信し、受信機5は受信された目標信号に対して帯域制限、位相検波、増幅を行いアナログ信号を出力する。また、受信機5により受信信号の受信機雑音が重畳される。A/D変換器6は受信機5からのアナログ信号をディジタル信号に変換し、復調回路7は変調回路1からの変調信号とA/D変換器6からのディジタル信号を入力し、コヒーレント積分を行ってレンジビン毎の復調信号を出力する。
【0007】
A/D変換器6のサンプリング周期をT/Nとすると、復調回路7は一般的に次の(2)式に従って復調処理を行う。
【数1】
上記(2)式において、Nは復調回路7において受信信号を復調するのに使用される信号点数、すなわちレンジビン数、iは信号sの番号、s* (t)はs(t)の複素共役、kは信号sと信号s* の時間的なずれを示し、コヒーレント積分を行った場合、k=0のときのみ1となり、k≠0のときには零となる。
【0008】
図18に示すレンジビン毎に構成された参照信号乗算回路25−1〜25−Nは、変調回路1から送信された変調信号と同一の変調信号を入力して、それぞれ所定時間ずつ遅延させた参照信号を生成し、A/D変換器6からのディジタル信号と生成した参照信号とを乗算して乗算結果を出力する。
【0009】
図19は復調回路7における復調処理を説明するタイミングチャートであり、ここでは、送信信号と受信信号に関するタイミングチャートと、送信信号と各参照信号の対応関係を示している。復調回路7における各参照信号乗算回路25−1〜25−Nは、図19に示すように、変調回路1からの変調信号により、それぞれサンプリング周期T/Nずつ遅らせた参照信号1〜Nを生成する。図19では目標がjレンジビンに存在し、電波を送信してからjT/N[s]だけ遅れて受信アンテナ4により受信されている状況を表している。
【0010】
このとき、jレンジビンのkサンプリング目の参照信号sr,j (k)は次の(3)式により示される。
【数2】
ここで、s0 (t)は時刻tの送信信号を示している。
【0011】
jレンジビンに存在する目標信号は、jレンジビンの参照信号により復調される。すなわち、参照信号乗算回路25−jは、次の(4)式によりA/D変換器6からのディジタル信号とjレンジビンの参照信号とを複素乗算する。
【数3】
ここで、s1 (t)は時刻tの受信信号、sr,j (k)はjレンジビンのkサンプリング目の参照信号、s2,j (k)は参照信号乗算回路25−jの出力信号を示している。
【0012】
そして、復調フィルタ26−jは、次の(5)式によりコヒーレント積分を行いjレンジビンの復調信号s3,j を出力する。
【数4】
【0013】
jレンジビンに目標が存在した場合には、jレンジビンの復調信号s3,j に目標信号成分が積分されている。復調回路7からのレンジビン毎の復調信号は検波回路8に出力される。検波回路8は入力したレンジビン毎の復調信号を検波し、レンジビン毎の検波信号が移動量補償積分回路24に出力される。移動量補償積分回路24は、検波回路8からのレンジビン毎の検波信号を入力し、初期レンジビンを基点とする目標が移動すると想定される積分路を設定し、設定された積分路におけるインコヒーレント積分処理を行う。
【0014】
移動量補償積分回路24は、以上の処理を1CPI(Coherent Processing Interval)単位として、あらかじめ設定しておいたM回行う。この時点で移動量補償積分回路24には、i(1≦i≦M)CPIにおけるj(1≦j≦N)レンジビンの復調信号si,j が2次元データとして蓄えられている。目標信号は2次元領域上を移動する。
【0015】
図20は移動量補償積分回路24の入力信号の状況を説明する図である。図20ではM=N=3のときに、1CPIあたり1レンジビン移動する目標信号成分の変化を表している。図20ではsi,j をiCPI目におけるjレンジビンのセルとして表している。斜線部分が目標信号成分を表している。
【0016】
図21は移動量補償積分処理を説明するフローチャートである。図21に示す処理手順では、レーダ方向に接近している、すなわち、CPIと伴にレンジビン番号が減少する目標の検出を対象としている。ステップST11において、CPI数M、レンジビン数Nをそれぞれ設定する。また、レンジビンの補償間隔δrを設定する。この場合、目標はレーダ方向に接近しているので、レンジビンの補償間隔δrは負である。
【0017】
ステップST12において、1CPI目において目標信号の存在を仮定する初期レンジビンrnをNに設定する。ステップST13において、目標のレンジ移動量drの初期値を0と設定し、レンジビンの補償間隔δrを単位とした移動レンジjaを1と設定する。
【0018】
ステップST14において、次の(6)式を用いてiCPIに目標信号成分の存在するレンジr0 (i)を計算する。
【数5】
ここで、rnは1CPI目において目標信号の存在を仮定する初期レンジビン、drはMCPIで移動する目標のレンジビン移動量を示している。また、レンジr0 (i)を含むレンジビンc(i)を選択する。この場合、計算されたレンジr0 (i)の値の小数点以下は、例えば四捨五入してレンジビンc(i)を選択する。
【0019】
ステップST15において、次の(7)式を用いて、MCPIでdrレンジビンだけ移動する目標に関するインコヒーレント積分PDI(Post Detection Integration)の値P(rn,ja)を計算する。
【数6】
ここで、pc(i)はc(i)レンジビンの検波出力値を示している。
【0020】
ステップST16において、MCPIあたりの目標のレンジビン移動量drをδrだけ増加する。ステップST17において、MCPIにおける目標のレンジ移動量drがrnより大きいときは、次のステップST19に移行する。また、そうでない場合はステップST18において、移動レンジjaをインクリメントして、上記ステップST14からST17までの処理を繰り返す。
【0021】
ステップST19において、目標の初期レンジビンrnを1だけ減少し、ステップST20において、目標の初期レンジビンrn<1の場合は処理を終了し、そうでない場合は上記ステップST13からST20までの処理を繰り返す。
【0022】
移動量補償積分回路24からはP(rn,ja)(1≦rn≦N,1≦ja≦N/δr)の2次元データが警報回路10に出力される。警報回路10では、M点のインコヒーレント積分処理を想定し、固定スレッショルドを設定している。すなわち、雑音についてM点のインコヒーレント積分処理を行った値が、固定スレッショルドを越える確率を誤警報確率とし、その誤警報確率が所望値となるように固定スレッショルドが設定されている。
【0023】
警報回路10は移動量補償積分回路24から出力されたP(rn,ja)の各成分と固定スレッショルドを比較し、固定スレッショルドを越えた成分を目標信号と判定する。P(rn0,j0)が固定スレッショルドを越えた場合、初期レンジビンrn0からMCPIの間にj0 δrレンジだけ移動する目標が存在すると判定される。
【0024】
以上の処理は各アジマス角、エレベーション角毎に行われ、警報回路10での検出目標については、アジマス角、エレベーション角、距離(検出目標の1CPI目における位置を距離とする)の3次元のデータが得られる。この処理はSPI(Signal Processing Interval)単位で行われ、1SPI毎にこの3次元の観測データは追尾回路11に出力されて追尾処理が行われる。
【0025】
追尾回路11は等速度運動を行っている目標を想定して追尾処理を行う。アジマス角、エレベーション角、距離の3次元データから、次の(8)式を用いてx,y,zの3次元データに変換する。
x=RcosElcosAz
y=RcosElsinAz
z=RsinEl (8)
ここで、Rは目標距離、Elはエレベーション角、Azはアジマス角を示している。
【0026】
また、等速度運動を想定した目標運動モデルは、次の(9)式の状態方程式により示される。
【数7】
【0027】
また、目標信号の観測モデルは次の(10)式により示される。観測雑音は、角度(アジマス角、エレベーション角)の観測誤差、距離の観測誤差をx−y−z座標系における観測誤差に変換して計算される。
【数8】
【0028】
上記(10)式では、各SPIにおいて目標の位置が観測されるモデルを表している。k−1SPI目の平滑値xハットk−1とその誤差共分散行列Pハットk−1が与えられたとき、次の(11)式を用いてkSPI目の予測値xチルダkを計算することで、kSPI目における目標の予測位置が計算される。
【数9】
【0029】
また、次の(12)式を用いて、xチルダkの誤差共分散行列Pチルダkを計算する。Pチルダkは6行6列の行列となっている。
【数10】
【0030】
次の(13)式を用いてkSPI目におけるカルマンゲインKk を計算する。
【数11】
ここで、Rvは観測雑音分散(E[vv* ])、Tは転置行列を示している。
【0031】
次の(14)式を用いて、kSPI目におけるxの平滑値xハットkを計算する。
【数12】
【0032】
そして、次の(15)式を用いて、xハットkの誤差共分散行列Pハットkを計算する。
【数13】
ここで、Iは単位行列を示している。
【0033】
また、k+1→kとして、(11)式を用いて、次のSPIの予測値を計算する。以降、同様の手順により各SPIにおける目標の予測位置とその誤差分散が計算される。実際の目標が(9)式で表される運動モデルに従って運動しており、(10)式で表される観測モデルに従って目標信号が観測されるとき、k−1SPIまでの観測データを用い、カルマンフィルタ処理を行った結果、kSPI目における目標観測位置zk の確率分布は、次の(16)式で表される確率密度関数ps (zk )の正規分布となる。
【数14】
ここで、Tは転置行列を示している。
【0034】
追尾回路11は、各SPI毎に目標位置と速度の平滑値、平滑誤差共分散行列、目標位置と速度の予測値、予測誤差共分散行列、目標信号電力を出力する。
【0035】
図22は特開平10−68771号公報に開示された従来の目標検出装置の構成を示すブロック図である。図において、27は入力信号についてコヒーレント積分処理を行い、時間領域の信号成分からドップラー周波数領域の信号成分に変換するコヒーレント積分回路、28は目標のドップラー周波数の移動(加速度)を考慮してインコヒーレント積分を行う加速度補償積分回路である。その他の構成は図17に示す構成と同等のものである。
【0036】
図22に示す信号処理系では、ドップラー周波数の変化する(加速度を有する)目標のドップラー速度を補償しながら、インコヒーレント積分を行っている。復調回路7からの復調信号はコヒーレント積分回路27により、ドップラー周波数領域の信号成分に変換される。検波回路8からの検波信号は時間とドップラー周波数の2次元データとなり、目標信号はこの2次元領域上を移動する。
【0037】
加速度補償積分回路28では、時間とドップラー周波数の2次元データに対して図21に示すフローチャートに従って(時間領域をドップラー周波数領域に置き換えて)インコヒーレント積分処理を行う。その結果、加速度を有する目標についてもインコヒーレント積分され目標検出性能が改善される。
【0038】
また、従来の目標検出装置として、特開平8−179037号公報に開示されたレーダ装置がある。これは、目標のレンジビン移動を考慮して、インコヒーレント処理を行うことで目標検出性能の改善を図るものである。
【0039】
【発明が解決しようとする課題】
一般に移動量補償積分回路24におけるインコヒーレント積分路数が増加すると誤警報確率が増加するが、従来のレーダ装置は以上のように構成され、移動量補償積分回路24におけるインコヒーレント積分路数に無関係に、警報回路10では積分路数を1としたときに発生する誤警報確率を基に固定スレッショルドを設定しているため、警報回路10において所望の誤警報確率とならないという課題があった。
【0040】
また、特開平10−68771号公報や特開平8−179037号公報に開示された目標検出装置では、積分路数の増加に伴う誤警報確率の増加を考慮した固定スレッショルドの設定については具体的に言及されていない。
【0041】
この発明は上記のような課題を解決するためになされたもので、インコヒーレント積分路数が増加しても、移動量補償積分回路において初期レンジビンに対応した固定スレッショルドを設定し、所望の誤警報確率を達成できる目標検出装置を得ることを目的とする。
【0042】
【課題を解決するための手段】
この発明に係る目標検出装置は、目標に反射して生成された目標信号をレンジビン毎に復調する復調回路と、この復調回路からのレンジビン毎の復調信号を検波する検波回路と、この検波回路からのレンジビン毎の検波信号を入力し、初期レンジビンを起点とする目標が移動すると想定される積分路を設定し、設定された積分路におけるインコヒーレント積分処理を行い、設定された積分路を基に所望の誤警報確率が得られるように、雑音の上記検波回路出力時の確率密度関数、積分路設定を行った際のセルに通じるパスを構成するセルの検波出力値の総和、セルに通じる全パス数及び初期セルの番号により、初期レンジビンに対応した固定スレッショルドを計算する移動量補償積分回路と、この移動量補償積分回路によるインコヒーレント積分処理結果と上記初期レンジビンに対応した固定スレッショルドを比較して目標信号を判定する警報回路とを備えたものである。
【0043】
この発明に係る目標検出装置は、移動量補償積分回路が初期レンジビンに対応した固定スレッショルドを、各積分路の積分結果が無相関であると近似させることで簡易的に計算するものである。
【0044】
この発明に係る目標検出装置は、警報回路により判定された目標信号の追尾処理を行い目標の位置を推定する追尾回路を備え、移動量補償積分回路が、上記追尾回路からの目標観測位置の確率分布に基づき、目標存在確率の高い領域を選択して積分路を設定するものである。
【0045】
この発明に係る目標検出装置は、移動量補償積分回路が、追尾回路からの目標観測位置の確率分布に基づき、目標存在確率の低い領域の積分路を棄却し、棄却した積分路の演算負荷を目標存在確率の高い領域の積分路における積分点数増加に転用することで、積分路数と積分点数との間の演算負荷を調整して、目標存在確率の高い領域の積分路と積分点数を設定するものである。
【0046】
この発明に係る目標検出装置は、目標に反射して生成された目標信号をレンジビン毎に復調する復調回路と、この復調回路からのレンジビン毎の復調信号を検波する検波回路と、この検波回路からのレンジビン毎の検波信号を入力し、所望の誤警報確率が得られるようにレンジビン数とレンジビンの補償間隔により積分路数を決定して、初期レンジビンを起点とする目標が移動すると想定される積分路を設定し、設定された積分路におけるインコヒーレント積分処理を行い、設定された積分路を基に初期レンジビンに対応した固定スレッショルドを計算する移動量補償積分回路と、この移動量補償積分回路によるインコヒーレント積分処理結果と上記初期レンジビンに対応した固定スレッショルドを比較して目標信号を判定する警報回路とを備えたものである。
【0047】
この発明に係る目標検出装置は、警報回路により判定された目標信号の追尾処理を行い目標の位置を推定する追尾回路を備え、移動量補償積分回路が、所望の誤警報確率が得られるように、かつ、上記追尾回路からの目標観測位置の確率分布に基づき、所望の目標検出確率が得られるように積分路数を決定するものである。
【0048】
【発明の実施の形態】
以下、この発明の実施の一形態を説明する。
実施の形態1.
図1はこの発明の実施の形態1による目標検出装置の構成を示すブロック図である。図において、9は検波回路8からのレンジビン毎の検波信号を入力し、初期レンジビンを起点とする目標が移動すると想定される積分路を設定し、設定された積分路におけるインコヒーレント積分処理を行い、設定された積分路を基に所望の誤警報確率が得られるように、初期レンジビンに対応した固定スレッショルドを計算する移動量補償積分回路、10は移動量補償積分回路9によるインコヒーレント積分処理結果と初期レンジビンに対応した固定スレッショルドを比較して目標信号を判定する警報回路である。その他の構成は従来の図17に示す構成と同等である。
【0049】
図2は移動量補償積分回路9の内部構成を示すブロック図である。図において、12は検波回路8からのレンジビン毎の検波信号を入力し、初期レンジビンを起点とする目標が移動すると想定される積分路を設定し、設定された積分路を基に所望の誤警報確率が得られるように、初期レンジビンに対応した固定スレッショルドを計算して警報回路10に出力する積分路設定回路、13は積分路設定回路12により設定された積分路におけるインコヒーレント積分処理を行い、そのインコヒーレント積分処理結果を警報回路10に出力する積分回路である。
【0050】
次に動作について説明する。
変調回路1、送信機2、送信アンテナ3、受信アンテナ4、受信機5、A/D変換器6、復調回路7、検波回路8は、図17に示す従来例と同様に動作し、検波回路8からのレンジビン毎の検波信号が移動量補償積分回路9に出力される。移動量補償積分回路9における積分路設定回路12は、検波回路8からのレンジビン毎の検波信号を入力し、初期レンジビンを起点とする目標が移動すると想定される積分路を設定して積分回路13に出力すると共に、設定された積分路を基に所望の誤警報確率が得られるように、初期レンジビンに対応した固定スレッショルドを計算して警報回路10に出力する。
【0051】
このとき、積分路設定回路12は、次の(17)式を用いて固定スレッショルドThを計算する。
【数15】
ここで、Nc は初期セルの番号、pn (x)は雑音の検波回路出力時の確率密度関数、si,j はセルci に通じるj番目のパス、p(si,j )はパスsi,j を構成するセルの検波出力値の総和、imax はセルci に通じる全パス数を示している。
【0052】
図3は移動量補償積分回路9におけるの入力信号の状況を説明する図であり、図3(a)はデータ例、図3(b)は積分路設定例、図3(c)は積分路セルの集合を示しており、3CPIの間、最大3レンジビン移動する目標を想定した場合のセルの設定例について示している。
【0053】
また、図4は誤警報確率計算に使用するパスの例を説明する図であり、図4(a)はセルc1 に通じるパス、図4(b)はセルc4 に通じるパス、図4(c)はセルc2 に通じるパスをそれぞれ示しており、図3の積分路設定を行った際のセルに通じるパスの例について示している。
【0054】
図3(a)では3CPIと3レンジビンの範囲におけるセルを示している。1CPI目における初期セルを1レンジビンとし、レンジビン移動補償間隔δr=1としたとき、設定される積分路は図3(b)に示す4種類となる。また、4種類の積分路に含まれる全てのセルの集合は図3(c)のようになる。図3(c)では3CPI目の左側のセルから順に番号付けをしている。
【0055】
図4(a)に示すように、セルc1 に通じるパスは{c6 ,c4 }の1種類となる。また、図4(b)に示すように、セルc4 に通じるパスは{c6 }の1種類となる。さらに、図4(c)に示すように、セルc2 に通じるパスは{c6 ,c4 }と{c6 ,c5 }の2種類となる。
【0056】
以上の要領で各セルに通じるパスを調べ、上記(17)式を用いて誤警報確率Pfaとなる固定スレッショルドThを計算する。実際に、全てのパスを調べると、次の(18)式のようになる。
s1,1 ={c6 ,c4 }
s2,1 ={c6 ,c4 }
s2,2 ={c6 ,c5 }
s3,1 ={c6 ,c5 }
s4,1 ={c6 }
s5,1 ={c6 } (18)
【0057】
また、各パスの検波回路出力値に関する積分値を比較したときの最大値は、次の(19)式のように計算される。
ここで、xi はセルci の検波回路出力値を示している。
【0058】
上記(19)式を上記(17)式に代入すると、誤警報確率は次の(20)式により計算される。
【数16】
ここで、pn (x1 )〜pn (x6 )は、各セルに含まれる雑音成分電力x1 〜x6 の確率密度関数を示している。
【0059】
移動量補償積分回路9における積分路設定回路12は、設定された積分路を積分回路13に出力すると共に、計算された初期レンジビンに対応した固定スレッショルドThを警報回路10に出力する。積分回路13は積分路設定回路12で設定された積分路に沿ってインコヒーレント積分処理を行い、そのインコヒーレント積分処理結果を警報回路10に出力する。警報回路10は積分回路13からのインコヒーレント積分処理結果を入力し、初期レンジビンに対応した固定スレッショルドThに基づき目標信号を判定する。
【0060】
以上のように、この実施の形態1によれば、移動量補償積分回路9が検波回路8からのレンジビン毎の検波信号を入力し、初期レンジビンを起点とする目標が移動すると想定される積分路を設定し、設定された積分路におけるインコヒーレント積分処理を行い、設定された積分路を基に所望の誤警報確率が得られるように、初期レンジビンに対応した固定スレッショルドを計算し、警報回路10がインコヒーレント積分処理結果と初期レンジビンに対応した固定スレッショルドを比較して目標信号を判定することにより、インコヒーレント積分路数が増加しても、所望の誤警報確率を達成できるという効果が得られる。
【0061】
実施の形態2.
図5はこの発明の実施の形態2による目標検出装置の構成を示すブロック図である。図において、14は検波回路8からのレンジビン毎の検波信号を入力し、追尾回路11が上記(16)式で求めた目標観測位置zk の確率分布ps (zk )に基づき目標存在確率の高い領域に限定して、初期レンジビンを起点とする目標が移動すると想定される積分路を設定し、設定された積分路におけるインコヒーレント積分処理を行い、設定された積分路を基に所望の誤警報確率が得られるように、初期レンジビンに対応した固定スレッショルドを計算する移動量補償積分回路である。その他の構成は実施の形態1の図1に示す構成と同等である。
【0062】
図6は移動量補償積分回路14の内部構成を示すブロック図である。図において、15は検波回路8からのレンジビン毎の検波信号を入力し、追尾回路11が上記(16)式で求めた目標観測位置zk の確率分布ps (zk )に基づき目標存在確率の高い領域に限定して、初期レンジビンを起点とする目標が移動すると想定される積分路を設定し、設定された積分路を基に所望の誤警報確率が得られるように、初期レンジビンに対応した固定スレッショルドを計算する選択型積分路設定回路であり、積分回路13は実施の形態1の図2に示すものと同等である。
【0063】
次に動作について説明する。
変調回路1、送信機2、送信アンテナ3、受信アンテナ4、受信機5、A/D変換器6、復調回路7、検波回路8は、図17に示す従来例と同様に動作して、検波回路8からのレンジビン毎の検波信号が移動量補償積分回路14に出力される。
【0064】
移動量補償積分回路14における選択型積分路設定回路15は、検波回路8からのレンジビン毎の検波信号を入力し、追尾回路11が上記(16)式で求めた目標観測位置zk の確率分布ps (zk )に基づき目標存在確率の高い領域に限定して、初期レンジビンを起点とする目標が移動すると想定される積分路を設定して積分回路13に出力すると共に、設定された積分路を基に所望の誤警報確率が得られるように、初期レンジビンに対応した固定スレッショルドを計算して警報回路10に出力する。
【0065】
図7は目標存在範囲を説明する図であり、選択型積分路設定回路15において設定される移動量補償範囲について示している。ここでは、便宜的に2次元の状況について示している。図7に示す目標観測位置の確率分布の等高線は、上記(16)式の正規分布の等高線より計算する。受信ビーム幅と目標存在確率の等高線の関係が図7となるようなアジマス角、エレベーション角を検出し、その角度におけるΔRの距離範囲に限定して移動量補償を行い積分路を選択する。図7ではkσの範囲で移動量補償を行う状況について表している。ここで、σは目標存在確率の標準偏差、kは適宜設定された係数を示している。
【0066】
移動量補償積分回路14における積分回路13と、警報回路10、追尾回路11は実施の形態1と同様に動作する。
【0067】
以上のように、この実施の形態2によれば、移動量補償積分回路14が検波回路8からのレンジビン毎の検波信号を入力し、追尾回路11が求めた目標観測位置の確率分布に基づき目標存在確率の高い領域に限定して、初期レンジビンを起点とする目標が移動すると想定される積分路を設定し、設定された積分路におけるインコヒーレント積分処理を行い、設定された積分路を基に所望の誤警報確率が得られるように、初期レンジビンに対応した固定スレッショルドを計算し、警報回路10がインコヒーレント積分処理結果と初期レンジビンに対応した固定スレッショルドを比較して目標信号を判定することにより、インコヒーレント積分路数が増加しても、所望の誤警報確率を達成できるという効果が得られる。
【0068】
また、この実施の形態2によれば、移動量補償積分回路14が目標存在確率の高い領域に限定して積分路を設定することにより、積分路数が少なくなり、積分回路13におけるインコヒーレント積分の演算量を少なくすることができるという効果が得られる。
【0069】
実施の形態3.
図8はこの発明の実施の形態3による目標検出装置の構成を示すブロック図である。図において、16は検波回路8からのレンジビン毎の検波信号を入力し、所望の誤警報確率が得られるように積分路数を決定して、初期レンジビンを起点とする目標が移動すると想定される積分路を設定し、設定された積分路におけるインコヒーレント積分処理を行い、設定された積分路を基に初期レンジビンに対応した固定スレッショルドを計算する移動量補償積分回路である。その他の構成は実施の形態1の図1に示す構成と同等である。
【0070】
図9は移動量補償積分回路16の内部構成を示すブロック図である。図において、17は検波回路8からのレンジビン毎の検波信号を入力し、所望の誤警報確率が得られるように積分路数を決定して、初期レンジビンを起点とする目標が移動すると想定される積分路を設定し、設定された積分路を基に初期レンジビンに対応した固定スレッショルドを計算する誤警報確率調整型積分路設定回路で、積分回路13は実施の形態1の図2に示すものと同等である。
【0071】
次に動作について説明する。
変調回路1、送信機2、送信アンテナ3、受信アンテナ4、受信機5、A/D変換器6、復調回路7、検波回路8は、図17に示す従来例と同様に動作して、検波回路8からのレンジビン毎の検波信号が移動量補償積分回路16に出力される。
【0072】
移動量補償積分回路16における誤警報確率調整型積分路設定回路17は、検波回路8からのレンジビン毎の検波信号を入力し、所望の誤警報確率が得られるように積分路数を決定して、初期レンジビンを起点とする目標が移動すると想定される積分路を設定して積分回路13に出力すると共に、設定された積分路を基に初期レンジビンに対応した固定スレッショルドを計算して警報回路10に出力する。
【0073】
一般に、積分路数を多くすると誤警報確率が大きくなり、積分路数を少なくすると誤警報確率が小さくなる。ここで、レンジビン数をN、レンジビンの補償間隔をδrとすると、積分路数はN/δrで示される。そのため、誤警報確率調整型積分路設定回路17は、所望の誤警報確率が得られるように、レンジビン数Nとレンジビンの補償間隔δrにより積分路数を決定する。
【0074】
移動量補償積分回路16における積分回路13と、警報回路10、追尾回路11は実施の形態1と同様に動作する。
【0075】
以上のように、この実施の形態3によれば、移動量補償積分回路16が検波回路8からのレンジビン毎の検波信号を入力し、所望の誤警報確率が得られるように積分路数を決定して、初期レンジビンを起点とする目標が移動すると想定される積分路を設定し、設定された積分路におけるインコヒーレント積分処理を行い、設定された積分路を基に初期レンジビンに対応した固定スレッショルドを計算し、警報回路10がインコヒーレント積分処理結果と初期レンジビンに対応した固定スレッショルドを比較して目標信号を判定することにより、インコヒーレント積分路数が増加しても、所望の誤警報確率を達成できるという効果が得られる。
【0076】
実施の形態4.
図10はこの発明の実施の形態4による目標検出装置の構成を示すブロック図である。図において、18は検波回路8からのレンジビン毎の検波信号を入力し、所望の誤警報確率が得られるように、かつ、追尾回路11が上記(16)式で求めた目標観測位置zk の確率分布ps (zk )に基づき、所望の目標検出確率が得られるように積分路数を決定して、初期レンジビンを起点とする目標が移動すると想定される積分路を設定し、設定された積分路におけるインコヒーレント積分処理を行い、設定された積分路を基に初期レンジビンに対応した固定スレッショルドを計算する移動量補償積分回路である。その他の構成は実施の形態1の図1に示す構成と同等である。
【0077】
図11は移動量補償積分回路18の内部構成を示すブロック図である。図において、19は検波回路8からのレンジビン毎の検波信号を入力し、所望の誤警報確率が得られるように、かつ、追尾回路11が上記(16)式で求めた目標観測位置zk の確率分布ps (zk )に基づき、所望の目標検出確率が得られるように積分路数を決定して、初期レンジビンを起点とする目標が移動すると想定される積分路を設定し、設定された積分路を基に初期レンジビンに対応した固定スレッショルドを計算する誤警報確率/目標検出確率調整型積分路設定回路であり、積分回路13は実施の形態1の図2に示すものと同等である。
【0078】
次に動作について説明する。
変調回路1、送信機2、送信アンテナ3、受信アンテナ4、受信機5、A/D変換器6、復調回路7、検波回路8は、図17に示す従来例と同様に動作して、検波回路8からのレンジビン毎の検波信号が移動量補償積分回路18に出力される。
【0079】
移動量補償積分回路18における誤警報確率/目標検出確率調整型積分路設定回路19は、検波回路8からのレンジビン毎の検波信号を入力し、所望の誤警報確率が得られるように、かつ、追尾回路11が上記(16)式で求めた目標観測位置zk の確率分布ps (zk )に基づき、所望の目標検出確率が得られるように積分路数を決定して、初期レンジビンを起点とする目標が移動すると想定される積分路を設定して積分回路13に出力すると共に、設定された積分路を基に初期レンジビンに対応した固定スレッショルドを計算して警報回路10に出力する。
【0080】
一般に、積分路数を多くすると目標検出確率が大きくなり、積分路数を少なくすると目標検出確率が小さくなる。そこで、誤警報確率と同様に、レンジビン数Nとレンジビンの補償間隔δrにより積分路数を決定することで、目標検出確率を制御することができる。そのため、誤警報確率/目標検出確率調整型積分路設定回路19は、所望の誤警報確率が得られるように、かつ、追尾回路11が上記(16)式で求めた目標観測位置zk の確率分布ps (zk )に基づき、所望の目標検出確率が得られるように、レンジビン数Nとレンジビンの補償間隔δrにより積分路数を決定する。
【0081】
ここで、目標検出確率Pdは次の(21)式により計算される。
【数17】
ここで、rd はMCPIの間に移動したレンジ数、Rrdは想定する移動レンジ数rd の全範囲、prd,i(xi )は移動レンジ数rd を仮定したときのi番目のセルにおける目標検出確率密度関数、zrdは移動レンジ数rd を仮定したときに定まるMCPI目での目標位置、p(si,j )はパスsi,j を構造するセルの検波出力値の総和、imax はセルci に通じる全パス数を示している。
【0082】
移動量補償積分回路18における積分回路13と、警報回路10、追尾回路11は実施の形態1と同様に動作する。
【0083】
以上のように、この実施の形態4によれば、移動量補償積分回路18が検波回路8からのレンジビン毎の検波信号を入力し、所望の誤警報確率が得られるように、かつ、追尾回路11が求めた目標観測位置の確率分布に基づき、所望の目標検出確率が得られるように積分路数を決定して、初期レンジビンを起点とする目標が移動すると想定される積分路を設定し、設定された積分路におけるインコヒーレント積分処理を行い、設定された積分路を基に初期レンジビンに対応した固定スレッショルドを計算し、警報回路10がインコヒーレント積分処理結果と初期レンジビンに対応した固定スレッショルドを比較して目標信号を判定することにより、所望の誤警報確率及び所望の目標検出確率を達成できるという効果が得られる。
【0084】
実施の形態5.
図12はこの発明の実施の形態5による目標検出装置の構成を示すブロック図である。20は検波回路8からのレンジビン毎の検波信号を入力し、初期レンジビンを起点とする目標が移動すると想定される積分路を設定し、設定された積分路におけるインコヒーレント積分処理を行い、設定された積分路を基に所望の誤警報確率が得られるように、初期レンジビンに対応した固定スレッショルドを、各積分路の積分結果が無相関であると近似させることで簡易的に計算する移動量補償積分回路である。その他の構成は実施の形態1の図1に示す構成と同等である。
【0085】
図13は移動量補償積分回路20の内部構成を示すブロック図である。図において、21は検波回路8からのレンジビン毎の検波信号を入力し、初期レンジビンを起点とする目標が移動すると想定される積分路を設定し、設定された積分路を基に所望の誤警報確率が得られるように、初期レンジビンに対応した固定スレッショルドを、各積分路の積分結果が無相関であると近似させることで簡易的に計算する誤警報簡易計算型積分路設定回路であり、積分回路13は実施の形態1の図2に示すものと同等である。
【0086】
上記実施の形態1では、移動量補償積分回路9における積分路設定回路12が上記(17)式を用いて固定スレッショルドを計算しているが、このとき、実用的には演算処理負荷が大きく、実際に目標を追尾する過程での実時間処理には向いていない。そこで、この実施の形態5は、実時間処理可能な簡易的な式を用いて固定スレッショルドを計算するものである。
【0087】
次に動作について説明する。
変調回路1、送信機2、送信アンテナ3、受信アンテナ4、受信機5、A/D変換器6、復調回路7、検波回路8は、図17に示す従来例と同様に動作して、検波回路8からのレンジビン毎の検波信号が移動量補償積分回路20に出力される。
【0088】
移動量補償積分回路20における誤警報簡易計算型積分路設定回路21は、検波回路8からのレンジビン毎の検波信号を入力し、初期レンジビンを起点とする目標が移動すると想定される積分路を設定して積分回路13に出力すると共に、設定された積分路を基に所望の誤警報確率が得られるように、初期レンジビンに対応した固定スレッショルドを、次の(22)式により各積分路の積分結果が無相関であると近似させることで簡易的に計算して警報回路10に出力する。
【数18】
ここで、Pfa0 は一つの積分路当たり発生する誤警報確率、Naは積分路数を示している。
【0089】
ここで、上記(22)式における一つの積分路当り発生する誤警報確率Pfa0 は次の(23)式により計算する。
【数19】
ここで、Thは固定スレッショルド、A(x)*B(x)は関数A(x)と関数B(x)の畳み込み演算を示している。
【0090】
上記(22)式と(23)式を用いて、誤警報簡易計算型積分路設定回路21は、所望の誤警報確率が得られるように固定スレッショルドを簡易的に計算して警報回路10に出力する。
【0091】
移動量補償積分回路20における積分回路13と、警報回路10、追尾回路11は実施の形態1と同様に動作する。
【0092】
以上のように、この実施の形態5によれば、移動量補償積分回路20が検波回路8からのレンジビン毎の検波信号を入力し、初期レンジビンを起点とする目標が移動すると想定される積分路を設定し、設定された積分路におけるインコヒーレント積分処理を行い、設定された積分路を基に所望の誤警報確率が得られるように、初期レンジビンに対応した固定スレッショルドを、各積分路の積分結果が無相関であると近似させることで簡易的に計算し、警報回路10がインコヒーレント積分処理結果と初期レンジビンに対応した固定スレッショルドを比較して目標信号を判定することにより、インコヒーレント積分路数が増加しても、所望の誤警報確率を達成できるという効果が得られる。
【0093】
また、この実施の形態5によれば、移動量補償積分回路20が(22)式、(23)式を用いて固定スレッショルドを簡易的に計算することにより、固定スレッショルドの計算を高速に行い実時間処理ができるという効果が得られる。
【0094】
実施の形態6.
図14はこの発明の実施の形態6による目標検出装置の構成を示すブロック図である。図において、22は検波回路8からのレンジビン毎の検波信号を入力し、追尾回路11が上記(16)式で求めた目標観測位置zk の確率分布ps (zk )に基づきMCPI目において目標存在確率の低い領域を含む積分路(以下、目標存在確率の低い領域の積分路)を棄却し、棄却した積分路の演算負荷をMCPI目において目標存在確率の高い領域を含む積分路(以下、目標存在確率の高い領域の積分路)目標存在確率の高い領域の積分路における積分点数増加に転用することで、積分路数と積分点数との間で負荷を調整して、初期レンジビンを起点とする目標が移動すると想定される目標存在確率の高い領域の積分路と積分点数を設定し、設定された積分路におけるインコヒーレント積分処理を行い、設定された積分路を基に所望の誤警報確率が得られるように、初期レンジビンに対応した固定スレッショルドを計算する移動量補償積分回路である。その他の構成は実施の形態1の図1に示す構成と同等である。
【0095】
図15は移動量補償積分回路22の内部構成を示すブロック図である。図において、23は検波回路8からのレンジビン毎の検波信号を入力し、追尾回路11が上記(16)式で求めた目標観測位置zk の確率分布ps (zk )に基づき目標存在確率の低い領域の積分路を棄却し、棄却した積分路の演算負荷を目標存在確率の高い領域の積分路における積分点数増加に転用することで、積分路数と積分点数との間で負荷を調整して、初期レンジビンを起点とする目標が移動すると想定される目標存在確率の高い領域の積分路と積分点数を設定し、設定された積分路を基に所望の誤警報確率が得られるように、初期レンジビンに対応した固定スレッショルドを計算する負荷配分型積分路設定回路であり、積分回路13は実施の形態1の図2に示すものと同等である。
【0096】
次に動作について説明する。
変調回路1、送信機2、送信アンテナ3、受信アンテナ4、受信機5、A/D変換器6、復調回路7、検波回路8は、図17に示す従来例と同様に動作して、検波回路8からのレンジビン毎の検波信号が移動量補償積分回路22に出力される。
【0097】
移動量補償積分回路22における負荷配分型積分路設定回路23は、予め演算負荷か決められており、検波回路8からのレンジビン毎の検波信号を入力し、追尾回路11が上記(16)式で求めた目標観測位置zk の確率分布ps (zk )に基づき目標存在確率の低い領域の積分路を棄却し、棄却した積分路の演算負荷を目標存在確率の高い領域の積分路における積分点数増加に転用することで、積分路数と積分点数との間で負荷を調整して、初期レンジビンを起点とする目標が移動すると想定される目標存在確率の高い領域の積分路と積分点数を設定して積分回路13に出力すると共に、設定された積分路を基に所望の誤警報確率が得られるように、初期レンジビンに対応した固定スレッショルドを計算して警報回路10に出力する。
【0098】
図16は演算負荷の配分例を説明する図であり、積分路数を制限した際の余剰リソースを積分点数に配分した例を示している。図16では、積分点数が各々2である積分路Aと積分路Bを、積分点数が3である1種類の積分路Cに限定している。積分路Aの和演算回数は1で、積分路Bの和演算回数も1であり、積分路Aと積分路Bの和演算回数の合計は2である。また、積分路Cの和演算回数も2であり、積分路数を制限する前と積分路数を制限した後の和演算回数は変化していない。
【0099】
移動量補償積分回路22における積分回路13と、警報回路10、追尾回路11は実施の形態1と同様に動作する。
【0100】
以上のように、この実施の形態6によれば、負荷配分型積分路設定回路23が検波回路8からのレンジビン毎の検波信号を入力し、追尾回路11が求めた目標観測位置の確率分布に基づき目標存在確率の低い領域の積分路を棄却し、棄却した積分路の演算負荷を目標存在確率の高い領域の積分路における積分点数増加に転用することで、積分路数と積分点数との間で負荷を調整して、初期レンジビンを起点とする目標が移動すると想定される目標存在確率の高い領域の積分路と積分点数を設定し、設定された積分路におけるインコヒーレント積分処理を行い、設定された積分路を基に所望の誤警報確率が得られるように、初期レンジビンに対応した固定スレッショルドを計算し、警報回路10がインコヒーレント積分処理結果と初期レンジビンに対応した固定スレッショルドを比較して目標信号を判定することにより、インコヒーレント積分路数が増加しても、所望の誤警報確率を達成できるという効果が得られる。
【0101】
また、この実施の形態6によれば、負荷配分型積分路設定回路23が棄却した積分路の演算負荷を目標存在確率の高い領域の積分路における積分点数増加に転用することで、目標検出確率を改善することができるという効果が得られる。
【0102】
【発明の効果】
以上のように、この発明によれば、検波回路からのレンジビン毎の検波信号を入力し、初期レンジビンを起点とする目標が移動すると想定される積分路を設定し、設定された積分路におけるインコヒーレント積分処理を行い、設定された積分路を基に所望の誤警報確率が得られるように、雑音の上記検波回路出力時の確率密度関数、積分路設定を行った際のセルに通じるパスを構成するセルの検波出力値の総和、セルに通じる全パス数及び初期セルの番号により、初期レンジビンに対応した固定スレッショルドを計算する移動量補償積分回路と、移動量補償積分回路によるインコヒーレント積分処理結果と初期レンジビンに対応した固定スレッショルドを比較して目標信号を判定する警報回路とを備えたことことにより、インコヒーレント積分路数が増加しても、所望の誤警報確率を達成できるという効果がある。
【0103】
この発明によれば、移動量補償積分回路が初期レンジビンに対応した固定スレッショルドを、各積分路の積分結果が無相関であると近似させることで簡易的に計算することにより、インコヒーレント積分路数が増加しても、所望の誤警報確率を達成できるという効果がある。
【0104】
この発明によれば、警報回路により判定された目標信号の追尾処理を行い目標の位置を推定する追尾回路を備え、移動量補償積分回路が、追尾回路からの目標観測位置の確率分布に基づき、目標存在確率の高い領域を選択して積分路を設定することにより、インコヒーレント積分路数が増加しても、所望の誤警報確率を達成できるという効果がある。
【0105】
この発明によれば、移動量補償積分回路が、追尾回路からの目標観測位置の確率分布に基づき、目標存在確率の低い領域の積分路を棄却し、棄却した積分路の演算負荷を目標存在確率の高い領域の積分路における積分点数増加に転用することで、積分路数と積分点数との間の演算負荷を調整して、目標存在確率の高い領域の積分路と積分点数を設定することにより、インコヒーレント積分路数が増加しても、所望の誤警報確率を達成できるという効果がある。
【0106】
この発明によれば、検波回路からのレンジビン毎の検波信号を入力し、所望の誤警報確率が得られるようにレンジビン数とレンジビンの補償間隔により積分路数を決定して、初期レンジビンを起点とする目標が移動すると想定される積分路を設定し、設定された積分路におけるインコヒーレント積分処理を行い、設定された積分路を基に初期レンジビンに対応した固定スレッショルドを計算する移動量補償積分回路と、移動量補償積分回路によるインコヒーレント積分処理結果と初期レンジビンに対応した固定スレッショルドを比較して目標信号を判定する警報回路とを備えたことにより、インコヒーレント積分路数が増加しても、所望の誤警報確率を達成できるという効果がある。
【0107】
この発明によれば、警報回路により判定された目標信号の追尾処理を行い目標の位置を推定する追尾回路を備え、移動量補償積分回路が、所望の誤警報確率が得られるように、かつ、追尾回路からの目標観測位置の確率分布に基づき、所望の目標検出確率が得られるように積分路数を決定することにより、所望の誤警報確率及び所望の目標検出確率を達成できるという効果がある。
【図面の簡単な説明】
【図1】 この発明の実施の形態1による目標検出装置の構成を示すブロック図である。
【図2】 この発明の実施の形態1の移動量補償積分回路の内部構成を示すブロック図である。
【図3】 この発明の実施の形態1の移動量補償積分回路の入力信号の状況を説明する図である。
【図4】 この発明の実施の形態1の誤警報確率計算に使用するパスの例を説明する図である。
【図5】 この発明の実施の形態2による目標検出装置の構成を示すブロック図である。
【図6】 この発明の実施の形態2の移動量補償積分回路の内部構成を示すブロック図である。
【図7】 この発明の実施の形態2の目標存在範囲を説明する図である。
【図8】 この発明の実施の形態3による目標検出装置の構成を示すブロック図である。
【図9】 この発明の実施の形態3の移動量補償積分回路の内部構成を示すブロック図である。
【図10】 この発明の実施の形態4による目標検出装置の構成を示すブロック図である。
【図11】 この発明の実施の形態4の移動量補償積分回路の内部構成を示すブロック図である。
【図12】 この発明の実施の形態5による目標検出装置の構成を示すブロック図である。
【図13】 この発明の実施の形態5の移動量補償積分回路の内部構成を示すブロック図である。
【図14】 この発明の実施の形態6による目標検出装置の構成を示すブロック図である。
【図15】 この発明の実施の形態6の移動量補償積分回路の内部構成を示すブロック図である。
【図16】 この発明の実施の形態6による演算負荷の配分例を説明する図である。
【図17】 従来の目標検出装置の構成を示すブロック図である。
【図18】 従来の復調回路の内部構成を示すブロック図である。
【図19】 従来の目標検出装置の復調処理を説明するタイミングチャートである。
【図20】 従来の目標検出装置による移動量補償積分回路の入力信号の状況を説明する図である。
【図21】 従来の目標検出装置の移動量補償積分処理を説明するフローチャートである。
【図22】 従来の目標検出装置の構成を示すブロック図である。
【符号の説明】
1 変調回路、2 送信機、3 送信アンテナ、4 受信アンテナ、5 受信機、6 A/D変換器、7 復調回路、8 検波回路、9 移動量補償積分回路、10 警報回路、11 追尾回路、12 積分路設定回路、13 積分回路、14 移動量補償積分回路、15 選択型積分路設定回路、16 移動量補償積分回路、17 誤警報確率調整型積分路設定回路、18 移動量補償積分回路、19 誤警報確率/目標検出確率調整型積分路設定回路、20 移動量補償積分回路、21 誤警報簡易計算型積分路設定回路、22 移動量補償回路、23 負荷配分型積分路設定回路。[0001]
BACKGROUND OF THE INVENTION
The present invention relates to a target detection apparatus that performs target detection processing from a radar reception signal.
[0002]
[Prior art]
FIG. 17 is a block diagram showing a configuration of a conventional target detection apparatus. In the figure, 1 is a modulation circuit that modulates a transmission signal and outputs a modulation signal, 2 is a transmitter that forms a transmission beam by the modulation signal from the
[0003]
In FIG. 17, 8 is a detection circuit that detects a demodulated signal for each range bin from the
[0004]
FIG. 18 is a block diagram showing the internal configuration of the
[0005]
Next, the operation will be described.
The
s 0 (T) = s 0 (T + T) (1)
[0006]
The
[0007]
Assuming that the sampling period of the A /
[Expression 1]
In the above equation (2), N is the number of signal points used for demodulating the received signal in the
[0008]
The reference signal multiplication circuits 25-1 to 25-N configured for each range bin shown in FIG. 18 are input with the same modulation signal as the modulation signal transmitted from the
[0009]
FIG. 19 is a timing chart for explaining demodulation processing in the
[0010]
At this time, the reference signal s of the kth sampling of the j-range bin r, j (K) is expressed by the following equation (3).
[Expression 2]
Where s 0 (T) indicates a transmission signal at time t.
[0011]
The target signal present in the j range bin is demodulated by the reference signal of the j range bin. That is, the reference signal multiplication circuit 25-j performs complex multiplication of the digital signal from the A /
[Equation 3]
Where s 1 (T) is a received signal at time t, and s r, j (K) is the reference signal of the kth sampling of the j-range bin, s 2, j (K) shows the output signal of the reference signal multiplication circuit 25-j.
[0012]
Then, the demodulation filter 26-j performs coherent integration according to the following equation (5), and performs a j-range bin demodulated signal s. 3, j Is output.
[Expression 4]
[0013]
When the target exists in the j range bin, the demodulated signal s of the j range bin 3, j The target signal component is integrated. The demodulated signal for each range bin from the demodulating
[0014]
The movement amount
[0015]
FIG. 20 is a diagram for explaining the state of the input signal of the movement amount
[0016]
FIG. 21 is a flowchart for explaining the movement amount compensation integration processing. The processing procedure shown in FIG. 21 is directed to the detection of a target that is approaching the radar direction, that is, the range bin number decreases with the CPI. In step ST11, the CPI number M and the range bin number N are set. Also, the range bin compensation interval δr is set. In this case, since the target is approaching the radar direction, the compensation interval δr of the range bin is negative.
[0017]
In step ST12, an initial range bin rn that assumes the presence of the target signal at the 1st CPI is set to N. In step ST13, the initial value of the target range movement amount dr is set to 0, and the movement range ja with the range bin compensation interval δr as a unit is set to 1.
[0018]
In step ST14, the range r where the target signal component exists in iCPI using the following equation (6): 0 (I) is calculated.
[Equation 5]
Here, rn indicates an initial range bin that assumes the presence of a target signal at the 1st CPI, and dr indicates a target range bin moving amount that moves by MCPI. Range r 0 The range bin c (i) including (i) is selected. In this case, the calculated range r 0 For example, the range bin c (i) is selected by rounding off after the decimal point of the value of (i).
[0019]
In step ST15, an incoherent integration PDI (Post Detection Integration) value P (rn, ja) for a target moving by the dr range bin by MCPI is calculated using the following equation (7).
[Formula 6]
Here, pc (i) represents the detection output value of the c (i) range bin.
[0020]
In step ST16, the target range bin moving amount dr per MCPI is increased by δr. In step ST17, when the target range shift amount dr in MCPI is larger than rn, the process proceeds to the next step ST19. Otherwise, in step ST18, the movement range ja is incremented, and the processing from steps ST14 to ST17 is repeated.
[0021]
In step ST19, the target initial range bin rn is decreased by 1. In step ST20, if the target initial range bin rn <1, the process ends. If not, the processes from step ST13 to ST20 are repeated.
[0022]
Two-dimensional data of P (rn, ja) (1 ≦ rn ≦ N, 1 ≦ ja ≦ N / δr) is output from the movement amount
[0023]
The
[0024]
The above processing is performed for each azimuth angle and elevation angle, and the detection target in the
[0025]
The
x = RcosElcosAz
y = RcosElsinAz
z = RsinEl (8)
Here, R represents a target distance, El represents an elevation angle, and Az represents an azimuth angle.
[0026]
A target motion model that assumes constant velocity motion is represented by the following equation (9).
[Expression 7]
[0027]
The target signal observation model is expressed by the following equation (10). The observation noise is calculated by converting an observation error of an angle (azimuth angle, elevation angle) and an observation error of a distance into an observation error in an xyz coordinate system.
[Equation 8]
[0028]
The above equation (10) represents a model in which the target position is observed in each SPI. When a smooth value x hat k-1 of k-1 SPI and its error covariance matrix P hat k-1 are given, a predicted value x tilde k of k SPI is calculated using the following equation (11). Thus, the predicted position of the target at the kSPI-th is calculated.
[Equation 9]
[0029]
Further, an error covariance matrix P tilde k of x tilde k is calculated using the following equation (12). The P tilde k is a 6 × 6 matrix.
[Expression 10]
[0030]
Kalman gain K at the kSPI-th using the following equation (13) k Calculate
[Expression 11]
Here, Rv is the observation noise variance (E [vv * ]), T represents a transposed matrix.
[0031]
Using the following equation (14), a smooth value x of hat x at the kSPI-th is calculated.
[Expression 12]
[0032]
Then, an error covariance matrix P hat k of x hat k is calculated using the following equation (15).
[Formula 13]
Here, I indicates a unit matrix.
[0033]
Further, as k + 1 → k, the predicted value of the next SPI is calculated using equation (11). Thereafter, the target predicted position and its error variance in each SPI are calculated by the same procedure. When the actual target is moving according to the motion model represented by the equation (9) and the target signal is observed according to the observation model represented by the equation (10), the observation data up to k-1SPI is used, and the Kalman filter is used. As a result of processing, the target observation position z at the kSPIth k Is a probability density function p expressed by the following equation (16): s (Z k ) Normal distribution.
[Expression 14]
Here, T represents a transposed matrix.
[0034]
The
[0035]
FIG. 22 is a block diagram showing a configuration of a conventional target detection apparatus disclosed in Japanese Patent Laid-Open No. 10-68771. In the figure, 27 is a coherent integration circuit that performs coherent integration processing on an input signal and converts the signal component in the time domain into a signal component in the Doppler frequency domain, and 28 is incoherent in consideration of movement (acceleration) of the target Doppler frequency. It is an acceleration compensation integration circuit that performs integration. Other configurations are the same as those shown in FIG.
[0036]
In the signal processing system shown in FIG. 22, incoherent integration is performed while compensating for a target Doppler velocity (having acceleration) whose Doppler frequency changes. The demodulated signal from the
[0037]
The acceleration
[0038]
As a conventional target detection apparatus, there is a radar apparatus disclosed in Japanese Patent Laid-Open No. 8-179037. This is intended to improve target detection performance by performing incoherent processing in consideration of target range bin movement.
[0039]
[Problems to be solved by the invention]
Generally, when the number of incoherent integration paths in the movement amount
[0040]
Further, in the target detection devices disclosed in Japanese Patent Application Laid-Open Nos. 10-68771 and 8-179037, the setting of the fixed threshold considering the increase in the false alarm probability accompanying the increase in the number of integration paths is specifically described. Not mentioned.
[0041]
The present invention has been made to solve the above-described problems. Even when the number of incoherent integration paths increases, a fixed threshold corresponding to the initial range bin is set in the movement amount compensation integration circuit, and a desired false alarm is set. An object of the present invention is to obtain a target detection apparatus that can achieve the probability.
[0042]
[Means for Solving the Problems]
A target detection apparatus according to the present invention includes a demodulation circuit that demodulates a target signal generated by reflection on a target for each range bin, a detection circuit that detects a demodulated signal for each range bin from the demodulation circuit, and a detection circuit The detection signal for each range bin is input, the integration path where the target starting from the initial range bin is assumed to move is set, incoherent integration processing is performed on the set integration path, and the integration path is set based on the set integration path. To get the desired false alarm probability, The probability density function at the time of the detection circuit output of noise, the sum of the detection output values of the cells constituting the path leading to the cell when setting the integration path, the total number of paths leading to the cell, and the initial cell number, A moving amount compensation integrating circuit for calculating a fixed threshold corresponding to the initial range bin, an alarm circuit for comparing the incoherent integration processing result by the moving amount compensating integrating circuit and the fixed threshold corresponding to the initial range bin to determine a target signal; It is equipped with.
[0043]
In the target detection apparatus according to the present invention, the movement compensation compensation integration circuit simply calculates the fixed threshold corresponding to the initial range bin by approximating that the integration result of each integration path is uncorrelated.
[0044]
A target detection apparatus according to the present invention includes a tracking circuit that performs tracking processing of a target signal determined by an alarm circuit and estimates a target position, and the movement amount compensation integration circuit has a probability of a target observation position from the tracking circuit. Based on the distribution, an integration path is set by selecting a region having a high target existence probability.
[0045]
In the target detection device according to the present invention, the movement compensation compensation circuit rejects the integration path in the region where the target existence probability is low based on the probability distribution of the target observation position from the tracking circuit, and calculates the calculation load of the rejected integration path. By diverting it to increase the number of integration points in the integration path in the region with high target existence probability, adjust the calculation load between the number of integration paths and the number of integration points, and set the integration path and integration point in the region with high target existence probability To do.
[0046]
A target detection apparatus according to the present invention includes a demodulation circuit that demodulates a target signal generated by reflection on a target for each range bin, a detection circuit that detects a demodulated signal for each range bin from the demodulation circuit, and a detection circuit The detection signal for each range bin is input so that the desired false alarm probability can be obtained. Depending on the number of range bins and the range bin compensation interval Determine the number of integration paths, set the integration path where the target starting from the initial range bin is supposed to move, perform incoherent integration processing on the set integration path, and set the initial range bin based on the set integration path The movement amount compensation integration circuit that calculates a fixed threshold corresponding to the above, and an alarm circuit that compares the incoherent integration processing result by the movement amount compensation integration circuit with the fixed threshold corresponding to the initial range bin to determine the target signal. It is a thing.
[0047]
A target detection apparatus according to the present invention includes a tracking circuit that performs tracking processing of a target signal determined by an alarm circuit and estimates a target position, so that a movement amount compensation integration circuit can obtain a desired false alarm probability. In addition, based on the probability distribution of the target observation position from the tracking circuit, the number of integration paths is determined so as to obtain a desired target detection probability.
[0048]
DETAILED DESCRIPTION OF THE INVENTION
An embodiment of the present invention will be described below.
1 is a block diagram showing a configuration of a target detection apparatus according to
[0049]
FIG. 2 is a block diagram showing an internal configuration of the movement amount
[0050]
Next, the operation will be described.
The
[0051]
At this time, the integration
[Expression 15]
Where N c Is the initial cell number, p n (X) is a probability density function at the time of noise detection circuit output, and s i, j Is cell c i Jth path leading to, p (s i, j ) Is the path s i, j Sum of detection output values of cells constituting i, i max Is cell c i The total number of paths leading to is shown.
[0052]
3A and 3B are diagrams for explaining the state of the input signal in the movement compensation
[0053]
FIG. 4 is a diagram for explaining an example of a path used for false alarm probability calculation, and FIG. 1 Fig. 4 (b) shows the path leading to cell c. Four FIG. 4 (c) shows a path leading to cell c. 2 3, and an example of the path leading to the cell when the integration path setting of FIG. 3 is performed.
[0054]
FIG. 3A shows cells in the range of 3CPI and 3 range bins. When the initial cell at the 1st CPI is 1 range bin and the range bin movement compensation interval δr = 1, there are four types of integration paths set as shown in FIG. A set of all cells included in the four types of integration paths is as shown in FIG. In FIG. 3C, numbering is performed sequentially from the left cell of the 3CPI.
[0055]
As shown in FIG. 4 (a), cell c 1 The path leading to is {c 6 , C Four }. In addition, as shown in FIG. Four The path leading to is {c 6 }. Furthermore, as shown in FIG. 2 The path leading to is {c 6 , C Four } And {c 6 , C Five }.
[0056]
The path leading to each cell is examined as described above, and the fixed threshold Th that becomes the false alarm probability Pfa is calculated using the above equation (17). Actually, when all the paths are examined, the following equation (18) is obtained.
s 1,1 = {C 6 , C Four }
s 2,1 = {C 6 , C Four }
s 2,2 = {C 6 , C Five }
s 3,1 = {C 6 , C Five }
s 4,1 = {C 6 }
s 5,1 = {C 6 } (18)
[0057]
In addition, the maximum value when the integrated values related to the detection circuit output value of each path are compared is calculated as the following equation (19).
Where x i Is cell c i This shows the output value of the detector circuit.
[0058]
Substituting the above equation (19) into the above equation (17), the false alarm probability is calculated by the following equation (20).
[Expression 16]
Where p n (X 1 ) ~ P n (X 6 ) Is the noise component power x included in each cell. 1 ~ X 6 The probability density function of is shown.
[0059]
The integration
[0060]
As described above, according to the first embodiment, the movement amount
[0061]
FIG. 5 is a block diagram showing a configuration of a target detection apparatus according to
[0062]
FIG. 6 is a block diagram showing an internal configuration of the movement amount
[0063]
Next, the operation will be described.
The
[0064]
The selection type integration
[0065]
FIG. 7 is a diagram for explaining the target existence range, and shows the movement amount compensation range set in the selection type integration
[0066]
The
[0067]
As described above, according to the second embodiment, the movement amount
[0068]
Further, according to the second embodiment, the movement amount
[0069]
FIG. 8 is a block diagram showing a configuration of a target detection apparatus according to
[0070]
FIG. 9 is a block diagram showing an internal configuration of the movement amount
[0071]
Next, the operation will be described.
The
[0072]
The false alarm probability adjustment type integration
[0073]
In general, increasing the number of integration paths increases the false alarm probability, and decreasing the number of integration paths decreases the false alarm probability. Here, assuming that the number of range bins is N and the compensation interval of the range bin is δr, the number of integration paths is represented by N / δr. Therefore, the false alarm probability adjustment type integration
[0074]
The
[0075]
As described above, according to the third embodiment, the movement compensation
[0076]
10 is a block diagram showing a configuration of a target detection apparatus according to
[0077]
FIG. 11 is a block diagram showing an internal configuration of the movement amount
[0078]
Next, the operation will be described.
The
[0079]
The false alarm probability / target detection probability adjustment type integration
[0080]
In general, increasing the number of integration paths increases the target detection probability, and decreasing the number of integration paths decreases the target detection probability. Therefore, similarly to the false alarm probability, the target detection probability can be controlled by determining the number of integration paths based on the number of range bins N and the range bin compensation interval δr. Therefore, the false alarm probability / target detection probability adjustment type integration
[0081]
Here, the target detection probability Pd is calculated by the following equation (21).
[Expression 17]
Where r d Is the number of ranges moved during MCPI, R rd Is the assumed number of movement ranges r d The full range of p rd, i (X i ) Is the number of movement ranges r d Z is the target detection probability density function in the i-th cell, rd Is the number of moving ranges r d The target position at the MCPI, which is determined by assuming that p (s i, j ) Is the path s i, j The sum of the detection output values of the cells constituting max Is cell c i The total number of paths leading to is shown.
[0082]
The
[0083]
As described above, according to the fourth embodiment, the movement amount
[0084]
12 is a block diagram showing a configuration of a target detection apparatus according to
[0085]
FIG. 13 is a block diagram showing an internal configuration of the movement amount
[0086]
In the first embodiment, the integration
[0087]
Next, the operation will be described.
The
[0088]
The false alarm simple calculation type integration
[Expression 18]
Where P fa0 Is the probability of false alarms occurring per integration path, and Na is the number of integration paths.
[0089]
Here, the false alarm probability P generated per integration path in the above equation (22). fa0 Is calculated by the following equation (23).
[Equation 19]
Here, Th represents a fixed threshold, and A (x) * B (x) represents a convolution operation of the function A (x) and the function B (x).
[0090]
Using the above equations (22) and (23), the false alarm simple calculation type integration
[0091]
The
[0092]
As described above, according to the fifth embodiment, the movement amount
[0093]
Further, according to the fifth embodiment, the movement amount
[0094]
FIG. 14 is a block diagram showing a configuration of a target detection apparatus according to
[0095]
FIG. 15 is a block diagram showing an internal configuration of the movement amount
[0096]
Next, the operation will be described.
The
[0097]
The load distribution type integration
[0098]
FIG. 16 is a diagram for explaining an example of distribution of calculation loads, and shows an example in which surplus resources when the number of integration paths is limited are distributed to the number of integration points. In FIG. 16, the integration path A and the integration path B each having two integration points are limited to one type of integration path C having three integration points. The number of sum operations of the integration path A is 1, the number of sum operations of the integration path B is 1, and the sum of the number of sum operations of the integration path A and the integration path B is 2. In addition, the number of sum operations for the integration path C is 2, and the number of sum operations before the number of integration paths is limited and after the number of integration paths is limited are not changed.
[0099]
The
[0100]
As described above, according to the sixth embodiment, the load distribution type integration
[0101]
Further, according to the sixth embodiment, the target detection probability is obtained by diverting the calculation load of the integration path rejected by the load distribution type integration
[0102]
【The invention's effect】
As described above, according to the present invention, the detection signal for each range bin from the detection circuit is input, the integration path where the target starting from the initial range bin is assumed to move is set, and the input in the set integration path is set. Perform a coherent integration process, so that the desired false alarm probability is obtained based on the set integration path, The probability density function at the time of the detection circuit output of noise, the sum of the detection output values of the cells constituting the path leading to the cell when setting the integration path, the total number of paths leading to the cell, and the initial cell number, A displacement compensation integration circuit that calculates a fixed threshold corresponding to the initial range bin, and an alarm circuit that compares the incoherent integration processing result by the displacement compensation integration circuit with the fixed threshold corresponding to the initial range bin to determine the target signal. As a result, even if the number of incoherent integration paths increases, the desired false alarm probability can be achieved.
[0103]
According to the present invention, the number of incoherent integration paths is calculated by simply calculating the fixed threshold corresponding to the initial range bin by the movement amount compensation integration circuit by approximating that the integration result of each integration path is uncorrelated. Even if this increases, the desired false alarm probability can be achieved.
[0104]
According to this invention, the tracking circuit that performs the tracking process of the target signal determined by the alarm circuit and estimates the target position is provided, and the movement amount compensation integration circuit is based on the probability distribution of the target observation position from the tracking circuit, By selecting an area having a high target existence probability and setting integration paths, a desired false alarm probability can be achieved even if the number of incoherent integration paths increases.
[0105]
According to the present invention, the moving amount compensation integration circuit rejects the integration path in the region having a low target existence probability based on the probability distribution of the target observation position from the tracking circuit, and sets the calculation load of the rejected integration path to the target existence probability. By adjusting the calculation load between the number of integration paths and the number of integration points, and by setting the integration path and the number of integration points in the area where the target existence probability is high Even if the number of incoherent integration paths increases, the desired false alarm probability can be achieved.
[0106]
According to the present invention, a detection signal for each range bin from the detection circuit is input so that a desired false alarm probability can be obtained. Depending on the number of range bins and the range bin compensation interval Determine the number of integration paths, set the integration path where the target starting from the initial range bin is supposed to move, perform incoherent integration processing on the set integration path, and set the initial range bin based on the set integration path The movement amount compensation integration circuit that calculates the fixed threshold corresponding to, and the alarm circuit that compares the incoherent integration processing result by the movement amount compensation integration circuit and the fixed threshold corresponding to the initial range bin to determine the target signal As a result, even if the number of incoherent integration paths is increased, a desired false alarm probability can be achieved.
[0107]
According to the present invention, the tracking circuit that performs the tracking process of the target signal determined by the warning circuit and estimates the position of the target is provided, and the movement amount compensation integration circuit can obtain a desired false alarm probability, and There is an effect that a desired false alarm probability and a desired target detection probability can be achieved by determining the number of integration paths so that a desired target detection probability can be obtained based on the probability distribution of the target observation position from the tracking circuit. .
[Brief description of the drawings]
FIG. 1 is a block diagram showing a configuration of a target detection apparatus according to
FIG. 2 is a block diagram showing an internal configuration of a movement amount compensation integrating circuit according to the first embodiment of the present invention.
FIG. 3 is a diagram for explaining a state of an input signal of a movement amount compensation integrating circuit according to the first embodiment of the present invention.
FIG. 4 is a diagram illustrating an example of a path used for false alarm probability calculation according to the first embodiment of the present invention.
FIG. 5 is a block diagram showing a configuration of a target detection apparatus according to
FIG. 6 is a block diagram showing an internal configuration of a movement amount compensation integrating circuit according to
FIG. 7 is a diagram illustrating a target existence range according to the second embodiment of the present invention.
FIG. 8 is a block diagram showing a configuration of a target detection apparatus according to
FIG. 9 is a block diagram showing an internal configuration of a movement amount compensation integrating circuit according to
FIG. 10 is a block diagram showing a configuration of a target detection apparatus according to
FIG. 11 is a block diagram showing an internal configuration of a movement amount compensation integration circuit according to a fourth embodiment of the present invention.
FIG. 12 is a block diagram showing a configuration of a target detection apparatus according to
FIG. 13 is a block diagram showing an internal configuration of a movement amount compensation integration circuit according to a fifth embodiment of the present invention.
FIG. 14 is a block diagram showing a configuration of a target detection apparatus according to
FIG. 15 is a block diagram showing an internal configuration of a movement amount compensation integrating circuit according to a sixth embodiment of the present invention.
FIG. 16 is a diagram for explaining an example of calculation load distribution according to
FIG. 17 is a block diagram showing a configuration of a conventional target detection device.
FIG. 18 is a block diagram showing an internal configuration of a conventional demodulation circuit.
FIG. 19 is a timing chart illustrating a demodulation process of a conventional target detection device.
FIG. 20 is a diagram for explaining a state of an input signal of a movement amount compensation integration circuit by a conventional target detection device.
FIG. 21 is a flowchart for explaining movement amount compensation integration processing of a conventional target detection device;
FIG. 22 is a block diagram showing a configuration of a conventional target detection device.
[Explanation of symbols]
DESCRIPTION OF
Claims (6)
この復調回路からのレンジビン毎の復調信号を検波する検波回路と、
この検波回路からのレンジビン毎の検波信号を入力し、初期レンジビンを起点とする目標が移動すると想定される積分路を設定し、設定された積分路におけるインコヒーレント積分処理を行い、設定された積分路を基に所望の誤警報確率が得られるように、雑音の上記検波回路出力時の確率密度関数、積分路設定を行った際のセルに通じるパスを構成するセルの検波出力値の総和、セルに通じる全パス数及び初期セルの番号により、初期レンジビンに対応した固定スレッショルドを計算する移動量補償積分回路と、
この移動量補償積分回路によるインコヒーレント積分処理結果と上記初期レンジビンに対応した固定スレッショルドを比較して目標信号を判定する警報回路とを備えたことを特徴とする目標検出装置。A demodulation circuit that demodulates the target signal generated by reflection on the target for each range bin;
A detection circuit for detecting a demodulated signal for each range bin from the demodulation circuit;
Input the detection signal for each range bin from this detection circuit, set the integration path where the target starting from the initial range bin is assumed to move, perform incoherent integration processing in the set integration path, and set the integration The probability density function at the time of the above detection circuit output of noise, the sum of the detection output values of the cells constituting the path leading to the cell when setting the integration path, so that the desired false alarm probability can be obtained based on the path, A moving amount compensation integration circuit that calculates a fixed threshold corresponding to the initial range bin according to the total number of paths leading to the cell and the number of the initial cell, and
A target detection apparatus comprising: an alarm circuit that determines a target signal by comparing a result of incoherent integration processing by the movement amount compensation integration circuit and a fixed threshold corresponding to the initial range bin.
移動量補償積分回路が、上記追尾回路からの目標観測位置の確率分布に基づき、目標存在確率の高い領域を選択して積分路を設定することを特徴とする請求項1記載の目標検出装置。A tracking circuit that performs tracking processing of the target signal determined by the alarm circuit and estimates the position of the target is provided.
2. The target detection apparatus according to claim 1, wherein the movement compensation compensation integration circuit selects an area having a high target existence probability based on the probability distribution of the target observation position from the tracking circuit and sets the integration path.
この復調回路からのレンジビン毎の復調信号を検波する検波回路と、
この検波回路からのレンジビン毎の検波信号を入力し、所望の誤警報確率が得られるようにレンジビン数とレンジビンの補償間隔により積分路数を決定して、初期レンジビンを起点とする目標が移動すると想定される積分路を設定し、設定された積分路におけるインコヒーレント積分処理を行い、設定された積分路を基に初期レンジビンに対応した固定スレッショルドを計算する移動量補償積分回路と、
この移動量補償積分回路によるインコヒーレント積分処理結果と上記初期レンジビンに対応した固定スレッショルドを比較して目標信号を判定する警報回路とを備えたことを特徴とする目標検出装置。A demodulation circuit that demodulates the target signal generated by reflection on the target for each range bin;
A detection circuit for detecting a demodulated signal for each range bin from the demodulation circuit;
When the detection signal for each range bin from this detection circuit is input, the number of integration paths is determined by the number of range bins and the compensation interval of the range bin so that the desired false alarm probability can be obtained, and the target starting from the initial range bin moves A movement compensation compensation integration circuit that sets an assumed integration path, performs incoherent integration processing in the set integration path, and calculates a fixed threshold corresponding to the initial range bin based on the set integration path;
A target detection apparatus comprising: an alarm circuit that determines a target signal by comparing a result of incoherent integration processing by the movement amount compensation integration circuit and a fixed threshold corresponding to the initial range bin.
移動量補償積分回路が、所望の誤警報確率が得られるように、かつ、上記追尾回路からの目標観測位置の確率分布に基づき、所望の目標検出確率が得られるように積分路数を決定することを特徴とする請求項5記載の目標検出装置。A tracking circuit that performs tracking processing of the target signal determined by the alarm circuit and estimates the position of the target is provided.
The movement amount compensation integration circuit determines the number of integration paths so that a desired false alarm probability can be obtained and a desired target detection probability can be obtained based on the probability distribution of the target observation position from the tracking circuit. The target detection apparatus according to claim 5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001397919A JP4001745B2 (en) | 2001-12-27 | 2001-12-27 | Target detection device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001397919A JP4001745B2 (en) | 2001-12-27 | 2001-12-27 | Target detection device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2003194929A JP2003194929A (en) | 2003-07-09 |
| JP4001745B2 true JP4001745B2 (en) | 2007-10-31 |
Family
ID=27603548
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2001397919A Expired - Fee Related JP4001745B2 (en) | 2001-12-27 | 2001-12-27 | Target detection device |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP4001745B2 (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3695650B2 (en) * | 2001-03-14 | 2005-09-14 | 東陶機器株式会社 | Object detection device |
| JP5315585B2 (en) * | 2005-09-05 | 2013-10-16 | 日本電気株式会社 | Wake generation system, error covariance matrix initial value setting device, wake generation method, and error covariance matrix initial value setting method |
| JP4912762B2 (en) * | 2006-06-14 | 2012-04-11 | 三菱電機株式会社 | Target detection device |
| JP5036392B2 (en) * | 2007-04-25 | 2012-09-26 | 三菱電機株式会社 | Target detection device |
| JP5137693B2 (en) * | 2008-05-30 | 2013-02-06 | 三菱電機株式会社 | Image processing device |
-
2001
- 2001-12-27 JP JP2001397919A patent/JP4001745B2/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JP2003194929A (en) | 2003-07-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7196656B2 (en) | Apparatus for estimating direction of arrival of signal | |
| AU2002340719B2 (en) | System and method for narrowband pre-detection signal processing for passive coherent location applications | |
| KR20030018076A (en) | Target Prediction and Collision Warning System | |
| CN104380135B (en) | Phase Monopulse Radar System and Target Detection Method | |
| JP2597027B2 (en) | Wave observation radar | |
| JPWO2020246002A1 (en) | Radar signal processing equipment, radar system and signal processing method | |
| JP4001745B2 (en) | Target detection device | |
| JP2022116702A (en) | controller and program | |
| JP3608001B2 (en) | Passive radar device | |
| JP3580270B2 (en) | Target classification device and radar device having the same | |
| JP2005221374A (en) | GPS positioning device | |
| CN108650009B (en) | Wireless link occlusion estimation method for passive detection of aircraft | |
| JP2008026035A (en) | Radar | |
| JP2002131423A (en) | Radar equipment | |
| JP2009036539A (en) | Radar signal processing apparatus and radar signal processing method | |
| JP2005308428A (en) | Position estimation device | |
| KR100800851B1 (en) | Radar device | |
| JP3287784B2 (en) | Radar equipment | |
| JP3590598B2 (en) | Distance measuring device and measuring distance processing device | |
| JP7519817B2 (en) | Direction of arrival estimation device, direction of arrival estimation method, and direction of arrival estimation program | |
| JP2865079B2 (en) | Pulse radar receiver | |
| JP2011174866A (en) | Target detecting and tracking device | |
| JP2007033129A (en) | Distributed network radar equipment | |
| JPH0367594B2 (en) | ||
| JP2003337167A (en) | Radar wave reflection reduction device |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20041115 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20061110 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20061114 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20061222 |
|
| 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: 20070717 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20070815 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100824 Year of fee payment: 3 |
|
| 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: 20110824 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110824 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120824 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120824 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130824 Year of fee payment: 6 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |