JPH0515971B2 - - Google Patents

Info

Publication number
JPH0515971B2
JPH0515971B2 JP24204488A JP24204488A JPH0515971B2 JP H0515971 B2 JPH0515971 B2 JP H0515971B2 JP 24204488 A JP24204488 A JP 24204488A JP 24204488 A JP24204488 A JP 24204488A JP H0515971 B2 JPH0515971 B2 JP H0515971B2
Authority
JP
Japan
Prior art keywords
noise
data
cross
correlation
aircraft
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.)
Expired - Lifetime
Application number
JP24204488A
Other languages
Japanese (ja)
Other versions
JPH0288928A (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

Description

【発明の詳細な説明】 〔産業上の利用分野〕 本発明は航空機騒音識別装置に関し、特に航空
機騒音に含まれる雑音成分の影響を軽減できるよ
うにしたものである。
DETAILED DESCRIPTION OF THE INVENTION [Industrial Field of Application] The present invention relates to an aircraft noise identification device, and in particular is capable of reducing the influence of noise components contained in aircraft noise.

〔発明の概要〕[Summary of the invention]

本発明は、相互相関演算方式の航空機騒音識別
装置において、騒音検出データを差分演算処理し
た後相互相関演算を実行するようにしたことによ
り、航空機騒音を表すピークを一段と急峻にし得
る。
According to the present invention, in an aircraft noise identification device using a cross-correlation calculation method, a cross-correlation calculation is executed after noise detection data is subjected to a difference calculation process, thereby making it possible to make the peak representing aircraft noise even steeper.

〔従来の技術〕[Conventional technology]

航空機騒音を評価する基準として、「航空機騒
音に係る環境基準(昭和48年12月27日環境庁告示
154)」があり、当該基準に従つて製造された航空
機騒音自動測定装置が従来から用いられている。
As a standard for evaluating aircraft noise, the ``Environmental Standards for Aircraft Noise'' (December 27, 1971, Environment Agency Notification)
154), and automatic aircraft noise measurement equipment manufactured in accordance with this standard has been used for some time.

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

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

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

因に第2の方法によつて2つの測定信号間の相
互相関値を求める方法は、騒音が白色雑音として
の性質をもつと考えた場合当該騒音のトータル音
量についての相互相関値のピークの位置の時間的
変化が航空機の飛行位置を表している点に着目し
て航空機騒音を識別するようになされている。
Incidentally, the second method for determining the cross-correlation value between two measurement signals is based on the position of the peak 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 temporal changes in the noise represent the flight position of the aircraft.

ところが実際上例えばプロペラを有する航空機
の場合には100〜200〔Hz〕程度の周波数範囲に周
期性が強く騒音レベルが大きい騒音成分を含んで
おり、また市街雑音の中にサイレン、クラクシヨ
ンなどが混入すると、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 [Hz], and street noises include sirens, horns, etc. As a result, 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 It changes based on movement, and at the same time changes with the period of the noise component.
This results in noise being mixed into the aircraft noise identification results.

このような雑音を含んだ識別結果に基づいて航
空機騒音であるか否かを判定することは、実用上
極めて困難であり、改善することが望ましい。
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
By making it possible to obtain cross-correlation values that are unaffected even when periodic noise components are mixed in the noise signal detected by the microphone, aircraft noise identification performance can be further improved. This paper attempts to propose a method for identifying aircraft noise.

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

