Distributed source coding

Short description: Problem in information theory and communication

Distributed source coding (DSC) is an important problem in information theory and communication. DSC problems regard the compression of multiple correlated information sources that do not communicate with each other. By modeling the correlation between multiple sources at the decoder side together with channel codes, DSC is able to shift the computational complexity from encoder side to decoder side, therefore provide appropriate frameworks for applications with complexity-constrained sender, such as sensor networks and video/multimedia compression (see distributed video coding). One of the main properties of distributed source coding is that the computational burden in encoders is shifted to the joint decoder.

History

In 1973, David Slepian and Jack Keil Wolf proposed the information theoretical lossless compression bound on distributed compression of two correlated i.i.d. sources X and Y. After that, this bound was extended to cases with more than two sources by Thomas M. Cover in 1975, while the theoretical results in the lossy compression case are presented by Aaron D. Wyner and Jacob Ziv in 1976.

Although the theorems on DSC were proposed on 1970s, it was after about 30 years that attempts were started for practical techniques, based on the idea that DSC is closely related to channel coding proposed in 1974 by Aaron D. Wyner. The asymmetric DSC problem was addressed by S. S. Pradhan and K. Ramchandran in 1999, which focused on statistically dependent binary and Gaussian sources and used scalar and trellis coset constructions to solve the problem. They further extended the work into the symmetric DSC case.

Syndrome decoding technology was first used in distributed source coding by the DISCUS system of SS Pradhan and K Ramachandran (Distributed Source Coding Using Syndromes). They compress binary block data from one source into syndromes and transmit data from the other source uncompressed as side information. This kind of DSC scheme achieves asymmetric compression rates per source and results in asymmetric DSC. This asymmetric DSC scheme can be easily extended to the case of more than two correlated information sources. There are also some DSC schemes that use parity bits rather than syndrome bits.

The correlation between two sources in DSC has been modeled as a virtual channel which is usually referred as a binary symmetric channel.

Starting from DISCUS, DSC has attracted significant research activity and more sophisticated channel coding techniques have been adopted into DSC frameworks, such as Turbo Code, LDPC Code, and so on.

Similar to the previous lossless coding framework based on Slepian–Wolf theorem, efforts have been taken on lossy cases based on the Wyner–Ziv theorem. Theoretical results on quantizer designs was provided by R. Zamir and S. Shamai, while different frameworks have been proposed based on this result, including a nested lattice quantizer and a trellis-coded quantizer.

Moreover, DSC has been used in video compression for applications which require low complexity video encoding, such as sensor networks, multiview video camcorders, and so on.

With deterministic and probabilistic discussions of correlation model of two correlated information sources, DSC schemes with more general compressed rates have been developed. In these non-asymmetric schemes, both of two correlated sources are compressed.

Under a certain deterministic assumption of correlation between information sources, a DSC framework in which any number of information sources can be compressed in a distributed way has been demonstrated by X. Cao and M. Kuijper. This method performs non-asymmetric compression with flexible rates for each source, achieving the same overall compression rate as repeatedly applying asymmetric DSC for more than two sources. Then, by investigating the unique connection between syndromes and complementary codewords of linear codes, they have translated the major steps of DSC joint decoding into syndrome decoding followed by channel encoding via a linear block code and also via its complement code, which theoretically illustrated a method of assembling a DSC joint decoder from linear code encoders and decoders.

Theoretical bounds

The information theoretical lossless compression bound on DSC (the Slepian–Wolf bound) was first purposed by David Slepian and Jack Keil Wolf in terms of entropies of correlated information sources in 1973. They also showed that two isolated sources can compress data as efficiently as if they were communicating with each other. This bound has been extended to the case of more than two correlated sources by Thomas M. Cover in 1975.

Similar results were obtained in 1976 by Aaron D. Wyner and Jacob Ziv with regard to lossy coding of joint Gaussian sources.

Slepian–Wolf bound

Distributed Coding is the coding of two or more dependent sources with separate encoders and joint decoder. Given two statistically dependent i.i.d. finite-alphabet random sequences X and Y, Slepian–Wolf theorem includes theoretical bound for the lossless coding rate for distributed coding of the two sources as below:

[math]\displaystyle{ R_X\geq H(X|Y), \, }[/math]

[math]\displaystyle{ R_Y\geq H(Y|X), \, }[/math]

[math]\displaystyle{ R_X+R_Y\geq H(X,Y). \, }[/math]

If both the encoder and decoder of the two sources are independent, the lowest rate we can achieve for lossless compression is [math]\displaystyle{ H(X) }[/math] and [math]\displaystyle{ H(Y) }[/math] for [math]\displaystyle{ X }[/math] and [math]\displaystyle{ Y }[/math] respectively, where [math]\displaystyle{ H(X) }[/math] and [math]\displaystyle{ H(Y) }[/math] are the entropies of [math]\displaystyle{ X }[/math] and [math]\displaystyle{ Y }[/math]. However, with joint decoding, if vanishing error probability for long sequences is accepted, the Slepian–Wolf theorem shows that much better compression rate can be achieved. As long as the total rate of [math]\displaystyle{ X }[/math] and [math]\displaystyle{ Y }[/math] is larger than their joint entropy [math]\displaystyle{ H(X,Y) }[/math] and none of the sources is encoded with a rate larger than its entropy, distributed coding can achieve arbitrarily small error probability for long sequences.

