noise2.c
上传用户:szhypcb168
上传日期:2007-01-06
资源大小:2187k
文件大小:2k
源码类别:

语音压缩

开发平台:

Unix_Linux

  1. /**************************************************************************
  2. *
  3. * ROUTINE
  4. *               noise2
  5. *
  6. * FUNCTION
  7. *
  8. *               generates gaussian noise using the polar method
  9. *
  10. * SYNOPSIS
  11. *               noise2(x1,x2)
  12. *   formal 
  13. *
  14. *                       data    I/O
  15. *       name            type    type    function
  16. *       -------------------------------------------------------------------
  17. *       x1              float     o       a sample of noise source
  18. *       x2              float     o       another sample of noise source
  19. *
  20. ***************************************************************************
  21. *       
  22. * USAGE
  23. *
  24. * noise2 generates two samples of a Gaussian noise
  25. * source for each call using the polar method.
  26. *
  27. ***************************************************************************
  28. *
  29. * CALLED FROM
  30. *
  31. * codebook
  32. *
  33. * CALLS
  34. *
  35. * random
  36. *
  37. ***************************************************************************
  38. *
  39. * REFERENCES
  40. *
  41. *       Knuth, The Art of Programming, Volume 2
  42. *
  43. **************************************************************************/
  44. #include <math.h>
  45. noise2(x1,x2)
  46. float *x1, *x2;
  47. {
  48.   float f1, f2, s, f[2];
  49.   int i, j;
  50.   /* f(i) are samples from a uniform distribution
  51. in the range of 0.0 inclusive to 1.0 inclusive.     */
  52.   do 
  53.   {
  54.     for (i=0; i < 2; i++)
  55.     {
  56.       for (j=1; j < 5; j++)
  57.         f[i]=(float)(random2()+32768)/65535.;
  58.     }
  59.     f1=2.*f[0]-1.;
  60.     f2=2.*f[1]-1.;
  61.     s=f1*f1 + f2*f2;
  62.   } while(s > 1.);
  63.   s=sqrt(-2.*log(s)/s);
  64.   *x1=f1*s;
  65.   *x2=f2*s;
  66. }