Systems and Methods for Performing Lossless Source Coding

Information

  • Patent Application
  • 20210255836
  • Publication Number
    20210255836
  • Date Filed
    January 21, 2021
    3 years ago
  • Date Published
    August 19, 2021
    3 years ago
Abstract
Systems and methods in accordance with various embodiments of the invention perform lossless source coding. In several embodiments, a nested code structure is utilized to perform Random Access Source Coding (RASC), where the number of active encoders is initially unknown. In several embodiments, the decoder can attempt to decode using a number of Slepian-Wolf decoders corresponding to an estimated number of sources. One embodiment includes multiple source encoders configured to receive a start message and transmit a portion of a codeword selected by encoding data from a source until an end of epoch message is received. A source decoder can transmit at least one start message, and receive codeword portions transmitted by the plurality of source encoders. When a decoding rule is satisfied, the source decoder can decode data from multiple source encoders based upon received codeword portions, and cause the broadcast transmitter to transmit an end of epoch message.
Description
FIELD OF THE INVENTION

The present invention relates generally to lossless source coding and more specifically to finite blocklength lossless source coding in applications involving one or more encoder.


BACKGROUND

The information-theoretic limit in any lossless source coding scenario is the set of code sizes or rates at which a desired level of reconstruction error is achievable. Shannon's theory analyzes this fundamental limit by allowing an arbitrarily long encoding blocklength in order to obtain a vanishing error probability.


Finite-blocklength limits, which are of particular interest in delay-sensitive and computationally-constrained coding environments, allow a non-vanishing error probability and enable study of refined asymptotics of the rates achievable withencoding blocklength n. Due to their non-vanishing error probability, the resulting codes are sometimes called “almost-lossless” source codes. The tem “source coding” can be used to refer to this almost-lossless coding paradigm.


In point-to-point source coding, non-asymptotic bounds and asymptotic expansions of the minimum achievable rate have been published. For example, a third-order characterization is known for an optimal code that can provide the minimum achievable rate R*(n, ϵ) at blocklength n and error probability ϵ. For a finite-alphabet, stationary, memoryless source with single-letter distribution PX, entropy H(X), and varentropy V(X)>0,












R
*



(

n
,
ϵ

)





H


(
X
)


+




V


(
X
)


n





Q

-
1




(
ϵ
)



-


log





n


2

n




,




(
1
)







where Q−1(⋅) is the inverse complementary Gaussian distribution function, and any higher-order term is bounded by







O


(

1
n

)


.




For a multiple access source coding (MASC) also known as a Slepian-Wolf (SW) source coding the fundamental limit is the set of achievable rate tuples known is the rate region. The first-order rate region for stationary, memoryless and general sources are known and Second-order asymptotic expansions of the MASC rate region for stationary, memoryless sources have also been published. However, optimal codes that can approach the fundamental limits identified by Shannon are only known for a small number of special communication scenarios.


SUMMARY OF THE INVENTION

Systems and methods in accordance with various embodiments of the invention perform lossless source coding. A variety of systems and methods are described including systems and methods that are capable of performing source coding in one or more of three scenarios: point-to-point communication, multiple access communication, and random access communication. It can be shown that for point-to-point coding on stationary, memoryless sources, random code design with maximum likelihood decoding can achieve the same coding rate up to the third-order as an optimal code.


SysteMS and methods in accordance with many embodiments of the invention axe capable of performing Multiple Access Source Coding (MASC) utilizing a decoder such as (but not limited to) a decoder formed by a number of stand alone Slepian Wolf decoders to decode k users by choosing the jointly most probable source realizations consistent with the received codewords. Random access source coding (RASC) systems are also introduced and various embodiments are described in which RASC is performed for all possible encoder activity patterns. In several embodiments, a nested structure of the code is utilized to perform RASC, and there is no need for the encoders to know the set of active encoders a priori. In several embodiments, the decoder can attempt to decode using a number of Slepian Wolf decoders corresponding to an estimated number of sources. In the event that the decoder is unable to successfully decode the received information, then the decoder can retry using a different estimate of the number of sources. The third-order-optimal MASC performance can be achievable even when the only information the encoders receive is the acknowledgment that tells them when to stop transmitting.


One embodiment of the random access source coding (RASC) system includes a plurality of source encoders that are each configured to receive data from one of a plurality of sources. When one of the plurality of source encoders receives data from one of the plurality of sources, then the source encoder is configured to: receive a start of epoch message; and transmit a portion of a codeword selected by encoding data from one of the plurality of sources until an end of epoch feedback message is received. In addition, the RASC system includes a receiver comprising a source decoder, where the source decoder is configured to: cause a broadcast transmitter to transmit at least one start of epoch message; receive at least one portion of a codeword transmitted by at least one of the plurality of source encoders; and when a decoding rule is satisfied, decode data from the at least one of the plurality of the source encoders based upon the received at least one portion of a codeword, and cause the broadcast transmitter to transmit an end of epoch message.


In a further elribodiunent, the source decoder is configured to deteimine whether the decoding rule is satisfied at each of a predetermined set of decode times and transmit a feedback message; each of the plurality of source encoders is configured to receive feedback messages at the predetermined set of potential decoding times; and each of the predetermined set of decode times corresponds to one of a plurality of blocklengths of the codeword selected by encoding data from one of the plurality of sources.


In another embodiment, a first of the plurality of blocklengths of the codeword selected by encoding data from one of the plurality of sources forms a first blocklength code; a second of the plurality of blocklengths of the codeword selected by encoding data from one of the plurality of sources forms a second blocklength code; and the first blocklength code forms a prefix for the second blocklength code.


In a still further embodiment, a third of the plurality of blocklengths of the codeword selected by encoding data from one of the plurality of sources forms a third blocklength code; and the second blocklength code forms a prefix for the third blocklength code.


In still another embodiment, each source encoder is configured to select codewords from a codebook.


In a yet another embodiment, the data received from the plurality of sources is dependent.


In yet another embodiment, the broadcast transmitter sends a negative acknowledgement feedback message when the decoding rule is lot satisfied at a predetermined decode time.


In a further embodiment again, each source encoder that is transmitting is configured to transmit an additional portion of the codeword selected by encoding data from one of the plurality of sources in response to receipt of a negative acknowledgement feedback message.


A source encoding system in accordance with another embodiment again includes: a receiver configured to receive a start of epoch message; and a source encoder configured to transmit a portion of a codeword selected by encoding data from a source until an end of epoch feedback message is received.


In a further additional embodiment, each of the plurality of source encoders is configured to receive feedback messages at each of a predetermined set of decode times; and each of the predetermined set of decode times corresponds to one of a plurality of blocklengths of the codeword selected by encoding data from the source.


In another additional embodiment, a first of the plurality of blocklengths of the codeword selected by encoding data from the source forms a first blocklength code; a second of the plurality of blocklengths of the codeword selected by encoding data from the source forms a second blocklength code; and the first blocklength code forms a prefix for the second blocklength code.


In a still yet further embodiment, a third of the plurality of blocklengths of the codeword selected by encoding data from the source forms a third blocklength code; and the second blocklength code forms a prefix for the third blocklength code.


In still yet another embodiment, the source encoder is configured to transmit an additional portion of the codeword selected by encoding data from the source in response to receipt of a negative acknowledgement feedback message.


In a still further embodiment again, the data received from the plurality of sources is dependent.


A receiver system in accordance with still another embodiment again, includes: a source decoder configured using a codebook; a transmitter configured to transmit at least one start of epoch message; and a receiver configured to receive at least one portion of a codeword. In addition, the source decoder is further configured to determine when a decoding rule is satisfied, and when the decoding rule is satisfied: decode at least one message based upon the received at least one portion of a codeword using the codebook; and cause the transmitter to transmit an end of epoch message.


In a still further additional embodiment, the source decoder is configured to determine whether the decoding rule is satisfied at each of a predetermined set of decode times and transmit a feedback message; and each of the predetermined set of decode times corresponds to one of a plurality of blocklengths of the codeword.


In a yet further embodiment again, a first of the plurality of blocklengths of the codeword forms a first blocklength code; a second of the plurality of blocklengths of the codeword forms a second blocklength code; and the first blocklength code forms a prefix for the second blocklength code.


In yet another embodiment again, a third of the plurality of blocklengths of the codeword forms a third blocklength code; and the second blocklength code forms a prefix for the third blocklength code.


In a yet further additional embodiment, the data received from the plurality of sources is dependent.


In yet another additional embodiment, the broadcast transmitter sends a negative acknowledgement feedback message when the decoding rule is not satisfied at a predetermined decode time.





BRIEF DESCRIPTION OF THE DRAWINGS


FIG. 1 is a schematic diagram illustrating a point-to-point coding system in accordance with an embodiment of the invention.



FIGS. 2A and 2B are plots that illustrate the rate-blocklength trade-offs at ϵ=10−1 (FIG. 2A) and ϵ=10−3 (FIG. 2B) for a pair of binary, stationary, memoryless sources.



FIG. 3 is a plot showing the third-order rate region for systems and methods that perform Multiple Access Source Coding (MASC) in accordance with various embodiments of the invention.



FIGS. 4A and 4B are schematic plots of the (n, ϵ)-rate region and the third-order-optimal sum rate when the pair of random variables X1, X2 are dependent (FIG. 4A), and when X1, X2 are independent (FIG. 4B).



FIG. 5A illustrates a MASC system incorporating feedback and FIG. 5B illustrates a MASC system that utilizes a cooperation facilitator.



