Signal-to-interference + noise ratio estimator and method, mobile terminal having this estimator

Information

  • Patent Grant
  • 8310944
  • Patent Number
    8,310,944
  • Date Filed
    Tuesday, October 10, 2006
    17 years ago
  • Date Issued
    Tuesday, November 13, 2012
    11 years ago
Abstract
The SINR (signal-to-interference+noise ratio) estimator comprises a low-pass filter (66) to filter a variance σκ computed from the received pilot symbol amplitude to obtain an estimated variance σκ with a reduced bias, wherein the low-pass filter has an adjustable coefficient and the estimator comprises an electronic controller able to adjust the value of the adjustable coefficient according to the number of pilot symbols received during a timeslot.
Description
FIELD OF THE INVENTION

The present invention relates to signal-to-interference+noise ratio estimation.


BACKGROUND OF THE INVENTION

Signal-to-Interference+Noise Ratio (SINR) is an important metric of communication link quality. SINR estimation is of particular importance for wireless data systems where resources are shared dynamically amongst users. Some applications of SINR estimates are: a) Power Control in CDMA (Code Division Multiple Access) Systems: the receiver estimates the SINR, compares it to a target and commands the transmitter to increase/decrease its transmitted power; and b) Rate Adaptation: the information bit-rate assigned to a user can be dynamically varied based on its link quality and the system load. While such adaptation has limited use in voice systems, it is extremely useful for wireless data systems. Consequently, inaccurate SINR estimates can severely degrade performance and resource utilization.


Typically, the received signal corresponding to the jth demodulated transmitted pilot symbol in a kth timeslot is defined as

Ykj=akjμk+Ekj  (1)


j=1, 2, . . . , NP


where μk represents the received signal amplitude (product of transmitted amplitude and channel gain), Ekj is a random variable that represents the noise+interference, akj represents the demodulated pilot symbol-value, and Np is the number of pilot symbols received during the timeslot. Pilot symbol values can be +1 or −1 (in BPSK—Binary Phase Shift Keying), while it is assumed (without any loss of generality) that demodulated pilot symbol values are always +1. It is also assumed that the distribution that characterizes the noise+interference is Gaussian with zero mean and variance σ2. The SINR in the kth timeslot is then defined as:










θ
k

=


μ
K
2


σ
2






(
2
)








and is the parameter to be estimated.


US 2003/0016740 to Jeske et al. proposes a SINR estimator that improves the SINR estimation accuracy. More precisely, Jeske et al. propose a SINR estimator that smoothes the variance σK of the pilot symbol amplitude received during a kth timeslot to obtain an estimated variance {circumflex over (σ)}k with a reduced bias. The use of the estimated variance {circumflex over (σ)}k for the computation of SINR instead of variance σK increases the SINR accuracy. However, Jeske's estimator can still be improved to further increase the accuracy of the estimated SINR.


SUMMARY OF THE INVENTION

Accordingly, it is an object of the invention to provide a SINR estimator with a better accuracy.


With the foregoing and other objectives in view there is provided, in accordance with the invention, an estimator to estimate a signal-to-interference+noise ratio (SINR) from the amplitude of known pilot symbols received during a timeslot, the estimator comprising a low-pass filter to filter a variance σK computed from the received pilot symbol amplitudes to obtain an estimated variance {circumflex over (σ)}K with a reduced bias, wherein the low-pass filter has an adjustable coefficient and the estimator comprises an electronic controller able to adjust the value of the adjustable coefficient according to the number of pilot symbols received during the timeslot.


It has been observed that to further reduce the bias of the estimated variance {circumflex over (σ)}k, it is suited to dynamically adjust at least one coefficient of the low-pass filter according to the number NP of the pilot symbols received during the timeslot. In particular, this improves the operation of such estimator when the number NP is as low as two pilots per timeslot.


The embodiments of the above estimator may comprise one or several of the following features:

    • The low-pass filter comprises a feedback loop so that the estimated variance {circumflex over (σ)}K for the kth timeslot depends also on an estimated variance {circumflex over (σ)}K−1 of the previous (k−1)th timeslot;
    • The low-pass filter is able to obtain the estimated variance {circumflex over (σ)}K which is the result of the following relation:

      {circumflex over (σ)}K=aσK+b{circumflex over (σ)}K−1

      Herein:
    • σK is the variance computed only from pilot symbols received during the kth timeslot,
    • {circumflex over (σ)}K−1 is the estimated variance for the (k−1)th timeslot,
    • a is a predetermined coefficient, and
    • b is the adjustable coefficient.


      Coefficient a may be equal to a predetermined value α and the coefficient b value adjusted such that:






b
=

1
-
α
+

α

N
P








Herein, NP is the number of pilot symbols received during the kth timeslot.

    • The low-pass filter is an IIR (Infinite Impulse Response) filter.
    • The controller is able to receive the number NP of pilot symbols from a wireless emitter of a wireless telecommunication system.


