Weighing matrix

Short description: Mathematical weight device

Weighing matrices are so called because of their use in optimally measuring the individual weights of multiple objects.

In mathematics, a weighing matrix of order [math]\displaystyle{ n }[/math] and weight [math]\displaystyle{ w }[/math] is a matrix [math]\displaystyle{ W }[/math] with entries from the set [math]\displaystyle{ \{0, 1, -1\} }[/math] such that:

[math]\displaystyle{ WW^\mathsf{T} = wI_n }[/math]

Where [math]\displaystyle{ W^\mathsf{T} }[/math] is the transpose of [math]\displaystyle{ W }[/math] and [math]\displaystyle{ I_n }[/math] is the identity matrix of order [math]\displaystyle{ n }[/math]. The weight [math]\displaystyle{ w }[/math] is also called the degree of the matrix. For convenience, a weighing matrix of order [math]\displaystyle{ n }[/math] and weight [math]\displaystyle{ w }[/math] is often denoted by [math]\displaystyle{ W(n,w) }[/math].

Weighing matrices are so called because of their use in optimally measuring the individual weights of multiple objects. When the weighing device is a balance scale, the statistical variance of the measurement can be minimized by weighing multiple objects at once, including some objects in the opposite pan of the scale where they subtract from the measurement.

Properties

Some properties are immediate from the definition. If [math]\displaystyle{ W }[/math] is a [math]\displaystyle{ W(n,w) }[/math], then:

• The rows of [math]\displaystyle{ W }[/math] are pairwise orthogonal. Similarly, the columns are pairwise orthogonal.
• Each row and each column of [math]\displaystyle{ W }[/math] has exactly [math]\displaystyle{ w }[/math] non-zero elements.
• [math]\displaystyle{ W^\mathsf{T}W = w I }[/math], since the definition means that [math]\displaystyle{ W^{-1} = w^{-1}W^\mathsf{T} }[/math], where [math]\displaystyle{ W^{-1} }[/math] is the inverse of [math]\displaystyle{ W }[/math].
• [math]\displaystyle{ \det W = \pm w^{n/2} }[/math] where [math]\displaystyle{ \det W }[/math] is the determinant of [math]\displaystyle{ W }[/math].

A weighing matrix is a generalization of Hadamard matrix, which does not allow zero entries. As two special cases, a [math]\displaystyle{ W(n,n) }[/math] is a Hadamard matrix and a [math]\displaystyle{ W(n,n-1) }[/math] is equivalent to a conference matrix.

Applications

Experiment design

Weighing matrices take their name from the problem of measuring the weight of multiple objects. If a measuring device has a statistical variance of [math]\displaystyle{ \sigma^2 }[/math], then measuring the weights of [math]\displaystyle{ N }[/math] objects and subtracting the (equally imprecise) tare weight will result in a final measurement with a variance of [math]\displaystyle{ 2\sigma^2 }[/math]. It is possible to increase the accuracy of the estimated weights by measuring different subsets of the objects, especially when using a balance scale where objects can be put on the opposite measuring pan where they subtract their weight from the measurement.

An order [math]\displaystyle{ n }[/math] matrix [math]\displaystyle{ W }[/math] can be used to represent the placement of [math]\displaystyle{ n }[/math] objects—including the tare weight—in [math]\displaystyle{ n }[/math] trials. Suppose the left pan of the balance scale adds to the measurement and the right pan subtracts from the measurement. Each element of this matrix [math]\displaystyle{ w_{ij} }[/math] will have:

[math]\displaystyle{ w_{ij} = \begin{cases} 0 & \text{if on the }i\text{th trial the }j\text{th object was not measured} \\ 1 & \text{if on the }i\text{th trial the }j\text{th object was placed in the left pan} \\ -1 & \text{if on the }i\text{th trial the }j\text{th object was placed in the right pan } \\ \end{cases} }[/math]

Let [math]\displaystyle{ \mathbf{x} }[/math] be a column vector of the measurements of each of the [math]\displaystyle{ n }[/math] trials, let [math]\displaystyle{ \mathbf{e} }[/math] be the errors to these measurements each independent and identically distributed with variance [math]\displaystyle{ \sigma^2 }[/math], and let [math]\displaystyle{ \mathbf{y} }[/math] be a column vector of the true weights of each of the [math]\displaystyle{ n }[/math] objects. Then we have:

[math]\displaystyle{ \mathbf{x} = W\mathbf{y} + \mathbf{e} }[/math]

Assuming that [math]\displaystyle{ W }[/math] is non-singular, we can use the method of least-squares to calculate an estimate of the true weights:

[math]\displaystyle{ \mathbf{y} = (W^T W)^{-1} W \mathbf{x} }[/math]

The variance of the estimated [math]\displaystyle{ \mathbf{y} }[/math] vector cannot be lower than [math]\displaystyle{ \sigma^2/n }[/math], and will be minimum if and only if [math]\displaystyle{ W }[/math] is a weighing matrix.

Optical measurement