FIG. 6 is a flow chart illustrating a process that can be utilized by an encoder to perform source coding within a Random Access Source Coding (RASC) system in accordance with an embodiment of the invention.



FIG. 7 conceptually illustrates an RASC system in accordance with an embodiment of the invention.



FIG. 8 conceptually illustrates a source coding scheme employed by encoders during a single time epoch in an RASC system in accordance with an embodiment of the invention.



FIG. 9 shows the relationship between decoding blocklength, code symbol alphabet sizes, and source coding rate vector.



FIG. 10 shows a number of sensors reading measurements of a common hidden state as an example of an application of RASC in accordance with an embodiment of the invention.



FIGS. 11A and 11B are schematic illustrations of the MASC rate regions for a less redundant source.





DETAILED DESCRIPTION

Turning now to the drawings, systems and methods for performing source coding using random codes in accordance with various embodiments of the invention are illustrated. In several embodiments, codes are utilized that can be characterized in that they can achieve performance similar to an optimal code in source coding applications involving point-to-point communication between an encoder and a decoder. In a number of embodiments, linear random codes are employed and achieve performance similar to an optimal code in source coding applications involving point-to-point communications. In several embodiments, any of a number of different types of decoders can be utilized to perform decoding including (but not limited to) maximum likelihood decoders and/or threshold decoders.


In point-to-point source coding, the encoder maps a discrete random variable X defined on finite or countably infinite alphabet X into a message from codebook [M]. The decoder reconstructs X from the compressed description. A random code is one in which each code word in the codebook [M] has a finite length and the association between the value of X and a particular code word from codebook [M] is determined randomly. In this way, the encoder need not possess any information regarding the distribution of X to encode the source. Random codes are typically sub-optimal, but random codes utilized in accordance with many embodiments of the invention can be shown to have near-optimal performance when utilized for point-to-point source coding. Furthermore, the set of random codes that have near-optimal performance includes codes that can be utilized by a linear compressor. Therefore, systems and methods in accordance with a number of embodiments utilize one of a set of linear random codes that are capable of achieving near-optimal performance. In certain embodiments, random code are utilized in combination with a maximum likelihood decoder. In a number of embodiments, random code are utilized in combination with a threshold decoder. As can readily be appreciated, the use of a threshold decoder will typically involve a requirement to operate at a higher rate than a system that relies upon use of a maximum likelihood decoder.


Systems and methods in accordance with certain embodiments of the invention perform multiple access source coding (MASC) using random codes. In several embodiments, the MASC process employs limited feedback and/or cooperation feedback. In several embodiments, MASC involves uses of a maximum likelihood decoder that chooses the jointly most probable source realizations consistent with received codewords. Where rate points converge to a non-corner point on the asymptotic sum-rate boundary, systems and methods in accordance with embodiments of the invention can be shown to not compromise the performance in lossless data compression up to the third-order term despite each encoder performing separate encoding. For independent sources, it is shown that there are no non-corner points, and MASC separate encoding incurs a positive penalty in the second-order term relative to joint encoding with a point-to-point code. When two sources have the same marginals, this penalty equals the penalty for using two independent blocklength-n codes rather than a single blocklength-2n point-to-point code for encoding 2n samples.


While a MASC assumes a fixed, known collection of encoders, the set of transmitters communicating with a given access point in applications like sensor networks, the internet of things, and random access communication may be unknown or time-varying. To address these scenarios, a new applicationof source coding is introduced which can be referred to as random access source coding (RASC). RASC extends MASC to scenarios where some encoders are inactive, and the decoder seeks to reliably reconstruct the sources associated with the active encoders assuming that the set of active encoders is unknown a priori.


Systems and methods in accordance with several embodiments of the invention perform RASC using rateless encoders that transmit codewords symbol by symbol until a receiver sends a message telling them to stop. Unlike typical rateless codes, which allow arbitrary decoding times the codes utilized in accordance with many embodiments of the invention possess a small set of decoding times. A single-bit feedback from the decoder to all encoders at each potential decoding time can provide, the encoders with feedback concerning whether or not to continue transmitting.


In a number of embodiments, a single random code is utilized to perform RASC that is characterized in that it is a single code that simultaneously achieves, for every possible set of active encoders, the third-order-optimal MASC performance for the active source set. Since traditional random coding arguments do not guarantee the existence of a single deterministic code that meets multiple independent constraints, prior code designs for multiple-constraint scenarios typically employ a family of codes indexed using common randomness shared by all communicators. Systems and methods in accordance with many embodiments of the invention employ an alternative approach that utilizes a single deterministic code. As can readily be appreciated, for stationary, memoryless, permutation-invariant sources, employing identical encoders at all transmitters can reduce RASC design complexity.


While mach of the discussion that follows relates to source coding, the hypothesis testing techniques described herein can be useful in many applications including (but not limited to) testing components (e.g., finding positive cases in pooled coronavirus tests, defective components in a factory). In the discussion below, hypothesis testing is used as a tool for analyzing the performance of various codes, but decoders in accordance with many embodiments of the invention can apply a type of composite hypothesis test that compares the possible decoder outputs as hypotheses of the form “the encoded message is ‘a,’” against the composite hypothesis “the encoded message is one of ‘b’, ‘c’, . . . , ‘z’.” Accordingly, the hypothesis testing techniques described herein can be utilized in any of a variety of applications including (but not limited to) the applications specifically described herein. Systems and methods for performing source coding in accordance with various embodiments of the invention are discussed further below.


Point-To-Point Random Source Coding

A compression system 100 that employs point-to-point random source coding in accordance with an embodiment of the invention is illustrated in FIG. 1. An (M, ϵ) code for a random variable X with discrete alphabet X can be implemented by an encoder 102 employing an encoding function f: X→[M] and a decoder 104 employing a decoding function g: [M]→X with error probability custom-character[g(f(X))≠X]≤ϵ.


In the illustrated embodiment, an (n, M, ϵ) block point-to-point source code is utilized, which is an (M, ϵ) code defined for a random vector Xn with discrete vector alphabet X2. The minimum code size M*(n, ϵ) and rate R*(n, ϵ) achievable at blocklength n and error probability ϵ for the code are











M
*



(

n
,
ϵ

)




=
Δ



min


{

M
:




(

n
,
M
,
ϵ

)


code



}






(
2
)








R
*



(

n
,
ϵ

)




=
Δ




1
n


log








M
*



(

n
,
ϵ

)


.






(
3
)







Shannon's source coding theorem describes the fundamental limit on the asymptotic performance for lossless source coding on a stationary, memoryless source, giving












lim

n
->






R
*



(

n
,
ϵ

)



=

H


(
X
)



,



ϵ



(

0
,
1

)

.







(
4
)







Optimal codes that can approach the fundamental limits identified by Shannon are only known for a small number of special communication scenarios. Consequently, systems and methods in accordance with many embodiments of the invention utilize random codes that have been determined to achieve comparable performance to optimal codes in point-to-point source coding applications. Rates that can be achieved in point-to-point source coding using random codes in accordance with various embodiments of the invention are discussed further below.


Rates Achievable in Point-To-Point Source Coding Using Random Codes

Achievability results that are based on Shannon's random coding argument are important because use of a random code does not require knowledge of the optimal code, which is available only in a few special communication scenarios. The following random coding achievability bound can be obtained by assigning source realizations to codewords independently and uniformly at random. The threshold decoder decodes to x∈X if and only if x is a unique source realization that (i) is compatible with the observed codeword under the given (random) code design, and (ii) has information i(x) below log M−γ.


It can be shown that there exists an (M, ϵ) code for discrete random variable X such that





ϵ≤custom-character[custom-character(X)>log M−γ]+exp(−γ), ∀γ>0.  (5)


Particularizing (5) to a stationary, memoryless source with single-letter distribution PX satisfying V(X)>0 and T(X)<∞, choosing log M and γ optimally, and applying the Berry-Esseen inequality gives











R
*



(

n
,
ϵ

)





H


(
X
)


+




V


(
X
)


n





Q

-
1




(
ϵ
)



+


log





n


2

n


+

O







(

1
n

)

.







(
6
)







Since (6) exceeds the rate bounds for optimal codes by






+


log





n

n





in the third-order term, a question arises as to whether random code design, threshold decoding, or both yield third-order performance penalties. Accordingly, a new random coding bound is derived using a maximum likelihood decoder; this result demonstrates that random coding suffices to achieve the third-order optimal performance for a stationary, memoryless source.


Random code design can be used to derive two new non-asymptotic achievability bounds for point-to-point source coding. These results can be referred to as the dependence testing (DT) bound and the random coding union (RCU) bound.


The Dependence Testing Bound

The DT bound states that, given a discrete random variable X, there exists an (M, ϵ) code with a threshold decoder for which





ϵ≤custom-character[exp{−|log M−custom-character(X)]|+}].  (7)


The DT bound provides a bound on the random coding performance of a threshold decoder with threshold log γ as










ϵ






[


i


(
X
)


>

log





γ


]


+


1
M



𝕌


[


i


(
X
)




log





γ


]





,




(
8
)







where custom-character[⋅] denotes a mass with respect to the counting measure UX on X, which assigns unit weight to each x∈X.


The right-hand side of (8) can be shown to equal







M
+
1

M




times the minimum measure of the error event in a Bayesian binary hypothesis test between PX with a priori probability






M

M
+
1





and UX with a priori probability







1

M
+
1


.




This error measure can be minimized by a test that compares the log likelihood ratio






log




U
X



(
X
)




P
X



(
X
)