The above embodiments of the estimator present the following advantages:

    • using a low-pass filter having a feedback loop further reduces the estimated variance bias because the variance σK of the pilot symbol amplitude cannot vary rapidly from one timeslot to the next one;
    • obtaining the estimated variance {circumflex over (σ)}k from the relation with coefficient a and b limits the number of computations required to obtain such an estimation; and
    • using an IIR filter further reduces the estimated variance bias.


The invention also relates to a mobile terminal comprising the above estimator.


The invention also relates to a method to estimate a signal-to-interference+noise ratio (SINR) from the amplitude of known pilot symbols received during a timeslot, the method comprising:

    • a step of filtering, using a low-pass filter, a variance σK computed from the received pilot symbol amplitudes to obtain an estimated variance {circumflex over (σ)}K with a reduced bias, and
    • a step of automatically adjusting a coefficient of the low-pass filter according to the number NP of pilot symbols received during the timeslot.


The embodiment of the above method may comprise the step of receiving the number NP of pilot symbols from a wireless emitter of a wireless telecommunication system.


These and other aspects of the invention will be apparent from the following description, drawings and claims.





BRIEF DESCRIPTION OF THE DRAWINGS


FIG. 1 is a schematic diagram of a part of a wireless communication system;



FIG. 2 is a schematic diagram of a SINR estimator used in the system of FIG. 1;



FIG. 3 is a schematic diagram of a low-pass filter used in the estimator of FIG. 2; and



FIG. 4 is a flowchart of a method to estimate a SINR.





DESCRIPTION OF THE PREFERRED EMBODIMENTS

To simplify the description and the understanding of the described estimator, the following description of embodiments will assume that:

    • only Binary Phase Shift Keying (BPSK) modulation is considered although the following disclosure and related analysis can be extended to other signalling schemes;
    • noise and interference are modelled together as Additive White Gaussian Noise (AWGN), but as will be appreciated from the following disclosure, this should not limit the application of the following estimation method;
    • transmission is organized into fixed duration timeslots, each containing known pilot symbols and data symbols;
    • the channel attenuation and phase shift are assumed fixed over a timeslot, and are treated as unknown constants rather than as random variables (no fading assumption).



FIG. 1 shows a part of a wireless communication system 2. For example, system 2 is a CDMA (Code Division Multiple Access) system or a TDMA (Time Division Multiple Access) system. For illustration purposes, system 2 is assumed to be a UMTS (Universal Mobile Telecommunications System).


In the following description, well-known functions or constructions known by a person of ordinary skill in the art are not described in detail.


System 2 has a base station 4 (also know as a node B) and a mobile terminal 6 such as a mobile phone.


Base station 4 has a wireless emitter 8 to transmit a radio signal 10.


Terminal 6 has a radio frequency receiver 12 connected to an antenna 14 to receive the radio signal 10 and to convert it into a baseband signal.


Terminal 6 has also a baseband processor 16 connected to receiver 12 to process the baseband signal.


A SINR estimator 18 is implemented as part of baseband processor 16, for example.


For example, radio signal 10 is divided into 10 ms frames, each frame being divided into 15 timeslots, each timeslot having a duration of 10/15 ms. Such timeslots are defined in the GSM (Global System for Mobile Communications) or UMTS (Universal Mobile Telecommunications System) standards. During each timeslot, Np pilot symbols are transmitted. Np may vary according to the transmission conditions. For example, Np is chosen from the following group of numbers: {2;4;8;16;32}.


Pilot symbols are predetermined symbols, the values of which are known by terminal 6 before they are received. More precisely, terminal 6 knows the expected value of each pilot symbol before this pilot symbol is received.



FIG. 2 shows in more details estimator 18.


Estimator 18 has an input 30 to receive demodulated pilot symbols. For example, each demodulated pilot symbol is computed by a multiplier 32 that multiplies the expected value of a pilot symbol by the received value for this pilot symbol.


Estimator 18 has also three outputs 34, 36 and 38 to output an estimation ŶK of the mean or average amplitude YK of the pilot symbol amplitude of the kth timeslot, the SINR estimation {circumflex over (θ)}k and the variance estimation {circumflex over (σ)}K, respectively.


Estimator 18 has an upper branch 40 that computes the mean amplitude estimation ŶK and a lower branch 42 that computes estimations {circumflex over (θ)}k and {circumflex over (σ)}K.


Branch 40 has an averaging module 44 that computes the average amplitude YK of the demodulated pilot symbol of the kth timeslot and a low-pass filter 46 that computes estimation ŶK from mean amplitude YK.


For example, averaging module 44 includes a summer 48 to sum up the amplitudes Ykj of all received demodulated pilot symbols received during the kth timeslot and a divider 50 to divide this sum by the number Np.


