home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 7 / 07.iso / c / c129 / 1.ddi / GRAPHFFT.C < prev    next >
Encoding:
C/C++ Source or Header  |  1989-04-11  |  1.6 KB  |  72 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. float xdata[512];
  12. float ydata[512];
  13. float tdata[512];
  14. int lower,upper,n,i,j,k;
  15. char c;
  16.  
  17.  
  18.  
  19. void main()
  20. {
  21.   n = 128;
  22.   for (i= 0; i <= n; i++ )
  23.   {
  24.     xdata[i] = 6 * cos(M_PI*i/(19.0 +frandom())) + 2.5 *cos(M_PI*i/32.0 );
  25.     ydata[i] = 0.0;
  26.     tdata[i] = i;
  27.   }
  28.  
  29.  /****************************************************
  30.     - IF YOU WANT TO CONFIGURE FOR CGA 4 COLORS    
  31.       CHOOSE InitSEGraphics(1);
  32.     - HERCULES SUPPORT IS NOW INCLUDED IN MICROSOFT
  33.       C VERSION 5.1.  IF YOU WANT TO USE THE         
  34.       HERCULES DRIVER YOU MUST FIRST INSTALL       
  35.       MSHERC.COM FROM DOS BY TYPING  
  36.                   msherc<enter>
  37.       THEN IN YOUR MAIN PROGRAM CALL
  38.           InitSEGraphics(12);
  39.     - REFER TO THE README.DOC FILE - PART 3
  40.       INCLUDED IN YOUR MICROSOFT C 5.1 COMPILER 
  41.       FOR DETAILS ON ADDED HERCULES SUPPORT.
  42.  **************************************************/
  43.  
  44.     InitSEGraphics(6);
  45.  
  46.     SetCurrentWindow(3);
  47.     ClearWindow();
  48.     SetAxesType(0,0);
  49.     SelectColor(15);
  50.     AutoAxes(tdata,xdata,n,0);
  51.     LinePlotData(tdata, xdata, n,3, 1);
  52.     SelectColor(15);
  53.     TitleWindow("Original data");
  54.  
  55.  
  56.  
  57.     FFTCalc(xdata,ydata,n);
  58.  
  59.  
  60.     SetCurrentWindow(4);
  61.     SelectColor(15);
  62.     SetAxesType(0,0);
  63.     AutoAxes(tdata,xdata,n/2,0);
  64.     LinePlotData(tdata, xdata, n/2,2, 1);
  65.     TitleWindow("FFT Data");
  66.  
  67.    c = getch();
  68.    CloseSEGraphics();
  69. }
  70.  
  71.  
  72.