to the log ratio of a priori probabilities







log



M
/

(

M
+
1

)



1
/

(

M
+
1

)




,




giving





H0:PX, selected if custom-character(X)≤log M





H1:UX, selected if custom-character(X)>log M.


Taking γ=M minimizes the right-hand side of (8), which implies that the DT bound is the tightest possible bound for random coding with threshold decoding.


Particularizing the DT bound to a stationary, memoryless source with a single-letter distribution PX satisfying V(X)>0 and T(X)<∞ and invoking the Berry-Esseen inequality, an asymptotic expansion can be obtained











R
*



(

n
,
ϵ

)





H


(
X
)


+




V


(
X
)


n





Q

-
1




(
ϵ
)



+

O







(

1
n

)

.







(
9
)







Unfortunately, (9) is sub-optimal in its third-order term. Thus, random code design with threshold-based decoding fails to achieve the optimal third-order performance. Despite the sub-optimal performance of utilizing random codes with a threshold-based decoder, the use of threshold-based decoding can significantly reduce computational complexity compared to several categories of decoders including (but not limited to maximum likelihood decoders. Furthermore, the result with respect to random codes suggest that many practically imple nentable codes including (but not limited to) linear codes can also achieve acceptable levels of performance and/or performance approaching the performance of optimal codes. Accordingly, systems and methods in accordance with various embodiments of the invention can utilize any decoder (including threshold-based decoders) as appropriate to the requirements of specific applications.


The Random Coding Union Bound

The RCU bound employs random code design and maximum likelihood decoding. The RCU bound states that, given a discrete random variable X, there exists an (M, ϵ) code with a maximum likelihood decoder for which










ϵ


𝔼


[

min


{

1
,


1
M



𝔼


[



exp


(

ι


(

X
_

)


)



1


{


ι


(

X
_

)




ι


(
X
)



}



X

]




}


]



,




(
10
)







where PXX(a, b)=PX(a)PX(b) for all a, b∈X.


The same RCU bound is obtained by randomizing only over linear encoding maps. Thus, there is no loss in performance when restricting to linear compressors.


The RCU bound can be shown to recover the first three terms of the achievability result for an optimal code. Thus,the sub-optimal third-order terms in (6) and (9) result from the sub-optimal decoder rather than the random encoder design. This is important since optimal codes are not available for scenarios like MASC, which is discussed below.


The RCU bound focuses on a stationary, memoryless source with single-letter distribution PX satisfying





V(X)>0  (11)





T(X)<∞  (12)


Define constants









B


=
Δ




C
0




T


(
X
)




V


(
X
)



3
/
2








(
13
)







C


=
Δ



2


(



log





2



2





π






V


(
X
)





+

2


B


(
X
)




)



,




(
14
)







where C0 is the absolute constant in the Berry-Esseen inequality for i.i.d. random variables.


The following demonstrates that third-order-optimal achievability is possible via random coding of a stationary, memoryless source satisfying the conditions in (11) and (12). For all 0<ϵ<1,












R
*



(

n
,
ϵ

)





H


(
X
)


+




V


(
X
)


n





Q

-
1




(
ϵ
)



-


log





n


2

n


+

ξ


(
n
)




,




(
15
)







where







ξ


(
n
)


=

O


(

1
n

)






is bounded more precisely as follows.


1) For all







0
<
ϵ


1
2












and











n
>


(


B
+
C

ϵ

)

2


,












ξ


(
n
)






1
n


log





C

+


1
n





B
+
C


ϕ


(


Φ

-
1




(


Φ


(


Q

-
1




(
ϵ
)


)


+


B
+
C


n



)


)



.







(
16
)







2) For all








1
2

<
ϵ

1











and











n
>


(


B
+
C


ϵ
-

1
2



)

2


,












ξ


(
η
)






1
n


log





C

+


1
n





B
+
C


ϕ


(


Q

-
1




(
ϵ
)


)



.







(
17
)







Accordingly, the RCU bound demonstrates that systems and methods in accordance with various embodiments of the invention can achieve third-order optimal performance using random codes. While various systems and methods that employ random codes to perform source coding are described above with reference to the system shown in FIG. 1, encoders and decoders in accordance with many embodiments of the invention can be utilized in MASC applications and/or in the newly described RASC applications detailed below. Systems and methods that can perform MASC and/or RASC in accordance with various embodiments of the invention are discussed further below.


Multiple Access Source Coding

To simplify notation when discussing systems that are utilized tea perform MASC, the following discussion focuses on MASC with two encoders. However, the definitions and results generalize to more than two encoders. Accordingly, it should be readily appreciated that systems and methods in accordance with many embodiments of the invention can utilize two or more sources and two or more decoders (e.g., Slepian Wolf decoders) as appropriate to the requirements of specific applications. In MASC, also known as a Slepian-Wolf source coding, independent encoders compress a pair of random variables (X1, X2) with discrete alphabets X1 and X2. Encoder i, i∈[2], observes only Xi, which it maps to a codeword in [Mi]; a single decoder jointly decodes the pair of codewords to reconstruct (X1, X2). Codes can be defined for abstract random objects and then particularized to random objects that live in an alphabet endowed with a Cartesian product structure.


By way of definition, an (M1, M2, ϵ) MASC for random variables (X1, X2) with discrete alphabets X1 and X2 includes two encoding functions f1: X1→[M1] and f2: X2→[M2] and a decoding function, g: [M1]×[M2]→X1×X2 with error probability






custom-character[g(f1(X1),f2(X2))≠(X1,X2)]≤ϵ.  (18)


In block coding, encoders individually observe X1n and X2n drawn from distribution PX1nX2n on X1n×X2n. An (n, M1, M2, ϵ) block MASC is an (M1, M2, ϵ) MASC for random vectors (X1n, X2n) on X1n×X2n. The code rate for the block MASC R=(R1, R2) is given by











R
1



=
Δ




1
n


log






M
1



,


R
2



=
Δ




1
n


log







M
2

.







(
19
)







Rate R=(R1, R2) is (n, ϵ)-achievable if there exists an (n, M1, M2, ϵ) MASC with







R
1




1
n


log






M
1







and






R
2




1
n


log







M
2

.






The (n, ϵ)-rate region custom-character*(n, ϵ) is the closure of the set of (n, ϵ)-achievable rate pairs.


While the above definitions can apply to arbitrary discrete random variables (X1i, X2i), i=1, 2, . . . , with transition probability kernels P(X1X2)i|(X1X2)i−1, the asymptotic analysis discussed below focuses on stationary, memoryless sources, where P(X1X2)i|(X1X2)i−1=PX1X2 for all i=1, 2, . . .


For stationary, memoryless sources with rate R=(R1, R2) and distribution PX1X2, the following can be defined








R
¯



=
Δ



[




R
1






R
2







R
1

+

R
2





]


,


H
¯



=
Δ




{




H


(


X
1

|

X
2


)







H


(


X
2

|

X
1


)







H


(


X
1

,

X
2


)





}

.






It can be proven if (X1n, X2n) are stationary and memoryless, then for every ϵ∈(0, 1),














lim

n








(

n
,
ϵ

)



=

{


(


R
1

,

R
2


)

:

R
1










H


(


X
1

|

X
2


)







R
2







H


(


X
2

|

X
1


)








R
1

+

R
2










H


(


X
1

,

X
2


)


}

,







(
21
)







(i.e., the strong converse holds). This region can be referred to as the asymptotic MASC rate region. While the above result is presented with respect to stationary and memoryless sources, systems and methods in accordance with various embodiments of the invention can be utilized with any of a variety of different sources as appropriate to the requirements of specific applications.


Development of systems and methods for use in the MASC in accordance with various embodiments of the invention involved development of a MASC RCU bound, extending the RCU bound introduced above to the multiple-encoder case.


Given discrete random variables (X1, X2) there exists an (M1, M2, ϵ) MASC with


















ϵ


𝔼




[

min






{

1
,


A
1

+

A
2

+

A

1

2




}


]















where











A
1



=
Δ




1

M
1





𝔼




