JP6087850B2 - Acoustic transfer characteristic measuring device, acoustic transfer characteristic measuring method, program - Google Patents

Acoustic transfer characteristic measuring device, acoustic transfer characteristic measuring method, program Download PDF

Info

Publication number
JP6087850B2
JP6087850B2 JP2014014046A JP2014014046A JP6087850B2 JP 6087850 B2 JP6087850 B2 JP 6087850B2 JP 2014014046 A JP2014014046 A JP 2014014046A JP 2014014046 A JP2014014046 A JP 2014014046A JP 6087850 B2 JP6087850 B2 JP 6087850B2
Authority
JP
Japan
Prior art keywords
signal
stsp
spatio
temporal
transfer characteristic
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2014014046A
Other languages
Japanese (ja)
Other versions
JP2015142248A (en
Inventor
翔一 小山
翔一 小山
島内　末廣
末廣 島内
仲 大室
仲 大室
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.)
Nippon Telegraph and Telephone Corp
Original Assignee
Nippon Telegraph and Telephone Corp
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 Nippon Telegraph and Telephone Corp filed Critical Nippon Telegraph and Telephone Corp
Priority to JP2014014046A priority Critical patent/JP6087850B2/en
Publication of JP2015142248A publication Critical patent/JP2015142248A/en
Application granted granted Critical
Publication of JP6087850B2 publication Critical patent/JP6087850B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Obtaining Desirable Characteristics In Audible-Bandwidth Transducers (AREA)
  • Circuit For Audible Band Transducer (AREA)

Description

この発明は、スピーカアレイとマイクロホンアレイを用いて、それらの間の音響伝達特性を測定する技術に関する。   The present invention relates to a technique for measuring a sound transfer characteristic between a speaker array and a microphone array.

音響伝達特性を測定する従来技術にTime Stretched Pulse（ＴＳＰ）信号を用いるものがある（非特許文献１）。ＴＳＰ信号とは、パルス信号を時間方向に伸長した信号である。非特許文献１の音響伝達特性の測定技術は、ＴＳＰ信号を１個のスピーカから出力し、１個のマイクロホンで収音することで、単一のスピーカと単一のマイクロホンの間の伝達特性を、高精度かつ頑健に測定することを可能にする技術である。この技術を多数のスピーカ・マイクロホン間で適用する場合には、各スピーカからＴＳＰ信号を出力し、マイクロホンアレイで観測した信号に対し、逆ＴＳＰ信号を畳み込むことにより、音響伝達特性を算出する。   A conventional technique for measuring sound transfer characteristics uses a Time Stretched Pulse (TSP) signal (Non-Patent Document 1). The TSP signal is a signal obtained by extending a pulse signal in the time direction. The measurement technique of acoustic transfer characteristics of Non-Patent Document 1 is that the TSP signal is output from one speaker and picked up by one microphone, so that the transfer characteristic between a single speaker and a single microphone is obtained. It is a technology that enables highly accurate and robust measurement. When this technology is applied between a large number of speakers and microphones, a TSP signal is output from each speaker, and an acoustic transfer characteristic is calculated by convolving an inverse TSP signal with a signal observed by a microphone array.

Y. Suzuki, F. Asano, H. Kim, and T. Sone. “An optimal computer-generated pulse signal suitable for the the measurement of very long impulse responses,” J. Acoust. Soc. Am., vol. 97, no. 2, pp.1119-1123, 1995.Y. Suzuki, F. Asano, H. Kim, and T. Sone. “An optimal computer-generated pulse signal suitable for the measurement of very long impulse responses,” J. Acoust. Soc. Am., Vol. 97, no. 2, pp.1119-1123, 1995.

従来技術では、多数のスピーカ、マイクロホン間の音響伝達特性を測定する場合において、各スピーカからＴＳＰ信号を出力する必要があるため、全スピーカからの特性を測定するためには長い時間を要するという課題があった。そこで本発明では、スピーカ・マイクロホンアレイ間の音響伝達特性の測定に要する時間を短縮することができる音響伝達特性測定装置を提供することを目的とする。   In the prior art, when measuring the acoustic transfer characteristics between a large number of speakers and microphones, it is necessary to output a TSP signal from each speaker. Therefore, it takes a long time to measure the characteristics from all the speakers. was there. Therefore, an object of the present invention is to provide an acoustic transfer characteristic measuring apparatus that can shorten the time required for measuring the acoustic transfer characteristic between a speaker and a microphone array.

本発明の音響伝達特性測定装置は、スピーカアレイと、マイクロホンアレイと、ＳＴＳＰ信号計算部と、逆ＳＴＳＰ信号計算部と、伝達特性推定値算出部を含む。   The acoustic transfer characteristic measuring apparatus of the present invention includes a speaker array, a microphone array, an STSP signal calculation unit, an inverse STSP signal calculation unit, and a transfer characteristic estimated value calculation unit.

スピーカアレイは、複数のスピーカを直線状に配置してなる。マイクロホンアレイは、複数のマイクロホンを直線状に配置してなる。ＳＴＳＰ信号計算部は、合成時の平面波の角度と周波数とが時間方向に変化する信号であるＳＴＳＰ信号の時空間領域の表現である時空間領域ＳＴＳＰ信号を計算してスピーカアレイから出力する。逆ＳＴＳＰ信号計算部は、時空間領域ＳＴＳＰ信号に基づいて時空間領域逆ＳＴＳＰ信号を計算する。伝達特性推定値算出部は、マイクロホンアレイに収音される信号である入力信号と、時空間領域逆ＳＴＳＰ信号に基づいて時空間音響伝達特性の推定値を算出する。   The speaker array is formed by arranging a plurality of speakers in a straight line. The microphone array is formed by arranging a plurality of microphones in a straight line. The STSP signal calculation unit calculates a spatio-temporal region STSP signal that is a representation of the spatio-temporal region of the STSP signal, which is a signal in which the angle and frequency of the plane wave at the time of synthesis change in the time direction, and outputs it from the speaker array. The inverse STSP signal calculation unit calculates a spatiotemporal domain inverse STSP signal based on the spatiotemporal domain STSP signal. The transfer characteristic estimated value calculation unit calculates an estimated value of the spatio-temporal acoustic transfer characteristic based on the input signal that is a signal collected by the microphone array and the spatio-temporal domain inverse STSP signal.

本発明の音響伝達特性測定装置によれば、スピーカ・マイクロホンアレイ間の音響伝達特性の測定に要する時間を短縮することができる。   According to the acoustic transfer characteristic measuring apparatus of the present invention, the time required for measuring the acoustic transfer characteristic between the speaker and the microphone array can be shortened.

本発明の実施例１の音響伝達特性測定装置の構成を示すブロック図。The block diagram which shows the structure of the acoustic transfer characteristic measuring apparatus of Example 1 of this invention. 本発明の実施例１の音響伝達特性測定装置の動作を示すフローチャート。The flowchart which shows operation | movement of the acoustic transfer characteristic measuring apparatus of Example 1 of this invention. 本発明の実施例２の音響伝達特性測定装置の構成を示すブロック図。The block diagram which shows the structure of the acoustic transfer characteristic measuring apparatus of Example 2 of this invention. 本発明の実施例２の音響伝達特性測定装置の動作を示すフローチャート。The flowchart which shows operation | movement of the acoustic transfer characteristic measuring apparatus of Example 2 of this invention.

