Method and apparatus for measuring displacement of object using multiple frequency signal
10785083 ยท 2020-09-22
Assignee
Inventors
Cpc classification
G01S13/536
PHYSICS
H04L27/3427
ELECTRICITY
H04L27/2273
ELECTRICITY
H04L27/366
ELECTRICITY
G01S13/538
PHYSICS
H04L27/389
ELECTRICITY
H04J1/06
ELECTRICITY
International classification
H04L27/34
ELECTRICITY
G01S13/538
PHYSICS
Abstract
A method and an apparatus for measuring a displacement of an object according to steps of: dividing a signal into an I signal and a Q signal according to a phase of the signal, wherein the signal is reflected by the object after a transmission signal having a plurality of frequencies is emitted toward the object by the radar measurement system; estimating a direct current (DC) component from an N-tuple information acquired from the I signal and the Q signal; removing the estimated DC component to correct the I signal and the Q signal; and measuring the displacement of the object based on the corrected I signal and Q signal are provided.
Claims
1. A method for measuring a displacement of an object by a radar measurement system, comprising: dividing a signal into an in-phase (I) signal and a quadrature-phase (Q) signal according to a phase of the signal, wherein the signal is reflected by the object after a transmission signal having a plurality of frequencies is emitted toward the object by the radar measurement system; estimating coordinates of a direct current (DC) component in an I/O coordinate system based on a polygon formed by N constellation points corresponding to N-tuple information acquired from the I signal and the Q signal, where N is a non-negative integer; removing the DC component based on the coordinates of the DC component to correct the I signal and the Q signal; and measuring the displacement of the object based on the corrected I signal and Q signal, wherein the estimating of the DC component comprises: selecting at least three constellation points among the N constellation points corresponding to the N-tuple information; and determining coordinates of a circumcenter of a triangle formed by the at least three selected constellation points as the coordinates of the DC component.
2. The method of claim 1, wherein the selecting of the at least three constellation points comprises: comparing an area of the triangle with a predetermined threshold; and determining that the at least three constellation points are valid when the area of the triangle is larger than the predetermined threshold.
3. The method of claim 2, wherein the selecting of the at least three constellation points further comprises: selecting at least three constellation points which form a triangle different from the previously selected constellation points among the N constellation points when the area of the triangle is smaller than the predetermined threshold.
4. The method of claim 1, further comprising emitting the transmission signal corresponding to each frequency of the plurality of frequencies toward the object.
5. The method of claim 1, wherein the transmission signal is a signal in which a plurality of signals having different frequencies is multiplexed.
6. An apparatus for measuring a displacement of an object, comprising: a processor and a memory, wherein the processor executes a program stored in the memory to perform: dividing a signal into an in-phase (I) signal and a quadrature-phase (Q) signal according to a phase of the signal, wherein the signal is reflected by the object after a transmission signal having a plurality of frequencies is emitted toward the object; estimating coordinates of a direct current (DC) component in an I/Q coordinate system based on a polygon formed by N constellation points corresponding to N-tuple information acquired from the I signal and the Q signal, where N is a non-negative integer; removing the DC component based on the coordinates of the DC component to correct the I signal and the Q signal; and measuring the displacement of the object based on the corrected I signal and Q signal, wherein when the processor performs the estimating of the DC component, the processor executes the program to perform: selecting at least three constellation points among the N constellation points corresponding to the N-tuple information; and determining coordinates of a circumcenter of a triangle formed by the at least three selected constellation points as the coordinates of the DC component.
7. The apparatus of claim 6, wherein when the processor performs the selecting of the at least three constellation points, the processor executes the program to perform: comparing an area of the triangle with a predetermined threshold; and determining that the at least three constellation points are valid when the area of the triangle is larger than the predetermined threshold.
8. The apparatus of claim 7, wherein when the processor performs the selecting of the at least three constellation points, the processor executes the program to perform selecting at least three constellation points which form a triangle different from the previously selected constellation points among the N constellation points when the area of the triangle is smaller than the predetermined threshold.
9. The apparatus of claim 6, wherein the processor executes the program to further perform emitting the transmission signal corresponding to each frequency of the plurality of frequencies toward the object.
10. The apparatus of claim 6, wherein the transmission signal is a signal in which a plurality of signals having different frequencies is multiplexed.
11. A radar system for measuring a displacement of an object, comprising: a transmitter configured to radiate a transmission signal having a plurality of frequencies toward the object; and a receiver configured to receive a signal reflected by the object and measure the displacement of the object, wherein the receiver comprises a processor configured to execute a program stored in a memory of the receiver to perform: dividing a signal into an in-phase (I) signal and a quadrature-phase (Q) signal according to a phase of the signal, wherein the signal is reflected by the object after a transmission signal having a plurality of frequencies is emitted toward the object by the radar measurement system; estimating coordinates of a direct current (DC) component in an I/O coordinate system based on a polygon formed by N constellation points corresponding to N-tuple information acquired from the I signal and the Q signal, where N is a non-negative integer; removing the DC component based on the coordinates of the DC component to correct the I signal and the Q signal; and measuring the displacement of the object based on the corrected I signal and Q signal, wherein when the processor performs the estimating of the coordinates of the DC component, the processor executes the program to perform: selecting at least three constellation points among the N constellation points corresponding to the N-tuple information; and determining coordinates of a circumcenter of a triangle formed by the at least three selected constellation points as the coordinates of the DC component.
12. The radar system of claim 11, wherein when the processor performs the selecting of the at least three constellation points, the processor executes the program to perform: comparing an area of the triangle with a predetermined threshold; and determining that the at least three constellation points are valid when the area of the triangle is larger than the predetermined threshold.
13. The radar system of claim 12, wherein when the processor performs the selecting of the at least three constellation points, the processor executes the program to perform selecting at least three constellation points which form a triangle different from the previously selected constellation points among the N constellation points when the area of the triangle is smaller than the predetermined threshold.
14. The radar system of claim 11, wherein the processor executes the program to further perform radiating the transmission signal corresponding to each frequency of the plurality of frequencies toward the object.
15. The radar system of claim 11, wherein the transmission signal is a signal in which a plurality of signals having different frequencies is multiplexed.
Description
BRIEF DESCRIPTION OF THE DRAWINGS
(1)
(2)
(3)
(4)
(5)
(6)
(7)
(8)
(9)
(10)
(11)
(12)
DETAILED DESCRIPTION OF THE EMBODIMENTS
(13) In the following detailed description, only certain exemplary embodiments of the present invention have been shown and described, simply by way of illustration. As those skilled in the art would realize, the described embodiments may be modified in various different ways, all without departing from the spirit or scope of the present invention. Accordingly, the drawings and description are to be regarded as illustrative in nature and not restrictive, and like reference numerals designate like elements throughout the specification.
(14)
(15) Referring to
(16) Referring to
(17)
(18) In Equation 1, is phase noise and .sub.CW is a wavelength of the radiated signal of single frequency radar. According to an exemplary embodiment, periodic or aperiodic motion (i.e., displacement) x.sub.p(n) of the object within one wavelength appears as a trajectory on an I/Q plot and the displacement of the object may be calculated from the phase information of the I signal and the Q signal. Equation 2 represents an arc-tangential method for measuring the displacement x.sub.p(n) of the object within the one wavelength.
arctan(Q(t)/I(t)).fwdarw.x.sub.p(n)[Equation 2]
(19) However, when the DCRM is applied, it is difficult to accurately estimate the displacement of the object because the phase information of the I signal and the Q signal is changed. Therefore, the DC component in the I signal and the Q signal need to be removed without the DCRM.
(20) The DC component of the I signal and the Q signal may be generated by circuit characteristics of the radar measurement system or by the reflection characteristics of an environment surrounding the radar system. The DC component generated by the circuit characteristic is not variable, but the DC component caused by the environment of the radar system may be changed according to the change of the surroundings. For example, conventionally, the DCRM is not applied in a water level sensor or the like in a water tank. The DC component is likely to be expected because a reflection environment of the water tank in which the radar sensor is located is static. The radar for detecting the sleep level may convert the I signal and the Q signal including the DC component into a digital signal, and then estimate and remove the DC component by an algorithm. Saturation of the ADC may be avoided through a DC tuning circuit.
(21) Even when the radar measurement system is used to detect physiological signals of a human body, the space range for the radar detection may be narrow (within a few meters) and therefore, the DCRM cannot be applied. The I signal and the Q signal converted into the digital signal without the DCRM are represented by Equation 3 below.
(22)
(23) The DC component in the signal of Equation 3 may be estimated as follows. First, time streams [I(1), I(2), . . . , I(N)] of the I signal and [Q(1), Q(2), . . . , Q(N)] of the Q signal are collected. From the collected time streams, DC.sub.I, DC.sub.Q, and R which satisfy Equation 4 below are computed based on the l.sub.1-norm minimization algorithm. The DC.sub.I and the DC.sub.Q are center coordinates of the DC component, that is, center of a circle having the trajectory as an arc. The R is the radius of the circle having the trajectory as the arc.
min[d.sub.1. . . d.sub.N].sub.1(d.sub.i=(I(i)DC.sub.I).sup.2+(Q(i)DC.sub.Q).sup.2R.sup.2)[Equation 4]
(24) The time streams of the I signal and the Q signal need to be collected for a sufficient time to derive DC.sub.I, DC.sub.Q, and R through equation 4. Therefore, if the object moves slowly, or the object moves hardly, the trajectory is not changed, so it takes a considerable amount of time to collect the time streams.
(25)
(26) The radar measurement system according to the exemplary embodiment includes a transmitter 200 configured to radiate a transmission signal having a plurality of frequencies as shown in
(27) The transmitter 200 of the radar measurement system may control the continuous wave source to sequentially radiate the transmission signal by changing the frequency of the transmission signal. Referring to
(28) The received signal (I signal and Q signal) corresponding to the transmission signal having the i-th frequency among the N frequencies of
(29)
(30) Referring to Equation 5, the I signal and the Q signal have a wavelength .sub.i corresponding to the i-th frequency, respectively. The I signal and the Q signal of Equation 5 are the signals output from the ADC of
(31)
(32) Referring to
(33) The N-tuple acquiring processor 110 acquires N number of I signals and N number of Q signals, that is, N-tuples, from the received signals divided according to a phase. When the distance between the radar measurement system and the object is as short as several meters or less, the N-tuple acquiring processor 110 may acquire the N-tuple within a very short period of time (e.g., for milliseconds) or less since the radar measurement system according to the exemplary embodiment may radiate the transmission signals by changing the frequency at a very short interval. Equation 6 below represents an N-tuple of the I signal and the Q signal obtained by the N-tuple acquiring processor 110.
[I.sub.1(n),I.sub.2(n), . . . ,I.sub.N(n)],[Q.sub.1(n),Q.sub.2(n), . . . ,Q.sub.N(n)][Equation 6]
(34) When a symbol in which a plurality of signals having different frequencies are multiplexed is transmitted by the transmitter 200, the N-tuple acquiring processor 110 may perform N-point fast Fourier transform (FFT) for the received signal from the ADC so as to acquire the N-tuple from the multiplexed symbols. Referring to
(35) The DC component estimating processor 120 estimates the DC component of the I signal and the Q signal by using the N-tuple. According to the exemplary embodiment, the DC component estimating processor 120 may select at least three constellation points among N constellation points corresponding to the N-tuple, and calculate the circumcenter of a triangle formed by the at least three constellation points as the I-coordinate (DC.sub.I) and the Q-coordinate (DC.sub.Q). The at least three constellation points may include the first constellation point (I.sub.1(n),Q.sub.1(n)), the last constellation point (I.sub.N(n),Q.sub.N(n)), and the middle constellation point (I.sub.N/2(n),Q.sub.N/2(n)) ((I.sub.N/2+1(n),Q.sub.N/2+1(n))) where N is odd) of the N constellation points. Three of the N constellation points may be selected in different ways, and the present invention is not limited thereto.
(36) The N-tuple correcting processor 130 corrects the N-tuple based on the estimated DC component.
(37) The displacement measuring processor 140 measures the displacement of the object by using the corrected N-tuple.
(38) Referring to
(39)
(40) When x.sub.p(n) is much shorter than the wavelength in Equation 7, a postulation of Equation 8 may be stand.
A.sub.I,1(n)A.sub.I,2(n) . . . A.sub.I,N(n)=A.sub.I(n)
A.sub.Q,1(n)A.sub.Q,2(n) . . . A.sub.Q,N(n)=A.sub.Q(n)[Equation 8]
(41) Then, when f.sub.1>>(f.sub.Nf.sub.1), since the difference between the DC components corresponding to each of the frequencies is not large, the following Equation 9 holds.
DC.sub.I,1(n)DC.sub.I,2(n) . . . DC.sub.I,N(n)=DC.sub.I(n)
DC.sub.Q,1(n)DC.sub.Q,2(n) . . . DC.sub.Q,N(n)=DC.sub.Q(n)[Equation 9]
(42) When phase noise is =0, the Equation 7 may be expressed as Equation 10 below based on Equations 8 and 9.
(43)
(44) I.sub.i(n) and Q.sub.i(n) according to Equation 10 may be expressed as constellation points which are far from the center (DC.sub.I(n),DC.sub.Q(n)) by a distance {square root over (A.sub.I.sup.2(n)+A.sub.Q.sup.2(n))} at a phase
(45)
The phase
(46)
is a linear function of which slope is
(47)
as shown in
(48)
(49) The DC component estimating processor 120 may then estimate the DC components DC.sub.I(n) and DC.sub.Q(n) at time n by using an N-tuple at time n. According to an exemplary embodiment, for estimating the DC component, the DC component estimating processor 120 selects at least three constellation points among N constellation points corresponding to the N-tuple (S120). The DC component estimating processor 120 may determine the circumcenter of the triangle formed by the at least three constellation points as the coordinates of the DC component on the I/Q plot. The DC component estimating processor 120 determines whether the selected three constellation points are valid (S130). For example, the DC component estimating processor 120 may determine that the selected constellation points are valid when the area S of the triangle formed by the three constellation points is larger than a predetermined threshold. In addition, the DC component estimating processor 120 may determine that the at least three constellation points are valid when the area S is smaller than the predetermined threshold. The predetermined threshold may be predetermined according to the frequency of the transmission signal, the surrounding environment of the radar measurement system, and the like.
(50)
(51) If it is determined that the selected constellation points (alternatively, a group of the constellation points, an ordered pair of the constellation points) are not valid, the DC component estimating processor 120 re-selects at least three constellation points from among the N constellation points. The newly selected constellation points may be different from the previously selected constellation points. For example, all the newly selected constellation points may be different from the previously selected constellation points.
(52) Alternatively, some of the three constellation points may be different from the previously selected constellation points.
(53) If it is determined that the selected constellation points are valid, the DC component estimating processor 120 calculates the circumcenter of the triangle formed by the selected three constellation points (S140). The DC component estimating processor 120 may determine the circumcenter corresponding to the three points as the coordinates of the DC component of the I signal and the Q signal.
(54) The N-tuple correcting processor 130 corrects the N-tuple by using the coordinates of the estimated DC component (S150). Equation 11 below represents the corrected N-tuple.
I.sub.comp(n)=[I.sub.1(n),I.sub.2(n), . . . I.sub.N(n)]DC.sub.I(n)
Q.sub.comp(n)=[Q.sub.1(n),Q.sub.2(n), . . . Q.sub.N(n)]DC.sub.I(n)[Equation 11]
(55) Then, the displacement measuring processor 140 measures the displacement of the object by using the corrected N-tuple. The displacement of the object {circumflex over (x)}.sub.p,i(n) measured by the signal having the i-th frequency may be expressed as Equation 12 below, and the total displacement of the object {circumflex over (x)}.sub.p(n) may be expressed as Equation 13 below.
(56)
(57) The displacement measuring processor 140 according to an exemplary embodiment may determine the total displacement {circumflex over (x)}.sub.p(n) of the object based on the measured displacement {circumflex over (x)}.sub.p,i(n) corresponding to each frequency. For example, the displacement measuring processor 140 may determine the average value of the N partial displacements {circumflex over (x)}.sub.p,i(n) as the total displacement. Alternatively, the displacement measuring processor 140 may randomly select one of the N partial displacements {circumflex over (x)}.sub.p,i(n) to determine the selected value as the total displacement.
(58) The displacement of the object can be accurately measured even when the object is hardly moved by using the reflection signal of the multi-frequency signal. Further, the DC component of the received signal for measuring the displacement of the object may be estimated quickly from the tuples collected in a very short time period. Consequently, the movement of objects can be efficiently measured with a small computing resource, and miniaturization of an apparatus for acquiring biometric information can be realized.
(59)
(60) The receiver according to an exemplary embodiment may be implemented as a computer system, for example a computer readable medium. Referring to
(61) Thus, embodiments of the present invention may be embodied as a computer-implemented method or as a non-volatile computer-readable medium having computer-executable instructions stored thereon. In the exemplary embodiment, when executed by a processor, the computer-readable instructions may perform the method according to at least one aspect of the present disclosure. The network connection 1290 may transmit or receive a wired signal or a wireless signal.
(62) On the contrary, the embodiments of the present invention are not implemented only by the apparatuses and/or methods described so far, but may be implemented through a program realizing the function corresponding to the configuration of the embodiment of the present disclosure or a recording medium on which the program is recorded. Such an embodiment can be easily implemented by those skilled in the art from the description of the embodiments described above. Specifically, methods (e.g., network management methods, data transmission methods, transmission schedule generation methods, etc.) according to embodiments of the present disclosure may be implemented in the form of program instructions that may be executed through various computer means, and be recorded in the computer-readable medium. The computer-readable medium may include program instructions, data files, data structures, and the like, alone or in combination. The program instructions to be recorded on the computer-readable medium may be those specially designed or constructed for the embodiments of the present disclosure or may be known and available to those of ordinary skill in the computer software arts. The computer-readable recording medium may include a hardware device configured to store and execute program instructions. For example, the computer-readable recording medium can be any type of storage media such as magnetic media like hard disks, floppy disks, and magnetic tapes, optical media like CD-ROMs, DVDs, magneto-optical media like floptical disks, and ROM, RAM, flash memory, and the like. Program instructions may include machine language code such as those produced by a compiler, as well as high-level language code that may be executed by a computer via an interpreter, or the like.
(63) While this invention has been described in connection with what is presently considered to be practical example embodiments, it is to be understood that the invention is not limited to the disclosed embodiments, but, on the contrary, is intended to cover various modifications and equivalent arrangements included within the spirit and scope of the appended claims.