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

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

语音压缩

开发平台：

C/C++

GAINPRED.C：源码内容
  1. /* Version 3.3    Last modified: December 26, 1995 */
  3. /*---------------------------------------------------------------------------*
  4.  *  Gain_predict()        : make gcode0(exp_gcode0)                          *
  5.  *  Gain_update()         : update table of past quantized energies.         *
  6.  *  Gain_update_erasure() : update table of past quantized energies.         *
  7.  *                                                        (frame erasure)    *
  8.  *    This function is used both Coder and Decoder.                          *
  9.  *---------------------------------------------------------------------------*/
  11. #include "typedef.h"
  12. #include "basic_op.h"
  13. #include "ld8k.h"
  14. #include "tab_ld8k.h"
  15. #include "oper_32b.h"
  17. /*---------------------------------------------------------------------------*
  18.  * Function  Gain_predict                                                    *
  19.  * ~~~~~~~~~~~~~~~~~~~~~~                                                    *
  20.  * MA prediction is performed on the innovation energy (in dB with mean      *
  21.  * removed).                                                                 *
  22.  *---------------------------------------------------------------------------*/
  23. void Gain_predict(
  24.    Word16 past_qua_en[],/* (i) Q10 :Past quantized energies        */
  25.    Word16 code[],       /* (i) Q13 :Innovative vector.             */
  26.    Word16 L_subfr,      /* (i)     :Subframe length.               */
  27.    Word16 *gcode0,      /* (o) Qxx :Predicted codebook gain        */
  28.    Word16 *exp_gcode0   /* (o)     :Q-Format(gcode0)               */
  29. )
  30. {
  31.    Word16  i, exp, frac;
  32.    Word32  L_tmp;
  34.   /*-------------------------------*
  35.    * Energy coming from code       *
  36.    *-------------------------------*/
  38.    L_tmp = 0;
  39.    for(i=0; i<L_subfr; i++)
  40.      L_tmp = L_mac(L_tmp, code[i], code[i]);
  42.   /*-----------------------------------------------------------------*
  43.    *  Compute: means_ener - 10log10(ener_code/ L_sufr)               *
  44.    *  Note: mean_ener change from 36 dB to 30 dB because input/2     *
  45.    *                                                                 *
  46.    * = 30.0 - 10 log10( ener_code / lcode)  + 10log10(2^27)          *
  47.    *                                          !!ener_code in Q27!!   *
  48.    * = 30.0 - 3.0103 * log2(ener_code) + 10log10(40) + 10log10(2^27) *
  49.    * = 30.0 - 3.0103 * log2(ener_code) + 16.02  + 81.278             *
  50.    * = 127.298 - 3.0103 * log2(ener_code)                            *
  51.    *-----------------------------------------------------------------*/
  53.    Log2(L_tmp, &exp, &frac);               /* Q27->Q0 ^Q0 ^Q15       */
  54.    L_tmp = Mpy_32_16(exp, frac, -24660);   /* Q0 Q15 Q13 -> ^Q14     */
  55.                                            /* hi:Q0+Q13+1            */
  56.                                            /* lo:Q15+Q13-15+1        */
  57.                                            /* -24660[Q13]=-3.0103    */
  58.    L_tmp = L_mac(L_tmp, 32588, 32);        /* 32588*32[Q14]=127.298  */
  60.   /*-----------------------------------------------------------------*
  61.    * Compute gcode0.                                                 *
  62.    *  = Sum(i=0,3) pred[i]*past_qua_en[i] - ener_code + mean_ener    *
  63.    *-----------------------------------------------------------------*/
  65.    L_tmp = L_shl(L_tmp, 10);                      /* From Q14 to Q24 */
  66.    for(i=0; i<4; i++)
  67.      L_tmp = L_mac(L_tmp, pred[i], past_qua_en[i]); /* Q13*Q10 ->Q24 */
  69.    *gcode0 = extract_h(L_tmp);                    /* From Q24 to Q8  */
  71.   /*-----------------------------------------------------------------*
  72.    * gcode0 = pow(10.0, gcode0/20)                                   *
  73.    *        = pow(2, 3.3219*gcode0/20)                               *
  74.    *        = pow(2, 0.166*gcode0)                                   *
  75.    *-----------------------------------------------------------------*/
  77.    L_tmp = L_mult(*gcode0, 5439);       /* *0.166 in Q15, result in Q24*/
  78.    L_tmp = L_shr(L_tmp, 8);             /* From Q24 to Q16             */
  79.    L_Extract(L_tmp, &exp, &frac);       /* Extract exponent of gcode0  */
  81.    *gcode0 = extract_l(Pow2(14, frac)); /* Put 14 as exponent so that  */
  82.                                         /* output of Pow2() will be:   */
  83.                                         /* 16768 < Pow2() <= 32767     */
  84.    *exp_gcode0 = sub(14,exp);
  85. }
  88. /*---------------------------------------------------------------------------*
  89.  * Function  Gain_update                                                     *
  90.  * ~~~~~~~~~~~~~~~~~~~~~~                                                    *
  91.  * update table of past quantized energies                                   *
  92.  *---------------------------------------------------------------------------*/
  93. void Gain_update(
  94.    Word16 past_qua_en[],/* (io) Q10 :Past quantized energies        */
  95.    Word32  L_gbk12   /* (i) Q13 : gbk1[indice1][1]+gbk2[indice2][1] */
  96. )
  97. {
  98.    Word16  i, tmp;
  99.    Word16  exp, frac;
  100.    Word32  L_acc;
  102.    for(i=3; i>0; i--){
  103.       past_qua_en[i] = past_qua_en[i-1];         /* Q10 */
  104.    }
  105.   /*----------------------------------------------------------------------*
  106.    * -- past_qua_en[0] = 20*log10(gbk1[index1][1]+gbk2[index2][1]); --    *
  107.    *    2 * 10 log10( gbk1[index1][1]+gbk2[index2][1] )                   *
  108.    *  = 2 * 3.0103 log2( gbk1[index1][1]+gbk2[index2][1] )                *
  109.    *  = 2 * 3.0103 log2( gbk1[index1][1]+gbk2[index2][1] )                *
  110.    *                                                 24660:Q12(6.0205)    *
  111.    *----------------------------------------------------------------------*/
  113.    Log2( L_gbk12, &exp, &frac );               /* L_gbk12:Q13       */
  114.    L_acc = L_Comp( sub(exp,13), frac);         /* L_acc:Q16           */
  115.    tmp = extract_h( L_shl( L_acc,13 ) );       /* tmp:Q13           */
  116.    past_qua_en[0] = mult( tmp, 24660 );        /* past_qua_en[]:Q10 */
  117. }
  120. /*---------------------------------------------------------------------------*
  121.  * Function  Gain_update_erasure                                             *
  122.  * ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~                                             *
  123.  * update table of past quantized energies (frame erasure)                   *
  124.  *---------------------------------------------------------------------------*
  125.  *     av_pred_en = 0.0;                                                     *
  126.  *     for (i = 0; i < 4; i++)                                               *
  127.  *        av_pred_en += past_qua_en[i];                                      *
  128.  *     av_pred_en = av_pred_en*0.25 - 4.0;                                   *
  129.  *     if (av_pred_en < -14.0) av_pred_en = -14.0;                           *
  130.  *---------------------------------------------------------------------------*/
  131. void Gain_update_erasure(
  132.    Word16 past_qua_en[]/* (i) Q10 :Past quantized energies        */
  133. )
  134. {
  135.    Word16  i, av_pred_en;
  136.    Word32  L_tmp;
  138.    L_tmp = 0;                                                    /* Q10 */
  139.    for(i=0; i<4; i++)
  140.       L_tmp = L_add( L_tmp, L_deposit_l( past_qua_en[i] ) );
  141.    av_pred_en = extract_l( L_shr( L_tmp, 2 ) );
  142.    av_pred_en = sub( av_pred_en, 4096 );                          /* Q10 */
  144.    if( sub(av_pred_en, -14336) < 0 ){
  145.       av_pred_en = -14336;                              /* 14336:14[Q10] */
  146.    }
  148.    for(i=3; i>0; i--){
  149.       past_qua_en[i] = past_qua_en[i-1];
  150.    }
  151.    past_qua_en[0] = av_pred_en;
  152. }