开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 通讯/手机编程 > 语音压缩 > 查看源码
pre_proc.c
资源名称：g729a_c_linux.rar　[点击查看]
上传用户：zhouyunkk
上传日期：2013-01-10
资源大小：59k
文件大小：3k
源码类别：

语音压缩

开发平台：

C/C++

pre_proc.c：源码内容
  1. /*
  2.    ITU-T G.729 Annex C - Reference C code for floating point
  3.                          implementation of G.729
  4.                          Version 1.01 of 15.September.98
  5. */
  7. /*
  8. ----------------------------------------------------------------------
  9.                     COPYRIGHT NOTICE
  10. ----------------------------------------------------------------------
  11.    ITU-T G.729 Annex C ANSI C source code
  12.    Copyright (C) 1998, AT&T, France Telecom, NTT, University of
  13.    Sherbrooke.  All rights reserved.
  15. ----------------------------------------------------------------------
  16. */
  18. /*
  19.  File : PRE_PROC.C
  20.  Used for the floating point version of both
  21.  G.729 main body and G.729A
  22. */
  25. /*------------------------------------------------------------------------*
  26.  * Function pre_process()                                                 *
  27.  *                                                                        *
  28.  * Preprocessing of input speech.                                         *
  29.  *   - 2nd order high pass filter with cut off frequency at 140 Hz.       *
  30.  *-----------------------------------------------------------------------*/
  32. #include "typedef.h"
  33. #include "version.h"
  34. #ifdef VER_G729A
  35.  #include "ld8a.h"
  36.  #include "tab_ld8a.h"
  37. #else
  38.  #include "ld8k.h"
  39.  #include "tab_ld8k.h"
  40. #endif
  42. /*------------------------------------------------------------------------*
  43.  * 2nd order high pass filter with cut off frequency at 140 Hz.           *
  44.  * Designed with SPPACK efi command -40 dB att, 0.25 ri.                  *
  45.  *                                                                        *
  46.  * Algorithm:                                                             *
  47.  *                                                                        *
  48.  *  y[i] = b[0]*x[i] + b[1]*x[i-1] + b[2]*x[i-2]                          *
  49.  *                   + a[1]*y[i-1] + a[2]*y[i-2];                         *
  50.  *                                                                        *
  51.  *     b[3] = {0.92727435E+00, -0.18544941E+01, 0.92727435E+00};          *
  52.  *     a[3] = {0.10000000E+01, 0.19059465E+01, -0.91140240E+00};          *
  53.  *-----------------------------------------------------------------------*/
  56. static FLOAT x0, x1;         /* high-pass fir memory          */
  57. static FLOAT y1, y2;         /* high-pass iir memory          */
  59. void init_pre_process( void
  60. )
  61. {
  62.   x0 = x1 = (F)0.0;
  63.   y2 = y1 = (F)0.0;
  65.   return;
  66. }
  68. void pre_process(
  69.    FLOAT signal[],      /* (i/o)  : signal                     */
  70.    int lg               /* (i)    : lenght of signal           */
  71. )
  72. {
  73.   int i;
  74.   FLOAT x2;
  75.   FLOAT y0;
  77.   for(i=0; i<lg; i++)
  78.   {
  79.     x2 = x1;
  80.     x1 = x0;
  81.     x0 = signal[i];
  83.     y0 = y1*a140[1] + y2*a140[2] + x0*b140[0] + x1*b140[1] + x2*b140[2];
  85.     signal[i] = y0;
  86.     y2 = y1;
  87.     y1 = y0;
  88.   }
  90.   return;
  91. }