Description of the program: mem_spec

This program estimates the power spectrum of a scalar data set on the basis of the maximum entropy principle. A theoretical description of the idea behind is given in the Numerical Recipes.


mem_spec [Options]

Everything not being a valid option will be interpreted as a potential datafile name. Given no datafile at all, means read stdin. Also - means stdin

Possible options are:

Option Description Default
-l# number of data to use whole file
-x# number of lines to be ignored 0
-c# column to be read 1
-p# number of poles 128
-P# number of frequences to print 2000
-f# sampling rate in Hz 1
-o# output file name without file name: 'datafile'.spec
(or stdin.spec if data were read from stdin)
Without -o data is written to stdout
-V# verbosity level
  0: only panic messages
  1: add input/output messages
  2: add the AR coefficients to the output
-h show these options none

Description of the Output:

The first line shows the average forecast error of the AR-model fitted. The following p lines contain the coefficients of the fitted AR-model and the last P lines contain the power spectrum.
View the C-source.
See also spectrum
Table of Contents * TISEAN home