JPH0288928A - Device for discriminating noise of aircraft - Google Patents

Device for discriminating noise of aircraft

Info

Publication number
JPH0288928A
JPH0288928A JP24204488A JP24204488A JPH0288928A JP H0288928 A JPH0288928 A JP H0288928A JP 24204488 A JP24204488 A JP 24204488A JP 24204488 A JP24204488 A JP 24204488A JP H0288928 A JPH0288928 A JP H0288928A
Authority
JP
Japan
Prior art keywords
noise
data
aircraft
cross
correlation
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
Application number
JP24204488A
Other languages
Japanese (ja)
Other versions
JPH0515971B2 (en
Inventor
Noriaki Hayashi
林 範章
Ichiro Yamada
一郎 山田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Rion Co Ltd
Kobayashi Institute of Physical Research
Original Assignee
Rion Co Ltd
Kobayashi Institute of Physical Research
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Rion Co Ltd, Kobayashi Institute of Physical Research filed Critical Rion Co Ltd
Priority to JP24204488A priority Critical patent/JPH0288928A/en
Publication of JPH0288928A publication Critical patent/JPH0288928A/en
Publication of JPH0515971B2 publication Critical patent/JPH0515971B2/ja
Granted legal-status Critical Current

Links

Abstract

PURPOSE:To surely discriminate the noise of aircraft without receiving any influences of periodic high-level noise components by performing white noising calculating processes on noise detecting signals obtained from a pair of microphones and at the same time, executing a mutually correlating process on the signals. CONSTITUTION:Difference data S5A and S5B from which periodic noise components contained in noise detecting data S4A and S4B are removed are found by performing difference calculating processes on the noise detecting data S4A and S4B obtained from two microphones 2A and 2B by means of white noising calculating means 15 and mutually correlating processes are performed on the difference data S5A and S5B. Therefore, wrong discrimination based on a periodic high-level noise component can be avoided and the aircraft noise discriminating ability of this device can be improved further.

Description

【発明の詳細な説明】 〔発明の概要〕 本発明は、相互相関演算方式の航空機騒音識別装置にお
いて、騒音検出データを差分演算処理した後相互相関演
算を実行するようにしたことにより、航空機騒音を表す
ピークを一段と急峻にし得る。
Detailed Description of the Invention [Summary of the Invention] The present invention provides an aircraft noise identification device using a cross-correlation calculation method, in which a cross-correlation calculation is executed after noise detection data is subjected to differential calculation processing. It is possible to make the peak representing the value even steeper.

〔従来の技術〕[Conventional technology]

航空機騒音を評価する基準として、[航空機騒音に係る
環境基準(昭和48年12月27日環境庁告示154)
 Jがあり、当該基準に従って製造された航空機騒音自
動測定装置が従来から用いられている。
As a standard for evaluating aircraft noise, [Environmental standards related to aircraft noise (Environment Agency Notification 154 of December 27, 1970)]
J, and automatic aircraft noise measurement devices manufactured in accordance with this standard have been used in the past.

ところがこの種の航空機騒音自動測定装置においては、
航空機騒音を例えば市街騒音のような周囲の騒音と区別
して航空機騒音だけについて評価をすることが実際上困
難で、この問題を解決するための手段として従来、第1
にマイク0ホンで受音された騒音が予め決められた設定
レベルを超えたとき当該騒音を航空機騒音として評価す
る方法(特公昭60−15143号公報)や、第2に互
いに離間して配置したマイクロホンの出力信号について
相互相関値を求め、当該相互相関値のピークの発生位置
に基づいて航空機騒音を識別する方法(特公昭61−1
3169号公報)などが提案されている。
However, in this type of automatic aircraft noise measurement device,
In practice, it is difficult to distinguish aircraft noise from ambient noise such as city noise and evaluate only aircraft noise.
The second method is to evaluate the noise as aircraft noise when the noise received by microphone 0 exceeds a predetermined set level (Japanese Patent Publication No. 15143/1983), and A method for determining the cross-correlation value of the microphone output signal and identifying aircraft noise based on the position of the peak of the cross-correlation value (Special Publication No. 61-1
3169), etc. have been proposed.

〔発明が解決しようとする問題点〕[Problem that the invention seeks to solve]

ところが第1の識別方法は、航空機騒音以外の騒音のレ
ベルが航空機騒音より大きい場合には当該航空機騒音以
外の騒音を誤って航空機騒音であると識別するおそれが
あり、この問題は第2の識別方法を併用することによっ
て改善するようになされている。しかし当該第2の識別
方法は、測定すべき騒音の周波数スペクトルに周期性が
強い騒音成分が含まれている場合には、当該周期性に基
づくノイズが相互相関値に混入することにより誤動作す
るおそれがある。
However, with the first identification method, if the level of noise other than aircraft noise is higher than aircraft noise, there is a risk that the noise other than aircraft noise may be mistakenly identified as aircraft noise.This problem can be solved by the second identification method. Improvements are being made by using methods in combination. However, if the frequency spectrum of the noise to be measured contains noise components with strong periodicity, the second identification method may malfunction due to the noise based on the periodicity being mixed into the cross-correlation value. There is.

囚に第2の方法によって2つの測定信号間の相互相関値
を求める方法は、騒音が白色雑音としての性質をもつと
考えた場合当該騒音のトータル音量についての相互相関
値のピークの位置の時間的変化が航空機の飛行位置を表
している点に着目して航空機騒音を識別するようになさ
れている。
The second method for calculating the cross-correlation value between two measurement signals is to calculate the time of the peak position of the cross-correlation value for the total volume of the noise, assuming that the noise has the characteristics of white noise. Aircraft noise is identified by focusing on the fact that the change in position represents the flight position of the aircraft.

ところが実際上例えばプロペラを有する航空機の場合に
は100〜200 ([2)程度の周波数範囲に周期性
が強くかつ騒音レベルが大きい騒音成分を含んでおり、
また市街雑音の中にサイレン、クラクションなどが混入
すると、250 (Hz) 〜1 (kHz)の間の周
波数を有する周期性が強くかつ騒音レベルが大きい騒音
成分が生ずる。
However, in reality, for example, aircraft with propellers contain noise components with strong periodicity and high noise levels in the frequency range of about 100 to 200 ([2)].
Furthermore, when sirens, horns, etc. are mixed into city noise, a noise component having a frequency between 250 (Hz) and 1 (kHz) with strong periodicity and a high noise level is generated.

このように周波数スペクトラムのうち一部の周波数につ
いて騒音レベルが大きい騒音成分が混入すると、当該騒
音成分の騒音レベルの変動周期が相互相関演算周期より
十分に長い場合には、相互相関値が航空機の移動に基づ
いて変化すると同時に、当該騒音成分の周期で変化する
ことになり、これが航空機騒音の識別結果に雑音として
混入する結果になる。
In this way, when a noise component with a high noise level for some frequencies in the frequency spectrum is mixed, if the fluctuation period of the noise level of the noise component is sufficiently longer than the cross-correlation calculation period, the cross-correlation value of the aircraft At the same time as it changes based on movement, it also changes with the period of the noise component, and this results in being mixed in as noise in the aircraft noise identification result.

このような雑音を含んだ識別結果に基づいて航空機騒音
であるか否かを判定することは、実用上極めて困難であ
り、改善することが望ましい。
It is extremely difficult in practice to determine whether or not it is aircraft noise based on identification results that include such noise, and it is desirable to improve this.

本発明は以上の点を考慮してなされたもので、マイクロ
ホンにおいて検出された騒音信号に周期性騒音成分が混
入している場合にも、その影響を受けない相互相関値を
得ることができるようにすることにより、航空機騒音の
識別性能を一段と向上させ得るようにした航空機騒音識
別方法を提案しようとするものである。
The present invention has been made in consideration of the above points, and it is possible to obtain a cross-correlation value that is not affected by periodic noise components even if the noise signal detected by the microphone is mixed with periodic noise components. This paper attempts to propose an aircraft noise identification method that can further improve aircraft noise identification performance.

〔問題点を解決するための手段〕[Means for solving problems]

かかる問題点を解決するため本発明においては、第1及
び第2のマイクロホン2A及び2Bから得られる騒音検
出信号SIA及びsl、をサンプリングして時間の経過
に従って順次配列するサンプリングデータでなる第1及
び第2の騒音検出データS4.及び34mを形成する騒
音検出データ形成手段9と、 第1及び第2の騒音検出
データS4A及びS4.についてそれぞれ、順次配列さ
れたサンプリングデータ間の差分を検出して第1及び第
2の差分データS5A及びS5++を形成する白色雑音
化演算手段15と、第1及び第2の差分データS5A及
びS5mの相互相関を演算し、当該相互相関演算信号S
llに生ずるピークに基づいて航空機騒音を識別する相
互相関判定手段2゜とを設けるようにする。
In order to solve this problem, in the present invention, the first and second microphones 2A and 2B have first and second noise detection signals SIA and sl obtained from the first and second microphones 2A and 2B, which are made up of sampling data that are sequentially arranged in accordance with the passage of time. Second noise detection data S4. and a noise detection data forming means 9 forming the first and second noise detection data S4A and S4. a white noise conversion calculation means 15 for detecting the difference between the sequentially arranged sampling data to form first and second difference data S5A and S5++; Calculate the cross-correlation and obtain the cross-correlation calculation signal S
A cross-correlation determining means 2° for identifying aircraft noise based on peaks occurring at ll is provided.

〔作用〕[Effect]

相互相関判定手段によって相互相関演算をするのに先立
って、騒音検出データS4A及びS4゜を白色雑音化演
算手段15によって差分演算処理する。
Prior to the cross-correlation calculation performed by the cross-correlation determining means, the noise detection data S4A and S4° are subjected to differential calculation processing by the white noise conversion calculation means 15.

その結果騒音検出データS4A及び34mに含まれてい
る周期性騒音成分を除去した差分データS5A及び55
11を得ることができ、これを相互相関演算処理するこ
とにより、周期性をもちがっ騒音レベルが高い騒音成分
に基づいて誤判定するおそれを有効に回避し得る航空機
騒音識別装置を実現し得る。
As a result, the difference data S5A and 55 from which periodic noise components included in the noise detection data S4A and 34m are removed
11 can be obtained, and by performing cross-correlation calculation processing on this, it is possible to realize an aircraft noise identification device that can effectively avoid the possibility of misjudgment based on noise components with different periodicity and a high noise level.

〔実施例〕〔Example〕

以下図面について、本発明の一実施例を詳述する。 An embodiment of the present invention will be described in detail below with reference to the drawings.

(1)第1実施例 第1図においてlは全体として航空機騒音識別装置を示
し、上側及び下側マイクロホン2A及び2Bから得られ
る騒音検出信号Sla及びSImを増幅回路3A及び3
Bを介してアナログ/ディジタル変換回路4A及び4B
に入力される。
(1) First Embodiment In FIG. 1, l indicates an aircraft noise identification device as a whole, and noise detection signals Sla and SIm obtained from upper and lower microphones 2A and 2B are transmitted to amplifier circuits 3A and 3.
Analog/digital conversion circuits 4A and 4B via B
is input.

この実施例の場合上側及び下側マイクロホン2A及び2
Bは、第2図に示すように、地表5上に鉛直方向に立接
されたスタンド6に、地表5からそれぞれ高さHA及び
H3の位置にアーム7A及び7Bによって互いに同一の
鉛直線Lv上に配列するように取り付けられ、かくして
航空機8の騒音が航空機8と上側及び下側マイクロホン
2A及び2Bとの間を結ぶ行路WA及びWBをそれぞれ
通って上側及び下側マイクロホン2A及び2Bに到来し
たとき、行路WA及びWBの行路差(従って上側及び下
側マイクロホン2A及び2B位置における行路WA及び
WBの仰角θ、及びθ8の差)に相当する時間遅れをも
ちながら時間の経過に従って変化する騒音検出信号S1
m及びS1mを騒音検出データ形成手段9に取り込むこ
とができるようになされている。
In this embodiment, upper and lower microphones 2A and 2
As shown in FIG. 2, B is placed on the same vertical line Lv by arms 7A and 7B at heights HA and H3 from the ground surface 5, respectively, on a stand 6 that stands vertically on the ground surface 5. Thus, when the noise from the aircraft 8 reaches the upper and lower microphones 2A and 2B through the paths WA and WB connecting the aircraft 8 and the upper and lower microphones 2A and 2B, respectively. , a noise detection signal that changes over time with a time delay corresponding to the path difference between the paths WA and WB (therefore, the difference in elevation angles θ and θ8 of the paths WA and WB at the upper and lower microphone positions 2A and 2B). S1
m and S1m can be taken into the noise detection data forming means 9.

アナログ/ディジタル変換回路4A及び4Bはシステム
コントローラ11から与えられるサンプリング信号S3
によって、そのサンプリング周期で増幅回路3A及び3
Bの出力S2.及び32mの瞬時値をサンプリングして
騒音検出データS4A及びS4.に変換した後、騒音検
出データ形成手段9の出力として差分検出回路12A及
び12Bに供給する。
The analog/digital conversion circuits 4A and 4B receive the sampling signal S3 from the system controller 11.
, the amplifier circuits 3A and 3 at that sampling period.
B's output S2. and 32m instantaneous values were sampled to obtain noise detection data S4A and S4. After converting the data into , it is supplied as the output of the noise detection data forming means 9 to the difference detection circuits 12A and 12B.

差分検出回路12A(及び12B)は直前データラッチ
回路13A(及び13B)と、減算回路14A(及び1
4B)とを有し、サンプリング信号S3が与えられるご
とにアナログ/ディジタル変換回路4A(及び4B)か
ら送出される騒音検出データ54m  (及びS4.)
を減算回路14A(及び14B)に受けると共に、直前
データラッチ回路13A(及び13B)に直前データS
4^8(及び3411X)としてラッチし、かくして白
色雑音化演算手段15を形成する。
The difference detection circuit 12A (and 12B) is connected to the immediately preceding data latch circuit 13A (and 13B) and the subtraction circuit 14A (and 1
4B) and noise detection data 54m (and S4.) sent from the analog/digital conversion circuit 4A (and 4B) every time the sampling signal S3 is given.
is received by the subtraction circuit 14A (and 14B), and the immediately preceding data S is received by the immediately preceding data latch circuit 13A (and 13B).
4^8 (and 3411X), thus forming the white noise conversion calculation means 15.

減算回路14A(及び14B)は、直前データラッチ回
路13A(及び13B)の直前データS4^X(及び3
4ax)を減算入力として受けて、第3図(A)(及び
第3図(B))に示すように、現在時点1fiの騒音検
出データ34A   (及び54R)から直前データ5
4ax(及び341X)を減算して得られる差分データ
S5A (及びS5.)を相互相関判定回路20に供給
する。
The subtraction circuit 14A (and 14B) inputs the immediately preceding data S4^X (and 3
4ax) as a subtraction input, and as shown in FIG. 3(A) (and FIG. 3(B)), the immediately preceding data 5 is obtained from the noise detection data 34A (and 54R) at the current time point 1fi.
The difference data S5A (and S5.) obtained by subtracting 4ax (and 341X) is supplied to the cross-correlation determination circuit 20.

かくして差分検出回路12A(及び12B)は、順次続
くサンプリング時点・・・・・・ta−1、tn、t7
゜、・・・・・・における騒音検出データ34m  (
及び54s)の変化分を1次差分演算によって抽出した
騒音変化分情報を差分データ35m   (及びS5、
)として相互相関判定回路20に供給することになる。
Thus, the difference detection circuit 12A (and 12B) sequentially detects the successive sampling points... ta-1, tn, t7.
Noise detection data 34m at ゜,...
The difference data 35m (and S5,
) is supplied to the cross-correlation determining circuit 20.

相互相関判定回路20は、相互相関演算回路部21及び
判定a踏部22でなる。
The cross-correlation determination circuit 20 includes a cross-correlation calculation circuit section 21 and a determination a-step section 22.

相互相関演算回路部21は、差分データS5A及びS5
mを順次取り込む時間データメモリ25A及び25Bを
有し、かくして時間データメモリ25A及び25Bに保
持した差分データ35m及びS5.をシステムコントロ
ーラ11から与えられる制御信号Sllの制御の下に相
関データ512A及び512Bとして順次読み出して掛
算回路26において乗算した後その掛算データ313を
平均値回路27において平均値演算し、かくして相互相
関値信号514を相互相関演算回路部21の出力として
送出する。
The cross-correlation calculation circuit section 21 calculates the difference data S5A and S5.
It has time data memories 25A and 25B which sequentially take in data 35m and S5. are sequentially read out as correlation data 512A and 512B under the control of the control signal Sll given from the system controller 11, multiplied in the multiplication circuit 26, and the average value of the multiplied data 313 is calculated in the average value circuit 27, thus obtaining the cross-correlation value. A signal 514 is sent out as the output of the cross-correlation calculation circuit section 21.

この実施例の場合、上側時間データメモリ25Aは第4
図(A)に示すように、所定数例えば1064のサンプ
リング時点において得られる差分データ35.を順次記
憶するシフトレジスタで構成され、 例えば1000個
のデータを記憶する同時刻デ−タエリアMs1.Aを中
央部に有すると共に、その前の時刻について例えば32
個分のデータを記憶する前時刻データエリアM、。、と
、後の時刻について32個分のデータを記憶する後時刻
データエリアMAF?Aとを有する。
In this embodiment, the upper time data memory 25A is the fourth
As shown in Figure (A), difference data 35. For example, a simultaneous data area Ms1.Ms1. A in the center and the previous time, for example, 32
a previous time data area M for storing data for each item; , and a later time data area MAF? that stores 32 pieces of data regarding later times. It has A.

下側時間データメモリ25Bも同様にして、同時刻デー
タエリアMs、14.と、前時刻データエリアM@E□
及び後時刻データエリアM @ y y @とを有する
Similarly, the lower time data memory 25B also stores the same time data areas Ms, 14. and previous time data area M@E□
and a later time data area M @ y y @.

上側時間データメモリ25A及び下側時間データメモリ
25Bの各データは、システムコントローラ11の制御
信号Sllによって互いに同期関係を維持しながら次の
ような手順で相互相関データ512A及び312Bとし
て読み出されて掛算回路26に出力される。
Each data in the upper time data memory 25A and the lower time data memory 25B is read out as cross-correlation data 512A and 312B and multiplied by the following procedure while maintaining a synchronized relationship with each other by the control signal Sll of the system controller 11. It is output to the circuit 26.

すなわち下側時間データメモリ25Bの同時刻データエ
リアMs1Mmに記憶されている第33、第34・・・
・・・第1032番目のデータが前時刻データに対応し
て繰り返し32回読み出され、続いて同時刻データに対
応して1回読み出され、続いて後時刻データに対応して
繰り返し32回読み出される。ここで当該同時刻データ
エリアM□。が1回読み出されている間に、上側時間デ
ータメモリ25Aから1000個分のデータが1データ
ずつシフトされながら順次読み出される。
That is, the 33rd, 34th, . . . stored in the same time data area Ms1Mm of the lower time data memory 25B.
...The 1032nd data is read out repeatedly 32 times in response to the previous time data, then read out once in response to the same time data, and then repeatedly read out 32 times in response to the later time data. Read out. Here, the same time data area M□. is read out once, 1000 pieces of data are sequentially read out from the upper time data memory 25A while being shifted one data at a time.

例えば前時刻データエリアMl!、、の第1〜第100
0番目のデータが順次読み出されるタイミングで下側時
間データメモリ25Bの第33〜第1032番目のデー
タが順次読み出され、掛算回路26において互いに掛算
された後掛算データS13として平均値回路27に出力
される。
For example, previous time data area Ml! , 1st to 100th
The 33rd to 1032nd data of the lower time data memory 25B are sequentially read out at the timing when the 0th data is sequentially read out, multiplied by each other in the multiplication circuit 26, and then outputted to the average value circuit 27 as multiplication data S13. be done.

平均値回路27は順次出力されて来る1000個の掛算
データS13から平均値を求めて第1番目のサンプリン
グ点の相関度を表す相互相関値信号S14として出力す
る。
The average value circuit 27 calculates an average value from the 1000 multiplication data S13 that are sequentially output, and outputs it as a cross-correlation value signal S14 representing the correlation degree of the first sampling point.

このようにして前時刻の第1番目のサンプリング点につ
いての相互相関演算が終了すると、続いて上側時間デー
タメモリ25Aから第2〜第1001番目のデータ(す
なわちlデー5分だけシフトしたデータ)が読み出され
て下側時間データメモリ25Bから繰り返して読み出さ
れる第33〜第1032番目のデータと順次掛算されて
掛算データS13として平均値回路27に出力される。
When the cross-correlation calculation for the first sampling point of the previous time is completed in this way, the 2nd to 1001st data (that is, the data shifted by 1 day and 5 minutes) are subsequently transferred from the upper time data memory 25A. The data is sequentially multiplied by the 33rd to 1032nd data read out and repeatedly read out from the lower time data memory 25B, and outputted to the average value circuit 27 as multiplied data S13.

このとき平均値回路27は順次出力されて来る1000
個の掛算データS13から平均値を求めて第2番目のサ
ンプリング点の相関度を表す相互相関値信号S14とし
て送出する。
At this time, the average value circuit 27 sequentially outputs 1000
An average value is calculated from the multiplied data S13 and sent as a cross-correlation value signal S14 representing the correlation degree of the second sampling point.

以下同様にして、上側時間データメモリ25Aの第32
〜第1031番目までのデータに基づいて第32番目の
サンプリング点についての相互相関値信号S14を送出
し、これにより前時刻データエリアM0、に関連する相
互相関演算を終了する。
Thereafter, in the same manner, the 32nd
The cross-correlation value signal S14 for the 32nd sampling point is sent out based on the data from the 1031st to the 1031st, thereby ending the cross-correlation calculation related to the previous time data area M0.

その後上側時間データメモリ25Aは、第4図(B)に
示すように、同時刻データエリアMs、M。
Thereafter, the upper time data memory 25A stores the same time data areas Ms and M, as shown in FIG. 4(B).

に記憶されている第33〜第1032番目のデータを読
み出して掛算回路26に供給する。このとき平均値回路
27は1000個の掛算データ313から平均値を求め
て第33番目のサンプリング点の相関度を表す相互相関
値信号S14を送出し、これにより同時刻データエリア
M84.に関連する相互相関演算を終了する。
The 33rd to 1032nd data stored in are read out and supplied to the multiplication circuit 26. At this time, the average value circuit 27 calculates the average value from the 1000 multiplication data 313 and sends out a cross-correlation value signal S14 representing the degree of correlation of the 33rd sampling point, thereby causing the same-time data area M84. Finish the cross-correlation operation related to .

その後続いて上側時間データメモリ25Aは、第4図(
C)に示すように、後時刻データエリアM A門Aに関
連して、前時刻データエリアM@、、Aについて上述し
たと同様にして、第34、第35・・・・・・第66番
目のサンプリング点の相関度を表す相互相関値信号S1
4を送出する。
Subsequently, the upper time data memory 25A is stored in the upper time data memory 25A as shown in FIG.
As shown in C), in relation to the later time data area M A gate A, the 34th, 35th, 66th, etc. Cross-correlation value signal S1 representing the correlation degree of the th sampling point
Send 4.

このようにして平均値回路27から得られる相互相関値
信号514は、第5図に示すように65個のサンプリン
グ点についての相関度を表すが、これらのデータのうち
、第33番目のデータは上側マイクロホン2A及び下側
マイクロホン2Bに航空機騒音が同時に到達した状態に
おける相関度を表すデータでなり、このデータにサンプ
リング点番号「0」を割り当てる。
The cross-correlation value signal 514 obtained from the average value circuit 27 in this way represents the correlation degree for 65 sampling points as shown in FIG. 5, but among these data, the 33rd data is This data represents the degree of correlation in a state where aircraft noise reaches the upper microphone 2A and the lower microphone 2B at the same time, and a sampling point number "0" is assigned to this data.

また第32、第31・・・・・・第1番目のデータは航
空機騒音が上側マイクロホン2Aに到達する時点が下側
マイクロホン2Bに到達する時点より前の時刻であるこ
とを表すデータでなり、このデータにサンプリング点番
号「+1」、「+2」・・・・・・「+32」を割り当
てる。
Further, the 32nd, 31st, etc., first data is data representing that the time when the aircraft noise reaches the upper microphone 2A is before the time when the aircraft noise reaches the lower microphone 2B, Sampling point numbers "+1", "+2", . . . "+32" are assigned to this data.

さらに第34、第35・・・・・・第65番目のデータ
は航空機騒音が上側マイクロホン2Aに到達する時点が
下側マイクロホン2Bに到達する時点より後の時刻であ
ることを表すデータでなり、このデータにサンプリング
点番号「−1」、「−2」・・・・・・「−32」を割
り当てる。
Furthermore, the 34th, 35th, . . . , 65th data are data representing that the time point at which the aircraft noise reaches the upper microphone 2A is later than the time point at which it reaches the lower microphone 2B, Sampling point numbers "-1", "-2", . . . "-32" are assigned to this data.

かくして第5図に示すように、サンプリング番号r−3
2」・・・・・・「−1」、「0」、「+1」・・・・
・・r+324のデータをその順序で配列するような相
互相関値信号514を形成してこれをモニタグラフに表
示すれば、相関度ピークがサンプリング点番号「0」か
ら右側のサンプリング点番号「+IJ、「+2」・・・
・・・r+32Jに移動して行くと、このことは航空機
が上側及び下側マイクロホン2A及び2Bから遠(離れ
た位置に有る状態(従って行路WA及びWBの長さがほ
ぼ等しい(仰角θ1及びθ、もほぼ等しい)状態)から
真上の位置に近づいて来つつある状態(従って行路WA
の長さが行路WBの長さより短かくなって行く状態)を
表すことになる。
Thus, as shown in FIG.
2"..."-1", "0", "+1"...
...If we form a cross-correlation value signal 514 in which r+324 data are arranged in that order and display this on a monitor graph, the correlation degree peak will change from the sampling point number "0" to the right sampling point number "+IJ, "+2"...
... As we move to r+32J, this means that the aircraft is far from the upper and lower microphones 2A and 2B (thus, the lengths of paths WA and WB are approximately equal (elevation angles θ1 and θ , are almost equal)) to a state that is approaching the position directly above (therefore, the path WA
This represents a state in which the length of the route WB becomes shorter than the length of the route WB.

これに対して相関度のピークがサンプリング点番号「0
」から左側のサンプリング点番号「−1」、「−2」・
・・・・・r−32Jに移動して行くとこのことは、航
空機が遠く離れた位置から上側及び下側マイクロホン2
A及び2Bの真上の位置に近づいて来た状態において航
空機騒音が地表5(第2図)で反射して下側マイクロホ
ン2B側から到来したことを表すことになる。
On the other hand, the peak of the correlation degree is the sampling point number “0”.
” to the left sampling point number “-1”, “-2”.
...Moving to r-32J, this means that the upper and lower microphones 2
This indicates that the aircraft noise is reflected from the ground surface 5 (FIG. 2) and arrives from the lower microphone 2B side when it approaches the position directly above A and 2B.

相互相関値信号314は、判定回路部22のピーク検出
回路28に与えられ、相互相関値信号S14が最大値に
なったときピーク検出回路28からピーク検出信号S1
5をピークアドレスメモリ29に送出する。
The cross-correlation value signal 314 is given to the peak detection circuit 28 of the determination circuit section 22, and when the cross-correlation value signal S14 reaches the maximum value, the peak detection circuit 28 outputs the peak detection signal S1.
5 to the peak address memory 29.

ピークアドレスメモリ29はピーク検出信号S15が得
られるごとにそのタイミングにおけるサンプリング点番
号を記憶更新して行き、かくしてサンプリング点番号「
+32」からr−32J  (第5図)に至るまでの相
互相関値信号S14について相関度が最大値のサンプリ
ング点番号JDGをピークアドレス検出信号S16とし
て判定回路30に送出する。
The peak address memory 29 stores and updates the sampling point number at that timing every time the peak detection signal S15 is obtained, and thus the sampling point number "
For the cross-correlation value signal S14 ranging from "+32" to r-32J (FIG. 5), the sampling point number JDG having the maximum correlation degree is sent to the determination circuit 30 as the peak address detection signal S16.

判定回路30はピーク値のサンプリング点番号JDGに
基づいて、航空機の飛行位置を表す航空機騒音識別信号
5174相互相関判定回路20の出力として送出する。
Based on the sampling point number JDG of the peak value, the determination circuit 30 sends out an aircraft noise identification signal 5174 representing the flight position of the aircraft as an output of the cross-correlation determination circuit 20.

この実施例の場合相互相関演算回路部21から送出され
る相互相関値信号S14は、ディジタル/アナログ変換
回路31においてアナログ信号に変換さた後モニタ信号
818として相互相関判定回路20から送出するように
なされ、このモニタ信号818を例えばオシュロスコー
プでなるモニタに与えることにより第5図に示すような
相関度曲線を表示させ、これによりそのピークの発生位
置によって航空機の飛行位置をオペレータが目視確認で
きるようになされている。
In this embodiment, the cross-correlation value signal S14 sent out from the cross-correlation calculation circuit section 21 is converted into an analog signal in the digital/analog conversion circuit 31, and then sent out from the cross-correlation determination circuit 20 as a monitor signal 818. By applying this monitor signal 818 to a monitor such as an oscilloscope, a correlation curve as shown in FIG. is being done.

以上の構成において、上側マイクロホン2A及び下側マ
イクロホン2Bから得られる騒音検出信号SIA及びS
1mをアナログ/ディジタル変換回路4A及び4Bにお
いてサンプリングして騒音検出データS4.及びS41
を得たとき、この騒音検出データS4A及びS41を差
分検出回路12A及び12Bにおいて差分データS5A
及びS5、に変換した後相互相関判定回路20において
相互相関の演算及び当該相互相関値信号S14のピーク
位置に対応するサンプリング点番号を表す航空機騒音識
別信号317を得ることができる。
In the above configuration, the noise detection signals SIA and S obtained from the upper microphone 2A and the lower microphone 2B are
1m is sampled in analog/digital conversion circuits 4A and 4B and noise detection data S4. and S41
When the noise detection data S4A and S41 are obtained, the difference detection circuits 12A and 12B convert the noise detection data S4A and S41 into difference data S5A.
and S5, the cross-correlation determination circuit 20 calculates the cross-correlation and obtains the aircraft noise identification signal 317 representing the sampling point number corresponding to the peak position of the cross-correlation value signal S14.

かくするにつき、差分検出回路12A及び12Bは、各
サンプリング時点における騒音検出データS4A及び3
4.の変化分を表す情報だけを差分データS5A及びs
5mとして取り出すことができることにより、たとえ騒
音検出信号SIA及びS1mに周期性がありかつ騒音レ
ベルが大きい騒音成分が含まれている場合にも、当該騒
音成分を除外できることによりほぼ−様な周波数スペク
トラムを有する白色雑音としての特性をもたせることが
でき(これを白色雑音化と呼ぶ)、かくして安定な信号
レベルを有する相互相関値信号S14を相互相関判定回
路20に供給できる。
Therefore, the difference detection circuits 12A and 12B detect the noise detection data S4A and 3 at each sampling point.
4. Only the information representing the change in difference data S5A and s
5m, even if the noise detection signals SIA and S1m contain noise components with periodicity and a large noise level, the noise components can be excluded and a nearly -like frequency spectrum can be obtained. In this way, the cross-correlation value signal S14 having a stable signal level can be supplied to the cross-correlation determination circuit 20.

因に例えばプロペラを持つ航空機の騒音として100〜
200 (Hz)程度の周期性騒音成分が混在している
場合には、 当該周期性騒音成分の周期が差分検出回路
12A及び12Bの差分演算周期より十分に長い周期を
もっているので、差分データS5A及びS5mに含まれ
る当該航空機騒音成分の信号レベルがほぼ等しくなるた
めに減算回路14A及び14Bにおける減算動作によっ
て相殺され、結局差分データ35A及びs5gには含ま
せないようにし得る。
For example, the noise of an aircraft with a propeller is 100~
If a periodic noise component of about 200 (Hz) is mixed, the period of the periodic noise component is sufficiently longer than the difference calculation period of the difference detection circuits 12A and 12B, so the difference data S5A and Since the signal levels of the aircraft noise components included in S5m are almost equal, they are canceled out by the subtraction operations in the subtraction circuits 14A and 14B, and can be made not to be included in the difference data 35A and s5g.

また同様にして航空機騒音以外の市街騒音例えば250
 (Hz) 〜1 (k)lx)程度のサイレン、クラ
クション等の騒音成分についても同様にして、当該騒音
成分の周期が差分検出回路12A及び12Bにおける差
分演算の周期より十分に長いことによりその影響を除去
することができる。
Similarly, city noise other than aircraft noise, for example, 250
Similarly, regarding noise components such as sirens and horns with a frequency of approximately (Hz) to 1 (k)lx), the influence of noise components is affected by the fact that the period of the noise component is sufficiently longer than the period of the difference calculation in the difference detection circuits 12A and 12B. can be removed.

従って上述の構成によれば、 たとえ上側及び下側マイ
クロホン2A及び2Bの騒音検出信号SIA及びSim
に周期性をもちかつ大きい騒音レベルの騒音成分が混入
しても、実質上その影響を受けずに航空機の飛行位置を
確実に検出することができる。
Therefore, according to the above configuration, even if the noise detection signals SIA and Sim of the upper and lower microphones 2A and 2B
Even if a noise component having periodicity and a high noise level is mixed into the noise component, the flight position of the aircraft can be reliably detected without being substantially affected by the noise component.

因に100〜200 (Hz)程度の周期性がある騒音
成分を有する航空機騒音について、差分検出回路12A
及び12Bの差分演算処理を実行せずに相互相関判定回
路20における相互相関の判定をすれば、第6図におい
て曲線314Xに示すように、航空機騒音に含まれてい
る周期性騒音成分に基づいて相関度が当該騒音成分の周
期に対応するように大きく変動することにより、ピーク
が複数個発生したり、ピーク波形が曖昧になったりする
のに対して、第1図の構成によれば、第5図の最大ピー
クサンプリング点番号JDGに示すように、航空機の飛
行位置を表す情報として、ピーク波形が急峻な相互相関
値信号を得ることができる。
Incidentally, regarding aircraft noise having noise components with periodicity of about 100 to 200 (Hz), the difference detection circuit 12A
If the cross-correlation determination circuit 20 determines the cross-correlation without executing the difference calculation processing of If the degree of correlation changes greatly to correspond to the period of the noise component, multiple peaks may occur or the peak waveform may become ambiguous. As shown by the maximum peak sampling point number JDG in FIG. 5, a cross-correlation value signal with a steep peak waveform can be obtained as information representing the flight position of the aircraft.

実験によれば、第7図に示すように、プロペラを存する
航空機が頭上を通過したときの航空機騒音について、そ
の相互相関値信号S14を測定したところ第7図(B)
に示すように、航空機の飛行位置を表すピーク波形がス
タート時点においてサンプリング点番号「0」近傍の位
置にあったのに対して、3.5秒後に航空機が頭上を通
過する際にはサンプリング点番号r+32Jの方向に大
きく移動して行き、その後7.5秒後までの間に航空機
が遠ざかるに従って再びサンプリング点番号「0」の方
向に戻って行くような波形変化を呈する測定情報を得る
ことができた。
According to experiments, as shown in Fig. 7, the cross-correlation value signal S14 of aircraft noise when an aircraft with a propeller passes overhead is measured, and the result is shown in Fig. 7 (B).
As shown in , the peak waveform representing the flight position of the aircraft was near sampling point number "0" at the start, but when the aircraft passed overhead 3.5 seconds later, the peak waveform representing the flight position of the aircraft was located near sampling point number "0". It is possible to obtain measurement information exhibiting a waveform change in which the aircraft moves largely in the direction of number r+32J and then returns to the direction of sampling point number "0" as the aircraft moves away from it until 7.5 seconds later. did it.

これに対して差分検出回路12A及び12Bを省略した
構成によって同じような測定をしたところ、第7図(A
)に示すように、相互相関値信号514Xに、スタート
時から7゜5秒経過するまでの間にプロペラの騒音に基
づく周期性のある騒音成分が相互相関値信号波形514
X上に集畳することにより、航空機の位置を表すはずの
ピーク値が急峻には出て来ないのみならず、航空機が頭
上を通過する時点(すなわち3.5秒後〜4.5秒後)
において2つのピーク波形が発生するような結果を得た
On the other hand, when similar measurements were made using a configuration in which the difference detection circuits 12A and 12B were omitted, the results were as shown in FIG.
), a periodic noise component based on propeller noise appears in the cross-correlation value signal 514
By converging on )
The result was that two peak waveforms were generated.

