METHOD AND DEVICE FOR MULTILEVEL POLAR-CODED MODULATION TRANSMITTING AND RECEIVING

Abstract

The present disclosure relates to a method (50) for transmitting a data stream from a transmitting device (30) to a receiving device (40), the data stream comprising m data sub-streams, said method (50) comprising: encoding (S51) the data sub-streams with respective polar codes, such as to produce m polar-encoded data sub-streams, modulating (S52) the m polar-encoded data sub-streams onto symbols of a multilevel modulation comprising m levels defining 2.sup.m different symbols, such as to produce a symbol stream, transmitting (S53) the symbol stream to the receiving device (40), wherein the 2.sup.m symbols of the multilevel modulation are distributed in the complex plane such that they are regularly spaced along the real axis and the complex axis by a spacing factor K, the predetermined labeling function Z being such that Φ (Z)<1.7.sup.m×(0.0425×m−0.06), expression in which:

Claims

1. A method for transmitting a data stream from a transmitting device to a receiving device, the data stream comprising m data sub-streams, m≥4 said method comprising: encoding the data sub-streams with respective polar codes, such as to produce m polar-encoded data sub-streams, modulating the m polar-encoded data sub-streams onto symbols of a multilevel modulation comprising m levels defining 2.sup.m different symbols in a complex plane, according to a predetermined labeling function Z which bijectively associates a symbol to a set c.sub.1, . . . , c.sub.m of m bits from respective polar-encoded data sub-streams, such that successive sets of m bits from respective polar-encoded data sub-streams are converted into successive symbols forming a symbol stream, transmitting the symbol stream to the receiving device, wherein the 2.sup.m symbols of the multilevel modulation are distributed in the complex plane such that they are regularly spaced along the real axis and the complex axis by a spacing factor K, wherein the predetermined labeling function Z is such that:
Φ(Z)<1.7.sup.m×(0.0425×m−0.06) expression in which: $\Phi \left(Z\right)=\underset{i=1}{\overset{m}{.Math.}}\underset{\left\{c_1,.Math.,c_m\right\}}{.Math.}\underset{\left\{c_i^{\prime },.Math.,c_m^{\prime }\right\}\text{(}{c}_i,.Math.,c_m)}{.Math.}{e}^{-\frac{4}{K^2}{.Math.Z\left(c_1,.Math.,c_m\right)-Z\left(c_1,.Math.,c_{i-1},c_i^{\prime },.Math.,c_m^{\prime }\right).Math.}^2}$

2. The method according to claim 1, wherein each of the m polar codes used is a polar code designed for a binary input additive white Gaussian noise, AWGN, channel.

3. The method according to claim 1, wherein the respective code rates of the m polar codes are adapted to propagation conditions.

4. The method according to claim 3, comprising: determining signal to noise ratios, SNRs, experienced by each data sub-stream during the propagation from the transmitting device to the receiving device, selecting the code rates of each of the m polar codes used based on the SNRs experienced by each data sub-stream.

5. The method according to claim 1, wherein the data sub-stream and/or each of the data sub-streams include an error detection code.

6. The method according to claim 1, wherein the multilevel modulation is a 16QAM modulation, and the labeling function is given by the following table, or by the following table modified by rotating the constellation of symbols in the complex plane by an angle multiple of 90°, and/or by applying an axial reflection to the constellation of symbols with respect to the real axis, and/or by applying an axial reflection to the constellation of symbols with respect to the complex axis, and/or by inverting all the bit values: TABLE-US-00005 (c.sub.1, . . . , c.sub.m) $Z\left(c_1,\mathrm{.Math.},c_m\right)\times \frac{2}{K}$ 0000 +3−1j 1000 +3−3j 0100 −3−3j 1100 −3−1j 0010 −1−1j 1010 −1−3j 0110 +1−3j 1110 +1−1j 0001 −3+1j 1001 −1+1j 0101 +1+1j 1101 +3+1j 0011 +3+3j 1011 +1+3j 0111 −1+3j 1111  −3+3j.

7. A computer program product comprising instructions which, when executed by a processor, configure said processor to carry out a transmitting method according to claim 1.

8. A computer-readable storage medium comprising instructions which, when executed by a processor, configure said processor to carry out a transmitting method according to claim 1.

9. A device for transmitting a data stream to a receiving device, comprising a processing circuit configured to carry out a transmitting method according to claim 1.

10. A method for receiving a data stream by a receiving device, said data stream comprising m data sub-streams, m≥4, comprising: receiving a symbol stream including the data stream from a transmitting device, the symbol stream including symbols of a multilevel modulation comprising m levels defining 2.sup.m different symbols in a complex plane, demodulating the symbols of the symbol stream in order to produce m polar-encoded data sub-streams, by applying an inverse function of a predetermined labeling function Z which bijectively associates a symbol to a set c.sub.1, . . . , c.sub.m of m bits from respective polar-encoded data sub-streams, decoding successively the m polar-encoded data sub-streams such as to produce m data sub-streams, by applying respective polar-code decoders, wherein the 2.sup.m symbols of the multilevel modulation used are distributed in the complex plane such that they are regularly spaced along the real axis and the complex axis by a spacing factor K, wherein the predetermined labeling function Z is such that:
Φ(Z)<1.7.sup.m×(0.0425×m−0.06) expression in which: $\Phi \left(Z\right)=\underset{i=1}{\overset{m}{.Math.}}\underset{\left\{c_1,.Math.,c_m\right\}}{.Math.}\underset{\left\{c_i^{\prime },.Math.,c_m^{\prime }\right\}\text{(}{c}_i,.Math.,c_m)}{.Math.}{e}^{-\frac{4}{K^2}{.Math.Z\left(c_1,.Math.,c_m\right)-Z\left(c_1,.Math.,c_{i-1},c_i^{\prime },.Math.,c_m^{\prime }\right).Math.}^2}$

11. The method according to claim 10, wherein each of the m polar code decoders used is designed for decoding a polar code designed for a binary input additive white Gaussian noise, AWGN, channel.

12. The method according to claim 10, wherein each polar-code decoder used is a list decoder.

13. A computer program product comprising instructions which, when executed by a processor, configure said processor to carry out a receiving method according to claim 10.

14. A computer-readable storage medium comprising instructions which, when executed by a processor, configure said processor to carry out a receiving method according to claim 10.

15. A device for receiving a data stream from a transmitting device, comprising a processing circuit configured to carry out a receiving method according to claim 10.

Description

BRIEF DESCRIPTION OF DRAWINGS

[0049] FIG. 1 is a schematic representation of an exemplary multilevel polar-coded modulation scheme.

[0050] FIG. 2 is a plot representing the mutual information for a 16QAM modulation using a Gray labeling function, and the distribution of the mutual information between the 4 bit levels.

