home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 10 / 10.iso / l / l292 / 1.ddi / FFT2DEMO.FOR < prev    next >
Encoding:
Text File  |  1989-10-10  |  1.0 KB  |  45 lines
  1.       PROGRAM FFTDEMO
  2.       INCLUDE 'STDHDR.FOR'
  3.  
  4.       REAL xr(0: maxr, 0: maxc)
  5.       REAL yi(0: maxr, 0: maxc)
  6.       INTEGER nd, i,j
  7.  
  8.  
  9.       !!!!!!!!!!!!!!!!!!!!!!!!!!!
  10.       !!!!!!!!!!!!!!!!!!!!!!!!!!!
  11.       ! {main program fft}
  12.       nd = 8
  13.       do i = 0, nd - 1
  14.         do j = 0, nd - 1
  15.           xr(i, j) = COS(pi * i / 4.0) + COS(pi * j / 2.0)
  16.           yi(i, j) = 0.0
  17.         end do
  18.       end do
  19.       CALL showdata(xr, nd)
  20.       CALL FFT2DCalc(xr, yi, nd, nd, 0)
  21.       CALL showdata(xr, nd)
  22.       CALL FFT2DCalc(xr, yi, nd, nd, 1)
  23.       CALL showdata(xr, nd)
  24.  
  25.       END
  26.  
  27.       !!!!!!!!!!!!!!!!!!!!!!!!!!!
  28. C     ****************************
  29.       !!!!!!!!!!!!!!!!!!!!!!!!!!!
  30.       SUBROUTINE showdata (xr, nd)
  31.       INCLUDE 'STDHDR.FOR'
  32.       REAL xr(0: maxr, 0: maxc)
  33.       INTEGER nd, i, j
  34.  
  35.       DO i = 0, nd - 1
  36.         DO j = 0, nd - 1
  37.           WRITE(*,10) xr(i, j)
  38.         END DO
  39.         WRITE (*,*)
  40.       END DO
  41. 10    FORMAT (F9.3,5X \)
  42.       END !! SUB ShowData
  43.  
  44.  
  45.