# Process for estimating the channel in a OFDM communication system, and receiver for doing the same

Imported: 17 Feb '17 | Published: 23 Sep '14

USPTO - Utility Patents

## Abstract

A process for estimating a channel in an OFDM communication system, wherein an OFDM sample is based on a time-frequency grid comprising a set of pilots may include the steps of receiving, by a receiver, OFDM sample based subframes comprising a set of pilot signals, performing an iterative loop for the purpose of computing (Mk, kk) respectively corresponding to the covariance matrix and the mean of the probability distribution of p(hk|yk, lk) at time instant k, yk, and lk respectively corresponding to the received signal and to a priori known information, the value of (Mk, kk) being computed as a function of (Mk−1, kk−1) as well as the values of the pilots received at instant k.

## Description

### TECHNICAL FIELD

The invention relates to telecommunications and more particularly to a method for estimating the transmission channel vector in an Orthogonal Frequency Division Multiplexing (OFDM) communication system, and receiver for doing the same.

### BACKGROUND ART

Orthogonal Frequency Division Multiplexing (OFDM) communication systems is a widely developed technique in the field of wireless digital communications thanks to the high possibilities offered by digital signal processing and Discrete Fourier Transform computations.

However, the efficiency of the process strongly depends on the estimation of the channel, which is achieved thanks to the communication of pilot signals within the subframe.

FIG. 1 illustrates the particular structure of the time-frequency OFDM grid, corresponding to a sub-frame, with the location of the pilots.

Many mechanisms are known in the art for computing an estimation of the channel vector hk (k corresponding to the discrete instant k), including iterative algorithms deriving the particular knowledge of hk from the knowledge of hk−1 and the new information brought be the pilots received at the instant k.

Generally speaking, the channel estimation is much more effective when one knows the law of variation of the channel—its time correlation—and when one may store the high volume of data derived from the received pilots so as to elaborate the most comprehensive knowledge of the channel estimation.

However, it is clear that the storage of a high volume of data increases the amount of storage required for the estimation process and moreover the amount of digital processing resources for processing such information.

For those reasons, many channel estimation processes are based on an iterative mechanism for computing the estimate ĥk of the vector channel only based on the preceding estimate ĥk−1 and the newly received information. Such iterative mechanisms show clear limitations particularly when the channel shows a large delay spread prohibiting accurate interpolation between the OFDM subcarriers.

In such an environment, but not exclusively, there is a need for an effective iterative mechanism, which takes advantage of all the information being available at a given instant k−1 for computing the estimated channel ĥk at the instant k.

Such is the technical problem to be solved by the present invention

### SUMMARY OF THE INVENTION

It is an object of the present invention to provide a process estimating a channel for an OFDM communication system which can be used in particular in a mobile handset.

It is a further object of the present invention to provide a process of channel estimation for an OFDM communication system which reduces the amount of data to be stored for computing an accurate estimate of the channel.

It is another object of the present invention to provide a channel estimation method, using the Maximum Entropy principle, which allows one to assign probability distributions for the channel, when the estimator only knows the last past inferred channel distribution and the new received pilot symbols.

Those and other objects are achieved by means of a process for estimating the channel in an OFDM communication system, with OFDM samples formed of a time-frequency grid comprising a set of pilots.

The process involves the steps of:

• receiving OFDM subframes comprising a set of pilot signals;
• performing an iterative loop for the purpose of computing (Mk, kk) respectively corresponding to the covariance matrix and the mean of the probability distribution of p(hk|yk, Ik) at time instant k, yk, and Ik respectively corresponding to the received signal vector and to a priori known information;
said value of (Mk, kk) being computed as a function of (Mk−1, kk−1), λ being the time correlation (left to be defined) between the true channels hk and hk−1 as well as the values of the pilots received at instant k.

More particularly, the process performs an estimation of the channel based on the mean kk and variance Mk of the distribution of the channel estimation in accordance with the following formulation:

${ M k = λ 2 ⁡ ( M k - 1 + 1 - λ 2 λ 2 ⁢ Q ) × ( λ 2 σ 2 ⁢ P k ⁡ ( M k - 1 + 1 - λ 2 λ 2 ⁢ Q ) + I N ) - 1 k k = λ ⁢ ⁢ k k - 1 + 1 σ 2 ⁢ M k ⁢ P k ⁡ ( h k ′ - λ ⁢ ⁢ k k - 1 )$
With

• N being the number of subcarriers in the OFDM symbol;
• kk being the estimated channel at instant k;
• Mk being an N×N matrix representative of the variance of the estimated channel at instant k;
• Q being a matrix initialized to a predetermined value;
• Pk is an N×N diagonal matrix comprising a “1” at the location corresponding to a subcarrier at pilot position and a “0” elsewhere;
• IN being the N×N identity matrix;
• σ2 being the variance of the noise process;
• h′k being equal to the received signal yk divided pointwise by the corresponding signal vector sk;
• λ corresponding to the channel time correlation between instant k−1 and k.

In one embodiment, the matrix Q is an a priori frequency covariance matrix for an L-tap uniform channel, i.e. Q=FLFLH, with FL the matrix formed by the L first columns of the N×N DFT matrix

Preferably, the time correlation factor (λ) corresponds to the correlation factor in the auto regressive model of the first order for the time series hk.

In one embodiment, the parameter Mk is forwarded to the base station as an information representative of the quality of the channel estimation.

Preferably, the process is used for the Long Term Evolution (LTE) standard.

A process is also provided, which can be used even without the knowledge of the time correlation variable λ.

In that case the process involves the following steps:

• performing an iterative loop for the purpose of computing (Mk(λ), kk(λ)) respectively corresponding to the variance and the mean of the probability distribution of p(hk|yk, Ik) at time instant k, yk, and Ik respectively corresponding to the received signal and to a priori known information;
said value of (Mk(λ), kk(λ)) being computed as a function of (Mk−1(λ), kk−1(λ)) for a predetermined set of values of time correlation λ n} with n=1 to M, for some positive integer M, between the true channel hk and hk−1, as well as the values of the pilots received at instant k,

