**APPROXIMATION OF LOG-LIKELIHOOD RATIOS FOR SOFT DECISION DECODING IN THE PRESENCE OF IMPULSE NOISE CHANNELS**

GELLE, Guillaume (3 rue Thiers, VILLE EN SELVE, 51500, FR)

GOUPIL, Alban (13 rue Goiot, REIMS, 51100, FR)

CLAVIER, Laurent (7 avenue Jacques Picavet, LAMBERSART, 59130, FR)

GUZZO BONIFACIO, Salvatore (102 rue des Stations, Bât. D2, LILLE, 59800, FR)

UNIVERSITE DE LILLE (42 rue Paul Duez, LILLE, 59800, FR)

CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE (3 rue Michel Ange, PARIS, 75016, FR)

*;*

**H03M13/45**

**H04L1/00**US20050058302A1 | 2005-03-17 |

YAZDANI R ET AL: "Linear LLR approximation for iterative decoding on wireless channels", IEEE TRANSACTIONS ON COMMUNICATIONS, IEEE SERVICE CENTER, PISCATAWAY, NJ. USA, vol. 57, no. 11, 1 November 2009 (2009-11-01), pages 3278 - 3287, XP011284838, ISSN: 0090-6778, DOI: 10.1109/TCOMM.2009.11.080038

DIMANCHE VINCENT ET AL: "On Detection Method for Soft Iterative Decoding in the Presence of Impulsive Interference", IEEE COMMUNICATIONS LETTERS, IEEE SERVICE CENTER, PISCATAWAY, NJ, US, vol. 18, no. 6, 1 June 2014 (2014-06-01), pages 945 - 948, XP011550325, ISSN: 1089-7798, [retrieved on 20140606], DOI: 10.1109/LCOMM.2014.2320921

HOU YI ET AL: "A non-linear LLR approximation for LDPC decoding over impulsive noise channels", PROC. 2014 IEEE/CIC INTERNATIONAL CONFERENCE ON COMMUNICATIONS IN CHINA (ICCC), IEEE, 13 October 2014 (2014-10-13), pages 86 - 90, XP032721806, DOI: 10.1109/ICCCHINA.2014.7008248

GAO CHENYU ET AL: "Gradient descent bit-flipping based on penalty factor for decoding LDPC codes over symmetric alpha-stable noise channels", PROC. 2017 IEEE/CIC INTERNATIONAL CONFERENCE ON COMMUNICATIONS IN CHINA (ICCC), IEEE, 22 October 2017 (2017-10-22), pages 1 - 4, XP033343368, DOI: 10.1109/ICCCHINA.2017.8330404

HASSEN BENMAAD ET AL.: "Clipping demapper for LDPC decoding in impulsive channel", IEEE COMMUNICATIONS LETTERS, INSTITUTE OF ELECTRICAL AND ELECTRONICS ENGINEERS, vol. 17, no. 5, 2013, pages 968 - 971, XP011511372, DOI: doi:10.1109/LCOMM.2013.031913.130060

POULLIAT C. ET AL.: "Performances des codes LDPCpour l'Ultra-Large Bande (ULB) impulsionnelle sous approximation gaussienne generalisee des interferences", THE PROCEEDINGS OF GRETSI CONFERENCE, 2007

CHUAN HSIAN: "Iterative SOVA decoding over alpha stable channels", JOURNAL OF ENGINEERING SCIENCE AND TECHNOLOGY, vol. 7, no. 3, June 2012 (2012-06-01)

CHUAN HISIAN: "SOVA decoding in Symmetric Alpha-Stable Noise", SYMPOSIUM ON INFORMATION & COMPUTER SCIENCES, 2011

M. JOHNSTON ET AL.: "Sum-Product Algorithm Utilizing Soft Distances on Additive impulsive Noise Channels", IEEE TRANSACTIONS ON COMMUNICATIONS, vol. 61, no. 6, June 2013 (2013-06-01), pages 2113 - 2116

WEI GUI: "Decoding Metric Study for Turbo Codes in Very impulsive Environment", IEEE COMMUNICATIONS LETTERS, vol. 16, no. 2, February 2012 (2012-02-01), pages 256 - 258, XP011409351, DOI: doi:10.1109/LCOMM.2011.112311.111504

YAZDANI ET AL., IEEE TRANSACTIONS ON COMMUNICATIONS, vol. 57, no. 11, November 2009 (2009-11-01)

