DISTANCE MEASUREMENT DEVICE AND DISTANCE MEASUREMENT METHOD

Abstract

A distance measurement device of an embodiment includes a first device including a first reference signal source and a first transmitter-receiver, a second device including a second reference signal source and a second transmitter-receiver, and a calculation unit configured to calculate a distance between the first device and the second device. One of a first distance measurement signal and a second distance measurement signal is transmitted once or more, and another is transmitted twice or more. The calculation unit calculates the distance based on a total of three or more pieces of first phase information and second phase information acquired through transmission of the distance measurement signals three times or more in total, a first sampling period based on a first reference signal, and a second sampling period based on a second reference signal.

Description

CROSS-REFERENCE TO RELATED APPLICATIONS

This application is based upon and claims the benefit of priority from the prior Japanese Patent Application No. 2020-159572 filed in Japan on Sep. 24, 2020; the entire contents of which are incorporated herein by reference.

FIELD

Embodiments described herein relate generally to a distance measurement device and a distance measurement method.

BACKGROUND

Recently, a system with which doors of an automobile can be locked and unlocked by pressing a button of a key-fob without inserting a key into a keyhole has been employed on a large number of automobiles. In the system, the doors are locked and unlocked by using wireless communication between the key-fob and the automobile.

In addition, recently, a system with which doors can be locked and unlocked and an engine can be started without touching a key-fob has been employed.

However, incidents in which an automobile is stolen by performing an attack called a “relay attack” on wireless communication between a key-fob and the automobile have been occurring. In a proposed countermeasure against the attack, a distance between the key-fob and the automobile is measured, and control of the automobile through wireless communication is prohibited when it is determined that the measured distance is equal to or longer than a predetermined distance.

The distance between the key-fob and the automobile is measured by, for example, a method in which a predetermined radio frequency (RF) signal (distance measurement signal) is transmitted from one of a distance measurement device mounted on the key-fob and a distance measurement device mounted on the automobile, the transmitted distance measurement signal is received by the other distance measurement device, a time difference between a transmission time point and a reception time point, in other words, a flight time of the distance measurement signal is measured or estimated, and the obtained flight time is multiplied by a light speed to obtain the distance.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a block diagram illustrating a configuration of a distance measurement device according to a first embodiment;

FIG. 2 is a graph illustrating a phase of a GMSK signal in the first embodiment;

FIG. 3 is a graph illustrating temporal change of the phase of the GMSK signal in the first embodiment;

FIG. 4 is a time chart illustrating an example in which a distance measurement signal makes one round trip between a device A and a device B in the first embodiment;

FIG. 5 is a time chart for description of error generated in a calculated flight time of the distance measurement signal as the distance measurement signal makes one round trip starting at the device A when a reference clock frequency of the device A has no error and a reference clock frequency of the device B has error in the first embodiment;

FIG. 6 is a time chart for description of error generated in a calculated flight time of the distance measurement signal as the distance measurement signal makes one round trip starting at the device B when the reference clock frequency of the device A has no error and the reference clock frequency of the device B has error in the first embodiment;

FIG. 7 is a time chart illustrating a sequence in which the distance measurement signal is transmitted in the following order of the device B to the device A and distance measurement is performed twice in the first embodiment;

FIG. 8 is a time chart illustrating an exemplary specific distance measurement sequence according to the first embodiment:

FIG. 9 is a time chart for description of exemplary decrease of an elapsed time in which the distance measurement signal is transmitted from the device A to the device B and detected at intervals of m sampling periods in a second embodiment;

FIG. 10 is a time chart for description of exemplary increase of an elapsed time in which the distance measurement signal is transmitted from the device B to the device A and detected at intervals of m sampling periods in the second embodiment;

FIG. 11 is timing charts for illustrating an example of averaging elapsed times obtained by performing phase detection K times at the device B and the device A in a sequence in which the distance measurement signal is transmitted from the device A to the device B and then the distance measurement signal is transmitted from the device B to the device A in the second embodiment;

FIG. 12 is a time chart illustrating an exemplary sequence in which phase detection is performed K times in accordance with a latter half part of FIG. 8 in comparison with the sequence of FIG. 8 in which the distance measurement signal is transmitted in the following order of the device A to the device B to the device B to the device A and phase detection is performed once in the second embodiment;

FIG. 13 is a time chart illustrating a 1.5 round trip sequence according to a third embodiment;

FIG. 14 is a time chart for description of a method of producing the distance measurement signal that reduces communication time in a distance measurement sequence in which the distance measurement signal is transmitted in the following order of the device A to the device B to the device B to the device A in a fourth embodiment;

FIG. 15 is a time chart illustrating a sequence employing a 1.5 round trip sequence and phase loop-back distance measurement sequence in a fifth embodiment:

FIG. 16 is a time chart illustrating an example in which an arrival time t_arvA of the distance measurement signal is between a sample time t₁ and a sample time t₂ in a sixth embodiment; and

FIG. 17 is a time chart for description of a method of calculating t_a by using four or more sample points in a seventh embodiment.

DETAILED DESCRIPTION

A distance measurement device of an embodiment includes a first device, a second device, and a calculation unit. The first device includes a first reference signal source configured to generate a first reference signal, and a first transmitter-receiver configured to transmit a modulated first distance measurement signal, receive a modulated second distance measurement signal, obtain a second demodulated signal through demodulation, and acquire second phase information of the second demodulated signal in a first sampling period based on the first reference signal. The second device includes a second reference signal source configured to operate independently from the first reference signal source and generate a second reference signal, and a second transmitter-receiver configured to transmit the second distance measurement signal, receive the first distance measurement signal, obtain a first demodulated signal through demodulation, and acquire first phase information of the first demodulated signal in a second sampling period based on the second reference signal. The calculation unit calculates a distance between the first device and the second device based on the first phase information and the second phase information. One of the first distance measurement signal and the second distance measurement signal is transmitted once or more, and another is transmitted twice or more. The calculation unit calculates the distance based on a total of three or more pieces of the first phase information and the second phase information acquired through transmission of the first distance measurement signal and the second distance measurement signal three times or more in total, the first sampling period, and the second sampling period.

First Embodiment

(Configuration)

FIG. 1 is a block diagram illustrating the configuration of a distance measurement device according to the present embodiment.

The distance measurement device includes a device A1 (hereinafter referred to as a device A) and a device B2 (hereinafter referred to as a device B) configured to perform wireless communication through a radio frequency (RF) signal. The device A and the device B are separated from each other by a certain distance, and the distance may time-sequentially change (at least one of the device A and the device B may be movable) but distance change in a time period in which distance measurement is performed once is negligible.

The device A includes a processor 11, a memory 12, a first reference signal source 13, and a first transmitter-receiver 14.

The device B includes a processor 21, a memory 22, a second reference signal source 23, and a second transmitter-receiver 24.

The first reference signal source 13 generates a first reference clock (first reference signal) and supplies the first reference clock to the processor 11, the memory 12, and the first transmitter-receiver 14.

The memory 12 stores a processing program to be executed by the processor 11, and various kinds of data used in processing at the processor 11.

The processor 11 controls transmission and reception operation of the first transmitter-receiver 14 in accordance with the processing program stored in the memory 12.

The first transmitter-receiver 14 generates a first sampling clock based on the first reference clock supplied from the first reference signal source 13. The first transmitter-receiver 14 transmits a modulated first distance measurement signal, receives a modulated second distance measurement signal transmitted from the second transmitter-receiver 24, and obtains a second demodulated signal through demodulation. The first transmitter-receiver 14 samples second phase information of the second demodulated signal at a timing of the first sampling clock.

In addition, the processor 11 calculates, based on a phase of a distance measurement signal sampled by the first transmitter-receiver 14, a flight time of the distance measurement signal in accordance with the processing program stored in the memory 12, and calculates a flight distance (distance between the device A and the device B) by multiplying the flight time by a light speed c.

The second reference signal source 23 operates independently from the first reference signal source 13, generates a second reference clock (second reference signal), and supplies the second reference clock to the processor 21, the memory 22, and the second transmitter-receiver 24.

The memory 22 stores a processing program to be executed by the processor 21, and various kinds of data used in processing at the processor 21.

The processor 21 controls transmission and reception operation of the second transmitter-receiver 24 in accordance with the processing program stored in the memory 22.

The second transmitter-receiver 24 generates a second sampling clock based on the second reference clock supplied from the second reference signal source 23. The second transmitter-receiver 24 transmits the second distance measurement signal, receives the first distance measurement signal, and obtains a first demodulated signal through demodulation. The second transmitter-receiver 24 samples first phase information of the first demodulated signal at a timing of the second sampling clock.

In addition, the processor 21 calculates, based on a phase of a distance measurement signal sampled by the second transmitter-receiver 24, the flight time of the distance measurement signal in accordance with the processing program stored in the memory 22, and calculates a flight distance (distance between the device A and the device B) by multiplying the flight time by the light speed c.

Note that the processor 11 and the processor 21 are a calculation unit configured to calculate the distance between the first device and the second device based on the first phase information and the second phase information. However, only one of the processor 11 and the processor 21 may be the calculation unit.

In the embodiment, a relatively narrow band signal using known modulation wave is used as each distance measurement signal. Specifically, a modulated signal is a gaussian filtered minimum-shift keying (GMSK) signal that repeats 0 and 1.

FIG. 2 is a graph illustrating a phase of a GMSK signal according to the present embodiment.

The first transmitter-receiver 14 and the second transmitter-receiver 24 each demodulate a received GMSK signal by separating the GMSK signal into an in-phase component (I signal) and an orthogonal component (Q signal). The phase of the GMSK signal can be calculated based on the I and Q signals.

As illustrated in FIG. 2, on a phase plane, the GMSK signal, which repeats 0 and 1, does not transition between (I, Q)=(1, 0) and (I, Q)=(0, 1) but transitions between (I, Q)=(a, b) and (I, Q)=(b, a) due to influence of a Gauss filter.