Filter 46 is an IIR filter with a feedback loop to take into account the estimation Ŷk−1.


Lower branch 42 has an energy calculator 52 to calculate the energy YK2 of the demodulated pilot symbols received during the kth timeslot and a subtracter 54 to establish variance σK.


Calculator 52 is able to calculate the energy of the demodulated pilot symbols according to the following relation:











Y
K
2

_

=


1

N
P




(




j
=
1


N
P




Y
Kj
2


)






(
3
)







For example, calculator 52 has the following units connected in series:

    • a squaring unit 56 to compute YKj2 from the demodulated pilot symbol amplitude Ykj;
    • a summer 58 to calculate the sum







(




j
=
1


N
P




Y
Kj
2


)

;





and

    • a divider 60 to divide the sum calculated by summer 58 by the number Np of pilot symbols received during the kth timeslot.


Subtraction unit 54 computes variance σK according to the following relation:

σK= YK2YK2  (4)


Herein, YK2 is the square of mean amplitude YK.


Estimator 18 has a squaring unit 62 to calculate YK2 from the mean amplitude YK output by averaging module 44.


Subtraction unit 54 has one input connected to an output of calculator 52, one input connected to an output of squaring unit 62 and one output to output variance σK.


Lower branch 42 also has an adaptive IIR low-pass filter 66 and a divider 68.


Filter 66 has an input connected to the output of subtraction unit 54 to receive variance σK and an output to output the estimated variance {circumflex over (σ)}K.


Filter 66 has at least one adjustable coefficient and an input 70 to receive instructions to adjust the value of the adjustable coefficient.


The output of filter 66 is connected to output 38 of estimator 18.


Filter 66 is described in more detail in FIG. 3.


Divider 68 is designed to calculate the SINR estimation {circumflex over (θ)}K according to the following relation:

{circumflex over (θ)}K=(ŶK)2/{circumflex over (σ)}K  (5)


To do so, divider 68 has two inputs: one of them is connected to the output of filter 66 and the other one is connected to the output of a squaring unit 74. Divider 68 has also an output connected to output 36 of estimator 18 to source the SINR estimation {circumflex over (θ)}K.


Squaring unit 74 is able to compute the square ŶK2 from the mean amplitude estimation ŶK outputted by filter 46.


Estimator 18 also has an electronic controller 78, which is able to adjust the value of the adjustable coefficient of filter 66 according to the value of number Np. Controller 78 has an output connected to input 70.



FIG. 3 shows filter 66 in more detail. Filter 66 has a feedback loop 80 to take into consideration the previously estimated variance {circumflex over (σ)}K−1 in the calculation of variance estimation {circumflex over (σ)}K.


For example, filter 66 obtains variance estimation {circumflex over (σ)}K from variance σk and estimation {circumflex over (σ)}K−1 according to the following relation:

{circumflex over (σ)}K=aσK+b{circumflex over (σ)}K−1  (6)

Herein:

    • a is a predetermined coefficient,
    • b is the adjustable coefficient, and
    • {circumflex over (σ)}K−1 is the variance estimation for the previous timeslot, i.e. the (K−1)th timeslot.


For example, a is equal to a constant value α. Coefficient b is given by the following relation:









b
=

1
-
α
+

α

N
P







(
7
)







In the illustrated example, filter 66 has:

    • a multiplier 82 to multiply variance σK with the predetermined constant value α,
    • a summer 84 to add the term α.σK to the term b.{circumflex over (σ)}k−1 of relation (6),
    • a multiplier 86 to multiply estimation {circumflex over (σ)}K−1 with the coefficient







1
-
α
+

α

N
P



,





and

    • a delay unit 88 to output estimation b.{circumflex over (σ)}k−1 to an input of summer 84.


The operation of estimator 18 will now be described with reference to FIG. 4.


Initially, during an initialization phase 100, base station 4 transmits, in step 102, the value of number NP to terminal 6. For example, this transmission takes place during the establishment of a communication link between terminal 6 and base station 4.


In step 104, controller 78 receives the value of number NP and automatically adjusts the value of coefficient b according to relation (7).


Once the value of coefficient b has been adjusted, during a phase 106, estimator 18 computes the variance estimation {circumflex over (θ)}K.


At the beginning of phase 106, multiplier 32, in step 108, demodulates the received pilot symbol and outputs to input 30 the demodulated pilot symbol Ykj.


In step 110, module 44 computes mean amplitude YK from the demodulated pilot symbols Ykj received during the kth timeslot.


In step 112, the mean amplitude YK is filtered by filter 46 to output the mean amplitude estimation ŶK. The mean amplitude estimation ŶK is outputted by estimator 18 in step 114.


In parallel to steps 110 to 114, in step 116, calculator 52 calculates the energy YK2.


