(57)【要約】  パルストレイン内の候補となる瞬間を選択する。かかる瞬間の各々について、その瞬間の前後両方において均一なインターバルでパルスが発生するかを判定する。発生していれば、関連するパルスを、特定のグループに属するものとして分類する。候補となる瞬間に対するパルスの到着時刻から繰り返しインターバルを算出し、重み付け角度平均を取ることにより、パルス間のインターバルが正確に計算される。
2010525315.tif 000002
【特許請求の範囲】
【請求項１】
パルストレインを含む受信信号を解析する方法であって、前記パルスはそれぞれの到着時刻において発生し、
前記方法は：
（ａ）パルスのグループを形成するステップであって、前記グループ内のパルスが、所定の最小値および最大値により定義される範囲内で実質的に均一なパルス繰り返しインターバルを有するように、形成するステップと；
（ｂ）平均グループパルス繰り返しインターバルを計算するステップであって：
（ｉ）前記グループにおける複数のパルスの各々について、それぞれの概算の繰り返しインターバルを決定するステップと；
（ｉｉ）前記所定の最小値および最大値の間の前記範囲が、π以下の長さの弧により表されるように線形マッピング処理を行うことにより、各前記概算の繰り返しインターバルからベクトルを決定するステップと；
（ｉｉｉ）前記決定されたベクトルのベクトル和の角度を決定するステップと；
（ｉｖ）前記マッピング処理の逆を行うことで、前記ベクトル和の前記角度から前記平均グループパルス繰り返しインターバルを導出するステップと
により計算するステップと；を含む、方法。
【請求項２】
各概算の繰り返しインターバルは、そのパルスの到着時刻と共通の選択された瞬間との差から決定される、請求項１に記載の方法。
【請求項３】
前記方法は、前記グループのパルスを選択するステップを含み、
前記共通の選択された瞬間は、前記選択されたパルスの到着時刻であり、
それにより、各概算の繰り返しインターバルは、それぞれのパルスの到着時刻と前記選択されたパルスの到着時刻との差から決定される、請求項２に記載の方法。
【請求項４】
前記共通の選択された瞬間は、パルスの到着時刻と一致しない、請求項２に記載の方法。
【請求項５】
前記共通の選択された瞬間は、前記グループにおける前記パルスの到着時刻の前記範囲の中央であるか、またはその付近である、請求項２〜４のいずれか一項に記載の方法。
【請求項６】
前記ベクトル和の前記角度は、前記パルスのそれぞれの到着時刻と前記選択された瞬間との間のインターバルの増加に伴ってそのパルスの影響が増加する重み付け処理を用いて計算される、請求項２〜５のいずれか一項に記載の方法。
【請求項７】
前記決定されたベクトルから、到着時刻間の差のばらつきを示すばらつき値を計算するステップを含む、請求項１〜６のいずれか一項に記載の方法。
【請求項８】
前記グループに分類された前記パルスを除外することにより、前記トレインから１組のパルスを導出するステップと、
その後、前記組におけるパルスを別のグループに割り当てるステップと
を含み、
前記除外ステップは、前記ばらつき値を条件として行われる、請求項７に記載の方法。
【請求項９】
ステップ（ａ）は：
（Ｉ）候補となる瞬間を選択し、少なくとも、前記候補となる瞬間の前の第１の所定の時間範囲を計算し、前記候補となる瞬間の後の第２の所定の時間範囲を計算するステップと；
（ＩＩ）到着時刻が前記所定の時間範囲内であるパルスについてカウント値を導出するステップであって、
前記カウント値は、それぞれの時間インターバルに各々関連付けられ、
前記カウント値は、それぞれのパルスの到着時刻と、候補となる瞬間との間の絶対差が、実質的に前記時間インターバルの整数倍に等しいようなパルスの数に各々依存する
カウント値を導出するステップと；
（ＩＩＩ）前記カウント値の１つが所定の閾値を超える場合、前記１つのカウント値に関連するパルスを前記パルスグループに割り当てるステップと；を含む、請求項１〜８のいずれか一項に記載の方法。
【請求項１０】
前記選択された瞬間は、前記候補となる瞬間である、請求項２に直接的または間接的に従属する請求項９に記載の方法。
【請求項１１】
パルストレインを含む受信信号を解析する方法であって、前記パルスはそれぞれの到着時刻において発生し、
前記方法は：
（Ｉ）候補となる瞬間を選択し、少なくとも、前記候補となる瞬間の前の第１の所定の時間範囲を計算し、前記候補となる瞬間の後の第２の所定の時間範囲を計算するステップと；
（ＩＩ）到着時刻が前記所定の時間範囲内であるパルスについてカウント値を導出するステップであって、
前記カウント値は、それぞれの時間インターバルに各々関連付けられ、
前記カウント値は、それぞれのパルスの到着時刻と、候補となる瞬間との間の絶対差が、実質的に前記時間インターバルの整数倍に等しいようなパルスの数に各々依存する
カウント値を導出するステップと；
（ＩＩＩ）前記カウント値の１つが所定の閾値を超える場合、前記１つのカウント値に関連付けられたパルスをパルスグループに割り当てるステップと；を含む方法。
【請求項１２】
前記グループに含まれる候補のパルスを選択するステップを含み、前記候補となる瞬間は、前記選択された候補のパルスの到着時刻である、請求項１１に記載の方法。
【請求項１３】
前記候補となる瞬間は、パルスの到着時刻と一致しない、請求項１１に記載の方法。
【請求項１４】
前記方法は、第１および第２の振幅閾値を前記受信信号に適用するステップを含み、
パルスを第１のパルスの種類と第２のパルスの種類とに分類するために、前記第２の閾値は前記第１の閾値よりも大きくなっており、
前記第２の種類のパルスは、前記第１の種類のパルスよりも大きい振幅を有し、
前記候補のパルスは、前記第２のタイプのパルスの組から選択される、請求項１２に記載の方法。
【請求項１５】
前記グループに分類された前記パルスを検討対象から除外した後にステップ（Ｉ）〜（ＩＩＩ）を繰り返すことにより、残りのパルスのいくつかを別のグループに割り当てるステップを含む、請求項９〜１４のいずれか一項に記載の方法。
【請求項１６】
前記第１および第２の範囲は、
前記候補となる瞬間から第１の有限所定量だけ異なるそれぞれの最小範囲値から、
前記候補となる瞬間から第２の有限所定量だけ異なるそれぞれの最大範囲値まで
及ぶ、請求項９〜１５のいずれか一項に記載の方法。
【請求項１７】
前記候補となる瞬間から前記第１の有限所定量の２倍だけ異なるそれぞれの最小範囲値から、
前記候補となる瞬間から前記第２の有限所定量の２倍だけ異なるそれぞれの最大範囲値まで
及ぶ、第３および第４の範囲を含む、請求項１６に記載の方法。
【請求項１８】
ステップ（Ｉ）〜（ＩＩＩ）は、第１および第２の所定の時間範囲を前記候補となる瞬間から漸進させながら連続的に繰り返される、請求項９〜１７のいずれか一項に記載の方法。
【請求項１９】
各カウント値は、パルスについての前記絶対差が前記時間インターバルを含む所定の範囲を定義する窓内に存在する値の整数倍に等しい、そのようなパルスの数を表し、
前記方法は、前記窓をシフトして他の時間インターバルについてのカウント値を得るステップを含む、請求項９〜１８のいずれか一項に記載の方法。
【請求項２０】
同じグループに分類されたパルスがバーストにおいて発生することを判定するステップを含み、
前記方法は、前記バーストの持続時間およびそれらの間のインターバルを示す情報を導出するステップを含む、請求項１〜１９のいずれか一項に記載の方法。
【請求項２１】
分類されたパルスが出現しないことが予想されるインターバルを計算し、そのインターバルを表す信号を生成するステップを含む、請求項１〜２０のいずれか一項に記載の方法。
【請求項２２】
請求項１〜２１のいずれか一項に記載の方法に従って動作するように編成された装置。
【請求項２３】
請求項２１に記載の方法に従って動作するように編成された請求項２２に記載の装置を含む送信システムであって、
電磁信号を送信するための送信手段と、
前記周期を表す信号に応答して前記送信手段の動作を可能にするための制御手段と
をさらに含む、送信システム。
【発明の詳細な説明】
【技術分野】
【０００１】
［発明の分野］
本発明は、特定の周期的または規則的なパターンを表す電磁パルスなどのイベントを解析および分類するための方法および装置に関し、特に、様々な無線測位サービスと同じ周波数帯域を共有する無線通信またはセンサネットワークにおいて発生し得る無線周波数パルスのバーストがない時間インターバルの予測に適用可能であるが、それに限定されるものではない。
【背景技術】
【０００２】
［従来技術の説明］
多くの事例において、注目対象の現象を表す観測信号には、未知数の異なるソースにより産生された複数のインターリーブされたパルストレイン（またはバースト）が含まれることがある。各観測パルストレインを、そのトレインを発する対応するソースに割り当てる必要があることがよくある。一般に、かかる関連付け手順は、異なるソースが、異なるパラメータ（パルスの振幅、パルスの形状および持続時間、ならびにパルス間のインターバル等）を有するパルストレインを産生する、という仮定に基づいている。電磁環境において、パルスの分類および関連付けは、観測された到着方向、キャリアの周波数および偏波、パルス内変調の種類などを含む他の特性に基づくこともある。
【０００３】
「パルスデインターリービング」と一般に呼ばれる上記のパルス関連付けには、地震探査、コンピュータ通信、ニューラルパルスネットワーク、電波天文学（例えば、パルサの検索）などを含む様々な分野における多くの周知の用途がある。最近では、初歩的な形態のパルスデインターリービングを、新興分野であるインテリジェント（または「コグニティブ」）無線通信システムに適用する試みがなされてきた。
【０００４】
インテリジェント無線の概念は、適切な規制当局によりその全体が様々な商用および軍用のユーザに割り当てられているスペクトル帯域の（空間、時間、または周波数における）遊休部分を、適応的に利用することに基づいている。スペクトルの適応的な利用により、様々な広帯域通信サービスおよびマルチメディアサービス、ならびに監視アプリケーションをサポートする周波数スペクトルへの高まる需要が、少なくとも部分的に満たされることが期待される。
【０００５】
スペクトルの利用を最適化するためにインテリジェントシステムが用いる一般的なアプローチは：
−周囲の電磁環境を感知して、スペクトル機会を識別するステップと；
−遊休スペクトルセグメントを好適に用いて、意図する動作（感知または情報送信等）を行うステップと；
−主たるまたは優先度が高いユーザが出現する前に、スペクトルセグメントを空けるステップと；を含む。
【０００６】
用いるアプローチに関わらず、インテリジェントシステムは、周囲の電磁環境を認識することと、自動的に反応する能力との両方を発揮する必要がある。
【０００７】
５ＧＨｚ周波数帯域（５ＧＨｚと６ＧＨｚとの間の周波数）において開発中のインテリジェント無線システムの一例には、各種の国際標準（電気電子学会（ＩＥＥＥ）および欧州電気通信標準化機構（ＥＴＳＩ）により策定されたもの等）に従って動作するように編成されたワイヤレスローカルエリアネットワーク（ＷＬＡＮ）の類がある。ＷＬＡＮがこのセグメントのスペクトルを利用できるようにした結果、主たるシステム、とりわけ５ＧＨｚ帯域において動作するライセンスを受けたレーダと、このスペクトルを共有する必要が生じている。
【０００８】
レーダシステムとＷＬＡＮワイヤレス装置との両方が同じ周波数帯域において動作することを許容されているため、ＷＬＡＮ業界は、相互干渉を排除するために周波数調整を行う必要を回避すべく、動的周波数選択（ＤＦＳ）のメカニズムを策定した。従って、特定のチャネル上でレーダ信号の存在が検出されると、ＷＬＡＮ装置は、自動的に他のチャネルに切り換えることで、レーダ信号との干渉を回避しなければならない。
【０００９】
国土防衛の支援におけるレーダの役割の拡張の結果、それらを配備する場所は、過去のものとは異なっている。よって、レーダおよび他のセンサネットワークは、同じ地理的領域において動作することが必要であり、同じ周波数帯域を共有することが有益である。ネットワーク構成において動作するレーダおよび他のセンサにより、重要インフラ（電力網および発電所、ガスおよびオイルのパイプライン、ならびに給水システムを含む）の連続的監視の改良が提供される。
【００１０】
他の重要な用途には、音響センサや他のセンサを用いる浮遊ブイのネットワークにより、注目対象の高速ボートおよび他の水上艦を検出および位置特定することが可能な沿岸警備隊または沿岸監視の用途がある。それらのセンサは、無線通信リンクを介して情報を交換する必要があり、固定や機上の監視レーダと同じ周波数帯域を共有することが有益である。
【００１１】
既存の規制基準が提供するレーダの検出および回避のためのガイドラインは単純なものであるが、複合電磁信号を傍受および解析してレーダ信号のパラメータを決定するための具体的な方法は与えられていない。
【００１２】
通信ネットワークにおけるレーダパルスの検出および回避に関する従来技術の例は、下記の特許および特許出願公開において提示されている：
１．欧州特許出願公開第１５０５７７２号明細書、Ｋｒｕｙｓ他、２００５年２月９日
２．米国特許第６８３１５８９号明細書、Ｓｈｅａｒｅｒ  ＩＩＩ、２００４年１２月１４日
３．米国特許出願公開第２００４／０１５６３３６号明細書、ＭｃＦａｒｌａｎｄ他、２００４年８月１２日
４．米国特許出願公開第２００４／０１５１１３７号明細書、ＭｃＦａｒｌａｎｄ他、２００４年８月５日
５．米国特許第６６９７０１３号明細書、ＭｃＦａｒｌａｎｄ他、２００４年２月２４日（国際公開第０３／０５０５６０号パンフレット、２００３年６月１９日）
６．米国特許出願公開第２００４／００３３７８９号明細書、Ｔｓｉｅｎ、２００４年２月１９日
７．米国特許出願公開第２００３／０２１４４３０号明細書、Ｈｕｓｔｅｄ他、２００３年１１月２０日
８．米国特許出願公開第２００３／０２０６１３０号明細書、Ｈｕｓｔｅｄ他、２００３年１１月６日
９．米国特許出願公開第２００２／０１５５８１１号明細書、Ｐｒｉｓｍａｎｔａｓ他、２００２年１０月２４日
１０．国際公開第０２／０８２８４４Ａ３号パンフレット、Ｚｉｍｍｅｒｍａｎｎ他、２００２年１０月１７日（欧州特許出願公開第１２４８４７７号明細書、２００２年１０月９日）
上記のリストは網羅的であることを意図したものではないが、引用文献において開示されている手法および実施形態は、ＷＬＡＮ装置およびレーダシステムの共存に関する最新の従来技術を表すものと思われる。
【発明の概要】
【発明が解決しようとする課題】
【００１３】
従来技術による手法のいずれにも、同じ周波数帯域において動作する複数のレーダにより発信されたインターリーブされたレーダパルスバーストの特性を決定する能力はない。特に、観測された複合パルスストリームにおいて、いくつかのレーダパルスが欠落しており、いくつかの「誤った」パルス（例えば、雑音により生成されたもの）が出現している場合において、そのような能力はない。そのため、インターリーブされたパルスバーストのパラメータ（各バースト内でのパルス間のインターバル、バースト持続時間、バースト繰り返し周波数、および異なるレーダシステムから発せられたバースト間の時間シフト等）を決定するためのロバストな手法を提供することが望ましい。
【００１４】
ほとんどの場合、監視レーダは、連続的なパルス送信ではなく、間欠的にパルスのバーストを（任意の特定の方向に）発信する。よって、送信に用いられる周波数帯域は、バーストとバーストの間では遊休状態のままである。このため、間欠的送信を行う複数の監視レーダが、感知または情報送信に広帯域信号または超広帯域（ＵＷＢ）信号を用いる他のＲＦ装置と同じ周波数帯域において動作する際に、改善されたスペクトル利用を提供する、適応的方法を提供することが望ましい。
【課題を解決するための手段】
【００１５】
［発明の要約］
本発明の態様は、付帯の請求項において述べる。
【００１６】
本発明のさらなる独立態様によれば、パルスの、１つ以上の別々のグループを発見する試行が行われる。各グループのパルスは、それらが、恐らくは共通のソースからの、まとまった組(coherent set)を形成することを示すパルス間インターバルを有する。各試行には、パルスの１つを候補のパルスとして選択するステップが含まれる。選択されたパルスに対して均一なパルス繰り返しインターバルで到着するパルスを位置合わせする（すなわち、それらのパルスを表す値を合成する）ように、他のパルスの到着時刻を処理する。（選択された候補のパルスの到着時刻を参照して他のパルスの到着時刻を処理するため、本明細書では候補のパルスを「ピボット」パルスとも呼ぶ。）十分な数のパルスが位置合わせされれば、これらのパルスおよび選択された候補のパルスがまとまったグループを形成するものと仮定する。
【００１７】
処理は、好ましくは、選択された候補のパルスの前に到着したパルスおよび後に到着したパルスの両方を考慮する。これにより、統計的により正確な結果を得ることが可能になる。
【００１８】
好ましくは、受信信号におけるパルスを、低振幅パルスと高振幅パルスとに分割する。後者は、雑音の影響であるよりは有効なパルスである可能性が高い。高振幅パルスのみを候補のパルスとして用いることにより、処理が著しく速められる。
【００１９】
本発明のさらなる独立態様（上記の独立態様と組み合わせて用いるのが好ましいが、別々に、当該技術において公知のものを含む異なるパルスグループ化編成とともに用いることも可能である）によれば、パルスのグループのパルス繰り返しインターバルを、パルス到着時刻を測定することにより決定される概算のパルス繰り返しインターバルの角度平均を決定することにより計算する。特に、測定されたパルス繰り返しインターバルを、π以下の角度範囲にわたって線形に分布する角度付きベクトルに変換する。範囲の端点は、所定の最小および最大の予想されるパルス繰り返しインターバルに対応する。平均化は、ベクトル和の角度を決定することにより果たされ、この角度は、実際のパルス繰り返しインターバルの正確な推定値を表す。
【００２０】
好ましくは、測定されたパルス間のインターバルは、各パルスの到着時刻と中間のパルス（好ましくは中央のパルス）の到着時刻との間のインターバルである。適切な場合は、この測定されたインターバルを到着時刻間のインターバルの数で除算して、概算のパルス繰り返しインターバルを得る。本発明の２つの上記態様を組み合わせる編成において、中間のパルスは、上述の候補のパルスであってもよい。
【００２１】
好ましくは、中央のパルスからより離れたパルスの到着時刻の影響を相対的に大きくして、重み付け角度平均を導出する。従って、測定されたインターバルが大きい（丸め誤差の影響はそれに比例して小さい）ほど、影響が大きい。
【００２２】
以下、添付の図面を参照して、本発明を実施する編成を例示により説明する。
【図面の簡単な説明】
【００２３】
【図１】センサ群Ｓ１、Ｓ２、．．．、ＳＫ、複数の監視レーダＲ１、Ｒ２、．．．、ＲＬ、情報融合センタＩＦＣ、およびレーダ信号解析器ＲＳＡを備える安全システムの例である。
【図２】本発明に従って動作するように編成されたレーダ信号解析器ＲＳＡのブロック図である。
【図３】図３ａは、３シーケンスのパルスバーストの例を示す。図３ｂは、部分的に重なったパルスバーストを含む複合信号のセグメントを示す。図３ｃは、２値の送信許可信号を示す。
【図４】本発明による２つの判定閾値と比較した時間変化信号包絡線の例を示す。
【図５】本発明に従って候補のパルス繰り返しインターバルを予備決定するために用いる方法を説明する図である。
【図６】本発明に従って単位半円の円周上にイベント５連体(event quintuple)をマッピングする例である。
【図７】最も不調和なイベント５連体を示す。
【図８】図２の解析器の繰り返しインターバル推定器のブロック図である。
【図９】ヒストグラムを産生するとともにバースト当たりパルス数を決定するため、複数の時間フレーム内で発生するイベントに対して本発明に従って行われる動作を概略的に示す。
【図１０】本発明の変形例による候補のパルス繰り返しインターバルを予備決定するために用いる方法を説明する図である。
【発明を実施するための形態】
【００２４】
図１は、注目対象の領域にわたり分散させた好適なセンサ群(a constellation of suitable sensors)Ｓ１、Ｓ２、．．．、ＳＫと、情報融合センタＩＦＣとを備える安全監視システムの例を示す。センサ群には、パッシブセンサもしくはアクティブセンサ、またはそれらの組み合わせを利用することができる。さらに、センサ群を形成するセンサは、異種のもの（例えば、音響センサ、赤外センサ、または電磁的センサ）であってもよい。
【００２５】
広帯域レーダ信号解析器ＲＳＡも、センサ群により利用される。
【００２６】
各センサおよび情報融合センタの間での高データレート情報交換に利用可能な広帯域通信チャネルが設けられるものと仮定する。また、センサ群の近傍において、広帯域通信チャネルを形成する周波数帯域と一致する周波数帯域（の少なくとも一部）を送信に利用する監視レーダＲ１、Ｒ２、．．．、ＲＬが動作しているものと仮定する。
【００２７】
また、監視レーダＲ１、Ｒ２、．．．、ＲＬの各々は、３６０度の探知範囲を提供する回転アンテナを用いるものと仮定する。その結果、センサ群は、レーダのアンテナがセンサ群の位置を直接指しているときにのみ、レーダパルス送信の影響を受ける。その位置では、各動作レーダが、センサ群に向けて送られる周期的パルスバーストにより自らの存在を明らかにする。
【００２８】
各パルスバーストの持続時間は、アンテナの回転速度およびアンテナビームの角度幅により決定される。各バーストにおいて生じるレーダパルスの数は、バースト持続時間とレーダにより用いられるパルス繰り返し周波数とに依存する。異なる特性を有するいくつかのレーダが同時に動作しているときは、センサ群は、インターリーブされたパルスバーストと、部分的または完全に重なったバースト内で発生するインターリーブされたパルスとを含む複合パルス信号を受信する。
【００２９】
本発明に従って構成された広帯域レーダ信号解析器ＲＳＡは、以下の機能および動作を行う：
−所定の検出閾値を超えるパルスの検出。各検出パルスは、「有効な」レーダパルスであるか、またはいずれの動作レーダにも関係しない雑音もしくは干渉により産生されたものである。
−検出パルスの到着時刻（ＴＯＡ）の登録された値(registered value)により決定されるイベントの複合ストリームを解析して、周期的パルスバーストの数、各バーストの周期、およびその持続時間を決定する。
−レーダパルスのバーストがない時間インターバルの予測。予測されたインターバルを用いて、センサ群により利用されるセンサおよび情報融合センタＩＦＣの間で高データレート情報交換が行われる。
−送信許可時間インターバルを示す好適な信号の生成。かかる信号を情報融合センタＩＦＣに送り、センサ群内での情報送信を調整することができる。
【００３０】
図２は、本発明に従って動作するように編成されたレーダ信号解析器（ＲＳＡ）のブロック図である。解析器は、広帯域受信機（ＷＢＲ）に接続された全方向性アンテナ（ＯＤＡ）を備え、広帯域受信機には包絡線検出器（ＥＤＲ）が続いている。解析器は、また、比較器（ＣＭＲ）、イベント記述子生成器（ＥＤＧ）、イベント記述子フィルタ（ＥＤＦ）、イベント記述子バッファ（ＥＤＢ）、繰り返しインターバル推定器（ＲＩＥ）、および送信許可信号ＣＴを産生する演算／制御ユニット（ＡＣＵ）を備え、これらのすべては、ハードウェアおよび／またはソフトウェアで実装することが可能である。
【００３１】
レーダ信号解析器ＲＳＡは、その動作中、広帯域雑音を傍受しているか、または、広帯域受信機（ＷＢＲ）によりカバーされる広周波数帯域内で動作しているかもしれない未知数のレーダ発信器により産生された複数のパルスバーストを傍受しているものと仮定する。各間欠的レーダ送信のパラメータ（パルス持続時間、パルスピークパワー、バースト持続時間（ＢＤ）、バースト繰り返しインターバル（ＢＲＩ）、およびバースト内でのパルス繰り返しインターバル（ＰＲＩ）等）は、先験的には未知であるものと仮定する。
【００３２】
図３ａは、３つのレーダ発信器によりそれぞれ産生されたパルスバーストの３つのシーケンスの例を示す。各シーケンスは、自らのバースト繰り返しインターバル（ＢＲＩ）、バースト内パルス繰り返しインターバル（ＰＲＩ）、およびバースト持続時間（ＢＤ）により特徴付けられる。レーダ信号解析器ＲＳＡは、それらのパルスバーストの加法的合成である複合信号を傍受している。
【００３３】
図３ｂは、部分的に重なったパルスバーストを含む複合信号のセグメントを示し、図３ｃは、演算／制御ユニットＡＣＵにより産生された所望の２値の送信許可信号ＣＴの例を示す。
【００３４】
包絡線検出器ＥＤＲにより供給される複合信号包絡線ＥＶのレベルは、比較器ＣＭＲにおいて、２つの好適に選定された所定の閾値、すなわち低閾値ＤＤおよび高閾値ＰＰと比較される。ある時点で、包絡線ＥＶが低閾値ＤＤを超えた場合は、パルス検出が宣言される。包絡線ＥＶが高閾値ＰＰも超えた場合は、観測信号が雑音スパイクでなく「有効な」レーダパルスを表す可能性がさらに高い。以下で説明するように、より高い確度で検出されたパルスは、パルスバーストおよびそれぞれのパラメータを効率的に決定するために重要な役割を果たす。
【００３５】
図４は、レーダ信号解析器ＲＳＡが傍受した信号の時間変化包絡線ＥＶの例を示す。本例において、包絡線は、３つの時点Ｔk1、Ｔk2、およびＴk3において低閾値ＤＤと交差している。これらの時点の各々は、イベントの到着時刻ＴＯＡ（これは閾値ＤＤと交差した時点に等しい）を決定している。図示のように、包絡線ＥＶは、Ｔk1において低閾値ＤＤと交差した直後に、高閾値ＰＰとも交差している。かかる場合を表すすべてのイベントは、他のイベント（包絡線が低閾値ＤＤのみと交差し、所定の時間インターバルの間、高閾値ＰＰを下回り続ける状況に対応する）から区別される。
【００３６】
以下では、低閾値ＤＤおよび高閾値ＰＰを、それぞれ「イベント検出」閾値および「ピボット判定」閾値と呼ぶ。よって、各「ピボットイベント」は、包絡線が、Ｔkで検出閾値ＤＤと交差し、時間インターバル［Ｔk，（Ｔk＋δ）］内でピボット判定閾値ＰＰと交差する状況を表す。ここで、δは所定の値（例えば、δ＝０．２μｓ）を有する。
【００３７】
イベント検出閾値ＤＤの値は、雑音のみに起因する正交差(upcrossings)の回数を許容可能なレベルに保つように選択される。例えば、中心周波数ｆ0および矩形周波数特性｜ｆ−ｆ0｜＜ΔＦ／２の帯域通過フィルタにガウス雑音を通過させる場合、フィルタの出力で得られる雑音信号の包絡線には、平均で、持続時間Ｔの時間インターバル内でＮ+回のレベルλの正交差が産生される。ここで、
【数１】
2010525315.tif 000003
であり、σΕは包絡線の平方二乗平均である。
【００３８】
例えば、ΔＦ＝１０ＭＨｚおよびλ＝３．５σΕであれば、平均で、持続時間１００ｍｓのインターバルにおいて１７回の正交差が登録される。しかし、交差レベルλをλ＝３σΕに低下させると、１００ｍｓのインターバルにおいて観測される正交差の平均回数は３８０に増加する。
【００３９】
ピボット判定閾値ＰＰの値は、ＰＰ＝ｇＤＤとして表現することが可能で、ここで定数ｇ＞１を経験的に、または理論的解析から決定することが可能である。
【００４０】
図２の比較器ＣＭＲの２つの２値出力ＸＬおよびＸＨは、それぞれイベント検出閾値ＤＤおよびピボット判定閾値ＰＰの交差を表す。イベント記述子生成器ＥＤＧは、比較器ＣＭＲの出力（ＸＬ、ＸＨ）を、演算／制御ユニットＡＣＵにより供給される「時間タグ」ＴＴと組み合わせる。演算／制御ユニットＡＣＵは、この目的のため、好適に選定された周波数を有するクロック回路を備える。イベント記述子生成器ＥＤＧにより産生されたイベント記述子ＥＤは、Ｌビットが検出されたイベントの時間タグＴＴ（すなわち、到着時刻Ｔk）を表し、単一の追加ビットがピボットイベントを特定する２進ワードの形態を有する。
【００４１】
イベント記述子フィルタＥＤＦの主要な仕事は、イベント記述子生成器ＥＤＧにより供給されているイベント記述子ＥＤのストリームがパルスバーストを構成するか否かを決定することである。そのため、イベント記述子フィルタＥＤＦは、以下の２つの基準を適用する：
Ａ．連続性基準：連続するイベントの間の時間インターバルは、所定の容認最大値ＩＭ未満でなければならない。すなわち、
Ｔk+1−Ｔk＜ＩＭ
Ｂ．完全性基準：基準Ａを満たす連続するイベントの数ＣＰは、所定の許容最小値ＭＮよりも大きくなければならない。すなわち、
ＣＰ＞ＭＮ
【００４２】
また、連続性基準Ａは、欠落している単一パルスを考慮に入れることもできる。例えば、最大パルス繰り返しインターバルＰＲＩが５ｍｓに等しく、単一パルスが欠落しているかもしれない場合、ＩＭの選択値は、１０ｍｓよりもわずかに大きいものになる。
【００４３】
完全性基準Ｂに関して、ＭＮの最小値は、単一のパルスバーストにおいて発生することが予想される最小パルス数から、許容される欠落パルスの総数を減じたものに等しいものになる。例えば、８つのパルスのみがバースト内で発生することが予想され、合計で３つのパルスの欠落が許容される場合、ＭＮの選択値は４（すなわち、ＣＰ＞４）となる。
【００４４】
検討中のイベント記述子ＥＤのストリームが有効なパルスバーストを構成するものであるか、またはいくつかの重なったパルスバーストを構成するものであると宣言するためには、上記の両方の基準が満たされなければならない。
【００４５】
イベント記述子フィルタＥＤＦは、Ｔk+1−Ｔk＜ＩＭである限り、イベント記述子ＥＤをイベント記述子バッファＥＤＢに連続的に転送する。Ｔk+1−Ｔk≧ＩＭとなり次第、またはバッファＥＤＢが満たされ次第、転送は終了する。バッファが満たされた状況は、例えば、固定レーダに加え、撮像機能を実行する機上システムが、注目対象の領域を照射する場合に発生し得る。
【００４６】
イベント記述子バッファＥＤＢは、ＣＰ＜ＭＮである場合に、入力ＣＬを介してクリアされる。そうでない場合は、バースト全体、またはバーストの組み合わせが、イベント記述子バッファＥＤＢに格納されている。かかる場合、データ利用可能信号ＤＡがイベント記述子フィルタＥＤＦから演算／制御ユニットＡＣＵに送られる。
【００４７】
データ利用可能信号ＤＡの受信により、バッファＥＤＢに格納されたイベント記述子のシーケンスの解析処理が開始される。解析は、演算／制御ユニットＡＣＵおよび繰り返しインターバル推定器ＲＩＥにより行われ、これらは、バッファＥＤＢに格納されたイベント記述子に対して共同で動作する。
【００４８】
イベント記述子の解析には、繰り返し行われる３つの基本的な動作が含まれる：
１．候補のパルス繰り返しインターバルＰＲＩの決定
２．仮定されたパルス繰り返しインターバルＰＲＩの高精度推定
３．仮定されたパルス繰り返しインターバルＰＲＩの検証、および解析中のバーストにおける発生パルス数の決定
【００４９】
それらの３つの動作の実施に成功すれば、ポジティブに識別されているパルスバーストを構成するイベント記述子がバッファＥＤＢから除去される（バッファ「間引き」）。次いで、十分な数のイベント記述子がバッファＥＤＢに残っていれば、これら３つの動作のシーケンスが繰り返される。
【００５０】
以下では、上記の３つの基本的な動作についてより詳細に説明する。
【００５１】
候補のパルス繰り返しインターバルＰＲＩの決定
候補のＰＲＩの予備決定は、イベントパターンの時間領域解析に基づく。好ましくは、用いられるイベントパターンは、少数のイベント（例えば、３または５個）を含む。以下では、かかるイベントパターンを、それぞれイベント３連体（event triple）およびイベント５連体（event quintuple）と呼ぶ。
【００５２】
用いられるイベントパターンの特定の形態に関わらず、まず、好適なピボットイベントを記述子バッファＥＤＢに格納されたイベントから選択する。「良好な」候補のピボットイベントは、格納されたイベントのストリームの中程で見付けられる、と仮定するのが合理的である。
【００５３】
さらに、イベントパターンの特定の形態に関わらず、候補のパルス繰り返しインターバルＰＲＩの探索は、ＰＲＩの最小の予想される値Ｄ（例えば、Ｄ＝ＰＲＩmin＝２００μｓ）から開始する。
【００５４】
イベント５連体
選択された候補のピボットイベントの到着時刻（ＴＯＡ）をＴ*と示す。イベントの４つのグループを以下のように定義する。
レベル１後続者（level-one successors）は、到着時刻Ｔm+が次の条件を満たすすべてのイベントである。
Ｄ＜Ｔm+−Ｔ*＜ｐＤ，ｐ＞１，ｍ＝１，２，．．．，Ｍ
レベル２後続者（level-two successors）は、到着時刻Ｔn++が次の条件を満たすすべてのイベントである。
２Ｄ＜Ｔn++−Ｔ*＜２ｐＤ，ｐ＞１，ｎ＝１，２，．．．，Ｎ
レベル１先行者（level-one predecessors）は、到着時刻Ｔj-が次の条件を満たすすべてのイベントである。
Ｄ＜Ｔ*−Ｔj-＜ｐＤ，ｐ＞１，ｊ＝１，２，．．．，Ｊ
レベル２先行者（level-two predecessors）は、到着時刻Ｔi--が次の条件を満たすすべてのイベントである。
２Ｄ＜Ｔ*−Ｔｉ--＜２ｐＤ，ｐ＞１，ｉ＝１，２，．．．，Ｉ
【００５５】
従って、レベル２のパルスを定義する範囲の上限値および下限値は、レベル１のパルスを定義するために用いられる限界値の２倍である。好ましくは、上記の条件において、パラメータｐの値は、１．５にわずかに満たない値を取る。
【００５６】
選択されたピボットイベントと、１つのレベル１後続者と、１つのレベル２後続者と、１つのレベル１先行者と、１つのレベル２先行者とにより、一意なイベント５連体が形成されることに留意されたい。このように形成されたすべてのイベント５連体の中に、Ｄ＜ＰＲＩ＜ｐＤであるＰＲＩを有するパルスバーストのセグメントを表す１つの５連体（プライミング５連体（priming quintuple）と呼ぶ）が存在し得る。
【００５７】
図５は、１つのピボットイベントと４つのイベントセット（レベル１後続者、レベル２後続者、レベル１先行者、レベル２先行者を含む）から、候補のパルス繰り返しインターバルＰＲＩを予備決定するために用いる方法を示す図である。
【００５８】
本発明の好適な実施形態では、レベル１およびレベル２のイベントを含む４つの時間インターバルを、以下の規則に従って、持続時間（ｐ−１）Ｄを有する単一の時間インターバルに落とし込む。
１．レベル１後続者を含むインターバルは、変更しない。
２．レベル２後続者を含むインターバルは、好適にシフトおよびスケーリングする。
３．レベル１先行者を含むインターバルの持続時間は、変更しない；ただし、ピボットイベントからの距離の絶対値（大きさ）のみを用いる（時間反転）。
４．レベル２先行者を含むインターバルの持続時間は、好適にスケーリングする；ピボット要素からの距離の絶対値は、好適にシフトおよびスケーリングする（時間反転後にスケーリング）。
【００５９】
これにより、ピボットイベントを除くすべてのイベントが単一のインターバル内に配置され、そのインターバルの端点は、想定される候補のＰＲＩの最小値および最大値を表す。イベントの複合ストリームが解析されたものにおいて、上記のインターバルからのある値を有するＰＲＩだけ隔てられたイベントが存在していれば、正確なＴＯＡ測定が行われかつ欠落しているイベントがない理想的な場合には、そのインターバルにおいて同じ位置を占めるイベントが正確に４つあることになる。よって、対応するＰＲＩの高精度な推定値が、図５ｂに示すような好適に構成されたヒストグラムにおける顕著なピークの位置から得られよう。
【００６０】
単一のレベル１イベントが欠落している場合も、かかるイベントセットをさらなる処理に用いることが可能である。
【００６１】
実際の場合、ヒストグラムにおいて観測されるピークは、以下の要因により低減され、平坦化される：
−時間量子化効果（すなわち、時間タグを表すビットは有限数である）
−傍受されたパルスの有限の立ち上がり時間
−雑音による到着時刻ＴＯＡの時間変動
【００６２】
このため、実際には、ピーク位置は、注目対象のＰＲＩの概算の推定値を提供するのみである。しかし、ヒストグラムを用いて、４つの特定のレベル１およびレベル２のイベントを識別することが可能で、それらを既に選択済みのピボットイベントとともに用いて、さらなる処理に要求されるプライミング５連体を構成する。よって、ヒストグラムを好適に解析することで、プライミング５連体に属するレベル１およびレベル２のイベントの４つの基礎となる到着時刻Ｔ0--、Ｔ0-、Ｔ0+、およびＴ0++が決定される。
【００６３】
ヒストグラムがプライミング５連体を識別できなければ、バッファＥＤＢに格納されたピボットのイベント記述子から新たなピボットイベントを選択し、上記の動作を繰り返す。すべての潜在的なピボットを用いたら、バッファＥＤＢに格納された他のイベント記述子の中から候補のピボットを選択することが可能である。このアプローチでも容認可能なプライミング５連体を提供できなければ、ＰＲＩ値の次の範囲、すなわちｐＤ＜ＰＲＩ＜ｐ2Ｄについて探索を繰り返す。例えば、注目対象の範囲全体において、ＰＲＩmin＝２００μｓかつＰＲＩmax＝４ｍｓであるとき、ｐ＝１．４６に対して、ＰＲＩ値のかかる範囲は、８つのパス（eight passes）でカバーされる。（代替的な実施形態では、容認可能なパルスグループが見付けられたかどうかに関わらず、後続のＰＲＩ値の範囲について探索を繰り返す。）最高の範囲からではなく最低の範囲から開始することで、パルス繰り返しインターバルの第２次またはより高次の高調波に到着するパルスのみをグループに割り当ててしまうことが回避される、という利点が得られる。
【００６４】
仮定されたパルス繰り返しインターバルＰＲＩの高精度推定
仮定されたパルス繰り返しインターバルＰＲＩについての情報は、５つのＴＯＡ値（すなわち、プライミング５連体を形成するピボットイベントのＴＯＡ値および他の４つのイベントのＴＯＡ値）に含まれる。それらの値は、次いで、繰り返しインターバル推定器ＲＩＥのブロックにより、仮定されたパルス繰り返しインターバルＰＲＩの値をより高精度に決定するために利用される。
【００６５】
イベント５連体の場合、レベル１およびレベル２のイベントを表す４つの点が、すべて、単位半円（０，π）の円周上にマッピングされる。単位半円は、その中心がイベントの１つ、好ましくはピボットイベントをシンボリックに表し、その角度域が範囲ＰＲＩmax−ＰＲＩmin＝Δを表すものである。（これらの点は、代わりに、他の任意の弧上にマッピングすることもできるが、範囲Δを表す角度域は、π以下とすべきである。）マッピングの結果、（ピボットパルスに対して決定された）パルスインターバルが、単位ベクトルの角度により表される。
【００６６】
検討中のイベント５連体は、ピボットイベント、レベル２先行者、レベル１先行者、レベル１後続者、およびレベル２後続者の到着時刻ＴＯＡの値、すなわち、それぞれＴ*、Ｔ0--、Ｔ0-、Ｔ0+、およびＴ0++により表されるものと仮定する。次いで、単位半円上における各イベントの角度位置を、次式から決定することが可能である：
レベル２およびレベル１の先行者
【数２】
2010525315.tif 000004
レベル１およびレベル２の後続者
【数３】
2010525315.tif 000005
【００６７】
図６は、単位半円の円周上にイベント５連体をマッピングする例である。半円（０，π）上の点は、ピボットイベントを始点とする単位ベクトルの終点を表す。各イベントの角度位置は、それぞれのパルスの到着時刻とピボットパルスの到着時刻との間のパルス間インターバルを測定することにより決定される概算のパルス繰り返しインターバルに対して、線形に関係している。
【００６８】
次に、最も可能性の高いパルス繰り返しインターバルＰＲＩmを、次式のように計算される重み付け角度平均とする：
【数４】
2010525315.tif 000006
式中、
【数５】
2010525315.tif 000007
はベクトルの和の角度であり、また、
【数６】
2010525315.tif 000008
【数７】
2010525315.tif 000009
が成立し、式中、ｗ1およびｗ2は所定の重み係数である。ｗ1＝ｗ2＝１であれば、ＰＲＩmは、各単位ベクトルの和の角度から導出される。しかし、好ましくはｗ1＜ｗ2であり、その場合、ＰＲＩmは、レベル２イベントを用いて測定されたインターバルを表すベクトルの長さが、比ｗ2／ｗ1において、レベル１イベントを用いて測定されたインターバルを表すベクトルの長さよりも大きくなるように調整されたベクトルの和の角度から導出される。この重み付けを行う理由は、狭い間隔のイベントのタイミングの測定に基づいてインターバルを計算する場合に、丸め誤差によりそれらの計算が不正確なものになってしまうからである。好ましくは、レベル２およびレベル１の重みの比（ｗ2／ｗ1）は、２〜４の値を呈する。
【００６９】
加えて、上記のＰＲＩ推定値の確度の測定値ＶPRIを、測定されたインターバルのばらつき(dispersion)から決定することが可能で、例えば、次式のように計算することが可能である：
【数８】
2010525315.tif 000010
ＰＲＩの高品質な推定値については、測定値ＶPRIは、１に非常に近い値を呈する。例えば、次式が成立する場合：
【数９】
2010525315.tif 000011
レベル１およびレベル２のイベントは、単一の点により表され、ＶPRIの値は１に等しい。しかし、「最も不調和なイベント」（図７に示すもの等）については、ＶPRIの値は０に等しい。
【００７０】
ＰＲＩ推定値の確度の測定値ＶPRIの低い値を独立に用いて、既に選択済のイベント５連体を却下することができる。かかる場合、他の５連体の探索を繰り返す。他方、測定値ＶPRIが所定の閾値を超える場合、イベントのグループのさらなる処理を後述のように行う。さらなる処理には、同じグループに属する他のイベントを位置特定することと、残りのパルスを他のグループに分類する目的での検討から、受信したパルストレインからのグループのイベントを排除することとが含まれる。
【００７１】
図８は、繰り返しインターバル推定器ＲＩＥのブロック図である。ブロックＲＩＥは、データセレクタＤＳＲ、ピボットレジスタＰＴＲおよびイベントレジスタＥＴＲ、合算／スケーリング回路ＳＳＲ、ルックアップテーブルＴＳＣ、重みレジスタＥＷＲ、２つの同一の乗算器Ｍ、３つのアキュムレータＡＣＣ、ＡＣＳ、およびＡＣＷ、ならびにピタゴラスプロセッサＰＹＴを備える。
【００７２】
演算／制御ユニットＡＣＵは、イベント５連体を選択し、対応する記述子のバッファＥＤＢからＲＩＥブロックへのポートＳＤを介する転送を制御する。また、イベントのレベルおよびタイプ（ピボットか否か）は、ユニットＡＣＵの出力ＥＬを介して、データセレクタＤＳＲ、合算／スケーリング回路ＳＳＲ、および重みレジスタＥＷＲへ送られる。その結果、ピボット記述子（すなわち、そのＴＯＡ）ＰＴはレジスタＰＴＲに格納され、イベント記述子ＥＴはレジスタＥＴＲに格納される。それらのレジスタの出力ＴＰおよびＴＥは、合算／スケーリング回路ＳＳＲに接続されている。合算／スケーリング回路ＳＳＲは、その出力ＤＴにおいて以下を産生する：
−レベル１イベントについては、ピボットＴＯＡとイベントＴＯＡとの間の差の絶対値
−レベル２イベントについては、ピボットＴＯＡとイベントＴＯＡとの間の差の絶対値の半分
【００７３】
値ＤＴは、ルックアップテーブルＴＳＣに格納された正弦関数および余弦関数の要求される値を得るために用いられる、好適な構造のアドレスである。要求される値Ｃ２およびＳ２は、適切な重みＥＷとともに２つの同一の乗算器Ｍにより用いられ、２つのアキュムレータＡＣＣおよびＡＣＳへの入力データを産生する。重みＥＷは、重みレジスタＥＷＲにより供給される。重みレジスタＥＷＲは、イベントレベルに関する情報ＥＬを演算／制御ユニットＡＣＵから受信する。イベント重みＥＷは、また、重みアキュムレータＡＣＷによっても用いられる。すべてのアキュムレータは、演算／制御ユニットＡＣＵにより供給される信号ＲＳにより同時にリセットされる。
【００７４】
完全なイベント５連体を処理している（かつ重みの和が既知である）ときには、アキュムレータＡＣＷは冗長なものであるが、イベント５連体から単一のレベル１イベントが欠落しているときには、単純な形態のアキュムレータＡＣＷが必要とされる。同様に、レーダ信号解析器ＲＳＡの動作中、イベント５連体をイベント３連体で置換するときは、重みの和は異なるものとなる。
【００７５】
アキュムレータＡＣＣおよびＡＣＳの出力ＳＣおよびＳＳは、ピタゴラスプロセッサに接続されている。ピタゴラスプロセッサは、その組み合わせ出力ＥＩ／ＥＣにおいて、推定されたパルス繰り返しインターバルＰＲＩmの値とその推定値に関連付けられた確度測定値ＶPRIとの両方を産生する。両方の値ＰＲＩmおよびＶPRIは、さらなる処理のため演算／制御ユニットＡＣＵにより利用される。
【００７６】
仮定されたパルス繰り返しインターバルＰＲＩの検証、およびバースト当たりパルス数の決定
ＲＩＥブロックにより供給される推定値ＰＲＩmは、演算／制御ユニットＡＣＵにより、この特定ＰＲＩにより特徴付けられるバースト内で発生するパルスの数を決定するために利用される。この目的のため、単純なインパルステンプレート（impulse template）を用いて、複数の連続フレーム、および各フレーム内でのイベントの相対位置を決定する。
【００７７】
図９は、演算／制御ユニットＡＣＵにより行われる動作を概略的に示す。ユニットＡＣＵは、受信信号を所定長の時間フレームに分割することにより動作を行う。各時間フレーム内で発生するイベントを合成して、ヒストグラムを産生する。ヒストグラムから、イベントカウント、すなわちバースト当たりパルス数を決定することが可能である。イベント統合に用いられる初期のフレームは、プライミング５連体（または３連体）からのイベントを含むものとする。そして、他のフレーム（先行するものと後続のものとの両方）について検討するとともに、得られたヒストグラムをパルスの連続性について監視する。例えば、２つの連続するフレームが特定のパルス繰り返しインターバルのパルスについてイベントカウントを増加させない場合、それぞれのバーストが終了したとの決定を行ってもよい。
【００７８】
このように、ＡＣＵは、受信トレインにおける他のパルスのどれが識別されたグループに属するかを判定することが可能で、次いで、そのグループにおけるパルスを、システムが他のグループの識別を試みるときに検討対象から排除することが可能である。
【００７９】
上記の方法は、同じＰＲＩを有する２つ（またはそれ以上）のインターリーブされたバーストを、開始時刻がシフトしていても識別する能力を有することを指摘しておく。
【００８０】
レーダパルスバーストがない時間インターバルの予測
演算／制御ユニットＡＣＵは、その動作の結果、情報を格納することとなり、その情報から、傍受されたレーダ送信の各々について、以下のパラメータを推測することが可能である：
１．パルス繰り返しインターバルＰＲＩ
２．バースト当たりパルス数
３．バースト持続時間
４．バースト到着時刻
５．バーストとバーストの間のインターバルの持続時間
【００８１】
上記の情報は、演算／制御ユニットＡＣＵが、当業者により直ちに理解される手法を用いてレーダパルスバーストがない時間インターバルを予測するために十分である。このように、ＡＣＵは、これらの時間インターバルを表す送信許可信号ＣＴを提供することが可能である。また、好適な予測手順を、傍受したパルス信号の連続的監視と組み合わせることで、予測品質を向上させ、新たな発信器が出現する場合も適時に反応することができる。
【００８２】
送信許可信号ＣＴは、無線送信機により、無線送信機からの送信タイミングを制御するために用いられる。図１の情報融合センタＩＦＣおよびセンサにより形成される通信システムにおいて、信号はセンタＩＦＣに送られ、センタＩＦＣは、それに応答して、センタＩＦＣからセンサへと、およびセンサからセンタＩＦＣへと、情報を送信するタイミングを制御する。
【００８３】
また、上で説明した手法は、欧州特許出願第０５２５６５７５．１号（２００５年１０月２４日出願）および国際出願第ＰＣＴ／ＧＢ２００６／００３９４２号（２００６年１０月２３日出願）においても記載および請求されている。上で説明した手法に行うことができる具体的な変形の例を以下で説明する。
【００８４】
上述の実施形態では、候補のパルスの選択後に、プライミング５連体（すなわち、実質的に等しいパルス繰り返しインターバルを有する５つのパルスであり、候補のパルスは中央のパルスとする）を位置特定する処理を行うように編成されている。示したように、５つのパルスを位置特定する代わりに、３つのパルスで構成される組を位置特定することが可能である。加えて、（ｉ）この処理は、候補のパルスが中央にあるものとしてパルスシーケンスを位置特定するように試みる必要はなく、また、（ｉｉ）この処理は、奇数個のパルスを含むシーケンスを位置特定する必要もない。例えば、用途によっては、４つのパルスを含むグループを解析することが有用であるかもしれない。従って、図５および６に関連して上で説明した処理を、可能性のあるレベル２後続者またはレベル２先行者を無視するように変更することで、プライミング４連体を形成する４つのパルスのグループを求め、それらのパルスに基づいて最も可能性の高いパルス繰り返しインターバルを計算することが可能である。
【００８５】
図５を参照して上で述べたように、ヒストグラムにおけるピークは、時間量子化効果や、到着時刻測定値における雑音を原因とする変動などの要因の影響を受けることがある。また、ピークは、パルスをそれらの公称位置（nominal positions）から故意に変位（ある種類のパルストレインにおいて生じる）させることによっても影響を受ける。これらの要因のいずれの１つまたはいずれの組み合わせにも対処するため、本発明の変形例によれば、候補のパルスに対する各パルスの到着時刻を決定する代わりに、パルスの相対的な到着時刻の平均値を決定する。これは、好ましくは、好適に選択されたスパンの移動窓を用い、窓内に入るイベントをカウントすることにより果たされる。
【００８６】
図１０は、この手法の例を概略的に説明する。この変形手法を用いる実施形態は、候補のパルス繰り返しインターバルを予備決定する方法が、図５でなく図１０を参照して説明するものであることを除き、上記の実施形態と同一である。図１０ａは、図５ａと同様に、適切なシフトおよびスケーリング操作後のパルスイベントの測定された到着時刻を表す。各イベントには、そのイベントの時刻と候補のパルスの到着時刻との間の絶対誤差を表す値（必要に応じ、所定のインターバルＰＲＩmin〜ＰＲＩmax内に配置されるように整数で除算する）が関連付けられている。図１０ａには、また、移動窓のスパンも示している。窓を漸次（窓幅未満の量だけ）シフトし、その各位置で窓内のイベントをカウントする。図１０ｂは、その結果、すなわち、移動窓をシフトさせながら窓内で観測されたイベントの数を示す。図１０は、パルスタイミングが故意にずらされたパルス４連体を検出するための例に関する。故に、合成されたストリームに一致するイベントはない。しかし、所定のスパンの窓を時間的にシフトさせると、ある窓位置において３つのイベントのクラスタが窓内で観測される。これにより、４連体の検出が宣言され、窓の対応する中央位置から平均ＰＲＩが決定される。
【００８７】
図１０ｂは、イベントカウントピークをイベントクラスタの真の位置に対して窓スパンの半分だけシフトさせた状態を示す。これは、図がイベント数を窓の前縁の関数として示すためである。単一の窓を用いる代わりに、２つの逆方向に伝播する窓（１つはＰＲＩminから開始し、１つはＰＲＩmaxから開始する）を用いることも可能である。
【００８８】
また、図９を参照して説明した手法も、１つ以上の移動窓を用いるように対応する方法で変形することができる。
【００８９】
上記の実施形態において、パルスのグループは、選定された候補のパルスの到着時刻に対するそれらの到着時刻により選択される。状況によっては、観測されたパルスの到着時刻の代わりに、ある候補となる瞬間を参照点として選択することが望ましいかもしれない。換言すれば、候補のパルスとして、実際の（観測された）パルスの代わりに、仮想の（観測不能な）パルスを用いることが可能である。これは、例えば、ずらされたパルスを有する短いパルスバーストの解析において有益であり得る。候補となる瞬間は、例えば、観測されたパルスの到着時刻（例えば、２つの到着時刻の間の中間点）によって選択してもよく、好ましくは、グループを形成する候補として考えられるパルスの到着時刻の範囲の中央であるか、またはその付近である。代替として、各々が直前の候補となる瞬間から所定の時間インターバルだけ隔てられた異なる候補である瞬間を、繰り返し選択するための反復的処理であってもよい。観測パルスの到着時刻の代わりに候補となる瞬間を用いる以外は、先に説明した実施形態に関して上で述べたものとまったく同じ動作であってもよい。
【００９０】
以上から、観測されたパルスの到着時刻と一致しない選択された候補となる瞬間を、候補のパルスの到着時刻の代わりに用いることで、（ｉ）パルスをグループに類別し、および／または（ｉｉ）概算の繰り返しインターバルを決定し、次いで、それらを用いて平均グループパルス繰り返しインターバルを計算することが可能である、ということが理解されよう。
【００９１】
本発明の好適な実施形態の上記の説明は、例示および説明の目的で提示したものである。網羅的であることや、本発明を開示したとおりの形態に限定することを意図するものではない。上記の説明に鑑み、多くの改変、修正、および変更を行うことで、当業者が、考慮される特定の用途に適した様々な実施形態において本発明を利用することが可能になる、ということは明白である。
【図１】
2010525315.tif 000012
【図２】
2010525315.tif 000013
【図３】
2010525315.tif 000014
【図４】
2010525315.tif 000015
【図５】
2010525315.tif 000016
【図６】
2010525315.tif 000017
【図７】
2010525315.tif 000018
【図８】
2010525315.tif 000019
【図９】
2010525315.tif 000020
【図１０】
2010525315.tif 000021
【国際調査報告】
2010525315.tif 000022
2010525315.tif 000023
2010525315.tif 000024
2010525315.tif 000025
2010525315.tif 000026