以下、本発明の実施の形態について、詳細に説明する。なお、同じ機能を有する構成部には同じ番号を付し、重複説明を省略する。   Hereinafter, embodiments of the present invention will be described in detail. In addition, the same number is attached | subjected to the structure part which has the same function, and duplication description is abbreviate | omitted.

＜Space-Time Stretched Pulse（ＳＴＳＰ）信号について＞
以下、本発明で用いられるSpace-Time Stretched Pulse（ＳＴＳＰ）信号について定義する。ＳＴＳＰ信号とは、簡単にはＴＳＰ信号を空間方向に伸長させた信号のことである。従って、ＳＴＳＰ信号は、パルス信号を時間方向および空間方向に伸長させた信号であり、合成時の平面波の角度と周波数とが時間方向に変化する信号であるとも表現できる。 <About Space-Time Stretched Pulse (STSP) signal>
Hereinafter, the Space-Time Stretched Pulse (STSP) signal used in the present invention will be defined. The STSP signal is simply a signal obtained by expanding the TSP signal in the spatial direction. Therefore, the STSP signal is a signal obtained by extending a pulse signal in the time direction and the spatial direction, and can also be expressed as a signal in which the angle and frequency of the plane wave at the time of synthesis change in the time direction.

以下、図１、図２を参照して実施例１の音響伝達特性測定装置について説明する。図１は、本実施例の音響伝達特性測定装置１の構成を示すブロック図である。図２は、本実施例の音響伝達特性測定装置１の動作を示すフローチャートである。   Hereinafter, the acoustic transfer characteristic measuring apparatus according to the first embodiment will be described with reference to FIGS. 1 and 2. FIG. 1 is a block diagram illustrating a configuration of an acoustic transfer characteristic measuring apparatus 1 according to the present embodiment. FIG. 2 is a flowchart showing the operation of the sound transfer characteristic measuring apparatus 1 of the present embodiment.

図１に示すように、本実施例の音響伝達特性測定装置１は、ＳＴＳＰ信号計算部１１と、スピーカアレイ１２と、逆ＳＴＳＰ信号計算部１３と、マイクロホンアレイ１４と、伝達特性推定値算出部１５を含む。ＳＴＳＰ信号計算部１１は、時空間周波数領域ＳＴＳＰ信号計算部１１１と、時空間領域ＳＴＳＰ信号計算部１１２を含む。逆ＳＴＳＰ信号計算部１３は、時空間周波数領域逆ＳＴＳＰ信号計算部１３１と、時空間領域逆ＳＴＳＰ信号計算部１３２を含む。   As shown in FIG. 1, the acoustic transfer characteristic measuring apparatus 1 of the present embodiment includes an STSP signal calculation unit 11, a speaker array 12, an inverse STSP signal calculation unit 13, a microphone array 14, and a transfer characteristic estimated value calculation unit. 15 is included. The STSP signal calculator 11 includes a spatio-temporal frequency domain STSP signal calculator 111 and a spatio-temporal domain STSP signal calculator 112. The inverse STSP signal calculator 13 includes a spatiotemporal frequency domain inverse STSP signal calculator 131 and a spatiotemporal domain inverse STSP signal calculator 132.

スピーカアレイ１２は、複数のスピーカ１２−１，１２−２，…，１２−Ｎｘを直線状に配置してなる。マイクロホンアレイ１４は、複数のマイクロホン１４−１，１４−２，…，１４−Ｎｘを直線状に配置してなる。ただし、Nxをスピーカ・マイクロホン数とし、Nxは２以上の任意の整数とする。 Speaker array 12 includes a plurality of speakers 12-1, 12-2, ..., formed by arranging in a straight line to 12-N x. Microphone array 14, a plurality of microphones 14-1 and 14-2, ..., formed by arranging in a straight line to 14-N x. However, N x is the number of speakers and microphones, and N x is an arbitrary integer of 2 or more.

＜ＳＴＳＰ信号計算部１１＞
以下、ＳＴＳＰ信号計算部１１について説明する。ＳＴＳＰ信号計算部１１の時空間周波数領域ＳＴＳＰ信号計算部１１１は、前述のＳＴＳＰ信号の時空間周波数領域の表現である時空間周波数領域ＳＴＳＰ信号を、時間周波数と空間周波数とが線形変化する条件で計算する（Ｓ１１１）。より詳細には、時空間周波数領域ＳＴＳＰ信号計算部１１１は、ＳＴＳＰ信号設計パラメータ(Δω,ω0,Δkx,kx0)を入力とし、時空間周波数領域ＳＴＳＰ信号F(kx,n,ωm)を、例えば以下の式（１）で計算する。 <STSP signal calculator 11>
Hereinafter, the STSP signal calculation unit 11 will be described. The spatio-temporal frequency domain STSP signal calculator 111 of the STSP signal calculator 11 converts the spatio-temporal frequency domain STSP signal, which is a representation of the above-mentioned spatio-temporal frequency domain of the STSP signal, under the condition that the temporal frequency and the spatial frequency change linearly. Calculate (S111). More specifically, the spatio-temporal frequency domain STSP signal calculation unit 111 receives the STSP signal design parameters (Δω, ω 0 , Δk x , k x0 ) as input, and the spatio-temporal frequency domain STSP signal F (k x, n , ω m ) is calculated by the following equation (1), for example.

Figure 0006087850
Figure 0006087850

ここで、kx,nは空間周波数、ωmは時間周波数、Δωは時間周波数の増加量を調整するパラメータ、Δkxは空間周波数の増加量を調整するパラメータ、ω0は時間周波数の初期値を調整するパラメータ、kx0は空間周波数の初期値を調整するパラメータ、n,mはそれぞれ空間周波数kx,n，時間周波数ωmのインデックスを示しており、1≦n≦Nkx,1≦m≦Nωである。NkxとNωはそれぞれ、空間周波数ビン数、時間周波数ビン数である。空間周波数kx,n，時間周波数ωmは例えば以下のように定義される。 Where k x, n is the spatial frequency, ω m is the time frequency, Δω is a parameter for adjusting the increase amount of the time frequency, Δk x is a parameter for adjusting the increase amount of the spatial frequency, and ω 0 is the initial value of the time frequency. K x0 is a parameter for adjusting the initial value of the spatial frequency, n and m are indices of the spatial frequency k x, n and the time frequency ω m , respectively, 1 ≦ n ≦ N kx , 1 ≦ m ≦ N ω . N kx and N ω are the spatial frequency bin number and the temporal frequency bin number, respectively. The spatial frequency k x, n and the temporal frequency ω m are defined as follows, for example.

Figure 0006087850
Figure 0006087850