FIG. 3 is a graph illustrating temporal change of the phase of the GMSK signal according to the present embodiment. In FIG. 3, a horizontal axis represents time, and a vertical axis represents the phase with respect to an origin.

The phase θ_GMSK of the GMSK signal can be approximated by Expression (1).

θ_GMSK≈A₀₁ sin(2πf₀₁t)+θ_OFF  (1)

In the expression, A₀₁ represents an amplitude of the phase of the GMSK signal, which repeats 0 and 1 (0 for a case in which the phase decreases, and 1 for a case in which the phase increases), f₀₁ represents a basic frequency (reciprocal of a basic period T₀₁) of the phase of the GMSK signal, which repeats 0 and 1, in other words, half of a transmission frequency of the GMSK signal, and θ_OFF represents an offset phase of the phase of the GMSK signal, which repeats 0 and 1, in other words, π/4. When a time is to be measured, for example, a time in which θ_GMSK indicated in Expression (1) becomes θ_OFF may be detected.

When frequencies (sampling frequencies) of the sampling clocks of the device A and the device B are set to be high, time resolution is high, and thus the time in which the phase θ_GMSK becomes θ_OFF is accurate. However, with taken into consideration low electric power consumption, it is desired that each sampling frequency is, for example, 16 times approximately higher than the basic frequency of the phase signal θ_GMSK. To clearly determine in which of four quadrants a phase of the sine wave exists, each sampling clock needs a frequency four times or more higher than the basic frequency of the phase signal θ_GMSK.

Under such a condition, accuracy of measurement of the flight time can be prevented from significantly decreasing by sampling the phase θ_GMSK based on the sampling clock of each of Devices A and B and applying a signal that can be approximated to the sine wave indicated in Expression (1). Note that, in the following description, a “signal that can be approximated to a sine wave” is simply referred to as a sinusoidal signal as long as no interpretation problem occurs. In addition, it is assumed that a sinusoidal signal such as the GMSK signal indicated in Expression (1) is used as a distance measurement signal.

FIG. 4 is a time chart illustrating an example in which a distance measurement signal makes one round trip between the device A and the device B in the present embodiment. FIG. 4 illustrates a case in which the first reference clock of the device A and the second reference clock of the device B have no error in frequency but have a reference time difference T_OFF.

FIG. 4 illustrates a time chart of the device A at an upper part and a time chart of the device B at a lower part. Since reference clock frequency is same between the device A and the device B, a sampling period T_S is same. However, a sampling timing of the device B is shifted temporally later from a sampling timing of the device A by the reference time difference T_OFF.

A sinusoidal signal (distance measurement signal) starting at the phase θ_OFF is transmitted at the sampling timing of the device A, and the distance measurement signal is received by the device B after a flight time t_x in accordance with the distance between Devices A and B has elapsed. In the time chart of FIG. 4 and later, a transmitted wave is illustrated with a dashed line, and received wave is illustrated with a solid line. In addition, phases θ_GMSK,b and θ_GMSK,a and the like sampled at the sampling clocks are illustrated with black circles.

A time point at which the phase θ_GMSK,b at the device B is sampled is a time point at which a time period to has elapsed since the distance measurement signal starting at the phase θ_OFF arrives at the device B. After a predetermined process time period T_pres in which various kinds of calculation and the like are performed has elapsed since reception of the distance measurement signal is detected, the device B transmits a distance measurement signal at the sampling timing of the device B.

The distance measurement signal transmitted from the device B is received by the device A after elapse of the flight time t_x. A time point at which the phase θ_GMSK,a is sampled at the device A is a time point at which a time period t_a has elapsed since the sinusoidal signal (distance measurement signal) starting at the phase θ_OFF arrives at the device A.

The elapsed time t_a is calculated by Expression (2) from the phase θ_GMSK,a sampled by the device A, and the elapsed time to is calculated by Expression (3) from the phase θ_GMSK,b sampled by the device B.

t _a≈{1/(2πf₀₁}sin⁻¹{(θ_GMSK,a˜θ_OFF)/A₀₁}  (2)

t _b≈{1/(2πf₀₁)}sin⁻¹{(θ_GMSK,b˜θ_OFF)/A₀₁)}  (3)

Note that, in FIG. 4, phases in which Devices A and B start transmission at the timings of the sampling clocks are θ_OFF. In each of Expressions (2) and (3), the elapsed times t_a and t_b are uniquely determined when the sampling frequency is four times or more higher than the frequency f₀₁ of the sinusoidal signal.

A round trip time RTT elapses until the device A receives a signal from the device B after the device A transmits a signal. As understood from FIG. 4, the round trip time RTT is expressed by Expression (4) as the sum of the flight time t_x from the device A to the device B, the elapsed time t_b until a signal is sampled since arrival at the device B, the process time period T_pres at the device B, the flight time t_x from the device B to the device A, and the elapsed time t_a until a signal is sampled since arrival at the device A.

RTT=2t_x+t_a+t_b+T_pres  (4)

In the expression, a difference between the round trip time RTT and the process time period T_pres is equal to the sampling period T_S, and thus Expression (5) is obtained.

T _S=2t_x+t_a+t_b  (5)

Expression (5) can be rewritten as Expression (6).

t _x=(T_S˜t_a˜t_b)/2  (6)

Thus, the flight time t_x can be accurately obtained through Expression (6) by calculating the elapsed time t_a through Expression (2), calculating the elapsed time t_b through Expression (3), and using the known sampling period T_S.

However, when there is error between the period (reciprocal of the frequency) of the first reference clock of the device A and the period of the second reference clock of the device B, error occurs to the flight time t_x calculated through Expression (6). The error occurs by a reason described below.

For example, assume that the reference clock period T_A of the device A is accurate and has no error and the reference clock period T_B of the device B has error. A ratio T_B/T_A of the reference clock period T_B of the device B relative to the reference clock period T_A of the device A is defined by Expression (7).

T _B /T _A=1˜r_err  (7)

In the expression, r_err represents an error ratio of the reference clock period T_B. A sample signal is generated by dividing or multiplying a reference clock signal, and each sampling frequency is proportional to the reference clock frequency. Thus, a ratio T_SB/T_S of the sampling period T_SB of the device B relative to the sampling period T_S of the device A is expressed by Expression (8) similar to Expression (7).

T _SB /T _S=1˜r_err  (8)

Note that the sampling period of the device A is not written as TSA because no error is assumed.

FIG. 5 is a time chart for description of error that occurs to the calculated flight time t_x of the distance measurement signal when the reference clock frequency of the device A has no error and the reference clock frequency of the device B has error in the present embodiment. Note that error between the reference clock frequencies of Devices A and B means that the reference clock period of the device B has error.

FIG. 5 illustrates a time chart of the device A at an upper part, illustrates a time chart of the device B when the reference clock frequency has no error at a middle part, and illustrates a time chart of the device B when the reference clock frequency has error at a lower part.

When a time period T_presB taken for a process at the device B is N times longer than the sampling period of the device B, Expression (9) is obtained.

T _presB =NT _SB  (9)

A numerator and a denominator on a left-hand side of Expression (8) are each multiplied by N to obtain Expression (10).

NT _SB /NT _S=1˜r_err  (10)

Both sides of Expression (10) are multiplied by NT_S to obtain Expression (11).

NT _SB =NT _S ˜r _err NT _S  (11)

Since NT_S is the process time period T_pres when the clock period has no error, Expression (12) is obtained by substituting Expression (9) into a left-hand side of Expression (11).

T _presB =T _pres ˜r _err T _pres  (12)

Expression (12) indicates that the process time period T_presB of the device B is shorter than the process time period T_pres when there is no error by r_errT_pres.

In FIG. 5, a time point at which the device B illustrated at the lower part starts transmission to the device A is earlier by r_errT_pres than a time point at which the device B illustrated at the middle part starts transmission to the device A. As a result, a reception phase θ′_GMSK,a detected at the sampling clock of the device A advances by an amount by which the time point of transmission from the device B is earlier.

Thus, an elapsed time t′_a since signal arrival at the device A, which is detected at the sampling timing of the device A, when the clock period of the device B has error is longer by r_errT_pres than the elapsed time t_a since signal arrival at the device A, which is detected at the sampling timing of the device A, when the clock period of the device B has no error, and accordingly, Expression (13) holds.

t′ _a =t _a +r _err T _pres  (13)

When t′_a to which error occurs is substituted into t_a in Expression (6) to calculate a flight time t′_x, Expression (14) is obtained.

t′ _x=(T_S˜t_a˜t_b˜r_errT_pres)/2  (14)

As understood from comparison of Expression (14) with Expression (6), error due to r_errT_pres occurs to the flight time t′_x. For example, when the process time period T_pres is 100 [μs] and r_err is 40 [ppm], r_errT_pres/2 is 2 [ns]. Conversion into distance by multiplying 2 [ns] by the light speed indicates that a distance error of 0.6 [m] occurs. The distance error further increases as the process time period and the clock period error further increase.

FIG. 5 illustrates that a distance measurement signal makes one round trip starting at the device A with a setting [Expression (7)] that the reference clock period of the device B is shorter by the error ratio r_err than the reference clock period of the device A, which is accurate, and as a result, the flight time is estimated to be shorter by r_errT_pres/2 as indicated in Expression (14).

Subsequently, discussion is made on error that occurs to the flight time when a distance measurement signal makes one round trip starting at the device B with the setting of Expression (7).

FIG. 6 is a time chart for description of error that occurs to the calculated flight time of a distance measurement signal as the distance measurement signal makes one round trip starting at the device B when the reference clock frequency of the device A has no error and the reference clock frequency of the device B has error in the present embodiment.

In FIG. 6, a time chart of the device A is illustrated at an upper part, a time chart of the device B when the reference clock frequency has no error is illustrated at a middle part, and a time chart of the device B when the reference clock frequency has error (the reference clock frequency is unique to the device B) is illustrated at a lower part. Note that FIG. 6 illustrates a case in which a time point at which the distance measurement signal is transmitted from the device B illustrated at the middle part to the device A and a time point at which the distance measurement signal is transmitted from the device B illustrated at the lower part to the device A are same.

