home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1992 #16 / NN_1992_16.iso / spool / comp / dsp / 1834 < prev    next >
Encoding:
Internet Message Format  |  1992-07-28  |  1.6 KB
  1. Path: sparky!uunet!newsstand.cit.cornell.edu!vax5.cit.cornell.edu!cpyy
  2. From: cpyy@vax5.cit.cornell.edu
  3. Newsgroups: comp.dsp
  4. Subject: Re: Numerical Recipes FFT routine bugs ???
  5. Message-ID: <1992Jul28.093232.14073@vax5.cit.cornell.edu>
  6. Date: 28 Jul 92 09:32:32 EDT
  7. References: <1992Jul24.203717.2527@bert.eecs.uic.edu> <$=c1Hsr*p8@atlantis.psu.edu>
  8. Distribution: comp
  9. Organization: Cornell University
  10. Lines: 27
  11.  
  12. In article <$=c1Hsr*p8@atlantis.psu.edu>,
  13. mek@guinan.psu.edu (Mark E. Kotanchek) writes: 
  14. > In article <1992Jul24.203717.2527@bert.eecs.uic.edu>  
  15. > eddins@bert.eecs.uic.edu (Dr. Steven L. Eddins) writes:
  16. >> croteau@ccwf.cc.utexas.edu (Ed Croteau) writes:
  17. >> 
  18. >> >Does anyone have any information on the "FOUR1.C" routine in the
  19. >> >Numerical Recipes in C book? I believe that the routine as printed
  20. >> >(and on a source disk) does not properly compute the FFT (try a 
  21. >> >simple real SIN wave and examine the imaginary).  Please let me know
  22. >> >if I am wrong about this but it seems like the imaginary part is
  23. >> >frequency reversed or something similar.
  24.  
  25. I also tried to implement the Numerical Recepie's algorithm for
  26. speech processing algorithms.  Strangely, if you do a forward
  27. transform followed by a reverse transform, you don't get back
  28. your original input.  I had to go through the Cooley-Tukey
  29. derivation to produce the correct algorithm.  It seems that
  30. the one in the book is bass-ackwards in such a way that it will
  31. produce a meaningful magnitude result, with a quasi-meaningless
  32. reversed phase result.
  33.  
  34. -- 
  35.  
  36. Jon C. Russo
  37. internet : cpyy@vax5.cit.cornell.edu
  38. phone    : 607.277.3295
  39.