A special case of distributed coding is compression with decoder side information, where source [math]\displaystyle{ Y }[/math] is available at the decoder side but not accessible at the encoder side. This can be treated as the condition that [math]\displaystyle{ R_Y=H(Y) }[/math] has already been used to encode [math]\displaystyle{ Y }[/math], while we intend to use [math]\displaystyle{ H(X|Y) }[/math] to encode [math]\displaystyle{ X }[/math]. The whole system is operating in an asymmetric way (compression rate for the two sources are asymmetric).

Wyner–Ziv bound

Shortly after Slepian–Wolf theorem on lossless distributed compression was published, the extension to lossy compression with decoder side information was proposed as Wyner–Ziv theorem. Similarly to lossless case, two statistically dependent i.i.d. sources [math]\displaystyle{ X }[/math] and [math]\displaystyle{ Y }[/math] are given, where [math]\displaystyle{ Y }[/math] is available at the decoder side but not accessible at the encoder side. Instead of lossless compression in Slepian–Wolf theorem, Wyner–Ziv theorem looked into the lossy compression case.

The Wyner–Ziv theorem presents the achievable lower bound for the bit rate of [math]\displaystyle{ X }[/math] at given distortion [math]\displaystyle{ D }[/math]. It was found that for Gaussian memoryless sources and mean-squared error distortion, the lower bound for the bit rate of [math]\displaystyle{ X }[/math] remain the same no matter whether side information is available at the encoder or not.

Virtual channel

Deterministic model

Probabilistic model

Asymmetric DSC vs. symmetric DSC

Asymmetric DSC means that, different bitrates are used in coding the input sources, while same bitrate is used in symmetric DSC. Taking a DSC design with two sources for example, in this example [math]\displaystyle{ X }[/math] and [math]\displaystyle{ Y }[/math] are two discrete, memoryless, uniformly distributed sources which generate set of variables [math]\displaystyle{ \mathbf{x} }[/math] and [math]\displaystyle{ \mathbf{y} }[/math] of length 7 bits and the Hamming distance between [math]\displaystyle{ \mathbf{x} }[/math] and [math]\displaystyle{ \mathbf{y} }[/math] is at most one. The Slepian–Wolf bound for them is:

[math]\displaystyle{ R_X+R_Y \geq 10 }[/math]

[math]\displaystyle{ R_X \geq 5 }[/math]

[math]\displaystyle{ R_Y \geq 5 }[/math]

This means, the theoretical bound is [math]\displaystyle{ R_X+R_Y=10 }[/math] and symmetric DSC means 5 bits for each source. Other pairs with [math]\displaystyle{ R_X+R_Y=10 }[/math] are asymmetric cases with different bit rate distributions between [math]\displaystyle{ X }[/math] and [math]\displaystyle{ Y }[/math], where [math]\displaystyle{ R_X=3 }[/math], [math]\displaystyle{ R_Y=7 }[/math] and [math]\displaystyle{ R_Y=3 }[/math], [math]\displaystyle{ R_X=7 }[/math] represent two extreme cases called decoding with side information.

Practical distributed source coding

Slepian–Wolf coding – lossless distributed coding

It was understood that Slepian–Wolf coding is closely related to channel coding in 1974, and after about 30 years, practical DSC started to be implemented by different channel codes. The motivation behind the use of channel codes is from two sources case, the correlation between input sources can be modeled as a virtual channel which has input as source [math]\displaystyle{ X }[/math] and output as source [math]\displaystyle{ Y }[/math]. The DISCUS system proposed by S. S. Pradhan and K. Ramchandran in 1999 implemented DSC with syndrome decoding, which worked for asymmetric case and was further extended to symmetric case.

The basic framework of syndrome based DSC is that, for each source, its input space is partitioned into several cosets according to the particular channel coding method used. Every input of each source gets an output indicating which coset the input belongs to, and the joint decoder can decode all inputs by received coset indices and dependence between sources. The design of channel codes should consider the correlation between input sources.

A group of codes can be used to generate coset partitions, such as trellis codes and lattice codes. Pradhan and Ramchandran designed rules for construction of sub-codes for each source, and presented result of trellis-based coset constructions in DSC, which is based on convolution code and set-partitioning rules as in Trellis modulation, as well as lattice code based DSC. After this, embedded trellis code was proposed for asymmetric coding as an improvement over their results.

After DISCUS system was proposed, more sophisticated channel codes have been adapted to the DSC system, such as Turbo Code, LDPC Code and Iterative Channel Code. The encoders of these codes are usually simple and easy to implement, while the decoders have much higher computational complexity and are able to get good performance by utilizing source statistics. With sophisticated channel codes which have performance approaching the capacity of the correlation channel, corresponding DSC system can approach the Slepian–Wolf bound.

