ATCコード(ユニークコード)は航空機AP1,AP2のフライトごとに管制官により指定され、各航空機AP1,AP2の操縦士が応答装置に入力する。示す情報はQNE(圧力補正無しの高度)であり、これは100ft単位で更新される。QNEは地上局のコンピュータによりQNH(圧力補正された高度)に変換される。各航空機AP1,AP2の応答装置は、各二次監視レーダSSR1,SSR2からの質問信号(1030MHz帯)を受信すると、予め決められた遅延時間をおいて応答信号(1090MHz帯)を送信する。このとき応答信号のパターンは質問信号のパターンに合致している。このため、一方の二次監視レーダSSR1からの質問信号に対する応答信号は以下のパターンとなる。(1)A,A,C,A,A,C,・・・また、他方の二次監視レーダSSR2からの質問信号に対する応答信号は以下のパターンとなる。マイクロホン−アンテナユニット10のレーダ受信アンテナ12を通じて各航空機AP1,AP2から発せられる応答信号を受信することができる。マイクロホン−アンテナユニット10の騒音計測用マイクロホン14、及び音の到来方向識別用マイクロホン16を用いて各航空機AP1,AP2に由来する航空騒音を観測することができる。上記の観測ユニット100は、例えば図示しない中央演算処理装置(CPU)や半導体メモリ(ROM、RAM)、ハードディスクドライブ(HDD)、入出力インタフェース、液晶ディスプレイ等を備えたコンピュータ機器で構成されている。観測ユニット100は、そのハードウェア資源をソフトウェア制御により複数の機能要素として機能させている。複数の機能要素は、主に応答信号の処理系統と騒音測定系統とに分かれている。以下、観測ユニット100内部で実現されている機能要素について説明する。レーダ受信アンテナ12に通じる受信部20や目的信号群抽出部22等の機能要素が含まれる。目的信号群抽出部22には、一時メモリ24及び時間間隔分析部26(時間間隔算出部)が含まれている。受信部20は、レーダ受信アンテナ12を通じて受信した応答信号(1090MHz帯)をパルス形式に復調し、所定ビット長(例えば13ビット)のフレームを生成する。生成されたフレームは、パルス形式の応答信号として一時メモリ24に蓄積される。一時メモリ24には、所定時間内(例えば1秒間)に受信された応答信号が時系列順に蓄積される。時間間隔分析部26は、蓄積された応答信号の全てについて相互の時間間隔を算出する。ヒストグラム作成部28(頻度算出部)、フィルタ部30等の機能要素が含まれている。ヒストグラム作成部28は、時間間隔分析部26で算出された複数通りの時間間隔を統計データとして、ヒストグラムを作成する。ヒストグラムは可視的なものである必要はなく、物理メモリ空間上でヒストグラムとして認識可能なものであればよい。【0043】フィルタ部30は、ヒストグラム上で最も頻度が高い時間間隔をタイムゲート(篩の目)として、一時メモリ24に蓄積された応答信号から当該時間間隔でペアを作れない応答信号を除去する。これにより、タイムゲートにマッチする特定の応答信号群だけが残り、その他の応答信号は全て篩い落とされる。32(情報算出部)、時刻データ作成部36(データ配列作成部)等の機能要素が含まれている。C」の内容別に電波強度の平均(平均強度)及び出現回数を算出する。C」の平均強度及び出現回数について、これらを受信時刻(例えば1秒間)ごとに配列した時刻データを作成する。騒音測定系統には、騒音計測用マイクロホン14に通じて騒音測定部38が含まれている。ここで騒音計測用マイクロホン14は、騒音計測用のマイクロホンであり、音の到来方向識別用マイクロホン16は4つのマイクロホンは3軸相関法による上空音識別用のものである。騒音測定部38は、3軸相関法により航空機AP1,AP2(音源)の仰角θ及び方位角δを算出し、観測点を基準とした3軸の観測空間(ベクトル空間)内で騒音の到来方向ベクトル(単位ベクトル)を算出する。騒音測定部38は、算出した到来方向ベクトルの値(軌跡)により、観測点を基準として航空機AP1,AP2(音源)の移動方向(どの方角からどの方角へ向かったか)を算出する。騒音測定部38は、観測点での騒音レベルがある閾値を超えた場合(騒音イベント発生時)、その時刻と仰角θ及び方位角δを関連付けて保存する。観測ユニット100には、整合部40の機能要素が含まれている。整合部40は、信号処理系統で作成された時刻データと騒音測定系統で作成された騒音データとの整合(関連付け)をとる。C」の情報を取得する。取得された情報はデータメモリ42に保存され、さらに出力部44にて出力される。出力部44は、上記の液晶ディスプレイや外部インタフェース等の構成を含む。外部インタフェースにコンピュータ等を接続することで、データメモリ42に保存された情報をコンピュータ等で利用することが可能となる。構成例と機能についての説明であるが、以下に観測ユニット100により行われる応答信号処理について具体例を挙げて説明する。second)でA,A,C,A,A,C,・・・のパターンの質問信号を放射しているとする。管制塔CT2の二次監視レーダSSR2は、474.3ppsでA,C,A,C,・・・のパターンの質問信号を放射しているとする。【0053】ここで各二次監視レーダSSR1,SSR2は、東西南北の全方位(360°)をスキャンするため、非常に指向性の鋭いアンテナを3〜5秒周期で1回転させている。したがって各航空機AP1,AP2は、各二次監視レーダSSR1,SSR2のアンテナと相対した瞬間に質問信号を受信し、かつ、予め決められた時間遅れ3.5μsで応答信号を発している。【0054】ここで、各航空機AP1,AP2が応答する時間は、各二次監視レーダSSR1,SSR2の指向性角度と回転周期から計算することができ、一般的には、概ね30〜50msの時間内で応答することができる。二次監視レーダSSR1に対する応答数は11〜18回程度となり、二次監視レーダSSR2に対する応答数は14〜24回程度となる。図3は、上記の運用条件で航空機AP1,AP2から発せられる応答信号のパターンを時系列で示した図である。図3中(A)に示されているように、航空機AP1は、ある時刻m1に二次監視レーダSSR1と相対し、その瞬間から質問信号(AACAAC・・・)を有効に受信できる時間内で11〜18パルスの応答信号(AACAAC・・・)を発している。次に航空機AP1が二次監視レーダSSR1と相対するのは、回転周期T1後の時刻m5であり、この瞬間から航空機AP1は再び11〜18パルスの応答信号(AACAAC・・・)を発することになる。一方で航空機AP1は、別の時刻m3に二次監視レーダSSR2とも相対し、その瞬間から質問信号(ACAC・・・)を有効に受信できる時間内で14〜24パルスの応答信号(ACAC・・・)を発する。次に航空機AP1が二次監視レーダSSR2と相対するのは回転周期T2後の時刻m7であり、この瞬間から航空機AP1は再び14〜24パルスの応答信号(ACAC・・・)を発することになる。別の航空機AP2については以下のようになる。図3中(B)に示されているように、航空機AP2は、ある時刻m2に二次監視レーダSSR2の方と相対し、その瞬間から質問信号(ACAC・・・)を有効に受信できる時間内で14〜24パルスの応答信号(ACAC・・・)を発している。次に航空機AP2が二次監視レーダSSR2と相対するのは、回転周期T2後の時刻m6であり、この瞬間から航空機AP2は再び14〜24パルスの応答信号(ACAC・・・)を発することになる。一方で航空機AP2は、別の時刻m4に二次監視レーダSSR1とも相対し、その瞬間から質問信号(AACAAC・・・)を有効に受信できる時間内で11〜18パルスの応答信号(AACAAC・・・)を発する。図3には示されていないが、次に航空機AP2が二次監視レーダSSR1と相対するのは回転周期T1後であり、その瞬間から航空機AP2は再び11〜18パルスの応答信号(AACAAC・・・)を発することになる。図4は、所定時間内(例えば1秒間)に受信された全ての応答信号を時系列に並べた例を示す図である。このような応答信号S1〜Spの配列は、上記の一時メモリ24に蓄積された1秒間ごとの応答信号を並べた結果(13ビットのコードと時刻、及び強度を含む)を表している。上記のように時間間隔分析部26は、全ての応答信号S1〜Spについて、相互の時間間隔tmnを算出する。図4中(A)に示されているように、先頭の応答信号S1であれば、他の全ての応答信号S2〜Spと(いわゆる「総当たり」)の時間間隔t12,t13,t14,・・・t1pを算出する。図4中(B)に示されているように、2番目の応答信号S2であれば、他の全ての応答信号S1,S3〜Spとの時間間隔t21,t23,t24,・・・t2(p-1)を算出する。図4中(C)に示されているように、3番目の応答信号S3であれば、他の全ての応答信号S1,S2,S4〜Spとの時間間隔t31,t32,t34,・・・t3(p-2)を算出する。図示していないが、以下同様に応答信号S4〜Spまで、全ての時間間隔tmnを算出する。図5は、算出した時間間隔から得られたヒストグラムの例を示す図である。ここでは便宜上、可視化して示しているが、上記のようにヒストグラム作成部28は、メモリ空間上で認識可能なヒストグラムを作成する。便宜上、ここでは個々の時間間隔を1秒間の受信パルス数(pps)で表している。次に、得られたヒストグラムから、最も頻度が高い(受信数が多い)時間間隔tgateを算出し、この算出した時間間隔tgateをフィルタリング用のタイムゲートとする。ここでは最頻値である「370pps」がタイムゲートの第1候補となり、次に頻度が高い「474.3pps」は第2候補となる。図6は、タイムゲートを用いて応答信号をフィルタリングした場合の一例を示した図である。上記のフィルタ部30は、最も頻度が高い時間間隔tgateをタイムゲートとして、目的とする応答信号群をフィルタリングする。仕様により第1候補と第2候補の両方を扱う場合もあるが、本実施形態では、第1候補のみを使用する。図6中(A)に示す1秒間に受信された全ての応答信号S1〜Spについて、図6中(B)〜(D)に示す時間間隔tgateをタイムゲートとしてペアとなり得る応答信号の組を抽出してフィルタリングすると、図6中(B),(D)に示すように、目的の応答信号群だけが残り(例えばS1,S4,S6)、図6中(C)に示すように、その他の応答信号は除去される。応答信号をペアにすることで、ノイズを除去することができる。観測ユニット100により実行される騒音監視処理の手順例を示すフローチャートである。用いられる各ステップの内容が明らかとなる。観測ユニット100は、騒音監視処理を1秒ごとの時間割り込みで実行することができる。【0067】ステップS10:観測ユニット100は、受信部20において応答信号受信処理を実行する。この処理では、上記のように1秒間ごとの応答信号を受信する。22においてデータ整列処理を実行する。この処理では、受信した応答信号に時刻と強度を関連付けて一時メモリ24に蓄積し、時系列にデータを整列させる(図4)。【0069】ステップS30:また観測ユニット100は、時間間隔分析部26において時間間隔分析処理を実行する。この処理では、上記のように全ての応答信号S1〜Spについて、相互の時間間隔tmnを算出する(図4)。【0070】ステップS40:次に観測ユニット100は、ヒストグラム作成部28においてヒストグラム処理を実行する。この処理では、上記のように複数通りの時間間隔tmnからヒストグラムを作成する。【0071】ステップS50:そして観測ユニット100は、フィルタ部30においてフィルタ処理を実行する。この処理では、上記のように最大の頻度である時間間隔tgateを算出し、これをタイムゲートとして目的の応答信号群をフィルタリングする(図6)。このときフィルタ部30は、全ての応答信号の中から、時間間隔tgateに一致する応答信号のペアを抽出する。32においてペア情報算出処理を実行する。この処理では、フィルタリング後の応答信号のうち、内容(コード)が同じであるペアの電波強度の平均及び出現回数を算出する。【0073】ステップS70:また観測ユニット100は、モード決定部34において、モード判定処理を実行する。C」のいずれに該当するかを判定する。モード判定処理の詳細については、別のフローチャートを用いてさらに後述する。【0074】ステップS80:また観測ユニット100は、時刻データ作成部36においてデータ配列処理を実行する。C」の別、電波強度及び出現回数を配列した時刻データを作成する。【0075】ステップS90:そして観測ユニット100は、整合部40において整合処理を実行する。C」の情報を取得する。図8は、モード判定処理の内容を詳細に示すフローチャートである。C」のいずれのコードに該当するかを判定する。判定済みの応答信号について、コードの違いにより「A1,A2,A3,・・・」又は「C1,C2,C3,・・・」等のラベリングを行う。以下、手順に沿って説明する。【0077】ステップS100:モード決定部34は、抽出したペアの応答信号が同じコード(例えば13ビットコード)を有するものであるかを確認する。抽出された応答信号のペアは理論上、「A,A」、「A,C」又は「C,A」のいずれかとなるため、ペアのコードが同じ場合、それは「A,A」のペアである(Yes→ステップS102へ)。ステップS104:またモード決定部34は、コード内容別(機体識別コード情報別)に「A1,A2,A3,・・・」としてラベリングする。ラベリングしたコードは適宜、図示しない内部メモリに保存され、必要に応じて参照することができる。【0079】ステップS106:この段階で、未だフィルタリング後のデータを全て抽出していなければ(No)、モード決定部34はステップS100に戻る。ステップS100:そしてモード決定部34は、次に抽出したペアのコードが同じであるかを確認する。ペアのコードが同じでない場合(No)、モード決定部34はステップS108に進む。【0080】ステップS108:コードが同じでないペアは、「A,C」又は「C,A」のいずれかである。C」である(Yes→ステップS110へ)。ステップS102:そしてモード決定部34は、コードの内容別(高度情報別)に「C1,C2,C3,・・・」としてラベリングする。ラベリングしたコードは、上記のように図示しない内部メモリに保存され、必要に応じて参照することができる。以上のようにして、1秒間の応答信号について「A1,A2,A3,・・・」又は「C1,C2,C3,・・・」のラベリングを行い、全データを抽出し終えると(ステップS106)、モード決定部34は呼び出し元の騒音監視処理(図7)に復帰する。C」も判定することができる。C」の判定を行うことができない。この場合、モード決定部34は以下の手順を実行する。【0084】ステップS100:すなわちモード決定部34は、新たに抽出したペアのコードが同じでないことを確認すると(No)、次にステップS108に進む。C」を判定することができない(No→ステップS114へ)。【0086】ステップS114:この場合、モード決定部34はペアのコードをひとまずバッファに記憶する。ステップS116:そしてモード決定部34は、記憶したコードの時間的な変化傾向を分析する。【0087】ステップS118:変化傾向の分析には相当数のサンプルが必要となるので、数組程度のサンプルでは変化傾向が顕著に表れない。したがってモード決定部34は、現状のサンプル数で充分に変化傾向を分析できない場合(No)、ステップS106に進んで全データを抽出済みか確認する。全データを抽出していなければ(ステップS106:No)、モード決定部34はステップS100に戻って処理を継続する。【0088】そして、次に抽出したペアのコードも同じでなく(ステップS100:No)、かつ、ペアの中に「A1,A2,A3,・・・」としてラベリングされているものがない場合(ステップS108:No)、モード決定部34は次のペアもバッファに記憶する。【0089】なお、ステップS106で1秒間の全データを抽出済みである場合(Yes)、ここで呼び出し元の騒音監視処理(図7)に復帰するが、ペアのコードが「A,C」又は「C,A」であり、ペアの中に「A1,A2,A3,・・・」としてラベリングされているものがない場合、次回の呼び出しでもステップS100(No)→ステップS108(No)→ステップS114の手順が実行される。【0090】このようにして、相当数のサンプルがバッファに記憶されてくると、コードの時間的な変化傾向が明らかになってくる。9は、コードの時間的な変化傾向の分析例を示す図である。C」であれば高度情報であるが、バッファに記憶した時点ではmodeが不明であるため、ひとまずペアのコードは「データ(13ビット値)」として認識される。C」であると容易に判定することができる。)に「An,An+1,An+2,・・・」、「Cn,Cn+1,Cn+2,・・・」としてラベリングする。このようにしてラベリングを行うと、モード決定部34は呼び出し元の騒音監視処理(図7)に復帰する。図10は、時刻データと騒音データとを整合させた場合の例を示す図である。図10中(A)は時刻データの構造例を示し、図10中(B)は騒音データの構造例を示している。配列した構造を有する。配列した構造を有している。騒音データには、「仰角」及び「方位角」が含まれるため、これらの情報をもとに観測点からみた騒音源の到来方向を知ることができる。C」の情報を取得すれば、「機体識別コード情報」や「高度情報」から騒音源(航空機AP1,AP2等の固有情報)との照合の精度向上を図ることができる。本発明は上述した実施形態に制約されることなく、種々に変形して実施することができる。飛行場以外でも利用することができる。【0098】また、観測ユニット100の構成は一例であり、複数の機能要素を1つに統合したり、あるいは1つの機能要素を複数に分割して観測ユニット100を構成してもよい。ekouhou.net。