JP2607861B2 - 信号形式判定方法 - Google Patents
信号形式判定方法Info
- Publication number
- JP2607861B2 JP2607861B2 JP7024928A JP2492895A JP2607861B2 JP 2607861 B2 JP2607861 B2 JP 2607861B2 JP 7024928 A JP7024928 A JP 7024928A JP 2492895 A JP2492895 A JP 2492895A JP 2607861 B2 JP2607861 B2 JP 2607861B2
- Authority
- JP
- Japan
- Prior art keywords
- frequency
- pulse
- signal
- residual
- peak
- 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 - Lifetime
Links
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Radar Systems Or Details Thereof (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Description
【０００１】
【産業上の利用分野】本発明は、例えば、ソナー等にお
いて、入力信号から形成したパルス信号の変調形式を判
定する信号形式判定方法に関するものである。
いて、入力信号から形成したパルス信号の変調形式を判
定する信号形式判定方法に関するものである。
【０００２】
【従来の技術】従来、このような分野の技術としては、
例えば次のようなものがあった。ソナー等が目標から放
射されるパルス音等を検知し、パルスの変調形式を推定
する場合、入力信号を一定の時間間隔で周波数分析して
強度が極大となるピーク周波数を測定し、このピークの
時間的連なりとしてパルス信号を形成して、このパルス
信号内での周波数の変化の仕方からパルス信号の変調形
式を判定するという方法が取られている。図２は、従来
の信号形式判定方法を実施するための信号形式判定装置
の機能ブロック図である。この図を参照しつつ従来の信
号形式判定方法を説明する。この信号形式判定装置は、
目標から放射されるパルス音を周囲雑音と共に入力する
入力端子１１を有している。入力端子１１は、周波数分
析器１２の入力側に接続されている。周波数分析器１２
は、入力端子１１から入力された信号を一定の時間間隔
で周波数分析する手段である。周波数分析器１２の出力
側は、周波数分析器１２で周波数分析した結果を各周波
数ビン毎に短時間の時定数で積分する短時間積分器１３
の入力側に接続されると共に、周波数分析器１２で周波
数分析した結果を各周波数ビン毎に長時間の時定数で積
分する長時間積分器１４の入力側に接続されている。
例えば次のようなものがあった。ソナー等が目標から放
射されるパルス音等を検知し、パルスの変調形式を推定
する場合、入力信号を一定の時間間隔で周波数分析して
強度が極大となるピーク周波数を測定し、このピークの
時間的連なりとしてパルス信号を形成して、このパルス
信号内での周波数の変化の仕方からパルス信号の変調形
式を判定するという方法が取られている。図２は、従来
の信号形式判定方法を実施するための信号形式判定装置
の機能ブロック図である。この図を参照しつつ従来の信
号形式判定方法を説明する。この信号形式判定装置は、
目標から放射されるパルス音を周囲雑音と共に入力する
入力端子１１を有している。入力端子１１は、周波数分
析器１２の入力側に接続されている。周波数分析器１２
は、入力端子１１から入力された信号を一定の時間間隔
で周波数分析する手段である。周波数分析器１２の出力
側は、周波数分析器１２で周波数分析した結果を各周波
数ビン毎に短時間の時定数で積分する短時間積分器１３
の入力側に接続されると共に、周波数分析器１２で周波
数分析した結果を各周波数ビン毎に長時間の時定数で積
分する長時間積分器１４の入力側に接続されている。
【０００３】短時間積分器１３及び長時間積分器１４の
各出力側は、ピーク検出器１５の入力側に接続されてい
る。ピーク検出器１５は、短時間積分器１３から送られ
た短時間積分結果及び長時間積分器１４から送られた長
時間積分結果に基づいて強度が極大となるピークを検出
する手段である。ピーク検出器１５の出力側は、入力端
子１１から入力する入力信号を分析してピーク検出器１
５で検出されたピークに対応する正確な周波数を測定す
る周波数精測器１６に接続されている。又、ピーク検出
器１５の出力側は、ピーク検出器１５から送られたピー
ク検出情報及び周波数精測器１６から送られた周波数情
報により、同一パルスの成分と判断されるピークの時間
的連なりからパルスを形成するパルス形成器１７の入力
側に接続されている。パルス形成器１７の出力側は、パ
ルス形成器１７で形成されたパルス信号の周波数変化の
仕方を線形回帰により推定する勾配推定器１８の入力側
に接続されている。勾配推定器１８の出力側は、勾配推
定器１８で推定されたパルス信号の周波数変化の勾配値
と、予め設定された閾値とからパルス信号の変調形式を
判定するパルス判定処理器１９の入力側に接続されてい
る。パルス判定処理器１９の出力側は、パルス信号の各
種パラメータの出力端子２０に接続されている。
各出力側は、ピーク検出器１５の入力側に接続されてい
る。ピーク検出器１５は、短時間積分器１３から送られ
た短時間積分結果及び長時間積分器１４から送られた長
時間積分結果に基づいて強度が極大となるピークを検出
する手段である。ピーク検出器１５の出力側は、入力端
子１１から入力する入力信号を分析してピーク検出器１
５で検出されたピークに対応する正確な周波数を測定す
る周波数精測器１６に接続されている。又、ピーク検出
器１５の出力側は、ピーク検出器１５から送られたピー
ク検出情報及び周波数精測器１６から送られた周波数情
報により、同一パルスの成分と判断されるピークの時間
的連なりからパルスを形成するパルス形成器１７の入力
側に接続されている。パルス形成器１７の出力側は、パ
ルス形成器１７で形成されたパルス信号の周波数変化の
仕方を線形回帰により推定する勾配推定器１８の入力側
に接続されている。勾配推定器１８の出力側は、勾配推
定器１８で推定されたパルス信号の周波数変化の勾配値
と、予め設定された閾値とからパルス信号の変調形式を
判定するパルス判定処理器１９の入力側に接続されてい
る。パルス判定処理器１９の出力側は、パルス信号の各
種パラメータの出力端子２０に接続されている。
【０００４】次に、この信号形式判定装置の動作を説明
する。入力端子１１から入力された広帯域成分を含んだ
入力信号は、周波数分析器１２で一定時間間隔ｄｔでサ
ンプリングされ、周波数分析される。周波数分析の方法
としては、例えば、ＦＦＴ（Fast Fourier Transform、
高速フーリエ変換）等の方法がある。即ち、入力端子１
１からの入力信号をｘ（ｔ）とし、その離散形を次の
（１）式で表す。 ｘ（ｔi ），ｔi ＝ｔ0 ＋ｉ・ｄｔ ・・・（１） 但し、 ｉ＝０，１，・・・ 一定時間間隔Δｔ毎に周波数分析結果を出力するものと
し、分析周期ｋにおける周波数サンプルｆn の周波数分
析出力信号を次の（２）項のように表す。 Ｘ（ｆn ；ｋ） ・・・（２） 但し、 ｎ＝０，１，・・・，Ｎ−１ 周波数分析結果Ｙ（ｆn ；ｋ）は、実部と虚部との二乗
和の形にして次の（３）式で表す。 Ｙ（ｆn ；ｋ）＝｜Ｘ（ｆn ；ｋ）｜2 ・・・（３） この周波数分析結果Ｙ（ｆn ；ｋ）は、短時間積分器１
３及び長時間積分器１４へ送られる。
する。入力端子１１から入力された広帯域成分を含んだ
入力信号は、周波数分析器１２で一定時間間隔ｄｔでサ
ンプリングされ、周波数分析される。周波数分析の方法
としては、例えば、ＦＦＴ（Fast Fourier Transform、
高速フーリエ変換）等の方法がある。即ち、入力端子１
１からの入力信号をｘ（ｔ）とし、その離散形を次の
（１）式で表す。 ｘ（ｔi ），ｔi ＝ｔ0 ＋ｉ・ｄｔ ・・・（１） 但し、 ｉ＝０，１，・・・ 一定時間間隔Δｔ毎に周波数分析結果を出力するものと
し、分析周期ｋにおける周波数サンプルｆn の周波数分
析出力信号を次の（２）項のように表す。 Ｘ（ｆn ；ｋ） ・・・（２） 但し、 ｎ＝０，１，・・・，Ｎ−１ 周波数分析結果Ｙ（ｆn ；ｋ）は、実部と虚部との二乗
和の形にして次の（３）式で表す。 Ｙ（ｆn ；ｋ）＝｜Ｘ（ｆn ；ｋ）｜2 ・・・（３） この周波数分析結果Ｙ（ｆn ；ｋ）は、短時間積分器１
３及び長時間積分器１４へ送られる。
【０００５】短時間積分器１３では、周波数分析器１２
から送られた周波数分析結果Ｙ（ｆn ；ｋ）を予め設定
された短時間の積分時定数τ1 で周波数ビン毎に積分
し、その結果をピーク検出器１５に送る。一方、長時間
積分器１４では、周波数分析器１２から送られた周波数
分析結果Ｙ（ｆn ；ｋ）を予め設定された長時間の積分
時定数τ2 （但し、τ2 ≫τ1 ）で周波数ビン毎に積分
し、その結果をピーク検出器１５に送る。短時間積分器
１３の積分時定数τ1 は、積分結果が入力信号の特徴を
失わないように短く設定され、長時間積分器１４の積分
時定数τ2 は、ピーク検出器１５のピーク検出における
背景の平均レベルを推定するために長く設定される。積
分は、有限の時間窓を用いてもよいし、瞬時値の影響が
指数関数的に減少する無限の時間窓を用いてもよい。例
えば、無限の時間窓を用いた場合、次の（４）式の形で
積分処理を実行する。 Ｙ1 （ｆn ；ｋ）＝α・Ｙ（ｆn ；ｋ） ＋（１−α）・Ｙ1 （ｆn ；ｋ−１） Ｙ2 （ｆn ；ｋ）＝β・Ｙ（ｆn ；ｋ） ＋（１−β）・Ｙ2 （ｆn ；ｋ−１）・・・（４） 但し、 Ｙ1 ；短時間積分器３の出力信号 α；積分時定数τ1 から定められる定数 Ｙ2 ；長時間積分器４の出力信号 β；積分時定数τ2 から定められる定数 ピーク検出器１５は、短時間積分器１３及び長時間積分
器１４から各積分結果を入力し、次の（５）式に示すよ
うに両者の出力信号の比を予め定めた閾値γと比較す
る。 Ｙ1 ／Ｙ2 ≧γ ・・・（５） この（５）式が成立し、かつ周波数軸上でピークになっ
ている場合には、ピーク検出器１５は、この周波数近傍
にパルス信号の周波数成分が存在すると判断し、周波数
ビンの中心周波数ｆnp（ｋ）又は周波数ビン番号ｎ
（ｋ）を周波数精測器１６へ送り、又、周波数ビン番号
ｎ（ｋ）をパルス形成器１７へ送る。
から送られた周波数分析結果Ｙ（ｆn ；ｋ）を予め設定
された短時間の積分時定数τ1 で周波数ビン毎に積分
し、その結果をピーク検出器１５に送る。一方、長時間
積分器１４では、周波数分析器１２から送られた周波数
分析結果Ｙ（ｆn ；ｋ）を予め設定された長時間の積分
時定数τ2 （但し、τ2 ≫τ1 ）で周波数ビン毎に積分
し、その結果をピーク検出器１５に送る。短時間積分器
１３の積分時定数τ1 は、積分結果が入力信号の特徴を
失わないように短く設定され、長時間積分器１４の積分
時定数τ2 は、ピーク検出器１５のピーク検出における
背景の平均レベルを推定するために長く設定される。積
分は、有限の時間窓を用いてもよいし、瞬時値の影響が
指数関数的に減少する無限の時間窓を用いてもよい。例
えば、無限の時間窓を用いた場合、次の（４）式の形で
積分処理を実行する。 Ｙ1 （ｆn ；ｋ）＝α・Ｙ（ｆn ；ｋ） ＋（１−α）・Ｙ1 （ｆn ；ｋ−１） Ｙ2 （ｆn ；ｋ）＝β・Ｙ（ｆn ；ｋ） ＋（１−β）・Ｙ2 （ｆn ；ｋ−１）・・・（４） 但し、 Ｙ1 ；短時間積分器３の出力信号 α；積分時定数τ1 から定められる定数 Ｙ2 ；長時間積分器４の出力信号 β；積分時定数τ2 から定められる定数 ピーク検出器１５は、短時間積分器１３及び長時間積分
器１４から各積分結果を入力し、次の（５）式に示すよ
うに両者の出力信号の比を予め定めた閾値γと比較す
る。 Ｙ1 ／Ｙ2 ≧γ ・・・（５） この（５）式が成立し、かつ周波数軸上でピークになっ
ている場合には、ピーク検出器１５は、この周波数近傍
にパルス信号の周波数成分が存在すると判断し、周波数
ビンの中心周波数ｆnp（ｋ）又は周波数ビン番号ｎ
（ｋ）を周波数精測器１６へ送り、又、周波数ビン番号
ｎ（ｋ）をパルス形成器１７へ送る。
【０００６】周波数精測器１６では、入力端子１１から
入力した入力信号を詳細に周波数分析し、ピーク検出器
１５から送られた中心周波数ｆnp（ｋ）又は周波数ビン
番号ｎ（ｋ）に基づいてピーク周波数を精測し、その精
測値Ｆ（ｋ）をパルス形成器１７へ送る。ピーク周波数
の精測方法としては、例えば、周波数精測器１６での周
波数分析結果の中から中心周波数ｆnp（ｋ）に最も近い
周波数をもつ極大点を抽出し、この極大点及びその両隣
の周波数ビンの周波数分析結果から放物線補間により算
出する方法がある。パルス形成器１７では、ピーク検出
器５からの周波数ビン番号ｎ（ｋ）及び周波数精測器１
６からの精測値Ｆ（ｋ）を入力し、周波数ビン番号ｎ
（ｋ）に基づき、時間的かつ周波数的に連続したピーク
の時間的連なりを１つのパルス信号として形成する。こ
の方法としては、例えば、連続した分析周期に生じたピ
ークの周波数ビン番号ｎ（ｋ）と周波数ビン番号ｎ（ｋ
−１）とを比較し、その差が予め定めた基準値以下であ
れば同一パルスの要素と見做す方法がある。又、（ｋ−
１）までの分析周期のデータから同一パルスと見做され
るピーク値の連なり（以下、トラックという）を形成
し、トラック内の周波数又は周波数ビンの時間的な変化
の仕方から分析周期ｋにおける精測値Ｆ（ｋ）又は周波
数ビン番号ｎ（ｋ）を予測し、予測値と実際に得られた
値との比較により同一パルスに属するピークか否かを判
断する方法もある。
入力した入力信号を詳細に周波数分析し、ピーク検出器
１５から送られた中心周波数ｆnp（ｋ）又は周波数ビン
番号ｎ（ｋ）に基づいてピーク周波数を精測し、その精
測値Ｆ（ｋ）をパルス形成器１７へ送る。ピーク周波数
の精測方法としては、例えば、周波数精測器１６での周
波数分析結果の中から中心周波数ｆnp（ｋ）に最も近い
周波数をもつ極大点を抽出し、この極大点及びその両隣
の周波数ビンの周波数分析結果から放物線補間により算
出する方法がある。パルス形成器１７では、ピーク検出
器５からの周波数ビン番号ｎ（ｋ）及び周波数精測器１
６からの精測値Ｆ（ｋ）を入力し、周波数ビン番号ｎ
（ｋ）に基づき、時間的かつ周波数的に連続したピーク
の時間的連なりを１つのパルス信号として形成する。こ
の方法としては、例えば、連続した分析周期に生じたピ
ークの周波数ビン番号ｎ（ｋ）と周波数ビン番号ｎ（ｋ
−１）とを比較し、その差が予め定めた基準値以下であ
れば同一パルスの要素と見做す方法がある。又、（ｋ−
１）までの分析周期のデータから同一パルスと見做され
るピーク値の連なり（以下、トラックという）を形成
し、トラック内の周波数又は周波数ビンの時間的な変化
の仕方から分析周期ｋにおける精測値Ｆ（ｋ）又は周波
数ビン番号ｎ（ｋ）を予測し、予測値と実際に得られた
値との比較により同一パルスに属するピークか否かを判
断する方法もある。
【０００７】パルス形成器１７は、同一パルスと見做せ
るピーク値の連なりを形成すると、その周波数精測値の
組Ｆ（ｋ），ｋ＝ｍ，・・・，ｍ＋Ｌ−１（ｍ；ピーク
値の連なりの頭の分析周期の番号、Ｌ；ピーク値の連な
りの個数）を勾配推定器１８へ出力する。勾配推定器１
８では、パルス内の周波数変化の仕方を次の（６）式で
表す。 Ｆ（ｔ）＝Ｆ0 ＋ｄＦ（ｔ−ｔm ） ・・・（６） 但し、 ｔm ；ピークの連なりの頭の分析周期ｍの中心時刻 Ｆ0 ；初期周波数 ｄＦ；周波数勾配 上記（６）式で表したときの初期周波数Ｆ0 と周波数勾
配ｄＦとを、パルス形成器１７から送られた周波数精測
値Ｆ（ｋ）の組から、次の（７）式及び（８）式により
推定する。
るピーク値の連なりを形成すると、その周波数精測値の
組Ｆ（ｋ），ｋ＝ｍ，・・・，ｍ＋Ｌ−１（ｍ；ピーク
値の連なりの頭の分析周期の番号、Ｌ；ピーク値の連な
りの個数）を勾配推定器１８へ出力する。勾配推定器１
８では、パルス内の周波数変化の仕方を次の（６）式で
表す。 Ｆ（ｔ）＝Ｆ0 ＋ｄＦ（ｔ−ｔm ） ・・・（６） 但し、 ｔm ；ピークの連なりの頭の分析周期ｍの中心時刻 Ｆ0 ；初期周波数 ｄＦ；周波数勾配 上記（６）式で表したときの初期周波数Ｆ0 と周波数勾
配ｄＦとを、パルス形成器１７から送られた周波数精測
値Ｆ（ｋ）の組から、次の（７）式及び（８）式により
推定する。
【数１】
【数２】 但し、 ｋ＝ｍ，・・・，ｍ＋Ｌ−１、 ｍ；ピーク値の連なりの頭の分析周期の番号 Ｌ；ピーク値の連なりの個数 （７）式を解くと次の（９）式が得られる。 Ｆ0 ＝（ＢＣ−ＡＤ）／（ＬＢ−Ａ2 ） ｄＦ＝（ＬＤ−ＡＣ）／（ＬＢ−Ａ2 ） ・・・（９） 尚、入力端子１１からパルス形成器１７までの動作がパ
ルス信号形成処理である。勾配推定器１８は、初期周波
数Ｆ0 と周波数勾配ｄＦとを算出すると、これらの値を
パルスの時間幅ΔＴ＝（ｔm+L-1 −ｔm ）とともにパル
ス判定処理器１９へ出力する。パルス判定処理におい
て、パルス判定処理器１９では、勾配推定器１８から入
力した各パラメータを用いて周波数勾配ｄＦ又はパルス
信号内の周波数変化幅ΔＦ＝ｄＦ・ΔＴを予め定めた閾
値と比較し、パルス信号が一定周波数の正弦波（Pulsed
Continuous Wave、以下、ＰＣＷという）であるか或い
は周波数が時間と共に線形に変化する周波数変調波（Li
near Frequency Modulation 、以下、ＬＦＭという）で
あるかを判定する。
ルス信号形成処理である。勾配推定器１８は、初期周波
数Ｆ0 と周波数勾配ｄＦとを算出すると、これらの値を
パルスの時間幅ΔＴ＝（ｔm+L-1 −ｔm ）とともにパル
ス判定処理器１９へ出力する。パルス判定処理におい
て、パルス判定処理器１９では、勾配推定器１８から入
力した各パラメータを用いて周波数勾配ｄＦ又はパルス
信号内の周波数変化幅ΔＦ＝ｄＦ・ΔＴを予め定めた閾
値と比較し、パルス信号が一定周波数の正弦波（Pulsed
Continuous Wave、以下、ＰＣＷという）であるか或い
は周波数が時間と共に線形に変化する周波数変調波（Li
near Frequency Modulation 、以下、ＬＦＭという）で
あるかを判定する。
【０００８】
【発明が解決しようとする課題】しかしながら、上記の
信号形式判定方法では、次のような問題点があった。即
ち、入力端子１１からＰＣＷやＬＦＭ以外の過渡音が入
力された場合にも、パルス判定処理器１９は、それをＰ
ＣＷ又はＬＦＭとして判定してしまうという問題点があ
った。特に、広い周波数帯域をもつ衝撃音や生物音等が
入力端子に入力された場合、ピーク周波数の時間的変化
が大きいので、周波数勾配ｄＦやパルス内の周波数変化
幅ΔＦが大きく推定されてＬＦＭと見做されることがあ
り、誤った判断となることがある。本発明は、以上述べ
たパルス判定処理器１９がＰＣＷやＬＦＭ以外の信号を
ＰＣＷやＬＦＭとして判定してしまうという問題を解決
するために、周波数勾配の推定値の信頼性に対する判断
機能を設け、広い周波数帯域をもつ過渡音をＰＣＷやＬ
ＦＭ以外の不明信号として判定し、衝撃音や生物音等の
雑音による誤認を除去する信号形式判定方法を提供する
ものである。
信号形式判定方法では、次のような問題点があった。即
ち、入力端子１１からＰＣＷやＬＦＭ以外の過渡音が入
力された場合にも、パルス判定処理器１９は、それをＰ
ＣＷ又はＬＦＭとして判定してしまうという問題点があ
った。特に、広い周波数帯域をもつ衝撃音や生物音等が
入力端子に入力された場合、ピーク周波数の時間的変化
が大きいので、周波数勾配ｄＦやパルス内の周波数変化
幅ΔＦが大きく推定されてＬＦＭと見做されることがあ
り、誤った判断となることがある。本発明は、以上述べ
たパルス判定処理器１９がＰＣＷやＬＦＭ以外の信号を
ＰＣＷやＬＦＭとして判定してしまうという問題を解決
するために、周波数勾配の推定値の信頼性に対する判断
機能を設け、広い周波数帯域をもつ過渡音をＰＣＷやＬ
ＦＭ以外の不明信号として判定し、衝撃音や生物音等の
雑音による誤認を除去する信号形式判定方法を提供する
ものである。
【０００９】
【課題を解決するための手段】本発明では、前記課題を
解決するために、入力信号を一定の時間間隔毎に周波数
分析して強度が極大となるピーク周波数を求め、かつ該
ピーク周波数の時間的連なりとしてパルス信号を形成す
るパルス信号形成処理と、前記パルス信号の周波数の時
間変化を１次関数で近似して該１次関数の時間勾配又は
変化幅に基づいて該パルス信号がＰＣＷか或いはＬＦＭ
かを判定するパルス判定処理とを、施す信号形式判定方
法において、次のような処理を行っている。即ち、前記
ピーク周波数の測定値の前記１次関数に対する残差の二
乗和を算出する推定残差算出処理と、残差の二乗和と、
予め定めた基準値と前記パルス信号形成処理により検出
されたピークの連なりの個数とから計算される残差の二
乗和の許容値とを比較してパルス判定処理の判定結果を
採用するか或いは棄却するかを判定する残差比較処理と
を、行うようにしている。
解決するために、入力信号を一定の時間間隔毎に周波数
分析して強度が極大となるピーク周波数を求め、かつ該
ピーク周波数の時間的連なりとしてパルス信号を形成す
るパルス信号形成処理と、前記パルス信号の周波数の時
間変化を１次関数で近似して該１次関数の時間勾配又は
変化幅に基づいて該パルス信号がＰＣＷか或いはＬＦＭ
かを判定するパルス判定処理とを、施す信号形式判定方
法において、次のような処理を行っている。即ち、前記
ピーク周波数の測定値の前記１次関数に対する残差の二
乗和を算出する推定残差算出処理と、残差の二乗和と、
予め定めた基準値と前記パルス信号形成処理により検出
されたピークの連なりの個数とから計算される残差の二
乗和の許容値とを比較してパルス判定処理の判定結果を
採用するか或いは棄却するかを判定する残差比較処理と
を、行うようにしている。
【００１０】
【作用】本発明によれば、以上のように信号形式判定方
法を構成したので、パルス信号形成処理において、入力
信号を一定の時間間隔毎に周波数分析して強度が極大と
なるピーク周波数を求め、かつ該ピーク周波数の時間的
変化を表すパルス信号を形成する。次に、パルス判定処
理において、前記パルス信号の周波数の時間的変化を１
次関数で近似して該１次関数の時間勾配又は変化幅に基
づいて該パルス信号がＰＣＷか或いはＬＦＭかを判定す
る。推定残差算出処理において、前記ピーク周波数の測
定値の前記１次関数に対する残差の二乗和を算出する。
更に、残差比較処理において、前記残差の二乗和と、予
め定めた基準値と前記パルス信号形成処理により検出さ
れたピークの連なりの個数とから計算される残差の二乗
和の許容値とを比較してパルス判定処理の判定結果の信
頼性について判定し、前記パルス判定処理の判定結果に
基づいてパルス信号がＰＣＷやＬＦＭであると見做すか
或いはＰＣＷやＬＦＭに属さない不明の信号と見做すか
を判断するようにしている。従って、前記課題を解決で
きるのである。
法を構成したので、パルス信号形成処理において、入力
信号を一定の時間間隔毎に周波数分析して強度が極大と
なるピーク周波数を求め、かつ該ピーク周波数の時間的
変化を表すパルス信号を形成する。次に、パルス判定処
理において、前記パルス信号の周波数の時間的変化を１
次関数で近似して該１次関数の時間勾配又は変化幅に基
づいて該パルス信号がＰＣＷか或いはＬＦＭかを判定す
る。推定残差算出処理において、前記ピーク周波数の測
定値の前記１次関数に対する残差の二乗和を算出する。
更に、残差比較処理において、前記残差の二乗和と、予
め定めた基準値と前記パルス信号形成処理により検出さ
れたピークの連なりの個数とから計算される残差の二乗
和の許容値とを比較してパルス判定処理の判定結果の信
頼性について判定し、前記パルス判定処理の判定結果に
基づいてパルス信号がＰＣＷやＬＦＭであると見做すか
或いはＰＣＷやＬＦＭに属さない不明の信号と見做すか
を判断するようにしている。従って、前記課題を解決で
きるのである。
【００１１】
【実施例】先ず、本実施例の基礎となっている原理を説
明する。パルス信号の周波数の時間的変化についての観
測値Ｆ（ｔk ）（但し、ｋ＝ｍ，・・・，ｍ＋Ｌ−１）
を（６）式で近似した場合、推定残差の二乗和Ｇは次の
（１０）式で表される。
明する。パルス信号の周波数の時間的変化についての観
測値Ｆ（ｔk ）（但し、ｋ＝ｍ，・・・，ｍ＋Ｌ−１）
を（６）式で近似した場合、推定残差の二乗和Ｇは次の
（１０）式で表される。
【数３】 （１０）式中の観測値Ｆ（ｔk ）の二乗和Ｈは次の（１
１）式で表される。
１）式で表される。
【数４】 パルス信号がＰＣＷかＬＦＭであったとすると、パルス
信号内での周波数の時間的変化は、従来の（６）式で表
わされ、（１０）式で表される推定残差の二乗和Ｇは非
常に小さくなる。ところが、入力信号が広帯域の衝撃音
や生物音の場合、ピーク周波数は時間に対して線形には
変化せず、（１０）式で表される推定残差の二乗和Ｇは
大きい値をもつ。
信号内での周波数の時間的変化は、従来の（６）式で表
わされ、（１０）式で表される推定残差の二乗和Ｇは非
常に小さくなる。ところが、入力信号が広帯域の衝撃音
や生物音の場合、ピーク周波数は時間に対して線形には
変化せず、（１０）式で表される推定残差の二乗和Ｇは
大きい値をもつ。
【００１２】ここで、パルス信号がＰＣＷかＬＦＭであ
る場合、（１０）式の推定残差の二乗和Ｇの期待値Ｅ
（Ｇ）は、次の（１２）式で表される。 Ｅ（Ｇ）＝Ｌ・σf 2 ・・・（１２） 但し、 Ｌ；パルス信号に属するピークデータの個数 σf 2 ；各ピークの周波数推定値の誤差分散 このσf 2 の値は、ピーク検出の閾値γ、周波数分析器
１２での周波数分析幅、周波数精測器での周波数分析等
により変化する。従って、これらの値を基にピーク周波
数推定誤差分散の許容値σmax 2 を定め、（１０）式の
推定残差の二乗和Ｇの許容値Ｇmax を次の（１３）式で
表す。 Ｇmax ＝Ｌ・σmax 2 ・・・（１３） 更に、次の（１４）式に示すように、推定残差の二乗和
Ｇが許容値Ｇmax よりも大きくなった場合にパルス信号
の変調形式を不明とすれば、入力信号中の広帯域の衝撃
音や生物音を誤ってＰＣＷやＬＦＭと判定する可能性が
大幅に低減する。 Ｇ＞Ｇmax ・・・（１４） 図１は、本発明の実施例の信号形式判定方法を実施する
ための信号形式判定装置の機能ブロック図であり、従来
の図２と共通の要素には共通の符号が付されている。
る場合、（１０）式の推定残差の二乗和Ｇの期待値Ｅ
（Ｇ）は、次の（１２）式で表される。 Ｅ（Ｇ）＝Ｌ・σf 2 ・・・（１２） 但し、 Ｌ；パルス信号に属するピークデータの個数 σf 2 ；各ピークの周波数推定値の誤差分散 このσf 2 の値は、ピーク検出の閾値γ、周波数分析器
１２での周波数分析幅、周波数精測器での周波数分析等
により変化する。従って、これらの値を基にピーク周波
数推定誤差分散の許容値σmax 2 を定め、（１０）式の
推定残差の二乗和Ｇの許容値Ｇmax を次の（１３）式で
表す。 Ｇmax ＝Ｌ・σmax 2 ・・・（１３） 更に、次の（１４）式に示すように、推定残差の二乗和
Ｇが許容値Ｇmax よりも大きくなった場合にパルス信号
の変調形式を不明とすれば、入力信号中の広帯域の衝撃
音や生物音を誤ってＰＣＷやＬＦＭと判定する可能性が
大幅に低減する。 Ｇ＞Ｇmax ・・・（１４） 図１は、本発明の実施例の信号形式判定方法を実施する
ための信号形式判定装置の機能ブロック図であり、従来
の図２と共通の要素には共通の符号が付されている。
【００１３】この信号形式判定装置では、従来の図２中
の勾配推定器１８の出力側に、パルス形成器１７の出力
信号に基づいて勾配推定の推定残差の二乗和を計算する
推定残差算出器２０が接続されている。又、パルス形成
器１７の出力側には、パルス形成器１７の出力信号に基
づいて残差判定の基準値を生成する基準値生成器２１が
接続されている。基準値生成器２１の出力側には、推定
残差算出器２０から送られた推定残差の二乗和と基準値
生成器２１から送られた基準値とを比較する残差比較器
２２が接続されている。残差比較器２２の出力側には、
パルスの各種パラメータの推定結果を出力する出力端子
２３が接続されている。尚、この信号形式判定装置中の
各構成要素２〜２２は、個別回路で構成されるか或いは
コンピュータを用いたプログラムで実現される。
の勾配推定器１８の出力側に、パルス形成器１７の出力
信号に基づいて勾配推定の推定残差の二乗和を計算する
推定残差算出器２０が接続されている。又、パルス形成
器１７の出力側には、パルス形成器１７の出力信号に基
づいて残差判定の基準値を生成する基準値生成器２１が
接続されている。基準値生成器２１の出力側には、推定
残差算出器２０から送られた推定残差の二乗和と基準値
生成器２１から送られた基準値とを比較する残差比較器
２２が接続されている。残差比較器２２の出力側には、
パルスの各種パラメータの推定結果を出力する出力端子
２３が接続されている。尚、この信号形式判定装置中の
各構成要素２〜２２は、個別回路で構成されるか或いは
コンピュータを用いたプログラムで実現される。
【００１４】次に、図１の動作を説明する。この信号形
式判定装置では、入力端子１１から周波数精測器１６ま
での動作は従来と同一であるが、パルス形成器１７以降
の動作が異なっている。即ち、パルス信号形成処理にお
いて、パルス形成器１７は、同一パルスと見做せるピー
ク値の連なりを形成すると、ピーク値の連なりの個数
Ｌ、分析周期ｋの中心時刻ｔk 、及び周波数精測値の組
Ｆ（ｋ），ｋ＝ｍ，・・・，ｍ＋Ｌ−１（ｍ；ピーク値
の連なりの頭の分析周期の番号、Ｌ；ピーク値の連なり
の個数）を勾配推定器１８及び推定残差算出器２０へ出
力すると共に、ピーク値の連なりの個数Ｌを基準値生成
器２１へ出力する。勾配推定器１８では、パルス信号内
の周波数変化の仕方を（６）式で表し、パルス形成器１
７からの精測値Ｆ（ｋ）から、（７）式及び（８）式に
より初期周波数Ｆ0 と周波数勾配ｄＦとを推定する。勾
配推定器１８は、初期周波数Ｆ0 と周波数勾配ｄＦとを
算出すると、これらの値をパルス信号の時間幅ΔＴ＝
（ｔm+L-1 −ｔm ）とともにパルス判定処理器１９へ出
力する。又、初期周波数Ｆ0 、周波数勾配ｄＦ、及び
（８）式中のＣを推定残差算出器２０へ出力する。パル
ス判定処理において、パルス判定処理器１９では、勾配
推定器１８から入力した初期周波数Ｆ0 、周波数勾配ｄ
Ｆ、及びパルス信号の時間幅ΔＴの各パラメータを用い
て周波数勾配ｄＦ又はパルス内の周波数変化幅ΔＦ＝ｄ
Ｆ・ΔＴを予め定めた閾値と比較し、パルス信号がＰＣ
ＷであるかＬＦＭであるかを判定する。この判定結果
は、初期周波数Ｆ0 及び周波数勾配ｄＦの推定結果と共
に残差比較器２２へ出力する。
式判定装置では、入力端子１１から周波数精測器１６ま
での動作は従来と同一であるが、パルス形成器１７以降
の動作が異なっている。即ち、パルス信号形成処理にお
いて、パルス形成器１７は、同一パルスと見做せるピー
ク値の連なりを形成すると、ピーク値の連なりの個数
Ｌ、分析周期ｋの中心時刻ｔk 、及び周波数精測値の組
Ｆ（ｋ），ｋ＝ｍ，・・・，ｍ＋Ｌ−１（ｍ；ピーク値
の連なりの頭の分析周期の番号、Ｌ；ピーク値の連なり
の個数）を勾配推定器１８及び推定残差算出器２０へ出
力すると共に、ピーク値の連なりの個数Ｌを基準値生成
器２１へ出力する。勾配推定器１８では、パルス信号内
の周波数変化の仕方を（６）式で表し、パルス形成器１
７からの精測値Ｆ（ｋ）から、（７）式及び（８）式に
より初期周波数Ｆ0 と周波数勾配ｄＦとを推定する。勾
配推定器１８は、初期周波数Ｆ0 と周波数勾配ｄＦとを
算出すると、これらの値をパルス信号の時間幅ΔＴ＝
（ｔm+L-1 −ｔm ）とともにパルス判定処理器１９へ出
力する。又、初期周波数Ｆ0 、周波数勾配ｄＦ、及び
（８）式中のＣを推定残差算出器２０へ出力する。パル
ス判定処理において、パルス判定処理器１９では、勾配
推定器１８から入力した初期周波数Ｆ0 、周波数勾配ｄ
Ｆ、及びパルス信号の時間幅ΔＴの各パラメータを用い
て周波数勾配ｄＦ又はパルス内の周波数変化幅ΔＦ＝ｄ
Ｆ・ΔＴを予め定めた閾値と比較し、パルス信号がＰＣ
ＷであるかＬＦＭであるかを判定する。この判定結果
は、初期周波数Ｆ0 及び周波数勾配ｄＦの推定結果と共
に残差比較器２２へ出力する。
【００１５】推定残差算出処理において、推定残差算出
器２０では、パルス形成器１７からピーク値の連なりの
個数Ｌ及び精測値Ｆ（ｋ）（但し、ｋ＝ｍ，・・・，ｍ
＋Ｌ−１）を入力し、（１１）式により観測値Ｆ
（ｔk ）の二乗和Ｈを計算し、この観測値Ｆ（ｔk ）の
二乗和Ｈと勾配推定器１８から入力した初期周波数
Ｆ0 、周波数勾配ｄＦ、及び（８）式中のＣとを用い
て、（１０）式により推定残差の二乗和Ｇを算出して残
差比較器２２へ出力する。基準値生成器２１では、パル
ス形成器１７からピーク値の連なりの個数Ｌを入力し、
予め設定されている周波数推定誤差分散の許容値σmax
2 から、（１３）式により推定残差の二乗和Ｇの許容値
Ｇmax を算出して残差比較器２２へ出力する。残差比較
処理において、残差比較器２２では、推定残差算出器２
０からの推定残差の二乗和Ｇと基準値生成器２１からの
許容値Ｇmax とを比較し、次の判定基準（ａ）及び
（ｂ）に基づいてパルス信号の形式を判定し、中心周波
数や勾配等の推定結果と共に出力端子２３から出力す
る。 （ａ）Ｇ＞Ｇmax のとき；パルス信号の形式を不明とす
る。 （ｂ）Ｇ≦Ｇmax のとき；パルス信号の形式をパルス判
定処理器１９で判定したとおりとする。
器２０では、パルス形成器１７からピーク値の連なりの
個数Ｌ及び精測値Ｆ（ｋ）（但し、ｋ＝ｍ，・・・，ｍ
＋Ｌ−１）を入力し、（１１）式により観測値Ｆ
（ｔk ）の二乗和Ｈを計算し、この観測値Ｆ（ｔk ）の
二乗和Ｈと勾配推定器１８から入力した初期周波数
Ｆ0 、周波数勾配ｄＦ、及び（８）式中のＣとを用い
て、（１０）式により推定残差の二乗和Ｇを算出して残
差比較器２２へ出力する。基準値生成器２１では、パル
ス形成器１７からピーク値の連なりの個数Ｌを入力し、
予め設定されている周波数推定誤差分散の許容値σmax
2 から、（１３）式により推定残差の二乗和Ｇの許容値
Ｇmax を算出して残差比較器２２へ出力する。残差比較
処理において、残差比較器２２では、推定残差算出器２
０からの推定残差の二乗和Ｇと基準値生成器２１からの
許容値Ｇmax とを比較し、次の判定基準（ａ）及び
（ｂ）に基づいてパルス信号の形式を判定し、中心周波
数や勾配等の推定結果と共に出力端子２３から出力す
る。 （ａ）Ｇ＞Ｇmax のとき；パルス信号の形式を不明とす
る。 （ｂ）Ｇ≦Ｇmax のとき；パルス信号の形式をパルス判
定処理器１９で判定したとおりとする。
【００１６】以上のように、本実施例では、パルス信号
内での周波数の時間変化特性の観測値Ｆ（ｔk ）に対し
て推定残差の二乗和Ｇを求め、推定残差の二乗和Ｇが許
容値Ｇmax を越えた場合にパルス信号の信号形式を不明
としてＰＣＷやＬＦＭと区別するようにしたので、パル
ス判定処理器１９がＰＣＷやＬＦＭ以外の信号をＰＣＷ
やＬＦＭとして判定して入力信号中の衝撃音や生物音等
の雑音がＰＣＷやＬＦＭ等の信号音に統合される可能性
が大幅に低減され、目標の信号音を正しく分析する追尾
システムを構成できる。尚、本発明は上記実施例に限定
されず、例えば、目標から放射されたレーダ信号を受信
するシステム等にも応用が可能である。
内での周波数の時間変化特性の観測値Ｆ（ｔk ）に対し
て推定残差の二乗和Ｇを求め、推定残差の二乗和Ｇが許
容値Ｇmax を越えた場合にパルス信号の信号形式を不明
としてＰＣＷやＬＦＭと区別するようにしたので、パル
ス判定処理器１９がＰＣＷやＬＦＭ以外の信号をＰＣＷ
やＬＦＭとして判定して入力信号中の衝撃音や生物音等
の雑音がＰＣＷやＬＦＭ等の信号音に統合される可能性
が大幅に低減され、目標の信号音を正しく分析する追尾
システムを構成できる。尚、本発明は上記実施例に限定
されず、例えば、目標から放射されたレーダ信号を受信
するシステム等にも応用が可能である。
【００１７】
【発明の効果】以上詳細に説明したように、本発明によ
れば、パルス信号内での周波数の時間的変化の推定値の
信頼性に対する判断機能を設けて周波数の時間変化特性
の観測値に対して推定残差の二乗和を求め、推定残差の
二乗和が許容値を越えた場合に信号特性を不明としてＰ
ＣＷやＬＦＭと区別するようにしたので、入力信号中の
衝撃音や生物音等の雑音がＰＣＷやＬＦＭ等の信号音に
統合される可能性が大幅に低減され、目標の信号音を正
しく分析する追尾システムを構成できる。
れば、パルス信号内での周波数の時間的変化の推定値の
信頼性に対する判断機能を設けて周波数の時間変化特性
の観測値に対して推定残差の二乗和を求め、推定残差の
二乗和が許容値を越えた場合に信号特性を不明としてＰ
ＣＷやＬＦＭと区別するようにしたので、入力信号中の
衝撃音や生物音等の雑音がＰＣＷやＬＦＭ等の信号音に
統合される可能性が大幅に低減され、目標の信号音を正
しく分析する追尾システムを構成できる。
【図１】本発明の実施例を示す信号形式判定方法を実施
するための信号形式判定装置の機能ブロック図である。
するための信号形式判定装置の機能ブロック図である。
【図２】従来の信号形式判定方法を実施するための信号
形式判定装置の機能ブロック図である。
形式判定装置の機能ブロック図である。
１１ 入力端
子 １２ 周波数
分析器 １３ 短時間
積分器 １４ 長時間
積分器 １５ ピーク
検出器 １６ 周波数
精測器 １７ パルス
形成器 １８ 勾配推
定器 １９ パルス
判定処理器 ２０ 推定残
差算出器 ２１ 基準値
生成器 ２２ 残差比
較器 ２３ 出力端
子
子 １２ 周波数
分析器 １３ 短時間
積分器 １４ 長時間
積分器 １５ ピーク
検出器 １６ 周波数
精測器 １７ パルス
形成器 １８ 勾配推
定器 １９ パルス
判定処理器 ２０ 推定残
差算出器 ２１ 基準値
生成器 ２２ 残差比
較器 ２３ 出力端
子
───────────────────────────────────────────────────── フロントページの続き (72)発明者 山下 正人 東京都港区虎ノ門１丁目７番12号 沖電 気工業株式会社内 (72)発明者 屋内 伸治 東京都港区虎ノ門１丁目７番12号 沖電 気工業株式会社内 (56)参考文献 特開 平４−64080（ＪＰ，Ａ) 特開 平５−134025（ＪＰ，Ａ) 特開 昭64−1996（ＪＰ，Ａ) 特開 昭61−91579（ＪＰ，Ａ)
Claims (1)
- 【請求項１】 入力信号を一定の時間間隔毎に周波数分
析して強度が極大となるピーク周波数を求め、かつ該ピ
ーク周波数の時間的連なりとしてパルス信号を形成する
パルス信号形成処理と、 前記パルス信号の周波数の時間変化を１次関数で近似し
て該１次関数の時間勾配又は周波数変化幅に基づいて該
パルス信号が一定周波数の正弦波か或いは周波数が時間
と共に線形に変化する周波数変調波かを判定するパルス
判定処理とを、 施す信号形式判定方法において、 前記ピーク周波数の測定値の前記１次関数に対する残差
の二乗和を算出する推定残差算出処理と、 前記残差の二乗和と、予め定めた基準値と前記パルス信
号形成処理により検出されたピークの連なりの個数とか
ら計算される残差の二乗和の許容値とを比較して前記パ
ルス判定処理の判定結果を採用するか或いは棄却するか
を判定する残差比較処理とを、 行うことを特徴とする信号形式判定方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP7024928A JP2607861B2 (ja) | 1995-02-14 | 1995-02-14 | 信号形式判定方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP7024928A JP2607861B2 (ja) | 1995-02-14 | 1995-02-14 | 信号形式判定方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH08220215A JPH08220215A (ja) | 1996-08-30 |
| JP2607861B2 true JP2607861B2 (ja) | 1997-05-07 |
Family
ID=12151773
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP7024928A Expired - Lifetime JP2607861B2 (ja) | 1995-02-14 | 1995-02-14 | 信号形式判定方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2607861B2 (ja) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2007050289A1 (en) * | 2005-10-21 | 2007-05-03 | Raytheon Company | Sonar system and method providing low probability of impact on marine mammals |
| JP6205685B2 (ja) * | 2012-08-02 | 2017-10-04 | 日本電気株式会社 | アクティブソーナー装置、アクティブソーナー信号処理方法、及び、その信号処理プログラム |
-
1995
- 1995-02-14 JP JP7024928A patent/JP2607861B2/ja not_active Expired - Lifetime
Also Published As
| Publication number | Publication date |
|---|---|
| JPH08220215A (ja) | 1996-08-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US5686919A (en) | Process for generating wind profiler data free of fixed ground clutter contamination | |
| US20080231496A1 (en) | Method for determining noise floor level and radar using the same | |
| JP2637439B2 (ja) | 流体分布検出装置 | |
| JP7499045B2 (ja) | 融合された閾値／位相変調検出を使用してパルスを検出するためのシステムおよび方法 | |
| CN110988834B (zh) | 一种基于自适应阈值双门限的脉冲到达时间测量方法 | |
| WO2019032371A1 (en) | SYSTEM AND METHOD FOR INTERFERENCE DETECTION IN AN RF RECEIVER | |
| AU2019229361B2 (en) | Input count rate estimation in radiation pulse detectors | |
| EP1788459B1 (en) | A time-delay discriminator | |
| EP1744173A2 (en) | Methods and apparatus for target radial extent determination using deconvolution | |
| JP2004508557A (ja) | バイナリインテグレーターの検波閾値の適応制御 | |
| JP2607861B2 (ja) | 信号形式判定方法 | |
| RU2376581C2 (ru) | Способ и устройство для получения информации о распределении макроскопических частиц в жидкости по размеру | |
| CN105223552B (zh) | 基于无线电谱纹识别的干扰辨识方法 | |
| JP6398689B2 (ja) | 信号処理装置、信号処理方法及びプログラム | |
| Chen et al. | Long-range channel prediction based on nonstationary parametric modeling | |
| US20200173904A1 (en) | Method and Apparatus for Detecting Particles | |
| JP2607854B2 (ja) | 信号追尾方法 | |
| JP3394412B2 (ja) | パルス音探知方法及びその装置 | |
| JP2002311131A (ja) | 車両の物体検知装置 | |
| JPH04318700A (ja) | 閾値決定手段を有する車間距離測定装置 | |
| JP2603758B2 (ja) | キャビテーション雑音検出方法及びその装置 | |
| Angrisani et al. | Performance assessment of a Kalman filter-based method for ultrasonic time-of-flight estimation | |
| RU2747577C1 (ru) | Способ неразрушающего обнаружения импульсных радиосигналов в условиях узкополосных помех | |
| JPH07128349A (ja) | 速度検出装置 | |
| Lucena et al. | An innovative ultrasonic time of flight method based on extended Kalman filter for wind speed measurement |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 19961224 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| EXPY | Cancellation because of completion of term |