Target sound enhancement device, noise estimation parameter learning device, target sound enhancement method, noise estimation parameter learning method, and program

Abstract

A noise estimation parameter learning device is provided according to which even in a large space causing a problem of the reverberation and the time frame difference, multiple microphones disposed at distant positions cooperate with each other, and a spectral subtraction method is executed, thereby allowing the target sound to be enhanced. A noise estimation parameter learning device for learning noise estimation parameters used to estimate noise included in observed signals through a plurality of microphones, the noise estimation parameter learning device comprising: a modeling part that models a probability distribution of observed signals of the predetermined microphone, models a probability distribution of time frame differences, and models a probability distribution of transfer function gains; a likelihood function setting part that sets a likelihood function pertaining to the time frame difference, and a likelihood function pertaining to the transfer function gain, based on the modeled probability distributions; and a parameter update part that alternately and repetitively updates two variables of two likelihood functions, and outputs the time frame difference and the transfer function gain that have converged, as the noise estimation parameters.

Claims

1. A target sound enhancement device, comprising: processing circuitry configured to implement an observed signal acquisition part that acquires observed signals from a plurality of microphones; a frequency transformation part that transforms the observed signals into frequency spectra using a time-frame shift with a predetermined shift width; a noise estimation part that associates (i) an observed signal from a predetermined microphone that is among the plurality of microphones and that is disposed closest to a target sound, (ii) a selected microphone that is among the plurality of microphones and that is different from the predetermined microphone, the selected microphone being disposed adjacent to a noise source (iii) a time frame difference that is caused according to an arrival-time difference between the arrival times of a noise from a noise source to the predetermined microphone and to the selected microphone, the arrival-time difference being equal to or more than the shift width and (iv) a transfer function gain caused according to the relative position difference between the predetermined microphone, the selected microphone and the noise source, with each other, and estimates noise included in observed signals through a plurality of the predetermined microphones; a filter generation part that generates a filter based at least on the estimated noise; and a filtering part that filters the observed signal obtained from the predetermined microphone through the filter.

2. The target sound enhancement device according to claim 1, wherein the observed signal of the predetermined microphone includes a target sound and noise, and the observed signal of the selected microphone includes noise.

3. The target sound enhancement device according to claim 2, wherein the observed signal is a signal obtained by frequency-transforming an acoustic signal collected by the microphone, and a difference of two arrival times is equal to or more than a shift width of the frequency transformation, the arrival times being an arrival time of the noise from the noise source to the predetermined microphone and an arrival time of the noise from the noise source to the selected microphone.

4. The target sound enhancement device according to claim 2, wherein the noise estimation part associates, with each other, a probability distribution of observed signals of the predetermined microphone, a probability distribution where a time frame difference caused according to a relative position difference between the predetermined microphone and the selected microphone and the noise source is modeled, and a probability distribution where a transfer function gain caused according to the relative position difference between the predetermined microphone and the selected microphone and the noise source is modeled, and estimates the noise included in the observed signals through the plurality of microphones.

5. The target sound enhancement device according to claim 4, wherein the noise estimation part associates two likelihood functions set with each other based on three probability distributions and estimates the noise included in the observed signals through the plurality of microphones, the three probability distributions being a probability distribution of observed signals of the predetermined microphone, a probability distribution where a time frame difference caused according to a relative position difference between the predetermined microphone and the selected microphone and the noise source is modeled, and a probability distribution where a transfer function gain caused according to the relative position difference between the predetermined microphone and the selected microphone and the noise source is modeled, a first likelihood function being based on at least the probability distribution where the time frame difference is modelled, a second likelihood function being based on at least the probability distribution where the transfer function gain is modeled.

6. The target sound enhancement device according to claim 5, wherein the noise estimation part alternately and repetitively updates a variable of the first likelihood function and a variable of the second likelihood function.

7. The target sound enhancement device according to claim 6, wherein the variable of the first likelihood function and the variable of the second likelihood function are updated with an assigned restriction that limits the transfer function gain to a nonnegative value.

8. The target sound enhancement device according to claim 7, wherein the probability distribution of the time frame difference is modeled with a Poisson distribution, and the probability distribution of the transfer function gain is modeled with an exponential distribution.

9. A noise estimation parameter learning device for learning noise estimation parameters used to estimate noise included in observed signals through a plurality of microphones, the noise estimation parameter learning device comprising: processing circuitry configured to implement a modeling part that models a probability distribution of observed signals of a predetermined microphone among the plurality of microphones, models a probability distribution of time frame differences caused according to a relative position difference between the predetermined microphone, a selected microphone and a noise source, and models a probability distribution of transfer function gains caused according to the relative position difference between the predetermined microphone, the selected microphone and the noise source; a likelihood function setting part that sets a likelihood function pertaining to the time frame difference, and a likelihood function pertaining to the transfer function gain, based on the modeled probability distributions; and a parameter update part that alternately and repetitively updates a variable of the likelihood function pertaining to the time frame difference and a variable of the likelihood function pertaining to the transfer function gain, and outputs the time frame difference and the transfer function gain that have been updated, as the noise estimation parameters.

10. The noise estimation parameter learning device according to claim 9, wherein the parameter update part comprises a transfer function gain update part that assigns a restriction for limiting the transfer function gain to a nonnegative value, and repetitively updates the variable of the likelihood function pertaining to the transfer function gain by a proximal gradient method.

11. The noise estimation parameter learning device according to claim 9, wherein the modeling part comprises: an observed signal modeling part that models the probability distribution of the observed signals with a Gaussian distribution; a time frame difference modeling part that models the probability distribution of the time frame differences with a Poisson distribution; and a transfer function gain modeling part that models the probability distribution of the transfer function gains with an exponential distribution.