ここで、Δx,Δtはそれぞれ空間（スピーカ・マイクロホン）および時間サンプリング間隔である。次に、ＳＴＳＰ信号計算部１１の時空間領域ＳＴＳＰ信号計算部１１２は、時空間周波数領域ＳＴＳＰ信号に基づいて、ＳＴＳＰ信号の時空間領域の表現である時空間領域ＳＴＳＰ信号を計算してスピーカアレイ１２から出力する（Ｓ１１２）。より詳細には、時空間領域ＳＴＳＰ信号計算部１１２は、前述した時空間周波数領域ＳＴＳＰ信号F(kx,n,ωm)に対して、例えば２次元逆離散フーリエ変換を適用することにより、時空間領域ＳＴＳＰ信号f(xi,tj)を計算する。時空間領域ＳＴＳＰ信号計算部１１２は、例えば以下の式（２）で定義される２次元逆離散フーリエ変換を実行して、時空間領域ＳＴＳＰ信号f(xi,tj)を計算する。 Here, Δx and Δt are space (speaker / microphone) and time sampling interval, respectively. Next, the spatio-temporal region STSP signal calculator 112 of the STSP signal calculator 11 calculates a spatio-temporal region STSP signal that is a representation of the spatio-temporal region of the STSP signal based on the spatio-temporal frequency region STSP signal. 12 (S112). More specifically, the spatio-temporal domain STSP signal calculation unit 112 applies, for example, a two-dimensional inverse discrete Fourier transform to the above-described spatio-temporal frequency domain STSP signal F (k x, n , ω m ). A space-time domain STSP signal f (x i , t j ) is calculated. The spatio-temporal domain STSP signal calculation unit 112 calculates a spatio-temporal domain STSP signal f (x i , t j ) by executing, for example, a two-dimensional inverse discrete Fourier transform defined by the following equation (2).

Figure 0006087850
Figure 0006087850

ただし、i,jは空間方向（スピーカ方向）xiおよび時間方向tjのインデックスを示す。前述したように、Nxをスピーカ・マイクロホン数とし、Ntを時間方向の信号長とすると、インデクスi,jは1≦i≦Nx，1≦j≦Ntである。時空間領域ＳＴＳＰ信号計算部１１２は、この時空間領域ＳＴＳＰ信号f(xi,tj)をスピーカアレイ１２を構成する各スピーカ１２−１，１２−２，…，１２−Ｎｘから出力する。上述した、ＳＴＳＰ信号設計パラメータ(Δω,ω0,Δkx,kx0)、及びその他のパラメータの具体的な値として、例えばΔω=400, ω0=0, Δkx=4000, kx0=0, Nx=512, Nt=512, Δx=0.06m, Δt=0.625ms（サンプリング周波数16kHz）を用いることができる。 Here, i and j indicate indexes in the spatial direction (speaker direction) x i and the time direction t j . As described above, the N x is the number of speaker microphones, when the N t and the time direction of the signal length is an index i, j is 1 ≦ i ≦ N x, 1 ≦ j ≦ N t. The spatio-temporal region STSP signal calculation unit 112 outputs the spatio-temporal region STSP signal f (x i , t j ) from each speaker 12-1, 12-2,..., 12-N x constituting the speaker array 12. . As specific values of the above-described STSP signal design parameters (Δω, ω 0 , Δk x , k x0 ) and other parameters, for example, Δω = 400, ω 0 = 0, Δk x = 4000, k x0 = 0 , N x = 512, N t = 512, Δx = 0.06 m, Δt = 0.625 ms (sampling frequency 16 kHz).

＜逆ＳＴＳＰ信号計算部１３＞
以下、逆ＳＴＳＰ信号計算部１３について説明する。逆ＳＴＳＰ信号計算部１３の時空間周波数領域逆ＳＴＳＰ信号計算部１３１は、時空間周波数領域ＳＴＳＰ信号に畳み込むことでパルス信号を生成する信号である時空間周波数領域逆ＳＴＳＰ信号を計算する（Ｓ１３１）。より詳細には、時空間周波数領域逆ＳＴＳＰ信号計算部１３１は、ＳＴＳＰ信号設計パラメータ(Δω,ω0,Δkx,kx0)を入力とし、例えば、以下の式（３）に基づき、時空間周波数領域逆ＳＴＳＰ信号G(kx,n,ωm)を計算する。 <Inverse STSP signal calculator 13>
Hereinafter, the inverse STSP signal calculation unit 13 will be described. The spatio-temporal frequency domain inverse STSP signal calculator 131 of the inverse STSP signal calculator 13 calculates a spatio-temporal frequency domain inverse STSP signal that is a signal for generating a pulse signal by convolution with the spatio-temporal frequency domain STSP signal (S131). . More specifically, the spatio-temporal frequency domain inverse STSP signal calculation unit 131 receives the STSP signal design parameters (Δω, ω 0 , Δk x , k x0 ) as an input, for example, based on the following equation (3), The frequency domain inverse STSP signal G (k x, n , ω m ) is calculated.

Figure 0006087850
Figure 0006087850

逆ＳＴＳＰ信号計算部１３の時空間領域逆ＳＴＳＰ信号計算部１３２は、時空間周波数領域逆ＳＴＳＰ信号に基づいて時空間領域逆ＳＴＳＰ信号を計算する（Ｓ１３２）。より詳細には、時空間領域逆ＳＴＳＰ信号計算部１３２は、時空間周波数領域逆ＳＴＳＰ信号G(kx,n,ωm)に対して、例えば２次元逆離散フーリエ変換を適用することにより、時空間領域逆ＳＴＳＰ信号g(xi,tj)を計算する。時空間領域逆ＳＴＳＰ信号計算部１３２は、例えば以下の式（４）で定義される２次元逆離散フーリエ変換を実行して、時空間領域逆ＳＴＳＰ信号g(xi,tj)を計算する。 The spatiotemporal domain inverse STSP signal calculator 132 of the inverse STSP signal calculator 13 calculates a spatiotemporal domain inverse STSP signal based on the spatiotemporal frequency domain inverse STSP signal (S132). More specifically, the spatio-temporal domain inverse STSP signal calculation unit 132 applies, for example, a two-dimensional inverse discrete Fourier transform to the spatio-temporal frequency domain inverse STSP signal G (k x, n , ω m ). A space-time domain inverse STSP signal g (x i , t j ) is calculated. The spatio-temporal domain inverse STSP signal calculation unit 132 performs, for example, a two-dimensional inverse discrete Fourier transform defined by the following equation (4) to calculate a spatio-temporal domain inverse STSP signal g (x i , t j ). .

Figure 0006087850
Figure 0006087850

