Detection of the number of speakers

図 1: An example of convergence to the true values

図 2: Input spectrum for Fig 1

It is generally known that ML estimates obtained by EM algorithm firmly depend on initial values and may often converge to undesirable values. To avoid this, we first prepare extra amount of harmonic-GMMs in the model in order to raise possibility of obtaining the true values. Then, obviously, the model may overfit the given observed specrum. If one Gaussian is enough for approximating the shape of one partial¹, the same number of underlying harmonic structures must be enough with the harmonic-GMMs. And this number can be detected by reducing harmonic-GMM one after another until they become the proper number on the basis of AIC. The specific operation is as follows:

1. Set initial values of in the limited frequency range.
2. Estimate the ML model parameters by EM algorithm. However, is constrained here as
   

   (11)

   
   This represents the degree of predominance of th harmonic-GMM. In Maximization-step, model parameters and should be updated to
   
   

   			(12)
   			(13)

   
   
   
   where is an integral of with respect to .
3. Calculate AIC with equation (10). Since there are two free parameters for each harmonic-GMM, the model has free parameters altogether. If the AIC increases, the number of harmonic-GMMs just before they are reduced in step4 will be the estimate of the number of harmonic structures.
4. Remove the harmonic-GMM(s) which conforms either of the two conditions as below and repeat from step 2.
   • The one whose is the minimum among all. Since the contribution to the maximum log-lik-
     elihood must be the least.
   • The one whose is smaller if the two adjacent representative means become closer than a certain distance (threshold). Since the two representative means are presumed to converge to the same optimal solution.

An example of how this process actually works is shown in Fig.1 where the observed spectrum used is depicted in Fig.2. The broken line represents the point where the model parameters were judged to be converged and the line graph indicates the value of AIC at each point. Since AIC takes minimum when harmonic-GMMs remain, the detected number here is .