[


exp






(

i


(



X
_

1


|

X
2


)


)


1


{





(


i


(



X
_

1


|

X
2


)




i


(


X
1

|

X
2


)



}













X
1


,

X
2


]














A
2



=
Δ




1

M
2




𝔼




[


exp






(

i


(



X
_

2


|

X
1


)


)


1


{





(


i


(



X
_

2


|

X
1


)




i


(


X
1

|

X
2


)



}













X
1


,

X
2


]










A
12



=
Δ




1


M
1



M
2





𝔼




[


exp






(

i


(



X
_

1


|

X
2


)


)


1


{





(


i


(



X
_

1


|

X
2


)




i


(


X
1

,

X
2


)



}













X
1


,

X
2


]











P


X
1



X
2




X
_

1




X
_

2




X
_

1





X
_

2






(

a
,
b
,

a
_

,

b
_

,


a
_









,


b
_










)


=



P


X
1



X
2





(

a
,
b

)





P


X
1



X
2





(


a
_

,

b
_


)









P



X
1





X
2



(


a
_












b



)




P



X
1





X
2



(


b
_












a

)

.




(
22
)








FIGS. 2A and 2B plots the point-to-point and MASC RCU bounds. Rate-blocklength trade-offs at ϵ=10−1 (FIG. 2A) and ϵ=10−3 (FIG. 2B) for a pair of binary, stationary, memoryless sources with joint distribution pX1X2(0, 0)=½, pX1X2(0, 1)=pX1X2(1, 0)=pX1X2(1, 1)=⅙ are shown. The plot labelled point-to-point (P2P) hypothesis testing (HT) converse coincides with the optimum R*(n, ϵ).


Asymptotics: Third-Order MACS Rate Region

The following third-order asymptotic characterization of the MASC rate region for stationary, memoryless sources closes the






O


(


log





n

n

)





gap between previous published bounds based upon second order characterization of the MASC rate-region for finite-alphabet stationary, memoryless sources in terms of the asymptotic rate region and the entropy dispersion matrix.


In order to develop a third-oreder asymptotic characterization, memoryless sources with single-letter joint distribution PX1X2 were considered for which





V(X1,X2)>0,custom-character[Vc(X1|X2)]>0,custom-character[V(X2|X1)]>0,  (27)





T(X1,X2)>∞,T(X1|X2)>∞,T(X2|X1)>∞,  (28)






custom-character[Tc2(X1|X2)]>∞,custom-character[Tc2(X2|X1)]>∞.  (29)


When (27) holds, rank (V)≥1 (in fact, the weaker condition V(X1|X2)>0, V(X2|X1)>0, V(X1, X2)>0 suffices). Assumption (28) is also satisfied automatically if the alphabets X1 and X2 are finite.


The following set can be defined












(

n
,
ϵ

)




=
Δ




{




R
_



3


:

R
_


=


H
_

+



inv



(

V
,
ϵ

)



n


-



log





n


2





n







1



}

.





(
30
)







where vector H is defined in (20), V is the entropy dispersion matrix for (X1, X2), and custom-characterinv(V, ϵ) is defined as the multidimensional counterpart of the function Q−1(⋅) as






custom-character
inv(V,ϵ)custom-character{z∈custom-characterd:custom-character[Z≤z]≤1−ϵ}.  (31)


Note that custom-character*(n, ϵ)⊂custom-character2 but custom-character*(n, ϵ)⊂custom-character3. The inner and outer bounding sets can be defined as follows











in
*



(

n
,
ϵ

)




=
Δ



{



R
_



2


:


R
_










(

n
,
ϵ

)


+


O


(

1
n

)



1




}





(
32
)








out
*



(

n
,
ϵ

)




=
Δ




{



R
_



2


:


R
_










(

n
,
ϵ

)


-


O


(

1
n

)



1




}

.





(
33
)







When a pair of stationary, memoryless sources is utilized with single-letter joint distribution PX1X2 satisfying (27)-(29), for any 0<ϵ<1, the (n, ϵ)-rate region custom-character*(n, ϵ) satisfies






custom-character*in(n,ϵ)⊆custom-character*(n,ϵ)⊆custom-character*out(n,ϵ).  (34)


Since the upper and lower bounds agree up to their third-order terms, we call custom-character(n, ϵ) the third-order MASC rate region.



FIG. 3 plots the boundaries of custom-character(n, ϵ) at different values of n for an example pair of sources. The third-order MASC rate regions custom-character(n, ϵ) at ϵ=10−3 are shown for stationary, memoryless sources (X1, X2) with pX1,X2(0, 0)=½, pX1,X2(1, 1)=pX1,X2(1, 1)=⅙.


For point-to-point source coding, zero varentropy means that the source is uniform; the






-


log





n


2

n






third-order term is absent in that case. While condition (27) limits the third-order MASC rate region considered above to sources with positive varentropies, the case where one or more varentropies are zero is discussed below. Roughly, each zero varentropy yields a zero dispersion, and the absence of a






-


log





n


2

n






third-order term, similar to the point-to-point case. Furthermore, if V(X1|X2)>0 but custom-character[Vc(X1|X2)]=0, the corresponding achievable third order term increases from






-


log





n


2

n






to 0.


In this context, the third-order MASC rate region described above can be generalized to any finite number of encoders. Let custom-charactercustom-character be a nonempty ordered set with a unique index for each encoder. For any vector custom-charactercustom-character, the (custom-character−1)-dimensional vector of its partial sums can be defined as










R
𝒯



=
Δ




(





i

A




R
i


,


𝒯
^



𝒫


(
𝒯
)




)

.





(
35
)







For any distribution custom-character defined on custom-character and any custom-charactercustom-character, (custom-character−1)-dimensional vectors can be defined as follows:






custom-character(custom-character)custom-character(custom-character(custom-character|custom-character),custom-charactercustom-character(custom-character))  (36)






custom-character
custom-character
custom-character[custom-character(custom-character)],  (37)


and (custom-character−1)×(custom-character−1) entropy dispersion matrix






custom-character
custom-characterCov[custom-character(custom-character)]  (38)


for random vector custom-character. The following set can be defined












(

n
,
ϵ

)




=
Δ





H
_

𝒯

+



inv



(


V
𝒯

,
ϵ

)



n


-



log





n


2





n







1.






(
39
)







Thus, custom-character(n, ϵ)⊂custom-character while custom-character(n, ϵ)⊂custom-character−1. Finally,












in
,
𝒯

*



(

n
,
ϵ

)




=
Δ



{



R
𝒯





𝒯




:



R
_

𝒯











𝒯



(

n
,
ϵ

)


+


O


(

1
n

)



1




}





(
40
)









out
,
𝒯

*



(

n
,
ϵ

)




=
Δ




{



R
𝒯





𝒯




:



R
_

𝒯











𝒯







(

n
,
ϵ

)


-


O


(

1
n

)



1




}

.





(
41
)







If every element of custom-character(custom-character) has a positive variance and a finite third center moment, then for any 0<ϵ<1,






custom-character(n,ϵ)⊆custom-character(n,ϵ)⊆custom-character(n,ϵ).  (42)


Comparison with Point-To-Point Sorce Coding


FIGS. 2A and 2B compare joint (point-to-point) compression of (X1n, X2n) to the MASC sum-rate at the symmetrical rate point, (R1=R2). The gap between the MASC and point-to-point HT converses captures a penalty due to separate encoding. For small n, the third-order Gaussian approximation (without the






O


(

1
n

)





term) is more accurate at ϵ=10−1 than at ϵ=10−3 since the






O


(

1
n

)





term blows up as ϵ approaches 0.


It can be shown that optimal MASC incurs no first-order penalty in achievable sum rate when compared to joint coding A quantity known as the local dspersion can be utilized to characterize the second-order speed of convergence to any asymptotic MASC rate point from any direction. For any non-corner point on the diagonal boundary of the asymptotic MASC rate region, the sum rate's second-order coefficient is optimal when approached either vertically or horizontally. Approaching corner points can incur a positive second-order penalty relative to point-to-point coding. The MASC penalty can be bounded by considering the achievable sum rate R1+R2 for different choices of R1 and R2. The cases where X1 and X2 are dependent and X1 and X2 are independent can be separately addressed, assuming throughout that (27) and (28) hold.


When X1 and X2 are dependent, H(X1)+H(X2)>H(X1, X2)>H(X1|X2)+H(X2|X1), and the asymptotic sum-rate boundary contains non-corner and corner points. It can be shown that a MASC incurs no first-, second-, or third-order performance penalty relative to joint coding at non-corner points (i.e., when R1<H(X1) and R2<H(X2)). In contrast, a MASC suffers a second-order performance penalty at corner points (i.e., when R1=H(X1) or R2=H(X2)). See FIG. 4A for an illustration. FIGS. 4A and 4B are schematic plots of the (n, ϵ)-rate region and the third-order-optimal sum rate when X1, X2 are dependent (FIG. A), X1, X2 are independent (FIG. 4B). In FIG. 4A, the boundary of custom-character*(n, ϵ) between H(X1) and H(X2) (excluding the end points) contains rate points that achieve the optimal point-to-point rate up to the third order, while the end points do not achieve that optimal rate. The value of r* in FIG. 4A is defined in (49). The values of r1*, r2* in FIG. 4B are defined as










(


r
1
*

,

r
2
*


)

=

arg







min



(


r
1

,

r
2


)

:




Φ


(

r
1

)




Φ


(

r
2

)





1
-
ϵ







(




V


(

X
1

)





r
1


+



V


(

X
2

)





r
2



)

.







(
43
)







Suppose that X1 and X2 are dependent.


1. The constants δ1, δ2, G>0 and ϵ∈(0, 1) can be fixed. Then there exists some constant n(δ1, δ2, G) such that if










R
1




H


(

X
1

)


-

δ
1






(
44
)







R
2




H


(

X
2

)


-

δ
2






(
45
)








R
1

+

R
2


=


H


(


X
1

,

X
2


)


+





V


(


X
1

,

X
2


)







n









Q

-
1




(

ϵ
-

G

n



)



-


log





n


2





n







(
46
)









    • then R=(R1, R2)∈custom-character(n, ϵ) for all n>n(δ1, δ2, G).





2. When ϵ∈(0, 1) is fixed. If










R
2




H


(


X
1



X
2


)


+


r
*


n


-


log





n


2





n


+

G
n






(
47
)









    • for some G>0, then R=(H(X1), R2)∈custom-character(n, ϵ). Conversely, if R=(H(X1), R2)∈custom-character(n, ϵ), then














R
2




H


(


X
2

|

X
1


)


+


r
*


n


-


log





n


2

n




,




(
48
)









    • where r* is the solution to equation








Φ(V2;r,r)=1−ϵ,  (49)

    • and V2 is the covariance matrix for random vector (custom-character(X2|X1), (custom-character(X1, X2)).


For independent sources, the asymptotic sum-rate boundary contains only the single (corner) point (R1, R2)=(H(X1), H(X2)), and the entropy dispersion matrix







[




V


(

X
1

)




0




0



V


(

X
2

)





]

.




is singular.


The next result concerns the third-order-opitimal sum rate







R

sum*(n,ϵ)custom-charactermin{R1+R2:custom-characterR=(R1,R2) s.t. Rcustom-character*(n,ϵ)}.  (50)


According to the third-order MASC rate region described above, Rsum*(n, ϵ) characterizes the best achievable sum rate in SW source coding up to an






O


(

1
n

)





gap.


For independent sources a unique (r1*, r2*) can capture the best MASC second-order sum-rate; the third-order term is achieved at all points on a segment of the rate region boundary. See FIG. 4B. Under assumption (27),













min


(


r
1

,

r
2


)

:





Φ


(

r
1

)




Φ


(

r
2

)





1
-
ϵ









(




V


(

X
1

)





r
1


+



V


(

X
1

)





r
2



)


>




V


(

X
1

)


+

V


(

X
2

)







Q

-
1




(
ϵ
)




,




(
51
)







where V(X1)+V(X2)=V(X1, X2) for (X1, X2) independent. Here (51) follows since its left-hand side solves











min

(


a
1

,

a
2


)








(


a
1

+

a
2


)











s
.
t
.

Φ


(


a
1



V


(

X
1

)




)





Φ


(


a
2



V


(

X
2

)




)





1
-
ϵ


,





(
52
)







and the constraint in (52) requires a1>√{square root over (V(X1))}Q−1(ϵ) and a2>√{square root over (V(X2))}Q−1(ϵ) , which gives the bound since





√{square root over (V(X1))}+√{square root over (V(X2))}>√{square root over (V(X1)+V(X2))}.  (53)


Therefore, when X1 and X2 are independent, a MASC approach incurs a positive second-order sum-rate penalty relative to joint coding. Closed-form expressions for this penalty are available in special cases. When V(X1)−V(X2), r1*=r2*=Q−1(1−√{square root over (1−ϵ)}), and the penalty is










2




V


(

X
1

)


n





Q

-
1




(

1
-


1
-
ϵ



)



-




2


V


(

X
1

)



n






Q

-
1




(
ϵ
)


.






(
54
)







When and X1 and X2 are i.i.d., the penalty equals the penalty fir coding a vector X2n of 2n i.i.d. outputs from PX by applying an independent (n, ϵ) (point-to-point) code with error probability 1−√{square root over (1−ϵ)} to each of (X1, . . . , Xn) and (Xn+1, . . . , X2n) instead of a single (2n, ϵ) code to vector X2n.


Limited Feedback and Cooperation

Systems and methods in accordance with various embodiments of the invention that are employed in MASC scenarios can utilize limited feedback. In several embodiments, the decoder broadcasts the same custom-character bits of feedback to all encoders. A bit sent at time i is a function of the encoder outputs received in time steps 1, . . . , i−1. (See FIG. 5A.) The impact of feedback can be bounded by studying a MASC with a cooperation facilitator (CF). The CF broadcasts the same custom-character-bit function of the sources to both encoders prior to their encoding operations. (See FIG. 5B.) Since the MASC network has no channel noise, feedback from the decoder cannot convey more information than feedback from the CF. As a result, the impact of feedback can be bound by bounding the impact of cooperation.


The CF-MASC and its rate region can be defined as follows. An (L, M1, M2, ϵ) CF-MASC for random variables (X1, X2) on X1×X2 includes a CF function L, two encoding functions f1 and f2, and a decoding function g given by





L: X1×X2→[L]





f1: [L]×X1→[M1]





f2: [L]×X2→[M2]





g: [M1]×[M2]→X1×X2,


with error probability






custom-character[g(f1(L(X1,X2),X1),f2(L(X1,X2),X2))≠(X1,X2)]≤ϵ.


An (n, L, M1, M2, ϵ) MASC is a CF MASC for random variables (X1n, X2n) on X1n×X2n. The code's finite blocklength rates are defined by








R
1

=


1
n


log






M
1



,






R
2

=


1
n


log







M
2

.







A rate pair (R1, R2) is (n, custom-character, ϵ)-CF achievable if there exists an (n, L, M1, M2, ϵ) CF-MASC with M1≤exp(nR1), M2≤exp(nR2), and L≤exp(custom-character). The (n, custom-character, ϵ)-CF rate region custom-characterCF*(n, custom-character, ϵ) can be defined as the closure of the set of all (n, custom-character, ϵ)-CF achievable rate pairs.



custom-character
FB*(n, custom-character, ϵ) can be used to denote the feedback-MASC (FB-MASC) rate region, which can be defined as the closure of the set of all (n, ϵ)-achievable rate pairs when the same custom-character bits of feedback from the decoder are available to both encoders.


Since the CF sees the source vectors while the decoder sees a coded description of those vectors (using a deterministic code), an custom-character-bit CF can implement any function used to determine the decoder's custom-character-bit feedback. As a result, any rate point that is achievable by an custom-character-bit FB-MASC is also achievable by an custom-character-bit CF-MASC. Therefore, for any 0<ϵ<1 and custom-character<∞,






custom-character
FB*(n,custom-character,ϵ)⊆custom-characterCF*(n,custom-character,ϵ).  (56)


Consider stationary, memoryless sources with single-letter distribution PX1X2 satisfying (27) and (28). For any 0<ϵ<1 and custom-character<∞,






custom-character
CF*(n,custom-character,ϵ)⊆custom-characterout*(n,ϵ).  (57)


Thus custom-characterCF*(n, custom-character, ϵ) and custom-character(n, ϵ) share the same outer bound.


Accordingly, the CF-MASC (and FB-MASC) performance is bounded, so that for any custom-character<∞, the third-order rate region for custom-character-bit CF-MASCs cannot exceed the corresponding MASC rate region. Hence finite feedback does not enlarge the third-order (n, ϵ) MASC rate region. This result generalizes to scenarios with more than two encoders.


While a variety of systems and methods for performing MASC in a manner that can achieve third-order optimal performance are described above with reference to FIGS. 2A-5B, any of a variety of systems and methods can be implemented in order to perform MASC with or without feedback and/or the use of a CF as appropriate to the requirements of specific applications in accordance with various embodiments of the invention. Furthermore, systems and methods in accordance with a number of embodiments of the invention can be utilized for the newly described RASC applications, which are discussed further below.


Random Access Source Coding

Systems and methods in accordnace with many embodiments are used to perform RASC. RASC is a generalization of MASC for networks where the set of participating encoders is unknown to both the encoders and the decoder a priori.


Definitions and Coding Strategy

In a RASC system, an assumption can be made that there is maximum number of active encoder K<∞. Accordingly, each encoder with a unique source from the set of sources can be indexed by [K]. In addition, each encoder can choose whether to be active or silent. Only sources associated with active encoders are compressed and reconstructed. By assumption, the decision to remain silent is independent of the observed source instance. Given the joint distribution PX[K] on countable alphabet X[K], when ordered set custom-charactercustom-character([K]) of [K] is active, the marginal on the transmitted sources is












P

X
𝒯




(

X
𝒯

)


=





ϰ


[
K
]



\


𝒯



χ


[
K
]


\

𝒯









P

X

[
K
]





(

X

[
K
]


)




,




X
𝒯







χ
𝒯

.







(
58
)







Thus, each encoder's state has no effect on the statistical relationship among sources observed by other encoders.


Systems and methods in accordance with many embodiments of the invention involve the use of a RASC process that organizes communication into epochs. At the beginning of each epoch, each encoder can independently deckle its activity state; that activity state remains unchanged until the end of the epoch. Thus, the active encoder set custom-character is fixed in each epoch.


A process 600 that can be implemented by an encoder within an RASC system in accordance with an embodiment of the invention is illustrated in FIG. 6. As illustrated in FIG. 6, each active encoder i∈custom-character observes source output Xi∈Xi. When a determination (602) is made that a source output is available, the encoder independently maps (604) the source data to a codeword that includes a sequence of code symbols from alphabet [Qi]. When a start of an epoch is detected (606), the |custom-character| codewords are sent (608) simultaneously to the decoder. Since set custom-character is unknown a priori, the encoder behavior cannot vary with custom-character. The decoder sees custom-character and decides a time custom-character, called the decoding blocklength, at which to jointly decode all received partial codewords. The set of potential decoding blocklengths custom-charactercustom-character{custom-character: custom-charactercustom-character([K])} is part of the code design; it is known to all encoders and to the decoder. When the decoder performs the joint decoding, the decoder provides feedback to the encoder indicating (610) that the epoch has ended. The encoding and transmission process 600 repeats as more data becomes available to each encoder and as each new epoch commences.



FIGS. 7 and 8 illustrate an RASC scheme employed in accordance with an embodiment of the invention during one epoch when custom-character=[k]. Each encoder i∈custom-character sends a single code symbol per time step. At each time m∈{m′∈custom-character: m′<custom-character}, the decoder sends a “0” to indicate that it is not yet ready to decode; at time m=custom-character, the decoder sends a “1,” ending one epoch and starting the next. The decoder then reconstructs source vector custom-character using the first custom-character code symbols from each active encoder. To avoid wasting time in an epoch with no active encoders, decoding time m=1 is included in set custom-character. The decoder sends at most 2K bits of feedback, and encoders need only listen for decoder feedback at the times in set custom-character.


To formalize the above strategy, K can be fixed to K≥1. In addition, the following vectors can be defined






ϵ
K
custom-character(custom-character,custom-charactercustom-character([K]))  (59)






m
K
custom-character(custom-character,custom-charactercustom-character([K])∪{∈})  (60)


with m=1 and mmaxcustom-charactermax{custom-character: custom-charactercustom-character([K])}.


An (mK, Q[K], ϵK) RASC for sources X[K] on source alphabet X[K] comprises a collection of encoding and decoding functions












f
i



:



χ
i





[

Q
i

]


m
max



,

i


[
K
]


,




(
61
)









g
𝒯







:




i

𝒯




[

Q
i

]



m
𝒯





χ
𝒯


,

𝒯


𝒫


(

[
K
]

)



,




(
62
)







where fi is the encoding function for source Xi and custom-character is the decoding function for active coder set custom-character. For each custom-charactercustom-character([K]), source vector custom-character is decoded at time custom-character with error probability custom-character[custom-character(fi(Xi)[custom-character], i∈custom-character)≠custom-character]≤custom-character, where fi(xi)[m] denotes the first m code symbols of fi(xi).


An (n, mK, Q[K], ϵK) RASC is an RASC for an m-block of source outcomes. The parameter n, called the encoding blocklength does not vary with custom-character. An (mK, Q[K], ϵK) RASC behaves, for each custom-character, like a ((custom-character, custom-charactercustom-character), custom-character) MASC with a finite number |{m∈custom-character: m≤custom-character}| of feedback bits. However, the RASC is one code. Its descriptions are nested so that the previously transmitted symbols are a prefix for the codeword of the next larger blocklength (i.e., for each custom-charactercustom-character, if custom-character′<custom-character, then fi(xi)custom-character is a prefix of fi(xi)[custom-character]). It simultaneously satisfies the error constraints for all custom-charactercustom-character([K]). And, since the code symbol alphabet sizes Q[K] are fixed, its rate vectors are coupled. See FIG. 9, which shows the relationship between decoding blocklength custom-character, code symbol alphabet sizes (Q1, Q2), and source coding rate vector custom-character (illustrated for custom-character={1, 2}).


The following definitions can be utilized to build a non-asymptotic fundamental limit of RASCs.


A collection custom-character of rate vectors is n-valid if ∃ (mK, Q[K]) s.t.











R
𝒯

=


1
n



(



m
𝒯


log






Q
i


,

i

𝒯


)



,



𝒯



𝒫


(

[
K
]

)


.







(
63
)







The set custom-charactervalid(n) is the set of n-valid rate collections. The collection is (n, ϵK)-achievable if there exists an (n, mK, Q[K], ϵK) RASC scheme. The (n, ϵK)-rate set custom-character* (n, ϵK) is the set of (n, ϵK)-achievable rate collections.


Asymptotics: Third-Order Performance of the RASC

The performance of an (n, mK, Q[K], ϵK) RASC for stationary, memoryless sources can be analyzed. Results include both achievability and converse characterizations of the (n, ϵK)-rate set custom-character* (n, ϵK) under the assumption that the single-letter joint source distribution PX[K] satisfies






custom-character[Vc(custom-character|custom-character)]>0 ∀custom-charactercustom-character⊆[K],custom-character,custom-character≠∅  (64)





T(custom-character|custom-character)>∞ ∀custom-charactercustom-character⊆[K],custom-character,custom-character≠∅  (65)






custom-character[Tc2(custom-character|custom-character)]>∞ ∀custom-charactercustom-character⊆[K],custom-character,custom-character≠∅  (66)


Constraints (64)-(66) enable the use of Berry-Esseen bounds in subsequent analysis. The resulting characterization is tight up to the third-order term. While the existence of an (n, mK, Q[K], ϵK) RASC implies the existence of an (n, (custom-character, i∈custom-character), custom-character) MASC for each custom-charactercustom-character([K]), the existence of individual MASCs does not imply the existence of a single RASC that simultaneously satisfies the error probability constraints for all possible configurations of active encoders. Indeed, the ability of systems and methods in accordance with various embodiments of the invention to employ a single RASC that simultaneously performs as well (up to the third-order term) as the optimal MASC for each possible nuinber of encoders (i.e. each custom-charactercustom-character([K])) is an important and surprising result.


The inner and outer bounding sets can be defined as follows:






custom-character
in*(n,ϵK)custom-character{(Rcustom-charactercustom-charactervalid(n): custom-charactercustom-character(n, ϵcustom-character)∀custom-charactercustom-character([K])}  (67)






custom-character
out*(n,ϵK)custom-character{(Rcustom-charactercustom-charactervalid(n): custom-charactercustom-character(n, ϵcustom-character)∀custom-charactercustom-character([K])},  (68)


where custom-character(n, ϵ) and custom-character(n, ϵ) are the third-order MASC bounding sets for distribution custom-character. (See (40) and (41).)


Third-Order RASC Performance

For any K<∞, stationary, memoryless sources specified by a single-letter joint distribution PX[K] satisfying (64)-(66) can be considered. For any 0<ϵK<1,






custom-character
in*(n,ϵK)⊆custom-character*(n,ϵK)≣Rout*(n,ϵK).  (69)


The achievability of third-order RASC performance provides a sufficient condition for the existence of a single RASC that is simultaneously good for all custom-charactercustom-character([K]). To prove this, an achievability result can be derived assuming that the encoders and decoder share the common randomness used to generate a random code. Unfortunately, the existence of a random code ensemble with expected error probability satisfying the error probability constraint for each custom-charactercustom-character([K]) does not guarantee the existence of a single deterministic code satisfying those constraints simultaneously. Therefore a different approach can be taken in accordance with many embodiments of the invention, which, unexpectedly, combines a converse bound on error probability and a random coding argument to show achievability.


The following refinement of the random coding argument provides a bound on the probability (with respect to the random code choice) that the error probability of a randomly chosen code exceeds a certain threshold. The code of interest here call be any type of source or channel code.


Given any RASC c, for each custom-charactercustom-character([K]) let custom-character(c) denote the error probability of code c under active encoder set custom-character. Before proceeding with the proof, a random code ensemble is defined in order to calculate its expected error probability.


For any K<∞, consider a source distribution PX[K] defined on countable alpha-bet X[K]. There exists a random code ensemble C defined on the set of all RASCs with decoding blocklengths mK and code alphabets Q[K] for which the following inequalities hold simultaneously for all custom-charactercustom-character([K]):
















𝔼


[


P

e
,
𝒯




(
C
)


]




𝔼


[

min


{

1
,





𝒯
^



ϵ𝒫


(
𝒯
)







exp


(


-

m
𝒯


·


Q
_



(

𝒯
^

)



)




A

𝒯
^





}


]



,








where





(
70
)













Q
_



(

𝒯
^

)




=
Δ






i





ϵ


𝒯
^





log






Q
i








(
71
)







A

𝒯
^




=
Δ



𝔼


[


exp
(


i


(



X
_


𝒯
^






X

𝒯

\


𝒯
^



)


)



.1


{


i


(



X
_


𝒯
^




X

𝒯

\


𝒯
^




)




i


(


X

𝒯
^




X

𝒯

\


𝒯
^




)



}






X
𝒯


]






(
72
)







and the expectation in (72) is with respect to the conditional distribution






custom-character=custom-character  (73)


When parameters (n, Q[K], ϵK) are fixed, increasing custom-character yields larger decoding blocklengths custom-character. Therefore, the choice of (custom-character to satisfy














𝒯ϵ𝒫


(

[
K
]

)





1


λ
𝒯

+
1



<
1

,




(
74
)







controls the RASC performance trade-off across different active encoder sets. This trade-off affects the performance of the RASC in the fourth- or higher-order terms.


RASC for Permutation-Invariant Sources

A permutation-invariant source is defined by the constraint






P
X

[K]
(x[K])=PX[K](xπ([K]))  (75)


for all permutations π on [K] and all x[K]∈X[K]. For example, given any PS and PX|S, the marginal PX[K] of PX[K]S=(PX|S)KPS satisfies (75). Such “hidden variable” models have applications in statistics, science, and economics, where latent variables (e.g., the health of the world economy or the state of the atmosphere) influence observabies (e.g., stock prices or climates). FIG. 10 shows an example with K sensors reading measmements of a common hidden state S.


Permutation-invariant source models can be of interest both because of their wide applicability and because they present an opportunity for code simplification through identical encoding, where all encoders employ the same encoding map. For any permutation-invariant source, (58) and (75) imply that Xi=X for all i∈[K] and, for any custom-charactercustom-character([K]) with |custom-character|=k,






custom-character=PX[k].  (76)


Thus, custom-character is permutation-invariant for every custom-character and the joint source distribution depends on the number of active encoders but not their identities. Assuming the same error probability ϵk is further employed for all custom-charactercustom-character([K]) with |custom-character|=k, a single decoding blocklength can be fixed for each number k∈[K] of active encoders. Furthermore, identical encoders can be utilized at all transmitters, allowing the accommodation of an arbitrarily large number of encoders without designing a unique encoder for each source.


In analyzing RASC performance with identical encoders on a permutation-invariant source, an assumption can be made in addition to (64) and (65) that no two sources are identical, i.e.,












[




i
,




[
K
]


,

i

j





{


X
i

=

X
j


}


]


<
1.




(
77
)







This can be important since using identical encoders on identical sources yields identical descriptions, in which case descriptions from multiple encoders are no better than descriptions from a single encoder. Under these assumptions, third-order RASC performance continues to hold. A number of embodiments of the invention utilize a modified decoder that outputs the most probable source vector custom-charactercustom-character that contains no repeated symbols, treating the case where custom-character contains repeated symbols as an error. In several embodiments, the system permits multiple identical sources and the decoder relies upon a system configuration in which multiple identical sources are active with low likelihood. In certain embodiments, a decoder can detect that multiple identical sources are active and can command all except one of the identical sources to cease activity. In several embodiments, a source encoder ceases activity for a specified and/or random period of time. In the asymptotic analysis for stationary, memoryless sources, the probability of this error event is bounded by














[




i
,

j


[
K
]


,

i

j





{


X
i
n

=

X
j
n


}


]





(




[




i
,

j


[
K
]


,

i

j





{


X
i

=

X
j


}


]


)

n


,




(
78
)







which decays exponentially in n by (77). Therefore, under the assumption in (77), identical encoding does not incur a first- second-, or third-order performance penalty.


MASC Using Sources with Less Redundancy

Referring again to systems and methods in accordance with various embodimennts of the invention that employ MASC, the asymptotic achievability result for the MASC schemes described above requires that all V(X1, X2), V(X1|X2), and V(X2|X1) are strictly positive (as an implication of assumption (27)). Thus, the analysis presented above can break down when any of these varentropies is equal to zero. Such a source can be referred to as being less redundant. Accordingly, the performance systems and methods that perform MASC in accordance with a number of embodiments of the invention can be analyzed when utilized in combination with less redundant sources. Specifically, a pair of stationary, memoryless sources can be analyzed with respect to the following three cases:


1) all three varentropies are equal to zero;


2) exactly two of the varentropies are equal to zero;


3) exactly one of the varentropies is equal to zero.


The assumption can be maintained that the joint distribution PX1X2 satisfies (28) and (29). For those cases in which V(X2|X1)>0 the assumption can be maintained that custom-character[Vc(X2|X1)]>0. Furthermore, if V(X2|X1)>0, the assumption can be maintained that custom-character[Vc(X1|X2)]>0.


When encoders are required to operate independently in a MASC, it is typically not a simple matter to find optimal codes. In Section A below, characterizations of the (n, ϵ)-rate region are provided with respect to the three general cases listed above. Then, in Section B, the case where PX1X2(x1, x2)>0 for every (x1, x2)∈X1×X2 is considered; under this condition, the optimal codes can be found and analyzed directly.


Section A: General Characterizations of the (n, ϵ)-Rate Region

The results in the three general cases are listed below.


Case 1): Suppose that V(X1|X2)=0, and V(X2|X1)=0, and V(X1, X2)=0. For any δ1, δ2, δ12>0, let














in

(
1
)




(

n
,

δ
1

,

δ
2

,

δ
12


)




=
Δ





{


(


R
1

,

R
2


)





2

:










R
1





H
(



X
1





X
2

)


+


1
n


log


1

δ
1












R
2





H
(



X
2





X
1

)


+


1
n


log


1

δ
2















R
1

+

R
2







H
(


X
1

,

X
2


)

+


1
n


log


1

δ
12





}

.







(
79
)







Define











in

(
1
)




(

n
,
ϵ

)




=
Δ








δ
1

,

δ
2

,


δ
12

>
0










δ
1

+

δ
2

+

δ
12


=
ϵ






in

(
1
)




(

n
,

δ
1

,

δ
2

,

δ
12


)







(
80
)











out

(
1
)




(

n
,
ϵ

)




=
Δ





{


(


R
1

,

R
2


)





2

:










R
1





H
(



X
1





X
2

)


+


1
n


log


1

1
-
ϵ












R
2





H
(



X
2





X
1

)


+


1
n


log


1

1
-
ϵ















R
1

+

R
2







H
(


X
1

,

X
2


)

+


1
n


log


1

1
-
ϵ





}

.







(
81
)







When V(X1|X2)=0, V(X2|X1)=0, and V(X1, X2)=0, the (n, ϵ)-rate region custom-character*(n, ϵ) satisfies






custom-character
in
(1)(n,ϵ)⊆custom-character*(n,ϵ)⊆custom-characterout(1)(n,ϵ).  (82)


As in the point-to-point scenario, there are no second-order dispersion terms or






-


log





n


2

n






third-order terms in the characterization of custom-character*(n, ϵ) in this case. For any n and ϵ, the achievable region custom-characterin(1)(n, ϵ) has a curved boundary due to the trade-off in the






O
(

1
n

)




fourth-order terms, while the converse region custom-characterout(1)(n, ϵ) has three linear boundaries.


Case 2): There are three possible cases where exactly two of the three varentropies are equal to zero. Here, we suppose that V(X1|X2)>0 while V(X2|X1)=V(X1, X2)=0. The other two cases can be analyzed in the same way.


When V(X1|X2)>0, V(X2|X1)=0, and V(X1, X2)=0, the (n, ϵ)-rate region custom-character*(n, ϵ) satisfies






custom-character
in
(2)(n,ϵ)⊆custom-character*(n,ϵ)⊆custom-characterout(2)(n,ϵ).  (83)


The achievable region custom-characterin(2)(n, ϵ) has a curved boundary due to the trade-off in δ1, ϵ2, and δ12. If











δ
1

=

ϵ
-

2

n




,


δ
2

=

1

n



,


δ
12

=

1

n



,




(
84
)







then it is apparent that the dispersion corresponding to R1 is V(X1|X2) with a






-


log





n


2

n






third-order term, while the dispersions of R2 and R1+R2 are zero.


Case 3): Similar to Case 2), there are three possible cases where exactly one of the three varentropies is equal to zero. Here, the case where V(X1|X2)=0 while V(X2|X1)>0 and V(X1, X2)>0 is considered. When V(X1|X2)=0, V(X2|X1)=0, and V(X1, X2)>0, the (n, ϵ)-rate region custom-character*(n, ϵ) satisfies






