JP2021183910A - Estimation method of acoustic power level - Google Patents

Estimation method of acoustic power level Download PDF

Info

Publication number
JP2021183910A
JP2021183910A JP2020088613A JP2020088613A JP2021183910A JP 2021183910 A JP2021183910 A JP 2021183910A JP 2020088613 A JP2020088613 A JP 2020088613A JP 2020088613 A JP2020088613 A JP 2020088613A JP 2021183910 A JP2021183910 A JP 2021183910A
Authority
JP
Japan
Prior art keywords
sound source
estimation
power level
sound
value
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2020088613A
Other languages
Japanese (ja)
Other versions
JP7387530B2 (en
Inventor
幹也 伊藤
Mikiya Ito
圭司 三神
Keiji Mikami
卓美 小曽戸
Takumi Kosodo
昌之 小西
Masayuki Konishi
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Ono Sokki Co Ltd
Original Assignee
Ono Sokki Co Ltd
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 Ono Sokki Co Ltd filed Critical Ono Sokki Co Ltd
Priority to JP2020088613A priority Critical patent/JP7387530B2/en
Publication of JP2021183910A publication Critical patent/JP2021183910A/en
Application granted granted Critical
Publication of JP7387530B2 publication Critical patent/JP7387530B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

To provide an estimation method of an acoustic power level that can easily calculate a correction value to be used when estimating the acoustic power level of an estimation target sound source.SOLUTION: An estimation method of an acoustic power level comprises: a positioning step (step S102) for determining setting positions of a plurality of microphones; a reproduction step (step S107) for reproducing a measurement environment of an estimation target sound source in a semi-anechoic space by using a reference sound source whose acoustic power level is known; a reproduction acquisition step (step S108) for acquiring a sound pressure level of the reference sound source that the plurality of microphones measure in the reproduced measurement environment; and a first calculation step (step S110) for calculating a correction value to be used when estimating the acoustic power level of the estimation target sound source based on the acquired sound pressure level.SELECTED DRAWING: Figure 3

Description

本発明は、音響パワーレベルの推定方法に関する。 The present invention relates to a method for estimating sound power level.

閉空間における音源の音響パワーレベルを推定する方法として、例えば特許文献1のように、音響パワーレベルが既知である基準音源を用いた方法が知られている。特許文献1に記載の方法では、推定対象音源が設置される閉空間内に複数のセンサを設置し、それらのセンサで計測した音圧レベルの平均値である第1の平均音圧レベルに補正値を加えることにより音響パワーレベルを算出している。 As a method of estimating the sound power level of a sound source in a closed space, for example, as in Patent Document 1, a method using a reference sound source whose sound power level is known is known. In the method described in Patent Document 1, a plurality of sensors are installed in a closed space where the sound source to be estimated is installed, and the sound pressure level is corrected to the first average sound pressure level which is the average value of the sound pressure levels measured by those sensors. The sound power level is calculated by adding a value.

補正値は、上記閉空間内に設定した複数の測定点に基準音源を順次設置し、各測定点で基準音源を作動させた場合に各センサが計測した音圧レベルに基づいて算出されている。具体的には、基準音源について各センサが計測した音圧レベルの平均値である第2の平均音圧レベルを測定点ごとに算出する。そして、それら第2の平均音圧レベルの最大値と最小値との中間値を算出し、その中間値と基準音源の音響パワーレベルとの差が補正値に設定されている。 The correction value is calculated based on the sound pressure level measured by each sensor when the reference sound source is sequentially installed at a plurality of measurement points set in the closed space and the reference sound source is operated at each measurement point. .. Specifically, the second average sound pressure level, which is the average value of the sound pressure levels measured by each sensor for the reference sound source, is calculated for each measurement point. Then, an intermediate value between the maximum value and the minimum value of the second average sound pressure level is calculated, and the difference between the intermediate value and the sound power level of the reference sound source is set as the correction value.

特開2002−365127号公報Japanese Unexamined Patent Publication No. 2002-365127

しかしながら、特許文献1に記載の方法では、閉空間が狭い場合に、その狭い閉空間への基準音源の設置とセンサによる計測とを繰り返し行う必要があるため、補正値を算出する工数が膨大なものとなる。また、各測定点に基準音源を設置することが困難である場合には、その閉空間から推定対象音源を取り出す必要がある。こうした問題は、狭い閉空間に限らず、無響空間や半無響空間以外の閉空間に推定対象音源が設置されている場合に共通する。 However, in the method described in Patent Document 1, when the closed space is narrow, it is necessary to repeatedly install the reference sound source in the narrow closed space and measure with the sensor, so that the man-hours for calculating the correction value are enormous. It becomes a thing. If it is difficult to install a reference sound source at each measurement point, it is necessary to take out the estimation target sound source from the closed space. These problems are not limited to a narrow closed space, but are common when the sound source to be estimated is installed in a closed space other than an anechoic space or a semi-anechoic space.

本発明は、推定対象音源の音響パワーレベルを推定する際の補正値を容易に算出することのできる音響パワーレベルの推定方法を提供することを目的とする。 An object of the present invention is to provide a method for estimating a sound power level that can easily calculate a correction value when estimating the sound power level of a sound source to be estimated.

上記課題を解決する音響パワーレベルの推定方法は、推定対象音源の周囲に設置された複数の電気音響変換器で前記推定対象音源の音圧レベルを測定し、前記測定した音圧レベルに基づいて前記推定対象音源の音響パワーレベルを推定する。この推定方法は、前記複数の電気音響変換器の設置位置を決定する位置決め工程と、音響パワーレベルが既知である基準音源を音源として半無響空間に前記推定対象音源の測定環境を再現する再現工程と、前記再現された測定環境において前記複数の電気音響変換器が測定した前記基準音源の音圧レベルを取得する取得工程と、前記取得した音圧レベルに基づいて、前記推定対象音源の音響パワーレベルを推定する際の補正値を算出する算出工程と、を備える。 The sound power level estimation method for solving the above problems is to measure the sound pressure level of the estimation target sound source with a plurality of electroacoustic converters installed around the estimation target sound source, and based on the measured sound pressure level. The sound power level of the estimation target sound source is estimated. This estimation method reproduces the positioning step of determining the installation position of the plurality of electroacoustic converters and the measurement environment of the estimation target sound source in a semi-insensitive space using a reference sound source having a known sound power level as a sound source. The sound of the estimation target sound source based on the step, the acquisition step of acquiring the sound pressure level of the reference sound source measured by the plurality of electroacoustic converters in the reproduced measurement environment, and the acquired sound pressure level. It includes a calculation step of calculating a correction value when estimating a power level.

上記構成によれば、再現された推定対象音源の測定環境において電気音響変換器が測定した基準音源の音圧レベルに基づいて、推定対象音源の音響パワーレベルを推定する際の補正値が算出される。これにより、基準音源の設置回数を少なくできるため、補正値を算出する工数を少なくすることができる。また、閉空間から推定対象音源を取り出す必要もない。その結果、推定対象音源の音響パワーレベルを推定する際の補正値を容易に算出することができる。 According to the above configuration, a correction value for estimating the acoustic power level of the estimation target sound source is calculated based on the sound pressure level of the reference sound source measured by the electroacoustic converter in the measured environment of the reproduced estimation target sound source. NS. As a result, the number of times the reference sound source is installed can be reduced, and the man-hours for calculating the correction value can be reduced. Moreover, it is not necessary to extract the sound source to be estimated from the closed space. As a result, the correction value when estimating the sound power level of the sound source to be estimated can be easily calculated.

上記構成の音響パワーレベルの推定方法において、前記推定対象音源の音響パワーレベルは、前記推定対象音源の形状に基づいて分割した複数の推定面要素ごとの部分音響パワーレベルの平均値により推定される値である。前記複数の推定面要素の各々には、前記複数の電気音響変換器のうち、少なくとも1つが対応付けられている。前記部分音響パワーレベルは、前記対応付けられた電気音響変換器が測定した音圧レベルの平均値である平均音圧レベルを補正した値である。前記取得工程では、前記推定面要素ごとに、前記対応付けられた電気音響変換器が測定した前記基準音源の音圧レベルを取得する。前記算出工程では、前記推定面要素ごとに、前記取得工程で取得した音圧レベルの平均値と前記基準音源の音響パワーレベルとの差を、前記平均音圧レベルを補正する補正値として算出する。 In the sound power level estimation method having the above configuration, the sound power level of the estimation target sound source is estimated by the average value of the partial sound power levels of each of the plurality of estimation surface elements divided based on the shape of the estimation target sound source. The value. At least one of the plurality of electroacoustic transducers is associated with each of the plurality of estimation surface elements. The partial acoustic power level is a value obtained by correcting the average sound pressure level, which is the average value of the sound pressure levels measured by the associated electroacoustic converter. In the acquisition step, the sound pressure level of the reference sound source measured by the associated electroacoustic transducer is acquired for each estimated surface element. In the calculation step, the difference between the average value of the sound pressure level acquired in the acquisition step and the sound power level of the reference sound source is calculated as a correction value for correcting the average sound pressure level for each estimation surface element. ..

上記構成によれば、部分音響パワーレベルの推定に必要な面積情報を取得できるので、算出する部分音響パワーレベルの精度を向上させることができる。
上記構成の音響パワーレベルの推定方法において、前記測定した音圧レベルは、前記電気音響変換器の測定値を補正した値である。また、この推定方法は、前記推定対象音源の測定環境で前記複数の電気音響変換器の各々が測定する前記基準音源の音圧レベルを取得する基準音源取得工程をさらに備える。前記算出工程では、前記複数の電気音響変換器ごとに、前記基準音源取得工程で取得した音圧レベルと前記再現取得工程で取得した音圧レベルとの差を、前記測定値を補正する補正値として算出する。
According to the above configuration, since the area information necessary for estimating the partial acoustic power level can be acquired, the accuracy of the calculated partial acoustic power level can be improved.
In the method for estimating the acoustic power level having the above configuration, the measured sound pressure level is a value obtained by correcting the measured value of the electroacoustic converter. Further, this estimation method further includes a reference sound source acquisition step of acquiring the sound pressure level of the reference sound source measured by each of the plurality of electroacoustic converters in the measurement environment of the estimation target sound source. In the calculation step, for each of the plurality of electroacoustic converters, the difference between the sound pressure level acquired in the reference sound source acquisition step and the sound pressure level acquired in the reproduction acquisition step is a correction value for correcting the measured value. Calculate as.

上記構成によれば、閉空間における推定対象音源の反射波や回折、透過の影響を補正値によって排除することができる。その結果、部分音響パワーレベルの精度、ひいては推定対象音源の音響パワーレベルの精度を向上することができる。 According to the above configuration, the influence of the reflected wave, diffraction, and transmission of the estimation target sound source in the closed space can be eliminated by the correction value. As a result, the accuracy of the partial acoustic power level and, by extension, the accuracy of the acoustic power level of the estimation target sound source can be improved.

上記構成の音響パワーレベルの推定方法は、前記複数の電気音響変換器が測定する音圧レベルの周波数を複数の周波数帯域に区分けする区分工程と、前記複数の周波数帯域の各々について、前記音響パワーレベルを算出する際の補正方法を選択する選択工程と、をさらに備えていてもよい。 The method for estimating the acoustic power level having the above configuration includes a division step of dividing the frequency of the sound pressure level measured by the plurality of electroacoustic converters into a plurality of frequency bands, and the acoustic power for each of the plurality of frequency bands. It may further include a selection step of selecting a correction method when calculating the level.

発明者らは、周波数帯域によって適正な補正方法が異なることを見出した。上記構成によれば、推定した音響パワーレベルと実際の音響パワーレベルとの誤差を小さくすることができる。 The inventors have found that the appropriate correction method differs depending on the frequency band. According to the above configuration, the error between the estimated sound power level and the actual sound power level can be reduced.

推定対象音源の設置状態を模式的に示す図。The figure which shows the installation state of the estimation target sound source schematically. 音響パワーレベルを推定する推定システムの概略構成を示す図。The figure which shows the schematic structure of the estimation system which estimates the sound power level. 第1実施形態における音響パワーレベルの推定方法を示すフローチャート。The flowchart which shows the estimation method of the acoustic power level in 1st Embodiment. 第1実施形態において基準音源取得工程を模式的に示す図。The figure which shows typically the reference sound source acquisition process in 1st Embodiment. 第1実施形態において、(a)半無響空間に基準音源を設置した状態を模式的に示す図、(b)マイクロホンによる測定環境を再現した状態を模式的に示す図。In the first embodiment, (a) a diagram schematically showing a state in which a reference sound source is installed in a semi-anechoic space, and (b) a diagram schematically showing a state in which a measurement environment using a microphone is reproduced. 第1実施形態において音響パワーレベルの比較結果の一例を示すグラフ。The graph which shows an example of the comparison result of the acoustic power level in 1st Embodiment. 第2実施形態において音響パワーレベルの比較結果の一例を示すグラフ。The graph which shows an example of the comparison result of the acoustic power level in 2nd Embodiment. 第3実施形態において音響パワーレベルの比較結果の一例を示すグラフ。The graph which shows an example of the comparison result of the acoustic power level in 3rd Embodiment. 変形例において音響パワーレベルの推定方法の一部を示すフローチャート。The flowchart which shows a part of the sound power level estimation method in the modification. 変形例において音響パワーレベルの比較結果の一例を示すグラフ。The graph which shows an example of the comparison result of the acoustic power level in the modification example.

(第1実施形態)
図1〜図6を参照して、音響パワーレベルの推定方法の第1実施形態について説明する。
(First Embodiment)
A first embodiment of the method for estimating the sound power level will be described with reference to FIGS. 1 to 6.

図1に示すように、音響パワーレベルの推定方法においては、推定対象音源10の音響パワーレベルLw([dB])を推定する。推定対象音源10は、閉空間12に設置されており、音響パワーレベルLwが未知の音源である。 As shown in FIG. 1, in the sound power level estimation method, the sound power level Lw ([dB]) of the estimation target sound source 10 is estimated. The estimation target sound source 10 is installed in the closed space 12, and the sound power level Lw is unknown.

推定対象音源10は、その形状に基づいて複数の推定面要素Zi(iは自然数。)に分割される。具体的には、推定対象音源10は、その形状に基づいて多面体に簡略化され、その多面体を構成する各面が推定面要素Ziに設定される。例えば、推定対象音源10の簡略化した多面体が直方体であれば、その直方体を構成する各面が推定面要素Ziに設定される。推定対象音源10の音響パワーレベルLwは、それら推定面要素Ziごとの部分音響パワーレベルLwi([dB])を平均することにより求められる。「i」は、推定面要素Ziを特定する際に用いられる記号である。例えば、部分音響パワーレベルLw1は、推定面要素Z1の部分音響パワーレベルを示している。 The estimation target sound source 10 is divided into a plurality of estimation surface elements Zi (i is a natural number) based on its shape. Specifically, the estimation target sound source 10 is simplified into a polyhedron based on its shape, and each surface constituting the polyhedron is set in the estimation surface element Zi. For example, if the simplified polyhedron of the estimation target sound source 10 is a rectangular parallelepiped, each surface constituting the rectangular parallelepiped is set in the estimation surface element Zi. The sound power level Lw of the estimation target sound source 10 is obtained by averaging the partial sound power level Lwi ([dB]) for each of the estimated surface elements Zi. "I" is a symbol used when specifying the estimation surface element Zi. For example, the partial sound power level Lw1 indicates the partial sound power level of the estimated surface element Z1.

閉空間12には、複数のマイクロホンMが設置される。マイクロホンMは、電気音響変換器である。複数の推定面要素Ziの各々には、複数のマイクロホンMのうちの少なくとも1つが対応付けられている。推定面要素Ziに対応付けられた第j番目(jは自然数)のマイクロホンMをマイクロホンMijと表記する。また、マイクロホンMijの測定値を測定値Lmijと表記する。部分音響パワーレベルLwiは、対応付けられたマイクロホンMijの測定値Lmijに基づいて算出される。推定面要素の総数はN、推定面要素ごとのマイクロホンの総数はXiである。なお、図1では、推定面要素Z1に対応付けられた4つのマイクロホンM11,M12,M13,M14のみを示している。 A plurality of microphones M are installed in the closed space 12. The microphone M is an electroacoustic transducer. At least one of the plurality of microphones M is associated with each of the plurality of estimation surface elements Zi. The j-th (j is a natural number) microphone M associated with the estimation surface element Zi is referred to as microphone Mij. Further, the measured value of the microphone Mij is referred to as a measured value Lmij. The partial sound power level Lwi is calculated based on the measured value Lmij of the associated microphone Mij. The total number of estimated surface elements is N, and the total number of microphones for each estimated surface element is Xi. Note that FIG. 1 shows only the four microphones M11, M12, M13, and M14 associated with the estimation surface element Z1.

閉空間12の一例は、自動車のエンジンルームであり、推定対象音源10の一例はエンジンルームに設置されたエンジンである。図1では閉空間12および推定対象音源10を簡略化して示しているが、実際のエンジンルームにはエンジン付近にボンネットや補機類などの障害物が設置されている。こうした場合、マイクロホンMijは、エンジンおよび障害物との位置関係を考慮したうえでエンジン周辺の狭い空間に設置される。 An example of the closed space 12 is an engine room of an automobile, and an example of the estimation target sound source 10 is an engine installed in the engine room. Although the closed space 12 and the estimation target sound source 10 are shown in a simplified manner in FIG. 1, obstacles such as a bonnet and auxiliary equipment are installed in the vicinity of the engine in the actual engine room. In such a case, the microphone Mij is installed in a narrow space around the engine in consideration of the positional relationship between the engine and obstacles.

図2を参照して、推定対象音源10の音響パワーレベルLwを推定する推定システム20について説明する。
推定システム20は、マイクロホンMij、入力部21、および、推定装置25によって構成される。
The estimation system 20 for estimating the sound power level Lw of the estimation target sound source 10 will be described with reference to FIG. 2.
The estimation system 20 includes a microphone Mij, an input unit 21, and an estimation device 25.

各マイクロホンMijは、各々の測定値Lmijを推定装置25に入力する。
入力部21は、例えばキーボードやマウスなどであり、オペレーターによって操作されることで推定装置25に各種情報が入力可能に構成されている。また入力部21は、推定装置25が実行する処理を指示可能に構成されている。入力部21は、マイクロホンMijと推定面要素Ziとの対応付け関係を示す対応付け情報、後述する基準音源34の音響パワーレベルLwを示す基準音源情報、推定装置25における測定値Lmijの取り扱いを切り替えるモード情報、などが入力可能に構成されている。モード情報は、再現環境モード、基準音源モード、測定モードのいずれかである。
Each microphone Mij inputs each measured value Lmij to the estimation device 25.
The input unit 21 is, for example, a keyboard or a mouse, and is configured to be able to input various information to the estimation device 25 by being operated by an operator. Further, the input unit 21 is configured to be able to instruct the processing to be executed by the estimation device 25. The input unit 21 switches the handling of the correspondence information indicating the correspondence relationship between the microphone Mij and the estimation surface element Zi, the reference sound source information indicating the acoustic power level Lw of the reference sound source 34 described later, and the measured value Lmij in the estimation device 25. Mode information, etc. can be input. The mode information is one of the reproduction environment mode, the reference sound source mode, and the measurement mode.

推定装置25は、マイクロホンMijおよび入力部21からの各種情報、入力部21からの指示、および、メモリーに記憶したプログラムや各種のデータに基づいて各種の処理を実行する。推定装置25は、ASIC等の1つ以上の専用のハードウェア回路、コンピュータプログラム(ソフトウェア)に従って動作する1つ以上のプロセッサ、或いは、それらの組み合わせ、を含む回路として構成し得る。プロセッサは、CPU並びに、RAM及びROM等のメモリーを含み、メモリーは、処理をCPUに実行させるように構成されたプログラムコードまたは指令を格納している。メモリーすなわちコンピュータ可読媒体は、汎用または専用のコンピュータでアクセスできるあらゆる利用可能な媒体を含む。 The estimation device 25 executes various processes based on various information from the microphone Mij and the input unit 21, instructions from the input unit 21, a program stored in the memory, and various data. The estimation device 25 may be configured as a circuit including one or more dedicated hardware circuits such as an ASIC, one or more processors operating according to a computer program (software), or a combination thereof. The processor includes a CPU and a memory such as a RAM and a ROM, and the memory stores a program code or a command configured to cause the CPU to execute a process. Memory or computer-readable media includes any available medium accessible by a general purpose or dedicated computer.

推定装置25は、各種機能部として、取得部26、第1補正値算出部27、第2補正値算出部28、推定部29、および、記憶部30を有している。
取得部26は、各種情報を取得し、その取得した各種情報を示すデータを記憶部30の所定領域に格納する。
The estimation device 25 has an acquisition unit 26, a first correction value calculation unit 27, a second correction value calculation unit 28, an estimation unit 29, and a storage unit 30 as various functional units.
The acquisition unit 26 acquires various information and stores the data indicating the acquired various information in a predetermined area of the storage unit 30.

第1補正値算出部27は、推定面要素Ziの部分音響パワーレベルLwiを算出する際の補正値である補正値(第1補正値)Ai([dB])を推定面要素Ziごとに算出する。 The first correction value calculation unit 27 calculates a correction value (first correction value) Ai ([dB]), which is a correction value when calculating the partial sound power level Lwi of the estimation surface element Zi, for each estimation surface element Zi. do.

第2補正値算出部28は、マイクロホンMijの測定値Lmijを補正する補正値(第2補正値)Bij([dB])をマイクロホンMijごとに算出する。
推定部29は、各推定面要素Ziの部分音響パワーレベルLwiを算出し、その算出した部分音響パワーレベルLwiに基づいて推定対象音源10の音響パワーレベルLwを推定する。
The second correction value calculation unit 28 calculates a correction value (second correction value) Bij ([dB]) for correcting the measured value Lmij of the microphone Mij for each microphone Mij.
The estimation unit 29 calculates the partial sound power level Lwi of each estimated surface element Zi, and estimates the sound power level Lw of the estimation target sound source 10 based on the calculated partial sound power level Lwi.

図2〜図5を参照して、上述した推定システム20を用いた音響パワーレベルLwの推定方法の一例について説明する。
図3に示すように、まず、面要素設定工程が行われる(ステップS101)。面要素設定工程においては、推定対象音源10の形状に基づき、その推定対象音源10に対して複数の推定面要素Ziが設定される。
An example of a method for estimating the sound power level Lw using the above-mentioned estimation system 20 will be described with reference to FIGS. 2 to 5.
As shown in FIG. 3, first, a surface element setting step is performed (step S101). In the surface element setting step, a plurality of estimation surface elements Zi are set for the estimation target sound source 10 based on the shape of the estimation target sound source 10.

次に、位置決め工程が行われる(ステップS102)。位置決め工程においては、閉空間12における各マイクロホンMの設置位置が決定され、閉空間12に設置される。各マイクロホンMの設置位置は、推定対象音源10の各推定面要素Ziとの位置関係や推定対象音源10周りの障害物との位置関係などを考慮して決定される。各マイクロホンMの設置位置は、例えば、閉空間12内の基準点を原点とする3次元の座標系にて示すことができる。基準点は、例えば、推定対象音源10における特定位置(例えば中心位置)や特定のマイクロホンMの位置などである。 Next, a positioning step is performed (step S102). In the positioning step, the installation position of each microphone M in the closed space 12 is determined and installed in the closed space 12. The installation position of each microphone M is determined in consideration of the positional relationship of the estimation target sound source 10 with each estimation surface element Zi, the positional relationship with obstacles around the estimation target sound source 10, and the like. The installation position of each microphone M can be indicated by, for example, a three-dimensional coordinate system having a reference point in the closed space 12 as an origin. The reference point is, for example, a specific position (for example, a center position) in the estimation target sound source 10 or a position of a specific microphone M.

また、位置決め工程では、マイクロホンMと推定面要素Ziとの対応付けも行われる。対応付けにおいては、各マイクロホンMが複数の推定面要素Ziのいずれかに対応付けられるとともに、各推定面要素Ziに少なくとも1つのマイクロホンMが対応付けられる。こうして推定面要素Ziに対応付けられた第j番目のマイクロホンMは、マイクロホンMijとして表記される。マイクロホンMijと推定面要素Ziとの対応付けは、対応付け情報として入力部21を通じて推定装置25に入力される。その対応付け情報は、取得部26によって取得され、対応付けデータ31として記憶部30の所定領域に格納される。 Further, in the positioning step, the microphone M and the estimation surface element Zi are also associated with each other. In the association, each microphone M is associated with any of the plurality of estimated surface elements Zi, and at least one microphone M is associated with each estimated surface element Zi. The j-th microphone M associated with the estimation surface element Zi in this way is expressed as microphone Mij. The association between the microphone Mij and the estimation surface element Zi is input to the estimation device 25 as the association information through the input unit 21. The correspondence information is acquired by the acquisition unit 26 and stored as the association data 31 in a predetermined area of the storage unit 30.

次に、測定工程(ステップS103)が行われる。測定工程においては、推定対象音源10の音圧レベルが各マイクロホンMijで測定される。具体的には、入力部21を通じて測定モードを示すモード情報が推定装置25に入力される。その後、推定対象音源10を作動させ、入力部21を通じて推定装置25に測定開始を指示する。測定モードにおいて、取得部26は、各マイクロホンMijの測定値Lmijを取得し、その取得した各マイクロホンMijの測定値Lmijに基づく推定用測定データ43を記憶部30の所定領域に格納する。以下では、推定用測定データ43を構成する測定値Lmijを推定用測定値Lrijという。 Next, the measurement step (step S103) is performed. In the measurement step, the sound pressure level of the estimation target sound source 10 is measured by each microphone Mij. Specifically, mode information indicating the measurement mode is input to the estimation device 25 through the input unit 21. After that, the estimation target sound source 10 is operated, and the estimation device 25 is instructed to start measurement through the input unit 21. In the measurement mode, the acquisition unit 26 acquires the measured value Lmij of each microphone Mij, and stores the estimation measurement data 43 based on the acquired measured value Lmij of each microphone Mij in a predetermined area of the storage unit 30. Hereinafter, the measured value Lmij constituting the estimation measurement data 43 is referred to as an estimation measurement value Lrig.

なお、推定用測定値Lrijは、閉空間12に全てのマイクロホンMijを設置してから測定工程を行うことで、全ての推定面要素Ziについて同時期に測定される構成であってもよい。また、推定用測定値Lrijは、面要素設定工程、位置決め工程、測定工程を繰り返すことで、推定面要素Ziごとに測定される構成であってもよい。 The estimated measured value Lrig may be configured to be measured at the same time for all the estimated surface elements Zi by performing the measuring step after installing all the microphones Mij in the closed space 12. Further, the estimated measured value Lrig may be configured to be measured for each estimated surface element Zi by repeating the surface element setting step, the positioning step, and the measuring step.

全ての推定面要素Ziについて推定用測定値Lrijが測定されると(ステップS104:YES)、基準音源取得工程(ステップS105)が行われる。
基準音源取得工程の開始時、各推定面要素Ziに対応付けられたマイクロホンMijは、閉空間12に設置された状態にある。基準音源取得工程においては、推定面要素Ziごとに、対応付けられたマイクロホンMijが測定する基準音源34の音圧レベルを取得する。基準音源34は、例えば、既知音響パワーレベルLwKを有する点音源である。基準音源34の既知音響パワーレベルLwKを示す基準音源情報は、入力部21を通じて推定装置25に入力される。推定装置25に入力された基準音源情報は、取得部26によって取得され、基準音源データ35として記憶部30の所定領域に格納される。
When the estimation measurement value Lrig is measured for all the estimation surface elements Zi (step S104: YES), the reference sound source acquisition step (step S105) is performed.
At the start of the reference sound source acquisition process, the microphone Mij associated with each estimation surface element Zi is in a state of being installed in the closed space 12. In the reference sound source acquisition step, the sound pressure level of the reference sound source 34 measured by the associated microphone Mij is acquired for each estimation surface element Zi. The reference sound source 34 is, for example, a point sound source having a known sound power level LwK. The reference sound source information indicating the known sound power level LwK of the reference sound source 34 is input to the estimation device 25 through the input unit 21. The reference sound source information input to the estimation device 25 is acquired by the acquisition unit 26 and stored in a predetermined area of the storage unit 30 as the reference sound source data 35.

図4を参照して、基準音源取得工程について推定面要素Z1を例に具体的に説明する。
図4に示すように、推定面要素Z1に対応付けられたマイクロホンM1jで基準音源34の音圧レベルを取得するときには、推定面要素Z1の重心位置などに基準音源34を設置する。
With reference to FIG. 4, the reference sound source acquisition process will be specifically described using the estimation surface element Z1 as an example.
As shown in FIG. 4, when the sound pressure level of the reference sound source 34 is acquired by the microphone M1j associated with the estimation surface element Z1, the reference sound source 34 is installed at the position of the center of gravity of the estimation surface element Z1 or the like.

次に、入力部21を通じて、基準音源モードを示すモード情報と、推定面要素Z1を示す情報とを推定装置25に入力する。その後、基準音源34を作動させ、入力部21を通じて推定装置25に測定開始を指示する。このときの基準音源モードにおいて、取得部26は、推定面要素Z1に対応付けられたマイクロホンM1jの測定値Lm1jを取得し、その取得したマイクロホンM1jの測定値Lm1jを示すデータを記憶部30の所定領域に格納する。 Next, the mode information indicating the reference sound source mode and the information indicating the estimation surface element Z1 are input to the estimation device 25 through the input unit 21. After that, the reference sound source 34 is operated, and the estimation device 25 is instructed to start measurement through the input unit 21. In the reference sound source mode at this time, the acquisition unit 26 acquires the measured value Lm1j of the microphone M1j associated with the estimation surface element Z1, and stores the data indicating the acquired measured value Lm1j of the microphone M1j in the predetermined storage unit 30. Store in the area.

基準音源取得工程では、上述した基準音源34の音圧レベルの測定を推定面要素Z1,Z2,…,Ziの各々に対して行う。そして、全ての推定面要素Ziについての測定が終了したとき(ステップS106:YES)、記憶部30には、基準音源取得工程で取得した測定値Lmijで構成される基準音源測定データ41が格納されている。基準音源測定データ41は、補正値Bijを算出する際に用いられる測定データである。以下では、基準音源測定データ41を構成する測定値Lmijを基準音源測定値Ltijという。 In the reference sound source acquisition step, the sound pressure level of the reference sound source 34 described above is measured for each of the estimation surface elements Z1, Z2, ..., Zi. When the measurement for all the estimated surface elements Zi is completed (step S106: YES), the storage unit 30 stores the reference sound source measurement data 41 composed of the measured values Lmij acquired in the reference sound source acquisition step. ing. The reference sound source measurement data 41 is measurement data used when calculating the correction value Bij. Hereinafter, the measured value Lmij constituting the reference sound source measurement data 41 is referred to as a reference sound source measurement value Ltij.

次に、推定面要素Ziごとに、再現工程(ステップS107)と再現取得工程(ステップS108)とが行われる。
図5を参照して、再現工程(ステップS107)と再現取得工程(ステップS108)とについて、推定面要素Z1を例に具体的に説明する。
Next, a reproduction step (step S107) and a reproduction acquisition step (step S108) are performed for each estimation surface element Zi.
With reference to FIG. 5, the reproduction step (step S107) and the reproduction acquisition step (step S108) will be specifically described by taking the estimation surface element Z1 as an example.

再現工程においては、床面を有する半無響空間に、基準音源34を測定対象として、推定面要素Z1に対応付けられたマイクロホンM1jによる推定対象音源10の測定環境が再現される。 In the reproduction step, the measurement environment of the estimation target sound source 10 by the microphone M1j associated with the estimation surface element Z1 is reproduced in the semi-anechoic space having the floor surface with the reference sound source 34 as the measurement target.

具体的には、まず、図5(a)に示すように、半無響空間33の床面32に基準音源34を設置し、その基準音源34を基準とした床面32の一部を推定面要素Z1と仮定する。基準音源34と推定面要素Z1との位置関係は、基準音源取得工程(ステップS105)における基準音源34と推定面要素Z1との位置関係を表現するものである。 Specifically, first, as shown in FIG. 5A, a reference sound source 34 is installed on the floor surface 32 of the semi-anechoic space 33, and a part of the floor surface 32 based on the reference sound source 34 is estimated. It is assumed that the surface element Z1. The positional relationship between the reference sound source 34 and the estimated surface element Z1 expresses the positional relationship between the reference sound source 34 and the estimated surface element Z1 in the reference sound source acquisition step (step S105).

次に、図5(b)に示すように、基準音源取得工程における基準音源34とマイクロホンM1jとの位置関係に基づいて、基準音源34を基準にマイクロホンM1jを配置する。これにより、推定面要素Z1について、マイクロホンM1jによる推定対象音源10の測定環境が半無響空間33に再現される。 Next, as shown in FIG. 5B, the microphone M1j is arranged with reference to the reference sound source 34 based on the positional relationship between the reference sound source 34 and the microphone M1j in the reference sound source acquisition process. As a result, the measurement environment of the estimation target sound source 10 by the microphone M1j is reproduced in the semi-anechoic space 33 for the estimation surface element Z1.

再現取得工程(ステップS108)においては、再現工程(ステップS107)で再現された測定環境でマイクロホンM1jが測定する基準音源34の音圧レベルを取得する。
具体的には、入力部21を通じて、再現環境モードを示すモード情報と、推定面要素Z1を示す情報とが推定装置25に入力される。その後、基準音源34を作動させ、入力部21を通じて推定装置25に測定開始を指示する。このときの再現環境モードにおいて、取得部26は、推定面要素Z1に対応付けられたマイクロホンM1jの測定値Lm1jを取得し、その取得したマイクロホンM1jの測定値Lm1jを示すデータを記憶部30の所定領域に格納する。
In the reproduction acquisition step (step S108), the sound pressure level of the reference sound source 34 measured by the microphone M1j in the measurement environment reproduced in the reproduction step (step S107) is acquired.
Specifically, the mode information indicating the reproduction environment mode and the information indicating the estimation surface element Z1 are input to the estimation device 25 through the input unit 21. After that, the reference sound source 34 is operated, and the estimation device 25 is instructed to start measurement through the input unit 21. In the reproduction environment mode at this time, the acquisition unit 26 acquires the measured value Lm1j of the microphone M1j associated with the estimation surface element Z1, and stores the data indicating the acquired measured value Lm1j of the microphone M1j in the predetermined storage unit 30. Store in the area.

再現工程(ステップS107)と再現取得工程(ステップS108)は、推定面要素Z1,Z2,…,Ziの各々に対して行われる。全ての推定面要素Ziについての測定が終了したとき(ステップS109:YES)、記憶部30には、再現取得工程で取得した測定値Lmijで構成される再現環境測定データ36が格納されている。以下では、再現環境測定データ36を構成する測定値Lmijを再現環境測定値Lcijという。再現環境測定データ36は、補正値Aiおよび補正値Bijを算出する際に用いられる測定データである。 The reproduction step (step S107) and the reproduction acquisition step (step S108) are performed for each of the estimation surface elements Z1, Z2, ..., Zi. When the measurement for all the estimated surface elements Zi is completed (step S109: YES), the storage unit 30 stores the reproduction environment measurement data 36 composed of the measured values Lmij acquired in the reproduction acquisition step. In the following, the measured value Lmij constituting the reproduction environment measurement data 36 is referred to as a reproduction environment measurement value Lcij. The reproduction environment measurement data 36 is measurement data used when calculating the correction value Ai and the correction value Bij.

次に、第1算出工程(ステップS110)が行われる。第1算出工程においては、部分音響パワーレベルLwiを算出する際の補正値である補正値Aiが推定面要素Ziごとに算出される。補正値Aiは、各推定面要素Ziの部分音響パワーレベルLwiの推定に必要な面積情報を示す補正値である。 Next, the first calculation step (step S110) is performed. In the first calculation step, the correction value Ai, which is the correction value when calculating the partial sound power level Lwi, is calculated for each estimated surface element Zi. The correction value Ai is a correction value indicating the area information necessary for estimating the partial sound power level Lwi of each estimated surface element Zi.

具体的には、入力部21を通じて補正値Aiの算出指示が推定装置25に入力される。補正値Aiの算出指示が入力されると、第1補正値算出部27は、算出対象となる推定面要素Ziに対応付けられたマイクロホンMijの再現環境測定値Lcijを再現環境測定データ36から抽出する。そして、第1補正値算出部27は、式(1)に示すように、その抽出した再現環境測定値Lcijの平均値を基準音源34の既知音響パワーレベルLwKから減算することにより補正値Aiを算出する。なお、Siは算出対象となる推定面要素Ziに対応付けられたマイクロホンMijの配置部分の面積である。 Specifically, the calculation instruction of the correction value Ai is input to the estimation device 25 through the input unit 21. When the calculation instruction of the correction value Ai is input, the first correction value calculation unit 27 extracts the reproduction environment measurement value Lcij of the microphone Mij associated with the estimation surface element Zi to be calculated from the reproduction environment measurement data 36. do. Then, as shown in the equation (1), the first correction value calculation unit 27 subtracts the average value of the extracted reproduction environment measurement value Lcij from the known sound power level LwK of the reference sound source 34 to obtain the correction value Ai. calculate. In addition, Si is the area of the arrangement portion of the microphone Mij associated with the estimation surface element Zi to be calculated.

Figure 2021183910
Figure 2021183910

第1補正値算出部27は、補正値Aiを推定面要素Ziごとに算出し、その算出した補正値Aiと推定面要素Ziとを対応付けた第1補正データ37を記憶部30の所定領域に格納する。 The first correction value calculation unit 27 calculates the correction value Ai for each estimated surface element Zi, and stores the first correction data 37 in which the calculated correction value Ai and the estimated surface element Zi are associated with each other in a predetermined area of the storage unit 30. Store in.

次に、第2算出工程(ステップS111)が行われる。第2算出工程においては、マイクロホンMijの測定値Lmijを補正する補正値BijがマイクロホンMijごとに算出される。補正値Bijは、測定値Lmijから閉空間における対象音源の反射波や回折、透過の影響を排除する補正値である。 Next, the second calculation step (step S111) is performed. In the second calculation step, a correction value Bij that corrects the measured value Lmij of the microphone Mij is calculated for each microphone Mij. The correction value Bij is a correction value that excludes the influence of the reflected wave, diffraction, and transmission of the target sound source in the closed space from the measured value Lmij.

具体的には、入力部21を通じて補正値Bijの算出指示が推定装置25に入力される。補正値Bijの算出指示が入力されると、第2補正値算出部28は、算出対象となるマイクロホンMijについて、再現環境測定データ36から再現環境測定値Lcij、基準音源測定データ41から基準音源測定値Ltijを抽出する。そして、第2補正値算出部28は、式(2)に示すように、その抽出した基準音源測定値Ltijから再現環境測定値Lcijを減算することにより補正値Bijを算出する。 Specifically, the calculation instruction of the correction value Bij is input to the estimation device 25 through the input unit 21. When the calculation instruction of the correction value Bij is input, the second correction value calculation unit 28 measures the reference sound source from the reproduction environment measurement data 36 to the reproduction environment measurement value Lcij and the reference sound source measurement data 41 for the microphone Mij to be calculated. Extract the value Ltij. Then, as shown in the equation (2), the second correction value calculation unit 28 calculates the correction value Bij by subtracting the reproduction environment measurement value Lcij from the extracted reference sound source measurement value Ltij.

Figure 2021183910
Figure 2021183910

第2補正値算出部28は、補正値BijをマイクロホンMijごとに算出し、その算出した補正値BijとマイクロホンMijとを対応付けた第2補正データ42を記憶部30の所定領域に格納する。 The second correction value calculation unit 28 calculates the correction value Bij for each microphone Mij, and stores the second correction data 42 in which the calculated correction value Bij and the microphone Mij are associated with each other in a predetermined area of the storage unit 30.

次に、推定工程(ステップS112)が行われる。推定工程においては、推定対象音源10の音響パワーレベルLwが推定される。
具体的には、入力部21を通じて推定対象音源10の音響パワーレベルLwの推定指示が推定装置25に入力される。推定指示が入力されると、推定部29は、推定用測定データ43と第2補正データ42とに基づいて、各マイクロホンMijについて、式(3)に示すように、推定用測定値Lrijから補正値Bijを減算した値を、マイクロホンMijが測定した音圧レベルである推定用音圧レベルLijとして算出する。そして、推定部29は、式(4)に示すように、推定面要素Ziごとに、対応付けられたマイクロホンMijの推定用音圧レベルLijの平均値を算出する。推定部29は、式(5)に示すように、その算出した平均値に対して補正値Aiを加算することにより推定面要素Ziごとに部分音響パワーレベルLwiを算出する。推定部29は、推定面要素Ziごとに算出された部分音響パワーレベルLwiの平均値を推定対象音源10の音響パワーレベルLwとして推定する。
Next, the estimation step (step S112) is performed. In the estimation process, the sound power level Lw of the estimation target sound source 10 is estimated.
Specifically, the estimation instruction of the sound power level Lw of the estimation target sound source 10 is input to the estimation device 25 through the input unit 21. When the estimation instruction is input, the estimation unit 29 corrects each microphone Mij from the estimation measurement value Lrig based on the estimation measurement data 43 and the second correction data 42, as shown in the equation (3). The value obtained by subtracting the value Bij is calculated as the estimated sound pressure level Lij, which is the sound pressure level measured by the microphone Mij. Then, as shown in the equation (4), the estimation unit 29 calculates the average value of the estimation sound pressure level Lij of the associated microphone Mij for each estimation surface element Zi. As shown in the equation (5), the estimation unit 29 calculates the partial sound power level Lwi for each estimated surface element Zi by adding the correction value Ai to the calculated average value. The estimation unit 29 estimates the average value of the partial sound power level Lwi calculated for each estimation surface element Zi as the sound power level Lw of the estimation target sound source 10.

Figure 2021183910
Figure 2021183910

Figure 2021183910
Figure 2021183910

Figure 2021183910
Figure 2021183910

図6は、上述した推定方法によって推定された推定対象音源10の音響パワーレベルLwについて行った比較実験結果の一例である。この推定対象音源10は、自動車のエンジンルームに設置されたエンジンである。 FIG. 6 is an example of the results of a comparative experiment performed on the sound power level Lw of the estimation target sound source 10 estimated by the above estimation method. The estimation target sound source 10 is an engine installed in the engine room of an automobile.

図6に示すように、推定された音響パワーレベルLwは、補正値Aiおよび補正値Bijによる補正を行わない場合に比べて音響パワーレベルLwの真値に対する誤差が小さくなることが確認された。すなわち、上述した推定方法によれば、エンジンルームに設置されたエンジンのような狭い閉空間12に設置された推定対象音源10であっても、推定対象音源10を閉空間12から取り出すことなく、その音響パワーレベルLwを精度よく推定できることが確認された。 As shown in FIG. 6, it was confirmed that the estimated sound power level Lw has a smaller error with respect to the true value of the sound power level Lw than in the case where the correction by the correction value Ai and the correction value Bij is not performed. That is, according to the estimation method described above, even if the estimation target sound source 10 is installed in a narrow closed space 12 such as an engine installed in an engine room, the estimation target sound source 10 is not taken out from the closed space 12. It was confirmed that the acoustic power level Lw can be estimated accurately.

第1実施形態の作用及び効果について説明する。
(1−1)上述した音響パワーレベルLwの推定方法においては、マイクロホンMijによる推定対象音源10の測定環境が、基準音源34を測定対象とした半無響空間33に再現される。そして、基準音源34の既知音響パワーレベルLwKと、その再現された測定環境における各マイクロホンMijの再現環境測定値Lcijとに基づいて、推定面要素Ziごとの補正値Aiが算出される。
The operation and effect of the first embodiment will be described.
(1-1) In the above-mentioned sound power level Lw estimation method, the measurement environment of the estimation target sound source 10 by the microphone Mij is reproduced in the semi-anechoic space 33 with the reference sound source 34 as the measurement target. Then, the correction value Ai for each estimation surface element Zi is calculated based on the known sound power level LwK of the reference sound source 34 and the reproduction environment measurement value Lcij of each microphone Mij in the reproduced measurement environment.

また、マイクロホンMijによる推定対象音源10の測定環境が、基準音源34を測定対象として再現される。そして、上述した再現環境測定値Lcijと、その再現された測定環境における各マイクロホンMijの基準音源測定値Ltijとに基づいて、補正値Bijが算出される。 Further, the measurement environment of the estimation target sound source 10 by the microphone Mij is reproduced with the reference sound source 34 as the measurement target. Then, the correction value Bij is calculated based on the above-mentioned reproduction environment measurement value Lcij and the reference sound source measurement value Ltij of each microphone Mij in the reproduced measurement environment.

そして、これら補正値Ai、補正値Bij、および、推定用測定値Lrijに基づいて、各推定対象音源10の音響パワーレベルLwが推定される。
この構成によれば、基準音源34の設置回数を少なくできるため、推定対象音源10の音響パワーレベルLwを推定する際の補正値をより少ない工数で算出できる。その結果、該補正値を容易に算出できる。
Then, the sound power level Lw of each estimation target sound source 10 is estimated based on the correction value Ai, the correction value Bij, and the estimation measurement value Lrig.
According to this configuration, since the number of times the reference sound source 34 is installed can be reduced, the correction value when estimating the sound power level Lw of the estimation target sound source 10 can be calculated with less man-hours. As a result, the correction value can be easily calculated.

(1−2)上記構成によれば、推定対象音源10が狭い閉空間12に配設されている場合であっても、その閉空間12から推定対象音源10を取り出して推定対象音源10の音圧レベルを測定する必要がない。そのため、推定対象音源10の音響パワーレベルLwを簡易な方法のもとで推定できる。 (1-2) According to the above configuration, even when the estimation target sound source 10 is arranged in the narrow closed space 12, the estimation target sound source 10 is taken out from the closed space 12 and the sound of the estimation target sound source 10 is taken out. There is no need to measure the pressure level. Therefore, the sound power level Lw of the estimation target sound source 10 can be estimated by a simple method.

(1−3)上記構成によれば、部分音響パワーレベルLwiの推定に必要な各推定面要素Ziの面積情報として補正値Aiを取得できるため、各推定面要素Ziの部分音響パワーレベルLwiの精度、ひいては推定対象音源10の音響パワーレベルLwの精度が向上できる。 (1-3) According to the above configuration, since the correction value Ai can be acquired as the area information of each estimated surface element Zi necessary for estimating the partial sound power level Lwi, the partial sound power level Lwi of each estimated surface element Zi can be acquired. The accuracy, and thus the accuracy of the sound power level Lw of the estimation target sound source 10, can be improved.

(1−4)上記構成によれば、推定用音圧レベルLijは、閉空間12における推定対象音源10の反射波や回折、透過の影響を補正値Bijによって推定用測定値Lrijから排除した値とすることができる。その結果、部分音響パワーレベルLwiの精度、ひいては推定対象音源10の音響パワーレベルLwの精度が向上できる。 (1-4) According to the above configuration, the estimation sound pressure level Lij is a value obtained by excluding the influence of the reflected wave, diffraction, and transmission of the estimation target sound source 10 in the closed space 12 from the estimation measurement value Lrig by the correction value Bij. Can be. As a result, the accuracy of the partial acoustic power level Lwi, and by extension, the accuracy of the acoustic power level Lw of the estimation target sound source 10 can be improved.

(第2実施形態)
図7を参照して、音響パワーレベルの推定方法の第2実施形態について説明する。なお、第2実施形態の音響パワーレベルの推定方法は、第1実施形態の音響パワーレベルの推定方法と補正値Aiの算出方法が異なる。そのため、第2実施形態においては、補正値Ai´の算出方法について詳細に説明する。
(Second Embodiment)
A second embodiment of the method for estimating the sound power level will be described with reference to FIG. 7. The method for estimating the sound power level of the second embodiment is different from the method for estimating the sound power level of the first embodiment and the method for calculating the correction value Ai. Therefore, in the second embodiment, the calculation method of the correction value Ai'will be described in detail.

第2実施形態の第1算出工程(ステップS110)においては、面要素設定工程(ステップS101)で設定された各推定面要素Ziの実際の面積Siが算出されるとともに、その面積Siに基づいて補正値Ai´が算出される。 In the first calculation step (step S110) of the second embodiment, the actual area Si of each estimated surface element Zi set in the surface element setting step (step S101) is calculated, and the actual area Si is calculated based on the area Si. The correction value Ai'is calculated.

具体的には、入力部21を通じて、推定対象音源10の形状に基づいて分割された推定面要素Ziの形状を示す面要素情報が推定装置25に入力される。面要素情報は、各推定面要素Ziについて、例えば、各頂点の座標を示す情報であってもよいし、各辺の長さや各内角の大きさを示す情報であってもよい。 Specifically, the surface element information indicating the shape of the estimated surface element Zi divided based on the shape of the estimation target sound source 10 is input to the estimation device 25 through the input unit 21. The surface element information may be, for example, information indicating the coordinates of each vertex or information indicating the length of each side or the size of each internal angle for each estimated surface element Zi.

面要素情報が入力されると、取得部26は、その入力された面要素情報を取得する。第1補正値算出部27は、取得部26が取得した面要素情報に基づいて推定面要素Ziの面積Siを算出し、その算出した面積Siを式(6)に代入して得られる値を推定面要素Ziの補正値Ai´として算出する。 When the surface element information is input, the acquisition unit 26 acquires the input surface element information. The first correction value calculation unit 27 calculates the area Si of the estimated surface element Zi based on the surface element information acquired by the acquisition unit 26, and substitutes the calculated area Si into the equation (6) to obtain a value. It is calculated as the correction value Ai'of the estimated surface element Zi.

Figure 2021183910
Figure 2021183910

図7は、こうして算出された補正値Ai´と前述の補正値Bijとを用いて推定した推定対象音源10の音響パワーレベルLwについて行った比較実験結果の一例である。推定対象音源10は、第1実施形態と同様、エンジンルームに設置されたエンジンである。 FIG. 7 is an example of the results of a comparative experiment performed on the sound power level Lw of the estimation target sound source 10 estimated using the correction value Ai'calculated in this way and the above-mentioned correction value Bij. The estimation target sound source 10 is an engine installed in the engine room as in the first embodiment.

図7に示すように、推定された音響パワーレベルLwは、補正値Ai´および補正値Bijによる補正を行わない場合に比べて音響パワーレベルLwの真値に対する誤差が小さくなることが確認された。特に、この推定方法においては、特定周波数F1(本比較実験では1000Hz)よりも低い周波数において誤差が小さいことが確認された。 As shown in FIG. 7, it was confirmed that the estimated sound power level Lw has a smaller error with respect to the true value of the sound power level Lw than when the correction by the correction value Ai'and the correction value Bij is not performed. .. In particular, in this estimation method, it was confirmed that the error is small at a frequency lower than the specific frequency F1 (1000 Hz in this comparative experiment).

第2実施形態の効果について説明する。
(2−1)上記構成であっても、第1実施形態の(1−1)〜(1−4)に記載した効果に準ずる効果を得ることができる。
The effect of the second embodiment will be described.
(2-1) Even with the above configuration, it is possible to obtain an effect similar to the effect described in (1-1) to (1-4) of the first embodiment.

(第3実施形態)
図8を参照して、音響パワーレベルの推定方法の第3実施形態について説明する。第3実施形態の音響パワーレベルの推定方法は、第1実施形態の音響パワーレベルの推定方法と基準音源測定値Ltijを用いて算出される補正値が異なる。そのため、第3実施形態については、その補正値の算出方法について詳細に説明する。
(Third Embodiment)
A third embodiment of the method for estimating the sound power level will be described with reference to FIG. The sound power level estimation method of the third embodiment differs from the sound power level estimation method of the first embodiment in the correction value calculated by using the reference sound source measurement value Ltij. Therefore, with respect to the third embodiment, a method of calculating the correction value will be described in detail.

第3実施形態の第2算出工程(ステップS111)においては、基準音源取得工程(ステップS106)で取得された基準音源測定値Ltijに基づいて、第2補正値算出部28によって推定面要素Ziごとに補正値(第3補正値)Ciが算出される。 In the second calculation step (step S111) of the third embodiment, for each estimated surface element Zi by the second correction value calculation unit 28 based on the reference sound source measurement value Ltij acquired in the reference sound source acquisition step (step S106). The correction value (third correction value) Ci is calculated.

具体的には、入力部21を通じて、推定対象音源10の形状に基づいて分割された推定面要素Ziの形状を示す面要素情報が推定装置25に入力される。面要素情報は、各推定面要素Ziについて、例えば、各頂点の座標を示す情報であってもよいし、各辺の長さや各内角の大きさを示す情報であってもよい。 Specifically, the surface element information indicating the shape of the estimated surface element Zi divided based on the shape of the estimation target sound source 10 is input to the estimation device 25 through the input unit 21. The surface element information may be, for example, information indicating the coordinates of each vertex or information indicating the length of each side or the size of each internal angle for each estimated surface element Zi.

面要素情報が入力されると、取得部26は、その入力された面要素情報を取得する。第2補正値算出部28は、推定面要素Ziごとに、式(7)に示すように補正値Ciを算出する。具体的には、第2補正値算出部28は、まず、式(7)における右辺第1項および第2項のように、基準音源測定値Ltijに基づく部分音響パワーレベルを算出する。Siは、取得部26が取得した面要素情報に基づいて算出される推定面要素Ziの面積である。そして第2補正値算出部28は、その部分音響パワーレベルから再現環境測定値Lcijに基づく部分音響パワーレベルを減算することにより、すなわち、その部分音響パワーレベルから既知音響パワーレベルLwKを減算することにより補正値Ciを算出する。 When the surface element information is input, the acquisition unit 26 acquires the input surface element information. The second correction value calculation unit 28 calculates the correction value Ci as shown in the equation (7) for each estimation surface element Zi. Specifically, the second correction value calculation unit 28 first calculates the partial sound power level based on the reference sound source measurement value Ltij, as in the first and second terms on the right side in the equation (7). Si is the area of the estimated surface element Zi calculated based on the surface element information acquired by the acquisition unit 26. Then, the second correction value calculation unit 28 subtracts the known acoustic power level LwK from the partial acoustic power level by subtracting the partial acoustic power level based on the reproduced environment measurement value Lcij from the partial acoustic power level. The correction value Ci is calculated by.

Figure 2021183910
Figure 2021183910

推定工程(ステップS112)において、推定部29は、式(8)に示すように、推定面要素Ziごとに、推定用測定値Lrijの平均値(右辺第1項)に推定面要素の面積に関する値(右辺第2項)を加算した値から補正値Ciを減算することで部分音響パワーレベルLwiを算出する。そして、推定部29は、式(9)に従って、推定面要素Ziごとに算出された部分音響パワーレベルLwiの平均値を推定対象音源10の音響パワーレベルLwとして推定する。 In the estimation step (step S112), as shown in the equation (8), the estimation unit 29 relates to the area of the estimation surface element in the average value (first term on the right side) of the estimation measurement value Lrig for each estimation surface element Zi. The partial sound power level Lwi is calculated by subtracting the correction value Ci from the value obtained by adding the value (second term on the right side). Then, the estimation unit 29 estimates the average value of the partial sound power level Lwi calculated for each estimation surface element Zi as the sound power level Lw of the estimation target sound source 10 according to the equation (9).

Figure 2021183910
Figure 2021183910

Figure 2021183910
Figure 2021183910

図8は、こうして算出された推定対象音源10の音響パワーレベルLwについて行った比較実験結果の一例である。図8に示すように、推定された音響パワーレベルLwは、補正値Ciによる補正を行わない場合に比べて音響パワーレベルLwの真値に対する誤差が小さくなることが確認された。特に、この推定方法においては、特定周波数F1(本比較実験では1000Hz)付近よりも高い周波数において誤差が小さいことが確認された。 FIG. 8 is an example of the results of a comparative experiment performed on the sound power level Lw of the estimation target sound source 10 calculated in this way. As shown in FIG. 8, it was confirmed that the estimated sound power level Lw has a smaller error with respect to the true value of the sound power level Lw than in the case where the correction by the correction value Ci is not performed. In particular, in this estimation method, it was confirmed that the error is small at a frequency higher than the specific frequency F1 (1000 Hz in this comparative experiment).

第3実施形態の効果について説明する。
(3−1)上記構成であっても、第1実施形態の(1−1)〜(1−3)に記載した効果に準ずる効果を得ることができる。
The effect of the third embodiment will be described.
(3-1) Even with the above configuration, it is possible to obtain an effect similar to the effect described in (1-1) to (1-3) of the first embodiment.

本実施形態は、以下のように変更して実施することができる。本実施形態及び以下の変更例は、技術的に矛盾しない範囲で互いに組み合わせて実施することができる。
・各種測定データ36,41,43を得るための測定は、マイクロホンMijの設置位置の決定後であれば、どのような順番で行われてもよい。また、各種補正値についても、その補正値の算出で必要な測定データが得られたあとに算出されればよく、どのような順番で算出されてもよい。
This embodiment can be modified and implemented as follows. The present embodiment and the following modified examples can be implemented in combination with each other within a technically consistent range.
-Measurements for obtaining various measurement data 36, 41, 43 may be performed in any order as long as the installation position of the microphone Mij is determined. Further, the various correction values may be calculated after the measurement data required for the calculation of the correction values is obtained, and may be calculated in any order.

・推定装置25は、各種補正値の少なくとも1つを算出し、その算出した補正値を用いて音響パワーレベルLwを推定する構成であってもよい。
・上述したように、本発明者らは、第2実施形態の推定方法においては特定周波数F1よりも低い周波数において真値に対する誤差が小さく、かつ、第3実施形態の推定方法においては特定周波数F1よりも高い周波数において真値に対する誤差が小さいことを見出した。すなわち、周波数帯域によって適正な補正方法が異なることを見出した。また、本発明者らは、特定周波数F1について、推定対象音源10の形状によって異なる値であり、基準音源34を用いて予め行った実験やシミュレーションの結果に基づいて、推定対象音源10の形状ごとに設定することが可能であることを見出した。このことから、推定対象音源10の音響パワーレベルLwは、推定対象音源10の形状に基づき、周波数帯域ごとに異なる補正方法で推定されてもよい。
The estimation device 25 may be configured to calculate at least one of various correction values and estimate the sound power level Lw using the calculated correction values.
-As described above, the present inventors have a small error with respect to the true value at a frequency lower than the specific frequency F1 in the estimation method of the second embodiment, and the specific frequency F1 in the estimation method of the third embodiment. We have found that the error with respect to the true value is small at higher frequencies. That is, it was found that the appropriate correction method differs depending on the frequency band. Further, the present inventors have different values for the specific frequency F1 depending on the shape of the estimation target sound source 10, and based on the results of experiments and simulations conducted in advance using the reference sound source 34, for each shape of the estimation target sound source 10. Found that it is possible to set to. From this, the sound power level Lw of the estimation target sound source 10 may be estimated by a correction method different for each frequency band based on the shape of the estimation target sound source 10.

上述した推定方法においては、図9に示すように、面要素設定工程(ステップS101)に先立って、周波数帯域を区分する区分工程(ステップS100−1)が行われる。また、面要素設定工程(ステップS101)と補正値を算出する算出工程(例えばステップS110)との間に、周波数帯域ごとの補正方法を選択する選択工程(ステップS100−2)が行われる。 In the estimation method described above, as shown in FIG. 9, a division step (step S100-1) for dividing the frequency band is performed prior to the surface element setting step (step S101). Further, a selection step (step S100-2) for selecting a correction method for each frequency band is performed between the surface element setting step (step S101) and the calculation step for calculating the correction value (for example, step S110).

区分工程(ステップS101−1)では、まず、基準音源34を仮想的な多面体に見立て、その多面体を構成する各面を推定面要素ZiとしてマイクロホンMijを設置する。次に、そのマイクロホンMijを用いて半無響空間や閉空間で基準音源34の音圧レベルを測定し、各種の補正方法を用いて音響パワーレベルLwを算出する。次に、算出した音響パワーレベルLwと基準音源34の既知音響パワーレベルLwKとを比較し、各種の補正方法に適した周波数帯域を特定する。 In the division step (step S101-1), first, the reference sound source 34 is regarded as a virtual polyhedron, and the microphone Mij is installed with each surface constituting the polyhedron as the estimated surface element Zi. Next, the sound pressure level of the reference sound source 34 is measured in a semi-anechoic space or a closed space using the microphone Mij, and the sound power level Lw is calculated using various correction methods. Next, the calculated sound power level Lw is compared with the known sound power level LwK of the reference sound source 34, and a frequency band suitable for various correction methods is specified.

基準音源34を様々な形状の多面体に見立てて周波数帯域の特定を行うことで、推定対象音源10の形状に対して、各種補正方法とその補正方法に適正な周波数帯域とを関連付けることができる。そして、その関連付けた情報である適正情報は、推定装置25の記憶部30に予め記憶される。なお、こうした適正情報は、ある程度のサンプルがとれたあと、多面体の形状を入力、補正方法およびその補正方法に適正な周波数帯域を出力とする機械学習によって構築されたモデルで構成されていてもよい。 By identifying the frequency band by simulating the reference sound source 34 as a polyhedron having various shapes, it is possible to associate various correction methods and an appropriate frequency band for the correction method with respect to the shape of the estimation target sound source 10. Then, the appropriate information, which is the associated information, is stored in advance in the storage unit 30 of the estimation device 25. It should be noted that such appropriate information may be composed of a model constructed by machine learning in which the shape of the polyhedron is input, the correction method and the frequency band appropriate for the correction method are output after a certain amount of sample is taken. ..

選択工程(ステップS100−2)においては、入力部21を通じて、推定対象音源10の形状を示す形状情報が推定装置25に入力される。推定部29は、その形状情報と記憶部30に記憶している適正情報とに基づいて、適正な補正方法、すなわち補正値算出部による適正な補正値の算出方法を周波数帯域ごとに選択する。 In the selection step (step S100-2), shape information indicating the shape of the estimation target sound source 10 is input to the estimation device 25 through the input unit 21. The estimation unit 29 selects an appropriate correction method, that is, an appropriate correction value calculation method by the correction value calculation unit for each frequency band, based on the shape information and the appropriate information stored in the storage unit 30.

この構成によれば、例えば、図10に示すように、特定周波数F1よりも低い周波数においては第2実施形態の推定方法を用い、特定周波数F1よりも高い周波数においては第3実施形態の推定方法を用いて音響パワーレベルLwを算出できる。これにより、音響パワーレベルLwの推定値と音響パワーレベルLwの真値との誤差をさらに小さくできる。なお、特定周波数F1における音響パワーレベルLwは、第2実施形態の推定方法を用いて算出された値、第3実施形態の推定方法を用いて算出された値、これらの値の中間値のいずれであってもよい。 According to this configuration, for example, as shown in FIG. 10, the estimation method of the second embodiment is used at a frequency lower than the specific frequency F1, and the estimation method of the third embodiment is used at a frequency higher than the specific frequency F1. Can be used to calculate the sound power level Lw. As a result, the error between the estimated value of the sound power level Lw and the true value of the sound power level Lw can be further reduced. The sound power level Lw at the specific frequency F1 is any of a value calculated by using the estimation method of the second embodiment, a value calculated by using the estimation method of the third embodiment, and an intermediate value between these values. May be.

・再現工程(ステップS107)および再現取得工程(ステップS108)は、コンピュータによるシミュレーションによって行われてもよい。この場合、再現工程では、位置決め工程において決定された設置位置の座標に基づいてマイクロホンMによる測定環境が再現される。また、再現取得工程は、例えば境界要素法や有限要素法などによる数値シミュレーションによって、各マイクロホンMが測定する基準音源の音圧レベルを推定する。この数値シミュレーションは推定面要素Ziごとに行われるため、コンピュータに対して多大な演算負荷がかかることもない。 The reproduction step (step S107) and the reproduction acquisition step (step S108) may be performed by computer simulation. In this case, in the reproduction step, the measurement environment by the microphone M is reproduced based on the coordinates of the installation position determined in the positioning step. Further, in the reproduction acquisition step, the sound pressure level of the reference sound source measured by each microphone M is estimated by numerical simulation by, for example, the boundary element method or the finite element method. Since this numerical simulation is performed for each estimation surface element Zi, a large computational load is not applied to the computer.

M…マイクロホン、10…推定対象音源、12…閉空間、20…推定システム、21…入力部、25…推定装置、26…取得部、27…第1補正値算出部、28…第2補正値算出部、29…推定部、30…記憶部、31…対応付けデータ、32…床面、33…半無響空間、34…基準音源、35…基準音源データ、36…再現環境測定データ、37…第1補正データ、41…基準音源測定データ、42…第2補正データ、43…推定用測定データ。 M ... Microphone, 10 ... Estimated sound source, 12 ... Closed space, 20 ... Estimate system, 21 ... Input unit, 25 ... Estimator, 26 ... Acquisition unit, 27 ... First correction value calculation unit, 28 ... Second correction value Calculation unit, 29 ... estimation unit, 30 ... storage unit, 31 ... associated data, 32 ... floor surface, 33 ... semi-insensitive space, 34 ... reference sound source, 35 ... reference sound source data, 36 ... reproduction environment measurement data, 37 ... 1st correction data, 41 ... Reference sound source measurement data, 42 ... 2nd correction data, 43 ... Measurement data for estimation.

Claims (4)

推定対象音源の周囲に設置された複数の電気音響変換器で前記推定対象音源の音圧レベルを測定し、前記測定した音圧レベルに基づいて前記推定対象音源の音響パワーレベルを推定する音響パワーレベルの推定方法であって、
前記複数の電気音響変換器の設置位置を決定する位置決め工程と、
音響パワーレベルが既知である基準音源を音源として半無響空間に前記推定対象音源の測定環境を再現する再現工程と、
前記再現された測定環境において前記複数の電気音響変換器が測定した前記基準音源の音圧レベルを取得する再現取得工程と、
前記取得した音圧レベルに基づいて、前記推定対象音源の音響パワーレベルを推定する際の補正値を算出する算出工程と、
を備える音響パワーレベルの推定方法。
Sound power of the estimation target sound source is measured by a plurality of electroacoustic converters installed around the estimation target sound source, and the sound power level of the estimation target sound source is estimated based on the measured sound pressure level. It ’s a level estimation method.
The positioning step of determining the installation position of the plurality of electroacoustic converters and
A reproduction process that reproduces the measurement environment of the estimation target sound source in a semi-anechoic space using a reference sound source with a known sound power level as a sound source.
A reproduction acquisition step of acquiring the sound pressure level of the reference sound source measured by the plurality of electroacoustic converters in the reproduced measurement environment, and a reproduction acquisition step.
A calculation step of calculating a correction value when estimating the sound power level of the estimation target sound source based on the acquired sound pressure level, and a calculation step.
A method of estimating sound power level.
前記推定対象音源の音響パワーレベルは、前記推定対象音源の形状に基づいて分割した複数の推定面要素ごとの部分音響パワーレベルの平均値により推定される値であり、
前記複数の推定面要素の各々には、前記複数の電気音響変換器のうちの少なくとも1つが対応付けられており、
前記部分音響パワーレベルは、前記対応付けられた電気音響変換器が測定した音圧レベルの平均値である平均音圧レベルを補正した値であり、
前記再現取得工程では、前記推定面要素ごとに、前記対応付けられた電気音響変換器が測定した前記基準音源の音圧レベルを取得し、
前記算出工程では、前記推定面要素ごとに、前記再現取得工程で取得した音圧レベルの平均値と前記基準音源の音響パワーレベルとの差を、前記平均音圧レベルを補正する補正値として算出する
請求項1に記載の音響パワーレベルの推定方法。
The sound power level of the estimation target sound source is a value estimated by the average value of the partial sound power levels of each of the plurality of estimation surface elements divided based on the shape of the estimation target sound source.
Each of the plurality of estimation surface elements is associated with at least one of the plurality of electroacoustic transducers.
The partial acoustic power level is a value obtained by correcting the average sound pressure level, which is the average value of the sound pressure levels measured by the associated electroacoustic converter.
In the reproduction acquisition step, the sound pressure level of the reference sound source measured by the associated electroacoustic transducer is acquired for each estimated surface element.
In the calculation step, the difference between the average value of the sound pressure level acquired in the reproduction acquisition step and the sound power level of the reference sound source is calculated as a correction value for correcting the average sound pressure level for each estimated surface element. The method for estimating the sound power level according to claim 1.
前記測定した音圧レベルは、前記電気音響変換器の測定値を補正した値であり、
前記推定対象音源の測定環境で前記複数の電気音響変換器の各々が測定する前記基準音源の音圧レベルを取得する基準音源取得工程をさらに備え、
前記算出工程では、前記複数の電気音響変換器ごとに、前記基準音源取得工程で取得した音圧レベルと前記再現取得工程で取得した音圧レベルとの差を、前記測定値を補正する補正値として算出する
請求項1または2に記載の音響パワーレベルの推定方法。
The measured sound pressure level is a value obtained by correcting the measured value of the electroacoustic converter.
Further provided with a reference sound source acquisition step of acquiring the sound pressure level of the reference sound source measured by each of the plurality of electroacoustic converters in the measurement environment of the estimation target sound source.
In the calculation step, for each of the plurality of electroacoustic converters, the difference between the sound pressure level acquired in the reference sound source acquisition step and the sound pressure level acquired in the reproduction acquisition step is a correction value for correcting the measured value. The method for estimating the acoustic power level according to claim 1 or 2.
前記複数の電気音響変換器が測定する音圧レベルの周波数を複数の周波数帯域に区分けする区分工程と、
前記複数の周波数帯域の各々について、前記音響パワーレベルを算出する際の補正方法を選択する選択工程と、をさらに備える
請求項1〜3のいずれか一項に記載の音響パワーレベルの推定方法。
A division step of dividing the frequency of the sound pressure level measured by the plurality of electroacoustic converters into a plurality of frequency bands, and a division step.
The method for estimating a sound power level according to any one of claims 1 to 3, further comprising a selection step of selecting a correction method for calculating the sound power level for each of the plurality of frequency bands.
JP2020088613A 2020-05-21 2020-05-21 How to estimate sound power level Active JP7387530B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2020088613A JP7387530B2 (en) 2020-05-21 2020-05-21 How to estimate sound power level

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2020088613A JP7387530B2 (en) 2020-05-21 2020-05-21 How to estimate sound power level

Publications (2)

Publication Number Publication Date
JP2021183910A true JP2021183910A (en) 2021-12-02
JP7387530B2 JP7387530B2 (en) 2023-11-28

Family

ID=78767280

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2020088613A Active JP7387530B2 (en) 2020-05-21 2020-05-21 How to estimate sound power level

Country Status (1)

Country Link
JP (1) JP7387530B2 (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS63153433A (en) * 1986-12-17 1988-06-25 Fuji Electric Co Ltd Acoustic power level measuring method
JPH0843189A (en) * 1994-07-27 1996-02-16 Ono Sokki Co Ltd Method and apparatus for estimating sound field
JP2002365127A (en) * 2001-06-06 2002-12-18 Mitsubishi Heavy Ind Ltd Method and device for estimating acoustic power level of sound source
JP2013511712A (en) * 2009-11-19 2013-04-04 ウニフェルジテイト・トゥウェンテ Method and apparatus for determining acoustic coefficients and power

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS63153433A (en) * 1986-12-17 1988-06-25 Fuji Electric Co Ltd Acoustic power level measuring method
JPH0843189A (en) * 1994-07-27 1996-02-16 Ono Sokki Co Ltd Method and apparatus for estimating sound field
JP2002365127A (en) * 2001-06-06 2002-12-18 Mitsubishi Heavy Ind Ltd Method and device for estimating acoustic power level of sound source
JP2013511712A (en) * 2009-11-19 2013-04-04 ウニフェルジテイト・トゥウェンテ Method and apparatus for determining acoustic coefficients and power

Also Published As

Publication number Publication date
JP7387530B2 (en) 2023-11-28

Similar Documents

Publication Publication Date Title
KR100938691B1 (en) Active noise suppressor
EP3631388B1 (en) A vehicle noise, vibration and harshness (nvh) testing method and system
JP2020511003A5 (en)
JP2010277577A (en) Real-time scheduling of linear model for control and estimation
JP6647992B2 (en) Design support equipment
JP2008014709A (en) Image processing apparatus
CN115278496A (en) Sparse sound source identification method and system for microphone array measurement
US9668075B2 (en) Estimating parameter values for a lumped parameter model of a loudspeaker
JP7387530B2 (en) How to estimate sound power level
US20170243111A1 (en) Neural network-based parameter estimation of loudspeakers
US10979833B2 (en) Acoustical performance evaluation method
JP5852935B2 (en) Transfer function estimation device, transfer function estimation method, and transfer function estimation program
CN115902773A (en) Sound source positioning method and system suitable for different microphone arrays
US9723419B2 (en) Systems and methods for determining metric for sound system evaluation
JP7488952B2 (en) Methods for reducing vehicle pass-by noise
JP2007309785A (en) Characteristic vibration mode extraction method, device, and program
JP2023024023A (en) Noise source contribution analysis device, noise source contribution analysis method and noise source contribution analysis program
JP2011139409A (en) Audio signal processor, audio signal processing method, and computer program
JP2002365127A (en) Method and device for estimating acoustic power level of sound source
JP4676920B2 (en) Signal separation device, signal separation method, signal separation program, and recording medium
JP7118626B2 (en) System, method and program
CN113903351A (en) Echo cancellation method, device, equipment and storage medium
CN113639850B (en) Sound source load test method, system, equipment and medium
JP6716933B2 (en) Noise estimation device, program and method, and voice processing device
JPH11148875A (en) Residual stress analyzing apparatus

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20230208

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20231027

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: 20231107

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20231115

R150 Certificate of patent or registration of utility model

Ref document number: 7387530

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150