12. A target sound enhancement method executed by a target sound enhancement device, the target sound enhancement method comprising: a step of acquiring observed signals from a plurality of microphones; a step of transforming the observed signals into frequency spectra using a time-frame shift with a predetermined shift width; a step of associating (i) an observed signal from a predetermined microphone that is among the plurality of microphones and that is disposed closest to a target sound, (ii) a selected microphone that is among the plurality of microphones and that is different from the predetermined microphone, the selected microphone being disposed adjacent to a noise source (iii) a time frame difference that is caused according to an arrival-time difference between the arrival times of a noise from a noise source to the predetermined microphone and to the selected microphone, the arrival-time difference being equal to or more than the shift width and (iv) a transfer function gain caused according to the relative position difference between the predetermined microphone, the selected microphone and the noise source, with each other, and of estimating noise included in observed signals through a plurality of the predetermined microphones; a step of generating a filter based at least on the estimated noise; and a step of filtering the observed signal obtained from the predetermined microphone through the filter.

13. A noise estimation parameter learning method executed by a noise estimation parameter learning device for learning noise estimation parameters used to estimate noise included in observed signals through a plurality of microphones, the noise estimation parameter learning method comprising: a step of modeling a probability distribution of observed signals of a predetermined microphone among the plurality of microphones, modeling a probability distribution of time frame differences caused according to a relative position difference between the predetermined microphone, a selected microphone and a noise source, and modeling a probability distribution of transfer function gains caused according to the relative position difference between the predetermined microphone, the selected microphone and the noise source; a step of setting a likelihood function pertaining to the time frame difference, and a likelihood function pertaining to the transfer function gain, based on the modeled probability distributions; and a step of alternately and repetitively updating a variable of the likelihood function pertaining to the time frame difference and a variable of the likelihood function pertaining to the transfer function gain, and of outputting the time frame difference and the transfer function gain that have been updated, as the noise estimation parameters.

14. A non-transitory computer readable medium that stores a program causing a computer to function as the target sound enhancement device according to claim 1.

15. A non-transitory computer readable medium that stores a program causing a computer to function as the noise estimation parameter learning device according to claim 9.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIG. 1 is a block diagram showing a configuration of a noise estimation parameter learning device of Embodiment 1;

(2) FIG. 2 is a flowchart showing an operation of the noise estimation parameter learning device of Embodiment 1;

(3) FIG. 3 is a flowchart showing an operation of a modeling part of Embodiment 1;

(4) FIG. 4 is a flowchart showing an operation of a likelihood function setting part of Embodiment 1;

(5) FIG. 5 is a flowchart showing an operation of a parameter update part of Embodiment 1;

(6) FIG. 6 is a block diagram showing a configuration of a target sound enhancement device of Embodiment 2;

(7) FIG. 7 is a flowchart showing an operation of the target sound enhancement device of Embodiment 2; and

(8) FIG. 8 is a block diagram showing a configuration of a target sound enhancement device of Modification 2.

DETAILED DESCRIPTION OF THE EMBODIMENTS

(9) Embodiments of the present invention are hereinafter described in detail. Components having the same functions are assigned the same numerals, and redundant description is omitted.

Embodiment 1

(10) Embodiment 1 solves the two problems. Embodiment 1 provides a technique of estimating the time frame difference and reverberation so as to cause microphones disposed at positions far apart in a large space to cooperate with each other for sound source enhancement. Specifically, the time frame difference and the reverberation (transfer function gain (Note *1)) are described in a statistical model, and are estimated with respect to a likelihood maximization reference for an observed signal. To model the reverberation that is caused by a distance sufficiently apart and cannot be described by instantaneous mixture, modeling is performed by convolution of the amplitude spectrum of the sound source and the transfer function gain in the time-frequency domain.

(11) (Note *1) The reverberation can be described as a transfer function in the frequency domain, and the gain thereof is called a transfer function gain.

(12) Hereinafter, referring to FIG. 1, a noise estimation parameter learning device in Embodiment 1 is described. As shown in FIG. 1, the noise estimation parameter learning device 1 in this embodiment includes a modeling part 11, a likelihood function setting part 12, and a parameter update part 13. In more detail, the modeling part 11 includes an observed signal modeling part 111, a time frame difference modeling part 112, and a transfer function gain modeling part 113. The likelihood function setting part 12 includes an objective function setting part 121, a logarithmic part 122, and a term factorization part 123. The parameter update part 13 includes a transfer function gain update part 131, a time frame difference update part 132, and a convergence determination part 133.

(13) Hereinafter, referring to FIG. 2, an overview of the operation of the noise estimation parameter learning device 1 in this embodiment is described.

(14) First, the modeling part 11 models the probability distribution of observed signals of a predetermined microphone (first microphone) among the plurality of microphones, models the probability distribution of time frame differences caused according to the relative position difference between the predetermined microphone, a freely selected microphone (m-th microphone) and a noise source, and models the probability distribution of transfer function gains caused according to the relative position difference between the predetermined microphone, the freely selected microphone and the noise source (S11).

(15) Next, the likelihood function setting part 12 sets a likelihood function pertaining to the time frame difference, and a likelihood function pertaining to the transfer function gain, based on the modeled probability distributions (S12).

(16) Next, the parameter update part 13 alternately and repetitively updates a variable of the likelihood function pertaining to the time frame difference and a variable of the likelihood function pertaining to the transfer function gain, and outputs the time frame difference and the transfer function gain that have converged, as the noise estimation parameters (S13).

(17) To describe the operation of the noise estimation parameter learning device 1 in further detail, required description is made in the following chapter <Preparation>.

(18) <Preparation>