Although most research focused on DSC with two dependent sources, Slepian–Wolf coding has been extended to more than two input sources case, and sub-codes generation methods from one channel code was proposed by V. Stankovic, A. D. Liveris, etc. given particular correlation models.

General theorem of Slepian–Wolf coding with syndromes for two sources

Theorem: Any pair of correlated uniformly distributed sources, [math]\displaystyle{ X, Y \in \left\{0,1\right\}^n }[/math], with [math]\displaystyle{ \mathbf{d_H}(X, Y) \leq t }[/math], can be compressed separately at a rate pair [math]\displaystyle{ (R_1, R_2) }[/math] such that [math]\displaystyle{ R_1, R_2 \geq n-k, R_1+R_2 \geq 2n-k }[/math], where [math]\displaystyle{ R_1 }[/math] and [math]\displaystyle{ R_2 }[/math] are integers, and [math]\displaystyle{ k \leq n-\log(\sum_{i=0}^t{n \choose i}) }[/math]. This can be achieved using an [math]\displaystyle{ (n,k,2t+1) }[/math] binary linear code.

Proof: The Hamming bound for an [math]\displaystyle{ (n,k,2t+1) }[/math] binary linear code is [math]\displaystyle{ k \leq n-\log(\sum_{i=0}^t{n \choose i}) }[/math], and we have Hamming code achieving this bound, therefore we have such a binary linear code [math]\displaystyle{ \mathbf{C} }[/math] with [math]\displaystyle{ k\times n }[/math] generator matrix [math]\displaystyle{ \mathbf{G} }[/math]. Next we will show how to construct syndrome encoding based on this linear code.

Let [math]\displaystyle{ R_1+R_2=2n-k }[/math] and [math]\displaystyle{ \mathbf{G_1} }[/math] be formed by taking first [math]\displaystyle{ (n-R_1) }[/math] rows from [math]\displaystyle{ \mathbf{G} }[/math], while [math]\displaystyle{ \mathbf{G_2} }[/math] is formed using the remaining [math]\displaystyle{ (n-R_2) }[/math] rows of [math]\displaystyle{ \mathbf{G} }[/math]. [math]\displaystyle{ \mathbf{C_1} }[/math] and [math]\displaystyle{ \mathbf{C_2} }[/math] are the subcodes of the Hamming code generated by [math]\displaystyle{ \mathbf{G_1} }[/math] and [math]\displaystyle{ \mathbf{G_2} }[/math] respectively, with [math]\displaystyle{ \mathbf{H_1} }[/math] and [math]\displaystyle{ \mathbf{H_2} }[/math] as their parity check matrices.

For a pair of input [math]\displaystyle{ \mathbf{(x, y)} }[/math], the encoder is given by [math]\displaystyle{ \mathbf{s_1}=\mathbf{H_1}\mathbf{x} }[/math] and [math]\displaystyle{ \mathbf{s_2}=\mathbf{H_2}\mathbf{y} }[/math]. That means, we can represent [math]\displaystyle{ \mathbf{x} }[/math] and [math]\displaystyle{ \mathbf{y} }[/math] as [math]\displaystyle{ \mathbf{x=u_1G_1+c_{s1}} }[/math], [math]\displaystyle{ \mathbf{y=u_2G_2+c_{s2}} }[/math], where [math]\displaystyle{ \mathbf{c_{s1}, c_{s2}} }[/math] are the representatives of the cosets of [math]\displaystyle{ \mathbf{s1, s2} }[/math] with regard to [math]\displaystyle{ \mathbf{C_1, C_2} }[/math] respectively. Since we have [math]\displaystyle{ \mathbf{y=x+e} }[/math] with [math]\displaystyle{ w(\mathbf{e}) \leq t }[/math]. We can get [math]\displaystyle{ \mathbf{x+y=uG+c_{s}=e} }[/math], where [math]\displaystyle{ \mathbf{u=\left[ u_1, u_2\right] } }[/math], [math]\displaystyle{ \mathbf{c_{s}=c_{s1}+c_{s2}} }[/math].

Suppose there are two different input pairs with the same syndromes, that means there are two different strings [math]\displaystyle{ \mathbf{u^1, u^2} \in \left\{ 0,1\right\}^k }[/math], such that [math]\displaystyle{ \mathbf{u^1G+c_{s}=e} }[/math] and [math]\displaystyle{ \mathbf{u^2G+c_{s}=e} }[/math]. Thus we will have [math]\displaystyle{ \mathbf{(u^1-u^2)G=0} }[/math]. Because minimum Hamming weight of the code [math]\displaystyle{ \mathbf{C} }[/math] is [math]\displaystyle{ 2t+1 }[/math], the distance between [math]\displaystyle{ \mathbf{u_1G} }[/math] and [math]\displaystyle{ \mathbf{u_2G} }[/math] is [math]\displaystyle{ \geq 2t+1 }[/math]. On the other hand, according to [math]\displaystyle{ w(\mathbf{e}) \leq t }[/math] together with [math]\displaystyle{ \mathbf{u^1G+c_{s}=e} }[/math] and [math]\displaystyle{ \mathbf{u^2G+c_{s}=e} }[/math], we will have [math]\displaystyle{ d_H(\mathbf{u^1G, c_{s}}) \leq t }[/math] and [math]\displaystyle{ d_H(\mathbf{u^2G, c_{s}}) \leq t }[/math], which contradict with [math]\displaystyle{ d_H(\mathbf{u^1G, u^2G}) \geq 2t+1 }[/math]. Therefore, we cannot have more than one input pairs with the same syndromes.