When the distance measurement signal makes one round trip in a reverse order, as well, the round trip time RTT is expressed by Expression (4) and the flight time t_x is expressed by Expression (6) as long as the reference clock periods of Devices A and B are both accurate.

However, when the distance measurement signal makes one round trip between the device B illustrated at the lower part and the device A illustrated at the upper part in FIG. 6, error due to the error ratio r_err occurs to an estimated flight time (that is, flight distance) as described below.

RTT_B represents the round trip time RTT detected based on an internal clock of the device B and is expressed by Expression (15).

RTT _B=2t_x+t_a+t′_b+T_pres  (15)

In the expression, t′_b is an elapsed time calculated by Expression (3) based on a detection phase θ′_GMSK,b of a GMSK signal received by the device B. Since RTT_B=(N+1)T_SB, Expression (16) is obtained by using Expression (8) and T_pres=NT_S.

RTT _B=(N+1)T_S˜r_err(N+1)T_S=T_pres+T_S˜r_err(N+1)T_S  (16)

Comparison of right hand sides of Expressions (15) and (16) obtains Expression (17).

T _S=2t_x+t_a+t′_b+r_err(N+1)T_S  (17)

Comparison of Expressions (5) and (17) obtains Expression (18).

t′ _b =t _b ˜r _err(N+1)T_S  (18)

Expression (18) means that a timing at which the distance measurement signal received from the device A is sampled by the device B at the lower part when the round trip time RTT_B has elapsed is earlier by r_err(N+1)T_S than a timing at which the distance measurement signal received from the device A is sampled by the device B at the middle part.

When t′_b to which error has occurred is substituted into Expression (6) to calculate a flight time t″_x and T_pres=NT_S is used, Expression (19) is obtained.

t″ _x=(T_S˜t_a˜t_b+r_errT_pres+r_errT_S)/2  (19)

Comparison between Expression (14), which is obtained as the distance measurement signal makes one round trip starting at the device A, and Expression (19), which is obtained as the distance measurement signal makes one round trip starting at the device B indicates that the error time period r_errT_pres has opposite signs between the expressions. In addition, the other error time period r_errT_S in Expression (19) is clearly small as compared to r_errT_pres. Thus, the terms r_errT_pres are eliminated by averaging the flight time t′_x estimated by Expression (14) and the flight time t″_x estimated by Expression (19), and the accuracy of the estimated time clearly improves.

For example, r_errT_pres/4=1.25 [μs] is obtained for r_err=40 [ppm] and T_S=125 [ns]. Conversion into distance by multiplying 1.25 [ps] by the light speed obtains a distance error of 0.4 [mm], which indicates significant improvement of measurement accuracy as compared to the distance error of 0.6 [m] when the distance measurement signal makes one round trip starting at the device A.

Thus, significant accuracy improvement of an estimated distance is achieved when the distance measurement signal is transmitted in the following order of the device A to the device B and distance measurement is performed, the distance measurement signal is transmitted in the following order of the device B to the device A and distance measurement is performed again, and results of the two distance measurements are averaged.

However, in FIG. 6, it is assumed that, when the distance measurement signal is transmitted by the device B, an elapsed time until the distance measurement signal is sampled since arrival at the device A is t_a, and the above-described calculation is performed based on the assumption.

The elapsed time t_a illustrated in FIG. 5 is a time that elapses until sampling since an arrival time when the distance measurement signal from the device B is sampled by the device A after a time period (N+1)T_S has elapsed since the distance measurement signal is transmitted by the device A.

Thus, the elapsed time t_a illustrated in FIG. 6 is different from the elapsed time t_a illustrated in FIG. 5. Despite of the difference between the elapsed times, the flight time of Expression (19) can be obtained as described below.

Specifically, it is described that, even when the distance measurement signal is transmitted by the device B with delay of a predetermined number of clocks, for example, L clocks, the flight time expressed by Expression (19) when the distance measurement signal is transmitted without delay holds, and the estimated flight distance is same irrespective of a time point at which the distance measurement signal is transmitted.

FIG. 7 is a time chart illustrating a sequence in which the distance measurement signal is transmitted in the following order of the device B to the device A and distance measurement is performed twice in the present embodiment. In the sequence illustrated in FIG. 7, the distance measurement signal is transmitted in the following order of the device B to the device A, and first-time distance measurement is performed. Then, after a time period LT_SB has elapsed since the first-time distance measurement signal is transmitted by the device B, the distance measurement signal is transmitted in the following order of the device B to the device A, and second-time distance measurement is performed.

A distance measurement result obtained by the first-time distance measurement is same as Expression (19).

In the second-time distance measurement, the distance measurement signal is transmitted from the device B after the time period LT_SB has elapsed since the first-time distance measurement signal is transmitted by the device B. A time at which the transmitted distance measurement signal is received by the device A is represented by t_aL. The time t_aL is a time calculated by substituting a phase θ_GMSK,aL observed (sampled) by the device A into Expression (2).

When the process time period T_pres=NT_S has elapsed since the phase θ_GMSK,aL of the distance measurement signal is detected, the distance measurement signal is transmitted by the device A. The device B observes (samples) a phase θ′_GMSK,bL of the distance measurement signal from the device A and substitutes the phase into Expression (3) to calculate an elapsed time t′_bL.

The time period LT_SB that has elapsed until the distance measurement signal is transmitted by the device B in the second-time distance measurement since the distance measurement signal is transmitted by the device B in the first-time distance measurement is expressed by Expression (20) by using Expression (8).

LT _SB =LT _S ˜r _err LT _S  (20)

A transmission time point of the distance measurement signal of the device B is earlier by r_errLT_S with respect to the device A. Thus, a time at which the distance measurement signal is received by the device A is earlier by the same amount, and accordingly, t_aL is given by Expression (21).

t _aL =t _a +r _err LT _S  (21)

However, since the round trip time RTT_B=(N+1)T_SB is constant, the round trip time in the first-time distance measurement and the round trip time in the second-time distance measurement are same. The round trip time RTT_B in the second-time distance measurement is given by Expression (22).

RTT _B=2t_x+t_aL+t′_bL+T_pres  (22)

Since Expression (22) is same as Expression (15), t′_bL is given by Expression (23).

t′ _bL =t _b ˜r _err(N+1)T_S˜r_errLT_S  (23)

The flight time t_x is calculated by substituting t_aL into t_a and substituting t′_bL into t_b in Expression (6) while assuming that the device A and the device B operate at accurate clocks. Accordingly. Expression (24) is obtained.