[0051] FIG. 3 is a plot representing the mutual information for a 16QAM modulation using a random labeling function, and the distribution of the mutual information between the 4 bit levels.

[0052] FIG. 4 is a plot representing the log likelihood ratio distribution at the input of the different polar code decoders of a multilevel polar-coded modulation scheme, when using a Gray labeling function.

[0053] FIG. 5 is a plot representing a preferred embodiment of a 16QAM labeling function.

[0054] FIG. 6 is a plot representing a preferred embodiment of a 64QAM labeling function.

[0055] FIG. 7 is a plot representing a preferred embodiment of a 256QAM labeling function.

[0056] FIG. 8 is a plot representing the log likelihood ratio distribution at the input of the different polar code decoders of a multilevel polar-coded modulation scheme, when using the 16QAM labeling function of FIG. 5.

[0057] FIG. 9 is a plot comparing the expected throughput performance for the 16QAM Gray labeling function and for the 16QAM labeling function of FIG. 5.

[0058] FIG. 10 is a diagram representing the main steps of a transmitting method according to a preferred embodiment.

[0059] FIG. 11 is a schematic representation of a transmitting device according to a preferred embodiment.

[0060] FIG. 12 is a diagram representing the main steps of a receiving method according to a preferred embodiment.

[0061] FIG. 13 is a schematic representation of a receiving device according to a preferred embodiment.

DESCRIPTION OF EMBODIMENTS

[0062] In these figures, references identical from one figure to another designate identical or analogous elements. For reasons of clarity, the elements shown are not to scale, unless explicitly stated otherwise.

[0063] As indicated above, the present disclosure relates to a multilevel polar-coded modulation scheme.

[0064] FIG. 1 represents schematically a multilevel polar-coded modulation scheme. It is assumed, in a non-limiting manner, that the multilevel modulation comprises m levels, implying that the corresponding constellation comprises 2.sup.m symbols in the complex plane. Preferably, m is higher than or equal to 4 (m≥4). Preferably, m is a power of 2, and is for instance equal to 4, equal to 6, equal to 8, etc.

[0065] As illustrated by FIG. 1, each level comprises a polar code encoder 20, such that m polar code encoders are used in parallel on the m different levels of the multilevel modulation. The m polar code encoders can be either identical or different. In preferred optional embodiments, that will be discussed hereinafter, the respective code rates of the m polar code encoders can be dynamically adapted to the propagation conditions, potentially resulting in the polar code encoders using different code rates.

[0066] As illustrated by FIG. 1, b.sub.i represents data bits to be encoded by the polar code encoder 20 of the i-th level, that produces polar-encoded data bits c.sub.i. Each level of polar-encoded data bits is then fed to a symbol modulator 21 that produces a stream of symbols. The constellation modulator takes one polar-encoded data bit per level, i.e. a total of m bits (c.sub.1, . . . , c.sub.m), and outputs one symbol X=Z(c.sub.1, . . . , c.sub.m) according to a predetermined labeling function Z. A symbol X is defined by a set of coordinates in real and complex dimensions and, as discussed above, the constellation of symbols comprises 2.sup.m different symbols in the complex plane. The labeling function Z(c.sub.1, . . . , c.sub.m) allows to make a one-to-one mapping between a set of m encoded data bits (c.sub.1, . . . , c.sub.m) and a symbol X=Z(c.sub.1, . . . , c.sub.m). Successive sets of m polar-encoded data bits (c.sub.1, . . . , c.sub.m) thus produce successive symbols, forming a symbol stream, i.e. a vector of symbols transmitted on respective resources.

[0067] In the sequel, c.sub.i denotes one polar-encoded data bit of the i-th level. If necessary to distinguish different polar-encoded data bits of a same level, a second subscript is used such that c.sub.i,t denotes the polar-encoded data bit of the i-th level that is used to produce the t-th symbol transmitted on the t-th resource. Also, b.sub.i denotes one data bit of the i-th level, and b.sub.i corresponds to a vector comprising several data bits of the i-th level. Similarly, X and Y denote respectively one symbol and one observation. If necessary to distinguish different symbols and different observations, a subscript is used such that X.sub.t and Y.sub.t denote respectively the t-th symbol and the t-th observation, transmitted and observed on the t-th resource. X and Y denote respectively a vector comprising T different symbols X.sub.t (X=(X.sub.1, . . . , X.sub.T)) and a vector comprising T different observations Y.sub.t (Y=(Y.sub.1, . . . , Y.sub.T)). It is noted that the T resources can correspond to different time slots, to different frequencies (e.g. different subcarriers of an orthogonal frequency division multiplex, OFDM, symbol), different spreading codes, etc., or any combination thereof. Each symbol X.sub.t is obtained based on a set of m polar-encoded data bits (c.sub.1,t, . . . , c.sub.m,t). The vector (c.sub.i,1, . . . , c.sub.i,T) comprising the T polar-encoded data bits for the i-th level is denoted c.sub.i. Each vector c.sub.i of polar-encoded data bits is obtained by encoding the vector b.sub.i of data bits with a polar code. Each vector b.sub.i comprises a number of data bits that depends on the code rate of polar code encoder of the i-th level, and that can vary from one level to another.

[0068] The vector X of symbols is then transmitted by a transmitting device 30 through a channel and an observation Y of the vector X of symbols is obtained at a receiving device 40. The channel modifies each symbol X=Z(c.sub.1, . . . , c.sub.m) into the observation Y according to a statistical behavior characterized by the transition probability p(Y|c.sub.1, . . . , c.sub.m) of the channel. For example, when considering an additive white Gaussian noise, AWGN, channel, the observation Y is the result of Z(c.sub.1, . . . , c.sub.m) plus a noise following a Gaussian distribution.

[0069] In order to properly design a multilevel polar-coded modulation, it is of interest to first characterize the model from an information theory standpoint. In general, the channel capacity is by definition the maximization of the channel mutual information I(X; Y), i.e. without taking into account any structure for the transmitting device 30 or the receiving device 40. When adding elements in the encoding/decoding structure, the mutual information decreases. For example, in an AWGN channel, by constraining the symbols to belong to a finite constellation, the mutual information I(X; Y) is lower than with a continuous Gaussian input which maximizes the channel mutual information.

[0070] Further to constraining each symbol X to belong to a finite constellation, we constrain the system to a multilevel modulation structure. The one-to-one mapping between X=Z(c.sub.1, . . . , c.sub.m) and (c.sub.1, . . . , c.sub.m) allows to now consider I(c.sub.1, . . . , c.sub.m; Y) without any loss. By using the chain rule, we obtain:

I(X;Y)=I(c.sub.1, . . . ,c.sub.m;Y)=Σ.sub.i=1.sup.mI(c.sub.i;Y|c.sub.1, . . . ,c.sub.i−1)

