# Full text of "Adaptive filtering techniques for interferometric data preparation: removal of long-term sinusoidal signals and oscillatory transients"

## See other formats

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 encouraging results. Keywords: Gravitational wave interferometric data, detector characterization, transient denoising, LMS method, adaptive line enhancement 1. Motivation 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 1 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 first data. 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 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. References 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 Cliffs, 1984. 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 750 '700 650 750 > 700 650 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 frequency [kHz] (c) power spec, output signal 0.4 0.6 0.8 1 frequency [kHz] 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.