Method, device and storage medium for frequency offset estimate

Abstract

The present invention provides a method, device and storage medium for frequency offset estimate, and the method for frequency offset estimate comprises: performing initial channel estimate by receiving an Sounding Reference Signal (SRS) or a Preamble code transmitted by a User Equipment (UE) to obtain a first frequency offset estimate value of a current sub-frame of the UE (S101); determining a frequency offset range in which the first frequency offset estimate value is located, and obtaining an initial frequency offset value corresponding to the first frequency offset estimate value according to the preset correspondence between the range of the frequency offset and the initial frequency offset value (S102); determining a phase difference of adjacent pilot positions based on a received pilot sequence included in the current sub-frame, and obtaining a dual pilot frequency offset estimate value according to the phase difference within the range limited by the initial frequency offset value corresponding to the first frequency offset estimate value (S103); and acquiring a frequency offset estimate value of the current sub-frame according to the initial frequency offset value corresponding to the first frequency offset estimate value and the dual pilot frequency offset estimate value (S104).

Claims

1. A method for frequency offset estimation comprising: performing initial channel estimation with a received Sounding Reference Signal (SRS) or Preamble code to obtain a first frequency offset estimation value of a current sub-frame; determining a frequency offset range in which the first frequency offset estimation value is located, and obtaining an initial frequency offset value corresponding to the first frequency offset estimation value according to a preset correspondence between the frequency offset range and the initial frequency offset value; determining a phase difference of adjacent pilot locations based on a received pilot sequence included in the current sub-frame, and obtaining, within the range defined by the initial frequency offset value corresponding to the first frequency offset estimation value, a dual pilot frequency offset estimation value according to the phase difference; and obtaining a frequency offset estimation value of the current sub-frame according to the initial frequency offset value corresponding to the first frequency offset estimation value and the dual pilot frequency offset estimation value.

2. The method according to claim 1, wherein the performing initial channel estimation with the received SRS or Preamble code to obtain the first frequency offset estimation value of the current sub-frame comprises: obtaining a channel estimation value H.sup.u and received data Y.sup.u by receiving the SRS or Preamble code of the current sub-frame; obtaining a time-domain sample value w.sup.u(n,ka.sub.Rx) of the current sub-frame on a receiving antenna ka.sub.Rx according to the channel estimation value H.sup.u and the received data Y.sup.u and in combination with a local pilot sequence P.sup.u and the following formula:
w.sup.u(n,ka.sub.Rx)=IDFT(Y.sup.u(k,ka.sub.Rx)H.sup.u*(k,ka.sub.Rx))(IDFT(P.sup.u(k)))*, where n represents a series number of a sample point, u represents a User Equipment (UE) corresponding to the current sub-frame, k represents a sub-carrier index corresponding to the UE, 0k12N.sub.RB.sup.u1, ka.sub.Rx represents an index of the receiving antenna, N.sub.RB.sup.u represents a size of a Resource Block (RB) occupied by the UE, * represents a conjugate operation, and IDFT represents an inverse Fourier transform; and obtaining the first frequency offset estimation value f.sub.0.sup.u(ka.sub.Rx) of the current sub-frame on the receiving antenna ka.sub.Rx according to phase information of the time-domain sample value w.sup.u(n,ka.sub.Rx) of the current sub-frame on the receiving antenna ka.sub.Rx.

3. The method according to claim 2, wherein the obtaining the first frequency offset estimation value f.sub.0.sup.u(ka.sub.Rx) of the current sub-frame on the receiving antenna ka.sub.Rx according to the phase information of the time-domain sample value w.sup.u(n,ka.sub.Rx) of the current sub-frame on the receiving antenna ka.sub.Rx comprises: averaging respectively phase values of a first half and a second half of the time-domain sample value w.sup.u(n,ka.sub.Rx) of the current sub-frame to obtain a phase difference .sub.1.sup.u(ka.sub.Rx) of the first half and the second half; and obtaining the first frequency offset estimation value f.sub.0.sup.u(ka.sub.Rx) of the current sub-frame on the receiving antenna ka.sub.Rx according to the following formula: $f_0^u\left({\mathrm{ka}}_{\mathrm{Rx}}\right)=\frac{{}_1^u\left({\mathrm{ka}}_{\mathrm{Rx}}\right)2f}{}.$

4. The method according to claim 2, wherein the method further comprises: performing a smoothing processing on the first frequency offset estimation value f.sub.0.sup.u(ka.sub.Rx) of the current sub-frame on the receiving antenna ka.sub.Rx to obtain a smoothed first frequency offset estimation value f.sub.1.sup.(u)(ka.sub.Rx) of the current sub-frame.

