JP5106936B2 - Sound field reproduction filter calculation device and sound field reproduction system - Google Patents

Description

  The present invention relates to a sound field reproduction filter calculating apparatus and a sound field reproduction system, and more specifically, based on the acoustic characteristics measured in the original sound field and the acoustic characteristics measured in the reproduced sound field, the original sound in the reproduced sound field. The present invention relates to a sound field reproduction filter calculation device that calculates a sound field reproduction filter for reproducing a sound field environment of a field, and a sound field reproduction system including the device.

  Conventionally, after recording (measuring) the sound field characteristics of a certain space (original sound field), a sound field that virtually reproduces the sound field environment of the measured original sound field in a different space (reproduced sound field). A reproduction method has been proposed (see, for example, Patent Document 1).

  In the sound field reproduction method described above, first, a sound field reproduction filter for performing sound field reproduction is calculated based on the acoustic characteristics of the original sound field and the acoustic characteristics of the reproduced sound field. Next, by performing acoustic processing (for example, convolution calculation processing using a FIR (Finite Impulse Response) filter) on the acoustic signal output in the reproduced sound field using the calculated sound field reproduction filter, The sound processing of the reproduced sound field is performed on the output sound output to the reproduced sound field. By outputting the acoustic signal subjected to the acoustic processing in the reproduced sound field in this way, the listener can obtain a feeling as if listening in the original sound field in the reproduced sound field.

  By the way, as described above, in order to reproduce the acoustic characteristics of the original sound field in the reproduced sound field, it is necessary to cancel the acoustic characteristics in the reproduced sound field and to effectively apply the acoustic characteristics of the original sound field in the reproduced sound field. . That is, the sound field reproduction filter has a role of canceling the acoustic characteristics of the reproduced sound field and making the acoustic characteristics in the original sound field obvious.

FIG. 7 shows the relationship between the above-described sound field reproduction filter, the acoustic characteristics of the original sound field, and the acoustic characteristics of the reproduced sound field. To reproduce the acoustic characteristics of the original sound field in the reproduced sound field, where D (Z) is the acoustic characteristic of the original sound field, A (Z) is the acoustic characteristic of the reproduced sound field, and Q (Z) is the sound field reproduction filter. ,
D (Z) = A (Z) Q (Z) Formula 1
It is necessary to obtain a sound field reproduction filter Q (Z) that satisfies the above relationship.

と表現することができる。 For example, as shown in FIG. 8A, the number of speakers (sound sources) c in the space A1 that is the target of the original sound field is C, and the control point m (output sound output from the sound source in the space A1 is collected. When the total number of points) is M, in the space A1, each speaker from each speaker c (c = 1, 2,..., C) to each control point m (m = 1, 2,..., M). The transfer function is in the form of Z-transform

It can be expressed as

と表現することができる。 Further, as shown in FIG. 8 (b), when the number of speakers s in the space A2 to be reproduced sound field is S and the total number of control points m is M, the speakers s (s = 1, 1) in the space A2. 2,..., S) to the control point m (m = 1, 2,..., M)

It can be expressed as

  Note that the transfer function of the reproduced sound field and the transfer function of the original sound field described above satisfy the causality (0 when t <0).

If the sound field reproduction system is a system that can reproduce the acoustic environment of the original sound field in the reproduced sound field, and the sound field reproduction system reproduces the acoustic characteristics of the original sound field in the reproduced sound field, the sound field reproduction filter For example, it is necessary to perform filter processing (acoustic processing) using an FIR filter. The sound field reproduction system used for this filtering process sets 0 ≦ t ≦ N _D −1 as the target of the time range of sound field reproduction as a finite time range of time response. Furthermore, since it is desirable to use an environment with few unnecessary reflections and a short reverberation time as a reproduction environment of the reproduced sound field, the time response of the space A2 can be handled with the reverberation time N cut off.

と表すことができ、
また、再現音場の伝達関数は

と表すことができる。 Under these conditions, the transfer function of the original sound field is

Can be expressed as
The transfer function of the reproduced sound field is

It can be expressed as.

で表すことができる。 When the sound reproduction filter to obtain a Q _{z (Z),} the relationship between D _z (Z) and A _z (Z) and Q _{z (Z),}
D _z (Z) = A _z (Z) Q _z (Z) Expression 6
A sound field reproduction filter (FIR filter) Q _z (Z) having a finite number of stages L satisfying this equation 6 can be expressed by:

Can be expressed as

である。 Where q _z (Z) is

It is.

  The stage number L is a constant given as a design parameter, and an appropriate value is set in consideration of the length of the transfer function to be reproduced, the arithmetic processing capability of the sound field reproduction system, and the like. The reverberation time of the acoustic signal to which the sound field reproduction filter is applied is determined according to the value of the stage number L.

As shown in FIG. 7, the sound field reproduction filter Q _z (Z) described above is subjected to acoustic processing using the sound field reproduction filter at the input stage of the speaker s to the acoustic signal, so that the transfer function of the original sound field is obtained. It is possible to reproduce in the reproduction sound field.
The above equation 6 expressed by Z conversion
D _z (Z) = A _z (Z) Q _z (Z) Expression 6
Is an equivalent linear algebraic equation D = AQ Equation 9
Can be expressed as

である。 Where D is

It is.

と表される。
ここで、０≦ｔ≦Ｎ_D−１の範囲外の時刻において、ｄt_i,j（t）の値を０とする。 d _{i, j} is a column vector in which N + L−1 sample values d t _{i, j} (t) at the time t of the time response of the transfer function dz _{i, j} (Z) from the sound source j to the control point i are arranged. Yes,

It is expressed.
Here, the value of dt _{i, j} (t) is set to 0 at a time outside the range of 0 ≦ t ≦ N _D −1.

と表すことができる。 Similarly for Q,

It can be expressed as.

で示される。 Here, q _{i, j} are vectors representing L-stage FIR filter coefficients, respectively.

Indicated by

A can be expressed as the following equation.

と表すことができ、この畳み込み行列を使用することによって、畳み込み演算処理を行列の積として表現することが可能となる。 Here, A _{i, j} which is a submatrix of A is a convolution matrix of N + L−1 rows and L columns shown in the following Expression 15, and is configured by a regular pattern in which elements in the diagonal direction have the same value.

By using this convolution matrix, the convolution operation processing can be expressed as a matrix product.

  As described above, by using the matrix D, the matrix A, and the matrix Q, Q is obtained as a least squares solution of D = AQ. The least squares solution is

