Complementary sequences

For complementary sequences in biology, see complementarity (molecular biology). For integer sequences with complementary sets of members see Lambek–Moser theorem.

In applied mathematics, complementary sequences (CS) are pairs of sequences with the useful property that their out-of-phase aperiodic autocorrelation coefficients sum to zero. Binary complementary sequences were first introduced by Marcel J. E. Golay in 1949. In 1961–1962 Golay gave several methods for constructing sequences of length 2^N and gave examples of complementary sequences of lengths 10 and 26. In 1974 R. J. Turyn gave a method for constructing sequences of length mn from sequences of lengths m and n which allows the construction of sequences of any length of the form 2^N10^K26^M.

Later the theory of complementary sequences was generalized by other authors to polyphase complementary sequences, multilevel complementary sequences, and arbitrary complex complementary sequences. Complementary sets have also been considered; these can contain more than two sequences.

Definition

Let (a₀, a₁, ..., a_N − 1) and (b₀, b₁, ..., b_N − 1) be a pair of bipolar sequences, meaning that a(k) and b(k) have values +1 or −1. Let the aperiodic autocorrelation function of the sequence x be defined by

R_{x}(k)=\sum _{j=0}^{N-k-1}x_{j}x_{j+k}.\,

Then the pair of sequences a and b is complementary if:

R_{a}(k)+R_{b}(k)=2N,\,

for k = 0, and

R_{a}(k)+R_{b}(k)=0,\,

for k = 1, ..., N − 1.

Or using Kronecker delta we can write:

R_{a}(k)+R_{b}(k)=2N\delta (k),\,

So we can say that the sum of autocorrelation functions of complementary sequences is a delta function, which is an ideal autocorrelation for many applications like radar pulse compression and spread spectrum telecommunications.

Examples

As the simplest example we have sequences of length 2: (+1, +1) and (+1, −1). Their autocorrelation functions are (2, 1) and (2, −1), which add up to (4, 0).
As the next example (sequences of length 4), we have (+1, +1, +1, −1) and (+1, +1, −1, +1). Their autocorrelation functions are (4, 1, 0, −1) and (4, −1, 0, 1), which add up to (8, 0, 0, 0).
One example of length 8 is (+1, +1, +1, −1, +1, +1, −1, +1) and (+1, +1, +1, −1, −1, −1, +1, −1). Their autocorrelation functions are (8, −1, 0, 3, 0, 1, 0, 1) and (8, 1, 0, −3, 0, −1, 0, −1).
An example of length 10 given by Golay is (+1, +1, −1, +1, −1, +1, −1, −1, +1, +1) and (+1, +1, −1, +1, +1, +1, +1, +1, −1, −1). Their autocorrelation functions are (10, −3, 0, −1, 0, 1,−2, −1, 2, 1) and (10, 3, 0, 1, 0, −1, 2, 1, −2, −1).

Properties of complementary pairs of sequences

Complementary sequences have complementary spectra. As the autocorrelation function and the power spectra form a Fourier pair, complementary sequences also have complementary spectra. But as the Fourier transform of a delta function is a constant, we can write

S_{a}+S_{b}=C_{S},

where C_S is a constant.

S_a and S_b are defined as a squared magnitude of the Fourier transform of the sequences. The Fourier transform can be a direct DFT of the sequences, it can be a DFT of zero padded sequences or it can be a continuous Fourier transform of the sequences which is equivalent to the Z transform for

Z = e j ω

.

CS spectra is upper bounded. As S_a and S_b are non-negative values we can write

S_{a}=C_{S}-S_{b}<C_{S},

also

S_{b}<C_{S}.

If either of the sequences of the CS pair is inverted (multiplied by −1) they remain complementary. More generally if any of the sequences is multiplied by e^jφ they remain complementary;
If either of the sequences is reversed they remain complementary;
If either of the sequences is delayed they remain complementary;
If the sequences are interchanged they remain complementary;
If both sequences are multiplied by the same constant (real or complex) they remain complementary;
If alternating bits of both sequences are inverted they remain complementary. In general for arbitrary complex sequences if both sequences are multiplied by e^jπkn/N (where k is a constant and n is the time index) they remain complementary;
A new pair of complementary sequences can be formed as [a b] and [a −b] where [..] denotes concatenation and a and b are a pair of CS;
A new pair of sequences can be formed as {a b} and {a −b} where {..} denotes interleaving of sequences.
A new pair of sequences can be formed as a + b and a − b.

Golay pair

A complementary pair a, b may be encoded as polynomials A(z) = a(0) + a(1)z + ... + a(N − 1)z^N−1 and similarly for B(z). The complementarity property of the sequences is equivalent to the condition

\vert A(z)\vert ^{2}+\vert B(z)\vert ^{2}=2N\,

for all z on the unit circle, that is, |z| = 1. If so, A and B form a Golay pair of polynomials. Examples include the Shapiro polynomials, which give rise to complementary sequences of length a power of two.

Applications of complementary sequences

Multislit spectrometry
Ultrasound measurements
Acoustic measurements
radar pulse compression
Wi-Fi networks,
3G CDMA wireless networks
OFDM communication systems
Train wheel detection systems^[1]^[2]
Non-destructive tests (NDT)
Communications
coded aperture masks are designed using a 2-dimensional generalization of complementary sequences.

References

^ Donato, P.G.; Ureña, J.; Mazo, M.; Alvarez, F. "Train wheel detection without electronic equipment near the rail line". 2004. doi:10.1109/IVS.2004.1336500
^ J.J. Garcia; A. Hernandez; J. Ureña; J.C. Garcia; M. Mazo; J.L. Lazaro; M.C. Perez; F. Alvarez. "Low cost obstacle detection for smart railway infrastructures". 2004.

Golay, Marcel J.E. (1949). "Multislit spectroscopy". J. Opt. Soc. Am. 39 (6): 437–444. doi:10.1364/JOSA.39.000437. PMID 18152021.
Golay, Marcel J.E. (April 1961). "Complementary series". IRE Trans. Inf. Theory. 7 (2): 82–87. doi:10.1109/TIT.1961.1057620.
Golay, Marcel J.E. (1962). "Note on "Complementary series"". Proc. IRE. 50: 84. doi:10.1109/JRPROC.1962.288278.
Turyn, R.J. (1974). "Hadamard matrices, Baumert-Hall units, four-symbol sequences, pulse compression, and surface wave encodings". J. Comb. Theory A. 16 (3): 313–333. doi:10.1016/0097-3165(74)90056-9.
Borwein, Peter (2002). Computational Excursions in Analysis and Number Theory. Springer. pp. 110–9. ISBN 978-0-387-95444-8.
Donato, P.G.; Ureña, J.; Mazo, M.; De Marziani, C.; Ochoa, A. (2006). "Design and signal processing of a magnetic sensor array for train wheel detection". Sensors and Actuators A: Physical. 132 (2): 516–525. Bibcode:2006SeAcA.132..516D. doi:10.1016/j.sna.2006.02.043.

[1] Donato, P.G.; Ureña, J.; Mazo, M.; Alvarez, F. "Train wheel detection without electronic equipment near the rail line". 2004. doi:10.1109/IVS.2004.1336500

[2] J.J. Garcia; A. Hernandez; J. Ureña; J.C. Garcia; M. Mazo; J.L. Lazaro; M.C. Perez; F. Alvarez. "Low cost obstacle detection for smart railway infrastructures". 2004.

[1]

[2]