(19) Now, an issue of estimating a target sound S.sup.(1).sub.ω,τ from observation through M microphones (M is an integer of two or more) is discussed. One or more of the microphones are assumed to be disposed (Note *2) at positions sufficiently apart from a microphone serving as a main one.

(20) (Note *2) a distance causing an arrival time difference equal to or more than the shift width of the short-time Fourier transform (STFT). That is, a distance causing the time frame difference in time-frequency analysis. For example, in a case where the microphone interval is 2 [m] or more with the sonic speed of C=340 [m/s], the sampling frequency of 48.0 [kHz] and the STFT shift width of 512, the time frame difference occurs. That is, this means that the observed signal is a signal obtained by frequency-transforming an acoustic signal collected by the microphone, and the difference of two arrival times is equal to or more than the shift width of the frequency transformation, the arrival times being the arrival time of the noise from the noise source to the predetermined microphone and the arrival time of the noise from the noise source to the freely selected microphone.

(21) The identification number of the predetermined microphone disposed closest to S.sup.(1).sub.ω,τ is assumed as one. Its observed signal X.sup.(1).sub.ω,τ is assumed to be obtained by Formula (1). It is assumed that in a space there are M−1 point noise sources (e.g., public-address announcement) or a group of point noise sources (e.g., the cheering by supporters)
[Formula 6]
S.sub.ω,τ.sup.(2, . . . ,M)

(22) It is also assumed that the m-th microphone is disposed adjacent to the m-th (m=2, . . . ,M) noise source. It is assumed that adjacent to the m-th microphone,
[Formula 7]
|S.sub.ω,τ.sup.(m)|>>|S.sub.ω,τ.sup.(1, . . . ,M,≠m)|
holds. It is also assumed that the observed signal X.sup.(m).sub.ω,τ can be approximately described as
[Formula 8]
X.sub.ω,τ.sup.(m)≈S.sub.ω,τ.sup.(m)  (7)
Formula (7) shows that the observed signal of the freely selected (m-th) microphone includes noise. It is assumed that the noise N.sub.ω,τ reaching the first microphone consists only of
[Formula 9]
S.sub.ω,τ.sup.(2, . . . ,M)

(23) The amplitude spectrum thereof can be approximately described as follows.

(24) $\left[\mathrm{Formula}10\right]$$\begin{array}{cc}.Math.N_{\omega ,\tau }.Math.\approx \underset{m=2}{\overset{M}{.Math.}}\underset{k=0}{\overset{K}{.Math.}}{a}_{\omega ,k}^{\left(m\right)}.Math.X_{\omega ,\tau -P_m-k}^{\left(m\right)}.Math.& \left(8\right)\end{array}$

(25) Here, P.sub.m∈N.sub.+ is the time frame difference in the time-frequency domain, the difference being caused according to the relative position difference between the first microphone, the m-th microphone and the noise source S(m).sub.ω,τ. Here, a.sup.(m).sub.ω,k∈R.sub.+ is the transfer function gain, which is caused according to the relative position difference between the first microphone, the m-th microphone and the noise source S.sup.(m).sub.ω,τ.

(26) Hereinafter, description of the reverberation due to convolution between the amplitude spectrum of the sound source
[Formula 11]
|X.sub.ω,τ−P.sub.m.sub.−k.sup.(m)|

(27) and the transfer function gain a.sup.(m).sub.ω,k in the time-frequency domain is illustrated in detail. In a case where the number of taps of impulse response is longer than the analysis width of short-time Fourier transform (STFT), the transfer characteristics cannot be described by instantaneous mixture in the time-frequency domain (Reference non-patent literature 1). For example, in a case where the sampling frequency is 48.0 [kHz] and the analysis width of STFT is 512, the time length of reverberation (impulse response) that can be described as instantaneous mixture is 10 [ms]. Typically, the reverberation time period in a sport field or a manufacturing factory is equal to or longer than this time length. Consequently, a simple instantaneous mixture model cannot be assumed. To describe a long reverberation approximately, the m-th sound source is assumed to arrive, with convolution of the amplitude spectrum of X.sup.(m).sub.ω,τ with the transfer function gain a.sup.(m).sub.ω,k in the time-frequency domain. Reference non-patent literature 1 describes this with complex spectral convolution. The present invention describes this with an amplitude spectrum for the sake of more simple description. (Reference non-patent literature 1: T. Higuchi and H. Kameoka, “Joint audio source separation and dereverberation based on multichannel factorial hidden Markov model”, in Proc MLSP 2014, 2014.)

(28) According to the above discussion, based on Formula (8), possible estimation of the time frame difference P.sub.2, . . . ,M of the noise sources and the transfer function gain
[Formula 12]
a.sub.1, . . . ,K.sup.(2, . . . ,M)
can, in turn, estimate the amplitude spectrum of noise. Consequently, the spectral subtraction method can be executed. That is, in this embodiment and Embodiment 2,

(29) $\left[\mathrm{Formula}13\right]$$\Theta =\left\{a_{1,\mathrm{.Math.},K}^{\left(2,\mathrm{.Math.},M\right)},P_{2,\mathrm{.Math.},M}\right\}$
is estimated, and the spectral subtraction method is executed, thereby allowing the target sound to be collected in the large space.

(30) First, it is assumed that Formula (1) holds even in the amplitude spectrum domain, and |X.sup.(1).sub.ω,τ| is approximately described as follows.
[Formula 14]
|X.sub.ω,τ.sup.(1)|=|S.sub.ω,τ.sup.(1)|+|N.sub.ω,τ|  (9)

(31) Here, to simplify the description, H.sub.ω.sup.(1) is omitted. To represent all frequency bins ω∈{1, . . . , Ω} and τ∈{1, . . . , T} at the same time, Formula (9) is represented with the following matrix operations.