^{Q = (A t A + lE} LS × LS) -1 A t D ··· formula 16
Can be obtained by solving
Here, A ^t represents the transposed matrix of A, E _{LS × LS} indicates a unit matrix of the LS line LS column.

Q obtained by Equation 16 is an optimal solution under the condition that the number of filter stages is L.
JP 2007-411164 A

  By the way, when the sound field reproduction filter Q is calculated by solving the above equation 16 in the time domain, the sound field reproduction system becomes large (when the number of speakers or the number of microphones is large). As the number of stages L of the filter (FIR filter) increases, there is a problem that the calculation amount (calculation scale) of the calculation increases. For this reason, the method of calculating the sound field reproduction filter in the time domain has not been a realistic method because of the burden of calculation amount and the like.

For example, Equation 16 above
^{Q = (A t A + lE} LS × LS) -1 A t D
Is calculated by numerical calculation,
R = (A ^t A + lE _{LS × LS} )
P = A ^t D
Can be reduced to the problem of solving the simultaneous linear equation P = RQ, that is, the problem of obtaining Q = R ⁻¹ P.

Consider the case where the simultaneous linear equations are solved by a direct solution using a computer or the like. As an example of the method of solving by the direct solution method, for example, Gaussian elimination method or the like corresponds to the example. When n is the number of unknowns obtained using the direct solution method, the memory capacity required in the calculation process for solving the above equation is on the order of n ² , and the number of operations is on the order of n ³ .

  When obtaining a sound field reproduction filter, it is necessary to solve the problem of n = L × S. For example, the filter stage number L is 4096 (L = 4096), and the number S of speakers s in the sound field reproduction system is 64 (S = 64), and the data accuracy of the filter calculation is a single precision real number (4 bytes), a computer having a memory capacity of about α × 256 Gbytes is required. There was a problem that it was difficult to do.

On the other hand, as a method for efficiently calculating the sound field reproduction filter,
D _z (Z) = A _z (Z) Q _z (Z) Expression 6
In general, a method is used in which a time domain filter coefficient is calculated by inverse transformation after Fourier transform of both sides is obtained in the frequency domain. In this method, the problem of convolution calculation processing can be replaced with the problem of multiplication for each frequency component, and the filter coefficient is calculated by repeatedly calculating a small matrix having a matrix size of S × S by the number of frequency components. Can do. For this reason, even when the scale of the reproduction sound field system becomes large, there is an advantage that calculation calculation of the sound field reproduction filter becomes easy.

Specifically, D _ω (ω) = A _ω (ω) Q _ω (ω) Expression 17 by Fourier transform of Expression 6
Can be expressed as

を示しており、Ａ_ω（ω）は

を示している。 Where D _ω (ω) is

A _ω (ω) is

Is shown.

Using this D _ω (ω) and A _ω (ω), the formula for obtaining Q _ω (ω) is:
Q _ω (ω) = (A _ω ^* (ω) A _ω (ω) + l E _{S × S} ) ⁻¹ (A _ω ^* (ω) D _ω (ω))
... Formula 20
It becomes.

で示される。なおconj(x)は、複素数ｘの複素共役を表す。 Where A _ω ^* (ω) represents the complex conjugate transpose matrix of A _ω (ω),

Indicated by Note that conj (x) represents a complex conjugate of the complex number x.

  However, the sound field reproduction filter calculated in the frequency domain does not agree with the optimal solution obtained in the time domain, and the accuracy of the time response tends to be insufficient. For this reason, there is a possibility that sufficient reproduction accuracy cannot be obtained when the acoustic characteristics of the original sound field are reproduced in the reproduced sound field by applying the sound field reproduction filter obtained in the frequency domain.

  The present invention has been made in view of the above problems, and by obtaining a sound field reproduction filter in the time domain, it is possible to avoid occurrence of a calculation error like the sound field reproduction filter obtained in the frequency domain, and in the time domain. It is an object of the present invention to provide a sound field reproduction filter calculation device and a sound field reproduction system that can reduce the calculation burden of the sound field reproduction filter.

  In order to solve the above problems, the sound field reproduction filter calculation device according to the present invention is based on the acoustic characteristics measured in the original sound field and the acoustic characteristics measured in the reproduction sound field, and The present invention relates to a sound field reproduction filter calculation device that calculates a sound field reproduction filter for reproducing a sound field environment, and the sound field reproduction filter calculation device calculates a sound field reproduction filter by applying an iterative method in the time domain. And performing only the correlation calculation process and the convolution calculation process that occur in the application of the iterative solution method by performing the calculation process by the conversion process to the frequency domain, and then using the calculation method that reversely converts the calculation result to the time domain. And

  In the sound field reproduction filter calculation device according to the present invention, the calculation process of the sound field reproduction filter is executed by applying an iterative solution in the time domain. For this reason, the calculated value does not contain an error as in the sound field reproduction filter calculated by converting the calculation result to the time domain after performing the calculation process after converting to the frequency domain. Therefore, it is possible to calculate a highly accurate sound field reproduction filter.

  In addition, when calculating the sound field reproduction filter using the iterative solution method, the solution in the time domain is obtained using the calculation result in the frequency domain only when the correlation calculation process and the convolution calculation process are performed. It is possible to reduce the memory load and the memory capacity used for arithmetic processing. For this reason, it becomes possible to aim at the speeding up of the whole process, and the efficiency of the calculation assets utilized for calculation processing.

  Furthermore, in the sound field reproduction filter calculation device described above, the correlation obtained by the sound field reproduction filter calculation device performing an arithmetic process by a conversion process to the frequency domain and then inversely converting the calculation result to the time domain. Of the calculation results of the calculation processing or the convolution calculation processing, only data in a range in which an error from the calculation result obtained when the correlation calculation processing or the convolution calculation processing is performed only in the time domain cannot be used. It is preferable that the calculation process of the sound field filter in the time domain is executed.

  In this way, there can be no error from the calculation result obtained when the calculation process is performed only in the time domain among the calculation results of the correlation calculation process or the convolution calculation process obtained by inversely transforming from the frequency domain to the time domain. Since the sound field filter calculation process is executed using only the range data, it is possible to effectively prevent errors and the like from being included in the finally calculated sound field reproduction filter value. Become.

  Further, the range in which the error does not occur may be determined according to the number of stages of the sound field reproduction filter set in advance according to the reverberation time of the output sound to which the sound field reproduction filter is applied. Good.

  Thus, when the range in which no error occurs is determined according to the number of stages of the sound field reproduction filter set in advance according to the reverberation time or the like, the performance of the calculated sound field reproduction filter is sufficiently high. While ensuring, it becomes possible to suppress an error etc. being included in the value of the sound field reproduction filter finally calculated.

  The sound field reproduction system according to the present invention includes a speaker installed in the reproduction sound field, a sound reproduction device that reproduces an acoustic signal to be output from the speaker, the above-described sound field reproduction filter calculation device, A sound processing device that performs sound processing on the sound signal calculated by the sound reproduction device using the sound field reproduction filter calculated by the sound field reproduction filter calculation device, and the sound processing device performs sound processing. And an acoustic output device for outputting an acoustic signal from the speaker.

  As described above, the sound field reproduction system includes the above-described sound field reproduction filter calculation device, so that the sound field calculated by performing the calculation process after converting to the frequency domain and inversely converting the calculation result to the time domain. Since the calculated value does not contain an error like the reproduction filter, it is possible to faithfully reproduce the acoustic environment of the original sound field in the reproduction sound field using a highly accurate sound field reproduction filter. .

  According to the sound field reproduction filter calculation device and the sound field reproduction system according to the present invention, it is possible to calculate a sound field reproduction filter with high accuracy, and further, to speed up and calculate the sound field reproduction filter. It is possible to improve the efficiency of the computing assets used.

  Hereinafter, a sound field reproduction filter calculation device and a sound field reproduction system according to the present invention will be described in detail with reference to the drawings. In this embodiment, a boundary sound field control method is used as the sound field reproduction method. The boundary sound field control method is that when listening to the output sound in the reproduced sound field, a plurality of control points are installed around the listener's head or in the vicinity of the auricle (earlobe). This is a method of perceiving the original sound field audibly by controlling the reproduced sound field so that the field characteristics are reproduced.

  First, a method for measuring the acoustic characteristics of the original sound field by the original sound field acoustic system will be described.

