Wiener filter

wiener1 [-f# -w# -o outfile -l# -x# -c# -V# -h] file

-f sampling rate (e.g. in Hz, default 1.)
-w frequency resolution (e.g. in Hz, default 1/N)>
-l number of values to be read (all)
-x number of values to be skipped (0)
-c column to be read (1 or file,#)
-o output file name, just -o means file_amp
-V verbosity level (0 = only fatal errors)
-h show this message
Now edit file_amp to give desired periodogram

wiener2 [-f# -w# -o outfile1 -O outfile -l# -x# -c# -V# -h] file


-f sampling rate (e.g. in Hz, default 1.)
-w frequency resolution (e.g. in Hz, default 1/N)>
-l number of values to be read (all)
-x number of values to be skipped (0)
-c column to be read (1 or file,#)
-o output file of wiener1, just -o means file_amp
-O final output file name, just -O means file_wc
-V verbosity level (0 = only fatal errors)
-h show this message
Wiener filter. The first call produces the original periodogram on stdout (or in file_amp with the -o option). This may then be edited to provide the desired periodogram. The second call uses file and the output of wiener1 to generate a filtered sequence. The final result is written to stdout, or file_wc with the -O option). The option setting of -o should be the same for both calls. Internally, the series is padded with zeroes in order to get a FFT-able number of points. A warning is issued if applicable. It is recommended to plot the spectral estimator computed by wiener1 in order to adjust the frequency resolution -w properly.

Table of Contents * TISEAN home