5. The method according to claim 1, wherein the step of determining the phase difference of adjacent pilot locations based on the received pilot sequence included in the current sub-frame and obtaining, within the range defined by the initial frequency offset value corresponding to the first frequency offset estimation value, the dual pilot frequency offset estimation value according to the phase difference comprises: obtaining a phase difference .sub.2.sup.u(ka.sub.Rx) of two pilot locations on the receiving antenna ka.sub.Rx based on two pilot channel estimation values H.sub.1.sup.u(k,ka.sub.Rx) and H.sub.2.sup.u(k,ka.sub.Rx) obtained from two pilot sequences included in the current sub-frame on the receiving antenna ka.sub.Rx and according to the following formula: ${}_2^u\left({\mathrm{ka}}_{\mathrm{Rx}}\right)=a\mathrm{ngle}\left(\underset{k=0}{\overset{12N_{\mathrm{RB}}^2-1}{.Math.}}{H}_1^u\left(k,{\mathrm{ka}}_{\mathrm{Rx}}\right)H_2^{u^{*}}\left(k,{\mathrm{ka}}_{\mathrm{Rx}}\right)\right),$ where * represents a conjugate operator, and angle( ) represents an operator for calculating an angle; obtaining a dual pilot frequency offset estimation result f.sub.2.sup.(u)(ka.sub.Rx) on the receiving antenna ka.sub.Rx based on a time interval t and the phase difference .sub.2.sup.u(ka.sub.Rx) of the two pilot sequences and according to the following formula: $f_2^{\left(u\right)}\left({\mathrm{ka}}_{\mathrm{Rx}}\right)=\frac{{}_2^u\left({\mathrm{ka}}_{\mathrm{Rx}}\right)}{2t};$ and obtaining, within the range defined by the initial frequency offset value f.sub.initial.sup.(u)(ka.sub.Rx), the dual pilot frequency offset estimation value f.sub.2(ka.sub.Rx) on the receiving antenna ka.sub.Rx based on the dual pilot frequency offset estimation result f.sub.2.sup.(u)(ka.sub.Rx) and according to the following formula:
f.sub.2(ka.sub.Rx)=f.sub.2.sup.(u)(ka.sub.Rx)f.sub.initial.sup.(u)(ka.sub.Rx)+lT, where T represents the maximum range values for the dual pilot frequency offset estimation value, and l represents a preset integer for defining the dual pilot frequency offset estimation value within a predetermined range.

6. The method according to claim 1, wherein after the obtaining the frequency offset estimation value of the current sub-frame according to the initial frequency offset value corresponding to the first frequency offset estimation value and the dual pilot frequency offset estimation value, the method further comprises: performing a smoothing processing on the frequency offset estimation value f.sup.u(ka.sub.Rx) of the current sub-frame by means of filtering to obtain a smoothed frequency offset estimation value f.sup.(u)(ka.sub.Rx) of the current sub-frame.

7. A device for frequency offset estimation comprising: a processor; and a memory for storing instructions executable by the processor; wherein the processor is configured to: perform initial channel estimation with a received Sounding Reference Signal (SRS) or Preamble code to obtain a first frequency offset estimation value of a current sub-frame; determine a frequency offset range in which the first frequency offset estimation value is located and obtain an initial frequency offset value corresponding to the first frequency offset estimation value according to a preset correspondence between the frequency offset range and the initial frequency offset value; determine a phase difference of adjacent pilot locations based on a received pilot sequence included in the current sub-frame and obtain, within the range defined by the initial frequency offset value corresponding to the first frequency offset estimation value, a duel pilot frequency offset estimation value according to the phase difference; and obtain a frequency offset estimation value of the current sub-frame according to the initial frequency offset value corresponding to the first frequency offset estimation value and the dual pilot frequency offset estimation value.

8. The device according to claim 7, wherein the processor is further configured to: obtain a channel estimation value H.sup.u and received data Y.sup.u by receiving the SRS or Preamble code of the current sub-frame; obtain a time-domain sample value w.sup.u(n,ka.sub.Rx) of the current sub-frame on a receiving antenna ka.sub.Rx according to the channel estimation value H.sup.u and the received data Y.sup.u and in combination with a local pilot sequence P.sup.u and the following formula:
w.sup.u(n,ka.sub.Rx)=IDFT(Y.sup.u(k,ka.sub.Rx)H.sup.u*(k,ka.sub.Rx))(IDFT(P.sup.u(k)))*, where n represents a series number of a sample point, u represents a User Equipment (UE) corresponding to the current sub-frame, k represents a sub-carrier index corresponding to the UE, 0k12N.sub.RB.sup.u1, ka.sub.Rx represents an index of the receiving antenna, N.sub.RB.sup.u represents a size of a Resource Block (RB) occupied by the UE, * represents a conjugate operation, and IDFT represents an inverse Fourier transform; and obtain the first frequency offset estimation value f.sub.0.sup.u(ka.sub.Rx) of the current sub-frame on the receiving antenna ka.sub.Rx according to phase information of the time-domain sample value w.sup.u(n,ka.sub.Rx) of the current sub-frame on the receiving antenna ka.sub.Rx.

9. The device according to claim 8, wherein the processor is further configured to: average respectively phase values of a first half and a second half of the time-domain sample value w.sup.u(n,ka.sub.Rx) of the current sub-frame to obtain a phase difference .sub.1.sup.u(ka.sub.Rx) of the first half and the second half; and obtain the first frequency offset estimation value f.sub.0.sup.u(ka.sub.Rx) of the current sub-frame on the receiving antenna ka.sub.Rx according to the following formula: $f_0^u\left({\mathrm{ka}}_{\mathrm{Rx}}\right)=\frac{{}_1^u\left({\mathrm{ka}}_{\mathrm{Rx}}\right)2f}{}.$

10. The device according to claim 8, wherein the processor is further configured to perform a smoothing processing on the first frequency offset estimation value f.sub.0.sup.u(ka.sub.Rx) of the current sub-frame on the receiving antenna ka.sub.Rx to obtain a smoothed first frequency offset estimation value f.sub.1.sup.(u)(ka.sub.Rx) of the current sub-frame.