(32) $\left[\mathrm{Formula}15\right]$$\begin{array}{cc}{X}_{\tau }^{\left(1\right)}\approx S_{\tau }^{\left(1\right)}+N_{\tau }& \left(10\right)\\ X_{\tau }^{\left(m\right)}\approx S_{\tau }^{\left(m\right)}& \left(11\right)\\ \begin{array}{c}{N}_{\tau }\approx \underset{m=2}{\overset{M}{.Math.}}\underset{k=0}{\overset{K}{.Math.}}{a}_k^{\left(m\right)}\circ X_{\tau -P_m-k}^{\left(m\right)}\\ \approx X_{\tau }a\end{array}& \left(12\right)\end{array}$
Note that ° is a Hadamard product. Here,
[Formula 16]
X.sub.τ.sup.(i)=(|X.sub.1,τ.sup.(i)|,|X.sub.2,τ.sup.(i)|, . . . ,|X.sub.Ω,τ.sup.(i)|).sup.T  (13)
S.sub.τ.sup.(i)=(|S.sub.1,τ.sup.(i)|,|S.sub.2,τ.sup.(i)|, . . . ,|S.sub.Ω,τ.sup.(i)|).sup.T  (14)
N.sub.τ=(|N.sub.1,τ|,|N.sub.2,τ|, . . . ,|N.sub.Ω,τ|).sup.T  (15)
a.sub.k.sup.(i)=(a.sub.1,k.sup.(i),a.sub.2,k.sup.(i), . . . ,a.sub.Ω,k.sup.(i)).sup.T  (16)
X.sub.τ=(X.sub.τ.sup.(2), . . . ,X.sub.τ.sup.(M))  (17)
X.sub.τ.sup.(m)=(diag(X.sub.τ−P.sub.m.sup.(m)), . . . ,diag(X.sub.τ−P.sub.m.sub.−K.sup.(m)))  (18)
a=(a.sup.(2), . . . ,a.sup.(M))  (19)
a.sup.(m)=(a.sub.0.sup.(m), . . . ,a.sub.K.sup.(m))  (20)

(33) diag(x) represents a diagonal matrix having a vector x as diagonal elements. Here, S.sup.(1).sub.ω,τ is often sparse in the time frame direction (the target sound is not present almost over the time period). In a specific example, it means that soccer ball kicking sounds and voices of referees are temporally short, and rarely occur. Consequently, over the most time period,
[Formula 17]
X.sub.τ.sup.(1)=N.sub.τ  (21)

(34) holds.

(35) <Detailed Operation of Modeling Part 11>

(36) Hereinafter, referring to FIG. 3, the details of the operation of the modeling part 11 are described. Data required for learning is input into the observed signal modeling part 111. Specifically, the observed signal
[Formula 18]
X.sub.τ.sup.(1)=N.sub.τ  (21)

(37) is input.

(38) The observed signal modeling part 111 models the probability distribution of the observed signal X.sup.(1).sub.τ of the predetermined microphone with a Gaussian distribution where NT is the average and a covariance matrix diag(G) is adopted
[Formula 19]
(N.sub.τ,diag(σ.sup.2))
(S111).

(39) $\left[\mathrm{Formula}20\right]$$\begin{array}{cc}{X}_{\tau }^{\left(1\right)}\sim \left(X_{\tau }^{\left(1\right)}.Math.N_{\tau },\mathrm{diag}\left(\sigma \right)\right)& \left(22\right)\\ \sim \frac{{.Math.\Lambda .Math.}^{1\text{/}2}}{{\left(2\pi \right)}^{\Omega \text{/}2}}\exp \left\{-\frac{1}{2}{\left(X_{\tau }^{\left(1\right)}-N_{\tau }\right)}^T\Lambda \left(X_{\tau }^{\left(1\right)}-N_{\tau }\right)\right\}& \left(23\right)\end{array}$
Here, A=(diag(σ)).sup.−1. σ=(σ.sub.1, . . . , σ.sub.Ω).sup.T is the power of X.sup.(1).sub.τ for each frequency, and is obtained by

(40) $\left[\mathrm{Formula}21\right]$$\begin{array}{cc}{\sigma }_{\omega }=\frac{1}{T}\underset{\tau =1}{\overset{T}{.Math.}}.Math.X_{\omega ,\tau }^{\left(1\right)}.Math.& \left(24\right)\end{array}$

(41) This is for the sake of correcting the difference of averages of amplitudes for the frequencies.

(42) The observed signal may be transformed from the time waveform into the complex spectrum using a method, such as STFT. As for the observed signal, in a case of batch learning, X.sup.(m).sub.ω,τ for M channels obtained by applying short-time Fourier transform to learning data is input. In a case of online learning, what is obtained by buffering data for T frames is input. Here, the buffer size is to be tuned according to the time frame difference and the reverberation length, and may be set to be about T=500.

(43) Microphone distance parameters, and signal processing parameters are input into the time frame difference modeling part 112. The microphone distance parameters include microphone distances ϕ.sub.2, . . . ,M, and the minimum value and the maximum value of the sound source distance estimated from the microphone distances ϕ.sub.2, . . . ,M
[Formula 22]
ϕ.sub.2, . . . ,M.sup.min,ϕ.sub.2, . . . ,M.sup.max

(44) The signal processing parameters include the number of frames K, the sampling frequency f.sub.s, the STFT analysis width, and the shift length f.sub.shift Here, K=15 and therearound are recommended. The signal processing parameters may be set in conformity with the recording environment. When the sampling frequency is 16.0 [kHz], the analysis width may be set to be about 512, and the shift length may be set to be about 256.

