Category

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

Voice Compress

Development Platform:

Unix_Linux

specdist.f：Code Content
  1. c==========================================================================
  2. c
  3. c ROUTINE
  4. c               specdist
  5. c
  6. c FUNCTION
  7. c
  8. c               Computes spectral distorion caused by quantization of
  9. c line spectral frequencies.
  10. c
  11. c SYNOPSIS
  12. c
  13. c subroutine specdist(unqfreq, newfreq, dm, sumdm, iframedm)
  14. c
  15. c   formal 
  16. c
  17. c                       data    I/O
  18. c       name            type    type    function
  19. c       -------------------------------------------------------------------
  20. c unqfreq real i unquantized line spectral frequencies
  21. c newfreq real i quantized line spectral frequencies
  22. c dm real o distortion array (subframe)
  23. c sumdm real o distortion array (current sum)
  24. c iframedm int o number of subframes  
  25. c
  26. c==========================================================================
  27. c
  28. c DESCRIPTION
  29. c
  30. c Calculate distortions/distances (log spectral error, etc.).
  31. c See the first reference below for a complete description.  A
  32. c "reference" system is compared against a "test" system.  Because
  33. c of the nonsymetric nature of the Itakura-Saito measure which some
  34. c of these distortion measures are based, poorer measures will be
  35. c obtained if the "reference" and "test" systems are reversed.
  36. c Because of gain uncertainties, a few measures are reported.
  37. c (Peter Kroon generally uses the measure DM(4).)
  38. c
  39. c==========================================================================
  40. c 
  41. c REFERENCES
  42. c
  43. c "Distance Measures for Speech Processing", A.H. Gray and J.D. Markel,
  44. c IEEE Trans. on ASSP, Vol. ASSP-24, no. 5, Oct. 1976
  45. c
  46. c "Quantization and Bit Allocation in Speech Processing",
  47. c A.H. Gray and J.D. Markel,IEEE Trans. on ASSP, Vol. ASSP-24
  48. c no. 6, Dec. 1976
  49. c
  50. c "A Note on Quantization and Bit Allocation in Speech Processing",
  51. c A.H. Gray and J.D. Markel,IEEE Trans. on ASSP, Vol. ASSP-25
  52. c no. 3, June 1977
  53. c
  54. c**************************************************************************
  55. c
  56. subroutine specdist(unqfreq, newfreq, dm, sumdm, iframedm)
  57. implicit undefined(a-z)
  58. include 'ccsub.com'
  59. convex #include "ccsub.com"
  60. real unqfreq(no), newfreq(no), dm(9), sumdm(10)
  61. integer iframedm
  62. c
  63. real unqpc(maxno+1), newpc(maxno+1), unqac(maxno+1), newac(maxno+1)
  64. real unqrc(maxno), newrc(maxno)
  65. c
  66. c Convert LSP's to autocorrelation coefficients for input to "dist"
  67. c
  68. call lsptopc(unqfreq, unqpc)
  69. call lsptopc(newfreq, newpc)
  70. call pctorc(unqpc, unqrc, no)
  71. call pctorc(newpc, newrc, no)
  72. call rctoac(unqrc, unqac, no)
  73. call rctoac(newrc, newac, no)
  74. c
  75. c *find distances between the ac sequences
  76. call dist(no,4*no,unqac,newac,dm,sumdm,iframedm)
  77. return
  78. end