JP6654106B2 - Sound wave monitoring device and airframe - Google Patents

Sound wave monitoring device and airframe Download PDF

Info

Publication number
JP6654106B2
JP6654106B2 JP2016122319A JP2016122319A JP6654106B2 JP 6654106 B2 JP6654106 B2 JP 6654106B2 JP 2016122319 A JP2016122319 A JP 2016122319A JP 2016122319 A JP2016122319 A JP 2016122319A JP 6654106 B2 JP6654106 B2 JP 6654106B2
Authority
JP
Japan
Prior art keywords
unit
sound wave
frequency spectrum
monitoring device
wave monitoring
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.)
Active
Application number
JP2016122319A
Other languages
Japanese (ja)
Other versions
JP2017227480A (en
Inventor
浩蔵 石前
浩蔵 石前
高行 黒川
高行 黒川
Original Assignee
ジェイ・アール・シー特機株式会社
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 ジェイ・アール・シー特機株式会社 filed Critical ジェイ・アール・シー特機株式会社
Priority to JP2016122319A priority Critical patent/JP6654106B2/en
Publication of JP2017227480A publication Critical patent/JP2017227480A/en
Application granted granted Critical
Publication of JP6654106B2 publication Critical patent/JP6654106B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Description

本発明は、音波監視装置に関し、特に、音波の周波数スペクトラムを求める装置に関する。   The present invention relates to a sound wave monitoring device, and more particularly, to a device for determining a frequency spectrum of a sound wave.

水中を探索するための航走体として、自律型無人潜水機(UUV:Unmanned Underwater Vehicle)が広く用いられている。UUVは、自らの移動距離および移動方向を測定することによって測位を行う測位装置を備えている。UUVは測位装置の測位結果を用い、予めプログラムされた航路に従って潜水しながら水中を観測する。観測の対象としては、海底、他の潜水艦、船舶等がある。   Autonomous underwater vehicles (UUV: Unmanned Underwater Vehicle) are widely used as sailboats for searching underwater. The UUV includes a positioning device that performs positioning by measuring its own moving distance and moving direction. UUV observes underwater while diving according to a pre-programmed route using the positioning result of the positioning device. Observation targets include the sea floor, other submarines, and ships.

一般にUUVには、水中を観測するための音波監視装置が搭載されている。海底の形状を観測する音波監視装置は、音波を送信し海底で反射した音波を受信し、受信信号に基づいて海底の形状データを取得する。他の潜水艦、船舶等の物標を観測する音波監視装置は、水中を伝搬する音波を受信し、受信信号に基づいて物標の位置データを取得する。   Generally, the UUV is equipped with a sound wave monitoring device for observing underwater. A sound wave monitoring device for observing the shape of the seabed transmits a sound wave, receives a sound wave reflected on the seabed, and acquires shape data of the seabed based on a received signal. A sound wave monitoring device for observing a target such as another submarine or ship receives a sound wave propagating in water and acquires position data of the target based on a received signal.

以下の特許文献1および2には水中の音波を観測する装置が記載されている。特許文献3には、このような装置に用いられる圧電素子が記載されている。   Patent Documents 1 and 2 below disclose devices for observing sound waves in water. Patent Literature 3 describes a piezoelectric element used for such an apparatus.

特開平8−61952号公報JP-A-8-61952 特開2013−160682号公報JP 2013-160682 A 特開平6−269091号公報JP-A-6-269091

音波監視装置の受信信号には観測に必要な所望音波成分の他、雑音成分が含まれる。雑音成分が大きい場合、所望音波成分に基づくデータを解析することが困難となることがある。   The received signal of the sound wave monitoring device includes a noise component in addition to a desired sound wave component required for observation. If the noise component is large, it may be difficult to analyze data based on the desired sound component.

本発明は、音波監視装置において観測される所望音波成分を雑音成分に対して十分に大きくすることを目的とする。   An object of the present invention is to make a desired sound component observed by a sound wave monitoring device sufficiently large with respect to a noise component.

本発明は、それぞれが音波を検出する複数のセンサ部と、複数の前記センサ部のそれぞれに対して設けられ、前記センサ部による検出信号に対して周波数変換処理を施す変換部と、各前記変換部から出力された周波数スペクトラムデータに対する合成処理を実行する合成処理部と、を備え、前記合成処理部による合成処理によって得られた信号のピーク周波数成分を検出する周波数検出部と、複数の前記センサ部による複数の検出信号のそれぞれに含まれる前記ピーク周波数成分の位相を揃える位相同期部と、位相が揃えられた上で前記位相同期部から出力される複数の同期化信号を合成する同期合成処理部と、前記同期合成処理部による合成信号に対して周波数変換処理を施す第2変換部と、を備えることを特徴とする。 The present invention provides a plurality of sensor units each for detecting a sound wave, a conversion unit provided for each of the plurality of sensor units, and performing a frequency conversion process on a detection signal from the sensor unit; A synthesis processing unit that performs synthesis processing on frequency spectrum data output from the unit, a frequency detection unit that detects a peak frequency component of a signal obtained by the synthesis processing by the synthesis processing unit, and a plurality of the sensors. A phase synchronizing unit that aligns the phases of the peak frequency components included in each of the plurality of detection signals, and a synthesizing process that synthesizes a plurality of synchronizing signals output from the phase synchronizing unit after the phases are aligned. and parts, and wherein Rukoto and a second converter for performing frequency conversion processing to the composite signal by the synchronization combination processing unit.

望ましくは、前記合成処理部は、各前記変換部から出力された周波数スペクトラムデータに基づいて、各前記検出信号に含まれていた同一周波数成分を掛け合わせる乗算器を含む。   Preferably, the synthesis processing unit includes a multiplier that multiplies the same frequency component included in each of the detection signals based on the frequency spectrum data output from each of the conversion units.

望ましくは、前記同期合成処理部は、複数の前記同期化信号を加算合計する。   Preferably, the synchronization synthesis processing unit adds and sums a plurality of the synchronization signals.

本発明の関連技術は、音波を検出するセンサ部と、前記センサ部による検出信号を、時系列順に所定時間長の信号に区切って分割する分割部と、前記検出信号を分割することで得られた複数の分割信号のそれぞれに対して設けられ、前記分割信号に対して周波数変換処理を施す変換部と、各前記変換部から出力された周波数スペクトラムデータに対する合成処理を実行する合成処理部と、を備え、前記合成処理部は、各前記変換部から出力された周波数スペクトラムデータに基づいて、各前記分割信号に含まれていた同一周波数成分を掛け合わせる乗算器を含む。 The related art of the present invention is obtained by dividing the detection signal by a sensor unit that detects a sound wave, a division unit that divides a detection signal by the sensor unit into signals of a predetermined time length in chronological order, and divides the detection signal. A conversion unit that is provided for each of the plurality of divided signals, and performs a frequency conversion process on the divided signal, and a synthesis processing unit that performs a synthesis process on the frequency spectrum data output from each of the conversion units. wherein the combining processing unit on the basis of the frequency spectrum data output from each of said conversion unit, including a multiplier for multiplying the same frequency components included in each of the divided signals.

本発明の関連技術は、音波を検出するセンサ部と、前記センサ部による検出信号を、時系列順に所定時間長の信号に区切って分割する分割部と、前記検出信号を分割することで得られた複数の分割信号のそれぞれに対して設けられ、前記分割信号に対して周波数変換処理を施す変換部と、各前記変換部から出力された周波数スペクトラムデータに対する合成処理を実行する合成処理部と、を備え、前記合成処理部は、各前記変換部から出力された周波数スペクトラムデータに基づいて、各前記分割信号に含まれていた同一周波数成分を加算合計する加算器を含む。 The related art of the present invention is obtained by dividing the detection signal by a sensor unit that detects a sound wave, a division unit that divides a detection signal by the sensor unit into signals of a predetermined time length in chronological order, and divides the detection signal. A conversion unit that is provided for each of the plurality of divided signals, and performs a frequency conversion process on the divided signal, and a synthesis processing unit that performs a synthesis process on the frequency spectrum data output from each of the conversion units. The synthesis processing unit includes an adder that adds and sums the same frequency components included in each of the divided signals based on the frequency spectrum data output from each of the conversion units.

望ましくは、各前記センサ部は、ピエゾフィルムセンサを備える。   Preferably, each of the sensor units includes a piezo film sensor.

望ましくは、各前記センサ部は、複数の振動子と、前記複数の振動子に基づく指向性を適応的に変化させる指向性制御部と、を備える。   Preferably, each of the sensor units includes a plurality of transducers and a directivity control unit that adaptively changes directivity based on the plurality of transducers.

本発明に係る航走体は、望ましくは、前記音波監視装置を搭載する。   The navigation body according to the present invention is desirably equipped with the sound wave monitoring device.

本発明によれば、音波監視装置において観測される所望音波成分を雑音成分に対して十分に大きくすることができる。   ADVANTAGE OF THE INVENTION According to this invention, the desired sound wave component observed by a sound wave monitoring apparatus can be made sufficiently large with respect to a noise component.

本発明の実施形態に係る音波監視装置の構成を示す図である。It is a figure showing composition of a sound wave monitoring device concerning an embodiment of the present invention. UUVに固定された音響センサおよび振動センサを模式的に示す図である。It is a figure which shows the acoustic sensor and the vibration sensor fixed to UUV typically. 第1信号処理部の構成を示す図である。FIG. 3 is a diagram illustrating a configuration of a first signal processing unit. 第2信号処理部の構成を示す図である。FIG. 3 is a diagram illustrating a configuration of a second signal processing unit. 第3信号処理部の構成を示す図である。FIG. 9 is a diagram illustrating a configuration of a third signal processing unit. 画像生成部および表示部を備える音波監視装置を示す図である。FIG. 2 is a diagram illustrating a sound wave monitoring device including an image generation unit and a display unit.

図1には、本発明の実施形態に係る音波監視装置の構成が示されている。音波監視装置は、例えばUUVに搭載され、他の潜水艦、船舶等から発せられる音波を検出する。音波監視装置は、音響センサ10−1〜10−n、音響受信部12−1〜12−n、S/N改善処理部14、信号処理ユニット16、音波データ記憶部18、振動センサ20、振動受信部22、振動高速フーリエ変換部(振動FFT)24、振動データ記憶部26、操作部28、制御部30および測位部31を備えている。   FIG. 1 shows a configuration of a sound wave monitoring device according to an embodiment of the present invention. The sound wave monitoring device is mounted on, for example, a UUV and detects sound waves emitted from other submarines, ships, and the like. The sound wave monitoring device includes sound sensors 10-1 to 10-n, sound receiving units 12-1 to 12-n, S / N improvement processing unit 14, signal processing unit 16, sound wave data storage unit 18, vibration sensor 20, vibration A receiving unit 22, a vibration fast Fourier transform unit (vibration FFT) 24, a vibration data storage unit 26, an operation unit 28, a control unit 30, and a positioning unit 31 are provided.

音響センサ10−1〜10−nのそれぞれは、音波を検出するセンサ部としての機能を有し、例えば、UUVのボデーに固定される。音響センサ10−1〜10−nには、例えば、ピエゾフィルムセンサ等の圧電素子が用いられる。音響センサ10−1〜10−nは、UUVに到来した音波を電気信号に変換し、受信信号としてそれぞれ音響受信部12−1〜12−nに出力する。音響受信部12−q(q=1〜n)は、音響センサ10−qから出力された受信信号を増幅し、さらに、受信信号を時間軸上でサンプリングされたディジタル信号に変換し、S/N改善処理部14に出力する。   Each of the acoustic sensors 10-1 to 10-n has a function as a sensor unit for detecting a sound wave, and is fixed to, for example, a UUV body. For the acoustic sensors 10-1 to 10-n, for example, a piezoelectric element such as a piezo film sensor is used. The acoustic sensors 10-1 to 10-n convert sound waves arriving at the UUV into electric signals and output the electric signals to the sound receiving units 12-1 to 12-n, respectively, as reception signals. The acoustic receiver 12-q (q = 1 to n) amplifies the received signal output from the acoustic sensor 10-q, further converts the received signal into a digital signal sampled on the time axis, and performs S / Output to the N improvement processing unit 14.

振動センサ20は、例えば、UUVのボデーに固定される。振動センサ20には、例えば、ピエゾフィルムセンサ等の圧電素子が用いられる。振動センサ20は、UUVのボデーの振動を電気信号に変換し、振動信号として振動受信部22に出力する。振動受信部22は、振動センサ20から出力された振動信号を増幅し、さらに、振動信号を時間軸上でサンプリングされたディジタル信号に変換し、そのディジタル信号をS/N改善処理部14および振動高速フーリエ変換部24に出力する。   The vibration sensor 20 is fixed to a UUV body, for example. As the vibration sensor 20, for example, a piezoelectric element such as a piezo film sensor is used. The vibration sensor 20 converts the vibration of the UUV body into an electric signal and outputs the electric signal to the vibration receiving unit 22 as a vibration signal. The vibration receiving unit 22 amplifies the vibration signal output from the vibration sensor 20, further converts the vibration signal into a digital signal sampled on a time axis, and converts the digital signal into an S / N improvement processing unit 14 and a vibration signal. Output to the fast Fourier transform unit 24.

振動高速フーリエ変換部24は、振動信号に対して高速フーリエ変換処理を施して振動周波数スペクトラムデータを生成し振動データ記憶部26に出力する。なお、振動信号は、フーリエ変換処理が施されることなく振動データ記憶部26に記憶されてもよい。   The vibration fast Fourier transform unit 24 performs a fast Fourier transform process on the vibration signal, generates vibration frequency spectrum data, and outputs it to the vibration data storage unit 26. Note that the vibration signal may be stored in the vibration data storage unit 26 without being subjected to Fourier transform processing.

測位部31は、複数方向のそれぞれについて加速度センサを備えており、UUVが航行を開始したときの位置情報と、各方向について検出された加速度に基づいて、UUVの位置情報を求める。測位部31は振動データ記憶部26に位置情報を出力する。振動データ記憶部26は、振動周波数スペクトラムデータまたは時間領域での振動信号と、測位部31から出力された位置情報とを対応付けて記憶する。   The positioning unit 31 includes an acceleration sensor for each of a plurality of directions, and obtains UUV position information based on position information when the UUV starts traveling and acceleration detected for each direction. The positioning unit 31 outputs the position information to the vibration data storage unit 26. The vibration data storage unit 26 stores the vibration frequency spectrum data or the vibration signal in the time domain and the position information output from the positioning unit 31 in association with each other.

S/N改善処理部14は、音響受信部12−1〜12−nから出力された各受信信号に対し、振動受信部22から出力された振動信号を用いて、各受信信号の信号対雑音比(S/N)を大きくする処理、すなわちS/Nを改善する処理を施す。このS/N改善処理には、時間軸上で離散的に連なる受信信号のサンプル値と、時間軸上で離散的に連なる振動信号のサンプル値に対し、サンプル時間ごとに1回の演算が実行される逐次計算を行う適応アルゴリズムがあり、例えばLMSアルゴリズム等がある。   The S / N improvement processing unit 14 uses the vibration signal output from the vibration reception unit 22 for each reception signal output from the sound reception units 12-1 to 12-n, and A process for increasing the ratio (S / N), that is, a process for improving the S / N is performed. In the S / N improvement processing, a single operation is performed for each sample time for a sample value of a received signal discretely connected on a time axis and a sample value of a vibration signal discretely connected on a time axis. There is an adaptive algorithm for performing the sequential calculation performed, for example, an LMS algorithm.

また、S/N改善処理部14は、振動信号を用いたS/N改善処理が施された各受信信号に対し、さらに、適応線スペクトル強調処理を施してもよい。適応線スペクトル強調処理は、処理対象の信号に含まれる周期的な成分を増大させるフィルタ処理である。適応線スペクトル強調処理には、時間軸上で離散的に連なる受信信号のサンプル値に対し、サンプル時間ごとに1回の演算が実行される逐次計算を行う適応アルゴリズムがある。適応線スペクトル強調処理は、例えば、特開平5−268105号公報、特開平5−75391号公報等に開示されている。   In addition, the S / N improvement processing unit 14 may further perform adaptive line spectrum enhancement processing on each of the received signals on which the S / N improvement processing using the vibration signal has been performed. The adaptive line spectrum enhancement process is a filter process for increasing a periodic component included in a signal to be processed. In the adaptive line spectrum enhancement processing, there is an adaptive algorithm that performs a sequential calculation in which a single calculation is performed for each sample time for sample values of a received signal that is discretely continuous on a time axis. The adaptive line spectrum enhancement processing is disclosed in, for example, JP-A-5-268105 and JP-A-5-75391.

後述のように、信号処理ユニット16は、1つまたは複数の音響センサによって観測された音波について、異なる処理によって生成された複数種の周波数スペクトラムデータを生成し、音波データ記憶部18に記憶する。測位部31は振動データ記憶部26の他、音波データ記憶部18にも位置情報を出力する。音波データ記憶部18は、信号処理ユニット16から出力された各周波数スペクトラムデータと、測位部31から出力された位置情報とを対応付けて記憶する。   As described later, the signal processing unit 16 generates a plurality of types of frequency spectrum data generated by different processes for the sound waves observed by one or a plurality of acoustic sensors, and stores the generated frequency spectrum data in the sound wave data storage unit 18. The positioning unit 31 outputs position information to the sound wave data storage unit 18 in addition to the vibration data storage unit 26. The sound wave data storage unit 18 stores the frequency spectrum data output from the signal processing unit 16 and the position information output from the positioning unit 31 in association with each other.

操作部28は、ユーザの操作に基づいて音波監視装置に対する制御情報を生成する。制御部30は、制御情報に基づいて音波監視装置に対する制御を行う。この制御には、例えば、電源電力のオンオフ、音波データ記憶部18および振動データ記憶部26に記憶されたデータの読み出し等がある。   The operation unit 28 generates control information for the sound wave monitoring device based on a user operation. The control unit 30 controls the sound wave monitoring device based on the control information. This control includes, for example, turning on / off the power supply, reading data stored in the sound wave data storage unit 18 and the vibration data storage unit 26, and the like.

図2(a)には、UUV32に固定された音響センサ10−1〜10−nが模式的に示されている。また、図2(b)には、音響センサ10−1〜10−nのうちの1つである音響センサ10−q、振動センサ20、およびボデー34の断面が模式的に示されている。この図では、振動センサ20はボデー34の表面に固定されており、音響センサ10−qは振動センサ20上に重ねて固定されている。これによって、振動センサ20は、UUV32のボデー34の振動を検出し、ボデー34の振動に基づく振動信号を出力する。音響センサ10−qは、UUV32に到来した音波に基づく受信信号を出力する。   FIG. 2A schematically shows the acoustic sensors 10-1 to 10-n fixed to the UUV32. FIG. 2B schematically shows a cross section of the acoustic sensor 10-q, the vibration sensor 20, and the body 34, which are one of the acoustic sensors 10-1 to 10-n. In this figure, the vibration sensor 20 is fixed to the surface of the body 34, and the acoustic sensor 10-q is fixed on the vibration sensor 20 so as to overlap. Thus, the vibration sensor 20 detects the vibration of the body 34 of the UUV 32 and outputs a vibration signal based on the vibration of the body 34. The acoustic sensor 10-q outputs a reception signal based on the sound wave arriving at the UUV32.

図1の信号処理ユニット16は、第1信号処理部16−1、第2信号処理部16−2および第3信号処理部16−3を備える。これらの信号処理部は、1つまたは複数の音響センサによって観測された音波について、異なる処理によって周波数スペクトラムデータを生成する。   The signal processing unit 16 in FIG. 1 includes a first signal processing unit 16-1, a second signal processing unit 16-2, and a third signal processing unit 16-3. These signal processing units generate frequency spectrum data by different processes for sound waves observed by one or more acoustic sensors.

図3には第1信号処理部16−1の構成が示されている。図3に示されている受信信号x(t)〜x(t)は、それぞれ、図1に示されているS/N改善処理部14から出力されたものである。「(t)」は時間領域の値であることを示す。受信信号x(t)〜x(t)は、高速フーリエ変換部(FFT)36−1〜36−nにそれぞれ入力される。受信信号x(t)〜x(t)は位相同期部42にも入力される。 FIG. 3 shows the configuration of the first signal processing unit 16-1. Received signal is shown in FIG. 3 x 1 (t) ~x n (t) are those respectively, output from the S / N improving module 14 shown in FIG. “(T)” indicates a value in the time domain. Received signal x 1 (t) ~x n ( t) are inputted fast Fourier transform unit (FFT) in 36-1 to 36-n. Received signal x 1 (t) ~x n ( t) is also input to the phase synchronization unit 42.

高速フーリエ変換部36−q(q=1〜n)は、受信信号x(t)に対して高速フーリエ変換処理を施し、受信信号x(t)についての周波数スペクトラムデータを求め、パワー演算部38−qに出力する。 The fast Fourier transform unit 36-q (q = 1 to n) performs a fast Fourier transform process on the received signal x q (t), obtains frequency spectrum data for the received signal x q (t), and performs power calculation. Output to the unit 38-q.

高速フーリエ変換処理は、時間軸上に連なるkサンプル(時間軸上の離散値の個数)の受信信号に対して実行される。高速フーリエ変換処理によって得られる1セットの周波数スペクトラムデータは、複数m個の周波数成分を表すm個の複素数を含む。   The fast Fourier transform process is performed on a received signal of k samples (the number of discrete values on the time axis) connected on the time axis. One set of frequency spectrum data obtained by the fast Fourier transform process includes m complex numbers representing a plurality of m frequency components.

高速フーリエ変換処理が施される受信信号x(t)は(数1)のように表される。 The received signal x q (t) subjected to the fast Fourier transform processing is represented as (Equation 1).

(数1)x(t)=[x(t),x(t),・・・・・,x(t)] (Equation 1) x q (t) = [x q (t 1 ), x q (t 2 ),..., X q (t k )]

受信信号x(t)に対し高速フーリエ変換を施すことにより得られる各周波数成分の実数部および虚数部は、それぞれ(数2)および(数3)のように表される。 The real part and the imaginary part of each frequency component obtained by performing the fast Fourier transform on the received signal x q (t) are expressed as (Equation 2) and (Equation 3), respectively.

(数2)R=[R(f),R(f),・・・・・,R(f)]
(数3)I=[I(f),I(f),・・・・・,I(f)]
(Equation 2) R q = [R q (f 1 ), R q (f 2 ),..., R q (f m )]
(Equation 3) I q = [I q (f 1 ), I q (f 2 ),..., I q (f m )]

受信信号x(t)の周波数f、f、・・・・・、fの各成分は、jを虚数単位とする複素数によって、R(f)+j・I(f)、R(f)+j・I(f)、R(f)+j・I(f)、・・・・・、R(f)+j・I(f)のように表される。 Frequency f 1, f 2 of the received signal x q (t), ·····, each component of f m is the complex number of the j and imaginary unit, R q (f 1) + j · I q (f 1 ), R q (f 2 ) + j · I q (f 2 ), R q (f 3 ) + j · I q (f 3 ),..., R q (f m ) + j · I q (f m ).

パワー演算部38−qは、各周波数成分についてパワー値を求める。パワー値は周波数成分を表す複素数の絶対値の自乗として定義される。すなわち、パワー演算部38−qが出力するパワー値は(数4)のように表される。   The power calculator 38-q obtains a power value for each frequency component. The power value is defined as the square of the absolute value of a complex number representing a frequency component. That is, the power value output by the power calculation unit 38-q is expressed as (Equation 4).

(数4)P=[L(f),L(f),・・・・・,L(f)]
ただし、1からmのうちいずれかの整数をhとして、L(f)=R(f+I(fである。
(Equation 4) P q = [L q (f 1 ), L q (f 2 ),..., L q (f m )]
However, as h any integer of from 1 m, which is L q (f h) = R q (f h) 2 + I q (f h) 2.

乗算器40は、複数の周波数スペクトラムデータに対する合成処理を実行する合成処理部として機能する。すなわち、乗算器40は、パワー演算部38−1〜38−nから出力されたパワー値の同一周波数成分を掛け合わせ、m個の新たなパワー値を含む積算周波数スペクトラムデータPαを生成し、図1の音波データ記憶部18に出力する。積算周波数スペクトラムデータPαは、(数5)のように表される。ただし、記号Πは、q=1〜nについて総てを掛け合わせることを意味する。   The multiplier 40 functions as a synthesis processing unit that executes a synthesis process on a plurality of frequency spectrum data. That is, the multiplier 40 multiplies the same frequency components of the power values output from the power calculators 38-1 to 38-n to generate integrated frequency spectrum data Pα including m new power values. The data is output to the first sound wave data storage unit 18. The integrated frequency spectrum data Pα is expressed as (Equation 5). Here, the symbol Π means that all of q = 1 to n are multiplied.

(数5)Pα=[ΠL(f),ΠL(f),・・・・・,ΠL(f)] (Equation 5) Pα = [ΠL q (f 1 ), ΠL q (f 2 ),..., ΠL q (f m )]

ピーク周波数検出部50は、乗算器40から出力される積算周波数スペクトラムデータPαに基づいて、積算周波数スペクトラムのピーク周波数を検出する。ここで、ピーク周波数とは、積算周波数スペクトラムデータPαが示す周波数成分のうち、最も大きい周波数成分の周波数をいう。ピーク周波数検出部50は、ピーク周波数を位相同期部42に出力する。   The peak frequency detector 50 detects the peak frequency of the integrated frequency spectrum based on the integrated frequency spectrum data Pα output from the multiplier 40. Here, the peak frequency refers to the frequency of the largest frequency component among the frequency components indicated by the integrated frequency spectrum data Pα. The peak frequency detection section 50 outputs the peak frequency to the phase synchronization section 42.

位相同期部42は、受信信号x(t)〜x(t)に対し、ピーク周波数検出部50によって求められたピーク周波数の成分について位相を揃える処理を実行する。位相同期部42は、このような位相同期処理が施された受信信号x(t)〜x(t)を加算器44に出力する。 The phase synchronization unit 42 executes a process of aligning the phases of the received signals x 1 (t) to x n (t) with respect to the components of the peak frequency determined by the peak frequency detection unit 50. The phase synchronizer 42 outputs the received signals x 1 (t) to x n (t) on which the phase synchronization processing has been performed to the adder 44.

加算器44は、同期化信号としての受信信号x(t)〜x(t)を合成する合成処理部としての機能を有する。すなわち、加算器44は、位相同期化処理が施された受信信号x(t)〜x(t)を加算合計し、それによって得られた加算合計・受信信号(合成信号)を高速フーリエ変換部46に出力する。高速フーリエ変換部46は、加算合計・受信信号に対して高速フーリエ変換処理を施し、加算合計・受信信号についての周波数スペクトラムデータを求め、パワー演算部48に出力する。パワー演算部48は、各周波数成分についてパワー値を求め、m個のパワー値を含む位相同期周波数スペクトラムデータPβを図1の音波データ記憶部18に出力する。 The adder 44 has a function as a combining processing unit for combining the received signal as a synchronization signal x 1 (t) ~x n ( t). That is, the adder 44 adds and totals the received signals x 1 (t) to x n (t) that have been subjected to the phase synchronization processing, and outputs the added total / received signal (synthesized signal) obtained by the fast Fourier transform. Output to the conversion unit 46. The fast Fourier transform unit 46 performs a fast Fourier transform process on the summation / reception signal, obtains frequency spectrum data on the summation / reception signal, and outputs the frequency spectrum data to the power calculation unit 48. The power calculation unit 48 obtains a power value for each frequency component and outputs phase-locked frequency spectrum data Pβ including m power values to the sound wave data storage unit 18 in FIG.

高速フーリエ変換部46およびパワー演算部48が実行する具体的な処理は、それぞれ、高速フーリエ変換部36−qおよびパワー演算部38−qが実行する上述の処理と同様である。   Specific processes executed by the fast Fourier transform unit 46 and the power operation unit 48 are the same as the above-described processes executed by the fast Fourier transform unit 36-q and the power operation unit 38-q, respectively.

なお、ここでは、ピーク周波数検出部50がピーク周波数を検出し、位相同期部42に出力する処理について説明した。このような処理を実行する代わりに、ユーザの操作によってピーク周波数が位相同期部42に入力されてもよい。   Here, the processing in which the peak frequency detection unit 50 detects the peak frequency and outputs the peak frequency to the phase synchronization unit 42 has been described. Instead of performing such processing, the peak frequency may be input to the phase synchronization unit 42 by a user operation.

第1信号処理部16−1が実行する処理によって、積算周波数スペクトラムデータPα、および位相同期周波数スペクトラムデータPβが図1の音波データ記憶部18に出力される。音波データ記憶部18は、積算周波数スペクトラムデータPαおよび位相同期周波数スペクトラムデータPβを、測位部31から出力された位置情報に対応付けて記憶する。   By the processing executed by the first signal processing unit 16-1, the integrated frequency spectrum data Pα and the phase synchronization frequency spectrum data Pβ are output to the sound wave data storage unit 18 in FIG. The sound wave data storage unit 18 stores the integrated frequency spectrum data Pα and the phase synchronization frequency spectrum data Pβ in association with the position information output from the positioning unit 31.

図4には第2信号処理部16−2の構成が示されている。図4に示されている受信信号x(t)〜x(t)は、それぞれ、図1に示されているS/N改善処理部14から出力されたものである。受信信号x(t)〜x(t)は、それぞれ、高速フーリエ変換部52−1〜52−nに入力される。高速フーリエ変換部52−q(q=1〜n)は、受信信号x(t)に対して高速フーリエ変換処理を施し、受信信号x(t)についての周波数スペクトラムデータを求め、複素加算部54に出力する。 FIG. 4 shows the configuration of the second signal processing unit 16-2. Received signal is shown in FIG. 4 x 1 (t) ~x n (t) are those respectively, output from the S / N improving module 14 shown in FIG. Received signal x 1 (t) ~x n ( t) are respectively input to the fast Fourier transform unit 52-1 to 52-n. The fast Fourier transform unit 52-q (q = 1 to n) performs a fast Fourier transform process on the received signal x q (t), obtains frequency spectrum data on the received signal x q (t), and performs complex addition. Output to the unit 54.

高速フーリエ変換部52−1〜52−nが実行する高速フーリエ変換処理は、図3の高速フーリエ変換部36−1〜36−nが実行する高速フーリエ変換処理と同様である。すなわち、受信信号x(t)に対し高速フーリエ変換が施されることにより、m個の周波数成分が得られる。各周波数成分の実数部および虚数部は、それぞれ(数2)および(数3)のように表され、受信信号x(t)の周波数f、f、・・・・・、fの各成分は、複素数によって、R(f)+j・I(f)、R(f)+j・I(f)、R(f)+j・I(f)、・・・・・、R(f)+j・Iq(f)のように表される。 The fast Fourier transform performed by the fast Fourier transforms 52-1 to 52-n is the same as the fast Fourier transform performed by the fast Fourier transforms 36-1 to 36-n in FIG. That is, by performing fast Fourier transform on the received signal x q (t), m frequency components are obtained. The real and imaginary parts of each frequency component is represented as the respective (number 2) and (Equation 3), the frequency f 1, f 2 of the received signal x q (t), ·····, f m Are represented by R q (f 1 ) + j · I q (f 1 ), R q (f 2 ) + j · I q (f 2 ), and R q (f 3 ) + j · I q (f 3), ... it is expressed as R q (f m) + j · Iq (f m).

複素加算部54は、受信信号x(t)の周波数f、f、・・・・・、fの成分それぞれにつき、実数部および虚数部のそれぞれを加算合計する。実数部の加算合計値は(数6)のように表され、虚数部の加算合計値は(数7)のように表される。ただし、記号Σは、q=1〜nについて総てを加算合計することを意味する。 Complex addition unit 54, the frequency f 1, f 2 of the received signal x q (t), ·····, per each component of f m, adds the sum of each of the real and imaginary parts. The addition sum of the real part is represented as (Equation 6), and the addition sum of the imaginary part is represented as (Equation 7). However, the symbol Σ means that all of q = 1 to n are added and summed.

(数6)R=[ΣR(f),ΣR(f),・・・・・,ΣR(fm)]
(数7)I=[ΣI(f),ΣI(f),・・・・・,ΣI(fm)]
(Equation 6) R = [ΣR q (f 1 ), ΣR q (f 2 ),..., ΣR q (fm)]
(Equation 7) I = [ΣI q (f 1 ), ΣI q (f 2 ),..., ΣI q (fm)]

複素加算部54が実行する処理は、横軸に実数部をとり縦軸に虚数部をとった複素平面において、各周波数成分の複素ベクトルを加算合計する処理であるといえる。複素加算部54は、受信信号x(t)の周波数f、f、・・・・・、fの成分のそれぞれについて、実数部の加算合計値、および、虚数部の加算合計値をパワー演算部56に出力する。パワー演算部56は、各周波数成分について、実数部の自乗および虚数部の自乗を加算して、加算合計周波数スペクトラムデータPγを求める。加算合計周波数スペクトラムデータPγは、(数8)のように表される。 The processing executed by the complex addition unit 54 can be said to be the processing of adding and summing the complex vectors of the respective frequency components on a complex plane having the real part on the horizontal axis and the imaginary part on the vertical axis. The complex addition unit 54, a frequency f 1, f 2 of the received signal x q (t), ·····, for each of the components of f m, sum total value of the real part, and, adding the total value of the imaginary part Is output to the power calculation unit 56. The power calculator 56 adds the square of the real part and the square of the imaginary part for each frequency component to obtain the sum total frequency spectrum data Pγ. The addition total frequency spectrum data Pγ is represented as (Equation 8).

(数8)Pγ=[{ΣR(f)}+{ΣI(f)},{ΣR(f)}+{ΣI(f)},・・・・・,{ΣR(f)}+{ΣI(f)}(Number 8) Pγ = [{ΣR q (f 1)} 2 + {ΣI q (f 1)} 2, {ΣR q (f 2)} 2 + {ΣI q (f 2)} 2, ··· ··, {ΣR q (f m )} 2 + {ΣI q (f m)} 2]

パワー演算部56は、加算合計周波数スペクトラムデータPγを図1の音波データ記憶部18に出力する。   The power calculation unit 56 outputs the added total frequency spectrum data Pγ to the sound wave data storage unit 18 in FIG.

第2信号処理部16−2が実行する処理によって、加算合計周波数スペクトラムデータPγが図1の音波データ記憶部18に出力される。音波データ記憶部18は、加算合計周波数スペクトラムデータPγを、測位部31から出力された位置情報に対応付けて記憶する。   By the processing executed by the second signal processing unit 16-2, the sum total frequency spectrum data Pγ is output to the sound wave data storage unit 18 in FIG. The sound wave data storage unit 18 stores the added total frequency spectrum data Pγ in association with the position information output from the positioning unit 31.

積算周波数スペクトラムデータPα、位相同期周波数スペクトラムデータPβ、および加算合計周波数スペクトラムデータPγには次のような特徴がある。上記のように、積算周波数スペクトラムデータPαは、パワー演算部38−1〜38−nから出力された周波数スペクトラムデータの同一周波数成分値を掛け合わせる処理を経て得られる。また、位相同期周波数スペクトラムデータPβは、ピーク周波数の成分について位相を揃えた上で受信信号x(t)〜x(t)を加算合計し、これによって得られた信号に対して高速フーリエ変換処理を施す処理を経て得られる。そして、加算合計周波数スペクトラムデータPγは、高速フーリエ変換部52−1〜52−nから出力された周波数スペクトラムデータの同一周波数成分値を、実数部および虚数部のそれぞれについて加算合計する処理を経て得られる。 The integrated frequency spectrum data Pα, the phase-locked frequency spectrum data Pβ, and the sum total frequency spectrum data Pγ have the following features. As described above, the integrated frequency spectrum data Pα is obtained through a process of multiplying the same frequency component value of the frequency spectrum data output from the power calculators 38-1 to 38-n. Further, the phase-locked frequency spectrum data Pβ is obtained by adding the received signals x q (t) to x n (t) after aligning the phases of the components of the peak frequency, and applying a fast Fourier transform to the obtained signal. It is obtained through a process of performing a conversion process. Then, the sum total frequency spectrum data Pγ is obtained through a process of summing the same frequency component values of the frequency spectrum data output from the fast Fourier transform units 52-1 to 52-n for each of the real part and the imaginary part. Can be

一般に、雑音成分には観測位置と大きさとの間に規則性がなく、雑音成分は観測位置によって大きさが異なることが通常である。これに対し、特定の周波数成分を有する音波は、観測位置に対する大きさの違いは雑音成分ほど顕著ではない。   Generally, the noise component has no regularity between the observation position and the size, and the noise component usually has a different size depending on the observation position. On the other hand, a sound wave having a specific frequency component has a smaller difference in magnitude with respect to the observation position than a noise component.

したがって、異なる位置に固定された複数の音響センサで検出された雑音成分の値を掛け合わせた値と、その複数の音響センサで検出された音波の周波数成分値を掛け合わせた値との比率であるS/Nは、音響センサの数が多い程大きくなる。   Therefore, the ratio of the value obtained by multiplying the value of the noise component detected by the plurality of acoustic sensors fixed at different positions to the value obtained by multiplying the value of the frequency component of the sound wave detected by the plurality of acoustic sensors is obtained. A certain S / N increases as the number of acoustic sensors increases.

同様の理由により、異なる位置に固定された複数の音響センサで検出された雑音成分の値を加算合計した値と、その複数の音響センサで検出された音波の周波数成分値を加算合計した値との比率であるS/Nは、音響センサの数が多い程大きくなる。   For the same reason, a value obtained by adding and adding values of noise components detected by a plurality of acoustic sensors fixed at different positions, and a value obtained by adding and adding frequency component values of sound waves detected by the plurality of acoustic sensors are used. S / N, which is the ratio of, increases as the number of acoustic sensors increases.

積算周波数スペクトラムデータPαは、このような原理に基づき、異なる位置に配置された複数の音響センサによって観測された周波数スペクトラムデータの各周波数成分を掛け合わせることで、1つの位置で観測される周波数スペクトラムに対しS/Nを改善するものである。   Based on such a principle, the integrated frequency spectrum data Pα is obtained by multiplying each frequency component of the frequency spectrum data observed by a plurality of acoustic sensors arranged at different positions to obtain a frequency spectrum observed at one position. To improve S / N.

1つの音響センサによって観測される周波数スペクトラムのS/Nをα[dB]とすると、積算周波数スペクトラムのS/Nは、αn=α+α(n−1)[dB]となる。これは、観測される音波の大きさがn乗されることに基づく。したがって、デシベルスケールでS/Nがα(n−1)だけ増加しS/Nが改善される。   Assuming that the S / N of the frequency spectrum observed by one acoustic sensor is α [dB], the S / N of the integrated frequency spectrum is αn = α + α (n−1) [dB]. This is based on the magnitude of the observed sound wave being raised to the nth power. Therefore, the S / N is increased by α (n−1) on the decibel scale, and the S / N is improved.

位相同期周波数スペクトラムデータPβおよび加算合計周波数スペクトラムデータPγは、同様の原理に基づき、異なる位置に配置された複数の音響センサによって観測された周波数スペクトラムの同一周波数成分を加算合計することで、1つの位置で観測される周波数スペクトラムよりもS/Nを良好とするものである。位相同期周波数スペクトラムデータPβが、ある特定の周波数成分の位相を揃えて加算合計するものであるのに対し、加算合計周波数スペクトラムデータPγは、各周波数成分の位相はそのままとした上で、各周波数成分の複素ベクトルを複素平面上で加算合計するものである点で異なる。しかし、これらの周波数スペクトラムデータにおいて得られる効果は同様である。   The phase-locked frequency spectrum data Pβ and the sum total frequency spectrum data Pγ are obtained by adding and summing the same frequency components of the frequency spectrum observed by a plurality of acoustic sensors arranged at different positions based on the same principle. This is to make the S / N better than the frequency spectrum observed at the position. While the phase-synchronous frequency spectrum data Pβ is to add and sum the phases of certain specific frequency components, the sum total frequency spectrum data Pγ is obtained by keeping the phase of each frequency component as it is, The difference is that the complex vectors of the components are added and summed on a complex plane. However, the effect obtained with these frequency spectrum data is the same.

1つの音響センサによって観測される周波数スペクトラムのS/Nをα[dB]とすると、位相同期周波数スペクトラムおよび加算合計周波数スペクトラムのそれぞれのS/Nは、α+10log(n)[dB]となる。これは、観測される音波の大きさがn倍となることに基づく。したがって、デシベルスケールで10log(n)だけS/Nが増加し、S/Nが改善される。   Assuming that the S / N of the frequency spectrum observed by one acoustic sensor is α [dB], the S / N of each of the phase locked frequency spectrum and the sum total frequency spectrum is α + 10 log (n) [dB]. This is based on the fact that the magnitude of the observed sound wave becomes n times. Therefore, the S / N is increased by 10 log (n) on the decibel scale, and the S / N is improved.

図5には第3信号処理部16−3の構成が示されている。図5に示されている受信信号x(t)は、図1に示されている音響センサ10−1の検出値に基づきS/N改善処理部14から出力されたものであり、分割部58に入力される。 FIG. 5 shows the configuration of the third signal processing unit 16-3. The received signal x 1 (t) shown in FIG. 5 is output from the S / N improvement processing unit 14 based on the detection value of the acoustic sensor 10-1 shown in FIG. 58 is input.

分割部58は、音響センサ10−1による検出信号としての受信信号を、時系列順に所定時間長の信号に区切って分割する機能を有する。すなわち、分割部58は、時間軸上で連なる受信信号x(t)をpサンプルごとに区切り、x11(t)、x12(t)、x13(t)、・・・・・、x1s(t)の順に、pサンプルの区切りでx11(t)〜x1s(t)を、それぞれ、高速フーリエ変換部60−1〜60−sに出力する。すなわち、分割部58は、時間軸上で連なるs・p個のサンプル値をs個に分割し、分割信号x11(t)〜x1s(t)を出力する。x11(t)〜x1s(t)のそれぞれには、時間軸上で連なるp個のサンプル値が含まれる。 The dividing unit 58 has a function of dividing a reception signal as a detection signal from the acoustic sensor 10-1 into a signal of a predetermined time length in a time series manner. That is, the dividing unit 58 divides the received signal x 1 (t) that continues on the time axis into p samples, and obtains x 11 (t), x 12 (t), x 13 (t),. In the order of x 1s (t), x 11 (t) to x 1s (t) are output to the fast Fourier transform units 60-1 to 60-s at p-sample breaks, respectively. That is, the dividing unit 58 divides the sp sample values continuous on the time axis into s, and outputs the divided signals x 11 (t) to x 1s (t). Each of x 11 (t) to x 1s (t) includes p sample values connected on the time axis.

高速フーリエ変換部60−u(u=1〜s)は、分割信号x1u(t)に対して高速フーリエ変換処理を施し、分割信号x1u(t)についての周波数スペクトラムデータを求め、パワー演算部62−uに出力する。 The fast Fourier transform unit 60-u (u = 1 to s ) performs a fast Fourier transform process on the divided signal x 1u (t), obtains frequency spectrum data on the divided signal x 1u (t), and performs power calculation. Output to the unit 62-u.

高速フーリエ変換部60−uが実行する高速フーリエ変換処理は、時間軸上に連なるpサンプルの受信信号に対して実行される。高速フーリエ変換処理によって得られる1セットの周波数スペクトラムデータは、複数m個の周波数成分を表すm個の複素数を含む。   The fast Fourier transform process performed by the fast Fourier transform unit 60-u is performed on a received signal of p samples connected on the time axis. One set of frequency spectrum data obtained by the fast Fourier transform process includes m complex numbers representing a plurality of m frequency components.

ここで、高速フーリエ変換部60−uにおいて高速フーリエ変換処理が施される分割信号x1u(t)は(数9)のように表される。 Here, the divided signal x 1u (t) to be subjected to the fast Fourier transform processing in the fast Fourier transform unit 60-u is expressed as in ( Equation 9).

(数9)x1u(t)=[x1u(t),x1u(t),・・・・・,x1u(t)] ( Equation 9) x 1u (t) = [x 1u (t 1 ), x 1u (t 2 ),..., X 1u (t p )]

分割信号x1u(t)に対し高速フーリエ変換を施すことにより得られる各周波数成分の実数部および虚数部は、それぞれ(数10)および(数11)のように表される。 The real part and the imaginary part of each frequency component obtained by performing the fast Fourier transform on the divided signal x 1u (t) are expressed as (Equation 10) and (Equation 11), respectively.

(数10)R1u=[R1u(f),R1u(f),・・・・・,R1u(f)]
(数11)I1u=[I1u(f),I1u(f),・・・・・,I1u(f)]
( Equation 10) R 1u = [R 1u (f 1 ), R 1u (f 2 ),..., R 1u (f m )]
( Equation 11) I 1u = [I 1u (f 1 ), I 1u (f 2 ),..., I 1u (f m )]

分割信号x1u(t)の周波数f、f、・・・・・、fの各成分は、jを虚数単位とする複素数によって、R1u(f)+j・I1u(f)、R1u(f)+j・I1u(f)、R1u(f)+j・I1u(f)、・・・・・、R1u(f)+j・I1u(f)のように表される。 Each component of the frequency f 1, f 2, ·····, f m split signal x 1u (t) is the complex number of the j and imaginary unit, R 1u (f 1) + j · I 1u (f 1 ), R 1u (f 2 ) + j · I 1u (f 2 ), R 1u (f 3 ) + j · I 1u (f 3 ),..., R 1u (f m ) + j · I 1u (f m ).

パワー演算部62−uは、各周波数成分についてパワー値を求める。パワー値は周波数成分を表す複素数の絶対値の自乗として定義される。すなわち、パワー演算部62−uが出力するパワー値は(数12)のように表される。   The power calculation unit 62-u calculates a power value for each frequency component. The power value is defined as the square of the absolute value of a complex number representing a frequency component. That is, the power value output from the power calculation unit 62-u is represented as (Equation 12).

(数12)P1u=[L1u(f),L1u(f),・・・・・,L1u(f)]
ただし、1からmのうちいずれかの整数をhとして、L1u(f)=R1u(f+I1u(fである。
( Equation 12) P 1u = [L 1u (f 1 ), L 1u (f 2 ),..., L 1u (f m )]
However, as h any integer of from 1 m, which is L 1u (f h) = R 1u (f h) 2 + I 1u (f h) 2.

乗算器64および加算器66のそれぞれは、複数の周波数スペクトラムデータに対する合成処理を実行する合成処理部として機能する。すなわち、乗算器64は、パワー演算部62−1〜62−sから出力されたパワー値の同一周波数成分を掛け合わせ、m個の新たなパワー値を含む分割・積算周波数スペクトラムデータPδを生成し、図1の音波データ記憶部18に出力する。分割・積算周波数スペクトラムデータPδは、(数13)のように表される。ただし、記号Πは、u=1〜sについて総てを掛け合わせることを意味する。   Each of the multiplier 64 and the adder 66 functions as a synthesis processing unit that executes synthesis processing on a plurality of frequency spectrum data. That is, the multiplier 64 multiplies the same frequency components of the power values output from the power calculation units 62-1 to 62-s to generate divided / integrated frequency spectrum data Pδ including m new power values. Are output to the sound wave data storage unit 18 in FIG. The divided / integrated frequency spectrum data Pδ is expressed as (Equation 13). Here, the symbol Π means that all of u = 1 to s are multiplied.

(数13)Pδ=[ΠL1u(f),ΠL1u(f),・・・・・,ΠL1u(f)] (Equation 13) Pδ = [ΠL 1u (f 1 ), ΠL 1u (f 2 ),..., ΠL 1u (f m )]

加算器66は、パワー演算部62−1〜62−sから出力されたパワー値の同一周波数成分を加算合計し、m個の新たなパワー値を含む分割・加算合計周波数スペクトラムデータPεを生成し、図1の音波データ記憶部18に出力する。分割・加算合計周波数スペクトラムデータPεは、(数14)のように表される。ただし、記号Σは、u=1〜sについて加算合計することを意味する。   The adder 66 adds and sums the same frequency components of the power values output from the power calculation units 62-1 to 62-s, and generates division / sum total frequency spectrum data Pε including m new power values. Are output to the sound wave data storage unit 18 in FIG. The division / addition total frequency spectrum data Pε is expressed as (Equation 14). However, the symbol Σ means that addition and sum are performed for u = 1 to s.

(数14)Pε=[ΣL1u(f),ΣL1u(f),・・・・・,ΣL1u(f)] (Equation 14) Pε = [ΣL 1u (f 1 ), ΣL 1u (f 2 ),..., ΣL 1u (f m )]

第3信号処理部16−3が実行する処理によって、分割・積算周波数スペクトラムデータPδ、および分割・加算合計周波数スペクトラムデータPεが図1の音波データ記憶部18に出力される。音波データ記憶部18は、分割・積算周波数スペクトラムデータPδおよび分割・加算合計周波数スペクトラムデータPεを、測位部31から出力された位置情報に対応付けて記憶する。   By the processing executed by the third signal processing unit 16-3, the divided / integrated frequency spectrum data Pδ and the divided / added total frequency spectrum data Pε are output to the sound wave data storage unit 18 in FIG. The sound wave data storage unit 18 stores the divided / integrated frequency spectrum data Pδ and the divided / added total frequency spectrum data Pε in association with the position information output from the positioning unit 31.

なお、ここでは、音響センサ10−1に基づく受信信号x(t)に基づいて、分割・積算周波数スペクトラムデータPδおよび分割・加算合計周波数スペクトラムデータPεを生成する構成について説明した。第3信号処理部16−3は、同様の構成によって、その他の音響センサ10−1〜10−nのいずれかに基づく受信信号に基づいて、分割・積算周波数スペクトラムデータPδおよび分割・加算合計周波数スペクトラムデータPεを生成するものであってもよい。 Here, the configuration has been described in which the divided / integrated frequency spectrum data Pδ and the divided / added total frequency spectrum data Pε are generated based on the reception signal x 1 (t) based on the acoustic sensor 10-1. The third signal processing unit 16-3 has the same configuration, and based on the received signal based on any of the other acoustic sensors 10-1 to 10-n, the divided / integrated frequency spectrum data Pδ and the divided / added total frequency The spectrum data Pε may be generated.

分割・積算周波数スペクトラムデータPδ、および分割・加算合計周波数スペクトラムデータPεには次のような特徴がある。上記のように、分割・積算周波数スペクトラムデータPδは、時間軸上で連なるx(t)をpサンプルごとに区切り、x11(t)、x12(t)、x13(t)、・・・・・、x1s(t)の順に、pサンプルの区切りでx11(t)〜x1s(t)を、それぞれ、高速フーリエ変換部60−1〜60−sに出力し、高速フーリエ変換部60−1〜60−nから出力された周波数スペクトラムデータの同一周波数成分値を掛け合わせる処理を経て得られる。そして、分割・加算合計周波数スペクトラムデータは、高速フーリエ変換部60−1〜60−nから出力された周波数スペクトラムデータの同一周波数成分値を加算合計する処理を経て得られる。 The divided / integrated frequency spectrum data Pδ and the divided / added total frequency spectrum data Pε have the following features. As described above, the divided / integrated frequency spectrum data Pδ is obtained by dividing x 1 (t), which is continuous on the time axis, into p samples, and obtaining x 11 (t), x 12 (t), x 13 (t),. ..., in the order of x 1s (t), x 11 in separator p samples (t) ~x 1s (t) , respectively, and outputs the fast Fourier transform unit 60-1 to 60-s, fast Fourier It is obtained through a process of multiplying the same frequency component values of the frequency spectrum data output from the conversion units 60-1 to 60-n. The divided / added total frequency spectrum data is obtained through a process of adding and adding the same frequency component values of the frequency spectrum data output from the fast Fourier transform units 60-1 to 60-n.

一般に、雑音成分には時間についての規則性がなく、雑音成分は時間経過に対してランダムに変化することが通常である。一方、特定の周波数成分を有する音波は、時間経過に対する大きさの違いは雑音成分ほど顕著ではない。   Generally, the noise component does not have a regularity with respect to time, and the noise component usually changes randomly with the passage of time. On the other hand, a sound wave having a specific frequency component has a smaller difference in magnitude over time than the noise component.

したがって、異なる複数の分割時間帯で検出された雑音成分の値を掛け合わせた値と、その複数の分割時間帯で検出された音波の周波数成分値を掛け合わせた値との比率であるS/Nは、分割時間帯の数sが多い程大きくなる。   Accordingly, S / S is a ratio of a value obtained by multiplying the value of the noise component detected in a plurality of different time periods and a value obtained by multiplying the frequency component value of the sound wave detected in the plurality of time periods. N increases as the number s of divided time zones increases.

同様の理由により、異なる複数の分割時間帯で検出された雑音成分の値を加算合計した値と、その複数の分割時間帯で検出された音波の周波数成分値を加算合計した値との比率であるS/Nは、分割時間帯の数sが多い程大きくなる。   For the same reason, the ratio of the value obtained by adding and summing the values of the noise components detected in a plurality of different divided time periods and the value obtained by adding and summing the frequency component values of the sound waves detected in the plurality of divided time periods is used. A certain S / N increases as the number s of divided time zones increases.

分割・積算周波数スペクトラムデータPδは、このような原理に基づき、異なる複数の分割時間帯のそれぞれで観測された周波数スペクトラムデータの各周波数成分を掛け合わせることで、1つの分割時間帯で観測される周波数スペクトラムに対しS/Nを良好とするものである。   The divided / integrated frequency spectrum data Pδ is observed in one divided time zone by multiplying each frequency component of the frequency spectrum data observed in each of a plurality of different divided time zones based on such a principle. This is to improve the S / N with respect to the frequency spectrum.

1つの分割時間帯によって観測される周波数スペクトラムのS/Nをα[dB]とすると、分割・積算周波数スペクトラムのS/Nは、αs=α+α(s−1)[dB]となる。これは、観測される音波の大きさがs乗されることに基づく。したがって、デシベルスケールでS/Nがα(s−1)だけ増加しS/Nが改善される。   Assuming that the S / N of the frequency spectrum observed in one divided time zone is α [dB], the S / N of the divided / integrated frequency spectrum is αs = α + α (s−1) [dB]. This is based on the magnitude of the observed sound wave being raised to the s power. Therefore, the S / N is increased by α (s−1) on the decibel scale, and the S / N is improved.

分割・加算合計周波数スペクトラムPεは、同様の原理に基づき、異なる複数の分割時間帯で観測された周波数スペクトラムの同一周波数成分を加算合計することで、1つの分割時間帯で観測される周波数スペクトラムよりもS/Nを良好とするものである。   Based on the same principle, the division / addition total frequency spectrum Pε is obtained by adding and adding the same frequency components of the frequency spectrum observed in a plurality of different division time zones to obtain a frequency spectrum observed in one division time zone. Also improve S / N.

1つの分割時間帯で観測される周波数スペクトラムのS/Nをα[dB]とすると、分割・加算合計周波数スペクトラムのS/Nは、α+10log(s)[dB]となる。これは、観測される音波の大きさがs倍となることに基づく。したがって、デシベルスケールで10log(s)だけS/Nが増加し、S/Nが改善される。   Assuming that the S / N of the frequency spectrum observed in one divided time zone is α [dB], the S / N of the divided and added total frequency spectrum is α + 10 log (s) [dB]. This is based on the fact that the magnitude of the observed sound wave becomes s times. Therefore, the S / N is increased by 10 log (s) on the decibel scale, and the S / N is improved.

図1に戻って音波監視装置について説明する。音波監視装置を搭載したUUVは、予めプログラムされた航路を航行した後、ユーザによって回収される。ユーザは、音波監視用のコンピュータを用いて音波データ記憶部18および振動データ記憶部26に記憶されたデータを回収する。   Returning to FIG. 1, the sound wave monitoring device will be described. A UUV equipped with an acoustic wave monitoring device is collected by a user after navigating a pre-programmed route. The user collects the data stored in the sound wave data storage unit 18 and the vibration data storage unit 26 using a computer for sound wave monitoring.

音波データ記憶部18に記憶された積算周波数スペクトラムデータPα、位相同期周波数スペクトラムデータPβ、加算合計周波数スペクトラムデータPγ、分割・積算周波数スペクトラムデータPδ、および分割・加算合計周波数スペクトラムデータPεは位置情報に対応付けられている。音波データ記憶部18に記憶されたデータによって、UUVの各航行位置についてUUVに到来した音波についての解析が行われる。例えば、音波監視用のコンピュータは、これらの周波数スペクトラムデータのいずれか1つ、または、組み合わせに基づく音波監視用画像データを生成し、音波監視用画像データに基づく画像を表示してもよい。   The integrated frequency spectrum data Pα, the phase-locked frequency spectrum data Pβ, the sum total frequency spectrum data Pγ, the divided / integrated frequency spectrum data Pδ, and the divided / added total frequency spectrum data Pε stored in the sound wave data storage unit 18 are included in the position information. Are associated. Based on the data stored in the sound wave data storage unit 18, the sound wave arriving at the UUV is analyzed for each UUV navigation position. For example, the computer for sound wave monitoring may generate sound wave monitoring image data based on one or a combination of these frequency spectrum data, and display an image based on the sound wave monitoring image data.

また、振動データ記憶部26に記憶された振動周波数スペクトラムデータは位置情報に対応付けられている。振動データ記憶部26に記憶されたデータによって、UUVの各航行位置についてUUVのボデーの振動についての解析が行われる。例えば、音波監視用のコンピュータは、振動周波数スペクトラムデータに基づく振動監視用画像データを生成し、振動監視用画像データに基づく画像を表示してもよい。   Further, the vibration frequency spectrum data stored in the vibration data storage unit 26 is associated with the position information. Based on the data stored in the vibration data storage unit 26, the analysis of the vibration of the UUV body is performed for each UUV navigation position. For example, the sound wave monitoring computer may generate vibration monitoring image data based on the vibration frequency spectrum data, and display an image based on the vibration monitoring image data.

音波監視装置は、航走体として有人の潜水艦や船舶に搭載されてもよい。この場合、音波監視用画像データおよび振動監視用画像データを生成する画像生成部と、音波監視用画像データまたは振動監視用画像データに基づく画像を表示する表示部を音波監視装置が備える。音波監視装置が船舶に搭載される場合、測位部31はGPS測位装置等、測位衛星から信号を受信する装置であってもよい。   The sound wave monitoring device may be mounted on a manned submarine or ship as a vehicle. In this case, the sound wave monitoring device includes an image generation unit that generates the sound wave monitoring image data and the vibration monitoring image data, and a display unit that displays an image based on the sound wave monitoring image data or the vibration monitoring image data. When the sound wave monitoring device is mounted on a ship, the positioning unit 31 may be a device that receives a signal from a positioning satellite, such as a GPS positioning device.

図6には、画像生成部68および表示部70を備える音波監視装置が示されている。画像生成部68は、信号処理ユニット16から出力された積算周波数スペクトラムデータPα、位相同期周波数スペクトラムデータPβ、加算合計周波数スペクトラムデータPγ、分割・積算周波数スペクトラムデータPδ、および分割・加算合計周波数スペクトラムデータPεのうちいずれか1つ、または、いずれかの組み合わせに基づく音波監視用画像データを生成する。また、画像生成部68は、振動高速フーリエ変換部24から出力された振動周波数スペクトラムデータに基づく振動監視用画像データを生成する。画像生成部68は、音波監視用画像データおよび振動監視用画像データに対し、各画像データが取得された位置を示す情報を含ませる。画像生成部68は、音波監視用画像データおよび振動監視用画像データに基づいて表示用画像データを生成し表示部70に出力する。表示用画像データは、音波監視用画像データ、振動監視用画像データまたはこれらの組み合わせに基づく画像を表す画像データである。表示部70は、表示用画像データに基づく画像を表示する。また、表示部70は、表示用画像データに基づく画像と共に、その画像が取得された位置を示す情報を表示してもよい。   FIG. 6 shows a sound wave monitoring device including an image generation unit 68 and a display unit 70. The image generation unit 68 includes the integrated frequency spectrum data Pα, the phase-locked frequency spectrum data Pβ, the sum total frequency spectrum data Pγ, the division / summation frequency spectrum data Pδ, and the division / summation total frequency spectrum data output from the signal processing unit 16. The sound wave monitoring image data is generated based on any one of Pε or any combination. Further, the image generation unit 68 generates vibration monitoring image data based on the vibration frequency spectrum data output from the vibration fast Fourier transform unit 24. The image generation unit 68 causes the sound wave monitoring image data and the vibration monitoring image data to include information indicating a position where each image data is obtained. The image generation unit 68 generates display image data based on the sound wave monitoring image data and the vibration monitoring image data, and outputs the display image data to the display unit 70. The display image data is image data representing an image based on sound wave monitoring image data, vibration monitoring image data, or a combination thereof. The display unit 70 displays an image based on the display image data. Further, the display unit 70 may display an image based on the display image data and information indicating a position where the image is obtained.

このような構成によれば、航走体に搭乗するユーザは、航走体が航行しているときに表示部70を参照することで、航走体に到来する音波および航走体のボデーの振動を観測することができる。   According to such a configuration, the user boarding the navigation body refers to the display unit 70 when the navigation body is navigating, so that the sound waves arriving at the navigation body and the body of the navigation body can be displayed. Vibration can be observed.

上記では、n個の音響センサを用いた音波監視装置について説明した。1つの音響センサは、音響センサアレイに置き換えられてもよい。音響センサアレイは、複数の振動子と指向性制御部とを備える。指向性制御部は、各振動子から出力される信号の大きさおよび遅延時間の少なくとも一方を調整して合成する。これによって、受信指向性が音波の到来方向に応じて変化するセンサ部が構成される。   In the above, the sound wave monitoring device using n acoustic sensors has been described. One acoustic sensor may be replaced by an acoustic sensor array. The acoustic sensor array includes a plurality of transducers and a directivity control unit. The directivity control unit adjusts and synthesizes at least one of the magnitude and the delay time of the signal output from each transducer. Thus, a sensor unit whose reception directivity changes according to the arrival direction of the sound wave is configured.

なお、図1に示されているS/N改善処理部14、信号処理ユニット16、振動FFT24、制御部30、および図6に示されている画像生成部68は、プロセッサ等の演算処理デバイスによって構成してもよい。この場合、予め読み込まれたプログラムによって、各機能を実行してもよい。また、各機能はハードウエアによって個別に構成してもよい。   In addition, the S / N improvement processing unit 14, the signal processing unit 16, the vibration FFT 24, the control unit 30, and the image generation unit 68 illustrated in FIG. 6 illustrated in FIG. You may comprise. In this case, each function may be executed by a program read in advance. Further, each function may be individually configured by hardware.

10−1〜10−n 音響センサ、12−1〜12−n 音響受信部、14 S/N改善処理部、16 信号処理ユニット、16−1 第1信号処理部、16−2 第2信号処理部、16−3 第3信号処理部、18 音波データ記憶部、20 振動センサ、22 振動受信部、24 振動高速フーリエ変換部、26 振動データ記憶部、28 操作部、30 制御部、31 測位部、32 UUV、34 ボデー、36−1〜36−n,46,52−1〜52−n,60−1〜60−s 高速フーリエ変換部、38−1〜38−n,48,56,62−1〜62−s パワー演算部、40,64 乗算器、42 位相同期部、44 加算器、50 ピーク周波数検出部、54 複素加算部、58 分割部、66 加算器、68 画像生成部、70 表示部。
10-1 to 10-n sound sensor, 12-1 to 12-n sound receiving unit, 14 S / N improvement processing unit, 16 signal processing unit, 16-1 first signal processing unit, 16-2 second signal processing Unit, 16-3 third signal processing unit, 18 sound wave data storage unit, 20 vibration sensor, 22 vibration reception unit, 24 vibration fast Fourier transform unit, 26 vibration data storage unit, 28 operation unit, 30 control unit, 31 positioning unit , 32 UUV, 34 bodies, 36-1 to 36-n, 46, 52-1 to 52-n, 60-1 to 60-s Fast Fourier transform units, 38-1 to 38-n, 48, 56, 62 -1 to 62-s power operation unit, 40, 64 multiplier, 42 phase synchronization unit, 44 adder, 50 peak frequency detection unit, 54 complex addition unit, 58 division unit, 66 adder, 68 image generation unit, 70 Display section.

Claims (6)

それぞれが音波を検出する複数のセンサ部と、
複数の前記センサ部のそれぞれに対して設けられ、前記センサ部による検出信号に対して周波数変換処理を施す変換部と、
各前記変換部から出力された周波数スペクトラムデータに対する合成処理を実行する合成処理部と、
を備え、
前記合成処理部による合成処理によって得られた信号のピーク周波数成分を検出する周波数検出部と、
複数の前記センサ部による複数の検出信号のそれぞれに含まれる前記ピーク周波数成分の位相を揃える位相同期部と、
位相が揃えられた上で前記位相同期部から出力される複数の同期化信号を合成する同期合成処理部と、
前記同期合成処理部による合成信号に対して周波数変換処理を施す第2変換部と、
を備えることを特徴とする音波監視装置。
A plurality of sensor units each detecting a sound wave,
A conversion unit that is provided for each of the plurality of sensor units and performs a frequency conversion process on a detection signal from the sensor unit;
A synthesis processing unit that performs synthesis processing on the frequency spectrum data output from each of the conversion units,
With
A frequency detection unit that detects a peak frequency component of a signal obtained by the synthesis processing by the synthesis processing unit,
A phase synchronization unit that aligns the phase of the peak frequency component included in each of the plurality of detection signals by the plurality of sensor units,
A synchronization synthesis processing unit that synthesizes a plurality of synchronization signals output from the phase synchronization unit after the phases are aligned,
A second conversion unit that performs a frequency conversion process on a synthesized signal by the synchronous synthesis processing unit;
A sound wave monitoring device comprising:
請求項1に記載の音波監視装置において、
前記合成処理部は、
各前記変換部から出力された周波数スペクトラムデータに基づいて、各前記検出信号に含まれていた同一周波数成分を掛け合わせる乗算器を含むことを特徴とする音波監視装置。
The sound wave monitoring device according to claim 1,
The synthesis processing unit includes:
A sound wave monitoring apparatus comprising: a multiplier that multiplies the same frequency component included in each of the detection signals based on frequency spectrum data output from each of the conversion units.
請求項1または請求項2に記載の音波監視装置において、
前記同期合成処理部は、
複数の前記同期化信号を加算合計することを特徴とする音波監視装置。
In the sound wave monitoring device according to claim 1 or 2,
The synchronous synthesis processing unit,
A sound wave monitoring device, wherein a plurality of the synchronization signals are added and totalized.
請求項1から請求項のいずれか1項に記載の音波監視装置において、
各前記センサ部は、ピエゾフィルムセンサを備えることを特徴とする音波監視装置。
In the sound wave monitoring device according to any one of claims 1 to 3 ,
The sound wave monitoring device, wherein each of the sensor units includes a piezo film sensor.
請求項1から請求項のいずれか1項に記載の音波監視装置において、
各前記センサ部は、
複数の振動子と、
前記複数の振動子に基づく指向性を適応的に変化させる指向性制御部と、
を備えることを特徴とする音波監視装置。
In the sound wave monitoring device according to any one of claims 1 to 4 ,
Each of the sensor units,
A plurality of transducers,
A directivity control unit that adaptively changes the directivity based on the plurality of transducers,
A sound wave monitoring device comprising:
請求項1から請求項のいずれか1項に記載の音波監視装置を搭載した航走体。 A marine vehicle equipped with the sound wave monitoring device according to any one of claims 1 to 5 .
JP2016122319A 2016-06-21 2016-06-21 Sound wave monitoring device and airframe Active JP6654106B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2016122319A JP6654106B2 (en) 2016-06-21 2016-06-21 Sound wave monitoring device and airframe

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2016122319A JP6654106B2 (en) 2016-06-21 2016-06-21 Sound wave monitoring device and airframe

Publications (2)

Publication Number Publication Date
JP2017227480A JP2017227480A (en) 2017-12-28
JP6654106B2 true JP6654106B2 (en) 2020-02-26

Family

ID=60891638

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2016122319A Active JP6654106B2 (en) 2016-06-21 2016-06-21 Sound wave monitoring device and airframe

Country Status (1)

Country Link
JP (1) JP6654106B2 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111401120B (en) * 2019-11-19 2023-09-05 浙江机电职业技术学院 River course ship flow analysis system and method

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6010190A (en) * 1983-06-30 1985-01-19 Tech Res & Dev Inst Of Japan Def Agency Mutual signal processing system for passive sonar
US5515300A (en) * 1993-09-30 1996-05-07 The United States Of America As Represented By The Secretary Of The Navy Coherent signal power detector using higher-order statistics
JPH11347035A (en) * 1998-06-11 1999-12-21 Hitachi Medical Corp Ultrasonic diagnostic device
JP2001051039A (en) * 1999-08-09 2001-02-23 Nec Corp Signal processing system
JP2001349940A (en) * 2000-06-06 2001-12-21 Oki Electric Ind Co Ltd Pulse sound detecting system
FR2836076B1 (en) * 2002-02-15 2005-10-14 Thales Sa SURFACE ACOUSTIC ANTENNA FOR SUBMARINES

Also Published As

Publication number Publication date
JP2017227480A (en) 2017-12-28

Similar Documents

Publication Publication Date Title
JP6536910B2 (en) Target extraction system, target extraction method, information processing apparatus, and control method and control program therefor
US20160131759A1 (en) Ctfm detection apparatus and underwater detection apparatus
JP2008236077A (en) Target sound extracting apparatus, target sound extracting program
JP2017227515A (en) Active sonar and control method for the same
JP6654106B2 (en) Sound wave monitoring device and airframe
JP2010127771A (en) Synthetic aperture sonar, and method and program for correcting phase error of synthetic aperture sonar
JP2018091647A (en) Signal processing apparatus, orientation calculation method and orientation calculation program
KR20170124279A (en) Method and Apparatus for DEMON Processing in order that Removal of External Target Noise When measuring underwater radiated noise
JP7180447B2 (en) Azimuth Estimation Device, Azimuth Estimation System, Azimuth Estimation Method and Program
US8045419B1 (en) Method for mitigating spatial aliasing
CN107957571B (en) Hydrophone direction finding method and device, computer readable storage medium and computer equipment
US11129595B2 (en) Ultrasonic diagnostic apparatus, interpolation processing unit, and interpolation processing method
Rudnicki Implementation of DIFAR processing in ASW dipping sonar
Tan et al. Fast Fourier transform overlap approach for underwater acoustic positioning system
JP5603355B2 (en) Ultrasonic measuring device
JP5599200B2 (en) Target detection / tracking device
US8750077B2 (en) Acoustic ranging system for multi-line towed acoustic arrays
JP6610224B2 (en) Bistatic active sonar device and its receiver
US20220260711A1 (en) Target detection device and target detection method
Cheng et al. Noise cancellation for an autonomous underwater vehicle‐towed thin line array through recursive adaptive filtering
JP2014020906A (en) Signal processor, underwater searching device, signal processing method, and program
JP2019074373A (en) Marine engine speed estimating device, marine engine speed estimation method, and marine engine speed estimating program
JP2865090B2 (en) Passive sonar device
JP5713933B2 (en) Sound source distance measuring device, acoustic direct ratio estimating device, noise removing device, method and program thereof
JP7267825B2 (en) Time difference measuring device and direction of arrival estimating device

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20180301

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20181226

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20190122

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20190320

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20190702

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20190827

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20200114

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20200129

R150 Certificate of patent or registration of utility model

Ref document number: 6654106

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250