11. The device according to claim 7, wherein the processor is further configured to: obtain a phase difference .sub.2.sup.u(ka.sub.Rx) of two pilot locations on the receiving antenna ka.sub.Rx based on two pilot channel estimation values H.sub.1.sup.u(k,ka.sub.Rx) and H.sub.2.sup.u(k,ka.sub.Rx) obtained from two pilot sequences included in the current sub-frame on the receiving antenna ka.sub.Rx and according to the following formula: ${}_2^u\left({\mathrm{ka}}_{\mathrm{Rx}}\right)=a\mathrm{ngle}\left(\underset{k=0}{\overset{12N_{\mathrm{RB}}^2-1}{.Math.}}{H}_1^u\left(k,{\mathrm{ka}}_{\mathrm{Rx}}\right)H_2^{u^{*}}\left(k,{\mathrm{ka}}_{\mathrm{Rx}}\right)\right),$ where * represents a conjugate operator, and angle( ) represents an operator for calculating an angle; obtain a dual pilot frequency offset estimation result f.sub.2.sup.(u)(ka.sub.Rx) on the receiving antenna ka.sub.Rx based on a time interval t and a phase difference .sub.2.sup.u(ka.sub.Rx) of the two pilot sequences and according to the following formula: $f_2^{\left(u\right)}\left({\mathrm{ka}}_{\mathrm{Rx}}\right)=\frac{{}_2^u\left({\mathrm{ka}}_{\mathrm{Rx}}\right)}{2t};$ and obtain, within the range defined by the initial frequency offset value f.sub.initial.sup.(u)(ka.sub.Rx), a dual pilot frequency offset estimation value f.sub.2(ka.sub.Rx) on the receiving antenna ka.sub.Rx based on the dual pilot frequency offset estimation result f.sub.2.sup.(u)(ka.sub.Rx) and according to the following formula:
f.sub.2(ka.sub.Rx)=f.sub.2.sup.(u)(ka.sub.Rx)f.sub.initial.sup.(u)(ka.sub.Rx)+lT, where T represents the maximum range values for the dual pilot frequency offset estimation value, and l represents a preset integer for defining the dual pilot frequency offset estimation value within a predetermined range.

12. The device according claim 7, wherein the processor is further configured to perform a smoothing processing on the frequency offset estimation value f.sup.u(ka.sub.Rx) of the current sub-frame by means of filtering to obtain a smoothed frequency offset estimation value f.sup.(u)(ka.sub.Rx) of the current sub-frame.

13. A non-transitory computer readable storage medium in which a computer program is stored, wherein the computer program comprising executable instructions that cause a computer to: perform initial channel estimation with a received Sounding Reference Signal (SRS) or Preamble code to obtain a first frequency offset estimation value of a current sub-frame; determine a frequency offset range in which the first frequency offset estimation value is located, and obtain an initial frequency offset value corresponding to the first frequency offset estimation value according to a preset correspondence between the frequency offset range and the initial frequency offset value; determine a phase difference of adjacent pilot locations based on a received pilot sequence included in the current sub-frame, and obtain, within the range defined by the initial frequency offset value corresponding to the first frequency offset estimation value, a dual pilot frequency offset estimation value according to the phase difference; and obtain a frequency offset estimation value of the current sub-frame according to the initial frequency offset value corresponding to the first frequency offset estimation value and the dual pilot frequency offset estimation value.

14. The non-transitory computer readable storage medium according to claim 13, wherein performing initial channel estimation with the received SRS or Preamble code to obtain the first frequency offset estimation value of the current sub-frame comprises: obtaining a channel estimation value H.sup.u and received data Y.sup.u by receiving the SRS or Preamble code of the current sub-frame; obtaining a time-domain sample value w.sup.u(n,ka.sub.Rx) of the current sub-frame on a receiving antenna ka.sub.Rx according to the channel estimation value H.sup.u and the received data Y.sup.u and in combination with a local pilot sequence P.sup.u and the following formula:
w.sup.u(n,ka.sub.Rx)=IDFT(Y.sup.u(k,ka.sub.Rx)H.sup.u*(k,ka.sub.Rx))(IDFT(P.sup.u(k)))*, where n represents a series number of a sample point, u represents a User Equipment (UE) corresponding to the current sub-frame, k represents a sub-carrier index corresponding to the UE, 0k12N.sub.RB.sup.u1, ka.sub.Rx represents an index of the receiving antenna, N.sub.RB.sup.u represents a size of a Resource Block (RB) occupied by the UE, * represents a conjugate operation, and IDFT represents an inverse Fourier transform; and obtaining the first frequency offset estimation value f.sub.0.sup.u(ka.sub.Rx) of the current sub-frame on the receiving antenna ka.sub.Rx according to phase information of the time-domain sample value w.sup.u(n,ka.sub.Rx) of the current sub-frame on the receiving antenna ka.sub.Rx.

15. The non-transitory computer readable storage medium according to claim 14, wherein obtaining the first frequency offset estimation value f.sub.0.sup.u(ka.sub.Rx) of the current sub-frame on the receiving antenna ka.sub.Rx according to the phase information of the time-domain sample value w.sup.u(n,ka.sub.Rx) of the current sub-frame on the receiving antenna ka.sub.Rx comprises: averaging respectively phase values of a first half and a second half of the time-domain sample value w.sup.u(n,ka.sub.Rx) of the current sub-frame to obtain a phase difference .sub.1.sup.u(ka.sub.Rx) of the first half and the second half; and obtaining the first frequency offset estimation value f.sub.0.sup.u(ka.sub.Rx) of the current sub-frame on the receiving antenna ka.sub.Rx according to the following formula: $f_0^u\left({\mathrm{ka}}_{\mathrm{Rx}}\right)=\frac{{}_1^u\left({\mathrm{ka}}_{\mathrm{Rx}}\right)2f}{}.$