[Measurement of acoustic characteristics in the original sound field]
As shown in FIG. 1, an original sound field acoustic system 2 includes a speaker unit 3 installed in an acoustic space (original sound field) 6 desired to be reproduced, a sound output device 4 that outputs sound to the speaker unit 3, and a sound A computer 5 that outputs an evaluation sound to the output device 4, a microphone array 7 installed in the original sound field 6, and a voice input device 8 that inputs the recorded sound collected by the microphone array 7. ing.

  As shown in FIG. 1, the speaker unit 3 includes four speakers that are equally installed at the four corners of the space B1 of the original sound field 6, and each speaker is a speaker FL that is installed at the left front position of the space B1. And a speaker FR installed at the front right position, a speaker RL installed at the rear left position, and a speaker RR installed at the rear right position.

  The audio output device 4 has a role of outputting the evaluation sound output from the computer 5 to each speaker FR, FL, RR, RL of the speaker unit 3. Note that the audio output device 4 can individually output the output destination of the evaluation sound to the speakers FR, FL, RR, RL according to the control command of the computer 5.

  Here, the evaluation sound is an output sound used for calculating the acoustic characteristics of the original sound field 6 based on the recorded sound collected via the microphone array 7 and the sound input device 8, and specifically, A sound suitable for measuring the collected recorded sound as an impulse response (impulse input sound) is applicable.

  The microphone array 7 is composed of a group of microphones in which a plurality of (for example, M) microphones 12 are arranged in a circle. The microphone array 7 is installed at the center position of the space B1 that is equidistant from each speaker FR, FL, RR, RL.

  The voice input device 8 has a role of outputting the output sound of the speaker unit 3 measured by the microphone array 7 to the computer 5. Here, since a plurality of microphones 12 are installed in the microphone array 7, each speaker FR, FL, RR, RL and each microphone are measured by measuring an acoustic signal of the microphone array 7 input in the voice input device 8. It is possible to measure the acoustic characteristics between 12.

  The computer 5 has a function of generating the above-described evaluation sound and calculating the acoustic characteristic of the original sound field 6 based on the recorded sound collected by the microphone array 7 in order to measure the acoustic characteristic of the space B1. ing.

  The computer 5 is roughly configured by an evaluation sound generation unit 15, an impulse response calculation unit 16, an impulse response recording unit 17, and a control unit 18. The evaluation sound generator 15 has a function of generating an evaluation sound for collecting the impulse response described above. The evaluation sound generation unit 15 generates an evaluation sound in accordance with a control command from the control unit 18, and outputs the evaluation sound to the audio output device 4 via the control unit 18.

  The impulse response calculation unit 16 calculates an impulse response to the evaluation sound based on the recorded sound collected via the microphone array 7 and the voice input device 8.

  The impulse response recording unit 17 records the calculated impulse response as acoustic characteristics. The impulse response recording unit 17 may be a general hard disk or the like, or may be composed of a RAM (Random Access Memory), a flash memory, or the like.

  The control unit 18 has a role of performing processing necessary for obtaining the acoustic characteristics in the original sound field 6. A CPU (Central Processing Unit) that mainly performs arithmetic control processing, and measures the acoustic characteristics in the original sound field 6. These are configured by a general arithmetic processor such as a ROM (Read Only Memory) in which the above program is recorded, a RAM (Random Access Memory) used in a work area, and the like.

  As described above, the control unit 18 causes the evaluation sound generation unit 15 to generate an evaluation sound and has a role of outputting the generated evaluation sound to the sound output device 4. In addition, the control unit 18 causes the impulse response calculation unit 16 to calculate an impulse response based on the recorded sound input from the voice input device 8 and records the calculated impulse response as an acoustic characteristic in the original sound field 6 as an impulse response recording. The unit 17 has a role of recording.

  In the original sound field acoustic system 2 configured as described above, the control unit 18 causes the evaluation sound generation unit 15 to generate an evaluation sound, and the generated evaluation sound is output via the audio output device 4 to the speaker unit. 3 from each speaker FR, RL, RR, RL. When outputting the evaluation sound from the speaker, the audio output device 4 outputs the evaluation sound individually for each speaker in accordance with the control command of the control unit 18.

  In the microphone array 7, output sounds output from any of the speakers are picked up by a plurality (M) of microphones 12 and output to the control unit 18 via the voice input device 8. The control unit 18 causes the impulse response calculation unit 16 to perform impulse response calculation processing based on the recorded sound acquired from the voice input device 8, and causes the impulse response recording unit 17 to record the calculated impulse response. The control unit 18 performs the impulse response calculation process for each speaker, causes the impulse response recording unit 17 to record the acoustic characteristics corresponding to each speaker, and ends the acoustic characteristic measurement process in the original sound field 6. .

  The acoustic characteristics measured in the original sound field acoustic system 2 are a number of impulse responses between the four speakers FR, RL, RR, RL of the speaker unit 3 and the M microphones 12 of the microphone array 7. The acoustic characteristics obtained based on all the impulse responses are recorded in the impulse response recording unit 17.