＜伝達特性推定値算出部１５＞
以下、伝達特性推定値算出部１５について説明する。伝達特性推定値算出部１５は、マイクロホンアレイに収音される信号である入力信号と、時空間領域逆ＳＴＳＰ信号に基づいて時空間音響伝達特性の推定値を算出する（Ｓ１５）。より詳細には、マイクロホンアレイ１４に収音される信号である入力信号（マイクロホンアレイ入力信号ともいう）であるf~(xi,tj)と時空間領域逆ＳＴＳＰ信号g(xi,tj)を入力とし、時空間音響伝達特性の推定値h~(xi,tj)（または、時空間周波数領域の表現であるH~(kx,n,ωm)）を出力とする。伝達特性推定値算出部１５が出力する時空間音響伝達特性の推定値は、時空間領域の表現であってもよいし、時空間周波数領域の表現であってもよい。例えば、伝達特性推定値算出部１５は、時空間音響伝達特性の推定値をマイクロホンアレイ入力信号f~(xi,tj)と時空間領域逆ＳＴＳＰ信号g(xi,tj)の畳み込みによって算出する。伝達特性推定値算出部１５は、例えば二次元の円状畳み込みを用いる方法、または二次元の線状畳み込みを用いる方法、または円状畳み込みと線状畳み込みを組み合わせる方法によって、時空間音響伝達特性の推定値を算出できる。以下に、これら３つの算出方法の一例を示す。 <Transfer characteristic estimated value calculation unit 15>
Hereinafter, the transfer characteristic estimated value calculation unit 15 will be described. The transfer characteristic estimated value calculation unit 15 calculates an estimated value of the spatiotemporal acoustic transfer characteristic based on the input signal that is a signal picked up by the microphone array and the spatiotemporal domain inverse STSP signal (S15). More specifically, the input signals (also referred to as microphone array input signals) f ~ (x i , t j ) that are signals collected by the microphone array 14 and the spatio-temporal domain inverse STSP signal g (x i , t j ) as input and spatio-temporal acoustic transfer characteristic estimation value h ~ (x i , t j ) (or H ~ (k x, n , ω m ) representing spatio-temporal frequency domain) as output . The estimated value of the spatio-temporal acoustic transfer characteristic output by the transfer characteristic estimated value calculation unit 15 may be a spatio-temporal domain expression or a spatio-temporal frequency domain expression. For example, the transfer characteristic estimated value calculation unit 15 convolves the estimated value of the spatio-temporal acoustic transfer characteristic with the microphone array input signal f˜ (x i , t j ) and the spatio-temporal domain inverse STSP signal g (x i , t j ). Calculated by The transfer characteristic estimation value calculation unit 15 performs the spatio-temporal acoustic transfer characteristic by using, for example, a method using two-dimensional circular convolution, a method using two-dimensional linear convolution, or a method combining circular convolution and linear convolution. An estimated value can be calculated. An example of these three calculation methods is shown below.

＜二次元円状畳み込みを用いる方法＞
二次元円状畳み込みを用いる場合、伝達特性推定値算出部１５は、マイクロホンアレイ入力信号f~(xi,tj)の２次元離散フーリエ変換であるF~(kx,n,ωm)と、時空間領域逆ＳＴＳＰ信号g(xi,tj)の２次元離散フーリエ変換であるG(kx,n,ωm)を計算する。伝達特性推定値算出部１５は、例えば以下の式（５）で定義される２次元離散フーリエ変換を実行してF~(kx,n,ωm)を計算する。 <Method using two-dimensional circular convolution>
When using a two-dimensional circular convolution, the transfer characteristic estimation value calculation unit 15 performs F˜ (k x, n , ω m ) which is a two-dimensional discrete Fourier transform of the microphone array input signal f˜ (x i , t j ). And G (k x, n , ω m ) , which is a two-dimensional discrete Fourier transform of the spatio-temporal domain inverse STSP signal g (x i , t j ). The transfer characteristic estimated value calculation unit 15 calculates F˜ (k x, n , ω m ) by executing a two-dimensional discrete Fourier transform defined by the following equation (5), for example.

Figure 0006087850
Figure 0006087850

これらを乗算した信号が、時空間周波数領域での時空間音響伝達特性の推定値H~(kx,n,ωm)となる。従って、伝達特性推定値算出部１５は、時空間周波数領域での時空間音響伝達特性の推定値H~(kx,n,ωm)を以下の式（６）で計算する。 A signal obtained by multiplying these becomes the estimated value H˜ (k x, n , ω m ) of the spatiotemporal acoustic transfer characteristic in the spatiotemporal frequency domain. Accordingly, the transfer characteristic estimated value calculation unit 15 calculates the estimated value H˜ (k x, n , ω m ) of the spatio-temporal acoustic transfer characteristic in the spatio-temporal frequency domain by the following equation (6).

Figure 0006087850
Figure 0006087850

伝達特性推定値算出部１５は、時空間領域での時空間音響伝達特性の推定値h~(xi,tj)を、H~(kx,n,ωm)の２次元逆離散フーリエ変換として算出できる。なお、伝達特性推定値算出部１５は、H~(kx,n,ωm)をそのまま出力しても構わない。 The transfer characteristic estimation value calculation unit 15 converts the spatio-temporal acoustic transfer characteristic estimation values h˜ (x i , t j ) in the spatio-temporal region into two-dimensional inverse discrete Fourier transforms H˜ (k x, n , ω m ). It can be calculated as a conversion. The transfer characteristic estimated value calculation unit 15 may output H˜ (k x, n , ω m ) as it is.

＜２次元線状畳み込みを用いる方法＞
２次元線状畳み込みを用いる場合、伝達特性推定値算出部１５は、マイクロホンアレイ入力信号f~(xi,tj)と時空間領域逆ＳＴＳＰ信号g(xi,tj)に対し、以下のように２次元線状畳み込みを行うことで、直接、時空間領域での時空間音響伝達特性の推定値h~(xi,tj)を算出する。従って、伝達特性推定値算出部１５は、以下の式（７）により、時空間領域での時空間音響伝達特性の推定値h~(xi,tj)を算出する。 <Method using two-dimensional linear convolution>
When using the two-dimensional linear convolution, the transfer characteristic estimation value calculation unit 15 performs the following for the microphone array input signal f˜ (x i , t j ) and the spatio-temporal domain inverse STSP signal g (x i , t j ). By performing the two-dimensional linear convolution as described above, the estimated values h˜ (x i , t j ) of the spatiotemporal acoustic transfer characteristics in the spatiotemporal region are directly calculated. Therefore, the transfer characteristic estimated value calculation unit 15 calculates the spatio-temporal acoustic transfer characteristic estimated values h˜ (x i , t j ) in the spatio-temporal region by the following equation (7).

Figure 0006087850
Figure 0006087850

χ、τはそれぞれ、空間方向、時間方向の遅延を表す。 χ and τ represent delays in the spatial direction and the time direction, respectively.

＜円状畳み込みと線状畳み込みの組み合わせ＞
ここでは、例えば空間方向xiで線状畳み込みを行い、時間方向tjで円状畳み込みを行う場合を考える。この場合、伝達特性推定値算出部１５は、マイクロホンアレイ入力信号f~(xi,tj)と時空間領域逆ＳＴＳＰ信号g(xi,tj)に対して、xi方向の両端にNx/2サンプルだけゼロづめをした信号を、それぞれf−(xi,tj),g−(xi,tj)と設定する。ただしここでは、1≦i≦2Nx,1≦j≦Ntである。伝達特性推定値算出部１５は、f−(xi,tj),g−(xi,tj)の２次元離散フーリエ変換であるF−(kx,n,ωm),G−(kx,n,ωm)を計算する。これらを乗算した信号が、時空間周波数領域での時空間音響伝達特性の推定値H~(kx,n,ωm)となる。従って、伝達特性推定値算出部１５は、以下の式（８）により、時空間周波数領域での時空間音響伝達特性の推定値H~(kx,n,ωm)を計算する。 <Combination of circular convolution and linear convolution>
Here, for example, consider a case where linear convolution is performed in the spatial direction x i and circular convolution is performed in the time direction t j . In this case, the transfer characteristic estimated value calculation unit 15 applies the microphone array input signal f˜ (x i , t j ) and the spatio-temporal domain inverse STSP signal g (x i , t j ) to both ends in the x i direction. Signals zeroed by N x / 2 samples are set as f − (x i , t j ) and g − (x i , t j ), respectively. However, here, 1 ≦ i ≦ 2N x and 1 ≦ j ≦ N t . Transfer characteristic estimate calculation unit 15, f - (x i, t j), g - (x i, t j) is a two-dimensional discrete Fourier transform of F - (k x, n, ω m), G - Calculate (k x, n , ω m ). A signal obtained by multiplying these becomes the estimated value H˜ (k x, n , ω m ) of the spatiotemporal acoustic transfer characteristic in the spatiotemporal frequency domain. Accordingly, the transfer characteristic estimated value calculation unit 15 calculates the spatio-temporal acoustic transfer characteristic estimated values H˜ (k x, n , ω m ) in the spatio-temporal frequency domain by the following equation (8).