Therefore, we can successfully compress the two dependent sources with constructed subcodes from an [math]\displaystyle{ (n,k,2t+1) }[/math] binary linear code, with rate pair [math]\displaystyle{ (R_1, R_2) }[/math] such that [math]\displaystyle{ R_1, R_2 \geq n-k, R_1+R_2 \geq 2n-k }[/math], where [math]\displaystyle{ R_1 }[/math] and [math]\displaystyle{ R_2 }[/math] are integers, and [math]\displaystyle{ k \leq n-\log(\sum_{i=0}^t{n \choose i}) }[/math]. Log indicates Log₂.

Slepian–Wolf coding example

Take the same example as in the previous Asymmetric DSC vs. Symmetric DSC part, this part presents the corresponding DSC schemes with coset codes and syndromes including asymmetric case and symmetric case. The Slepian–Wolf bound for DSC design is shown in the previous part.

Asymmetric case

In the case where [math]\displaystyle{ R_X=3 }[/math] and [math]\displaystyle{ R_Y=7 }[/math], the length of an input variable [math]\displaystyle{ \mathbf{y} }[/math] from source [math]\displaystyle{ Y }[/math] is 7 bits, therefore it can be sent lossless with 7 bits independent of any other bits. Based on the knowledge that [math]\displaystyle{ \mathbf{x} }[/math] and [math]\displaystyle{ \mathbf{y} }[/math] have Hamming distance at most one, for input [math]\displaystyle{ \mathbf{x} }[/math] from source [math]\displaystyle{ X }[/math], since the receiver already has [math]\displaystyle{ \mathbf{y} }[/math], the only possible [math]\displaystyle{ \mathbf{x} }[/math] are those with at most 1 distance from [math]\displaystyle{ \mathbf{y} }[/math]. If we model the correlation between two sources as a virtual channel, which has input [math]\displaystyle{ \mathbf{x} }[/math] and output [math]\displaystyle{ \mathbf{y} }[/math], as long as we get [math]\displaystyle{ \mathbf{y} }[/math], all we need to successfully "decode" [math]\displaystyle{ \mathbf{x} }[/math] is "parity bits" with particular error correction ability, taking the difference between [math]\displaystyle{ \mathbf{x} }[/math] and [math]\displaystyle{ \mathbf{y} }[/math] as channel error. We can also model the problem with cosets partition. That is, we want to find a channel code, which is able to partition the space of input [math]\displaystyle{ X }[/math] into several cosets, where each coset has a unique syndrome associated with it. With a given coset and [math]\displaystyle{ \mathbf{y} }[/math], there is only one [math]\displaystyle{ \mathbf{x} }[/math] that is possible to be the input given the correlation between two sources.

In this example, we can use the [math]\displaystyle{ (7,4, 3) }[/math] binary Hamming Code [math]\displaystyle{ \mathbf{C} }[/math], with parity check matrix [math]\displaystyle{ \mathbf{H} }[/math]. For an input [math]\displaystyle{ \mathbf{x} }[/math] from source [math]\displaystyle{ X }[/math], only the syndrome given by [math]\displaystyle{ \mathbf{s}=\mathbf{H}\mathbf{x} }[/math] is transmitted, which is 3 bits. With received [math]\displaystyle{ \mathbf{y} }[/math] and [math]\displaystyle{ \mathbf{s} }[/math], suppose there are two inputs [math]\displaystyle{ \mathbf{x_1} }[/math] and [math]\displaystyle{ \mathbf{x_2} }[/math] with same syndrome [math]\displaystyle{ \mathbf{s} }[/math]. That means [math]\displaystyle{ \mathbf{H}\mathbf{x_1}=\mathbf{H}\mathbf{x_2} }[/math], which is [math]\displaystyle{ \mathbf{H}(\mathbf{x_1}-\mathbf{x_2})=0 }[/math]. Since the minimum Hamming weight of [math]\displaystyle{ (7,4,3) }[/math] Hamming Code is 3, [math]\displaystyle{ d_H(\mathbf{x_1}, \mathbf{x_2})\geq 3 }[/math]. Therefore, the input [math]\displaystyle{ \mathbf{x} }[/math] can be recovered since [math]\displaystyle{ d_H(\mathbf{x}, \mathbf{y})\leq 1 }[/math].