$h ^ k = ∑ n = 1 N ⁢ ⅇ - x ⁡ ( λ n ) ⁢ det ⁡ [ X ⁡ ( λ n ) ] ⁢ k k ( λ n ) 1 N ⁢ ∑ n = 1 N ⁢ ⁢ ⅇ - x ⁡ ( λ n ) ⁢ det ⁡ [ X ⁢ ( λ n ) ] ⁢ ⁢ With ⁢  ⁢ { M k ( λ ) = λ 2 ⁡ ( M k - 1 ( λ ) + 1 - λ 2 λ 2 ⁢ Q ) × ( λ 2 σ 2 ⁢ P k ⁡ ( M k - 1 ( λ ) + 1 - λ 2 λ 2 ⁢ Q ) + I N ) - 1 k k ( λ ) = λ ⁢ ⁢ k k - 1 ( λ ) + 1 σ 2 ⁢ M k ( λ ) ⁢ P k ⁡ ( h k ′ - λ ⁢ ⁢ k k - 1 ( λ ) ) And ⁢  ⁢ { X ⁡ ( λ ) = ( I + P k σ 2 ⁢ ( λ 2 ⁢ M k - 1 ( λ ) + ( 1 - λ 2 ) ⁢ Q ) ) - 1 X ⁡ ( λ ) = ( λ ⁢ ⁢ k k - 1 ( λ ) - h k ′ ) H ⁢ X ⁡ ( λ ) ⁢ P k σ 2 ⁢ ( λ ⁢ ⁢ k k - 1 ( λ ) - h k ′ )$
With kk(λ) being the estimated channel at instant k depending on λ,

• Mk(λ) being a matrix representative of the variance of the estimated channel at instant k; depending on λ,
• Q being a matrix initialized to a predetermined value;
• Pk is a diagonal N×N matrix comprising a “1” at the location of a pilot symbol;
• IN being the N×N identity matrix;
• σ2 being the variance of the noise;
• h′k being equal to the received signal yk pointwise divided by the corresponding pilot sk.

At last, the invention encompasses a receiver for an OFDM communication system comprising:

• a means for receiving a OFDM subframes comprising a set of pilots signals;
• a means for performing an iterative loop for the purpose of computing (Mk , kk) respectively corresponding to the variance and the mean of the probability distribution of p(hk|yk, Ik) at time instant k, yk, and Ik respectively corresponding to the received signal and to a priori known information; said value of (Mk, kk) being computed as a function of (Mk−1, kk−1), λ being the known time correlation between the true channel estimate hk and hk−1, as well as the values of the pilots received at instant k.

### DESCRIPTION OF THE PREFERRED EMBODIMENT

Some preliminary theoretical considerations (I) show to be appropriate before describing one embodiment of a channel estimation process (II)

I. Preliminary Theoretical Considerations

Channel estimation, along with most synchronization procedures, is an intricate information theoretic problem. Indeed, while Shannon (C. E. Shannon, “A mathematical theory of communications”, Bell System Technical Journal, vol. 27, no. 7, pp. 379-423, 1948) allows us to derive the capacity of a channel for which all synchronization parameters, plus the noise variance, are perfectly known, no such theory exists when the knowledge of some of these parameters is missing.

More precisely, for a scalar communication y=hx+n, x, for some codebook , if h is unknown, then the maximum a posteriori estimate for x is

$x ^ = arg ⁢ ⁢ max x ∈ X ⁢ ∫ h ⁢ ∫ n ⁢ p ⁡ ( x ❘ h , n ) ⁢ p ⁡ ( h ) ⁢ p ⁡ ( n ) ⁢ ⅆ h ⁢ ⅆ n ( 1 )$

which requires to have an a priori p(h) for h. But this a priori is too impractical to obtain and would require to know all possible channel realizations and their respective probability. As a consequence to this strong difficulty, most contributions in the synchronization field have provided various empirical models based on field observations in order either to give an expression to p(h) or, more practically, to propose good channel estimators ĥ to h.

The difficulty of handling estimation problems when little side information is available is treated by Jaynes, through the Bayesian probability field, thanks to the maximum entropy principle (MaxEnt) (E. T. Jaynes, “Information Theory and Statistical Mechanics”, Physical Review, APS, vol. 106, no. 4, pp. 620-630, 1957).

However MaxEnt does not allow to perform updates of probability when new information, such as new pilots in the channel estimation problem, is available. In this case, the complete set of past symbols along with a prior distribution for the channel h0 at time t=0.

The inventors have discussed such question fro the OFDM framework in the article “A maximum entropy approach to OFDM channel estimation”, arxiv Preprint http://arxiv.org/pdf/0811.0778 when the channel delay spread, the channel time correlation and the signal-to noise ratio (SNR) are alternatively known or unknown. When these parameters are not perfectly known, MaxEnt provides channel estimates minimizing the estimate mean square error (MMSE estimates) that outperform classical estimates which use empirical (often erroneous) models. Recent contributions in the Bayesian probability field enable one to perform probability updates, in particular based minimum cross entropy considerations [See A. Caticha, “Lectures on Probability, Entropy and Statistical Physics”, arXiv:0808.0012v1 [physics.data-an], 2008].

1. OFDM System Model

Preliminary notations: in the following boldface lower case represent vectors, capital boldface characters denote matrices (IN is the N×N identity matrix). The transposition operation is denoted (·)T. The Hermitian transpose is denoted (·)H. The operator diag(x) turns the vector x into a diagonal matrix. The symbol det(X) is the determinant of matrix X. The symbol E[·] denotes expectation. The Kronecker delta function is denoted δx that equals 1 if x=0 and equals 0 otherwise.

Consider a single cell OFDM system with N subcarriers. The cyclic prefix (CP) length is NCP samples. In the timefrequency OFDM symbol grid, pilots are found in the symbol positions indexed by the function φt(n){0, 1} which equals 1 if a pilot symbol is present at subcarrier n, at symbol time index t, and 0 otherwise. The time-frequency grid is depicted in FIG. 1. Both data and pilots at time t are modeled by the frequency-domain vector St CN with pilot entries of amplitude |st,k|2=1. The transmission channel is denoted ht CN in the frequency-domain with entries of variance E[|ht,k|2]=1. The additive noise is denoted nt CN with entries of variance E[|nt,k|2]=σ2. The noise process is modeled as Gaussian with independent and identically distributed (i.i.d.) entries, n˜(0, σ2IN)

The time-domain representation of ht is denoted vt ε {dot over (L)}{dot over ( )}with L the channel length, i.e. the channel delay spread expressed in OFDM-sample unit. The frequency-domain received signal yt CN is then:
yt=diag(ht)st+nt   (2)

We will also denote,
∀k ε {1, . . . , N}
h′k=yk/sk=hk+nk/sk

and h′=(h′1, . . . , h′N)T (here, the time index t is implicit).

The channel ht evolves in time with coherence time function λ(τ) such that, independently of the channel delay spread index

$E ⁡ [ v i , t ⁢ v i , t + τ * ] = λ ⁡ ( τ ) L ( 3 )$

This model is justified by maximum entropy considerations (see e.g. the reference R. Couillet, M. Debbah, “A maximum entropy approach to OFDM channel estimation”, arxiv Preprint http://arxiv.orq/pdf/0811.0778) and leads for p(ht|ht+τ) to be assigned a Gaussian distribution of mean λ(τ)ht+τ and variance

$E ⁡ [ h t ⁢ h t + τ H ] = ( 1 - λ ⁡ ( τ ) 2 ) ⁢ Q ⁢ ⁢ With ( 4 ) Q nm = E ⁡ [ ∑ k = 0 L - 1 ⁢ ∑ l = 0 L - 1 ⁢ v k ⁢ v l * ⁢ ⅇ - 2 ⁢ πⅈ ⁢ kn - lm N ] = 1 L ⁢ ∑ k = 0 L - 1 ⁢ ⅇ - 2 ⁢ πⅈ ⁢ ⁢ k ⁢ n - m N ( 5 )$

This work aims at providing MMSE estimates ĥtCN of the vector ht when the receiver, at time t is only aware of:

(i) the last stream of pilots received at some time tp≦t,

(ii) the last density function p(ht0) assigned to ht0. All other system parameters, e.g. channel delay spread, channel coherence time, SNR, are either known or unknown. This will be clearly stated in the following and will be denoted I.

2. ME Principle and Application to Channel Estimation.

The fundamental difference between orthodox and Bayesian probabilities lies in the existence of prior distributions in the Bayesian realm. In the present situation, when estimating ht from the observation yt, orthodox probabilists decide for a maximum likelihood estimator ĥt=arg maxht p(ht|h′t), following a given model linking h′t to ht. Bayesian probabilists would note that
p(ht|h′t, I)=p(h′t|ht, I)p(ht|I)/p(h′t|I).
p(ht|h′t, I)=p(h′t|ht, I)p(ht|I)/p(h′t|I),

where I is the side information of the estimator. The probability p(ht|I) is the subjective prior distribution for ht, which is immaterial in orthodox probability.

To derive p(ht|I) based on prior statistical knowledge on ht, e.g. first order moments, Jaynes introduces the maximum entropy principle in the reference E. T. Jaynes, “Information Theory and Statistical Mechanics”, Physical Review, APS, vol. 106, no. 4, pp. 620-630, 1957.

The latter allows to assign a distribution to p(ht|I) which satisfies the statistical constraints given by I and which minimizes the risk of making unwanted assumptions, by maximizing the entropy of p(ht|I).

When it comes to update probability assignments, Caticha proposes an extension of the maximum entropy principle, namely the minimum cross entropy principle (ME) in A. Caticha, “Lectures on Probability, Entropy and Statistical Physics”, arXiv:0808.0012v1 [physics.data-an], 2008.

When p(ht|I1) has been assigned for some side information I1, and new cogent information I2 is later available, then the ME principle consists in assigning to p(ht|I2) the distribution

$p ⁡ ( h t ❘ I 2 ) = arg ⁢ ⁢ min q ⁢ ⁢ S ⁡ [ q , p ⁡ ( h t ❘ I 1 ) ] ⁢ ⁢ Where ( 6 ) S ⁡ [ q , p ] = ∫ ⁢ q ⁡ ( x ) ⁢ log ⁡ ( p ⁡ ( x ) q ⁡ ( x ) ) ⁢ ⅆ x ( 7 )$

It can therefore be derived a method based on a minimal update requirement, which in essence assigns to p(ht|I2) the unique distribution which minimizes the changes brought to p(ht|I1) while satisfying the new constraints given by I2.

In the following, additional side information on ht (which possibly varies over time) will come from new available pilots at later time positions.

3. dynamic Optimal Channel Estimation for OFDM with Perfect System Parameter Knowledge

Channel estimation is performed at different time instants tk, k ε .

Denote Ik the knowledge at time tk. Since memory restrictions impose to discard past, we decide here only to consider at time tk the last received pilot data symbols, the last assigned probability p(hk|Ik−1), hk=htk and the supposedly known time correlation λ between time t and the past channel hk−1. Assume prior assigned distribution p(hk−1|Ik−1) at time index k−1. We have in general

$p ⁡ ( h k ❘ y k , I k ) ( 8 ) = p ⁡ ( y k ❘ h k , I k ) · p ⁡ ( h k ❘ I k ) p ⁡ ( y k ❘ I k ) ( 9 ) = p ⁡ ( y k ❘ h k , I k ) · ∫ p ⁡ ( h k ❘ h k - 1 , I k ) · p ⁡ ( h k - 1 ❘ I k ) ⁢ ⅆ h k - 1 p ⁡ ( y k ❘ I k ) ( 10 )$

Using Caticha's ME principle and set p(hk−1|Ik) to the previous p(hk−1|yk−1, Ik−1). The reason lies in the minimal update principle: if no additional information is given in Ik, compared to Ik−1, then p(hk−1|yk−1, Ik−1) is the distribution q that minimizes the cross-entropy S[q, p(hk−1|yk−1, Ik−1)]2

Let us now perform a recursive reasoning over the channel estimates at time indexes k ε . Assume that p(hk−1|yk−1, Ik−1) is Gaussianz, (kk−1, Mk−1).

One may show that this implies p(hk|yk, Ik) is still Gaussian. This will therefore be denoted (kk, Mm−1)

We have

$⁢ p ⁡ ( h k ❘ y k , I k ) ( 11 ) ⁢ = α 1 ⁢ p ⁡ ( y k ❘ h k , I k ) · ∫ p ⁡ ( h k ❘ h k - 1 , I k ) · p ⁡ ( h k - 1 , I k ) ⁢ ⅆ h k - 1 ( 12 ) = lim Q ~ ⟶ Q ⁢ ⅇ ( h k - h k ′ ) H ⁢ P k σ 2 ⁢ ( h k - h k ′ ) ⁢ ∫ ⅇ ( h k - λ ⁢ ⁢ h k - 1 ) H ⁢ Q ~ - 1 1 - λ 2 ⁢ ( h k - λ ⁢ ⁢ h k - 1 ) × α 2 ⁢ ⅇ ( h k - 1 - k k - 1 ) H ⁢ M k - 1 - 1 ⁡ ( h k - 1 - k k - 1 ) ⁢ ⅆ h k - 1 ( 13 )$

where the {tilde over (Q)}'s are taken from a set of invertible matrices in the neighborhood of Q, and the αi's are constants.

First we need to write the exponents of the Gaussian products in the integrand in a single Gaussian exponent form of the vector hk−1 times a constant independent of hk−1. By expansion and simplification, this is

$( h k - λ ⁢ ⁢ h k - 1 ) H ⁢ Q ~ - 1 1 - λ 2 ⁢ ( h k - λ ⁢ ⁢ h k - 1 ) + ( h k - 1 - k k - 1 ) H ⁢ M k - 1 - 1 ⁡ ( h k - 1 - k k - 1 ) = ( h k - 1 - 1 ) H ⁢ N ⁡ ( h k - 1 - 1 ) + C ⁡ ( h k ) ⁢ ⁢ With ⁢ ⁢ { N = λ 2 ⁢ Q _ - 1 1 - λ 2 + M k - 1 - 1 1 = N - 1 ⁡ ( λ 1 - λ 2 ⁢ Q ~ - 1 ⁢ h k + M k - 1 - 1 ⁢ k k - 1 ) C = h k H ⁢ Q _ - 1 1 - λ 2 ⁢ h k + k k - 1 H ⁢ M k - 1 - 1 ⁢ k k - 1 - ( h k + 1 - λ 2 λ ⁢ Q ~ ⁢ M k - 1 - 1 ) H ⁢ λ 2 ( 1 - λ 2 ) 2 ⁢ Q ~ - 1 ( λ 2 ⁢ Q _ - 1 1 - λ 2 + M k - 1 - 1 ) - 1 ⁢ Q ~ - 1 ( h k + 1 - λ 2 λ ⁢ Q ~ ⁢ M k - 1 - 1 ⁢ k k - 1 ) ( 14 )$

The integral (13) is then a constant times eC, which depends on hk. The term C must then be written again into a quadratic expression of hk. This is

$C = ( h k - j ) H ⁢ R ⁡ ( h k - j ) + B ⁢ ⁢ with ( 15 ) { R = ( λ 2 ⁢ M k - 1 + ( 1 - λ 2 ) ⁢ Q ~ ) - 1 j = λ ⁢ ⁢ k k - 1 B = 0 ( 16 )$

Together with the term outside the integral (13), this is
p(hk|yk, Ik)=α·e(hk−kk)HMk−1(hk−kk)   (17)

With α=(∫p(hk|yk, Ik)dhk)−1.

Finally, after some arithmetic derivation, in the limit {tilde over (Q)}→Q,

${ M k = λ 2 ( M k - 1 + 1 - λ 2 λ 2 ⁢ Q ) × ( λ 2 σ 2 ⁢ P k ( M k - 1 + 1 - λ 2 λ 2 ⁢ Q ) + I N ) - 1 k k = λ ⁢ ⁢ k k - 1 + 1 σ 2 ⁢ M k ⁢ P k ⁡ ( h k ′ - λ ⁢ ⁢ k k - 1 ) ( 18 )$

And the MMSE estimator ĥk for the channel at time index k is the first order moment of a Gaussian distribution centered in kk which is ĥk=kk. At initial time instant t0, if nothing but the channel delay spread L is known, M0=Q from MaxEnt, and k0=0. Therefore, we prove by the above recursion that, under this state of initial knowledge, for all k ε p(hk|Ik) is Gaussian with mean kk and variance Mk, and ĥk=kk. Note that, while regularized inverses of Q were used along the derivations, the final formulas are properly conditioned with respect to Q

The process which will be described hereafter will take advantage of this particular observation, and will achieve an iterative mechanism of channel estimation, based on an iterative update of both the mean kk and variance Mk of the distribution of the channel estimation.

Before entering into the particular detail of the process, the inventors have further discovered that the channel estimation can even be achieved without perfect system parameter knowledge 4. Dynamic Optimal Channel Estimation for OFDM with Perfect System Parameter Knowledge

In practical applications, contrary to what was stated above, the different parameters λ, σ2 and L especially, might not be perfectly known. For simplicity we assume those parameters are constant over the duration of the channel estimation process. Following the MaxEnt principle, these parameters must be assigned an a priori distribution. Let us focus on the time correlation λ, which is typically the most difficult parameter to track. In this respect, one has
p(hk|yk, Ik)=∫p(hk|yk, λIk)p(λ|yk, Ik)  (19)

Since yk cannot bring alone any cogent information on λ, p(λ|yk, Ik)=p(λ|Ik). The probability p(hk|yk, λ, Ik) was computed above and is given by the right-hand side of Equation (17), in which a depends on A and must therefore be made explicit.

p(hk|yk, Ik)=∫p(λ|Ik)α(λ)e(hk−kk(λ))H(Mk(λ))−1(hk−kk(λ))  (20)

with Mk(λ) and kk(λ) given by Equation (18) above for the λ in question and

$α ⁡ ( λ ) = βⅇ - x ⁡ ( λ ) ⁢ det ⁡ [ X ⁡ ( λ ) ] ⁢ ⁢ With ( 21 ) { X ⁡ ( λ ) = ( I + P k σ 2 ⁢ ( λ 2 ⁢ M k - 1 ( λ ) + ( 1 - λ 2 ) ⁢ Q ) ) - 1 x ⁡ ( λ ) = ( λ ⁢ ⁢ k k - 1 ( λ ) - h k ′ ) H ⁢ X ⁡ ( λ ) ⁢ P k σ 2 ⁢ ( λ ⁢ ⁢ k k - 1 ( λ ) - h k ′ ) ⁢ ⁢ And ⁢ ⁢ β = ( ∫ p ⁡ ( h k ❘ y k , I k ) ⁢ ⅆ ⋀ ⁢ h k ) - 1 , independent ⁢ ⁢ of ⁢ ⁢ λ . ( 22 )$

The MMSE estimate ĥk is then the weighted sum:

$h ^ k = ∫ p ⁡ ( λ ❘ I k ) ⁢ ⅇ - x ⁡ ( λ ) ⁢ det ⁡ [ X ⁡ ( λ ) ] ⁢ k k ( λ ) ⁢ ⅆ λ ∫ p ⁡ ( λ ❘ I k ) ⁢ ⅇ - x ⁡ ( λ ) ⁢ det ⁡ [ X ⁡ ( λ ) ] ⁢ ⅆ λ ( 23 )$

This integral is however very involved. In practice, it must be broken into a finite sum over a set of potential values for λ. Denoting S this set and |S| its cardinality, the recursive algorithm that provides the successive estimates ĥnk k=1, . . . , K, requires that at every step, the values for Mk(λ) and kk(λ), λS are kept in memory.

One may then describe the preferred embodiment of a channel estimation process which is based on the underlying above described principle.

II. Preferred Embodiment of the Channel Estimation Process for a OFDM Communication System

It has been demonstrated how the knowledge of the variance and estimate Mk(λ) and kk(λ) at k−1 can be advantageously used for achieving an iterative estimation method based on the update mechanism shown in formulas (18). With the hypothesis of a Gaussian channel, the information collected by the estimate at instant k−1 and the variance shows to be minimal and, therefore, allows to reduce the quantity of data to be stored within the memory since the variance will only need to be represented by an N×N matrix.

Furthermore, while the processes described below will be discussed in the frequency domain, where a matrix shows to be an N×N matrix which might have a significant size, the computations below may also be performed in the time domain what will result in a significant reduction of the size of the matrix.

With respect to FIG. 2, there is now described one embodiment of a process which assumes the knowledge of the value of λ

In a step 21, the process starts with an initialization step. At time t=0, the values of M0 and k0 are respectively set to Q and 0.

Then, in a step 22, the process enters into a loop;

Then, in a step 23, the process receives an OFDM symbol and may extract the received signal and the pilots therein contained.

Then, in a step 24, the process computes an estimate value of the variance of the probability distribution of p(hk|yk, Ik) at time instant k, yk, and Ik respectively corresponding to the received signal and to a priori known information, in accordance with the following formula:

$M k = λ 2 ( M k - 1 + 1 - λ 2 λ 2 ⁢ Q ) × ( λ 2 σ 2 ⁢ P k ( M k - 1 + 1 - λ 2 λ 2 ⁢ Q ) + I N ) - 1$
With

• N being the number of sub carriers in the OFDM symbol;
• kk being the estimated channel at instant k;
• Mk being a N×N matrix representative of the variance of the estimated channel at instant k;
• Q being a matrix initialized to a predetermined value;
• Pk is an N×N diagonal matrix comprising a “1” at the location of a pilot symbol and a “0” elsewhere;
• IN being the N×N identity matrix;
• λ corresponding to the known channel time correlation between instant k−1 and k.

Then, in a step 25the process computes an estimate value of the channel kk in accordance with the formula:

$k k = λ ⁢ ⁢ k k - 1 + 1 σ 2 ⁢ M k ⁢ P k ⁡ ( h k ′ - λ ⁢ ⁢ k k - 1 )$
With

• σ2 being the variance of the noise;
• h′k being equal to the received signal yk pointwise divided by the corresponding pilot sk;

Then, in a step 26, the process issues the estimate value ĥk=kk

Then, in a step 27, the process loops back to step 22 in order to process the next OFDM symbol.

With respect to FIG. 3, there is illustrated a second embodiment of a process of channel estimation which takes no assumption on the value of λ.

In a step 31, the process starts with an initialization step similar to step 21 above.

In a step 32 the process enters a loop.

In a step 33, the process receives an OFDM sample.

Then, in step 34, the process enters a new loop for processing a predetermined set of possible coefficients number of values (λ1 λ2 λ3 . . . λn . . . λM) with n=1 to M. To each value λn one associates a couple (Mk(λ), kk(λ), which will be calculated in accordance with the same formula as above.

Indeed, in steps 35 and 36 the process successively computes Mk(λ) and kk(λ) in accordance with the formulas:

${ M k ( λ ) = λ 2 ( M k - 1 ( λ ) + 1 - λ 2 λ 2 ⁢ Q ) × ( λ 2 σ 2 ⁢ P k ( M k - 1 ( λ ) + 1 - λ 2 λ 2 ⁢ Q ) + I N ) - 1 k k ( λ ) = λ ⁢ ⁢ k k - 1 ( λ ) + 1 σ 2 ⁢ M k ( λ ) ⁢ P k ⁡ ( h k ′ - λ ⁢ ⁢ k k - 1 ( λ ) )$
With

• N being the number of sub carriers in the OFDM symbol;
• kk being the estimated channel at instant k;
• Mk being a N×N matrix representative of the variance of the estimated channel at instant k;
• Q being a matrix initialized to a predetermined value;
• Pk is an N×N diagonal matrix comprising a “1” at the location of a pilot symbol and a “0” elsewhere;
• IN being the N×N identity matrix;
• λ corresponding to the known channel time correlation between instant k−1 and k.

Then, in a step 37, a test is considered in order to determine whether a subsequent value λn+1 is to be processed.

If so, the process returns to step 34 above for the purpose of computing that new value and, conversely, if all the values of (λ1 λ2 λ3 . . . λn . . . λM) have been processed, then the process proceeds to step 38 where it computes the estimate ĥk of the channel from a weighted sum of the latter, as follows:

$h ^ k = ∑ n = 1 N ⁢ ⅇ - x ⁡ ( λ n ) ⁢ det ⁡ [ X ⁡ ( λ n ) ] ⁢ k k ( λ n ) 1 N ⁢ ∑ n = 1 N ⁢ ⅇ - x ⁡ ( λ n ) ⁢ det ⁡ [ X ⁡ ( λ n ) ]$ $With$ ${ X ⁡ ( λ ) = ( I + P k σ 2 ⁢ ( λ 2 ⁢ M k - 1 ( λ ) + ( 1 - λ 2 ) ⁢ Q ) ) - 1 x ⁡ ( λ ) = ( λ ⁢ ⁢ k k - 1 ( λ ) - h k ′ ) H ⁢ X ⁡ ( λ ) ⁢ P k σ 2 ⁢ ( λ ⁢ ⁢ k k - 1 ( λ ) - h k ′ )$
With kk(λ) being the estimated channel at instant k depending on λ,

Mk(λ) being a matrix representative of the variance of the estimated channel at instant k assuming λ time correlation.

The process then proceeds to step 39 which loops back to step 32 for the purpose of processing a new OFDM symbol.

Simulations and Results

In this section, we provide simulation plots to compare, at time t, the minimal channel estimation update method against (i) the one-dimensional MMSE, taken at time t, which takes only into account the last past pilot symbols and uses a fixed empirical covariance matrix, (ii) an typical optimum two-dimensional MMSE, with K=4 pilot time indexes. The OFDM DFT size is N=64, the channel length L=6 is known to the receiver, the vehicular speed is v=120 km/h, pilot sequences are transmitted every 0.29 ms (as in 3GPPLTE [7]), and the induced Jake's time correlation λ between t and the past pilot sequence arrival time is known to the receiver. In scenario (ii), all K past received pilot sequences and time correlations are perfectly known. The channel time correlation model is a K-order autoregressive model following the reference K. E. Baddour and, N. C. Beaulieu, “Autoregressive modeling for fading channel simulation,” IEEE Transaction on Wireless Communications, July 2005.

A performance comparison is proposed in FIG. 5. We notice here that the minimal update algorithm does not show significant performance decay compared to the optimal two-dimensional MMSE estimator, while the one-dimensional MMSE estimator, also relying on the last past pilot sequence, shows large performance impairment.

In FIG. 3, with the same assumptions as previously, we consider the hypothesis where the vehicular speed v is a priori known to be (with equal probabilities) either 5, 50, 120 km/h. The performance is compared against the optimal 2-D algorithm where v is known but erroneously estimated (v=5, 50, 100 km/h). It is observed that, again, even when λ, or equivalently v, is a priori unknown, the Bayesian minimal update framework manages to ideally recover the channel with no performance decay. On the opposite, when λ is erroneously estimated, the performance decay of the optimal estimator might be dramatic.

## Claims

1. A process for estimating a channel in an OFDM communication system, wherein an OFDM sample is based on a time-frequency grid comprising a set of pilots, the process comprising the steps of:
receiving, by a receiver, OFDM sample based subframes comprising a set of pilot signals;
performing an iterative loop for the purpose of computing (Mk, kk) respectively corresponding to the covariance matrix and the mean of the probability distribution of p(hk|yk, lk) at time instant k, yk, and lk respectively corresponding to the received signal and to a priori known information and hk being a true channel at time instant k; and
said value of (Mk, kk) being computed as a function of (Mk−1, kk−1) as well as the values of the pilots received at instant k.
receiving, by a receiver, OFDM sample based subframes comprising a set of pilot signals;
performing an iterative loop for the purpose of computing (Mk, kk) respectively corresponding to the covariance matrix and the mean of the probability distribution of p(hk|yk, lk) at time instant k, yk, and lk respectively corresponding to the received signal and to a priori known information and hk being a true channel at time instant k; and
said value of (Mk, kk) being computed as a function of (Mk−1, kk−1) as well as the values of the pilots received at instant k.
2. The process for estimating the channel in an OFDM communication system according to claim 1 further comprising the steps of:
receiving OFDM sample based subframes comprising a set of pilots signals;
performing an estimation of the channel based on the mean kk and variance Mk of the distribution of the true channel hk in accordance with the following formulation:
${ M k = λ 2 ( M k - 1 + 1 - λ 3 λ 2 ⁢ Q ) × ( λ 2 σ 2 ⁢ P k ( M k - 1 + 1 - λ 2 λ 2 ⁢ Q ) + I N ) - 1 k k = λ ⁢ ⁢ k k - 1 + 1 σ 2 ⁢ M k ⁢ P k ⁡ ( h k ′ - λ ⁢ ⁢ k k - 1 )$
with
N being the number of sub carriers in the OFDM symbol;
kk being the mean of the inferred probability distribution p(hk|yk, lk) at time instant k;
Mk being a N×N matrix representative of the variance of the probability distribution of p(hk|yk, lk) at time instant k;
Q being a matrix initialized to a predetermined value;
Pk is an N×N diagonal matrix comprising a “1” at the location of a pilot symbol and a “0”elsewhere;
IN being the N×N identity matrix;
λ corresponding to the known channel time correlation between instant k−1 and k;
σ2 being the variance of the noise; and
h′k being equal to the received signal yk pointwise divided by the corresponding pilot sk.
receiving OFDM sample based subframes comprising a set of pilots signals;
performing an estimation of the channel based on the mean kk and variance Mk of the distribution of the true channel hk in accordance with the following formulation:
N being the number of sub carriers in the OFDM symbol;
kk being the mean of the inferred probability distribution p(hk|yk, lk) at time instant k;
Mk being a N×N matrix representative of the variance of the probability distribution of p(hk|yk, lk) at time instant k;
Q being a matrix initialized to a predetermined value;
Pk is an N×N diagonal matrix comprising a “1” at the location of a pilot symbol and a “0”elsewhere;
IN being the N×N identity matrix;
λ corresponding to the known channel time correlation between instant k−1 and k;
σ2 being the variance of the noise; and
h′k being equal to the received signal yk pointwise divided by the corresponding pilot sk.
3. The process according to claim 2 wherein Q is an a priori frequency covariance matrix for an L-tap uniform channel.
4. The process according to claim 2 wherein the time correlation factor (λ) corresponds to the correlation factor in the auto regressive model of the first order of hk.
5. The process according to claim 1 wherein said parameter Mk is forwarded to the base station as an information representative of the quality of the channel estimation.
6. The process according to claim 1 adapted to the 3GPP Long Term Evolution (LTE) standard.
7. A process for estimating a channel in an OFDM communication system, wherein an OFDM sample is based on a time-frequency grid comprising a set of pilots, the process comprising the steps of:
receiving, by a receiver, OFDM sample based subframes comprising a set of pilots signals;
performing an iterative loop for the purpose of computing (Mk(λ), kk(λ)) respectively corresponding to the variance and the mean of the probability distribution of p(hk|yk, lk) at time instant k, yk, and lk respectively corresponding to the received signal and to a priori known information;
said value of (Mk(λ), kk(λ)) being computed as a function of (Mk−1(λ), kk−1(λ)) for λspanning a predetermined set of values {λn} with n=1 to M between true channel hk and hk−1 respectively at time instant k and k−1, as well as the values of the pilots received at instant k,
$h ^ k = ∑ n = 1 N ⁢ ⅇ - x ⁡ ( λ n ) ⁢ det ⁡ [ X ⁡ ( λ n ) ] ⁢ k k ( λ n ) 1 N ⁢ ∑ n = 1 N ⁢ ⅇ - x ⁡ ( λ n ) ⁢ det ⁡ [ X ⁡ ( λ n ) ] ⁢ ⁢ With ⁢ ⁢ { M k ( λ ) = λ 2 ( M k - 1 ( λ ) + 1 - λ 2 λ 2 ⁢ Q ) × ( λ 2 σ 2 ⁢ P k ( M k - 1 ( λ ) + 1 - λ 2 λ 2 ⁢ Q ) + I N ) - 1 k k ( λ ) = λ ⁢ ⁢ k k - 1 ( λ ) + 1 σ 2 ⁢ M k ( λ ) ⁢ P k ⁡ ( h k ′ - λ ⁢ ⁢ k k - 1 ( λ ) ) And ⁢ ⁢ { X ⁡ ( λ ) = ( I + P k σ 2 ⁢ ( λ 2 ⁢ M k - 1 ( λ ) + ( 1 - λ 2 ) ⁢ Q ) ) - 1 x ⁡ ( λ ) = ( λ ⁢ ⁢ k k - 1 ( λ ) - h k ′ ) H ⁢ X ⁡ ( λ ) ⁢ P k σ 2 ⁢ ( λ ⁢ ⁢ k k - 1 ( λ ) - h k ′ )$
with kk(λ) being the mean of the probability distribution of p(hk|yk, lk) at instant k depending on λ, . . .
Mk(λ) being a matrix representative of the variance of the estimated channel at instant k; depending on λ,
Q being a matrix initialized to a predetermined value;
Pk is a diagonal N×N matrix comprising a “1” at the location of a pilot symbol;
IN being the N×N identity matrix;
σ2 being the variance of the noise;
h′k being equal to the received signal yk pointwise divided by the corresponding pilot sk;
λ corresponding to the known channel time correlation between instant k−1 and k;
H is the Hermitian transpose; and
det[X] is the determinant of the matrix X.
receiving, by a receiver, OFDM sample based subframes comprising a set of pilots signals;
performing an iterative loop for the purpose of computing (Mk(λ), kk(λ)) respectively corresponding to the variance and the mean of the probability distribution of p(hk|yk, lk) at time instant k, yk, and lk respectively corresponding to the received signal and to a priori known information;
said value of (Mk(λ), kk(λ)) being computed as a function of (Mk−1(λ), kk−1(λ)) for λspanning a predetermined set of values {λn} with n=1 to M between true channel hk and hk−1 respectively at time instant k and k−1, as well as the values of the pilots received at instant k,
Mk(λ) being a matrix representative of the variance of the estimated channel at instant k; depending on λ,
Q being a matrix initialized to a predetermined value;
Pk is a diagonal N×N matrix comprising a “1” at the location of a pilot symbol;
IN being the N×N identity matrix;
σ2 being the variance of the noise;
h′k being equal to the received signal yk pointwise divided by the corresponding pilot sk;
λ corresponding to the known channel time correlation between instant k−1 and k;
H is the Hermitian transpose; and
det[X] is the determinant of the matrix X.
8. The process according to claim 7 wherein Q is an a priori frequency correlation matrix for an L-tap uniform channel.
9. The process according to claim 7 wherein the set of time correlation factors {λn} corresponds to hypothetical values for the time correlation parameter λ.
10. The process according to claim 7 wherein said parameter Mk is forwarded to the base station as an information representative of the quality of the channel estimation.
11. The process according to claim 7 wherein it is adapted to the Long Term Evolution (LTE) standard.
12. A receiver for an OFDM communication system comprising:
a means for receiving a OFDM sample based subframes comprising a set of pilots signals;
a means for performing an iterative loop for the purpose of computing (Mk , kk) respectively corresponding to the variance and the mean of the probability distribution of p(hk|yk, lk) at time instant k, yk, and lk respectively corresponding to the received signal and to a priori known information; and
said value of (Mk, kk) being computed as a function of (Mk−1, kk−1) as well as the values of the pilots received at instant k, with hk corresponding to the channel at the kth sampled time instant.
a means for receiving a OFDM sample based subframes comprising a set of pilots signals;
a means for performing an iterative loop for the purpose of computing (Mk , kk) respectively corresponding to the variance and the mean of the probability distribution of p(hk|yk, lk) at time instant k, yk, and lk respectively corresponding to the received signal and to a priori known information; and
said value of (Mk, kk) being computed as a function of (Mk−1, kk−1) as well as the values of the pilots received at instant k, with hk corresponding to the channel at the kth sampled time instant.
13. A receiver according to claim 12, said receiver being incorporated to a mobile communication device for an LTE communication system.
14. The receiver according to claim 12, further comprising:
a means for performing an estimation of the channel based on the mean kk and variance Mk of the distribution of the channel estimation in accordance with the following formulation:
${ M k = λ 2 ( M k - 1 + 1 - λ 2 λ 2 ⁢ Q ) × ( λ 2 σ 2 ⁢ P k ( M k - 1 + 1 - λ 2 λ 2 ⁢ Q ) + I N ) - 1 k k = λ ⁢ ⁢ k k - 1 + 1 σ 2 ⁢ M k ⁢ P k ⁡ ( h k ′ - λ ⁢ ⁢ k k - 1 )$
with
N being the number of sub carriers in the OFDM symbol;
kk being the mean of the probability distribution of p(hk|yk, lk) at time instant k;
Mk being a N×N matrix representative of the variance of the probability distribution of p(hk|yk, lk) at instant k;
Q being a matrix initialized to a predetermined value;
Pk is an N×N diagonal matrix comprising a “1” at the location of a pilot symbol and a “0”elsewhere;
IN being the N×N identity matrix; and
λ corresponding to the known channel time correlation between instant k−1 and k;
σ2 being the variance of the noise; and
h′k being equal to the received signal yk pointwise divided by the corresponding pilot sk.
a means for performing an estimation of the channel based on the mean kk and variance Mk of the distribution of the channel estimation in accordance with the following formulation:
N being the number of sub carriers in the OFDM symbol;
kk being the mean of the probability distribution of p(hk|yk, lk) at time instant k;
Mk being a N×N matrix representative of the variance of the probability distribution of p(hk|yk, lk) at instant k;
Q being a matrix initialized to a predetermined value;
Pk is an N×N diagonal matrix comprising a “1” at the location of a pilot symbol and a “0”elsewhere;
IN being the N×N identity matrix; and
λ corresponding to the known channel time correlation between instant k−1 and k;
σ2 being the variance of the noise; and
h′k being equal to the received signal yk pointwise divided by the corresponding pilot sk.
15. A receiver for an OFDM communication system comprising a means for estimating the channel and characterized in that it further comprises:
a means for receiving a OFDM sample based subframes comprising a set of pilots signals;
a means for performing an iterative loop for the purpose of computing (Mk(λ), kk(λ)) respectively corresponding to the variance and the mean of the probability distribution of p(hk|yk, lk) at time instant k, yk, and lk respectively corresponding to the received signal and to a priori known information;
said value of (Mk(λ), kk(λ)) being computed as a function of (Mk−1, kk−1) for a predetermined set of values of time correlation (λn) with n=1 to M between the true channels hk and hk−1, with hk corresponding to the channel at the kth sampled time instant,as well as the values of the pilots received at instant k,
$h ^ k = ∑ n = 1 N ⁢ ⅇ - x ⁡ ( λ n ) ⁢ det ⁡ [ X ⁡ ( λ n ) ] ⁢ k k ( λ n ) 1 N ⁢ ∑ n = 1 N ⁢ ⅇ - x ⁡ ( λ n ) ⁢ det ⁡ [ X ⁡ ( λ n ) ] ⁢ ⁢ with ⁢ ⁢ { M k ( λ ) = λ 2 ( M k - 1 ( λ ) + 1 - λ 2 λ 2 ⁢ Q ) × ( λ 2 σ 2 ⁢ P k ( M k - 1 ( λ ) + 1 - λ 2 λ 2 ⁢ Q ) + I N ) - 1 k k ( λ ) = λ ⁢ ⁢ k k - 1 ( λ ) + 1 σ 2 ⁢ M k ( λ ) ⁢ P k ⁡ ( h k ′ - λ ⁢ ⁢ k k - 1 ( λ ) ) and ⁢ ⁢ { X ⁡ ( λ ) = ( I + P k σ 2 ⁢ ( λ 2 ⁢ M k - 1 ( λ ) + ( 1 - λ 2 ) ⁢ Q ) ) - 1 x ⁡ ( λ ) = ( λ ⁢ ⁢ k k - 1 ( λ ) - h k ′ ) H ⁢ X ⁡ ( λ ) ⁢ P k σ 2 ⁢ ( λ ⁢ ⁢ k k - 1 ( λ ) - h k ′ )$
With kk(λ) being the mean of the probability distribution of p(hk|yk, lk) at instant k depending on λ,
Mk(λ) being a matrix representative of the variance of the estimated channel at instant k; depending on λ,
Q being a matrix initialized to a predetermined value;
Pk is a diagonal N×N matrix comprising a “1” at the location of a pilot symbol;
IN being the N×N identity matrix;
σ2 being the variance of the noise; and
h′k being equal to the received signal yk pointwise divided by the corresponding pilot sk;
λ corresponding to the known channel time correlation between instant k−1 and k;
H is the Hermitian transpose; and
det[X] is the determinant of the matrix X.
a means for receiving a OFDM sample based subframes comprising a set of pilots signals;
a means for performing an iterative loop for the purpose of computing (Mk(λ), kk(λ)) respectively corresponding to the variance and the mean of the probability distribution of p(hk|yk, lk) at time instant k, yk, and lk respectively corresponding to the received signal and to a priori known information;
said value of (Mk(λ), kk(λ)) being computed as a function of (Mk−1, kk−1) for a predetermined set of values of time correlation (λn) with n=1 to M between the true channels hk and hk−1, with hk corresponding to the channel at the kth sampled time instant,as well as the values of the pilots received at instant k,
Mk(λ) being a matrix representative of the variance of the estimated channel at instant k; depending on λ,
Q being a matrix initialized to a predetermined value;
Pk is a diagonal N×N matrix comprising a “1” at the location of a pilot symbol;
IN being the N×N identity matrix;
σ2 being the variance of the noise; and
h′k being equal to the received signal yk pointwise divided by the corresponding pilot sk;
λ corresponding to the known channel time correlation between instant k−1 and k;
H is the Hermitian transpose; and
det[X] is the determinant of the matrix X.