[Measurement of acoustic characteristics in reproduced sound field]
Next, a reproduced sound field acoustic system for measuring acoustic characteristics in a reproduced sound field will be described. In the present embodiment, the laboratory space B2 is used as an example of a reproduced sound field. By using the laboratory as a reproduction sound field in this way, the acoustic environment of the original sound field 6 described above can be faithfully reproduced in the laboratory.

  As shown in FIG. 2, the reproduced sound field acoustic system 20 includes a speaker unit 22 installed in a laboratory space B <b> 2 that is a reproduced sound field 21, and an audio output device 4 that outputs sound to the speaker unit 22. The computer 5 that outputs the evaluation sound to the sound output device 4, the microphone array 7 installed in the original sound field 6, and the sound input device 8 that inputs the recorded sound collected by the microphone array 7. It is configured.

  Here, the audio output device 4, the computer 5, the microphone array 7, and the audio input device 8 used in the reproduced sound field acoustic system 20 are the audio output device 4 used in the original sound field acoustic system 2 described above, Since the computer 5, the microphone array 7, and the voice input device 8 have the same configuration and the same functions, detailed descriptions thereof are omitted.

  The evaluation sound generator 15, the impulse response calculator 16, the impulse response recording unit 17, and the controller 18 that roughly configure the computer 5 are evaluated in accordance with the computer 5 of the original sound field acoustic system 2 described above. Since the sound generating unit 15, the impulse response calculating unit 16, the impulse response recording unit 17, and the control unit 18 have the same configuration and the same functions, detailed description thereof is omitted.

  As shown in FIG. 2, the speaker unit 22 includes a plurality of (for example, M) speakers 23 arranged in a circle so as to maintain an equal distance from the microphone array 7. The speaker 23 can output the evaluation sound output from the audio output device 4.

  On the other hand, the audio output device 4 can individually output the output destination of the evaluation sound to each speaker 23 according to the control command of the computer 5.

  In the reproduced sound field acoustic system 20 configured as described above, when the control unit 18 performs the acoustic characteristic measurement process of the reproduced sound field 21, the control unit 18 causes the evaluation sound generation unit 15 to generate an evaluation sound, The generated evaluation sound is output from each speaker 23 of the speaker unit 22 via the audio output device 4. When outputting the evaluation sound from the speaker 23, the audio output device 4 outputs the evaluation sound individually for each speaker 23 according to the control command of the control unit 18.

  In the microphone array 7, output sounds output from any of the speakers 23 are collected by a plurality (M) of microphones 12 and output to the control unit 18 via the audio input device 8. The control unit 18 causes the impulse response calculation unit 16 to perform impulse response calculation processing based on the recorded sound acquired from the voice input device 8 and causes the impulse response recording unit 17 to record the calculated impulse response. The control unit 18 executes the impulse response calculation process for each speaker 23, records the acoustic characteristics corresponding to each speaker 23 in the impulse response recording unit 17, and performs the acoustic characteristic measurement process in the reproduced sound field 6 Exit.

  The acoustic characteristics measured in the reproduced sound field acoustic system 20 are obtained based on a number of impulse responses between the M speakers 23 of the speaker unit 22 and the M microphones 12 of the microphone array 7. The acoustic characteristics obtained based on all the impulse responses are recorded in the impulse response recording unit 17.

[Calculation of sound field reproduction filter]
Next, the sound for realizing the acoustic environment of the original sound field 6 in the reproduced sound field 21 based on the acoustic characteristics calculated in the original sound field acoustic system 2 and the acoustic characteristics calculated in the reproduced sound field acoustic system 20. Processing for calculating the field reproduction filter will be described.

  FIG. 3 is a block diagram showing a sound field reproduction filter calculation device (sound field reproduction filter calculation device) 30 for calculating a sound field reproduction filter.

  The sound field reproduction filter calculation device 30 is roughly configured by a filter calculation processing unit 31, an original sound field characteristic recording unit 32, a reproduction sound field characteristic recording unit 33, and a filter information recording unit 34.

  In the original sound field characteristic recording unit 32, acoustic characteristics measured in the original sound field 6 are recorded. Like the impulse response recording unit 17, a general hard disk, RAM (Random Access Memory), flash memory, or the like is used. It is constituted by a simple recording device. Similarly, the reproduced sound field characteristic recording unit 33 records acoustic characteristics measured in the reproduced sound field 21, and is configured by a general recording device such as a hard disk.

  The filter information recording unit 34 is a recording medium on which the sound field reproduction filter calculated by the filter calculation processing unit 31 is recorded. Similarly to the original sound field characteristic recording unit 32 and the reproduced sound field characteristic recording unit 33, the filter information recording unit 34 is also configured by a general recording device.

  The filter calculation processing unit 31 performs a process of calculating a sound field reproduction filter based on the acoustic characteristics recorded in the original sound field characteristic recording unit 32 and the acoustic characteristics recorded in the reproduction sound field characteristic recording unit 33. The filter calculation processing unit 31 includes a CPU (Central Processing Unit) mainly performing calculation processing of the sound field reproduction filter, a ROM (Read Only Memory) in which a program for calculating the sound field reproduction filter is recorded, a work area, and the like. It is constituted by a general arithmetic processor such as a RAM (Random Access Memory) used for the above.

As described above, the filter calculation processing unit 31
^{Q = (A t A + lE} LS × LS) -1 A t D ··· formula 16
R = (A ^t A + lE _{LS × LS} )
P = A ^t D
By solving the simultaneous linear equation P = RQ, the sound field reproduction filter Q (Q = R ⁻¹ P) in the time domain is obtained.

  Hereinafter, a process in which the filter calculation processing unit 31 obtains a sound field reproduction filter in the time domain will be described. The filter calculation processing unit 31 uses an iterative solution to obtain an optimal solution of Equation 16 in the time domain. In particular, in this embodiment, the filter calculation processing unit 31 performs an operation using a conjugate gradient method as a solution example of the iterative solution. To do.

