CN101599939B - Method and device for estimating reference signal received power of orthogonal frequency division multiplexing system - Google Patents

Method and device for estimating reference signal received power of orthogonal frequency division multiplexing system Download PDF

Info

Publication number
CN101599939B
CN101599939B CN2009100887531A CN200910088753A CN101599939B CN 101599939 B CN101599939 B CN 101599939B CN 2009100887531 A CN2009100887531 A CN 2009100887531A CN 200910088753 A CN200910088753 A CN 200910088753A CN 101599939 B CN101599939 B CN 101599939B
Authority
CN
China
Prior art keywords
mrow
rsrp
msub
math
munder
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN2009100887531A
Other languages
Chinese (zh)
Other versions
CN101599939A (en
Inventor
许百成
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing T3G Technology Co Ltd
Original Assignee
Beijing T3G Technology Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing T3G Technology Co Ltd filed Critical Beijing T3G Technology Co Ltd
Priority to CN2009100887531A priority Critical patent/CN101599939B/en
Publication of CN101599939A publication Critical patent/CN101599939A/en
Application granted granted Critical
Publication of CN101599939B publication Critical patent/CN101599939B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Monitoring And Testing Of Transmission In General (AREA)

Abstract

The invention provides a method and a device for estimating the reference signal received power (RSRP) of an orthogonal frequency division multiplexing (OFDM) system. The method comprises the following steps: A. carrying out the channel estimation of a received RS (reference signal); B. obtaining the corresponding RSRP through calculation according to the channel estimation results of all the RSsin the coherent block; C. perform the weighted average to the RSRPs corresponding to all the coherent blocks within the ranges of measuring frequency bandwidth and measuring time; D. carrying out thenoise power elimination on the RSRP obtained through the weighted average. The invention realizes the unbiased estimation of the RSRP.

Description

Method and device for estimating reference signal received power of orthogonal frequency division multiplexing system
Technical Field
The present invention relates to the field of wireless communications technologies, and in particular, to a method and an apparatus for estimating Reference Signal Received Power (RSRP) of an Orthogonal Frequency Division Multiplexing (OFDM) system.
Background
In a Long Term Evolution (LTE) wireless communication system, OFDM is the physical layer technology at its core. The data stream is modulated on a plurality of orthogonal subcarriers, and the frequency spectrums of the orthogonal subcarriers can be overlapped, so that the frequency spectrum utilization rate is greatly improved.
Referring to fig. 1, in an ofdm system, transmitted information may be described by a Resource grid, which represents resources of an entire time domain and a frequency domain in the form of Resource Elements (REs), i.e., a minimum square in fig. 1, which corresponds to a subcarrier in the frequency domain and a time length of a symbol (symbol) in the time domain. All information that needs to be transmitted is carried over resource elements. A plurality of resource elements constitute one resource block, and specifically, in the case of a normal cyclic prefix, 12 (the number of subcarriers) × 7 (the number of symbols) constitutes one resource block; in the case of an extended cyclic prefix, 12 (the number of subcarriers) × 6 (the number of symbols) constitutes one resource block.
In order to demodulate data correctly, the channel conditions at each RE location need to be known. In the LTE system, the channel estimation at each RE position is obtained as follows: and inserting a known Reference Signal (RS) at a specific resource element position of each resource block to calculate the channel estimation at the RE position, and then obtaining the channel estimation at all other resource element positions by using an interpolation method. Wherein the reference signal is a random sequence generated according to a specific rule.
RSRP is one of the standard measurement items in LTE systems, which is defined and required in the standard as follows: RSRP is the linear power average of the reference signal over the considered measurement frequency bandwidth and measurement time range. Reference signal R0(indicating that the RS signal sent by antenna port 0, carried by the RE labeled 0 in fig. 1) can be used to calculate RSRP if the terminal can accurately detect R1Is valid, then R1(watch)The RS signal sent by antenna port 1, carried by the RE labeled 1 in fig. 1) may also be used to calculate RSRP. If the terminal adopts the receiving diversity, the reported measurement value is not lower than the value of the RSRP corresponding to any branch in the diversity. However, no specific method of measuring RSRP is given in the standard.
Disclosure of Invention
The technical problem to be solved by the present invention is to provide a method and a device for estimating the reference signal received power of an orthogonal frequency division multiplexing system, so as to achieve unbiased estimation of the reference signal received power.
In order to solve the technical problems, the invention provides the following technical scheme:
a method for estimating reference signal received power of an orthogonal frequency division multiplexing system comprises the following steps:
A. performing channel estimation on the received RS;
B. calculating according to the channel estimation results of all RSs in the coherent block to obtain corresponding RSRP;
C. carrying out weighted average on the RSRP corresponding to all coherent blocks in the measurement frequency bandwidth and the measurement time range;
D. and carrying out noise power elimination on the RSRP obtained by the weighted average.
In the estimation method, in step B, RSRP is calculated by using the following formula:
<math><mrow><msup><mi>P</mi><mrow><mo>&prime;</mo><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow></mrow></msup><mo>=</mo><msup><mrow><mo>|</mo><mfrac><mn>1</mn><mi>KL</mi></mfrac><munder><mi>&Sigma;</mi><mi>k</mi></munder><munder><mi>&Sigma;</mi><mi>l</mi></munder><msubsup><mover><mi>h</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>.</mo><mi>l</mi></mrow><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow></msubsup><mo>|</mo></mrow><mn>2</mn></msup></mrow></math>
wherein, P'(i,j,n)The corresponding RSRP of the nth coherent block transmitted by the ith antenna port and received by the jth antenna is shown, K represents the number of subcarriers occupied by the RS in the coherent block, L represents the number of OFDM symbols occupied by the RS in the coherent block,
Figure G2009100887531D00022
represents the channel estimation of the kth RS signal in the nth coherent block, the kth RS signal in the frequency domain and the ith RS signal in the time domain, which are transmitted by the ith antenna port and received by the jth antenna, | · | represents the modulo operator.
In the estimation method, in step C, the weighted average is performed by using the following formula:
<math><mrow><mover><mi>P</mi><mo>&OverBar;</mo></mover><mo>=</mo><munder><mi>&Sigma;</mi><mi>i</mi></munder><msub><mi>&alpha;</mi><mi>i</mi></msub><munder><mi>&Sigma;</mi><mi>j</mi></munder><msub><mi>&beta;</mi><mi>j</mi></msub><mrow><mo>(</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munder><mi>&Sigma;</mi><mi>n</mi></munder><msup><mi>P</mi><mrow><mo>&prime;</mo><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow></mrow></msup><mo>)</mo></mrow><mo>,</mo><mn>1</mn><mo>&le;</mo><mi>i</mi><mo>&le;</mo><mi>I</mi><mo>,</mo><mn>1</mn><mo>&le;</mo><mi>j</mi><mo>&le;</mo><mi>J</mi></mrow></math>
wherein P represents P 'to all'(i,j,n)RSRP after weighted averaging, I denotes the total number of transmit antenna ports used to calculate RSRP, J denotes the total number of receive antennas used to calculate RSRP, and N denotes coherence within the measurement frequency bandwidth and measurement time rangeTotal number of blocks, αiAnd betajWeighting coefficients for the transmit antenna ports and the receive antennas respectively, <math><mrow><mn>0</mn><mo>&lt;</mo><msub><mi>&alpha;</mi><mi>i</mi></msub><mo>&lt;</mo><mn>1</mn><mo>,</mo><munder><mi>&Sigma;</mi><mi>i</mi></munder><msub><mi>&alpha;</mi><mi>i</mi></msub><mo>=</mo><mn>1</mn><mo>,</mo></mrow></math> <math><mrow><mn>0</mn><mo>&lt;</mo><msub><mi>&beta;</mi><mi>j</mi></msub><mo>&lt;</mo><mn>1</mn><mo>,</mo><munder><mi>&Sigma;</mi><mi>j</mi></munder><msub><mi>&beta;</mi><mi>j</mi></msub><mo>=</mo><mn>1</mn><mo>.</mo></mrow></math>
the above-mentioned estimation method, wherein, <math><mrow><msub><mi>&alpha;</mi><mi>i</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>I</mi></mfrac><mo>,</mo></mrow></math> <math><mrow><msub><mi>&beta;</mi><mi>j</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>J</mi></mfrac><mo>.</mo></mrow></math>
in the estimation method, in step D, the following formula is adopted to perform noise power cancellation:
<math><mrow><mover><mi>P</mi><mo>^</mo></mover><mo>=</mo><mover><mi>P</mi><mo>&OverBar;</mo></mover><mo>-</mo><mfrac><msup><mi>&sigma;</mi><mn>2</mn></msup><mi>KL</mi></mfrac></mrow></math>
wherein,
Figure G2009100887531D00034
represents RSRP, σ after eliminating noise power component contained in P2Is the noise power.
The estimation method further comprises the following steps: E. and smoothing the RSRP obtained by eliminating the noise power.
In the above estimation method, in step E, the following formula is adopted for smoothing:
<math><mrow><msub><mi>P</mi><mi>m</mi></msub><mo>=</mo><mi>&alpha;</mi><msub><mover><mi>P</mi><mo>^</mo></mover><mi>m</mi></msub><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&alpha;</mi><mo>)</mo></mrow><msub><mi>P</mi><mrow><mi>m</mi><mo>-</mo><mn>1</mn></mrow></msub></mrow></math>
wherein, PmFor the smoothed RSRP current value,
Figure G2009100887531D00036
is the current value of RSRP without smoothing, Pm-1Alpha is a smoothing factor, and alpha is more than 0 and less than or equal to 1.
The above estimation method, wherein the coherent block is a block formed by all REs in the coherence time range and the coherence bandwidth range.
In the above estimation method, the coherent block is a resource block.
An apparatus for estimating reference signal received power of an orthogonal frequency division multiplexing system, comprising:
an RS channel estimator for performing channel estimation on the received RS;
the RSRP calculator is used for calculating according to the channel estimation results of all RSs in the coherent block to obtain corresponding RSRP;
the RSRP averager is used for carrying out weighted average on the RSRPs corresponding to all coherent blocks in the measurement frequency bandwidth and the measurement time range;
and the noise power eliminator is used for eliminating the noise power of the RSRP obtained by the weighted average.
The above estimation device further includes: and the RSRP smoother is used for smoothing the RSRP obtained by eliminating the noise power.
The embodiment of the invention divides the channel estimation of the reference signal into a plurality of coherent blocks, carries out RSRP estimation on each coherent block, then carries out weighted average on the estimation results of the coherent blocks on all the transceiving antennas, and finally carries out noise power elimination on the average result, thereby realizing unbiased estimation on the RSRP.
Drawings
Fig. 1 is a schematic diagram of a distribution situation of a reference signal in a resource block under a condition of a normal cyclic prefix and a number of transmit antennas less than or equal to 2;
fig. 2 is a schematic structural diagram of an apparatus for estimating reference signal received power of an ofdm system according to an embodiment of the present invention;
fig. 3 is a flowchart of a method for estimating a reference signal received power of an ofdm system according to an embodiment of the present invention.
Detailed Description
As can be seen from the theory of wireless communication, reference signals within the coherence bandwidth and coherence time range may have the same or similar channel estimates. The invention is based on this principle to accomplish the estimation of RSRP. For simplicity of the following description, a concept of a coherent block is proposed, which is a block composed of all REs within a coherence time range and a coherence bandwidth range.
Referring to fig. 2, the apparatus for estimating reference signal received power of an ofdm system according to an embodiment of the present invention mainly includes: RS channel estimator, RSRP calculator, RSRP averager, noise power canceller, and RSRP smoother (optional).
Is as followsIn the description, i represents a serial number of a transmitting antenna port; j represents the receiving antenna serial number; n represents a coherent block number; k represents the sequence number of the RS in the frequency domain in the coherent block; l represents the sequence number of the RS in the time domain in the coherent block;
Figure G2009100887531D00041
the system comprises a first antenna port, a second antenna port, a first coherent block, a second coherent block and a third coherent block, wherein the first coherent block is transmitted by the first antenna port and received by the second antenna port;
Figure G2009100887531D00042
to represent
Figure G2009100887531D00043
Estimating a channel of (1); p'(i,j,n)Representing the estimated RSRP of all RSs in the nth coherent block transmitted by the ith antenna port and received by the jth antenna; p represents P'(i,j,n)Performing weighted average processing on the RSRP;
Figure G2009100887531D00044
an RSRP indicating the noise power component contained in the P is removed; p represents a pair
Figure G2009100887531D00045
RSRP after smoothing.
RS channel estimator
The RS channel estimator is configured to perform channel estimation on the received RS, and specifically includes:
h ^ k , l ( i , j , n ) = r ^ k , l ( i , j , n ) r k , l ( i , j , n ) *
wherein r isk,l (i,j,n)Indicating the transmitted reference signal, (.)*Representing the conjugate operator.
RSRP calculator
The RSRP calculator is configured to calculate according to channel estimation results of all RSs in the coherent block to obtain a corresponding RSRP, and specifically includes:
<math><mrow><msup><mi>P</mi><mrow><mo>&prime;</mo><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow></mrow></msup><mo>=</mo><msup><mrow><mo>|</mo><mfrac><mn>1</mn><mi>KL</mi></mfrac><munder><mi>&Sigma;</mi><mi>k</mi></munder><munder><mi>&Sigma;</mi><mi>l</mi></munder><msubsup><mover><mi>h</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>.</mo><mi>l</mi></mrow><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow></msubsup><mo>|</mo></mrow><mn>2</mn></msup></mrow></math>
wherein, K represents the number of sub-carriers occupied by the RS in the coherent block, L represents the number of OFDM symbols occupied by the RS in the coherent block, and | - | represents the modulus operator.
The coherent block can be divided by using a coherent bandwidth estimation module and a coherent time estimation module in the receiver, and as a conservative implementation method, if the coherent bandwidth estimation module and the coherent time estimation module do not exist in the receiver, one resource block can be regarded as one coherent block.
RSRP averager
The RSRP averager is configured to perform weighted averaging on RSRPs corresponding to all coherent blocks in a measurement frequency bandwidth and a measurement time range to reduce an estimated error, and specifically includes:
<math><mrow><mover><mi>P</mi><mo>&OverBar;</mo></mover><mo>=</mo><munder><mi>&Sigma;</mi><mi>i</mi></munder><msub><mi>&alpha;</mi><mi>i</mi></msub><munder><mi>&Sigma;</mi><mi>j</mi></munder><msub><mi>&beta;</mi><mi>j</mi></msub><mrow><mo>(</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munder><mi>&Sigma;</mi><mi>n</mi></munder><msup><mi>P</mi><mrow><mo>&prime;</mo><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow></mrow></msup><mo>)</mo></mrow><mo>,</mo><mn>1</mn><mo>&le;</mo><mi>i</mi><mo>&le;</mo><mi>I</mi><mo>,</mo><mn>1</mn><mo>&le;</mo><mi>j</mi><mo>&le;</mo><mi>J</mi></mrow></math>
where I denotes the total number of transmit antenna ports used to calculate RSRP (one antenna may have multiple ports, where the total number is the sum of the number of ports of all transmit antennas), J denotes the total number of receive antennas used to calculate RSRP, N denotes the total number of coherent blocks within the measurement frequency bandwidth and measurement time range considered, αiAnd betajWeighting coefficients of the transmitting antenna port and the receiving antenna respectively, satisfy <math><mrow><mn>0</mn><mo>&lt;</mo><msub><mi>&alpha;</mi><mi>i</mi></msub><mo>&lt;</mo><mn>1</mn><mo>,</mo><munder><mi>&Sigma;</mi><mi>i</mi></munder><msub><mi>&alpha;</mi><mi>i</mi></msub><mo>=</mo><mn>1</mn><mo>,</mo></mrow></math> <math><mrow><mn>0</mn><mo>&lt;</mo><msub><mi>&beta;</mi><mi>j</mi></msub><mo>&lt;</mo><mn>1</mn><mo>,</mo><munder><mi>&Sigma;</mi><mi>j</mi></munder><msub><mi>&beta;</mi><mi>j</mi></msub><mo>=</mo><mn>1</mn><mo>.</mo></mrow></math> As a simple embodiment, it is preferable <math><mrow><msub><mi>&alpha;</mi><mi>i</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>I</mi></mfrac><mo>,</mo></mrow></math> <math><mrow><msub><mi>&beta;</mi><mi>j</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>J</mi></mfrac><mo>.</mo></mrow></math>
Noise power canceller
And the noise power eliminator is used for eliminating the noise power of the RSRP obtained by the weighted average. The RSRP value calculated by the RSRP averager is not an unbiased estimation result, and includes a noise power with a certain component, and in order to implement the unbiased estimation of RSRP, the component of the noise power needs to be eliminated, specifically:
<math><mrow><mover><mi>P</mi><mo>^</mo></mover><mo>=</mo><mover><mi>P</mi><mo>&OverBar;</mo></mover><mo>-</mo><mfrac><msup><mi>&sigma;</mi><mn>2</mn></msup><mi>KL</mi></mfrac></mrow></math>
wherein σ2Is the noise power.
RSRP smoother
And the RSRP smoother is used for smoothing the RSRP obtained by eliminating the noise power. To improve the stability of RSRP estimation, as an optimization, the current RSRP estimate and the historical values may be smoothed, i.e.,
<math><mrow><msub><mi>P</mi><mi>m</mi></msub><mo>=</mo><mi>&alpha;</mi><msub><mover><mi>P</mi><mo>^</mo></mover><mi>m</mi></msub><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&alpha;</mi><mo>)</mo></mrow><msub><mi>P</mi><mrow><mi>m</mi><mo>-</mo><mn>1</mn></mrow></msub></mrow></math>
wherein, PmFor the smoothed RSRP current value,
Figure G2009100887531D00062
is the current value of RSRP without smoothing, Pm-1Alpha is a smoothing factor (can be taken as a value according to the change condition of a channel) for the smoothed RSRP historical value, and alpha is more than 0 and less than or equal to 1.
How to divide the resources within the measurement frequency bandwidth and the measurement time range into a plurality of coherent blocks is the prior art, and as an implementation manner of the prior art, reference may be made to the chinese patent application document, which is entitled "a channel estimation method and apparatus for an orthogonal frequency division multiplexing system" and has an application number of 200910082322.4 and an application date of 2009, 4.14, which is provided by the present applicant, and only the principle thereof is briefly introduced here. As described above, coherent blocks may be divided by using a coherent bandwidth estimation module and a coherent time estimation module in a receiver, and as a conservative implementation method, if the coherent bandwidth estimation module and the coherent time estimation module do not exist in the receiver, a resource block may be regarded as a coherent block.
By coherent block is meant a block consisting of all REs within the coherence bandwidth and within the coherence time. By coherence time is meant a time interval during which the amplitude of the received signal has a strong correlation, i.e. the channel impulse response remains substantially constant during the coherence time. By coherence bandwidth is meant a frequency domain bandwidth within which the amplitude of the received signal has a strong correlation, i.e. the amplitude-frequency response of the channel remains substantially constant within the coherence bandwidth. The channel delay spread and the coherent bandwidth have a strong corresponding relation, and the larger the channel delay spread is, the narrower the coherent bandwidth is; conversely, the smaller the channel delay spread, the wider the coherence bandwidth. The term "channel delay spread" refers to the difference of propagation paths experienced by each path component of a channel, and thus has different time delays, which causes the energy of a signal to be spread in time.
The coherent bandwidth estimation module firstly utilizes the frequency domain channel estimation of the reference signal to calculate the time delay expansion tau in the wireless channelnThen, the corresponding relation between the time delay expansion and the coherent bandwidth is carried outCalculating to obtain the coherent bandwidth <math><mrow><mi>B</mi><mo>=</mo><mfrac><mn>1</mn><mrow><mi>&zeta;</mi><msub><mi>&tau;</mi><mi>n</mi></msub></mrow></mfrac><mo>,</mo></mrow></math> Where ζ is a mapping factor of delay spread and coherence bandwidth. The delay spread has various expression forms, and the maximum delay spread, the average delay spread, the root mean square delay spread and the like are commonly used. In either representation, the effective paths of the time domain channel impulse response are first determined. The time delay expansion can be obtained by performing inverse Fourier transform on the frequency domain channel estimation of the position of the reference signal to obtain time domain channel estimation, then obtaining an effective diameter of the time domain channel estimation, and calculating according to the effective diameter to obtain the time delay expansion of the channel.
And the coherent time estimation module calculates according to the frequency domain channel estimation at different time domain positions to obtain the coherent time of the channel. Specifically, the coherence time estimation module may include: a difference factor calculation unit for calculating a difference factor between the frequency domain channel estimation of the reference signal and an available frequency domain channel estimation adjacent in the time domain; and the coherence time calculation unit is used for determining the coherence time of the channel according to the difference factor.
Referring to fig. 3, the method for estimating the reference signal received power of the ofdm system according to the embodiment of the present invention mainly includes the following steps:
step 301: performing channel estimation on the received RS, specifically:
h ^ k , l ( i , j , n ) = r ^ k , l ( i , j , n ) r k , l ( i , j , n ) *
wherein r isk,l (i,j,n)Indicating the transmitted reference signal, (.)*Representing the conjugate operator.
Step 302: calculating according to the channel estimation results of all the RSs in the coherent block to obtain the corresponding RSRP, specifically:
<math><mrow><msup><mi>P</mi><mrow><mo>&prime;</mo><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow></mrow></msup><mo>=</mo><msup><mrow><mo>|</mo><mfrac><mn>1</mn><mi>KL</mi></mfrac><munder><mi>&Sigma;</mi><mi>k</mi></munder><munder><mi>&Sigma;</mi><mi>l</mi></munder><msubsup><mover><mi>h</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>.</mo><mi>l</mi></mrow><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow></msubsup><mo>|</mo></mrow><mn>2</mn></msup></mrow></math>
wherein, K represents the number of sub-carriers occupied by the RS in the coherent block, L represents the number of OFDM symbols occupied by the RS in the coherent block, and | - | represents the modulus operator.
The coherent block can be divided by using a coherent bandwidth estimation module and a coherent time estimation module in the receiver, and as a conservative implementation method, if the coherent bandwidth estimation module and the coherent time estimation module do not exist in the receiver, one resource block can be regarded as one coherent block.
Step 303: performing weighted average on RSRP corresponding to all coherent blocks within a measurement frequency bandwidth and a measurement time range to reduce an estimation error, specifically:
<math><mrow><mover><mi>P</mi><mo>&OverBar;</mo></mover><mo>=</mo><munder><mi>&Sigma;</mi><mi>i</mi></munder><msub><mi>&alpha;</mi><mi>i</mi></msub><munder><mi>&Sigma;</mi><mi>j</mi></munder><msub><mi>&beta;</mi><mi>j</mi></msub><mrow><mo>(</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munder><mi>&Sigma;</mi><mi>n</mi></munder><msup><mi>P</mi><mrow><mo>&prime;</mo><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow></mrow></msup><mo>)</mo></mrow><mo>,</mo><mn>1</mn><mo>&le;</mo><mi>i</mi><mo>&le;</mo><mi>I</mi><mo>,</mo><mn>1</mn><mo>&le;</mo><mi>j</mi><mo>&le;</mo><mi>J</mi></mrow></math>
where I denotes the total number of transmit antenna ports used to calculate RSRP (one antenna may have multiple ports, where the total number is the sum of the number of ports of all transmit antennas), J denotes the total number of receive antennas used to calculate RSRP, N denotes the total number of coherent blocks within the measurement frequency bandwidth and measurement time range considered, αiAnd betajWeighting coefficients of the transmitting antenna port and the receiving antenna respectively, satisfy <math><mrow><mn>0</mn><mo>&lt;</mo><msub><mi>&alpha;</mi><mi>i</mi></msub><mo>&lt;</mo><mn>1</mn><mo>,</mo><munder><mi>&Sigma;</mi><mi>i</mi></munder><msub><mi>&alpha;</mi><mi>i</mi></msub><mo>=</mo><mn>1</mn><mo>,</mo></mrow></math> <math><mrow><mn>0</mn><mo>&lt;</mo><msub><mi>&beta;</mi><mi>j</mi></msub><mo>&lt;</mo><mn>1</mn><mo>,</mo><munder><mi>&Sigma;</mi><mi>j</mi></munder><msub><mi>&beta;</mi><mi>j</mi></msub><mo>=</mo><mn>1</mn><mo>.</mo></mrow></math> As a simple embodiment, it is preferable <math><mrow><msub><mi>&alpha;</mi><mi>i</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>I</mi></mfrac><mo>,</mo></mrow></math> <math><mrow><msub><mi>&beta;</mi><mi>j</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>J</mi></mfrac><mo>.</mo></mrow></math>
Step 304: and carrying out noise power elimination on the RSRP obtained by the weighted average. The RSRP value calculated in step 303 is not a result of unbiased estimation, and includes a noise power component, and to implement the unbiased estimation of RSRP, the noise power component needs to be eliminated, specifically:
<math><mrow><mover><mi>P</mi><mo>^</mo></mover><mo>=</mo><mover><mi>P</mi><mo>&OverBar;</mo></mover><mo>-</mo><mfrac><msup><mi>&sigma;</mi><mn>2</mn></msup><mi>KL</mi></mfrac></mrow></math>
wherein σ2Is the noise power.
Step 305: and smoothing the RSRP obtained by eliminating the noise power. To improve the stability of RSRP estimation, as an optimization, the current RSRP estimate and the historical values may be smoothed, i.e.,
<math><mrow><msub><mi>P</mi><mi>m</mi></msub><mo>=</mo><mi>&alpha;</mi><msub><mover><mi>P</mi><mo>^</mo></mover><mi>m</mi></msub><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&alpha;</mi><mo>)</mo></mrow><msub><mi>P</mi><mrow><mi>m</mi><mo>-</mo><mn>1</mn></mrow></msub></mrow></math>
wherein, PmFor the smoothed RSRP current value,
Figure G2009100887531D00086
is the current value of RSRP without smoothing, Pm-1Alpha is a smoothing factor (can be taken as a value according to the change condition of a channel) for the smoothed RSRP historical value, and alpha is more than 0 and less than or equal to 1.
In summary, the embodiments of the present invention implement unbiased estimation of RSRP by dividing the channel estimation of the reference signal into a plurality of coherent blocks, performing RSRP estimation on each coherent block, then performing weighted averaging on the estimation results of the coherent blocks on all the transceiving antennas, and finally performing noise power cancellation on the averaged result.
Finally, it should be noted that the above embodiments are only for illustrating the technical solutions of the present invention and not for limiting, and those skilled in the art should understand that the technical solutions of the present invention can be modified or substituted with equivalents without departing from the spirit scope of the technical solutions of the present invention, which should be covered by the scope of the claims of the present invention.

