JPH10206523A - 目標検出装置 - Google Patents
目標検出装置Info
- Publication number
- JPH10206523A JPH10206523A JP9011608A JP1160897A JPH10206523A JP H10206523 A JPH10206523 A JP H10206523A JP 9011608 A JP9011608 A JP 9011608A JP 1160897 A JP1160897 A JP 1160897A JP H10206523 A JPH10206523 A JP H10206523A
- Authority
- JP
- Japan
- Prior art keywords
- target
- signal
- ratio
- likelihood ratio
- received signal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Radar Systems Or Details Thereof (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
(57)【要約】
【課題】 受信信号列における目標の位置とそのＳＮ比
をユーザに容易に認識させること。 【解決手段】 音波等の送信信号に対する応答としての
受信信号列各々からは、領域設定器71，分散算出器72お
よび尤度比算出器60により、同一位置に定常的に存在し
ている目標からの反射信号が目標信号として検出された
上、その目標の位置が受信信号列画像上に表示されてい
るが、その際に、尤度比算出器60では、その目標の位置
に併せてその目標のＳＮ比が特定された上、受信信号列
画像上にその位置とともに表示されているものである。
をユーザに容易に認識させること。 【解決手段】 音波等の送信信号に対する応答としての
受信信号列各々からは、領域設定器71，分散算出器72お
よび尤度比算出器60により、同一位置に定常的に存在し
ている目標からの反射信号が目標信号として検出された
上、その目標の位置が受信信号列画像上に表示されてい
るが、その際に、尤度比算出器60では、その目標の位置
に併せてその目標のＳＮ比が特定された上、受信信号列
画像上にその位置とともに表示されているものである。
Description
【０００１】
【発明の属する技術分野】本発明は、ソーナー、レーダ
ー等、時間的に、または周波数的に振幅値が変動してい
る複数の受信信号列から、各受信信号列の同位置に定常
的に存在している物体を目標物体として検出するための
目標検出装置に関するものである。
ー等、時間的に、または周波数的に振幅値が変動してい
る複数の受信信号列から、各受信信号列の同位置に定常
的に存在している物体を目標物体として検出するための
目標検出装置に関するものである。
【０００２】
【従来の技術】一般に、ソーナー、レーダー等、一定の
領域に音波、または電波を複数回に亘って送信し、この
送信に対する応答としての反射音波、または反射電波と
しての受信信号列中から、各受信信号列中の同位置に定
常的に存在している物体を目標物体として検出するため
の目標検出装置がこれまでに知られている。一般に、受
信信号は、目標とする物体（以下、目標と称す）からの
反射信号（以下、目標信号と称す）の他に、目標以外
（海面、海底および海中の散乱体等）からの、各受信信
号毎にランダムに発生する反射信号、更には、雑音から
構成されたものとなっている。その目標検出装置では、
目標信号以外の反射信号や雑音が効果的に抑圧された状
態として、目標信号が状態良好に検出された上、可視表
示されているわけであるが、このような目標信号の検出
・表示技術に関しては、重要な技術として従来より研究
が進められているのが実情である。例えば「科学計測の
ための波形データ処理」（南茂夫編著、１９９５年８月
１日 ＣＱ出版社発行：pp105〜106）に記載された積算
平均化法が知られている。積算平均化法とは、同一条件
の下で、繰り返し測定される多数の波形を対象として、
各波形の同一測定位置に対する集合平均化（単純平均
化）を行う手法であり、これを、簡単ながら説明すれば
以下のようである。
領域に音波、または電波を複数回に亘って送信し、この
送信に対する応答としての反射音波、または反射電波と
しての受信信号列中から、各受信信号列中の同位置に定
常的に存在している物体を目標物体として検出するため
の目標検出装置がこれまでに知られている。一般に、受
信信号は、目標とする物体（以下、目標と称す）からの
反射信号（以下、目標信号と称す）の他に、目標以外
（海面、海底および海中の散乱体等）からの、各受信信
号毎にランダムに発生する反射信号、更には、雑音から
構成されたものとなっている。その目標検出装置では、
目標信号以外の反射信号や雑音が効果的に抑圧された状
態として、目標信号が状態良好に検出された上、可視表
示されているわけであるが、このような目標信号の検出
・表示技術に関しては、重要な技術として従来より研究
が進められているのが実情である。例えば「科学計測の
ための波形データ処理」（南茂夫編著、１９９５年８月
１日 ＣＱ出版社発行：pp105〜106）に記載された積算
平均化法が知られている。積算平均化法とは、同一条件
の下で、繰り返し測定される多数の波形を対象として、
各波形の同一測定位置に対する集合平均化（単純平均
化）を行う手法であり、これを、簡単ながら説明すれば
以下のようである。
【０００３】即ち、図２６に積算平均化法の原理を示
す。各グラフの横軸は時間（あるいは音波、または電波
の単位当たりの伝搬距離を乗じることにより求まる距
離）、縦軸は振幅値の大きさを表す。図示のように、１
回目の測定波形Ｘ1(ｔ)、２回目の測定波形Ｘ2(ｔ)、…
… …、Ｋ回目の測定波形ＸK(ｔ)、積算平均化結果
ｚ(ｔ)がそれぞれ示されているが、それら測定波形Ｘ
1(ｔ)〜ＸK(ｔ)において、目標信号は時刻ｔ2 に存在す
る場合が想定されたものとなっている。１回目の測定波
形Ｘ1(ｔ)には、時刻ｔ1 に目標信号以外の信号が偶発
的に混入されているが、他の測定波形Ｘ2(ｔ)〜ＸK(ｔ)
各々には、その時刻ｔ1 に目標信号以外の信号が何等混
入されていないことから、測定波形Ｘ1(ｔ)〜ＸK(ｔ)に
対する積算平均化により目標信号以外の信号は、その
分、抑圧される結果として、積算平均化結果ｚ(ｔ)で
は、時刻ｔ2 における目標信号のみが顕在化された状態
の波形として検出され得るものであり、また、その際
に、積算平均後の振幅値からは、時刻ｔ2 における目標
信号についての平均的なＳＮ比（＝ＳＮＲ）が併せて計
算可とされているものである。
す。各グラフの横軸は時間（あるいは音波、または電波
の単位当たりの伝搬距離を乗じることにより求まる距
離）、縦軸は振幅値の大きさを表す。図示のように、１
回目の測定波形Ｘ1(ｔ)、２回目の測定波形Ｘ2(ｔ)、…
… …、Ｋ回目の測定波形ＸK(ｔ)、積算平均化結果
ｚ(ｔ)がそれぞれ示されているが、それら測定波形Ｘ
1(ｔ)〜ＸK(ｔ)において、目標信号は時刻ｔ2 に存在す
る場合が想定されたものとなっている。１回目の測定波
形Ｘ1(ｔ)には、時刻ｔ1 に目標信号以外の信号が偶発
的に混入されているが、他の測定波形Ｘ2(ｔ)〜ＸK(ｔ)
各々には、その時刻ｔ1 に目標信号以外の信号が何等混
入されていないことから、測定波形Ｘ1(ｔ)〜ＸK(ｔ)に
対する積算平均化により目標信号以外の信号は、その
分、抑圧される結果として、積算平均化結果ｚ(ｔ)で
は、時刻ｔ2 における目標信号のみが顕在化された状態
の波形として検出され得るものであり、また、その際
に、積算平均後の振幅値からは、時刻ｔ2 における目標
信号についての平均的なＳＮ比（＝ＳＮＲ）が併せて計
算可とされているものである。
【０００４】
【発明が解決しようとする課題】ところで、積算平均化
法による場合、目標信号についての平均的なＳＮ比が併
せて計算可とされているも、その位置に目標信号が定常
的に存在しているか否かについての判定は困難となって
いる。以下の表１に受信信号列のある２点Ａ，Ｂ各々に
ついての、１０回分に亘る振幅値の積算平均結果（平均
ＳＮＲ：単位はdB）を示す。
法による場合、目標信号についての平均的なＳＮ比が併
せて計算可とされているも、その位置に目標信号が定常
的に存在しているか否かについての判定は困難となって
いる。以下の表１に受信信号列のある２点Ａ，Ｂ各々に
ついての、１０回分に亘る振幅値の積算平均結果（平均
ＳＮＲ：単位はdB）を示す。
【０００５】
【表１】
【０００６】その表１から判るように、両地点Ａ，Ｂと
もに積算平均結果として、平均ＳＮ比が６dBであるとの
結果が得られているが、地点Ａについて、目標信号が存
在しているのは１回目の受信信号列のみであり、従っ
て、地点Ａに定常的に目標物体が存在しているとはいえ
ないものとなっている。
もに積算平均結果として、平均ＳＮ比が６dBであるとの
結果が得られているが、地点Ａについて、目標信号が存
在しているのは１回目の受信信号列のみであり、従っ
て、地点Ａに定常的に目標物体が存在しているとはいえ
ないものとなっている。
【０００７】よって、本発明の目的は、時間的に、また
は周波数的に振幅値が変動している複数の受信信号列か
ら、各受信信号列中の同一位置に定常的に存在している
目標を検出するに際して、その目標の位置に併せてその
ＳＮ比が容易に認識可とされた目標検出装置を供するに
ある。
は周波数的に振幅値が変動している複数の受信信号列か
ら、各受信信号列中の同一位置に定常的に存在している
目標を検出するに際して、その目標の位置に併せてその
ＳＮ比が容易に認識可とされた目標検出装置を供するに
ある。
【０００８】
【課題を解決するための手段】上記目的は、目標検出装
置が、ソーナー、レーダー等、音波、または電波が複数
回に亘って送信信号として送信される度に、その送信信
号各々に対する応答としての反射信号を受信信号列とし
て受信した上、その受信信号列各々は映像信号に変換さ
れた状態で画像として表示されつつ、上記受信信号列各
々からは、同一位置に定常的に存在している物体からの
反射信号が目標信号として検出された上、その目標の位
置が上記画像上に併せて表示されるべく構成されている
場合に、目標の位置とその目標のＳＮ比が画像上に同時
に表示されるべく、受信信号列各々からその目標の位置
とその目標のＳＮ比を特定する処理を行う目標位置・Ｓ
Ｎ比特定手段を少なくとも具備せしめることで達成され
る。
置が、ソーナー、レーダー等、音波、または電波が複数
回に亘って送信信号として送信される度に、その送信信
号各々に対する応答としての反射信号を受信信号列とし
て受信した上、その受信信号列各々は映像信号に変換さ
れた状態で画像として表示されつつ、上記受信信号列各
々からは、同一位置に定常的に存在している物体からの
反射信号が目標信号として検出された上、その目標の位
置が上記画像上に併せて表示されるべく構成されている
場合に、目標の位置とその目標のＳＮ比が画像上に同時
に表示されるべく、受信信号列各々からその目標の位置
とその目標のＳＮ比を特定する処理を行う目標位置・Ｓ
Ｎ比特定手段を少なくとも具備せしめることで達成され
る。
【０００９】
【発明の実施の形態】以下、本発明の実施形態を図１か
ら図２５により説明する。先ず本発明の具体的説明に先
立って、その理論的な背景について説明しておく。既
に、これまでに、目標信号を非目標信号や雑音から効果
的に弁別した上、その目標信号の検出性能をより向上さ
せる方法としては、例えば、モラード―バーケット著に
よる“シグナル―デティクション アンド エスティメ
ーション 頁１１５〜１７４”（Mourad-Barket；SIGNA
L-Detection & Estimation pp115〜174）に記載された
もの、即ち、尤度比（likelihood ratio）を用いた信号
検出法が知られている。この信号検出法による場合、各
地点の振幅値および周辺区間の分散値から尤度比が算出
された上、予め設定された閾値との比較により目標信号
が検出可とされているが、その原理は以下のようであ
る。
ら図２５により説明する。先ず本発明の具体的説明に先
立って、その理論的な背景について説明しておく。既
に、これまでに、目標信号を非目標信号や雑音から効果
的に弁別した上、その目標信号の検出性能をより向上さ
せる方法としては、例えば、モラード―バーケット著に
よる“シグナル―デティクション アンド エスティメ
ーション 頁１１５〜１７４”（Mourad-Barket；SIGNA
L-Detection & Estimation pp115〜174）に記載された
もの、即ち、尤度比（likelihood ratio）を用いた信号
検出法が知られている。この信号検出法による場合、各
地点の振幅値および周辺区間の分散値から尤度比が算出
された上、予め設定された閾値との比較により目標信号
が検出可とされているが、その原理は以下のようであ
る。
【００１０】即ち、その信号検出法は統計的検定法の１
つであり、仮説Ｈ0 ：母集団分布の確率密度関数はｆ0
(ｘ)である、を対立仮説Ｈ1 ：母集団分布の確率密度関
数はｆ1 (ｘ)である、に対して検定する方式であり、そ
の手順は観測値ｘ1，ｘ2，…を独立に１つずつ取り、観
測値を取る度に、次の３通りの決定のうち、何れか１通
りの決定を行うというものである。具体的には、ｍ番目
の観測値（ｍ=１,２,…）に対し、以下の数式１で定義
されている尤度比Λ(x)が求められる。
つであり、仮説Ｈ0 ：母集団分布の確率密度関数はｆ0
(ｘ)である、を対立仮説Ｈ1 ：母集団分布の確率密度関
数はｆ1 (ｘ)である、に対して検定する方式であり、そ
の手順は観測値ｘ1，ｘ2，…を独立に１つずつ取り、観
測値を取る度に、次の３通りの決定のうち、何れか１通
りの決定を行うというものである。具体的には、ｍ番目
の観測値（ｍ=１,２,…）に対し、以下の数式１で定義
されている尤度比Λ(x)が求められる。
【００１１】
【数１】
【００１２】条件１：もしも、Λ(x)が以下の数式２を
満たすならば、(ｍ+１)番目の観測を行う。
満たすならば、(ｍ+１)番目の観測を行う。
【００１３】
【数２】
【００１４】条件２：もしも、Λ(x)が以下の数式３を
満たすならば、観測値を取るのを止めＨ0 を採択する。
満たすならば、観測値を取るのを止めＨ0 を採択する。
【００１５】
【数３】
【００１６】条件３：もしも、Λ(x)が以下の数式４を
満たすならば、観測を止めＨ1 を採択する。
満たすならば、観測を止めＨ1 を採択する。
【００１７】
【数４】
【００１８】但し、それら数式２〜４中における閾値η
1 ，η0 は第１種の誤り（Ｈ0 が正しいのにＨ0 を棄却
する）、第２種の誤り（Ｈ1 が正しいのにＨ1 を棄却す
る）の確率がそれぞれα、βになるように選ばれる。尤
も、α、βが小さい値であれば、閾値η1 ，η0 はそれ
ぞれ近似的に以下の数式５として表現され得る。
1 ，η0 は第１種の誤り（Ｈ0 が正しいのにＨ0 を棄却
する）、第２種の誤り（Ｈ1 が正しいのにＨ1 を棄却す
る）の確率がそれぞれα、βになるように選ばれる。尤
も、α、βが小さい値であれば、閾値η1 ，η0 はそれ
ぞれ近似的に以下の数式５として表現され得る。
【００１９】
【数５】
【００２０】さて、以上の信号検出法を目標信号の検出
に適用する場合、 ｆ0(x)：目標信号を含まない受信信号列の分布が従うべ
き確率密度関数 ｆ1(x)：目標信号を含む受信信号列の分布が従うべき確
率密度関数 観測値：複数の受信信号列中の同一時刻（あるいは位
置）における振幅値ｘ1,ｘ2,… として設定した上、以下の判定を行う。
に適用する場合、 ｆ0(x)：目標信号を含まない受信信号列の分布が従うべ
き確率密度関数 ｆ1(x)：目標信号を含む受信信号列の分布が従うべき確
率密度関数 観測値：複数の受信信号列中の同一時刻（あるいは位
置）における振幅値ｘ1,ｘ2,… として設定した上、以下の判定を行う。
【００２１】即ち、上記条件１を満たすならば、(ｍ+
１)番目の受信信号列の観測を行い、また、もしも、上
記条件２を満たすならば、振幅値ｘ1,ｘ2,…の示す時刻
（あるいは位置）に目標が存在しないと判定し、更に、
上記条件３を満たすならば、振幅値ｘ1,ｘ2,…の示す時
刻（あるいは位置）に目標が存在すると判定する。
１)番目の受信信号列の観測を行い、また、もしも、上
記条件２を満たすならば、振幅値ｘ1,ｘ2,…の示す時刻
（あるいは位置）に目標が存在しないと判定し、更に、
上記条件３を満たすならば、振幅値ｘ1,ｘ2,…の示す時
刻（あるいは位置）に目標が存在すると判定する。
【００２２】以下、ソーナー、レーダー等、送受波器が
受信する受信信号の狭帯域フィルタ出力がレイリー分布
に従う場合を例に採り、数式１で定義されている尤度比
Λ(x)の具体的な設定方法を説明すれば、先ずレイリー
分布は以下の数式６で表現される。
受信する受信信号の狭帯域フィルタ出力がレイリー分布
に従う場合を例に採り、数式１で定義されている尤度比
Λ(x)の具体的な設定方法を説明すれば、先ずレイリー
分布は以下の数式６で表現される。
【００２３】
【数６】
【００２４】ところで、そのレイリー分布に従う受信信
号列に対し、振幅値Ａの正弦波信号が重畳された場合、
出力ｘは以下の数式７として示すライス分布に従うこと
が知られている。
号列に対し、振幅値Ａの正弦波信号が重畳された場合、
出力ｘは以下の数式７として示すライス分布に従うこと
が知られている。
【００２５】
【数７】
【００２６】従って、目標信号を含まない受信信号列の
分布が従うべき確率密度関数としてレイリー分布を、ま
た、目標信号を含む受信信号列の分布が従うべき確率密
度関数としてライス分布を代入することにより、尤度比
Λ(x)は以下の数式８として具体的に定義され得るもの
である。
分布が従うべき確率密度関数としてレイリー分布を、ま
た、目標信号を含む受信信号列の分布が従うべき確率密
度関数としてライス分布を代入することにより、尤度比
Λ(x)は以下の数式８として具体的に定義され得るもの
である。
【００２７】
【数８】
【００２８】その数式８より、尤度比Λ(x)の値は３つ
のパラメータ、即ち、振幅値ｘ、分散値σ2 、および受
信信号列に重畳される信号振幅値Ａから決定される。こ
のうち、信号振幅値Ａは、重畳される信号のＳＮ比から
求まる定数である。従って、受信信号が数式６として示
すレイリー分布に従う場合、受信信号列中の任意の領域
内の受信信号振幅値ｘおよびその周辺区間に含まれる受
信信号振幅値の分散値σ2 から尤度比Λ(x)が求まり、
数式５により予め定められた閾値η1 ，η0 との比較に
より、その領域に目標信号が存在するか否かの判定が可
能となる。以下、検出対象としての受信信号振幅値が含
まれている領域を目標捜索領域、分散値計算対象として
の受信信号振幅値が含まれている領域を分散計算領域と
記す。
のパラメータ、即ち、振幅値ｘ、分散値σ2 、および受
信信号列に重畳される信号振幅値Ａから決定される。こ
のうち、信号振幅値Ａは、重畳される信号のＳＮ比から
求まる定数である。従って、受信信号が数式６として示
すレイリー分布に従う場合、受信信号列中の任意の領域
内の受信信号振幅値ｘおよびその周辺区間に含まれる受
信信号振幅値の分散値σ2 から尤度比Λ(x)が求まり、
数式５により予め定められた閾値η1 ，η0 との比較に
より、その領域に目標信号が存在するか否かの判定が可
能となる。以下、検出対象としての受信信号振幅値が含
まれている領域を目標捜索領域、分散値計算対象として
の受信信号振幅値が含まれている領域を分散計算領域と
記す。
【００２９】さて、図２には、レイリー分布形状、ライ
ス分布形状およびこれら両者の比としての尤度比Λ(ｘ)
の変化が、横軸は振幅値、左縦軸は確率密度、右縦軸は
尤度比として示されているが、本例では、ライス分布に
おける信号振幅値Ａとして、Ａ＝２.４５(×σ：以下省
略)が想定されており、また、図面の簡単化上、尤度比
は対数尤度比lnΛ(ｘ)として示されたものとなっている
（以下、「尤度比Λ(ｘ)」は対数尤度比lnΛ(ｘ)を意味
するものとする）。
ス分布形状およびこれら両者の比としての尤度比Λ(ｘ)
の変化が、横軸は振幅値、左縦軸は確率密度、右縦軸は
尤度比として示されているが、本例では、ライス分布に
おける信号振幅値Ａとして、Ａ＝２.４５(×σ：以下省
略)が想定されており、また、図面の簡単化上、尤度比
は対数尤度比lnΛ(ｘ)として示されたものとなっている
（以下、「尤度比Λ(ｘ)」は対数尤度比lnΛ(ｘ)を意味
するものとする）。
【００３０】図３にはまた、尤度比Λ(ｘ)とＳＮ比の関
係が、横軸はＳＮ比(dB)、縦軸は尤度比Λ(ｘ)として示
されているが、これより、尤度比Λ(ｘ)はＳＮ比に対し
て単調増加傾向にあり、しかも、本例での場合（Ａ＝
２.４５）、約４dB以上のＳＮ比に対して正の値をとる
ことが判る。例えば閾値η1 ，η0 の算出式（数式５）
にα＝10-3(％)、β＝10(％)を代入した場合、閾値η
1 ，η0 は以下の数式９として得られる。
係が、横軸はＳＮ比(dB)、縦軸は尤度比Λ(ｘ)として示
されているが、これより、尤度比Λ(ｘ)はＳＮ比に対し
て単調増加傾向にあり、しかも、本例での場合（Ａ＝
２.４５）、約４dB以上のＳＮ比に対して正の値をとる
ことが判る。例えば閾値η1 ，η0 の算出式（数式５）
にα＝10-3(％)、β＝10(％)を代入した場合、閾値η
1 ，η0 は以下の数式９として得られる。
【００３１】
【数９】
【００３２】従って、定常的に４dB以上の信号が存在す
る領域では、複数の受信信号列に亘って大きい値の尤度
比Λ(ｘ)が算出され、複数受信信号の累積の結果、閾値
η1より大きくなった時点で、その領域に目標が存在し
ていると判定し得るものである。これとは逆に、図２６
における時刻ｔ1 のように、ある受信信号列のみに信号
（非目標信号）が出現している領域については、その受
信信号では一時的に大きい値の尤度比Λ(ｘ)が算出され
るが、他の受信信号列では信号が存在しないために尤度
比Λ(ｘ)の値が小さくなる結果として、その領域には目
標が存在しないと判定し得るものである。換言すれば、
累積尤度比Λ(ｘ)（以下、「累積尤度比Λ(ｘ)」は累加
対数尤度比ΣlnΛ(ｘ)を意味する）が閾値η0 より小さ
くなった時点で、その領域には目標が存在しないと判定
し得るものである。
る領域では、複数の受信信号列に亘って大きい値の尤度
比Λ(ｘ)が算出され、複数受信信号の累積の結果、閾値
η1より大きくなった時点で、その領域に目標が存在し
ていると判定し得るものである。これとは逆に、図２６
における時刻ｔ1 のように、ある受信信号列のみに信号
（非目標信号）が出現している領域については、その受
信信号では一時的に大きい値の尤度比Λ(ｘ)が算出され
るが、他の受信信号列では信号が存在しないために尤度
比Λ(ｘ)の値が小さくなる結果として、その領域には目
標が存在しないと判定し得るものである。換言すれば、
累積尤度比Λ(ｘ)（以下、「累積尤度比Λ(ｘ)」は累加
対数尤度比ΣlnΛ(ｘ)を意味する）が閾値η0 より小さ
くなった時点で、その領域には目標が存在しないと判定
し得るものである。
【００３３】以上の判定論理により、複数の受信信号列
に亘り定常的に存在する目標の検出が可能となるもので
ある。ただ、尤度比の原理をそのまま目標検出に適用し
ただけでは、その目標検出と同時に、その目標について
の平均ＳＮ比算出は困難である。以下の表２には、受信
信号列のある３点について各受信信号列におけるＳＮ比
と累積尤度比Λ(ｘ)の変化とともに、各ＳＮＲに対する
尤度比の算出値が参考として併せて示されているが、目
標物体が存在すると判定された領域について、平均的な
ＳＮ比を直接算出することは不可となっている。
に亘り定常的に存在する目標の検出が可能となるもので
ある。ただ、尤度比の原理をそのまま目標検出に適用し
ただけでは、その目標検出と同時に、その目標について
の平均ＳＮ比算出は困難である。以下の表２には、受信
信号列のある３点について各受信信号列におけるＳＮ比
と累積尤度比Λ(ｘ)の変化とともに、各ＳＮＲに対する
尤度比の算出値が参考として併せて示されているが、目
標物体が存在すると判定された領域について、平均的な
ＳＮ比を直接算出することは不可となっている。
【００３４】
【表２】
【００３５】なお、用いた尤度比算出式は数式８におい
て、Ａ＝２.４５を代入したものであり、尤度比Λ(ｘ)
は図３の対ＳＮ比特性を示す。
て、Ａ＝２.４５を代入したものであり、尤度比Λ(ｘ)
は図３の対ＳＮ比特性を示す。
【００３６】その表２より判るように、各地点Ｃ，Ｄ，
Ｅは定常的に５dB以上の信号が存在する領域であり、10
回以内の受信信号列の累積でその領域に目標が存在する
と判定されている。しかしながら、本検出法では、定常
的にある一定値以上（本例では４dB以上）のＳＮ比の信
号が存在する領域は全て目標が存在すると判定されてい
るため、目標が存在すると判定された領域について、平
均的なＳＮ比を直接算出することは不可となっている。
Ｅは定常的に５dB以上の信号が存在する領域であり、10
回以内の受信信号列の累積でその領域に目標が存在する
と判定されている。しかしながら、本検出法では、定常
的にある一定値以上（本例では４dB以上）のＳＮ比の信
号が存在する領域は全て目標が存在すると判定されてい
るため、目標が存在すると判定された領域について、平
均的なＳＮ比を直接算出することは不可となっている。
【００３７】そこで、以上の事情を踏まえ、本願発明者
等は尤度比の原理を拡張し、上記尤度比を用いた目標検
出処理に必要とされる手段に加え、目標信号を含む受信
信号列の分布が従うべきライス分布の確率密度関数を目
標信号のＳＮ比に応じて複数設定することにより、複数
の尤度比算出式を定義する手段と、検出された受信信号
振幅値および算出された受信信号振幅値の分散値を各々
の尤度比算出式に代入し、目標捜索領域に対する尤度比
を算出する手段と、算出された尤度比を予め定めた閾値
と比較し、閾値を上回る尤度比が算出された目標捜索領
域に対して目標が存在すると判定する手段と、目標が存
在すると判定された領域について、用いた尤度比算出式
毎に分類する手段とを追加することにより、各受信信号
列中の同一位置に定常的に存在している目標を検出する
と同時に、その目標のＳＮ比を算出する機能を有する目
標検出装置を考案するに至ったものである。ここで、目
標が存在すると判定された領域の平均的なＳＮ比を算出
する方法として、信号振幅値Ａが異なるライス分布を用
いた複数の尤度比算出式を設定する理由について述べれ
ば以下のようである。
等は尤度比の原理を拡張し、上記尤度比を用いた目標検
出処理に必要とされる手段に加え、目標信号を含む受信
信号列の分布が従うべきライス分布の確率密度関数を目
標信号のＳＮ比に応じて複数設定することにより、複数
の尤度比算出式を定義する手段と、検出された受信信号
振幅値および算出された受信信号振幅値の分散値を各々
の尤度比算出式に代入し、目標捜索領域に対する尤度比
を算出する手段と、算出された尤度比を予め定めた閾値
と比較し、閾値を上回る尤度比が算出された目標捜索領
域に対して目標が存在すると判定する手段と、目標が存
在すると判定された領域について、用いた尤度比算出式
毎に分類する手段とを追加することにより、各受信信号
列中の同一位置に定常的に存在している目標を検出する
と同時に、その目標のＳＮ比を算出する機能を有する目
標検出装置を考案するに至ったものである。ここで、目
標が存在すると判定された領域の平均的なＳＮ比を算出
する方法として、信号振幅値Ａが異なるライス分布を用
いた複数の尤度比算出式を設定する理由について述べれ
ば以下のようである。
【００３８】即ち、一般に、レイリー分布に従う受信信
号列に信号振幅値Ａの信号が重畳された場合、その際で
のＳＮ比は以下の数式１０として定義される。
号列に信号振幅値Ａの信号が重畳された場合、その際で
のＳＮ比は以下の数式１０として定義される。
【００３９】
【数１０】
【００４０】例えば、図２では、信号振幅値Ａとして、
Ａ＝２.４５が代入されているが、これは、数式１０に
よりＳＮ比＝６dBの信号が重畳された場合に相当する。
図２における横軸がＳＮ比に変換されたグラフを図４に
示す。そのグラフより６dB付近の確率密度が最も大きく
なっていることが判る。また、図３に示す尤度比とＳＮ
比の関係より、Ａ＝２.４５の場合は、ＳＮ比＝４dB以
上の信号に対して正の値の尤度比Λ(ｘ)が算出されるの
が判る。尤度比を用いた目標検出法では、各受信信号に
おいて算出される尤度比Λ(ｘ)の累積結果と閾値との比
較により目標が検出されており、この場合は、定常的に
４dB以上の信号が存在する領域では常に正の値の尤度比
Λ(ｘ)が算出されるので、複数受信信号列の累積により
目標の検出が可能となる。ここで、図５に、信号振幅値
Ａの値が異なる場合でのＳＮ比と確率密度関数および尤
度比Λ(ｘ)の関係を示す。図中では、Ａ＝４.２４（数
式１０よりＳＮ比10dBに相当）を用いている。Ａ＝４.
２４の場合では、確率密度関数のピークはＳＮ比10dB付
近に位置し、また、尤度比Λ(ｘ)はＳＮ比＝６dB以上の
信号に対して正の値が算出される。
Ａ＝２.４５が代入されているが、これは、数式１０に
よりＳＮ比＝６dBの信号が重畳された場合に相当する。
図２における横軸がＳＮ比に変換されたグラフを図４に
示す。そのグラフより６dB付近の確率密度が最も大きく
なっていることが判る。また、図３に示す尤度比とＳＮ
比の関係より、Ａ＝２.４５の場合は、ＳＮ比＝４dB以
上の信号に対して正の値の尤度比Λ(ｘ)が算出されるの
が判る。尤度比を用いた目標検出法では、各受信信号に
おいて算出される尤度比Λ(ｘ)の累積結果と閾値との比
較により目標が検出されており、この場合は、定常的に
４dB以上の信号が存在する領域では常に正の値の尤度比
Λ(ｘ)が算出されるので、複数受信信号列の累積により
目標の検出が可能となる。ここで、図５に、信号振幅値
Ａの値が異なる場合でのＳＮ比と確率密度関数および尤
度比Λ(ｘ)の関係を示す。図中では、Ａ＝４.２４（数
式１０よりＳＮ比10dBに相当）を用いている。Ａ＝４.
２４の場合では、確率密度関数のピークはＳＮ比10dB付
近に位置し、また、尤度比Λ(ｘ)はＳＮ比＝６dB以上の
信号に対して正の値が算出される。
【００４１】以上より、ライス分布の分布形状を決定す
るパラメータ、即ち、信号振幅値Ａは信号のＳＮ比を表
すと同時に、尤度比を用いて検出可能な信号のＳＮ比の
最小値を決定する。しかしながら、課題の項でも述べた
通り、１つの信号振幅値Ａでは、検出可能なＳＮ比の下
限は決定できるが、上限は決定できず、それ故、検出さ
れた目標信号の平均ＳＮ比の算出は困難になっている。
そこで、信号振幅値Ａの値を１つではなく、２つ以上設
定することを考える。例えば、図５に示すように、Ａ＝
２.４５とＡ＝４.２４の２通り設定すれば、図５より、
Ａ＝４.２４の場合は、ＳＮ比＝６dB以上の信号に対し
て、また、Ａ＝２.４５の場合には、ＳＮ比＝４dB以上
の信号に対して正の値の尤度比が算出される。従って、
４dB〜６dBの信号が定常的に存在する領域に対しては、
Ａ＝２.４５のライス分布の確率密度関数を尤度比算出
式に用いた場合は、常に正の値の尤度比Λ(ｘ)が算出さ
れるため、複数受信信号の累積により目標信号として検
出可能であるが、Ａ＝４.２４の場合は、４dB〜６dBの
信号が定常的に存在しても、算出される尤度比Λ(ｘ)は
負の値なので、閾値η1 を超えることはなく、目標信号
としての検出は不可となる。一方、６dB以上の信号が定
常的に存在する領域に対しては、Ａの値が何れであって
も正の値の尤度比Λ(ｘ)が算出されるため、目標信号と
して検出可能である。換言すれば、尤度比算出式がＡ＝
２.４５とＡ＝４.２４の２通りに設定された状態で目標
検出処理が行われた場合、Ａ＝２.４５の尤度比算出式
では、目標信号として検出され得るも、Ａ＝４.２４の
尤度比算出式では、目標信号として検出され得なかった
領域には、定常的に４〜６dBの信号が存在していると判
断し得るものである。
るパラメータ、即ち、信号振幅値Ａは信号のＳＮ比を表
すと同時に、尤度比を用いて検出可能な信号のＳＮ比の
最小値を決定する。しかしながら、課題の項でも述べた
通り、１つの信号振幅値Ａでは、検出可能なＳＮ比の下
限は決定できるが、上限は決定できず、それ故、検出さ
れた目標信号の平均ＳＮ比の算出は困難になっている。
そこで、信号振幅値Ａの値を１つではなく、２つ以上設
定することを考える。例えば、図５に示すように、Ａ＝
２.４５とＡ＝４.２４の２通り設定すれば、図５より、
Ａ＝４.２４の場合は、ＳＮ比＝６dB以上の信号に対し
て、また、Ａ＝２.４５の場合には、ＳＮ比＝４dB以上
の信号に対して正の値の尤度比が算出される。従って、
４dB〜６dBの信号が定常的に存在する領域に対しては、
Ａ＝２.４５のライス分布の確率密度関数を尤度比算出
式に用いた場合は、常に正の値の尤度比Λ(ｘ)が算出さ
れるため、複数受信信号の累積により目標信号として検
出可能であるが、Ａ＝４.２４の場合は、４dB〜６dBの
信号が定常的に存在しても、算出される尤度比Λ(ｘ)は
負の値なので、閾値η1 を超えることはなく、目標信号
としての検出は不可となる。一方、６dB以上の信号が定
常的に存在する領域に対しては、Ａの値が何れであって
も正の値の尤度比Λ(ｘ)が算出されるため、目標信号と
して検出可能である。換言すれば、尤度比算出式がＡ＝
２.４５とＡ＝４.２４の２通りに設定された状態で目標
検出処理が行われた場合、Ａ＝２.４５の尤度比算出式
では、目標信号として検出され得るも、Ａ＝４.２４の
尤度比算出式では、目標信号として検出され得なかった
領域には、定常的に４〜６dBの信号が存在していると判
断し得るものである。
【００４２】ところで、以上の例では、説明の簡単化
上、信号振幅値Ａは２通り分、設定されているが、一般
に、信号振幅値Ａが多数設定されることで、信号はその
ＳＮ比の大きさ毎に分類され得、また、信号振幅値Ａ間
の間隔が小さく設定される程に、信号はそのＳＮ比の大
きさが細かく分類され得るものである。
上、信号振幅値Ａは２通り分、設定されているが、一般
に、信号振幅値Ａが多数設定されることで、信号はその
ＳＮ比の大きさ毎に分類され得、また、信号振幅値Ａ間
の間隔が小さく設定される程に、信号はそのＳＮ比の大
きさが細かく分類され得るものである。
【００４３】以上、本発明の理論的な背景について説明
した。その背景を踏まえ、以下、本発明を具体的に説明
すれば、図１に本発明による目標検出装置の一例での要
部構成を示す。これによる場合、その要部としての尤度
比算出器６０は本発明に直接係わるものとされた上、計
算器６１、乗算器６２および判定器６３から構成されて
いるが、このうち、計算器６１は複数の計算器６１ｉ
(ｉ＝１,…,ｅ：以下省略)から、また、乗算器６２は複
数の乗算器６２ｉから、更に、判定器６３は複数の判定
器６３ｉと比較器６４とからそれぞれ構成されたものと
なっている。計算器６１ｉには目標捜索領域内に含まれ
る受信信号振幅値ｘおよび分散計算領域内に含まれる受
信信号振幅値の分散値σ2 が入力される一方、比較器６
４からは目標が存在すると判定された領域について、用
いた尤度比算出器毎の分類結果が出力されているもので
あるが、尤度比算出器６０での動作をより詳細に説明す
れば以下のようである。
した。その背景を踏まえ、以下、本発明を具体的に説明
すれば、図１に本発明による目標検出装置の一例での要
部構成を示す。これによる場合、その要部としての尤度
比算出器６０は本発明に直接係わるものとされた上、計
算器６１、乗算器６２および判定器６３から構成されて
いるが、このうち、計算器６１は複数の計算器６１ｉ
(ｉ＝１,…,ｅ：以下省略)から、また、乗算器６２は複
数の乗算器６２ｉから、更に、判定器６３は複数の判定
器６３ｉと比較器６４とからそれぞれ構成されたものと
なっている。計算器６１ｉには目標捜索領域内に含まれ
る受信信号振幅値ｘおよび分散計算領域内に含まれる受
信信号振幅値の分散値σ2 が入力される一方、比較器６
４からは目標が存在すると判定された領域について、用
いた尤度比算出器毎の分類結果が出力されているもので
あるが、尤度比算出器６０での動作をより詳細に説明す
れば以下のようである。
【００４４】即ち、計算器６１ｉには、時間的に変動し
ている複数の受信信号列Ｓ1(ｔ)、Ｓ2(ｔ)、…の時間ｔ
における目標捜索領域内に含まれる受信信号振幅値ｘt
および該目標捜索領域に対応する分散計算領域内に含ま
れる受信信号振幅値の分散値σt 2が順次入力されるもの
となっている。以下、ｍ番目に入力される受信信号列を
Ｓm(ｔ)として記すものとして、計算器６１ｉ各々で
は、受信信号振幅値ｘtおよび分散値σt 2が所定に処理
されるものとなっている。計算器６１ｉ各々にはそれぞ
れ異なる信号振幅値Ａのライス分布に基づく尤度比算出
式が事前設定されているが、これら尤度比算出式に受信
信号振幅値ｘtおよび分散値σt 2が代入されることによ
って、ｍ番目の受信信号列Ｓm(ｔ)における尤度比Λ
m(ｘt；ｉ)が算出された上、乗算器６２ｉへ出力されて
いるものである。
ている複数の受信信号列Ｓ1(ｔ)、Ｓ2(ｔ)、…の時間ｔ
における目標捜索領域内に含まれる受信信号振幅値ｘt
および該目標捜索領域に対応する分散計算領域内に含ま
れる受信信号振幅値の分散値σt 2が順次入力されるもの
となっている。以下、ｍ番目に入力される受信信号列を
Ｓm(ｔ)として記すものとして、計算器６１ｉ各々で
は、受信信号振幅値ｘtおよび分散値σt 2が所定に処理
されるものとなっている。計算器６１ｉ各々にはそれぞ
れ異なる信号振幅値Ａのライス分布に基づく尤度比算出
式が事前設定されているが、これら尤度比算出式に受信
信号振幅値ｘtおよび分散値σt 2が代入されることによ
って、ｍ番目の受信信号列Ｓm(ｔ)における尤度比Λ
m(ｘt；ｉ)が算出された上、乗算器６２ｉへ出力されて
いるものである。
【００４５】その乗算器６２ｉではまた、各受信信号列
Ｓ1(ｔ)、Ｓ2(ｔ)、…、Ｓm(ｔ)の同一位置ｔ(＝Ｔとお
く)における尤度比Λ1(ｘT；ｉ)、Λ2(ｘT；ｉ)、…、
Λm(ｘT；ｉ)が順次累積され累積尤度比Λ(ｘT,ｍ；ｉ)
が算出されるものとなっている。例えば累積尤度比Λ
(ｘT,ｍ；ｉ)が算出されるに際しては、乗算器６２ｉに
は受信信号列Ｓ(m-1)(ｔ)までの位置ｔにおける累積尤
度比Λ(ｘt,ｍ-１；ｉ)がそれまでの累積尤度比として
保持されているが、計算器６１ｉから新たにｍ番目の受
信信号列Ｓm(ｔ)における尤度比Λm(ｘt；ｉ)が入力さ
れれば、位置ｔが一致するΛ(ｘt,ｍ-１；ｉ)が取り出
された上、その新たなΛm(ｘt；ｉ)と乗算されることに
より新たな累積尤度比Λ(ｘt,ｍ；ｉ)が算出されている
ものである。このようにして算出された累積尤度比Λ
(ｘt,ｍ；ｉ)は、判定器６３ｉへ出力される一方では、
乗算器６２ｉにそれまでに保持されていた(ｍ-１)番目
までの累積尤度比Λ(ｘt,ｍ-１；ｉ)の値と置換されて
いるものである。
Ｓ1(ｔ)、Ｓ2(ｔ)、…、Ｓm(ｔ)の同一位置ｔ(＝Ｔとお
く)における尤度比Λ1(ｘT；ｉ)、Λ2(ｘT；ｉ)、…、
Λm(ｘT；ｉ)が順次累積され累積尤度比Λ(ｘT,ｍ；ｉ)
が算出されるものとなっている。例えば累積尤度比Λ
(ｘT,ｍ；ｉ)が算出されるに際しては、乗算器６２ｉに
は受信信号列Ｓ(m-1)(ｔ)までの位置ｔにおける累積尤
度比Λ(ｘt,ｍ-１；ｉ)がそれまでの累積尤度比として
保持されているが、計算器６１ｉから新たにｍ番目の受
信信号列Ｓm(ｔ)における尤度比Λm(ｘt；ｉ)が入力さ
れれば、位置ｔが一致するΛ(ｘt,ｍ-１；ｉ)が取り出
された上、その新たなΛm(ｘt；ｉ)と乗算されることに
より新たな累積尤度比Λ(ｘt,ｍ；ｉ)が算出されている
ものである。このようにして算出された累積尤度比Λ
(ｘt,ｍ；ｉ)は、判定器６３ｉへ出力される一方では、
乗算器６２ｉにそれまでに保持されていた(ｍ-１)番目
までの累積尤度比Λ(ｘt,ｍ-１；ｉ)の値と置換されて
いるものである。
【００４６】一方、判定器６３ｉでは、乗算器６２ｉか
らの累積尤度比Λ(ｘt,ｍ；ｉ)は、数式５により定めた
閾値η0 ，η1 各々との比較により、目標の存否が判定
されたものとなっている。即ち、η0 ＜Λ(ｘt,ｍ；ｉ)
＜η1 ならば、目標の存否判定は不可とされているも、
Λ(ｘt,ｍ；ｉ)≧η1 ならば、位置ｔに目標が存在する
として、また、Λ(ｘt,ｍ；ｉ)≦η0 ならば、位置ｔに
目標が存在しないとしてそれぞれ判定されているもので
ある。
らの累積尤度比Λ(ｘt,ｍ；ｉ)は、数式５により定めた
閾値η0 ，η1 各々との比較により、目標の存否が判定
されたものとなっている。即ち、η0 ＜Λ(ｘt,ｍ；ｉ)
＜η1 ならば、目標の存否判定は不可とされているも、
Λ(ｘt,ｍ；ｉ)≧η1 ならば、位置ｔに目標が存在する
として、また、Λ(ｘt,ｍ；ｉ)≦η0 ならば、位置ｔに
目標が存在しないとしてそれぞれ判定されているもので
ある。
【００４７】上記判定により目標が存在すると判定され
た場合には、目標捜索領域の位置ｔは比較器６４に出力
されているが、比較器６４では、判定器６３ｉからの、
各々の尤度比算出式を用いて目標が存在すると判定され
た目標捜索領域の位置ｔについて、用いた尤度比算出式
毎の分類が行われ、分類結果が出力されているものであ
る。例えば、計算器６１１にＡ＝２.４５のライス分布
に基づく尤度比算出式が、また、計算器６１２にＡ＝
４.２４のライス分布に基づく尤度比算出式がそれぞれ
設定されている場合、判定器６３ｉ(ｉ＝１,２)各々か
らの、目標が存在する旨の判定分類パターンとしては、
２通りの分類パターンが考えられるものとなっている。
先ず分類パターン１としては、判定器６３１（Ａ＝２.
４５）からのみ、目標が存在する旨の判定結果が得られ
る場合であり、他の分類パターン２としては、判定器６
３１、６３２各々から、ともに目標が存在する旨の判定
結果が得られる場合である。
た場合には、目標捜索領域の位置ｔは比較器６４に出力
されているが、比較器６４では、判定器６３ｉからの、
各々の尤度比算出式を用いて目標が存在すると判定され
た目標捜索領域の位置ｔについて、用いた尤度比算出式
毎の分類が行われ、分類結果が出力されているものであ
る。例えば、計算器６１１にＡ＝２.４５のライス分布
に基づく尤度比算出式が、また、計算器６１２にＡ＝
４.２４のライス分布に基づく尤度比算出式がそれぞれ
設定されている場合、判定器６３ｉ(ｉ＝１,２)各々か
らの、目標が存在する旨の判定分類パターンとしては、
２通りの分類パターンが考えられるものとなっている。
先ず分類パターン１としては、判定器６３１（Ａ＝２.
４５）からのみ、目標が存在する旨の判定結果が得られ
る場合であり、他の分類パターン２としては、判定器６
３１、６３２各々から、ともに目標が存在する旨の判定
結果が得られる場合である。
【００４８】さて、以上の如くにしてなる尤度比算出器
６０は、尤度比による目標検出処理が行われている目標
検出装置、あるいは尤度比演算器に適用されるが、図６
にその尤度比算出器６０を含む尤度比演算器の一例での
概要構成を示す。図示のように、尤度比演算器７０は、
領域設定器７１、分散算出器７２および尤度比算出器６
０から構成されたものとなっている。尤度比演算器７０
への入力は受信信号列Ｓm(ｔ)とされ、その出力は、目
標が存在すると判定された目標捜索領域の位置とされて
いるわけであるが、その動作を、図７に示す時間（ある
いは音波、または電波の単位時間当たりの伝搬距離を乗
じることにより求まる距離）方向に変動する受信信号列
を例に採り説明すれば以下のようである。
６０は、尤度比による目標検出処理が行われている目標
検出装置、あるいは尤度比演算器に適用されるが、図６
にその尤度比算出器６０を含む尤度比演算器の一例での
概要構成を示す。図示のように、尤度比演算器７０は、
領域設定器７１、分散算出器７２および尤度比算出器６
０から構成されたものとなっている。尤度比演算器７０
への入力は受信信号列Ｓm(ｔ)とされ、その出力は、目
標が存在すると判定された目標捜索領域の位置とされて
いるわけであるが、その動作を、図７に示す時間（ある
いは音波、または電波の単位時間当たりの伝搬距離を乗
じることにより求まる距離）方向に変動する受信信号列
を例に採り説明すれば以下のようである。
【００４９】即ち、領域設定器７１では、外部から入力
される受信信号列Ｓm(ｔ)に対し、図７に示すように、
時間（距離）方向に一定幅をもつ目標捜索領域とこの目
標捜索領域に対応する分散計算領域とが対として設定さ
れるが、この場合、目標捜索領域はただ１つの受信信号
振幅値ｘを含む幅として設定されるものとなっている。
換言すれば、受信信号列Ｓm(ｔ)を構成している受信信
号振幅値ｘ各々が目標捜索領域として設定されているも
のである。受信信号列Ｓm(ｔ)の全領域がカバーされる
べく、対としての目標捜索領域および分散計算領域は時
間（距離）方向に順次シフト設定されているわけである
が、これら両領域がシフト設定される毎に、目標捜索領
域に含まれている受信信号振幅値ｘは尤度比算出器６０
に、また、分散計算領域に含まれている複数の受信信号
振幅値ｘは分散算出器７２にそれぞれ出力されているも
のである。その分散算出器７２では、領域設定器７１か
らの、分散計算領域内に含まれている複数の受信信号振
幅値に基づき、以下の数式１１によりその分散計算領域
上での受信信号振幅値ｘの分散値σ2 が計算された上、
尤度比算出器６０に出力されているものである。
される受信信号列Ｓm(ｔ)に対し、図７に示すように、
時間（距離）方向に一定幅をもつ目標捜索領域とこの目
標捜索領域に対応する分散計算領域とが対として設定さ
れるが、この場合、目標捜索領域はただ１つの受信信号
振幅値ｘを含む幅として設定されるものとなっている。
換言すれば、受信信号列Ｓm(ｔ)を構成している受信信
号振幅値ｘ各々が目標捜索領域として設定されているも
のである。受信信号列Ｓm(ｔ)の全領域がカバーされる
べく、対としての目標捜索領域および分散計算領域は時
間（距離）方向に順次シフト設定されているわけである
が、これら両領域がシフト設定される毎に、目標捜索領
域に含まれている受信信号振幅値ｘは尤度比算出器６０
に、また、分散計算領域に含まれている複数の受信信号
振幅値ｘは分散算出器７２にそれぞれ出力されているも
のである。その分散算出器７２では、領域設定器７１か
らの、分散計算領域内に含まれている複数の受信信号振
幅値に基づき、以下の数式１１によりその分散計算領域
上での受信信号振幅値ｘの分散値σ2 が計算された上、
尤度比算出器６０に出力されているものである。
【００５０】
【数１１】
【００５１】その尤度比演算器６０では、既述のよう
に、領域設定器７１からの受信信号振幅値ｘと分散計算
器７２からの分散値σ2 とに基づき、尤度比を用いた目
標検出が行われているものである。
に、領域設定器７１からの受信信号振幅値ｘと分散計算
器７２からの分散値σ2 とに基づき、尤度比を用いた目
標検出が行われているものである。
【００５２】また、図８に尤度比演算器の他の例での概
要構成を示す。図示のように、尤度比演算器９０は、領
域設定器９１、分散算出器９２、平均値算出器９３およ
び尤度比算出器９４から構成されているが、図６に示す
ものとの大きな相違は、領域設定器７１と尤度比算出器
６０との間に平均値算出器９３が新たに介在設置されて
いることである。領域設定器９１では、外部から入力さ
れる受信信号列Ｓm(ｔ)に対し、図６の場合と同様、時
間（距離）方向に一定幅をもつ目標捜索領域とこの目標
捜索領域に対する分散計算領域とが対として設定されて
いるわけであるが、図６の場合とは異なり、目標捜索領
域各々はＮ（≠１）個の受信信号振幅値ｘを含む幅とし
て設定されたものとなっている。図６の場合と同様、受
信信号列Ｓm(ｔ)の全領域がカバーされるべく、対とし
ての目標捜索領域および分散計算領域は時間（距離）方
向に順次シフト設定されているが、これら両領域がシフ
ト設定される毎に、その目標捜索領域に含まれているＮ
個の受信信号振幅値ｘはその２乗平均値ｘ2 が平均値算
出器９３で計算された上、尤度比算出器９４に出力され
ている一方、分散計算領域に含まれている複数の受信信
号振幅値ｘはその分散値σ2 が分散算出器９２で計算さ
れた上、尤度比算出器９４に出力されているものであ
る。尤度比算出器９４では、平均値算出器９３からの２
乗平均値ｘ2 と分散算出器９２からの分散値σ2 とに基
づき目標検出処理が行われているわけであるが、この尤
度比算出器９４自体の構成は、図１に示す尤度比算出器
６０とほぼ同様ながらも、計算器６１ｉ各々に設定され
る尤度比算出式は、図１の場合とは異なるもの（後述の
数式１８として定義される尤度比算出式）として設定さ
れたものとなっている。ここで、レイリー分布およびラ
イス分布に従う確率変数の２乗平均値が従う確率密度関
数の導出過程および尤度比算出式の設定方法を述べれ
ば、以下のようである。
要構成を示す。図示のように、尤度比演算器９０は、領
域設定器９１、分散算出器９２、平均値算出器９３およ
び尤度比算出器９４から構成されているが、図６に示す
ものとの大きな相違は、領域設定器７１と尤度比算出器
６０との間に平均値算出器９３が新たに介在設置されて
いることである。領域設定器９１では、外部から入力さ
れる受信信号列Ｓm(ｔ)に対し、図６の場合と同様、時
間（距離）方向に一定幅をもつ目標捜索領域とこの目標
捜索領域に対する分散計算領域とが対として設定されて
いるわけであるが、図６の場合とは異なり、目標捜索領
域各々はＮ（≠１）個の受信信号振幅値ｘを含む幅とし
て設定されたものとなっている。図６の場合と同様、受
信信号列Ｓm(ｔ)の全領域がカバーされるべく、対とし
ての目標捜索領域および分散計算領域は時間（距離）方
向に順次シフト設定されているが、これら両領域がシフ
ト設定される毎に、その目標捜索領域に含まれているＮ
個の受信信号振幅値ｘはその２乗平均値ｘ2 が平均値算
出器９３で計算された上、尤度比算出器９４に出力され
ている一方、分散計算領域に含まれている複数の受信信
号振幅値ｘはその分散値σ2 が分散算出器９２で計算さ
れた上、尤度比算出器９４に出力されているものであ
る。尤度比算出器９４では、平均値算出器９３からの２
乗平均値ｘ2 と分散算出器９２からの分散値σ2 とに基
づき目標検出処理が行われているわけであるが、この尤
度比算出器９４自体の構成は、図１に示す尤度比算出器
６０とほぼ同様ながらも、計算器６１ｉ各々に設定され
る尤度比算出式は、図１の場合とは異なるもの（後述の
数式１８として定義される尤度比算出式）として設定さ
れたものとなっている。ここで、レイリー分布およびラ
イス分布に従う確率変数の２乗平均値が従う確率密度関
数の導出過程および尤度比算出式の設定方法を述べれ
ば、以下のようである。
【００５３】即ち、例えば領域内にＮ個の受信信号振幅
値を含む大きさの目標捜索領域が設定される場合、その
領域内の各振幅値の２乗平均値をｘ2 とおくと、各振幅
値の２乗和は(Ｎ*ｘ2 )で表される。ところで、レイリ
ー分布に従う確率変数Ｒは、平均値０、分散値σ2 の正
規分布に従う確率変数Ｘ，Ｙを用いて以下の数式１２で
表される。
値を含む大きさの目標捜索領域が設定される場合、その
領域内の各振幅値の２乗平均値をｘ2 とおくと、各振幅
値の２乗和は(Ｎ*ｘ2 )で表される。ところで、レイリ
ー分布に従う確率変数Ｒは、平均値０、分散値σ2 の正
規分布に従う確率変数Ｘ，Ｙを用いて以下の数式１２で
表される。
【００５４】
【数１２】
【００５５】数式１２より、Ｎ個の振幅値の２乗和ΣＲ
i 2は、以下の数式１３として表現され得る。
i 2は、以下の数式１３として表現され得る。
【００５６】
【数１３】
【００５７】従って、２乗和Ｎｘ2 を分散値σ2 で正規
化することで得られる確率変数χ2は、以下に数式１４
として示す自由度２Ｎのカイ２乗確率密度関数に従う。
化することで得られる確率変数χ2は、以下に数式１４
として示す自由度２Ｎのカイ２乗確率密度関数に従う。
【００５８】
【数１４】
【００５９】確率変数Ｐが数式７として示すライス分布
に従う場合、確率変数Ｐは、平均値２-1/2Ａ、分散値σ
2 の正規分布に従う確率変数Ｗ，Ｚを用いて以下の数式
１５で表される。
に従う場合、確率変数Ｐは、平均値２-1/2Ａ、分散値σ
2 の正規分布に従う確率変数Ｗ，Ｚを用いて以下の数式
１５で表される。
【００６０】
【数１５】
【００６１】数式１５より、Ｎ個の振幅値の２乗和ΣＰ
i 2は、以下の数式１６として表し得るものである。
i 2は、以下の数式１６として表し得るものである。
【００６２】
【数１６】
【００６３】従って、２乗和Ｎｘ2 を分散値σ2 で正規
化することで得られる確率変数κ2は、以下に数式１７
として示す自由度２Ｎの非心カイ２乗確率密度関数(non
center chi-square p.d.f)に従う。
化することで得られる確率変数κ2は、以下に数式１７
として示す自由度２Ｎの非心カイ２乗確率密度関数(non
center chi-square p.d.f)に従う。
【００６４】
【数１７】
【００６５】故に、尤度比算出器９４に設定される尤度
比算出式は、以下の数式１８として定義される。
比算出式は、以下の数式１８として定義される。
【００６６】
【数１８】
【００６７】ところで、数式１８で定義されている尤度
比Λ(ｘ)の値は、４つのパラメータ、即ち、目標捜索領
域に含まれる受信信号振幅値の２乗平均値ｘ2 、分散計
算領域に含まれる受信信号振幅値の分散値σ2 、信号振
幅値Ａおよび目標捜索領域内に含まれる受信信号振幅値
の数Ｎによって決定されているが、信号振幅値Ａおよび
受信信号振幅値数Ｎは事前設定定数であることから、結
局、２乗平均値ｘ2 および分散値σ2 から、その目標捜
索領域における尤度比Λ(ｘ)が求まり、数式５により予
め定められている閾値η1 ，η0 各々との比較により目
標信号の検出が可能とされているものである。
比Λ(ｘ)の値は、４つのパラメータ、即ち、目標捜索領
域に含まれる受信信号振幅値の２乗平均値ｘ2 、分散計
算領域に含まれる受信信号振幅値の分散値σ2 、信号振
幅値Ａおよび目標捜索領域内に含まれる受信信号振幅値
の数Ｎによって決定されているが、信号振幅値Ａおよび
受信信号振幅値数Ｎは事前設定定数であることから、結
局、２乗平均値ｘ2 および分散値σ2 から、その目標捜
索領域における尤度比Λ(ｘ)が求まり、数式５により予
め定められている閾値η1 ，η0 各々との比較により目
標信号の検出が可能とされているものである。
【００６８】その目標信号の検出についてより詳細に説
明すれば、図９（ａ）に数式１４、数式１６それぞれで
定義されているカイ２乗分布、非心カイ２乗分布（Ａ＝
２.４５、４.２４の場合）の分布形状を示す。図示のよ
うに、２乗平均値ｘ2 と確率密度関数の関係が、横軸は
２乗平均値ｘ2 、縦軸は確率密度として示されている。
また、図９（ｂ）には、数式１８により求まる尤度比Λ
(ｘ)とＳＮ比の関係が、横軸はＳＮ比(dB)、縦軸は尤度
比Λ(ｘ)として示されているが、図９（ｂ）より、尤度
比算出式に数式１８を用いた場合でも、尤度比を用いて
検出可能な信号の最小ＳＮ比（正の値の尤度比Λ(ｘ)を
算出し得るＳＮ比の最小値）は、パラメータとしての信
号振幅値Ａによって変化するものであることが判る。従
って、信号振幅値Ａが異なる非心カイ２乗分布が尤度比
算出式として複数設定される場合であっても、目標が存
在すると判定された領域での平均的なＳＮ比が算出可能
となるものである。
明すれば、図９（ａ）に数式１４、数式１６それぞれで
定義されているカイ２乗分布、非心カイ２乗分布（Ａ＝
２.４５、４.２４の場合）の分布形状を示す。図示のよ
うに、２乗平均値ｘ2 と確率密度関数の関係が、横軸は
２乗平均値ｘ2 、縦軸は確率密度として示されている。
また、図９（ｂ）には、数式１８により求まる尤度比Λ
(ｘ)とＳＮ比の関係が、横軸はＳＮ比(dB)、縦軸は尤度
比Λ(ｘ)として示されているが、図９（ｂ）より、尤度
比算出式に数式１８を用いた場合でも、尤度比を用いて
検出可能な信号の最小ＳＮ比（正の値の尤度比Λ(ｘ)を
算出し得るＳＮ比の最小値）は、パラメータとしての信
号振幅値Ａによって変化するものであることが判る。従
って、信号振幅値Ａが異なる非心カイ２乗分布が尤度比
算出式として複数設定される場合であっても、目標が存
在すると判定された領域での平均的なＳＮ比が算出可能
となるものである。
【００６９】次に、本発明による目標検出装置の一例で
の全体構成を図１０に示すが、その構成と動作概要を説
明すれば以下のようである。即ち、目標検出装置１１０
自体は船舶に搭載された状態で、水中に存在する目標を
検出すべく機能したものとなっている。図１０に示すよ
うに、送信装置１１１で周期的に生成されている送信信
号は送受波装置１１２で超音波に変換された上、水中内
に送波される一方、送信超音波に対する反射超音波はそ
の送受波装置１１２で受波され、電気信号に変換された
状態として受信装置１１３で所定に前処理されるものと
なっている。受信装置１１３では、送受波装置１１２か
らの受信信号は電力増幅・ディジタル変換された上、複
数方位対応の音響ビームとして形成されているが、これ
ら音響ビーム各々が信号処理装置１１４で所定に信号処
理されているものである。信号処理装置１１４での信号
処理結果は表示装置１１５上に可視画像として表示され
ることによって、その表示より複数方位の何れに目標が
存在しているかや、その目標の位置・ＳＮ比がほぼ知れ
るように構成されているものである。また、目標検出装
置１１０とは別に、船舶には速度測定装置１１６および
水深測定装置１１７が搭載され、これら装置１１６，１
１７各々からの測定速度、測定深度は信号処理装置１１
４に取込まれた上、信号処理に供されているが、本発明
に直接係る尤度比演算器はこの信号処理装置１１４内に
含まれたものとなっている。
の全体構成を図１０に示すが、その構成と動作概要を説
明すれば以下のようである。即ち、目標検出装置１１０
自体は船舶に搭載された状態で、水中に存在する目標を
検出すべく機能したものとなっている。図１０に示すよ
うに、送信装置１１１で周期的に生成されている送信信
号は送受波装置１１２で超音波に変換された上、水中内
に送波される一方、送信超音波に対する反射超音波はそ
の送受波装置１１２で受波され、電気信号に変換された
状態として受信装置１１３で所定に前処理されるものと
なっている。受信装置１１３では、送受波装置１１２か
らの受信信号は電力増幅・ディジタル変換された上、複
数方位対応の音響ビームとして形成されているが、これ
ら音響ビーム各々が信号処理装置１１４で所定に信号処
理されているものである。信号処理装置１１４での信号
処理結果は表示装置１１５上に可視画像として表示され
ることによって、その表示より複数方位の何れに目標が
存在しているかや、その目標の位置・ＳＮ比がほぼ知れ
るように構成されているものである。また、目標検出装
置１１０とは別に、船舶には速度測定装置１１６および
水深測定装置１１７が搭載され、これら装置１１６，１
１７各々からの測定速度、測定深度は信号処理装置１１
４に取込まれた上、信号処理に供されているが、本発明
に直接係る尤度比演算器はこの信号処理装置１１４内に
含まれたものとなっている。
【００７０】ここで、受信装置１１３、信号処理装置１
１４、表示装置１１５各々についてより詳細に説明す
る。先ず受信装置１１３であるが、受信装置１１３で
は、送受波装置１１２からの受信信号は電力増幅された
上、図１１（ａ），（ｂ）にディジタル変換前、ディジ
タル変換後の状態として示すように、事前設定されたサ
ンプリング周期ｔs でディジタル信号に変換されるもの
となっている。更に、ディジタル変換後の受信信号に対
する整相処理により、図１２に示すように、一定の方位
範囲内に一定数のビームが、方位方向対応受信信号波形
として放射状に形成されるものとなっている。図１３に
はまた、送信装置１１１からの送信信号波形と受信装置
１１３で形成されている方位方向対応受信信号波形との
関係が、縦軸は受信信号の振幅値、横軸は時間（あるい
は音波の単位時間当たりの水中伝搬距離を乗じることに
より求まる距離）として示されているが、これを、時間
（距離）方向および方位方向に２次元状に配列したもの
を図１４として示す。以下、ｍ番目の送信信号に対する
受信信号列をＳm(ｎ,ｔ)として表すことにする。但し、
ｎは方位を、ｔは時間を表す。
１４、表示装置１１５各々についてより詳細に説明す
る。先ず受信装置１１３であるが、受信装置１１３で
は、送受波装置１１２からの受信信号は電力増幅された
上、図１１（ａ），（ｂ）にディジタル変換前、ディジ
タル変換後の状態として示すように、事前設定されたサ
ンプリング周期ｔs でディジタル信号に変換されるもの
となっている。更に、ディジタル変換後の受信信号に対
する整相処理により、図１２に示すように、一定の方位
範囲内に一定数のビームが、方位方向対応受信信号波形
として放射状に形成されるものとなっている。図１３に
はまた、送信装置１１１からの送信信号波形と受信装置
１１３で形成されている方位方向対応受信信号波形との
関係が、縦軸は受信信号の振幅値、横軸は時間（あるい
は音波の単位時間当たりの水中伝搬距離を乗じることに
より求まる距離）として示されているが、これを、時間
（距離）方向および方位方向に２次元状に配列したもの
を図１４として示す。以下、ｍ番目の送信信号に対する
受信信号列をＳm(ｎ,ｔ)として表すことにする。但し、
ｎは方位を、ｔは時間を表す。
【００７１】図１５にはまた、信号処理装置１１４の一
例での構成が示されているが、受信装置１１３からの受
信信号列Ｓm(ｎ,ｔ)は映像信号処理器１６１で所定に映
像信号に変換処理された上、表示装置１１５上に可視画
像として表示されている他、受信信号列Ｓm(ｎ,ｔ)中に
存在している目標のその位置を特定するための処理が尤
度比演算器７０で行われているものである。具体的に、
尤度比演算器７０では、受信装置１１３からの受信信号
列Ｓm(ｎ,ｔ)に対して、尤度比を用いた目標検出処理が
行われており、目標が存在すると判定された目標捜索領
域の位置（方位ｎ、距離ｔ）は表示装置１１５上で可視
表示されているものである。なお、速度測定器１１６、
水深測定器１１７各々からの測定速度、測定深度は、後
述のように、乗算器６２ｉに取り込まれた上、所定に処
理されるものとなっている。
例での構成が示されているが、受信装置１１３からの受
信信号列Ｓm(ｎ,ｔ)は映像信号処理器１６１で所定に映
像信号に変換処理された上、表示装置１１５上に可視画
像として表示されている他、受信信号列Ｓm(ｎ,ｔ)中に
存在している目標のその位置を特定するための処理が尤
度比演算器７０で行われているものである。具体的に、
尤度比演算器７０では、受信装置１１３からの受信信号
列Ｓm(ｎ,ｔ)に対して、尤度比を用いた目標検出処理が
行われており、目標が存在すると判定された目標捜索領
域の位置（方位ｎ、距離ｔ）は表示装置１１５上で可視
表示されているものである。なお、速度測定器１１６、
水深測定器１１７各々からの測定速度、測定深度は、後
述のように、乗算器６２ｉに取り込まれた上、所定に処
理されるものとなっている。
【００７２】その尤度比演算器７０での目標検出処理で
は、具体的に、図１４に示す受信信号列Ｓm(ｎ,ｔ)に対
し目標捜索領域および分散計算領域が設定されるが、そ
れら領域の設定態様を図１６に示す。その際に、目標捜
索領域は、ただ１つの振幅値を含む大きさとして設定さ
れつつ、目標捜索領域および分散計算領域が時間および
方位の方向に順次シフト設定されることで、受信信号列
Ｓm(ｎ,ｔ)はその全領域がカバーされているものである
（なお、目標捜索領域および分散計算領域のシフト設定
や処理等に関しては、本願出願人による特願平８ー３３
０８８２号に詳細に記載されている）。尤も、船舶自体
の移動を伴い超音波が水中に送波される場合には、送信
信号対応に得られる受信信号列の座標系は異なるため、
各受信信号列中の同一位置（時間ｔ、方位ｎ）の目標捜
索領域における尤度比同士を乗算するためには、各受信
信号間で位置修正を行うべく、座標変換が行われる必要
がある。図１７にその際でのモデルを示す。図示のよう
に、船舶が右側に位置している際に、ｍ番目の送信信号
が送波され、それに対する受信信号列Ｓm(ｎ,ｔ)では、
目標物体は距離ｒ、方位φに位置しているが、その後、
船舶が左方向に距離Ｘ0 分移動された状態で、(ｍ+１)
番目の送信信号が送波されたとすれば、それに対する受
信信号列Ｓ(m+1)(ｎ,ｔ)では、船舶からすれば、目標物
体は距離ｒ'、方位φ'に位置していることになるという
わけである。よって、受信信号列Ｓm(ｎ,ｔ)の受信信号
列Ｓ(m+1)(ｎ,ｔ)への位置修正が行われる必要があり、
その位置修正に際しては、ｍ番目の受信信号列Ｓm(ｎ,
ｔ)における距離ｔおよび方位ｎは、以下の数式１９に
代入されることによって、(ｍ+１）番目の受信信号列Ｓ
(m+1)(ｎ,ｔ)における距離ｔ'および方位ｎ'に補正され
ることで、船舶が移動状態にある間での目標検出処理が
可能とされているものである。
は、具体的に、図１４に示す受信信号列Ｓm(ｎ,ｔ)に対
し目標捜索領域および分散計算領域が設定されるが、そ
れら領域の設定態様を図１６に示す。その際に、目標捜
索領域は、ただ１つの振幅値を含む大きさとして設定さ
れつつ、目標捜索領域および分散計算領域が時間および
方位の方向に順次シフト設定されることで、受信信号列
Ｓm(ｎ,ｔ)はその全領域がカバーされているものである
（なお、目標捜索領域および分散計算領域のシフト設定
や処理等に関しては、本願出願人による特願平８ー３３
０８８２号に詳細に記載されている）。尤も、船舶自体
の移動を伴い超音波が水中に送波される場合には、送信
信号対応に得られる受信信号列の座標系は異なるため、
各受信信号列中の同一位置（時間ｔ、方位ｎ）の目標捜
索領域における尤度比同士を乗算するためには、各受信
信号間で位置修正を行うべく、座標変換が行われる必要
がある。図１７にその際でのモデルを示す。図示のよう
に、船舶が右側に位置している際に、ｍ番目の送信信号
が送波され、それに対する受信信号列Ｓm(ｎ,ｔ)では、
目標物体は距離ｒ、方位φに位置しているが、その後、
船舶が左方向に距離Ｘ0 分移動された状態で、(ｍ+１)
番目の送信信号が送波されたとすれば、それに対する受
信信号列Ｓ(m+1)(ｎ,ｔ)では、船舶からすれば、目標物
体は距離ｒ'、方位φ'に位置していることになるという
わけである。よって、受信信号列Ｓm(ｎ,ｔ)の受信信号
列Ｓ(m+1)(ｎ,ｔ)への位置修正が行われる必要があり、
その位置修正に際しては、ｍ番目の受信信号列Ｓm(ｎ,
ｔ)における距離ｔおよび方位ｎは、以下の数式１９に
代入されることによって、(ｍ+１）番目の受信信号列Ｓ
(m+1)(ｎ,ｔ)における距離ｔ'および方位ｎ'に補正され
ることで、船舶が移動状態にある間での目標検出処理が
可能とされているものである。
【００７３】
【数１９】
【００７４】なお、数式１９中、受信信号間の移動距離
Ｘ0 の値としては、速度測定装置１１６から得られる速
度に経過時間を乗じて得られる値が、また、水深ｚの値
としては、水深測定装置１１７から得られる深度値がそ
れぞれ用いられる。
Ｘ0 の値としては、速度測定装置１１６から得られる速
度に経過時間を乗じて得られる値が、また、水深ｚの値
としては、水深測定装置１１７から得られる深度値がそ
れぞれ用いられる。
【００７５】ここで、サフィックス（ｎ，ｔ）はサフィ
ックス（ｎｔ）に、また、サフィックス（ｎ'，ｔ'）は
サフィックス（ＮＴ）にそれぞれ置換されているとし
て、乗算器６２ｉ各々での計算処理方法を説明すれば以
下のようである。即ち、それまでに保持されているｍ番
目までの受信信号列Ｓm(ｎ,ｔ)中の距離ｔ、方位ｎにお
ける累積尤度比Λ(ｘ(nt),ｍ；ｉ)のｎおよびｔは数式
１９によりｎ'およびｔ'に変換された上、計算器６１ｉ
からの新たな(ｍ+１)番目の受信信号列Ｓ(m+1)(ｎ,ｔ)
おける尤度比Λ(m+1)(ｘ(nt)；ｉ)のうち、ｎ＝ｎ'およ
びｔ＝ｔ'である尤度比Λ(m+1)(ｘ(NT)；ｉ)と乗算され
ることで、(ｍ+１)番目までの受信信号列Ｓ(m+1)(ｎ,
ｔ)中の距離ｔ＝ｔ'、方位ｎ＝ｎ'における累積尤度比
Λ(ｘ(NT),ｍ+１；ｉ)が算出されているものである。こ
の計算処理結果は判定器６３ｉに出力されると同時に、
乗算器６２ｉにそれまで保持されていたΛ(ｘ( NT),ｍ；
ｉ)（ｍ番目の受信信号列Ｓm(ｎ,ｔ)の方位ｎ'、位置
ｔ'における累積尤度比）と置換されているものであ
る。
ックス（ｎｔ）に、また、サフィックス（ｎ'，ｔ'）は
サフィックス（ＮＴ）にそれぞれ置換されているとし
て、乗算器６２ｉ各々での計算処理方法を説明すれば以
下のようである。即ち、それまでに保持されているｍ番
目までの受信信号列Ｓm(ｎ,ｔ)中の距離ｔ、方位ｎにお
ける累積尤度比Λ(ｘ(nt),ｍ；ｉ)のｎおよびｔは数式
１９によりｎ'およびｔ'に変換された上、計算器６１ｉ
からの新たな(ｍ+１)番目の受信信号列Ｓ(m+1)(ｎ,ｔ)
おける尤度比Λ(m+1)(ｘ(nt)；ｉ)のうち、ｎ＝ｎ'およ
びｔ＝ｔ'である尤度比Λ(m+1)(ｘ(NT)；ｉ)と乗算され
ることで、(ｍ+１)番目までの受信信号列Ｓ(m+1)(ｎ,
ｔ)中の距離ｔ＝ｔ'、方位ｎ＝ｎ'における累積尤度比
Λ(ｘ(NT),ｍ+１；ｉ)が算出されているものである。こ
の計算処理結果は判定器６３ｉに出力されると同時に、
乗算器６２ｉにそれまで保持されていたΛ(ｘ( NT),ｍ；
ｉ)（ｍ番目の受信信号列Ｓm(ｎ,ｔ)の方位ｎ'、位置
ｔ'における累積尤度比）と置換されているものであ
る。
【００７６】表示装置１１５ではまた、映像信号処理器
１６１からの、映像信号に変換された状態の受信信号列
Ｓm(ｎ,ｔ)は表示画面上に可視画像として表示されてい
るが、この表示に併せて、比較器６４からの、目標が存
在すると判定された目標捜索領域の位置（ｎ,ｔ）とそ
の分類結果が同一表示画面上に表示されているものであ
るが、その表示画面上での画像表示例を図１８に示す。
図示のように、半径方向は距離を表しており、これは、
受信信号列Ｓm(ｎ,ｔ)の距離ｔに相当する。また、円周
方向は方位を表しており、受信信号列Ｓm(ｎ,ｔ)の方位
ｎに相当する。その表示画面上では、受信信号列Ｓ
m(ｎ,ｔ)の振幅値は輝度表示されているが、その際に、
比較器６４からの、目標捜索領域の位置（ｎ,ｔ）の分
類結果が併せて表示されているものである。既に、比較
器６４での機能に関連して説明したように、計算器６１
１にはＡ＝２.４５のライス分布に基づく尤度比算出式
が、また、計算器６１２にはＡ＝４.２４のライス分布
に基づく尤度比算出式が設定されている場合、比較器６
４からの判定分類パターンとして、２通りの分類パター
ン１，２が考えられるものとなっている。図１８に示す
表示例では、分類パターン１として分類された目標捜索
領域の位置はマーカー○で、また、分類パターン２とし
て分類された目標捜索領域の位置はマーカー□で表示さ
れたものとなっている。なお、目標捜索領域の表示方法
としては、上記例のように、マーカー形状で区別する方
法の他、マーカー色で区別する方法、あるいはマーカー
を用いずに、その目標捜索領域の振幅値をカラーで区別
された状態として表示する方法等が容易に考えられるも
のとなっている。また、ある一定条件を満たす特定分類
パターンとして分類された目標捜索領域のみをマーカー
で表示することで、特定のＳＮ比の信号が存在する領域
のみを検出・表示することも可能である。例えば上記分
類パターン１が特定分類パターンとして設定されている
場合を想定すれば、該当目標捜索領域のみ、即ち、平均
的なＳＮ比が４〜６dBの信号が存在している目標捜索領
域のみがマーカー等で表示され得るものである。
１６１からの、映像信号に変換された状態の受信信号列
Ｓm(ｎ,ｔ)は表示画面上に可視画像として表示されてい
るが、この表示に併せて、比較器６４からの、目標が存
在すると判定された目標捜索領域の位置（ｎ,ｔ）とそ
の分類結果が同一表示画面上に表示されているものであ
るが、その表示画面上での画像表示例を図１８に示す。
図示のように、半径方向は距離を表しており、これは、
受信信号列Ｓm(ｎ,ｔ)の距離ｔに相当する。また、円周
方向は方位を表しており、受信信号列Ｓm(ｎ,ｔ)の方位
ｎに相当する。その表示画面上では、受信信号列Ｓ
m(ｎ,ｔ)の振幅値は輝度表示されているが、その際に、
比較器６４からの、目標捜索領域の位置（ｎ,ｔ）の分
類結果が併せて表示されているものである。既に、比較
器６４での機能に関連して説明したように、計算器６１
１にはＡ＝２.４５のライス分布に基づく尤度比算出式
が、また、計算器６１２にはＡ＝４.２４のライス分布
に基づく尤度比算出式が設定されている場合、比較器６
４からの判定分類パターンとして、２通りの分類パター
ン１，２が考えられるものとなっている。図１８に示す
表示例では、分類パターン１として分類された目標捜索
領域の位置はマーカー○で、また、分類パターン２とし
て分類された目標捜索領域の位置はマーカー□で表示さ
れたものとなっている。なお、目標捜索領域の表示方法
としては、上記例のように、マーカー形状で区別する方
法の他、マーカー色で区別する方法、あるいはマーカー
を用いずに、その目標捜索領域の振幅値をカラーで区別
された状態として表示する方法等が容易に考えられるも
のとなっている。また、ある一定条件を満たす特定分類
パターンとして分類された目標捜索領域のみをマーカー
で表示することで、特定のＳＮ比の信号が存在する領域
のみを検出・表示することも可能である。例えば上記分
類パターン１が特定分類パターンとして設定されている
場合を想定すれば、該当目標捜索領域のみ、即ち、平均
的なＳＮ比が４〜６dBの信号が存在している目標捜索領
域のみがマーカー等で表示され得るものである。
【００７７】以上のように、図１０に示す目標検出装置
による場合は、受信信号列に尤度比を用いた目標検出法
が適用された上、目標が存在すると判定された領域の位
置が受信信号列の画像上に表示されるに際して、その目
標のＳＮ比も併せて表示され得るものとなっている。
による場合は、受信信号列に尤度比を用いた目標検出法
が適用された上、目標が存在すると判定された領域の位
置が受信信号列の画像上に表示されるに際して、その目
標のＳＮ比も併せて表示され得るものとなっている。
【００７８】次に、本発明による目標検出装置の他の例
での要部としての構成、即ち、その要部としての信号処
理装置１１４の構成について説明すれば、図１９にその
構成を示す。図示のように、図１５に示すものと比し、
尤度比演算器７０は尤度比演算器９０に置換されている
こと以外は、事情は図１０や図１５に示すものに同様と
なっている。既に、図１１（ａ），（ｂ）で説明したよ
うに、受信装置１１３では、受信信号はサンプリング周
期ｔs でディジタル信号に変換された上、信号処理装置
１１４へ出力されているが、その際に、図１１（ａ），
（ｂ）に示すように、受信信号に含まれる目標信号の時
間幅がサンプリング周期ｔs 程度と短い場合には、目標
信号はディジタル変換後の受信信号列中にただ１つの振
幅値として出現するものとなっている。これに対し、図
２０（ａ），（ｂ）に示すように、目標信号の時間幅が
サンプリング周期ｔs より長い場合には、目標信号はデ
ィジタル変換後の受信信号列中に複数の振幅値として連
続的に出現するものとなっている。これと同様に、受信
信号列が２次元（時間および方位）に分布している状態
で、目標信号の方位幅が受信信号列のビーム幅（図１２
を参照のこと）より大きい場合にも、目標信号は複数の
振幅値に跨がって現れるものとなっている。
での要部としての構成、即ち、その要部としての信号処
理装置１１４の構成について説明すれば、図１９にその
構成を示す。図示のように、図１５に示すものと比し、
尤度比演算器７０は尤度比演算器９０に置換されている
こと以外は、事情は図１０や図１５に示すものに同様と
なっている。既に、図１１（ａ），（ｂ）で説明したよ
うに、受信装置１１３では、受信信号はサンプリング周
期ｔs でディジタル信号に変換された上、信号処理装置
１１４へ出力されているが、その際に、図１１（ａ），
（ｂ）に示すように、受信信号に含まれる目標信号の時
間幅がサンプリング周期ｔs 程度と短い場合には、目標
信号はディジタル変換後の受信信号列中にただ１つの振
幅値として出現するものとなっている。これに対し、図
２０（ａ），（ｂ）に示すように、目標信号の時間幅が
サンプリング周期ｔs より長い場合には、目標信号はデ
ィジタル変換後の受信信号列中に複数の振幅値として連
続的に出現するものとなっている。これと同様に、受信
信号列が２次元（時間および方位）に分布している状態
で、目標信号の方位幅が受信信号列のビーム幅（図１２
を参照のこと）より大きい場合にも、目標信号は複数の
振幅値に跨がって現れるものとなっている。
【００７９】さて、図１５に示す信号処理装置１１４で
は、受信信号列に対して、ただ１つの振幅値を含む大き
さの目標捜索領域が設定されており、その目標捜索領域
に含まれる振幅値とその周辺に設定されている分散計算
領域上での振幅値の分散値から、数式８により尤度比が
算出されたものとなっている。しかしながら、上述のよ
うに、目標信号が１振幅値ではなく、複数の振幅値に連
続的に跨がっている場合には、１振幅値を用いるより
も、複数の振幅値を含む目標捜索領域が設定された上、
その領域内の各振幅値をパラメータとして用いるのがよ
り実状に即しており、目標信号がより容易に検出可能に
なると考えられる。そこで、図１９に示す信号処理装置
１１４では、尤度比演算器７０の代りに尤度比演算器９
０が設置されたことによって、数式１８で定義されてい
る、複数の受信信号振幅値の２乗平均値ｘ2 および周辺
区間の振幅値の分散値σ2 をパラメータとした尤度比算
出式を用いた目標検出処理が可能とされているものであ
る。受信信号が２次元（時間および方位）に分布してい
る場合にも、目標捜索領域および分散計算領域は、図１
６で説明したように、設定されればよいものである。尤
も、目標捜索領域は、その領域内にＮ（≠１）個の振幅
値を含む大きさとして設定される。
は、受信信号列に対して、ただ１つの振幅値を含む大き
さの目標捜索領域が設定されており、その目標捜索領域
に含まれる振幅値とその周辺に設定されている分散計算
領域上での振幅値の分散値から、数式８により尤度比が
算出されたものとなっている。しかしながら、上述のよ
うに、目標信号が１振幅値ではなく、複数の振幅値に連
続的に跨がっている場合には、１振幅値を用いるより
も、複数の振幅値を含む目標捜索領域が設定された上、
その領域内の各振幅値をパラメータとして用いるのがよ
り実状に即しており、目標信号がより容易に検出可能に
なると考えられる。そこで、図１９に示す信号処理装置
１１４では、尤度比演算器７０の代りに尤度比演算器９
０が設置されたことによって、数式１８で定義されてい
る、複数の受信信号振幅値の２乗平均値ｘ2 および周辺
区間の振幅値の分散値σ2 をパラメータとした尤度比算
出式を用いた目標検出処理が可能とされているものであ
る。受信信号が２次元（時間および方位）に分布してい
る場合にも、目標捜索領域および分散計算領域は、図１
６で説明したように、設定されればよいものである。尤
も、目標捜索領域は、その領域内にＮ（≠１）個の振幅
値を含む大きさとして設定される。
【００８０】更に、本発明による目標検出装置の更に異
なる他の例での要部としての構成、即ち、その要部とし
ての信号処理装置１１４の構成について説明すれば、図
２１にその構成を示す。図示のように、図１５に示すも
のに比し、映像信号処理器２２２の前段には距離圧縮器
２２１が新たに設置され、また、尤度比演算器７０は尤
度比演算器２２３に置換されたものとなっている。さ
て、受信装置１１３で受信信号列がディジタル変換され
るに際して、目標信号の検出確率を上げるためには、図
２２（ａ）に示すように、サンプリング周期ｔs は小さ
く設定される必要があるが、その半面、サンプリング周
期ｔs が小さく設定される程に、ディジタル変換後の受
信信号列の振幅値数が増大する結果として、表示装置１
１５の表示画面上に表示しきれないという不具合が新た
に生じるものとなっている。これとは逆に、その表示画
面上に収まるデータ数にすべく、サンプリング周期ｔs
が大きく設定される場合にはまた、図２２（ｂ）に示す
ように、目標信号の検出確率が低下することは否めない
ものとなっている。
なる他の例での要部としての構成、即ち、その要部とし
ての信号処理装置１１４の構成について説明すれば、図
２１にその構成を示す。図示のように、図１５に示すも
のに比し、映像信号処理器２２２の前段には距離圧縮器
２２１が新たに設置され、また、尤度比演算器７０は尤
度比演算器２２３に置換されたものとなっている。さ
て、受信装置１１３で受信信号列がディジタル変換され
るに際して、目標信号の検出確率を上げるためには、図
２２（ａ）に示すように、サンプリング周期ｔs は小さ
く設定される必要があるが、その半面、サンプリング周
期ｔs が小さく設定される程に、ディジタル変換後の受
信信号列の振幅値数が増大する結果として、表示装置１
１５の表示画面上に表示しきれないという不具合が新た
に生じるものとなっている。これとは逆に、その表示画
面上に収まるデータ数にすべく、サンプリング周期ｔs
が大きく設定される場合にはまた、図２２（ｂ）に示す
ように、目標信号の検出確率が低下することは否めない
ものとなっている。
【００８１】以上の不具合を解決すべく、図２３
（ａ），（ｂ）に示すように、サンプリング周期ｔs で
ディジタル変換された受信信号列を、Ｒ個の振幅値を含
む区間毎に領域分割した上、それら区間各々における最
大振幅値を検出出力する距離圧縮処理が距離圧縮器２２
１で行われているものである。この処理により表示上で
のデータ量は１／Ｒに圧縮されつつも、目標信号の検出
性能は従来通り維持されているものである。図２１に示
す信号処理装置１１４ではまた、距離圧縮後の受信信号
列に対して尤度比を用いた目標検出処理が可能とされた
ものとなっている。
（ａ），（ｂ）に示すように、サンプリング周期ｔs で
ディジタル変換された受信信号列を、Ｒ個の振幅値を含
む区間毎に領域分割した上、それら区間各々における最
大振幅値を検出出力する距離圧縮処理が距離圧縮器２２
１で行われているものである。この処理により表示上で
のデータ量は１／Ｒに圧縮されつつも、目標信号の検出
性能は従来通り維持されているものである。図２１に示
す信号処理装置１１４ではまた、距離圧縮後の受信信号
列に対して尤度比を用いた目標検出処理が可能とされた
ものとなっている。
【００８２】より具体的に、図２１に示す信号処理装置
１１４での処理について説明すれば、受信装置１１３か
らの受信信号列Ｓm(ｎ,ｔ)は距離圧縮器２２１および尤
度比演算器２２３へ出力されているが、このうち、距離
圧縮器２２１では、受信装置１１３からの受信信号列Ｓ
m(ｎ,ｔ)はＲ個の振幅値を含む区間毎に領域分割された
上、それら区間各々における最大振幅値が検出される、
といった距離圧縮処理が行われており、距離圧縮処理後
の受信信号列ＳPm(ｎ,ｔ,Ｒ)は映像信号処理器２２２お
よび尤度比演算器２２３に出力されているものである。
その映像信号処理器２２２では、距離圧縮器２２１から
の受信信号列ＳPm(ｎ,ｔ,Ｒ)は映像信号に変換された
上、表示装置１１５上で画像として表示されているもの
である。一方、尤度比演算器２２３では、受信信号列Ｓ
m(ｎ,ｔ)と受信信号列ＳPm(ｎ,ｔ,Ｒ)とに基づく処理が
所定に行われているが、その処理を図２４を参照しつ
つ、説明すれば以下のようである。
１１４での処理について説明すれば、受信装置１１３か
らの受信信号列Ｓm(ｎ,ｔ)は距離圧縮器２２１および尤
度比演算器２２３へ出力されているが、このうち、距離
圧縮器２２１では、受信装置１１３からの受信信号列Ｓ
m(ｎ,ｔ)はＲ個の振幅値を含む区間毎に領域分割された
上、それら区間各々における最大振幅値が検出される、
といった距離圧縮処理が行われており、距離圧縮処理後
の受信信号列ＳPm(ｎ,ｔ,Ｒ)は映像信号処理器２２２お
よび尤度比演算器２２３に出力されているものである。
その映像信号処理器２２２では、距離圧縮器２２１から
の受信信号列ＳPm(ｎ,ｔ,Ｒ)は映像信号に変換された
上、表示装置１１５上で画像として表示されているもの
である。一方、尤度比演算器２２３では、受信信号列Ｓ
m(ｎ,ｔ)と受信信号列ＳPm(ｎ,ｔ,Ｒ)とに基づく処理が
所定に行われているが、その処理を図２４を参照しつ
つ、説明すれば以下のようである。
【００８３】即ち、その尤度比演算器２２３自体は、領
域設定器２５１、振幅値検出器２５２、分散算出器２５
３および尤度比算出器２５４から構成されているが、図
８に示す尤度比演算器９０に比し、平均値算出器９３が
振幅値検出器２５２に置換されており、また、尤度比算
出器２５４は図１に示す尤度比演算器６０に比し、計算
器６１ｉに設定される尤度比算出式が異なるだけで、こ
れら以外の構成は同様となっている。その尤度比演算器
２２３での処理であるが、先ず領域設定器２５１では、
受信装置１１３からの受信信号列Ｓm(ｎ,ｔ)に対し、時
間（距離）および方位方向に一定幅をもつ目標捜索領域
および分散計算領域が設定されているが、その際、目標
捜索領域は、距離圧縮器２２１で設定されている領域と
同一として領域設定器２５１で設定されるものとなって
いる。受信信号列Ｓm(ｎ,ｔ)の全領域がカバーされるべ
く、目標捜索領域および分散計算領域は時間（距離）お
よび方位方向に順次シフト設定されているわけである
が、そのシフト設定の度に、目標捜索領域間境界信号が
振幅値保持器２５２へ、分散計算領域に含まれる受信信
号振幅値が分散算出器２５３に出力されているものであ
るが、この分散算出器２５３の機能は尤度比演算器７０
における分散算出器７２と同様である。一方、振幅値保
持器２５２では、距離圧縮器２２１からの、距離圧縮処
理後の受信信号列ＳPm(ｎ,ｔ,Ｒ)が目標捜索領域間境界
信号により、所望の目標捜索領域に相当する領域の最大
振幅値ｘPHとして一時的に保持された上、尤度比算出器
２５４へ出力され、また、尤度比算出器２５４では、振
幅値保持器２５２からの最大振幅値ｘPHと分散算出器２
５３からの受信信号振幅値の分散値σ2 とから尤度比Λ
m(ｘ(nt))が算出されているものである。
域設定器２５１、振幅値検出器２５２、分散算出器２５
３および尤度比算出器２５４から構成されているが、図
８に示す尤度比演算器９０に比し、平均値算出器９３が
振幅値検出器２５２に置換されており、また、尤度比算
出器２５４は図１に示す尤度比演算器６０に比し、計算
器６１ｉに設定される尤度比算出式が異なるだけで、こ
れら以外の構成は同様となっている。その尤度比演算器
２２３での処理であるが、先ず領域設定器２５１では、
受信装置１１３からの受信信号列Ｓm(ｎ,ｔ)に対し、時
間（距離）および方位方向に一定幅をもつ目標捜索領域
および分散計算領域が設定されているが、その際、目標
捜索領域は、距離圧縮器２２１で設定されている領域と
同一として領域設定器２５１で設定されるものとなって
いる。受信信号列Ｓm(ｎ,ｔ)の全領域がカバーされるべ
く、目標捜索領域および分散計算領域は時間（距離）お
よび方位方向に順次シフト設定されているわけである
が、そのシフト設定の度に、目標捜索領域間境界信号が
振幅値保持器２５２へ、分散計算領域に含まれる受信信
号振幅値が分散算出器２５３に出力されているものであ
るが、この分散算出器２５３の機能は尤度比演算器７０
における分散算出器７２と同様である。一方、振幅値保
持器２５２では、距離圧縮器２２１からの、距離圧縮処
理後の受信信号列ＳPm(ｎ,ｔ,Ｒ)が目標捜索領域間境界
信号により、所望の目標捜索領域に相当する領域の最大
振幅値ｘPHとして一時的に保持された上、尤度比算出器
２５４へ出力され、また、尤度比算出器２５４では、振
幅値保持器２５２からの最大振幅値ｘPHと分散算出器２
５３からの受信信号振幅値の分散値σ2 とから尤度比Λ
m(ｘ(nt))が算出されているものである。
【００８４】ここで、レイリー分布およびライス分布に
従う確率変数の一定区間の最大値が従う確率密度関数の
導出過程および尤度比算出式の設定方法を以下に示す。
即ち、一般に、確率変数ｘが確率密度関数ｐ(ｘ)に従う
場合、Ｒ個の出力の最大値ｘPHが従う確率密度関数ｐPH
(ｘPH)は以下の数式２０で表される。
従う確率変数の一定区間の最大値が従う確率密度関数の
導出過程および尤度比算出式の設定方法を以下に示す。
即ち、一般に、確率変数ｘが確率密度関数ｐ(ｘ)に従う
場合、Ｒ個の出力の最大値ｘPHが従う確率密度関数ｐPH
(ｘPH)は以下の数式２０で表される。
【００８５】
【数２０】
【００８６】従って、ｐ(ｘ)にレイリー分布の確率密度
関数（数式６）を代入することにより、目標信号が含ま
れない受信信号列の、一定区間内の最大振幅値の分布が
従うべき確率密度関数ｐ0(ｘPH)は以下の数式２１とし
て得られる。
関数（数式６）を代入することにより、目標信号が含ま
れない受信信号列の、一定区間内の最大振幅値の分布が
従うべき確率密度関数ｐ0(ｘPH)は以下の数式２１とし
て得られる。
【００８７】
【数２１】
【００８８】同様に、ｐ(ｘ)にライス分布の確率密度関
数（数式７）を代入することにより、目標信号を含む受
信信号列の、一定区間内の最大振幅値の分布が従うべき
確率密度関数ｐ1(ｘPH)は以下の数式２２として得られ
る。
数（数式７）を代入することにより、目標信号を含む受
信信号列の、一定区間内の最大振幅値の分布が従うべき
確率密度関数ｐ1(ｘPH)は以下の数式２２として得られ
る。
【００８９】
【数２２】
【００９０】故に、計算器６１ｉにおける尤度比算出式
は以下の数式２３として定義される。
は以下の数式２３として定義される。
【００９１】
【数２３】
【００９２】その数式２３より、尤度比Λ(ｘ)の値を決
定する４つのパラメータ、即ち、目標捜索領域に含まれ
る受信信号振幅値の最大値ｘPH、分散計算領域に含まれ
る受信信号振幅値の分散値σ2 、信号振幅値Ａ、および
一定区間内（目標捜索領域内）に含まれる受信信号振幅
値の数Ｒのうち、信号振幅値Ａおよび振幅値の数Ｒは事
前設定される定数であることから、結局、最大値ｘPHお
よび分散値σ2 からその目標捜索領域における尤度比Λ
(ｘ)が求まり、数式５により予め設定されている閾値η
1 ，η0 各々との比較により目標信号の検出が可能とさ
れているものである。
定する４つのパラメータ、即ち、目標捜索領域に含まれ
る受信信号振幅値の最大値ｘPH、分散計算領域に含まれ
る受信信号振幅値の分散値σ2 、信号振幅値Ａ、および
一定区間内（目標捜索領域内）に含まれる受信信号振幅
値の数Ｒのうち、信号振幅値Ａおよび振幅値の数Ｒは事
前設定される定数であることから、結局、最大値ｘPHお
よび分散値σ2 からその目標捜索領域における尤度比Λ
(ｘ)が求まり、数式５により予め設定されている閾値η
1 ，η0 各々との比較により目標信号の検出が可能とさ
れているものである。
【００９３】因みに、図２５（ａ）に一定区間内の最大
振幅値が従うべき確率密度関数の分布形状と最大値ｘPH
との関係を示す。図中、レイリー分布と併記されている
分布形状は、数式２０で定義されている確率密度関数
に、また、ライス分布と併記されている分布形状（Ａ＝
２.４５、Ａ＝４.２４の場合）は、数式２２で定義され
ている確率密度関数にそれぞれ基づいている。図２５
（ｂ）にはまた、数式２３により求まる尤度比Λ(ｘ)と
ＳＮ比の関係を示す。これより、尤度比算出式として数
式２３が用いられた場合でも、尤度比を用いて検出可能
な信号の最小ＳＮ比（正の値の尤度比Λ(ｘ)を算出し得
るＳＮ比の最小値）は、パラメータとしての信号振幅値
Ａによって変化するものであることが判る。従って、信
号振幅値Ａが異なる尤度比算出式が複数設定される場合
は、目標物体が存在すると判定された領域の平均的なＳ
Ｎ比が算出可能となる。
振幅値が従うべき確率密度関数の分布形状と最大値ｘPH
との関係を示す。図中、レイリー分布と併記されている
分布形状は、数式２０で定義されている確率密度関数
に、また、ライス分布と併記されている分布形状（Ａ＝
２.４５、Ａ＝４.２４の場合）は、数式２２で定義され
ている確率密度関数にそれぞれ基づいている。図２５
（ｂ）にはまた、数式２３により求まる尤度比Λ(ｘ)と
ＳＮ比の関係を示す。これより、尤度比算出式として数
式２３が用いられた場合でも、尤度比を用いて検出可能
な信号の最小ＳＮ比（正の値の尤度比Λ(ｘ)を算出し得
るＳＮ比の最小値）は、パラメータとしての信号振幅値
Ａによって変化するものであることが判る。従って、信
号振幅値Ａが異なる尤度比算出式が複数設定される場合
は、目標物体が存在すると判定された領域の平均的なＳ
Ｎ比が算出可能となる。
【００９４】以上、説明したように、本発明に係る各実
施の形態での目標検出装置では、尤度比を用いた目標検
出法が適用されたことにより、受信信号列中の目標信号
が検出され得るばかりか、検出された目標信号はそのＳ
Ｎ比毎に分類されていることから、ユーザに目標の位置
と同時に、その目標のＳＮ比をも容易に認識させ得るも
のとなっている。また、目標信号、あるいは受信信号列
の特徴に応じて、目標捜索領域および尤度比算出式が変
更される場合には、目標信号の検出性能は高性能に維持
されるものとなっている。更に、以上の説明では、目標
検出装置は船舶に搭載された状態で、水中に超音波を送
波し、反射超音波を信号処理して目標信号が検出される
構成とされているが、これに限定されることなく、例え
ば超音波の代わりに電波を用い、大気中に存在している
目標物体を検出する場合にも適用可とされたものとなっ
ている。
施の形態での目標検出装置では、尤度比を用いた目標検
出法が適用されたことにより、受信信号列中の目標信号
が検出され得るばかりか、検出された目標信号はそのＳ
Ｎ比毎に分類されていることから、ユーザに目標の位置
と同時に、その目標のＳＮ比をも容易に認識させ得るも
のとなっている。また、目標信号、あるいは受信信号列
の特徴に応じて、目標捜索領域および尤度比算出式が変
更される場合には、目標信号の検出性能は高性能に維持
されるものとなっている。更に、以上の説明では、目標
検出装置は船舶に搭載された状態で、水中に超音波を送
波し、反射超音波を信号処理して目標信号が検出される
構成とされているが、これに限定されることなく、例え
ば超音波の代わりに電波を用い、大気中に存在している
目標物体を検出する場合にも適用可とされたものとなっ
ている。
【００９５】
【発明の効果】以上、説明したように、請求項１〜７に
よれば、時間的に、または周波数的に振幅値が変動して
いる複数の受信信号列から、各受信信号列中の同一位置
に定常的に存在している物体を目標として検出するに際
して、その目標の位置に併せてそのＳＮ比が容易に認識
され得るものとなっている。
よれば、時間的に、または周波数的に振幅値が変動して
いる複数の受信信号列から、各受信信号列中の同一位置
に定常的に存在している物体を目標として検出するに際
して、その目標の位置に併せてそのＳＮ比が容易に認識
され得るものとなっている。
【図１】図１は、本発明による目標検出装置の一例での
要部構成を示す図
要部構成を示す図
【図２】図２は、レイリー分布形状、ライス分布形状お
よびこれら両者の比としての尤度比の変化を示す図
よびこれら両者の比としての尤度比の変化を示す図
【図３】図３は、尤度比とＳＮ比の関係を示す図
【図４】図４は、ＳＮ比と確率密度関数の関係を示す図
【図５】図５は、パラメータ（信号振幅値）Ａと尤度比
の関係を示す図
の関係を示す図
【図６】図６は、尤度比による目標検出処理が行われて
いる、尤度比算出器を含む尤度比演算器の一例での概要
構成を示す図
いる、尤度比算出器を含む尤度比演算器の一例での概要
構成を示す図
【図７】図７は、時間方向に変動する受信信号列への、
目標捜索領域および分散計算領域の設定態様を示す図
目標捜索領域および分散計算領域の設定態様を示す図
【図８】図８は、尤度比演算器の他の例での概要構成を
示す図
示す図
【図９】図９（ａ）は、カイ２乗確率密度関数の分布形
状を、また、図９（ｂ）は、ＳＮ比（ＳＮＲ）と尤度比
の関係を示す図
状を、また、図９（ｂ）は、ＳＮ比（ＳＮＲ）と尤度比
の関係を示す図
【図１０】図１０は、本発明による目標検出装置の一例
での全体構成を示す図
での全体構成を示す図
【図１１】図１１（ａ），（ｂ）は、それぞれディジタ
ル変換前の受信信号、ディジタル変換後の受信信号を示
す図
ル変換前の受信信号、ディジタル変換後の受信信号を示
す図
【図１２】図１２は、整相処理による複数ビーム形成を
説明するための図
説明するための図
【図１３】図１３は、送信信号波形に対する、方位方向
対応受信信号波形との関係を示す図
対応受信信号波形との関係を示す図
【図１４】図１４は、その方位方向対応受信信号波形を
２次元状に配列した状態として示す図
２次元状に配列した状態として示す図
【図１５】図１５は、本発明による目標検出装置におけ
る映像信号・尤度比演算処理に係る一例での構成を示す
図
る映像信号・尤度比演算処理に係る一例での構成を示す
図
【図１６】図１６は、２次元状に配列された受信信号列
への、目標捜索領域および分散計算領域の設定態様を示
す図
への、目標捜索領域および分散計算領域の設定態様を示
す図
【図１７】図１７は、受信信号間で座標変換が行われる
際でのモデルを示す図
際でのモデルを示す図
【図１８】図１８は、受信信号列の画像としての表示例
を示す図
を示す図
【図１９】図１９は、本発明による目標検出装置の他の
例での要部としての、映像信号・尤度比演算処理に係る
構成を示す図
例での要部としての、映像信号・尤度比演算処理に係る
構成を示す図
【図２０】図２０（ａ），（ｂ）は、目標信号の時間幅
がサンプリング周期より長い場合での、ディジタル変換
前の受信信号、ディジタル変換後の受信信号をそれぞれ
示す図
がサンプリング周期より長い場合での、ディジタル変換
前の受信信号、ディジタル変換後の受信信号をそれぞれ
示す図
【図２１】図２１は、本発明による目標検出装置の更に
異なる他の例での要部としての、映像信号・尤度比演算
処理に係る構成を示す図
異なる他の例での要部としての、映像信号・尤度比演算
処理に係る構成を示す図
【図２２】図２２（ａ），（ｂ）は、ディジタル変換前
の受信信号がサンプリング周期の相違によって、如何に
ディジタル変換されるかを示す図
の受信信号がサンプリング周期の相違によって、如何に
ディジタル変換されるかを示す図
【図２３】図２３（ａ），（ｂ）は、それぞれ距離圧縮
処理前の受信信号、距離圧縮処理後の受信信号をそれぞ
れ示す図
処理前の受信信号、距離圧縮処理後の受信信号をそれぞ
れ示す図
【図２４】図２４は、その距離圧縮処理後の受信信号が
考慮された、尤度比演算処理に係る構成を示す図
考慮された、尤度比演算処理に係る構成を示す図
【図２５】図２５（ａ）は、一定区間内の最大振幅値が
従うべき確率密度関数の分布形状を、図２５（ｂ）は、
尤度比とＳＮ比の関係をそれぞれ示す図
従うべき確率密度関数の分布形状を、図２５（ｂ）は、
尤度比とＳＮ比の関係をそれぞれ示す図
【図２６】図２６は、積算平均化法の原理を説明するた
めの図
めの図
６０，９４，２５４…尤度比算出器、６１，６１ｉ…計
算器、６２，６２ｉ…乗算器、６３，６３ｉ…判定器、
６４…比較器、７０，９０，２２３…尤度比演算器、７
１，９１，２５１…領域設定器、９２，２５３…分散算
出器、９３…平均値算出器、１１１…送信装置、１１２
…送受波装置、１１３…受信装置、１１４…信号処理装
置、１１５…表示装置、１１６…速度測定装置、１１７
…水深測定装置、１６２、２２２…映像信号処理器、２
２１…距離圧縮器、２５２…振幅値保持器
算器、６２，６２ｉ…乗算器、６３，６３ｉ…判定器、
６４…比較器、７０，９０，２２３…尤度比演算器、７
１，９１，２５１…領域設定器、９２，２５３…分散算
出器、９３…平均値算出器、１１１…送信装置、１１２
…送受波装置、１１３…受信装置、１１４…信号処理装
置、１１５…表示装置、１１６…速度測定装置、１１７
…水深測定装置、１６２、２２２…映像信号処理器、２
２１…距離圧縮器、２５２…振幅値保持器
───────────────────────────────────────────────────── フロントページの続き (72)発明者 佐々木 潤 神奈川県横浜市戸塚区戸塚町216番地 株 式会社日立製作所情報通信事業部内 (72)発明者 藤井 正司 神奈川県横浜市戸塚区戸塚町216番地 株 式会社日立製作所情報通信事業部内
Claims (7)
- 【請求項１】 ソーナー、レーダー等、音波、または電
波が複数回に亘って送信信号として送信される度に、該
送信信号各々に対する応答としての反射信号を受信信号
列として受信した上、該受信信号列各々は映像信号に変
換された状態で画像として表示されつつ、上記受信信号
列各々からは、同一位置に定常的に存在している目標か
らの反射信号が目標信号として検出された上、該目標の
位置が上記画像上に併せて表示されるようにした目標検
出装置であって、目標の位置と該目標のＳＮ比が画像上
に同時に表示されるべく、受信信号列各々から該目標の
位置と該目標のＳＮ比を特定する処理を行う目標位置・
ＳＮ比特定手段を少なくとも含む構成の目標検出装置。 - 【請求項２】 ソーナー、レーダー等、音波、または電
波が複数回に亘って送信信号として送信される度に、該
送信信号各々に対する応答としての反射信号を受信信号
列として受信した上、該受信信号列各々は映像信号に変
換された状態で画像として表示されつつ、上記受信信号
列各々からは、同一位置に定常的に存在している目標か
らの反射信号が目標信号として検出された上、該目標の
位置が上記画像上に併せて表示されるようにした目標検
出装置であって、目標の位置と該目標のＳＮ比が画像上
に同時に表示されるべく、受信信号列各々から該目標の
位置と該目標のＳＮ比を尤度比により特定する処理を行
う目標位置・ＳＮ比特定手段を少なくとも含む構成の目
標検出装置。 - 【請求項３】 ソーナー、レーダー等、音波、または電
波が複数回に亘って送信信号として送信される度に、該
送信信号各々に対する応答としての反射信号を受信信号
列として受信した上、該受信信号列各々は映像信号に変
換された状態で画像として表示されつつ、上記受信信号
列各々からは、同一位置に定常的に存在している目標か
らの反射信号が目標信号として検出された上、該目標の
位置が上記画像上に併せて表示されるようにした目標検
出装置であって、目標の位置と該目標のＳＮ比が画像上
に同時に表示されるべく、受信信号列各々から該目標の
位置と該目標のＳＮ比を尤度比により特定する処理を行
う目標位置・ＳＮ比特定手段を少なくとも含む場合に、
該目標位置・ＳＮ比特定手段は、受信信号列が得られる
度に、該受信信号列中に任意の目標捜索領域および該目
標捜索領域に対応する分散計算領域を順次更新設定する
手段と、該目標捜索領域から該領域内に含まれる受信信
号振幅値を検出する手段と、該分散計算領域から該領域
内に含まれる受信信号振幅値の分散値を算出する手段
と、目標信号を含まない受信信号列の分布が従うべきレ
イリー分布の確率密度関数と目標信号を含む受信信号列
の分布が従うべきライス分布の、目標信号のＳＮ比に応
じた複数の確率密度関数各々との比として定義される尤
度比算出式が複数設定された状態で、上記受信信号列が
得られる度に、該目標捜索領域各々から検出された受信
信号振幅値と上記分散計算領域各々から算出された分散
値とにもとづき、上記尤度比算出式各々から目標捜索領
域毎に尤度比を算出する手段と、目標捜索領域毎に算出
された尤度比各々の事前設定閾値との比較によって、該
閾値を上回る累積尤度比に係る目標捜索領域に目標が存
在すると判定する手段と、目標が存在すると判定された
目標捜索領域を該目標の位置として、用いられた尤度比
算出式毎に分類する手段と、目標が存在すると判定され
た目標捜索領域を該目標の位置として、該位置を分類結
果とともに画像上に表示せしめる手段とから構成されて
いる目標検出装置。 - 【請求項４】 ソーナー、レーダー等、音波、または電
波が複数回に亘って送信信号として送信される度に、該
送信信号各々に対する応答としての反射信号を受信信号
列として受信した上、該受信信号列各々は映像信号に変
換された状態で画像として表示されつつ、上記受信信号
列各々からは、同一位置に定常的に存在している目標か
らの反射信号が目標信号として検出された上、該目標の
位置が上記画像上に併せて表示されるようにした目標検
出装置であって、目標の位置と該目標のＳＮ比が画像上
に同時に表示されるべく、受信信号列各々から該目標の
位置と該目標のＳＮ比を尤度比により特定する処理を行
う目標位置・ＳＮ比特定手段を少なくとも含む場合に、
該目標位置・ＳＮ比特定手段は、受信信号列が得られる
度に、該受信信号列中に任意の目標捜索領域および該目
標捜索領域に対応する分散計算領域を順次更新設定する
手段と、該目標捜索領域から該領域内に含まれる受信信
号振幅値の２乗平均値を算出する手段と、該分散計算領
域から該領域内に含まれる受信信号振幅値の分散値を算
出する手段と、目標信号を含まない受信信号列の２乗平
均値の分布が従うべきカイ２乗分布の確率密度関数と目
標信号を含む受信信号列の２乗平均値の分布が従うべき
非心カイ２乗分布の、目標信号のＳＮ比に応じた複数の
確率密度関数各々との比として定義される尤度比算出式
が複数設定された状態で、上記受信信号列が得られる度
に、該目標捜索領域各々から検出された受信信号振幅値
の２乗平均値と上記分散計算領域各々から算出された分
散値とにもとづき、上記尤度比算出式各々から目標捜索
領域毎に尤度比を算出する手段と、目標捜索領域毎に算
出された尤度比各々の事前設定閾値との比較によって、
該閾値を上回る累積尤度比に係る目標捜索領域に目標が
存在すると判定する手段と、目標が存在すると判定され
た目標捜索領域を該目標の位置として、用いられた尤度
比算出式毎に分類する手段と、目標が存在すると判定さ
れた目標捜索領域を該目標の位置として、該位置を分類
結果とともに画像上に表示せしめる手段とから構成され
ている目標検出装置。 - 【請求項５】 ソーナー、レーダー等、音波、または電
波が複数回に亘って送信信号として送信される度に、該
送信信号各々に対する応答としての反射信号を受信信号
列として受信した上、該受信信号列各々は映像信号に変
換された状態で画像として表示されつつ、上記受信信号
列各々からは、尤度比を用いた処理により同一位置に定
常的に存在している目標からの反射信号が目標信号とし
て検出された上、該目標の位置が上記画像上に併せて表
示されるようにした目標検出装置であって、目標のＳＮ
比を算出する手段と、該ＳＮ比の算出結果を該目標の位
置とともに画像上に表示せしめる手段とを少なくとも含
む構成の目標検出装置。 - 【請求項６】 ソーナー、レーダー等、音波、または電
波が複数回に亘って送信信号として送信される度に、該
送信信号各々に対する応答としての反射信号を受信信号
列として受信した上、該受信信号列各々は映像信号に変
換された状態で画像として表示されつつ、上記受信信号
列各々に更新設定された目標捜索領域及び分散計算領域
からは、尤度比を用いた処理により同一位置に定常的に
存在している目標からの反射信号が目標信号として検出
された上、該目標の位置が上記画像上に併せて表示され
るようにした目標検出装置であって、目標のＳＮ比を算
出するＳＮ比算出手段と、該ＳＮ比の算出結果を該目標
の位置とともに画像上に表示せしめる手段とを少なくと
も含む場合に、ＳＮ比算出手段は、目標信号を含む受信
信号列の分布が従うべきライス分布の、目標信号のＳＮ
比に応じた確率密度関数として定義される尤度比算出式
が複数設定された状態で、受信信号列が得られる度に、
目標捜索領域各々から検出された受信信号振幅値と上記
分散計算領域各々から算出された分散値とにもとづき、
上記尤度比算出式各々から目標捜索領域毎に尤度比を算
出する手段と、目標捜索領域毎に算出された尤度比各々
の事前設定閾値との比較によって、該閾値を上回る累積
尤度比に係る目標捜索領域に目標が存在すると判定する
手段と、目標が存在すると判定された目標捜索領域を該
目標の位置として、用いられた尤度比算出式毎に分類す
る手段と、目標が存在すると判定された目標捜索領域を
該目標の位置として、該位置を分類結果とともに画像上
に表示せしめる手段とから構成されている目標検出装
置。 - 【請求項７】 ソーナー、レーダー等、音波、または電
波が複数回に亘って送信信号として送信される度に、該
送信信号各々に対する応答としての反射信号を受信信号
列として受信した上、該受信信号列各々は映像信号に変
換された状態で画像として表示されつつ、上記受信信号
列各々に更新設定された目標捜索領域及び分散計算領域
からは、尤度比を用いた処理により同一位置に定常的に
存在している目標からの反射信号が目標信号として検出
された上、該目標の位置が上記画像上に併せて表示され
るようにした目標検出装置であって、目標のＳＮ比を算
出するＳＮ比算出手段と、該ＳＮ比の算出結果を該目標
の位置とともに画像上に表示せしめる手段とを少なくと
も含む場合に、ＳＮ比算出手段は、目標信号を含む受信
信号列の２乗平均値の分布が従うべき非心カイ２乗分布
分布の、目標信号のＳＮ比に応じた確率密度関数として
定義される尤度比算出式が複数設定された状態で、受信
信号列が得られる度に、目標捜索領域各々から算出され
た受信信号振幅値の２乗平均値と上記分散計算領域各々
から算出された分散値とにもとづき、上記尤度比算出式
各々から目標捜索領域毎に尤度比を算出する手段と、目
標捜索領域毎に算出された尤度比各々の事前設定閾値と
の比較によって、該閾値を上回る累積尤度比に係る目標
捜索領域に目標が存在すると判定する手段と、目標が存
在すると判定された目標捜索領域を該目標の位置とし
て、用いられた尤度比算出式毎に分類する手段と、目標
が存在すると判定された目標捜索領域を該目標の位置と
して、該位置を分類結果とともに画像上に表示せしめる
手段とから構成されている目標検出装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP01160897A JP3785715B2 (ja) | 1997-01-24 | 1997-01-24 | 目標検出装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP01160897A JP3785715B2 (ja) | 1997-01-24 | 1997-01-24 | 目標検出装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH10206523A true JPH10206523A (ja) | 1998-08-07 |
| JP3785715B2 JP3785715B2 (ja) | 2006-06-14 |
Family
ID=11782632
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP01160897A Expired - Lifetime JP3785715B2 (ja) | 1997-01-24 | 1997-01-24 | 目標検出装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3785715B2 (ja) |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006194887A (ja) * | 2005-01-13 | 2006-07-27 | Ma Com Inc | 車両用センサシステム及び方法 |
| JP2006343252A (ja) * | 2005-06-10 | 2006-12-21 | Nec Corp | 信号検出方法及び信号検出装置 |
| JP2011209238A (ja) * | 2010-03-30 | 2011-10-20 | Fujitsu Ltd | 検知装置および検知プログラム |
| JP2012093286A (ja) * | 2010-10-28 | 2012-05-17 | Mitsubishi Electric Corp | 目標検出装置および目標検出方法 |
| JP2012163350A (ja) * | 2011-02-03 | 2012-08-30 | Toshiba Corp | レーダ装置、レーダ受信装置及び目標検出方法 |
| JP2017078694A (ja) * | 2015-10-22 | 2017-04-27 | 株式会社東芝 | 目標検出装置、目標検出方法及び目標検出プログラム |
| JP2017078693A (ja) * | 2015-10-22 | 2017-04-27 | 株式会社東芝 | 信号検出装置、信号検出方法及び信号検出プログラム |
| JP2017156148A (ja) * | 2016-02-29 | 2017-09-07 | 株式会社東芝 | 信号検出装置、信号検出方法及び信号検出プログラム |
| JP2017156147A (ja) * | 2016-02-29 | 2017-09-07 | 株式会社東芝 | 目標検出装置、目標検出方法及び目標検出プログラム |
| JP2017156146A (ja) * | 2016-02-29 | 2017-09-07 | 株式会社東芝 | 目標検出装置、目標検出方法及び目標検出プログラム |
-
1997
- 1997-01-24 JP JP01160897A patent/JP3785715B2/ja not_active Expired - Lifetime
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006194887A (ja) * | 2005-01-13 | 2006-07-27 | Ma Com Inc | 車両用センサシステム及び方法 |
| JP2006343252A (ja) * | 2005-06-10 | 2006-12-21 | Nec Corp | 信号検出方法及び信号検出装置 |
| JP2011209238A (ja) * | 2010-03-30 | 2011-10-20 | Fujitsu Ltd | 検知装置および検知プログラム |
| JP2012093286A (ja) * | 2010-10-28 | 2012-05-17 | Mitsubishi Electric Corp | 目標検出装置および目標検出方法 |
| JP2012163350A (ja) * | 2011-02-03 | 2012-08-30 | Toshiba Corp | レーダ装置、レーダ受信装置及び目標検出方法 |
| US8810446B2 (en) | 2011-02-03 | 2014-08-19 | Kabushiki Kaisha Toshiba | Radar device, radar receiver, and target detection method |
| JP2017078694A (ja) * | 2015-10-22 | 2017-04-27 | 株式会社東芝 | 目標検出装置、目標検出方法及び目標検出プログラム |
| JP2017078693A (ja) * | 2015-10-22 | 2017-04-27 | 株式会社東芝 | 信号検出装置、信号検出方法及び信号検出プログラム |
| JP2017156148A (ja) * | 2016-02-29 | 2017-09-07 | 株式会社東芝 | 信号検出装置、信号検出方法及び信号検出プログラム |
| JP2017156147A (ja) * | 2016-02-29 | 2017-09-07 | 株式会社東芝 | 目標検出装置、目標検出方法及び目標検出プログラム |
| JP2017156146A (ja) * | 2016-02-29 | 2017-09-07 | 株式会社東芝 | 目標検出装置、目標検出方法及び目標検出プログラム |
Also Published As
| Publication number | Publication date |
|---|---|
| JP3785715B2 (ja) | 2006-06-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7852709B1 (en) | Sonar system and process | |
| JP5398195B2 (ja) | レーダ装置 | |
| RU2653189C1 (ru) | Способ обнаружения шумящих объектов в мелком и глубоком море | |
| JP2014002085A (ja) | 信号処理装置、レーダ装置、信号処理方法、およびプログラム | |
| CN113884986B (zh) | 波束聚焦增强的强冲击信号空时域联合检测方法及系统 | |
| JP2016206153A (ja) | 信号処理装置、及びレーダ装置 | |
| JPH10206523A (ja) | 目標検出装置 | |
| EP4307003A1 (en) | Precipitation clutter suppression algorithm for marine radar | |
| RU2711406C1 (ru) | Способ классификации гидроакустических сигналов шумоизлучения морских объектов | |
| JP2014089056A (ja) | 追尾処理装置、及び追尾処理方法 | |
| JP2000230972A (ja) | レーダ信号処理装置 | |
| EP0393741B1 (en) | Acoustic detection device | |
| JP7337575B2 (ja) | レーダ監視装置及び方法 | |
| JP6347660B2 (ja) | 魚群探知機、単体魚検出方法、および単体魚検出プログラム | |
| CN110703261A (zh) | 一种基于多峰值评分制的侧扫声呐底跟踪方法 | |
| RU2654335C1 (ru) | Способ обнаружения шумящих в море объектов с помощью комбинированного приемника | |
| CN111965629B (zh) | 一种主动声呐非均匀背景抑制恒虚警检测方法及装置 | |
| JP3871875B2 (ja) | 目標類別方法及び装置 | |
| JP2000111646A (ja) | ドップラ検出表示方法及びシステム | |
| De Roos et al. | The detection and classification of objects lying on the seafloor | |
| JPH10206524A (ja) | 目標検出装置 | |
| KR20170082080A (ko) | 수중 고정표적과 클러터 식별을 위한 가변 데이터 맵 및 가변 통계적 특징정보를 적용하는 방법 및 장치 | |
| RU2305297C2 (ru) | Устройство с корреляционным формирователем характеристики направленности для обнаружения сигналов и определения направления на их источник | |
| RU2736188C9 (ru) | Способ отображения гидроакустической информации | |
| JPH10253758A (ja) | 周波数推定方法、周波数推定装置、ドップラソナーおよび潮流計 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20050802 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20051122 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20060123 |
|
| 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: 20060228 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20060313 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090331 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100331 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110331 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110331 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120331 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130331 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130331 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140331 Year of fee payment: 8 |
|
| EXPY | Cancellation because of completion of term |