home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 7 / 07.iso / c / c129 / 1.ddi / FFTPOWER.C < prev    next >
Encoding:
C/C++ Source or Header  |  1989-04-11  |  1.8 KB  |  77 lines
  1. # include <math.h>
  2. # include <graph.h>
  3. # include "stdlib.h"
  4. # include "stdio.h"
  5. # include "miscio.h"
  6. # include "fft.h"
  7. # include "worlddr.h"
  8. # include "segraph.h"
  9.  
  10.  
  11.  
  12. float xdata[512], ydata[512], tdata[512];
  13. int lower,upper,n,i,j,k;
  14. char c;
  15.  
  16.  
  17.  
  18. void main()   /* main program fft */
  19. {
  20.   n = 128;
  21.   for ( i= 0; i <= n-1; i++ )
  22.   {
  23.     xdata[i] = 6.0 * cos(M_PI*i/(19.0 + frandom())) + 3.5 *cos(M_PI*i/32.0 ) + 3.0 * sin(M_PI*i/4.0) ;
  24.     ydata[i] = 0.0;
  25.  
  26.   tdata[i] = i;
  27.   }
  28.  /****************************************************
  29.     - IF YOU WANT TO CONFIGURE FOR CGA 4 COLORS    
  30.       CHOOSE InitSEGraphics(1);
  31.     - HERCULES SUPPORT IS NOW INCLUDED IN MICROSOFT
  32.       C VERSION 5.1.  IF YOU WANT TO USE THE         
  33.       HERCULES DRIVER YOU MUST FIRST INSTALL       
  34.       MSHERC.COM FROM DOS BY TYPING  
  35.                   msherc<enter>
  36.       THEN IN YOUR MAIN PROGRAM CALL
  37.           InitSEGraphics(12);
  38.     - REFER TO THE README.DOC FILE - PART 3
  39.       INCLUDED IN YOUR MICROSOFT C 5.1 COMPILER 
  40.       FOR DETAILS ON ADDED HERCULES SUPPORT.
  41.  **************************************************/
  42.     InitSEGraphics(6);
  43.  
  44.     SetCurrentWindow(3);
  45.     ClearWindow();
  46.     SetAxesType(0,0);
  47.     BorderCurrentWindow(15);
  48.     SelectColor(1);
  49.     AutoAxes(tdata,xdata,n,0);
  50.     LinePlotData(tdata, xdata, n,2, 1);
  51.     SelectColor(1);
  52.     TitleWindow("Original data");
  53.  
  54.  
  55.  
  56.     PowerSpectrumCalc(xdata,ydata,n,0.01);
  57.     SetCurrentWindow(4);
  58.     ClearWindow();
  59.     BorderCurrentWindow(15);
  60.     SelectColor(1);
  61.     SetAxesType(0,0);
  62.     AutoAxes(ydata,xdata,n / 2,0);
  63.     LinePlotData(ydata, xdata, n / 2,4, 1);
  64.  
  65.  
  66.     TitleWindow("Power Spectrum");
  67.     TitleXAxis("Frequency");
  68.     TitleYAxis("Power Density");
  69.  
  70.     scanf( "%c",&c );
  71.  
  72.     CloseSEGraphics();
  73.  
  74. }
  75.  
  76.  
  77.