An LFMCW detector with new structure and FRFT based differential distance estimation method

This paper describes a linear frequency modulated continuous wave (LFMCW) detector which is designed for a collision avoidance radar. This detector can estimate distance between the detector and pedestrians or vehicles, thereby it will help to reduce the likelihood of traffic accidents. The detector consists of a transceiver and a signal processor. A novel structure based on the intermediate frequency signal (IFS) is designed for the transceiver which is different from the traditional LFMCW transceiver using the beat frequency signal (BFS) based structure. In the signal processor, a novel fractional Fourier transform (FRFT) based differential distance estimation (DDE) method is used to detect the distance. The new IFS based structure is beneficial for the FRFT based DDE method to reduce the computation complexity, because it does not need the scan of the optimal FRFT order. Low computation complexity ensures the feasibility of practical applications. Simulations are carried out and results demonstrate the efficiency of the detector designed in this paper.

and fractional Fourier transform (FRFT) (Grennberg and Sandell 1994;Lin and Chern 1998;Sharma and Joshi 2007;Zhou et al. 2014). Because of the advantage of the FRFT in terms of processing linear frequency modulation (LFM) signals, FRFT are attracting more and more attention in the field of LFM signal processing. However, high computation complexity coming from the scan of the optimal FRFT order restricts its practical applications. Therefore, to utilize FRFT solving FM signal problems in practical applications, a novel intermediate frequency signal (IFS) based structure is designed for the transceiver and a related differential distance estimation (DDE) method based on FRFT is implemented in the signal processor. The LFMCW detector made up of this transceiver and signal processor can estimate the distance effectively.
The rest of this paper is organized as follows. System model is given in "Preliminaries" section. "Proposed LFMCW detector" section describes the proposed IFS based structure and DDE method based on FRFT in detail. In "Simulation and results" section, the simulation results are presented. Finally, the conclusion is drew in the fifth section.

Preliminaries
As shown in Fig. 1, a collision avoidance radar can estimate the distance between a vehicle and itself to avoid that the two cars collide. Normally, an FM radar or a pulse radar is used as the collision avoidance radar. Figure 2 describes an FM radar traditional LFMCW detector with the traditional beat frequency signal (BFS) based structure. Usually, correlation, FFT or other methods in time or frequency domain are utilized in this detector to estimate the distance. However, the distance estimation accuracy of these methods is inversely proportional to the modulation band of the FM signal. It means that better distance estimation accuracy needs wider modulation band of the FM signal. Nevertheless, wider modulation band will bring other problems, such as frequency modulation nonlinear, which can affect the distance estimation performance. Therefore we hope a distance estimation method whose accuracy has nothing to do with the modulation band. Then a novel distance estimation method using FRFT is proposed in this paper. However, if we using the method in the LFMCW detector shown as Fig. 2, we have to search for the optimal rotation angle and this makes the computation complexity enhanced obviously. R C o ll is io n a v o id a n c e r a d a r

Fig. 1 System model
For avoiding the search for the optimal rotation angle, the LFMCW detector with new IFS based structure is proposed in paper and it is shown in Fig. 3. Compared with the traditional one, the decrease of computation complexity is based on the additional hardware cost. The additional hardware consists of a LPF and a mixer. Either the structure burden or the financial burden of the additional hardware can be ignored under the current level of technology. However, the decrease of computation complexity is multiple, because each search needs carrying out once FRFT.