The iterative solution method has a feature that the number of iterations of processing necessary for solution convergence is smaller than that of the direct solution method, and calculation resources (memory use amount and calculation amount) necessary for each iteration are reduced. When the number of unknowns is n and the matrix of equations for calculation is a sparse matrix (a matrix in which most components of the elements are 0), the solution is generally obtained by an iterative method. The memory usage required for the coefficient matrix and the number of calculations required for one iteration are on the order of n, and the calculation efficiency is increased. On the other hand, when the matrix of equations is a dense matrix, the order of the number of operations is on the order of n ² , and when n is large, the calculation is difficult. As shown in the present embodiment, when a sound field reproduction filter is obtained by an iterative solution, it is necessary to handle a large-scale dense matrix instead of a sparse matrix as a coefficient matrix.

In order to calculate the sound field reproduction filter Q, as described above, the equation 16
^{Q = (A t A + lE} LS × LS) -1 A t D ··· formula 16
Need to be solved. Here, when the sound field reproduction filter Q is calculated using an iterative solution method, the matrix form required in the iterative calculation is characterized by a regular convolution matrix. For this reason, in the sound field reproduction filter calculation device 30 according to the present invention, the memory amount necessary for the calculation process is reduced to an order proportional to the filter stage number L by performing a characteristic calculation process using this convolution matrix. It becomes possible.

  Specifically, the main matrix operation used in the process of assembling the equation and the process of the iterative calculation is calculated by a correlation operation process using fast Fourier transform (FFT) and a convolution operation process, thereby reducing the amount of calculation. It can be greatly reduced (specifically, it can be reduced to a level proportional to the calculation order by the fast Fourier transform (calculation order proportional to LlogL)). For this reason, even if the scale of the sound field reproduction system is increased, it is possible to perform arithmetic processing using a general computer.

  By using the arithmetic processing according to the present invention, specifically, the sound field when the number of filter stages L is 4096 (L = 4096) and the number S of speakers s in the sound field reproduction system is 64 (S = 64). The calculation of the reproduction filter can be executed by a computer having a memory of about 1 Gbyte.

FIG. 4 is a flowchart showing a process for calculating the sound field reproduction filter using the conjugate gradient method. When the calculation process is performed in the time domain according to the flowchart, there is a problem that the processing load increases in the calculation processes (1) to (4) described below.
(1) Matrix R: R = (A ^t A + L E _{LS × LS} ) arithmetic processing (2) Matrix P: P = A ^t D arithmetic processing (3) Rq ₀ product of R and vector in initial value calculation Arithmetic processing (4) Arithmetic processing of the product Rσ _μ of R and vector in the iterative calculation In addition, there is a problem that the storage capacity required for the calculation of the matrix R increases.

  As described above, the filter calculation processing unit 31 according to the present embodiment performs arithmetic processing in the frequency domain using fast Fourier transform (FFT) when performing the calculations (1) to (4). The processed calculation result is inversely converted into a time domain calculation result by inverse fast Fourier transform (IFFT) (returned to the time domain calculation result), thereby reducing the calculation processing load.

で示される。 First, calculation of the matrix R and the matrix P will be described as an initial calculation.
The matrix R is

Indicated by

で示される。 Here, each element of the matrix R is

Indicated by

で示される。 Each element of the matrix R _{i, j} is

Indicated by

である。 However,

It is.

The sub-matrix R _{i, j} of R is a correlation matrix, has a regular pattern, and has the characteristic that the values of the elements in the diagonal direction are the same value. The partial matrix R _{i, j} described above has 2L−1 different values. For this reason, when calculating the coefficient matrix R, the overlapping values are not stored in a storage area such as a memory, but only different values are stored as vectors, so that they are used in the calculation processing of the filter calculation processing unit 31. The amount of memory used can be reduced on the order of L.

において{}内は伝達関数同士の相関演算の形となっている。 Further, Expression 24 indicating the value of each element

In {}, the inside of {} is the form of correlation calculation between transfer functions.

For this reason, when performing the correlation calculation process, the fast Fourier transform (FFT) is used to convert the value of rt _{i, j} (τ) to the frequency domain value to obtain a solution. The inverse transform is performed by the transformation (IFFT), and the solution of the correlation calculation processing of rt _{i, j} (τ) in the time domain is obtained. By performing the arithmetic processing in the frequency domain in this way, the processing load can be reduced.

で表すと、両辺をフーリエ変換することによって、
Ｓ_ω（ω）＝conj（ｇ_ω（ω））ｈ_ω（ω） ・・・・式２６
となる。ここでconj（ｘ）はｘの複素共役を示している。 However, when performing correlation calculation processing using fast Fourier transform, caution is required in the following points.
A general expression for performing correlation calculation processing

In terms of the Fourier transform of both sides,
S _ω (ω) = conj (g _ω (ω)) h _ω (ω)...
It becomes. Here, conj (x) represents a complex conjugate of x.

As described above, the value of S _ω (ω) subjected to Fourier transform is obtained using frequency component multiplication, and then the obtained solution is subjected to inverse Fourier transform to convert it to time-series data in the time domain. A solution can be obtained quickly and easily as compared with the case where computation is performed in a region.

となる。ここで、Ｎ_FFTは、高速フーリエ変換を実施するときのデータ長であり、mod（ｔ＋τ，Ｎ_FFT）は、ｔ＋τをＮ_FFTで割ったときの余りを示している。 However, when performing the correlation calculation process using the fast Fourier transform, the calculation result S ′ _t (τ) obtained by obtaining the solution of S _ω (ω) using frequency domain multiplication and then returning it to the time domain by inverse Fourier transform. )

It becomes. Here, N _FFT is a data length when the fast Fourier transform is performed, and mod (t + τ, N _FFT ) indicates a remainder when t + τ is divided by N _FFT .

Thus, when S ′ _t (τ) is obtained by inverse Fourier transform, S ′ _t (τ) = S _t (τ) for all τ: Equation 28
However, if a solution that does not satisfy Expression 28 is used as it is in the calculation process of the conjugate gradient method shown in FIG. 4, there is a problem that an error occurs in the obtained sound field reproduction filter.

In order to avoid such a problem and obtain a solution satisfying the relational expression of Expression 28, the filter calculation processing unit 31 according to the present embodiment adds an appropriate number of 0 data to g and h, thereby increasing the speed. by increasing the data length N _FFT Fourier transform, S in the range of tau to be obtained 'value of _t (tau) is _{S' t (τ) = S} t (τ) ··· equation 28
Make adjustments so that