custom-character
in
(3)(n,ϵ)⊆custom-character*(n,ϵ)⊆custom-characterout(3)(n,ϵ).  (85)


For any n and ϵ the achievable region custom-characterin(3)(n, ϵ) has a curved boundary that is characterized by the trade-off between a separate bound on R1 and a region in custom-character2 that bounds (R2, R1+R2) jointly. The converse region custom-characterout(3)(n, ϵ) is the intersection of a region with a linear boundary that bounds R1 only and a region with a curved boundary that bounds (R2, R1+R2) jointly. Letting










δ
=

1

n



,




(
86
)







then it is apparent that the dispersion corresponding to R2 and R1+R2 is given by V2 with a






-


log





n


2

n






third-order term, while the dispersion of R1 is zero.


A less redundant stationary memoryless source has some useful properties. When V(X1, X2)=0,





PX1nX2n(x1n,x2n)∈{0, exp(−nH(X1,X2))},  (87)


for every (x1n, x2n)∈X1n×X2n; in other words, (X1, X2) is uniformly distributed over its support in X1×X2. When V(X1|X2)=0,










P


X
1
n





X
2
n




(



x
1
n





x
2
n

)


=

{




exp
(



-
n







H


(


X
1





X
2

)


)



,








if







P


X
1
n



X
2
n



(


x
1
n

,

x
2
n


)


>
0






0
,

otherwise
;











(
88
)







in other words, X1 is uniformly distributed over its conditional support for each x2∈X2. When V(X2|X1)=0, a result analogous to (88) holds. These properties do not reduce the difficulty of characterizing the optimal MASCs in general. As a result, the random coding techniques described above are utilized. For the achievability argument, the MASC bound described above can be invoked.


Section B: Two Special Cases

The analysis above applies to any stationary, memoryless source with single-letter distribution PX1X2 that satisfies (28). In such a general setting, it is hard to find an optimal code. However, there are some special cases in which the optimal codes for a less redundant source can be characterized.


To enable the following analysis on these special cases, it is assumed that PX1X2(x1, x2)>0 for every (x1, x2)∈X1×X2. Under this assumption, V(X1, X2)=0 if and only if V(X1|X2)=V(X2|X1)=0. As a result, the three cases discussed in Section A reduce to only two possible scenarios:


1) V(X1, X2)=V(X1|X2)=V(X2|X1)=0;


