International Journal of Modern Physics D,
© World Scientific Publishing Company
Adaptive filtering techniques for interferometric data preparation :
removal of long-term sinusoidal signals and oscillatory transients
E. CHASSANDE-MOTTIN (1) and S. DHURANDHAR (1) & (2)
(1) Albert- Einstein- Institut, Am Muhlenberg, 1
(2) IUCAA, Postbag 4, Ganeshkhind
Pune 411 007, India
Received (received date)
Revised (revised date)
We propose an adaptive denoising scheme for poorly modeled non-Gaussian features
in the gravitational wave interferometric data. Preliminary tests on real data show
Keywords: Gravitational wave interferometric data, detector characterization, transient
denoising, LMS method, adaptive line enhancement
Several large-scale interferometric gravitational wave detectors will come on-line
soon, such as L1GO in the U.S., the French/Italian Virgo project, GEO600 the
German/British interferometer and TAMA in Japan. Gravitational wave detectors
produce an enormous volume of output. Data analysis techniques will have to be
developed to optimally extract the weak signature of a gravitational wave from
these data. Many of the techniques developed so far are based on matched filtering
and assume stationary Gaussian noise.
However, the real data stream from the detectors is not expected to satisfy the
stationary and Gaussian assumptions. This disparity between standard Gaussian
assumptions and real data characteristics poses a major problem to the direct appli-
cation of matched filtering techninues in particular for burst sources such as black
hole binary quasinormal ringings □ or inspiral waveforms i. In fact, the data from
the Caltech 40 meter proto-type interferometer has the expected broadband noise
spectrum, but superposed on this are several other noise features: such as long-
term sinusoidal disturbances coming from suspensions and electric main harmonics
and also ringdown transients occurring occasionally, typically due to servo-controls
instabilities or mechanical relaxation in suspension system etc. While no precise
model can be given for this noise until the detector is completed and fully tested,
matched filtering techniques cannot be used to locate/remove these noisy signals.
We propose a denoising method based on adaptive linear prediction techniques
2 Adaptive filtering techniques for interferometric data
which does not require any precise a priori information about the noise character-
istics. Although our method does not pretend to optimality, we believe that its
simplicity makes it useful for data preparation and for the understanding of the
In the following, we present the structure of the proposed algorithm and some re-
sults obtained with the data from the Caltech 40 meter proto-type interferometer a.
For a more detailed presentation, we refer the reader to □.
2. Adaptive linear prediction
The idea is to predict the current signal sample x k with a collection of past samples
Xk = (xk—d-n, n = 0, 1, . . . , N — 1)*, the delay d > 1 being fixed arbitrarily. This is
possible, only if the target sample shares enough information with (i.e., is sufficiently
correlated to) the previous ones. In other words, the only predictable part of the sig-
nal is the one whose correlation length is sufficiently large (i.e., long-term sinusoids
or ringdowns). On the other hand, the broad band noise cannot be predicted, as it is
not possible to guess the next value in this way. The prediction y k of x k is obtained
through a linear combining of these data samples weighed by the corresponding
coefficients w n , forming the tap-weight vector W — (w n ,n = 0, 1, ...,N— 1)',
therefore leading to y k — W t X k .
The optimal tap-weight vector W* which leads to the smallest prediction error
efc = yk — Xk in the mean square sense can be proved to minimize a convex cost
function. This minimization can be done using an approach similar to the steepest
descent method with the following evolution equation for the tap-weights referred
to as Least Mean Square or LMS algorithm i :
W k+1 =W k + f ie k X kl (1)
where the step gain parameter fi is an adjustable parameter. The tap-weight coef-
ficients are renewed iteratively so that to converge and stabilize in a neighborhood
of W* whose size is defined by /z.
Once the filter has converged, we reject the predicted part of the signal (i.e., y k )
corresponding to the long-term sinusoids and the ringdown noise and we send the
rest of the signal (i.e., e k ) for further analysis for detection.
This LMS based prediction method is referred to as adaptive line enhancer
(ALE) In this context, the term "adaptive" has two different meanings. First,
it means that it will auto-adjust to reach for the best setup for a problem which
is not initially precisely defined. Second, it is also able to follow changes in the
characteristics of the data being processed in case they occur.
Convergence time, frequency tracking ability and frequency resolution are con-
trolled by the three adjustable parameters : the number of tap-weight coefficients
N, the step gain parameter [i and the prediction depth d. One can take advantage
of certain settings of the ALE to select a family of signals intead of another.
Adaptive filtering techniques for interferometric data.
3. The ALE in practice
Structure of the algorithm — We decompose the frequency axis in p disjoint
frequency subbands of the same size. In each subbands, we apply twice the ALE
with different sets of parameters. In the first step, the adjustable parameters are
tuned to best remove long-term sinusoidal components of the noise ;_ whereas in
the second one, the target is the short-time oscillatory transients (see □ for details
about parameter adjustment).
The ALE needs to be applied only in the parts of the signal which appears
non-Gaussian. Some refinements are adjunct for this purpose : the acceptation
(dismissal) of the first cleaning step relies on the detection of a long-term sinusoids
of sufficient amplitude in the data. In the second step, the ALE is applied only
if the filtered output ul deviates from Gaussianity. Details about these additional
vetos can be found in □.
Results on Caltech 40m proto-type data — We have applied the algorithm
to the Caltech 40meter proto-type data taken in October 1994. Figure |l| illustrates
how the algorithm is operating in the fifth frequency subband (from 617 Hz to 771
Hz) among the p = 32 ones being processed. Figure || shows comparisons between
the power spectra and histograms of the signal before and after denoising.
4. Concluding remarks
The originality of the proposed approach lies in the fact that it is possible to treat
oscillatory transients. However, it remains that a comparison of the performances
achieved by this algorithm on noise features of longer duration with other existing
methods need to be done. Finally, although undertaken in □, this method suffers
from the lack of a complete statistical characterization.
Acknowledgments — We would like to thank the LIGO collaboration for provid-
ing us the Caltech 40meter proto-type data.
1. J. D. E. Creighton. Listening for ringing black holes, jr-qc 9712044, 1997.
2. B. Allen et al. Observational limit on gravitational waves from binary neutron star in
the Galaxy. Phys. Rev. Lett, 83(8):1498-1501, 1999.
3. A. Abramovici et al. Improved sensitivity in a gravitational wave interferometer and
implication for LIGO. Phys. Letters A, 218:157-163, 1996.
4. E. Chassande-Mottin and S. Dhurandhar. Adaptive filtering techniques for the removal
of noisy oscillatory transients from interferometric data. In preparation, 2000.
5. B. Widrow and S. D. Stearns. Adaptive Signal Processing. Prentice Hall, Englewoods
6. The Caltech signal has been downloaded and calibrated with the GRASP package
(version 1.9.3, http://www.lsc-group.phys.uwm.edu).
4 Adaptive filtering techniques for interferometric data .
(a) input signal, band #5
(b) output step #1
(c) output step #2
Fig. 1. Illustration of the denoising procedure «i Caltech proto-type data. In the
subband #5 (between 617 Hz and 771 Hz), the signal Q (the data were taken on the October,
14th 1994, frame #2) contains two power line harmonics (at 660 Hz and 720 Hz), as we see on
the spectrogram (a) (this is a time-frequency representation of the signal energy. Dark regions
are associated to large values of the energy density). We apply the ALE a first time to suppress
long-term components (see spectrogram (b)) and a second run (c) eliminates artefacts of shorter
duration (such as fast fluctuations in the harmonic envelope).
(a) power spec, input sig.[300 Hz-1kHz]
(b) histogram, input signal
0.4 0.6 0.8 1
(c) power spec, output signal
0.4 0.6 0.8 1
Fig. 2. Comparison between power spectra and histograms of input/output signals.
The figure depicts power spectra (left column) and histograms (right column) of the Caltech 40
meter signal (top row) in the operating frequency band, between 300 Hz and 1kHz and the same
signal after denoising (bottom row). The histograms are displayed in a graph with special axes
where a Gaussian bell curve should appeared as a straight line.