Therefore, even when rt _{i, j} (τ) described above is obtained, an error is not generated in the solution obtained by the calculation in the time domain so that the solution obtained by the inverse Fourier transform becomes an accurate value. ), The calculation is performed using only the data in the range of −L + 1 ≦ t ≦ L−1. That is, the range in which no error occurs is determined according to the number L of stages of the sound field reproduction filter set in advance according to the reverberation time of the output sound to which the sound field reproduction filter is applied.

で示される。 Similarly, each element of the matrix P can be calculated as a result of the correlation calculation process. The matrix P is

Indicated by

で示される。 Here, each element of the matrix P is

Indicated by

で示される。 Each element of the matrix P _{i, j} is

Indicated by

Since the inside of {} of this pt _{i, j} (τ) is a form of correlation calculation between transfer functions, the fast Fourier transform is applied to pt _{i, j} (τ) _, similarly to the matrix R, A value can be obtained by converting a value of pt _{i, j} (τ) to a frequency domain value to obtain a solution, and inverse Fourier transforming the obtained solution to inversely transform to a time domain solution. Thus, by using the process of inversely transforming the frequency domain solution into the time domain solution, it is possible to reduce the burden of the arithmetic processing.

Also when calculating pt _{i, j} (τ), the solution is obtained using only data in the range of 0 ≦ t ≦ L−1 so that the value obtained by the inverse Fourier transform becomes an accurate value.

  In this way, it is necessary to obtain the calculation results of the matrix R and the matrix P in advance and store them in a memory or the like, or to perform the calculation using the matrices R and P in the calculation process of the flowchart shown in FIG. The calculation process of the sound field reproduction filter is executed by calculating R and P each time.

  Next, specific processing of the filter calculation processing unit 31 will be described based on the flowchart shown in FIG.

  As shown in FIG. 4, the filter calculation processing unit 31 assigns 1 to j as an initial value and 0 to μ (step S.1), where j is the number of speakers in the original sound field 6 Is shown. Μ is a parameter indicating the number of operations (number of iterations) executed to calculate the value of the j-th column of the sound field reproduction filter Q. 3-step S.E. 5 shows the number of times the process of 5 has been repeatedly executed.

の初期値であるε₀の値
ε₀＝ｐ−Ｒｑ_{0 ・・・}式３３
を求める。なお、ε_μは、反復回数μでの暫定的な解をｑ_μとしたときの、誤差ベクトルを示したものである。 Next, the filter calculation unit 31 performs initial setting (step S.2). In the initial setting, the values in the jth column of P and Q are set as p and q, and ε _μ is set using an appropriate initial value q _0.

The value of ε ₀ , which is the initial value of ε ₀ = p−Rq _0, Equation 33
Ask for. Note that ε _μ represents an error vector when a provisional solution at the number of iterations _μ is q _μ .

で示される。 here. A calculation process (calculation of Rq ₀ ) between the matrix R and the vector used for calculating ε ₀ will be described.
The coefficient matrix R is

Indicated by

Here arbitrary vector x

とする。）
とし、行列Ｒと任意のベクトルｘとの積ｙを

(However, each element of x is x _{i, 1} is

And )
And the product y of the matrix R and an arbitrary vector x is

とする。）とすると、 (However, each element of y is y _{i, 1} is

And )

となり、 y is

And

となり、

として計算される。 y _{i, 1} which is a submatrix of y is

And

Is calculated as

Here, {} in yt _{i, 1} (τ) indicates a convolution operation. When performing convolution calculation processing, a solution is obtained in the frequency domain by converting the value of yt _{i, 1} (τ) into a frequency domain value using fast Fourier transform, as in the case of performing correlation calculation processing. Then, a value in the time domain can be obtained by performing an inverse Fourier transform (IFFT). Thus, by inversely transforming the solution obtained in the frequency domain into a time domain solution, it is possible to reduce the burden of arithmetic processing.

で表すと、両辺をフーリエ変換することによって、
Ｓ_ω（ω）＝ｇ_ω（ω）ｈ_ω（ω） ・・式４３
となる。 However, when performing the convolution calculation process using the fast Fourier transform, attention should be paid to the following points.
A general expression for performing correlation calculation processing

In terms of the Fourier transform of both sides,
S _ω (ω) = g _ω (ω) h _ω (ω)... 43
It becomes.

As described above, the value of S _ω (ω) subjected to Fourier transform is obtained using frequency component multiplication, and then the obtained solution is subjected to inverse Fourier transform to convert it to time-series data in the time domain. A solution can be obtained quickly and easily as compared with the case where computation is performed in a region.

となり、相関演算処理と同様に、求められたＳ’_t（τ）の全てのτについて
Ｓ’_t（τ）＝Ｓ_t（τ） ・・・式４５
になるとは限らないという問題がある。ここで、Ｎ_FFTは、高速フーリエ変換を実施するときのデータ長であり、mod（τ−ｔ，Ｎ_FFT）は、τ−ｔをＮ_FFTで割ったときの余りを示している。 However, also in the case of performing the convolution operation processing using the fast Fourier transform, the result S ′ _t obtained by obtaining the solution of S _ω (ω) using frequency domain multiplication and then returning to the time domain by inverse Fourier transform. (Τ) is

Next, similar to the correlation operation, 'for all tau of _t (τ) S' obtained _{S t (τ) = S t} (τ) ··· Equation 45
There is a problem that it does not always become. Here, N _FFT is a data length when the fast Fourier transform is performed, and mod (τ−t, N _FFT ) indicates a remainder when τ−t is divided by N _FFT .

In order to avoid the above-described problem and obtain a solution satisfying the relational expression of Expression 45, the filter calculation processing unit 31 according to the present embodiment adds an appropriate number of 0 data to g and h, by increasing the data length N _FFT fast Fourier transform, 'value S of _t (tau)' _t S in the range of tau to be obtained _{(τ) = S t (τ} ) ··· equation 45
Make adjustments so that

Accordingly, when Rq ₀ is obtained, ₀ ≦≦ so that the solution obtained by the inverse Fourier transform becomes an accurate value (so that an error does not occur in the solution obtained by the calculation in the time domain only). ε _μ is calculated using only data in the range of t ≦ L−1.

の初期値であるσ₀の値

の値を算出する。ここで、σ_μは、後述するステップＳ．４におけるＲσ_μの計算に用いられる。 In addition, the filter calculation processing unit 31 uses σ _μ

The initial value of σ ₀

