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

源码开发语言/平台
当前位置：首页 > 源码/资料 > 通讯/手机编程 > 语音压缩 > 查看源码
gainncd2.f
资源名称：celp_3.2a.tar.gz　[点击查看]
上传用户：szhypcb168
上传日期：2007-01-06
资源大小：2187k
文件大小：2k
源码类别：

语音压缩

开发平台：

Unix_Linux

gainncd2.f：源码内容
  1. c==========================================================================
  2. c
  3. c ROUTINE
  4. c gainencode2
  5. c
  6. c FUNCTION
  7. c  
  8. c encode and quantize code book gain
  9. c
  10. c SYNOPSIS
  11. c subroutine gainencode2(numer, denom, bits, type, index)
  12. c
  13. c   formal 
  14. c
  15. c                       data    I/O
  16. c       name            type    type    function
  17. c       -------------------------------------------------------------------
  18. c numer i i code book gain numerator
  19. c denom i i code book gain denominator
  20. c bits i i # bits for encode
  21. c type c i type of quantization
  22. c index r o encoded code book gain ZERO BASED index
  23. c gainencode2 r f encoded code book gain
  24. c
  25. c==========================================================================
  26. c
  27. c DESCRIPTION
  28. c
  29. c   Fast nonuniform division for code book gain quantization to allow
  30. c practical quantization inside the code book search loop.  A binary
  31. c tree search with cross multiply quantization is implemented below.
  32. c
  33. C==========================================================================
  34. C 
  35. C INPUT FILES
  36. C cbgain.tbl code book gain coding
  37. C
  38. c**************************************************************************
  39. c
  40. real function gainencode2(numer, denom, bits, type, index)
  41. implicit undefined(a-z)
  42. c
  43. real numer, denom
  44. integer bits, index
  45. character*10 type
  46. c
  47. c Hard coded for 5 bit quantization to achieve high speed
  48. c
  49. integer bit, bitm1
  50. parameter (bit   = 5)
  51. parameter (bitm1 = bit-1)
  52. real gainlog5(2**bit), midpoints(2**bit-1)
  53. integer i, shift(bitm1)
  54. data shift/8, 4, 2, 1/
  55. c
  56. include './cbgain.tbl'
  57. convex #include "./cbgain.tbl"
  58. c
  59. ccc if (bits .ne. bit) stop ' gainencode2:  bad bits'
  60. c
  61. c *Binary tree search for closest gain
  62. c
  63. index = 2**bitm1
  64. do 69 i = 1, bitm1
  65.    if (numer .gt. denom*midpoints(index)) then
  66.       index = index + shift(i)
  67.    else
  68.       index = index - shift(i)
  69.    end if
  70. 69 continue
  71. if (numer .gt. denom*midpoints(index)) index = index + 1
  72. c
  73. c *Return quantized gain and ZERO based index
  74. c
  75. gainencode2= gainlog5(index)
  76. index      = index - 1
  78. return
  79. end