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

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

语音压缩

开发平台：

Unix_Linux

lsptopc.f：源码内容
  1. c==========================================================================
  2. c
  3. c NAME
  4. c lsptopc 
  5. c
  6. c FUNCTION
  7. c
  8. c convert lsp frequencies to predictor coefficients
  9. c
  10. c SYNOPSIS
  11. c
  12. c subroutine lsptopc(f, pc)
  13. c 
  14. c   formal
  15. c data I/O
  16. c name type type function
  17. c -------------------------------------------------------------------
  18. c f real i lsp frequencies
  19. c pc real o LPC predictor coefficients
  20. c
  21. c   global 
  22. c                       data    I/O
  23. c       name            type    type    function
  24. c       -------------------------------------------------------------------
  25. c  /ccsub/       see description include file
  26. c
  27. c==========================================================================
  28. c
  29. c DESCRIPTION
  30. c
  31. c LSPTOPC converts line spectral frequencies to LPC predictor 
  32. c coefficients.  
  33. c
  34. c The analysis filter may be reconstructed:
  35. c
  36. c A(z) = 1/2 [ P(z) + Q(z) ]
  37. c
  38. c CELP's LPC predictor coefficient convention is:
  39. c              p+1         -(i-1)
  40. c       A(z) = SUM   a   z          where a  = +1.0
  41. c              i=1    i                    1
  42. c
  43. c**************************************************************************
  44. c-
  45. subroutine lsptopc(f, pc)
  46. implicit undefined(a-z)
  47. include 'ccsub.com'
  48. convex #include "ccsub.com"
  49. c
  50. real f(maxno),pc(maxno+1)
  51. real freq(maxno),p(maxno/2),q(maxno/2)
  52. real a(maxno/2+1),a1(maxno/2+1),a2(maxno/2+1)
  53. real b(maxno/2+1),b1(maxno/2+1),b2(maxno/2+1)
  54. real pi,xx,xf
  55. integer i,j,k,noh
  56. logical lspflag
  57. c
  58. c *** check input for ill-conditioned cases
  59. c
  60. if (f(1).le.0.0 .or. f(1).ge.0.5) print *,' lsptopc:  ',
  61.      & 'LSPs out of bounds; f(1)=',f(1),' at frame', frame
  62. lspflag = .false.
  63. do 69 i = 2, no
  64.    if (f(i) .le. f(i-1)) lspflag = .true.
  65.    if (f(i).le.0.0 .or. f(i).ge.0.5) print *,' lsptopc:  ',
  66.      &    'LSPs out of bounds; f(',i,')=',f(i),' at frame', frame
  67. 69 continue
  68. if (lspflag) print *,' lsptopc:  nonmonotonic LSPs at frame', frame
  69. c
  70. c *** initialization
  71. c
  72. pi=4.*atan(1.)
  73. noh=no/2
  74. do 10 j=1,no
  75.    freq(j)=f(j)
  76. 10 continue
  77. do 100 i=1,noh+1
  78.    a(i) =0.
  79.    a1(i)=0.
  80.    a2(i)=0.
  81.    b(i) =0.
  82.    b1(i)=0.
  83.    b2(i)=0.
  84. 100 continue
  85. c
  86. c *** lsp filter parameters
  87. c
  88. do 110 i=1,noh
  89.    p(i)=-2.*cos(2.*pi*freq(2*i-1))
  90.    q(i)=-2.*cos(2.*pi*freq(2*i))
  91. 110 continue
  92. c
  93. c *** impulse response of analysis filter
  94. c
  95. xf = 0.
  96. do 120 k=1,no+1
  97.    xx=0.
  98.    if (k.eq.1) xx=1.
  99.    a(1)=xx+xf
  100.    b(1)=xx-xf
  101.    xf=xx
  102.    do 130 i=1,noh
  103.       a(i+1)=a(i)+p(i)*a1(i)+a2(i)
  104.       b(i+1)=b(i)+q(i)*b1(i)+b2(i)
  105.       a2(i)=a1(i)
  106.       a1(i)=a(i)
  107.       b2(i)=b1(i)
  108.       b1(i)=b(i)
  109. 130    continue
  110.    if (k.eq.1) go to 120
  111.    pc(k-1)=-.5*(a(noh+1)+b(noh+1))
  112. 120 continue
  113. c
  114. c *** convert to CELP's predictor coefficient array configuration
  115. c
  116. do 140 i = no, 1, -1
  117.    pc(i+1) = -pc(i)
  118. 140 continue
  119. pc(1)=1.
  120. return
  121. end