# Patent

## Title
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).

## Background
TECHNICAL FIELD
The disclosure relates to the technology of communication signal processing, and in particular to a method, device and storage medium for frequency offset estimation.
BACKGROUND
A technology of Orthogonal Frequency Division Multiplexing (OFDM) is a kind of multicarrier transmission technology with subcarriers orthogonal to each other, which requires the subcarriers to be orthogonal to each other and is sensitive to the carrier frequency offset. Therefore, the problem of carrier frequency offset becomes one of the principal problems of the current OFDM technology.
At present, to solve the problem of a small frequency offset, there are many methods of a high precision and with simple implementation. However, for the problem of a large frequency offset, the precision for these solutions is relatively low, even though the precision of frequency offset estimation for some methods is improved, the complexity of the implementation are increased.
SUMMARY
To solve the above technical problem, embodiments of the disclosure intend to provide a method, device and storage medium for frequency offset estimation, which, for the problem of large frequency offset, ensure the precision of frequency offset estimation and reduce the complexity of estimation of large frequency offset.
The solutions claimed in the embodiments of the disclosure are implemented as follows.
In the first aspect, a method for frequency offset estimation is provided in the embodiments of the disclosure, the method may include:
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.
As one of the embodiments, 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 includes:
obtaining a channel estimation value Hu and received data Yu by receiving the SRS or Preamble code of the current sub-frame;
obtaining a time-domain sample value wu(n,kaRx) of the current sub-frame on a receiving antenna kaRx according to the channel estimation value Hu and the received data Yu and in combination with a local pilot sequence Pu and the following formula:

wu(n,kaRx)=IDFT(Yu(k,kaRx)×Hu*(k,kaRx))×(IDFT(Pu(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, 0≦k≦12NRBu−1, kaRx represents an index of the receiving antenna, NRBu 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 f0u(kaRx) of the current sub-frame on the receiving antenna kaRx according to phase information of the time-domain sample value wu(n,kaRx) of the current sub-frame on the receiving antenna kaRx.
As one of the embodiments, the obtaining the first frequency offset estimation value f0u(kaRx) of the current sub-frame on the receiving antenna kaRx according to the phase information of the time-domain sample value wu(n,kaRx) of the current sub-frame on the receiving antenna kaRx includes:
averaging respectively phase values of a first half and a second half of the time-domain sample value (n,kaRx) of the current sub-frame to obtain a phase difference Δφ1u(kaRx) of the first half and the second half; and
obtaining the first frequency offset estimation value f0u(kaRx) of the current sub-frame on the receiving antenna kaRx according to the following formula:




  
    
      f
      0
      u
    
    ⁡
    
      (
      
        ka
        Rx
      
      )
    
  
  =
  
    
      
        
          
            Δφ
            1
            u
          
          ⁡
          
            (
            
              ka
              Rx
            
            )
          
        
        ×
        2
        ×
        Δ
        ⁢
        
          
        
        ⁢
        f
      
      π
    
    .
  




As one of the embodiments, the method further includes:
performing a smoothing processing on the first frequency offset estimation value f0u(kaRx) of the current sub-frame on the receiving antenna kaRx to obtain a smoothed first frequency offset estimation value f1(u)(kaRx) of the current sub-frame.
As one of the embodiments, 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 includes:
obtaining a phase difference Δφ2u(kaRx) of two pilot locations on the receiving antenna kaRx based on two pilot channel estimation values H1u(k,kaRx) and H2u(k,kaRx) obtained from two pilot sequences included in the current sub-frame on the receiving antenna kaRx and according to the following formula:




  
    
      
        Δφ
        2
        u
      
      ⁡
      
        (
        
          ka
          Rx
        
        )
      
    
    =
    
      angle
      ⁢
      
        
      
      ⁢
      
        (
        
          
            ∑
            
              k
              =
              0
            
            
              
                12
                ⁢
                
                  N
                  RB
                  2
                
              
              -
              1
            
          
          ⁢
          
            
              
                H
                1
                u
              
              ⁡
              
                (
                
                  k
                  ,
                  
                    ka
                    Rx
                  
                
                )
              
            
            ×
            
              
                H
                2
                
                  u
                  *
                
              
              ⁡
              
                (
                
                  k
                  ,
                  
                    ka
                    Rx
                  
                
                )
              
            
          
        
        )
      
    
  
  ,




where * represents a conjugate operator, and angle( ) represents an operator for calculating an angle;
obtaining a dual pilot frequency offset estimation result Δf2(u)(kaRx) on the receiving antenna kaRx based on a time interval Δt and the phase difference Δφ2u(kaRx) of the two pilot sequences and according to the following formula:




  
    
      Δ
      ⁢
      
        
      
      ⁢
      
        
          f
          2
          
            (
            u
            )
          
        
        ⁡
        
          (
          
            ka
            Rx
          
          )
        
      
    
    =
    
      
        
          Δφ
          2
          u
        
        ⁡
        
          (
          
            ka
            Rx
          
          )
        
      
      
        2
        ⁢
        π
        ×
        Δ
        ⁢
        
          
        
        ⁢
        t
      
    
  
  ;




and

obtaining, within the range defined by the initial frequency offset value Δfinitial(u)(kaRx), the dual pilot frequency offset estimation value ΔΔf2(kaRx) on the receiving antenna kaRx based on the dual pilot frequency offset estimation result Δf2(u)(kaRx) and according to the following formula:

ΔΔf2(kaRx)=Δf2(u)(kaRx)−Δfinitial(u)(kaRx)+l×T, 

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.
As one of the embodiments, 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 includes:
performing a smoothing processing on the frequency offset estimation value fu(kaRx) of the current sub-frame by means of filtering to obtain a smoothed frequency offset estimation value f(u)(kaRx) of the current sub-frame.
In the second aspect, a device for frequency offset estimation is provided in the embodiments of the disclosure, and the device includes a first acquisition unit, a second acquisition unit, a third acquisition unit and a fourth acquisition unit,
the first acquisition unit 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;
the second acquisition unit is configured to 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;
the third acquisition unit is configured to 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
the fourth acquisition unit is configured to 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.
As one of the embodiments, the first acquisition unit is configured to:
obtain a channel estimation value Hu and received data Yu by receiving the SRS or Preamble code of the current sub-frame;
obtain a time-domain sample value wu(n,kaRx) of the current sub-frame on a receiving antenna kaRx according to the channel estimation value Hu and the received data Yu and in combination with a local pilot sequence Pu and the following formula:

wu(n,kaRx)=IDFT(Yu(k,kaRx)×Hu*(k,kaRx))×(IDFT(Pu(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, 0≦k≦12NRBu−1, kaRx represents an index of the receiving antenna, NRBu 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 f0u(kaRx) of the current sub-frame on the receiving antenna kaRx according to phase information of the time-domain sample value wu(n,kaRx) of the current sub-frame on the receiving antenna kaRx.
As one of embodiments, the first acquisition unit is configured to:
average respectively phase values of a first half and a second half of the time-domain sample value wu(n,kaRx) of the current sub-frame to obtain a phase difference Δφ1u(kaRx) of the first half and the second half; and
obtain the first frequency offset estimation value f0u(kaRx) of the current sub-frame on the receiving antenna kaRx according to the following formula:




  
    
      f
      0
      u
    
    ⁡
    
      (
      
        ka
        Rx
      
      )
    
  
  =
  
    
      
        
          
            Δφ
            1
            u
          
          ⁡
          
            (
            
              ka
              Rx
            
            )
          
        
        ×
        2
        ×
        Δ
        ⁢
        
          
        
        ⁢
        f
      
      π
    
    .
  




As one of the embodiments, the first acquisition unit is further configured to perform a smoothing processing on the first frequency offset estimation value f0u(kaRx) of the current sub-frame on the receiving antenna kaRx to obtain a smoothed first frequency offset estimation value f1(u)(kaRx) of the current sub-frame.
As one of the embodiments, the third acquisition unit is configured to:
obtain a phase difference Δφ2u(kaRx) of two pilot locations on the receiving antenna kaRx based on two pilot channel estimation values H1u(k,kaRx) and H2u(k,kaRx) obtained from two pilot sequences included in the current sub-frame on the receiving antenna kaRx and according to the following formula:




  
    
      
        Δφ
        2
        u
      
      ⁡
      
        (
        
          ka
          Rx
        
        )
      
    
    =
    
      angle
      ⁢
      
        
      
      ⁢
      
        (
        
          
            ∑
            
              k
              =
              0
            
            
              
                12
                ⁢
                
                  N
                  RB
                  2
                
              
              -
              1
            
          
          ⁢
          
            
              
                H
                1
                u
              
              ⁡
              
                (
                
                  k
                  ,
                  
                    ka
                    Rx
                  
                
                )
              
            
            ×
            
              
                H
                2
                
                  u
                  *
                
              
              ⁡
              
                (
                
                  k
                  ,
                  
                    ka
                    Rx
                  
                
                )
              
            
          
        
        )
      
    
  
  ,




where * represents a conjugate operator, and angle( ) represents an operator for calculating an angle;
obtain a dual pilot frequency offset estimation result Δf2(u)(kaRx) on the receiving antenna kaRx based on a time interval Δt and a phase difference Δφ2u(kaRx) of the two pilot sequences and according to the following formula:




  
    
      Δ
      ⁢
      
        
      
      ⁢
      
        
          f
          2
          
            (
            u
            )
          
        
        ⁡
        
          (
          
            ka
            Rx
          
          )
        
      
    
    =
    
      
        
          Δφ
          2
          u
        
        ⁡
        
          (
          
            ka
            Rx
          
          )
        
      
      
        2
        ⁢
        π
        ×
        Δ
        ⁢
        
          
        
        ⁢
        t
      
    
  
  ;




and

obtain, within the range defined by the initial frequency offset value Δfinitial(u)(kaRx), a dual pilot frequency offset estimation value ΔΔf2(kaRx) on the receiving antenna kaRx based on the dual pilot frequency offset estimation result Δf1(u)(kaRx) and according to the following formula:

ΔΔf2(kaRx)=Δf2(u)(kaRx)−Δfinitial(u)(kaRx)+l×T, 

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.
As one of the embodiments, the fourth acquisition unit is further configured to perform a smoothing processing on the frequency offset estimation value fu(kaRx) of the current sub-frame by means of filtering to obtain a smoothed frequency offset estimation value f(u)(kaRx) of the current sub-frame.
A storage medium, in which a computer program is stored, is provided, and the computer program is configured to execute the above method for frequency offset estimation.
In 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 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.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a flowchart of a method for frequency offset estimation according to an embodiment of the disclosure; and
FIG. 2 is a structural diagram of a device for frequency offset estimation according to an embodiment of the disclosure.



DETAILED DESCRIPTION
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.
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.
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.
Exemplarily, step S101 may specifically include:
firstly, obtaining a channel estimation value Hu and received data Yu 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;
secondly, obtaining a time-domain sample value wu(n,kaRx) of the current sub-frame of the UE on a receiving antenna kaRx based on the channel estimation value Hu and the received data Yu of the UE and in combination with a local pilot sequence Pu and Formula 1:

wu(n,kaRx)=IDFT(Yu(k,kaRx)×Hu*(k,kaRx))×(IDFT(Pu(k)))*  (1)

where n represents a series number of a sample point, u represents the UE, k represents a sub-carrier index corresponding to the UE 0≦k≦12NRBu−1, kaRx represents an index of the receiving antenna, NRBu represents a size of a Resource Block (RB) occupied by the UE, * represents a conjugate operation, and IDFT represents an inverse Fourier transform; and
thirdly, obtaining a first frequency offset estimation value f0u(kaRx) of the current sub-frame of the UE on the receiving antenna kaRx according to phase information of a time-domain sample value wu(n,kaRx) of the current sub-frame of the UE on the receiving antenna kaRx. In detail, in the present embodiment, phase values of a first half and a second half of the time-domain sample value wu(n,kaRx) of the current sub-frame can be averaged respectively to obtain a phase difference Δφ1u(kaRx) of the first half and the second half, and Δφ1u(kaRx) can be also obtained by averaging a phase difference of adjacent sample points of the time-domain sample value wu(n,kaRx) of the current sub-frame. Then, a first frequency offset estimation value f0u(kaRx) of the current sub-frame of the UE on the receiving antenna kaRx is obtained according to Formula 2:




  
    
      
        
          
            f
            0
            u
          
          ⁡
          
            (
            
              ka
              Rx
            
            )
          
        
        =
        
          
            
              
                
                  Δφ
                  1
                  u
                
                ⁡
                
                  (
                  
                    ka
                    Rx
                  
                  )
                
              
              ×
              2
              ×
              Δ
              ⁢
              
                
              
              ⁢
              f
            
            
          
          .
        
      
    
    
      
        (
        2
        )
      
    
  




As one of the embodiments, a smoothing processing can be performed on f0u(kaRx) obtained by Formula 2, so as to obtain a smoothed first frequency offset estimation value f0(u)(kaRx) of the current sub-frame. Specifically, a weighted summation can be made on the first frequency offset estimation value f0u(kaRx) 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.
In the embodiments, a smoothing processing can be performed on f0u(kaRx) by making the weighted summation on the basis of a first frequency offset estimation value fn_subframe-1(u) corresponding to the previous sub-frame before the current sub-frame and f0u(kaRx) according to Formula 3, to obtain a smoothed first frequency offset estimation value f0(u)(kaRx) of the current sub-frame:

f1(u)(kaRx)=(1−p1)fn_subframe-1(u)(kaRx)+p1f0u(kaRx)  (3),

where p1 is a filtering factor, and 0≦p1≦1.
Understandably, f1(u)(kaRx) is smoother than f0u(kaRx), but f0u(kaRx) can be also used for subsequent processing, and taking f1(u)(kaRx) 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 f0u(kaRx) and f1(u)(kaRx) may be selected. However, there is no specific limitation in the embodiments of the disclosure.
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.
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.










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









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, Δfinitial(u)(kaRx) 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 kaRx.
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.
As one of the embodiments, in this embodiment, S103 may specifically include:
firstly, obtaining a phase difference Δφ2u(kaRx) of two pilot locations of the UE on the receiving antenna kaRx based on two pilot channel estimation values H1u(k,kaRx) and H2u(k,kaRx) obtained from two pilot sequences included in the current sub-frame of the UE on the receiving antenna kaRx and according to Formula 3:




  
    
      
        
          
            
              Δφ
              2
              u
            
            ⁡
            
              (
              
                ka
                Rx
              
              )
            
          
          =
          
            
              a
              ⁢
              ngle
            
            (
            
              
                ∑
                
                  k
                  =
                  0
                
                
                  
                    12
                    ⁢
                    
                      N
                      RB
                      2
                    
                  
                  -
                  1
                
              
              ⁢
              
                
              
              ⁢
              
                
                  
                    H
                    1
                    u
                  
                  ⁡
                  
                    (
                    
                      k
                      ,
                      
                        ka
                        Rx
                      
                    
                    )
                  
                
                ×
                
                  
                    H
                    2
                    
                      u
                      *
                    
                  
                  ⁡
                  
                    (
                    
                      k
                      ,
                      
                        ka
                        Rx
                      
                    
                    )
                  
                
              
            
            )
          
        
        ,
      
    
    
      
        (
        3
        )
      
    
  




where * represents a conjugate operator, and angle( ) represents an operator for calculating an angle;
secondly, obtaining a dual pilot frequency offset estimation result Δf1(u)(kaRx) of the UE on the receiving antenna kaRx based on a time interval Δt and a phase difference Δφ2u(kaRx) of two pilot sequences and according to Formula 4:




  
    
      
        
          
            Δ
            ⁢
            
              
            
            ⁢
            
              
                f
                2
                
                  (
                  u
                  )
                
              
              ⁡
              
                (
                
                  ka
                  Rx
                
                )
              
            
          
          =
          
            
              
                Δφ
                2
                u
              
              ⁡
              
                (
                
                  ka
                  Rx
                
                )
              
            
            
              2
              ⁢
              
              ×
              Δt
            
          
        
        ;
      
    
    
      
        (
        4
        )
      
    
  




and

thirdly, obtaining, within the range defined by the initial frequency offset value Δfinitial(u)(kaRx), a dual pilot frequency offset estimation value ΔΔf2(kaRx) of the UE on the receiving antenna kaRx according to the dual pilot frequency offset estimation result Δf2(u)(kaRx) and Formula 5:

ΔΔf2(kaRx)=Δf2(u)(kaRx)−Δfinitial(u)(kaRx)+l×T  (5),

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.
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.
Exemplarily, in the embodiments, a frequency offset estimation value fu(kaRx) of the current sub-frame of the UE on the receiving antenna kaRx may be obtained by adding the initial frequency offset value Δfinitial(u)(kaRx) and the dual pilot frequency offset estimation value ΔΔf2(kaRx).
As one of the embodiments, after obtaining the frequency offset estimation value fu(kaRx) of the current sub-frame of the UE on the receiving antenna kaRx, a smoothing processing may be performed on fu(kaRx) by means of filtering to obtain a smoothed frequency offset estimation value f(u)(kaRx) of the current sub-frame.
As one of the embodiments, in this embodiment, a smoothed frequency offset estimation value f(u)(kaRx) of the current sub-frame can be obtained by performing a weighed summation on a frequency offset estimation value fn_subframe-1u(kaRx) of the previous sub-frame of the UE on the receiving antenna kaRx and f(u)(kaRx) according to Formula 6:

f(u)(kaRx)=(1−p2)fn_subframe-1u(kaRx)+p2•fu(kaRx)  (6),

where p2 is a filtering factor, and 0≦p2≦1.
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.
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:
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;
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;
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
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.
Exemplarily, the first acquisition unit 201 is configured to:
obtain a channel estimation value Hu and received data Yu of the UE by receiving a SRS or Preamble code of the current sub-frame sent by the UE;
obtain a time-domain sample value wu(n,kaRx) of the current sub-frame of the UE on the receiving antenna kaRx based on the channel estimation value Hu and the received data Yu of the UE and in combination with a local pilot sequence Pu and according to the following formula:

wu(n,kaRx)=IDFT(Yu(k,kaRx)×Hu*(k,kaRx))×(IDFT(Pu(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, 0≦k≦12NRBu−1, kaRx represents an index of the receiving antenna, NRBu 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 a first frequency offset estimation value f0u(kaRx) of the current sub-frame of the UE on the receiving antenna kaRx according to the phase information of the time-domain sample value wu(n,kaRx) of the current sub-frame of the UE on the receiving antenna kaRx.
As one of the embodiments, the first acquisition unit 201 is configured to:
average respectively phase values of a first half and s second half of a time-domain sample value wu(n,kaRx) of the current sub-frame to obtain a phase difference Δφ1u(kaRx) of the first half and the second half; and
obtain a first frequency offset estimation value f0u(kaRx) of the current sub-frame of the UE on the receiving antenna kaRx according to the following formula:




  
    
      f
      0
      u
    
    ⁡
    
      (
      
        ka
        Rx
      
      )
    
  
  =
  
    
      
        
          
            Δφ
            1
            u
          
          ⁡
          
            (
            
              ka
              Rx
            
            )
          
        
        ×
        2
        ×
        Δ
        ⁢
        
          
        
        ⁢
        f
      
      
    
    .
  




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 f0u(kaRx) of the current sub-frame of the UE on the receiving antenna kaRx to obtain a smoothed first frequency offset estimation value f1(u)(kaRx) of the current sub-frame. Specifically, the first acquisition unit 201 can perform a weighted summation on a first frequency offset estimation value f0(u)(kaRx) 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.
In the embodiments, the first acquisition unit 201 can perform a smoothing processing on f0u(kaRx) by performing a weighted summation on the basis of a first frequency offset estimation value fn_subframe-1(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 f1(u)(kaRx) of the current sub-frame:

f1(u)(kaRx)=(1−p1)fn_subframe-1(u)(kaRx)+p1f0u(kaRx),

where p1 is a filtering factor, and 0≦p1≦1.
Understandably, f1(u)(kaRx) is smoother than f0u(kaRx), but f0u(kaRx) can be also used for processing by a subsequent module, and taking f1(u)(kaRx) 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 f0u(kaRx) and f1(u)(kaRx) may be selected. However, there is no specific limitation in the embodiments of the disclosure.
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.
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, Δfinitial(u)(kaRx) 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 kaRx.
Exemplarily, the third acquisition unit 203 is configured to:
obtain a phase difference Δφ2u(kaRx) of two pilot locations of the UE on the receiving antenna kaRx based on two pilot channel estimation values, H1u(k,kaRx) and H2u(k,kaRx), obtained from two pilot sequences included in the current sub-frame of the UE on the receiving antenna kaRx and according to the following formula:




  
    
      
        Δφ
        2
        u
      
      ⁡
      
        (
        
          ka
          Rx
        
        )
      
    
    =
    
      
        a
        ⁢
        ngle
      
      (
      
        
          ∑
          
            k
            =
            0
          
          
            
              12
              ⁢
              
                N
                RB
                2
              
            
            -
            1
          
        
        ⁢
        
          
        
        ⁢
        
          
            
              H
              1
              u
            
            ⁡
            
              (
              
                k
                ,
                
                  ka
                  Rx
                
              
              )
            
          
          ×
          
            
              H
              2
              
                u
                *
              
            
            ⁡
            
              (
              
                k
                ,
                
                  ka
                  Rx
                
              
              )
            
          
        
      
      )
    
  
  ,




where * represents a conjugate operator, and angle( ) represents an operator for calculating an angle;
obtain a dual pilot frequency offset estimation result Δf2(u)(kaRx) of the UE on the receiving antenna kaRx based on a time interval Δt and a phase difference Δφ2u(kaRx) of the two pilot sequences and according to the following formula:




  
    
      Δ
      ⁢
      
        
      
      ⁢
      
        
          f
          2
          
            (
            u
            )
          
        
        ⁡
        
          (
          
            ka
            Rx
          
          )
        
      
    
    =
    
      
        
          Δφ
          2
          u
        
        ⁡
        
          (
          
            ka
            Rx
          
          )
        
      
      
        2
        ⁢
        
        ×
        Δt
      
    
  
  ;




and

obtain, within the range defined by the initial frequency offset value, a dual pilot frequency offset estimation value ΔΔf2(kaRx) of the UE on the receiving antenna kaRx according to the dual pilot frequency offset estimation result Δf2(u)(kaRx) and the following formula:

ΔΔf2(kaRx)=Δf2(u)(kaRx)−Δfinitial(u)(kaRx)+l×T, 

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.

As one of the embodiments, in this embodiment, the fourth acquisition unit 204 can obtain a frequency offset estimation value fu(kaRx) of the current sub-frame of the UE on the receiving antenna kaRx by adding the initial frequency offset value Δfinitial(u)(kaRx) and the dual pilot frequency offset estimation value ΔΔf2(kaRx).
As one of the embodiments, after obtaining a frequency offset estimation value Δfinitial(u)(kaRx) of the current sub-frame of the UE on the receiving antenna kaRx, the fourth acquisition unit 204 is further configured to perform a smoothing processing on the frequency offset estimation value fu(kaRx) of the current sub-frame by means of filtering to obtain a smoothed frequency offset estimation value f(u)(kaRx) of the current sub-frame.
As one of the embodiments, in this embodiment, the fourth acquisition unit 204 can obtain a smoothed frequency offset estimation value f(u)(kaRx) of the current sub-frame by performing a weighed summation on a frequency offset estimation value fn_subframe-1u(kaRx) of a previous sub-frame of the UE on the receiving antenna kaRx and fu(kaRx) according to the following formula:

f(u)(kaRx)=(1−p2)fn_subframe-1u(kaRx)+p2•fu(kaRx),

where p2 is a filtering factor, and 0≦p2≦1.
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.
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.
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.
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.
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.
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.
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
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.

## 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 Hu and received data Yu by receiving the SRS or Preamble code of the current sub-frame;
obtaining a time-domain sample value wu(n,kaRx) of the current sub-frame on a receiving antenna kaRx according to the channel estimation value Hu and the received data Yu and in combination with a local pilot sequence Pu and the following formula:


wu(n,kaRx)=IDFT(Yu(k,kaRx)×Hu*(k,kaRx))×(IDFT(Pu(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, 0≦k≦12NRBu−1, kaRx represents an index of the receiving antenna, NRBu 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 f0u(kaRx) of the current sub-frame on the receiving antenna kaRx according to phase information of the time-domain sample value wu(n,kaRx) of the current sub-frame on the receiving antenna kaRx.

3. The method according to claim 2, wherein the obtaining the first frequency offset estimation value f0u(kaRx) of the current sub-frame on the receiving antenna kaRx according to the phase information of the time-domain sample value wu(n,kaRx) of the current sub-frame on the receiving antenna kaRx comprises:
averaging respectively phase values of a first half and a second half of the time-domain sample value wu(n,kaRx) of the current sub-frame to obtain a phase difference Δφ1u(kaRx) of the first half and the second half; and
obtaining the first frequency offset estimation value f0u(kaRx) of the current sub-frame on the receiving antenna kaRx according to the following formula:





  
    
      f
      0
      u
    
    ⁡
    
      (
      
        ka
        Rx
      
      )
    
  
  =
  
    
      
        
          
            Δφ
            1
            u
          
          ⁡
          
            (
            
              ka
              Rx
            
            )
          
        
        ×
        2
        ×
        Δ
        ⁢
        
          
        
        ⁢
        f
      
      
    
    .

4. The method according to claim 2, wherein the method further comprises:
performing a smoothing processing on the first frequency offset estimation value f0u(kaRx) of the current sub-frame on the receiving antenna kaRx to obtain a smoothed first frequency offset estimation value f1(u)(kaRx) 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 Δφ2u(kaRx) of two pilot locations on the receiving antenna kaRx based on two pilot channel estimation values H1u(k,kaRx) and H2u(k,kaRx) obtained from two pilot sequences included in the current sub-frame on the receiving antenna kaRx and according to the following formula:





  
    
      
        Δφ
        2
        u
      
      ⁡
      
        (
        
          ka
          Rx
        
        )
      
    
    =
    
      
        a
        ⁢
        ngle
      
      (
      
        
          ∑
          
            k
            =
            0
          
          
            
              12
              ⁢
              
                N
                RB
                2
              
            
            -
            1
          
        
        ⁢
        
          
        
        ⁢
        
          
            
              H
              1
              u
            
            ⁡
            
              (
              
                k
                ,
                
                  ka
                  Rx
                
              
              )
            
          
          ×
          
            
              H
              2
              
                u
                *
              
            
            ⁡
            
              (
              
                k
                ,
                
                  ka
                  Rx
                
              
              )
            
          
        
      
      )
    
  
  ,



where * represents a conjugate operator, and angle( ) represents an operator for calculating an angle;
obtaining a dual pilot frequency offset estimation result Δf2(u)(kaRx) on the receiving antenna kaRx based on a time interval Δt and the phase difference Δφ2u(kaRx) of the two pilot sequences and according to the following formula:





  
    
      Δ
      ⁢
      
        
      
      ⁢
      
        
          f
          2
          
            (
            u
            )
          
        
        ⁡
        
          (
          
            ka
            Rx
          
          )
        
      
    
    =
    
      
        
          Δφ
          2
          u
        
        ⁡
        
          (
          
            ka
            Rx
          
          )
        
      
      
        2
        ⁢
        
        ×
        Δt
      
    
  
  ;




and
obtaining, within the range defined by the initial frequency offset value Δfinitial(u)(kaRx), the dual pilot frequency offset estimation value ΔΔf2(kaRx) on the receiving antenna kaRx based on the dual pilot frequency offset estimation result Δf2(u)(kaRx) and according to the following formula:


ΔΔf2(kaRx)=Δf2(u)(kaRx)−Δfinitial(u)(kaRx)+l×T, 


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 fu(kaRx) of the current sub-frame by means of filtering to obtain a smoothed frequency offset estimation value f(u)(kaRx) 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 Hu and received data Yu by receiving the SRS or Preamble code of the current sub-frame;
obtain a time-domain sample value wu(n,kaRx) of the current sub-frame on a receiving antenna kaRx according to the channel estimation value Hu and the received data Yu and in combination with a local pilot sequence Pu and the following formula:


wu(n,kaRx)=IDFT(Yu(k,kaRx)×Hu*(k,kaRx))×(IDFT(Pu(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, 0≦k≦12NRBu−1, kaRx represents an index of the receiving antenna, NRBu 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 f0u(kaRx) of the current sub-frame on the receiving antenna kaRx according to phase information of the time-domain sample value wu(n,kaRx) of the current sub-frame on the receiving antenna kaRx.

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 wu(n,kaRx) of the current sub-frame to obtain a phase difference Δφ1u(kaRx) of the first half and the second half; and
obtain the first frequency offset estimation value f0u(kaRx) of the current sub-frame on the receiving antenna kaRx according to the following formula:





  
    
      f
      0
      u
    
    ⁡
    
      (
      
        ka
        Rx
      
      )
    
  
  =
  
    
      
        
          
            Δφ
            1
            u
          
          ⁡
          
            (
            
              ka
              Rx
            
            )
          
        
        ×
        2
        ×
        Δ
        ⁢
        
          
        
        ⁢
        f
      
      
    
    .

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 f0u(kaRx) of the current sub-frame on the receiving antenna kaRx to obtain a smoothed first frequency offset estimation value f1(u)(kaRx) of the current sub-frame.

11. The device according to claim 7, wherein the processor is further configured to:
obtain a phase difference Δφ2u(kaRx) of two pilot locations on the receiving antenna kaRx based on two pilot channel estimation values H1u(k,kaRx) and H2u(k,kaRx) obtained from two pilot sequences included in the current sub-frame on the receiving antenna kaRx and according to the following formula:





  
    
      
        Δφ
        2
        u
      
      ⁡
      
        (
        
          ka
          Rx
        
        )
      
    
    =
    
      
        a
        ⁢
        ngle
      
      (
      
        
          ∑
          
            k
            =
            0
          
          
            
              12
              ⁢
              
                N
                RB
                2
              
            
            -
            1
          
        
        ⁢
        
          
        
        ⁢
        
          
            
              H
              1
              u
            
            ⁡
            
              (
              
                k
                ,
                
                  ka
                  Rx
                
              
              )
            
          
          ×
          
            
              H
              2
              
                u
                *
              
            
            ⁡
            
              (
              
                k
                ,
                
                  ka
                  Rx
                
              
              )
            
          
        
      
      )
    
  
  ,



where * represents a conjugate operator, and angle( ) represents an operator for calculating an angle;
obtain a dual pilot frequency offset estimation result Δf2(u)(kaRx) on the receiving antenna kaRx based on a time interval Δt and a phase difference Δφ2u(kaRx) of the two pilot sequences and according to the following formula:





  
    
      Δ
      ⁢
      
        
      
      ⁢
      
        
          f
          2
          
            (
            u
            )
          
        
        ⁡
        
          (
          
            ka
            Rx
          
          )
        
      
    
    =
    
      
        
          Δφ
          2
          u
        
        ⁡
        
          (
          
            ka
            Rx
          
          )
        
      
      
        2
        ⁢
        
        ×
        Δt
      
    
  
  ;




and
obtain, within the range defined by the initial frequency offset value Δfinitial(u)(kaRx), a dual pilot frequency offset estimation value ΔΔf2(kaRx) on the receiving antenna kaRx based on the dual pilot frequency offset estimation result Δf2(u)(kaRx) and according to the following formula:


ΔΔf2(kaRx)=Δf2(u)(kaRx)−Δfinitial(u)(kaRx)+l×T, 


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 fu(kaRx) of the current sub-frame by means of filtering to obtain a smoothed frequency offset estimation value f(u)(kaRx) 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 Hu and received data Yu by receiving the SRS or Preamble code of the current sub-frame;
obtaining a time-domain sample value wu(n,kaRx) of the current sub-frame on a receiving antenna kaRx according to the channel estimation value Hu and the received data Yu and in combination with a local pilot sequence Pu and the following formula:


wu(n,kaRx)=IDFT(Yu(k,kaRx)×Hu*(k,kaRx))×(IDFT(Pu(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, 0≦k≦12NRBu−1, kaRx represents an index of the receiving antenna, NRBu 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 f0u(kaRx) of the current sub-frame on the receiving antenna kaRx according to phase information of the time-domain sample value wu(n,kaRx) of the current sub-frame on the receiving antenna kaRx.

15. The non-transitory computer readable storage medium according to claim 14, wherein obtaining the first frequency offset estimation value f0u(kaRx) of the current sub-frame on the receiving antenna kaRx according to the phase information of the time-domain sample value wu(n,kaRx) of the current sub-frame on the receiving antenna kaRx comprises:
averaging respectively phase values of a first half and a second half of the time-domain sample value wu(n,kaRx) of the current sub-frame to obtain a phase difference Δφ1u(kaRx) of the first half and the second half; and
obtaining the first frequency offset estimation value f0u(kaRx) of the current sub-frame on the receiving antenna kaRx according to the following formula:





  
    
      f
      0
      u
    
    ⁡
    
      (
      
        ka
        Rx
      
      )
    
  
  =
  
    
      
        
          
            Δφ
            1
            u
          
          ⁡
          
            (
            
              ka
              Rx
            
            )
          
        
        ×
        2
        ×
        Δ
        ⁢
        
          
        
        ⁢
        f
      
      
    
    .

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 f0u(kaRx) of the current sub-frame on the receiving antenna kaRx to obtain a smoothed first frequency offset estimation value f1(u)(kaRx) 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 Δφ2u(kaRx) of two pilot locations on the receiving antenna kaRx based on two pilot channel estimation values H1u(k,kaRx) and H2u(k,kaRx) obtained from two pilot sequences included in the current sub-frame on the receiving antenna kaRx and according to the following formula:





  
    
      
        Δφ
        2
        u
      
      ⁡
      
        (
        
          ka
          Rx
        
        )
      
    
    =
    
      
        a
        ⁢
        ngle
      
      (
      
        
          ∑
          
            k
            =
            0
          
          
            
              12
              ⁢
              
                N
                RB
                2
              
            
            -
            1
          
        
        ⁢
        
          
        
        ⁢
        
          
            
              H
              1
              u
            
            ⁡
            
              (
              
                k
                ,
                
                  ka
                  Rx
                
              
              )
            
          
          ×
          
            
              H
              2
              
                u
                *
              
            
            ⁡
            
              (
              
                k
                ,
                
                  ka
                  Rx
                
              
              )
            
          
        
      
      )
    
  
  ,



where * represents a conjugate operator, and angle( ) represents an operator for calculating an angle;
obtaining a dual pilot frequency offset estimation result Δf2(u)(kaRx) on the receiving antenna kaRx based on a time interval Δt and the phase difference Δφ2u(kaRx) of the two pilot sequences and according to the following formula:





  
    
      Δ
      ⁢
      
        
      
      ⁢
      
        
          f
          2
          
            (
            u
            )
          
        
        ⁡
        
          (
          
            ka
            Rx
          
          )
        
      
    
    =
    
      
        
          Δφ
          2
          u
        
        ⁡
        
          (
          
            ka
            Rx
          
          )
        
      
      
        2
        ⁢
        
        ×
        Δt
      
    
  
  ;




and
obtaining, within the range defined by the initial frequency offset value Δfinitial(u)(kaRx), the dual pilot frequency offset estimation value ΔΔf2(kaRx) on the receiving antenna kaRx based on the dual pilot frequency offset estimation result Δf2(u)(kaRx) and according to the following formula:


ΔΔf2(kaRx)=Δf2(u)(kaRx)−Δfinitial(u)(kaRx)+l×T, 


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 fu(kaRx) of the current sub-frame by means of filtering to obtain a smoothed frequency offset estimation value f(u)(kaRx) 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.