NAME

dwiener - perform Wiener deconvolution of two signals

SYNOPSIS

SUBROUTINE DWIENER(NPTS, ACOR, XCOR, FILTER, PRDERR, ISW, IERR)

INTEGER NPTS, ISW, IERR

DOUBLE PRECISION ACOR(∗), XCOR(∗), FILTER(∗), PRDERR(∗)

#include <sunperf.h>

void dwiener(int npts, double ∗acor, double ∗xcor, double ∗filter, double ∗prderr, int isw, int ∗ierr) ;

PURPOSE

DWIENER performs Wiener deconvolution of two signals. 
 

PARAMETERS

NPTS - INTEGER
On entry, the number of points in the input correlations. Unchanged on exit.

ACOR - DOUBLE PRECISION array of DIMENSION (NPTS)
On entry, autocorrelation coefficients. Unchanged on exit.

XCOR - DOUBLE PRECISION array of DIMENSION (NPTS)
On entry, cross-correlation coefficients. Unchanged on exit.

FILTER - DOUBLE PRECISION array of DIMENSION (NPTS)
On exit, filter coefficients.  Unchanged on exit.

PRDERR - DOUBLE PRECISION array of DIMENSION (NPTS)
On exit, the prediction error.

ISW - INTEGER
On entry, if ISW .EQ. 0 then perform spiking deconvolution, otherwise perform general deconvolution.  Unchanged on exit.

IERR - INTEGER
On exit, the deconvolution was successful iff IERR .EQ. 0, otherwise there was an error.

SunOS 5.0  —  Last change: 10 Dec 1998