home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1992 #16 / NN_1992_16.iso / spool / comp / dsp / 1829 < prev    next >
Encoding:
Text File  |  1992-07-26  |  2.0 KB  |  47 lines
  1. Newsgroups: comp.dsp
  2. Path: sparky!uunet!destroyer!gatech!psuvax1!atlantis.psu.edu!news
  3. From: mek@guinan.psu.edu (Mark E. Kotanchek)
  4. Subject: Re: Numerical Recipes FFT routine bugs ???
  5. Message-ID: <$=c1Hsr*p8@atlantis.psu.edu>
  6. Sender: news@atlantis.psu.edu (Usenet)
  7. Organization: Penn State Center for Academic Computing
  8. References: <1992Jul24.203717.2527@bert.eecs.uic.edu>
  9. Date: Sun, 26 Jul 92 18:35:36 GMT
  10. Lines: 35
  11.  
  12. In article <1992Jul24.203717.2527@bert.eecs.uic.edu>  
  13. eddins@bert.eecs.uic.edu (Dr. Steven L. Eddins) writes:
  14. > croteau@ccwf.cc.utexas.edu (Ed Croteau) writes:
  15. > 
  16. > >Does anyone have any information on the "FOUR1.C" routine in the
  17. > >Numerical Recipes in C book? I believe that the routine as printed
  18. > >(and on a source disk) does not properly compute the FFT (try a 
  19. > >simple real SIN wave and examine the imaginary).  Please let me know
  20. > >if I am wrong about this but it seems like the imaginary part is
  21. > >frequency reversed or something similar.
  22. > 
  23. > Perhaps your difficulty is in their definition of the DFT (Discrete
  24. > Fourier Transform).  See equations (12.1.7) and (12.1.9), and note
  25. > that the signs in the complex exponential terms are reversed from the
  26. > engineering norm.  Or maybe I shouldn't say "norm"---it's reversed
  27. > from the definition *I* use.  ;-)
  28.  
  29. Bingo! The problem is that mathematicians view the Fourier transform as a  
  30. transform pair and "a tranform is a transform is a transform" while us  
  31. engineers have this perverted idea that a physical meaning should be  
  32. attached to the domain representations.
  33.  
  34. Anyhow, to solve the problem, simply use the inverse FFT whenever you want  
  35. an FFT and vice versa. Unfortunately, Numerical Recipes is not an isolated  
  36. case of this potentially dangerous definition of the FFT--the Extend  
  37. simulation language on the Mac uses the mathematician's definition and (I  
  38. think) MathCAD also is "in error".
  39.  
  40. --
  41. Mark Kotanchek
  42. Guidance & Control Dept - 363 ASB
  43. Applied Research Lab/Penn State
  44. P.O. Box 30
  45. State College, PA 16804
  46.  
  47.