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

源码开发语言/平台
当前位置：首页 > 源码/资料 > 通讯/手机编程 > 语音压缩 > 查看源码
PRE_PROC.C
资源名称：G711-729.rar　[点击查看]
上传用户：meifeng08
上传日期：2013-06-18
资源大小：5304k
文件大小：3k
源码类别：

语音压缩

开发平台：

C/C++

PRE_PROC.C：源码内容
  1. /*
  2.    ITU-T G.729A Speech Coder    ANSI-C Source Code
  3.    Version 1.1    Last modified: September 1996
  5.    Copyright (c) 1996,
  6.    AT&T, France Telecom, NTT, Universite de Sherbrooke
  7.    All rights reserved.
  8. */
  10. /*------------------------------------------------------------------------*
  11.  * Function Pre_Process()                                                 *
  12.  *                                                                        *
  13.  * Preprocessing of input speech.                                         *
  14.  *   - 2nd order high pass filter with cut off frequency at 140 Hz.       *
  15.  *   - Divide input by two.                                               *
  16.  *-----------------------------------------------------------------------*/
  18. #include "typedef.h"
  19. #include "basic_op.h"
  20. #include "oper_32b.h"
  22. #include "ld8a.h"
  23. #include "tab_ld8a.h"
  25. /*------------------------------------------------------------------------*
  26.  * 2nd order high pass filter with cut off frequency at 140 Hz.           *
  27.  * Designed with SPPACK efi command -40 dB att, 0.25 ri.                  *
  28.  *                                                                        *
  29.  * Algorithm:                                                             *
  30.  *                                                                        *
  31.  *  y[i] = b[0]*x[i]/2 + b[1]*x[i-1]/2 + b[2]*x[i-2]/2                    *
  32.  *                     + a[1]*y[i-1]   + a[2]*y[i-2];                     *
  33.  *                                                                        *
  34.  *     b[3] = {0.92727435E+00, -0.18544941E+01, 0.92727435E+00};          *
  35.  *     a[3] = {0.10000000E+01, 0.19059465E+01, -0.91140240E+00};          *
  36.  *                                                                        *
  37.  *  Input are divided by two in the filtering process.                    *
  38.  *-----------------------------------------------------------------------*/
  41. /* Static values to be preserved between calls */
  42. /* y[] values is keep in double precision      */
  44. static Word16 y2_hi, y2_lo, y1_hi, y1_lo, x0, x1;
  47. /* Initialization of static values */
  49. void Init_Pre_Process(void)
  50. {
  51.   y2_hi = 0;
  52.   y2_lo = 0;
  53.   y1_hi = 0;
  54.   y1_lo = 0;
  55.   x0   = 0;
  56.   x1   = 0;
  57. }
  60. void Pre_Process(
  61.   Word16 signal[],    /* input/output signal */
  62.   Word16 lg)          /* length of signal    */
  63. {
  64.   Word16 i, x2;
  65.   Word32 L_tmp;
  67.   for(i=0; i<lg; i++)
  68.   {
  69.      x2 = x1;
  70.      x1 = x0;
  71.      x0 = signal[i];
  73.      /*  y[i] = b[0]*x[i]/2 + b[1]*x[i-1]/2 + b140[2]*x[i-2]/2  */
  74.      /*                     + a[1]*y[i-1] + a[2] * y[i-2];      */
  76.      L_tmp     = Mpy_32_16(y1_hi, y1_lo, a140[1]);
  77.      L_tmp     = L_add(L_tmp, Mpy_32_16(y2_hi, y2_lo, a140[2]));
  78.      L_tmp     = L_mac(L_tmp, x0, b140[0]);
  79.      L_tmp     = L_mac(L_tmp, x1, b140[1]);
  80.      L_tmp     = L_mac(L_tmp, x2, b140[2]);
  81.      L_tmp     = L_shl(L_tmp, 3);      /* Q28 --> Q31 (Q12 --> Q15) */
  82.      signal[i] = round(L_tmp);
  84.      y2_hi = y1_hi;
  85.      y2_lo = y1_lo;
  86.      L_Extract(L_tmp, &y1_hi, &y1_lo);
  87.   }
  88.   return;
  89. }