expression in which I(c.sub.i; Y|c.sub.1, . . . , c.sub.i−1) corresponds to the conditional mutual information.

[0071] We can observe that by associating a code rate I(c.sub.i; Y|c.sub.1, . . . , c.sub.i−1) to the i-th polar-coded level, no loss of information is observed if successive cancellation decoding is performed at the receiving device (which is suggested by the knowledge of c.sub.1, . . . , c.sub.i−1 when decoding c.sub.i). By definition of the mutual information, no error occurs if the code rate is adapted to the equivalent level channel. An equivalent level channel corresponds to a virtual channel experienced by a single polar-encoded data sub-stream, the m polar-encoded data sub-streams being transmitted on m parallel equivalent level channels. In practice, when using a finite length coding scheme, errors occur and propagate from one level to another when using a successive decoding scheme. This error propagation should be as limited as possible, by design.

[0072] We assume that the bits c.sub.i are equiprobable, i.e., that the probability p(c.sub.i=0)=p(c.sub.i=1)=0.5, such that the conditional mutual information can be expressed as:

[00003]$I\left(c_i;Y|c_1,.Math.,c_{i-1}\right)=E_{c_1,.Math.,c_m}{E}_{Y.Math.c_1,.Math.,c_m}\left[\text{⁠}1-{\log }_2\left(1+\frac{{.Math.}_{\left\{c_{i+1}^{\prime },.Math.,c_m^{\prime }\right\}}p\left(Y.Math.c_1,.Math.,c_{i-1},1-c_i,c_{i+1}^{\prime },.Math.,c_m^{\prime }\right)}{{.Math.}_{\left\{c_{i+1}^{\prime },.Math.,c_m^{\prime }\right\}}p\left(Y.Math.c_1,.Math.,c_{i-1},c_i,c_{i+1}^{\prime },.Math.,c_m^{\prime }\right)}\right)\right]$

expression in which:

[0073] E.sub.c.sub.1.sub., . . . , c.sub.m is an expectation over all possible polar-encoded data bits {c.sub.1, . . . , c.sub.m}, or equivalently on all possible constellation symbols {X};

[0074] E.sub.Y|c.sub.1.sub., . . . , c.sub.m is an expectation over all possible observations Y when the transmitted symbol Z(c.sub.1, . . . , c.sub.m) is known; it is fully characterized by the transition probability p(Y|Z(c.sub.1, . . . , c.sub.m));

[0075] the sum Σ.sub.{c′.sub.1+1.sub., . . . , c′.sub.m} is performed over all the 2.sup.m−i sets of polar-encoded data bits {c′.sub.i+1, . . . , c′.sub.m}; this marginalization makes appear the transition probabilities p(Y|Z(c.sub.1, . . . , c.sub.i−1, c.sub.i, c′.sub.i+1, . . . , c′.sub.m)) or p (Y|Z(c.sub.1, . . . , c.sub.i−1, 1−c.sub.i, c′.sub.i+1, . . . , c′.sub.m));

[0076] Z(c.sub.1, . . . , c.sub.i−1, c.sub.i, c′.sub.i+1, . . . , c′.sub.m) is the symbol with the i first polar-encoded data bits in their labeling fixed to (c.sub.1, . . . , c.sub.i), the same values as the i first polar-encoded data bits of the labeling of Z(c.sub.1, . . . , c.sub.m) selected in the expectation E.sub.c.sub.1.sub., . . . , c.sub.m; and the last (m−i) polar-encoded data bits of the labeling of Z(c.sub.1, . . . , c.sub.i−1, c.sub.i, c′.sub.i+1, . . . , c′.sub.m) are set by the configuration (c′.sub.i+1, . . . , c′.sub.m) selected in the sum Σ.sub.{c′.sub.i+1.sub., . . . , c′.sub.m.sub.}; thus, Z(c.sub.1, . . . , c.sub.i−1, c.sub.i, c′.sub.i+1, . . . , c′.sub.m) is one of the multilevel modulation symbols, and p(Y|Z(c.sub.1, . . . , c.sub.i−1, c.sub.i, c′.sub.i+1, . . . , c′.sub.m)) can be easily computed by the knowledge of the observation Y and of the channel transition probability;

[0077] similarly, Z(c.sub.1, . . . , c.sub.i−1, 1−c.sub.i, c′.sub.i+1, . . . , c′.sub.m) is the constellation symbol with labeling having the (i−1) first polar-encoded data bits equal to (c.sub.1, . . . , c.sub.i−1), the i-th polar-encoded data bit value is (1−c.sub.i), and the last (m−i) polar-encoded data bits values are (c′.sub.i+1, . . . , c′.sub.m).

[0078] The chain rule decomposition of the Gray labeling for a 16QAM (16-level QAM) is shown in FIG. 2. The sum of the 4 plots associated each to one of the 4 bit levels is indeed equal to the mutual information of the 16QAM constellation without constraint on the polar-coding scheme. We can observe that, since the Gray labeling is independent and the same on the real and imaginary parts, the mutual information of bits 1 and 3 (i.e. of the levels 1 and 3 of the multilevel modulation) and bits 2 and 4 (i.e. of the levels 2 and 4) are equal. The chain rule decomposition of a random labeling for a 16QAM is shown in FIG. 3. We can observe that, for a given signal to noise ratio, SNR, the rate splitting is different from the one of Gray labeling. Yet, the sum of the rates is still equal to the mutual information of the 16QAM constellation.

[0079] Hence, it is important to note that, from the mutual information perspective, only the vector X of symbols has an impact on the performance. The labeling function itself does not improve or degrade the overall theoretical performance. It just changes the repartition of the rate between the parallel levels of the multilevel polar-coded modulation. This remark is important as it highlights that any labeling function performs theoretically the same in view of information theory.

[0080] However, we will see later that, when combined with a practical polar code, this remark is not true anymore. In other words, when considering practical polar codes in a multilevel polar-coded modulation context, the labeling can influence the overall actual performance, such that optimization of the labeling function can be made by taking into account the polar coding strategy of each parallel level of the multilevel polar-coded modulation.

[0081] In the sequel, we assume that the channel is memory-less. Hence, the channel effect is independent from one transmission resource to another (e.g., additive white Gaussian noise is usually independent in time).

[0082] Thus, the likelihood p(Y|b.sub.1, . . . , b.sub.m) can be decomposed through the T resources used for transmitting the symbol vector X=(X.sub.1, . . . , X.sub.T), as follows:

[00004]$p\left(\stackrel{\_}{Y}|{\overline{b}}_1,.Math.,{\overline{b}}_m\right)=\underset{t=1}{\overset{T}{.Math.}}p\left(Y_t|{\overline{b}}_1,.Math.,{\overline{b}}_m\right)$