CLAIMS 1. Method for processing an encoded signal (S) received from a noisy channel, in particular with impulsive noise in wireless transmission, prior to the decoding of the signal by a decoder, the method comprising: - applying to the encoded signal (S) a first function ( f said at least two parameters (a, b) being estimated beforehand by minimizing a second function (J) defined as a cost function of the first function (f 2. The method according to claim 1, wherein the encoded signal (S) comprising a plurality of symbols, the first function ( f 3. The method according to claim 2, wherein the first function (f wherein x is the amplitude of the symbol under analysis, and a and b are the beforehand- estimated parameters (a, b). 4. The method according to the preceding claim, wherein all input amplitudes (x) in absolute value that are below a threshold (sm) are processed using a linear part of the first function (f 5. The method according to the preceding claim, wherein the threshold ( s 6. The method according to any one of the preceding claims, wherein the encoded signal (S) is sampled with a number of n samples ( S , S 7. The method according to claim 6, wherein the second function (J) takes as input a signal (AS) obtained from the encoded signal (S) modified by a processing function using a random binary sequence increasing diversity of an estimated noise sequence. 8. The method according to the claim 1, wherein the cost function (J) is defined by wherein (f 9. The method according to any one of claims 2 to 8, comprising normalizing the amplitudes of the encoded signal (5) before estimation of the parameters (a, b) and application of the first function ( f 10. The method according to claim 8, wherein the input signal (AS) for the estimation of the parameters (a, b) is generated by the following expression: AS N 11. The method according to any one of the proceeding claims, wherein the estimation of the parameters (a, b) is performed by a non-linear optimization algorithm, preferably a Nelder-Mead algorithm. 12. The method according to the preceding claim, wherein an initial guess of the parameters (a, b) is defined before the first iteration of the minimization of the second 13. The method according to claim 6, comprising normalization of the samples of the modified signal (S’), in particular by computing a RMS value, the normalized modified signal (S”) being defined 14. System for processing an encoded signal (S) received from a noisy channel, in particular with impulsive noise in wireless transmission, prior to the decoding of the signal by a decoder, the system comprising: - a signal processing module adapted for applying to the encoded signal (S) a first function (f - a parameter estimation module adapted for estimating said at least two parameters (a, b) beforehand by minimization of a second function (J 15. System according to claim 14, comprising a module adapted for normalizing the amplitudes of the encoded signal (5) before estimation of the parameters (a, b) and application of the first function f 16. System according to claim 14 or 15, comprising a module adapted for normalizing the samples of the modified signal (S’), in particular by computing a RMS value, the normalized modified signal (S”) being defined by S" 17. A computer program product for processing an encoded signal (S) received from a noisy channel, in particular with impulsive noise in wireless transmission, prior to the decoding of the signal by a decoder, the computer program product comprising a support and stored on this support instructions that can be read by a processor, these instructions being configured for: - applying to the encoded signal (S) a first function ( f - estimating beforehand said at least two parameters (a, b) by minimizing a second function (J |

PRESENCE OF IMPULSE NOISE CHANNELS

The invention is related to methods and systems for processing an encoded signal received from a noisy channel, in particular before the decoding of the encoded signal.

Background of the invention

Sensor networks and ad hoc networks are wireless communication networks capable of organizing themselves without previously defined infrastructures. These networks are experiencing developments in many areas, from home applications to military applications. The implementation of these networks particularly in the case where they have a large number of communicating entities, such as several sensor nodes in the case of sensor networks, raises new issues such as the appearance of disturbances in transmission, notably interferences, presenting the peculiarity of being of impulsive nature in sufficiently dense networks.

To mitigate the disturbances of the system, it is necessary to code the information transmitted in a robust manner. It is known to code the information according to an error correction code which adds a redundancy.

Modem codes such as turbo codes or Low Density Parity Check (LDPC) codes can be used, with iterative decoding algorithms, trying to improve the quality of information being decoded iteration after iteration. The decoding of these codes uses as input the signal exiting from the channel, which may be in the form of log likelihood ratio (LLR).

The calculation of the LLR is based on channel modeling and depends largely on the distribution of the additive noise and interference. To do so, it is known to use a Gaussian model, which quickly reaches its limits and poorly matches the multiple communications in a network, i.e. the multiple access interferences generated by the nodes of a sufficiently dense network, which have an impulsive character and strongly influence the quality of the reception of information at the destination nodes. These interferences do not present Gaussian distributions but can be modeled by stable distributions which are characterized by a slower decay of the tails than in the Gaussian case, as described in the article of Hassen Benmaad et al. “ Clipping demapper for LDPC decoding in impulsive channel”, IEEE Communications Letters, Institute of Electrical and Electronics Engineers, 2013, 17 (5), pp.968-971.

The use of known iterative-based decoding methods is not efficient anymore to transmit information in an impulsive environment and a significant degradation of the performances of these methods appears and affects the reliability of the decoding of the information. The mismatch between the assumption made to design the receiver and the true interference distribution makes the decoding inefficient.

However, the main difficulty of changing the metric lies in the fact that the distribution of the interference is not necessarily a priori known. It depends on the network configuration, for instance the density of users and the design of the layers of Media Access Control (MAC). If stable distribution is quite a general model, there is not however an explicit form of the density, making the LLR extremely complex to compute.

The article of Poulliat C. et al, " Performances des codes LDPC pour VUltra-Large Bande (ULB) impulsionneUe sous approximation gaussienne genera/isee des interferences », in the proceedings of GRETSI conference, 2007, shows the use of a generalized Gaussian approximation of multi-users interferences. The type of considered impulsivity is quite specific and renders the approximation more restrictive.

The article of Chuan Hsian,“ Iterative SOVA decoding over alpha stable channels”, Journal of Engineering Science and Technology, June 2012, vol. 7(3) describes the use of a Cauchy-Bayesan metric. The article of Chuan Hisian, « SOVA decoding in Symmetric Alpha- Stable Noise », Symposium on Information & Computer Sciences, 2011, describes the use of Cauchy function in the path-metric. The Cauchy density is a special case of stable distribution but assumes a high impulsivity and is not adapted to low levels of impulsivity.

The article of M. Johnston et al, « Sum-Product Algorithm Utilizing Soft Distances on Additive impulsive Noise Channels », in IEEE Transactions on Communications, vol. 61, no. 6, pp. 2113-2116, June 2013., describes an adapted Sum-Product Algorithm (SPA) increasing the performances with respect to already-known SPA. It modifies the decoding algorithm in a deep manner, which implies a great complexity.

The known p-norm is an advantageously chosen solution because it only requires a rough estimate of the stable distribution, and not the estimation of the noise parameters

The article of Wei Gui,“Decoding Metric Study for Turbo Codes in Very impulsive Environment in IEEE Communications Letters, vol. 16, no. 2, pp. 256-258, February 2012, describes the utilization of the p-norm. This norm leads to good performances, but is quite complex to implement.

From a simplified point of view, a transmission chain is constituted of three main parts:

- A transmitter (emitter), wherein the information bits are generated by a source, then coded (for example using LDPC code), and after a modulation step, a signal X(t) is sent. For example, when the modulation is of BPSK (Binary phase-shift keying) type, the bits whom the value is zero are associated to X = -1 and the bits whom the value is one are associated to X = +l .

- A transmission channel, which propagates the emitted signal, and renders it noisy through the influence of a noise N.

The channel capacity C is the maximum average information that can be sent per channel use, between a source A sending a signal X and a receiver B receiving the signal S. This channel capacity can also be defined as the maximum mutual information I achieved for a given transition probability matrix (that is to say, for a fixed channel set of characteristics):

C = max / ( X , S )

- A receptor, wherein the signal received (also called“channel output”) is under the general form S = X + N, and is processed in order to obtain a modified signal S’ to be decoded, the modified signal S’ being a“LLR- like” signal.

From the receptor point of view, the decoding module, which is intended to decode the sent coded message, requires as input the signal exiting from the channel expressed under the form of a log likelihood ratio (LLR).

Thus, in order to obtain a proper working of the decoder, one has to transform the channel output S (the received value) into LLR (likelihood), which is thereafter called “demapping” or“demapper”.

Classically, this LLR is given as:

LLR (s) = log

where g _{N } is the probability density function of the noise N.

As the LLR expression involves the probability density function of the noise, which is usually not accessible in closed- form, it is highly desirable to find some sets of functions able to approximate at once and the most adequately the whole panel of LLR known and interesting curves for the targeted application(s).

Summary

There thus remains a need to further increase decoding performances in the case of transmission over a noisy wireless channel. More specifically, there is a need to reduce the impact that high intensity and low duration noises, also called impulsive noise, can have on the decoding process. Reduction of such an impact can improve the amount of data packets successfully decoded by a receiver, or, equivalently, reducing the packet error rate. The invention aims at satisfying this need by providing a method for processing an encoded signal S received from a noisy channel, in particular with impulsive noise in wireless transmission, prior to the decoding of the signal by a decoder, the method comprising:

- applying to the encoded signal S a first function f _{a b } based on at least two parameters a, b, depending on the channel noise, for obtaining a modified signal S’ to be decoded, the first function being chosen so as to mitigate the influence of the channel noise during the decoding of the modified signal S’,

wherein said at least two parameters a, b, are estimated beforehand by minimizing a second function J taking as input a signal AS obtained from the encoded signal S modified by a processing function using a random binary sequence.

The invention also provides a method for processing an encoded signal S received from a noisy channel, in particular with impulsive noise in wireless transmission, prior to the decoding of the signal by a decoder, the method comprising:

- applying to the encoded signal S a first function f _{a b } based on at least two parameters a, b depending on the channel noise, for obtaining a modified signal S’, S” to be decoded, the first function being chosen so as to mitigate the influence of the channel noise during the decoding of the modified signal S’, S”,

said at least two parameters a, b being estimated beforehand by minimizing a second function J _{a.b } defined as a cost function of the first function f _{a b } such that the (a, b ) = argmin ] (a, b).

The signal processing of the invention provides a generic method that can be implemented on many different types of receivers. It does not necessitate modification in the transmitter and can simply be added to existing protocols at the receiver side. It also has the advantage of not degrading communications in case of traditional communication conditions with a simple thermal Gaussian noise.

Furthermore, the method is self-adaptive and can be used in many different settings, with different noise distributions. Indeed, the set of f _{a b } functions, through the“on-demand” determination of the at least two parameters a and b, is able to provide a specific function f _{a,b } which approximate quickly and at lower cost the optimal LLR expected. Thus, through the use of self-adaptive analytical expressions it is possible to obtain high performances.

Accordingly, the demappers proposed according to the invention are able to self- adapt to a large panel of environments, and are easy to implement. The processing method according to the invention provides better results than state- of-the-art methods, allowing obtaining much more correctly decoded packets. The improvement is particularly significant if noise has an impulsive behavior.

Synchronization

Transmitters and receivers used for implementing the method according to the invention are usually driven by separate hardware clocks, the carrier frequency set up at the transmitter and the receiver sides may not be exactly the same. It is thus advantageous to perform a step of synchronization, for estimating and correcting carrier shift that may be introduced, to generate the encoded signal S.

The synchronization included in state-of-the-art systems is capable to apply the necessary correction in the case of noiseless signals or with a low level of noise. In the case of noisy signals, especially with impulsive noise, the default correction mechanism is not precise enough. An example of a BPSK modulated signal, with the presence of carrier offset, is illustrated in Fig. 7. The signal is composed of around 40000 samples, with a carrier offset of approximately l40Hz. The presence of carrier impairment produces the identified effect in the image: the amplitude, affected by noise, is not regular, and seems almost modulated, with a sinusoidal shape for the envelop. Carrier shift effect is also visible in constellation diagrams, with the spread of the samples in a circle.

In order to compensate the carrier offset, a synchronization metric M may be introduced, defined as: M = mean

represents the corrected sampled sequence from the encoded signal S, using a D/ offset, 5 indicates the imaginary part of the complex number, whereas abs and mean represent absolute value and average operations respectively. Corrected sampled sequence may be obtained by shifting every sample of a D/ offset, in frequency, and defined as: where T _{s } represents the sampling time. For the estimation and correction operation, an iterative process is advantageously defined, performed over a range of frequencies Nf ^{in } £ f £ N ^{ax } , in order to obtain frequency offset correction.

The aforementioned loop may be composed of the following steps: apply a temporary offset correction D^-, as defined above, estimate phase shift of the constellation, using for example linear approximation, compensate estimated phase shift, and compute synchronization metric M as defined above. Once all the synchronization metrics are collected, the frequency offset that allows minimizing the average dispersion over the imaginary coordinate of the complex space may be selected to be applied to the sampled signal under analysis, following the expression below:

Offset = argmin( f).

D/

By applying this kind of synchronization, it is possible to obtain corrective measures that are effective in a wider noise scenario, compared with the known solutions.

The step of symbol synchronization may be used to detect the synchronization sequence, which depends on the used protocol and is included in the encoded signal in order to identify the beginning of the useful data, also called“payload”.

Symbols and amplitudes

Preferably, the encoded signal S comprises a plurality of symbols a _{k } , the first function being applied to an amplitude x associated to each symbol. The encoded signal S can be expressed as: S(t) = å _{k } a _{k }h(t— kT _{symb } ) + n(t), where a _{k } are the sent symbols at rate T _{symb }, h(t ) is a filter comprising a shaping filter and a channel filter, and n(t) is the noisy part.

Thus, the encoded signal S is under the form S = X’ + N, where X’ is the sent message modified by h(t) and N the noisy part added during the transmission via the channel.

The encoded signal S may thus be sampled with a number of n samples

( S , S _{2 }, . , S _{n }) per symbol, n being an integer, for example n= 8. This samples can be viewed as {Si} = CLi. hi + rii (N per symbol).

A matched filtering may also be applied, for example after a synchronization step, in order to minimize Inter-Symbol Interference. Preferably all the processing steps to which the signal is subjected, may be followed by a decimation operation for extracting, from the oversampled signal, one sample per symbol.

A step of symbol synchronization may be used to detect the synchronization sequence, which depends on the used protocol and is included in the encoded signal in order to identify the beginning of the useful data, also called“payload”. In state-of-art systems, a PSK Detect step is responsible for the analysis of the signal, in order to associate one of the admissible point in the constellation to each symbol of the sequence.

The method may comprise a step of normalizing the amplitudes of the encoded signal S before estimation of the parameters, in particular at least the parameters a, b among others, and application of the first function f _{a b } . The amplitude normalization may be obtained through division of the amplitudes of the sampled symbols by an attenuation factor, preferably defined by A = ^ ^{SynCR } wherein

\Sync _{D } \

\Sync _{R } | represents a received synchronization sequence and \Sync _{D } \ represents a default synchronization sequence, for example defined at the receiver and by the receiver, for example known or chosen, and |n|being, by definition, an average absolute amplitude defined by |W| = sequence W of n samples W _{ΐ } W _{2 }, . , W _{h }.

Such amplitude normalization is useful to compensate path loss attenuation, due to the propagation of the signal over the channel, in particular in case of a wireless channel.

The samples can then be viewed, after normalization and decimation, as: { S _{d } } = a _{d }. (^) + ^ where (h _{d } / A) is a value very near to 1 (which is the goal of the normalization), or equivalently:

{¾} = a _{d } + n _{d } = x _{d } (1 per symbol, which is the goal of the decimation and n being equivalent to an impulsive noise due to this mathematical shaping), i being the index initial and d being the index after the decimation step.

Case of a BPSK modulation

In the case of a BPSK modulation, only two symbols are possible and a positive value will be associated to 1 when a negative value will be associated to 0. A step of symbol demapping then transforms the symbol sequence to a binary sequence, which corresponds to the output for the decoding. When a Forward Error Correction code is used, such a binary decision is not optimal, and a soft decision, that is to say based on real values, is preferable.

The quantized decimated values, optionally along with some information received from the PSK detect step to correct the phase shift, may be used for extract amplitudes, as can be seen in Fig. 9a and 9b.

When the method comprises a step of amplitude normalization for the encoded signal S, \Sync _{D } | is a default synchronization sequence, represented by BPSK symbol +1 and -1 (a BPSK modulation having been applied to the encoded signal).

First and second functions

The first function is preferably a non-linear function providing a parametrized approximation of log likelihood ratio. It is known how to find, write and choose a non-linear function providing a parametrized approximation of the log likelihood ratio defined previously as FFR (s).

The first function may be defined by:

wherein .¾ _{/ } is the amplitude of the received symbol under analysis, and a and b are the beforehand-estimated parameters, {¾} = a _{d } + n _{d } = x _{d } (1 per symbol ) as mentioned above.

Other non-limiting examples of the first function are:

• , illustrated in Fig.17,

• , illustrated in Fig.18,

• , illustrated in

Fig.19, and

• fa,b Od) = min _{a b } (a. x _{d } / (1 + b. x _{d } ^{2 })), illustrated in Fig.20.

An example of the first function f _{a b } (x _{d }) = min _{a b } (a. \x _{d }\ , . sign(x _{d } ) i is

represented in Fig. 12. Three examples, depending on the values of a and b are represented in Fig. 16. As visible, all input amplitudes .¾ _{/ } in absolute value that are below a threshold Sm are preferably processed using a linear part of the first function f _{a b }, and all amplitudes . _{¾/ } that exceed said threshold Sm are processed using an hyperbola type part of the first function f _{a b }.

The threshold may be defined by s _{x } =

To optimize the LLR approximating function and obtain good Bit Error Rate (BER) performance, different methods can be applied. For example:

• direct estimation of a and g if the noise distribution follows an alpha-stable law (as a reminder, an alpha-stable distribution is said symmetrical in zero when it has as characteristic function f _{a }(ί) = exp(— |yt|“)) ,

• kernel density estimation followed by the fitting of the estimated LLR with the chosen approximation,

• maximization of the mutual information between the source and the output of the approximated LLR.

According to the invention, the optimization of the LLR approximating function is done via the third method previously mentioned (maximization of the mutual information).

The second function ./ may be defined by a cost function / (a, b ) of the first function f _{a b }, the parameters a, b being estimated such that (a, b ) = argmin ] (a, b ). The second function J, defined as a cost function of the first function ( f _{a b } ), may take as input a signal obtained from the encoded signal modified by a processing function increasing diversity of an estimated noise sequence.

Preferably, the cost function J is defined by , wherein

f _{a b } is the first function, and AS) are the amplitudes of the input signal AS obtained from the encoded signal S modified by a processing function using a random binary sequence.

Estimation of a, b parameters

As explained above, the non-linear processing requires the definition of at least two parameters: a and b. Those parameters must be found in order to configure the non-linear block of the process. The goal is to estimate the values of a and b that can optimize the mutual information based on the channel condition estimated for the encoded signal S.

The receiver parameters a and b may be given by maximizing of the mutual information between the source and the output of the approximated LLR (before the decoding step), which means that one has to find a solution to the following optimization problem, related to the channel capacity C:

a.rg max _{a b }C = arg max _{a b }I (x , LLR _{a b } (S _{d })^j = arg max _{a b }l(X , / _{a,fc } (¾))

It is known how to solve this problem in the general case of memoryless binary input symmetric output channel (MBISO). Indeed, it is known (Yazdani et al. IEEE Transactions on communications, vol. 57, no. 11, november 2009) that then the capacity of the channel is defined by:

where:

• E _{L } is the expected value with respect of L

• L is the LLR, a random variable

• l is a realization of L

• f _{L } (0 is the probability density function of L

which can be also expressed, for approximated LLR l, as:

The preceding formula may be modified to be able to use the approximation function f _{a b } as: where U is a function depending on the f _{a b } function

The problem is then to maximize the expressed capacity between the source and the output of the approximated LLR, in order to determine the optimal receiver parameters a and b. In other terms, the goal is to find (estimate) a couple (a,b) such that:

(a, b) = arg

This can also be viewed as a need to minimize a second function J _{a b } .

The second function J _{a b } may be defined as a cost function for the optimization problem, allowing the maximization of channel capacity C. In other terms, the couple (a,b) is found through:

(a, b) = arg min _{a b } J _{U b }

The coefficient a and b may be estimated using an information theory criterion in order to maximize the mutual information between the encoded signal and the sequence after the non-linear mapping, based on a learning sequence and/or minimize a second function J _{a b } which is defined as a cost function for the optimization problem, allowing the maximization of channel capacity C.

Preferably, the estimation of the parameters a, b is performed by a non-linear optimization algorithm, preferably a Nelder-Mead algorithm. An initial guess of the parameters a, b may be defined before the first iteration of the minimization of the second function (./), for example by setting a = 4 and = 6.

Other algorithms may be used, for example

• a search on a predefined grid,

• Branch & bound method (if the problem is not convex).

• Newton method (Gradient descent).

Two methods may be used as input for the second function J _{a b } :

• Supervised way: X = {a , ... , a _{D }} is the transmitted sequence comprising at least in part a learning sequence and S _{d } = {x , ... , x _{D }} is the output of the channel resulting of the input X,

• Unsupervised way: only the received sequence S _{d } = {c , ... , x _{D }] is known by the receiver, and a step of rebuilding (from the sequence S _{d }) a“corresponding” input sequence AS is necessary.

In the case of supervised way, the optimization is triggered after receiving the channel output, by maximizing the mutual information according to the formula:

where:

• U is a function depending of the f _{a b } function

• X = {a , ... , a _{D } } is the transmitted sequence and a _{d } is one emitted symbol

• S _{d } = {x _{l } ... , x _{D }} is the received sequence and x _{d } is one received symbol (as a reminder {¾} = a _{d } + n _{d } = x _{d })

This method requires the use of a learning sequence of important length, which also means that non-informative transmitted bits must be added to the transmitted sequence. The learning sequence has to be long enough to guarantee the presence of impulsive events (else, the estimation of the parameters a and b will not be good enough to provide an efficient LLR estimation).

In the case of unsupervised way, the optimization is triggered after receiving the channel output, by maximizing the mutual information according to the formula:

where:

• U is a function depending of the f _{a b } function

• f _{a.b } is the first function,

• AS i are the amplitudes of an input signal AS (described thereafter) obtained from the encoded signal S,

• The second function J _{a b } is defined as a cost function of the first function f _{a b } .

Input signal AS

The input signal AS for the estimation of the parameters a, b used in the unsupervised way is called“ autodidact signal”, and is used to feed the second function J _{a b }. This autodidact signal corresponds to an approximated isolation of the estimated noise component introduced by the channel.

It has to be underlined that the unsupervised way of feeding the second function J _{a b } allows to the invention to be adapted to any transmission scheme, without modifying any standard nor adding non-informative transmitted bits (contrary to the supervised way). The input signal AS for the estimation of the parameters a, b may be generated by the following steps:

• determining an estimated noise sequence n _{d }

• optionally, increasing diversity in the estimated noise sequence n _{d } to obtain a sequence n' _{d }, by the use of a processing function,

• generate (at the receiver) a source sequence t _{s } .

• generate the received generated sequence AS _{d } , this step being equivalent to simulate a new channel for which the input sequence t _{s } is known (as it is generated artificially at the receiver), AS _{d } = t _{s } + N _{d } , where N _{d } corresponds to n _{d } or n' _{d }.

Preferentially, the estimated noise sequence (estimated noise samples) is obtained directly from the received samples x _{d } , using a sign detector, as : n _{d } = x _{d }— sign(x _{d }).

The processing function can be written as:

• n' _{d } = [ n _{d },— n _{d }] where the estimated noise sequence and the opposite of the same sequence are concatenated, or

• n' _{d } = r _{2d } [ n _{d },— n _{d }] where r _{2d } follows a Rademacher’s law of length 2d, which is the same length as the concatenated vector [ n _{d },— n _{d }] , or

• preferentially, the processing function consists in a independent and identically distributed sequence r _{d } where r _{d } follows a Rademacher’s law of length d (i.e. r _{d } is a random binary sequence). In that case, one has: n' _{d } = n _{d }. r _{d }

The source sequence t _{s } is a binary sequence produced at the receiver. As an example, the source sequence t _{s } can be an i.i.d. Rademacher’s sequence. More preferentially, as an alternative, the source sequence t _{s } contains only bits which values are 1.

Cumulatively, all the preferred embodiments of the invention lead to the following particularly preferred generated sequence: AS _{d } = 1 + n' _{d } = 1 + n _{d }. r _{d }.

Since the input sequence t _{s } is known and generated artificially at the receiver, the unsupervised way will not suffer from throughput loss. This is an important advantage of the invention compared to the other methods known from the state of the art.

For example, the input signal AS for the estimation of the parameters a, b may be generated by the following expression (k being an arbitrary index in these formulas): AS _{k } = 1 + n' _{k } wherein n' _{k } is a training sequence defined by n' _{k } = n _{k } r _{k }, with n _{k } = y _{k }— sign(y _{k } ), r _{k } being a random binary sequence distributed in {-1, 1 }, and y _{k } being the k ^{th } normalized amplitude of the encoded signal S.

Once the parameters estimation step is done, the normalized amplitudes x _{d } and the parameters a and b are known and used to feed the non-linear first function f _{a b } . The first function f _{a b } then provides a modified signal S’ which can be decoded.

The method, supervised or unsupervised, may comprise a step of normalization of the samples of the modified signal S’, in particular by computing a RMS value, the normalized modified signal S’ being defined by S" _{d } = Normalized f _{a b } (x _{d }) = =— ^{5 d }

^{6 6 J a l a }’ ^{b a } RMS(S' _{d }y with the general formula RMS(OL>) = wherein m is the number of samples available

after processing by the first function f _{a b } . In this case, a modified signal S”, which is to be decoded, is obtained after the processing of the signal S’.

The linear part of the first function f _{a b } , as can be seen in Fig. 12, may be determined by a slope that can introduce a non-negligible amplification effect on the output signal. This amplification may introduce amplitude values out of the admissible range of the decoding algorithm. This normalization operation on the modified signal S’ is useful to resize overall amplitudes and re-establish the limits within tolerable boundaries. The normalized samples keep the original shape of the non-normalized samples but with an adapted dynamic range.

System

The invention also provides a system for processing an encoded signal S received from a noisy channel, in particular with impulsive noise in wireless transmission, prior to the decoding of the signal by a decoder, the system comprising:

- a signal processing module adapted for applying to the encoded signal S a first function based on at least two parameters for obtaining a modified signal S’, the first function being chosen so as to mitigate the influence of the channel noise during the decoding of the modified signal, and

- a parameter estimation module adapted for estimating said at least two parameters beforehand by minimization of a second function taking as input a signal AS obtained from the encoded signal S modified by a processing function using a random binary sequence.

The invention also provides a system for processing an encoded signal S received from a noisy channel, in particular with impulsive noise in wireless transmission, prior to the decoding of the signal by a decoder, the system comprising: - a signal processing module adapted for applying to the encoded signal S a first function based on at least two parameters for obtaining a modified signal S’, the first function being chosen so as to mitigate the influence of the channel noise during the decoding of the modified signal S’, S”, and

- a parameter estimation module adapted for estimating said at least two parameters a, b beforehand by minimization of a second function , defined as a cost function of the first function f _{a b }, such that (a, b) = argmin J(a, b).

The system may further comprise a module adapted for normalizing the amplitudes of the encoded signal S before estimation of the parameters a, b and application of the first function f _{a b }.

The system may comprise a module adapted for normalizing the samples of the modified signal (S’), in particular by computing a RMS value, the normalized modified signal

(S’) being defined by S" _{d } = Normalized f _{a b } (x _{d }) with the g

^{6 }eneral formula RMS (w) = wherein m is the number of samples available after processing

by the first function f _{a b } .

The features defined above for the method apply to the system and vice-versa.

Computer program product

The invention also provides a computer program product for processing an encoded signal S received from a noisy channel, in particular with impulsive noise in wireless transmission, prior to the decoding of the signal by a decoder, the computer program product comprising a support and stored on this support instructions that can be read by a processor, these instructions being configured for:

- applying to the encoded signal S a first function based on at least two parameters for obtaining a modified signal S’ to be decoded, the first function being chosen so as to mitigate the influence of the channel noise during the decoding of the modified signal S’, and

- estimating beforehand said at least two parameters by minimizing a second function taking as input a signal AS obtained from the encoded signal S modified by a processing function using a random binary sequence.

The invention also provides a computer program product for processing an encoded signal S received from a noisy channel, in particular with impulsive noise in wireless transmission, prior to the decoding of the signal by a decoder, the computer program product comprising a support and, stored on this support, instructions that can be read by a processor, these instructions being configured for:

- applying to the encoded signal S a first function based on at least two parameters for obtaining a modified signal S’ to be decoded, the first function being chosen so as to mitigate the influence of the channel noise during the decoding of the modified signal S’, S”, and

- estimating beforehand said at least two parameters a,b by minimizing a second function J _{a },b, defined as a cost function of the first function f _{a b } , such that (a, b) = ar grain ] a, b ).

The features defined above for the method and the system apply to the computer program product.

Brief description of the drawings

The invention will be better understood on reading the following detailed description of non-limiting exemplary embodiments thereof and on examining the appended drawings in which:

- Figure 1 illustrates schematically an example of hardware system components for the transmission and reception of packets over a wireless network,

- Figures 2 and 3 are respective block diagrams of the implementation steps of a universal software radio peripheral (USRP) at the transmitter and receiver sides,

- Figure 4 illustrates an example of a receiver block diagram,

- Figure 5 shows an example of parallel processing for signal extraction,

- Figure 6 is a block diagram showing some steps for extracting amplitudes according to the invention,

- Figure 7, previously described, illustrates an example of a received baseband signal with carrier offset according to the state-of-the-art,

- Figure 8 shows the result, in time domain, of the synchronization algorithm according to the invention,

- Figures 9a and 9b, previously described, show a demodulation and amplitude extraction block diagram according to the invention,

- Figure 10 is a block diagram showing some steps of the method for processing an encoded signal according to the invention,

- Figure 11 is a block diagram showing some steps of the estimation of parameters depending on the channel noise according to the invention,

- Figure 12, previously described, illustrates an example of a non-linear function first function according to the invention, - Figure 13 illustrates the unsupervised way enabling the estimation of the parameters a and b,

- Figure 14 illustrates the LLR in case of an Alpha-stable (a = 1.5; y = 0.1) plus Gaussian noise (s = 0,2),

- Figure 15 illustrates the LLR in case of an e-contaminated (e = 0.1; k = 50) noise, and

- Figures 16 to 20 illustrate examples of the first function f _{a b } as previously described.

Detailed description

As illustrated in Fig. l, an example of a system used in the invention for the transmission and reception of packets over a wireless network comprises an universal software radio peripheral (USRP) and an associated driving computer (PC) at the transmitter side, and a USRP and a PC at the receiver side. The packets are transmitted over a noisy wireless channel. The invention is not limited to wireless transmissions, and can apply to non-wireless transmissions, for example in Power-Line communications or acoustic transmissions. The USRPs in use are for example NI-USRP 2952R. The invention is not limited to the use of a USRP, and a system-on-chip, a DSP or FPGA may also be used. Preferably, a single USRP equipped with two independent radio interfaces is used to send and receive the signal. In this case, synchronization between transmitter and receiver is automatically achieved due to the share of the same clock at transmitter and receiver sides.

The operations performed by the hardware inside the USRP can be abstracted both for transmission and for reception of signals as illustrated in Fig. 2 and Fig. 3 respectively. These operations are necessary for the conversion of the signal between baseband and passband and vice versa. In the present invention, bands centered in 915 MHz and 2.45 GHz may be used, corresponding to bands of Industrial Scientific and Medical bands (ISM). The invention is however not limited to specific carrier bands.

In order to transmit a signal over a communication medium, it is necessary to transpose it from the baseband, which is centered around 0Hz, to the desired passband. To do so, the USRP receives the waveform to transmit from the computer, filters and transposes, in a step of up-conversion, the spectrum of the signal. The signal is then encoded and modulated. The signal in passband is amplified prior to be directed toward the antenna to be transmitted over the wireless medium.

The operations to be performed at the reception are meant to revert the process performed by the transmitter. In this example, in the receiving chain of the USRP, the signal captured by the antenna is amplified using a low noise amplifier. In a step of down-conversion, the received signal is transposed back to baseband applying the same carrier frequency used by the transmitter. A filtering of the down converted signal is advantageously performed in order to retain only the frequencies of interest. The filtered signal is then sampled and transferred to the computer to be further processed.

The general structure of the receiver application which can be used for the implementation of the invention is represented by the block diagram of Fig. 4.

A step 2 for establishing the setup of the receiver is performed once at startup, similarly to the one at the transmitter, by initializing the USRP with parameters such as carrier frequency, sampling frequency, or amplification factor. During such step, it is also possible to qualify some of the properties that are specific to the coding technique in use.

The next step 3 is to acquire the baseband sampled signal produced by the USRP. In this example, the data type is an array of complex numbers containing the In-Phase and Quadrature (I/Q) components of each sample, obtained by the down-conversion, as the real and imaginary part of a complex number. The structure of the received signal can be expressed as:

where T _{s } represents the sampling time.

Due to its structure, the USRP is continuously down-converting and sampling what is detected. Consequently, the received sampled signal may or may not contain an actual transmission. Hence, a step 4 of extracting useful signal may be performed, in order to estimate the signal and noise level, for extracting from all received samples only those that are likely to carry out a message. This is meant to reduce the execution time avoiding to process samples containing merely noise. In order to improve efficiency in signal processing and, consequently, being able to achieve real time decoding, a parallelized signal analysis may be implemented, as depicted in Fig. 5. Each“ Process signal” component of parallel processing is in charge to detect boundaries for the presence of a packet. All results are recombined to produce boundaries, and used to create set of samples to be further analyzed. The remaining have been estimated to contain just channel noise.

The signal may then be subjected to another sampling process 5, as shown in Fig. 4. This is useful to have the correct number of samples per symbol necessary for the further demodulation step. Even though the sampling frequency is often a configurable parameter at the receiverside, there can be a difference between the requested sampling frequency and the actual sampling frequency, due to hardware limitations. Thus, resampling may be performed to eliminate possible impairments.

A further step 10 of demodulation and amplitude extraction shown in Fig. 4, which can be decomposed as illustrated in Fig. 6, is performed in the illustrated example. The resampled noisy signal is preferably low-pass filtered 11 prior to continue further in the demodulation chain, to form the encoded signal S.

As previously described, it is advantageous to perform a step of synchronization 12, for estimating and correcting carrier shift that may be introduced.

An example of a BPSK modulated signal, with the presence of carrier offset, is illustrated in Fig. 7.

Fig. 8 shows the result of this synchronization step, performed on the encoded signal S, in time domain.

A step of demodulation with amplitude extraction 20, for example BPSK demodulation with amplitude extraction, as represented in Fig.6 and previously described, and detailed in Fig. 9a and 9b, is performed in order to extract amplitude associated to every symbol of the encoded signal S. This amplitude sequence x _{d } and parameters a and b will be the inputs for non-linear operation performed in the invention.

The method for processing the encoded signal according to the invention, marked by the black arrow of Fig. 4, will be further detailed below.

Once the encoded signal has been processed, the transmitted message is retrieved in a step of signal decoding.

The validity of the decoded sequence may be verified in a step of“Message display & statistics update”. In this example, the sequence is subjected to a verification of the cyclic redundancy check (CRC) code added as tail of the payload at the transmitter side. If such verification is successfully passed, the payload, the essential data that is being carried within a packet, is extracted, converted in the string format and displayed for the user. The reception statistics may be updated to keep count of the current received packet as successfully or unsuccessfully decoded.

An example of detailed steps for processing the encoded signal S according to the invention is shown in Fig. 10.

As previously described, a first step 30, detailed in Fig.l l, is the beforehand estimation of parameters a, b depending on the channel noise, which are necessary to properly configure the computation of the non-linear first function /,, _{/ },, followed by an operation of RMS normalization, before the obtaining of the decoded sequence by a decoding step. This step uses the second function / that is a cost function / (a, b) of the first function f _{a b }, the parameters a, b being estimated such that (a, b ) = argmin ] (a, b ) , by using, preferably and as shown in Fig. l 1, a particular type of simplex method known as Nelder-Mead algorithm. It makes use of three points, in a bi-dimensional space, that are relocated during the iterations of the algorithm in order to obtain the estimated optimized pair a, b.

Preferably, as previously described and as shown on Fig. 11, for achieving the estimation of the parameters a, b, a preliminary operation of amplitude normalization is realized, leading to { S _{d } } = x _{d } as this step is realized after the decimation step.

Once obtained the normalized amplitudes x _{d } , as previously described, the autodidact Signal AS is generated, for example by the following expression (k being an arbitrary index in these formulas): AS _{k } = 1 + n' _{k }, wherein n' _{k } is a training sequence defined by n' _{k } = n _{k } ^{r } _{/ } with ^{n }k = k— signfyk , ^{r }k being a random binary sequence distributed in (-1, 1 }, and y _{k } being the k ^{th } normalized amplitude of the encoded signal S. This signal AS is used as input for the cost function / (a, b ) , by defining an initial guess for parameters a and b, for example a = 4 and b = 6.

Preferably and as previously described, to every amplitude x of every symbols under analysis, is applied a non-linear function f _{a,b } i ^{x } _{d } ) ^{= m m }a,b ( ^{a } l _{d } sign{x _{d }), as

illustrated in Fig. 12, for obtaining a modified signal S’ to be decoded. In this example, all input amplitudes .¾ _{/ } in absolute value that are below a threshold sm are processed using a linear part of the first function f _{a b }, and all amplitudes .¾ _{/ } that exceed said threshold sm are processed using an hyperbola type part of the first function f _{a b }. The threshold may be defined by: s _{TH } =

In the illustrated example, a normalization of the samples of the modified signal S’ is performed, by computing a RMS value, the normalized modified signal S” being defined by the general formula RMS(x) =

^{2 } åml ^{x }ml wherein m is the number of samples available after processing by the first function

( fa,b )·

The processing of the encoded signal S according to the invention is then concluded, the produced output S” may be adopted to perform decoding, as above-written.

Figure 13 illustrates the unsupervised way enabling the estimation of the parameters a and b, where only the received sequence S _{d } = {x _{l } ... , x _{D }] is known by the receiver. In this case, as can be seen on Fig. 13, a step of building a new channel from: • an artificially generated input source sequence t _{s } , and

• an estimated noise sequence (diversified ( n' _{d } ) or not (n _{d })), obtained directly from the received samples x _{d } = {¾},

with a“corresponding” input sequence AS _{d } = t _{s } + N _{d } is necessary. Once the parameters a and b obtained, by minimizing the second function ] _{a b } , those parameters are used as input of the first function f _{a b } in order to obtain the signal S’ or S” to be decoded.

Figure 14 illustrates the exact LLR obtained in case of an Alpha-stable (a = 1,5 ; g = 0,1) plus Gaussian noise (s = 0,2).

Figure 15 illustrates the exact LLR obtained in case of an e-contaminated (e = 0,1 ; K = 50) noise.

Figures 16 to 20 illustrate examples of the first function f _{a b } as previously described, chosen:

• so as to mitigate the influence of the channel noise during the decoding of the modified signal S’, S”,

· so as to provide a parametrized approximation of the log likelihood ratio.

• preferably under the form of a non-linear function.

Decoding performance tests are now going to be presented. In the following test examples, a single USRP equipped with two independent radio interfaces is used, and received signals have been added with a-stable random noise aiming to simulate an impulsive channel noise. Such a noise depends on several parameters a is the characteristic exponent of the probability distribution function. For an a -stable distributed noise, a lower a implies a higher impulsivity of the noise b indicates the skewness, g is used as a scale parameter, allowing to adjust the produced noise in amplitude, and 5 is defined as a location parameter introducing a shifting offset. All test presented hereafter have been produced using b = 0 and 5 = 0.

Tests have been performed using three different coding techniques, more precisely

Low Density Parity Check (LDPC) coding, Convolutional coding and Turbo codes.

The results, presented in table form indicating main configuration parameters, are grouped based on the type of coding technique under evaluation. The tables represent the percentage of correctly decoded packets; that is to say which have a CRC check with positive result. Results for pure decoding, using amplitudes without applying any further processing, for classic LLR decoding, amplitudes processed using a standard equation for the computation of the LLR, LLR (s) = log X the expected bit and s the channel output, assuming an additive white gaussian noise with standard deviation s = 0.8, and for the method according to the invention are presented and compared.

Table 1 : LDPC coding

Table 2: Convolutional coding - short payload

Table 3: Convolutional coding - medium payload

Table 4: Convolutional coding - long payload

Table 5: Turbo coding

From these results, it is obvious that the processing method according to the invention provides better results than state-of-the-art methods. From 13% up to 80% of additional correctly decoded packets can be obtained.

Performances of the automatic estimation of the channel configuration parameters according to the invention are also shown, in order to determine the improvement in terms of decoding capabilities and to validate the implemented channel estimation technique.

All tests have been performed using convolutional coding technique, with b = 0 and 5 = 0, and the results are presented in table form, indicating main configuration parameters.

In this example, transmitter and receiver are played by two completely separated couples of USRP and controlling unit. In such configuration, the two systems do not share the same clock, hence the synchronization operations are important to allow a proper analysis of the processed sampled signal S’, S”. The tables represent the percentage of correctly decoded packets, for pure decoding, for classic LLR decoding, for a processing according to the 2 invention but with theoretical a and b parameters mathematically calculated by a =— and b = 2 Va + 1, and for the parameters estimation method according to the invention.

Table 6: Convolutional coding - a=1.7

Table 7: Convolutional coding - a=1.8

Table 8: Convolutional coding - a=1.9

Table 9: Convolutional coding - Natural thermal noise

From these results, it is obvious that the estimated parameters according to the invention do not degrade the decoding performances, and even improve these performances most of the time.

**Previous Patent:**METHOD AND SYSTEM FOR MEASURING AN OBJECT BY MEANS OF STEREOSCOPY

**Next Patent: HEAT-INSULATED PRECAST CONCRETE PART**