Claims (16)

1. A method for estimating reference signal received power of an orthogonal frequency division multiplexing system, comprising the steps of:
A. performing channel estimation on a received Reference Signal (RS);
B. calculating according to the channel estimation results of all RSs in the coherent block to obtain corresponding Reference Signal Received Power (RSRP);
C. carrying out weighted average on the RSRP corresponding to all coherent blocks in the measurement frequency bandwidth and the measurement time range;
D. carrying out noise power elimination on the RSRP obtained by weighted average;
in step B, RSRP is calculated using the following formula:
<math><mrow><msup><mi>P</mi><mrow><mo>&prime;</mo><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow></mrow></msup><mo>=</mo><msup><mrow><mo>|</mo><mfrac><mn>1</mn><mi>KL</mi></mfrac><munder><mi>&Sigma;</mi><mi>k</mi></munder><munder><mi>&Sigma;</mi><mi>l</mi></munder><msubsup><mover><mi>h</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>,</mo><mi>l</mi></mrow><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow></msubsup><mo>|</mo></mrow><mn>2</mn></msup></mrow></math>
wherein, P'(i,j,n)The corresponding RSRP of the nth coherent block transmitted by the ith antenna port and received by the jth antenna is shown, K represents the number of subcarriers occupied by the RS in the coherent block, L represents the number of orthogonal frequency division multiplexing OFDM symbols occupied by the RS in the coherent block,
Figure FSB00000525382500012
represents the channel estimation of the kth RS signal in the nth coherent block, the kth RS signal in the frequency domain and the ith RS signal in the time domain, which are transmitted by the ith antenna port and received by the jth antenna, | · | represents the modulo operator.
2. The estimation method according to claim 1, wherein in step C, the weighted average is performed using the following formula:
<math><mrow><mover><mi>P</mi><mo>&OverBar;</mo></mover><mo>=</mo><munder><mi>&Sigma;</mi><mi>i</mi></munder><msub><mi>&alpha;</mi><mi>i</mi></msub><munder><mi>&Sigma;</mi><mi>j</mi></munder><msub><mi>&beta;</mi><mi>j</mi></msub><mrow><mo>(</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munder><mi>&Sigma;</mi><mi>n</mi></munder><msup><mi>P</mi><mrow><mo>&prime;</mo><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow></mrow></msup><mo>)</mo></mrow><mo>,</mo><mn>1</mn><mo>&le;</mo><mi>i</mi><mo>&le;</mo><mi>I</mi><mo>,</mo><mn>1</mn><mo>&le;</mo><mi>j</mi><mo>&le;</mo><mi>J</mi></mrow></math>
wherein,
Figure FSB00000525382500014
represents to all P'(i,j,n)RSRP after weighted averaging, I denotes the total number of transmit antenna ports used to calculate RSRP, J denotes the total number of receive antennas used to calculate RSRP, N denotes the total number of coherent blocks within the measurement frequency bandwidth and measurement time range, αiAnd betajWeighting coefficients for the transmit antenna ports and the receive antennas respectively,
Figure FSB00000525382500015
3. the estimation method according to claim 2, characterized in that:
<math><mrow><msub><mi>&alpha;</mi><mi>i</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>I</mi></mfrac><mo>,</mo></mrow></math> <math><mrow><msub><mi>&beta;</mi><mi>j</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>J</mi></mfrac><mo>.</mo></mrow></math>
4. the estimation method according to claim 2, wherein in step D, the noise power cancellation is performed using the following formula:
<math><mrow><mover><mi>P</mi><mo>^</mo></mover><mo>=</mo><mover><mi>P</mi><mo>&OverBar;</mo></mover><mo>-</mo><mfrac><msup><mi>&sigma;</mi><mn>2</mn></msup><mi>KL</mi></mfrac></mrow></math>
wherein,
Figure FSB00000525382500024
representation elimination
Figure FSB00000525382500025
RSRP, σ after the noise power component contained in (1)2Is the noise power.
5. The estimation method according to claim 1, further comprising the steps of:
E. and smoothing the RSRP obtained by eliminating the noise power.
6. The estimation method according to claim 5, wherein in step E, the smoothing is performed using the following formula:
<math><mrow><msub><mi>P</mi><mi>m</mi></msub><mo>=</mo><mi>&alpha;</mi><msub><mover><mi>P</mi><mo>^</mo></mover><mi>m</mi></msub><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&alpha;</mi><mo>)</mo></mrow><msub><mi>P</mi><mrow><mi>m</mi><mo>-</mo><mn>1</mn></mrow></msub></mrow></math>
wherein, PmFor the smoothed RSRP current value,
Figure FSB00000525382500027
is the current value of RSRP without smoothing, Pm-1Alpha is a smoothing factor, and alpha is more than 0 and less than or equal to 1.
7. The estimation method according to claim 1, characterized in that:
the coherent block is a block composed of all resource elements REs within a coherence time range and a coherence bandwidth range.
8. The estimation method according to claim 1, characterized in that:
the coherent block is a resource block.
9. An apparatus for estimating a reference signal received power of an orthogonal frequency division multiplexing system, comprising:
an RS channel estimator for performing channel estimation on the received RS;
the RSRP calculator is used for calculating according to the channel estimation results of all RSs in the coherent block to obtain corresponding RSRP;
the RSRP averager is used for carrying out weighted average on the RSRPs corresponding to all coherent blocks in the measurement frequency bandwidth and the measurement time range;
the noise power eliminator is used for eliminating noise power of the RSRP obtained by the weighted average;
wherein the RSRP calculator calculates the RSRP using the following formula:
<math><mrow><msup><mi>P</mi><mrow><mo>&prime;</mo><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow></mrow></msup><mo>=</mo><msup><mrow><mo>|</mo><mfrac><mn>1</mn><mi>KL</mi></mfrac><munder><mi>&Sigma;</mi><mi>k</mi></munder><munder><mi>&Sigma;</mi><mi>l</mi></munder><msubsup><mover><mi>h</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>,</mo><mi>l</mi></mrow><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow></msubsup><mo>|</mo></mrow><mn>2</mn></msup></mrow></math>
wherein, P'(i,j,n)The RSRP corresponding to the nth coherent block transmitted by the ith antenna port and received by the jth antenna is represented, and K represents the RS occupied in the coherent blockThe number of subcarriers, L represents the number of OFDM symbols occupied by RS in the coherent block,
Figure FSB00000525382500032
represents the channel estimation of the kth RS signal in the nth coherent block, the kth RS signal in the frequency domain and the ith RS signal in the time domain, which are transmitted by the ith antenna port and received by the jth antenna, | · | represents the modulo operator.
10. The estimation apparatus of claim 9 wherein the RSRP averager performs a weighted average using the formula:
<math><mrow><mover><mi>P</mi><mo>&OverBar;</mo></mover><mo>=</mo><munder><mi>&Sigma;</mi><mi>i</mi></munder><msub><mi>&alpha;</mi><mi>i</mi></msub><munder><mi>&Sigma;</mi><mi>j</mi></munder><msub><mi>&beta;</mi><mi>j</mi></msub><mrow><mo>(</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munder><mi>&Sigma;</mi><mi>n</mi></munder><msup><mi>P</mi><mrow><mo>&prime;</mo><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow></mrow></msup><mo>)</mo></mrow><mo>,</mo><mn>1</mn><mo>&le;</mo><mi>i</mi><mo>&le;</mo><mi>I</mi><mo>,</mo><mn>1</mn><mo>&le;</mo><mi>j</mi><mo>&le;</mo><mi>J</mi></mrow></math>
wherein,
Figure FSB00000525382500034
represents to all P'(i,j,n)RSRP after weighted averaging, I denotes the total number of transmit antenna ports used to calculate RSRP, J denotes the total number of receive antennas used to calculate RSRP, N denotes the total number of coherent blocks within the measurement frequency bandwidth and measurement time range, αiAnd betajWeighting coefficients for the transmit antenna ports and the receive antennas respectively,
Figure FSB00000525382500035
Figure FSB00000525382500036
11. the estimation apparatus of claim 10, wherein:
<math><mrow><msub><mi>&alpha;</mi><mi>i</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>I</mi></mfrac><mo>,</mo></mrow></math> <math><mrow><msub><mi>&beta;</mi><mi>j</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>J</mi></mfrac><mo>.</mo></mrow></math>
12. the estimation apparatus as set forth in claim 10, wherein the noise power canceller performs noise power cancellation using the following equation:
<math><mrow><mover><mi>P</mi><mo>^</mo></mover><mo>=</mo><mover><mi>P</mi><mo>&OverBar;</mo></mover><mo>-</mo><mfrac><msup><mi>&sigma;</mi><mn>2</mn></msup><mi>KL</mi></mfrac></mrow></math>
wherein,representation eliminationRSRP, σ after the noise power component contained in (1)2Is the noise power.
13. The estimation apparatus of claim 9, further comprising:
and the RSRP smoother is used for smoothing the RSRP obtained by eliminating the noise power.
14. The estimation apparatus of claim 13 wherein the RSRP smoother smoothes using the formula:
<math><mrow><msub><mi>P</mi><mi>m</mi></msub><mo>=</mo><mi>&alpha;</mi><msub><mover><mi>P</mi><mo>^</mo></mover><mi>m</mi></msub><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&alpha;</mi><mo>)</mo></mrow><msub><mi>P</mi><mrow><mi>m</mi><mo>-</mo><mn>1</mn></mrow></msub></mrow></math>
wherein, PmFor the smoothed RSRP current value,
Figure FSB00000525382500041
is the current value of RSRP without smoothing, Pm-1Alpha is a smoothing factor, and alpha is more than 0 and less than or equal to 1.
15. The estimation apparatus of claim 9, wherein:
the coherent block is a block composed of all REs within a coherence time range and a coherence bandwidth range.
16. The estimation apparatus of claim 9, wherein:
the coherent block is a resource block.
CN2009100887531A 2009-07-10 2009-07-10 Method and device for estimating reference signal received power of orthogonal frequency division multiplexing system Expired - Fee Related CN101599939B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2009100887531A CN101599939B (en) 2009-07-10 2009-07-10 Method and device for estimating reference signal received power of orthogonal frequency division multiplexing system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009100887531A CN101599939B (en) 2009-07-10 2009-07-10 Method and device for estimating reference signal received power of orthogonal frequency division multiplexing system