Figure 0006087850
Figure 0006087850

伝達特性推定値算出部１５は、時空間領域での時空間音響伝達特性の推定値h~(xi,tj)は、H~(kx,n,ωm)の２次元逆離散フーリエ変換として算出できる。なお、伝達特性推定値算出部１５は、H~(kx,n,ωm)をそのまま出力しても構わない。ただし、この方法ではゼロづめした部分に相当する信号は、測定値としては扱わないものとする。 The transfer characteristic estimation value calculation unit 15 calculates the spatio-temporal acoustic transfer characteristic estimation value h˜ (x i , t j ) in the spatio-temporal region as a two-dimensional inverse discrete Fourier of H˜ (k x, n , ω m ). It can be calculated as a conversion. The transfer characteristic estimated value calculation unit 15 may output H˜ (k x, n , ω m ) as it is. However, in this method, the signal corresponding to the zeroed portion is not treated as a measurement value.

なお、上述では、周波数領域−時間領域の変換の具体的方法として離散フーリエ変換（および逆変換）を用いた。例えばステップＳ１１２，Ｓ１３２，Ｓ１５である。しかしながら、本発明において周波数領域−時間領域の変換は、離散フーリエ変換（逆変換）には限定されず、ステップＳ１１２，Ｓ１３２，Ｓ１５では、他の方法により、周波数領域−時間領域の変換を実行してもよい。   In the above description, the discrete Fourier transform (and inverse transform) is used as a specific method of the frequency domain-time domain transformation. For example, steps S112, S132, and S15. However, in the present invention, the frequency domain-time domain transformation is not limited to the discrete Fourier transformation (inverse transformation). In steps S112, S132, and S15, the frequency domain-time domain transformation is executed by another method. May be.

本実施例の音響伝達特性測定装置１によれば、ＳＴＳＰ信号計算部１１がパルス信号を時間方向だけでなく空間方向にも伸長させたＳＴＳＰ信号を直線状に配置されたスピーカアレイ１２から出力し、この信号を直線状に配置されたマイクロホンアレイ１４で収音し、伝達特性推定値算出部１５が時空間音響伝達特性の推定値をマイクロホンアレイ入力信号と時空間領域逆ＳＴＳＰ信号に基づいて算出するため、一度の測定で伝達特性を求めることが可能となり、時空間音響伝達特性の測定に要する時間を大幅に短縮することができる。   According to the acoustic transfer characteristic measuring apparatus 1 of the present embodiment, the STSP signal calculation unit 11 outputs the STSP signal obtained by extending the pulse signal not only in the time direction but also in the spatial direction from the linearly arranged speaker array 12. The signal is picked up by the microphone array 14 arranged in a straight line, and the transfer characteristic estimated value calculation unit 15 calculates the estimated value of the spatio-temporal acoustic transfer characteristic based on the microphone array input signal and the spatio-temporal domain inverse STSP signal. Therefore, the transfer characteristic can be obtained by one measurement, and the time required for measuring the spatio-temporal acoustic transfer characteristic can be greatly reduced.

なお、本実施例の音響伝達特性測定装置１は、波面合成法を用い、マイクロホンアレイで収音された音を遠隔地のスピーカアレイで忠実に再生することによって臨場感を実現する音声コミュニケーションにシステムに利用することができる。本実施例の音響伝達特性測定装置１は、波面合成法によるフィルタを生成する際に必要となる伝達特性を求めるのに有用である。   The sound transfer characteristic measuring apparatus 1 of the present embodiment uses a wavefront synthesis method, and is a system for voice communication that realizes a sense of reality by faithfully reproducing sound collected by a microphone array with a remote speaker array. Can be used. The acoustic transfer characteristic measuring apparatus 1 according to the present embodiment is useful for obtaining a transfer characteristic necessary for generating a filter by the wavefront synthesis method.

以下、図３、図４を参照して実施例２の音響伝達特性測定装置について説明する。図３は、本実施例の音響伝達特性測定装置２の構成を示すブロック図である。図４は、本実施例の音響伝達特性測定装置２の動作を示すフローチャートである。   Hereinafter, the acoustic transfer characteristic measuring apparatus according to the second embodiment will be described with reference to FIGS. 3 and 4. FIG. 3 is a block diagram showing the configuration of the acoustic transfer characteristic measuring apparatus 2 of the present embodiment. FIG. 4 is a flowchart showing the operation of the sound transfer characteristic measuring apparatus 2 of the present embodiment.

図３に示すように、本実施例の音響伝達特性測定装置２は、ＳＴＳＰ信号計算部２１と、スピーカアレイ１２と、逆ＳＴＳＰ信号計算部２３と、マイクロホンアレイ１４と、伝達特性推定値算出部１５を含む。   As shown in FIG. 3, the acoustic transfer characteristic measuring apparatus 2 of the present embodiment includes an STSP signal calculation unit 21, a speaker array 12, an inverse STSP signal calculation unit 23, a microphone array 14, and a transfer characteristic estimated value calculation unit. 15 is included.

実施例１のＳＴＳＰ信号計算部１１、逆ＳＴＳＰ信号計算部１３は、本実施例においてＳＴＳＰ信号計算部２１、逆ＳＴＳＰ信号計算部２３に変更されている。このほかの構成は、実施例１と実施例２において共通するため説明を割愛する。   The STSP signal calculator 11 and the inverse STSP signal calculator 13 of the first embodiment are changed to an STSP signal calculator 21 and an inverse STSP signal calculator 23 in the present embodiment. Since other configurations are common in the first embodiment and the second embodiment, description thereof is omitted.

＜ＳＴＳＰ信号計算部２１＞
本実施例のＳＴＳＰ信号計算部２１は、時空間周波数領域での計算を行わずに、時空間領域ＳＴＳＰ信号を直接計算してスピーカアレイから出力する（Ｓ２１）。具体的には、ＳＴＳＰ信号計算部２１は、以下の式（９）により、時空間領域ＳＴＳＰ信号を直接計算する。 <STSP signal calculator 21>
The STSP signal calculation unit 21 of this embodiment directly calculates the spatio-temporal region STSP signal without performing the calculation in the spatio-temporal frequency region and outputs it from the speaker array (S21). Specifically, the STSP signal calculation unit 21 directly calculates a spatio-temporal region STSP signal by the following equation (9).

Figure 0006087850
Figure 0006087850

