Costas array
In mathematics, a Costas array can be regarded geometrically as a set of n points lying on the squares of a n×n checkerboard, such that each row or column contains only one point, and that all of the n(n − 1)/2 displacement vectors between each pair of dots are distinct. This results in an ideal 'thumbtack' auto-ambiguity function, making the arrays useful in applications such as sonar and radar. Costas arrays can be regarded as two-dimensional cousins of the one-dimensional Golomb ruler construction, and, as well as being of mathematical interest, have similar applications in experimental design and phased array radar engineering.
Costas arrays are named after John P. Costas, who first wrote about them in a 1965 technical report. Independently, Edgar Gilbert also wrote about them in the same year, publishing what is now known as the logarithmic Welch method of constructing Costas arrays.[1]
Numerical representation
A Costas array may be represented numerically as an n×n array of numbers, where each entry is either 1, for a point, or 0, for the absence of a point. When interpreted as binary matrices, these arrays of numbers have the property that, since each row and column has the constraint that it only has one point on it, they are therefore also permutation matrices. Thus, the Costas arrays for any given n are a subset of the permutation matrices of order n.
Arrays are usually described as a series of indices specifying the column for any row. Since it is given that any column has only one point, it is possible to represent an array one-dimensionally. For instance, the following is a valid Costas array of order N = 4:
| 0 | 0 | 0 | 1 |
| 0 | 0 | 1 | 0 |
| 1 | 0 | 0 | 0 |
| 0 | 1 | 0 | 0 |
There are dots at coordinates: (1,2), (2,1), (3,3), (4,4)
Since the x-coordinate increases linearly, we can write this in shorthand as the set of all y-coordinates. The position in the set would then be the x-coordinate. Observe: {2,1,3,4} would describe the aforementioned array. This makes it easy to communicate the arrays for a given order of N.
Known arrays
All Costas arrays of size up to and including 29×29 are known. Costas arrays are known for infinitely many sizes because there exist two methods of generating them, both of which work near primes (of which there are infinitely many) and powers of primes. These are known as the Welch and Lempel-Golomb generation methods, and come from a branch of mathematics known as finite field theory.
The following table describes all possible Costas arrays of size up to six 6×6:
N = 1
{1}
N = 2
{1,2}
{2,1}
N = 3
{1,3,2}
{2,1,3}
{2,3,1}
{3,1,2}
N = 4
{1,2,4,3}
{1,3,4,2}
{1,4,2,3}
{2,1,3,4}
{2,3,1,4}
{2,4,3,1}
{3,1,2,4}
{3,2,4,1}
{3,4,2,1}
{4,1,3,2}
{4,2,1,3}
{4,3,1,2}
N = 5
{1,3,4,2,5}
{1,4,2,3,5}
{1,4,3,5,2}
{1,4,5,3,2}
{1,5,3,2,4}
{1,5,4,2,3}
{2,1,4,5,3}
{2,1,5,3,4}
{2,3,1,5,4}
{2,3,5,1,4}
{2,3,5,4,1}
{2,4,1,5,3}
{2,4,3,1,5}
{2,5,1,3,4}
{2,5,3,4,1}
{2,5,4,1,3}
{3,1,2,5,4}
{3,1,4,5,2}
{3,1,5,2,4}
{3,2,4,5,1}
{3,4,2,1,5}
{3,5,1,4,2}
{3,5,2,1,4}
{3,5,4,1,2}
{4,1,2,5,3}
{4,1,3,2,5}
{4,1,5,3,2}
{4,2,3,5,1}
{4,2,5,1,3}
{4,3,1,2,5}
{4,3,1,5,2}
{4,3,5,1,2}
{4,5,1,3,2}
{4,5,2,1,3}
{5,1,2,4,3}
{5,1,3,4,2}
{5,2,1,3,4}
{5,2,3,1,4}
{5,2,4,3,1}
{5,3,2,4,1}
N = 6
{1,2,5,4,6,3}
{1,2,6,4,3,5}
{1,3,2,5,6,4}
{1,3,2,6,4,5}
{1,3,6,4,5,2}
{1,4,3,5,6,2}
{1,4,5,3,2,6}
{1,4,6,5,2,3}
{1,5,3,4,6,2}
{1,5,3,6,2,4}
{1,5,4,2,3,6}
{1,5,4,6,2,3}
{1,5,6,2,4,3}
{1,5,6,3,2,4}
{1,6,2,4,5,3}
{1,6,3,2,4,5}
{1,6,3,4,2,5}
{1,6,3,5,4,2}
{1,6,4,3,5,2}
{2,3,1,5,4,6}
{2,3,5,4,1,6}
{2,3,6,1,5,4}
{2,4,1,6,5,3}
{2,4,3,1,5,6}
{2,4,3,6,1,5}
{2,4,5,1,6,3}
{2,4,5,3,6,1}
{2,5,1,6,3,4}
{2,5,1,6,4,3}
{2,5,3,4,1,6}
{2,5,3,4,6,1}
{2,5,4,6,3,1}
{2,6,1,4,3,5}
{2,6,4,3,5,1}
{2,6,4,5,1,3}
{2,6,5,3,4,1}
{3,1,2,5,4,6}
{3,1,5,4,6,2}
{3,1,5,6,2,4}
{3,1,6,2,5,4}
{3,1,6,5,2,4}
{3,2,5,1,6,4}
{3,2,5,6,4,1}
{3,2,6,1,4,5}
{3,2,6,4,5,1}
{3,4,1,6,2,5}
{3,4,2,6,5,1}
{3,4,6,1,5,2}
{3,5,1,2,6,4}
{3,5,1,4,2,6}
{3,5,2,1,6,4}
{3,5,4,1,2,6}
{3,5,4,2,6,1}
{3,5,6,1,4,2}
{3,5,6,2,1,4}
{3,6,1,5,4,2}
{3,6,4,5,2,1}
{3,6,5,1,2,4}
{4,1,2,6,5,3}
{4,1,3,2,5,6}
{4,1,6,2,3,5}
{4,2,1,5,6,3}
{4,2,1,6,3,5}
{4,2,3,5,1,6}
{4,2,3,6,5,1}
{4,2,5,6,1,3}
{4,2,6,3,5,1}
{4,2,6,5,1,3}
{4,3,1,6,2,5}
{4,3,5,1,2,6}
{4,3,6,1,5,2}
{4,5,1,3,2,6}
{4,5,1,6,3,2}
{4,5,2,1,3,6}
{4,5,2,6,1,3}
{4,6,1,2,5,3}
{4,6,1,5,2,3}
{4,6,2,1,5,3}
{4,6,2,3,1,5}
{4,6,5,2,3,1}
{5,1,2,4,3,6}
{5,1,3,2,6,4}
{5,1,3,4,2,6}
{5,1,6,3,4,2}
{5,2,3,1,4,6}
{5,2,4,3,1,6}
{5,2,4,3,6,1}
{5,2,6,1,3,4}
{5,2,6,1,4,3}
{5,3,2,4,1,6}
{5,3,2,6,1,4}
{5,3,4,1,6,2}
{5,3,4,6,2,1}
{5,3,6,1,2,4}
{5,4,1,6,2,3}
{5,4,2,3,6,1}
{5,4,6,2,3,1}
{6,1,3,4,2,5}
{6,1,4,2,3,5}
{6,1,4,3,5,2}
{6,1,4,5,3,2}
{6,1,5,3,2,4}
{6,2,1,4,5,3}
{6,2,1,5,3,4}
{6,2,3,1,5,4}
{6,2,3,5,4,1}
{6,2,4,1,5,3}
{6,2,4,3,1,5}
{6,3,1,2,5,4}
{6,3,2,4,5,1}
{6,3,4,2,1,5}
{6,4,1,3,2,5}
{6,4,5,1,3,2}
{6,4,5,2,1,3}
{6,5,1,3,4,2}
{6,5,2,3,1,4}
A full database of the arrays for all sizes which have been exhaustively checked is available It is not known whether Costas arrays exist for all sizes. Currently, the smallest sizes for which no arrays are known are 32×32 and 33×33.
Constructions
Welch
A Welch–Costas array, or just Welch array, is a Costas array generated using the following method, first discovered by Edgar Gilbert in 1965 and rediscovered in 1982 by Lloyd R. Welch.
The Welch–Costas array is constructed by taking a primitive root g of a prime number p and defining the array A by 
if 
, otherwise 0. The result is a Costas array of size p − 1.
Example:
3 is a primitive element modulo 5.
- 31 = 3
- 32 = 9 = 4 (mod 5)
- 33 = 27 = 2 (mod 5)
- 34 = 81 = 1 (mod 5)
Therefore [3 4 2 1] is a Costas permutation. More specifically, this is an exponential Welch array. The transposition of the array is a logarithmic Welch array.
The number of Welch–Costas arrays which exist for a given size depends on the totient function.
Lempel–Golomb
The Lempel–Golomb construction takes α and β to be primitive elements of the finite field GF(q) and similarly defines if 
, otherwise 0. The result is a Costas array of size q − 2. If α + β = 1 then the first row and column may be deleted to form another Costas array of size q − 3: such a pair of primitive elements exists for every prime power q>2.
See also
Notes
- ↑ An independent discovery of Costas arrays, Aaron Sterling, October 9, 2011.
References
- Barker, L.; Drakakis, K.; Rickard, S. (2009), "On the complexity of the verification of the Costas property" (PDF), Proceedings of the IEEE 97 (3): 586–593, doi:10.1109/JPROC.2008.2011947.
- Beard, J.; Russo, J.; Erickson, K.; Monteleone, M.; Wright, M. (2004), "Combinatoric collaboration on Costas arrays and radar applications", IEEE Radar Conference, Philadelphia, Pennsylvania (PDF), pp. 260–265, doi:10.1109/NRC.2004.1316432.
- Costas, J. P. (1965), Medium constraints on sonar design and performance, Class 1 Report R65EMH33, G.E. Corporation
- Costas, J. P. (1984), "A study of a class of detection waveforms having nearly ideal range-Doppler ambiguity properties" (PDF), Proceedings of the IEEE 72 (8): 996–1009, doi:10.1109/PROC.1984.12967.
- Gilbert, E. N. (1965), "Latin squares which contain no repeated digrams", SIAM Review 7: 189–198, doi:10.1137/1007035, MR 0179095.
- Golomb, S. W.; Taylor, H. (1984), "Construction and properties of Costas arrays" (PDF), Proceedings of the IEEE 72 (9): 1143–1163, doi:10.1109/PROC.1984.12994.
- Guy, Richard K. (2004), "Sections C18 and F9", Unsolved Problems in Number Theory (3rd ed.), Springer Verlag, ISBN 0-387-20860-7.
- Moreno, Oscar (1999), "Survey of results on signal patterns for locating one or multiple targets", in Pott, Alexander; Kumar, P. Vijay; Helleseth, Tor; et al., Difference Sets, Sequences and Their Correlation Properties, NATO Advanced Science Institutes Series 542, Kluwer, p. 353, ISBN 0-7923-5958-5.
External links
- MacTech 1999 Programmer's challenge: Costas arrays
- On-Line Encyclopedia of Integer Sequences: