Complementary sequences

For complementary sequences in biology, see complementarity (molecular biology).

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_jx_{j+k}.\,$

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

$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) = C\delta(k),\,$

where C is a constant.

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 both sequences are decimated in time by K they remain complementary. More precisely if from a complementary pair (a(k), b(k)) we form a new pair (a(Nk), b(Nk)) with skipped samples discarded then the new sequences are 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 2.

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.

See also

• Pseudorandom binary sequences (also called maximum length sequences or M-sequences)
• Gold sequences
• Kasami sequences
• Polyphase sequence
• Walsh–Hadamard sequences
• Zadoff–Chu sequence
• Binary Golay code (Error-correcting code)
• Ternary Golay code (Error-correcting code)

References

• Golay, Marcel J.E. (1949). "Multislit spectroscopy". J. Opt. Soc. Am. 39: 437–444. doi:10.1364/JOSA.39.000437.
• Golay, Marcel J.E. (April 1961). "Complementary series". IRE Trans. Inform. 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. Combin. 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. 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.