Similarly, the bits distribution with [math]\displaystyle{ R_X=7 }[/math], [math]\displaystyle{ R_Y=3 }[/math] can be achieved by reversing the roles of [math]\displaystyle{ X }[/math] and [math]\displaystyle{ Y }[/math].

Symmetric case

In symmetric case, what we want is equal bitrate for the two sources: 5 bits each with separate encoder and joint decoder. We still use linear codes for this system, as we used for asymmetric case. The basic idea is similar, but in this case, we need to do coset partition for both sources, while for a pair of received syndromes (corresponds to one coset), only one pair of input variables are possible given the correlation between two sources.

Suppose we have a pair of linear code [math]\displaystyle{ \mathbf{C_1} }[/math] and [math]\displaystyle{ \mathbf{C_2} }[/math] and an encoder-decoder pair based on linear codes which can achieve symmetric coding. The encoder output is given by: [math]\displaystyle{ \mathbf{s_1}=\mathbf{H_1}\mathbf{x} }[/math] and [math]\displaystyle{ \mathbf{s_2}=\mathbf{H_2}\mathbf{y} }[/math]. If there exists two pair of valid inputs [math]\displaystyle{ \mathbf{x_1}, \mathbf{y_1} }[/math] and [math]\displaystyle{ \mathbf{x_2}, \mathbf{y_2} }[/math] generating the same syndromes, i.e. [math]\displaystyle{ \mathbf{H_1}\mathbf{x_1} = \mathbf{H_1}\mathbf{x_2} }[/math] and [math]\displaystyle{ \mathbf{H_1}\mathbf{y_1} = \mathbf{H_1}\mathbf{y_2} }[/math], we can get following([math]\displaystyle{ w() }[/math] represents Hamming weight):

[math]\displaystyle{ \mathbf{y_1}=\mathbf{x_1}+\mathbf{e_1} }[/math], where [math]\displaystyle{ w(\mathbf{e_1}) \leq 1 }[/math]

[math]\displaystyle{ \mathbf{y_2}=\mathbf{x_2}+\mathbf{e_2} }[/math], where [math]\displaystyle{ w(\mathbf{e_2}) \leq 1 }[/math]

Thus: [math]\displaystyle{ \mathbf{x_1}+\mathbf{x_2} \in \mathbf{C_1} }[/math]

[math]\displaystyle{ \mathbf{y_1}+\mathbf{y_2}=\mathbf{x_1}+\mathbf{x_2}+\mathbf{e_3} \in \mathbf{C_2} }[/math]

where [math]\displaystyle{ \mathbf{e_3}=\mathbf{e_2}+\mathbf{e_1} }[/math] and [math]\displaystyle{ w(\mathbf{e_3}) \leq 2 }[/math]. That means, as long as we have the minimum distance between the two codes larger than [math]\displaystyle{ 3 }[/math], we can achieve error-free decoding.

The two codes [math]\displaystyle{ \mathbf{C_1} }[/math] and [math]\displaystyle{ \mathbf{C_2} }[/math] can be constructed as subcodes of the [math]\displaystyle{ (7, 4, 3) }[/math] Hamming code and thus has minimum distance of [math]\displaystyle{ 3 }[/math]. Given the generator matrix [math]\displaystyle{ \mathbf{G} }[/math] of the original Hamming code, the generator matrix [math]\displaystyle{ \mathbf{G_1} }[/math] for [math]\displaystyle{ \mathbf{C_1} }[/math] is constructed by taking any two rows from [math]\displaystyle{ \mathbf{G} }[/math], and [math]\displaystyle{ \mathbf{G_2} }[/math] is constructed by the remaining two rows of [math]\displaystyle{ \mathbf{G} }[/math]. The corresponding [math]\displaystyle{ (5\times7) }[/math] parity-check matrix for each sub-code can be generated according to the generator matrix and used to generate syndrome bits.

Wyner–Ziv coding – lossy distributed coding

In general, a Wyner–Ziv coding scheme is obtained by adding a quantizer and a de-quantizer to the Slepian–Wolf coding scheme. Therefore, a Wyner–Ziv coder design could focus on the quantizer and corresponding reconstruction method design. Several quantizer designs have been proposed, such as a nested lattice quantizer, trellis code quantizer and Lloyd quantization method.

Large scale distributed quantization

Unfortunately, the above approaches do not scale (in design or operational complexity requirements) to sensor networks of large sizes, a scenario where distributed compression is most helpful. If there are N sources transmitting at R bits each (with some distributed coding scheme), the number of possible reconstructions scales [math]\displaystyle{ 2^{NR} }[/math]. Even for moderate values of N and R (say N=10, R = 2), prior design schemes become impractical. Recently, an approach, using ideas borrowed from Fusion Coding of Correlated Sources, has been proposed where design and operational complexity are traded against decoder performance. This has allowed distributed quantizer design for network sizes reaching 60 sources, with substantial gains over traditional approaches.

The central idea is the presence of a bit-subset selector which maintains a certain subset of the received (NR bits, in the above example) bits for each source. Let [math]\displaystyle{ \mathcal{B} }[/math] be the set of all subsets of the NR bits i.e.