(45) The time frame difference modeling part 112 models the probability distribution of the time frame differences with a Poisson distribution (S112). In a case where the m-th microphone is disposed adjacent to the m-th noise source, P.sub.m can be approximately estimated by the distances between the first microphone and the m-th microphone. That is, provided that the distance between the first microphone and the m-th microphone is ϕ.sub.m, the sonic speed is C, the sampling frequency is f and the STFT shift width is f.sub.shift, the time frame difference D.sub.m is approximately obtained by

(46) $\left[\mathrm{Formula}23\right]$$\begin{array}{cc}{D}_m=\mathrm{round}\left\{\frac{{\varphi }_m}{C}.Math.\frac{f_s}{f_{\mathrm{shift}}}\right\}& \left(25\right)\end{array}$

(47) Here, round{•} indicates rounding off to an integer. However, in actuality, the distance between the m-th microphone and the m-th noise source is not zero. Consequently, P.sub.m may stochastically fluctuate in proximity to D.sub.m. To model this, the time frame difference modeling part 112 models the probability distribution of the time frame difference with a Poisson distribution having the average value D.sub.m (S112).

(48) 0$\left[\mathrm{Formula}24\right]$$\begin{array}{cc}\begin{array}{c}{P}_m\sim \mathrm{Poisson}\left(P_m.Math.D_m\right)\\ \sim \frac{D_m^{P_m}}{P_m!}\exp \left\{-D_m\right\}\end{array}& \left(26\right)\end{array}$

(49) Transfer function gain parameters are input into the transfer function gain modeling part 113. The transfer function gain parameters include the initial value of the transfer function gain,
[Formula 25]
a.sub.1, . . . ,Ω,1, . . . ,K.sup.(2, . . . ,M)

(50) the average value α.sub.k of the transfer function gain, the time attenuation weight β of the transfer function gain, and the step size λ. If there is any knowledge, the initial value of the transfer function gain may be set accordingly. On the contrary, without any knowledge, the value may be set to
[Formula 26]
a.sub.1, . . . ,Ω,1, . . . ,K.sup.(2, . . . ,M)=1.0

(51) Likewise, if there is any knowledge, α.sub.k may be set accordingly. Without any knowledge, to reduce α.sub.k according to frame passage, α.sub.k may be set as follows.
[Formula 27]
α.sub.k=max(α−βk,ε)  (27)

(52) Here, α is the value of α.sub.0, β is the attenuation weight according to frame passage, and ε is a small coefficient for preventing division by zero. As various parameters, α=1.0 or therearound, β=0.05, and λ=10.sup.−3 or therearound are recommended.

(53) The transfer function gain modeling part 113 models the probability distribution of the transfer function gains with an exponential distribution (S113). a.sup.(m).sub.ω,k is a positive real number. In general, the value of the transfer function gain increases with increase in time k. To model this, the transfer function gain modeling part 113 models the probability distribution of the transfer function gains with an exponential distribution having the average value α.sub.k (S113).

(54) $\left[\mathrm{Formula}28\right]$$\begin{array}{cc}\begin{array}{c}{a}_{\omega ,k}^{\left(m\right)}\sim \mathrm{Exponential}\left(a_{\omega ,k}^{\left(m\right)}.Math.{\alpha }_k\right)\\ \sim \frac{1}{{\alpha }_k}\exp \left\{-\frac{a_{\omega ,k}^{\left(m\right)}}{{\alpha }_k}\right\}\end{array}& \left(28\right)\end{array}$

(55) As described above, the probability distributions for the observed signal and each parameter can be defined. In this embodiment, the parameters are estimated by maximizing the likelihood.

(56) <Detailed Operation of Likelihood Function Setting Part 12>

(57) Hereinafter, referring to FIG. 4, the details of the operation of the likelihood function setting part 12 are described. Specifically, the objective function setting part 121 sets the objective function as follows, on the basis of the modeled probability distribution (S121).

(58) $\left[\mathrm{Formula}29\right]$$\begin{array}{cc}L=p\left(X_{1,\mathrm{.Math.},T},\Theta \right)& \left(29\right)\\ =p\left(X_{1,\mathrm{.Math.},T}.Math.\Theta \right)p\left(a_{1,\mathrm{.Math.},K}^{\left(2,\mathrm{.Math.},M\right)}\right)p\left(P_{2,\mathrm{.Math.},M}\right)& \left(30\right)\\ p\left(X_{1,\mathrm{.Math.},T}.Math.\Theta \right)=\underset{\tau =1}{\overset{T}{.Math.}}\left(X_{\tau }^{\left(1\right)}.Math.N_{\tau },\mathrm{diag}\left(\sigma \right)\right)& \left(31\right)\\ p\left(a_{1,\mathrm{.Math.},K}^{\left(2,\mathrm{.Math.},M\right)}\right)=\underset{\omega =1}{\overset{Q}{.Math.}}\underset{m=2}{\overset{M}{.Math.}}\underset{k=1}{\overset{K}{.Math.}}\mathrm{Exponential}\left(a_{\omega ,k}^{\left(m\right)}.Math.{\alpha }_k\right)& \left(32\right)\\ p\left(P_{2,\mathrm{.Math.},M}\right)=\underset{m=2}{\overset{M}{.Math.}}\mathrm{Poisson}\left(P_m.Math.D_m\right)& \left(33\right)\end{array}$

(59) Here,
[Formula 30]
a.sub.1, . . . ,K.sup.(2, . . . ,M)

(60) is required to have a nonnegative value. Consequently, this optimization is a multivariable maximization problem with a limitation of L as follows.

(61) $\left[\mathrm{Formula}31\right]$$\begin{array}{cc}\Theta \leftarrow \arg \underset{\Theta }{\max }L\mathrm{subject}\mathrm{to}0\le a_{1,\mathrm{.Math.},\Omega ,1,\mathrm{.Math.},K}^{\left(2,\mathrm{.Math.},M\right)}& \left(34\right)\end{array}$

