METHOD FOR ESTIMATING AND CORRECTING A DRIFT BETWEEN CLOCKS OF A RECEIVING TRANSCEIVER AND A CORRESPONDING EMITTING TRANSCEIVER, AND RECEIVER FOR IMPLEMENTING SAID METHOD

• US 20110051847A1
• Filed: 09/01/2010
• Published: 03/03/2011
• Est. Priority Date: 09/01/2009
• Status: Active Grant

• Current Assignee: Ecole polytechnique fãdãrale de lausanne epfl (Schweizerische Eidgenossenschaft)
• Original Assignee: Ecole Polytechnique Federal de Lausanne
• Inventors: Merz, Ruben, FLURY, Manuel, Le Boudec, Jean-Yves
• Granted Patent: US 8,396,175 B2
• Time in Patent Office: Days
• US Class Current: 375/316
• CPC Class Codes: H04L 7/0054 Detection of the synchronisation error by features other than the received signal transition by means of signal transition H04L7/033

Abstract

Method for estimating and correcting a frequency offset between a local clock of a receiving transceiver and a remote clock of an emitting transceiver comprises receiving by said receiving transceiver, a signal emitted by the emitting transceiver, said signal lacking a continuous sinusoidal carrier and having at least one symbol, each symbol being emitted within a time frame having a first duration; sampling the signal by taking for each time frame a set of energy measurements of said signal; forming a matrix representative of said energy measurements; determining from said energy matrix, points of maximum of energy within each time frame; determining a line connecting a plurality of said points within each of said time frames; determining a slope of said line; calculating a frequency offset between the clocks of the emitting and receiving transceivers using the slope; and adjusting the receiving transceiver'"'"'s clock frequency according to the frequency offset.

Claims

1. Method for estimating and correcting a frequency offset between a local clock of a receiving transceiver and a remote clock of a corresponding emitting transceiver, said receiving transceiver receiving from the emitting transceiver electromagnetic signals that lack a continuous sinusoidal carrier, said method comprising the steps of:

receiving by said receiving transceiver, a signal emitted by the emitting transceiver, said signal lacking a continuous sinusoidal carrier and being formed of at least one symbol, each symbol being emitted within a time frame having a first duration T_f^pre;

sampling the said received signal by taking for each frame a set of energy measurements of said received signal;

forming a matrix representative of said energy measurements;

determining from said energy matrix, points of maximum energy within each of said time frames;

determining a straight line connecting a plurality of said points of maximum energy;

determining a slope of said straight line;

calculating a frequency offset between the clock of the emitting transceiver and the clock of the receiving transceiver by using the slope of said straight line;

adjusting the frequency of the receiving transceiver'"'"'s clock according to the calculated frequency offset.

2. Method according to claim 1, wherein the step of forming said energy matrix comprises the step of introducing in two subsequent entries of a column of said energy matrix, two subsequent energy measurements obtained by sampling the received signal.

3. Method according to claim 1, wherein the step of forming said energy matrix comprises the step of introducing in two subsequent entries of a row of said energy matrix, two energy measurements obtained by sampling the received signal at two times differing by a multiple of said first duration T_f^pre.

4. Method according to claim 1, wherein the step of determining from said energy matrix a straight line connecting a plurality of points of maximum of energy comprises the step of using one of:

a maximum likelihood estimation technique;

a least square technique;

a general fitting technique;

a line detection technique used in image processing.

5. Method according to claim 4, wherein the step of determining from said energy matrix a straight line connecting a plurality of points of maximum of energy uses a line detection technique and wherein said line detection technique comprises the step of:

transforming lines of the energy matrix in points of a Radon space, wherein the coordinates (ρ

,θ

) of a point in said Radon space corresponds to a line of slope θ and

of distance to origin ρ

in said energy matrix;

integrating the energy along each line of the energy matrix and attributing an energy value corresponding to the integrated energy to the corresponding point of the Radon space and obtaining a Radon matrix;

determining at least one point of maximal energy value in the Radon matrix;

determining a column of the Radon matrix corresponding to the estimated clock drift and containing at least one of said at least one point of maximal energy value.

6. Method according to claim 5, wherein the values of the Radon matrix are convolved with a vector of predetermined length.

7. Method according to claim 5, wherein the step of determining said at least one point of maximal value in the Radon space comprises a step of squaring the intensity values of the Radon matrix.

