CN102457452B - Method for calculating frequency offset estimated value and channel estimated value and LTE base station - Google Patents

Abstract

The invention discloses a method for calculating a frequency offset estimated value and a channel estimated value and an LTE base station. Compared with the prior art, in the invention, through establishing function operation relation of a user between different channel estimation of each sampling period of a reference symbol upper frequency domain channel, or establishing function operation relation of the user between different theoretical signals and practical signals of each sampling period of the reference symbol upper frequency domain channel, when a target function value is minimal, a frequency offset estimated value and a channel estimated value of the user corresponding to a minimum value are taken as a frequency offset estimated value and a channel estimated value of the user. According to the invention, the LTE base station can estimate a frequency offset estimated value and a channel estimated value of each user at Normal CP format 2 in an uplink control channel in various channel scenes.

Description

A kind of method and LTE base station calculating frequency deviation estimated value and channel estimation value

Technical field

The present invention relates to 3GPP technical field, be specifically related to a kind of method that LTE base station calculates frequency deviation estimated value and channel estimation value.

Background technology

LTE (Long Term Evolution, Long Term Evolution) in ascending control channel be divided into form (format) 1/1a/1b and form 2/2a/2b, wherein form series 1 is for the transmission of dispatch request (SR) and feedback signaling (ACK/NACK), form series 2 is for channel quality instruction (CQI) and ACK/NACK transmission, and subscriber equipment (UE) selects corresponding form according to the category combinations of signal to be sent.Control channel is the frequency-domain and time-domain resource that multiple user's concurrent multiplexing is identical with an important difference of Uplink Shared Channel (PUSCH), causes existing method to be unable to estimate the frequency deviation of each user, and then affects the accuracy of channel estimating.

Summary of the invention

The main technical problem to be solved in the present invention is to provide a kind of LTE base station ascending control channel calculates frequency deviation estimation and channel estimating method when general cyclic prefix (Normal CP), estimates frequency deviation exactly and estimates and channel estimating.

For solving the problems of the technologies described above, the invention provides a kind of method calculating frequency deviation estimated value and channel estimation value, comprising step:

Predetermined frequency deviation sampling period q;

To user i at the first reference symbol l ₁the signal of upper reception frequency deviation f _iby described predetermined sampling period sampling, obtain the frequency deviation influence matrix M (f of user i _{i, q}), according to user i at the first reference symbol l ₁the signal of upper reception frequency deviation influence matrix M (f _{i, q}) and the first reference symbol l ₁the reference sequences that upper user i sends estimate, obtain user i at the first reference symbol l ₁on first channel estimation value in each sampling period of frequency domain channel

According to user i at the second reference symbol l ₂the signal of upper reception is estimated, obtain user i at the second reference symbol l ₂on the second channel estimated value in each sampling period of frequency domain channel estimate according to described second channel calculate user i at the first reference symbol l ₁3rd channel estimation value in each sampling period of upper frequency domain channel

First channel estimation value in each sampling period is calculated by preset function computing with the 3rd channel estimation value functional operation value between the two, and using frequency deviation corresponding for the minimum value in described functional operation value and channel estimation value as the frequency deviation estimated value of user and channel estimation value.

In a kind of embodiment, according to following formula:

${\hat{h}}_{l_1,i,q}=\frac{1}{12}{p}_{l_1,i}^{H}M(-f_{i,q})y_{l_1}$

Estimation user i is at the first reference symbol l ₁on first channel estimation value in each sampling period of frequency domain channel wherein for described reference sequences p _{l, i}complex-conjugate transpose, frequency deviation influence matrix M (-f _{i, q}) be described frequency deviation influence matrix M (f _{i, q}) inverse matrix.

In a kind of embodiment, described 3rd channel estimation value draws as follows:

In a kind of embodiment, described functional operation comprises: the first channel estimation value with the 3rd channel estimation value absolute value, first channel estimation value of both differences with the 3rd channel estimation value the absolute value of both difference square.

In a kind of embodiment, described minimum value adopts golden section search or binary search method or three search method search to draw.