[math]\displaystyle{ \mathcal{B} = 2^{\{1,...,NR\}} }[/math]

Then, we define the bit-subset selector mapping to be

[math]\displaystyle{ \mathcal{S} : \{1,...,N\} \rightarrow \mathcal{B} }[/math]

Note that each choice of the bit-subset selector imposes a storage requirement (C) that is exponential in the cardinality of the set of chosen bits.

[math]\displaystyle{ C = \sum_{n=1}^N 2^{|\mathcal{S}(n)|} }[/math]

This allows a judicious choice of bits that minimize the distortion, given the constraints on decoder storage. Additional limitations on the set of allowable subsets are still needed. The effective cost function that needs to be minimized is a weighted sum of distortion and decoder storage

[math]\displaystyle{ J = D + \lambda C }[/math]

The system design is performed by iteratively (and incrementally) optimizing the encoders, decoder and bit-subset selector till convergence.

Non-asymmetric DSC for more than two sources

The syndrome approach can still be used for more than two sources. Consider [math]\displaystyle{ a }[/math] binary sources of length-[math]\displaystyle{ n }[/math] [math]\displaystyle{ \mathbf{x}_1,\mathbf{x}_2,\cdots, \mathbf{x}_a \in \{0,1\}^n }[/math]. Let [math]\displaystyle{ \mathbf{H}_1, \mathbf{H}_2, \cdots, \mathbf{H}_s }[/math] be the corresponding coding matrices of sizes [math]\displaystyle{ m_1 \times n, m_2 \times n, \cdots, m_a \times n }[/math]. Then the input binary sources are compressed into [math]\displaystyle{ \mathbf{s}_1 = \mathbf{H}_1 \mathbf{x}_1, \mathbf{s}_2 = \mathbf{H}_2 \mathbf{x}_2, \cdots, \mathbf{s}_a = \mathbf{H}_a \mathbf{x}_a }[/math] of total [math]\displaystyle{ m= m_1 + m_2 + \cdots m_a }[/math] bits. Apparently, two source tuples cannot be recovered at the same time if they share the same syndrome. In other words, if all source tuples of interest have different syndromes, then one can recover them losslessly.

General theoretical result does not seem to exist. However, for a restricted kind of source so-called Hamming source that only has at most one source different from the rest and at most one bit location not all identical, practical lossless DSC is shown to exist in some cases. For the case when there are more than two sources, the number of source tuple in a Hamming source is [math]\displaystyle{ 2^n (a n + 1) }[/math]. Therefore, a packing bound that [math]\displaystyle{ 2^m \ge 2^n (a n + 1) }[/math] obviously has to satisfy. When the packing bound is satisfied with equality, we may call such code to be perfect (an analogous of perfect code in error correcting code).

A simplest set of [math]\displaystyle{ a, n, m }[/math] to satisfy the packing bound with equality is [math]\displaystyle{ a=3, n=5, m=9 }[/math]. However, it turns out that such syndrome code does not exist. The simplest (perfect) syndrome code with more than two sources have [math]\displaystyle{ n = 21 }[/math] and [math]\displaystyle{ m = 27 }[/math]. Let