かかる問題点を解決するため本発明において
は、第1及び第2のマイクロホン2A及び2Bか
ら得られる騒音検出信号S1A及びS1Bをサンプ
リングして時間の経過に従つて順次配列するサン
プリングデータでなる第1及び第2の騒音検出デ
ータS4A及びS4Bを形成する騒音検出データ形
成手段9と、第1及び第2の騒音検出データS4
及びS4Bについてそれぞれ、順次配列されたサ
ンプリングデータ間の差分を検出して第1及び第
2の差分データS5A及びS5Bを形成する白色雑
音化演算手段15と、第1及び第2の差分データ
S5A及びS5Bの相互相関を演算し、当該相互相
関演算信号S11に生ずるピークに基づいて航空
機騒音を識別する相互相関判定手段20とを設け
るようにする。
In order to solve this problem, in the present invention, the noise detection signals S1 A and S1 B obtained from the first and second microphones 2A and 2B are sampled and are made up of sampling data that is sequentially arranged as time passes. Noise detection data forming means 9 for forming first and second noise detection data S4 A and S4 B , and first and second noise detection data S4
A and S4B , respectively, detecting the difference between the sequentially arranged sampling data to form first and second difference data S5A and S5B ; A cross-correlation determination means 20 is provided which calculates the cross-correlation of the difference data S5A and S5B and identifies aircraft noise based on the peak that occurs in the cross-correlation calculation signal S11.

〔作用〕[Effect]

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

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

〔実施例〕〔Example〕

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

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

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

アナログ/デイジタル変換回路4A及び4Bは
システムコントローラ11から与えられるサンプ
リング信号S3によつて、そのサンプリング周期
で増幅回路3A及び3Bの出力S2A及びS2B
瞬時値をサンプリングして騒音検出データS4A
及びS4Bに変換した後、騒音検出データ形成手
段9の出力として差分検出回路12A及び12B
に供給する。
The analog/digital conversion circuits 4A and 4B sample the instantaneous values of the outputs S2 A and S2 B of the amplifier circuits 3A and 3B at the sampling period according to the sampling signal S3 given from the system controller 11, and generate noise detection data S4 A.
and S4B , the difference detection circuits 12A and 12B are output as the output of the noise detection data forming means 9.
supply to.

差分検出回路12A(及び12B)は直前デー
タラツチ回路13A(及び13B)と、減算回路
14A(及び14B)とを有し、サンプリング信
号S3が与えられるごとにアナログ/デイジタル
変換回路4A(及び4B)から送出される騒音検
出データS4A(及びS4B)を減算回路14A
(及び14B)に受けると共に、直前データラツ
チ回路13A(及び13B)に直前データS4AX
(及びS4BX)としてラツチし、かくして白色雑
音化演算手段15を形成する。
The difference detection circuit 12A (and 12B) has an immediately preceding data latch circuit 13A (and 13B) and a subtraction circuit 14A (and 14B). The noise detection data S4 A (and S4 B ) to be sent is subtracted by a subtraction circuit 14A.
(and 14B), and the immediately preceding data latch circuit 13A (and 13B) receives the immediately preceding data S4 AX.
(and S4 BX ), thus forming the white noise conversion calculation means 15.

減算回路14A(及び14B)は、直前データ
ラツチ回路13A(及び13B)の直前データS
AX(及びS4BX)を減算入力として受けて、第
3図A(及び第3図B)に示すように、現在時点
toの騒音検出データS4A(及びS4B)から直前
データS4AX(及びS4BX)を減算して得られる
差分データS5A(及びS5B)を相互相関判定回
路20に供給する。
The subtracting circuit 14A (and 14B) receives the immediately preceding data S of the immediately preceding data latch circuit 13A (and 13B).
4 AX (and S4 BX ) as subtraction input, the current point is calculated as shown in Figure 3A (and Figure 3B).
Difference data S5 A (and S5 B ) obtained by subtracting immediately preceding data S4 AX (and S4 BX ) from noise detection data S4 A (and S4 B ) of t o is supplied to the cross-correlation determination circuit 20 .

かくして差分検出回路12A(及び12B)は、
順次続くサンプリング時点……to-1、to、to+1
…における騒音検出データS4A(及びS4B)の
変化分を1次差分演算によつて抽出した騒音変化
分情報を差分データS5A(及びS5B)として相
互相関判定回路20に供給することになる。
Thus, the difference detection circuit 12A (and 12B)
Sequential sampling points... t o-1 , t o , t o+1 ...
The noise change information extracted from the change in noise detection data S4 A (and S4 B ) in ... by a first-order difference calculation is supplied to the cross-correlation determination circuit 20 as difference data S5 A (and S5 B ). Become.

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

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