A kind of LTE base station, comprises sampling unit and computing unit, wherein:

Sampling unit, for user i at the first reference symbol l ₁the signal of upper reception frequency deviation f _isample by the predetermined sampling period, obtain the frequency deviation influence matrix M (f of user i _{i, q}), wherein q is the sampling period;

Computing unit, is connected with sampling unit, for according to user i at the first reference symbol l ₁the signal of upper reception frequency deviation influence matrix M (f _{i, q}) and the first reference symbol l ₁the reference sequences that upper user i sends estimate, obtain user i at the first reference symbol l ₁on first channel estimation value in each sampling period of frequency domain channel according to user i at the second reference symbol l ₂the signal of upper reception is estimated, obtain user i at the second reference symbol l ₂on the second channel estimated value in each sampling period of frequency domain channel estimate according to described second channel calculate user i at the first reference symbol l ₁3rd channel estimation value in each sampling period of upper frequency domain channel first channel estimation value in each sampling period is calculated by preset function computing with the 3rd channel estimation value functional operation value between the two, and using frequency deviation corresponding for the minimum value in described functional operation value and channel estimation value as the frequency deviation estimated value of user and channel estimation value.

Calculate a method for frequency deviation estimation and channel estimating, comprise step:

Predetermined frequency deviation sampling period q;

According to user i at the second reference symbol l ₂upper actual Received signal strength carries out estimation and obtains user i at the second reference symbol l ₂on the second channel estimated value in each sampling period of frequency domain channel estimate according to described second channel calculate user i at the first reference symbol l ₁3rd channel estimation value in each sampling period of upper frequency domain channel

According to the first reference symbol l ₁the reference sequences p that upper user i sends _{l, i}, the 3rd channel estimation value with frequency deviation influence matrix M (f _{i, q}) construct user i at the first reference symbol l ₁the theoretical Received signal strength in each sampling period upper;

Theoretical Received signal strength and the actual Received signal strength functional operation value between the two in each sampling period is calculated by preset function computing, and using frequency deviation corresponding for the minimum value in described functional operation value and channel estimation value as the frequency deviation estimated value of user and channel estimation value.

In a kind of embodiment, described 3rd channel estimation value draws as follows:

In a kind of embodiment, according to following formula:

$y_{l_1}=\underset{i}{\mathrm{\Σ}}M\left(f_{i,q_i}\right)p_{l_1,i}{\tilde{h}}_{l_1,i,q_i}$

Construct the first reference symbol l ₁the theoretical Received signal strength in each sampling period upper, wherein for described reference sequences p _{l, i}complex-conjugate transpose, q _ifor the frequency deviation sampling period of user i.

In a kind of embodiment, described functional operation comprises: the absolute value of the difference of the absolute value of the difference of theoretical Received signal strength and actual Received signal strength, theoretical Received signal strength and actual Received signal strength square.

A kind of LTE base station, comprises sampling unit and computing unit, wherein:

Sampling unit, for user i at the first reference symbol l ₁the signal of upper reception frequency deviation f _isample by the predetermined sampling period, obtain the frequency deviation influence matrix M (f of user i _{i, q}), wherein q is the sampling period;

Computing unit, is connected with sampling unit, for according to user i at the second reference symbol l ₂upper actual Received signal strength carries out estimation and obtains user i at the second reference symbol l ₂on the second channel estimated value in each sampling period of frequency domain channel estimate according to described second channel calculate user i at the first reference symbol l ₁3rd channel estimation value in each sampling period of upper frequency domain channel according to the first reference symbol l ₁the reference sequences p that upper user i sends _{l, i}, the 3rd channel estimation value with frequency deviation influence matrix M (f _{i, q}) construct user i at the first reference symbol l ₁the theoretical Received signal strength in each sampling period upper; Theoretical Received signal strength and the actual Received signal strength functional operation value between the two in each sampling period is calculated by preset function computing, and using frequency deviation corresponding for the minimum value in described functional operation value and channel estimation value as the frequency deviation estimated value of user and channel estimation value.

