next previous
Up: Spectral detector for

2. Generalized chi-square test (tex2html_wrap_inline1391) theory

This test is based on the evaluation of the quadratic error between sample mean and ensemble averaged values of Q non-linear functions of the observations:

Let tex2html_wrap_inline1395 be the discrete version of the measured noise s(t). From a set of non-linear Lk-variate functions Fk and from appropriate choice of Lk time-lags tex2html_wrap_inline1405, a set of Q sample mean generalized moments, wk, are computed on N successive samples. The general form of such a sample mean value can be written as:


 equation243
where tex2html_wrap_inline1413.

Let tex2html_wrap_inline1415 denote the Q-dimensional vector formed by the juxtaposition of the sample mean wk:


 equation258
where tex2html_wrap_inline1423 symbolizes the transpose operator.

When the noise s(t) is free of RFI (termed the tex2html_wrap_inline1427 hypothesis throughout this paper), this vector tex2html_wrap_inline1415 converges to a multivariate Gaussian variable with ensemble average mean vector tex2html_wrap_inline1431 and ensemble average covariance matrix tex2html_wrap_inline1433 (Moulines et al. 1993). tex2html_wrap_inline1431 and tex2html_wrap_inline1433 depend on the statistical properties of the noise under the tex2html_wrap_inline1427 hypothesis. In the present case, under the tex2html_wrap_inline1427 hypothesis, the noise is assumed to be Gaussian. Thus, only the second order statistical properties are involved.

The test function tex2html_wrap_inline1443 consists in computing the quadratic error between the measured tex2html_wrap_inline1415 and the expected tex2html_wrap_inline1431. To normalize and take the statistical links between the wk into account, this quadratic error is weighted by the inverse of the covariance matrix tex2html_wrap_inline1433. Namely,


 equation288

Under the tex2html_wrap_inline1427 hypothesis, the test function tex2html_wrap_inline1443 is distributed asymptotically as a central chi-square variable whose degree of freedom is related to vector size Q. When no parameter of the s(t) statistics needs to be estimated, the degree of freedom is exactly Q (Moulines et al. 1993). Knowing the test distribution under the tex2html_wrap_inline1427 hypothesis, it is easy to establish a detection level tex2html_wrap_inline1467 as a function of the desired false alarm probability.


next previous
Up: Spectral detector for

Copyright by the European Southern Observatory (ESO)
web@ed-phys.fr