(62) Here, L has a form of a product of probability value. Consequently, there is a possibility that underflow occurs during calculation. Accordingly, the fact that a logarithmic function is a monotonically increasing function is used, and the logarithms of both sides are taken. Specifically, the logarithmic part 122 takes logarithms of both sides of the objective function, and transforms Formulae (34) and (33) as follows (S122).

(63) $\left[\mathrm{Formula}32\right]$$\begin{array}{cc}\Theta \leftarrow \arg \underset{\Theta }{\max }ℒ\mathrm{subject}\mathrm{to}0\le a_{1,\mathrm{.Math.},\Omega ,1,\mathrm{.Math.},K}^{\left(2,\mathrm{.Math.},M\right)}& \left(35\right)\\ ℒ=\ln p\left(X_{1,\mathrm{.Math.},T}.Math.\Theta \right)+\ln p\left(a_{1,\mathrm{.Math.},K}^{\left(2,\mathrm{.Math.},M\right)}\right)+\ln p\left(P_{2,\mathrm{.Math.},M}\right)& \left(36\right)\end{array}$

(64) Here,
[Formula 33]
=ln(L)

(65) Each element can be described as follows.

(66) $\left[\mathrm{Formula}34\right]$$\begin{array}{cc}\ln p\left(X_{1,\mathrm{.Math.},T}.Math.\Theta \right)\propto -\frac{1}{2}\underset{\tau =1}{\overset{T}{.Math.}}{\left(X_{\tau }^{\left(1\right)}-X_{\tau }a\right)}^T\Lambda \left(X_{\tau }^{\left(1\right)}-X_{\tau }a\right)& \left(37\right)\\ \ln p\left(a_{1,\mathrm{.Math.},K}^{\left(2,\mathrm{.Math.},M\right)}\right)\propto \underset{\omega =1}{\overset{\Omega }{.Math.}}\underset{m=2}{\overset{M}{.Math.}}\underset{k=1}{\overset{K}{.Math.}}-\ln {\alpha }_k-\frac{a_k^{\left(m\right)}}{{\alpha }_k}& \left(38\right)\\ \ln p\left(P_{2,\mathrm{.Math.},M}\right)\propto \underset{m=2}{\overset{M}{.Math.}}-\ln \left(P_m!\right)+P_m\ln \left(D_m\right)-D_m& \left(39\right)\end{array}$

(67) The above transformation facilitates maximization of each likelihood function constituting
[Formula 35]

(68) Formula (35) achieves maximization using the coordinate descent (CD) method. Specifically, the term factorization part 123 factorizes the likelihood function (logarithmic objective function) to a term related to a (a term related to the transfer function gain), and a term related to P (a term related to the time frame difference) (S123).
[Formula 36]
.sub.a=ln p(X.sub.1, . . . ,T|Θ)+ln p(a.sub.1, . . . ,K.sup.(2, . . . ,M))  (40)
.sub.P=ln p(X.sub.1, . . . ,T|Θ)+ln p(P.sub.(2, . . . ,M))  (41)

(69) Alternate optimization of each variable (repetitive update) approximately maximizes
[Formula 37]

(70) $\left[\mathrm{Formula}38\right]$$\begin{array}{cc}{a}_{1,\mathrm{.Math.},K}^{\left(2,\mathrm{.Math.},M\right)}\leftarrow \arg \underset{\Theta }{\max }{ℒ}_a\mathrm{subject}\mathrm{to}0\le a_{1,\mathrm{.Math.},\Omega ,1,\mathrm{.Math.},K}^{\left(2,\mathrm{.Math.},M\right)}& \left(42\right)\\ P_{2,\mathrm{.Math.},M}\leftarrow \arg \underset{\Theta }{\max }{ℒ}_P& \left(43\right)\end{array}$

(71) Formula (42) is optimization with the limitation. Accordingly, the optimization is achieved using the proximal gradient method.

(72) <Detailed Operation of Parameter Update Part 13>

(73) Hereinafter, referring to FIG. 5, the details of the operation of the parameter update part 13 are described. The transfer function gain update part 131 assigns a restriction that limits the transfer function gain to a nonnegative value, and repetitively updates the variable of the likelihood function pertaining to the transfer function gain by the proximal gradient method (S131).

(74) In more detail, the transfer function gain update part 131 obtains the gradient vector of
[Formula 39]
.sub.a with respect to a

(75) by the following formula.

(76) $\left[\mathrm{Formula}40\right]$$\begin{array}{cc}\frac{\partial {ℒ}_a}{\partial a}=\frac{1}{T}\underset{\tau =1}{\overset{T}{.Math.}}{X}_{\tau }^T\Lambda \left(-X_{\tau }^{\left(1\right)}+X_{\tau }a\right)-\alpha & \left(44\right)\\ \alpha =\left(\underset{\underset{M-1}{︸}}{\tilde{\alpha },\tilde{\alpha },\mathrm{.Math.},\tilde{\alpha }}\right)& \left(45\right)\\ \tilde{\alpha }=\left(\underset{\underset{\Omega }{︸}}{\frac{1}{{\alpha }_0},\mathrm{.Math.},\frac{1}{{\alpha }_0}},\underset{\underset{\Omega }{︸}}{\frac{1}{{\alpha }_1},\mathrm{.Math.},\frac{1}{{\alpha }_1}},\mathrm{.Math.},\underset{\underset{\Omega }{︸}}{\frac{1}{{\alpha }_K},\mathrm{.Math.},\frac{1}{{\alpha }_K}}\right)& \left(46\right)\end{array}$

(77) Execution is made by repetitive optimization of alternately performing the gradient method of Formula (47) and flooring of Formula (48).