Compared with prior art, functional operation relation between the present invention is estimated by the different channels in each sampling period of frequency domain channel of setting up user in reference symbol, or the functional operation relation set up between the different theories signal in each sampling period of user's frequency domain channel in reference symbol and actual signal, when target function value is minimum, the frequency deviation of the user that minimum value is corresponding and channel estimation value are as the frequency deviation estimated value of this user and channel estimation value, the present invention can make LTE base station estimate comparatively accurately under various channel scenario in ascending control channel at the frequency deviation estimated value of each user of general cyclic prefix (Normal CP) form 2 and channel estimation value.

Accompanying drawing explanation

A kind of method flow diagram calculating frequency deviation estimated value and channel estimation value that Fig. 1 provides for the embodiment of the present invention one;

A kind of LTE base station block diagram that Fig. 2 provides for the embodiment of the present invention two;

A kind of method flow diagram calculating frequency deviation estimated value and channel estimation value that Fig. 3 provides for the embodiment of the present invention three;

A kind of LTE base station block diagram that Fig. 4 provides for the embodiment of the present invention four.

Embodiment

According to LTE protocol (TS36.211) when general cyclic prefix (Normal CP), a subframe comprises two time slots, and each time slot comprises 7 symbols, and numbering 1 is 0,1,2,3,4,5,6,0,1,2,3,4,5,6.In the embodiment of the present invention in units of each time slot, therefore reference symbol numbering 1 is 0,1,2,3,4,5,6.

According to LTE protocol (TS36.211), in each time slot with l=1 and 1=5 for reference symbol, the reference sequences that reference symbol sends transmitting terminal and receiving terminal all known.The sequences y that receiving terminal receives in reference symbol _lbe the mixed signal of multi-user, be expressed as follows:

$y_l=\underset{i}{\mathrm{\Σ}}M\left(f_i\right)p_{l,i}{h}_{l,i}+n---\left(1\right)$

Wherein i is Customs Assigned Number, 1≤i≤12.F _iit is the frequency deviation of i-th user.H _{l, i}be the frequency domain channel of i-th user on l symbol, n is white noise.P _{l, i}be the reference sequences that on l symbol, i-th user sends, because ascending control channel only takies 1 Resource Block (RB), therefore p _{l, i}to be length be 12 column vector, and the reference sequences of two different users is orthogonal, namely as i ≠ j

M (f _i)=(a _{k, m}) be 12 rank Top profit thatch (Toeplitz) matrixes, its each element definition is as follows:

$a_{k,m}=\left\{\begin{array}{c}\left\{\frac{\sin \left({\mathrm{\π\γ}}_i\right)}{N\sin \left(\frac{{\mathrm{\π\γ}}_i}{N}\right)}\right\}\exp \left(j{\mathrm{\π\γ}}_i\frac{N-1}{N}\right),k=m\\ \left\{\frac{\sin \left({\mathrm{\π\γ}}_i\right)}{N\sin \left(\frac{\mathrm{\π}(m-k+{\mathrm{\γ}}_i)}{N}\right)}\right\}\exp (-\mathrm{j\π}\frac{m-k}{N})\exp \left(j{\mathrm{\π\γ}}_i\frac{N-1}{N}\right),k\≠m\end{array}\right.$

Wherein make γ _i=f _i/ 15000 (15000 is the sub-spectrum of carrier intervals of LTE base station one, and unit is hertz), N is that the IFFT (N ties up inverse fast Fourier transform) that OFDM (OFDM) modulates counts.K represents row k, and m represents m row.

By reference to the accompanying drawings the present invention is described in further detail below by embodiment.

Embodiment one

As shown in Figure 1, Fig. 1 is a kind of method flow diagram calculating frequency deviation estimated value and channel estimation value, comprises step:

Step S101: predetermined frequency deviation sampling period q.

Step S102: to user i at the first reference symbol l ₁the signal of upper reception frequency deviation f _iby described predetermined sampling period sampling, obtain the frequency deviation influence matrix M (f of user i _{i, q}), q=1,2 ..., Q _i(Q _ifor sampling sum).

In practice, f _iexcursion generally in [-1300Hz, 1300Hz], be that interval samples with 30Hz, then f _i=[-1300Hz ,-1270Hz ..., 1250Hz, 1280Hz], then there are 87 sampled points altogether.

Step S103: according to user i at the first reference symbol l ₁the signal of upper reception frequency deviation influence matrix M (f _{i, q}) and the first reference symbol l ₁the reference sequences that upper user i sends estimate, obtain user i at the first reference symbol l ₁on first channel estimation value in each sampling period of frequency domain channel

To (1) formula both sides premultiplication simultaneously

$p_{l,i}^{H}M(-f_{i,q})y_l=p_{l,i}^{H}M(f_1-f_{i,q})p_{l,1}{h}_{l,1}+p_{l,i}^{H}M(f_2-f_{i,q})p_{l,2}{h}_{l,2}---\left(2\right)$

$+p_{l,i}^{H}M(f_i-f_{i,q})p_{l,i}{h}_{l,i}+\·\·\·+p_{l,i}^{H}M(f_N-f_{i,q})p_{l,N}{h}_{l,N}+\tilde{n}$

If f _{i, q}equal actual f _i, then (2) formula becomes:

$p_{l,i}^{H}M(-f_{i,q})y_l$

$=12\·h_{l,i}+[p_{l,i}^{H}M(f_1-f_{i,q})p_{l,1}{h}_{l,1}+p_{l,i}^{H}M(f_2-f_{i,q})p_{l,2}{h}_{l,2}+\·\·\·+p_{l,i}^{H}M(f_N-f_{i,q})p_{l,N}{h}_{l,N}+\tilde{n}]$

$=12\·h_{l,i}+\mathrm{IN}$

Wherein IN is interference noise item, for p _{l, i}complex-conjugate transpose, frequency deviation influence matrix M (-f _{i, q}) be described frequency deviation influence matrix M (f _{i, q}) inverse matrix.

In the present embodiment, definition for user i samples f in frequency deviation _{i, q}under channel estimating on symbol l (=1 and 5).

Step S104: according to user i at the second reference symbol l ₂the signal of upper reception is estimated, obtain user i at the second reference symbol l ₂on the second channel estimated value in each sampling period of frequency domain channel estimate according to described second channel calculate user i at the first reference symbol l ₁3rd channel estimation value in each sampling period of upper frequency domain channel

The duration of a LTE system time slot is 0.5 millisecond (ms), when general cyclic prefix (NormalCP) time slot be divided into 7 symbols (l=0,1 ..., 6), the time difference therefore between l=1 and l=5 is in the so short time, channel variation is very little, but the phase rotating that the larger frequency deviation of existence is brought, at f _{i, q}lower twiddle factor is in LTE system, suppose the first reference symbol l ₁=1, second reference symbol l ₂=5, then the channel estimating of reference symbol 5 is obtained by reference to symbol 1 conversely, if the channel estimating obtaining reference symbol 1 by reference to symbol 5 is then

Step S105: the first channel estimation value calculating each sampling period by preset function computing with the 3rd channel estimation value functional operation value between the two, and using frequency deviation corresponding for the minimum value in described functional operation value and channel estimation value as the frequency deviation estimated value of user and channel estimation value.

In one embodiment, described functional operation comprises: the first channel estimation value with the 3rd channel estimation value the absolute value of both differences namely first channel estimation value with the 3rd channel estimation value the absolute value of both difference square certainly these two kinds are also not limited to, as long as other modes thought according to the invention is also within protection scope of the present invention.

Δ in frequency deviation region-1300 hertz ~+1300 hertz _{i, q}for single peak function, and peak value is minimum.Therefore minimum delta is searched for _{i, q}method can adopt golden section search (Golden Section Search), binary search method (Quadratic Search), three search methods (Cubic Search) and at equal intervals frequency deviation adopt trial-and-error method one by one to find Δ from 87 sampled points _{i, q}minimum value, then the f of its correspondence _{i, q}with for the frequency deviation of user i is estimated and channel estimating.Work as f _{i, q}equal actual f _itime, Δ _{i, q}obtain minimum value.Therefore minimum Δ _{i, q}corresponding f _{i, q}can estimate as the frequency deviation of i-th user formula is expressed as: simultaneously corresponding (l=1 and l=5) is the channel estimation value in corresponding reference symbol.

Embodiment two

According to embodiment one, embodiment two provides a kind of LTE base station, and as shown in Figure 2, LTE base station comprises sampling unit 201 and computing unit 202, wherein:

Sampling unit 201, for user i at the first reference symbol l ₁the signal of upper reception frequency deviation f _isample by the predetermined sampling period, obtain the frequency deviation influence matrix M (f of user i _{i, q}), wherein q is the sampling period.

Computing unit 202, is connected with sampling unit 201, for according to user i at the first reference symbol l ₁the signal of upper reception frequency deviation influence matrix M (f _{i, q}) and the first reference symbol l ₁the reference sequences that upper user i sends estimate, obtain user i at the first reference symbol l ₁on first channel estimation value in each sampling period of frequency domain channel according to user i at the second reference symbol l ₂the signal of upper reception is estimated, obtain user i at the second reference symbol l ₂on the second channel estimated value in each sampling period of frequency domain channel estimate according to described second channel calculate user i at the first reference symbol l ₁3rd channel estimation value in each sampling period of upper frequency domain channel first channel estimation value in each sampling period is calculated by preset function computing with the 3rd channel estimation value functional operation value between the two, and using frequency deviation corresponding for the minimum value in described functional operation value and channel estimation value as the frequency deviation estimated value of user and channel estimation value.

Embodiment three

Refer to shown in Fig. 3, a kind of method calculating frequency deviation estimated value and channel estimation value, comprises step:

Step S301: predetermined frequency deviation sampling period q.

Step S302: in frequency deviation region to user i at the first reference symbol l ₁the signal of upper reception frequency deviation f _icarry out predetermined period sampling, obtain the frequency deviation influence matrix M (f of user i _{i, q}), wherein q is the sampling period.

Step S303: according to user i at the second reference symbol l ₂upper actual Received signal strength carries out estimation and obtains user i at the second reference symbol l ₂on the second channel estimated value in each sampling period of frequency domain channel estimate according to described second channel calculate user i at the first reference symbol l ₁3rd channel estimation value in each sampling period of upper frequency domain channel

Step S304: according to the first reference symbol l ₁the reference sequences p that upper user i sends _{l, i}, the 3rd channel estimation value with frequency deviation influence matrix M (f _{i, q}) construct user i at the first reference symbol l ₁the theoretical Received signal strength in each sampling period upper.

In the present embodiment, according to following formula:

$y_{l_1}=\underset{i}{\mathrm{\Σ}}M\left(f_{i,q_i}\right)p_{l_1,i}{\tilde{h}}_{l_1,i,q_i}$

Reconstructing user i is at the first reference symbol l ₁on theoretical Received signal strength, wherein for p _{l, i}complex-conjugate transpose, q _ifor the frequency deviation sampling period of user i.

Step S305: the theoretical Received signal strength and the actual Received signal strength functional operation value between the two that calculate each sampling period by preset function computing, and using frequency deviation corresponding for the minimum value in described functional operation value and channel estimation value as the frequency deviation estimated value of user and channel estimation value.

In one embodiment, described target function can be the absolute value of the difference of theoretical Received signal strength and actual Received signal strength, also can be the absolute value of the difference of theoretical Received signal strength and actual Received signal strength square.Certainly these two kinds are also not limited to, as long as other modes thought according to the invention is also within protection scope of the present invention.

Embodiment four

According to embodiment three, embodiment four provides a kind of LTE base station, and as shown in Figure 4, LTE base station comprises sampling unit 401 and computing unit 402, wherein:

Sampling unit 401, for user i at the first reference symbol l ₁the signal of upper reception frequency deviation f _isample by the predetermined sampling period, obtain the frequency deviation influence matrix M (f of user i _{i, q}), wherein q is the sampling period.

Computing unit 402, is connected with sampling unit 401, for according to user i at the second reference symbol l ₂upper actual Received signal strength carries out estimation and obtains user i at the second reference symbol l ₂on the second channel estimated value in each sampling period of frequency domain channel estimate according to described second channel calculate user i at the first reference symbol l ₁3rd channel estimation value in each sampling period of upper frequency domain channel according to the first reference symbol l ₁the reference sequences p that upper user i sends _{l, i}, the 3rd channel estimation value with frequency deviation influence matrix M (f _{i, q}) construct user i at the first reference symbol l ₁the theoretical Received signal strength in each sampling period upper; Theoretical Received signal strength and the actual Received signal strength functional operation value between the two in each sampling period is calculated by preset function computing, and using frequency deviation corresponding for the minimum value in described functional operation value and channel estimation value as the frequency deviation estimated value of user and channel estimation value.

In sum, the frequency deviation that the present invention can make LTE base station estimate each user of form 2 in ascending control channel comparatively accurately under various channel scenario is estimated and channel estimating.

Above content is in conjunction with concrete execution mode further description made for the present invention, can not assert that specific embodiment of the invention is confined to these explanations.For general technical staff of the technical field of the invention, without departing from the inventive concept of the premise, some simple deduction or replace can also be made, all should be considered as belonging to protection scope of the present invention.

Claims (7)

1. calculate a method for frequency deviation estimated value and channel estimation value, it is characterized in that, comprise step:

Predetermined frequency deviation sampling period q;

To user i in the first reference symbol l ₁the signal y of upper reception _l1frequency deviation f _iby described predetermined sampling period sampling, obtain the frequency deviation influence matrix M (f of user i _{i, q}), f _i,qit is the frequency deviation value in q sampling period of i-th user; According to user i in the first reference symbol l ₁the signal y of upper reception _l1, frequency deviation influence matrix M (f _{i, q}) and the first reference symbol l ₁the reference sequences p that upper user i sends _{l1, i}estimate, obtain user i in the first reference symbol l ₁on first channel estimation value in each sampling period of frequency domain channel

According to user i in the second reference symbol l ₂the signal of upper reception is estimated, obtain user i in the second reference symbol l ₂on the second channel estimated value in each sampling period of frequency domain channel estimate according to described second channel calculate user i in the first reference symbol l ₁3rd channel estimation value in each sampling period of upper frequency domain channel described 3rd channel estimation value draws as follows:

First channel estimation value in each sampling period is calculated by preset function computing with the 3rd channel estimation value functional operation value between the two, and using frequency deviation corresponding for the minimum value in described functional operation value and channel estimation value as the frequency deviation estimated value of user and channel estimation value; Described functional operation comprises: the first channel estimation value with the 3rd channel estimation value the absolute value of both differences, or the first channel estimation value with the 3rd channel estimation value the absolute value of both difference square.

2. the method for claim 1, is characterized in that, according to following formula:

Estimation user i is in the first reference symbol l ₁on first channel estimation value in each sampling period of frequency domain channel wherein for described reference sequences p _{l, i}complex-conjugate transpose, frequency deviation influence matrix M (-f _{i, q}) be described frequency deviation influence matrix M (f _{i, q}) inverse matrix.

3. method as claimed in claim 1 or 2, is characterized in that, described minimum value adopts golden section search or binary search method or three search method search to draw.

4. a LTE base station, is characterized in that, comprises sampling unit and computing unit, wherein:

Sampling unit, for user i in the first reference symbol l ₁the signal y of upper reception _l1frequency deviation f _isample by the predetermined sampling period, obtain the frequency deviation influence matrix M (f of user i _{i, q}), f _i,qit is the frequency deviation value in q sampling period of i-th user; Wherein q is the sampling period;

Computing unit, is connected with sampling unit, for according to user i in the first reference symbol l ₁the signal y of upper reception _l1, frequency deviation influence matrix M (f _{i, q}) and the first reference symbol l ₁the reference sequences p that upper user i sends _{l1, i}estimate, obtain user i in the first reference symbol l ₁on first channel estimation value in each sampling period of frequency domain channel according to user i in the second reference symbol l ₂the signal of upper reception is estimated, obtain user i in the second reference symbol l ₂on the second channel estimated value in each sampling period of frequency domain channel estimate according to described second channel calculate user i in the first reference symbol l ₁3rd channel estimation value in each sampling period of upper frequency domain channel described 3rd channel estimation value draws as follows: first channel estimation value in each sampling period is calculated by preset function computing with the 3rd channel estimation value functional operation value between the two, and using frequency deviation corresponding for the minimum value in described functional operation value and channel estimation value as the frequency deviation estimated value of user and channel estimation value; Described functional operation comprises: the first channel estimation value with the 3rd channel estimation value the absolute value of both differences, or the first channel estimation value with the 3rd channel estimation value the absolute value of both difference square.

5. calculate a method for frequency deviation estimated value and channel estimation value, it is characterized in that, comprise step:

Predetermined frequency deviation sampling period q;

According to user i in the second reference symbol l ₂upper actual Received signal strength carry out estimating obtain user i in the second reference symbol l ₂on the second channel estimated value in each sampling period of frequency domain channel estimate according to described second channel calculate user i in the first reference symbol l ₁3rd channel estimation value in each sampling period of upper frequency domain channel described 3rd channel estimation value draws as follows:

According to user first reference symbol l ₁the reference sequences p that upper user i sends _{l, i}, the 3rd channel estimation with frequency deviation influence matrix M (f _{i, q}) construct user i in the first reference symbol l ₁the theoretical Received signal strength in each sampling period upper; f _i,qit is the frequency deviation value in q sampling period of i-th user;

Theoretical Received signal strength and the actual Received signal strength functional operation value between the two in each sampling period is calculated by preset function computing, and using frequency deviation corresponding for the minimum value in described functional operation value and channel estimation value as the frequency deviation estimated value of user and channel estimation value; Described functional operation comprises: the first channel estimation value with the 3rd channel estimation value the absolute value of both differences, or the first channel estimation value with the 3rd channel estimation value the absolute value of both difference square.

6. method as claimed in claim 5, is characterized in that, according to following formula:

Construct the first reference symbol l ₁the theoretical Received signal strength in each sampling period upper, wherein for described reference sequences p _{l, i}complex-conjugate transpose, q _ifor the frequency deviation sampling period of user i.

7. a LTE base station, is characterized in that, comprises sampling unit and computing unit, wherein:

Sampling unit, for user i in the first reference symbol l ₁the signal y of upper reception _l1frequency deviation f _isample by the predetermined sampling period, obtain the frequency deviation influence matrix M (f of user i _{i, q}), f _i,qit is the frequency deviation value in q sampling period of i-th user; Wherein q is the sampling period;

Computing unit, is connected with sampling unit, for according to user i in the second reference symbol l ₂upper actual Received signal strength carry out estimating obtain user i in the second reference symbol l ₂on the second channel estimated value in each sampling period of frequency domain channel estimate according to described second channel calculate user i in the first reference symbol l ₁3rd channel estimation value in each sampling period of upper frequency domain channel described 3rd channel estimation value draws as follows: according to user first reference symbol l ₁the reference sequences p that upper user i sends _{l, i}, the 3rd channel estimation with frequency deviation influence matrix M (f _{i, q}) construct user i in the first reference symbol l ₁the theoretical Received signal strength in each sampling period upper; Theoretical Received signal strength and the actual Received signal strength functional operation value between the two in each sampling period is calculated by preset function computing, and using frequency deviation corresponding for the minimum value in described functional operation value and channel estimation value as the frequency deviation estimated value of user and channel estimation value; Described functional operation comprises: the first channel estimation value with the 3rd channel estimation value the absolute value of both differences, or the first channel estimation value with the 3rd channel estimation value the absolute value of both difference square.