t _x=(T_S˜t_aL˜t′_bL=(T_S˜t_a˜t_b+r_errT_pres+r_errT_S)/2  (24)

Since Expressions (19) and (24) are same, a distance measurement result when the distance measurement signal is transmitted in the following order of the device B to the device A is same irrespective of a distance measurement start time. Thus, after the distance measurement signal is transmitted in the following order of the device A to the device B and the first-time distance measurement is performed, the distance measurement signal is transmitted in the following order of the device B to the device A and the second-time distance measurement is performed, and results of the two distance measurements are averaged to obtain the flight time and hence the flight distance at improved accuracy.

FIG. 8 is a time chart illustrating an exemplary specific distance measurement sequence according to the present embodiment. In the distance measurement sequence illustrated in FIG. 8, the distance measurement signal is transmitted in the following order of the device A to the device B to the device B to the device A.

In FIG. 8, since the numbers of times of transmission and reception of the distance measurement signal by the device A and the device B are two, respectively, elapsed times from the first reception of the distance measurement signal to sampling by the device B and the device A are represented by t_b1 and t_a1, respectively, and elapsed times from the second reception of the distance measurement signal to sampling by the device B and the device A are represented by t_b2 and t_a2, respectively.

When a flight time estimated from t_b1 and t_a1 is represented by t_x1 and a flight time estimated from t_b2 and t_a2 is represented by t_x2, Expressions (25) and (26) are obtained from Expressions (14) and (24), respectively.

t _x1=(T_S˜t_a1˜t_b1)/2=(T_S˜t_a˜t_b˜r_errT_pres)/2  (25)

t _x2=(T_S˜t_a2˜t_b2)/2=(T_S˜t_a˜t_b+r_errT_pres+r_errT_S)/2  (26)

Thus, the flight time t_x obtained by averaging t_x1 and t_x2 is given by Expression (27).

t _x=(t_x1+t_x2)/2={T_S˜t_a˜t_b+(r_errT_S)/2}/2  (27)

According to the first embodiment, round trip sequences with different starting points, namely, a distance measurement sequence in which the distance measurement signal makes a round trip starting at the device A and a distance measurement sequence in which the distance measurement signal makes a round trip starting at the device B are performed, and flight times obtained by the respective distance measurement sequences are averaged. Accordingly, when the reference clock periods of the device A and the device B have error, error of an estimated flight time can be significantly reduced, and thus, error of the flight distance obtained by multiplying the flight time by the light speed can be significantly reduced.

Specifically, in the present embodiment, a time offset (T_OFF illustrated in FIG. 6) can be corrected through one round trip of the distance measurement signal, and a frequency offset can be corrected by adding one round trip of the distance measurement signal with a different starting point and performing averaging. When a reference clock frequency corresponding to the reference clock period T_A of the device A is represented by f_A and a reference clock frequency corresponding to the reference clock period T_B of the device B is represented by f_B, Expression (7) becomes Expression (46) to be described later, and a corrected frequency offset is r_errf_B obtained from Expression (46).

Note that although flight times are averaged in the above description, flight distances may be calculated from the flight times and then averaging may be performed, or sampled phases may be averaged and then a flight time and a flight distance may be calculated (refer to Expression (72) in a seventh embodiment, for example).

Second Embodiment

In a second embodiment, any part same as a part in the first embodiment is denoted by the same reference sign, description of the part is omitted, and difference will be mainly described.

In the above description, elapsed times until phase signals are sampled since crossing the phase offset θ_OFF corresponding to a zero cross point are measured by Devices A and B. In description of the first embodiment, the phase θ_GMSK,a or θ_GMSK,b at a sample point right after reception is started is measured, and one phase data is targeted.

However, when one reception duration of the distance measurement signal is long as compared to a period 1/f₀₁ of the distance measurement signal, and for example, the distance measurement signal of K periods (K is a natural number) can be received, an elapsed time until sampling since crossing of a plurality of phase offsets θ_OFF can be measured in one reception duration of the distance measurement signal. In this case, measurement variance of phase detection due to noise and the like can be prevented by detecting the phase a plurality of times and averaging the phases, and thus the flight time t_x can be accurately estimated.

However, it is needed to prove that, when there is a frequency offset between the reference clocks of Devices A and B, the flight time t_x calculated by one phase detection is equal to the flight time t_x calculated by detecting the phase a plurality of times and performing averaging. The following proves that these flight times are equal.

Assume that, when there is no frequency error, the period 1/f₀₁ of the distance measurement signal is m times longer than the sampling period T_S, and a duration in which the distance measurement signal is received is K/f₀₁.

First, when the distance measurement signal is transmitted from the device A to the device B, a sum t_sumB of detection elapsed times sampled K times (elapsed times until sampling since crossing of θ_OFF) is calculated with reference to FIG. 9. FIG. 9 is a time chart for description of exemplary decrease of an elapsed time in which the distance measurement signal is transmitted from the device A to the device B and detected at intervals of m sampling periods in the present embodiment. In other words, the sampling period of the device B is shorter than the sampling period of the device A.

A time at which the distance measurement signal is sampled by the device B for the first time is referred to as a time t_b1 based on a detection phase θ_GMSK,b1. An elapsed time until next phase detection since phase detection for the first time is represented by mT_SB.

A detection time for the second time is earlier than a correct time by r_errmT_S by using Expression (8) because mT_SB has elapsed since a detection time for the first time, and thus an observed time is t_b1˜r_errmT_S.

Similarly, by analogy, a k-th detection time is t_b1˜(k˜1)r_errmT_S. Thus, the sum t_sumB of the first to K-th detection times is given by Expression (28).

t _sumB=Σ₁^K{t_b1˜(k˜1)r_errmT_S}=K(t_b1+r_errmT_S)˜{K(K+1)/2}r_errmT_S  (28)

Thus, an average of t_sumB when detection is performed K times is given by Expression (29).

t _sumB /K=t _b1˜{(K˜1)/2}r_errmT_S  (29)

When the distance measurement signal is transmitted from the device A to the device B, an average of elapsed times detected K times at intervals of m sampling periods is shorter than an elapsed time detected only once on the first time by an elapsed time {(K˜1)/2}r_errmT_S.

Subsequently, an elapsed time based on phases detected K times by the device A when the distance measurement signal is transmitted from the device B to the device A after the time period NT_SB has elapsed is calculated with reference to FIG. 10. FIG. 10 is a time chart for description of exemplary increase of an elapsed time in which the distance measurement signal is transmitted from the device B to the device A and detected at intervals of m sampling periods according to the present embodiment.

An elapsed time calculated based on initial phase detection is represented by t_a1. A time difference between first phase detection and second phase detection is mT_S, time advances by r_errmT_S at the device B in mT_S, and the distance measurement signal is transmitted earlier accordingly. Thus, an elapsed time calculated based on the second phase detection is t_a1+r_errmT_S.

Similarly, by analogy, a k-th detection elapsed time is t_a1+(k˜1)r_errmT_S. Thus, a sum t_sumA of the first to K detection elapsed times is given by Expression (30).

t _sumA=Σ₁^K{t_a1+(k˜1)r_errmT_S}=K(t_a1˜r_errmT_S)+{K(K+1)/2}r_errmT_S  (30)

Thus, an average of t_sumA when detection is performed K times is given by Expression (31).

t _sumA /K=t _a1+{(K˜1)/2}r_errmT_S  (31)

When the distance measurement signal is transmitted from the device B to the device A, an average of elapsed times detected K times at intervals of m sampling periods is longer than an elapsed time detected only once on the first time by the elapsed time {(K˜1)/2}r_errmT_S.

However, when t_sumA/K and t_sumB/K are used in place of t_a and t_b in Expression (6), a flight time t_{x1_avg} of an initial round trip is estimated as in Expression (32).

t _{x1_avg} ={T _S˜(t_sumA/K)˜(t_sumB/K)}/2=(T_S˜t_a1˜t_b1)/2=(T_S˜t_a˜t_b˜r_errT_pres)/2=t_x1  (32)

This is same as a flight time without averaging.

FIG. 11 is timing charts for illustrating an example of averaging of elapsed times obtained by performing phase detection K times at the device B and the device A in a sequence in which the distance measurement signal is transmitted from the device A to the device B and then the distance measurement signal is transmitted from the device B to the device A in the present embodiment.

In the first embodiment, when a flight time is calculated through one phase detection, error due to a frequency difference included in the estimated flight time t_x1 is reduced by averaging elapsed times calculated by using a sequence in which the distance measurement signal is transmitted from the device A to the device B, from the device B to the device A, from the device B to the device A, and from the device A to the device B as illustrated in FIG. 8.

Similarly, in the present embodiment, error of flight time estimation is reduced by using the sequence illustrated in FIG. 8 when elapsed times obtained from phases detected K times at intervals of m sampling periods are averaged.

FIG. 12 is a time chart illustrating an exemplary sequence in which phase detection is performed K times in accordance with a latter half part of FIG. 8 in comparison with the sequence of FIG. 8 in which the distance measurement signal is transmitted in the following order of the device A to the device B to the device B to the device A and phase detection is performed once per transmission in the present embodiment.

A sequence on first and second rows in FIG. 12 is same as the sequence in FIG. 8. A sequence on third and fourth rows in FIG. 12 illustrates an image of a sequence corresponding to the latter half part in FIG. 8, in which the distance measurement signal is transmitted from the device B to the device A and from the device A to the device B and phase detection is performed K times. Note that an image of a sequence corresponding to a first half part of FIG. 8, in which the distance measurement signal is transmitted from the device A to the device B and from the device B to the device A respectively and phase detection is performed K times is illustrated in FIG. 11 and thus omitted in FIG. 12.

As illustrated on the third and fourth rows in FIG. 12, when the distance measurement signal is transmitted from the device B to the device A and the distance measurement signal is received by the device A, a first detection phase is θ_GMSK,a2, and an elapsed time corresponding to the detection phase is t_a2.

Thereafter, when phases are detected K times, a sum t_sumA2 of elapsed times calculated from the detected phases is equal to an expression obtained by replacing t_a1 with t_a2 in Expression (30). Thus, an average t_sumA2/K of t_sumA2 is given by Expression (33).

t _sumA2 /K=t _a2+{(K˜1)/2}r_errmT_S  (33)

Subsequently, when the distance measurement signal is transmitted from the device A to the device B and the distance measurement signal is received by the device B, a first detection phase is θ_GMSK,b2, and an elapsed time corresponding to the detection phase is t_b2.

Thereafter, when phases are detected K times, a sum t_sumB2 of elapsed times calculated from the detected phases is equal to an expression obtained by replacing t_b1 with t_b2 in Expression (29). Thus, an average t_sumB2/K of t_sumB2 is given by Expression (34).

t _sumB2 /K=t _b2˜{(K˜1)/2}r_errmT_S  (34)

When t_sumA2/K and t_sumB2/K are used in place of t_a and t_b in Expression (6), a flight time of a second round trip is estimated as in Expression (35), which is same as an estimated flight time without averaging.

t _{x2_avg} ={T _S˜(t_sumA2/K)˜(t_sumB2/K)}/2=(T_S˜t_a˜t_b2)/2=(T_S˜t_a˜t_b+r_errT_pres+r_errT_S)/2=t_x2  (35)

A flight time t_{x_avg} is calculated by averaging t_{x1_avg} and t_{x2_avg} as in Expression (36).

t _{x_avg}=(t_{x1_sum}+t_{x2_sum})/2={T_S˜t_a˜t_b+(r_errT_S/2)}/2=t_x  (36)

A flight time obtained by performing phase detection K times in a sequence in which the distance measurement signal is transmitted from the device A to the device B, from the device B to the device A, from the device B to the device A, and from the device A to the device B in the stated order and averaging elapsed times obtained from detected phases is same as a flight time obtained by performing phase detection once. However, the same result is obtained by an expectation value, and a result of the averaging has reduced error attributable to measurement variance.

According to the second embodiment, effects substantially same as effects of the first embodiment described above are achieved, and in addition, distance measurement can be performed at higher accuracy by increasing the amount of data sampled through the distance measurement sequence and averaging the obtained data.

Third Embodiment

FIG. 13 is a time chart illustrating a 1.5 round trip sequence according to the present embodiment. In a third embodiment, any part same as a part in the first and second embodiments is denoted by the same reference sign, description of the part is omitted, and difference will be mainly described.

In the distance measurement sequence illustrated in FIG. 8, a time period from transmission of the distance measurement signal by the device A or the device B once to transmission of the next distance measurement signal by the device A or the device B corresponds to N sampling clocks, and thus a time period corresponding to 3N clocks is needed when the distance measurement signal is transmitted in the following order of the device A to the device B to the device B to the device A. The time period is 3NT_S when converted with the reference clock of the device A. The present embodiment employs a sequence (hereinafter referred to as 1.5 round trip sequence) illustrated in FIG. 13, thereby reducing a distance measurement time period from 3NT_S to 2NT_S approximately.

In the sequence in FIG. 8, phase measurement is performed twice by the device A to obtain the elapsed times t_a1 and t_a2. In the 1.5 round trip sequence of the present embodiment, as illustrated in FIG. 13, phase measurement is performed once by the device A, and the distance measurement signal is transmitted in the following order of the device A to the device B to the device A. The 1.5 round trip sequence corresponds to a case in which L=N is set in the sequence in FIG. 8, and thus t_a1=t_a2 is obtained. This corresponds to a time shift of one round trip sequence (the device B to the device A) in the latter half by NT_S.

As indicated by comparison between Expressions (19) and (24), one-round-trip distance measurement sequence of the device B to the device A obtains a unique flight time t_x2 irrespective of a starting point of time. Thus, an elapsed time observed in the sequence of the device B to the device A illustrated in FIG. 13 is t_x2.

Note that although not illustrated, it is clear that the method of performing phase detection K times for each m samples and averaging elapsed times obtained from detected phases according to the second embodiment is also applicable to the present embodiment.

According to the third embodiment, it is possible to reduce the distance measurement time period by using the 1.5-round-trip distance measurement sequence in effect in place of the two-round-trip distance measurement sequence in the first and second embodiments described above. In a specific example, a time period reduced by the sequence in FIG. 13 is 33% approximately of a time period by the sequence in FIG. 8.

Fourth Embodiment

In a fourth embodiment, any part same as a part in the first to third embodiments is denoted by the same reference sign, description of the part is omitted, and difference will be mainly described.

Calculation of the flight time t_x needs data of the elapsed times t_a1 and t_a2, t_b1, and t_b2 until sampling since crossing of the phase offset θ_OFF, the data being obtained through phase detection by the device A and the device B.

Data of t_b1 and t_b2 needs to be sent from the device B to the device A to calculate the flight time t_x at the device A, and data of t_a1 and t_a2 needs to be sent from the device A to the device B to calculate the flight time t_x at the device B.

This communication of distance measurement data is performed by using typical digital communication after the distance measurement data is detected. A time taken for communication of distance measurement data depends on requested accuracy, transmission rate, and the like of the distance measurement data and is, for example, a time period of several hundreds [μs] approximately in some cases. Thus, when communication of distance measurement data is unnecessary, a time taken for a distance measurement sequence is reduced. The present embodiment relates to a technology of producing the distance measurement signal that reduces or eliminates the time of communication of distance measurement data.

FIG. 14 is a time chart for description of a method of producing the distance measurement signal that reduces communication time in a distance measurement sequence in which transmission of the distance measurement signal is performed in the following order of the device A to the device B to the device B to the device A in the present embodiment.

FIG. 14 is different from FIG. 8 of the first embodiment as follows. When the distance measurement signal is transmitted for a first time from the device B to the device A, the phase at the device B at the timing NT_SB is the phase θ_GMSK,b1, at which the distance measurement signal is received from the device A, instead of the offset phase θ_OFF. In addition, when the distance measurement signal is transmitted for a second time from the device A to the device B, the phase at the timing NT_S of the device A is the phase θ_GMSK,a2, at which the distance measurement signal is received from the device B, instead of the offset phase θ_OFF. Difference of these processes is that loop-back transmission is performed at a reception phase.

Note that the distance measurement time period is reduced by performing at least one of transmission of the distance measurement signal of the phase θ_GMSK,b1 by the device B at the timing NT_SB and transmission of the distance measurement signal of the phase θ_GMSK,a2 by the device A at the timing NT_S.

The following describes a reason why the distance measurement time period is reduced by transmitting the distance measurement signal of the phase θ_GMSK,b1 from the device B at the timing NT_SB. Assume that distance measurement calculation is performed by the device A.

When the distance measurement signal of the phase θ_GMSK,b1 is transmitted by the device B at the timing NT_SB, a timing of the offset phase θ_OFF is earlier than the timing NT_SB by t_b1. When a time point of the offset phase θ_OFF is a time point of transmission start of the distance measurement signal, the time point of transmission start is NT_SB˜t_b1.

Subsequently, a case in which the distance measurement signal is received by the device A will be described. First, when the device B starts transmitting the distance measurement signal of the phase θ_OFF at the timing NT_SB, the device A samples a phase θ_GMSK,a1 at the sampling timing and converts the phase θ_GMSK,a1 into the detection elapsed time t_a1.

On the other hand, when the device B starts transmitting the distance measurement signal of the phase θ_OFF at the timing “NT_SB˜t_b1”, the device A determines that the distance measurement signal arrives earlier than the elapsed time t_a1 of the phase θ_OFF by the elapsed time t_b1. Accordingly, the phase sampled by the device A at the sampling timing is not the phase θ_GMSK,a1 but the phase θ_GMSK,ab1 corresponding to the reception earlier by the elapsed time t_b1. Thus, a relation of Expression (37) holds.

t _a1 +t _b1≈{1/(2πf₀₁)}sin⁻¹{(θ_GMSK,ab1˜θ_OFF)/A₀₁}  (37)

Thus, data of the elapsed time t_b1 detected by the device B can be acquired at the device A by detecting the phase θ_GMSK,ab1. When the distance measurement signal of the phase θ_OFF is not transmitted from the device B at the timing “NT_SB˜t_b1”, the device B needs to transmit data of the detected elapsed time t_b1 to the device A.

Similarly, when the distance measurement signal of the phase θ_GMSK,a2 is transmitted from the device A at the timing NT_S, the device B detects a phase θ_GMSK,ab2, and this phase and an elapsed time “t_a2+t_b2” to be detected through a second round trip sequence have a relation of Expression (38).

t _a2 +t _b2≈{1/(2πf₀₁)}sin⁻¹{(θ_GMSK,ab2˜θ_OFF)/A₀₁}  (38)

However, this result needs to be transmitted to the device A when distance calculation is to be performed by the device A, and thus the phase of the distance measurement signal to be transmitted does not need to be shifted in a second round trip sequence. The result leads to a reduced distance measurement time period when distance calculation is performed by the device B. A result detected by Expression (28) needs to be transmitted to the device B when distance calculation is to be performed by the device B. Thus, phase shift of the distance measurement signal may be performed in accordance with whether to perform distance calculation at the device A or the device B.

For verification, improvement of distance measurement accuracy by using Expressions (37) and (38) will be described below. The averaged flight time t_x can be rewritten as Expression (39) by using Expressions (25) to (27).

t _x=(t_x1+t_x2)/2={2T_S˜(t_a1˜t_b1)˜(t_a2+t_b2)}/4  (39)

The term (t_a1+t_b1) in Expression (39) can be calculated by Expression (37), and the term (t_a2+t_b2) in Expression (39) can be calculated by Expression (38). Thus, any one of Expressions (37) and (38) can be used depending on whether to perform distance calculation at the device A or the device B. and the amount of data communication in a detection time can be reduced.

Note that although not illustrated, it is clear that the method of performing phase detection K times for each m samples and averaging elapsed times obtained from detected phases in the second embodiment is also applicable to the present embodiment.

According to the fourth embodiment, effects substantially same as effects of the first to third embodiments described above are achieved, and in addition, a data communication time is reduced by transmitting a loop-back distance measurement signal at the phase of a received distance measurement signal, and the distance measurement time period can be further reduced.

Fifth Embodiment

In a fifth embodiment, any part same as a part in the first to fourth embodiments is denoted by the same reference sign, description of the part is omitted, and difference will be mainly described.

In the present embodiment, the distance measurement time period is further reduced by employing 1.5-round-trip distance measurement in the third embodiment and loop-back distance measurement at a received phase in the fourth embodiment. FIG. 15 is a time chart illustrating a sequence employing the 1.5 round trip sequence and the phase loop-back distance measurement sequence in the present embodiment. In the sequence in FIG. 15, the distance measurement signal is transmitted in the following order of the device A to the device B to the device A.

Difference between the sequence in FIG. 15 and the 1.5 round trip sequence in FIG. 13 is the phase of the distance measurement signal transmitted last by the device A. In FIG. 13, the distance measurement signal of the offset phase GOFF is transmitted at the timing NT_S of the device A.

In contrast, in FIG. 15, a phase transmitted at the timing NT_S is a phase θ_GMSK,TXa corresponding to an arrival time twice longer than the arrival time tai estimated from the phase θ_GMSK,a1 received by the device A, in other words, θ_GMSK,TXa satisfying Expression (40).

2t_a1≈{1/(2πf₀₁)}sin⁻¹{(θ_GMSK,TXa˜θ_OFF)/A₀₁}  (40)

When the distance measurement signal is transmitted from the device A to the device B at the phase θ_GMSK,TXa, an estimated arrival time at the device B is a time longer than the elapsed time t_b2 estimated with θ_GMSK,TXa=θ_OFF by 2t_a1, and accordingly, the received phase is θ_GMSK,ab2. Thus, Expression (41) holds.

2t_a1+t_b2≈{1/(2πf₀₁)}sin⁻¹{(θ_GMSK,ab2˜θ_OFF)/A₀₁}  (41)

Alternatively, in a case of the 1.5 round trip, a left-hand side of Expression (43) is rewritten as a right hand side with interpretation of Expression (42) as described above.

2t_a1=t_a2  (42)

2t_a1+t_b2=t_a1+t_a2+t_b2  (43)

A phase detected by the device B having received the initial distance measurement signal from the device A is obtained by changing b to b1 in Expression (3), and thus the elapsed time t_b1 obtained from the detected phase is given by Expression (44).

t _b1≈{1/(2πf₀₁)}sin⁻¹{(θ_GMSK,b1˜θ_OFF)/A₀₁}  (44)

At the device B, t_b1 and (t_a1+t_a2+t_b2) can be estimated from Expression (41), Expression (43), and Expression (44). Expression (39) can be rewritten to obtain Expression (45).

t _x=(t_x1+t_x2)/2={2T_S˜(t_a1+t_b1+t_a2+t_b2)}/4  (45)

Thus, the flight distance can be estimated by using Expression (45) with addition of t_b1 and (t_a1+t_a2+t_b2) estimated at the device B.

Note that although not illustrated, it is clear that the method of performing phase detection K times for each m samples and averaging elapsed times obtained from detected phases in the second embodiment is also applicable to the present embodiment.

According to the fifth embodiment, 1.5 round trip distance measurement in the third embodiment and phase loop-back distance measurement in the fourth embodiment can be combined to estimate t_a1+t_b1+t_a2+t_b2 at one of the device A and the device B (in the above description, the device B). Accordingly, the data communication time can be further reduced, and the distance measurement time period can be further reduced.

Sixth Embodiment

In a sixth embodiment, any part same as a part in the first to fifth embodiments is denoted by the same reference sign, description of the part is omitted, and difference will be mainly described.

The elapsed times t_a1, t_a2, t_b1, and t_b2 until sampling since crossing of the offset phase θ_OFF are calculated by using, for example, Expressions (2) and (3) based on phases θ_GMSK,ai and θ_GMSK,bi (i=1, 2) sampled first after the offset phase θ_OFF.

Assume that the device A has an accurate clock and the device B has frequency error as defined in Expressions (7) and (8). Expression (7) is given as Expression (46) by using the reference clock frequency f_A of the device A and the reference clock frequency f_B of the device B.

f _A /f _B=1˜r_err  (46)

The frequency of the RF signal (distance measurement signal) generated at each of the device A and the device B is a frequency obtained by multiplying the reference clock frequency by a predetermined integer multiple, and thus error occurs to the frequency of the RF signal at a ratio same as a ratio for the reference clock frequency. Specifically, when an RF frequency of the device A is represented by f_RFA and an RF frequency of the device B is represented by f_RFB, Expression (47) is obtained.

f _RFA /f _RFB=1˜r_err  (47)

Thus, the RF frequency f_RFB of the device B is higher than the RF frequency f_RFA of the device A by r_errf_RFB. Thus, when the RF frequencies of the device A and the device B are set to be equal, the RF signal of the device B is received by the device A not as a direct current (DC) signal but as a signal of the frequency r_errf_RFB.

When a frequency corresponding to this difference is represented by f_A=r_errf_RFB and a difference between initial phases of the device B and the device A is represented by Δθ_BA, the phase of the distance measurement signal received by the device A is given by Expression (48).

θ_GMSK,a≈A₀₁ sin(2πf₀₁t)+2πf_Δt+θ_OFF+Δθ_BA  (48)

A time needed to be calculated is a time of crossing the offset phase θ_OFF with no frequency offset between the device A and the device B. Phase change is maximum at the time of crossing the offset phase θ_OFF. The maximum phase change at the time of crossing the offset phase θ_OFF occurs irrespective of the frequency offset f_Δ and the initial phase difference Δθ_BA. Thus, a time at which a value obtained by differentiating Expression (48) with respect to time is maximum is a time at which the distance measurement signal arrives at the device A. The time differential of Expression (48) is given by Expression (49).

dθ _GMSK,a /dt≈2πf₀₁A₀₁ cos(2πf₀₁t)+2πf_Δ  (49)

Expression (49) has a maximum value of 2πf₀₁A₀₁+2πf_Δ at t=0 [s].

When there is no frequency offset because the reference clocks of the device A and the device B are completely identical to each other and when initial phases are identical to each other, the distance measurement signal of Expression (1) is received by the device A. A time differential of Expression (1) has a maximum value of 2πf₀₁A₀₁ at a time of t=0 [s]. Although the maximum value has a difference due to a frequency difference, a time t_arvA of arrival at the device A can be detected by detecting a time at which a time differential of the phase has a maximum value.

However, since the sampling period cannot be infinitesimally small but is finite, the time at which the differential of the phase has a maximum value needs to be estimated from phases acquired at sample points. When, for example, phases at three points are sampled per quadrant to increase accuracy of the estimation, the sampling frequency is 12f₀₁. The following describes processing of estimating the arrival time t_arvA from sampled phases.

The estimation of the arrival time t_arvA is performed through a difference operation replacing a differential operation based on Expression (48). The offset phase θ_OFF+Δθ_BA in the phase signal of Expression (48) is replaced with θ_FFA=θ_OFF˜Δθ_BA for simplification.

First, the device A observes the phase θ_GMSK,a of a received phase signal for a duration equal to or longer than the modulation period T₀₁=1/f₀₁. The modulation period is known through beforehand communication between the device A and the device B, and thus the number of samples is determined accordingly. In the present assumption, a maximum value of a phase change amount can be detected through phase observation at 12 samples or more in the sampling period T_S. The 12 samples are represented by t_i (i=0 to 11). Sampling is performed a larger number of times in reality, in which phase observation is performed at 12 samples or more, but for description, the arrival time t_arvA estimated by the device A based on particular 12 samples is calculated.

The maximum value of the change amount of the received phase θ_GMSK,a can be calculated by calculating a value of Expression (50) for i=0 to 10 and determining a maximum value among the values.

Δθ_GMSK,a(t_i)=θ_GMSK,a(t_i+1)˜θ_GMSK,a(t_i) for i=0 to 10  (50)

Three samples are extracted based on the sample time t_i that provides the maximum value of the phase change amount. Note that estimation accuracy of the arrival time t_arvA can be further increased with samples in the number of three or larger, but the present description is made on a process of estimating the arrival time t_sumA with three samples.

FIG. 16 illustrates a relation between θ_GMSK,a(t_i) and the sample point t_i when the maximum value of Δθ_GMSK,a(t_i) is provided at i=1, in other words, when the arrival time t_arvA of the distance measurement signal from the device B, which is to be calculated, is between the sample time t₁ and the sample time t₂. FIG. 16 is a time chart illustrating an example in which the arrival time t_arvA of the distance measurement signal is between the sample time t₁ and the sample time t₂ in the present embodiment.

In this case, three sample points used to estimate the arrival time t_arvA are t₁, t₂, and t₃. Note that the same result can be obtained when the sample points are t₀, t₁, and t₂, but description of such a case is omitted.

Formulae that associate the detection phase θ_GMSK,a(t_i) with the sample points of times t₁, t₂, and t₃ are Expressions (51) to (53) below.

θ_GMSK,a(t₁)≈A₀₁ sin(2πf₀₁t₁)+2πf_Δt₁+θ_OFFA  (51)

θ_GMSK,a(t₂)≈A₀₁ sin(2πf₀₁t₂)+2πf_Δt₂+θ_OFFA  (52)

θ_GMSK,a(t₃)≈A₀₁ sin(2πf₀₁t₃)+2πf_Δt₃+θ_OFFA  (53)

Since T_S=t₂˜t₁=t₃˜t₂ holds, Expression (54) is obtained by subtracting, from a result obtained by subtracting Expression (52) from Expression (53), a result subtracting Expression (51) from Expression (52).

θ_GMSK,a(t₃)˜2θ_GMSK,a(t₂)+θ_GMSK,a(t₁)=A₀₁{sin(2πf₀₁t₃)˜2 sin(2πf₀₁t₂)+sin(2πf₀₁t₁)}  (54)

Expression (57) is obtained when the sample times t₁ and t₃ in Expression (54) are expressed as Expressions (55) and (56), respectively, by using t₂ and substituted into Expression (54).

t ₁ =˜T _S +t ₂  (55)

t ₃ =T _S +t ₂  (56)

θ_GMSK,a(t₃)˜2θ_GMSK,a(t₂)+θ_GMSK,a(t₁)=2A₀₁ sin(2πf₀₁t₂){cos(2πf₀₁T_S)˜1}  (57)

Since an elapsed time from the arrival time t_arvA between the sample time t₁ and the sample time t₂ to the sample time t₂ is t_a, t₂ is expressed by Expression (58).

t ₂ =t _arvA +t _a  (58)

Since the arrival time t_arvA is a time at which the phase of the sinusoidal signal as the distance measurement signal becomes the offset phase θ_OFF, Expression (59) is obtained.

sin(2πf₀₁t_arvA)=0  (59)

When sin(2πf₀₁t₂) in Expression (57) is rewritten by using Expression (58), Expression (60) is obtained.

sin(2πf₀₁t₂)=sin{2πf₀₁(t_arvA+t_a)}=sin(2πf₀₁t_a)  (60)

Expression (61) is obtained from Expressions (57) and (60).

sin(2πf₀₁t_a)=[θ_GMSK,a(t₃)˜2θ_GMSK,a(t₂)+θ_GMSK,a(t₁)]/[2A₀₁{cos(2πf₀₁T_S)˜1}]  (61)

Expression (61) can be rewritten as Expression (62), which is an expression that calculates t_a.

t _a={1/(2πf₀₁)}sin⁻¹[{θ_GMSK,a(t₃)˜2θ_GMSK,a(t₂)+θ_GMSK,a(t₁)}/(2A₀₁{cos(2πf₀₁T_S)˜1})]  (62)

Expression (62) indicates that t_a can be calculated from a result of the phase detection at three points, the sampling period T_S, the frequency f₀₁ of the GMSK signal of 0 and 1, and the amplitude A₀₁. From Expression (58), the arrival time t_arvA of the distance measurement signal is earlier than the sample time t₂ by the elapsed time t_a. Note that, as described up to the fifth embodiment, t_a calculated by Expression (62) is needed to calculate the flight time t_x.

It is trivial that the same processing can be performed to calculate the elapsed time t_b from an arrival time point of the distance measurement signal detected by the device B to the sample point, and t_a1, t_b1, t_a2, and t_b2 detected in an alternation sequence.

According to the sixth embodiment, the times t_arvA and t_arvB at which the distance measurement signal arrives can be accurately calculated at Devices A and B described in the first to fifth embodiments.

Seventh Embodiment

In the seventh embodiment, any part same as a part in the first to sixth embodiments is denoted by the same reference sign, description of the part is omitted, and difference will be mainly described.

In the sixth embodiment, the elapsed time t_a from the arrival time point of the distance measurement signal to the sample point is calculated by extracting only three points from among 12 sample points. The present embodiment describes a method to which the method of the sixth embodiment is extended and that calculates t_a by using four or more sample points. Note that the description assumes that conditions on the sample points are identical to conditions in the sixth embodiment.

FIG. 17 is a time chart for description of a method of calculating t_a by using four or more sample points in the present embodiment.

As illustrated in FIG. 17, times of 12 samples in total, six samples before and six samples after the arrival time t_arvA of the distance measurement signal, are sequentially referred to as t₈, t₉, t₁₀, t₁₁, t₀, t₁, t₂, t₃, t₄, t₅, t₆, and t₇.

The sixth embodiment uses samples at t₁ and t₃, one sample time T_S before and after the sample time t₂ when the elapsed time t_a has elapsed since the arrival time t_arvA. In contrast, the present embodiment uses samples at t_mod(2˜i) and t_2+i, i sample times iT_S before and after the sample time t₂, where mod represents a residue of 12, and for example, mod(˜1)=11. Exemplary calculation in a case of i=2 will be described below.

In a case of i=2, sample times other than the sample time t₂ are t₀ and t₄, and phases are given by Expressions (63) and (64). Note that Expression (52) at the sample time t₂ is listed below again.

θ_GMSK,a(t₀)≈A₀₁ sin(2πf₀₁t₀)+27πf_Δt₀+θ_OFFA  (63)

θ_GMSK,a(t₂)≈A₀₁ sin(2πf₀₁t₂)+2πf_Δt₂+θ_OFFA  (52)

θ_GMSK,a(t₄)≈A₀₁ sin(2πf₀₁t₄)+2πf_Δt₄+θ_OFFA  (64)

Similarly to derivation of Expression (54) by using Expressions (51) to (53), Expression (65) can be obtained by using Expressions (63). (52), and (64).

θ_GMSK,a(t₄)˜2θ_GMSK,a(t₂)+θ_GMSK,a(t₀)=A₀₁{sin(2πf₀₁t₄)˜2 sin(2πf₀₁t₂)+sin(2πf₁₀t₀)}  (65)

In addition, Expressions (66) and (67) below hold.

t ₀=˜2T_S+t₂  (66)

t ₄=2T_S+t₂  (67)

Expression (65) can be rewritten as Expression (68) by using Expressions (66) and (67).

θ_GMSK,a(t₄)˜2θ_GMSK,a(t₂)+θ_GMSK,a(t₀)=2A₀₁ sin(2πf₀₁t₂){cos(2πf₀₁2T_S)˜1}  (68)

Processing same as processing for obtaining a result of Expression (62) is performed to obtain the elapsed time t_a as in Expression (69).

t _a={1/(2πf₀₁)}sin⁻¹[{(θ_GMSK,a(t₃)˜2θ_GMSK,a(t₂)+θ_GMSK,a(t₁)}/(2A₀₁{cos(2πf₀₁2T_S)˜1})]  (69)

When calculated through the same processing for i≤5 in the present assumption, the elapsed time t_a is given by Expression (70).

t _a={1/(2πf₀₁)}sin⁻¹[{θ_GMSK,a(t_2+i)˜2θ_GMSK,a(t₂)+θ_GMSK,a(t_{mod (2˜i)})}/(2A₀₁{cos(2πf₀₁iT_S)˜1})]  (70)

When the elapsed time t_a calculated by Expression (70) for i=1, 2, 3, 4, and 5 is newly written an elapsed time t_a,i, the elapsed time t_a can be calculated by Expression (71) as an average value over a large number of sample points.

t _a=(⅕)×Σ_i=1⁵t_a,i  (71)

However, when the number of phase sample points for the frequency f₀₁ of the GMSK signal of 0 and 1 is 12 as in the present assumption, a sample at t₈, which is shifted from the sample time t₂ by half the period, is not used in averaging. To avoid this, the number of sample points in one period needs to be an odd number.

Note that the averaging of the elapsed time t_a may be performed by averaging the phase term of sin⁻¹ in Expression (70) and then converting the average into time. In this case, Expression (72) is obtained.

t _a={1/(2πf₀₁)}sin⁻¹[(⅕)×Σ_i=1⁵{θ_GMSK,a(t_2+i)˜2θ_GMSK,a(t₂)+θ_GMSK,a(t_mod(2˜i))}/(2A₀₁{cos(2πf₀₁iT_S)˜1})]  (72)

Thus, averaging over a large number of sample points is possible through phase averaging as well.

Note that it is trivial that the same processing can be performed to calculate the elapsed time t_b from the arrival time point of the distance measurement signal detected by the device B to the sample point, and t_a1, t_b1, t_a2, and t_b2 detected through an alternation sequence.

According to the seventh embodiment, effects substantially same as effects of the first to sixth embodiments described above are achieved, and in addition, it is possible to easily obtain a time until sampling since the distance measurement signal arrives at each device, which is necessary for calculating the distance. The same result can be obtained through phase averaging in place of time averaging.

While certain embodiments have been described, these embodiments have been presented by way of example only, and are not intended to limit the scope of the inventions. Indeed, the novel devices described herein may be embodied in a variety of other forms; furthermore, various omissions, substitutions and changes in the form of the devices described herein may be made without departing from the spirit of the inventions. The accompanying claims and their equivalents are intended to cover such forms or modifications as would fall within the scope and spirit of the inventions.

Claims

1. A distance measurement device comprising: a first device including a first reference signal source configured to generate a first reference signal, anda first transmitter-receiver configured to transmit a modulated first distance measurement signal, receive a modulated second distance measurement signal, obtain a second demodulated signal through demodulation, and acquire second phase information of the second demodulated signal in a first sampling period based on the first reference signal;a second device including a second reference signal source configured to operate independently from the first reference signal source and generate a second reference signal, anda second transmitter-receiver configured to transmit the second distance measurement signal, receive the first distance measurement signal, obtain a first demodulated signal through demodulation, and acquire first phase information of the first demodulated signal in a second sampling period based on the second reference signal; anda calculation unit configured to calculate a distance between the first device and the second device based on the first phase information and the second phase information, whereinone of the first distance measurement signal and the second distance measurement signal is transmitted once or more, and another of the first distance measurement signal and the second distance measurement signal is transmitted twice or more, andthe calculation unit calculates the distance based on a total of three or more pieces of the first phase information and the second phase information acquired through transmission of the first distance measurement signal and the second distance measurement signal three times or more in total, the first sampling period, and the second sampling period.

2. The distance measurement device according to claim 1, wherein when calculating the distance, the calculation unit corrects a time offset between the first sampling period and the second sampling period based on first phase information and second phase information obtained through loop-back transmission of the first distance measurement signal and the second distance measurement signal, andcorrects a frequency offset between the first sampling period and the second sampling period based on the first phase information and the second phase information obtained through the loop-back transmission and at least one of first phase information and second phase information obtained through transmission one or more times of at least one of the first distance measurement signal and the second distance measurement signal other than the loop-back transmission.

3. The distance measurement device according to claim 2, wherein the calculation unit calculates the distance based on two or more pieces of the first phase information and one or more pieces of the second phase information obtained through a distance measurement sequence in which the first transmitter-receiver transmits the first distance measurement signal, and the second transmitter-receiver receives the first distance measurement signal,the second transmitter-receiver transmits the second distance measurement signal, and the first transmitter-receiver receives the second distance measurement signal, andthe first transmitter-receiver transmits the first distance measurement signal, and the second transmitter-receiver receives the first distance measurement signal.

4. The distance measurement device according to claim 2, wherein the calculation unit calculates the distance based on two or more pieces of the first phase information and two or more pieces of the second phase information obtained through a distance measurement sequence in which the first transmitter-receiver transmits the first distance measurement signal, and the second transmitter-receiver receives the first distance measurement signal,the second transmitter-receiver transmits the second distance measurement signal, and the first transmitter-receiver receives the second distance measurement signal,the second transmitter-receiver transmits the second distance measurement signal, and the first transmitter-receiver receives the second distance measurement signal, andthe first transmitter-receiver transmits the first distance measurement signal, and the second transmitter-receiver receives the first distance measurement signal.

5. The distance measurement device according to claim 2, wherein the calculation unit calculates the distance based on two or more pieces of the first phase information and one or more pieces of the second phase information obtained through a distance measurement sequence in which the first transmitter-receiver transmits the first distance measurement signal, and the second transmitter-receiver receives the first distance measurement signal,the second transmitter-receiver transmits the second distance measurement signal, and the first transmitter-receiver receives the second distance measurement signal, andthe first transmitter-receiver sets a phase of the received second distance measurement signal as a phase at transmission start and transmits the first distance measurement signal, and the second transmitter-receiver receives the first distance measurement signal.

6. The distance measurement device according to claim 1, wherein the first transmitter-receiver samples the second phase information for one period of the second demodulated signal or longer,the second transmitter-receiver performs, twice or more, sampling of the first phase information for one period of the first demodulated signal or longer, andthe calculation unit obtains, for each of the second demodulated signal and the first demodulated signal, two samples, phase change of which is maximized through subtraction of a temporally previous sample of a sample,calculates, by using three continuous samples including the two samples of the first demodulated signal for a first time, an elapsed time tb1 of the second device from a time point at which the reception of the first demodulated signal is started, to a time point at which the first demodulated signal is sampled is obtained,calculates, by using three continuous samples including the two samples of the second demodulated signal, an elapsed time ta1 of the first device from a time point at which the reception of the second demodulated signal is started, to a time point at which the second demodulated signal is sampled is obtained,calculates, by using three continuous samples including the two samples of the first demodulated signal for a second time, an elapsed time tb2 of the second device from a time point at which the reception of the first demodulated signal is started, to a time point at which the first demodulated signal is sampled is obtained,calculates a first flight time tx,err1 including error by tx,err1=(TS˜ta1˜tb1)/2where TS represents the first or second sampling period,calculates a second flight time tx,err2 including error by tx,err2=(TS˜ta1˜tb2)/2calculates an average flight time by averaging the first flight time tx,err1 and the second flight time tx,err2 andcalculates the distance based on the average flight time.

7. The distance measurement device according to claim 6, wherein the first transmitter-receiver performs, twice or more, sampling of the second phase information for one period of the second demodulated signal or longer, andthe calculation unit calculates, by using three continuous samples including the two samples of the second demodulated signal for a second time, an elapsed time ta2 of the first device from a time point at which the reception of the second demodulated signal is started, to a time point at which the second demodulated signal is sampled is obtained,calculates the second flight time tx,err2 by tx,err2=(TS˜ta2˜tb2)/2in place of the calculation expression,calculates an average flight time by averaging the first flight time tx,err1 and the second flight time tx,err2, andcalculates the distance based on the average flight time.

8. The distance measurement device according to claim 6, wherein the calculation unit uses, as the elapsed time ta1, an average elapsed time obtained by averaging a plurality of elapsed times of the first device obtained based on the second distance measurement signal for two periods or longer, anduses, as the elapsed time tb1 and the elapsed time tb2, an average elapsed time obtained by averaging a plurality of elapsed times of the second device obtained based on the first distance measurement signal for two periods or longer.

9. The distance measurement device according to claim 1, wherein the first transmitter-receiver samples the second phase information for one period of the second demodulated signal or longer,the second transmitter-receiver performs, twice or more, sampling of the first phase information for one period of the first demodulated signal or longer,when k is a natural number, the calculation unit obtains, for each of the second demodulated signal and the first demodulated signal, two samples, phase change of which is maximized through subtraction of a temporally previous sample of a sample, and sets a reference sample to be a temporally later sample among the two samples, phase change of which is maximized, andfor three samples of the reference sample and two samples temporally k samples before and k samples after the reference sample, the calculation unit calculates, by using the three samples of the first demodulated signal for a first time, an elapsed time tb1 of the second device from a time point at which the reception of the first demodulated signal is started, to a time point at which the first demodulated signal is sampled is obtained,calculates, by using the three samples of the second demodulated signal, an elapsed time ta1 of the first device from a time point at which the reception of the second demodulated signal is started, to a time point at which the second demodulated signal is sampled is obtained,calculates, by using the three samples of the first demodulated signal for a second time, an elapsed time tb2 of the second device from a time point at which the reception of the first demodulated signal is started, to a time point at which the first demodulated signal is sampled is obtained,calculates a first flight time tx,err1 including error by tx,err1=(TS˜ta1˜tb1)/2where TS represents the first or second sampling period,calculates a second flight time tx,err2 including error by tx,err2=(TS˜ta1˜tb2)/2calculates an average flight time by averaging the first flight time tx,err1 and the second flight time tx,err2, andcalculates the distance based on the average flight time.

10. The distance measurement device according to claim 9, wherein when kmax is a maximum value of the k, the calculation unit uses, as the elapsed time ta1, an average elapsed time obtained by averaging kmax or less elapsed times of the first device obtained based on the second distance measurement signal for two periods or longer, anduses, as the elapsed time tb1 and the elapsed time tb2, an average elapsed time obtained by averaging kmax or less elapsed times of the second device obtained based on the first distance measurement signal for two periods or longer.

11. A distance measurement method of calculating a distance between a first device and a second device, wherein the first device generates a first reference signal,the first device transmits a modulated first distance measurement signal, receives a modulated second distance measurement signal, obtains a second demodulated signal through demodulation, and acquires second phase information of the second demodulated signal in a first sampling period based on the first reference signal,the second device generates a second reference signal independent from the first reference signal,the second device transmits the second distance measurement signal, receives the first distance measurement signal, obtains a first demodulated signal through demodulation, and acquires first phase information of the first demodulated signal in a second sampling period based on the second reference signal,a distance between the first device and the second device is calculated based on the first phase information and the second phase information,one of the first distance measurement signal and the second distance measurement signal is transmitted once or more, and another of the first distance measurement signal and the second distance measurement signal is transmitted twice or more, andthe distance is calculated based on a total of three or more pieces of the first phase information and the second phase information acquired through transmission of the first distance measurement signal and the second distance measurement signal three times or more in total, the first sampling period, and the second sampling period.

12. The distance measurement method according to claim 11, wherein when the distance is calculated, a time offset between the first sampling period and the second sampling period is corrected based on first phase information and second phase information obtained through loop-back transmission of the first distance measurement signal and the second distance measurement signal, anda frequency offset between the first sampling period and the second sampling period is corrected based on the first phase information and the second phase information obtained through the loop-back transmission and at least one of first phase information and second phase information obtained through transmission one or more times of at least one of the first distance measurement signal and the second distance measurement signal other than the loop-back transmission.

13. The distance measurement method according to claim 12, wherein the distance is calculated based on two or more pieces of the first phase information and one or more pieces of the second phase information obtained through a distance measurement sequence in which the first device transmits the first distance measurement signal, and the second device receives the first distance measurement signal,the second device transmits the second distance measurement signal, and the first device receives the second distance measurement signal, andthe first device transmits the first distance measurement signal, and the second device receives the first distance measurement signal.

14. The distance measurement method according to claim 12, wherein the distance is calculated based on two or more pieces of the first phase information and two or more pieces of the second phase information obtained through a distance measurement sequence in which the first device transmits the first distance measurement signal, and the second device receives the first distance measurement signal,the second device transmits the second distance measurement signal, and the first device receives the second distance measurement signal,the second device transmits the second distance measurement signal, and the first device receives the second distance measurement signal, andthe first device transmits the first distance measurement signal, and the second device receives the first distance measurement signal.

15. The distance measurement method according to claim 12, wherein the distance is calculated based on two or more pieces of the first phase information and one or more pieces of the second phase information obtained through a distance measurement sequence in which the first device transmits the first distance measurement signal, and the second device receives the first distance measurement signal,the second device transmits the second distance measurement signal, and the first device receives the second distance measurement signal, andthe first device sets a phase of the received second distance measurement signal as a phase at transmission start and transmits the first distance measurement signal, and the second device receives the first distance measurement signal.

16. The distance measurement method according to claim 11, wherein the first device samples the second phase information for one period of the second demodulated signal or longer,the second device performs, twice or more, sampling of the first phase information for one period of the first demodulated signal or longer,two samples, phase change of which is maximized through subtraction of a temporally previous sample of a sample are obtained for each of the second demodulated signal and the first demodulated signal,an elapsed time tb1 of the second device from a time point at which the reception of the first demodulated signal is started, to a time point at which the reception of the first demodulated signal is sampled is calculated by using three continuous samples including the two samples of the first demodulated signal for a first time,an elapsed time ta1 of the first device from a time point at which the reception of the second demodulated signal is started, to a time point at which the second demodulated signal is sampled is calculated by using three continuous samples including the two samples of the second demodulated signal,an elapsed time tb2 of the second device from a time point at which the reception of the first demodulated signal is started, to a time point at which the first demodulated signal is sampled is calculated by using three continuous samples including the two samples of the first demodulated signal for a second time,a first flight time tx,err1 including error is calculated by tx,err1=(TS˜ta1˜tb1)/2where TS represents the first or second sampling period,a second flight time tx,err2 including error is calculated by tx,err2=(TS˜ta1˜tb2)/2an average flight time is calculated by averaging the first flight time tx,err1 and the second flight time tx,err2, andthe distance is calculated based on the average flight time.

17. The distance measurement method according to claim 16, wherein the first device performs, twice or more, sampling of the second phase information for one period of the second demodulated signal or longer,an elapsed time ta2 of the first device from a time point at which the reception of the second demodulated signal is started, to a time point at which the second demodulated signal is sampled is calculated by using three continuous samples including the two samples of the second demodulated signal for a second time,the second flight time tx,err2 is calculated by tx,err2=(TS˜ta2˜tb2)/2in place of the calculation expression,an average flight time is calculated by averaging the first flight time tx,err1 and the second flight time tx,err2, andthe distance is calculated based on the average flight time.

18. The distance measurement method according to claim 16, wherein an average elapsed time obtained by averaging a plurality of elapsed times of the first device obtained based on the second distance measurement signal for two periods or longer is used as the elapsed time ta1, andan average elapsed time obtained by averaging a plurality of elapsed times of the second device obtained based on the first distance measurement signal for two periods or longer is used as the elapsed time tb1 and the elapsed time tb2.

19. The distance measurement method according to claim 11, wherein the first device samples the second phase information for one period of the second demodulated signal or longer,the second device performs, twice or more, sampling of the first phase information for one period of the first demodulated signal or longer,when k is a natural number, two samples, phase change of which is maximized through subtraction of a temporally previous sample of a sample are obtained for each of the second demodulated signal and the first demodulated signal, a reference sample is set to be a temporally later sample among the two samples, phase change of which is maximized,for three samples of the reference sample and two samples temporally k samples before and k samples after the reference sample, an elapsed time tb1 of the second device from a time point at which the reception of the first demodulated signal is started, to a time point at which the first demodulated signal is sampled is calculated by using the three samples of the first demodulated signal for a first time,an elapsed time ta1 of the first device from a time point at which the reception of the second demodulated signal started, to a time point at which the second demodulated signal is sampled is calculated by using the three samples of the second demodulated signal,an elapsed time tb2 of the second device from a time point at which the reception of the first demodulated signal is started, to a time point at which the first demodulated signal is sampled is calculated by using the three samples of the first demodulated signal for a second time,a first flight time tx,err1 including error is calculated by tx,err1=(TS˜ta1˜tb1)/2where TS represents the first or second sampling period,a second flight time tx,err2 including error is calculated by tx,err2=(TS˜ta1˜tb2)/2an average flight time is calculated by averaging the first flight time tx,err1 and the second flight time tx,err2, andthe distance is calculated based on the average flight time.

20. The distance measurement method according to claim 19, wherein when kmax is a maximum value of the k, an average elapsed time obtained by averaging kmax or less elapsed times of the first device obtained based on the second distance measurement signal for two periods or longer is used as the elapsed time ta1, andan average elapsed time obtained by averaging kmax or less elapsed times of the second device obtained based on the first distance measurement signal for two periods or longer is used as the elapsed time tb1 and the elapsed time tb2.