Nothing Special   »   [go: up one dir, main page]

WO2019215930A1 - 伝搬パス推定装置、伝搬パス推定方法及びプログラム - Google Patents

伝搬パス推定装置、伝搬パス推定方法及びプログラム Download PDF

Info

Publication number
WO2019215930A1
WO2019215930A1 PCT/JP2018/018414 JP2018018414W WO2019215930A1 WO 2019215930 A1 WO2019215930 A1 WO 2019215930A1 JP 2018018414 W JP2018018414 W JP 2018018414W WO 2019215930 A1 WO2019215930 A1 WO 2019215930A1
Authority
WO
WIPO (PCT)
Prior art keywords
propagation path
feature amount
solid
path estimation
sensors
Prior art date
Application number
PCT/JP2018/018414
Other languages
English (en)
French (fr)
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 US17/040,266 priority Critical patent/US20210010855A1/en
Priority to JP2020517749A priority patent/JP6984744B2/ja
Priority to PCT/JP2018/018414 priority patent/WO2019215930A1/ja
Publication of WO2019215930A1 publication Critical patent/WO2019215930A1/ja

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H3/00—Measuring characteristics of vibrations by using a detector in a fluid
    • G01H3/04—Frequency
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H3/00—Measuring characteristics of vibrations by using a detector in a fluid
    • G01H3/10—Amplitude; Power
    • G01H3/12—Amplitude; Power by electric means
    • G01H3/125—Amplitude; Power by electric means for representing acoustic field distribution
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H3/00—Measuring characteristics of vibrations by using a detector in a fluid
    • G01H3/10—Amplitude; Power
    • G01H3/12—Amplitude; Power by electric means