[0083] Then, by using the fact that the encoder converts the vector of data bits b.sub.i into a vector of polar-encoded data bits c.sub.i=(c.sub.i,1, . . . , c.sub.i,T), we can write:

[00005]$p\left(\overline{Y}|{\overline{b}}_1,.Math.,{\overline{b}}_m\right)=\underset{t=1}{\overset{T}{.Math.}}p\left(Y_t|c_{1,t},.Math.,c_{m,t}\right)$

[0084] This involves that the relationship between Y and b.sub.1, . . . , b.sub.m can be analyzed resource per resource by working on Y.sub.t and (c.sub.1,t, . . . , c.sub.m,t), and the associated symbol X.sub.t=Z(c.sub.1,t, . . . , c.sub.m,t). It is recalled that the role of the labeling function Z is to determine which symbol is uniquely associated with a polar-encoded data bit vector (c.sub.1,t, . . . , c.sub.m,t). This allows reducing the complexity of the decoder by considering a symbol-to-bit demodulation (i.e. inverse function of the labeling function Z) performed independently on each resource.

[0085] It is assumed that the decoder successively decodes the first level (i.e. the level with index 1), then the second level (i.e. the level with index 2), etc., until the last level (i.e. the level with index m). For instance, the decoder of the multilevel polar-coded modulation with list decoding can typically involve the following steps, performed successively for each level, from the level with index 1 to the level with index m, which include for the i-th level:

[0086] obtain the data bit vectors b.sub.1, . . . , b.sub.i−1 associated to each surviving path of the list decoding procedure;

[0087] re-encode and obtain the polar-encoded data bits according to the polar code encoder of each level (this can also be done during the decoding process of b.sub.1, . . . , b.sub.i−1);

[0088] for each transmission resource with index t with observation Y.sub.t, compute the likelihoods p(Y.sub.t|c.sub.1,t, . . . , c.sub.i−1,t, 1, c′.sub.i+1, . . . , c′.sub.m) and p(Y.sub.t|c.sub.1,t, . . . , c.sub.i−1,t, 0, c′.sub.i+1, . . . , c′.sub.m) for the obtained re-encoded data bits (c.sub.1,t, . . . , c.sub.i−1,t) and for all the 2.sup.m−i possible configurations {c′.sub.i+1, . . . , c′.sub.m},

[0089] compute the following log-likelihood ratios, for all indexes t:

[00006]$LL{R}_{i,t}=\log \left(\frac{{.Math.}_{\left\{c_{i+1}^{\prime },.Math.,c_m^{\prime }\right\}}p\left(Y_t.Math.c_{1,t},.Math.,c_{i-1,t},0,c_{i+1}^{\prime },.Math.,c_m^{\prime }\right)}{{.Math.}_{\left\{c_{i+1}^{″},.Math.,c_m^{″}\right\}}p\left(Y_t.Math.c_{1,t},.Math.,c_{i-1,t},1,c_{i+1}^{″},.Math.,c_m^{″}\right)}\right)$

expression in which Σ.sub.{c′.sub.i+1.sub., . . . , c′.sub.m.sub.} corresponds to the sum over the 2.sup.m−i possible configurations {c′.sub.i+1, . . . , c′.sub.m}, and Σ.sub.{c″.sub.i+1.sub., . . . , c″.sub.m.sub.} corresponds to the sum over the 2.sup.m−i possible configurations {c″.sub.i+1, . . . , c″.sub.m};

[0090] feed the input of the decoder of the i-th level with the last cumulated metric resulting from the decoding of the (i−1)-th level for each list decoding survivor, and perform the decoding to obtain the data bit vector b.sub.i for each surviving path and the associated cumulated metric.

[0091] It is noted that the above steps are only given for exemplary purposes, and that the decoding can be performed differently. In particular, it is possible, in other embodiments, to apply a simple successive decoding scheme without list decoding, i.e. by computing only a single survivor for each level. More generally speaking, the present disclosure is not limited to a specific decoding process, and the choice of a specific decoding process merely constitutes an exemplary embodiment of the present disclosure.

[0092] However, the present disclosure relies on the assumption that the different levels are decoded successively at the decoder, with the level with index 1 being decoded first and the level with index m being decoded last.

[0093] Also, the above exemplary embodiment of a successive decoding scheme has been given in order to highlight the structure of the log-likelihood ratios LLR.sub.i,t, which is not expected to vary significantly from a successive decoding scheme to another, and which help understanding why the labeling function can be optimized in order to increase the overall actual performance of the system, when considering practical implementations.

[0094] As previously said, practical implementations of polar code encoders and decoders correspond to polar codes designed for binary input AWGN channels. Such polar codes expect to be used on an AWGN channel with binary input, and consequently expect to have log likelihood ratios LLR.sub.i,t at the input of the polar code decoder that follow a Gaussian distribution.

[0095] However, when using a high order multilevel coded modulation (m≥4) with a Gray labeling, at least one of the m polar code decoders experiences an input log likelihood ratio distribution which is not Gaussian, such as illustrated in FIG. 4. This leads to degraded performance.

[0096] The inventors have noticed that the labeling function influences the input log likelihood ratio distribution of the different levels, and that some labeling functions outperform the conventional Gray labeling when considering a multilevel polar-coded modulation using an off-the-shelf polar codes designed for binary input AWGN channels.

[0097] In particular, the inventors have discovered that it is possible to find labeling functions which provide input log likelihood ratio distributions that are substantially Gaussian for all the levels of the multilevel polar-coded modulation. Hence, despite not influencing the mutual information, the labeling function does influence the overall actual performance if it enables providing the inputs of the different polar code decoders with log likelihood ratios that tend to follow the expected Gaussian distributions.

[0098] For that purpose, the inventors have discovered that the following function could be used to identify suitable labeling functions:

[00007]$\begin{array}{cc}\Phi \left(Z\right)={.Math.}_{i=1}^m{.Math.}_{\left\{c_1,.Math.,c_m\right\}}{.Math.}_{\left\{c_i^{\prime },.Math.,c_m^{\prime }\right\}\setminus \left(c_i,.Math.,c_m\right)}{e}^{-\frac{4}{K^2}{.Math.Z\left(c_1,.Math.,c_m\right)-Z\left(c_1,.Math.,c_{i-1},c_i^{\prime },.Math.,c_m^{\prime }\right).Math.}^2}& \left(1\right)\end{array}$

[0099] In the expression (1) above:

[0100] e.sup.x corresponds to the exponential function applied to the real number x;

[0101] |z| corresponds to the modulus (or magnitude) of the complex number z (i.e. |z|=√{square root over (z.sub.1.sup.2+z.sub.2.sup.2)} if z=z.sub.1+j×z.sub.2, with z.sub.1, z.sub.2 being real numbers and j being the imaginary unit);

