Interesting analogy Mal ! may
be we should start a CFRD
group , Campaign For Real Data ? you can
the 'KEY' person : )
G..
----- Original Message -----
Sent: Monday, August 27, 2012 10:22
PM
Subject: LF: WD2XES - Opera detected
and identified by correlation
Recently I have been pondering wether it is
possible to detect QRSS or Opera transmissions by signal correlation against
a known waveform. This should far be more sensitive than the conventional
incoherent decoding process. John's stable and phase-coherent Opera
transmission last night's provided a welcome opportunity to test this
scheme.
Using SndInput from DL4YHF, I recorded a long
IQ file at 2 samples / second, ie. 2 Hz wide centered on 137561 Hz. The
audio was taken straight from the Rubidium-locked receiver, with no noise
blanking inserted. The data was then postprocessed using a MathCad
spreadsheet. Some results can be viewed in
http://df6nm.bplaced.net/opera/xes/
,
along with TA grabber screenshots showing both XGJ and a weak trace
from XES in 21 mHz FFT (testTA.jpg).
First a high resolution spectrogram was
generated, at 1.9 mHz per bin (xes_spectrogram.png). The central carrier
component of the transmission can clearly be seen. Best SNR occured between
3:30 and 4:30 UT (as indicated by marker ticks), when the peak was 9.0 dB
above the noise (xes_spectrum.png). Scaling noise bandwidth from 2.9 mHz to
2.5 kHz (-49.5 dB), and adding 6 dB for 50% duty cycle, we get a peak-power
SNR of -44.5 dB. This corresponds to -48.5 dB on the Opera SNR scale, about
9 dB below the current decoding threshold for Op-32.
The peak appeared about 0.2 Hz off-center
because the 12 kHz samplerate had not been not calibrated. Once the peak
frequency is accurately identified, the received signal can be correlated
against a "prototype" waveform, which contains the Opera sequence for
WD2XES, 16-fold oversampled. The correlation is efficiently implemented as a
multiplication in Fourier space.
The result (xes_correlation_wd2xes.png) shows
four distinct peaks in time domain at 2:15, 2:48, 3:21 and 3:54 UT, which
should correspond to the a-priori unknown start times of John's Opera
sequences. The repetition period was apparently 32.92 minutes. As the DC
component in the reference waveform had not been removed, the peaks are
riding on a pedestal caused by the self-correlation of the carrier
component.
To check the ability to identify an unknown
station, the correlation to a different callsign was also plotted (WD2XGJ
just as an arbitrary example, see xes_correlation_wrongcode). In spite of
the weak cross-correlation peaks, we find that a correct selection from a
list of potential candidates would certainly be feasible.
It would not be too difficult to automate this
process and create an "Opera deep search" software, which should be able to
detect and identify signals reliably down to about 12 dB below the threshold
of the current Opera decoder. This means at same sensitivity we could
go 16 times faster!
To reap this benefit, the following
prerequisits need to be fulfilled:
- As we need a carrier component,
transmit keying has to be phase coherent. Thus simple keying schemes which
interrupt the oscillator or divider would not work.
- as we look at
phase-sensitive integration over the whole sequence rather than a single
symbol duration, the frequency stability has to be much tighter.
- As
there is no bit-wise decoding involved, we will need to supply a list of
potential candidate callsigns (similar to deep-search in K1JT EME
modes).
Thanks again to John and Warren for the
signals!
Best 73,
Markus
(DF6NM in
JN59NJ)
