home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 10 / 10.iso / l / l292 / 1.ddi / FFTPOWER.FOR < prev    next >
Encoding:
Text File  |  1989-11-21  |  1.3 KB  |  46 lines
  1.  
  2.       PROGRAM FFTPOWER
  3.       INCLUDE 'stdhdr.for'
  4.  
  5.       REAL xdata(0: maxv), ydata(0: maxv)
  6.       REAL tdata(0: maxv)
  7.       INTEGER  n, i
  8.       LOGICAL BW
  9.  
  10.        n = 2048
  11.        DO i = 0, n - 1
  12.          xdata(i) = 6 * COS(pi * i / (19.0)) + 3.5 *
  13.      +      COS(pi * i / 32.0) + 3.0 * SIN(pi * i / 4.0)
  14.          ydata(i) = 0.01
  15.          tdata(i) = REAL(i) + 1.0
  16.        END DO
  17.        CALL InitSEGraphics(-1, 'C:\FOR\LIB\*.FON')
  18.        CALL BlackAndWhite(BW)
  19.        CALL SetCurrentWindow(3)
  20.        CALL ClearWindow
  21.        IF (.NOT. BW)  CALL SetViewBackground (15)
  22.        CALL SetAxesType(0, 0)
  23.        CALL SelectColor(1)
  24.        CALL AutoAxes(tdata, xdata, n, 0)
  25.        CALL LinePlotData(tdata, xdata, n, 2, 0)
  26.        CALL TitleWindow('Original data')
  27.  
  28.        CALL PowerSpectrumCalc(xdata, ydata, n, .01)
  29.        ydata(0) = 0.01
  30.        CALL SetCurrentWindow(4)
  31.        IF (.NOT. BW)  CALL SetViewBackground (15)
  32.        CALL SelectColor(1)
  33.        CALL SetAxesType(0, 0)
  34.        CALL AutoAxes(ydata, xdata, n / 2, 0)
  35.        CALL LinePlotData(ydata, xdata, n / 2, 2, 1)
  36.        CALL TitleWindow('Power Spectrum')
  37.        CALL TitleXAxis('Frequency')
  38.        CALL TitleYAxis('Power')
  39.        READ (*,*)
  40.       CALL CloseSEGraphics
  41.       END
  42.  
  43.  
  44.  
  45.  
  46.