Definitions

  • the present invention relates to a propagation path estimation device, a propagation path estimation method, and a program.
  • Non-Patent Document 1 describes an example of a technique for estimating the direction of a transmission point of an acoustic signal received at a reception point.
  • a direction of a corresponding signal is estimated using a plurality of acoustic sensors and using a difference in arrival times to the plurality of sensors.
  • Patent Document 1 discloses an analysis method capable of accurately analyzing an acoustic transmission path in a system between a sound source and a sound receiving side.
  • the acoustic transfer analysis method includes a step of identifying all sound paths from a sound source to a sound receiving side, a step of measuring a transfer function from each path represented by a particle velocity vector to the sound receiving side, and each sound path. Measuring the operating load of the method. Further, the acoustic transfer analysis method includes a step of multiplying the operation load of the path by a transfer function from the path to the sound receiving side to represent the contribution of the path to the particle velocity vector on the sound receiving side. And determining a dominant sound path represented by a particle velocity vector from the calculated contribution of each path.
  • the frequency characteristics of the signal obtained at the receiving point may vary greatly. For example, when recording is performed at a later stage of the transmission path, it is considered that the sound quality can be improved and unnecessary sound can be removed by correcting the frequency characteristics according to the transmission path. In addition, for example, when detection is performed at a later stage of the transmission path, similarly, it is considered that the detection performance can be improved by correcting the frequency characteristics and switching the detection method according to the transmission path.
  • An object of the present invention is to provide a propagation path estimation device, a propagation path estimation method, and a program that can contribute to the enrichment of the estimation technique of a route through which a signal has been transmitted.
  • a propagation path estimation device According to a first aspect, a plurality of sensors, a feature amount calculation unit that calculates a feature amount from sensor output signals of the plurality of sensors, a propagation path determination unit that determines a propagation path corresponding to the feature amount, , A propagation path estimation device is provided.
  • a propagation path estimation method including a step of calculating a feature amount from sensor output signals of a plurality of sensors and a step of determining a propagation path corresponding to the feature amount.
  • This method is linked to a specific machine called a propagation path estimation device that determines a propagation path based on output signals from a plurality of sensors.
  • a program for causing a computer to execute a process for calculating a feature amount from sensor output signals of a plurality of sensors and a process for determining a propagation path corresponding to the corresponding feature amount can be recorded on a computer-readable (non-transitory) storage medium. That is, the present invention can be embodied as a computer program product.
  • FIG. 1st Embodiment of this invention It is a block diagram which shows the structure of the 1st Embodiment of this invention.
  • the 1st Embodiment of this invention it is a figure showing the relationship with the sound source used in the air as a propagation path.
  • the 1st Embodiment of this invention it is a figure showing the relationship with the sound source used as a propagation path in solid.
  • connection lines between blocks such as drawings referred to in the following description include both bidirectional and unidirectional directions.
  • the unidirectional arrow schematically shows the main signal (data) flow and does not exclude bidirectionality.
  • ports and interfaces at input / output connection points of each block in the figure they are not shown.
  • FIG. 1 is a block diagram showing the configuration of the first embodiment of the present invention.
  • a propagation path estimation apparatus 100 includes a plurality of sensors 101 and 102 (hereinafter also referred to as sensors 1 and 2), a feature amount calculation unit 201, a propagation path. And a determination unit 202.
  • the plurality of sensors 101 and 102 are fixed at regular intervals, and each receives information on the received signal.
  • An example of the fixing method is to arrange the sensors 101 and 102 along the wall as shown in FIG.
  • the outputs of the sensors 101 and 102 can be handled as digital time series signals.
  • the feature amount calculation unit 201 calculates a feature amount representing a spatial path required for arrival and transmission of a signal for each fixed time period from the signals received by the sensors 101 and 102.
  • the propagation path determination unit 202 determines whether the signal input to the sensors 101 and 102 has propagated through the air or the solid by using the feature amount, and the determination result. 401 is output.
  • the sound propagating through the air means that air is interposed between the sound source 300 and the sensors 101 and 102 as shown in FIG.
  • the acoustic path reaching the sensor 1101 from the sound source 300 in addition to the acoustic path 1-1 that is a direct sound, the acoustic path 1-2 that is a reflected sound, the acoustic path 1- 3.
  • There are a plurality of acoustic paths such as acoustic paths 1-4.
  • the senor 1 serving as a reception point is usually observed in a form in which sounds passing through these acoustic paths are mixed.
  • the acoustic path reaching the sensor 2 102 from the sound source 300 is also observed in the same manner.
  • the sound propagating through the solid is exemplified by a solid, in this case a wall, as shown in FIG. 3, but the sound source 300 is in contact with the solid.
  • a solid in this case a wall, as shown in FIG. 3, but the sound source 300 is in contact with the solid.
  • FIG. 6 only an acoustic path that is a direct sound exists as an acoustic path that reaches the sensor 1101 from the sound source 300, and there is no reflected sound.
  • the acoustic path reaching the sensor 2 102 from the sound source 300 is also observed in the same manner.
  • This cross spectrum itself is the reverse of the similarity of the propagation function of the path from the sound source 300 to the sensor 1 101 and the path to the sensor 2 102, that is, the difference. To express.
  • the cross-correlation function between the plurality of sensors 101 and 102 is obtained by inverse Fourier transforming this difference.
  • this cross-correlation function is output as a feature amount.
  • the propagation path determination unit 202 determines that the sound has propagated through the solid, and outputs the determination result 401.
  • the propagation path determination unit 202 determines that the sound has propagated through the air, and outputs the determination result 401.
  • the number of sensors has been described as two. However, the number of sensors may be three or more, and a determination may be made between them, and a determination may be made by taking a majority vote, a logical sum, or a logical product, thereby improving the estimation accuracy. It is possible to increase.
  • the propagation path determination unit 202 may operate only when the received signal has a certain level of power. As a result, it is also possible to reduce errors that occur under low power signals and thus under low S / N ratio conditions.
  • a plurality of sensors are exposed and fixed on the surface of the wall 301, which is solid, but may be installed behind the wall 301 or embedded in the wall 301 as shown in FIG. is there.
  • air holes 501 and 502 are preferably opened in order to operate the sensor appropriately.
  • microphones are used for the plurality of sensors 101 and 102.
  • the application range of the present invention is not limited to acoustic signals in the audible range, and ultrasonic microphones, ultra-low frequency sensors, and the like are used. May be.
  • the sound source 300 is in contact with the wall 301, but it is also conceivable that sound propagating from the sound source 300 in the air shakes the wall 301, and as a result, propagates through the wall 301.
  • sound life noise
  • the daily noise propagates through the air, but even when the windows are closed, it propagates through a solid such as a housing of a housing complex.
  • the propagation through the air and the wall corresponds to an example of the living noise of an apartment house, and it can be easily seen that the frequency characteristics are different from each other.
  • the propagation path is typically in the air and in the solid, but in the air is a typical path where there is reflection, and in the solid is a typical path where there is no reflection. It doesn't matter.
  • the air can be replaced with a gas such as nitrogen or a liquid such as water.
  • the solid can be replaced with a sufficiently viscous gel-like object.
  • the present embodiment since a signal transmission path is estimated using a plurality of sensors, it is possible to estimate a propagation path if there is information on only the reception point, and information on the transmission side or transmission path is not necessary. Absent.
  • the present embodiment has an advantage that the installation cost can be reduced because the range can be calculated in a normal operation without requiring special in-space calibration.
  • FIG. 9 is a block diagram showing the configuration of the second exemplary embodiment of the present invention.
  • a propagation path estimation apparatus 100A according to the second embodiment of the present invention includes a storage unit 203, a range determination unit 204, and a determination model 205 in addition to the configuration of the first embodiment. Since the other configuration is the same as that of the first embodiment, the difference will be mainly described below.
  • the accumulation unit 203 accumulates the feature amounts calculated in the past by the feature amount calculation unit 201 for a certain period.
  • the accumulation period of the feature amount in the accumulation unit 203 may be a long period of time since, for example, the installation of the plurality of sensors 101 and 102.
  • the feature amount of the past one day may be accumulated.
  • the storage unit 203 will be described assuming that the length of one frame is one second and 86400 frames for one day without overlap are stored.
  • the range determination unit 204 maps the feature amounts of all frames accumulated by the accumulation unit 203 to the feature amount space.
  • FIG. 10 shows an example of feature amount mapping by the range determination unit 204. In the example of FIG. 10, only 44 frames are shown, but the same applies even if the number of accumulated frames changes.
  • FIG. 10 is a two-dimensional histogram (also referred to as a heat map) using two cross-correlation functions as variables, which expresses how many frames correspond to the corresponding feature value by numerical values.
  • 37 frames have almost the same feature amount, and the remaining frames have different feature amounts.
  • the range indicated by the dotted circle in FIG. It can be determined that the range of the feature amount is a propagation path.
  • the other range is the feature amount range in which the propagation path is in the air.
  • a condition of taking all points exceeding a predetermined threshold D can be employed.
  • the feature-value which takes the maximum value instead of the conditions that the threshold value D was exceeded can be employed.
  • the radius ⁇ of the circle in the range can be determined in advance with a small value assuming a range affected by noise.
  • Judgment model 205 is obtained by storing the information of the range obtained in this way as a judgment model. Therefore, the above-described range determination unit 204 can be rephrased as a determination model creation unit.
  • the storage unit 203 described above can be configured using a storage device of a computer that configures the propagation path estimation device.
  • the judgment model 205 can also be stored in a storage device of a computer constituting the propagation path estimation device.
  • the propagation path determination unit 202 compares the value of the feature amount output from the feature amount calculation unit 201 with the range information stored in the determination model 205 to determine whether the feature amount is a propagation path in the air. It is determined whether it is an internal propagation path, and the corresponding determination result is output as a determination output 401.
  • the effect of this embodiment since it is possible to make a determination using past information, it is possible to make an optimal determination according to the installed environment and improve the accuracy.
  • the feature amount of the past day is accumulated and the determination model 205 is created.
  • the accumulated feature amount is stratified from various viewpoints, and a plurality of decision models are created. May be. For example, if the environment where the propagation path estimation device is installed changes depending on the time or season, create a judgment model using the time zone that includes the time and the feature values obtained in the relevant season. The determination can also be performed using the determination model.
  • propagation path estimation apparatus of the present invention it is possible to configure a recorder that performs appropriate processing according to the propagation path and performs easy-to-hear recording.
  • an acoustic event detector having high detection performance can be configured.
  • the procedures shown in the first and second embodiments described above can be realized by a program that causes a computer (9000 in FIG. 11) functioning as a propagation path estimation device to realize the function as the propagation path estimation device.
  • a computer is exemplified by a configuration including a CPU (Central Processing Unit) 9010, a communication interface 9020, a memory 9030, and an auxiliary storage device 9040 in FIG. That is, the CPU 9010 shown in FIG. 11 may execute a feature amount calculation program or a propagation path determination program to update each calculation parameter held in the auxiliary storage device 9040 or the like.
  • a CPU Central Processing Unit
  • each unit (processing means, function) of the propagation path estimation apparatus described in the first and second embodiments described above uses the hardware of the processor installed in the propagation path estimation apparatus. It can be realized by a computer program for executing processing.
  • the propagation path estimation apparatus described above is An aerial transmission path has multiple reflections, and a transmission path through a wall can be used to determine whether the propagation path is in the air or solid by utilizing only direct sound transmission.
  • the propagation path estimation apparatus described above is It is possible to adopt a configuration in which it is determined that the propagation path is a solid-state propagation path by utilizing the fact that the difference in the propagation path consists only of a term corresponding to a delay.
  • the propagation path estimation apparatus described above is It is possible to adopt a configuration in which the difference between the propagation paths is accumulated, a range for identifying whether the propagation path is in the air or in the solid is determined, and whether the propagation path is in the air or in the solid is determined based on the range. it can.
  • the propagation path estimation apparatus described above further includes: An accumulation unit for accumulating the feature amount calculated by the feature amount calculation unit; A range determination unit that creates a determination model for determining a range in which the solid is a propagation path among the accumulated feature amount; and
  • the propagation path determination unit may adopt a configuration in which a propagation path is determined using the feature amount calculated by the feature amount calculation unit and the determination model.
  • the propagation path estimation apparatus described above further includes: As the plurality of sensors, it is possible to adopt a configuration in which two or more of the same kind of a microphone, an ultrasonic microphone, and an ultra-low frequency sensor are connected. [Seventh form] (See the propagation path estimation method from the second viewpoint above) [Eighth form] (Refer to the program from the third viewpoint) Note that the seventh to eighth embodiments can be developed into the second to sixth embodiments in the same manner as the first embodiment.

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