[0102] K is a spacing factor of the symbols in the complex plane;

[0103] Σ.sub.{c.sub.1.sub., . . . , c.sub.m.sub.} is a sum over the 2.sup.m possible values of the vector (c.sub.1, . . . , c.sub.m), bijectively associated to the 2.sup.m symbols of the constellation by the labeling function Z;

[0104] Σ.sub.{c.sub.i.sub.′, . . . , c.sub.m.sub.′}\(c.sub.i.sub., . . . , c.sub.m.sub.) is a sum over the (2.sup.m−i+1−1) possible values of the vector (c.sub.i′, . . . , c.sub.m′) when excluding the specific value (c.sub.i, . . . , c.sub.m) as currently set by the sum Σ.sub.{c.sub.1.sub., . . . , c.sub.m}, the value of i being set by the sum Σ.sub.i=1.sup.m; for completeness, it is emphasized that, for i=1, Z(c.sub.1, . . . , c.sub.i−1, c′.sub.i, . . . , c′.sub.m)=Z(c′.sub.1, . . . , c′.sub.m) in said sum and that, for i=2, Z(c.sub.1, . . . , c.sub.i−1, c′.sub.i, . . . , c′.sub.m)=Z(c.sub.1, c′.sub.2, . . . , c′.sub.m) in said sum, etc.

[0105] More specifically, regarding the spacing factor K, defined in expression (1) above, it implies that the 2.sup.m symbols of the constellation of the multilevel modulation are distributed in the complex plane such that they are regularly spaced along the real axis and the complex axis by said spacing factor K. For instance, if the sum of the 2.sup.m symbols is equal to zero (i.e. the 2.sup.m symbols are centered around (0,0) the center of the complex plane), then the 2.sup.m symbols are in the following set:

[00008]$\left\{\left(\pm K\times \left(k-\frac{1}{2}\right),\pm K\times \left(k^{\prime }-\frac{1}{2}\right)\right),k,k^{\prime }\in .Math.\left[1,2^{m/2-1}\right]\right\}.$

expression in which corresponds to the set of natural integers, such that:

[0106] for m=4 (16QAM modulation): k, k′∈{1, 2};

[0107] for m=6 (64QAM modulation): k, k′∈{1, 2, 3, 4};

[0108] for m=8 (256QAM modulation): k, k′∈{1, 2, 3, 4, 5, 6, 7, 8};

[0109] etc.

[0110] In general, the function Φ is a function that produces a value representative of fairness between the elementary distances |Z(c.sub.1, . . . , c.sub.m)−Z(c.sub.1, . . . , c.sub.i−1, c′.sub.i, . . . , c′.sub.m)|.sup.2 playing a role in each log likelihood ratio LLR.sub.i,t production for each level of the multilevel coded modulation. By optimizing the function Φ, i.e. by minimizing the function Φ, it is possible to obtain a labeling function Z for which the log likelihood ratios LLR.sub.i,t will tend to follow a Gaussian distribution for all the levels of the multilevel coded modulation.

[0111] As a remark, this function Φ does not take into account the channel transition probabilities. It helps determining labeling functions mimicking the log likelihood ratio distribution of a binary input AWGN channel regardless the actual channel. This is a key feature of our approach, since the labeling can be used in a context where the channel is unknown at the transmitting device 30.

[0112] As indicated above, this function Φ assumes that the different levels are decoded successively at the receiving device 40, with the level with index 1 being decoded first and the level with index m being decoded last.

[0113] In order to optimize the function Φ, it is for instance possible to execute a genetic-like algorithm, for example by executing the following steps:

[0114] initialize a labeling function Z; for example, it is possible to start from the Gray labeling or any random labeling;

compute BestMetric=Φ(Z);

[0115] repeat the following steps until a stopping criterion is verified:

TABLE-US-00001  for all j = 1 ••• m   for all {c.sub.1, ••• , c.sub.m}    copy Z into Z′    find the symbols with the labeling (c.sub.1, ••• , c.sub.j =   0, ... , c.sub.m) and (c.sub.1, ... , c.sub.j = 1, ... , c.sub.m) and switch their labeling   in Z′    compute Metric = Φ(Z′)    if Metric < BestMetric     save Z′ into Zsav     set BestMetric = Metric   end  end  set Z = Zsav end

[0116] For instance, the stopping criterion can be considered to be verified when no improvement of BestMetric is observed between two iterations. Of course other optimization algorithms can be considered and other stopping criteria can be used for optimizing the function Φ.

[0117] Without necessarily optimizing the function Φ, i.e. without necessarily finding the specific labeling function that minimizes the function Φ, the inventors have observed that the labeling functions satisfying the following criterion are good in the sense that, for such labeling functions, the log likelihood ratios LLR.sub.i,t tend to follow a Gaussian distribution for all the levels of the multilevel coded modulation:

Φ(Z)<1.7.sup.m×(0.0425×m−0.06)  (2)

expression in which m, as defined previously, corresponds to the number of levels of the multilevel modulation. Of course, the conventional Gray labeling function, inter alia, does not satisfy the expression (2) since the log likelihood ratios it produces do not follow a Gaussian distribution for all levels.

[0118] Hence: [0119] for m=4 (16QAM modulation): g(4)≈0.9187; [0120] for m=6 (64QAM modulation): g(6)≈4.7068; [0121] for m=8 (256QAM modulation): g(8)≈19.5321; [0122] etc.

[0123] Table 1 gives an example of labeling function Z that verifies the expression (2) above, in the case of m=4 (16QAM modulation). In table 1, j corresponds to the imaginary unit (j.sup.2=1).

TABLE-US-00002 TABLE 1 Example of 16 QAM labeling function (c.sub.1, . . . , c.sub.m) [00009]$Z\left(c_1,\mathrm{.Math.},c_m\right)\times \frac{2}{K}$ 0000 +3−1j 1000 +3−3j 0100 −3−3j 1100 −3−1j 0010 −1−1j 1010 −1−3j 0110 +1−3j 1110 +1−1j 0001 −3+1j 1001 −1+1j 0101 +1+1j 1101 +3+1j 0011 +3+3j 1011 +1+3j 0111 −1+3j 1111 −3+3j

[0124] FIG. 5 represents the labeling function Z given by table 1 in the complex plane. In FIG. 5, K is assumed to be equal to 2. It should be noted that a rotation of the constellation represented in FIG. 5 by a multiple of 90°, or an axial reflection of said constellation with respect to the real axis, or an axial reflection of said constellation with respect to the complex axis, or an inversion of all the bit values, or any combination thereof, does not modify the value of Φ(Z) and therefore yields a labeling function that is equally good.