16. The non-transitory computer readable storage medium according to claim 14, wherein the computer program comprising executable instructions that further cause a computer to: perform a smoothing processing on the first frequency offset estimation value f.sub.0.sup.u(ka.sub.Rx) of the current sub-frame on the receiving antenna ka.sub.Rx to obtain a smoothed first frequency offset estimation value f.sub.1.sup.(u)(ka.sub.Rx) of the current sub-frame.

17. The non-transitory computer readable storage medium according to claim 13, wherein the step of determining the phase difference of adjacent pilot locations based on the received pilot sequence included in the current sub-frame and obtaining, within the range defined by the initial frequency offset value corresponding to the first frequency offset estimation value, the dual pilot frequency offset estimation value according to the phase difference comprises: obtaining a phase difference .sub.2.sup.u(ka.sub.Rx) of two pilot locations on the receiving antenna ka.sub.Rx based on two pilot channel estimation values H.sub.1.sup.u(k,ka.sub.Rx) and H.sub.2.sup.u(k,ka.sub.Rx) obtained from two pilot sequences included in the current sub-frame on the receiving antenna ka.sub.Rx and according to the following formula: ${}_2^u\left({\mathrm{ka}}_{\mathrm{Rx}}\right)=a\mathrm{ngle}\left(\underset{k=0}{\overset{12N_{\mathrm{RB}}^2-1}{.Math.}}{H}_1^u\left(k,{\mathrm{ka}}_{\mathrm{Rx}}\right)H_2^{u^{*}}\left(k,{\mathrm{ka}}_{\mathrm{Rx}}\right)\right),$ where * represents a conjugate operator, and angle( ) represents an operator for calculating an angle; obtaining a dual pilot frequency offset estimation result f.sub.2.sup.(u)(ka.sub.Rx) on the receiving antenna ka.sub.Rx based on a time interval t and the phase difference .sub.2.sup.u(ka.sub.Rx) of the two pilot sequences and according to the following formula: $f_2^{\left(u\right)}\left({\mathrm{ka}}_{\mathrm{Rx}}\right)=\frac{{}_2^u\left({\mathrm{ka}}_{\mathrm{Rx}}\right)}{2t};$ and obtaining, within the range defined by the initial frequency offset value f.sub.initial.sup.(u)(ka.sub.Rx), the dual pilot frequency offset estimation value f.sub.2(ka.sub.Rx) on the receiving antenna ka.sub.Rx based on the dual pilot frequency offset estimation result f.sub.2.sup.(u)(ka.sub.Rx) and according to the following formula:
f.sub.2(ka.sub.Rx)=f.sub.2.sup.(u)(ka.sub.Rx)f.sub.initial.sup.(u)(ka.sub.Rx)+lT, where T represents the maximum range values for the dual pilot frequency offset estimation value, and l represents a preset integer for defining the dual pilot frequency offset estimation value within a predetermined range.

18. The non-transitory computer readable storage medium according to claim 13, wherein the computer program comprising executable instructions that further cause a computer to perform a smoothing processing on the frequency offset estimation value f.sup.u(ka.sub.Rx) of the current sub-frame by means of filtering to obtain a smoothed frequency offset estimation value f.sup.(u)(ka.sub.Rx) of the current sub-frame, after obtaining the frequency offset estimation value of the current sub-frame according to the initial frequency offset value corresponding to the first frequency offset estimation value and the dual pilot frequency offset estimation value.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIG. 1 is a flowchart of a method for frequency offset estimation according to an embodiment of the disclosure; and

(2) FIG. 2 is a structural diagram of a device for frequency offset estimation according to an embodiment of the disclosure.

DETAILED DESCRIPTION

(3) The technical solutions in the embodiments of the disclosure will be described clearly and completely in connection with the drawings of the embodiments of the disclosure.

(4) Referring to FIG. 1, which shows a flow of a method for frequency offset estimation provided by the embodiments of the disclosure. It should be noted that in order to clearly describe the technical solutions of the embodiments of the disclosure, the embodiments take one User Equipment (UE) as an example. It will be understood for those skilled in the art to apply the embodiments taking one UE as the example to embodiments with a plurality of UEs without creative labor. The method may include the following steps.

(5) At S101: initial channel estimation is performed by receiving a Sounding Reference Signal (SRS) or a Preamble code sent by a User Equipment (UE) to obtain a first frequency offset estimation value of a current sub-frame of the UE.

(6) Exemplarily, step S101 may specifically include:

(7) firstly, obtaining a channel estimation value H.sup.u and received data Y.sup.u of the UE by receiving a SRS or a Preamble code of a current sub-frame sent by the UE, and the specific process for the obtaining is conventional implementation means for those skilled in the art, which will not be repeated here;

(8) secondly, obtaining a time-domain sample value w.sup.u(n,ka.sub.Rx) of the current sub-frame of the UE on a receiving antenna ka.sub.Rx based on the channel estimation value H.sup.u and the received data Y.sup.u of the UE and in combination with a local pilot sequence P.sup.u and Formula 1:
w.sup.u(n,ka.sub.Rx)=IDFT(Y.sup.u(k,ka.sub.Rx)H.sup.u*(k,ka.sub.Rx))(IDFT(P.sup.u(k)))*(1)

(9) where n represents a series number of a sample point, u represents the UE, k represents a sub-carrier index corresponding to the UE 0k12N.sub.RB.sup.u1, ka.sub.Rx represents an index of the receiving antenna, N.sub.RB.sup.u represents a size of a Resource Block (RB) occupied by the UE, * represents a conjugate operation, and IDFT represents an inverse Fourier transform; and

