home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 7 / 07.iso / c / c129 / 1.ddi / FFT2DEMO.C < prev    next >
Encoding:
C/C++ Source or Header  |  1988-02-02  |  738 b   |  45 lines
  1.  
  2. # include <math.h>
  3. # include "miscio.h"
  4. # include <stdio.h>
  5. # include "fft.h"
  6.  
  7. float xr[8][8];
  8. float yi[8][8];
  9. int lower;
  10. int upper;
  11. int nd;
  12. int i;
  13. int j;
  14. int k;
  15.  
  16.  
  17. void printdata()
  18. {
  19.  
  20.    printf("\n");
  21.    for ( i = 0; i <= nd - 1; ++i ) {
  22.       for ( j = 0; j <= nd - 1; ++j ) {
  23.          printf("%8.3f  ", xr[i][j]);
  24.       }
  25.       printf("\n");
  26.    }
  27. }
  28.  
  29.  
  30. void main()
  31. {
  32.    nd = 8;
  33.    for ( i = 0; i <= nd - 1; ++i ) {
  34.       for ( j = 0; j <= nd - 1; ++j ) {
  35.          xr[i][j] = cos(M_PI * i / 4.0) + cos(M_PI * j / 2.0);
  36.          yi[i][j] = 0.0;
  37.       }
  38.    }
  39.    printdata();
  40.    FFT2DCalc(&xr[0][0],&yi[0][0],nd,nd,0);
  41.    printdata();
  42.    FFT2DCalc(&xr[0][0],&yi[0][0],nd,nd,1);
  43.    printdata();
  44. }
  45.