Is calculated. Here, σ _μ is step S. described later. Used to calculate Arushiguma _mu in 4.

  Further, the filter calculation processing unit 31 sets a parameter used for the convergence determination of the conjugate gradient method as η. η is a constant as a threshold parameter for determining convergence, and when the error becomes a value satisfying the convergence determination condition based on this parameter (when the following equation 48 is satisfied), the filter calculation processing unit 31 It is determined that the solution has converged, and the iterative calculation is terminated.

で表される。フィルタ算出処理部３１は、この判断式に基づいて解の収束を判断する。 Then, the filter calculation processing unit 31 performs a convergence condition determination process (step S.3). The judgment formula for the convergence condition is

It is represented by The filter calculation processing unit 31 determines the convergence of the solution based on this determination formula.

When the convergence condition is not satisfied (No in step S.3), the filter calculation processing unit 31 calculates Rσ _μ (step S.4).

Even when the calculation of Rσ _μ, it is possible to perform calculation by the convolution arithmetic processing falls into the calculation for determining the product of the above-mentioned matrix R and vector. Therefore, the filter calculation processing unit 31 obtains a solution of Rσ _{μ in} the frequency domain by fast Fourier transform, and then obtains a time domain value by performing inverse Fourier transform. Even in this case, the solution is obtained using only data in the range of 0 ≦ t ≦ L−1 so that no error occurs in the solution in the time domain.

を算出し、さらに、μに１を加える（μ＝μ＋１）処理を行う（ステップＳ．５）。 Then, the filter calculation processing unit 31 uses the calculated value of Rσ _μ and ε _μ and σ _μ ,

Further, the process of adding 1 to μ (μ = μ + 1) is performed (step S.5).

Then, the filter calculation processing unit 31 performs step S.1. Using the ε _μ and P calculated in step 5, the convergence condition determination process (step S.3) is repeatedly executed.

When the convergence condition is satisfied (Yes in step S.3), the filter calculation processing unit 31 performs step S.3. The value of q _μ calculated in step 5 is taken as the solution of the j-th column in the sound field reproduction filter Q (step S.6).

  Thereafter, the filter calculation processing unit 31 determines whether or not j = C, that is, whether or not the solution of all the columns (1 to C columns) in the sound field reproduction filter has been obtained (step S.7). If j is not C (No in step S.7), 1 is added to j and 0 is substituted for μ (step S.8), and then the process is performed in step S.7. The process proceeds to 2 and the above-described processing is repeatedly executed.

  When j is C (j = C) (Yes in step S.7), the filter calculation processing unit 31 determines that the values of the sound field reproduction filters of all the columns have been obtained, and the sound field The reproduction filter calculation process is terminated. The calculated sound field reproduction filter is recorded in the filter information recording unit 34 by the filter calculation processing unit 31.

[Sound environment reproduction processing of original sound field in reproduced sound field]
Next, a process for reproducing the acoustic environment of the laboratory that is the original sound field 6 in the laboratory space B2 that is the reproduced sound field 21 will be described. The reproduction of the acoustic environment of the original sound field 6 in the reproduced sound field 21 is realized by a sound field reproduction system.

  FIG. 5 is a block diagram showing a schematic configuration of the sound field reproduction system 1. As shown in FIG. 5, the sound field reproduction system 1 is installed in a laboratory space B <b> 2 that is a reproduction sound field 21, and a sound output device that outputs sound to the speaker unit 22. (Sound output device) 4, a computer (sound processing device) 40 that outputs an evaluation sound to the sound output device 4, an audio device (sound reproduction device) 41 that can output an acoustic signal, and an audio device And an audio input device 8 that outputs an acoustic signal reproduced by the computer 41 to the computer 40.

  Here, the audio output device 4 used in the sound field reproduction system 1, the audio input device 8, and the speaker unit 22 are the audio output device 4 used in the reproduction sound field acoustic system 20, the audio input device 8, Since it is the same structure as the speaker part 22 and is provided with the same function, the detailed description is abbreviate | omitted.

  The audio device 41 is a device that reads out music information recorded on a recording medium such as a CD, a DVD, and an MD and outputs it as an acoustic signal. The acoustic signal reproduced in the audio device 41 is input to the computer 40 via the audio input device 8.

  The computer 40 acquires an acoustic signal reproduced by the audio device 41 via the audio input device 8, and performs an acoustic process on the acquired acoustic signal, thereby generating an original sound from the acoustic signal output from the reproduction sound field 21. It has a role to reproduce the acoustic environment of the field 6.

  The computer 40 is roughly configured by a control unit 42, a convolution operation processing unit 43, and a sound field reproduction filter recording unit 44. The sound field reproduction filter recording unit 44 records the sound field reproduction filter calculated by the sound field reproduction filter calculation device 30 described above. The sound field reproduction filter recording unit 44 is configured by a general recording device.

  The convolution operation processing unit 43 performs acoustic processing on the acoustic signal input by the voice input device 8 under the control of the control unit 42. More specifically, acoustic processing is performed on the acoustic signal output from the speaker 23 of the speaker unit 22 by performing convolution calculation processing on the acoustic signal using the sound field reproduction filter recorded in the sound field reproduction filter recording unit 44. Apply.

  For acoustic signals that have been subjected to acoustic processing by the sound field reproduction filter by the convolution operation processing unit 43, the acoustic characteristics in the reproduced sound field 21 are canceled and the acoustic characteristics of the original sound field 6 are effectively applied in the reproduced sound field 21. The sound effect for making it add is added. Therefore, in the sound field reproduction system 1, by outputting an acoustic signal to which the sound field reproduction filter is applied from the speaker 22, it is possible to reproduce the acoustic characteristics of the original sound field 6 in the reproduced sound field 21. A feeling as if listening in the original sound field 6 can be obtained in the reproduced sound field 21.

  As described above, in the filter calculation processing unit 31 according to the embodiment, the calculation process of the sound field reproduction filter obtained from the acoustic characteristics of the original sound field 6 and the acoustic characteristics of the reproduced sound field 21 is performed in all calculation processes. Rather than using a calculation method that reversely converts the calculation result in the frequency domain to the time domain, the calculation method that uses the iterative solution method to find the sound field reproduction filter in the time domain is adopted in principle, so accuracy without errors It is possible to calculate a high sound field reproduction filter.

  Furthermore, when calculating the sound field reproduction filter using the iterative solution method, the solution in the time domain is obtained using the calculation result in the frequency domain only when the correlation calculation process and the convolution calculation process are performed. It is possible to reduce the memory load and the memory capacity used for arithmetic processing. For this reason, it becomes possible to aim at the speeding up of the whole process, and the efficiency of the calculation assets utilized for calculation processing.

  In particular, when using the result of the correlation calculation process and the result of the convolution calculation process calculated in the frequency domain in the time domain, the calculation is performed using only the solution obtained by converting to the frequency domain and the time domain. In order to prevent an error from occurring in the solution, the sound field reproduction filter calculation process is performed using only data in a predetermined range in which the solution values do not become mismatched among the solutions subjected to inverse Fourier transform into the time domain. For this reason, it is possible to effectively prevent an error or the like from being included in the finally calculated value of the sound field reproduction filter.

  The sound field reproduction filter calculation device and the sound field reproduction system according to the present invention have been described in detail with reference to the drawings, but the sound field reproduction filter calculation device and the sound field reproduction system according to the present invention have been described above. The form is not limited. It will be apparent to those skilled in the art that various changes and modifications can be made within the scope of the claims, and these are naturally within the technical scope of the present invention. Understood.

  For example, in the above-described embodiment, after the sound field reproduction filter is calculated by the sound field reproduction filter calculation device 30, the calculated sound field reproduction filter is recorded in the sound field reproduction filter recording unit 44 of the sound field reproduction system 1. Although the sound signal is subjected to sound processing, the sound field reproduction filter may be calculated in the sound field reproduction system 1.

  FIG. 6 shows a sound field reproduction system capable of performing a sound field reproduction filter calculation process by causing the control unit 56 of the computer 55 to process the function of the filter calculation processing unit 31 in the sound field reproduction filter calculation device 30. 5 is a block diagram showing 57. FIG.

  In the sound field reproduction system 57 shown in FIG. 6, compared to the sound field reproduction system 1 (see FIG. 5) described in the embodiment, an original sound field characteristic recording unit 32 and a reproduction sound field characteristic recording unit are included in the computer 55. 33 in that it is provided.

  In the sound field reproduction system 57, the control unit 56 calculates a sound field reproduction filter based on the acoustic characteristics recorded in the original sound field characteristic recording unit 32 and the acoustic characteristics recorded in the reproduced sound field characteristic recording unit 33. The calculated sound field reproduction filter is recorded in the sound field reproduction filter recording unit 44. And the control part 56 performs the convolution calculation process with respect to the input acoustic signal with respect to the convolution calculation process part 43 using the sound field reproduction filter recorded on the sound field reproduction filter recording part 44. FIG.

  As described above, in the sound field reproduction system 57, by performing the calculation process of the sound field reproduction filter, it is not necessary to separately provide the sound field reproduction filter calculation device 30, so that the system configuration can be simplified. Become. In addition, by causing the sound field reproduction system 57 to perform the calculation process of the sound field reproduction filter, it is possible to facilitate and speed up the process.