[math]\displaystyle{ \mathbf{Q}_1 = \begin{pmatrix} 1 \; 0 \; 0 \; 0 \; 0 \; 0 \; 1 \; 0 \; 0 \; 0 \; 0 \; 1 \; 1 \; 1 \; 0 \; 1 \; 1 \; 0 \; 0 \; 0 \; 0 \\ 0 \; 1 \; 0 \; 0 \; 0 \; 0 \; 1 \; 1 \; 0 \; 0 \; 0 \; 0 \; 1 \; 0 \; 0 \; 0 \; 0 \; 0 \; 1 \; 1 \; 1 \\ 0 \; 0 \; 1 \; 0 \; 0 \; 0 \; 0 \; 1 \; 1 \; 0 \; 0 \; 0 \; 0 \; 1 \; 1 \; 1 \; 0 \; 1 \; 0 \; 1 \; 1 \\ 0 \; 0 \; 0 \; 1 \; 0 \; 0 \; 0 \; 0 \; 1 \; 1 \; 0 \; 0 \; 0 \; 1 \; 0 \; 0 \; 1 \; 1 \; 1 \; 1 \; 0 \\ 0 \; 0 \; 0 \; 0 \; 1 \; 0 \; 0 \; 0 \; 0 \; 1 \; 1 \; 0 \; 1 \; 0 \; 1 \; 1 \; 0 \; 1 \; 1 \; 1 \; 1 \\ 0 \; 0 \; 0 \; 0 \; 0 \; 1 \; 0 \; 0 \; 0 \; 0 \; 1 \; 1 \; 0 \; 0 \; 1 \; 0 \; 0 \; 1 \; 1 \; 0 \; 1 \end{pmatrix}, }[/math] [math]\displaystyle{ \mathbf{Q}_2= \begin{pmatrix} 0 \; 0 \; 0 \; 1 \; 0 \; 1 \; 1 \; 0 \; 1 \; 1 \; 1 \; 1 \; 0 \; 1 \; 0 \; 0 \; 0 \; 1 \; 1 \; 1 \; 1 \\ 1 \; 0 \; 0 \; 0 \; 1 \; 0 \; 1 \; 1 \; 0 \; 1 \; 1 \; 1 \; 1 \; 0 \; 1 \; 1 \; 1 \; 1 \; 0 \; 0 \; 0 \\ 0 \; 1 \; 0 \; 0 \; 0 \; 1 \; 1 \; 1 \; 1 \; 0 \; 1 \; 1 \; 1 \; 0 \; 0 \; 0 \; 0 \; 0 \; 1 \; 0 \; 1 \\ 1 \; 0 \; 1 \; 0 \; 0 \; 0 \; 1 \; 1 \; 1 \; 1 \; 0 \; 1 \; 0 \; 1 \; 1 \; 1 \; 0 \; 0 \; 1 \; 1 \; 1 \\ 0 \; 1 \; 0 \; 1 \; 0 \; 0 \; 1 \; 1 \; 1 \; 1 \; 1 \; 0 \; 0 \; 0 \; 1 \; 0 \; 1 \; 1 \; 0 \; 1 \; 1 \\ 0 \; 0 \; 1 \; 0 \; 1 \; 0 \; 0 \; 1 \; 1 \; 1 \; 1 \; 1 \; 1 \; 1 \; 0 \; 1 \; 0 \; 1 \; 1 \; 1 \; 0 \end{pmatrix}, }[/math] [math]\displaystyle{ \mathbf{Q}_3= \begin{pmatrix} 1 \; 0 \; 0 \; 1 \; 0 \; 1 \; 0 \; 0 \; 1 \; 1 \; 1 \; 0 \; 1 \; 0 \; 0 \; 1 \; 1 \; 1 \; 1 \; 1 \; 1 \\ 1 \; 1 \; 0 \; 0 \; 1 \; 0 \; 0 \; 0 \; 0 \; 1 \; 1 \; 1 \; 0 \; 0 \; 1 \; 1 \; 1 \; 1 \; 1 \; 1 \; 1 \\ 0 \; 1 \; 1 \; 0 \; 0 \; 1 \; 1 \; 0 \; 0 \; 0 \; 1 \; 1 \; 1 \; 1 \; 1 \; 1 \; 0 \; 1 \; 1 \; 1 \; 0 \\ 1 \; 0 \; 1 \; 1 \; 0 \; 0 \; 1 \; 1 \; 0 \; 0 \; 0 \; 1 \; 0 \; 0 \; 1 \; 1 \; 1 \; 1 \; 0 \; 0 \; 1 \\ 0 \; 1 \; 0 \; 1 \; 1 \; 0 \; 1 \; 1 \; 1 \; 0 \; 0 \; 0 \; 1 \; 0 \; 0 \; 1 \; 1 \; 0 \; 1 \; 0 \; 0 \\ 0 \; 0 \; 1 \; 0 \; 1 \; 1 \; 0 \; 1 \; 1 \; 1 \; 0 \; 0 \; 1 \; 1 \; 1 \; 1 \; 0 \; 0 \; 0 \; 1 \; 1 \end{pmatrix}, }[/math] [math]\displaystyle{ \mathbf{G} = [ \mathbf{0} | \mathbf{I}_9] }[/math], and [math]\displaystyle{ \mathbf{G}=\begin{pmatrix} \mathbf{G}_1 \\ \mathbf{G}_2 \\ \mathbf{G}_3 \end{pmatrix} }[/math] such that [math]\displaystyle{ \mathbf{G}_1, \mathbf{G}_2, \mathbf{G}_3 }[/math] are any partition of [math]\displaystyle{ \mathbf{G} }[/math].