音信号が伝達されてきた経路を推定する。伝搬パス推定装置は、音を検出する複数のセンサのセンサ出力信号から相互相関関数である特徴量を計算する特徴量算出部と、該当特徴量に対応して前記音が伝搬する媒質である伝搬パスを判断する伝搬パス判断部と、を含む。

Description

伝搬パス推定装置、伝搬パス推定方法及びプログラム
　本発明は、伝搬パス推定装置、伝搬パス推定方法及びプログラムに関する。
　受信点で受信された音響信号が近傍でどの経路を伝わって受信されたかを推定する推定装置が知られている。非特許文献１に、受信点において受信された音響信号の、発信点の方向を推定する技術の例が記載されている。この文献記載の方法では、複数の音響センサを用いて、複数センサへの到来時間の差を用いて、該当信号の方向を推定する。
　また、特許文献１に音源と受音側との間の系における音響伝達経路の正確な解析ができるという解析方法が開示されている。この音響伝達解析方法は、音源から受音側までの全ての音経路を識別する段階と、粒子速度ベクトルで表した各経路から前記受音側までの伝達関数を測定する段階と、各音経路の動作負荷を測定する段階と、を含む。さらに、この音響伝達解析方法は、前記経路の前記動作負荷に前記経路から前記受音側までの伝達関数を掛けて、前記受音側での前記粒子速度ベクトルへの前記経路の寄与を表わす段階と、前記計算した各経路の寄与から、粒子速度ベクトルで表した支配的音経路を決定する段階と、を含む。
特開２０１３－７９９５３号公報
Masanori Kato, Yuzo Senda and Reishi Kondo, "TDOA Estimation Based on Phase-Voting Cross Correlation and Circular Standard Deviation," AASP-P1, EUSIPCO 2017 (2017.08)
　以下の分析は、本発明によって与えられたものである。伝達経路に依って受信点で得られる信号の周波数特性が大きく変わる場合がある。例えば、伝達経路の後段で録音を行う場合に、伝達経路に応じた周波数特性の補正を行うことで、音質の向上や不要音の除去を行なうことができると考えられる。また例えば、伝達経路の後段で検知を行う場合も同様に、伝達経路に応じた周波数特性の補正や検知方式の切り替えを行うことで、検知性能を向上することができると考えられる。
　本発明の目的は、信号が伝達されてきた経路の推定技術の豊富化に貢献できる伝搬パス推定装置、伝搬パス推定方法及びプログラムを提供することにある。
　第１の視点によれば、複数のセンサと、前記複数のセンサのセンサ出力信号から特徴量を計算する特徴量算出部と、該当特徴量に対応して伝搬パスを判断する伝搬パス判断部と、を有する伝搬パス推定装置が提供される。
　第２の視点によれば、複数のセンサのセンサ出力信号から特徴量を計算するステップと、該当特徴量に対応して伝搬パスを判断するステップと、を含む、伝搬パス推定方法が提供される。本方法は、複数のセンサからの出力信号に基づいて伝搬パスを判断する伝搬パス推定装置という、特定の機械に結びつけられている。
　第３の視点によれば、コンピュータに、複数のセンサのセンサ出力信号から特徴量を計算する処理と、該当特徴量に対応して伝搬パスを判断する処理と、を実行させるプログラムが提供される。なお、このプログラムは、コンピュータが読み取り可能な（非トランジトリーな）記憶媒体に記録することができる。即ち、本発明は、コンピュータプログラム製品として具現することも可能である。
　本発明によれば、信号が伝達されてきた経路を推定することが可能となる。
本発明の第１の実施形態の構成を示すブロック図である。 本発明の第１の実施形態において、空気中が伝搬経路となる音源との関係を表す図である。 本発明の第１の実施形態において、固体中が伝搬経路となる音源との関係を表す図である。 本発明の第１の実施形態において、空気中が伝搬経路となる音源から、センサ１に到達する音の経路を説明する図である。 本発明の第１の実施形態において、空気中が伝搬経路となる音源から、センサ２に到達する音の経路を説明する図である。 本発明の第１の実施形態において、固体中が伝搬経路となる音源から、センサ１に到達する音の経路を説明する図である。 本発明の第１の実施形態において、固体中が伝搬経路となる音源から、センサ２に到達する音の経路を説明する図である。 本発明の第１の実施形態において、別のセンサの配置方法の例を説明する図である。 本発明の第２の実施形態を実施するための形態を示すブロック図である。 本発明の第２の実施形態の動作を説明する図である。 本発明の伝搬パス推定装置を構成するコンピュータの構成を示す図である。
［第１の実施の形態］
　続いて、本発明の第１の実施形態について図面を参照して詳細に説明する。なお、この概要に付記した図面参照符号は、理解を助けるための一例として各要素に便宜上付記したものであり、本発明を図示の態様に限定することを意図するものではない。また、以降の説明で参照する図面等のブロック間の接続線は、双方向及び単方向の双方を含む。一方向矢印については、主たる信号（データ）の流れを模式的に示すものであり、双方向性を排除するものではない。また、図中の各ブロックの入出力の接続点には、ポート乃至インタフェースがあるが図示省略する。
　図１は、本発明の第１の実施形態の構成を示すブロック図である。図１を参照すると、本発明の第１の実施形態の伝搬パス推定装置１００は、複数のセンサ１０１と１０２（以降、センサ１、センサ２とも記す）と、特徴量算出部２０１と、伝搬パス判断部２０２とを含む。
　これらの手段はそれぞれ概略すると次のように動作する。複数のセンサ１０１と１０２は、一定の間隔をあけて固定され、それぞれ受信した信号の情報を得る。固定の方法の一例は、図２に示すように壁に沿ってセンサ１０１、１０２を配置するものである。本実施形態では、センサ１０１、１０２の出力はディジタル量の時系列信号として扱うことができるものとする。
　特徴量算出部２０１は、センサ１０１、１０２で受信した信号から、一定時間周期ごとに、信号の到来、伝達に要した空間パスを表現する特徴量を算出する。
　伝搬パス判断部２０２は、前記特徴量を用いて、センサ１０１、１０２に入力されている信号が、空気中を伝搬してきたものか、あるいは固体中を伝搬してきたものかを判定し、判断結果４０１を出力する。
　次に、図１～８を参照して、本実施形態の全体の動作について詳細に説明する。はじめに、本実施形態の伝搬パス推定装置が識別する対象とする２つの事象、空気中を伝搬する音と固体中を伝搬する音について、図２から図７を用いて説明する。空気中を伝搬してくる音とは、図２に示すように、音源３００とセンサ１０１およびセンサ１０２の間に空気が介在しているものである。この場合、図４に示すように、音源３００からセンサ１　１０１に到達する音響パスとしては直接音である音響パス１－１の他に、反射音である音響パス１－２、音響パス１－３、音響パス１－４など複数の音響パスが存在する。この場合、受信点となるセンサ１ではこれらの音響パスを通った音が混じった形で観測されることが通例である。同様に、図５に示すように、音源３００からセンサ２　１０２に到達する音響パスも、同様に観測される。
　これに対して、固体を伝搬してくる音とは、図３に示すように固体、この場合は壁を例としているが、この固体に接して音源３００が存在する。この場合、図６に示すように、音源３００からセンサ１　１０１に到達する音響パスとしては直接音である音響パスだけが存在し、反射音は存在しない。同様に、図７に示すように、音源３００からセンサ２　１０２に到達する音響パスも、同様に観測される。
　ここではセンサ１　１０１とセンサ２　１０２には、マイクロフォンを用いる。特徴量算出部２０１は、センサ１０１とセンサ１０２から入力される信号のクロススペクトルを逐次計算する。すなわち、センサ１０１の信号系列ｘ１（ｔ）とセンサ１０２の信号系列ｘ２（ｔ）について、それぞれのフーリエ変換をＸ１（ｆ），Ｘ２（ｆ）と置き、Ｘ２（ｆ）の複素共役をＸ２＊（ｆ）とすれば、ある時刻におけるクロススペクトルＷ（ｆ）は、Ｗ（ｆ）＝Ｘ１（ｆ）　Ｘ２＊（ｆ）として計算することができる。
　このクロススペクトルそのもの、あるいはクロススペクトルの形状を適切な形状のフィルタで切り出したものが、音源３００からセンサ１　１０１へのパスとセンサ２　１０２へのパスの伝搬関数の類似度の逆、すなわち差異を表す。
　ここでクロススペクトルを計算する際にノルム正規化を行うことにより、音の大きさへの依存性を除去することも可能である。
　この差異を逆フーリエ変換することで、複数のセンサ１０１、１０２の間の相互相関関数を得る。ここでは、この相互相関関数を特徴量として出力する。
　次に、伝搬パス判断部２０２の動作について説明する。もしも特徴量算出部２０１が生成した相互相関関数が単一のピークを持つ場合は、複数のセンサ１０１、１０２の間に時間遅れの関係しか存在しないことは明らかである。この場合、反射波による影響は存在しないので、伝搬パス判断部２０２は固体中を伝搬してきた音であると判断し、判断結果４０１として出力する。
　一方、特徴量算出部２０１が生成した相互相関関数が複数のピークを持つ場合は、複数のセンサ１０１、１０２の間に時間遅れ以外の関係が存在することから、反射波による影響が存在するので、伝搬パス判断部２０２は空気中を伝搬してきた音であると判断し、判断結果４０１として出力する。
　ここではセンサ数を２であるとして説明したが、センサ数を３以上とし、それぞれの間で判断を行い、多数決または論理和あるいは論理積を取って判断を行っても良く、それによって推定精度を高めることが可能である。
　また、受信した信号がある一定以上のパワーを持っている時だけ、伝搬パス判断部２０２が動作しても良い。それによってパワーの小さい信号ひいてはＳ／Ｎ比の小さい条件下で発生する誤りを低減することも可能となる。
　なお、図２では複数のセンサを固体である壁３０１の表面に露出して固定して設置したが、図８に示すように壁３０１の裏側、あるいは壁３０１に埋め込んで設置しても同様である。この場合、センサを適切に動作させるためにそれぞれ空気穴５０１、５０２を空けると良い。
　また、本実施形態では、複数のセンサ１０１、１０２にマイクロフォンを用いたが、本発明の適用範囲は、可聴域の音響信号に限定する必要はなく、超音波マイクロフォンや超低周波センサなどを用いても良い。
　また、図３では音源３００は壁３０１に接しているとしたが、音源３００から空気中を伝搬してきた音が壁３０１を震わせ、その結果壁３０１を伝搬することも考えられる。例えば、集合住宅で隣家や階の異なる家の中で発する音（生活騒音）が聞こえる場合には次のようなケースがある。例えば、窓が全て開いている場合、当該生活騒音は、空気中を伝搬して届くが、窓を閉めていても集合住宅の躯体などの固体を伝搬して届く。上記空気と壁を経由する伝搬は、集合住宅の生活騒音の事例に相当し、それぞれで、周波数特性が異なることも容易に分かる。
　ここでは典型的に伝搬経路を空気中と固体中としたが、空気中は反射の存在する典型的な経路であり、固体中は反射の無い典型的な経路である範囲で、他の媒質としてもかまわない。例えば、空気中は、窒素などの気体、水などの液体で代えることができる。また固体中は、十分に粘性の高いゲル状の物体などで代えることができる。
　次に、本実施形態の効果について説明する。本実施形態によれば、複数のセンサを用いて信号の伝達するパスを推定するため、受信点だけの情報があれば伝搬パスの推定が可能であり、送信側あるいは伝達経路に関する情報は必要ではない。換言すると、本実施形態は、特別な空間内キャリブレーションを必要とせず通常の運用で範囲が計算できるため、設置コストを低減することができるという利点もある。
［第２の実施形態］
　続いて、より判定精度を向上できる第２の実施形態について図面を参照して詳細に説明する。図９は、本発明の第２の実施形態の構成を示すブロック図である。図９を参照すると、本発明の第２の実施形態の伝搬パス推定装置１００Ａは、第１の実施形態の構成に加えて、蓄積部２０３、範囲決定部２０４及び判断モデル２０５を有する。その他の構成は第１の実施形態と同様であるので、以下、その相違点を中心に説明する。
　蓄積部２０３は、特徴量算出部２０１が過去において算出した特徴量を一定期間蓄積する。蓄積部２０３における特徴量の蓄積期間は、例えば複数のセンサ１０１、１０２を設置して以来といった長期間のものでも良いが、典型的には例えば過去１日間の特徴量を蓄積するものでもよい。以下、本実施形態では、蓄積部２０３は、１フレームの長さを１秒間とし、重なりの無い１日分の８６４００フレームを蓄積するものとして説明する。
　範囲決定部２０４は、蓄積部２０３が蓄積した全フレームの特徴量を、特徴量空間にマップする。図１０は、範囲決定部２０４による特徴量のマッピングの例を示す。図１０の例では、４４フレームのみを記載しているが、蓄積されたフレーム数が変わっても同様である。
　ここで、図１０は、該当特徴量に幾つのフレームが該当するかを数値で表現した、２つの相互相関関数を変数とした２次元ヒストグラム（ヒートマップとも言う。）である。図１０の例では、３７フレームがほぼ同一の特徴量を持ち、残るフレームはそれぞれ異なる特徴量を持つ。ここで、多数のフレームが同一の特徴量を持つということは、特徴量の分散が小さく遅延項のみで構成されていることであるから、図１０で点線の円で示した範囲を、固体中を伝搬経路とする特徴量の範囲であると決定することができる。一方、その他の領域は特徴量の分散が大きいということであるから、それ以外の範囲を空気中を伝搬経路とする特徴量の範囲であると決定することができる。ここで多数のフレームである条件としては、予め定められた閾値Ｄを超えた点全てを取るという条件を採用することができる。もちろん、閾値Ｄを超えたという条件の代わりに最大値を取る特徴量としても良い。また範囲の円の半径εは、雑音により影響を受ける範囲を想定して小さな値を予め決めることができる。
　判断モデル２０５は、このようにして得られた範囲の情報を判断モデルとして保存したものである。従って、上記した範囲決定部２０４は、判断モデル作成部と言い換えることもできる。
　なお、上記した蓄積部２０３は、伝搬パス推定装置を構成するコンピュータの記憶装置を用いて構成することができる。同様に、判断モデル２０５も、伝搬パス推定装置を構成するコンピュータの記憶装置に保存させることができる。
　伝搬パス判断部２０２は、特徴量算出部２０１が出力した特徴量の値を、前記判断モデル２０５に保存された範囲の情報と比較して、該当特徴量が空気中の伝搬経路であるか固体中の伝搬経路であるかを決定し、該当判断結果を判断出力４０１として出力する。
　次に、本実施形態の効果について説明する。本実施形態では、過去の情報を用いて判断を行うことができるため、設置した環境に応じて最適な判定を行い、その精度を向上させることが可能である。なお、上記した実施形態では、過去１日間の特徴量を蓄積し、判断モデル２０５を作成するものとして説明したが、蓄積した特徴量をさまざまな観点で層別し、複数の判断モデルを作成してもよい。例えば、伝搬パス推定装置が設置される環境が時刻や季節によって環境が変わるというような場合には、当該時刻を含む時間帯、該当する季節に得られた特徴量を用いて判断モデルを作成し、その判断モデルを用いて判定を行うこともできる。
　また、本発明の伝搬パス推定装置を組み込むことで、伝搬経路に応じた適切な処理を行い、聞き取りやすい録音を行う録音機を構成することが可能となる。また、本発明の伝搬パス推定装置を組み込むことで高い検知性能を持つ音響イベント検知器を構成することも可能となる。
　以上、本発明の各実施形態を説明したが、本発明は、上記した実施形態に限定されるものではなく、本発明の基本的技術的思想を逸脱しない範囲で、更なる変形・置換・調整を加えることができる。例えば、各図面に示した装置構成、各要素の構成、データの表現形態は、本発明の理解を助けるための一例であり、これらの図面に示した構成に限定されるものではない。また、以下の説明において、「Ａ及び／又はＢ」は、Ａ及びＢの少なくともいずれかという意味で用いる。
　また、上記した第１～第２の実施形態に示した手順は、伝搬パス推定装置として機能するコンピュータ（図１１の９０００）に、伝搬パス推定装置としての機能を実現させるプログラムにより実現可能である。このようなコンピュータは、図１１のＣＰＵ（Ｃｅｎｔｒａｌ　Ｐｒｏｃｅｓｓｉｎｇ　Ｕｎｉｔ）９０１０、通信インタフェース９０２０、メモリ９０３０、補助記憶装置９０４０を備える構成に例示される。すなわち、図１１のＣＰＵ９０１０にて、特徴量算出プログラムや伝搬パス判断プログラムを実行し、その補助記憶装置９０４０等に保持された各計算パラメーターの更新処理を実施させればよい。
　即ち、上記した第１、第２の実施形態に示した伝搬パス推定装置の各部（処理手段、機能）は、伝搬パス推定装置に搭載されたプロセッサに、そのハードウェアを用いて、上記した各処理を実行させるコンピュータプログラムにより実現することができる。
　最後に、本発明の好ましい形態を要約する。
［第１の形態］
　（上記第１の視点による伝搬パス推定装置参照）
［第２の形態］
　上記した伝搬パス推定装置は、
　空中伝達パスは多重反射が有り、壁越し伝達パスは直接音のみが伝達することを利用して、前記伝搬パスが空気中か固体中かを求める構成を採ることができる。
［第３の形態］
　上記した伝搬パス推定装置は、
　前記伝搬パスの差異が遅延に相当する項だけからなることを利用して固体中の伝搬経路であると判断する構成を採ることができる。
［第４の形態］
　上記した伝搬パス推定装置は、
　前記伝搬パスの差異を蓄積し、前記伝搬パスが空気中か固体中かを識別する範囲を決定し、その範囲に基づいて前記伝搬パスが空気中か固体中かを判断する構成を採ることができる。
［第５の形態］
　上記した伝搬パス推定装置は、さらに、
　前記特徴量算出部により計算された特徴量を蓄積する蓄積部と、
　前記蓄積した特徴量のうち、固体を伝搬パスとする範囲を定める判断モデルを作成する範囲決定部と、を備え、
　前記伝搬パス判断部は、前記特徴量算出部により計算された特徴量と、前記判断モデルとを用いて伝搬パスを判断する構成を採ることができる。
［第６の形態］
　上記した伝搬パス推定装置は、さらに、
　前記複数のセンサとして、マイクロフォン、超音波マイクロフォン及び超低周波センサのうちのいずれか同種のもの２以上が接続される構成を採ることができる。
［第７の形態］
　（上記第２の視点による伝搬パス推定方法参照）
［第８の形態］
　（上記第３の視点によるプログラム参照）
　なお、上記第７～第８の形態は、第１の形態と同様に、第２～第６の形態に展開することが可能である。
　なお、上記の特許文献および非特許文献の各開示を、本書に引用をもって繰り込むものとする。本発明の全開示（請求の範囲を含む）の枠内において、さらにその基本的技術思想に基づいて、実施形態ないし実施例の変更・調整が可能である。また、本発明の開示の枠内において種々の開示要素（各請求項の各要素、各実施形態ないし実施例の各要素、各図面の各要素等を含む）の多様な組み合わせ、ないし選択（又は削除）が可能である。すなわち、本発明は、請求の範囲を含む全開示、技術的思想にしたがって当業者であればなし得るであろう各種変形、修正を含むことは勿論である。特に、本書に記載した数値範囲については、当該範囲内に含まれる任意の数値ないし小範囲が、別段の記載のない場合でも具体的に記載されているものと解釈されるべきである。
　１００、１００Ａ　伝搬パス推定装置
　１０１、１０２　センサ
　２０１　特徴量算出部
　２０２　伝搬パス判断部
　２０３　蓄積部
　２０４　範囲決定部
　２０５　判断モデル
　３００　音源
　３０１　壁
　４０１　判断結果
　９０００　コンピュータ
　９０１０　ＣＰＵ
　９０２０　通信インタフェース
　９０３０　メモリ
　９０４０　補助記憶装置