＜逆ＳＴＳＰ信号計算部２３＞
本実施例の逆ＳＴＳＰ信号計算部２３は、時空間周波数領域での計算を行わずに、時空間領域ＳＴＳＰ信号に基づいて時空間領域逆ＳＴＳＰ信号を直接計算する（Ｓ２３）。より詳細には、逆ＳＴＳＰ信号計算部２３は、式（９）により計算した時空間領域ＳＴＳＰ信号f(xi,tj)を時空間方向に反転したものを時空間領域の逆ＳＴＳＰ信号とする。具体的には、逆ＳＴＳＰ信号計算部２３は、以下の式（１０）により、時空間領域逆ＳＴＳＰ信号を直接計算する。 <Inverse STSP signal calculator 23>
The inverse STSP signal calculation unit 23 of the present embodiment directly calculates the spatio-temporal domain inverse STSP signal based on the spatio-temporal domain STSP signal without performing the calculation in the spatio-temporal frequency domain (S23). More specifically, the inverse STSP signal calculation unit 23 converts the spatio-temporal domain STSP signal f (x i , t j ) calculated by Equation (9) into the spatio-temporal domain inverse STSP signal. To do. Specifically, the inverse STSP signal calculation unit 23 directly calculates a spatio-temporal domain inverse STSP signal by the following equation (10).

Figure 0006087850
Figure 0006087850

本実施例の音響伝達特性測定装置２によれば、時空間周波数領域の計算を行わずに、実施例１と同様の効果を得ることができる。 上述の各種の処理は、記載に従って時系列に実行されるのみならず、処理を実行する装置の処理能力あるいは必要に応じて並列的にあるいは個別に実行されてもよい。その他、本発明の趣旨を逸脱しない範囲で適宜変更が可能であることはいうまでもない。   According to the acoustic transfer characteristic measuring apparatus 2 of the present embodiment, the same effect as that of the first embodiment can be obtained without calculating the spatio-temporal frequency region. The various processes described above are not only executed in time series according to the description, but may also be executed in parallel or individually as required by the processing capability of the apparatus that executes the processes. Needless to say, other modifications are possible without departing from the spirit of the present invention.

また、上述の構成をコンピュータによって実現する場合、各装置が有すべき機能の処理内容はプログラムによって記述される。そして、このプログラムをコンピュータで実行することにより、上記処理機能がコンピュータ上で実現される。   Further, when the above-described configuration is realized by a computer, processing contents of functions that each device should have are described by a program. The processing functions are realized on the computer by executing the program on the computer.

この処理内容を記述したプログラムは、コンピュータで読み取り可能な記録媒体に記録しておくことができる。コンピュータで読み取り可能な記録媒体としては、例えば、磁気記録装置、光ディスク、光磁気記録媒体、半導体メモリ等どのようなものでもよい。   The program describing the processing contents can be recorded on a computer-readable recording medium. As the computer-readable recording medium, for example, any recording medium such as a magnetic recording device, an optical disk, a magneto-optical recording medium, and a semiconductor memory may be used.

また、このプログラムの流通は、例えば、そのプログラムを記録したＤＶＤ、ＣＤ−ＲＯＭ等の可搬型記録媒体を販売、譲渡、貸与等することによって行う。さらに、このプログラムをサーバコンピュータの記憶装置に格納しておき、ネットワークを介して、サーバコンピュータから他のコンピュータにそのプログラムを転送することにより、このプログラムを流通させる構成としてもよい。   The program is distributed by selling, transferring, or lending a portable recording medium such as a DVD or CD-ROM in which the program is recorded. Furthermore, the program may be distributed by storing the program in a storage device of the server computer and transferring the program from the server computer to another computer via a network.

このようなプログラムを実行するコンピュータは、例えば、まず、可搬型記録媒体に記録されたプログラムもしくはサーバコンピュータから転送されたプログラムを、一旦、自己の記憶装置に格納する。そして、処理の実行時、このコンピュータは、自己の記録媒体に格納されたプログラムを読み取り、読み取ったプログラムに従った処理を実行する。また、このプログラムの別の実行形態として、コンピュータが可搬型記録媒体から直接プログラムを読み取り、そのプログラムに従った処理を実行することとしてもよく、さらに、このコンピュータにサーバコンピュータからプログラムが転送されるたびに、逐次、受け取ったプログラムに従った処理を実行することとしてもよい。また、サーバコンピュータから、このコンピュータへのプログラムの転送は行わず、その実行指示と結果取得のみによって処理機能を実現する、いわゆるＡＳＰ（Application Service Provider）型のサービスによって、上述の処理を実行する構成としてもよい。なお、本形態におけるプログラムには、電子計算機による処理の用に供する情報であってプログラムに準ずるもの（コンピュータに対する直接の指令ではないがコンピュータの処理を規定する性質を有するデータ等）を含むものとする。   A computer that executes such a program first stores, for example, a program recorded on a portable recording medium or a program transferred from a server computer in its own storage device. When executing the process, the computer reads a program stored in its own recording medium and executes a process according to the read program. As another execution form of the program, the computer may directly read the program from a portable recording medium and execute processing according to the program, and the program is transferred from the server computer to the computer. Each time, the processing according to the received program may be executed sequentially. Also, the program is not transferred from the server computer to the computer, and the above-described processing is executed by a so-called ASP (Application Service Provider) type service that realizes the processing function only by the execution instruction and result acquisition. It is good. Note that the program in this embodiment includes information that is used for processing by an electronic computer and that conforms to the program (data that is not a direct command to the computer but has a property that defines the processing of the computer).

また、この形態では、コンピュータ上で所定のプログラムを実行させることにより、本装置を構成することとしたが、これらの処理内容の少なくとも一部をハードウェア的に実現することとしてもよい。   In this embodiment, the present apparatus is configured by executing a predetermined program on a computer. However, at least a part of these processing contents may be realized by hardware.

Claims (8)

