Category

Development Platform
Home > SourceCode/Document > Communication-Mobile > Voice Compress > View Code
polefilt.f
Package: celp_3.2a.tar.gz　[view]
Upload User: szhypcb168
Upload Date: 2007-01-06
Package Size: 2187k
Code Size: 2k
Category:

Voice Compress

Development Platform:

Unix_Linux

polefilt.f：Code Content
  1. C==========================================================================
  2. C
  3. C ROUTINE
  4. C polefilt
  5. C
  6. C FUNCTION
  7. C Direct form all-pole filter
  8. C
  9. C SYNOPSIS
  10. C subroutine polefilt(a, n, z, xy, len)
  11. C
  12. C   formal 
  13. C
  14. C                       data    I/O
  15. C       name            type    type    function
  16. C       -------------------------------------------------------------------
  17. C a real i N+1 filter coefficients
  18. C n int i Filter order 
  19. C z real i N+1 filter delay elements
  20. C (maintained by the user)
  21. C (z(0) is a dummy delay)
  22. C xy real i/o Input/Output data array 
  23. C len int i Number of samples to filter
  24. C
  25. C==========================================================================
  26. C       
  27. C DESCRIPTION
  28. C
  29. C Recursive all-pole in-place time-domain filter.
  30. C The transfer function 1/A(z) is implemented
  31. C in a direct form realisation.
  32. C 
  33. C                   N       -i
  34. C       H(z) = 1 / SUM a(i)z         with a(0) = +1.0
  35. C                  i=0
  36. C
  37. C NOTE:  a(0) is not used for the actual computing,
  38. C as can easily be seen in the following flow graph.
  39. C
  40. C       x(t) ->---+------->--------(z0)-----> y(t)
  41. C                 |                  |
  42. C                 +-------< -a1 ----z1
  43. C                 |                  |
  44. C                 +-------< -a2 ----z2
  45. C                 |                  |
  46. C                 :                  :
  47. C                 |                  |
  48. C                 +-------< -aN ----zN
  49. C
  50. C==========================================================================
  51. C
  52. C CALLED BY
  53. C
  54. C celp confg impulse postfilt
  55. C
  56. C CALLS
  57. C
  58. C
  59. C
  60. C==========================================================================
  61. C
  62. C REFERENCES
  63. C
  64. C Oppenheim & Schafer, Digital Signal Processing, PH, 1975, p. 149.
  65. C
  66. C**************************************************************************
  67. C*-
  68. subroutine polefilt(a, n, z, xy, len)
  69. implicit undefined(a-z)
  70. integer n, len
  71. real a(0:n), z(0:n), xy(len)
  72. integer t, j
  73. c
  74. if (a(0) .ne. 1.0) stop 'polefilt:  bad coefficients'
  75. c
  76. do 69 t = 1, len
  77.    z(0) = xy(t)
  78.    do 10 j = n, 1, -1
  79.       z(0) = z(0) - a(j)*z(j)
  80.       z(j) = z(j-1)
  81. 10    continue
  82.    xy(t) = z(0)
  83. 69 continue
  84. return
  85. end