Claims (10)

  1. 　複数のセンサのセンサ出力信号から特徴量を計算する特徴量算出部と、
    　該当特徴量に対応して伝搬パスを判断する伝搬パス判断部と、
    　を含むことを特徴とする、伝搬パス推定装置。
  2. 　空中伝達パスは多重反射が有り、壁越し伝達パスは直接音のみが伝達することを利用して、前記伝搬パスが空気中か固体中かを求める、請求項１記載の伝搬パス推定装置。
  3. 　前記伝搬パスの差異が遅延に相当する項だけからなることを利用して固体中の伝搬経路であると判断する、請求項２記載の伝搬パス推定装置。
  4. 　前記伝搬パスの差異を蓄積し、前記伝搬パスが空気中か固体中かを識別する範囲を決定し、その範囲に基づいて前記伝搬パスが空気中か固体中かを判断する、請求項２記載の伝搬パス推定装置。
  5. 　さらに、
    　前記特徴量算出部により計算された特徴量を蓄積する蓄積部と、
    　前記蓄積した特徴量のうち、固体を伝搬パスとする範囲を定める判断モデルを作成する範囲決定部と、を備え、
    　前記伝搬パス判断部は、前記特徴量算出部により計算された特徴量と、前記判断モデルとを用いて伝搬パスを判断する請求項１から４いずれか一に記載の伝搬パス推定装置。
  6. 　さらに、前記複数のセンサとして、マイクロフォン、超音波マイクロフォン及び超低周波センサのうちのいずれか同種のもの２以上備える請求項１から５いずれか一に記載の伝搬パス推定装置。
  7. 　複数のセンサのセンサ出力信号から特徴量を計算するステップと、
    　該当特徴量に対応して伝搬パスを判断するステップと、
    　を含む、伝搬パス推定方法。
  8. 　空中伝達パスは多重反射が有り、壁越し伝達パスは直接音のみが伝達することを利用して、前記伝搬パスが空気中か固体中かを求める、請求項７記載の伝搬パス推定方法。
  9. 　伝搬パスの差異を蓄積し、前記伝搬パスが空気中か固体中かを識別する範囲を決定し、その範囲に基づいて前記伝搬パスが空気中か固体中かを判断する、請求項７記載の伝搬パス推定方法。
  10. 　コンピュータに、
    　複数のセンサのセンサ出力信号から特徴量を計算する処理と、
    　該当特徴量に対応して伝搬パスを判断する処理と、
    　を実行させるプログラム。