8. Method according to claim 5, wherein the integration of the energy is performed while assuming that the energy is distributed in each cell of the energy matrix according to a specific distribution.

9. Method according to claim 8, wherein said specific distribution corresponds to one of:

all the energy in one cell is in the center of the cell of the energy matrix;

the energy is uniformly distributed in the cell of the energy matrix;

or

there is a Gaussian distribution of the energy in the cell of the energy matrix.

10. Method according to claim 5, wherein energy measurements below a threshold are dismissed prior to calculating the Radon transform.

11. Method according to claim 5, wherein energy measurements above a threshold are dismissed prior to calculating the Radon transform.

12. Method according to claim 5, wherein several energy measurements are combined using a moving average operation prior to calculating the Radon transform.

13. Method according to claim 12, wherein the moving average operation uses the median.

14. Method according to claim 12, wherein energy measurements above a threshold are set to a common value prior to averaging.

15. Method according to claim 5, wherein the column of the Radon matrix corresponding to the estimated clock drift is used to determine an estimate of the channel energy-delay profile.

16. Method according to claim 5, wherein the point of the Radon matrix with maximal value is used as a coarse estimate of the time of arrival of the signal.

17. Method according to claim 16, wherein the estimated time of arrival is refined by a search-back procedure starting from the coarse estimate of the time of arrival.

31. A receiver comprising:

an energy detection receiver structure;

a memory; and

a processor connected to the memory, the processor being configured to perform the steps of claim 1.

18. Method for estimating and correcting a frequency offset between a local clock of a receiving transceiver and a remote clock of a corresponding emitting transceiver, said receiving transceiver receiving from the emitting transceiver electromagnetic signals that lack a continuous sinusoidal carrier, said method comprising the steps of:

receiving by said receiving transceiver, a signal emitted by the emitting transceiver, said signal lacking a continuous sinusoidal carrier and being formed of at least one symbol, each symbol being emitted within a time frame having a first duration T_f^pre;

sampling the said received signal by taking for each frame a set of energy measurements of said received signal;

defining an initial Radon matrix;

for each energy measurement, calculating a contribution of this energy measurement to an entry of the Radon matrix;

updating the Radon matrix by adding said contribution;

determining at least one point of maximal energy value in the Radon space;

determining a column of the Radon matrix containing at least one of said at least one point of maximal energy value;

calculating a frequency offset between the clock of the emitting transceiver and the clock of the receiving transceiver by using the determined column of the Radon matrix;

adjusting the frequency of the receiving transceiver'"'"'s clock according to the calculated frequency offset.

19. Method according to claim 18, wherein the values of the Radon space are convolved with a vector of predetermined length.

20. Method according to claim 19, wherein the step of determining said at least one point of maximal value in the Radon space comprises a step of squaring the intensity values of the Radon matrix.

21. Method according to claim 18, wherein the integration of the energy is performed while assuming that the energy is distributed in each cell of the energy matrix according to a specific distribution.

22. Method according to claim 21, wherein said specific distribution corresponds to one of:

all the energy in one cell is in the center of the cell of the energy matrix;

the energy is uniformly distributed in the cell of the energy matrix;

or

there is a Gaussian distribution of the energy in the cell of the energy matrix.

23. Method according to claim 18, wherein energy measurement below a threshold are dismissed prior to calculating the Radon transform.

24. Method according to claim 18, wherein energy measurement above a threshold are dismissed prior to calculating the Radon transform.

25. Method according to claim 18, wherein several energy measurement are combined using a moving average operation prior to calculating the Radon transform.

26. Method according to claim 25, wherein the moving average operation uses the median.

27. Method according to claim 25, wherein energy measurement above a threshold are set to a common value prior to averaging.

28. Method according to claim 18, wherein the column of the Radon matrix corresponding to the estimated clock drift is used to determine an estimate of the channel energy-delay profile.

29. Method according to claim 18, wherein the point of the Radon transform with maximal value is used as a coarse estimate of the time of arrival of the signal.

30. Method according to claim 29, wherein the estimated time of arrival is refined by a search-back procedure starting from the coarse estimate of the time of arrival.

32. A receiver comprising:

an energy detection receiver structure;

a memory; and

a processor connected to the memory, the processor being configured to perform the steps of claim 18.