かかる実験結果から差分検出回路12A及び12Bが航
空機騒音に含まれている周期性騒音成分の影響を受けな
いような相互相関値信号314を形成し得ることを確認
し得た。
From these experimental results, it has been confirmed that the difference detection circuits 12A and 12B can form a cross-correlation value signal 314 that is not affected by periodic noise components included in aircraft noise.

〔2〕他の実施例 (1)上述の実施例においては白色雑音化演算手段15
として1次差分演算回路を用いた実施例について述べた
がこれに代え、例えばn次差分演算回路、線形予測フィ
ルタ、偏自己相関演算回路及びこれらの組合せなど種々
の構成のものを適用し得、要は白色雑音化処理手段を適
用すれば良い。
[2] Other embodiments (1) In the above embodiment, the white noise conversion calculation means 15
Although an embodiment using a first-order difference calculation circuit has been described, instead of this, various configurations such as an n-order difference calculation circuit, a linear prediction filter, a partial autocorrelation calculation circuit, and a combination thereof may be applied. In short, it is sufficient to apply a white noise processing means.

(2)上述の実施例においては2つのマイクロホン2A
及び2Bを鉛直方向に配列した場合について述べたが、
配列方向はこれに限らず例えば水平方向など種々の方向
に配列しても上述の場合と同様の効果を得ることができ
る。
(2) In the above embodiment, two microphones 2A
We have described the case where 2B and 2B are arranged vertically,
The arrangement direction is not limited to this, and the same effect as described above can be obtained by arranging in various directions, such as the horizontal direction.