この実施例の場合、上側時間データメモリ25
Aは第4図Aに示すように、所定数例えば1064の
サンプリング時点において得られる差分データS
Aを順次記憶するシフトレジスタで構成され、
例えば1000個のデータを記憶する同時刻データエ
リアMSIMAを中央部に有すると共に、その前の時
刻について例えば32個分のデータを記憶する前時
刻データエリアMBEFAと、後の時刻について32個
分のデータを記憶する後時刻データエリアMAFTA
とを有する。
In this embodiment, the upper time data memory 25
As shown in FIG. 4A, A is the difference data S obtained at a predetermined number of sampling points, for example, 1064.
Consists of a shift register that sequentially stores 5 A ,
For example, there is a simultaneous data area M SIMA in the center that stores 1000 pieces of data, and a previous time data area M BEFA that stores, for example, 32 pieces of data for the previous time, and 32 pieces of data for the later time. The rear time data area M AFTA stores the data of
and has.

下側時間データメモリ25Bも同様にして、同
時刻データエリアMSIMBと、前時刻データエリア
MBEFA及び後時刻データエリアMAFTBとを有する。
Similarly, the lower time data memory 25B has a same time data area M SIMB and a previous time data area.
It has M BEFA and a later time data area M AFTB .

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

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

例えば前時刻データエリアMBEFAの第1〜第
1000番目のデータが順次読み出されるタイミング
で下側時間データメモリ25Bの第33〜第1032番
目のデータが順次読み出され、掛算回路26にお
いて互いに掛算された後掛算データS13として
平均値回路27に出力される。
For example, the 1st to 1st of the previous time data area M BEFA
At the timing when the 1000th data is read out in sequence, the 33rd to 1032nd data of the lower time data memory 25B are read out in sequence, 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番目の
サンプリング点の相関度を表す相互相関値信号S
14として出力する。
The average value circuit 27 calculates the average value from the 1000 multiplication data S13 that are sequentially output, and generates a cross-correlation value signal S representing the correlation degree of the first sampling point.
Output as 14.

このようにして前時刻の第1番目のサンプリン
グ点についての相互相関演算が終了すると、続い
て上側時間データメモリ25Aから第2〜第1001
番目のデータ(すなわち1データ分だけシフトし
たデータ)が読み出されて下側時間データメモリ
25Bから繰り返して読み出される第33〜第1032
番目のデータと順次掛算されて掛算データS13
として平均値回路27に出力される。このとき平
均値回路27は順次出力されて来る1000個の掛算
データS13から平均値を求めて第2番目のサン
プリング点の相関度を表す相互相関値信号S14
として送出する。
When the cross-correlation calculation for the first sampling point of the previous time is completed in this way, the second to 1001th sampling points from the upper time data memory 25A are
The 33rd to 1032nd data (that is, data shifted by one data) are read out and repeatedly read out from the lower time data memory 25B.
The multiplied data S13 is sequentially multiplied with the th data.
It is output to the average value circuit 27 as a value. At this time, the average value circuit 27 calculates an average value from the 1000 multiplication data S13 that are sequentially output, and generates a cross-correlation value signal S14 representing the correlation degree of the second sampling point.
Send as.

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

その後上側時間データメモリ25Aは、第4図
Bに示すように、同時刻データエリアASIMAに記
憶されている第33〜第1032番目のデータを読み出
して掛算回路26に供給する。このとき平均値回
路27は1000個の掛算データS13から平均値を
求めて第33番目のサンプリング点の相関度を表す
相互相関値信号S14を送出し、これにより同時
刻データエリアMSIMAに関連する相互相関演算を
終了する。
Thereafter, the upper time data memory 25A reads out the 33rd to 1032nd data stored in the simultaneous data area A SIMA and supplies it to the multiplication circuit 26, as shown in FIG. 4B. At this time, the average value circuit 27 calculates the average value from the 1000 multiplication data S13 and sends out a cross-correlation value signal S14 representing the degree of correlation of the 33rd sampling point, and thereby Finish the cross-correlation calculation.