2) V(X1, X2)>0, and either V(X1|X2)=0 or V(X2|X1)=0.


Note that X1 and X2 are independent in both of these scenarios. The results are summarized below.


Special Case 1): Suppose that V(X1|X2)=0, V(X2|X1)=0, and V(X1, X2)=0. If PX1X2 satisfies PX1X2(x1, x2)>0 for every (x1, x2)∈X1×X2, then






custom-character*(n,ϵ)=custom-characterout(1)(n,ϵ),  (89)


where custom-characterout(1)(n, ϵ) is defined in (81).


This scenario is a special example of Case 1) discussed in Section A above. The (n, ϵ)-rate region here coincides with the converse region custom-characterout(1)(n, ϵ) presented in (81) for general source distributions. See FIGS. 11A and 11B for a comparison among custom-characterin(1)(n, ϵ), custom-characterout(1)(n, ϵ), and custom-character*(n, ϵ) in this special case. FIGS. 11A and 11B are schematic illustrations of the MASC rate regions for a less redundant source. The drawing in FIG. 11A illustrates the achievable and converse regions in Case 1) (V(X1, X2)=V(X1|X2)=V(X2|X1)=0) and the (n, ϵ)-rate region custom-character*(n, ϵ) when PX1X2 is assumed to have no zeros (Special Case 1)). The drawing in FIG. 11B illustrates the achievable and converse regions in Case 3) (V(X1|X2)=0, V(X1, X2), V(X2|X1)>0) and the (n, ϵ)-rate region custom-character*(n, ϵ) when PX1X2 is assumed to have no zeros (Special Case 2)).