(10) thirdly, obtaining a first frequency offset estimation value f.sub.0.sup.u(ka.sub.Rx) of the current sub-frame of the UE on the receiving antenna ka.sub.Rx according to phase information of a time-domain sample value w.sup.u(n,ka.sub.Rx) of the current sub-frame of the UE on the receiving antenna ka.sub.Rx. In detail, in the present embodiment, phase values of a first half and a second half of the time-domain sample value w.sup.u(n,ka.sub.Rx) of the current sub-frame can be averaged respectively to obtain a phase difference .sub.1.sup.u(ka.sub.Rx) of the first half and the second half, and .sub.1.sup.u(ka.sub.Rx) can be also obtained by averaging a phase difference of adjacent sample points of the time-domain sample value w.sup.u(n,ka.sub.Rx) of the current sub-frame. Then, a first frequency offset estimation value f.sub.0.sup.u(ka.sub.Rx) of the current sub-frame of the UE on the receiving antenna ka.sub.Rx is obtained according to Formula 2:

(11) $\begin{array}{cc}{f}_0^u\left({\mathrm{ka}}_{\mathrm{Rx}}\right)=\frac{{}_1^u\left({\mathrm{ka}}_{\mathrm{Rx}}\right)2f}{}.& \left(2\right)\end{array}$

(12) As one of the embodiments, a smoothing processing can be performed on f.sub.0.sup.u(ka.sub.Rx) obtained by Formula 2, so as to obtain a smoothed first frequency offset estimation value f.sub.0.sup.(u)(ka.sub.Rx) of the current sub-frame. Specifically, a weighted summation can be made on the first frequency offset estimation value f.sub.0.sup.u(ka.sub.Rx) corresponding to the sub-frame before the current sub-frame. It will be understood the weight becomes larger as it is closer to the current sub-frame in time, and the weight becomes smaller as it is further away from the current sub-frame in time.

(13) In the embodiments, a smoothing processing can be performed on f.sub.0.sup.u(ka.sub.Rx) by making the weighted summation on the basis of a first frequency offset estimation value f.sub.n.sub._.sub.subframe-1.sup.(u) corresponding to the previous sub-frame before the current sub-frame and f.sub.0.sup.u(ka.sub.Rx) according to Formula 3, to obtain a smoothed first frequency offset estimation value f.sub.0.sup.(u)(ka.sub.Rx) of the current sub-frame:
f.sub.1.sup.(u)(ka.sub.Rx)=(1p.sub.1)f.sub.n.sub._.sub.subframe-1.sup.(u)(ka.sub.Rx)+p.sub.1f.sub.0.sup.u(ka.sub.Rx)(3),

(14) where p.sub.1 is a filtering factor, and 0p.sub.11.

(15) Understandably, f.sub.1.sup.(u)(ka.sub.Rx) is smoother than f.sub.0.sup.u(ka.sub.Rx), but f.sub.0.sup.u(ka.sub.Rx) can be also used for subsequent processing, and taking f.sub.1.sup.(u)(ka.sub.Rx) for subsequent processing is only a particular solution in the embodiments; therefore, for a first frequency offset estimation value used during the subsequent processing, either one of f.sub.0.sup.u(ka.sub.Rx) and f.sub.1.sup.(u)(ka.sub.Rx) may be selected. However, there is no specific limitation in the embodiments of the disclosure.

(16) At S102: a frequency offset range in which the first frequency offset estimation value is located is determined, and an initial frequency offset value corresponding to the first frequency offset estimation value is obtained according to a preset correspondence between the frequency offset range and the initial frequency offset value.

(17) Exemplarily, the preset correspondence between the frequency offset range and the initial frequency offset value may be a specific table of mapping an upper and lower limit of the frequency offset range with an initial frequency offset value, as shown in Table 1.

(18) TABLE-US-00001 TABLE 1 Initial frequency Range number Upper limit: Hz Lower limit: Hz offset value: Hz 0 T1 T2 f0 1 T3 T4 f1 2 T5 T6 f2 3 T7 T8 f3 4 T9 T10 f4 5 T11 T12 f5

(19) It can be seen from Table 1, the first frequency offset estimation value can be compared with the upper and lower limit of each frequency offset range in Table 1 to determine the frequency offset range in which the first frequency offset estimation value is located, and then correspondingly an initial frequency offset value corresponding to the first frequency offset estimation value may be obtained. In the embodiments, f.sub.initial.sup.(u)(ka.sub.Rx) can be used to express the initial frequency offset value corresponding to the first frequency offset estimation value of the UE on the receiving antenna ka.sub.Rx.

(20) At S103: a phase difference of adjacent pilot locations is determined based on a received pilot sequence included in the current sub-frame, and a dual pilot frequency offset estimation value is obtain according to the phase difference within the range defined by the initial frequency offset value corresponding to the first frequency offset estimation value.

(21) As one of the embodiments, in this embodiment, S103 may specifically include:

(22) firstly, obtaining a phase difference .sub.2.sup.u(ka.sub.Rx) of two pilot locations of the UE on the receiving antenna ka.sub.Rx based on two pilot channel estimation values H.sub.1.sup.u(k,ka.sub.Rx) and H.sub.2.sup.u(k,ka.sub.Rx) obtained from two pilot sequences included in the current sub-frame of the UE on the receiving antenna ka.sub.Rx and according to Formula 3:

(23) $\begin{array}{cc}{}_2^u\left({\mathrm{ka}}_{\mathrm{Rx}}\right)=a\mathrm{ngle}\left(\underset{k=0}{\overset{12N_{\mathrm{RB}}^2-1}{.Math.}}{H}_1^u\left(k,{\mathrm{ka}}_{\mathrm{Rx}}\right)H_2^{u^{*}}\left(k,{\mathrm{ka}}_{\mathrm{Rx}}\right)\right),& \left(3\right)\end{array}$