その後続いて上側時間データメモリ25Aは、
第4図Cに示すように、後時刻データエリア
MAFTAに関連して、前時刻データエリアMBEFA
ついて上述したと同様にして、第34、第32……第
66番目のサンプリング点の相関度を表す相互相関
値信号S14を送出する。
Subsequently, the upper time data memory 25A
As shown in Figure 4C, the later time data area
In relation to M AFTA , the 34th, 32nd...
A cross-correlation value signal S14 representing the degree of correlation at the 66th sampling point is sent out.

このようにして平均値回路27から得られる相
互相関値信号S14は、第5図に示すように65個
のサンプリング点についての相関度を表すが、こ
れらのデータのうち、第33番目のデータは上側マ
イクロホン2A及び下側マイクロホン2Bに航空
機騒音が同時に到達した状態における相関度を表
すデータでなり、このデータにサンプリング点番
号「0」を割り当てる。
The cross-correlation value signal S14 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」
を割り当てる。
In addition, the 32nd, 31st...first data is data indicating that the time when the aircraft noise reaches the upper microphone 2A is before the time when it reaches the lower microphone 2B, and this data is sampled. Point number "+1", "+2"..."+32"
Assign.

さらに第34、第35……第65番目のデータは航空
機騒音が上側マイクロホン2Aに到達する時点が
下側マイクロホン2Bに到達する時点より後の時
刻であることを表すデータでなり、このデータに
サンプリング点番号「−1」、「−2」……「−
32」を割り当てる。
Further, the 34th, 35th...65th data is data indicating that the time when the aircraft noise reaches the upper microphone 2A is later than the time when the aircraft noise reaches the lower microphone 2B, and this data is sampled. Point number "-1", "-2"... "-
32”.

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

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

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

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

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

この実施例の場合相互相関演算回路部21から
送出される相互相関値信号S14は、デイジタ
ル/アナログ変換回路31においてアナログ信号
に変換された後モニタ信号S18として相互相関
判定回路20から送出するようになされ、このモ
ニタ信号S18を例えばオシユロスコープでなる
モニタに与えることにより第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 S18. By applying this monitor signal S18 to a monitor such as an oscilloscope, a correlation curve as shown in FIG. It is possible to confirm.

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