Special Case 2): With V(X1, X2)>0, it can be assumed that V(X1|X2)=0 and V(X2|X1)>0. The other case can be analyzed similarly. For any δ∈[0, ϵ), the following can be defined
















in
s



(

n
,
δ

)




=
Δ



{


(


R
1

,

R
2


)





2

:










R
1




H


(

X
1

)


-


1
n


log


1

1
-
δ












R
2






H


(

X
2

)


+




V


(

X
2

)


n





Q

-
1




(


ϵ
-
δ


1
-
δ


)



-


log





n


2

n


+


ξ
in



(

ϵ
,
δ
,
n

)






}






(
90
)














out
s



(

n
,
δ

)




=
Δ



{


(


R
1

,

R
2


)





2

:










R
1




H


(

X
1

)


-


1
n


log


1

1
-
δ












R
2






H


(

X
2

)


+




V


(

X
2

)


n





Q

-
1




(


ϵ
-
δ


1
-
δ


)



-


log





n


2

n


-


ξ
out



(

ϵ
,
δ
,
n

)






}


,





(
91
)







where the functions ξin(ϵ, δ, n) and ξout(ϵ, δ, n) are characterized as follows. For any fixed δ,








ξ
out



(

ϵ
,
δ
,
n

)


=



O


(

1
n

)







