alaw_test.c
上传用户:sy_wanhua
上传日期:2013-07-25
资源大小:3048k
文件大小:6k
源码类别:

流媒体/Mpeg4/MP4

开发平台:

C/C++

  1. /*
  2. ** Copyright (C) 1999-2000 Erik de Castro Lopo <erikd@zip.com.au>
  3. **  
  4. ** This program is free software; you can redistribute it and/or modify
  5. ** it under the terms of the GNU General Public License as published by
  6. ** the Free Software Foundation; either version 2 of the License, or
  7. ** (at your option) any later version.
  8. ** 
  9. ** This program is distributed in the hope that it will be useful,
  10. ** but WITHOUT ANY WARRANTY; without even the implied warranty of
  11. ** MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  12. ** GNU General Public License for more details.
  13. ** 
  14. ** You should have received a copy of the GNU General Public License
  15. ** along with this program; if not, write to the Free Software 
  16. ** Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
  17. */
  18. #include <stdio.h>
  19. #include <unistd.h>
  20. #include <sndfile.h>
  21. #define BUFFER_SIZE (65536)
  22. static unsigned char alaw_encode (int sample) ;
  23. static int alaw_decode (unsigned int alawbyte) ;
  24. static short short_buffer [BUFFER_SIZE] ;
  25. static unsigned char alaw_buffer [BUFFER_SIZE] ;
  26. int main (int argc, char *argv[])
  27. { SNDFILE *file ;
  28. SF_INFO sfinfo ;
  29. char *filename ;
  30. int k ;
  31. filename = "test.au" ;
  32. sfinfo.format      = SF_FORMAT_AU | SF_FORMAT_ALAW ;
  33. sfinfo.samplerate  = 44100 ;
  34. sfinfo.samples     = 123456789 ; /* Wrong length. Library should correct this on sf_close. */
  35. sfinfo.channels    = 1 ;
  36. sfinfo.pcmbitwidth = 16 ;
  37. if (! (file = sf_open_write (filename, &sfinfo)))
  38. { printf ("sf_open_write failed with error : ") ;
  39. sf_perror (NULL) ;
  40. exit (1) ;
  41. } ;
  42. /* Generate a file containing all possible 16 bit sample values 
  43. ** and write it to disk as alaw encoded samples. 
  44. */
  45. for (k = 0 ; k < 0x10000 ; k++)
  46. short_buffer [k] = k & 0xFFFF ;
  47. sf_write_short (file, short_buffer, BUFFER_SIZE) ;
  48. sf_close (file) ;
  49. /* Now open that file and compare the alaw encoded sample values 
  50. ** with what they should be.
  51. */
  52. if (! (file = sf_open_read (filename, &sfinfo)))
  53. { printf ("sf_open_write failed with error : ") ;
  54. sf_perror (NULL) ;
  55. exit (1) ;
  56. } ;
  57. if (sf_read_raw (file, alaw_buffer, BUFFER_SIZE) != BUFFER_SIZE)
  58. { printf ("sf_read_raw : ") ;
  59. sf_perror (file) ;
  60. exit (1) ;
  61. } ;
  62. for (k = 0 ; k < 0x10000 ; k++)
  63. if (alaw_encode (short_buffer [k]) != alaw_buffer [k])
  64. { printf ("Encoder error : sample #%d (0x%02X should be 0x%02X)n", k, alaw_buffer [k], alaw_encode (short_buffer [k])) ;
  65. exit (1) ;
  66. } ;
  67. sf_close (file) ;
  68. printf ("    alaw_test : encoder ... okn") ;
  69. /* Now generate a file containing all possible 8 bit encoded
  70. ** sample values and write it to disk as alaw encoded samples. 
  71. */
  72. if (! (file = sf_open_write (filename, &sfinfo)))
  73. { printf ("sf_open_write failed with error : ") ;
  74. sf_perror (NULL) ;
  75. exit (1) ;
  76. } ;
  77. for (k = 0 ; k < 256 ; k++)
  78. alaw_buffer [k] = k & 0xFF ;
  79. sf_write_raw (file, alaw_buffer, 256) ;
  80. sf_close (file) ;
  81. /* Now open that file and compare the alaw decoded sample values 
  82. ** with what they should be.
  83. */
  84. if (! (file = sf_open_read (filename, &sfinfo)))
  85. { printf ("sf_open_write failed with error : ") ;
  86. sf_perror (NULL) ;
  87. exit (1) ;
  88. } ;
  89. if (sf_read_short (file, short_buffer, 256) != 256)
  90. { printf ("sf_read_short : ") ;
  91. sf_perror (file) ;
  92. exit (1) ;
  93. } ;
  94. for (k = 0 ; k < 256 ; k++)
  95. if (short_buffer [k] != alaw_decode (alaw_buffer [k]))
  96. { printf ("Decoder error : sample #%d (0x%02X should be 0x%02X)n", k, short_buffer [k], alaw_decode (alaw_buffer [k])) ;
  97. exit (1) ;
  98. } ;
  99. sf_close (file) ;
  100. printf ("    alaw_test : decoder ... okn") ;
  101. unlink (filename) ;
  102. return 0 ;
  103. } /* main */
  104. /*=================================================================================
  105. ** The following routines came from the sox-12.15 (Sound eXcahcnge) distribution.
  106. **
  107. ** This code is not compiled into libsndfile. It is only used to test the 
  108. ** libsndfile lookup tables for correctness.
  109. **
  110. ** I have included the original authors comments.
  111. */
  112. /*
  113. ** A-law routines by Graeme W. Gill.
  114. ** Date: 93/5/7
  115. **
  116. ** References:
  117. ** 1) CCITT Recommendation G.711
  118. **
  119. */
  120. #define ACLIP 31744
  121. static
  122. unsigned char alaw_encode (int sample)
  123. { static int exp_lut[128] = {1,1,2,2,3,3,3,3,
  124.                                4,4,4,4,4,4,4,4,
  125.                                5,5,5,5,5,5,5,5,
  126.                                5,5,5,5,5,5,5,5,
  127.                                6,6,6,6,6,6,6,6,
  128.                                6,6,6,6,6,6,6,6,
  129.                                6,6,6,6,6,6,6,6,
  130.                                6,6,6,6,6,6,6,6,
  131.                                7,7,7,7,7,7,7,7,
  132.                                7,7,7,7,7,7,7,7,
  133.                                7,7,7,7,7,7,7,7,
  134.                                7,7,7,7,7,7,7,7,
  135.                                7,7,7,7,7,7,7,7,
  136.                                7,7,7,7,7,7,7,7,
  137.                                7,7,7,7,7,7,7,7,
  138.                                7,7,7,7,7,7,7,7};
  139.     int sign, exponent, mantissa;
  140.     unsigned char Alawbyte;
  141.     /* Get the sample into sign-magnitude. */
  142.     sign = ((~sample) >> 8) & 0x80 ; /* set aside the sign */
  143.     if (sign == 0) 
  144. sample = -sample ; /* get magnitude */
  145.     if (sample > ACLIP) 
  146. sample = ACLIP ; /* clip the magnitude */
  147.     /* Convert from 16 bit linear to ulaw. */
  148.     if (sample >= 256)
  149. { exponent = exp_lut [(sample >> 8) & 0x7F] ;
  150. mantissa = ( sample >> ( exponent + 3 ) ) & 0x0F ;
  151. Alawbyte = ((exponent << 4) | mantissa) ;
  152. }
  153.     else
  154. Alawbyte = (sample >> 4) ;
  155.     
  156. Alawbyte ^= (sign ^ 0x55) ;
  157.     return Alawbyte ;
  158. } /* alaw_encode */
  159. static
  160. int alaw_decode (unsigned int Alawbyte)
  161. { static int exp_lut[8] = { 0, 264, 528, 1056, 2112, 4224, 8448, 16896 } ;
  162.     int sign, exponent, mantissa, sample ;
  163.     Alawbyte ^= 0x55 ;
  164.     sign = ( Alawbyte & 0x80 ) ;
  165.     Alawbyte &= 0x7f ; /* get magnitude */
  166.     if (Alawbyte >= 16)
  167. { exponent = (Alawbyte >> 4 ) & 0x07 ;
  168. mantissa = Alawbyte & 0x0F ;
  169. sample = exp_lut[exponent] + (mantissa << ( exponent + 3 )) ;
  170. }
  171.     else
  172. sample = (Alawbyte << 4) + 8 ;
  173.     if (sign == 0) 
  174. sample = -sample ;
  175.     return sample ;
  176. } /* alaw_decode */