In 1916, Albert Einstein predicted the existence of gravitational waves  .
The discovery of the binary pulsar system PSR B1913 + 16 by Hulse and Taylor  and subsequent observations of its energy loss by Taylor and Weisberg  demonstrated the existence of gravitational waves.
Experiments to detect gravitational waves began with Weber  .
It was announced the observation of gravitational waves from merging black holes thanks to the two signals registered according to observatories LIGO Hanford (H1) and the LIGO Livingston (L1) in  .
This paper presents the results of processing of registration data posted on the website LIGO ( https://losc.ligo.org/events/GW150914)  .
2. Registration Data of Two Observatories
The time of arrival of gravitational waves:
−Tdetect = 1,126,259,462.39 = September 14, 2015, 09:50:45.39 UTC, respective signals detected in the frequency band [35 ¸ 350 Hz], the events have a combined signal-to-noise ratio SNRannounced ≈ 24   .
These data were presented in February-May, 2016 on the website  , later, in July 2016 the time of the gravitational wave arrival was corrected and been removed 0.39 seconds, this is roughening of accuracy the time of arrival.
Because of this it is now the time of arrival of gravitational waves is specified without hundredths of a second on the website  :
Tdetect = 1,126,259,462 = September 14 2015, 09:50:45 UTC.
It was announced in  about the delay between the observatories which is ≈ 6.9 ms (arrived first at L1 and 6.9 ms later at H1).
It is strange that removed hundredths of a second, no time of signal arrivals at the L1 and H1, but announced about the delay between the arrival of the GW- wave at different observatories to within a few thousandths of a second.
Time series (fragments) strains , are contained in the files:
The corresponding random processes have a bandwidth [10 ¸ 2048 Hz].
Fragments , have a duration of 32 seconds, strain time series centered at
GPS 1126259462 = September 14 2015, 09:50:45 UTC.
Sampling frequency .
Start Time of the fragment Tstart = 1,126,259,446 = September 14 2015, 09:50:29 UTC, it is 16 seconds before the time of signal arrival Tdetect = 16.39 sec.
Further seconds and its share are relative or .
On the page “Data release for event GW150914” ( https://losc.ligo.org/events/GW150914/) are given:
-results of processing carried out by the members of LIGO Collaboration (Figure 1);
-signals , (Figure 1) according to observation 14.09.2015.
Experimental templates, shown in Figure 1, can be described as a first approximation by the relation:
where ―the signal amplitude;
f0―the carrier frequency;
φ0―the initial phase;
β―the speed of change of frequency;
Fmax, Fmin―minimum and maximum value of the signal frequency;
These signals are called linear frequency modulated (LFM) signals.
Figure 1. Top part. The gravitational-wave event GW150914 observed by the LIGO Hanford (H1, left column panels) and Livingston (L1, right column panels) detectors. Times are shown relative to September 14, 2015 at 09:50:45 UTC. For visualization, all time series are filtered with a 35 - 350 Hz band-pass filter to suppress large fluctuations outside the detectors’ most sensitive frequency band, and band-reject filters to remove the strong instrumental spectral lines  . Declared time of the entry signal (0.39 s) indicated by blue circles  . Numbers and frequencies of wave phases (left): The true signal entry into 0.33 s (green circle). Bottom part: Experimental templates on observatories (Hanford), (Livingston)  . The horizontal axis―number in the data array, on the ordinate axis scaling factor k = 10−21 is sometimes omitted.
3. Invariants of Signals before and after Filtering
1) Linear-frequency modulated signal (LFM-signal), LFM-signal = Chirp signal.
Chirp signal stays the chirp signal after passing through the band filter.
Filtering is defined by the convolution:
where ―impulse response of a filter; ―input process; ―the filtered process; is the sign of the convolution.
2) The integral of the square of the signal remains practically unchanged if the signal in the filter band. This is a consequence of Parseval’s theorem.
where ―the transfer function of the filter, ―a spectrum filter’s impulse response, ―amplitude-frequency characteristic(AFC) of the filter.
3) The spectrum of the signal after filtering remains practically unchanged if the signal is filtered in band, .
4) If the registered signals and and they coincide with the delay and inversion (Figure 1), as stated in   , their spectra ( ) are almost equal after filtering.
Amplitude spectra ( ) are almost equal after filtering in the filter band ( ).
4. The Procedure for the Detection of Known Signals
Detection of gravitational chirp signal is based on the selection of one from two alternative hypotheses:
H0: ―GW-wave chirp signal is absent;
H1: ―GW-wave chirp signal is present.
From the mathematical theory of statistics follows  that the optimal receiver represents matched filtering procedure and then comparing to a threshold.
The matched filter has an impulse response equal to the inverted templates (Figure 1):
Ttemplate―template duration (duration of impulse response).
The matched filter of input records in the band [35 ¸ 350 Hz] with the signals, represented in Figure 1, is written by the relations:
where ―filtered processes after bandpass filtering in the band [35 ¸ 350 Hz].
The procedure for testing the hypothesis of the presence of chirp signals in the records:
1) Produced filtering and and find and .
2) The convolution is produced in accordance with (Equation (7)) and find and .
3) Modules and are compared with thresholds (P1, P2).
The signal/noise ratio is determined by the formula:
The results of the matched filter in accordance with the Equations (7) are shown in Figure 2, the signal/noise ratio SNRreal ≈ 3.
Such a small value of the leads to the inevitable conclusion about doubtfulness of the detected signals, which are stated in   .
For example, a mixtures of signal and noise were formed to demonstrate the matched filtering efficiency:
where ―impulse response of a band filter [35 ¸ 350 Hz];
, ―time series of strains;
―rms signal ;
―rms noise ;
, kSNR―factor depending on SNR;
―time of signal arrival.
Mixtures of signal and noise are calculated for the following parameters: time of model signal arrival .
Matched filtering results with model signals and shown in Figure 2, the lower part.
The signal/noise ratio at the output SNRout is determined by the SNRin and gain at the expense of processing (BSNR):
For formed processes signal/noise ratio at the input in the simulation:
It is clearly distinguished model signals to the 15-th second after matched filtering on Figure 2, when the input signal/noise ratio ≈ 2.323.
Input signal/noise ratio for gravitational signals in Figure 1.
At such values SNRin by the Equation (10)
Figure 2. (a) The results of the matched filter , ; (b) The results matched filtering the mixture of signals and noises , .
that more than doubled when the model signals.
A small rise seen near 16-th second, corresponding SNRreal ≈ 3, as a result of matched filtering (the upper part of Figure 2), which is much smaller than it should be (Equation (11)).
Thus, the hypothesis H1 is not confirmed.
This means: the chirp signals are absent in the original records (Figure 1).
5. Definition of the Signals and Their Parameters
If the chirp signals are absent in registration data (Figure 1), it is necessary to determine the type of signals that are in LIGO14.09.2015 records.
5.1. Evaluation of the Times of Signal Arrivals
For this purpose, the optimum detector is used, which calculates the likelihood ratio:
From (Equation (12)) it follows that detection of the unknown signal is based on the computation of the functional in the first approximation   :
where ―energy process spectra before and after the current time t, ―lower and upper frequency range of the signal.
Looking at Figure 3, you can see:
Figure 3. The results of the optimal detector Hanford (H1)―blue line, Livingston (L1)― red line.
- Confident detection signals according to the observatory L1, H1;
- signal arrived first at L1 and ≤ 10 ms later at H1;
- SNR ≥ 7;
- arrival time at L1 earlier 16.4 sec.
5.2. Definition of Signal Waveforms Using Filtration
This is done using a two-stage causal filtering:
1st step. Butterworth filter, 35 - 350 Hz band, the filter order = 8.
2nd step. Butterworth filter, 60 - 450 Hz band, the filter order = 4 (to notch components at frequencies near 32 Hz and 60 Hz).
According to   , the amplitude of the wave phases №4 - 7 (Figure 1), and their frequency is greater than 80 Hz and below 250 Hz, so after a two-stage filtration ( ), these wave phases have not practically change.
However, these wave phases are not observed in the time interval 16.39 sec - 16.45 sec in Figure 4, and there are pulse signals (amplitude signals from to ).
This apparent contradiction, and thereby violated №1 invariance property. Similarly, we can be sure of violating the invariance property №2.
In more detail the detected signals are represented in Figure 4 (Regime “Magnifier”).
Signals (Figure 4, Regime “Magnifier”) like wavelets “Mexican hat”, “Sombrero”.
They have the property of solitons:
the higher the frequency, the greater the amplitude,
Spectra of these signals are shown in Figure 5.
It follows from Figure 5:
Spectrum of Signal L1 is a higher frequency than spectrum of signal H1 and invariants (№3, 4) are not performed.
This is a clear contradiction to the declared data about detected signals   , and contradicts the invariance properties of the spectra of signals at filtering (Invariants№№3, 4).
It means: !
6. Discussion of Processing Results
Different types of filtration are used in detection of signal with a priori unknown form: bandpass, whitening(rectify), etc.   .
Physically realizable filter whitening is based on the use of optimal filtering Wiener-Kolmogorov theory   and Levinson-Durbin procedures   .
In   abandoned this approbated way of whitening of noises.
(a)Inverted signal L1 = - sL(t)(b)
Figure 4. (a) Recordings of filtered information, the time interval 13.2 sec - 16.8 sec. (b) Recordings of filtered information, the time interval 16.2 sec - 16.5 sec. Recordings of filtered information, the time interval 16.425 sec - 16.445 sec. Regime “Magnifier”. For a visual comparison the L1 data are inverted.
Figure 5. Signal spectra on observatories L1 (red line) and H1 (blue line).
Whitening (rectify) is proposed to conduct in the frequency domain in   , using the equation:
where ―Fourier transform of ;
―transfer function of the whitening filter;
―inverse Fourier transform.
The transfer function of the whitening filter is in   :
This function has a zero phase.
The filter (Equation (14), Equation (15)) is physical be unrealizable, because any physically realizable filter performs a phase shift   .
For example, the Butterworth filter is physically realizable (causal) filter.
If the real (true) signals shown in Figure 1, the effect of change in frequency observed well by 0.33 seconds by which you can visually detect signals. With such a difference between the spectra of processes that can be seen before and after 0.33 seconds in Figure 1, obviously the optimum detector would detect these signals (Equation (13)).
However, the detection of the signal is observed on the 0.395 sec on the results of optimal detector (Figure 2, red line) using the input data and .
Also according to   (February-June 2016) Signal Detection is indicated on the 0.39 seconds.
The difference between the times of the signal arrivals indicates a contradiction with the assumption of trueness of image of waveform in Figure 1.
The difference, which is equal to 60 milliseconds, is the result of violation of the principle of causality in the processing of registration data   .
The reason is that the transfer function (Equation (15)) defines a physically unrealizable (not causal) filter, its impulse response contains non-zero values for t < 0, and the use of such a filter leads to false detection of signals, shifts the time of their arrival and distortion of their forms (Figure 6)   .
1) Chirp signals are absent in data LIGO 14.09.2015.
This is established by means of a matched filter and two-stage filtration, therefore signals of the merger of black holes are absent in these data.
Figure 6. The results of the processing of the same registration data. Left part―causal filtering, the right part―not causal filtering.
2) Signals of type “soliton” are distinguished by means of bandpass filtering (two steps).
Two pulse signals of different polarity were found in these records.
The corresponding wavelet is: “Mexican hat” or “Sombrero”.
3) Optimal detector detects these signals, SNR ≥ 7.
4) Signal arrived first at L1and ≤10 ms later at H1.
5) The signal spectra have a maximum of:
The spectra differ significantly as a result of a violation of invariance of spectra of signals at a filtration (properties №3, 4).
6) It is found that the detected signals at L1 and H1:
7) Physical unrealizability (not causal) filter is used to whitening the registration data in   and it leads to a false detection signals.
8) It is necessary to be based on the principle of causality when processing signals in physics experiments.
P.S. This article describes a repeatable physical experiment.
Initial data are known: parameters of used Butterworth filters, time intervals of registration and also file names.
Data of registration can be downloaded from website:
This experiment can be repeated to check the main conclusions.
 Taylor, J. H. and Weisberg, J.M. (1982) A New Test of General Relativity— Gravitational Radiation and the Binary Pulsar PSR 1913+16. Astrophysical Journal, 253, 908-920.