[math]\displaystyle{ \mathbf{H}_1= \begin{pmatrix} \mathbf{G}_1 \\ \mathbf{Q}_1 \end{pmatrix}, \mathbf{H}_2= \begin{pmatrix} \mathbf{G}_2 \\ \mathbf{Q}_2 \end{pmatrix}, \mathbf{H}_3= \begin{pmatrix} \mathbf{G}_3 \\ \mathbf{Q}_3 \end{pmatrix} }[/math] can compress a Hamming source (i.e., sources that have no more than one bit different will all have different syndromes). For example, for the symmetric case, a possible set of coding matrices are [math]\displaystyle{ \mathbf{H}_1 = \begin{pmatrix} 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 1 \; 0 \; 0 \\ 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 1 \; 0 \\ 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 1 \\ 1 \; 0 \; 0 \; 0 \; 0 \; 0 \; 1 \; 0 \; 0 \; 0 \; 0 \; 1 \; 1 \; 1 \; 0 \; 1 \; 1 \; 0 \; 0 \; 0 \; 0 \\ 0 \; 1 \; 0 \; 0 \; 0 \; 0 \; 1 \; 1 \; 0 \; 0 \; 0 \; 0 \; 1 \; 0 \; 0 \; 0 \; 0 \; 0 \; 1 \; 1 \; 1 \\ 0 \; 0 \; 1 \; 0 \; 0 \; 0 \; 0 \; 1 \; 1 \; 0 \; 0 \; 0 \; 0 \; 1 \; 1 \; 1 \; 0 \; 1 \; 0 \; 1 \; 1 \\ 0 \; 0 \; 0 \; 1 \; 0 \; 0 \; 0 \; 0 \; 1 \; 1 \; 0 \; 0 \; 0 \; 1 \; 0 \; 0 \; 1 \; 1 \; 1 \; 1 \; 0 \\ 0 \; 0 \; 0 \; 0 \; 1 \; 0 \; 0 \; 0 \; 0 \; 1 \; 1 \; 0 \; 1 \; 0 \; 1 \; 1 \; 0 \; 1 \; 1 \; 1 \; 1 \\ 0 \; 0 \; 0 \; 0 \; 0 \; 1 \; 0 \; 0 \; 0 \; 0 \; 1 \; 1 \; 0 \; 0 \; 1 \; 0 \; 0 \; 1 \; 1 \; 0 \; 1 \end{pmatrix}, }[/math] [math]\displaystyle{ \mathbf{H}_2= \begin{pmatrix} 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 1 \; 0 \; 0 \; 0 \; 0 \; 0 \\ 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 1 \; 0 \; 0 \; 0 \; 0 \\ 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 1 \; 0 \; 0 \; 0 \\ 0 \; 0 \; 0 \; 1 \; 0 \; 1 \; 1 \; 0 \; 1 \; 1 \; 1 \; 1 \; 0 \; 1 \; 0 \; 0 \; 0 \; 1 \; 1 \; 1 \; 1 \\ 1 \; 0 \; 0 \; 0 \; 1 \; 0 \; 1 \; 1 \; 0 \; 1 \; 1 \; 1 \; 1 \; 0 \; 1 \; 1 \; 1 \; 1 \; 0 \; 0 \; 0 \\ 0 \; 1 \; 0 \; 0 \; 0 \; 1 \; 1 \; 1 \; 1 \; 0 \; 1 \; 1 \; 1 \; 0 \; 0 \; 0 \; 0 \; 0 \; 1 \; 0 \; 1 \\ 1 \; 0 \; 1 \; 0 \; 0 \; 0 \; 1 \; 1 \; 1 \; 1 \; 0 \; 1 \; 0 \; 1 \; 1 \; 1 \; 0 \; 0 \; 1 \; 1 \; 1 \\ 0 \; 1 \; 0 \; 1 \; 0 \; 0 \; 1 \; 1 \; 1 \; 1 \; 1 \; 0 \; 0 \; 0 \; 1 \; 0 \; 1 \; 1 \; 0 \; 1 \; 1 \\ 0 \; 0 \; 1 \; 0 \; 1 \; 0 \; 0 \; 1 \; 1 \; 1 \; 1 \; 1 \; 1 \; 1 \; 0 \; 1 \; 0 \; 1 \; 1 \; 1 \; 0 \end{pmatrix}, }[/math] [math]\displaystyle{ \mathbf{H}_3= \begin{pmatrix} 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 1 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \\ 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 1 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \\ 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \; 1 \; 0 \; 0 \; 0 \; 0 \; 0 \; 0 \\ 1 \; 0 \; 0 \; 1 \; 0 \; 1 \; 0 \; 0 \; 1 \; 1 \; 1 \; 0 \; 1 \; 0 \; 0 \; 1 \; 1 \; 1 \; 1 \; 1 \; 1 \\ 1 \; 1 \; 0 \; 0 \; 1 \; 0 \; 0 \; 0 \; 0 \; 1 \; 1 \; 1 \; 0 \; 0 \; 1 \; 1 \; 1 \; 1 \; 1 \; 1 \; 1 \\ 0 \; 1 \; 1 \; 0 \; 0 \; 1 \; 1 \; 0 \; 0 \; 0 \; 1 \; 1 \; 1 \; 1 \; 1 \; 1 \; 0 \; 1 \; 1 \; 1 \; 0 \\ 1 \; 0 \; 1 \; 1 \; 0 \; 0 \; 1 \; 1 \; 0 \; 0 \; 0 \; 1 \; 0 \; 0 \; 1 \; 1 \; 1 \; 1 \; 0 \; 0 \; 1 \\ 0 \; 1 \; 0 \; 1 \; 1 \; 0 \; 1 \; 1 \; 1 \; 0 \; 0 \; 0 \; 1 \; 0 \; 0 \; 1 \; 1 \; 0 \; 1 \; 0 \; 0 \\ 0 \; 0 \; 1 \; 0 \; 1 \; 1 \; 0 \; 1 \; 1 \; 1 \; 0 \; 0 \; 1 \; 1 \; 1 \; 1 \; 0 \; 0 \; 0 \; 1 \; 1 \end{pmatrix}. }[/math]

See also

• Linear code
• Syndrome decoding
• Low-density parity-check code
• Turbo Code