複数のスピーカを直線状に配置してなるスピーカアレイと、
複数のマイクロホンを直線状に配置してなるマイクロホンアレイと、
合成時の平面波の角度と周波数とが時間方向に変化する信号であるＳＴＳＰ信号の時空間領域の表現である時空間領域ＳＴＳＰ信号を計算して前記スピーカアレイから出力するＳＴＳＰ信号計算部と、
前記時空間領域ＳＴＳＰ信号に基づいて時空間領域逆ＳＴＳＰ信号を計算する逆ＳＴＳＰ信号計算部と、
前記マイクロホンアレイに収音される信号である入力信号と、前記時空間領域逆ＳＴＳＰ信号に基づいて時空間音響伝達特性の推定値を算出する伝達特性推定値算出部と、
を含む音響伝達特性測定装置。 A speaker array in which a plurality of speakers are arranged in a straight line;
A microphone array in which a plurality of microphones are linearly arranged;
An STSP signal calculation unit that calculates a spatio-temporal region STSP signal that is a representation of the spatio-temporal region of the STSP signal that is a signal in which the angle and frequency of the plane wave at the time of synthesis change in the time direction, and outputs from the speaker array;
An inverse STSP signal calculator for calculating a spatiotemporal domain inverse STSP signal based on the spatiotemporal domain STSP signal;
An input signal that is a signal picked up by the microphone array, and a transfer characteristic estimated value calculation unit that calculates an estimated value of the spatiotemporal acoustic transfer characteristic based on the spatiotemporal domain inverse STSP signal;
An acoustic transfer characteristic measuring apparatus including: 請求項１に記載の音響伝達特性測定装置であって、
Δωを時間周波数の増加量を調整するパラメータとし、Δkxを空間周波数の増加量を調整するパラメータとし、ω0を時間周波数の初期値を調整するパラメータとし、kx0を空間周波数の初期値を調整するパラメータとし、i,jは空間方向（スピーカ方向）xiおよび時間方向tjのインデックスを示すものとし、Nxをスピーカ・マイクロホン数とし、Ntを時間方向の信号長とし、
前記ＳＴＳＰ信号計算部が、時空間領域ＳＴＳＰ信号を
Figure 0006087850
として計算し、
前記逆ＳＴＳＰ信号計算部が、時空間領域逆ＳＴＳＰ信号を
Figure 0006087850
として計算する
音響伝達特性測定装置。 The acoustic transfer characteristic measuring device according to claim 1,
Δω is a parameter for adjusting the increase amount of the temporal frequency, Δk x is a parameter for adjusting the increase amount of the spatial frequency, ω 0 is a parameter for adjusting the initial value of the temporal frequency, and k x0 is an initial value of the spatial frequency. As parameters to be adjusted, i and j indicate indexes in the spatial direction (speaker direction) x i and time direction t j , N x is the number of speakers and microphones, N t is the signal length in the time direction,
The STSP signal calculator calculates a spatio-temporal region STSP signal.
Figure 0006087850
Calculate as
The inverse STSP signal calculation unit calculates a space-time domain inverse STSP signal.
Figure 0006087850
Sound transfer characteristic measuring device that calculates as 複数のスピーカを直線状に配置してなるスピーカアレイと、
複数のマイクロホンを直線状に配置してなるマイクロホンアレイと、
合成時の平面波の角度と周波数とが時間方向に変化する信号であるＳＴＳＰ信号の時空間周波数領域の表現である時空間周波数領域ＳＴＳＰ信号を、時間周波数と空間周波数とが線形変化する条件で計算する時空間周波数領域ＳＴＳＰ信号計算部と、
前記時空間周波数領域ＳＴＳＰ信号に基づいて、前記ＳＴＳＰ信号の時空間領域の表現である時空間領域ＳＴＳＰ信号を計算して前記スピーカアレイから出力する時空間領域ＳＴＳＰ信号計算部と、
前記時空間周波数領域ＳＴＳＰ信号に畳み込むことでパルス信号を生成する信号である時空間周波数領域逆ＳＴＳＰ信号を計算する時空間周波数領域逆ＳＴＳＰ信号計算部と、
前記時空間周波数領域逆ＳＴＳＰ信号に基づいて時空間領域逆ＳＴＳＰ信号を計算する時空間領域逆ＳＴＳＰ信号計算部と、
前記マイクロホンアレイに収音される信号である入力信号と、前記時空間領域逆ＳＴＳＰ信号に基づいて時空間音響伝達特性の推定値を算出する伝達特性推定値算出部と、
を含む音響伝達特性測定装置。 A speaker array in which a plurality of speakers are arranged in a straight line;
A microphone array in which a plurality of microphones are linearly arranged;
A spatio-temporal frequency domain STSP signal, which is a representation of the spatio-temporal frequency domain of the STSP signal, which is a signal in which the angle and frequency of the plane wave at the time of synthesis change in the time direction, is calculated under the condition that the temporal frequency and the spatial frequency change linearly. A spatio-temporal frequency domain STSP signal calculator,
A spatio-temporal region STSP signal calculator that calculates a spatio-temporal region STSP signal that is a representation of the spatio-temporal region of the STSP signal based on the spatio-temporal frequency region STSP signal, and outputs the spatio-temporal region STSP signal from the speaker array;
A spatio-temporal frequency domain inverse STSP signal calculation unit for calculating a spatio-temporal frequency domain inverse STSP signal, which is a signal for generating a pulse signal by convolution with the spatio-temporal frequency domain STSP signal;
A spatio-temporal domain inverse STSP signal calculator that calculates a spatio-temporal domain inverse STSP signal based on the spatio-temporal frequency domain inverse STSP signal;
An input signal that is a signal picked up by the microphone array, and a transfer characteristic estimated value calculation unit that calculates an estimated value of the spatiotemporal acoustic transfer characteristic based on the spatiotemporal domain inverse STSP signal;
An acoustic transfer characteristic measuring apparatus including: 請求項３に記載の音響伝達特性測定装置であって、
前記時空間周波数領域ＳＴＳＰ信号計算部が、
kx,nを空間周波数とし、ωmを時間周波数とし、nを空間周波数kx,nのインデクスとし、mを時間周波数ωmのインデクスとし、Δωを時間周波数の増加量を調整するパラメータとし、Δkxを空間周波数の増加量を調整するパラメータとし、ω0を時間周波数の初期値を調整するパラメータとし、kx0を空間周波数の初期値を調整するパラメータとし、前記時空間周波数領域ＳＴＳＰ信号を、
Figure 0006087850
として計算する
音響伝達特性測定装置。 The acoustic transfer characteristic measuring device according to claim 3,
The spatio-temporal frequency domain STSP signal calculator is
k x, n is the spatial frequency, ω m is the time frequency , n is the index of the spatial frequency k x, n , m is the index of the time frequency ω m , and Δω is a parameter that adjusts the amount of increase of the time frequency. , Δk x is a parameter for adjusting the increase amount of the spatial frequency, ω 0 is a parameter for adjusting the initial value of the temporal frequency, k x0 is a parameter for adjusting the initial value of the spatial frequency, and the spatio-temporal frequency domain STSP signal The
Figure 0006087850
Sound transfer characteristic measuring device that calculates as 請求項１から４の何れかに記載の音響伝達特性測定装置であって、
前記伝達特性推定値算出部が、
２次元円状畳み込みを用いる方法、２次元線状畳み込みを用いる方法、円状畳み込みと線状畳み込みを組み合わせる方法、のいずれかにより時空間音響伝達特性の推定値を算出する
音響伝達特性測定装置。 The acoustic transfer characteristic measuring device according to any one of claims 1 to 4,
The transfer characteristic estimated value calculation unit
An acoustic transfer characteristic measuring apparatus that calculates an estimated value of a spatio-temporal acoustic transfer characteristic by any one of a method using a two-dimensional circular convolution, a method using a two-dimensional linear convolution, and a method combining a circular convolution and a linear convolution. 複数のスピーカを直線状に配置してなるスピーカアレイと、複数のマイクロホンを直線状に配置してなるマイクロホンアレイとを用いる音響伝達特性測定方法であって、
合成時の平面波の角度と周波数とが時間方向に変化する信号であるＳＴＳＰ信号の時空間領域の表現である時空間領域ＳＴＳＰ信号を計算して前記スピーカアレイから出力するＳＴＳＰ信号計算ステップと、
前記時空間領域ＳＴＳＰ信号に基づいて時空間領域逆ＳＴＳＰ信号を計算する逆ＳＴＳＰ信号計算ステップと、
前記マイクロホンアレイに収音される信号である入力信号と、前記時空間領域逆ＳＴＳＰ信号に基づいて時空間音響伝達特性の推定値を算出する伝達特性推定値算出ステップと、
を含む音響伝達特性測定方法。 An acoustic transfer characteristic measurement method using a speaker array in which a plurality of speakers are linearly arranged and a microphone array in which a plurality of microphones are linearly arranged,
An STSP signal calculation step of calculating a spatio-temporal region STSP signal that is a representation of the spatio-temporal region of the STSP signal that is a signal in which the angle and frequency of the plane wave at the time of synthesis change in the time direction, and outputting from the speaker array;
An inverse STSP signal calculation step of calculating a spatiotemporal domain inverse STSP signal based on the spatiotemporal domain STSP signal;
A transfer characteristic estimated value calculating step for calculating an estimated value of a spatiotemporal acoustic transfer characteristic based on the input signal that is a signal collected by the microphone array and the spatiotemporal domain inverse STSP signal;
A method for measuring acoustic transfer characteristics including: 複数のスピーカを直線状に配置してなるスピーカアレイと、複数のマイクロホンを直線状に配置してなるマイクロホンアレイとを用いる音響伝達特性測定方法であって、
合成時の平面波の角度と周波数とが時間方向に変化する信号であるＳＴＳＰ信号の時空間周波数領域の表現である時空間周波数領域ＳＴＳＰ信号を、時間周波数と空間周波数とが線形変化する条件で計算する時空間周波数領域ＳＴＳＰ信号計算ステップと、
前記時空間周波数領域ＳＴＳＰ信号に基づいて、前記ＳＴＳＰ信号の時空間領域の表現である時空間領域ＳＴＳＰ信号を計算して前記スピーカアレイから出力する時空間領域ＳＴＳＰ信号計算ステップと、
前記時空間周波数領域ＳＴＳＰ信号に畳み込むことでパルス信号を生成する信号である時空間周波数領域逆ＳＴＳＰ信号を計算する時空間周波数領域逆ＳＴＳＰ信号計算ステップと、
前記時空間周波数領域逆ＳＴＳＰ信号に基づいて時空間領域逆ＳＴＳＰ信号を計算する時空間領域逆ＳＴＳＰ信号計算ステップと、
前記マイクロホンアレイに収音される信号である入力信号と、前記時空間領域逆ＳＴＳＰ信号に基づいて時空間音響伝達特性の推定値を算出する伝達特性推定値算出ステップと、
を含む音響伝達特性測定方法。 An acoustic transfer characteristic measurement method using a speaker array in which a plurality of speakers are linearly arranged and a microphone array in which a plurality of microphones are linearly arranged,
A spatio-temporal frequency domain STSP signal, which is a representation of the spatio-temporal frequency domain of the STSP signal, which is a signal in which the angle and frequency of the plane wave at the time of synthesis change in the time direction, is calculated under the condition that the temporal frequency and the spatial frequency change linearly. A spatio-temporal frequency domain STSP signal calculation step;
A spatio-temporal region STSP signal calculation step of calculating a spatio-temporal region STSP signal, which is a representation of the spatio-temporal region of the STSP signal, based on the spatio-temporal frequency region STSP signal and outputting from the speaker array;
A spatio-temporal frequency domain inverse STSP signal calculation step of calculating a spatio-temporal frequency domain inverse STSP signal that is a signal that generates a pulse signal by convolution with the spatio-temporal frequency domain STSP signal;
A spatio-temporal domain inverse STSP signal calculation step of calculating a spatio-temporal domain inverse STSP signal based on the spatio-temporal frequency domain inverse STSP signal;
A transfer characteristic estimated value calculating step for calculating an estimated value of a spatiotemporal acoustic transfer characteristic based on the input signal that is a signal collected by the microphone array and the spatiotemporal domain inverse STSP signal;
A method for measuring acoustic transfer characteristics including: コンピュータを、請求項１から５の何れかに記載の音響伝達特性測定装置として機能させるためのプログラム。   A program for causing a computer to function as the acoustic transfer characteristic measuring device according to any one of claims 1 to 5.
JP2014014046A 2014-01-29 2014-01-29 Acoustic transfer characteristic measuring device, acoustic transfer characteristic measuring method, program Active JP6087850B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2014014046A JP6087850B2 (en) 2014-01-29 2014-01-29 Acoustic transfer characteristic measuring device, acoustic transfer characteristic measuring method, program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2014014046A JP6087850B2 (en) 2014-01-29 2014-01-29 Acoustic transfer characteristic measuring device, acoustic transfer characteristic measuring method, program

