Как я могу адаптировать speex echo canceller для обработки выборок с плавающей запятой?

#echo #speex

#эхо #speex

Вопрос:

Добрый день! как я могу использовать выборки с плавающей точкой для обработки подавления эха? Я попытался изменить интерфейс и тело центральной функции:

От

 void speex_echo_cancellation(SpeexEchoState *st, const spx_int16_t *rec, const spx_int16_t *play, spx_int16_t *out);
  

Для

 void float_speex_echo_cancellation(SpeexEchoState *st, const float rec[], const float play[], float out[]);
  

и из

 ...
       for (i=0;i<st->frame_size;i  )
      {
         spx_word32_t tmp_out;
         tmp_out = SUB32(EXTEND32(st->input[chan*st->frame_size i]), EXTEND32(st->e[chan*N i st->frame_size]));
         tmp_out = ADD32(tmp_out, EXTEND32(MULT16_16_P15(st->preemph, st->memE[chan])));

         if (in[i*C chan] <= -32000 || in[i*C chan] >= 32000)
         {
         if (st->saturated == 0)
            st->saturated = 1;
         }

         **out[i*C chan] = (spx_int16_t)WORD2INT(tmp_out);**

         st->memE[chan] = tmp_out;
      }
...
  

Для

   ...

     for (i=0;i<st->frame_size;i  )
          {
             spx_word32_t tmp_out;
             tmp_out = SUB32(EXTEND32(st->input[chan*st->frame_size i]), EXTEND32(st->e[chan*N i st->frame_size]));
             tmp_out = ADD32(tmp_out, EXTEND32(MULT16_16_P15(st->preemph, st->memE[chan])));

             if (in[i*C chan] <= -32000 || in[i*C chan] >= 32000)
             {
             if (st->saturated == 0)
                st->saturated = 1;
             }

             **out[i*C chan] = /*(spx_int16_t)WORD2INT(*/tmp_out*/)*/;**

             st->memE[chan] = tmp_out;
          }
    ...
  

и из

 static inline void filter_dc_notch16(const spx_int16_t *in, spx_word16_t radius, spx_word16_t *out, int len, spx_mem_t *mem, int stride)
{
   int i;
   spx_word16_t den2;
   den2 = (spx_word16_t)(radius*radius   .7f*(1.f-radius)*(1.f-radius));
   for (i=0;i<len;i  )
   {
      spx_int16_t vin = in[i*stride];
      spx_word32_t vout = mem[0]   SHL32(EXTEND32(vin),15);
      mem[0] = mem[1]   2*(-vin   radius*vout);
      mem[1] = SHL32(EXTEND32(vin),15) - MULT16_32_Q15(den2,vout);
      out[i] = SATURATE32(PSHR32(MULT16_32_Q15(radius,vout),15),32767);
   }
}
  

Для

  static inline void float_filter_dc_notch16(const /*spx_int16_t*/spx_word16_t *in, spx_word16_t radius, spx_word16_t *out, int len, spx_mem_t *mem, int stride)
{
   int i;
   spx_word16_t den2;
   den2 = /*(spx_word16_t)*/(radius*radius   .7f*(1.f-radius)*(1.f-radius));
   for (i=0;i<len;i  )
   {
      /*spx_int16_t*/spx_word16_t vin = in[i*stride];
      spx_word32_t vout = mem[0]   SHL32(EXTEND32(vin),15);
      mem[0] = mem[1]   2*(-vin   radius*vout);
      mem[1] = SHL32(EXTEND32(vin),15) - MULT16_32_Q15(den2,vout);
      out[i] = /*SATURATE32(*/PSHR32(MULT16_32_Q15(radius,vout),15)/*,32767)*/;
   }
}
  

Итак, я предотвратил преобразование выходного результата с плавающей запятой в short int, но теперь я получаю предупреждение:

speex_warning («Эхоподавитель начал вести себя странно и получил пощечину (сброс). Он клянется, что теперь будет вести себя нормально «.);

это указывает на параметр st-> screwed_up, имеющий 50 значений, и указывает на обнуление всех выборок:

 ...

 if (!(Syy>=0 amp;amp; Sxx>=0 amp;amp; See >= 0)
   || !(Sff < N*1e9 amp;amp; Syy < N*1e9 amp;amp; Sxx < N*1e9)
  )

{ st->screwed_up  = 50; for (i=0;iframe_size*C;i  ) out[i] = 0; }

...
  

Что я могу сделать?

 enter code here
  

Ответ №1:

Почему вы хотите использовать выборки с плавающей запятой? Стандартный линейный звук PCM представлен в виде целых выборок в соответствии с выбранным битрейтом — 8 бит, 16 бит и так далее.

Откуда вы получаете эти входные данные?

На вашем месте я бы просто преобразовал все, что у вас есть, в shorts и предоставил это Speex, чтобы он мог с этим работать.