Проблема реализации алгоритма порога энергии для обнаружения голосовой активности - PullRequest
3 голосов
/ 20 мая 2011

Я пытаюсь реализовать алгоритм пороговых значений энергии для обнаружения голосовой активности, но не получаю значимых значений энергии для кадров размером wL.

wL = 1784 // about 40 ms (
const double decay_constant = 0.90 // some optimal value between 0 and 1
double prevrms = 1.0  // avoid DivideByZero
double threshold = some optimal value after some experimentation

for (int i = 0; i < noSamples ; i += wL)
{
   for (int j = 0; j < wL; j++)
   {
     // Exponential decay
     total = total * decay_constant;
     total += (audioSample[j] * audioSample[j]); // sum of squares
   }

   double mean = total / wL;
   double rms = Math.Round(Math.Sqrt(mean),2); // root mean sqare
   double prevrms = 1.0;

   if(rms/prevrms > threshold)
   {
  // voice detected
   }

   prevrms = rms;
   rms = 0.0;
}

Что не так с вышеуказанной реализацией?rms вычисляется для каждого кадра как 0,19.

Другая проблема - это скорость, так как выполнение описанного выше заняло около 30 минут.В настоящее время реализация O (n 2 ).Я работаю с автономными данными, так что это не так уж важно - точность является главной целью, но любые предложения по повышению эффективности будут высоко оценены.

Кроме того, я должен использовать другие факторы, такие как автокорреляция и скорость пересечения нуля, или достаточно одной энергии?

Ниже приводится краткое описание файла WAV (только с учетом чистой разговорной речи), который я использую:

// WAV file information
Sampling Frequency: 44100     Bits Per Sample:  16 
Channels: 2    nBlockAlign: 4   wavdata size: 557941248 bytes
Duration: 3162.932 sec    Samples: 139485312    Time between samples: 0.0227 ms
Byte position at start of samples: 44 bytes  (0x2C)

Chosen first sample to display:  1   (0.000 ms)
Chosen end  sample to display:  1784   (40.431 ms)

16 bit max possible value is:  32767  (0x7FFF)
16 bit min possible value is: -32768  (0x8000)

1 Ответ

1 голос
/ 24 мая 2011

Я нашел проблему. Мой второй цикл for был настроен неправильно. По сути, второй цикл for должен выглядеть примерно так:

for(j = i; j <= i + wL ;j++)

Вместо:

for(j = 0; j < wL; j++)

Который снова и снова повторял одни и те же значения выборки.

...