Publications (2)

Publication Number Publication Date
JP2015142248A JP2015142248A (en) 2015-08-03
JP6087850B2 true JP6087850B2 (en) 2017-03-01

Family

ID=53772342

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2014014046A Active JP6087850B2 (en) 2014-01-29 2014-01-29 Acoustic transfer characteristic measuring device, acoustic transfer characteristic measuring method, program

Country Status (1)

Country Link
JP (1) JP6087850B2 (en)

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE10254470B4 (en) * 2002-11-21 2006-01-26 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Apparatus and method for determining an impulse response and apparatus and method for presenting an audio piece
JP2005341267A (en) * 2004-05-27 2005-12-08 Victor Co Of Japan Ltd Sound field correcting device

Also Published As

Publication number Publication date
JP2015142248A (en) 2015-08-03

Similar Documents

Publication Publication Date Title
EP3232686B1 (en) Neural network-based loudspeaker modeling with a deconvolution filter
JP6458738B2 (en) Sound field reproduction apparatus and method, and program
KR102009274B1 (en) Fir filter coefficient calculation for beam forming filters
US8787501B2 (en) Distributed sensing of signals linked by sparse filtering
JP6485711B2 (en) Sound field reproduction apparatus and method, and program
JPWO2006011356A1 (en) Impulse response measuring method and apparatus
CN108293170A (en) Adaptive no phase distortion amplitude response in beam forming application is balanced
JP2019078864A (en) Musical sound emphasis device, convolution auto encoder learning device, musical sound emphasis method, and program
JP2015079080A (en) Sound source position estimation device, method, and program
Dietrich Uncertainties in acoustical transfer functions: modeling, measurement and derivation of parameters for airborne and structure-borne sound
Antonello et al. Source localization and signal reconstruction in a reverberant field using the FDTD method
JP5986966B2 (en) Sound field recording / reproducing apparatus, method, and program
JP2018077139A (en) Sound field estimation device, sound field estimation method and program
JP6087850B2 (en) Acoustic transfer characteristic measuring device, acoustic transfer characteristic measuring method, program
JP6087856B2 (en) Sound field recording and reproducing apparatus, system, method and program
Katzberg et al. Spatial interpolation of room impulse responses using compressed sensing
JP2017118189A (en) Sound collection signal estimating device, sound collection signal estimating method and program
JP5583181B2 (en) Cascade connection type transmission system parameter estimation method, cascade connection type transmission system parameter estimation device, program
JP5713964B2 (en) Sound field recording / reproducing apparatus, method, and program
JP2018120129A (en) Sound field estimation device, method and program
JP2018191255A (en) Sound collecting device, method thereof, and program
Liang et al. Sound field reconstruction using neural processes with dynamic kernels
JP7375905B2 (en) Filter coefficient optimization device, filter coefficient optimization method, program
JP6044043B2 (en) Plane wave expansion method, apparatus and program for sound field
JP5325134B2 (en) Echo canceling method, echo canceling apparatus, program thereof, and recording medium

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20160302

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20170120

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20170202

R150 Certificate of patent or registration of utility model

Ref document number: 6087850

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150