An optical mask (3) based on a weighing matrix is used in the measurement of the spectrum of incoming light (4). Depending on the corresponding element of the matrix, the light is either absorbed, or passed to one of two intensity detectors (1,2).

Weighing matrices appear in the engineering of spectrometers, image scanners, and optical multiplexing systems. The design of these instruments involve an optical mask and two detectors that measure the intensity of light. The mask can either transmit light to the first detector, absorb it, or reflect it toward the second detector. The measurement of the second detector is subtracted from the first, and so these three cases correspond to weighing matrix elements of 1, 0, and −1 respectively. As this is essentially the same measurement problem as in the previous section, the usefulness of weighing matrices also applies.

Examples

Note that when weighing matrices are displayed, the symbol [math]\displaystyle{ - }[/math] is used to represent −1. Here are some examples:

This is a [math]\displaystyle{ W(2,2) }[/math]:

[math]\displaystyle{ \begin{pmatrix}1 & 1 \\ 1 & -\end{pmatrix} }[/math]

This is a [math]\displaystyle{ W(4,3) }[/math]:

[math]\displaystyle{ \begin{pmatrix} 1 & 1 & 1 & 0 \\ 1 & - & 0 & 1 \\ 1 & 0 & - & - \\ 0 & 1 & - & 1 \end{pmatrix} }[/math]

This is a [math]\displaystyle{ W(7,4) }[/math]:

[math]\displaystyle{ \begin{pmatrix} 1 & 1 & 1 & 1 & 0 & 0 & 0 \\ 1 & - & 0 & 0 & 1 & 1 & 0 \\ 1 & 0 & - & 0 & - & 0 & 1 \\ 1 & 0 & 0 & - & 0 & - & - \\ 0 & 1 & - & 0 & 0 & 1 & - \\ 0 & 1 & 0 & - & 1 & 0 & 1 \\ 0 & 0 & 1 & - & - & 1 & 0 \end{pmatrix} }[/math]

Another [math]\displaystyle{ W(7,4) }[/math]:

[math]\displaystyle{ \begin{pmatrix} - & 1 & 1 & 0 & 1 & 0 & 0 \\ 0 & - & 1 & 1 & 0 & 1 & 0 \\ 0 & 0 & - & 1 & 1 & 0 & 1 \\ 1 & 0 & 0 & - & 1 & 1 & 0 \\ 0 & 1 & 0 & 0 & - & 1 & 1 \\ 1 & 0 & 1 & 0 & 0 & - & 1 \\ 1 & 1 & 0 & 1 & 0 & 0 & - \end{pmatrix} }[/math]

Which is circulant, i.e. each row is a cyclic shift of the previous row. Such a matrix is called a [math]\displaystyle{ CW(n,k) }[/math] and is determined by its first row. Circulant weighing matrices are of special interest since their algebraic structure makes them easier for classification. Indeed, we know that a circulant weighing matrix of order [math]\displaystyle{ n }[/math] and weight [math]\displaystyle{ k }[/math] must be of square weight. So, weights [math]\displaystyle{ 1,4,9,16,... }[/math] are permissible and weights [math]\displaystyle{ k \leq 25 }[/math] have been completely classified. Two special (and actually, extreme) cases of circulant weighing matrices are (A) circulant Hadamard matrices which are conjectured not to exist unless their order is less than 5. This conjecture, the circulant Hadamard conjecture first raised by Ryser, is known to be true for many orders but is still open. (B) [math]\displaystyle{ CW(n,k) }[/math] of weight [math]\displaystyle{ k=s^2 }[/math] and minimal order [math]\displaystyle{ n }[/math] exist if [math]\displaystyle{ s }[/math] is a prime power and such a circulant weighing matrix can be obtained by signing the complement of a finite projective plane. Since all [math]\displaystyle{ CW(n,k) }[/math] for [math]\displaystyle{ k \leq 25 }[/math] have been classified, the first open case is [math]\displaystyle{ CW(105,36) }[/math]. The first open case for a general weighing matrix (certainly not a circulant) is [math]\displaystyle{ W(35,25) }[/math].

Equivalence

Two weighing matrices are considered to be equivalent if one can be obtained from the other by a series of permutations and negations of the rows and columns of the matrix. The classification of weighing matrices is complete for cases where [math]\displaystyle{ w \leq 5 }[/math] as well as all cases where [math]\displaystyle{ n \leq 15 }[/math] are also completed. However, very little has been done beyond this with exception to classifying circulant weighing matrices.

Open questions

There are many open questions about weighing matrices. The main question about weighing matrices is their existence: for which values of [math]\displaystyle{ n }[/math] and [math]\displaystyle{ w }[/math] does there exist a [math]\displaystyle{ W(n,w) }[/math]? A great deal about this is unknown. An equally important but often overlooked question about weighing matrices is their enumeration: for a given [math]\displaystyle{ n }[/math] and [math]\displaystyle{ w }[/math], how many [math]\displaystyle{ W(n,w) }[/math]'s are there?

This question has two different meanings. Enumerating up to equivalence and enumerating different matrices with same n,k parameters. Some papers were published on the first question but none were published on the second important question.