(3)上述の実施例においては本発明による航空機騒音
識別装置をハード的回路手段を用いて構成した場合につ
いて述べたがこれに代え、ソフト的演算手段を用いて構
成するようにしても上述の場合と同様の効果を得ること
ができる。
(3) In the above-described embodiments, the case where the aircraft noise identification device according to the present invention is configured using hardware circuit means has been described, but instead of this, it may be configured using software calculation means as well. You can get the same effect as in the case.

〔発明の効果〕〔Effect of the invention〕

上述のように本発明によれば、一対のマイクロホンから
得た騒音ネ★出信号をディジタルデータに変換した後白
色雑音化演算処理をすると共に相互相関処理を実行する
ようにしたことにより、航空機騒音に周期性をもつ強い
騒音レベルの騒音成分が混入していてもその影響を受け
ることなく確実に航空機騒音を識別することができる。
As described above, according to the present invention, the noise output signals obtained from a pair of microphones are converted into digital data, and then subjected to white noise calculation processing and cross-correlation processing. Even if a noise component with a strong noise level having periodicity is mixed in, aircraft noise can be reliably identified without being affected by the noise component.

段、20・・・・・・相互相関判定回路、21・・・・
・・相互相関演算回路部、22・・・・・・判定回路部
Stage, 20... Cross correlation determination circuit, 21...
. . . Cross-correlation calculation circuit section, 22 . . . Determination circuit section.