It is a block diagram which shows schematic structure of the original sound field acoustic system which concerns on embodiment. It is a block diagram which shows schematic structure of the reproduction sound field acoustic system which concerns on embodiment. It is a block diagram which shows schematic structure of the sound field reproduction filter calculation apparatus which concerns on embodiment. It is the flowchart which showed the calculation process of the sound field reproduction filter by the filter calculation process part which concerns on embodiment. It is a block diagram which shows schematic structure of the sound field reproduction system which concerns on embodiment. It is a block diagram which shows the other schematic structural example of the sound field reproduction system which concerns on embodiment. It is the figure which showed the relationship between the sound field reproduction filter in the general sound field reproduction system, the acoustic characteristic of an original sound field, and the acoustic characteristic of a reproduction sound field. (A) It is the figure which showed the structure of the general original sound field, (b) is the figure which showed the structure of the general reproduction sound field.

Explanation of symbols

DESCRIPTION OF SYMBOLS 1, 57 ... Sound field reproduction system 2 ... Original sound field sound system 3 ... Speaker part 4 (original sound field sound system) ... Sound output device (sound output device)
5 ... computer 6 (of the original sound field acoustic system) ... original sound field 7 ... microphone array 8 ... voice input device 12 ... microphone 15 (of microphone array) ... evaluation sound generating unit 16 ... impulse response calculating unit 17 ... impulse response recording unit 18 , 42, 56 ... control unit 20 ... reproduced sound field acoustic system 21 ... reproduced sound field 22 ... (reproduced sound field) speaker unit 23 ... (speaker) speaker 30 ... sound field reproduction filter calculation device 31 ... filter calculation processing unit 32 ... Original sound field characteristic recording unit 33 ... Reproduced sound field characteristic recording unit 34 ... Filter information recording units 40 and 55 ... Computer (acoustic processing device)
41. Audio equipment (sound reproduction device)
43 ... Convolution calculation processing unit 44 ... Sound field reproduction filter recording unit A1 ... Space A2 (conventional original sound field) ... Space B1 (conventional reproduction sound field) ... Space B2 (of the original sound field according to this embodiment) ... space c, s (of the reproduced sound field according to this embodiment) ... speakers FL, FR, RL, RR (installed in the conventional original sound field and reproduced sound field) ... (original sound field according to this embodiment) Loudspeaker m ... control point

Claims (2)

Sound field reproduction filter calculation that calculates the sound field reproduction filter to reproduce the sound field environment of the original sound field in the reproduction sound field based on the sound characteristic measured in the original sound field and the sound characteristic measured in the reproduction sound field Regarding the device,
The sound field reproduction filter calculation device includes:
By applying the iterative method in the time domain, the sound field reproduction filter calculation process is executed,
Only the correlation calculation process and the convolution calculation process that occur in the application of the iterative solution method, using the calculation method that reversely converts the calculation result to the time domain after performing the calculation process by the conversion process to the frequency domain,
In the calculation method, an error may not occur between the calculation result of the correlation calculation process or the convolution calculation process and the calculation result obtained when the correlation calculation process or the convolution calculation process is performed only in the time domain. Using only the data of the range, execute the calculation process of the sound field filter in the time domain,
The range in which the error does not occur is determined according to the number of stages of the sound field reproduction filter set in advance according to the reverberation time of the output sound to which the sound field reproduction filter is applied. Reproduction filter calculation device. A speaker installed in the reproduction sound field;
An acoustic reproduction device for reproducing an acoustic signal to be output from the speaker;
The sound field reproduction filter calculation device according to claim 1,
An acoustic processing device that performs acoustic processing on the acoustic signal reproduced by the acoustic reproduction device using the sound field reproduction filter calculated by the sound field reproduction filter calculation device;
A sound field reproduction system comprising: an acoustic output device that outputs an acoustic signal acoustically processed by the acoustic processing device from the speaker .

Images