PCT/JP2018/018414 2018-05-11 2018-05-11 伝搬パス推定装置、伝搬パス推定方法及びプログラム WO2019215930A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
US17/040,266 US20210010855A1 (en) 2018-05-11 2018-05-11 Propagation path estimation apparatus, method, and program
JP2020517749A JP6984744B2 (ja) 2018-05-11 2018-05-11 伝搬パス推定装置、伝搬パス推定方法及びプログラム
PCT/JP2018/018414 WO2019215930A1 (ja) 2018-05-11 2018-05-11 伝搬パス推定装置、伝搬パス推定方法及びプログラム

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2018/018414 WO2019215930A1 (ja) 2018-05-11 2018-05-11 伝搬パス推定装置、伝搬パス推定方法及びプログラム

Publications (1)

Publication Number Publication Date
WO2019215930A1 true WO2019215930A1 (ja) 2019-11-14

Family

ID=68467371

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2018/018414 WO2019215930A1 (ja) 2018-05-11 2018-05-11 伝搬パス推定装置、伝搬パス推定方法及びプログラム

Country Status (3)

Country Link
US (1) US20210010855A1 (ja)
JP (1) JP6984744B2 (ja)
WO (1) WO2019215930A1 (ja)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US12044810B2 (en) 2021-12-28 2024-07-23 Samsung Electronics Co., Ltd. On-device user presence detection using low power acoustics in the presence of multi-path sound propagation

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007057460A (ja) * 2005-08-26 2007-03-08 Honda Motor Co Ltd 振動・音圧伝達特性解析装置及び方法
JP2009186481A (ja) * 2008-02-08 2009-08-20 Avl List Gmbh 伝達経路毎の成分を算定するための方法
JP2012021969A (ja) * 2010-06-18 2012-02-02 Panasonic Electric Works Co Ltd 部屋、配線器具
WO2017025248A1 (en) * 2015-08-10 2017-02-16 Harman Becker Automotive Systems Gmbh Noise and vibration sensing

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE3015323A1 (de) * 1979-04-23 1980-11-06 Antoni Marjan Ziolkowski Verfahren und vorrichtung zur lagebestimmung unterirdischer grenzflaechen
DE10064756A1 (de) * 2000-12-22 2002-07-04 Daimler Chrysler Ag Verfahren und Anordnung zur Verarbeitung von Geräuschsignalen einer Geräuschquelle
US8050142B2 (en) * 2007-12-06 2011-11-01 Sanyo Electric Co., Ltd. Sound collection environment deciding device, sound processing device, electronic appliance, sound collection environment deciding method and sound processing method
JP6012621B2 (ja) * 2010-12-15 2016-10-25 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． リモートノイズ検知器を使用したノイズ削減システム
US20130083929A1 (en) * 2011-09-30 2013-04-04 Hitachi, Ltd. Method for analyzing sound transmission paths in a system
US20150362465A1 (en) * 2014-06-13 2015-12-17 Reece Innovation Centre Limited System for monitoring and/or surveying conduits
GB2545441B (en) * 2015-12-15 2022-09-14 Water Intelligence Int Ltd System for monitoring and/or surveying conduits
WO2017125613A1 (de) * 2016-01-24 2017-07-27 Sensaction Ag Verfahren zur akustischen bestimmung von eigenschaften eines mediums und vorrichtung zur akustischen bestimmung von eigenschaften eines mediums mit reflexionselement
US10171638B2 (en) * 2016-02-01 2019-01-01 The Regents Of The University Of Michigan Force sensing based on structure-borne sound propagation
US10598543B1 (en) * 2017-12-04 2020-03-24 Amazon Technologies, Inc. Multi microphone wall detection and location estimation

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007057460A (ja) * 2005-08-26 2007-03-08 Honda Motor Co Ltd 振動・音圧伝達特性解析装置及び方法
JP2009186481A (ja) * 2008-02-08 2009-08-20 Avl List Gmbh 伝達経路毎の成分を算定するための方法
JP2012021969A (ja) * 2010-06-18 2012-02-02 Panasonic Electric Works Co Ltd 部屋、配線器具
WO2017025248A1 (en) * 2015-08-10 2017-02-16 Harman Becker Automotive Systems Gmbh Noise and vibration sensing

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
KATO, MASANORI ET AL.: "TDOA Estimation Based on Phase-Voting Cross Correlation and Circular Standard Deviation", 2017 25TH EUROPEAN SIGNAL PROCESSING CONFERENCE, August 2017 (2017-08-01), pages 1230 - 1234, XP033236133, DOI: 10.23919/EUSIPCO.2017.8081404 *