[0125] Table 2 gives an example of labeling function Z that verifies the expression (2) above, in the case of m=6 (64QAM modulation).

TABLE-US-00003 TABLE 2 Example of 64 QAM labeling function (c.sub.1, . . . , c.sub.m) [00010]$Z\left(c_1,\mathrm{.Math.},c_m\right)\times \frac{2}{K}$ 000000 +5−7j 100000 +5−5j 010000 +3−5j 110000 +3−7j 001000 −7−7j 101000 −3−5j 011000 +1−3j 111000 −5−7j 000100 −3−7j 100100 −7−5j 010100 −1−5j 110100 −1−7j 001100 +1−7j 101100 +1−5j 011100 +7−5j 111100 +7−7j 000010 −3−3j 100010 −1−3j 010010 −7−3j 110010 −5−3j 001010 +5−3j 101010 +7−3j 011010 −5−5j 111010 +3−3j 000110 −7−1j 100110 +3−1j 010110 +5−1j 110110 −1−1j 001110 +1−1j 101110 −5−1j 011110 −3−1j 111110 +7−1j 000001 −3+3j 100001 +5+1j 010001 −1+1j 110001 −1+3j 001001 +1+3j 101001 +1+1j 011001 −5+1j 111001 +3+3j 000101 +3+1j 100101 −5+3j 010101 −7+3j 110101 −7+1j 001101 +7+1j 101101 +7+3j 011101 +5+3j 111101 −3+1j 000011 +7+5j 100011 −5+7j 010011 −7+7j 110011 +5+5j 001011 −5+5j 101011 +7+7j 011011 +5+7j 111011 −7+5j 000111 −3+7j 100111 +1+5j 010111 +3+5j 110111 −1+7j 001111 +1+7j 101111 −3+5j 011111 −1+5j 111111 +3+7j

[0126] FIG. 6 represents the labeling function Z given by table 2 in the complex plane. In FIG. 6, K is assumed to be equal to 2. It should be noted that a rotation of the constellation represented in FIG. 6 by a multiple of 90°, or an axial reflection of said constellation with respect to the real axis, or an axial reflection of said constellation with respect to the complex axis, or an inversion of all the bit values, or any combination thereof, does not modify the value of Φ(Z) and therefore yields a labeling function that is equally good.

[0127] Table 3 gives an example of labeling function Z that verifies the expression (2) above, in the case of m=8 (256QAM modulation).

TABLE-US-00004 TABLE 3 Example of 256 QAM labeling function (c.sub.1, . . . , c.sub.m) [00011]$Z\left(c_1,\mathrm{.Math.},c_m\right)\times \frac{2}{K}$ 00000000  −1−13j 10000000 −13−15j 01000000 +13−15j 11000000  −9−15j 00100000  −9−13j 10100000 +11−15j 01100000 −11−15j 11100000  +7−15j 00010000  +9−15j 10010000  +3−15j 01010000  −7−15j 11010000 +15−15j 00110000  +5−15j 10110000  −5−15j 01110000  +1−15j 11110000  −1−15j 00001000 −15−15j 10001000  +9−13j 01001000 −13−13j 11001000  +5−13j 00101000  −5−13j 10101000 −11−13j 01101000  −3−15j 11101000 +13−13j 00011000  +3−13j 10011000  −7−13j 01011000  +7−13j 11011000  −3−13j 00111000 +15−13j 10111000  +1−13j 01111000 +11−13j 11111000 −15−13j 00000100  +9−11j 10000100  +3−11j 01000100  +1−11j 11000100 +13−9j  00100100  +5−11j 10100100  −5−11j 01100100  −7−11j 11100100  −1−11j 00010100 −11−11j 10010100 +15−11j 01010100 +13−11j 11010100  −9−11j 00110100 −1−9j 10110100 −3−9j 01110100 −13−9j  11110100 −7−9j 00001100 +3−9j 10001100 −11−9j  01001100 +7−9j 11001100 −13−11j 00101100 +15−9j  10101100 +1−9j 01101100 +11−9j  11101100 −15−9j  00011100 −5−9j 10011100 +9−9j 01011100 −9−9j 11011100 +5−9j 00111100 −15−11j 10111100 +11−11j 01111100  −3−11j 11111100  +7−11j 00000010 −5−5j 10000010 +1−5j 01000010 +13−7j  11000010 −9−7j 00100010 −1−5j 10100010 +15−7j  01100010 +1−7j 11100010 +7−7j 00010010 +11−5j  10010010 −11−5j  01010010 +7−5j 11010010 +13−5j  00110010 +9−7j 10110010 −13−7j  01110010 −7−7j 11110010 +11−7j  00001010 +15−5j  10001010 −5−7j 01001010 −9−5j 11001010 +5−5j 00101010 −11−7j  10101010 −7−5j 01101010 −13−5j  11101010 −3−5j 00011010 −15−7j  10011010 +3−7j 01011010 −3−7j 11011010 −1−7j 00111010 +3−5j 10111010 +9−5j 01111010 +5−7j 11111010 −15−5j  00000110 −11−3j  10000110 +7−3j 01000110 −15−3j  11000110 +13−1j  00100110 +9−3j 10100110 −1−3j 01100110 +13−3j  11100110 +15−3j  00010110 +5−3j 10010110 +11−3j  01010110 +1−3j 11010110 −5−3j 00110110 −7−3j 10110110 −11−1j  01110110 −9−1j 11110110 −9−3j 00001110 −1−1j 10001110 −3−1j 01001110 +11−1j  11001110 −13−3j  00101110 +3−1j 10101110 −15−1j  01101110 +7−1j 11101110 +5−1j 00011110 +15−1j  10011110 −7−1j 01011110 −5−1j 11011110 +1−1j 00111110 −13−1j  10111110 +3−3j 01111110 −3−3j 11111110 +9−1j 00000001 +9+1j 10000001 −9+3j 01000001 +5+1j 11000001 −9+1j 00100001 −7+1j 10100001 +11+1j  01100001 −7+3j 11100001 −1+1j 00010001 +5+3j 10010001 +15+3j  01010001 −11+1j  11010001 +7+3j 00110001 +1+1j 10110001 −5+1j 01110001 +9+3j 11110001 +15+1  00001001 −15+1j  10001001 +3+1j 01001001 −15+3j  11001001 +3+3j 00101001 +13+3j  10101001 −1+3j 01101001 −3+1j 11101001 +11+3j  00011001 −11+3j  10011001 +7+1j 01011001 +1+3j 11011001 −13+1j  00111001 −3+3j 10111001 −13+3j  01111001 +13+1j  11111001 −5+3j 00000101 −5+5j 10000101 −15+5j  01000101 −9+5j 11000101 +7+7j 00100101 +3+5j 10100101 −7+5j 01100101 +15+5j  11100101 −11+5j  00010101 +11+5j  10010101 +1+5j 01010101 +7+5j 11010101 +13+5j  00110101 −13+5j  10110101 +9+5j 01110101 −1+5j 11110101 +5+5j 00001101 −15+7j  10001101 −5+7j 01001101 −3+7j 11001101 −3+5j 00101101 −7+7j 10101101 −13+7j  01101101 −11+7j  11101101 −1+7j 00011101 +1+7j 10011101 +11+7j  01011101 +13+7j  11011101 −9+7j 00111101 +9+7j 10111101 +3+7j 01111101 +5+7j 11111101 +15+7j  00000011 +5+9j 10000011 +7+9j 01000011 −15+9j  11000011 −13+9j  00100011 +13+9j  10100011 −1+9j 01100011 +9+9j 11100011 +11+9j  00010011 −11+9j  10010011 −9+9j 01010011 +1+9j 11010011 +3+9j 00110011 −3+9j 10110011 +15+9j  01110011 −7+9j 11110011 −5+9j 00001011 +11+11j 10001011  +1+11j 01001011  +7+11j 11001011 +13+11j 00101011  +3+11j 10101011  +9+11j 01101011 +15+11j 11101011  +5+11j 00011011  −5+11j 10011011 −15+11j 01011011  −9+11j 11011011  −3+11j 00111011 −13+11j 10111011  −7+11j 01111011  −1+11j 11111011 −11+11j 00000111 −15+13j 10000111  −5+13j 01000111  −3+13j 11000111  +7+13j 00100111  −7+13j 10100111 −13+13j 01100111 −11+13j 11100111  −1+13j 00010111  +1+13j 10010111 +11+13j 01010111 +13+13j 11010111  −9+13j 00110111  +9+13j 10110111  +3+13j 01110111  +5+13j 11110111 +15+13j 00001111  −5+15j 10001111 −15+15j 01001111  −9+15j 11001111  −3+15j 00101111 −13+15j 10101111  −7+15j 01101111  −1+15j 11101111 −11+15j 00011111 +11+15j 10011111  +1+15j 01011111  +7+15j 11011111 +13+15j 00111111  +3+15j 10111111  +9+15j 01111111 +15+15j 11111111  +5+15j