and







ξ
in



(

ϵ
,
δ
,
n

)



=


O


(

1
n

)


.






For any fixed n, both ξout(ϵ, δ, n) and ξin(ϵ, δ, n) blow up as δ approaches ϵ.


Also the following can be defined
















in
s



(

n
,
ϵ

)




=
Δ






δ




0
,
t

)










in
s



(

n
,
δ

)








(
92
)








out
s



(

n
,
ϵ

)




=
Δ






δ




0
,
t

)










out
s




(

n
,
δ

)

.







(
93
)







Suppose that V(X1|X2)=0, V(X2|X1)>0, and V(X1, X2)>0. If PX1X2 satisfies PX1X2(x1, x2)>0 for every (x1, x2)∈X1×X2, then






custom-character
in
s(n,ϵ)⊆custom-character*(n,ϵ)⊆custom-characterouts(n,ϵ).  (94)


This scenario is a special example of Case 3) discussed in Section A above. The achievable region is

















in

(
3
)




(

n
,
ϵ

)


=




δ


(

0
,
ϵ

)





{


(


R
1

,

R
2


)





2

:










R
1




H


(

X
1

)


+


1
n


log


1
δ







R
2






H


(

X
2

)


+




V


(

X
2

)


n





Q

-
1




(

ϵ
-
δ
-


C
in


n



)



-


log





n


2

n


+


1
n


log


1

1
-
δ







}



,





(
95
)







and the converse region is



















out

(
3
)




(

n
,
ϵ

)


=



{


(


R
1

,

R
2


)





2

:











R
1






H


(

X
1

)


+


1
n


log


1

1
-
ϵ













R
2






H


(

X
2

)


+




V


(

X
2

)


n





Q

-
1




(
ϵ
)



-


log





n


2

n


-

O


(

1
n

)




}

.







(
96
)







As δ approaches ϵ, the boundary of the (n, ϵ)-rate region given in (92) approaches the line








R
1

=


H


(

X
1

)


-


1
n


log


1

1
-
ϵ





,




which matches the vertical segment of the boundary of the converse region custom-characterout(3)(n, ϵ). See FIG. 11B for a comparison of custom-characterin(3)(n, ϵ), custom-characterout(3)(n, ϵ), and custom-character*(n, ϵ) in this case.


Although the present invention has been described in certain specific aspects, many additional modifications and variations would be apparent to those skilled in the art. It is therefore to be understood that the present invention can be practiced otherwise than specifically described including using any of a variety of different decoders and/or sources without departing from the scope and spirit of the present invention. Thus, embodiments of the present invention should be considered in all respects as illustrative and not restrictive. Accordingly, the scope of the invention should be determined not by the embodiments illustrated, but by the appended claims and their equivalents.

Claims
  • 1. A random access source coding (RASC) system, comprising: a plurality of source encoders that are each configured to receive data from one of a plurality of sources;when one of the plurality of source encoders receives data from one of the plurality of sources, then the source encoder is configured to: receive a start of epoch message; andtransmit a portion of a codeword selected by encoding data from one of the plurality of sources until an end of epoch feedback message is received; anda receiver comprising a source decoder, where the source decoder is configured to: cause a broadcast transmitter to transmit at least one start of epoch message;receive at least one portion of a codeword transmitted by at least one of the plurality of source encoders; andwhen a decoding rule is satisfied, decode data from the at least one of the plurality of the source encoders based upon the received at least one portion of a codeword, and cause the broadcast transmitter to transmit an end of epoch message.
  • 2. The RASC system of claim 1, wherein: the source decoder is configured to determine whether the decoding rule is satisfied at each of a predetermined set of decode times and transmit a feedback message;each of the plurality of source encoders is configured to receive feedback messages at the predetermined set of potential decoding times; andeach of the predetermined set of decode times corresponds to one of a plurality of blocklengths of the codeword selected by encoding data from one of the plurality of sources.
  • 3. The RASC system of claim 2, wherein: a first of the plurality of blocklengths of the codeword selected by encoding data from one of the plurality of sources forms a first blocklength code;a second of the plurality of blocklengths of the codeword selected by encoding data from one of the plurality of sources forms a second blocklength code; andthe first blocklength code forms a prefix for the second blocklength code.
  • 4. The RASC system of claim 3, wherein: a third of the plurality of blocklengths of the codeword selected by encoding data from one of the plurality of sources forms a third blocklength code; andthe second blocklength code forms a prefix for the third blocklength code.
  • 5. The RASC system of claim 1, wherein each source encoder is configured to select codewords from a codebook.
  • 6. The RASC system of claim 1, wherein the data received from the plurality of sources is dependent.
  • 7. The RASC system of claim 1, wherein the broadcast transmitter sends a negative acknowledgement feedback message when the decoding rule is not satisfied at a predetermined decode time.
  • 8. The RASC system of claim 1, wherein each source encoder that is transmitting is configured to transmit an additional portion of the codeword selected by encoding data from one of the plurality of sources in response to receipt of a negative acknowledgement feedback message.
  • 9. A source encoding system, comprising: a receiver configured to receive a start of epoch message; anda source encoder configured to transmit a portion of a codeword selected by encoding data from a source until an end of epoch feedback message is received.
  • 10. The source encoding system of claim 9, wherein: each of the plurality of source encoders is configured to receive feedback messages at each of a predetermined set of decode times; andeach of the predetermined set of decode times corresponds to one of a plurality of blocklengths of the codeword selected by encoding data from the source.
  • 11. The source encoding system of claim 10, wherein: a first of the plurality blocklengths of the codeword selected by encoding data from the source forms a first blocklength code;a second of the plurality of blocklengths of the codeword selected by encoding data from the source forms a second blocklength code; andthe first blocklength code forms a prefix for the second blocklength code.
  • 12. The source encoding system of claim 11, wherein: a third of the plurality of blocklengths of the codeword selected by encoding data from the source forms a third blocklength code; andthe second blocklength code forms a prefix for the third blocklength code.
  • 13. The source encoding system of claim 9, wherein the source encoder is configured to transmit an additional portion of the codeword selected by encoding data from the source in response to receipt of a negative acknowledgement feedback message.
  • 14. The source encoding system of claim 9, wherein the data received from the plurality of sources is dependent.
  • 15. A receiver system, comprising: a source decoder configured using a codebook;a transmitter configured to transmit at least one start of epoch message; anda receiver configured to receive at least one portion of a codeword;wherein the source decoder is further configured to determine when a decoding rule is satisfied, and when the decoding rule is satisfied: decode at least one message based upon the received at least one portion of a codeword using the codebook; andcause the transmitter to transmit an end of epoch message.
  • 16. The receiver system of claim 15, wherein: the source decoder is configured to determine whether the decoding rule is satisfied at each of a predetermined set of decode times and transmit a feedback message; andeach of the predetermined set of decode times corresponds to one of a plurality of blocklengths of the codeword.
  • 17. The receiver system of claim 16, wherein: a first of the plurality of blocklengths of the codeword forms a first blocklength code;a second of the plurality of blocklengths of the codeword forms a second blocklength code; andthe first blocklength code forms a prefix for the second blocklength code.
  • 18. The receiver system of claim 17, wherein: a third of the plurality of blocklengths of the codeword forms a third blocklength code; andthe second blocklength code forms a prefix for the third blocklength code.
  • 19. The source encoding system of claim 17, wherein the data received from the plurality of sources is dependent.
  • 20. The receiver system of claim 15, wherein the broadcast transmitter sends a negative acknowledgement feedback message when the decoding rule is not satisfied at a predetermined decode time.
STATEMENT OF FEDERALLY SPONSORED RESEARCH

This invention was made with government support under Grant No. CCF1817241 awarded by the National Science Foundation. The government has certain rights in the invention. The present invention claims priority to U.S. Provisional Patent Application Ser. No. 62/963,693 entitled “Lossless Source Coding In The Point-to-Point, Multiple Access, And Random Access Regimes” to Chen et al., filed Jan. 21, 2020, the disclosures of which is herein incorporated by reference in its entirety.

Provisional Applications (1)
Number Date Country
62963693 Jan 2020 US