Publications (2)

Publication Number Publication Date
CN101599939A CN101599939A (en) 2009-12-09
CN101599939B true CN101599939B (en) 2011-09-28

Family

ID=41421192

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009100887531A Expired - Fee Related CN101599939B (en) 2009-07-10 2009-07-10 Method and device for estimating reference signal received power of orthogonal frequency division multiplexing system

Country Status (1)

Country Link
CN (1) CN101599939B (en)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101688551B1 (en) 2010-02-11 2016-12-22 삼성전자주식회사 Method for indicating user specific dmrs antenna port in wireless communication systems
US9130719B2 (en) 2010-02-11 2015-09-08 Samsung Electronics Co., Ltd Method for indicating a DM-RS antenna port in a wireless communication system
CN102196486B (en) * 2010-12-31 2012-07-04 重庆重邮信科通信技术有限公司 Reference signal received power measurement method and device for orthogonal frequency-division multiplexing system
CN102098265B (en) * 2011-03-08 2013-11-06 华为技术有限公司 Method and device for determining reference signal received power (RSRP)
CN102111181B (en) 2011-03-18 2013-09-11 华为技术有限公司 Method and device for estimating energy per chip (Ec)/total transmitting power spectral density (Ior)
CN103929772B (en) * 2013-01-16 2017-06-16 展讯通信(上海)有限公司 Reference signal received power measurement method and device, the method for gravity treatment, user terminal
CN104283820A (en) * 2013-07-03 2015-01-14 普天信息技术研究院有限公司 Method for determining signal to noise ratio (SNR) of communication system
CN104469837A (en) * 2013-09-24 2015-03-25 富士通株式会社 Method and terminal for determining reference signal receiving power and system
US20160135066A1 (en) * 2014-03-07 2016-05-12 Telefonaktiebolaget L M Ericsson (Publ) Method and Device for Calculating Reference Signal Received Power
CN113595702A (en) * 2015-05-14 2021-11-02 英特尔公司 Enhanced radio resource management reporting in cellular systems
CN107787011B (en) * 2016-08-31 2021-02-12 普天信息技术有限公司 Simplified method for reference signal receiving power smoothing processing
CN109150335B (en) * 2018-09-30 2021-04-16 武汉虹信科技发展有限责任公司 Method and device for measuring uplink signal strength of target terminal in complex environment

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1823485A (en) * 2003-05-14 2006-08-23 高通股份有限公司 Interference and noise estimation in an OFDM system

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1823485A (en) * 2003-05-14 2006-08-23 高通股份有限公司 Interference and noise estimation in an OFDM system