【図面の簡単な説明】[Brief explanation of drawings]

第1図は本発明による航空機騒音識別装置の一実施例を
示すブロック図、第2図はそのマイクロホン装置の構成
を示す路線的側面図、第3図は差分演算の説明に供する
曲線図、第4図は相互相関演算の説明に供する路線図、
第5図及び第6図は相互相関演算結果の説明に供する信
号波形図、第7図は実験結果を示す信号波形図である。
FIG. 1 is a block diagram showing an embodiment of the aircraft noise identification device according to the present invention, FIG. 2 is a line side view showing the configuration of the microphone device, FIG. 3 is a curve diagram for explaining the difference calculation, and FIG. Figure 4 is a route map for explaining the cross-correlation calculation.
5 and 6 are signal waveform diagrams for explaining the cross-correlation calculation results, and FIG. 7 is a signal waveform diagram showing the experimental results.

Claims (1)

【特許請求の範囲】 第1及び第2のマイクロホンから得られる騒音検出信号
をサンプリングして時間の経過に従つて順次配列するサ
ンプリングデータでなる第1及び第2の騒音検出データ
を形成する騒音検出データ形成手段と、 上記第1及び第2の騒音検出データについてそれぞれ、
順次配列されたサンプリングデータ間の差分を検出して
第1及び第2の差分データを形成する白色雑音化演算手
段と、 上記第1及び第2の差分データの相互相関を演算し、当
該相互相関演算信号に生ずるピークに基づいて航空機騒
音を識別する相互相関判定手段とを具えることを特徴と
する航空機騒音識別装置。
[Scope of Claims] Noise detection that forms first and second noise detection data consisting of sampling data that samples noise detection signals obtained from first and second microphones and sequentially arranges them over time. a data forming means, and each of the first and second noise detection data;
white noise conversion calculation means for detecting the difference between the sequentially arranged sampling data to form first and second difference data; and calculating a cross-correlation between the first and second difference data; 1. An aircraft noise identification device comprising: cross-correlation determination means for identifying aircraft noise based on peaks occurring in a calculated signal.
JP24204488A 1988-09-27 1988-09-27 Device for discriminating noise of aircraft Granted JPH0288928A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP24204488A JPH0288928A (en) 1988-09-27 1988-09-27 Device for discriminating noise of aircraft

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP24204488A JPH0288928A (en) 1988-09-27 1988-09-27 Device for discriminating noise of aircraft