(78) $\left[\mathrm{Formula}41\right]$$\begin{array}{cc}a\leftarrow a+\lambda \frac{\partial {ℒ}_a}{\partial a}& \left(47\right)\\ a_{1,\mathrm{.Math.},\Omega ,1,\mathrm{.Math.},K}^{\left(2,\mathrm{.Math.},M\right)}\leftarrow \max \left(0,a_{1,\mathrm{.Math.},\Omega ,1,\mathrm{.Math.},K}^{\left(2,\mathrm{.Math.},M\right)}\right)& \left(48\right)\end{array}$

(79) Here, λ is an update step size. The number of repetitions of the gradient method, i.e., Formulae (47) and (48), is about 30 in the case of the batch learning, and about one in the case of the online learning. The gradient of Formula (44) may be adjusted using an inertial term (Reference non-patent literature 2) or the like. (Reference non-patent literature 2: Hideki Asoh and other 7 authors, “ShinSo GakuShu, Deep Learning”, Kindai kagaku sha Co., Ltd., Nov. 2015).

(80) Formula (43) is combinatorial optimization of discrete variables. Accordingly, update is performed by grid searching. Specifically, the time frame difference update part 132 defines the possible maximum value and minimum value of P.sub.m for every in, evaluates, for every combination of the minimum and maximum for P.sub.m, the likelihood function related to the time frame difference
[Formula 42]
.sub.P

(81) and updates P.sub.m with the combination of maximizing the function (S132). For practical use, the minimum value
[Formula 43]
ϕ.sub.2, . . . ,M.sup.min
and the maximum value
[Formula 44]
ϕ.sub.2, . . . ,M.sup.max

(82) estimated from each microphone distance ϕ.sub.2, . . . ,M are input, and the possible maximum value and minimum value for P.sub.m may be calculated therefrom. The maximum value and the minimum value of the sound source distance is to be set in conformity with the environment, and may be set to about ϕ.sub.m.sup.min=ϕ.sub.m−20, and ϕ.sub.m.sup.max=ϕ.sub.m+20.

(83) The above update can be executed by a batch process of preliminarily estimating Θ using the learning data. In a case where an online process is intended, the observed signal may be buffered for a certain time period, and estimation of Θ may then be executed using the buffer.

(84) After Θ is successfully estimated by the above update, noise may be estimated by Formula (8), and the target sound may be enhanced by Formulae (4) and (5).

(85) The convergence determination part 133 determines whether the algorithm has converged or not (S133). As for the convergence condition, in the case of the batch learning, the determination method may be, for example, the sum of absolute values of the update amount of a.sup.(m).sub.ω,k, whether the learning times are equal to or more than a predetermined number (e.g., 1000 times) or the like. In the case of the online learning, dependent on the frequency of learning, the learning may be finished after a certain number of repetitions of learning (e.g., 1 to 5).

(86) When the algorithm converges (S133Y), the convergence determination part 133 outputs the converged time frame difference and transfer function gain as noise estimation parameter Θ.

(87) As described above, according to the noise estimation parameter learning device 1 of this embodiment, even in a large space causing a problem of the reverberation and the time frame difference, multiple microphones disposed at distant positions cooperate with each other, and the spectral subtraction method is executed, thereby allowing the target sound to be enhanced.

Embodiment 2

(88) In Embodiment 2, a target sound enhancement device that is a device of enhancing the target sound on the basis of the noise estimation parameter Θ obtained in Embodiment 1 is described. Referring to FIG. 6, the configuration of the target sound enhancement device 2 of this embodiment is described. As shown in FIG. 6, the target sound enhancement device 2 of this embodiment includes a noise estimation part 21, a time-frequency mask generation part 22, and a filtering part 23. Hereinafter, referring to FIG. 7, the operation of the target sound enhancement device 2 of this embodiment is described.

(89) Data required for enhancement is input into the noise estimation part 21. Specifically, the observed signal
[Formula 45]
X.sub.1, . . . ,Ω,τ.sup.(1, . . . ,M)

(90) and the noise estimation parameter Θ are input. The observed signal may be transformed from the time waveform into the complex spectrum using a method, such as STFT. Note that, for m=M, the spectrum
[Formula 46]
X.sub.1, . . . ,Ω,τ−P.sub.m.sub.−−K, . . . ,τ−P.sub.m.sup.(2, . . . ,M)

(91) buffered according to the time frame difference P, and the number of frames K of the transfer function gain are input.

(92) The noise estimation part 21 estimates noise included in the observed signals through M (multiple) microphones on the basis of the observed signals and the noise estimation parameter Θ by Formula (8) (S21).

(93) The noise estimation parameter Θ and Formula (8) may be construed as a parameter and formula where an observed signal from the predetermined microphone among the plurality of microphones, the time frame difference caused according to the relative position difference between the predetermined microphone, the freely selected microphone that is among the plurality of microphones and is different from the predetermined microphone and the noise source, and the transfer function gain caused according to the relative position difference between the predetermined microphone, the freely selected microphone and the noise source, are associated with each other.

(94) The target sound enhancement device 2 may have a configuration independent of the noise estimation parameter learning device 1. That is, independent of the noise estimation parameter Co, according to Formula (8), the noise estimation part 21 may associate the observed signal from the predetermined microphone among the plurality of microphones, the time frame difference caused according to the relative position difference between the predetermined microphone, the freely selected microphone that is among the plurality of microphones and is different from the predetermined microphone and the noise source, and the transfer function gain caused according to the relative position difference between the predetermined microphone, the freely selected microphone and the noise source, with each other, and estimate noise included in observed signals through a plurality of the predetermined microphones.

