JP2001077677A - パルス列分離装置及びパルス列分離方法 - Google Patents
パルス列分離装置及びパルス列分離方法Info
- Publication number
- JP2001077677A JP2001077677A JP25162399A JP25162399A JP2001077677A JP 2001077677 A JP2001077677 A JP 2001077677A JP 25162399 A JP25162399 A JP 25162399A JP 25162399 A JP25162399 A JP 25162399A JP 2001077677 A JP2001077677 A JP 2001077677A
- Authority
- JP
- Japan
- Prior art keywords
- pulse train
- pulse
- pri
- separating
- staggered
- 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
- Manipulation Of Pulses (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
(57)【要約】
【課題】 信号源の数が多い場合や、固定ＰＲＩ以外の
複雑なＰＲＩのパターンを有する信号源が存在する場合
にも、信号源のＰＲＩを正確に推定し、精度よくパルス
列を分離させるパルス列分離装置を得る。 【解決手段】 ＰＲＩの種類ごとにそれに応じた検出手
段を用いて、段階的に検出・分離していくようにしたも
のであり、ＰＲＩ変換を用いる固定ＰＲＩパルス列検出
・分離手段３００と、３連パルスのヒストグラムを用い
るスタガーパルス列検出・分離手段４００と、変形ＰＲ
Ｉ変換を用いるジッターパルス列検出・分離手段を備え
ている。
複雑なＰＲＩのパターンを有する信号源が存在する場合
にも、信号源のＰＲＩを正確に推定し、精度よくパルス
列を分離させるパルス列分離装置を得る。 【解決手段】 ＰＲＩの種類ごとにそれに応じた検出手
段を用いて、段階的に検出・分離していくようにしたも
のであり、ＰＲＩ変換を用いる固定ＰＲＩパルス列検出
・分離手段３００と、３連パルスのヒストグラムを用い
るスタガーパルス列検出・分離手段４００と、変形ＰＲ
Ｉ変換を用いるジッターパルス列検出・分離手段を備え
ている。
Description
【０００１】
【発明の属する技術分野】この発明は、複数の信号源か
らのパルス列が重なり合った複合信号を、信号源ごとの
パルス列に分離するパルス列分離装置及びパルス列分離
方法に関するものである。
らのパルス列が重なり合った複合信号を、信号源ごとの
パルス列に分離するパルス列分離装置及びパルス列分離
方法に関するものである。
【０００２】
【従来の技術】従来から、複数の信号源からのパルス列
が重なり合った複合信号を、信号源ごとのパルス列に分
離するために、パルスの到着時間（以下、ＴＯＡ（Time
of Arrival）とする）データに基づくパルス列分離装
置が提案されてきた。
が重なり合った複合信号を、信号源ごとのパルス列に分
離するために、パルスの到着時間（以下、ＴＯＡ（Time
of Arrival）とする）データに基づくパルス列分離装
置が提案されてきた。
【０００３】図１６は、例えば米国特許第５、５８３，
５０５号公報に示された、パルスのＴＯＡデータからパ
ルス列を検出し分離するための従来のパルス列分離装置
における処理を示したフローチャートである。
５０５号公報に示された、パルスのＴＯＡデータからパ
ルス列を検出し分離するための従来のパルス列分離装置
における処理を示したフローチャートである。
【０００４】図１６において、Ｓ６０は自己相関を求め
るステップ、Ｓ６２はパルス繰り返し間隔（以下、ＰＲ
Ｉ（Pulse Repetition Interval）とする）の候補を求
めるステップ、Ｓ６４は積分器にメモリーを割り当てる
ステップ、Ｓ６６は各パルスをすべての積分器のビンに
積み重ねるステップ、Ｓ６８は最小ＰＲＩの積分器を分
析するステップ、Ｓ７０は積分器のピークとしきい値を
比較するステップ、Ｓ７２は積分器のピークがしきい値
を越えているとき該当するＰＲＩが存在すると判定する
ステップ、Ｓ７４はそのＰＲＩを有するパルス列を時間
窓から除去するステップ、Ｓ７６は自己相関を再度求め
るステップ、Ｓ７７は次のＰＲＩの候補の積分器に進む
ステップ、Ｓ７８は分析すべき積分器がまだ残っている
か否か判定するステップ、Ｓ７９は残ったパルスが単独
のものであるか否かを判定するステップである。
るステップ、Ｓ６２はパルス繰り返し間隔（以下、ＰＲ
Ｉ（Pulse Repetition Interval）とする）の候補を求
めるステップ、Ｓ６４は積分器にメモリーを割り当てる
ステップ、Ｓ６６は各パルスをすべての積分器のビンに
積み重ねるステップ、Ｓ６８は最小ＰＲＩの積分器を分
析するステップ、Ｓ７０は積分器のピークとしきい値を
比較するステップ、Ｓ７２は積分器のピークがしきい値
を越えているとき該当するＰＲＩが存在すると判定する
ステップ、Ｓ７４はそのＰＲＩを有するパルス列を時間
窓から除去するステップ、Ｓ７６は自己相関を再度求め
るステップ、Ｓ７７は次のＰＲＩの候補の積分器に進む
ステップ、Ｓ７８は分析すべき積分器がまだ残っている
か否か判定するステップ、Ｓ７９は残ったパルスが単独
のものであるか否かを判定するステップである。
【０００５】次に、前述した従来のパルス列分離装置の
動作について図面を参照しながら説明する。
動作について図面を参照しながら説明する。
【０００６】図１６において、まず時間窓内にあるパル
スのＴＯＡデータの自己相関関数を求める（ステップＳ
６０）。得られた自己相関関数のピークの判定を行って
ピークの面積（ピークの高さと幅の積）がしきい値を超
えていれば、そのピークの位置をＰＲＩの候補として抽
出する（ステップＳ６２）。次に、これらＰＲＩの候補
ごとに積分器を割り当てる（ステップＳ６４）。積分器
はＰＲＩ分の時間幅をもつテーブルから構成される。積
分器では、入力した各パルスに対してそのＴＯＡをＰＲ
Ｉで割った余りの時間に対応するテーブル内の位置を求
め、テーブルの値をカウントアップしていく。時間窓内
のすべてのパルスに対して、各積分器で積分を行う（ス
テップＳ６６）。
スのＴＯＡデータの自己相関関数を求める（ステップＳ
６０）。得られた自己相関関数のピークの判定を行って
ピークの面積（ピークの高さと幅の積）がしきい値を超
えていれば、そのピークの位置をＰＲＩの候補として抽
出する（ステップＳ６２）。次に、これらＰＲＩの候補
ごとに積分器を割り当てる（ステップＳ６４）。積分器
はＰＲＩ分の時間幅をもつテーブルから構成される。積
分器では、入力した各パルスに対してそのＴＯＡをＰＲ
Ｉで割った余りの時間に対応するテーブル内の位置を求
め、テーブルの値をカウントアップしていく。時間窓内
のすべてのパルスに対して、各積分器で積分を行う（ス
テップＳ６６）。
【０００７】次いで信号源の検出を行うために、小さい
ＰＲＩの積分器から順番に分析していく（ステップＳ６
８）。積分器のテーブルのピークの値がしきい値を超え
ているか否かを調べ（ステップＳ７０）、超えていれば
そのＰＲＩを持つ信号源が存在すると判定する（ステッ
プＳ７２）。そのあとそのＰＲＩをもつパルス列をもと
のパルス列から除去し（ステップＳ７４）、残ったパル
ス列に対して再度自己相関関数を求めて（ステップＳ７
６）、上記の処理を繰り返す。
ＰＲＩの積分器から順番に分析していく（ステップＳ６
８）。積分器のテーブルのピークの値がしきい値を超え
ているか否かを調べ（ステップＳ７０）、超えていれば
そのＰＲＩを持つ信号源が存在すると判定する（ステッ
プＳ７２）。そのあとそのＰＲＩをもつパルス列をもと
のパルス列から除去し（ステップＳ７４）、残ったパル
ス列に対して再度自己相関関数を求めて（ステップＳ７
６）、上記の処理を繰り返す。
【０００８】一方、ステップＳ７０において、積分器の
テーブルのピークの値がしきい値を超えていなければ、
次に小さいＰＲＩに対する積分器の分析に進む（ステッ
プＳ７７）。これは分析する積分器が無くなるまで続け
られる（ステップＳ７８）。最後に、残ったパルスが単
独のものであるか否かを調べ、単独であれば処理を終了
する。単独のものでなくグループを構成していればレー
ダスキャンによって生じたものとみなし、スキャン周期
の分析を行う（ステップＳ７９）。
テーブルのピークの値がしきい値を超えていなければ、
次に小さいＰＲＩに対する積分器の分析に進む（ステッ
プＳ７７）。これは分析する積分器が無くなるまで続け
られる（ステップＳ７８）。最後に、残ったパルスが単
独のものであるか否かを調べ、単独であれば処理を終了
する。単独のものでなくグループを構成していればレー
ダスキャンによって生じたものとみなし、スキャン周期
の分析を行う（ステップＳ７９）。
【０００９】
【発明が解決しようとする課題】上記のような、従来の
パルス列分離装置では、まず自己相関関数でＰＲＩの候
補を検出することが処理の前提になっているが、信号源
の数が多い場合や、固定ＰＲＩ以外の複雑なＰＲＩのパ
ターンを有する信号源が存在する場合には、自己相関関
数に明確なピークが現われにくく、また、自己相関関数
にはＰＲＩの整数倍のところにできるハーモニックスの
問題があり、ＰＲＩの候補が適切に抽出できず、パルス
列を正確に検出し分離できないという問題点があった。
パルス列分離装置では、まず自己相関関数でＰＲＩの候
補を検出することが処理の前提になっているが、信号源
の数が多い場合や、固定ＰＲＩ以外の複雑なＰＲＩのパ
ターンを有する信号源が存在する場合には、自己相関関
数に明確なピークが現われにくく、また、自己相関関数
にはＰＲＩの整数倍のところにできるハーモニックスの
問題があり、ＰＲＩの候補が適切に抽出できず、パルス
列を正確に検出し分離できないという問題点があった。
【００１０】本発明は上記のような課題を解決するため
になされたもので、複雑なＰＲＩのパターンを有する信
号源も含めて信号源のＰＲＩを正確に推定し、パルス列
の分離精度を向上させることが可能なパルス列分離装置
及びパルス列分離方法を得ることを目的とする。
になされたもので、複雑なＰＲＩのパターンを有する信
号源も含めて信号源のＰＲＩを正確に推定し、パルス列
の分離精度を向上させることが可能なパルス列分離装置
及びパルス列分離方法を得ることを目的とする。
【００１１】
【課題を解決するための手段】この発明は、複数の信号
源からのパルス列が重なり合った複合パルス列の各パル
スの到着時間（ＴＯＡ）データに基づき、固定パルス繰
り返し間隔（ＰＲＩ）を有する固定ＰＲＩパルス列を検
出する固定ＰＲＩパルス列検出手段と、検出された固定
ＰＲＩパルス列をもとの複合パルス列から分離するとと
もに、分離後の残りのパルス列を出力する固定ＰＲＩパ
ルス列分離手段と、固定ＰＲＩパルス列分離手段から出
力される残りのパルス列から、複合パルス列の各パルス
の到着時間（ＴＯＡ）データに基づき、スタガーパルス
列を検出するスタガーパルス列検出手段と、検出された
スタガーパルス列を残りのパルス列から分離するととも
に、分離後の固定ＰＲＩパルス列及びスタガーパルス列
を含まない第二の残りのパルス列を出力するスタガーパ
ルス列分離手段と、スタガーパルス列分離手段から出力
される第二の残りのパルス列から、複合パルス列の各パ
ルスの到着時間（ＴＯＡ）データに基づき、ジッターパ
ルス列を検出するジッターパルス列検出手段と、検出さ
れたジッターパルス列を第二の残りのパルス列から分離
するジッターパルス列分離手段と、を備えたパルス列分
離装置である。
源からのパルス列が重なり合った複合パルス列の各パル
スの到着時間（ＴＯＡ）データに基づき、固定パルス繰
り返し間隔（ＰＲＩ）を有する固定ＰＲＩパルス列を検
出する固定ＰＲＩパルス列検出手段と、検出された固定
ＰＲＩパルス列をもとの複合パルス列から分離するとと
もに、分離後の残りのパルス列を出力する固定ＰＲＩパ
ルス列分離手段と、固定ＰＲＩパルス列分離手段から出
力される残りのパルス列から、複合パルス列の各パルス
の到着時間（ＴＯＡ）データに基づき、スタガーパルス
列を検出するスタガーパルス列検出手段と、検出された
スタガーパルス列を残りのパルス列から分離するととも
に、分離後の固定ＰＲＩパルス列及びスタガーパルス列
を含まない第二の残りのパルス列を出力するスタガーパ
ルス列分離手段と、スタガーパルス列分離手段から出力
される第二の残りのパルス列から、複合パルス列の各パ
ルスの到着時間（ＴＯＡ）データに基づき、ジッターパ
ルス列を検出するジッターパルス列検出手段と、検出さ
れたジッターパルス列を第二の残りのパルス列から分離
するジッターパルス列分離手段と、を備えたパルス列分
離装置である。
【００１２】また、固定ＰＲＩパルス列検出手段が、Ｐ
ＲＩ変換を用いて検出を行う。
ＲＩ変換を用いて検出を行う。
【００１３】また、スタガーパルス列検出手段が、等間
隔に並ぶ３つのパルスからなる３連パルスの集合に対し
てそのパルス繰り返し間隔のヒストグラムを求める３連
パルスヒストグラム演算部と、得られたヒストグラムか
ら複数のピークを検出し、各ピークがパルス列の基本周
期であるか否かを判定する判定部と、を備えている。
隔に並ぶ３つのパルスからなる３連パルスの集合に対し
てそのパルス繰り返し間隔のヒストグラムを求める３連
パルスヒストグラム演算部と、得られたヒストグラムか
ら複数のピークを検出し、各ピークがパルス列の基本周
期であるか否かを判定する判定部と、を備えている。
【００１４】また、ジッターパルス列検出手段が、オー
バーラップしたＰＲＩビンを持ち、時間原点の移動を行
う、変形ＰＲＩ変換を用いる。
バーラップしたＰＲＩビンを持ち、時間原点の移動を行
う、変形ＰＲＩ変換を用いる。
【００１５】また、この発明は、複数の信号源からのパ
ルス列が重なり合った複合パルス列の各パルスの到着時
間（ＴＯＡ）データに基づき、固定パルス繰り返し間隔
（ＰＲＩ）を有する固定ＰＲＩパルス列を検出する固定
ＰＲＩパルス列検出工程と、検出された固定ＰＲＩパル
ス列をもとの複合パルス列から分離するとともに、分離
後の残りのパルス列を出力する固定ＰＲＩパルス列分離
工程と、固定ＰＲＩパルス列分離手段から出力される残
りのパルス列から、複合パルス列の各パルスの到着時間
（ＴＯＡ）データに基づき、スタガーパルス列を検出す
るスタガーパルス列検出工程と、検出されたスタガーパ
ルス列を残りのパルス列から分離するとともに、分離後
の固定ＰＲＩパルス列及びスタガーパルス列を含まない
第二の残りのパルス列を出力するスタガーパルス列分離
工程と、スタガーパルス列分離手段から出力される第二
の残りのパルス列から、複合パルス列の各パルスの到着
時間（ＴＯＡ）データに基づき、ジッターパルス列を検
出するジッターパルス列検出工程と、検出されたジッタ
ーパルス列を第二の残りのパルス列から分離するジッタ
ーパルス列分離工程と、を備えたパルス列分離方法であ
る。
ルス列が重なり合った複合パルス列の各パルスの到着時
間（ＴＯＡ）データに基づき、固定パルス繰り返し間隔
（ＰＲＩ）を有する固定ＰＲＩパルス列を検出する固定
ＰＲＩパルス列検出工程と、検出された固定ＰＲＩパル
ス列をもとの複合パルス列から分離するとともに、分離
後の残りのパルス列を出力する固定ＰＲＩパルス列分離
工程と、固定ＰＲＩパルス列分離手段から出力される残
りのパルス列から、複合パルス列の各パルスの到着時間
（ＴＯＡ）データに基づき、スタガーパルス列を検出す
るスタガーパルス列検出工程と、検出されたスタガーパ
ルス列を残りのパルス列から分離するとともに、分離後
の固定ＰＲＩパルス列及びスタガーパルス列を含まない
第二の残りのパルス列を出力するスタガーパルス列分離
工程と、スタガーパルス列分離手段から出力される第二
の残りのパルス列から、複合パルス列の各パルスの到着
時間（ＴＯＡ）データに基づき、ジッターパルス列を検
出するジッターパルス列検出工程と、検出されたジッタ
ーパルス列を第二の残りのパルス列から分離するジッタ
ーパルス列分離工程と、を備えたパルス列分離方法であ
る。
【００１６】また、固定ＰＲＩパルス列検出工程におい
て、ＰＲＩ変換を用いて検出を行う。
て、ＰＲＩ変換を用いて検出を行う。
【００１７】また、スタガーパルス列検出工程が、等間
隔に並ぶ３つのパルスからなる３連パルスの集合に対し
てそのパルス繰り返し間隔のヒストグラムを求める３連
パルスヒストグラム演算ステップと、得られたヒストグ
ラムから複数のピークを検出し、各ピークがパルス列の
基本周期であるか否かを判定する判定ステップと、を備
えている。
隔に並ぶ３つのパルスからなる３連パルスの集合に対し
てそのパルス繰り返し間隔のヒストグラムを求める３連
パルスヒストグラム演算ステップと、得られたヒストグ
ラムから複数のピークを検出し、各ピークがパルス列の
基本周期であるか否かを判定する判定ステップと、を備
えている。
【００１８】また、ジッターパルス列検出工程におい
て、オーバーラップしたＰＲＩビンを持ち、時間原点の
移動を行う、変形ＰＲＩ変換を用いる。
て、オーバーラップしたＰＲＩビンを持ち、時間原点の
移動を行う、変形ＰＲＩ変換を用いる。
【００１９】
【発明の実施の形態】実施の形態１．以下、この発明の
一実施の形態を図１乃至図１５について説明する。
一実施の形態を図１乃至図１５について説明する。
【００２０】図１は本発明のパルス列分離装置が処理の
対象とするパルス列の種類を説明するための説明図であ
り、図２は本発明のパルス列分離装置の機能ブロック
図、図３は固定ＰＲＩパルス列検出・分離手段のブロッ
ク図、図４は固定ＰＲＩパルス列検出・分離手段の中で
用いるＰＲＩ変換に関するＰＲＩビンの説明図、図５は
固定ＰＲＩパルス列検出・分離手段の中で用いるＰＲＩ
変換のフローチャート、図６はスタガーパルス列検出・
分離手段のブロック図、図７はスタガーパルス列検出・
分離手段の中で用いる３連パルスのヒストグラムのフロ
ーチャート、図８はスタガーパルス列検出・分離手段の
中で用いる３連パルスのヒストグラムからの基本周期検
出のフローチャート、図９はジッターパルス列検出・分
離手段のブロック図、図１０はジッターパルス列検出・
分離手段の中で用いる変形ＰＲＩ変換に関するＰＲＩビ
ンの説明図、図１１はジッターパルス列検出・分離手段
の中で用いる変形 ＰＲＩ変換のフローチャート、図１
２は入力パルス列の説明図、図１３は固定ＰＲＩパルス
列検出・分離の説明図、図１４はスタガーパルス列検出
・分離の説明図、図１５はジッターパルス列検出・分離
の説明図である。
対象とするパルス列の種類を説明するための説明図であ
り、図２は本発明のパルス列分離装置の機能ブロック
図、図３は固定ＰＲＩパルス列検出・分離手段のブロッ
ク図、図４は固定ＰＲＩパルス列検出・分離手段の中で
用いるＰＲＩ変換に関するＰＲＩビンの説明図、図５は
固定ＰＲＩパルス列検出・分離手段の中で用いるＰＲＩ
変換のフローチャート、図６はスタガーパルス列検出・
分離手段のブロック図、図７はスタガーパルス列検出・
分離手段の中で用いる３連パルスのヒストグラムのフロ
ーチャート、図８はスタガーパルス列検出・分離手段の
中で用いる３連パルスのヒストグラムからの基本周期検
出のフローチャート、図９はジッターパルス列検出・分
離手段のブロック図、図１０はジッターパルス列検出・
分離手段の中で用いる変形ＰＲＩ変換に関するＰＲＩビ
ンの説明図、図１１はジッターパルス列検出・分離手段
の中で用いる変形 ＰＲＩ変換のフローチャート、図１
２は入力パルス列の説明図、図１３は固定ＰＲＩパルス
列検出・分離の説明図、図１４はスタガーパルス列検出
・分離の説明図、図１５はジッターパルス列検出・分離
の説明図である。
【００２１】まず、本発明のパルス列分離装置が対象と
するパルス列の種類を図１を用いて説明する。図１にお
いて、１１はＰＲＩが一定の固定ＰＲＩパルス列であ
る。パルス列には固定ＰＲＩパルス列１１以外に、１２
に示すスタガーパルス列と１３に示すジッターパルス列
がある。スタガーパルス列１２は、一定のフレーム周期
１２ａの中に、複数の不等間隔（図１では、ＰＲＩ１，
ＰＲＩ２，ＰＲＩ３の３つ）のパルスが存在し、それが
周期１２ａごとに繰り返すパルス列である。ジッターパ
ルス列１３は、ＰＲＩが平均ＰＲＩのまわりにランダム
に揺らぐパルス列である。本発明のパルス列分離装置へ
の入力パルス列は、固定ＰＲＩパルス列１１、スタガー
パルス列１２、ジッターパルス列１３がそれぞれ複数重
なり合った複合パルス列である。
するパルス列の種類を図１を用いて説明する。図１にお
いて、１１はＰＲＩが一定の固定ＰＲＩパルス列であ
る。パルス列には固定ＰＲＩパルス列１１以外に、１２
に示すスタガーパルス列と１３に示すジッターパルス列
がある。スタガーパルス列１２は、一定のフレーム周期
１２ａの中に、複数の不等間隔（図１では、ＰＲＩ１，
ＰＲＩ２，ＰＲＩ３の３つ）のパルスが存在し、それが
周期１２ａごとに繰り返すパルス列である。ジッターパ
ルス列１３は、ＰＲＩが平均ＰＲＩのまわりにランダム
に揺らぐパルス列である。本発明のパルス列分離装置へ
の入力パルス列は、固定ＰＲＩパルス列１１、スタガー
パルス列１２、ジッターパルス列１３がそれぞれ複数重
なり合った複合パルス列である。
【００２２】次に、本発明の構成及び動作について説明
する。
する。
【００２３】図２において、１００は、信号源から送信
されたパルスを受信し、当該パルスのＴＯＡデータを出
力する広帯域受信機、２００は本発明のパルス列分離装
置に設けられた入力端、３００は固定ＰＲＩパルス列の
検出・分離を行うための固定ＰＲＩパルス列検出・分離
手段、４００は、スタガーパルス列の検出・分離を行う
ためのスタガーパルス列検出・分離手段、５００は、ジ
ッターパルス列の検出・分離を行うためのジッターパル
ス列検出・分離手段、６００は本発明のパルス列分離装
置に設けられた出力端である。
されたパルスを受信し、当該パルスのＴＯＡデータを出
力する広帯域受信機、２００は本発明のパルス列分離装
置に設けられた入力端、３００は固定ＰＲＩパルス列の
検出・分離を行うための固定ＰＲＩパルス列検出・分離
手段、４００は、スタガーパルス列の検出・分離を行う
ためのスタガーパルス列検出・分離手段、５００は、ジ
ッターパルス列の検出・分離を行うためのジッターパル
ス列検出・分離手段、６００は本発明のパルス列分離装
置に設けられた出力端である。
【００２４】動作について説明する。広帯域受信機１０
０からのパルスのＴＯＡデータは入力端２００へ入力さ
れる。これらＴＯＡデータはまず固定ＰＲＩパルス列検
出・分離手段３００で検出・分離処理が行われる。
０からのパルスのＴＯＡデータは入力端２００へ入力さ
れる。これらＴＯＡデータはまず固定ＰＲＩパルス列検
出・分離手段３００で検出・分離処理が行われる。
【００２５】図３には、図２の固定ＰＲＩパルス列検出
・分離手段３００の内部構成が示されている。図におい
て、３１０は入力端、３３０はＰＲＩ変換を行うＰＲＩ
変換演算部、３５０は適当なしきい値の値を設定すると
ともに検出可能範囲を決定するしきい値演算部、３６０
は固定ＰＲＩパルス列のＰＲＩを抽出するＰＲＩ抽出
部、３７０は抽出されたＰＲＩの値を利用して、もとの
パルス列から該当する部分パルス列を分離する固定ＰＲ
Ｉパルス列分離演算部、３８０は固定ＰＲＩパルス列が
残存するか否かを判定する残存パルス列判定部、３９０
は出力端である。
・分離手段３００の内部構成が示されている。図におい
て、３１０は入力端、３３０はＰＲＩ変換を行うＰＲＩ
変換演算部、３５０は適当なしきい値の値を設定すると
ともに検出可能範囲を決定するしきい値演算部、３６０
は固定ＰＲＩパルス列のＰＲＩを抽出するＰＲＩ抽出
部、３７０は抽出されたＰＲＩの値を利用して、もとの
パルス列から該当する部分パルス列を分離する固定ＰＲ
Ｉパルス列分離演算部、３８０は固定ＰＲＩパルス列が
残存するか否かを判定する残存パルス列判定部、３９０
は出力端である。
【００２６】動作について説明する。まず、入力端３１
０へ入力されたパルス列のＴＯＡデータにはＰＲＩ変換
演算部３３０においてＰＲＩ変換が施される。このＰＲ
Ｉ変換は、例えば、西口（三菱電機（株））、''パルス
繰り返し間隔推定のための新手法''、昭和５８年度電子
情報通信学会情報・システム部門全国大会、に示されて
いる方法を用いる。この方法は、まず、入力データのＴ
ＯＡを t1、t2、・・・、tN とし、入力パルス列をデル
タ関数の和で次式のように表わす。
０へ入力されたパルス列のＴＯＡデータにはＰＲＩ変換
演算部３３０においてＰＲＩ変換が施される。このＰＲ
Ｉ変換は、例えば、西口（三菱電機（株））、''パルス
繰り返し間隔推定のための新手法''、昭和５８年度電子
情報通信学会情報・システム部門全国大会、に示されて
いる方法を用いる。この方法は、まず、入力データのＴ
ＯＡを t1、t2、・・・、tN とし、入力パルス列をデル
タ関数の和で次式のように表わす。
【００２７】
【数１】
【００２８】ＰＲＩ変換は関数 g(t) に対して次式で定
義されるような変換である。
義されるような変換である。
【００２９】
【数２】
【００３０】これは上式（１）のパルス列に対しては以
下のような形になる。
下のような形になる。
【００３１】
【数３】
【００３２】このＰＲＩ変換は、連続時間での表現であ
るが、実際の計算機で行う離散時間でのＰＲＩ変換は後
でフローチャートを用いて詳述する。このＰＲＩ変換を
入力パルス列に施すと、入力パルス列に含まれている部
分パルス列のＰＲＩのところにだけ鋭いピークがたつよ
うな変換結果が得られる。さらに、しきい値演算部３５
０でＰＲＩ検出のための適当なしきい値を定める。しき
い値演算部３５０では、同時に検出可能なＰＲＩの範囲
も求める。次に、ＰＲＩ抽出部３６０で、そのしきい値
を越えるピーク値を検出するようにすれば、固定ＰＲＩ
パルス列の検出可能な範囲にあるＰＲＩがすべて検出さ
れる。このようにして検出されたＰＲＩの値を利用し
て、固定ＰＲＩパルス列分離演算部３７０で、もとのパ
ルス列から該当する部分パルス列を分離する。この分離
のためには、例えば、特公平１−４７９３６号公報に詳
しく示されているパルス繰り返し間隔フィルターを用い
る。残存パルス列判定部３８０では、調べたいＰＲＩの
範囲と検出可能なＰＲＩの範囲を比較し、もし検出可能
なＰＲＩの範囲の方が狭ければ、すなわち、固定ＰＲＩ
パルス列が残存する可能性があれば、残ったパルス列に
再びＰＲＩ変換演算部３３０でＰＲＩ変換を施し処理を
繰り返す。そうでなければ、残ったパルス列と共に出力
端３９０へ進む。
るが、実際の計算機で行う離散時間でのＰＲＩ変換は後
でフローチャートを用いて詳述する。このＰＲＩ変換を
入力パルス列に施すと、入力パルス列に含まれている部
分パルス列のＰＲＩのところにだけ鋭いピークがたつよ
うな変換結果が得られる。さらに、しきい値演算部３５
０でＰＲＩ検出のための適当なしきい値を定める。しき
い値演算部３５０では、同時に検出可能なＰＲＩの範囲
も求める。次に、ＰＲＩ抽出部３６０で、そのしきい値
を越えるピーク値を検出するようにすれば、固定ＰＲＩ
パルス列の検出可能な範囲にあるＰＲＩがすべて検出さ
れる。このようにして検出されたＰＲＩの値を利用し
て、固定ＰＲＩパルス列分離演算部３７０で、もとのパ
ルス列から該当する部分パルス列を分離する。この分離
のためには、例えば、特公平１−４７９３６号公報に詳
しく示されているパルス繰り返し間隔フィルターを用い
る。残存パルス列判定部３８０では、調べたいＰＲＩの
範囲と検出可能なＰＲＩの範囲を比較し、もし検出可能
なＰＲＩの範囲の方が狭ければ、すなわち、固定ＰＲＩ
パルス列が残存する可能性があれば、残ったパルス列に
再びＰＲＩ変換演算部３３０でＰＲＩ変換を施し処理を
繰り返す。そうでなければ、残ったパルス列と共に出力
端３９０へ進む。
【００３３】次に、実際の計算機で行う離散時間でのＰ
ＲＩ変換について説明する。図４は、ＰＲＩ変換で用い
るＰＲＩビンの説明図である。調べたいＰＲＩの下限と
上限をそれぞれ、τmin、τmax と置く。ＰＲＩ軸での
区間[τmin，τmax] を K 個に等分割し、分割された細
分区間をＰＲＩビンと呼ぶ。各ビンの幅は次式（４）で
与えられる。
ＲＩ変換について説明する。図４は、ＰＲＩ変換で用い
るＰＲＩビンの説明図である。調べたいＰＲＩの下限と
上限をそれぞれ、τmin、τmax と置く。ＰＲＩ軸での
区間[τmin，τmax] を K 個に等分割し、分割された細
分区間をＰＲＩビンと呼ぶ。各ビンの幅は次式（４）で
与えられる。
【００３４】
【数４】
【００３５】また、中心ＰＲＩは次式（５）で与えられ
る。
る。
【００３６】
【数５】
【００３７】図５にはＰＲＩ変換の処理のフローチャー
トを示す。入力データのパルスのＴＯＡを t1、 t2、
…、 tN とする。ＰＲＩ変換の処理は以下の手順で行
う。
トを示す。入力データのパルスのＴＯＡを t1、 t2、
…、 tN とする。ＰＲＩ変換の処理は以下の手順で行
う。
【００３８】（１）初期設定を行う（ステップＳ３３
１）。すなわち、ＰＲＩ変換 Dk（k =1， 2， … ，
K）をゼロクリアするとともに、パルスペアの右パルス
のパルス番号を m = 1 に設定する。
１）。すなわち、ＰＲＩ変換 Dk（k =1， 2， … ，
K）をゼロクリアするとともに、パルスペアの右パルス
のパルス番号を m = 1 に設定する。
【００３９】（２）パルスペアの右パルスのパルス番号
を更新する（ステップＳ３３２）。すなわち、m = m +
1を行う。
を更新する（ステップＳ３３２）。すなわち、m = m +
1を行う。
【００４０】（３）終了判定を行う（ステップＳ３３
３）。すなわち、m≦ NならばステップＳ３３４へ。そ
うでなければ終了する。
３）。すなわち、m≦ NならばステップＳ３３４へ。そ
うでなければ終了する。
【００４１】（４）パルスペアの左パルスのパルス番号
をn ＝ m に設定する（ステップＳ３３４）。
をn ＝ m に設定する（ステップＳ３３４）。
【００４２】（５）パルスペアの左パルスのパルス番号
を更新する（ステップＳ３３５）。すなわち、n ＝n−1
を行う。
を更新する（ステップＳ３３５）。すなわち、n ＝n−1
を行う。
【００４３】（６）パルスペアの左パルスの終了判定を
行う（ステップＳ３３６）。すなわち、n ≧ 1 なら
ば、ステップＳ３３７へ。そうでなければステップＳ３
３２へ。
行う（ステップＳ３３６）。すなわち、n ≧ 1 なら
ば、ステップＳ３３７へ。そうでなければステップＳ３
３２へ。
【００４４】（７）パルスペアの時間間隔を求める（ス
テップＳ３３７）。すなわち、τ=tm−tnを行う。
テップＳ３３７）。すなわち、τ=tm−tnを行う。
【００４５】（８）時間間隔の大きさ判定１を行う（ス
テップＳ３３８）。すなわち、τ≧τmin ならばステッ
プＳ３３９へ。そうでなければステップＳ３３５へ。
テップＳ３３８）。すなわち、τ≧τmin ならばステッ
プＳ３３９へ。そうでなければステップＳ３３５へ。
【００４６】（９）時間間隔の大きさ判定２を行う（ス
テップＳ３３９）。すなわち、τ≦τmax ならばステッ
プＳ３４０へ。そうでなければＳ３３２へ。
テップＳ３３９）。すなわち、τ≦τmax ならばステッ
プＳ３４０へ。そうでなければＳ３３２へ。
【００４７】（１０）τ軸上でのＰＲＩビンの選出（ス
テップＳ３４０）。すなわち、ＰＲＩビンの番号kを次
式（６）で計算する。
テップＳ３４０）。すなわち、ＰＲＩビンの番号kを次
式（６）で計算する。
【００４８】
【数６】
【００４９】ここで、[・] は Gaussの整数、すなわ
ち、(・) を越えない最大の整数を表わす。
ち、(・) を越えない最大の整数を表わす。
【００５０】（１１）ＰＲＩ変換の更新を行う（ステッ
プＳ３４１）。すなわち、次式（７）でＰＲＩ変換の値
を更新する。このあと、ステップＳ３３５へ戻る。
プＳ３４１）。すなわち、次式（７）でＰＲＩ変換の値
を更新する。このあと、ステップＳ３３５へ戻る。
【００５１】
【数７】
【００５２】図３のしきい値演算部３５０で求めるしき
い値は、例えば次のような式にする。
い値は、例えば次のような式にする。
【００５３】
【数８】
【００５４】ここで、T は観測時間で、最初と最後のパ
ルスのＴＯＡの差から求める。αは調整可能なパラメー
タで、例えば 0．3 くらいの値にする。上のしきい値の
式（８）は、もし、全観測時間にわたってＰＲＩがτk
のパルス列が存在すればＰＲＩ変換の価は理想的には、
次式（９）となることに基づく。
ルスのＴＯＡの差から求める。αは調整可能なパラメー
タで、例えば 0．3 くらいの値にする。上のしきい値の
式（８）は、もし、全観測時間にわたってＰＲＩがτk
のパルス列が存在すればＰＲＩ変換の価は理想的には、
次式（９）となることに基づく。
【００５５】
【数９】
【００５６】しきい値はτkの減少関数であるが、一
方、ＰＲＩ変換のノイズフロアーはτkと共に減少はし
ない。したがって、τk がある大きさ（これをτlim と
おく）になるとそこからはしきい値よりもノイズフロア
ーの値の方が大きな値をとってＰＲＩの検出は不能にな
る。すなわち、[τmin，τlim] が検出可能範囲にな
る。
方、ＰＲＩ変換のノイズフロアーはτkと共に減少はし
ない。したがって、τk がある大きさ（これをτlim と
おく）になるとそこからはしきい値よりもノイズフロア
ーの値の方が大きな値をとってＰＲＩの検出は不能にな
る。すなわち、[τmin，τlim] が検出可能範囲にな
る。
【００５７】以上のようにして図２の固定ＰＲＩパルス
列検出・分離手段３００で固定ＰＲＩパルス列がもとの
パルス列から除去されると、次に、残ったパルス列に対
して、スタガーパルス列検出・分離手段４００でスタガ
ーパルス列の検出・分離処理が行われる。
列検出・分離手段３００で固定ＰＲＩパルス列がもとの
パルス列から除去されると、次に、残ったパルス列に対
して、スタガーパルス列検出・分離手段４００でスタガ
ーパルス列の検出・分離処理が行われる。
【００５８】図６にスタガーパルス列検出・分離手段の
内部構成を示したブロック図を示す。図において、４１
０は入力端、４３０は３連パルスのヒストグラムを求め
る３連パルスヒストグラム演算部、４５０は適当なしき
い値の値を設定するとともに検出可能範囲を決定するし
きい値演算部、４６０はスタガーパルス列のＰＲＩを抽
出するＰＲＩ抽出部、４７０は抽出されたＰＲＩの値を
利用して、もとのパルス列から該当する部分パルス列を
分離するスタガーパルス列分離演算部、４８０はスタガ
ーパルス列が残存するか否かを判定する残存パルス列判
定部、４９０は出力端である。
内部構成を示したブロック図を示す。図において、４１
０は入力端、４３０は３連パルスのヒストグラムを求め
る３連パルスヒストグラム演算部、４５０は適当なしき
い値の値を設定するとともに検出可能範囲を決定するし
きい値演算部、４６０はスタガーパルス列のＰＲＩを抽
出するＰＲＩ抽出部、４７０は抽出されたＰＲＩの値を
利用して、もとのパルス列から該当する部分パルス列を
分離するスタガーパルス列分離演算部、４８０はスタガ
ーパルス列が残存するか否かを判定する残存パルス列判
定部、４９０は出力端である。
【００５９】なお、しきい値演算部３５０及びＰＲＩ抽
出部３６０は、３連パルスヒストグラム演算部３３０に
おいて得られたヒストグラムから複数のピークを検出
し、各ピークがパルス列の基本周期であるか否かを判定
するための判定手段を構成している。
出部３６０は、３連パルスヒストグラム演算部３３０に
おいて得られたヒストグラムから複数のピークを検出
し、各ピークがパルス列の基本周期であるか否かを判定
するための判定手段を構成している。
【００６０】次に、動作について説明する。入力端４１
０に入力したパルス列は、３連パルスヒストグラム演算
部４３０で３連パルスのヒストグラムが計算される。３
連パルスのヒストグラムは、図１で説明したようなほぼ
等間隔のパルスの３つ組の個数をカウントするヒストグ
ラムである。このヒストグラムを作成する処理はあとで
詳述するが、得られるヒストグラムでは、スタガーパル
ス列の基本周期（フレーム周期１２ａ）とその整数倍の
ところにするどいピークがたつ。しきい値演算部４５０
では スタガーのＰＲＩ（基本周期）検出のための適当
なしきい値を定める。しきい値演算部４５０では同時
に、検出可能なＰＲＩの範囲も求める。次に、ＰＲＩ抽
出部４６０で、そのしきい値を超える値を検出するよう
にすれば、スタガーパルス列の検出可能な範囲にあるＰ
ＲＩがすべて検出される。ＰＲＩ抽出部４６０で行う３
連パルスのヒストグラムからのＰＲＩ（基本周期）の検
出についてはあとでフローチャートを用いて詳述する。
このようにして検出されたＰＲＩの値を利用して、スタ
ガーパルス列分離演算部４７０で、もとのパルス列から
該当する部分パルス列を分離する。この分離のために
は、例えば特公平１−４７９３６号公報に詳しく示され
ているパルス繰り返し間隔フィルターを用いる。残存パ
ルス列判定部４８０では、調べたいＰＲＩの範囲と検出
可能なＰＲＩの範囲を比較し、もし検出可能なＰＲＩの
範囲の方が狭ければ、すなわち、スタガーパルス列が残
存する可能性があれば、残ったパルス列に再び３連パル
スヒストグラム演算部４３０でヒストグラムを作り、処
理を繰り返す。そうでなければ、残ったパルス列と共に
出力端４９０へ進む。
０に入力したパルス列は、３連パルスヒストグラム演算
部４３０で３連パルスのヒストグラムが計算される。３
連パルスのヒストグラムは、図１で説明したようなほぼ
等間隔のパルスの３つ組の個数をカウントするヒストグ
ラムである。このヒストグラムを作成する処理はあとで
詳述するが、得られるヒストグラムでは、スタガーパル
ス列の基本周期（フレーム周期１２ａ）とその整数倍の
ところにするどいピークがたつ。しきい値演算部４５０
では スタガーのＰＲＩ（基本周期）検出のための適当
なしきい値を定める。しきい値演算部４５０では同時
に、検出可能なＰＲＩの範囲も求める。次に、ＰＲＩ抽
出部４６０で、そのしきい値を超える値を検出するよう
にすれば、スタガーパルス列の検出可能な範囲にあるＰ
ＲＩがすべて検出される。ＰＲＩ抽出部４６０で行う３
連パルスのヒストグラムからのＰＲＩ（基本周期）の検
出についてはあとでフローチャートを用いて詳述する。
このようにして検出されたＰＲＩの値を利用して、スタ
ガーパルス列分離演算部４７０で、もとのパルス列から
該当する部分パルス列を分離する。この分離のために
は、例えば特公平１−４７９３６号公報に詳しく示され
ているパルス繰り返し間隔フィルターを用いる。残存パ
ルス列判定部４８０では、調べたいＰＲＩの範囲と検出
可能なＰＲＩの範囲を比較し、もし検出可能なＰＲＩの
範囲の方が狭ければ、すなわち、スタガーパルス列が残
存する可能性があれば、残ったパルス列に再び３連パル
スヒストグラム演算部４３０でヒストグラムを作り、処
理を繰り返す。そうでなければ、残ったパルス列と共に
出力端４９０へ進む。
【００６１】図７に、３連パルスヒストグラム演算部４
３０で行う、３連パルスのヒストグラム演算の処理の流
れを示すフローチャートを示す。このヒストグラムで用
いるＰＲＩビンはＰＲＩ変換の場合の図４と同一であ
る。入力データのパルスのＴＯＡを t1、 t2、 …、 tN
とする。３連パルスのヒストグラムの作成は以下の手
順で行う。
３０で行う、３連パルスのヒストグラム演算の処理の流
れを示すフローチャートを示す。このヒストグラムで用
いるＰＲＩビンはＰＲＩ変換の場合の図４と同一であ
る。入力データのパルスのＴＯＡを t1、 t2、 …、 tN
とする。３連パルスのヒストグラムの作成は以下の手
順で行う。
【００６２】（１）まずはじめに、初期設定（ステップ
Ｓ４３１）を行う。すなわち、ヒストグラム Ek（ｋ=
1， 2， … ， K）をゼロクリアする。また、３連パル
スの中心パルスのパルス番号を m = 1 に設定する。
Ｓ４３１）を行う。すなわち、ヒストグラム Ek（ｋ=
1， 2， … ， K）をゼロクリアする。また、３連パル
スの中心パルスのパルス番号を m = 1 に設定する。
【００６３】（２）３連パルスの中心パルスのパルス番
号を更新する（ステップＳ４３２）。すなわち、m = m
+ 1を行う。
号を更新する（ステップＳ４３２）。すなわち、m = m
+ 1を行う。
【００６４】（３）終了判定を行う（ステップＳ４３
３）。すなわち、m≦N−1ならばステップＳ４３４へ。
そうでなければ終了。
３）。すなわち、m≦N−1ならばステップＳ４３４へ。
そうでなければ終了。
【００６５】（４）３連パルスの右パルスのパルス番号
を F = m + 1 に、左パルスのパルス番号をB = m に設
定する（ステップＳ４３４）。
を F = m + 1 に、左パルスのパルス番号をB = m に設
定する（ステップＳ４３４）。
【００６６】（５）３連パルスの左パルスのパルス番号
を更新する（ステップＳ４３５）。すなわち、B = B−1
を行う。
を更新する（ステップＳ４３５）。すなわち、B = B−1
を行う。
【００６７】（６）３連パルスの左パルスの終了判定を
行う（ステップＳ４３６）。すなわち、B ≧1 ならば、
ステップＳ４３７へ。そうでなければステップＳ４３２
へ。
行う（ステップＳ４３６）。すなわち、B ≧1 ならば、
ステップＳ４３７へ。そうでなければステップＳ４３２
へ。
【００６８】（７）３連パルスの左ペア（後ろ向きパル
スペア）の時間間隔を求める（ステップＳ４３７）。す
なわち、τ= tm−tB。
スペア）の時間間隔を求める（ステップＳ４３７）。す
なわち、τ= tm−tB。
【００６９】（８）時間間隔の大きさ判定１（ステップ
Ｓ４３８）を行う。すなわち、τ≧τmin ならばステッ
プＳ４３９へ。そうでなければステップＳ４３５へ。
Ｓ４３８）を行う。すなわち、τ≧τmin ならばステッ
プＳ４３９へ。そうでなければステップＳ４３５へ。
【００７０】（９）時間間隔の大きさ判定２（ステップ
Ｓ４３９）を行う。すなわち、τ≦τmax ならばステッ
プＳ４４０へ。そうでなければＳ４３２へ。
Ｓ４３９）を行う。すなわち、τ≦τmax ならばステッ
プＳ４４０へ。そうでなければＳ４３２へ。
【００７１】（１０）３連パルスの右ペア（前向きパル
スペア）の時間間隔を求める（ステップＳ４４０）。す
なわち、τ’= tF−tmを行う。
スペア）の時間間隔を求める（ステップＳ４４０）。す
なわち、τ’= tF−tmを行う。
【００７２】（１１）次に、時間間隔の大きさ判定３
（ステップＳ４４１）を行う。すなわち、τ’≧τ−Δ
τ/2ならばステップＳ４４５へ。そうでなければステッ
プＳ４４２へ。
（ステップＳ４４１）を行う。すなわち、τ’≧τ−Δ
τ/2ならばステップＳ４４５へ。そうでなければステッ
プＳ４４２へ。
【００７３】（１２）３連パルスの右パルスの終了判定
（ステップＳ４４２）を行う。すなわち、F = N ならば
ステップＳ４３５へ戻る。そうでなければ右パルスを F
= F+ 1 と更新（ステップＳ４４３）してステップＳ４
４０へ。
（ステップＳ４４２）を行う。すなわち、F = N ならば
ステップＳ４３５へ戻る。そうでなければ右パルスを F
= F+ 1 と更新（ステップＳ４４３）してステップＳ４
４０へ。
【００７４】（１３）次に、時間間隔の大きさ判定４
（ステップＳ４４５）を行う。すなわち、τ’≦τ＋Δ
τ/2ならばステップＳ４４６へ。そうでなければＳ４３
５へ。
（ステップＳ４４５）を行う。すなわち、τ’≦τ＋Δ
τ/2ならばステップＳ４４６へ。そうでなければＳ４３
５へ。
【００７５】（１４）τ軸上でのＰＲＩビンの選出を行
う（ステップＳ４４６）。すなわち、ＰＲＩビンの番号
k を次式（１０）で計算する。
う（ステップＳ４４６）。すなわち、ＰＲＩビンの番号
k を次式（１０）で計算する。
【００７６】
【数１０】
【００７７】ここで、[・] は Gauss の整数、すなわ
ち、(・) を越えない最大の整数を表わす。
ち、(・) を越えない最大の整数を表わす。
【００７８】（１５）ヒストグラムの更新を行う（ステ
ップＳ４４７）。すなわち、次式（１１）で３連パルス
のヒストグラムの値を更新する。そのあと、ステップＳ
４３５へ。
ップＳ４４７）。すなわち、次式（１１）で３連パルス
のヒストグラムの値を更新する。そのあと、ステップＳ
４３５へ。
【００７９】
【数１１】
【００８０】次に、図８に、上述のようにして求めた３
連パルスのヒストグラムからＰＲＩ（基本周期）を検出
する処理のフローチャートを示す。検出は以下のような
手順で行う。
連パルスのヒストグラムからＰＲＩ（基本周期）を検出
する処理のフローチャートを示す。検出は以下のような
手順で行う。
【００８１】（１）まず、しきい値の初期設定を行う
（ステップＳ４６１）。すなわち、次式（１２）でしき
い値を設定する。
（ステップＳ４６１）。すなわち、次式（１２）でしき
い値を設定する。
【００８２】
【数１２】
【００８３】ここで、αはパラメータで、1．6 程度の
値を用いる。また、T は観測時間であり、最初のパルス
と最後のパルスのＴＯＡの差から求める。
値を用いる。また、T は観測時間であり、最初のパルス
と最後のパルスのＴＯＡの差から求める。
【００８４】（２）ピークの検出を行う（ステップＳ４
６２）。すなわち、Ek = 1， …，K のうち、Ek ≧A
k 、かつ 、Ek ≧max {Ek-1，Ek+1}（局所最大値）とな
るものを見つけ、それらを Ekp（p = 1， …， pmax）
とおく。
６２）。すなわち、Ek = 1， …，K のうち、Ek ≧A
k 、かつ 、Ek ≧max {Ek-1，Ek+1}（局所最大値）とな
るものを見つけ、それらを Ekp（p = 1， …， pmax）
とおく。
【００８５】（３）ピーク番号の初期設定を行う（ステ
ップＳ４６３）。すなわち、 p = 0とする。
ップＳ４６３）。すなわち、 p = 0とする。
【００８６】（４）ピーク番号の更新を行う（ステップ
Ｓ４６４）。すなわち、p = p + 1とする。
Ｓ４６４）。すなわち、p = p + 1とする。
【００８７】（５）終了判定を行う（ステップＳ４６
５）。すなわち、p ≦ pmax ならばステップＳ４６６
へ。そうでなければ終了。
５）。すなわち、p ≦ pmax ならばステップＳ４６６
へ。そうでなければ終了。
【００８８】（６）しきい値の更新を行う（ステップＳ
４６６）。すなわち、次式（１３）でしきい値 Ak の値
を更新する。
４６６）。すなわち、次式（１３）でしきい値 Ak の値
を更新する。
【００８９】
【数１３】
【００９０】（７）ピークの判定を行う（ステップＳ４
６７）。すなわち、p番目のピーク値 Ekp が、更新され
たしきい値Akp を上回っているか否かを判定する。Ekp
≧Akpならば、ステップＳ４６８へ。そうでなければス
テップＳ４６４へ。
６７）。すなわち、p番目のピーク値 Ekp が、更新され
たしきい値Akp を上回っているか否かを判定する。Ekp
≧Akpならば、ステップＳ４６８へ。そうでなければス
テップＳ４６４へ。
【００９１】（８）次に、基本周期か否かの判定を行う
（ステップＳ４６８）。すなわち、すでに検出されてい
る基本周期のすべてについて、τkp がその整数倍にな
っているか否かの判定を行う。ある基本周期のほぼ整数
倍（ｊ倍）になっていれば、τをτkp/j で置き換える
ことにより、基本周期の精度を上げる。もし、どの基本
周期の整数倍にもなっていなければ、τkp を新たな基
本周期として付け加える。このあとステップＳ４６４へ
戻る。
（ステップＳ４６８）。すなわち、すでに検出されてい
る基本周期のすべてについて、τkp がその整数倍にな
っているか否かの判定を行う。ある基本周期のほぼ整数
倍（ｊ倍）になっていれば、τをτkp/j で置き換える
ことにより、基本周期の精度を上げる。もし、どの基本
周期の整数倍にもなっていなければ、τkp を新たな基
本周期として付け加える。このあとステップＳ４６４へ
戻る。
【００９２】図２のスタガーパルス列検出・分離手段４
００でスタガーパルス列がもとのパルス列から除去され
ると、残ったパルス列に対してジッターパルス列検出・
分離手段５００でジッターパルス列の検出・分離処理を
行う。
００でスタガーパルス列がもとのパルス列から除去され
ると、残ったパルス列に対してジッターパルス列検出・
分離手段５００でジッターパルス列の検出・分離処理を
行う。
【００９３】図９にジッターパルス列検出・分離手段５
００の内部構成を示したブロック図を示す。図におい
て、５１０は入力端、５３０は変形ＰＲＩ変換を行う変
形ＰＲＩ変換演算部、５５０は適当なしきい値の値を設
定するとともに検出可能範囲を決定するしきい値演算
部、５６０はジッターパルス列のＰＲＩを抽出するＰＲ
Ｉ抽出部、５７０は抽出されたＰＲＩの値を利用して、
もとのパルス列から該当する部分パルス列を分離するジ
ッターパルス列分離演算部、５８０はジッターパルス列
が残存するか否かを判定する残存パルス列判定部、５９
０は出力端である。
００の内部構成を示したブロック図を示す。図におい
て、５１０は入力端、５３０は変形ＰＲＩ変換を行う変
形ＰＲＩ変換演算部、５５０は適当なしきい値の値を設
定するとともに検出可能範囲を決定するしきい値演算
部、５６０はジッターパルス列のＰＲＩを抽出するＰＲ
Ｉ抽出部、５７０は抽出されたＰＲＩの値を利用して、
もとのパルス列から該当する部分パルス列を分離するジ
ッターパルス列分離演算部、５８０はジッターパルス列
が残存するか否かを判定する残存パルス列判定部、５９
０は出力端である。
【００９４】動作について説明する。入力端５１０に入
力したパルス列は、変形ＰＲＩ変換演算部５３０で変形
ＰＲＩ変換が計算される。変形ＰＲＩ変換についてはあ
とで詳述するが、これから求められるＰＲＩスペクトラ
ムには、ジッターパルス列の平均ＰＲＩのところにする
どいピークがたつ。しきい値演算部５５０ではジッター
パルス列のＰＲＩ検出のための適当なしきい値を定め
る。しきい値演算部５５０では同時に、検出可能なＰＲ
Ｉの範囲も求める。次に、ＰＲＩ抽出部５６０で、その
しきい値を超える値を検出するようにすれば、ジッター
パルス列の検出可能な範囲にあるＰＲＩがすべて検出さ
れる。このようにして検出されたＰＲＩの値を利用し
て、ジッターパルス列分離演算部５７０で、もとのパル
ス列から該当する部分パルス列を分離する。この分離の
ためには、例えば特公平１−４７９３６号公報に詳しく
示されているパルス繰り返し間隔フィルターを用いる。
残存パルス列判定部５８０では、調べたいＰＲＩの範囲
と検出可能なＰＲＩの範囲を比較し、もし検出可能なＰ
ＲＩの範囲の方が狭ければ、すなわち、ジッターパルス
列が残存する可能性があれば、残ったパルス列に再び変
形ＰＲＩ変換演算部５３０でＰＲＩスペクトルを作り、
処理を繰り返す。そうでなければ、出力端５９０へ進
む。
力したパルス列は、変形ＰＲＩ変換演算部５３０で変形
ＰＲＩ変換が計算される。変形ＰＲＩ変換についてはあ
とで詳述するが、これから求められるＰＲＩスペクトラ
ムには、ジッターパルス列の平均ＰＲＩのところにする
どいピークがたつ。しきい値演算部５５０ではジッター
パルス列のＰＲＩ検出のための適当なしきい値を定め
る。しきい値演算部５５０では同時に、検出可能なＰＲ
Ｉの範囲も求める。次に、ＰＲＩ抽出部５６０で、その
しきい値を超える値を検出するようにすれば、ジッター
パルス列の検出可能な範囲にあるＰＲＩがすべて検出さ
れる。このようにして検出されたＰＲＩの値を利用し
て、ジッターパルス列分離演算部５７０で、もとのパル
ス列から該当する部分パルス列を分離する。この分離の
ためには、例えば特公平１−４７９３６号公報に詳しく
示されているパルス繰り返し間隔フィルターを用いる。
残存パルス列判定部５８０では、調べたいＰＲＩの範囲
と検出可能なＰＲＩの範囲を比較し、もし検出可能なＰ
ＲＩの範囲の方が狭ければ、すなわち、ジッターパルス
列が残存する可能性があれば、残ったパルス列に再び変
形ＰＲＩ変換演算部５３０でＰＲＩスペクトルを作り、
処理を繰り返す。そうでなければ、出力端５９０へ進
む。
【００９５】ジッターパルス列に対しては、先に述べた
ＰＲＩ変換ではうまくＰＲＩのところにピークをたたせ
ることができないが、これから述べる変形ＰＲＩ変換で
はジッターパルス列であってもその平均ＰＲＩのところ
にピークがたつようにできる。
ＰＲＩ変換ではうまくＰＲＩのところにピークをたたせ
ることができないが、これから述べる変形ＰＲＩ変換で
はジッターパルス列であってもその平均ＰＲＩのところ
にピークがたつようにできる。
【００９６】変形ＰＲＩ変換は、ＰＲＩ変換のＰＲＩビ
ンをオーバーラップさせ、時間原点の移動を行うように
したものである。
ンをオーバーラップさせ、時間原点の移動を行うように
したものである。
【００９７】図１０は変形ＰＲＩ変換で用いるＰＲＩビ
ンの説明図である。調べたいＰＲＩの下限と上限をそれ
ぞれ、τmin、τmax と置く。 ＰＲＩ軸での区間 [τ
min，τmax] 内のＰＲＩビンの個数をKとし、ＰＲＩの
揺らぎの大きさの上限をεとする（εは片側で表現す
る。例えばＰＲＩの揺らぎの上限がピーク・ツー・ピー
クで３０％であったとするとε= 0．15 となる）。各Ｐ
ＲＩビンの中心ＰＲＩの値τk とビンの幅 bk は次式
（１４）のように決める。
ンの説明図である。調べたいＰＲＩの下限と上限をそれ
ぞれ、τmin、τmax と置く。 ＰＲＩ軸での区間 [τ
min，τmax] 内のＰＲＩビンの個数をKとし、ＰＲＩの
揺らぎの大きさの上限をεとする（εは片側で表現す
る。例えばＰＲＩの揺らぎの上限がピーク・ツー・ピー
クで３０％であったとするとε= 0．15 となる）。各Ｐ
ＲＩビンの中心ＰＲＩの値τk とビンの幅 bk は次式
（１４）のように決める。
【００９８】
【数１４】
【００９９】
【数１５】
【０１００】図１１には変形ＰＲＩ変換の処理のフロー
チャートを示す。入力データのパルスのＴＯＡを t1、
t2、 …、 tNとする。変形ＰＲＩ変換の処理の手順は以
下のように行う。
チャートを示す。入力データのパルスのＴＯＡを t1、
t2、 …、 tNとする。変形ＰＲＩ変換の処理の手順は以
下のように行う。
【０１０１】（１）初期設定を行う（ステップＳ５３
１）。すなわち、ＰＲＩ変換 Dk（ｋ= 1、 2、 … 、
K）をゼロクリアする。また、パルスペアの右パルスの
パルス番号を m = 1 に設定する。
１）。すなわち、ＰＲＩ変換 Dk（ｋ= 1、 2、 … 、
K）をゼロクリアする。また、パルスペアの右パルスの
パルス番号を m = 1 に設定する。
【０１０２】（２）パルスペアの右パルスのパルス番号
を更新する（ステップＳ５３２）。すなわち、m = m +
1とする。
を更新する（ステップＳ５３２）。すなわち、m = m +
1とする。
【０１０３】（３）終了判定を行う（ステップＳ５３
３）。すなわち、m≦Nならばステップ５３４ へ。そう
でなければ終了。
３）。すなわち、m≦Nならばステップ５３４ へ。そう
でなければ終了。
【０１０４】（４）パルスペアの左パルスのパルス番号
を n = m に設定する。（ステップＳ５３４）。
を n = m に設定する。（ステップＳ５３４）。
【０１０５】（５）パルスペアの左パルスのパルス番号
を更新する（ステップＳ５３５）。すなわち、n = n−
1。
を更新する（ステップＳ５３５）。すなわち、n = n−
1。
【０１０６】（６）パルスペアの左パルスの終了判定を
行う（ステップ ５３６）。すなわち、n ≧ 1 ならば、
ステップＳ５３７へ。そうでなければステップＳ５３２
へ。
行う（ステップ ５３６）。すなわち、n ≧ 1 ならば、
ステップＳ５３７へ。そうでなければステップＳ５３２
へ。
【０１０７】（７）パルスペアの時間間隔を求める（ス
テップ ５３７）。すなわち、τ=tm−tnを行う。
テップ ５３７）。すなわち、τ=tm−tnを行う。
【０１０８】（８）時間間隔の大きさ判定１を行う（ス
テップＳ５３８）。すなわち、τ≧τmin ならば、ステ
ップＳ５３９へ。そうでなければステップＳ５３５へ。
テップＳ５３８）。すなわち、τ≧τmin ならば、ステ
ップＳ５３９へ。そうでなければステップＳ５３５へ。
【０１０９】（９）時間間隔の大きさ判定２を行う（ス
テップＳ５３９）。すなわち、τ≦τmax ならば、ステ
ップＳ５４０へ。そうでなければステップＳ５３２へ。
テップＳ５３９）。すなわち、τ≦τmax ならば、ステ
ップＳ５４０へ。そうでなければステップＳ５３２へ。
【０１１０】（１０）τ軸上でのＰＲＩビンの選定を行
う（ステップＳ５４０）。すなわち、ＰＲＩビンの下限
k1 と上限 k2 を次式で計算する。
う（ステップＳ５４０）。すなわち、ＰＲＩビンの下限
k1 と上限 k2 を次式で計算する。
【０１１１】
【数１６】
【０１１２】ここで、[・] は Gauss の整数、すなわ
ち、(・) を越えない最大の整数を表わす。
ち、(・) を越えない最大の整数を表わす。
【０１１３】（１１）ＰＲＩ変換の更新を行う（ステッ
プＳ５４１）。すなわち、k = k1，…， k2 について以
下の処理を行う。
プＳ５４１）。すなわち、k = k1，…， k2 について以
下の処理を行う。
【０１１４】１）時間原点の初期設定を行う。すなわ
ち、ＰＲＩビンｋが最初に使われる場合に限ってOk = t
m とする。
ち、ＰＲＩビンｋが最初に使われる場合に限ってOk = t
m とする。
【０１１５】２）位相の予備計算と分解を行う。位相η
0を次式で求め、それを整数部分νと小数部分ζに分解
する。
0を次式で求め、それを整数部分νと小数部分ζに分解
する。
【０１１６】
【数１７】
【０１１７】３）時間原点の移動を行う。すなわち、次
の条件のいずれかが満たされるとき、時間原点を移動し
て、tm をあらたな時間原点とする。
の条件のいずれかが満たされるとき、時間原点を移動し
て、tm をあらたな時間原点とする。
【０１１８】
【数１８】
【０１１９】４）位相の計算を行う。次式（１９）で新
たな時間原点のもとでの位相を計算する。
たな時間原点のもとでの位相を計算する。
【０１２０】
【数１９】
【０１２１】５）ＰＲＩ変換の価の更新を行う。次式
（２０）でＰＲＩ変換の価を更新する。このあと、ステ
ップ５３２へ。
（２０）でＰＲＩ変換の価を更新する。このあと、ステ
ップ５３２へ。
【０１２２】
【数２０】
【０１２３】図１２乃至図１５にはこの実施の形態によ
る発明の効果を示す。
る発明の効果を示す。
【０１２４】図１２は、本発明のパルス列分離装置への
入力パルスデータの説明図である。図１２において、６
１（出力１）と６２（出力２）はスタガーパルス列、６
３（出力３）と６５（出力５）はジッターパルス列、６
４（出力４）と６６（出力６）は固定ＰＲＩパルス列、
６７は６１（出力１）から６６（出力６）までの６つの
パルス列が重なりあった複合パルス列である。本発明の
パルス列分離装置へは、複合パルス列６７が入力する。
入力パルスデータの説明図である。図１２において、６
１（出力１）と６２（出力２）はスタガーパルス列、６
３（出力３）と６５（出力５）はジッターパルス列、６
４（出力４）と６６（出力６）は固定ＰＲＩパルス列、
６７は６１（出力１）から６６（出力６）までの６つの
パルス列が重なりあった複合パルス列である。本発明の
パルス列分離装置へは、複合パルス列６７が入力する。
【０１２５】図１３は、図２の固定ＰＲＩパルス列検出
・分離手段３００で固定ＰＲＩパルス列が検出され分離
される様子を示した図である。図１３において、７１は
固定ＰＲＩパルス列検出・分離手段３００への入力パル
ス列、７２はＰＲＩ変換と検出結果を示すグラフ、７７
と７８は分離されたパルス列を表わす。７１の入力パル
ス列が図３のＰＲＩ変換演算部３３０へ入力されると、
図５にフローチャートを示したＰＲＩ変換により、７３
のようなＰＲＩスペクトルが得られる。しきい値演算部
３５０で求めたしきい値７４とＰＲＩスペクトル７３と
の比較をＰＲＩ抽出部３６０で行うことにより、７５
（ピーク１）及び７６（ピーク２）の２つのＰＲＩが検
出される。検出されたＰＲＩをもとに、固定ＰＲＩパル
ス列分離演算部３７０により、もとのパルス列７１から
抜き出したパルス列が７７と７８である。７５のピーク
１は７７のパルス列に対応し、７６のピーク２は７８の
パルス列に対応する。
・分離手段３００で固定ＰＲＩパルス列が検出され分離
される様子を示した図である。図１３において、７１は
固定ＰＲＩパルス列検出・分離手段３００への入力パル
ス列、７２はＰＲＩ変換と検出結果を示すグラフ、７７
と７８は分離されたパルス列を表わす。７１の入力パル
ス列が図３のＰＲＩ変換演算部３３０へ入力されると、
図５にフローチャートを示したＰＲＩ変換により、７３
のようなＰＲＩスペクトルが得られる。しきい値演算部
３５０で求めたしきい値７４とＰＲＩスペクトル７３と
の比較をＰＲＩ抽出部３６０で行うことにより、７５
（ピーク１）及び７６（ピーク２）の２つのＰＲＩが検
出される。検出されたＰＲＩをもとに、固定ＰＲＩパル
ス列分離演算部３７０により、もとのパルス列７１から
抜き出したパルス列が７７と７８である。７５のピーク
１は７７のパルス列に対応し、７６のピーク２は７８の
パルス列に対応する。
【０１２６】図１４は、図２のスタガーパルス列検出・
分離手段４００でスタガーパルス列が検出され分離され
る様子を示した図である。図１４において、８１はスタ
ガーパルス列検出・分離手段４００への入力パルス列、
８２は３連パルスのヒストグラムと検出結果を示すグラ
フ、８７と８８は分離されたパルス列を表わす。８１の
入力パルス列が、図６の３連パルスヒストグラム演算部
４３０へ入力されると、図７にフローチャートを示した
３連パルスのヒストグラムの処理により、８３のような
ヒストグラムが得られる。このヒストグラム８３としき
い値８４との比較により、８５（ピーク１）及び８６
（ピーク２）の２つのＰＲＩが検出される。検出された
ＰＲＩをもとに、もとのパルス列８１から抜き出したパ
ルス列が８７と８８である。８５のピーク１は８７のパ
ルス列に対応し、８６のピーク２は８８のパルス列に対
応する。
分離手段４００でスタガーパルス列が検出され分離され
る様子を示した図である。図１４において、８１はスタ
ガーパルス列検出・分離手段４００への入力パルス列、
８２は３連パルスのヒストグラムと検出結果を示すグラ
フ、８７と８８は分離されたパルス列を表わす。８１の
入力パルス列が、図６の３連パルスヒストグラム演算部
４３０へ入力されると、図７にフローチャートを示した
３連パルスのヒストグラムの処理により、８３のような
ヒストグラムが得られる。このヒストグラム８３としき
い値８４との比較により、８５（ピーク１）及び８６
（ピーク２）の２つのＰＲＩが検出される。検出された
ＰＲＩをもとに、もとのパルス列８１から抜き出したパ
ルス列が８７と８８である。８５のピーク１は８７のパ
ルス列に対応し、８６のピーク２は８８のパルス列に対
応する。
【０１２７】図１５は、図２のジッターパルス列検出・
分離手段５００でジッターパルス列が検出され分離され
る様子を示した図である。図１５において、９１はジッ
ターパルス列検出・分離手段５００への入力パルス列、
９２は変形ＰＲＩ変換と検出結果を示すグラフ、９７と
９８は分離されたパルス列を表わす。９１の入力パルス
列が図９の変形ＰＲＩ変換演算部５３０へ入力される
と、図１１にフローチャートを示した変形ＰＲＩ変換の
処理により、９３のようなＰＲＩスペクトルが得られ
る。このＰＲＩスペクトル９３としきい値９４との比較
により、９５（ピーク１）及び９６（ピーク２）の２つ
のＰＲＩが検出される。検出されたＰＲＩをもとにもと
のパルス列９１から抜き出したパルス列が９７と９８で
ある。９５のピーク１は９７のパルス列に対応し、９６
のピーク２は９８のパルス列に対応する。
分離手段５００でジッターパルス列が検出され分離され
る様子を示した図である。図１５において、９１はジッ
ターパルス列検出・分離手段５００への入力パルス列、
９２は変形ＰＲＩ変換と検出結果を示すグラフ、９７と
９８は分離されたパルス列を表わす。９１の入力パルス
列が図９の変形ＰＲＩ変換演算部５３０へ入力される
と、図１１にフローチャートを示した変形ＰＲＩ変換の
処理により、９３のようなＰＲＩスペクトルが得られ
る。このＰＲＩスペクトル９３としきい値９４との比較
により、９５（ピーク１）及び９６（ピーク２）の２つ
のＰＲＩが検出される。検出されたＰＲＩをもとにもと
のパルス列９１から抜き出したパルス列が９７と９８で
ある。９５のピーク１は９７のパルス列に対応し、９６
のピーク２は９８のパルス列に対応する。
【０１２８】以上のように、本発明のパルス列分離装置
によれば、ＰＲＩ変換を用いる固定ＰＲＩパルス列検出
・分離手段３００と、３連パルスのヒストグラムを用い
るスタガーパルス列検出・分離手段４００と、変形ＰＲ
Ｉ変換を用いるジッターパルス列検出・分離手段を備え
て、ＰＲＩの種類ごとに該当するパルス列を段階的に検
出・分離するようにしたので、信号源の数が多い場合
や、固定ＰＲＩ以外の複雑なＰＲＩのパターンを有する
信号源まで含めて、信号源のＰＲＩを正確に推定し、精
度良く、各パルス列を検出・分離できるという効果があ
る。
によれば、ＰＲＩ変換を用いる固定ＰＲＩパルス列検出
・分離手段３００と、３連パルスのヒストグラムを用い
るスタガーパルス列検出・分離手段４００と、変形ＰＲ
Ｉ変換を用いるジッターパルス列検出・分離手段を備え
て、ＰＲＩの種類ごとに該当するパルス列を段階的に検
出・分離するようにしたので、信号源の数が多い場合
や、固定ＰＲＩ以外の複雑なＰＲＩのパターンを有する
信号源まで含めて、信号源のＰＲＩを正確に推定し、精
度良く、各パルス列を検出・分離できるという効果があ
る。
【０１２９】
【発明の効果】この発明は、複数の信号源からのパルス
列が重なり合った複合パルス列の各パルスの到着時間
（ＴＯＡ）データに基づき、固定パルス繰り返し間隔
（ＰＲＩ）を有する固定ＰＲＩパルス列を検出する固定
ＰＲＩパルス列検出手段と、検出された固定ＰＲＩパル
ス列をもとの複合パルス列から分離するとともに、分離
後の残りのパルス列を出力する固定ＰＲＩパルス列分離
手段と、固定ＰＲＩパルス列分離手段から出力される残
りのパルス列から、複合パルス列の各パルスの到着時間
（ＴＯＡ）データに基づき、スタガーパルス列を検出す
るスタガーパルス列検出手段と、検出されたスタガーパ
ルス列を残りのパルス列から分離するとともに、分離後
の固定ＰＲＩパルス列及びスタガーパルス列を含まない
第二の残りのパルス列を出力するスタガーパルス列分離
手段と、スタガーパルス列分離手段から出力される第二
の残りのパルス列から、複合パルス列の各パルスの到着
時間（ＴＯＡ）データに基づき、ジッターパルス列を検
出するジッターパルス列検出手段と、検出されたジッタ
ーパルス列を第二の残りのパルス列から分離するジッタ
ーパルス列分離手段と、を備えたパルス列分離装置であ
るので、ＰＲＩの種類ごとに該当するパルス列を段階的
に検出・分離するようにしたので、信号源の数が多い場
合や、固定ＰＲＩ以外の複雑なＰＲＩのパターンを有す
る信号源まで含めて、信号源のＰＲＩを正確に推定し、
精度良く、各パルス列を検出・分離できるという効果が
ある。
列が重なり合った複合パルス列の各パルスの到着時間
（ＴＯＡ）データに基づき、固定パルス繰り返し間隔
（ＰＲＩ）を有する固定ＰＲＩパルス列を検出する固定
ＰＲＩパルス列検出手段と、検出された固定ＰＲＩパル
ス列をもとの複合パルス列から分離するとともに、分離
後の残りのパルス列を出力する固定ＰＲＩパルス列分離
手段と、固定ＰＲＩパルス列分離手段から出力される残
りのパルス列から、複合パルス列の各パルスの到着時間
（ＴＯＡ）データに基づき、スタガーパルス列を検出す
るスタガーパルス列検出手段と、検出されたスタガーパ
ルス列を残りのパルス列から分離するとともに、分離後
の固定ＰＲＩパルス列及びスタガーパルス列を含まない
第二の残りのパルス列を出力するスタガーパルス列分離
手段と、スタガーパルス列分離手段から出力される第二
の残りのパルス列から、複合パルス列の各パルスの到着
時間（ＴＯＡ）データに基づき、ジッターパルス列を検
出するジッターパルス列検出手段と、検出されたジッタ
ーパルス列を第二の残りのパルス列から分離するジッタ
ーパルス列分離手段と、を備えたパルス列分離装置であ
るので、ＰＲＩの種類ごとに該当するパルス列を段階的
に検出・分離するようにしたので、信号源の数が多い場
合や、固定ＰＲＩ以外の複雑なＰＲＩのパターンを有す
る信号源まで含めて、信号源のＰＲＩを正確に推定し、
精度良く、各パルス列を検出・分離できるという効果が
ある。
【０１３０】また、固定ＰＲＩパルス列検出手段が、Ｐ
ＲＩ変換を用いて検出を行うようにしたので、複合パル
ス列に含まれている部分パルス列のＰＲＩだけをピーク
として抽出できるので、ＰＲＩを明確に識別することが
できるという効果がある。
ＲＩ変換を用いて検出を行うようにしたので、複合パル
ス列に含まれている部分パルス列のＰＲＩだけをピーク
として抽出できるので、ＰＲＩを明確に識別することが
できるという効果がある。
【０１３１】また、スタガーパルス列検出手段が、等間
隔に並ぶ３つのパルスからなる３連パルスの集合に対し
てそのパルス繰り返し間隔のヒストグラムを求める３連
パルスヒストグラム演算部と、得られたヒストグラムか
ら複数のピークを検出し、各ピークがパルス列の基本周
期であるか否かを判定する判定部と、を備えているの
で、得られるヒストグラムでは、スタガーパルス列の基
本周期とその整数倍のところにするどいピークがたつの
で、スタガーパルスの基本周期を明確に識別することが
できるという効果がある。
隔に並ぶ３つのパルスからなる３連パルスの集合に対し
てそのパルス繰り返し間隔のヒストグラムを求める３連
パルスヒストグラム演算部と、得られたヒストグラムか
ら複数のピークを検出し、各ピークがパルス列の基本周
期であるか否かを判定する判定部と、を備えているの
で、得られるヒストグラムでは、スタガーパルス列の基
本周期とその整数倍のところにするどいピークがたつの
で、スタガーパルスの基本周期を明確に識別することが
できるという効果がある。
【０１３２】また、ジッターパルス列検出手段が、オー
バーラップしたＰＲＩビンを持ち、時間原点の移動を行
う、変形ＰＲＩ変換を用いるようにしたので、ジッター
パルス列の平均ＰＲＩのところにするどいピークがたつ
ので、ジッターパルス列の平均ＰＲＩを明確に識別する
ことができるという効果がある。
バーラップしたＰＲＩビンを持ち、時間原点の移動を行
う、変形ＰＲＩ変換を用いるようにしたので、ジッター
パルス列の平均ＰＲＩのところにするどいピークがたつ
ので、ジッターパルス列の平均ＰＲＩを明確に識別する
ことができるという効果がある。
【０１３３】また、この発明は、複数の信号源からのパ
ルス列が重なり合った複合パルス列の各パルスの到着時
間（ＴＯＡ）データに基づき、固定パルス繰り返し間隔
（ＰＲＩ）を有する固定ＰＲＩパルス列を検出する固定
ＰＲＩパルス列検出工程と、検出された固定ＰＲＩパル
ス列をもとの複合パルス列から分離するとともに、分離
後の残りのパルス列を出力する固定ＰＲＩパルス列分離
工程と、固定ＰＲＩパルス列分離手段から出力される残
りのパルス列から、複合パルス列の各パルスの到着時間
（ＴＯＡ）データに基づき、スタガーパルス列を検出す
るスタガーパルス列検出工程と、検出されたスタガーパ
ルス列を残りのパルス列から分離するとともに、分離後
の固定ＰＲＩパルス列及びスタガーパルス列を含まない
第二の残りのパルス列を出力するスタガーパルス列分離
工程と、スタガーパルス列分離手段から出力される第二
の残りのパルス列から、複合パルス列の各パルスの到着
時間（ＴＯＡ）データに基づき、ジッターパルス列を検
出するジッターパルス列検出工程と、検出されたジッタ
ーパルス列を第二の残りのパルス列から分離するジッタ
ーパルス列分離工程と、を備えたパルス列分離方法であ
るので、ＰＲＩの種類ごとに該当するパルス列を段階的
に検出・分離するようにしたので、信号源の数が多い場
合や、固定ＰＲＩ以外の複雑なＰＲＩのパターンを有す
る信号源まで含めて、信号源のＰＲＩを正確に推定し、
精度良く、各パルス列を検出・分離できるという効果が
ある。
ルス列が重なり合った複合パルス列の各パルスの到着時
間（ＴＯＡ）データに基づき、固定パルス繰り返し間隔
（ＰＲＩ）を有する固定ＰＲＩパルス列を検出する固定
ＰＲＩパルス列検出工程と、検出された固定ＰＲＩパル
ス列をもとの複合パルス列から分離するとともに、分離
後の残りのパルス列を出力する固定ＰＲＩパルス列分離
工程と、固定ＰＲＩパルス列分離手段から出力される残
りのパルス列から、複合パルス列の各パルスの到着時間
（ＴＯＡ）データに基づき、スタガーパルス列を検出す
るスタガーパルス列検出工程と、検出されたスタガーパ
ルス列を残りのパルス列から分離するとともに、分離後
の固定ＰＲＩパルス列及びスタガーパルス列を含まない
第二の残りのパルス列を出力するスタガーパルス列分離
工程と、スタガーパルス列分離手段から出力される第二
の残りのパルス列から、複合パルス列の各パルスの到着
時間（ＴＯＡ）データに基づき、ジッターパルス列を検
出するジッターパルス列検出工程と、検出されたジッタ
ーパルス列を第二の残りのパルス列から分離するジッタ
ーパルス列分離工程と、を備えたパルス列分離方法であ
るので、ＰＲＩの種類ごとに該当するパルス列を段階的
に検出・分離するようにしたので、信号源の数が多い場
合や、固定ＰＲＩ以外の複雑なＰＲＩのパターンを有す
る信号源まで含めて、信号源のＰＲＩを正確に推定し、
精度良く、各パルス列を検出・分離できるという効果が
ある。
【０１３４】また、固定ＰＲＩパルス列検出工程におい
て、ＰＲＩ変換を用いて検出を行うようにしたので、複
合パルス列に含まれている部分パルス列のＰＲＩだけを
ピークとして抽出できるので、ＰＲＩを明確に識別する
ことができるという効果がある。
て、ＰＲＩ変換を用いて検出を行うようにしたので、複
合パルス列に含まれている部分パルス列のＰＲＩだけを
ピークとして抽出できるので、ＰＲＩを明確に識別する
ことができるという効果がある。
【０１３５】また、スタガーパルス列検出工程が、等間
隔に並ぶ３つのパルスからなる３連パルスの集合に対し
てそのパルス繰り返し間隔のヒストグラムを求める３連
パルスヒストグラム演算ステップと、得られたヒストグ
ラムから複数のピークを検出し、各ピークがパルス列の
基本周期であるか否かを判定する判定ステップと、を備
えているので、得られるヒストグラムでは、スタガーパ
ルス列の基本周期とその整数倍のところにするどいピー
クがたつので、スタガーパルスの基本周期を明確に識別
することができるという効果がある。
隔に並ぶ３つのパルスからなる３連パルスの集合に対し
てそのパルス繰り返し間隔のヒストグラムを求める３連
パルスヒストグラム演算ステップと、得られたヒストグ
ラムから複数のピークを検出し、各ピークがパルス列の
基本周期であるか否かを判定する判定ステップと、を備
えているので、得られるヒストグラムでは、スタガーパ
ルス列の基本周期とその整数倍のところにするどいピー
クがたつので、スタガーパルスの基本周期を明確に識別
することができるという効果がある。
【０１３６】また、ジッターパルス列検出工程におい
て、オーバーラップしたＰＲＩビンを持ち、時間原点の
移動を行う、変形ＰＲＩ変換を用いるようにしたので、
ジッターパルス列の平均ＰＲＩのところにするどいピー
クがたつので、ジッターパルス列の平均ＰＲＩを明確に
識別することができるという効果がある。
て、オーバーラップしたＰＲＩビンを持ち、時間原点の
移動を行う、変形ＰＲＩ変換を用いるようにしたので、
ジッターパルス列の平均ＰＲＩのところにするどいピー
クがたつので、ジッターパルス列の平均ＰＲＩを明確に
識別することができるという効果がある。
【図１】 ＰＲＩの種類を説明した説明図である。。
【図２】 本発明の実施の形態１にかかわるパルス列分
離装置の構成を示すブロック図である。
離装置の構成を示すブロック図である。
【図３】 この発明の実施の形態１にかかわるパルス列
分離装置における固定ＰＲＩパルス列検出・分離手段の
内部構成を示すブロック図である。
分離装置における固定ＰＲＩパルス列検出・分離手段の
内部構成を示すブロック図である。
【図４】 この発明の実施の形態１にかかわるパルス列
分離装置における固定ＰＲＩパルス列検出・分離手段の
中で用いるＰＲＩ変換に関するＰＲＩビンの説明図であ
る。
分離装置における固定ＰＲＩパルス列検出・分離手段の
中で用いるＰＲＩ変換に関するＰＲＩビンの説明図であ
る。
【図５】 この発明の実施の形態１にかかわるパルス列
分離装置における固定ＰＲＩパルス列検出・分離手段の
中で用いるＰＲＩ変換のフローチャートである。
分離装置における固定ＰＲＩパルス列検出・分離手段の
中で用いるＰＲＩ変換のフローチャートである。
【図６】 この発明の実施の形態１にかかわるパルス列
分離装置におけるスタガーパルス列検出・分離手段の内
部構成を示すブロック図である。
分離装置におけるスタガーパルス列検出・分離手段の内
部構成を示すブロック図である。
【図７】 この発明の実施の形態１にかかわるパルス列
分離装置におけるスタガーパルス列検出・分離手段の中
で用いる３連パルスのヒストグラムのフローチャートで
ある。
分離装置におけるスタガーパルス列検出・分離手段の中
で用いる３連パルスのヒストグラムのフローチャートで
ある。
【図８】 この発明の実施の形態１にかかわるパルス分
離装置におけるスタガーパルス列検出・分離手段の中で
用いる３連パルスのヒストグラムからの基本周期検出の
フローチャートである。
離装置におけるスタガーパルス列検出・分離手段の中で
用いる３連パルスのヒストグラムからの基本周期検出の
フローチャートである。
【図９】 この発明の実施の形態１にかかわるパルス列
分離装置におけるジッターパルス列検出・分離手段の内
部構成を示すブロック図である。
分離装置におけるジッターパルス列検出・分離手段の内
部構成を示すブロック図である。
【図１０】 この発明の実施の形態１にかかわるパルス
列分離装置におけるジッターパルス列検出・分離手段の
中で用いる変形ＰＲＩ変換に関するＰＲＩビンの説明図
である。
列分離装置におけるジッターパルス列検出・分離手段の
中で用いる変形ＰＲＩ変換に関するＰＲＩビンの説明図
である。
【図１１】 この発明の実施の形態１にかかわるパルス
列分離装置におけるジッターパルス列検出・分離手段の
中で用いる変形ＰＲＩ変換のフローチャートである。
列分離装置におけるジッターパルス列検出・分離手段の
中で用いる変形ＰＲＩ変換のフローチャートである。
【図１２】 この発明の実施の形態１にかかわるパルス
列分離装置への入力パルス列の説明図である。
列分離装置への入力パルス列の説明図である。
【図１３】 この発明の実施の形態１にかかわるパルス
列分離装置における固定ＰＲＩパルス列検出・分離手段
での処理結果を示す図である。
列分離装置における固定ＰＲＩパルス列検出・分離手段
での処理結果を示す図である。
【図１４】 この発明の実施の形態１にかかわるパルス
列分離装置におけるスタガーパルス列検出・分離手段で
の処理結果を示す図である。
列分離装置におけるスタガーパルス列検出・分離手段で
の処理結果を示す図である。
【図１５】 この発明の実施の形態１にかかわるパルス
列分離装置におけるジッターパルス列検出・分離手段で
の処理結果を示す図である。
列分離装置におけるジッターパルス列検出・分離手段で
の処理結果を示す図である。
【図１６】 従来のパルス列分離装置の動作を示すフロ
ーチャートである。
ーチャートである。
１１ 固定ＰＲＩパルス列、１２ スタガーパルス列、
１３ ジッターパルス列、６１，６２ スタガーパルス
列、６３，６５ ジッターパルス列、６４，６６ 固定
ＰＲＩパルス列、６７ 複合パルス列、７１ 固定ＰＲ
Ｉパルス列検出・分離手段への入力パルス列、７２ Ｐ
ＲＩ変換と検出結果、７３ ＰＲＩスペクトル、７４
ＰＲＩスペクトルのしきい値、７５，７６ 検出された
固定ＰＲＩパルス列のＰＲＩ、７７，７８ 分離された
固定ＰＲＩパルス列、８１ スタガーパルス列検出・分
離手段への入力パルス列、８２ ３連パルスのヒストグ
ラムと検出結果、８３ ３連パルスのヒストグラム、８
４ ヒストグラムのしきい値、８５，８６ 検出された
スタガーパルス列のＰＲＩ、８７，８８ 分離されたス
タガーパルス列、９１ ジッターパルス列検出・分離手
段への入力パルス列、９２ 変形ＰＲＩ変換と検出結
果、９３ ＰＲＩスペクトル、９４ ＰＲＩスペクトル
のしきい値、９５，９６ 検出されたジッターパルス列
のＰＲＩ、９７，９８ 分離されたジッターパルス列、
１００ 広帯域受信機、２００ パルスＴＯＡデータの
入力端、３００ 固定ＰＲＩパルス列検出・分離手段、
４００スタガーパルス列検出・分離手段、５００ ジッ
ターパルス列検出・分離手段、６００ 出力端、３１０
固定ＰＲＩパルス列検出・分離手段への入力端、３３
０ ＰＲＩ変換演算部、３５０ しきい値演算部、３６
０ 固定ＰＲＩパルス列のＰＲＩ抽出部、３７０ 固定
ＰＲＩパルス列分離演算部、３８０ 残存パルス列判定
部、３９０ 固定ＰＲＩパルス列検出・分離手段からの
出力端、４１０スタガーパルス列検出・分離手段への入
力端、４３０ ３連パルスヒストグラム演算部、４５０
しきい値演算部、４６０ スタガーパルス列のＰＲＩ
抽出部、４７０ スタガーパルス列分離演算部、４８０
残存パルス列判定部、４９０スタガーパルス列検出・
分離手段からの出力端、５１０ ジッターパルス列検出
・分離手段への入力端、５３０ 変形ＰＲＩ変換演算
部、５５０ しきい値演算部、５６０ ジッターパルス
列のＰＲＩ抽出部、５７０ ジッターパルス列分離演算
部、５８０ 残存パルス列判定部、５９０ ジッターパ
ルス列検出・分離手段からの出力端。
１３ ジッターパルス列、６１，６２ スタガーパルス
列、６３，６５ ジッターパルス列、６４，６６ 固定
ＰＲＩパルス列、６７ 複合パルス列、７１ 固定ＰＲ
Ｉパルス列検出・分離手段への入力パルス列、７２ Ｐ
ＲＩ変換と検出結果、７３ ＰＲＩスペクトル、７４
ＰＲＩスペクトルのしきい値、７５，７６ 検出された
固定ＰＲＩパルス列のＰＲＩ、７７，７８ 分離された
固定ＰＲＩパルス列、８１ スタガーパルス列検出・分
離手段への入力パルス列、８２ ３連パルスのヒストグ
ラムと検出結果、８３ ３連パルスのヒストグラム、８
４ ヒストグラムのしきい値、８５，８６ 検出された
スタガーパルス列のＰＲＩ、８７，８８ 分離されたス
タガーパルス列、９１ ジッターパルス列検出・分離手
段への入力パルス列、９２ 変形ＰＲＩ変換と検出結
果、９３ ＰＲＩスペクトル、９４ ＰＲＩスペクトル
のしきい値、９５，９６ 検出されたジッターパルス列
のＰＲＩ、９７，９８ 分離されたジッターパルス列、
１００ 広帯域受信機、２００ パルスＴＯＡデータの
入力端、３００ 固定ＰＲＩパルス列検出・分離手段、
４００スタガーパルス列検出・分離手段、５００ ジッ
ターパルス列検出・分離手段、６００ 出力端、３１０
固定ＰＲＩパルス列検出・分離手段への入力端、３３
０ ＰＲＩ変換演算部、３５０ しきい値演算部、３６
０ 固定ＰＲＩパルス列のＰＲＩ抽出部、３７０ 固定
ＰＲＩパルス列分離演算部、３８０ 残存パルス列判定
部、３９０ 固定ＰＲＩパルス列検出・分離手段からの
出力端、４１０スタガーパルス列検出・分離手段への入
力端、４３０ ３連パルスヒストグラム演算部、４５０
しきい値演算部、４６０ スタガーパルス列のＰＲＩ
抽出部、４７０ スタガーパルス列分離演算部、４８０
残存パルス列判定部、４９０スタガーパルス列検出・
分離手段からの出力端、５１０ ジッターパルス列検出
・分離手段への入力端、５３０ 変形ＰＲＩ変換演算
部、５５０ しきい値演算部、５６０ ジッターパルス
列のＰＲＩ抽出部、５７０ ジッターパルス列分離演算
部、５８０ 残存パルス列判定部、５９０ ジッターパ
ルス列検出・分離手段からの出力端。
Claims (8)
- 【請求項１】 複数の信号源からのパルス列が重なり合
った複合パルス列の各パルスの到着時間（ＴＯＡ）デー
タに基づき、固定パルス繰り返し間隔（ＰＲＩ）を有す
る固定ＰＲＩパルス列を検出する固定ＰＲＩパルス列検
出手段と、 検出された上記固定ＰＲＩパルス列をもとの上記複合パ
ルス列から分離するとともに、分離後の残りのパルス列
を出力する固定ＰＲＩパルス列分離手段と、 上記固定ＰＲＩパルス列分離手段から出力される上記残
りのパルス列から、上記複合パルス列の各パルスの到着
時間（ＴＯＡ）データに基づき、スタガーパルス列を検
出するスタガーパルス列検出手段と、 検出されたスタガーパルス列を上記残りのパルス列から
分離するとともに、分離後の固定ＰＲＩパルス列及びス
タガーパルス列を含まない第二の残りのパルス列を出力
するスタガーパルス列分離手段と、 上記スタガーパルス列分離手段から出力される上記第二
の残りのパルス列から、上記複合パルス列の各パルスの
到着時間（ＴＯＡ）データに基づき、ジッターパルス列
を検出するジッターパルス列検出手段と、 検出された上記ジッターパルス列を上記第二の残りのパ
ルス列から分離するジッターパルス列分離手段と、 を備えたことを特徴とするパルス列分離装置。 - 【請求項２】 上記固定ＰＲＩパルス列検出手段が、Ｐ
ＲＩ変換を用いて検出を行うことを特徴とする請求項１
記載のパルス列分離装置。 - 【請求項３】 上記スタガーパルス列検出手段が、 等間隔に並ぶ３つのパルスからなる３連パルスの集合に
対してそのパルス繰り返し間隔のヒストグラムを求める
３連パルスヒストグラム演算部と、 得られた上記ヒストグラムから複数のピークを検出し、
各ピークがパルス列の基本周期であるか否かを判定する
判定部と、 を備えたことを特徴とする請求項１または２に記載のパ
ルス列分離装置。 - 【請求項４】 上記ジッターパルス列検出手段が、オー
バーラップしたＰＲＩビンを持ち、時間原点の移動を行
う、変形ＰＲＩ変換を用いることを特徴とする請求項１
ないし３のいずれかに記載のパルス列分離装置。 - 【請求項５】 複数の信号源からのパルス列が重なり合
った複合パルス列の各パルスの到着時間（ＴＯＡ）デー
タに基づき、固定パルス繰り返し間隔（ＰＲＩ）を有す
る固定ＰＲＩパルス列を検出する固定ＰＲＩパルス列検
出工程と、 検出された上記固定ＰＲＩパルス列をもとの上記複合パ
ルス列から分離するとともに、分離後の残りのパルス列
を出力する固定ＰＲＩパルス列分離工程と、 上記固定ＰＲＩパルス列分離手段から出力される上記残
りのパルス列から、上記複合パルス列の各パルスの到着
時間（ＴＯＡ）データに基づき、スタガーパルス列を検
出するスタガーパルス列検出工程と、 検出された上記スタガーパルス列を上記残りのパルス列
から分離するとともに、分離後の固定ＰＲＩパルス列及
びスタガーパルス列を含まない第二の残りのパルス列を
出力するスタガーパルス列分離工程と、 上記スタガーパルス列分離手段から出力される上記第二
の残りのパルス列から、上記複合パルス列の各パルスの
到着時間（ＴＯＡ）データに基づき、ジッターパルス列
を検出するジッターパルス列検出工程と、 検出された上記ジッターパルス列を上記第二の残りのパ
ルス列から分離するジッターパルス列分離工程と、を備
えたことを特徴とするパルス列分離方法。 - 【請求項６】 上記固定ＰＲＩパルス列検出工程におい
て、ＰＲＩ変換を用いて検出を行うことを特徴とする請
求項５記載のパルス列分離方法。 - 【請求項７】 上記スタガーパルス列検出工程が、等間
隔に並ぶ３つのパルスからなる３連パルスの集合に対し
てそのパルス繰り返し間隔のヒストグラムを求める３連
パルスヒストグラム演算ステップと、 得られた上記ヒストグラムから複数のピークを検出し、
各ピークがパルス列の基本周期であるか否かを判定する
判定ステップと、 を備えたことを特徴とする請求項５または６に記載のパ
ルス列分離方法。 - 【請求項８】 上記ジッターパルス列検出工程におい
て、オーバーラップしたＰＲＩビンを持ち、時間原点の
移動を行う、変形ＰＲＩ変換を用いることを特徴とする
請求項５ないし７のいずれかに記載のパルス列分離方
法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP25162399A JP3559205B2 (ja) | 1999-09-06 | 1999-09-06 | パルス列分離装置及びパルス列分離方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP25162399A JP3559205B2 (ja) | 1999-09-06 | 1999-09-06 | パルス列分離装置及びパルス列分離方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2001077677A true JP2001077677A (ja) | 2001-03-23 |
| JP3559205B2 JP3559205B2 (ja) | 2004-08-25 |
Family
ID=17225589
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP25162399A Expired - Fee Related JP3559205B2 (ja) | 1999-09-06 | 1999-09-06 | パルス列分離装置及びパルス列分離方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3559205B2 (ja) |
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007163406A (ja) * | 2005-12-16 | 2007-06-28 | Mitsubishi Electric Corp | レーダ受信パルス分離装置 |
| JP2008281517A (ja) * | 2007-05-14 | 2008-11-20 | Mitsubishi Electric Corp | パルス変調信号列の同定装置 |
| JP4855473B2 (ja) * | 2005-10-24 | 2012-01-18 | ミツビシ・エレクトリック・アールアンドディー・センター・ヨーロッパ・ビーヴィ | パルス列の解析 |
| JP2012107925A (ja) * | 2010-11-16 | 2012-06-07 | Mitsubishi Electric Corp | レーダ装置 |
| JP2013029420A (ja) * | 2011-07-28 | 2013-02-07 | Mitsubishi Electric Corp | パッシブレーダ装置 |
| JP2013213757A (ja) * | 2012-04-03 | 2013-10-17 | Mitsubishi Electric Corp | パルス内変調分析装置およびパルス内変調分析方法 |
| KR101357689B1 (ko) | 2012-09-04 | 2014-02-04 | 국방과학연구소 | 레이더 신호 펄스 열 분리 방법 |
| JP2015187575A (ja) * | 2014-03-27 | 2015-10-29 | 三菱電機株式会社 | レーダ電波識別装置、レーダ電波識別方法及びプログラム |
| KR101620453B1 (ko) * | 2014-11-18 | 2016-05-12 | 국방과학연구소 | 레이더 신호 분석 방법 |
| KR101640790B1 (ko) * | 2016-02-05 | 2016-07-19 | 국방과학연구소 | 모노펄스 안테나의 표적 신호 수신 장치 및 그 수신 방법 |
| JP2020024113A (ja) * | 2018-08-06 | 2020-02-13 | 株式会社東芝 | パルス間隔分析装置、パルス間隔分析方法、およびプログラム |
| CN111781564A (zh) * | 2020-06-12 | 2020-10-16 | 中国船舶重工集团公司第七二四研究所 | 一种抖动信号脉冲序列快速预判方法 |
| JP7426717B2 (ja) | 2019-09-12 | 2024-02-02 | セルスター工業株式会社 | 探知機 |
| JP7470955B2 (ja) | 2019-09-12 | 2024-04-19 | セルスター工業株式会社 | 探知機 |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105403863B (zh) * | 2015-12-09 | 2018-03-30 | 哈尔滨工业大学 | 基于多特征提取的雷达脉冲重复间隔调制识别方法 |
-
1999
- 1999-09-06 JP JP25162399A patent/JP3559205B2/ja not_active Expired - Fee Related
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4855473B2 (ja) * | 2005-10-24 | 2012-01-18 | ミツビシ・エレクトリック・アールアンドディー・センター・ヨーロッパ・ビーヴィ | パルス列の解析 |
| JP2007163406A (ja) * | 2005-12-16 | 2007-06-28 | Mitsubishi Electric Corp | レーダ受信パルス分離装置 |
| JP2008281517A (ja) * | 2007-05-14 | 2008-11-20 | Mitsubishi Electric Corp | パルス変調信号列の同定装置 |
| JP2012107925A (ja) * | 2010-11-16 | 2012-06-07 | Mitsubishi Electric Corp | レーダ装置 |
| JP2013029420A (ja) * | 2011-07-28 | 2013-02-07 | Mitsubishi Electric Corp | パッシブレーダ装置 |
| JP2013213757A (ja) * | 2012-04-03 | 2013-10-17 | Mitsubishi Electric Corp | パルス内変調分析装置およびパルス内変調分析方法 |
| KR101357689B1 (ko) | 2012-09-04 | 2014-02-04 | 국방과학연구소 | 레이더 신호 펄스 열 분리 방법 |
| JP2015187575A (ja) * | 2014-03-27 | 2015-10-29 | 三菱電機株式会社 | レーダ電波識別装置、レーダ電波識別方法及びプログラム |
| KR101620453B1 (ko) * | 2014-11-18 | 2016-05-12 | 국방과학연구소 | 레이더 신호 분석 방법 |
| KR101640790B1 (ko) * | 2016-02-05 | 2016-07-19 | 국방과학연구소 | 모노펄스 안테나의 표적 신호 수신 장치 및 그 수신 방법 |
| JP2020024113A (ja) * | 2018-08-06 | 2020-02-13 | 株式会社東芝 | パルス間隔分析装置、パルス間隔分析方法、およびプログラム |
| JP7426717B2 (ja) | 2019-09-12 | 2024-02-02 | セルスター工業株式会社 | 探知機 |
| JP7470955B2 (ja) | 2019-09-12 | 2024-04-19 | セルスター工業株式会社 | 探知機 |
| CN111781564A (zh) * | 2020-06-12 | 2020-10-16 | 中国船舶重工集团公司第七二四研究所 | 一种抖动信号脉冲序列快速预判方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP3559205B2 (ja) | 2004-08-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP2001077677A (ja) | パルス列分離装置及びパルス列分離方法 | |
| JP5148457B2 (ja) | 異常判定装置、方法、及びプログラム | |
| KR100677574B1 (ko) | 비디오 시퀀스로부터 노이즈 분포를 추정하는 방법 | |
| EP2081326A1 (en) | Statistical processing apparatus capable of reducing storage space for storing statistical occurence frequency data and a processing method therefor | |
| US8050910B2 (en) | Pitch extraction with inhibition of harmonics and sub-harmonics of the fundamental frequency | |
| US9192311B2 (en) | Pulse wave detection device and pulse wave detection method | |
| JP2009210593A (ja) | 特徴量抽出装置、方法及びプログラム | |
| EP2021979B1 (en) | Pattern matching | |
| CN116150592A (zh) | 基于人工智能的轮椅智能监测预警方法及系统 | |
| JP3733057B2 (ja) | 製造プロセスにおける操業解析装置、その方法、及びコンピュータ読み取り可能な記憶媒体 | |
| KR102461702B1 (ko) | 심전도 신호 처리 방법 | |
| KR20080069392A (ko) | 심전도 신호의 알 피크 검출 방법 | |
| JP5094937B2 (ja) | 周波数解析装置 | |
| KR101081050B1 (ko) | 비음수 행렬 인수분해에 기초한 목표 신호 검출 방법 및 시스템 | |
| JP5094281B2 (ja) | 信号分離装置 | |
| KR102367639B1 (ko) | 심전도 신호 처리 방법 | |
| JP2015096831A (ja) | 情報処理装置、情報処理方法及びプログラム | |
| CN113608190A (zh) | 基于奇异空间三特征的海面目标检测方法及系统 | |
| JP2000046893A (ja) | 異常診断装置および異常診断方法 | |
| JP2003133922A (ja) | パルス列検出装置及びパルス列検出方法 | |
| CN112948745A (zh) | 一种基于Jerk模型的自校正Kalman高机动目标跟踪滤波器 | |
| JP2002358526A (ja) | 映像オブジェクト検出・追跡装置 | |
| JP2004020427A (ja) | ノイズ除去方法及びノイズ除去フィルタ | |
| JP6932898B1 (ja) | 信号判定装置及びプログラム | |
| JPH0643893A (ja) | 音声認識方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20040302 |
|
| 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: 20040518 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20040520 |
|
| 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: 20080528 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090528 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100528 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100528 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110528 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110528 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120528 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120528 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130528 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140528 Year of fee payment: 10 |
|
| LAPS | Cancellation because of no payment of annual fees |