[0128] FIG. 7 represents the labeling function Z given by table 3 in the complex plane. In FIG. 7, K is assumed to be equal to 2. It should be noted that a rotation of the constellation represented in FIG. 7 by a multiple of 90°, or an axial reflection of said constellation with respect to the real axis, or an axial reflection of said constellation with respect to the complex axis, or an inversion of all the bit values, or any combination thereof, does not modify the value of Φ(Z) and therefore yields a labeling function that is equally good.

[0129] FIG. 8 represents the input log likelihood ratio distribution for the m polar code decoders when using the 16QAM labeling function given in table 1 above, which to be compared with the input log likelihood ratio distribution for the m polar code decoders when using the Gray labeling function (FIG. 4). As can be seen in FIG. 8, thanks to the particular labeling function identified by the inventors, the input log likelihood ratio distribution is substantially Gaussian for all the m polar code decoders. This makes this specific labeling function suitable for being used with any polar-code expecting input log likelihood ratios following a Gaussian distribution at the decoder, such as those polar-codes designed for binary input AWGN channels.

[0130] Another advantage of labeling functions verifying the expression (2) is that they enable performing simple code rate adaptation for all the polar codes used in the multilevel polar-coded modulation.

[0131] Indeed, since these labeling functions make each level of the system behave as if it were a binary input AWGN channel, and since the performance for such a binary input AWGN channel can be easily predicted, it is possible to adapt dynamically the respective code rates of the m polar codes to the propagation conditions experienced on each level. This can be done by reusing, on each level, the existing tools for predicting the performance on the binary input AWGN channel. For instance, it is possible to compute a function γ which gives the binary input mutual information as a function of the SNR of an AWGN channel. This formula for computing the function γ is considered to be known to the skilled person for the binary input AWGN channel. Then, it is possible to compute the inverse function of γ, denoted γ.sup.−1 in the sequel. It should be noted that γ.sup.−1 needs to be computed only once, and can then be memorized in storage means of the transmitting device 30. In order to adapt the code rate of the polar code of the i-th level, it is possible to compute the conditional mutual information I(b.sub.i; Y|b.sub.1, . . . , b.sub.i−1). Again, the formula for computing the conditional mutual information is considered to be known to the skilled person for the binary input AWGN channel. The SNR experienced on the i-th equivalent level channel, denoted {circumflex over (μ)}.sub.i, can then be estimated as {circumflex over (μ)}.sub.i=.sup.−1 (I(b.sub.i; Y|b.sub.1, . . . , b.sub.i−1)). Once {circumflex over (μ)}.sub.i is estimated, the code rate to be used for the polar code of the i-th level can be chosen to be ({circumflex over (μ)}.sub.i). Of course, it also possible to consider a backoff BO, for instance (−1) dB in logarithmic scale, in which case the code rate to be used for the polar code of the i-th level can be chosen to be (BO×{circumflex over (μ)}.sub.i).

[0132] FIG. 9 represents schematically the simulated expected throughput performance for both the 16QAM Gray labeling function and the 16QAM labeling function given in table 1 above, both in the context of a multilevel polar-coded modulation scheme as represented by FIG. 1, and both using polar-codes designed for binary input AWGN channels. Packets of 1024 16QAM symbols are sent. An automatic repeat request (ARQ) scheme is considered for evaluation, which involves that a packet is retransmitted when it is not correctly received. For a given SNR value, the rate splitting for each SNR is the one provided in FIG. 3 and defined according to the above-described procedure. The frozen bits of each of the 4 polar codes are chosen according to its rate. FIG. 9 represents also the mutual information for the 16QAM multilevel polar-coded modulation, i.e. the maximum achievable throughput as defined by information theory. As can be seen in FIG. 9, the optimized 16QAM labeling function introduces a gain of around 1 dB over the conventional 16QAM Gray labeling function.

[0133] FIG. 10 represents schematically the main steps of a method 50 for transmitting a data stream from a transmitting device 30 to a receiving device 40. Said data stream corresponds to one or several data packets comprising binary information to be transmitted from the transmitting device 30 to the receiving device 40.

[0134] In the present disclosure, the transmitting device 30 can be e.g. a user equipment (UE), a base station, a laptop, a tablet, a mobile phone, or any communicating object that can transmit a data stream to a receiving device 40. Similarly, the receiving device 40 can be e.g. a UE, a base station, a laptop, a tablet, a mobile phone, or any communicating object that can receive a data stream from a transmitting device 30.