Also Published As

Publication number Publication date
JP6984744B2 (ja) 2021-12-22
JPWO2019215930A1 (ja) 2021-05-27
US20210010855A1 (en) 2021-01-14

Similar Documents

Publication Publication Date Title
KR102340999B1 (ko) 시간 지연 추정을 기반으로 하는 에코 제거 방법 및 장치
KR101910679B1 (ko) 마이크로폰 어레이를 위한 잡음 적응형 빔 형성 기법
KR101529516B1 (ko) 음원 위치 추정 장치 및 음원 위치 추정 방법
KR100899836B1 (ko) 실내 충격응답 모델링 방법 및 장치
KR101562904B1 (ko) 도래각 측정 장치 및 방법
EP3078210B1 (en) Estimating a room impulse response for acoustic echo cancelling
CN203811527U (zh) 用于声学应用的检测装置
US20160171965A1 (en) Vibration source estimation device, vibration source estimation method, and vibration source estimation program
Saqib et al. Estimation of acoustic echoes using expectation-maximization methods
KR101086304B1 (ko) 로봇 플랫폼에 의해 발생한 반사파 제거 신호처리 장치 및 방법
Canclini et al. Exact localization of acoustic reflectors from quadratic constraints
WO2019215930A1 (ja) 伝搬パス推定装置、伝搬パス推定方法及びプログラム
KR101483513B1 (ko) 음원위치추적장치 및 음원위치추적방법
Rathsam et al. Analysis of absorption in situ with a spherical microphone array
Padois et al. On the use of geometric and harmonic means with the generalized cross-correlation in the time domain to improve noise source maps
Chen et al. Bayesian estimations of dissipation, sound speed, and microphone positions in impedance tubes
JP5946039B2 (ja) 騒音監視システム
Wen et al. Robust time delay estimation for speech signals using information theory: A comparison study
Noël et al. A new temporal method for the identification of source directions in a reverberant hall
CN101645701A (zh) 一种基于滤波器组的时延估计方法及其系统
JP7211419B2 (ja) パターン認識装置、パターン認識方法及びパターン認識プログラム
US11835625B2 (en) Acoustic-environment mismatch and proximity detection with a novel set of acoustic relative features and adaptive filtering
Vieira et al. A sound localizer robust to reverberation
Fernandez et al. Indoor simulation of amplitude modulated wind turbine noise
Dankers et al. Blind acoustic source separation via system identification for leak detection in pipelines

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 18917816

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2020517749

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 18917816

Country of ref document: EP

Kind code of ref document: A1