(24) where * represents a conjugate operator, and angle( ) represents an operator for calculating an angle;

(25) secondly, obtaining a dual pilot frequency offset estimation result f.sub.1.sup.(u)(ka.sub.Rx) of the UE on the receiving antenna ka.sub.Rx based on a time interval t and a phase difference .sub.2.sup.u(ka.sub.Rx) of two pilot sequences and according to Formula 4:

(26) $\begin{array}{cc}{f}_2^{\left(u\right)}\left({\mathrm{ka}}_{\mathrm{Rx}}\right)=\frac{{}_2^u\left({\mathrm{ka}}_{\mathrm{Rx}}\right)}{2t};& \left(4\right)\end{array}$
and

(27) thirdly, obtaining, within the range defined by the initial frequency offset value f.sub.initial.sup.(u)(ka.sub.Rx), a dual pilot frequency offset estimation value f.sub.2(ka.sub.Rx) of the UE on the receiving antenna ka.sub.Rx according to the dual pilot frequency offset estimation result f.sub.2.sup.(u)(ka.sub.Rx) and Formula 5:
f.sub.2(ka.sub.Rx)=f.sub.2.sup.(u)(ka.sub.Rx)f.sub.initial.sup.(u)(ka.sub.Rx)+lT(5),

(28) where T represents the maximum range values for the dual pilot frequency offset estimation value, and l represents a preset integer for defining a dual pilot frequency offset estimation value in a predetermined range.

(29) At S104: a frequency offset estimation value of the current sub-frame is obtained according to the initial frequency offset value corresponding to the first frequency offset estimation value and the dual pilot frequency offset estimation value.

(30) Exemplarily, in the embodiments, a frequency offset estimation value f.sup.u(ka.sub.Rx) of the current sub-frame of the UE on the receiving antenna ka.sub.Rx may be obtained by adding the initial frequency offset value f.sub.initial.sup.(u)(ka.sub.Rx) and the dual pilot frequency offset estimation value f.sub.2(ka.sub.Rx).

(31) As one of the embodiments, after obtaining the frequency offset estimation value f.sup.u(ka.sub.Rx) of the current sub-frame of the UE on the receiving antenna ka.sub.Rx, a smoothing processing may be performed on f.sup.u(ka.sub.Rx) by means of filtering to obtain a smoothed frequency offset estimation value f.sup.(u)(ka.sub.Rx) of the current sub-frame.

(32) As one of the embodiments, in this embodiment, a smoothed frequency offset estimation value f.sup.(u)(ka.sub.Rx) of the current sub-frame can be obtained by performing a weighed summation on a frequency offset estimation value f.sub.n.sub._.sub.subframe-1.sup.u(ka.sub.Rx) of the previous sub-frame of the UE on the receiving antenna ka.sub.Rx and f.sup.(u)(ka.sub.Rx) according to Formula 6:
f.sup.(u)(ka.sub.Rx)=(1p.sub.2)f.sub.n.sub._.sub.subframe-1.sup.u(ka.sub.Rx)+p.sub.2f.sup.u(ka.sub.Rx)(6),

(33) where p.sub.2 is a filtering factor, and 0p.sub.21.

(34) The above scheme provides a method of frequency offset estimation for the embodiments of the disclosure, an initial frequency offset value is obtained according to a first frequency estimation value determined with a SRS or Preamble code sent by the UE, and a frequency offset estimation value of the current sub-frame is obtained within a range of the initial frequency offset estimation value and in combination with a dual pilot frequency offset estimation value. Thus, the precision of frequency offset estimation is ensured, and the complexity of large frequency offset estimation is reduced.

(35) Based on the same technical idea as that of the embodiment shown in FIG. 1, and referring to FIG. 2 which shows a device for frequency offset estimation 20 provided by the embodiments of the disclosure, the device 20 may include a first acquisition unit 201, a second acquisition unit 202, a third acquisition unit 203 and a fourth acquisition unit 204, in this device:

(36) the first acquisition unit 201 is configured to perform an initial channel estimation by receiving a SRS or Preamble code sent by a UE to obtain a first frequency offset estimation value of a current sub-frame of the UE;

(37) the second acquisition unit 202 is configured to determine a frequency range in which the first frequency offset estimation value is located and obtain an initial frequency offset value corresponding to the first frequency offset estimation value according to the frequency range and a preset correspondence between the frequency range and an initial frequency offset value;

(38) the third acquisition unit 203 is configured to calculate a phase difference of two pilot locations based on two received pilot sequences included in the current sub-frame of the UE and obtain, within a range defined by the initial frequency offset value, a duel pilot frequency offset estimation value according to the phase difference; and

(39) the fourth acquisition unit 204 is configured to obtain a frequency offset estimation value of the current sub-frame according to the initial frequency offset value and the dual pilot frequency offset estimation value.

(40) Exemplarily, the first acquisition unit 201 is configured to:

(41) obtain a channel estimation value H.sup.u and received data Y.sup.u of the UE by receiving a SRS or Preamble code of the current sub-frame sent by the UE;

(42) obtain a time-domain sample value w.sup.u(n,ka.sub.Rx) of the current sub-frame of the UE on the receiving antenna ka.sub.Rx based on the channel estimation value H.sup.u and the received data Y.sup.u of the UE and in combination with a local pilot sequence P.sup.u and according to the following formula:
w.sup.u(n,ka.sub.Rx)=IDFT(Y.sup.u(k,ka.sub.Rx)H.sup.u*(k,ka.sub.Rx))(IDFT(P.sup.u(k)))*,