Publications (2)

Publication Number Publication Date
JPH0288928A true JPH0288928A (en) 1990-03-29
JPH0515971B2 JPH0515971B2 (en) 1993-03-03

Family

ID=17083434

Family Applications (1)

Application Number Title Priority Date Filing Date
JP24204488A Granted JPH0288928A (en) 1988-09-27 1988-09-27 Device for discriminating noise of aircraft

Country Status (1)

Country Link
JP (1) JPH0288928A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011237323A (en) * 2010-05-12 2011-11-24 Rion Co Ltd Noise measuring device

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011237323A (en) * 2010-05-12 2011-11-24 Rion Co Ltd Noise measuring device

Also Published As

Publication number Publication date
JPH0515971B2 (en) 1993-03-03

Similar Documents

Publication Publication Date Title
US10245724B2 (en) System, method and product for utilizing prediction models of an environment
EA007881B1 (en) Estimating the time of arrival of a seismic wave
US20190162812A1 (en) System and method for determining the source location of a firearm discharge
JP2005354676A (en) Signal processing circuit and method or time delay detection apparatus and object position specification apparatus using the method
WO2008048365A3 (en) Process and system of energy signal detection
US10969464B2 (en) Method for registering presence of a stationary object exhibiting motion at a location in a scene monitored by a radar detector
CN114492629A (en) Abnormality detection method, abnormality detection device, electronic apparatus, and storage medium
JPH0288928A (en) Device for discriminating noise of aircraft
GB2530374A (en) Device and method for sound-based environment detection
KR20170129361A (en) System and method for tracking target using asynchronous sensors
KR102146893B1 (en) Vehicle opening and closing system using sound wave communication and control method thereof
KR20190079183A (en) Apparatus for recognizing pedestrian and vehicle using doppler information and method thereof
US20170343644A1 (en) Detection of acoustic events
KR940004332A (en) Circuit test method and delay defect detector
Tejedor et al. Towards detection of pipeline integrity threats using a SmarT fiber-OPtic surveillance system: PIT-STOP project blind field test results
JPH05134032A (en) Method and circuit apparatus for identifying signal accompanying noise
CN106100771B (en) A kind of two-way time delay detecting method and device
JPH05346464A (en) Uniting equipment of target
JP2002054986A (en) Sound source contribution analyzing method and device having background noise separating function
CN115309618B (en) Input and output delay detection method and system based on RT-Linux system
JP2590999B2 (en) Image signal motion vector detector
SU1122641A1 (en) Device for determining properties of random process
JPS6035284A (en) Varying time difference detecting system
SU696483A2 (en) Adaptive time discretization device
JPH01172789A (en) Appearance detector of underwater acoustic signal

Legal Events

Date Code Title Description
R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080303

Year of fee payment: 15

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090303

Year of fee payment: 16

EXPY Cancellation because of completion of term
FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090303

Year of fee payment: 16