(57)【特許請求の範囲】
【請求項１】
複数の信号源からのパルス列が重なり合った複合パルス列の各パルスの到着時間（ＴＯＡ）データに基づき、固定パルス繰り返し間隔（ＰＲＩ）を有する固定ＰＲＩパルス列を検出する固定ＰＲＩパルス列検出手段と、
検出された上記固定ＰＲＩパルス列をもとの上記複合パルス列から分離するとともに、分離後の残りのパルス列を出力する固定ＰＲＩパルス列分離手段と、
上記固定ＰＲＩパルス列分離手段から出力される上記残りのパルス列から、上記複合パルス列の各パルスの到着時間（ＴＯＡ）データに基づき、スタガーパルス列を検出するスタガーパルス列検出手段と、
検出されたスタガーパルス列を上記残りのパルス列から分離するとともに、分離後の固定ＰＲＩパルス列及びスタガーパルス列を含まない第二の残りのパルス列を出力するスタガーパルス列分離手段と、
上記スタガーパルス列分離手段から出力される上記第二の残りのパルス列から、上記複合パルス列の各パルスの到着時間（ＴＯＡ）データに基づき、ジッターパルス列を検出するジッターパルス列検出手段と、
検出された上記ジッターパルス列を上記第二の残りのパルス列から分離するジッターパルス列分離手段と、
を備えたことを特徴とするパルス列分離装置。
【請求項２】
上記固定ＰＲＩパルス列検出手段が、ＰＲＩ変換を用いて検出を行うことを特徴とする請求項１記載のパルス列分離装置。
【請求項３】
上記スタガーパルス列検出手段が、
等間隔に並ぶ３つのパルスからなる３連パルスの集合に対してそのパルス繰り返し間隔のヒストグラムを求める３連パルスヒストグラム演算部と、
得られた上記ヒストグラムから複数のピークを検出し、各ピークがパルス列の基本周期であるか否かを判定する判定部と、
を備えたことを特徴とする請求項１または２に記載のパルス列分離装置。
【請求項４】
上記ジッターパルス列検出手段が、オーバーラップしたＰＲＩビンを持ち、時間原点の移動を行う、変形ＰＲＩ変換を用いることを特徴とする請求項１ないし３のいずれかに記載のパルス列分離装置。
【請求項５】
複数の信号源からのパルス列が重なり合った複合パルス列の各パルスの到着時間（ＴＯＡ）データに基づき、固定パルス繰り返し間隔（ＰＲＩ）を有する固定ＰＲＩパルス列を検出する固定ＰＲＩパルス列検出工程と、
検出された上記固定ＰＲＩパルス列をもとの上記複合パルス列から分離するとともに、分離後の残りのパルス列を出力する固定ＰＲＩパルス列分離工程と、
上記固定ＰＲＩパルス列分離手段から出力される上記残りのパルス列から、上記複合パルス列の各パルスの到着時間（ＴＯＡ）データに基づき、スタガーパルス列を検出するスタガーパルス列検出工程と、
検出された上記スタガーパルス列を上記残りのパルス列から分離するとともに、分離後の固定ＰＲＩパルス列及びスタガーパルス列を含まない第二の残りのパルス列を出力するスタガーパルス列分離工程と、
上記スタガーパルス列分離手段から出力される上記第二の残りのパルス列から、上記複合パルス列の各パルスの到着時間（ＴＯＡ）データに基づき、ジッターパルス列を検出するジッターパルス列検出工程と、
検出された上記ジッターパルス列を上記第二の残りのパルス列から分離するジッターパルス列分離工程と、
を備えたことを特徴とするパルス列分離方法。
【請求項６】
上記固定ＰＲＩパルス列検出工程において、ＰＲＩ変換を用いて検出を行うことを特徴とする請求項５記載のパルス列分離方法。
【請求項７】
上記スタガーパルス列検出工程が、
等間隔に並ぶ３つのパルスからなる３連パルスの集合に対してそのパルス繰り返し間隔のヒストグラムを求める３連パルスヒストグラム演算ステップと、
得られた上記ヒストグラムから複数のピークを検出し、各ピークがパルス列の基本周期であるか否かを判定する判定ステップと、
を備えたことを特徴とする請求項５または６に記載のパルス列分離方法。
【請求項８】
上記ジッターパルス列検出工程において、オーバーラップしたＰＲＩビンを持ち、時間原点の移動を行う、変形ＰＲＩ変換を用いることを特徴とする請求項５ないし７のいずれかに記載のパルス列分離方法。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
この発明は、複数の信号源からのパルス列が重なり合った複合信号を、信号源ごとのパルス列に分離するパルス列分離装置及びパルス列分離方法に関するものである。
【０００２】
【従来の技術】
従来から、複数の信号源からのパルス列が重なり合った複合信号を、信号源ごとのパルス列に分離するために、パルスの到着時間（以下、ＴＯＡ（Ｔｉｍｅ  ｏｆ  Ａｒｒｉｖａｌ）とする）データに基づくパルス列分離装置が提案されてきた。
【０００３】
図１６は、例えば米国特許第５、５８３，５０５号公報に示された、パルスのＴＯＡデータからパルス列を検出し分離するための従来のパルス列分離装置における処理を示したフローチャートである。
【０００４】
図１６において、Ｓ６０は自己相関を求めるステップ、Ｓ６２はパルス繰り返し間隔（以下、ＰＲＩ（Ｐｕｌｓｅ  Ｒｅｐｅｔｉｔｉｏｎ  Ｉｎｔｅｒｖａｌ）とする）の候補を求めるステップ、Ｓ６４は積分器にメモリーを割り当てるステップ、Ｓ６６は各パルスをすべての積分器のビンに積み重ねるステップ、Ｓ６８は最小ＰＲＩの積分器を分析するステップ、Ｓ７０は積分器のピークとしきい値を比較するステップ、Ｓ７２は積分器のピークがしきい値を越えているとき該当するＰＲＩが存在すると判定するステップ、Ｓ７４はそのＰＲＩを有するパルス列を時間窓から除去するステップ、Ｓ７６は自己相関を再度求めるステップ、Ｓ７７は次のＰＲＩの候補の積分器に進むステップ、Ｓ７８は分析すべき積分器がまだ残っているか否か判定するステップ、Ｓ７９は残ったパルスが単独のものであるか否かを判定するステップである。
【０００５】
次に、前述した従来のパルス列分離装置の動作について図面を参照しながら説明する。
【０００６】
図１６において、まず時間窓内にあるパルスのＴＯＡデータの自己相関関数を求める（ステップＳ６０）。得られた自己相関関数のピークの判定を行ってピークの面積（ピークの高さと幅の積）がしきい値を超えていれば、そのピークの位置をＰＲＩの候補として抽出する（ステップＳ６２）。次に、これらＰＲＩの候補ごとに積分器を割り当てる（ステップＳ６４）。積分器はＰＲＩ分の時間幅をもつテーブルから構成される。積分器では、入力した各パルスに対してそのＴＯＡをＰＲＩで割った余りの時間に対応するテーブル内の位置を求め、テーブルの値をカウントアップしていく。時間窓内のすべてのパルスに対して、各積分器で積分を行う（ステップＳ６６）。
【０００７】
次いで信号源の検出を行うために、小さいＰＲＩの積分器から順番に分析していく（ステップＳ６８）。積分器のテーブルのピークの値がしきい値を超えているか否かを調べ（ステップＳ７０）、超えていればそのＰＲＩを持つ信号源が存在すると判定する（ステップＳ７２）。そのあとそのＰＲＩをもつパルス列をもとのパルス列から除去し（ステップＳ７４）、残ったパルス列に対して再度自己相関関数を求めて（ステップＳ７６）、上記の処理を繰り返す。
【０００８】
一方、ステップＳ７０において、積分器のテーブルのピークの値がしきい値を超えていなければ、次に小さいＰＲＩに対する積分器の分析に進む（ステップＳ７７）。これは分析する積分器が無くなるまで続けられる（ステップＳ７８）。最後に、残ったパルスが単独のものであるか否かを調べ、単独であれば処理を終了する。単独のものでなくグループを構成していればレーダスキャンによって生じたものとみなし、スキャン周期の分析を行う（ステップＳ７９）。
【０００９】
【発明が解決しようとする課題】
上記のような、従来のパルス列分離装置では、まず自己相関関数でＰＲＩの候補を検出することが処理の前提になっているが、信号源の数が多い場合や、固定ＰＲＩ以外の複雑なＰＲＩのパターンを有する信号源が存在する場合には、自己相関関数に明確なピークが現われにくく、また、自己相関関数にはＰＲＩの整数倍のところにできるハーモニックスの問題があり、ＰＲＩの候補が適切に抽出できず、パルス列を正確に検出し分離できないという問題点があった。
【００１０】
本発明は上記のような課題を解決するためになされたもので、複雑なＰＲＩのパターンを有する信号源も含めて信号源のＰＲＩを正確に推定し、パルス列の分離精度を向上させることが可能なパルス列分離装置及びパルス列分離方法を得ることを目的とする。
【００１１】
【課題を解決するための手段】
この発明は、複数の信号源からのパルス列が重なり合った複合パルス列の各パルスの到着時間（ＴＯＡ）データに基づき、固定パルス繰り返し間隔（ＰＲＩ）を有する固定ＰＲＩパルス列を検出する固定ＰＲＩパルス列検出手段と、検出された固定ＰＲＩパルス列をもとの複合パルス列から分離するとともに、分離後の残りのパルス列を出力する固定ＰＲＩパルス列分離手段と、固定ＰＲＩパルス列分離手段から出力される残りのパルス列から、複合パルス列の各パルスの到着時間（ＴＯＡ）データに基づき、スタガーパルス列を検出するスタガーパルス列検出手段と、検出されたスタガーパルス列を残りのパルス列から分離するとともに、分離後の固定ＰＲＩパルス列及びスタガーパルス列を含まない第二の残りのパルス列を出力するスタガーパルス列分離手段と、スタガーパルス列分離手段から出力される第二の残りのパルス列から、複合パルス列の各パルスの到着時間（ＴＯＡ）データに基づき、ジッターパルス列を検出するジッターパルス列検出手段と、検出されたジッターパルス列を第二の残りのパルス列から分離するジッターパルス列分離手段と、を備えたパルス列分離装置である。
【００１２】
また、固定ＰＲＩパルス列検出手段が、ＰＲＩ変換を用いて検出を行う。
【００１３】
また、スタガーパルス列検出手段が、等間隔に並ぶ３つのパルスからなる３連パルスの集合に対してそのパルス繰り返し間隔のヒストグラムを求める３連パルスヒストグラム演算部と、得られたヒストグラムから複数のピークを検出し、各ピークがパルス列の基本周期であるか否かを判定する判定部と、を備えている。
【００１４】
また、ジッターパルス列検出手段が、オーバーラップしたＰＲＩビンを持ち、時間原点の移動を行う、変形ＰＲＩ変換を用いる。
【００１５】
また、この発明は、複数の信号源からのパルス列が重なり合った複合パルス列の各パルスの到着時間（ＴＯＡ）データに基づき、固定パルス繰り返し間隔（ＰＲＩ）を有する固定ＰＲＩパルス列を検出する固定ＰＲＩパルス列検出工程と、検出された固定ＰＲＩパルス列をもとの複合パルス列から分離するとともに、分離後の残りのパルス列を出力する固定ＰＲＩパルス列分離工程と、固定ＰＲＩパルス列分離手段から出力される残りのパルス列から、複合パルス列の各パルスの到着時間（ＴＯＡ）データに基づき、スタガーパルス列を検出するスタガーパルス列検出工程と、検出されたスタガーパルス列を残りのパルス列から分離するとともに、分離後の固定ＰＲＩパルス列及びスタガーパルス列を含まない第二の残りのパルス列を出力するスタガーパルス列分離工程と、スタガーパルス列分離手段から出力される第二の残りのパルス列から、複合パルス列の各パルスの到着時間（ＴＯＡ）データに基づき、ジッターパルス列を検出するジッターパルス列検出工程と、検出されたジッターパルス列を第二の残りのパルス列から分離するジッターパルス列分離工程と、を備えたパルス列分離方法である。
【００１６】
また、固定ＰＲＩパルス列検出工程において、ＰＲＩ変換を用いて検出を行う。
【００１７】
また、スタガーパルス列検出工程が、等間隔に並ぶ３つのパルスからなる３連パルスの集合に対してそのパルス繰り返し間隔のヒストグラムを求める３連パルスヒストグラム演算ステップと、得られたヒストグラムから複数のピークを検出し、各ピークがパルス列の基本周期であるか否かを判定する判定ステップと、を備えている。
【００１８】
また、ジッターパルス列検出工程において、オーバーラップしたＰＲＩビンを持ち、時間原点の移動を行う、変形ＰＲＩ変換を用いる。
【００１９】
【発明の実施の形態】
実施の形態１．
以下、この発明の一実施の形態を図１乃至図１５について説明する。
【００２０】
図１は本発明のパルス列分離装置が処理の対象とするパルス列の種類を説明するための説明図であり、図２は本発明のパルス列分離装置の機能ブロック図、図３は固定ＰＲＩパルス列検出・分離手段のブロック図、図４は固定ＰＲＩパルス列検出・分離手段の中で用いるＰＲＩ変換に関するＰＲＩビンの説明図、図５は固定ＰＲＩパルス列検出・分離手段の中で用いるＰＲＩ変換のフローチャート、図６はスタガーパルス列検出・分離手段のブロック図、図７はスタガーパルス列検出・分離手段の中で用いる３連パルスのヒストグラムのフローチャート、図８はスタガーパルス列検出・分離手段の中で用いる３連パルスのヒストグラムからの基本周期検出のフローチャート、図９はジッターパルス列検出・分離手段のブロック図、図１０はジッターパルス列検出・分離手段の中で用いる変形ＰＲＩ変換に関するＰＲＩビンの説明図、図１１はジッターパルス列検出・分離手段の中で用いる変形  ＰＲＩ変換のフローチャート、図１２は入力パルス列の説明図、図１３は固定ＰＲＩパルス列検出・分離の説明図、図１４はスタガーパルス列検出・分離の説明図、図１５はジッターパルス列検出・分離の説明図である。
【００２１】
まず、本発明のパルス列分離装置が対象とするパルス列の種類を図１を用いて説明する。図１において、１１はＰＲＩが一定の固定ＰＲＩパルス列である。パルス列には固定ＰＲＩパルス列１１以外に、１２に示すスタガーパルス列と１３に示すジッターパルス列がある。スタガーパルス列１２は、一定のフレーム周期１２ａの中に、複数の不等間隔（図１では、ＰＲＩ１，ＰＲＩ２，ＰＲＩ３の３つ）のパルスが存在し、それが周期１２ａごとに繰り返すパルス列である。ジッターパルス列１３は、ＰＲＩが平均ＰＲＩのまわりにランダムに揺らぐパルス列である。本発明のパルス列分離装置への入力パルス列は、固定ＰＲＩパルス列１１、スタガーパルス列１２、ジッターパルス列１３がそれぞれ複数重なり合った複合パルス列である。
【００２２】
次に、本発明の構成及び動作について説明する。
【００２３】
図２において、１００は、信号源から送信されたパルスを受信し、当該パルスのＴＯＡデータを出力する広帯域受信機、２００は本発明のパルス列分離装置に設けられた入力端、３００は固定ＰＲＩパルス列の検出・分離を行うための固定ＰＲＩパルス列検出・分離手段、４００は、スタガーパルス列の検出・分離を行うためのスタガーパルス列検出・分離手段、５００は、ジッターパルス列の検出・分離を行うためのジッターパルス列検出・分離手段、６００は本発明のパルス列分離装置に設けられた出力端である。
【００２４】
動作について説明する。広帯域受信機１００からのパルスのＴＯＡデータは入力端２００へ入力される。これらＴＯＡデータはまず固定ＰＲＩパルス列検出・分離手段３００で検出・分離処理が行われる。
【００２５】
図３には、図２の固定ＰＲＩパルス列検出・分離手段３００の内部構成が示されている。図において、３１０は入力端、３３０はＰＲＩ変換を行うＰＲＩ変換演算部、３５０は適当なしきい値の値を設定するとともに検出可能範囲を決定するしきい値演算部、３６０は固定ＰＲＩパルス列のＰＲＩを抽出するＰＲＩ抽出部、３７０は抽出されたＰＲＩの値を利用して、もとのパルス列から該当する部分パルス列を分離する固定ＰＲＩパルス列分離演算部、３８０は固定ＰＲＩパルス列が残存するか否かを判定する残存パルス列判定部、３９０は出力端である。
【００２６】
動作について説明する。まず、入力端３１０へ入力されたパルス列のＴＯＡデータにはＰＲＩ変換演算部３３０においてＰＲＩ変換が施される。このＰＲＩ変換は、例えば、西口（三菱電機（株））、’’パルス繰り返し間隔推定のための新手法’’、昭和５８年度電子情報通信学会情報・システム部門全国大会、に示されている方法を用いる。この方法は、まず、入力データのＴＯＡを  ｔ１、ｔ２、・・・、ｔＮ  とし、入力パルス列をデルタ関数の和で次式のように表わす。
【００２７】
【数１】
0003559205.tif 000002
【００２８】
ＰＲＩ変換は関数  ｇ（ｔ）  に対して次式で定義されるような変換である。
【００２９】
【数２】
0003559205.tif 000003
【００３０】
これは上式（１）のパルス列に対しては以下のような形になる。
【００３１】
【数３】
0003559205.tif 000004
【００３２】
このＰＲＩ変換は、連続時間での表現であるが、実際の計算機で行う離散時間でのＰＲＩ変換は後でフローチャートを用いて詳述する。このＰＲＩ変換を入力パルス列に施すと、入力パルス列に含まれている部分パルス列のＰＲＩのところにだけ鋭いピークがたつような変換結果が得られる。さらに、しきい値演算部  ３５０でＰＲＩ検出のための適当なしきい値を定める。しきい値演算部３５０では、同時に検出可能なＰＲＩの範囲も求める。次に、ＰＲＩ抽出部３６０で、そのしきい値を越えるピーク値を検出するようにすれば、固定ＰＲＩパルス列の検出可能な範囲にあるＰＲＩがすべて検出される。このようにして検出されたＰＲＩの値を利用して、固定ＰＲＩパルス列分離演算部３７０で、もとのパルス列から該当する部分パルス列を分離する。この分離のためには、例えば、特公平１−４７９３６号公報に詳しく示されているパルス繰り返し間隔フィルターを用いる。残存パルス列判定部３８０では、調べたいＰＲＩの範囲と検出可能なＰＲＩの範囲を比較し、もし検出可能なＰＲＩの範囲の方が狭ければ、すなわち、固定ＰＲＩパルス列が残存する可能性があれば、残ったパルス列に再びＰＲＩ変換演算部３３０でＰＲＩ変換を施し処理を繰り返す。そうでなければ、残ったパルス列と共に出力端３９０へ進む。
【００３３】
次に、実際の計算機で行う離散時間でのＰＲＩ変換について説明する。図４は、ＰＲＩ変換で用いるＰＲＩビンの説明図である。調べたいＰＲＩの下限と上限をそれぞれ、τｍｉｎ、τｍａｘ  と置く。ＰＲＩ軸での区間［τｍｉｎ，τｍａｘ］  を  Ｋ  個に等分割し、分割された細分区間をＰＲＩビンと呼ぶ。各ビンの幅は次式（４）で与えられる。
【００３４】
【数４】
0003559205.tif 000005
【００３５】
また、中心ＰＲＩは次式（５）で与えられる。
【００３６】
【数５】
0003559205.tif 000006
【００３７】
図５にはＰＲＩ変換の処理のフローチャートを示す。入力データのパルスの  ＴＯＡを  ｔ１、  ｔ２、  …、  ｔＮ  とする。ＰＲＩ変換の処理は以下の手順で行う。
【００３８】
（１）初期設定を行う（ステップＳ３３１）。すなわち、ＰＲＩ変換  Ｄｋ（ｋ  ＝１，  ２，  …  ，  Ｋ）をゼロクリアするとともに、パルスペアの右パルスのパルス番号を  ｍ  ＝  １  に設定する。
【００３９】
（２）パルスペアの右パルスのパルス番号を更新する（ステップＳ３３２）。すなわち、ｍ  ＝  ｍ  ＋  １を行う。
【００４０】
（３）終了判定を行う（ステップＳ３３３）。すなわち、ｍ≦  ＮならばステップＳ３３４へ。そうでなければ終了する。
【００４１】
（４）パルスペアの左パルスのパルス番号をｎ  ＝  ｍ  に設定する（ステップＳ３３４）。
【００４２】
（５）パルスペアの左パルスのパルス番号を更新する（ステップＳ３３５）。すなわち、ｎ  ＝ｎ−１を行う。
【００４３】
（６）パルスペアの左パルスの終了判定を行う（ステップＳ３３６）。すなわち、ｎ  ≧  １  ならば、ステップＳ３３７へ。そうでなければステップＳ３３２へ。
【００４４】
（７）パルスペアの時間間隔を求める（ステップＳ３３７）。すなわち、τ＝ｔｍ−ｔｎを行う。
【００４５】
（８）時間間隔の大きさ判定１を行う（ステップＳ３３８）。すなわち、τ≧τｍｉｎ  ならばステップＳ３３９へ。そうでなければステップＳ３３５へ。
【００４６】
（９）時間間隔の大きさ判定２を行う（ステップＳ３３９）。すなわち、τ≦τｍａｘ  ならばステップＳ３４０へ。そうでなければＳ３３２へ。
【００４７】
（１０）τ軸上でのＰＲＩビンの選出（ステップＳ３４０）。すなわち、ＰＲＩビンの番号ｋを次式（６）で計算する。
【００４８】
【数６】
0003559205.tif 000007
【００４９】
ここで、［・］  は  Ｇａｕｓｓの整数、すなわち、（・）  を越えない最大の整数を表わす。
【００５０】
（１１）ＰＲＩ変換の更新を行う（ステップＳ３４１）。すなわち、次式（７）でＰＲＩ変換の値を更新する。このあと、ステップＳ３３５へ戻る。
【００５１】
【数７】
0003559205.tif 000008
【００５２】
図３のしきい値演算部３５０で求めるしきい値は、例えば次のような式にする。
【００５３】
【数８】
0003559205.tif 000009
【００５４】
ここで、Ｔ  は観測時間で、最初と最後のパルスのＴＯＡの差から求める。αは調整可能なパラメータで、例えば  ０．３  くらいの値にする。上のしきい値の式（８）は、もし、全観測時間にわたってＰＲＩがτｋのパルス列が存在すればＰＲＩ変換の価は理想的には、次式（９）となることに基づく。
【００５５】
【数９】
0003559205.tif 000010
【００５６】
しきい値はτｋの減少関数であるが、一方、ＰＲＩ変換のノイズフロアーはτｋと共に減少はしない。したがって、τｋ  がある大きさ（これをτｌｉｍ  とおく）になるとそこからはしきい値よりもノイズフロアーの値の方が大きな値をとってＰＲＩの検出は不能になる。すなわち、［τｍｉｎ，τｌｉｍ］  が検出可能範囲になる。
【００５７】
以上のようにして図２の固定ＰＲＩパルス列検出・分離手段３００で固定ＰＲＩパルス列がもとのパルス列から除去されると、次に、残ったパルス列に対して、スタガーパルス列検出・分離手段４００でスタガーパルス列の検出・分離処理が行われる。
【００５８】
図６にスタガーパルス列検出・分離手段の内部構成を示したブロック図を示す。図において、４１０は入力端、４３０は３連パルスのヒストグラムを求める３連パルスヒストグラム演算部、４５０は適当なしきい値の値を設定するとともに検出可能範囲を決定するしきい値演算部、４６０はスタガーパルス列のＰＲＩを抽出するＰＲＩ抽出部、４７０は抽出されたＰＲＩの値を利用して、もとのパルス列から該当する部分パルス列を分離するスタガーパルス列分離演算部、４８０はスタガーパルス列が残存するか否かを判定する残存パルス列判定部、４９０は出力端である。
【００５９】
なお、しきい値演算部３５０及びＰＲＩ抽出部３６０は、３連パルスヒストグラム演算部３３０において得られたヒストグラムから複数のピークを検出し、各ピークがパルス列の基本周期であるか否かを判定するための判定手段を構成している。
【００６０】
次に、動作について説明する。入力端４１０に入力したパルス列は、３連パルスヒストグラム演算部４３０で３連パルスのヒストグラムが計算される。３連パルスのヒストグラムは、図１で説明したようなほぼ等間隔のパルスの３つ組の個数をカウントするヒストグラムである。このヒストグラムを作成する処理はあとで詳述するが、得られるヒストグラムでは、スタガーパルス列の基本周期（フレーム周期１２ａ）とその整数倍のところにするどいピークがたつ。しきい値演算部４５０では  スタガーのＰＲＩ（基本周期）検出のための適当なしきい値を定める。しきい値演算部４５０では同時に、検出可能なＰＲＩの範囲も求める。次に、ＰＲＩ抽出部４６０で、そのしきい値を超える値を検出するようにすれば、スタガーパルス列の検出可能な範囲にあるＰＲＩがすべて検出される。ＰＲＩ抽出部４６０で行う３連パルスのヒストグラムからのＰＲＩ（基本周期）の検出についてはあとでフローチャートを用いて詳述する。このようにして検出されたＰＲＩの値を利用して、スタガーパルス列分離演算部４７０で、もとのパルス列から該当する部分パルス列を分離する。この分離のためには、例えば特公平１−４７９３６号公報に詳しく示されているパルス繰り返し間隔フィルターを用いる。残存パルス列判定部４８０では、調べたいＰＲＩの範囲と検出可能なＰＲＩの範囲を比較し、もし検出可能なＰＲＩの範囲の方が狭ければ、すなわち、スタガーパルス列が残存する可能性があれば、残ったパルス列に再び３連パルスヒストグラム演算部４３０でヒストグラムを作り、処理を繰り返す。そうでなければ、残ったパルス列と共に出力端４９０へ進む。
【００６１】
図７に、３連パルスヒストグラム演算部４３０で行う、３連パルスのヒストグラム演算の処理の流れを示すフローチャートを示す。このヒストグラムで用いるＰＲＩビンはＰＲＩ変換の場合の図４と同一である。入力データのパルスのＴＯＡを  ｔ１、  ｔ２、  …、  ｔＮ  とする。３連パルスのヒストグラムの作成は以下の手順で行う。
【００６２】
（１）まずはじめに、初期設定（ステップＳ４３１）を行う。すなわち、ヒストグラム  Ｅｋ（ｋ＝  １，  ２，  …  ，  Ｋ）をゼロクリアする。また、３連パルスの中心パルスのパルス番号を  ｍ  ＝  １  に設定する。
【００６３】
（２）３連パルスの中心パルスのパルス番号を更新する（ステップＳ４３２）。すなわち、ｍ  ＝  ｍ  ＋  １を行う。
【００６４】
（３）終了判定を行う（ステップＳ４３３）。すなわち、ｍ≦Ｎ−１ならばステップＳ４３４へ。そうでなければ終了。
【００６５】
（４）３連パルスの右パルスのパルス番号を  Ｆ  ＝  ｍ  ＋  １  に、左パルスのパルス番号をＢ  ＝  ｍ  に設定する（ステップＳ４３４）。
【００６６】
（５）３連パルスの左パルスのパルス番号を更新する（ステップＳ４３５）。すなわち、Ｂ  ＝  Ｂ−１を行う。
【００６７】
（６）３連パルスの左パルスの終了判定を行う（ステップＳ４３６）。すなわち、Ｂ  ≧１  ならば、ステップＳ４３７へ。そうでなければステップＳ４３２へ。
【００６８】
（７）３連パルスの左ペア（後ろ向きパルスペア）の時間間隔を求める（ステップＳ４３７）。すなわち、τ＝  ｔｍ−ｔＢ。
【００６９】
（８）時間間隔の大きさ判定１（ステップＳ４３８）を行う。すなわち、τ≧τｍｉｎ  ならばステップＳ４３９へ。そうでなければステップＳ４３５へ。
【００７０】
（９）時間間隔の大きさ判定２（ステップＳ４３９）を行う。すなわち、τ≦τｍａｘ  ならばステップＳ４４０へ。そうでなければＳ４３２へ。
【００７１】
（１０）３連パルスの右ペア（前向きパルスペア）の時間間隔を求める（ステップＳ４４０）。すなわち、τ’＝  ｔＦ−ｔｍを行う。
【００７２】
（１１）次に、時間間隔の大きさ判定３（ステップＳ４４１）を行う。すなわち、τ’≧τ−Δτ／２ならばステップＳ４４５へ。そうでなければステップＳ４４２へ。
【００７３】
（１２）３連パルスの右パルスの終了判定（ステップＳ４４２）を行う。すなわち、Ｆ  ＝  Ｎ  ならばステップＳ４３５へ戻る。そうでなければ右パルスを  Ｆ  ＝  Ｆ
＋  １  と更新（ステップＳ４４３）してステップＳ４４０へ。
【００７４】
（１３）次に、時間間隔の大きさ判定４（ステップＳ４４５）を行う。すなわち、τ’≦τ＋Δτ／２ならばステップＳ４４６へ。そうでなければＳ４３５へ。
【００７５】
（１４）τ軸上でのＰＲＩビンの選出を行う（ステップＳ４４６）。すなわち、ＰＲＩビンの番号  ｋ  を次式（１０）で計算する。
【００７６】
【数１０】
0003559205.tif 000011
【００７７】
ここで、［・］  は  Ｇａｕｓｓ  の整数、すなわち、（・）  を越えない最大の整数を表わす。
【００７８】
（１５）ヒストグラムの更新を行う（ステップＳ４４７）。すなわち、次式（１１）で３連パルスのヒストグラムの値を更新する。そのあと、ステップＳ４３５へ。
【００７９】
【数１１】
0003559205.tif 000012
【００８０】
次に、図８に、上述のようにして求めた３連パルスのヒストグラムからＰＲＩ（基本周期）を検出する処理のフローチャートを示す。検出は以下のような手順で行う。
【００８１】
（１）まず、しきい値の初期設定を行う（ステップＳ４６１）。すなわち、次式（１２）でしきい値を設定する。
【００８２】
【数１２】
0003559205.tif 000013
【００８３】
ここで、αはパラメータで、１．６  程度の値を用いる。また、Ｔ  は観測時間であり、最初のパルスと最後のパルスのＴＯＡの差から求める。
【００８４】
（２）ピークの検出を行う（ステップＳ４６２）。すなわち、Ｅｋ  ＝  １，  …，  Ｋ  のうち、Ｅｋ  ≧Ａｋ  、かつ  、Ｅｋ  ≧ｍａｘ  ｛Ｅｋ−１，Ｅｋ＋１｝（局所最大値）となるものを見つけ、それらを  Ｅｋｐ（ｐ  ＝  １，  …，  ｐｍａｘ）とおく。
【００８５】
（３）ピーク番号の初期設定を行う（ステップＳ４６３）。すなわち、  ｐ  ＝  ０とする。
【００８６】
（４）ピーク番号の更新を行う（ステップＳ４６４）。すなわち、ｐ  ＝  ｐ  ＋  １とする。
【００８７】
（５）終了判定を行う（ステップＳ４６５）。すなわち、ｐ  ≦  ｐｍａｘ  ならばステップＳ４６６へ。そうでなければ終了。
【００８８】
（６）しきい値の更新を行う（ステップＳ４６６）。すなわち、次式（１３）でしきい値  Ａｋ  の値を更新する。
【００８９】
【数１３】
0003559205.tif 000014
【００９０】
（７）ピークの判定を行う（ステップＳ４６７）。すなわち、ｐ番目のピーク値  Ｅｋｐ  が、更新されたしきい値Ａｋｐ  を上回っているか否かを判定する。Ｅｋｐ  ≧  Ａｋｐならば、ステップＳ４６８へ。そうでなければステップＳ４６４へ。
【００９１】
（８）次に、基本周期か否かの判定を行う（ステップＳ４６８）。すなわち、すでに検出されている基本周期のすべてについて、τｋｐ  がその整数倍になっているか否かの判定を行う。ある基本周期のほぼ整数倍（ｊ倍）になっていれば、τをτｋｐ／ｊ  で置き換えることにより、基本周期の精度を上げる。もし、どの基本周期の整数倍にもなっていなければ、τｋｐ  を新たな基本周期として付け加える。このあとステップＳ４６４へ戻る。
【００９２】
図２のスタガーパルス列検出・分離手段４００でスタガーパルス列がもとのパルス列から除去されると、残ったパルス列に対してジッターパルス列検出・分離手段５００でジッターパルス列の検出・分離処理を行う。
【００９３】
図９にジッターパルス列検出・分離手段５００の内部構成を示したブロック図を示す。図において、５１０は入力端、５３０は変形ＰＲＩ変換を行う変形ＰＲＩ変換演算部、５５０は適当なしきい値の値を設定するとともに検出可能範囲を決定するしきい値演算部、５６０はジッターパルス列のＰＲＩを抽出するＰＲＩ抽出部、５７０は抽出されたＰＲＩの値を利用して、もとのパルス列から該当する部分パルス列を分離するジッターパルス列分離演算部、５８０はジッターパルス列が残存するか否かを判定する残存パルス列判定部、５９０は出力端である。
【００９４】
動作について説明する。入力端５１０に入力したパルス列は、変形ＰＲＩ変換演算部５３０で変形ＰＲＩ変換が計算される。変形ＰＲＩ変換についてはあとで詳述するが、これから求められるＰＲＩスペクトラムには、ジッターパルス列の平均ＰＲＩのところにするどいピークがたつ。しきい値演算部５５０ではジッターパルス列のＰＲＩ検出のための適当なしきい値を定める。しきい値演算部５５０では同時に、検出可能なＰＲＩの範囲も求める。次に、ＰＲＩ抽出部５６０で、そのしきい値を超える値を検出するようにすれば、ジッターパルス列の検出可能な範囲にあるＰＲＩがすべて検出される。このようにして検出されたＰＲＩの値を利用して、ジッターパルス列分離演算部５７０で、もとのパルス列から該当する部分パルス列を分離する。この分離のためには、例えば特公平１−４７９３６号公報に詳しく示されているパルス繰り返し間隔フィルターを用いる。残存パルス列判定部５８０では、調べたいＰＲＩの範囲と検出可能なＰＲＩの範囲を比較し、もし検出可能なＰＲＩの範囲の方が狭ければ、すなわち、ジッターパルス列が残存する可能性があれば、残ったパルス列に再び変形ＰＲＩ変換演算部５３０でＰＲＩスペクトルを作り、処理を繰り返す。そうでなければ、出力端５９０へ進む。
【００９５】
ジッターパルス列に対しては、先に述べたＰＲＩ変換ではうまくＰＲＩのところにピークをたたせることができないが、これから述べる変形ＰＲＩ変換ではジッターパルス列であってもその平均ＰＲＩのところにピークがたつようにできる。
【００９６】
変形ＰＲＩ変換は、ＰＲＩ変換のＰＲＩビンをオーバーラップさせ、時間原点の移動を行うようにしたものである。
【００９７】
図１０は変形ＰＲＩ変換で用いるＰＲＩビンの説明図である。調べたいＰＲＩの下限と上限をそれぞれ、τｍｉｎ、τｍａｘ  と置く。  ＰＲＩ軸での区間  ［τｍｉｎ，τｍａｘ］  内のＰＲＩビンの個数をＫとし、ＰＲＩの揺らぎの大きさの上限をεとする（εは片側で表現する。例えばＰＲＩの揺らぎの上限がピーク・ツー・ピークで３０％であったとするとε＝  ０．１５  となる）。各ＰＲＩビンの中心ＰＲＩの値τｋ  とビンの幅  ｂｋ  は次式（１４）のように決める。
【００９８】
【数１４】
0003559205.tif 000015
【００９９】
【数１５】
0003559205.tif 000016
【０１００】
図１１には変形ＰＲＩ変換の処理のフローチャートを示す。入力データのパルスのＴＯＡを  ｔ１、  ｔ２、  …、  ｔＮとする。変形ＰＲＩ変換の処理の手順は以下のように行う。
【０１０１】
（１）初期設定を行う（ステップＳ５３１）。すなわち、ＰＲＩ変換  Ｄｋ（ｋ  ＝  １、  ２、  …  、  Ｋ）をゼロクリアする。また、パルスペアの右パルスのパルス番号を  ｍ  ＝  １  に設定する。
【０１０２】
（２）パルスペアの右パルスのパルス番号を更新する（ステップＳ５３２）。すなわち、ｍ  ＝  ｍ  ＋  １とする。
【０１０３】
（３）終了判定を行う（ステップＳ５３３）。すなわち、ｍ≦Ｎならばステップ５３４  へ。そうでなければ終了。
【０１０４】
（４）パルスペアの左パルスのパルス番号を  ｎ  ＝  ｍ  に設定する。（ステップＳ５３４）。
【０１０５】
（５）パルスペアの左パルスのパルス番号を更新する（ステップＳ５３５）。すなわち、ｎ  ＝  ｎ−１。
【０１０６】
（６）パルスペアの左パルスの終了判定を行う（ステップ  ５３６）。すなわち、ｎ  ≧  １  ならば、ステップＳ５３７へ。そうでなければステップＳ５３２へ。
【０１０７】
（７）パルスペアの時間間隔を求める（ステップ  ５３７）。すなわち、τ＝  ｔｍ−ｔｎを行う。
【０１０８】
（８）時間間隔の大きさ判定１を行う（ステップＳ５３８）。すなわち、τ≧τｍｉｎ  ならば、ステップＳ５３９へ。そうでなければステップＳ５３５へ。
【０１０９】
（９）時間間隔の大きさ判定２を行う（ステップＳ５３９）。すなわち、τ≦τｍａｘ  ならば、ステップＳ５４０へ。そうでなければステップＳ５３２へ。
【０１１０】
（１０）τ軸上でのＰＲＩビンの選定を行う（ステップＳ５４０）。すなわち、ＰＲＩビンの下限  ｋ１  と上限  ｋ２  を次式で計算する。
【０１１１】
【数１６】
0003559205.tif 000017
【０１１２】
ここで、［・］  は  Ｇａｕｓｓ  の整数、すなわち、（・）  を越えない最大の整数を表わす。
【０１１３】
（１１）ＰＲＩ変換の更新を行う（ステップＳ５４１）。すなわち、ｋ  ＝  ｋ１，…，  ｋ２  について以下の処理を行う。
【０１１４】
１）時間原点の初期設定を行う。すなわち、ＰＲＩビンｋが最初に使われる場合に限ってＯｋ  ＝  ｔｍ  とする。
【０１１５】
２）位相の予備計算と分解を行う。位相η０を次式で求め、それを整数部分νと小数部分ζに分解する。
【０１１６】
【数１７】
0003559205.tif 000018
【０１１７】
３）時間原点の移動を行う。すなわち、次の条件のいずれかが満たされるとき、時間原点を移動して、ｔｍ  をあらたな時間原点とする。
【０１１８】
【数１８】
0003559205.tif 000019
【０１１９】
４）位相の計算を行う。次式（１９）で新たな時間原点のもとでの位相を計算する。
【０１２０】
【数１９】
0003559205.tif 000020
【０１２１】
５）ＰＲＩ変換の価の更新を行う。次式（２０）でＰＲＩ変換の価を更新する。このあと、ステップ５３２へ。
【０１２２】
【数２０】
0003559205.tif 000021
【０１２３】
図１２乃至図１５にはこの実施の形態による発明の効果を示す。
【０１２４】
図１２は、本発明のパルス列分離装置への入力パルスデータの説明図である。図１２において、６１（出力１）と６２（出力２）はスタガーパルス列、６３（出力３）と６５（出力５）はジッターパルス列、６４（出力４）と６６（出力６）は固定ＰＲＩパルス列、６７は６１（出力１）から６６（出力６）までの６つのパルス列が重なりあった複合パルス列である。本発明のパルス列分離装置へは、複合パルス列６７が入力する。
【０１２５】
図１３は、図２の固定ＰＲＩパルス列検出・分離手段３００で固定ＰＲＩパルス列が検出され分離される様子を示した図である。図１３において、７１は固定ＰＲＩパルス列検出・分離手段３００への入力パルス列、７２はＰＲＩ変換と検出結果を示すグラフ、７７と７８は分離されたパルス列を表わす。７１の入力パルス列が図３のＰＲＩ変換演算部３３０へ入力されると、図５にフローチャートを示したＰＲＩ変換により、７３のようなＰＲＩスペクトルが得られる。しきい値演算部３５０で求めたしきい値７４とＰＲＩスペクトル７３との比較をＰＲＩ抽出部３６０で行うことにより、７５（ピーク１）及び７６（ピーク２）の２つのＰＲＩが検出される。検出されたＰＲＩをもとに、固定ＰＲＩパルス列分離演算部３７０により、もとのパルス列７１から抜き出したパルス列が７７と７８である。７５のピーク１は７７のパルス列に対応し、７６のピーク２は７８のパルス列に対応する。
【０１２６】
図１４は、図２のスタガーパルス列検出・分離手段４００でスタガーパルス列が検出され分離される様子を示した図である。図１４において、８１はスタガーパルス列検出・分離手段４００への入力パルス列、８２は３連パルスのヒストグラムと検出結果を示すグラフ、８７と８８は分離されたパルス列を表わす。８１の入力パルス列が、図６の３連パルスヒストグラム演算部４３０へ入力されると、図７にフローチャートを示した３連パルスのヒストグラムの処理により、８３のようなヒストグラムが得られる。このヒストグラム８３としきい値８４との比較により、８５（ピーク１）及び８６（ピーク２）の２つのＰＲＩが検出される。検出されたＰＲＩをもとに、もとのパルス列８１から抜き出したパルス列が８７と８８である。８５のピーク１は８７のパルス列に対応し、８６のピーク２は８８のパルス列に対応する。
【０１２７】
図１５は、図２のジッターパルス列検出・分離手段５００でジッターパルス列が検出され分離される様子を示した図である。図１５において、９１はジッターパルス列検出・分離手段５００への入力パルス列、９２は変形ＰＲＩ変換と検出結果を示すグラフ、９７と９８は分離されたパルス列を表わす。９１の入力パルス列が図９の変形ＰＲＩ変換演算部５３０へ入力されると、図１１にフローチャートを示した変形ＰＲＩ変換の処理により、９３のようなＰＲＩスペクトルが得られる。このＰＲＩスペクトル９３としきい値９４との比較により、９５（ピーク１）及び９６（ピーク２）の２つのＰＲＩが検出される。検出されたＰＲＩをもとにもとのパルス列９１から抜き出したパルス列が９７と９８である。９５のピーク１は９７のパルス列に対応し、９６のピーク２は９８のパルス列に対応する。
【０１２８】
以上のように、本発明のパルス列分離装置によれば、ＰＲＩ変換を用いる固定ＰＲＩパルス列検出・分離手段３００と、３連パルスのヒストグラムを用いるスタガーパルス列検出・分離手段４００と、変形ＰＲＩ変換を用いるジッターパルス列検出・分離手段を備えて、ＰＲＩの種類ごとに該当するパルス列を段階的に検出・分離するようにしたので、信号源の数が多い場合や、固定ＰＲＩ以外の複雑なＰＲＩのパターンを有する信号源まで含めて、信号源のＰＲＩを正確に推定し、精度良く、各パルス列を検出・分離できるという効果がある。
【０１２９】
【発明の効果】
この発明は、複数の信号源からのパルス列が重なり合った複合パルス列の各パルスの到着時間（ＴＯＡ）データに基づき、固定パルス繰り返し間隔（ＰＲＩ）を有する固定ＰＲＩパルス列を検出する固定ＰＲＩパルス列検出手段と、検出された固定ＰＲＩパルス列をもとの複合パルス列から分離するとともに、分離後の残りのパルス列を出力する固定ＰＲＩパルス列分離手段と、固定ＰＲＩパルス列分離手段から出力される残りのパルス列から、複合パルス列の各パルスの到着時間（ＴＯＡ）データに基づき、スタガーパルス列を検出するスタガーパルス列検出手段と、検出されたスタガーパルス列を残りのパルス列から分離するとともに、分離後の固定ＰＲＩパルス列及びスタガーパルス列を含まない第二の残りのパルス列を出力するスタガーパルス列分離手段と、スタガーパルス列分離手段から出力される第二の残りのパルス列から、複合パルス列の各パルスの到着時間（ＴＯＡ）データに基づき、ジッターパルス列を検出するジッターパルス列検出手段と、検出されたジッターパルス列を第二の残りのパルス列から分離するジッターパルス列分離手段と、を備えたパルス列分離装置であるので、ＰＲＩの種類ごとに該当するパルス列を段階的に検出・分離するようにしたので、信号源の数が多い場合や、固定ＰＲＩ以外の複雑なＰＲＩのパターンを有する信号源まで含めて、信号源のＰＲＩを正確に推定し、精度良く、各パルス列を検出・分離できるという効果がある。
【０１３０】
また、固定ＰＲＩパルス列検出手段が、ＰＲＩ変換を用いて検出を行うようにしたので、複合パルス列に含まれている部分パルス列のＰＲＩだけをピークとして抽出できるので、ＰＲＩを明確に識別することができるという効果がある。
【０１３１】
また、スタガーパルス列検出手段が、等間隔に並ぶ３つのパルスからなる３連パルスの集合に対してそのパルス繰り返し間隔のヒストグラムを求める３連パルスヒストグラム演算部と、得られたヒストグラムから複数のピークを検出し、各ピークがパルス列の基本周期であるか否かを判定する判定部と、を備えているので、得られるヒストグラムでは、スタガーパルス列の基本周期とその整数倍のところにするどいピークがたつので、スタガーパルスの基本周期を明確に識別することができるという効果がある。
【０１３２】
また、ジッターパルス列検出手段が、オーバーラップしたＰＲＩビンを持ち、時間原点の移動を行う、変形ＰＲＩ変換を用いるようにしたので、ジッターパルス列の平均ＰＲＩのところにするどいピークがたつので、ジッターパルス列の平均ＰＲＩを明確に識別することができるという効果がある。
【０１３３】
また、この発明は、複数の信号源からのパルス列が重なり合った複合パルス列の各パルスの到着時間（ＴＯＡ）データに基づき、固定パルス繰り返し間隔（ＰＲＩ）を有する固定ＰＲＩパルス列を検出する固定ＰＲＩパルス列検出工程と、検出された固定ＰＲＩパルス列をもとの複合パルス列から分離するとともに、分離後の残りのパルス列を出力する固定ＰＲＩパルス列分離工程と、固定ＰＲＩパルス列分離手段から出力される残りのパルス列から、複合パルス列の各パルスの到着時間（ＴＯＡ）データに基づき、スタガーパルス列を検出するスタガーパルス列検出工程と、検出されたスタガーパルス列を残りのパルス列から分離するとともに、分離後の固定ＰＲＩパルス列及びスタガーパルス列を含まない第二の残りのパルス列を出力するスタガーパルス列分離工程と、スタガーパルス列分離手段から出力される第二の残りのパルス列から、複合パルス列の各パルスの到着時間（ＴＯＡ）データに基づき、ジッターパルス列を検出するジッターパルス列検出工程と、検出されたジッターパルス列を第二の残りのパルス列から分離するジッターパルス列分離工程と、を備えたパルス列分離方法であるので、ＰＲＩの種類ごとに該当するパルス列を段階的に検出・分離するようにしたので、信号源の数が多い場合や、固定ＰＲＩ以外の複雑なＰＲＩのパターンを有する信号源まで含めて、信号源のＰＲＩを正確に推定し、精度良く、各パルス列を検出・分離できるという効果がある。
【０１３４】
また、固定ＰＲＩパルス列検出工程において、ＰＲＩ変換を用いて検出を行うようにしたので、複合パルス列に含まれている部分パルス列のＰＲＩだけをピークとして抽出できるので、ＰＲＩを明確に識別することができるという効果がある。
【０１３５】
また、スタガーパルス列検出工程が、等間隔に並ぶ３つのパルスからなる３連パルスの集合に対してそのパルス繰り返し間隔のヒストグラムを求める３連パルスヒストグラム演算ステップと、得られたヒストグラムから複数のピークを検出し、各ピークがパルス列の基本周期であるか否かを判定する判定ステップと、を備えているので、得られるヒストグラムでは、スタガーパルス列の基本周期とその整数倍のところにするどいピークがたつので、スタガーパルスの基本周期を明確に識別することができるという効果がある。
【０１３６】
また、ジッターパルス列検出工程において、オーバーラップしたＰＲＩビンを持ち、時間原点の移動を行う、変形ＰＲＩ変換を用いるようにしたので、ジッターパルス列の平均ＰＲＩのところにするどいピークがたつので、ジッターパルス列の平均ＰＲＩを明確に識別することができるという効果がある。
【図面の簡単な説明】
【図１】ＰＲＩの種類を説明した説明図である。。
【図２】本発明の実施の形態１にかかわるパルス列分離装置の構成を示すブロック図である。
【図３】この発明の実施の形態１にかかわるパルス列分離装置における固定ＰＲＩパルス列検出・分離手段の内部構成を示すブロック図である。
【図４】この発明の実施の形態１にかかわるパルス列分離装置における固定ＰＲＩパルス列検出・分離手段の中で用いるＰＲＩ変換に関するＰＲＩビンの説明図である。
【図５】この発明の実施の形態１にかかわるパルス列分離装置における固定ＰＲＩパルス列検出・分離手段の中で用いるＰＲＩ変換のフローチャートである。
【図６】この発明の実施の形態１にかかわるパルス列分離装置におけるスタガーパルス列検出・分離手段の内部構成を示すブロック図である。
【図７】この発明の実施の形態１にかかわるパルス列分離装置におけるスタガーパルス列検出・分離手段の中で用いる３連パルスのヒストグラムのフローチャートである。
【図８】この発明の実施の形態１にかかわるパルス分離装置におけるスタガーパルス列検出・分離手段の中で用いる３連パルスのヒストグラムからの基本周期検出のフローチャートである。
【図９】この発明の実施の形態１にかかわるパルス列分離装置におけるジッターパルス列検出・分離手段の内部構成を示すブロック図である。
【図１０】この発明の実施の形態１にかかわるパルス列分離装置におけるジッターパルス列検出・分離手段の中で用いる変形ＰＲＩ変換に関するＰＲＩビンの説明図である。
【図１１】この発明の実施の形態１にかかわるパルス列分離装置におけるジッターパルス列検出・分離手段の中で用いる変形ＰＲＩ変換のフローチャートである。
【図１２】この発明の実施の形態１にかかわるパルス列分離装置への入力パルス列の説明図である。
【図１３】この発明の実施の形態１にかかわるパルス列分離装置における固定ＰＲＩパルス列検出・分離手段での処理結果を示す図である。
【図１４】この発明の実施の形態１にかかわるパルス列分離装置におけるスタガーパルス列検出・分離手段での処理結果を示す図である。
【図１５】この発明の実施の形態１にかかわるパルス列分離装置におけるジッターパルス列検出・分離手段での処理結果を示す図である。
【図１６】従来のパルス列分離装置の動作を示すフローチャートである。
【符号の説明】
１１  固定ＰＲＩパルス列、１２  スタガーパルス列、１３  ジッターパルス列、６１，６２  スタガーパルス列、６３，６５  ジッターパルス列、６４，６６  固定ＰＲＩパルス列、６７  複合パルス列、７１  固定ＰＲＩパルス列検出・分離手段への入力パルス列、７２  ＰＲＩ変換と検出結果、７３  ＰＲＩスペクトル、７４  ＰＲＩスペクトルのしきい値、７５，７６  検出された固定ＰＲＩパルス列のＰＲＩ、７７，７８  分離された固定ＰＲＩパルス列、８１  スタガーパルス列検出・分離手段への入力パルス列、８２  ３連パルスのヒストグラムと検出結果、８３  ３連パルスのヒストグラム、８４  ヒストグラムのしきい値、８５，８６  検出されたスタガーパルス列のＰＲＩ、８７，８８  分離されたスタガーパルス列、９１  ジッターパルス列検出・分離手段への入力パルス列、９２  変形ＰＲＩ変換と検出結果、９３  ＰＲＩスペクトル、９４  ＰＲＩスペクトルのしきい値、９５，９６  検出されたジッターパルス列のＰＲＩ、９７，９８  分離されたジッターパルス列、１００  広帯域受信機、２００  パルスＴＯＡデータの入力端、３００  固定ＰＲＩパルス列検出・分離手段、４００  スタガーパルス列検出・分離手段、５００  ジッターパルス列検出・分離手段、６００  出力端、３１０  固定ＰＲＩパルス列検出・分離手段への入力端、３３０  ＰＲＩ変換演算部、３５０  しきい値演算部、３６０  固定ＰＲＩパルス列のＰＲＩ抽出部、３７０  固定ＰＲＩパルス列分離演算部、３８０  残存パルス列判定部、３９０  固定ＰＲＩパルス列検出・分離手段からの出力端、４１０  スタガーパルス列検出・分離手段への入力端、４３０  ３連パルスヒストグラム演算部、４５０  しきい値演算部、４６０  スタガーパルス列のＰＲＩ抽出部、４７０  スタガーパルス列分離演算部、４８０  残存パルス列判定部、４９０  スタガーパルス列検出・分離手段からの出力端、５１０  ジッターパルス列検出・分離手段への入力端、５３０  変形ＰＲＩ変換演算部、５５０  しきい値演算部、５６０  ジッターパルス列のＰＲＩ抽出部、５７０  ジッターパルス列分離演算部、５８０  残存パルス列判定部、５９０  ジッターパルス列検出・分離手段からの出力端。
【図１】
0003559205.tif 000022
【図２】
0003559205.tif 000023
【図３】
0003559205.tif 000024
【図４】
0003559205.tif 000025
【図５】
0003559205.tif 000026
【図６】
0003559205.tif 000027
【図７】
0003559205.tif 000028
【図８】
0003559205.tif 000029
【図９】
0003559205.tif 000030
【図１０】
0003559205.tif 000031
【図１１】
0003559205.tif 000032
【図１２】
0003559205.tif 000033
【図１３】
0003559205.tif 000034
【図１４】
0003559205.tif 000035
【図１５】
0003559205.tif 000036
【図１６】
0003559205.tif 000037