(95) The time-frequency mask generation part 22 generates the time-frequency mask G.sub.ω,τ based on the spectral subtraction method by Formula (4), on the basis of the observed signal |X.sup.(1).sub.ω,τ| of the predetermined microphone and the estimated noise |N.sub.ω,τ| (S22). The time-frequency mask generation part 22 may be called a filter generation part. The filter generation part generates a filter, based at least on the estimated noise by Formula (4) or the like.

(96) The filtering part 23 filters the observed signal |X.sup.(1).sub.ω,τ| of the predetermined microphone on the basis of the generated time-frequency mask G.sub.ω,τ (Formula (5)), and obtains and outputs an acoustic signal (complex spectrum Y.sub.ω,τ) where the sound (target sound) present adjacent to the predetermined microphone is enhanced (S23). To return the complex spectrum Y.sub.ω,τ to the waveform, inverse short-time Fourier transform (ISTFT) or the like may be used, or the function of ISTFT may be implemented in the filtering part 23.

(97) [Modification 1]

(98) Embodiment 2 has the configuration where the noise estimation part 21 receives (accepts) the noise estimation parameter Θ from another device (noise estimation parameter learning device 1) as required. It is a matter of course that another mode of the target sound enhancement device can be considered. For example, as a target sound enhancement device 2a of Modification 1 shown in FIG. 8, the noise estimation parameter Θ may be preliminarily received from the other device (noise estimation parameter learning device 1), and preliminarily stored in a parameter storage part 20.

(99) In this case, the parameter storage part 20 preliminarily stores and holds the time frame difference and transfer function gain having been converged by alternately and repetitively updating the variables of the two likelihood functions set based on the three probability distributions described above, as the noise estimation parameter Θ.

(100) As described above, according to the target sound enhancement devices 2 and 2a of this embodiment and this modification, even in the large space causing the problem of the reverberation and the time frame difference, the multiple microphones disposed at distant positions cooperate with each other, and the spectral subtraction method is executed, thereby allowing the target sound to be enhanced.

(101) <Supplement>

(102) The device of the present invention includes, as a single hardware entity, for example: an input part to which a keyboard and the like can be connected; an output part to which a liquid crystal display and the like can be connected; a communication part to which a communication device (e.g., a communication cable) communicable with the outside of the hardware entity can be connected; a CPU (Central Processing Unit, which may include a cache memory and a register); a RAM and a ROM, which are memories; an external storage device that is a hard disk; and a bus that connects these input part, output part, communication part, CPU, RAM, ROM and external storing device to each other in a manner allowing data to be exchanged therebetween. The hardware entity may be provided with a device (drive) capable of reading and writing from and to a recording medium, such as CD-ROM, as required. A physical entity including such a hardware resource may be a general-purpose computer or the like.

(103) The external storage device of the hardware entity stores programs required to achieve the functions described above and data required for the processes of the programs (not limited to the external storage device; for example, programs may be stored in a ROM, which is a storage device dedicated for reading, for example). Data and the like obtained by the processes of the programs are appropriately stored in the RANI or the external storage device.

(104) In the hardware entity, each program stored in the external storage device (or a ROM etc.), and data required for the process of each program are read into the memory, as required, and are appropriately subjected to analysis, execution and processing by the CPU. As a result, the CPU achieves predetermined functions (each component represented as . . . part, . . . portion, etc. described above).

(105) The present invention is not limited to the embodiments described above, and can be appropriately changed in a range without departing from the spirit of the present invention. The processes described in the above embodiments may be executed in a time series manner according to the described order. Alternatively, the processes may be executed in parallel or separately, according to the processing capability of the device that executes the processes, or as required.

(106) As described above, in a case where the processing functions of the hardware entity (the device of the present invention) described in the embodiments are achieved by a computer, the processing details of the functions to be held by the hardware entity are described in a program. The program is executed by the computer, thereby achieving the processing functions in the hardware entity on the computer.

(107) The program that describes the processing details can be recorded in a computer-readable recording medium. The computer-readable recording medium may be, for example, any of a magnetic recording device, an optical disk, a magneto-optical recording medium, a semiconductor memory and the like. Specifically, for example, a hard disk device, a flexible disk, a magnetic tape and the like may be used as the magnetic recording device. A DVD (Digital Versatile Disc), a DVD-RAM (Random Access Memory), a CD-ROM (Compact Disc Read Only Memory), CD-R (Recordable)/RW (ReWritable) and the like may be used as the optical disk. An MO (Magneto-Optical disc) and the like may be used as the magneto-optical recording medium. An EEP-ROM (Electronically Erasable and Programmable-Read Only Memory) and the like may be used as the semiconductor memory.

(108) For example, the program may be distributed by selling, assigning, lending and the like of portable recording media, such as a DVD and a CD-ROM, which record the program. Alternatively, a configuration may be adopted that distributes the program by storing the program in the storage device of the server computer and then transferring the program from the server computer to another computer via a network.

(109) For example, the computer that executes such a program temporarily stores, in the own storage device, the program stored in the portable recording medium or the program transferred from the server computer. During execution of the process, the computer reads the program stored in the own recording medium, and executes the process according to the read program. Alternatively, according to another execution mode of the program, the computer may directly read the program from the portable recording medium, and execute the process according to the program. Further alternatively, every time the program is transferred to this computer from the server computer, the process according to the received program may be sequentially executed. Alternatively, a configuration may be adopted that does not transfer the program to this computer from the server computer but executes the processes described above by what is called an ASP (Application Service Provider) service that achieves the processing functions only through execution instructions and result acquisition. It is assumed that the program of this mode includes information that is to be provided for the processes by a computer and is equivalent to the program (data and the like having characteristics that are not direct instructions to the computer but define the processes of the computer).

(110) In this mode, the hardware entity can be configured by executing a predetermined program on the computer. Alternatively, at least one or some of the processing details may be achieved by hardware.