[0135] The data stream to be transmitted is assumed to comprise m data sub-streams (e.g. data bit vectors b.sub.1, . . . , b.sub.m). Said data sub-streams can be independent from one another, or can correspond to a same data stream that has been split into m data sub-streams, e.g. by using a serial to parallel converter. As indicated above, in some embodiments, the data stream and/or the data sub-streams can comprise error correcting codes such as CRCs.

[0136] As illustrated by FIG. 10, said transmitting method 50 comprises the following steps, executed by the transmitting device 30:

[0137] a step S51 of encoding the data sub-streams with respective polar codes, such as to produce m polar-encoded data sub-streams (e.g. polar encoded data bit vectors c.sub.1, . . . , c.sub.m);

[0138] a step S52 of modulating the m polar-encoded data sub-streams onto symbols of a multilevel modulation comprising m levels defining 2.sup.m different symbols in the complex plane, according to a predetermined labeling function Z which bijectively associates a symbol to a set (c.sub.1, . . . , c.sub.m) of m bits from respective polar-encoded data sub-streams, such that successive different sets of m bits from respective polar-encoded data sub-streams are converted into successive symbols forming a symbol stream (e.g. vector X);

[0139] a step S53 of transmitting the symbol stream to the receiving device.

[0140] As discussed above, the 2.sup.m symbols of the multilevel modulation are distributed in the complex plane such that they are regularly spaced along the real axis and the complex axis by a spacing factor K, and the labeling function Z is such that the expression (2) above is verified. By using such a labeling function, it is possible to use off-the-shelf polar codes designed for binary input AWGN channels. Accordingly, each of the m polar codes used can be, in preferred optional embodiments, a polar code designed for a binary input AWGN channel.

[0141] As indicated above, especially when using polar codes designed for a binary input AWGN channel, code rate adaptation is simple. In such a case, the transmitting method 50 can comprise, in preferred optional embodiments:

[0142] determining SNRs experienced on each level by each data sub-stream during the propagation from the transmitting device 30 to the receiving device 40, on the equivalent level channel,

[0143] selecting the code rates of each of the m polar codes used based on the SNRs experienced by each data sub-stream.

[0144] For instance, the SNR determination and code rate selection can be performed as discussed above. However, other SNR determination algorithms and/or code rate selection algorithms can also be considered in other embodiments of the present disclosure.

[0145] FIG. 11 represents schematically an exemplary embodiment of a transmitting device 30 suitable for implementing the transmitting method 50.

[0146] In this exemplary embodiment, the transmitting device 30 comprises a processing circuit 31 comprising one or more processors and storage means (magnetic hard disk, solid-state disk, optical disk, electronic memory, etc.) in which a computer program product is stored, in the form of a set of program-code instructions to be executed in order to implement all or a part of the steps of the transmitting method 50. Alternatively, or in combination thereof, the processing circuit 31 can comprise one or more programmable logic circuits (FPGA, PLD, etc.), and/or one or more specialized integrated circuits (ASIC), and/or a set of discrete electronic components, etc., adapted for implementing all or part of said steps of the transmitting method 50.

[0147] The transmitting device 30 comprises also a communication unit 32, coupled to the processing circuit 31, allowing said transmitting device 30 to transmit the symbol stream. The communication unit 32 is preferably a wireless communication unit, in which case it corresponds to a radiofrequency circuit comprising components (antenna(s), amplifier(s), local oscillator(s), mixer(s), analog and/or digital filter(s), etc.) considered known to the skilled person.

[0148] In other words, the processing circuit 31 and the communication unit 32 of the transmitting device 30 form a set of means configured by software (specific computer program product) and/or by hardware (processor, FPGA, PLD, ASIC, discrete electronic components, radiofrequency circuit, etc.) to implement all or part of the steps of the transmitting method 50.

[0149] FIG. 12 represents schematically the main steps of a method 60 for receiving a data stream by a receiving device 40, the data stream comprising m data sub-streams transmitted by using the transmitting method 60.

[0150] As illustrated by FIG. 12, said receiving method 60 comprises the following steps, executed by the receiving device 40:

[0151] a step S61 of receiving the symbol stream (e.g. vector Y) including the data stream,

[0152] a step S62 of demodulating the symbols of the symbol stream in order to produce m polar-encoded data sub-streams, by applying an inverse function of the labeling function Z,

[0153] a step S63 of decoding successively the m polar-encoded data sub-streams such as to produce m data sub-streams, by applying respective polar-code decoders.

[0154] It should be noted that the demodulating S62 and decoding S63 steps can be performed simultaneously, by performing joint demodulation/decoding.

[0155] Of course, the polar code decoders used are adapted to the polar code encoders used by the transmitting device 30. Hence, in preferred embodiments, each of the m polar code decoders used is designed for decoding a polar code designed for a binary input AWGN channel. In preferred embodiments, each polar-code decoder used in the decoding step S63 is a list decoder. The list decoder is considered to be known to the skilled person.

[0156] FIG. 13 represents schematically an exemplary embodiment of a receiving device 40 suitable for implementing the receiving method 60.

[0157] In this exemplary embodiment, the receiving device 40 comprises a processing circuit 41 comprising one or more processors and storage means (magnetic hard disk, solid-state disk, optical disk, electronic memory, etc.) in which a computer program product is stored, in the form of a set of program-code instructions to be executed in order to implement all or a part of the steps of the receiving method 60. Alternatively, or in combination thereof, the processing circuit 41 can comprise one or more programmable logic circuits (FPGA, PLD, etc.), and/or one or more specialized integrated circuits (ASIC), and/or a set of discrete electronic components, etc., adapted for implementing all or part of said steps of the receiving method 60.

[0158] The receiving device 40 comprises also a communication unit 42, coupled to the processing circuit 41, allowing said receiving device 40 to receive the symbol stream. The communication unit 42 is preferably a wireless communication unit, in which case it corresponds to a radiofrequency circuit comprising components (antenna(s), amplifier(s), local oscillator(s), mixer(s), analog and/or digital filter(s), etc.) considered known to the skilled person.

[0159] In other words, the processing circuit 41 and the communication unit 42 of the receiving device 40 form a set of means configured by software (specific computer program product) and/or by hardware (processor, FPGA, PLD, ASIC, discrete electronic components, radiofrequency circuit, etc.) to implement all or part of the steps of the receiving method 60.

[0160] The above description clearly illustrates that by its various features and their respective advantages, the present disclosure reaches the goals set for it. In particular, by optimizing the labeling function such as to produce log likelihood ratios which follow a Gaussian distribution for all the levels of the multilevel polar-coded modulation, it is possible to have a simple and efficient multilevel polar-coded modulation scheme that can use off-the-shelf polar codes designed for a binary input AWGN channel. Also, such a multilevel polar-coded modulation scheme can provide predictable performance on each of the levels, which is suitable for allowing simple code rate adaptation, if necessary.