(43) where n represents a series number of a sample point, u represents a User Equipment (UE) corresponding to the current sub-frame, k represents a sub-carrier index corresponding to the UE, 0k12N.sub.RB.sup.u1, ka.sub.Rx represents an index of the receiving antenna, N.sub.RB.sup.u represents a size of a Resource Block (RB) occupied by the UE, * represents a conjugate operation, and IDFT represents an inverse Fourier transform; and

(44) obtain a first frequency offset estimation value f.sub.0.sup.u(ka.sub.Rx) of the current sub-frame of the UE on the receiving antenna ka.sub.Rx according to the phase information of the time-domain sample value w.sup.u(n,ka.sub.Rx) of the current sub-frame of the UE on the receiving antenna ka.sub.Rx.

(45) As one of the embodiments, the first acquisition unit 201 is configured to:

(46) average respectively phase values of a first half and s second half of a time-domain sample value w.sup.u(n,ka.sub.Rx) of the current sub-frame to obtain a phase difference .sub.1.sup.u(ka.sub.Rx) of the first half and the second half; and

(47) obtain a first frequency offset estimation value f.sub.0.sup.u(ka.sub.Rx) of the current sub-frame of the UE on the receiving antenna ka.sub.Rx according to the following formula:

(48) 0$f_0^u\left({\mathrm{ka}}_{\mathrm{Rx}}\right)=\frac{{}_1^u\left({\mathrm{ka}}_{\mathrm{Rx}}\right)2f}{}.$

(49) As one of the embodiments, the first acquisition unit 201 is further configured to perform a smoothing processing on the first frequency offset estimation value f.sub.0.sup.u(ka.sub.Rx) of the current sub-frame of the UE on the receiving antenna ka.sub.Rx to obtain a smoothed first frequency offset estimation value f.sub.1.sup.(u)(ka.sub.Rx) of the current sub-frame. Specifically, the first acquisition unit 201 can perform a weighted summation on a first frequency offset estimation value f.sub.0.sup.(u)(ka.sub.Rx) corresponding to the sub-frame before the current sub-frame. Understandably, the weight becomes larger as it is closer to the current sub-frame in time, and the weight becomes smaller when it is further away from the current sub-frame in time.

(50) In the embodiments, the first acquisition unit 201 can perform a smoothing processing on f.sub.0.sup.u(ka.sub.Rx) by performing a weighted summation on the basis of a first frequency offset estimation value f.sub.n.sub._.sub.subframe-1.sup.(u) corresponding to the previous sub-frame before the current sub-frame and according to the following formula, to obtain a smoothed first frequency offset estimation value f.sub.1.sup.(u)(ka.sub.Rx) of the current sub-frame:
f.sub.1.sup.(u)(ka.sub.Rx)=(1p.sub.1)f.sub.n.sub._.sub.subframe-1.sup.(u)(ka.sub.Rx)+p.sub.1f.sub.0.sup.u(ka.sub.Rx),

(51) where p.sub.1 is a filtering factor, and 0p.sub.11.

(52) Understandably, f.sub.1.sup.(u)(ka.sub.Rx) is smoother than f.sub.0.sup.u(ka.sub.Rx), but f.sub.0.sup.u(ka.sub.Rx) can be also used for processing by a subsequent module, and taking f.sub.1.sup.(u)(ka.sub.Rx) used for processing by a subsequent module is only a particular solution in the embodiments; therefore, for a first frequency offset estimation value used in the processing by a subsequent module, either one of f.sub.0.sup.u(ka.sub.Rx) and f.sub.1.sup.(u)(ka.sub.Rx) may be selected. However, there is no specific limitation in the embodiments of the disclosure.

(53) Exemplarily, the preset correspondence between the frequency range and the initial frequency offset value can be a specific table of mapping an upper and lower limit of a frequency range and an initial frequency offset value, as shown in Table 1.

(54) As one of the embodiments, it can be seen from Table 1 that the second acquisition unit 202 can be configured to compare the first frequency offset estimation value with an upper limit and a lower limit of each frequency range in Table 1 in order to determine a frequency range in which the first frequency offset estimation value is located, and also correspondingly obtain an initial frequency offset value corresponding to the first frequency offset estimation value. In the embodiments, f.sub.initial.sup.(u)(ka.sub.Rx) can be used to express an initial frequency offset value corresponding to the first frequency offset estimation value of the UE on the receiving antenna ka.sub.Rx.

(55) Exemplarily, the third acquisition unit 203 is configured to:

(56) obtain a phase difference .sub.2.sup.u(ka.sub.Rx) of two pilot locations of the UE on the receiving antenna ka.sub.Rx based on two pilot channel estimation values, H.sub.1.sup.u(k,ka.sub.Rx) and H.sub.2.sup.u(k,ka.sub.Rx), obtained from two pilot sequences included in the current sub-frame of the UE on the receiving antenna ka.sub.Rx and according to the following formula:

(57) ${}_2^u\left({\mathrm{ka}}_{\mathrm{Rx}}\right)=a\mathrm{ngle}\left(\underset{k=0}{\overset{12N_{\mathrm{RB}}^2-1}{.Math.}}{H}_1^u\left(k,{\mathrm{ka}}_{\mathrm{Rx}}\right)H_2^{u^{*}}\left(k,{\mathrm{ka}}_{\mathrm{Rx}}\right)\right),$

(58) where * represents a conjugate operator, and angle( ) represents an operator for calculating an angle;