Fractional Fourier transform
As a generalization of conventional Fourier transform (FT), the FRFT has been shown to be more effective than the FT in many signal processing areas. Several signal processing methods based on FRFT have been proposed and their properties also have been derived and discussed, such as convolution, filtering, correlation and so on (Ozaktas and Barshan 1994;Almeida 1994;Ozaktas et al. 2001). To make the FRFT come true, some fast discrete fractional Fourier transform (DFRFT) have been done (Ozaktas et al. 1996;Pei and Yeh 1999;Candan and Kutay 2000;Ran et al. 2001). Among them, the Ozaktas' DFRFT fast algorithm is widely regarded as an efficient one (Ozaktas et al. 1996), which has high precision, fast speed and the computation complexity matched with FFT.
The FRFT with a rotation angle α of signal f(t), denoted as F α (u), is defined as (Ozaktas et al. 2001 where the transform kernel function K α (u, t) is given as: in which t is time, and u is the coordinate in fractional Fourier domain (FRFD). Meanwhile j is the imaginary unit and the function exp() is the exponential function. The parameter α is the rotation angle of FRFT. When α increases from 0 to π/2, the FRFT transforms a continuous signal from its time domain to the Fourier image. If α or α + π is multiples of 2π, the kernel K α (u, t) is simplified as δ(t − u) or δ(t + u), respectively. FRFT with rotation angle α defined in the time-frequency domain is shown in Fig. 4. In this figure, t is time and f is frequency. Meanwhile u is the coordinate in FRFD and v is the amplitude in FRFD. f if1 (t) and f if2 (t) are instantaneous frequency of two signals with IF1 and IF2 as their initial frequency (IF). u p1 and u p2 are peak positions of the two signals in FRFD after FRFT. It is noticed that with the change of the rotation angle α, the axis u will perpendicular to the instantaneous frequency f if1 (t) and f if2 (t) when the rotation angle α = β + π/2. Under this condition, the relationship between the peak positions in FRFD and the initial frequency of the signal can be expressed as Ozaktas' DFRFT can be expressed as Ozaktas' et al. (1996) (1)

Fig. 4 Definition of FRFT in the time-frequency domain
where A α = exp (−jπsgn(sin α) / 4 + jα / 2) |sin α| 1/2 and γ = cot α· 1 2 x is the normalized sampling interval in both time domain and FRFD. This DFRFT can be realized by three steps: (1) Multiply a chirp signal in the time domain; (2) Carry out FFT to the signal after first step.
Multiply another chirp signal and we get the DFRFT of the original signal.
We defined that multiplication times of an algorithm is the computation complexity. As a result, the computation complexity of this DFRFT with N points are 2N + Nlog 2 N/2, which is not so more than FFT.

Proposed LFMCW detector
According to Figs. 2 and 3, it can be noticed that the structure of this detector is different from the traditional one. Traditional LFMCW detectors use the beat frequency signal (BFS) to gain the wanted information, whereas in the proposed detector information is acquired from the IFS. The mixer puts out BFS by mixing the echo signal with the transmit signal. However, IFS is obtained by mixing the echo signal with the reference carrier signal. Compared to BFS, the frequency modulation slope of the transmit signal is remained in the IFS. Because of the knowledge of the frequency modulation slope, FRFT can be carried out without hunting for the optimal order. As a result, the computation complexity of the FRFT based DDE method is significantly low.
Supposed that the LFM transmitting signal s 1 (t) with amplitude a 1 , carrier frequency f 1 and the frequency modulation slope K 1 is expressed as According to Fig. 2, this signal is mixed with an echo signal s 2 (t) which can be expressed as where τ(t) is the time delay and λ 1 is the amplitude factor. Then we get the BFS s bf (t) and it can be formulated as (R 0 is the initial distance and v 0 is the relative speed. c is the velocity of electromagnetic wave) and the s bf (t) can be expressed as where a bf is the amplitude. We can find that the frequency modulation slope of s bf (t) is and it is unknown. Consequently, we have to make search for the optimal rotation angle to determine it.
On the contrary, according to Fig. 3, the echo signal is mixed with a reference carrier signal s 0 (t) = cos2π(f 1 t), and we obtain the IFS s if (t) which can be expressed as where a if is the amplitude. Considering τ (t) = 2(R 0 −v 0 t) c , the equation can be rewritten as Because of the relationship v 0 ≪ c in practice, the frequency modulation slope is approximate K 1 . Therefore, the frequency modulation slope of s if (t) is known for us and the FRFT with the optimal rotation angle can be carried out without scan for the optimal rotation angle.
Additionally, when the relative speed is small enough, the time delay τ(t) can be regarded as a constant τ for each computation. So the equation can also be formulated as Normally, this equation is utilized for distance estimation. From this equation, we can also find that the frequency modulation slope of s if (t) is known for us and the FRFT with the optimal rotation angle can be carried out without scan for the optimal rotation angle.

Derivation of the DDE method
The principle block diagram of the FRFT based DDE method is shown in Fig. 5. Two IFSs are sampled by A/D and become discrete signals. Then DFRFT is carried out to them and locations of the two IFSs in fractional Fourier domain are estimated. By means of the locations, IFs of the two IFSs are obtained. Time delay is calculated through the difference of the two IFSs' IFs. Finally, distance is estimated according to the time delay.
Supposed that the modulating signal of the LFM signal is a saw tooth wave, then the instantaneous frequency f if (t) as shown in Fig. 6 can be denoted as: where f 0 is the carrier frequency and T is the period of the modulated signal. K = B/T is the frequency modulation slope and B is the modulation bandwidth. Parameter n is the cycle times. So the transmit signal within a modulation period can be expressed as: where Φ n 0 = n(n + 1)KT 2 /2 + Φ 0 (Φ 0 is the initial phase), and a is the amplitude of the transmit signal.
Then, the reflected echo signal can be described as: where Φ τ = (Kτ 2 − 2f 0 τ − 2nKTτ)/2 + n(n + 1)KT 2 /2 + Φ 1 , λ is the amplitude factor of the echo signal, and the parameter τ is the time delay. The transmit signal and the echo signal are mixed with the reference carrier signal x(t) = cos 2πf 0 t in two separate mixers respectively, then the IFSs can be derived as: It's noticed that the frequency modulation slope K is remained in the IFSs mentioned as above. So FRFT will be carried out to the IFSs with angle α = arctan(K) + π/2, without hunting for the optimal rotation angle.
Based on the Eq. (15), the IFs f it and f ir of the IFSs can be simply obtained as the Eq. (16).
(16) f it = −nKT f ir = −(nKT + K τ ). Obviously, the time delay τ can be estimated through the difference between the f it and f ir . According to Eq. (3), while the rotation angle (α = arctan(K) + π/2) is optimal, the relationship between f ir (f it ) and the peak position u pr (u pt ) in FRFD of the IFSs can be written as With the IFs obtained through the estimation of u pr (u pt ), the time delay estimation τ is achieved and expressed as: According to the Eq. (18) and τ = 2R/c (R is the distance and c is the velocity of electromagnetic wave), the distance estimation value can be formulated as:

Distance estimation precision and computation complexity
According to Eq. (19), the distance estimation precision of this method is proportional to the estimation accuracy of u pr (u pt ). The estimation accuracy of u pr (u pt ) is equal to the sampling interval in the fractional domain, and can be expressed as u = 1 √ N after normalization (Ozaktas et al. 1996). According to the Ozaktas et al. (1996), the real sampling interval Δu r can be expressed as where f s is the sampling frequency in time domain and N is the computation point for each processing. Therefore, the distance estimation precision ΔR is For engineering applications, f s = 10B is reasonable. Meanwhile, the frequency modulation slope K = B/T. Consequently, the distance estimation precision can be rewrite as Because parameter c and α are constants, the period of the modulated signal T and sampling points are the main influencing factors for the distance estimation precision. Improvement of the distance estimation precision depends on T and N, but not the modulation bandwidth. It is beneficial for engineering applications. Figure 7 shows the relationships between the distance estimation precision and the period of the modulated signal or the sampling points.
For the proposed DDE method in this paper, the multiplication times is used to measure its computation complexity. According to Fig. 4, the multiplication times mainly come from the DFRFT of the IFS. As a result, according to the computation complexity of DFRFT, the computation of this DDE method can be expressed as (17) f ir = u pr csc α f it = u pt csc α.
(20) u r = u f s T = f s N .
(22) �R = 5cT csc α N On the contrary, if searching the adaptive angle of DFRFT is needed, the computation complexity with M times scan can be written as Compared with Eq. (24), the computation complexity of the proposed DDE method described as Eq. (23) decreases obviously.

Simulation and results
In this section, simulations are carried out to validate the performance of the detector proposed in this paper. The carrier frequency is 3.0 GHz, and the modulation bandwidth is 30 MHz. The modulation period is 3.3 × 10 −6 s. When the target is stationary relative to the LFMCW detector, the time delay between the transmitting signal and the echo signal is a constant. Simulations were carried out for this scenario, and the results are shown in Fig. 8, where data1, data2 and data3 is the transmitting signal, echo signal with time delay of 1e−7 s and echo signal with time delay of 2e−7 s respectively. It can be seen that the locations in the fractional Fourier domain are different for the three signals with different initial frequencies. Moreover, to test the performance of the proposed DDE method with noise, a Monte Carlo simulation with 1000 points was also conducted. The results are shown in Fig. 9. The simulation results prove that the proposed method in this paper is feasible for a fixed-distance measurement when the signal-to-noise ratio (SNR) is no less than 0 dB. This method would certainly be feasible with lower SNR values if the noise was filtered out prior to applying the DFRFT.
If there is some relative velocity between the target and the LFMCW detector, the time delay value will vary. For example, when the target is gradually approaching the system, the time delay will gradually reduce. Simulations of this scenario were carried out, and the results of the instantaneous time delay estimation with SNR value 0 dB are shown in Fig. 10. This figure proves that the proposed DDE method in this paper is also feasible for the instantaneous distance measurements for a moving target.

Conclusion
In this paper, a novel FRFT based LFMCW detector is proposed to achieve distance measurement for a collision avoidance radar. On the one hand, an FRFT based DDE method is proposed for the detector to realize distance estimation. On the other hand, a new IFS based structure transceiver is designed for the detector to reduce computation complexity of FRFT by taking advantage of the knowledge of frequency modulation slope. Benefited from the FRFT based DDE method and structure, the detector can estimate the distance effectively and make the FRFT feasibly applied in practice rather than merely researched theoretically. Besides, practical FRFT applications in LFM systems have important significance for improving the performance of LFM systems.
Authors' contributions KY and XH made substantial contributions to conception and design of the advanced strategy of the algorithm proposed. PL revised the manuscript critically to meet the expected standards in scientific publishing. All authors read and approved the final manuscript.  Technology, in 1996Technology, in , 1999Technology, in and 2007 She is now an Associate Professor in Beijing Institute of Technology. Her main research interests include target detection theory of radio sensor and signal processing, real-time signal processing and time-frequency analysis.

Authors information
Ping Li received the B.S. and M.S. degree in Mechantronical Engineering from Dalian Jiaotong University, Dalian, China, in 1985 and1987, respectively, and the Ph.D. degrees from the Beijing Institute of Technology, China, in 1995. In Sept. 1996, she joined School of Mechatronical Engineering, Beijing Institute of Technology, China. Her research interests include information countermeasure in wireless systems, information security, terminal information resistance and information security in wireless communication.