Then, in step 118, subtraction unit 54 computes variance σK from the outputs of calculator 52 and squaring unit 62.


In step 120, filter 66 filters variance σK to obtain the variance estimation {circumflex over (σ)}K computed according to relation (6).


In step 122, the variance estimation {circumflex over (σ)}K is outputted by estimator 18.


In parallel to step 122, in step 124, divider 68 computes the SINR estimation {circumflex over (θ)}K according to relation (5).


In step 126, the SINR estimation {circumflex over (θ)}K is outputted by estimator 18.


Many additional embodiments are possible. For example, low-pass filter 66 may be replaced with a low-pass filter which is not an adaptive IIR filter but has at least one adjustable coefficient that can be adjusted by controller 78.


In the present specification and claims the word “a” or “an” preceding an element does not exclude the presence of a plurality of such elements. Further, the word “comprising” does not exclude the presence of other elements or steps than those listed.


The inclusion of reference signs in parentheses in the claims is intended to aid understanding and is not intended to be limiting.

Claims
  • 1. An estimator to estimate a signal-to-interference plus noise ratio SINR from the amplitude of known pilot symbols received from a base station within a wireless communication system, said pilot symbols being received during a timeslot according to {circumflex over (θ)}K=(ŶK)2/{circumflex over (σ)}K wherein {circumflex over (θ)}K is SINR estimation, and YK2 is a mean amplitude estimation of pilot symbols, the estimator comprising a low-pass filter to filter an interference plus noise variance σK computed from the received pilot symbol amplitudes to obtain an estimated variance {circumflex over (σ)}K with a reduced bias, wherein the low-pass filter has an adjustable coefficient and the estimator comprises an electronic controller able to adjust the value of the adjustable coefficient according to the number NP of received pilot symbols being received during the timeslot wherein the low-pass filter comprises a feedback loop so that the estimated variance {circumflex over (σ)}K for the kth timeslot depends also on an estimated variance {circumflex over (σ)}K−1 of the previous (k−1)th timeslot.
  • 2. The estimator according to claim 1, wherein the low-pass filter is able to obtain the estimated variance {circumflex over (σ)}K as the result of the following relation: {circumflex over (σ)}K=aσK+b{circumflex over (σ)}K−1
  • 3. The estimator according to claim 2, wherein coefficient a is equal to a predetermined value α and the coefficient b value is adjusted so that:
  • 4. The estimator according to claim 1, wherein the low-pass filter is an IIR (Infinite Impulse Response) filter.
  • 5. The estimator according to claim 1, wherein the controller is able to receive the number of Np of received pilot symbols from a wireless emitter of a wireless telecommunication system.
  • 6. Mobile terminal comprising an estimator according to claim 1 to estimate the SINR used in signal processing.
  • 7. A method to estimate a signal-to-interference plus noise ratio SINR from the amplitude of known pilot symbols received from a base station within a wireless communication system, said pilot symbols being received during a timeslot according to {circumflex over (θ)}K=(ŶK)2/{circumflex over (σ)}K wherein {circumflex over (σ)}K is SINR estimation, and YK2 is a mean amplitude estimation of pilot symbols, the method comprising: a step of filtering, using a low-pass filter, an interference plus noise variance σK computed from the amplitudes of the received pilot symbols to obtain an estimated variance {circumflex over (σ)}K with a reduced bias, anda step of automatically adjusting a coefficient of the low-pass filter according to the number Np of received pilot symbols being received during the timeslot wherein the low-pass filter comprises a feedback loop so that the estimated variance {circumflex over (σ)}K for the kth timeslot depends also on an estimated variance {circumflex over (σ)}K−1 of the previous (k−1)th timeslot.
  • 8. The method according to claim 7, wherein the method comprises the step of receiving the number Np of received pilot symbols from a wireless emitter of a wireless telecommunication system.
Priority Claims (2)
Number Date Country Kind
05300827 Oct 2005 EP regional
PCT/IB2006/053713 Oct 2006 WO international
PCT Information
Filing Document Filing Date Country Kind 371c Date
PCT/IB2006/053713 10/10/2006 WO 00 4/11/2008
Publishing Document Publishing Date Country Kind
WO2007/046028 4/26/2007 WO A
US Referenced Citations (8)
Number Name Date Kind
6335954 Bottomley et al. Jan 2002 B1
6970795 Burgett et al. Nov 2005 B1
7283510 Ito et al. Oct 2007 B2
20030016740 Jeske et al. Jan 2003 A1
20050163194 Gore et al. Jul 2005 A1
20050207585 Christoph Sep 2005 A1
20060023775 Rimini et al. Feb 2006 A1
20060133457 Wang et al. Jun 2006 A1
Foreign Referenced Citations (1)
Number Date Country
2005050865 Feb 2005 WO
Related Publications (1)
Number Date Country
20080247329 A1 Oct 2008 US