(59) obtain a dual pilot frequency offset estimation result f.sub.2.sup.(u)(ka.sub.Rx) of the UE on the receiving antenna ka.sub.Rx based on a time interval t and a phase difference .sub.2.sup.u(ka.sub.Rx) of the two pilot sequences and according to the following formula:

(60) $f_2^{\left(u\right)}\left({\mathrm{ka}}_{\mathrm{Rx}}\right)=\frac{{}_2^u\left({\mathrm{ka}}_{\mathrm{Rx}}\right)}{2t};$
and

(61) obtain, within the range defined by the initial frequency offset value, a dual pilot frequency offset estimation value f.sub.2(ka.sub.Rx) of the UE on the receiving antenna ka.sub.Rx according to the dual pilot frequency offset estimation result f.sub.2.sup.(u)(ka.sub.Rx) and the following formula:
f.sub.2(ka.sub.Rx)=f.sub.2.sup.(u)(ka.sub.Rx)f.sub.initial.sup.(u)(ka.sub.Rx)+lT,
where T represents the maximum range values of the dual pilot frequency offset estimation value, and l represents a preset integer for defining the dual pilot frequency offset estimation value within a predetermined range.

(62) As one of the embodiments, in this embodiment, the fourth acquisition unit 204 can obtain a frequency offset estimation value f.sup.u(ka.sub.Rx) of the current sub-frame of the UE on the receiving antenna ka.sub.Rx by adding the initial frequency offset value f.sub.initial.sup.(u)(ka.sub.Rx) and the dual pilot frequency offset estimation value f.sub.2(ka.sub.Rx).

(63) As one of the embodiments, after obtaining a frequency offset estimation value f.sub.initial.sup.(u)(ka.sub.Rx) of the current sub-frame of the UE on the receiving antenna ka.sub.Rx, the fourth acquisition unit 204 is further configured to perform a smoothing processing on the frequency offset estimation value f.sup.u(ka.sub.Rx) of the current sub-frame by means of filtering to obtain a smoothed frequency offset estimation value f.sup.(u)(ka.sub.Rx) of the current sub-frame.

(64) As one of the embodiments, in this embodiment, the fourth acquisition unit 204 can obtain a smoothed frequency offset estimation value f.sup.(u)(ka.sub.Rx) of the current sub-frame by performing a weighed summation on a frequency offset estimation value f.sub.n.sub._.sub.subframe-1.sup.u(ka.sub.Rx) of a previous sub-frame of the UE on the receiving antenna ka.sub.Rx and f.sup.u(ka.sub.Rx) according to the following formula:
f.sup.(u)(ka.sub.Rx)=(1p.sub.2)f.sub.n.sub._.sub.subframe-1.sup.u(ka.sub.Rx)+p.sub.2f.sup.u(ka.sub.Rx),

(65) where p.sub.2 is a filtering factor, and 0p.sub.21.

(66) The embodiments provide a device for frequency offset estimation 20, an initial frequency offset value is obtained according to a first frequency offset estimation value determined with a SRS or Preamble code sent by the UE, and a frequency offset estimation value of the current sub-frame is obtained within the range of the initial frequency offset value and in combination with a dual pilot frequency offset estimation value. Thus, the precision of frequency offset estimation is ensured, and the complexity of large frequency offset estimation is reduced.

(67) The embodiments of the disclosure further disclose a storage medium in which a computer program is stored, and the computer program is configured to execute the method for frequency offset estimation mentioned in above embodiments.

(68) Those skilled in the art should understand that the embodiments of the disclosure can provide a method, a system or a computer program product. Thus, forms of hardware embodiments, software embodiments or embodiments combining both software and hardware can be adopted in the disclosure. Moreover, a form of the computer program product implemented on one or more computer available storage media (including, without limitation to, a magnetic disk storage, an optical storage and the like) containing computer available program codes can be adopted in the disclosure.

(69) The disclosure is described with reference to flowcharts and/or block diagrams of the method, the device (system) and the computer program product according to the embodiments of the disclosure. It should be understood that each flow and/or block in the flowcharts and/or the block diagrams and a combination of the flows and/or the blocks in the flowcharts and/or the block diagrams can be realized by computer program instructions. These computer program instructions can be provided for a general computer, a dedicated computer, an embedded processor or processors of other programmable data processing devices to generate a machine, so that an apparatus for realizing functions assigned in one or more flows of the flowcharts and/or one or more blocks of the block diagrams is generated via instructions executed by the computers or the processors of the other programmable data processing devices.

(70) These computer program instructions can also be stored in a computer readable memory capable of guiding the computers or the other programmable data processing devices to work in a specific mode, so that a manufactured product including an instruction apparatus is generated via the instructions stored in the computer readable memory, and the instruction apparatus realizes the functions assigned in one or more flows of the flowcharts and/or one or more blocks of the block diagrams.

(71) These computer program instructions can also be loaded to the computers or the other programmable data processing devices, so that processing realized by the computers is generated by executing a series of operation steps on the computers or the other programmable devices, and therefore the instructions executed on the computers or the other programmable devices provide a step of realizing the functions assigned in one or more flows of the flowcharts and/or one or more blocks of the block diagrams.

(72) What are described above are only the particular embodiments of the disclosure, and are not used for limiting the protection scope of the disclosure.

INDUSTRIAL APPLICABILITY

(73) The disclosure can obtain an initial frequency offset value according to a first frequency offset estimation value determined with a SRS or Preamble code sent by a UE, and a frequency offset estimation value of the current sub-frame is obtained within a range of the initial frequency offset estimation value and in combination with a dual pilot frequency offset estimation value. Thus, the precision of frequency offset estimation is ensured, and the complexity of large frequency offset estimation is reduced.