かくするにつき、差分検出回路12A及び12
Bは、各サンプリング時点における騒音検出デー
タS4A及びS4Bの変化分を表す情報だけを差分
データS5A及びS5Bとして取り出すことができ
ることにより、たとえ騒音検出信号S1A及びS
Bに周期性がありかつ騒音レベルが大きい騒音
成分が含まれている場合にも、当該騒音成分を除
外できることによりほぼ一様な周波数スペクトラ
ムを有する白色雑音としての特性をもたせること
ができ(これを白色雑音化と呼ぶ)、かくして安
定な信号レベルを有する相互相関値信号S14を
相互相関判定回路20に供給できる。
Therefore, the difference detection circuits 12A and 12
B can extract only the information representing the changes in the noise detection data S4 A and S4 B at each sampling time as the difference data S5 A and S5 B , so that even if the noise detection signals S1 A and S
1 Even if B contains a noise component with periodicity and a high noise level, by excluding the noise component, it can be given the characteristics of white noise with a nearly uniform frequency spectrum (this In this way, the cross-correlation value signal S14 having a stable signal level can be supplied to the cross-correlation determining circuit 20.

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

また同様にして航空機騒音以外の市街騒音例え
ば250〔Hz〕〜1〔kHz〕程度のサイレン、クラク
シヨン等の騒音成分についても同様にして、当該
騒音成分の周期が差分検出回路12A及び12B
における差分演算の周期より十分に長いことによ
りその影響を除去することができる。
In addition, in the same manner, for noise components of city noise other than aircraft noise, such as sirens and horns of about 250 [Hz] to 1 [kHz], the period of the noise component is determined by the difference detection circuits 12A and 12B.
This effect can be removed by making the period sufficiently longer than the period of the difference calculation in .

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

因に100〜200〔Hz〕程度の周期性がある騒音成
分を有する航空機騒音について、差分検出回路1
2A及び12Bの差分演算処理を実行せずに相互
相関判定回路20における相互相関の判定をすれ
ば、第6図において曲線S14Xに示すように、
航空機騒音に含まれている周期性騒音成分に基づ
いて相関度が当該騒音成分の周期に対応するよう
に大きく変動することにより、ピークが複数個発
生したり、ピーク波形が曖昧になつたりするのに
対して、第1図の構成によれば、第5図の最大ピ
ークサンプリング点番号JDGに示すように、航
空機の飛行位置を表す情報として、ピーク波形が
急峻な相互相関値信号を得ることができる。
Incidentally, regarding aircraft noise that has a noise component with periodicity of about 100 to 200 [Hz], the difference detection circuit 1
If the cross-correlation is determined in the cross-correlation determining circuit 20 without executing the difference calculation process of 2A and 12B, as shown by the curve S14X in FIG.
Based on the periodic noise components contained in aircraft noise, the degree of correlation changes greatly to correspond to the period of the noise components, resulting in multiple peaks or an ambiguous peak waveform. On the other hand, according to the configuration shown in Fig. 1, it is possible to obtain a cross-correlation value signal with a steep peak waveform as information representing the flight position of the aircraft, as shown by the maximum peak sampling point number JDG in Fig. 5. can.

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

これに対して差分検出回路12A及び12Bを
省略した構成によつて同じような測定をしたとこ
ろ、第7図Aに示すように、相互相関値信号S1
4Xに、スタート時から7.5秒経過するまでの間
にプロペラの騒音に基づく周期性のある騒音成分
が相互相関値信号波形S14X上に集畳すること
により、航空機の位置を表すはずのピーク値が急
峻には出て来ないのみならず、航空機が頭上を通
過する時点(すなわち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, as shown in FIG. 7A, the cross-correlation value signal S1
4X, periodic noise components based on the propeller noise are concentrated on the cross-correlation value signal waveform S14X from the time of start until 7.5 seconds have elapsed, so that the peak value that should represent the aircraft position is We obtained results in which not only the waveforms do not appear steeply, but also two peak waveforms occur at the time when the aircraft passes overhead (ie, after 3.5 seconds to 4.5 seconds).

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

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

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

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

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

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

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

第1図は本発明による航空機騒音識別装置の一
実施例を示すブロツク図、第2図はそのマイクロ
ホン装置の構成を示す略線的側面図、第3図は差
分演算の説明に供する曲線図、第4図は相互相関
演算の説明に供する略線図、第5図及び第6図は
相互相関演算結果の説明に供する信号波形図、第
7図は実験結果を示す信号波形図である。 1……航空機騒音識別装置、2A,2B……上
側、下側マイクロホン、9……騒音検出データ形
成回路、12A,12B……差分検出回路、13
A,13B……直前データラツチ回路、14A,
14B……減算回路、15……白色雑音化演算手
段、20……相互相関判定回路、21……相互相
関演算回路部、22……判定回路部。
FIG. 1 is a block diagram showing an embodiment of the aircraft noise identification device according to the present invention, FIG. 2 is a schematic side view showing the configuration of the microphone device, and FIG. 3 is a curve diagram for explaining the difference calculation. FIG. 4 is a schematic diagram for explaining the cross-correlation calculation, FIGS. 5 and 6 are signal waveform diagrams for explaining the cross-correlation calculation results, and FIG. 7 is a signal waveform diagram for showing the experimental results. DESCRIPTION OF SYMBOLS 1...Aircraft noise identification device, 2A, 2B...Upper and lower microphones, 9...Noise detection data formation circuit, 12A, 12B...Difference detection circuit, 13
A, 13B...Immediate data latch circuit, 14A,
14B... Subtraction circuit, 15... White noise conversion calculation means, 20... Cross-correlation determination circuit, 21... Cross-correlation calculation circuit section, 22... Determination circuit section.

Claims (1)

【特許請求の範囲】 1 第1及び第2のマイクロホンから得られる騒
音検出信号をサンプリングして時間の経過に従つ
て順次配列するサンプリングデータでなる第1及
び第2の騒音検出データを形成する騒音検出デー
タ形成手段と、 上記第1及び第2の騒音検出データについてそ
れぞれ、順次配列されたサンプリングデータ間の
差分を検出して第1及び第2の差分データを形成
する白色雑音化演算手段と、 上記第1及び第2の差分データの相互相関を演
算し、当該相互相関演算信号に生ずるピークに基
づいて航空機騒音を識別する相互相関判定手段と を具えることを特徴とする航空機騒音識別装置。
[Claims] 1. Noise forming 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 detection data forming means; a white noise calculation means for detecting a difference between sequentially arranged sampling data of the first and second noise detection data to form first and second difference data; An aircraft noise identification device comprising: cross-correlation determining means for calculating a cross-correlation between the first and second difference data and identifying aircraft noise based on a peak occurring in the cross-correlation calculation 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 JPH0288928A (en) 1990-03-29
JPH0515971B2 true 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)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5480009B2 (en) * 2010-05-12 2014-04-23 リオン株式会社 Noise measurement device

Also Published As

Publication number Publication date
JPH0288928A (en) 1990-03-29

Similar Documents

Publication Publication Date Title
JP3224344B2 (en) Multi-branch optical line test equipment
KR950703144A (en) METHOD AND APPARATUS FOR PROCESSING SIGNALS OF ULTRASONIC FLAW DETECTOR
CA1218457A (en) Method and apparatus for determining the endpoints of a speech utterance
JP2009103672A (en) Analysis method for discriminating between earthquake and vibration caused by noise
FR2438301A1 (en) SIGNAL DETECTION CIRCUIT OF A MAGNETIC SENSOR
JPH0515971B2 (en)
US5131743A (en) Apparatus and method for inspecting optical fibers
CN112927720A (en) Audio anomaly detection method and device
JP2007520128A (en) Process signal storage method for technical equipment
JPH05134032A (en) Method and circuit apparatus for identifying signal accompanying noise
JP4826999B2 (en) Vibration frequency measuring device
JPS58195172A (en) Radar device
SU1122641A1 (en) Device for determining properties of random process
JPH05130512A (en) Picture element defect measuring instrument for solid-state image pickup element
JPH11118939A (en) Earthquake observation system
EP0566928A2 (en) Tone detection method for telephone subset
JP4011740B2 (en) Optical signal receiving position determination device and its evaluation device
SU1146696A1 (en) Device for determining behaviour of random process
JPS58174820A (en) Laser spot checking method
JP3264089B2 (en) Human body detection device
JPS63234173A (en) Partial electric discharge measuring instrument
JP2999108B2 (en) Method and apparatus for continuous detection of waveform peak of ultrasonic flaw detection signal
KR100263927B1 (en) Automatic target tracking apparatus and formance evaluation method thereof
SU696483A2 (en) Adaptive time discretization device
JPS5821584A (en) Specific signal detection system

Legal Events

Date Code Title Description
R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

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

Free format text: PAYMENT UNTIL: 20080303

Year of fee payment: 15

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

Year of fee payment: 16

Free format text: PAYMENT UNTIL: 20090303

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

Free format text: PAYMENT UNTIL: 20090303

Year of fee payment: 16