Also Published As

Publication number Publication date
CN101599939A (en) 2009-12-09

Similar Documents

Publication Publication Date Title
CN101599939B (en) Method and device for estimating reference signal received power of orthogonal frequency division multiplexing system
JP4832261B2 (en) Channel estimation device
CN101534266B (en) Channel estimation method for Orthogonal Frequency Division Multiplexing system and device
US9020050B2 (en) Accounting for inter-carrier interference in determining a response of an OFDM communication channel
CN105391671B (en) Method and apparatus for channel estimation and equalization
CN101064571B (en) Apparatus for enhancing channel evaluation in OFDM receiver and its method
CN102196486B (en) Reference signal received power measurement method and device for orthogonal frequency-division multiplexing system
US20110159831A1 (en) Channel estimation for communication systems with multiple transmit antennas
CN102804719B (en) Method and apparatus for enhanced channel estimation in wireless communication systems
CN101577692A (en) Channel estimating method of orthogonal frequency division multiplexing system and device thereof
EP3125483B1 (en) Communication system determining time of arrival using matching pursuit
KR100945419B1 (en) Apparatus and method for estimating channel using sliding windows in broadband wireless communication system
CN104486266B (en) A kind of channel estimation methods and device based on MIMO-OFDM systems
Rodriguez-Pineiro et al. Experimental validation of ICI-aware OFDM receivers under time-varying conditions
US20100266078A1 (en) Radio communication device, and reception quality estimation method
CN101635598A (en) Method and device for estimating noise power
KR100835164B1 (en) Apparatus and method for channel estimate in broadband wireless access communication system
CN101599940B (en) Method and device for estimating noise power of orthogonal frequency division multiplexing system
CN102386953A (en) Transmitting and receiving calibration method for antenna in broadband orthogonal frequency division multiplexing (OFDM) system
US20080240262A1 (en) Ofdm/ofdma channel estimation
CN102487364B (en) Channel estimation method and apparatus thereof
KR101017865B1 (en) Channel estimation apparatus of OFDM receiver
CN104468426A (en) Method and system for estimating LTE uplink channel
CN102404269A (en) Method and device for origination of ranging process in TDD-OFDMA (time division duplex-orthogonal frequency division multiple access) system
CN102223336B (en) Wireless communication method and equipment

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20110928

Termination date: 20180710

CF01 Termination of patent right due to non-payment of annual fee