Experimental Log: Roger Bland, December 26-27, 2001
Matched Filtering on Frequency-Shifted Signals
ifremer on RAFOS
I have more or less convinced myself that the optimal filtering for
this signal is carried out as follows:
- Multiply by a replica of the signal, with a time offset tau.
- Multiply by a replica of the signal shifted in phase by
- Square and add preceeding results.
- Vary time offset tau to find maximum.
I also think that almost the same result is obtained by
frequency-shifting the signal from a reference frequency, and averaging
to effect low-pass filtering (down-sampling). This seems mathematically
to simply shift the frequency by an amount equal to the frequency of the
carrier. Then quadrature replicas are constructed by subtracting the
frequency of the carrier, and the detection carried out in the same
I'm going to try this out now. Then think more later about its
December 27. Get to work. I will use data in file s1245b.dat, 90 sec
of 1000-Hz data: npt (long) = 90,000, yd(npt) (double), websig(npt) (float).
(See November 25.)
21:30 I have a nice first result. I multiplied the data of
s1245b.dat by sin(260 Hz * 2 p t) and averaged
in blocks of 128 points, plotted here.
Next two replicas are constructed. Replica q1raf is a swept
frequency, starting at f1 = 259.375 Hz and sweeping over 1.523 Hz in 80
seconds. A constant 260 Hz is subtracted from this frequency. The
quadrature replica, q2raf, is 90 degrees behind q1raf in phase. The two
replicas are shown here.
The detection now proceeds by multiplying the signal by the two
replicas, summing the arrays, and adding in quadrature. Here is the result.
That's all for today.
Friday December 28, 2001. To proceed with the analysis of the
Webb signals, I need to do some mass data analysis:
- Carry out the cross-correlation with the RAFOS quadrature
replicas for all of the beam-formed, detected and down-sampled files.
From each one, create a file of the cross-correlation output.
- Scan the files of cross-correlation output for peaks. See how
effectively the peaks are detected. Make a file of detected peaks.
- Identify triple detections and try to determine the current
speed. (We expect the resolution to be crappy.)