Next Article in Journal
NiO Nanoparticles for Electrochemical Insulin Detection
Next Article in Special Issue
Stochastic Memristive Interface for Neural Signal Processing
Previous Article in Journal
Voltammetric E-Tongue for Honey Adulteration Detection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Spectrum Correction Algorithm Based on Beat Signal of FMCW Laser Ranging System

1
Key Laboratory of Biomimetic Robots and Systems (Ministry of Education), Beijing Institute of Technology, Beijing 100081, China
2
Instrument of Science and Technology, Beijing Information Science and Technology University, Beijing 100192, China
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(15), 5057; https://doi.org/10.3390/s21155057
Submission received: 23 June 2021 / Revised: 21 July 2021 / Accepted: 23 July 2021 / Published: 26 July 2021
(This article belongs to the Special Issue Data, Signal and Image Processing and Applications in Sensors II)

Abstract

:
The accuracy of target distance obtained by a frequency modulated continuous wave (FMCW) laser ranging system is often affected by factors such as white Gaussian noise (WGN), spectrum leakage, and the picket fence effect. There are some traditional spectrum correction algorithms to solve the problem above, but the results are unsatisfactory. In this article, a decomposition filtering-based dual-window correction (DFBDWC) algorithm is proposed to alleviate the problem caused by these factors. This algorithm reduces the influence of these factors by utilizing a decomposition filtering, dual-window in time domain and two phase values of spectral peak in the frequency domain, respectively. With the comparison of DFBDWC and these traditional algorithms in simulation and experiment on a built platform, the results show a superior performance of DFBDWC based on this platform. The maximum absolute error of target distance calculated by this algorithm is reduced from 0.7937 m of discrete Fourier transform (DFT) algorithm to 0.0407 m, which is the best among all mentioned spectrum correction algorithms. A high performance FMCW laser ranging system can be realized with the proposed algorithm, which has attractive potential in a wide scope of applications.

1. Introduction

The frequency modulated continuous wave (FMCW) laser ranging system is a noncontact detecting and distance measurement system, which has a large detection range and high measurement accuracy and has been widely used in high precision ranging. The system utilizes the corresponding relationship between frequency and distance, which means that the accuracy of distance value relies on the frequency resolution of beat signal obtained by a series of processing with the emitted signal and echo signal. Therefore, the key of ranging lies in the frequency calculation of beat signal [1].
The frequency value of the beat signal can be computed by discrete Fourier transform (DFT) after the signal is sampled and digitized. Ideally, the frequency resolution of it is closely related to the number of sampling points. Too few sampling points will decrease the frequency resolution and lead to the picket fence effect, which affects ranging accuracy, while too many will increase the computing time and the complexity of signal processing. An optimization method is to add points whose value are zero after the sampled beat signal [2]. The accuracy of this method for calculating the frequency completely depends on the number of added points. However, this operation is equivalent to utilizing a rectangular window function on the beat signal in the time domain. This not only cannot change the width of the main lobe in the spectrum but also causes spectrum leakage to a certain extent. Consequently, spectrum correction algorithms to improve the accuracy and resolution of the beat signal frequency have become more significant.
The ratio algorithm is the interpolation-based correction method [3,4,5]. Agrez et al. [6] and Belega et al. [7,8,9] have conducted a further studies on it and proposed new methods based on it to reduce the influence of spectrum leakage on the accuracy of correction. However, the above methods are all at the expense of noise adaptability. The phase difference (PD) algorithm originally comes from a phase interpolation estimator of a single tone frequency in noise proposed by McMahon et al. [10]. Zhu et al. [11] and Kang et al. [12] have performed further research on it, which indicates that the PD algorithm provides superior accuracy in frequency estimates compared with the ratio algorithm and has good adaptability. Luo et al. [13] proposed a new PD method based on asymmetric windows, which can be used to correct the errors of frequency. The main advantages of this algorithm are its characteristics of simple application and strong anti-noise performance, but its reduction of spectrum leakage is unsatisfactory.
The concept of energy centrobaric correction (ECC) algorithm [14] is originally proposed by Offelli et al. [15]. Many researchers have investigated the interferences from spectral components, wideband noise, and other precision factors related to the estimated parameters [16,17]. This algorithm has fast speed and great accuracy of frequency calculation, so it has been applied to engineering after improvement [18]. However, the correction accuracy is too dependent on the symmetric window function and is easily affected by white Gaussian noise (WGN). The Chirp z-transform (CZT) algorithm is a z-transformation method [19]. Because of the low complexity of calculation and the high correction accuracy, there are a lot of CZT-based related applications [20,21,22,23,24,25]. Because this algorithm still analyzes the truncated signal, it only reduces the influence of the picket fence effect on the local spectrum, whereas it does not significantly solve the problem caused by spectrum leakage. The Zoom fast Fourier transform (ZFFT) algorithm achieves spectrum correction by reducing the sampling rate of the signal. It blends complex down-conversion, low-pass filtering, and sample-rate change by way of decimation, thereby improving the frequency resolution [26]. Al-Qudsi et al. [27] presented an implementation method of the ZFFT approach to estimate the spectral peak in the FMCW radar, utilizing a field programmable gate array (FPGA). This algorithm can decrease the complexity of calculations and alleviate the influence of picket fence effect. However, it is severely affected by spectrum leakage and WGN.
In view of the unsatisfactory accuracy and resolution of beat signal frequency affected by WGN, spectrum leakage, and picket fence effect, which cannot be solved by the traditional algorithms above, we propose a new spectrum correction algorithm called decomposition filtering-based dual-window correction (DFBDWC). The main contributions are as follows:
(1)
This algorithm reduces the influence of WGN, affecting the correction accuracy. In the decomposition and filter part, the beat signal is divided into several components, and each component has its characteristics in the frequency domain. Among them, the first few components possess the widest frequency coverage, and there are no obvious peaks in their power spectrum. The sum can be used as the input of the filter, and the WGN in the beat signal will be mostly removed with the weight parameter.
(2)
This algorithm minimizes the impact of spectrum leakage effectively. The Hann window has a narrow main lobe, low side lobe, and fast attenuation speed from the main lobe to the first side lobe. Using two Hann windows in the correction part can concentrate more energy of the signal, thereby making the spectral peak of the desired frequency more obvious.
(3)
This algorithm diminishes the picket fence effect that may decrease the frequency resolution of the beat signal. We utilize phase values and the delay value of two signals in the frequency domain after DFT processing. The phase values correspond to the spectral peaks that are at the same position in these signals. Therefore, the calculation error caused by broad adjacent spectral lines near the peak in only one used signal is avoided, and an accurate frequency value of the beat signal is obtained.
(4)
This algorithm is different from the traditional spectrum correction algorithm, which can reduce the influence caused by WGN, spectrum leakage, and the picket fence effect at the same time, so that the frequency value obtained by this algorithm is more accurate and the distance ranged by this system is more precise.
This article is organized as follows. In Section 2, the principle of the FMCW laser ranging system is firstly briefly introduced, and we explain the DFBDWC algorithm in detail. In Section 3, we built an experimental platform based on the principle of the FMCW laser-ranging system. The results are obtained via simulation and experiment on this platform. Afterwards, the discussion that evaluates the spectrum correction performance of this algorithm by comparing it with these traditional algorithms is conducted. Finally, Section 4 concludes the article.

2. Methods

Using the method shown in Figure 1, we can obtain the high-precision distance value of the target.
The FMCW laser-ranging system emits a modulated laser signal that is reflected by the target and received by the system. After the processing of the received laser signal, the system will output the sampled beat signal. In the software part, we can calculate the precise frequency value of the sampled beat signal with the DFBDWC algorithm and obtain the distance value of this target by taking the frequency value into the equation. In this section, we will introduce the principle of the FMCW laser ranging system and the DFBDWC algorithm, respectively, in detail.

2.1. FMCW Laser Ranging System

The FMCW laser ranging system can be mainly divided into seven parts. The schematic diagram of it is as shown in Figure 2. The signal processing part controls the signal emitting part to generate the FMCW emitted signal, and it drives the laser diode to emit a linear beam, which is the emitted laser signal. The avalanche photo diode (APD) receives the laser signal that is focused by the lens and outputs the echo signal, which is a FMCW signal with a certain delay of emitted signal. The echo signal and the local oscillator signal synchronized by the signal emitting part are mixed in the signal mixing part, and with a series of processing, the beat signal is obtained. In the signal processing part, the beat signal is digitized and transformed into data, which are stored and sent to the PC. Finally, the beat signal is analyzed and processed by the algorithm in the PC, and the distance is computed.
In this FMCW laser ranging system, the frequency of the emitted signal is modulated by a triangle wave, which is as shown in Figure 3. Because of the static ranging target, the effect of the Doppler shift does not have to be considered. Then, the emitted signal s T ( t ) is expressed by
s T ( t ) = A 0 cos ( 2 π f 0 t + π k t 2 + φ 0 ) ,
where A 0 is the amplitude of emitted signal, f 0 is the initial frequency, φ 0 is the initial phase, k = 2 B / T is the modulation slope, B is the modulation bandwidth, and T is the modulation period. With the delay τ = 2 R / c , the echo signal s R ( t ) can be obtained by
s R ( t ) = η A 0 cos ( 2 π f 0 ( t τ ) + π k ( t τ ) 2 + φ 0 ) ,
where η is the amplitude decay rate of echo signal, R is the distance, and c is the speed of light. With the mixing of the emitted and echo signal, the beat signal s B ( t ) can be calculated by
s B ( t ) = η A 0 2 2 cos ( 2 π f 0 τ + 2 π k τ t π k τ 2 ) .
Obviously, the frequency of the beat signal f B = k τ . Thus, the relationship between R and f B is
R = T c 4 B f B .
As shown in the above equation, the factors affecting the accuracy of FMCW laser ranging system are T , B , and f B . Because T and B are the inherent parameters of this system, and they have already reached the limit of system performance, they do not have a decisive impact on the ranging accuracy. Improving the frequency resolution of the beat signal and obtaining the accurate f B become the most significant work of this system.

2.2. DFBDWC Algorithm

This section depicts a new spectrum correction algorithm that is different from the six traditional algorithms introduced in Section 1. The key step of it, whose process chart is as shown in Figure 4, is as follows: Improving the signal to noise ratio (SNR) of s B ( t ) in decomposition and filter, extracting two sub-signals with a dual-window on each of them, and calculating accurate f B according to the phase values in correction and DFT spectrum analysis. All the parameters are in digital form as ( n ) .
WGN is a kind of noise whose probability density function satisfies the statistical characteristics of normal distribution and whose power spectral density function is constant. The most noteworthy feature of it is that the signal contains all frequency components from negative infinity to positive infinity, so it can be apparently distinguished from useful signals with a spectral peak in the spectrum. Accordingly, a similar sequence of WGN can be decomposed from the beat signal.
First, decomposition is a new method based on empirical mode decomposition (EMD) [28,29,30] that is a direct extraction of the energy associated with various intrinsic time scales and the most important parameters of the system. After processing with EMD, a signal can be expressed as a sum of amplitude- and frequency-modulated functions called modes. Each mode is intrinsic and has unique characteristics in the frequency domain, which means several of them enable the estimation of WGN. However, there are some phenomena in EMD, such as oscillations with very disparate scales in one mode or oscillations with similar scales in different modes. These phenomena will cause a problem called “mode mixing”, and some decomposed modes, strictly speaking, will not be a single component signal, so it is not conducive to estimate the noise component accurately.
In order to alleviate “mode mixing”, we take advantage of the dyadic filter bank behavior of EMD and the addition of WGN that populates the whole time–frequency space. Then, the K sub-signal of s B ( n ) can be expressed by
( s B ( n ) ) k = s B ( n ) + ( 1 ) k · β 1 · E 1 ( G k ( n ) ) ,
where G k ( n )   ( k = 1 , 2 , , K 1 , K ) is the k th added WGN signal of zero mean unit variance, K is the number of WGN signals, E j ( · )   ( j = 1 , 2 , , J 1 , J ) is the j th mode obtained by EMD, and J is the number of modes or components, β j   ( j = 1 , 2 , , J 1 , J ) is the j th parameter used to adjust SNR between added WGN signals and s B ( n ) . The main purpose of adding the WGN signal with known features operated by EMD is to generate new extreme points. Additionally, with the operation of plus and minus, s B ( n ) will be forced to focus on some specific values in the scale energy space.
After that, let M ( · ) be the operator that produces the mean envelope of each signal in parentheses, which is same as the procedure in EMD and will make use of these new extreme points in parentheses. Additionally, let A ( · ) be the operator that produces an average signal of all signals in parentheses. With the operation of EMD, we can obtain the first component s ^ 1 ( n ) of beat signal s B ( n ) , which is
s ^ 1 ( n ) = A ( E 1 ( ( s B ( n ) ) k ) ) = A ( ( s B ( n ) ) k M ( ( s B ( n ) ) k ) ) = s B ( n ) A ( M ( ( s B ( n ) ) k ) ) .
Averaging is meant to better estimate the mean envelope value, which reduces “mode mixing” and produces more distinct components. Because of the concept in EMD called residue, the first residue of the beat signal can be expressed by
R 1 ( n ) = s B ( n ) s ^ 1 ( n ) = A ( M ( ( s B ( n ) ) k ) ) .
With s ^ 1 ( n ) and R 1 ( n ) , we can estimate the second residue R 2 ( n ) and the second component s ^ 2 ( n ) of beat signal, respectively, by
R 2 ( n ) = A ( M ( R 1 ( n ) + ( 1 ) k · β 2 · E 2 ( G k ( n ) ) ) ) ,
s ^ 2 ( n ) = R 1 ( n ) R 2 ( n ) = R 1 ( n ) A ( M ( R 1 ( n ) + ( 1 ) k · β 2 · E 2 ( G k ( n ) ) ) ) .
Similarly, for the j th ( j = 3 , 4 , , J 1 , J ) residue R j ( n ) and the j th ( j = 3 , 4 , , J 1 , J ) component s ^ j ( n ) of the beat signal can be calculated by
R j ( n ) = A ( M ( R j 1 ( n ) + ( 1 ) k · β j · E j ( G k ( n ) ) ) ) ,
s ^ j ( n ) = R j 1 ( n ) R j ( n ) = R j 1 ( n ) A ( M ( R j 1 ( n ) + ( 1 ) k · β j · E j ( G k ( n ) ) ) ) .
In this way, we not only utilize the advantages of EMD to make the frequency distribution of components more obvious but also reduce the effect of “mode mixing” so as to estimate the noise components more accurately.
Next, these J components can be distinguished according to the characteristics of each component in the frequency domain. Among them, the first to j th components possess the widest frequency coverage, and there are no obvious peaks in their power spectrum. At the same time, they are scattered in the time domain. Therefore, the sum of the first to j th components is regarded as the evaluation signal e B ( n ) of s B ( n ) , which can be considered to be the estimation of WGN, and we can put it into the filter.
The filter in this algorithm has three inputs and one output [31], which is expressed by
s ˜ B ( n ) = s B ( n ) w ( n ) e B ( n ) ,
where w ( n ) is a coefficient of weight, s ˜ B ( n ) is reconstruction signal of s B ( n ) . In this filter, e B ( n ) is weighted by a parameter at the same instant, so we consider that it is the possible interference signal. If it is subtracted from s B ( n ) , the useful information can be saved as much as possible in s ˜ B ( n ) . The weight w ( n ) is not a fixed parameter, which needs to be updated with the input e B ( n ) and the output s ˜ B ( n ) at the same instant; then, its value of the next instant will be obtained. In order to ensure the best result of this filter, we consult the calculation method of w ( n ) in [31].
However, these parameters above are described in matrices or vectors according to [31], which will lead to the great cost of increased computational complexity and some stability problems. In order to reduce the complexity of computation and the cost of calculation, we have changed the order of w ( n ) into 1 without affecting the performance of the filter. The expression of the weight update can be expressed by
w ( n ) = w ( n 1 ) + r ( n ) s ˜ B ( n ) e B ( n ) ,
where r ( n ) is a relevant coefficient. According to [31], it can be obtained by
r ( n ) = r ( n 1 ) λ [ 1 g ( n ) e B ( n ) ] ,
where λ is the forgetting factor. It is introduced to give a greater forgetting effect to s ˜ B ( n ) of the latest moment, and give less forgetting effect to s ˜ B ( n ) of the farther moment, so as to ensure that s ˜ B ( n ) in the past period is “forgotten” well, so that the filter can work in a more stable state. g ( n ) is a coefficient of gain, which controls the effect of the output s ˜ B ( n ) . Referring to the steps in [31], g ( n ) is calculated by
g ( n ) = r ( n 1 ) e B ( n ) λ + r ( n 1 ) e B 2 ( n ) .
Afterwards, we utilize correction to process s ˜ B ( n ) without the interference of WGN. Correction is the processing of s ˜ B ( n ) in the time domain [32,33], which is used to reduce the picket fence effect and spectrum leakage by spectral peaks in only two sub-signals of s ˜ B ( n ) with the time delay and dual-window of each sub-signal, respectively. With the processing of this step in the DFBDWC algorithm, the calculated frequency value will be more accurate and precise. Firstly, we extract two series of sub-signals, s ˜ B ( 1 ) ( n ) and s ˜ B ( 2 ) ( n ) , from s ˜ B ( n ) . There are L points of delay between them, which means that s ˜ B ( 2 ) ( n ) is L points behind s ˜ B ( 1 ) ( n ) . By putting the first L points of s ˜ B ( 1 ) ( n ) and the whole points of s ˜ B ( 2 ) ( n ) together, we can acquire s ˜ B ( n ) .
After that, the first correction signal s ˜ 1 ( n ) can be obtained by
s ˜ 1 ( n ) = S ( s ˜ B ( 1 ) ( n ) · N ( W ( n ) W ( n ) ) ) ,
where W ( n ) is a window signal, which is usually the Hann window because of its excellent performance in side lobe suppression, is the operator of convolution, and N ( · ) is the operator of producing normalization signal. S ( · ) is the operator that produces a sum signal of the signal’s front and back halves in parentheses. Figuratively speaking, there is a signal whose length is 2 N in parentheses of S ( · ) ; this signal’s front half means the first to N th points and back half means the N + 1 th to 2 N th points. The sum signal produced by S ( · ) is in the length of N , which is formed by adding the values of the first and N + 1 th, the second and N + 2 th, …, the N th points and 2 N th points.
Similarly, the second correction signal s ˜ 2 ( n ) can be obtained by
s ˜ 2 ( n ) = S ( s ˜ B ( 2 ) ( n ) · N ( W ( n ) W ( n ) ) ) .
With a dual-window in the time domain, the influence of spectrum leakage can be decreased more than the one-window and none-window, that is, the energy is more concentrated in the main lobe of these signals, which is more conducive to the subsequent operation in the frequency domain. Moreover, there cannot be more than two windows applied to these sub-signals of s ˜ B ( n ) , because the mathematical model of correction processing is only in two dimensions.
Finally, in order to make signals turn from time domain into the frequency domain, we process s ˜ 1 ( n ) and s ˜ 2 ( n ) with DFT to obtain their spectrum signals S ˜ 1 ( q ) and S ˜ 2 ( q ) . Before DFT, s ˜ 1 ( n ) and s ˜ 2 ( n ) can be also expressed in exponential form as
s ˜ 1 ( n ) = A e j ( ω B n + θ ) ,
s ˜ 2 ( n ) = A e j [ ω B ( n L ) + θ ] ,
where A is the amplitude of s ˜ B ( n ) , and θ is the initial phase of s ˜ B ( n ) . ω B is the angular frequency of s ˜ B ( n ) , and it can be calculated by
ω B = 2 π f B f S ,
where f S is the sample rate.
After the processing of s ˜ B ( 1 ) ( n ) and s ˜ B ( 2 ) ( n ) with DFT, we can obtain their spectrum signals S ˜ 1 ( q ) and S ˜ 2 ( q ) , respectively, by
S ˜ 1 ( q ) = A e j θ F g 2 ( q Δ ω ω 0 ) ,
S ˜ 2 ( q ) = A e j ( θ ω B L ) F g 2 ( q Δ ω ω 0 ) ,
where F g is the amplitude spectrum of W ( n ) , q is the serial number of spectral lines, Δ ω is the angular frequency between each spectral lines, and ω 0 is the initial angular frequency.
In the amplitude spectrum of S ˜ 1 ( q ) , the corresponding serial number of its spectral peak is q . According to this, we can find the phase values φ 1 ( q ) and φ 2 ( q ) in the phase spectrum of S ˜ 1 ( q ) and S ˜ 2 ( q ) , respectively. These phase values are expressed by
φ 1 ( q ) = θ ,
φ 2 ( q ) = θ ω L .
With Equations (23) and (24), we can only obtain an estimation of ω B as
φ 1 ( q ) φ 2 ( q ) = ω ^ B L .
This is because φ 1 ( q ) φ 2 ( q ) is still different from the ideal value; a compensated value of the phase difference has to be introduced. The corresponding angular frequency at the spectral peak q is 2 π q / I , where I is the length of W ( n ) . After the delay of L , this angular frequency will lead to an additional phase shift of 2 π q L / I , which will increase with this delay. Since the position of spectral peak can be observed, 2 π q L / I is considered to be the compensated value of the phase difference. Then, we will calculate ω B by
φ 1 ( q ) φ 2 ( q ) + 2 π q L I = ω B L .
Eventually, according to Equations (20) and (26), the frequency value f B of s B ( n ) can be calculated by
f B = f S 2 π [ φ 1 ( q ) φ 2 ( q ) L + 2 π q I ] .
The relationship between s ˜ B ( 1 ) ( n ) and s ˜ B ( 2 ) ( n ) with L delay will overcome the error caused by two wide spectral lines. With 2 π q / I , we can compute f B more precisely, and the influence of the picket fence effect will be reduced well. Therefore, the frequency resolution of s B ( n ) can be improved, and the purpose of spectrum correction may be achieved.

3. Results and Discussion

In this section, the performance of DFBDWC algorithm is evaluated by both a simulation and an experiment. In the simulation part, an original signal is constructed with Equation (3). Furthermore, to reach the real situation as much as possible, a WGN signal with appropriate SNR value is added to it, which is regarded as a beat signal. In the experiment part, we built an experimental platform according to the scheme shown in Figure 2, and a beat signal obtained with it is analyzed and processed by this algorithm. Table 1 shows the parameters used in the simulation and the experiment.
Among them, η and A 0 are only used in the simulation. f 0 , B , and T are the key parameters of the FMCW laser ranging system, which influence the theoretical accuracy of distance resolution and are determined by the direct digital synthesizer (DDS) in the experimental platform. Moreover, f S influences the number of sampled points in beat signal and is determined by analog-to-digital converter (ADC) in the experimental platform. f U and f L are considered to be a band-pass filter applied to the beat signal, and they are determined by the performance of the low-pass filter (LPF) in the experimental platform and T , respectively. Because τ cannot be greater than T / 2 , the value of f L is 2 / T ; otherwise, it will violate the principle of the FMCW laser ranging system. According to Equation (4), the distance values that this system can obtain are from 1.5141 to 22.7115   m , which correspond to f L and f U , respectively. Moreover, because the minimum sample time is T / 2 , the maximum range resolution of the system is 1.5141   m , which is given by DFT. In order to ensure the best working states of this platform and verify the performance of this algorithm better, the test distance values are shown in Table 2.

3.1. Simulation

As shown in Figure 5, it is a comparison diagram of s ˜ B ( n ) with s B ( n ) at two specific distances in both the time and frequency domain.
It can be seen that there is a great interference in s B ( n ) . Compared with s B ( n ) , the WGN interference in s ˜ B ( n ) is reduced well, and the wave pattern of s ˜ B ( n ) becomes smoother. Additionally, s ˜ B ( n ) also retains the information of shape, amplitude, and frequency. In their power spectrum, the amplitude of the spectral peak in s ˜ B ( n ) is greater than the amplitude of the spectral peak in s B ( n ) . Moreover, the power spectrum of WGN has also been decreased. The above results indicate that the DFBDWC algorithm can suppress the interference of WGN with a large SNR value in the beat signal and retain the useful information in the signal.
In order to verify the comparison of these signals further, we apply S N R and noise power P n o i s e here. They can be expressed, respectively, by
S N R = 10 lg ( n = 1 N x 2 ( n ) n = 1 N [ x ( n ) s O ( n ) ] 2 ) ,
P n o i s e = 1 N n = 1 N [ x ( n ) s O ( n ) ] 2 ,
where x ( n ) represents s B ( n ) or s ˜ B ( n ) . S N R judges the macroscopic performance of this algorithm. The larger the S N R is, and the smaller P n o i s e is, the better the performance of this algorithm is, that is, the energy will be more concentrated in the position of the spectral peak and the probability of selecting a “fake” peak as the target position due to enormous WGN will be reduced better.
We calculate the value of these parameters at different distances and noise powers, and the results are as shown in Figure 6. Among the first two figures, S N R increases from 10   dB to more than 25   dB , and P n o i s e reduces from 0.13   W to the order of 10 3   W . In the last two figures, S N R increases from 1   dB to more than 11   dB , and P n o i s e reduces from more than 5   W to about 0.1   W . The simulation results above indicate that this algorithm has a great suppression effect on WGN interference, and it saves useful information in s ˜ B ( n ) as much as possible. This will improve the accuracy of f B .
Furthermore, we compare the performance of the DFBDWC algorithm with six other algorithms in terms of computing f B , that is, the computed distances, and the original S N R of WGN in s B ( n ) is 10   dB . To show the results of the comparison better, absolute error ( A E ) and root mean square error ( R M S E ) are applied here. They can be expressed, respectively, by
A E = | R C ( p ) R T ( p ) | ,
R M S E = 1 P p = 1 P [ R C ( p ) R T ( p ) ] 2 ,
where R C ( p ) is the computed distance of the p th test distance, R T ( p ) is the p th test distance, and P is the number of test distance. The smaller the R M S E is, the better the performance of the algorithm is.
The computed distance and calculation results of A E and R M S E can be seen from Table 3 and Figure 7, respectively. At the first test distance, A E of the PD algorithm has the maximum value, and at the last test distance, A E of DFT algorithm has the maximum value. Among these traditional spectrum algorithm, DFT algorithm has the biggest jump of A E , while the ECC algorithm has the smallest jump of A E . The A E of the Ratio, ECC, CZT, and ZFFT algorithm are relatively stable. At each test distance, the A E of the DFBDWC algorithm basically has the minimum value. R M S E macroscopically evaluates the deviation between R C ( p ) and R T ( p ) . It can be seen that the DFT algorithm has the maximum value, and the DFBDWC algorithm has the minimum value, which indicates that the DFBDWC algorithm performs the best compared with other traditional algorithms.
Additionally, in order to illustrate that this algorithm can reduce the influence of the picket fence effect and spectrum leakage further, we conduct a simulation with s B ( n ) without WGN, that is s O ( n ) , at different distances, whose results of A E are as shown in Table 4. The DFBDWC algorithm still basically has the most minimum value of A E among all the spectrum correction algorithms. According to [34,35], the picket fence effect and spectrum leakage significantly decrease the precision of DFT when applying asynchronous sampling in practical applications. Additionally, there are disadvantages of each traditional spectrum correction algorithm that are described in Section 1. Therefore, the distance calculation accuracy of the DFBDWC algorithm is much better than any other six algorithms when processing with s O ( n ) , which can prove our new spectrum correction algorithm not only decreases the influence of spectrum leakage, but also reduces the picket fence effect.
The simulation results demonstrate that the DFT algorithm cannot accurately obtain the distance value, since it cannot overcome the problems described in Section 1. Although the other five algorithms have improved the accuracy of R C ( p ) compared with the DFT algorithm and they have achieved a certain spectrum correction effect, they are still inferior to the performance of the DFBDWC algorithm. As a consequence, our algorithm will improve the accuracy of distance calculation.

3.2. Experiment

The experimental platform and scene in the experiment part are as shown in Figure 8.
The laser diode whose power is 500   mW and wavelength is 950   nm is driven by an emitted signal generated by an emitted signal (ES) DDS named AD9958. The laser beam is reflected on the target surface at R T ( p ) and focused by the lens; then, the echo signal is outputted by an APD with 16 linearly arrayed receiving units. In each signal mixing part, two series of echo signals can be processed. However, the echo signal is too weak and needs to be amplified to a certain amplitude by an amplifier named AD8001, after which a local oscillator (LO) DDS synchronized by ES DDS is mixed with it in a mixer named AD831 to form a mixed signal that contains two frequency values because of the working principle in the mixer. The large frequency value, which is an interference, needs to be filtered out by an LPF named MAX274 whose bandwidth is 150   kHz , and the small frequency value passed through an automatic gain control (AGC) named AD8367 is amplified. Then, the beat signal of appropriate amplitude can be obtained. In the signal processing part, the beat signal is sampled by an ADC named AD9253. Finally, the data are transferred to a FPGA named XC7Z100 for temporary storage and transmitted to the PC for the distance calculation using the DFBDWC algorithm. The experimental platform is placed in a corridor with a length of 15   m , and the distance between the target and APD is considered to be R T ( p ) , whose value is shown in Table 2. In order to place the target in a precise position, a tape measure with centimeter accuracy is used specially, and its starting position is the surface of APD. At the same time, two benchmarks are placed at 2 and 10   m , respectively, to indicate the placement range of the target.
Above all, we utilize only one of the units in APD to receive the laser signal from a plane target. With the analysis and processing of the DFBDWC algorithm, we obtain a comparison diagram of s ˜ B ( n ) and s B ( n ) at two specific distances in the time domain, which is as shown in Figure 9. It can be seen that they are different from s ˜ B ( n ) and s B ( n ) in the simulation part, but they still contain the frequency information of s B ( n ) . Compared with s B ( n ) , the WGN interference in s ˜ B ( n ) is reduced well, and the wave pattern of s ˜ B ( n ) becomes smoother. This indicates that the algorithm can suppress the interference of noise in the beat signal and retain the useful information in the signal.
Similarly, we compare the performance of the DFBDWC algorithm with other six algorithms in terms of computed distances in Table 5 and error analysis in Figure 10. We can note from Figure 10 that as for the other six algorithms, the A E of DFT and CZT algorithm has the maximum value, respectively, at the first and the last test distance. Overall, the DFT algorithm has the biggest jump of A E , while the ECC algorithm has the smallest jump of A E . At each test distance, the A E and R M S E of the DFBDWC algorithm basically have the minimum value. This indicates that the performance of this algorithm is the best in all these algorithms. Additionally, the maximum and minimum A E of each algorithm are listed in Table 6. The maximum A E is decreased from 0.7937 to 0.0407   m by using the DFBDWC algorithm. As expected, this algorithm overcomes the problem to a certain extent caused by spectrum leakage and the picket fence effect and improves the accuracy of distance calculation.
On the basis of experimental results above, we can conclude that because of decomposition and filter, WGN has been estimated accurately from the beat signal and reduced to a certain extent, which will improve the SNR of the beat signal. Additionally, due to the dual-window applied in correction, the energy is more concentrated and the influence of spectrum leakage has been decreased. Moreover, utilizing two main spectral lines at the peaks of these sub-signals with delay to calculate the frequency of the beat signal has alleviated the problem of poor accuracy of results caused by the picket fence effect, and it also avoid the interference that arises from using multiple spectral lines in some traditional spectrum correction algorithms. When the beat signal is not affected by WGN severely, correction will become the key step that makes the distance calculation more accurate in the DFBDWC algorithm.
In addition to the accuracy comparison of computed distance, we also calculate computation time consumed for these seven algorithms by processing the same group of sampled beat signals so as to judge the efficiency of each algorithm. Using an PC with CPU of Intel i7-7700 and RAM of 16 GB, we obtain the results that are as shown in Table 7. It can be seen that in different sample times, the computation time consumption of the ZFFT algorithm is the least, the DFBDWC algorithm is the most, and the other algorithms are almost the same. Additionally, with the increase in sample time, the computation time consumption of the DFBDWC algorithm is doubled, and there are not many rises in the other algorithms. This is because decomposition and filter in this algorithm have to process by iterative operation. The larger the sampled points of the beat signal, the more computation time consumption will be needed. In practical application, we only focus on the accuracy of the distance calculation, while we do not require any real-time computation. Therefore, we sacrifice the efficiency for the precision of our algorithm.
It can be found from the experimental results that the performance of each algorithm is consistent with simulation results. This indicates that simulation has achieved the real situation well, and the parameters for evaluating the performance of these algorithms is also reasonable. However, the values of each parameter obtained in experiment are larger than those in simulation, which is caused by errors from the experimental platform, the factors of the environment and the target placement, such as the sensitivity of APD, the response speed of the laser diode, the bandwidth of the emitted signal, the interference of ambient light, and the accuracy of distance and angle when we place the target. This can still verify that the DFBDWC algorithm reduces the influence of WGN, spectrum leakage, and the picket fence effect. Moreover, it performs the best among the existing spectrum correction algorithms, and the maximum A E of it is not more than 0.05   m .

4. Conclusions

In this article, we proposed a new spectrum correction algorithm named DFBDWC, and built an experimental platform based on the principle of the FMCW laser ranging system. The experimental platform outputs the data of the beat signal, which is analyzed and processed by the DFBDWC algorithm in the PC, and the target distance detected by the system is obtained. Comparing this algorithm to the traditional DFT algorithm and other spectrum correction algorithms in both the simulation and the experiment, including the PD, ECC, Ratio, CZT, and ECC algorithm, we achieve the performance evaluation of this algorithm. The results indicate that DFBDWC algorithm can reduce the influence of WGN, spectrum leakage, and the picket fence effect. Additionally, it can also improve the accuracy and frequency resolution of the beat signal. The maximum absolute error of the target distance calculated by this algorithm is reduced from 0.7937 m of the DFT algorithm to 0.0407 m, which is the best among all the spectrum correction algorithms. The most remarkable performance of our algorithm is because decomposition can estimate WGN accurately in the beat signal and the filter reduces it to a certain extent. The double Hann window applied in correction concentrates more energy in the spectrum, which minimizes the impact of spectrum leakage well. Utilizing two main spectral lines at the peaks of these sub-signals with a delay to calculate the frequency of the beat signal has alleviated the problem of poor accuracy of results caused by the picket fence effect, and it also avoids the interference that arises from using multiple spectral lines in some traditional spectrum correction algorithms. Therefore, the DFBDWC algorithm can improve the performance of the FMCW laser ranging system. In future work, it is necessary to upgrade this platform of the system, such as by choosing more sensitive APD, selecting a laser diode with a faster response speed, and increasing the bandwidth of the emitted signal, which makes it adapt to this algorithm better. In addition, we still need to optimize the structure and computational complexity of our algorithm so that the efficiency of distance calculation in engineering can be greatly raised while keeping the accuracy. Furthermore, we will carry out experiments by utilizing 16 units in APD to figure out the surface fitting uncertainty for different object shapes of this system so that it could make our algorithm and platform more valuable for 3D scanning.

Author Contributions

Y.H. and P.S. proposed the decomposition filtering based dual-window correction algorithm for FMCW laser ranging system. Y.H., X.W., and Z.P. built the experimental platform and analyzed the simulation and experimental results. All authors have read and agreed to the published version of the manuscript.

Funding

The research was funded by Beijing Advanced Innovation Center for Intelligent Robots and Systems, grant number 2019IRS13.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Amann, M.C.; Bosch, T.; Lescure, M.; Myllyla, R.; Rioux, M. Laser ranging: A critical review of usual techniques for distance measurement. Opt. Eng. 2001, 40, 10–19. [Google Scholar]
  2. Borkowski, J.; Mroczka, J. LIDFT method with classic data windows and zero padding in multifrequency signal analysis. Measurement 2010, 43, 1595–1602. [Google Scholar] [CrossRef]
  3. Rife, D.C.; Vincent, G.A. Use of the discrete Fourier transform in the measurement of frequencies and levels of tones. Bell Syst. Tech. J. 1970, 49, 197–228. [Google Scholar] [CrossRef]
  4. Grandke, T. Interpolation algorithms for discrete Fourier transforms of weighted signals. IEEE Trans. Instrum. Meas. 1983, 32, 350–355. [Google Scholar] [CrossRef]
  5. Ming, X.; Kang, D. Corrections for frequency, amplitude and phase in a fast Fourier transform of a harmonic signal. Mech. Syst. Signal Proc. 1996, 10, 211–221. [Google Scholar] [CrossRef]
  6. Agrez, D. Weighted multipoint interpolated DFT to improve amplitude estimation of multifrequency signal. IEEE Trans. Instrum. Meas. 2002, 51, 287–292. [Google Scholar] [CrossRef]
  7. Belega, D.; Dallet, D. Frequency estimation via weighted multipoint interpolated DFT. IET Sci. Meas. Technol. 2008, 2, 1–8. [Google Scholar] [CrossRef]
  8. Belega, D.; Dallet, D. Multifrequency signal analysis by interpolated DFT method with maximum sidelobe decay windows. Measurement 2009, 42, 420–426. [Google Scholar] [CrossRef]
  9. Belega, D.; Dallet, D.; Petri, D. Accuracy of sine wave frequency estimation by multipoint interpolated DFT approach. IEEE Trans. Instrum. Meas. 2010, 59, 2808–2815. [Google Scholar] [CrossRef]
  10. Mcmahon, D.R.A.; Barrett, R.F. An efficient method for the estimation of the frequency of a single tone in noise from the phases of discrete Fourier transforms. Signal Process. 1986, 11, 169–177. [Google Scholar] [CrossRef]
  11. Zhu, L.; Li, H.; Ding, H.; Xiong, Y. Noise influence on estimation of signal parameter from the phase difference of discrete Fourier transforms. Mech. Syst. Signal Proc. 2002, 16, 991–1004. [Google Scholar] [CrossRef]
  12. Kang, D.; Ming, X.; Xiaofei, Z. Phase difference correction method for phase and frequency in spectral analysis. Mech. Syst. Signal Proc. 2000, 14, 835–843. [Google Scholar] [CrossRef]
  13. Luo, J.; Xie, M. Phase difference methods based on asymmetric windows. Mech. Syst. Signal Proc. 2015, 54/55, 52–67. [Google Scholar] [CrossRef]
  14. Ding, K.; Cao, D.; Li, W. An approach to discrete spectrum correction based on energy centroid. Key Eng. Mater. 2006, 321-323, 1270–1273. [Google Scholar]
  15. Offelli, C.; Petri, D. A frequency-domain procedure for accurate real-time signal parameter measurement. IEEE Trans. Instrum. Meas. 1990, 39, 363–368. [Google Scholar] [CrossRef]
  16. Belega, D.; Dallet, D.; Petri, D. Accuracy of the normalized frequency estimation of a discrete-time sine-wave by the energy-based method. IEEE Trans. Instrum. Meas. 2012, 61, 111–121. [Google Scholar] [CrossRef]
  17. Lin, H.; Ding, K. Energy based signal parameter estimation method and a comparative study of different frequency estimators. Mech. Syst. Signal Proc. 2011, 25, 452–464. [Google Scholar]
  18. Zhang, Q.; Zhong, S.; Lin, J.; Huang, Y.; Nsengiyumva, W.; Chen, W.; Luo, M.; Zhong, J.; Yu, Y.; Peng, Z.; et al. Anti-noise frequency estimation performance of Hanning-windowed energy centrobaric method for optical coherence velocimeter. Opt. Lasers Eng. 2020, 134. [Google Scholar] [CrossRef]
  19. Rabiner, L.; Schafer, R.; Rader, C. The chirp z-transform algorithm. IEEE Trans. Audio Electroacoust. 1969, 17, 86–92. [Google Scholar] [CrossRef]
  20. Leng, J.; Shan, C. Application of chirp-z transformation in high accuracy measurement of radar. Appl. Mech. Mater. 2013, 392, 730–733. [Google Scholar] [CrossRef]
  21. Li, D.; Liu, H.; Liao, Y.; Gui, X. A novel helicopter-borne rotating SAR imaging model and algorithm based on inverse chirp-z transform using frequency-modulated continuous wave. IEEE Geosci. Remote Sens. Lett. 2015, 12, 1625–1629. [Google Scholar]
  22. Masso, E.; Bolognini, N. Dynamic multiple-image encryption based on chirp z-transform. J. Opt. 2019, 21, 035704. [Google Scholar] [CrossRef]
  23. Qin, M.; Li, D.; Tang, X.; Zeng, C.; Li, W.; Xu, L. A fast high-resolution imaging algorithm for helicopter-borne rotating array SAR based on 2-D chirp-z transform. Remote Sens. 2019, 11, 1669. [Google Scholar] [CrossRef] [Green Version]
  24. Shen, S.; Nie, X.; Tang, L.; Bai, Y.; Zhang, X.; Li, L.; Ben, D. An improved coherent integration method for wideband radar based on two-dimensional frequency correction. Electronics 2020, 9, 840. [Google Scholar] [CrossRef]
  25. Wei, D.; Nagata, Y.; Aketagawa, M. Partial phase reconstruction for zero optical path difference determination using a chirp z-transform-based algorithm. Opt. Commun. 2020, 463, 125456. [Google Scholar] [CrossRef]
  26. Lyons, R.G. The Zoom FFT. In Understanding Digital Signal Processing, 3rd ed.; Bernard, G., Ed.; Pearson Education: Boston, MA, USA, 2011; pp. 548–550. [Google Scholar]
  27. Al-Qudsi, B.; Joram, N.; Strobel, A.; Ellinger, F. Zoom FFT for precise spectrum calculation in FMCW radar using FPGA. In Proceedings of the 2013 9th Conference on Ph.D. Research in Microelectronics and Electronics, Villach, Austria, 24–27 June 2013. [Google Scholar]
  28. Huang, N.E.; Shen, Z.; Long, S.R.; Wu, M.C.; Shih, H.H.; Zheng, Q.; Yen, N.C.; Tung, C.C.; Liu, H.H. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis. Proc. R. Soc. A Math. Phys. Eng. Sci. 1998, 454, 903–995. [Google Scholar] [CrossRef]
  29. Torres, M.E.; Colominas, M.A.; Schlotthauer, G.; Flandrin, P. A complete ensemble empirical mode decomposition with adaptive noise. In Proceedings of the IEEE International Conference on Acoustics, Speech, and Signal Processing (ICASSP), Prague Congress Center, Prague, Czech Republic, 22–27 May 2011. [Google Scholar]
  30. Colominas, M.A.; Schlotthauer, G.; Torres, M.E. Improved complete ensemble EMD: A suitable tool for biomedical signal processing. Biomed. Signal Process. Control. 2014, 14, 19–29. [Google Scholar] [CrossRef]
  31. Diniz, P.S.R. Adaptive Filtering: Algorithms and Practical Implementation, 5th ed.; Springer Nature: Cham, Switzerland, 2020; pp. 157–160. [Google Scholar]
  32. Huang, X.; Wang, Z.; Hou, G. New method of estimation of phase, amplitude, and frequency based on all phase FFT spectrum analysis. In Proceedings of the IEEE 2007 International Symposium on Intelligent Signal Processing and Communication Systems, Xiamen, China, 28 November–1 December 2007. [Google Scholar]
  33. Su, T.; Yang, M.; Jin, T.; Flesch, R.C.C. Power harmonic and interharmonic detection method in renewable power based on Nuttall double-window all-phase FFT algorithm. IET Renew. Power Gener. 2018, 12, 953–961. [Google Scholar] [CrossRef]
  34. Zhang, D.; Sun, S.; Zhao, H.; Yang, J. Laser Doppler dignal processing based on trispectral interpolation of Nuttall window. Optik 2019, 205, 163364. [Google Scholar] [CrossRef]
  35. Lin, H.C. Power harmonics and interharmonics measurement using recursive group-harmonic power minimizing algorithm. IEEE Trans. Ind. Electron. 2012, 59, 1184–1193. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of distance obtained by the FMCW laser ranging system and DFBDWC algorithm.
Figure 1. Schematic diagram of distance obtained by the FMCW laser ranging system and DFBDWC algorithm.
Sensors 21 05057 g001
Figure 2. Scheme of FMCW laser ranging system.
Figure 2. Scheme of FMCW laser ranging system.
Sensors 21 05057 g002
Figure 3. Type of FMCW modulation. The blue and red line represent the emitted signal s T ( t ) and echo signal s R ( t ) , respectively.
Figure 3. Type of FMCW modulation. The blue and red line represent the emitted signal s T ( t ) and echo signal s R ( t ) , respectively.
Sensors 21 05057 g003
Figure 4. Process of DFBDWC algorithm. The input and output of this algorithm are the beat signal and frequency value of the beat signal, respectively.
Figure 4. Process of DFBDWC algorithm. The input and output of this algorithm are the beat signal and frequency value of the beat signal, respectively.
Sensors 21 05057 g004
Figure 5. Simulation results of signals comparison at two different distances. (a) Beat signals s B ( n ) obtained by original signals and WGN signals with SNR of 10   dB . (b) Reconstruction signals s ˜ B ( n ) . (c) The power spectrum of s B ( n ) . (d) The power spectrum of s ˜ B ( n ) .
Figure 5. Simulation results of signals comparison at two different distances. (a) Beat signals s B ( n ) obtained by original signals and WGN signals with SNR of 10   dB . (b) Reconstruction signals s ˜ B ( n ) . (c) The power spectrum of s B ( n ) . (d) The power spectrum of s ˜ B ( n ) .
Sensors 21 05057 g005
Figure 6. Simulation results of two parameters for evaluating signal comparison at different distances and noise powers. (a,c) are S N R , (b,d) are P n o i s e . Additionally, the original S N R of WGN in (a,b) is 10   dB , and the original S N R of WGN in (c,d) is 1   dB .
Figure 6. Simulation results of two parameters for evaluating signal comparison at different distances and noise powers. (a,c) are S N R , (b,d) are P n o i s e . Additionally, the original S N R of WGN in (a,b) is 10   dB , and the original S N R of WGN in (c,d) is 1   dB .
Sensors 21 05057 g006
Figure 7. Simulation results of seven algorithms comparison at different distances. There are WGN signals with S N R of 10   dB in s B ( n ) . (a) is A E . (b) is R M S E .
Figure 7. Simulation results of seven algorithms comparison at different distances. There are WGN signals with S N R of 10   dB in s B ( n ) . (a) is A E . (b) is R M S E .
Sensors 21 05057 g007
Figure 8. The experimental platform and scene in the experiment part. (a) Scheme of this experimental platform. (b) Experimental platform in the experimental scene. (c) Details of this experimental platform.
Figure 8. The experimental platform and scene in the experiment part. (a) Scheme of this experimental platform. (b) Experimental platform in the experimental scene. (c) Details of this experimental platform.
Sensors 21 05057 g008aSensors 21 05057 g008b
Figure 9. Experimental results of signals comparison at two different distances in time domain. (a) is beat signals s B ( n ) . (b) is reconstruction signals s ˜ B ( n ) .
Figure 9. Experimental results of signals comparison at two different distances in time domain. (a) is beat signals s B ( n ) . (b) is reconstruction signals s ˜ B ( n ) .
Sensors 21 05057 g009
Figure 10. Experimental results of seven algorithms comparison at different distances. (a) is A E . (b) is R M S E .
Figure 10. Experimental results of seven algorithms comparison at different distances. (a) is A E . (b) is R M S E .
Sensors 21 05057 g010
Table 1. The parameters used in the simulation and experiment.
Table 1. The parameters used in the simulation and experiment.
ParameterInterpretationValue
η The amplitude decay rate of echo signal0.8
A 0 The amplitude of emitted signal1
f 0 The initial frequency1 MHz
c The speed of light299,792,458 m/s
B The modulation bandwidth99 MHz
T The modulation period200 μs
f S The sample rate20 MHz
f U The upper limit frequency0.8
f L The lower limit frequency1
Table 2. The test distance value used in the simulation and experiment.
Table 2. The test distance value used in the simulation and experiment.
Lower Limit (m)Test Distance (m)Upper Limit (m)
1.5141234567891022.7115
2.53.54.55.56.57.58.59.5
Table 3. Simulation results of computed distances comparison between seven algorithms. There are WGN signals with S N R of 10   dB in s B ( n ) .
Table 3. Simulation results of computed distances comparison between seven algorithms. There are WGN signals with S N R of 10   dB in s B ( n ) .
Test Distance
(Real Distance) (m)
Computed Distance
(m)
DFTPDECCRatioCZTZFFTDFBDWC
21.84832.27122.00922.01622.01092.12302.0021
2.52.31032.50132.48982.51022.51362.51212.5011
33.23452.99973.00182.99773.01642.93543.0004
3.53.69653.51023.51633.51183.51913.67943.5038
44.15863.99893.99924.00743.99234.06854.0003
4.54.62074.49454.49674.49394.52464.49194.4981
55.08274.98365.01285.00644.99774.88104.9988
5.55.54485.49795.49775.50695.50055.62505.5021
66.00695.99385.99685.99386.00326.01416.0028
6.56.46896.50566.51196.50536.50596.43746.4969
76.93106.99306.99427.00277.00866.82657.0017
7.57.39317.49547.49877.49687.51147.57057.5032
87.85527.99568.02028.00798.01417.99398.0022
8.58.31728.50208.50158.50798.51688.38308.4985
98.77938.99489.00088.99439.01969.12709.0030
9.59.70349.50659.52069.50439.49279.51619.4996
1010.16559.999710.003810.00569.99549.939410.0027
Table 4. Simulation results of A E comparison between seven algorithms. There are no WGN signals in s B ( n ) .
Table 4. Simulation results of A E comparison between seven algorithms. There are no WGN signals in s B ( n ) .
Test Distance
(Real Distance) (m)
AE
(cm)
DFTPDECCRatioCZTZFFTDFBDWC
215.172927.11550.75931.60491.091812.29820.0064
2.518.96620.36120.51001.21721.36481.20920.0066
323.44740.12850.00440.38791.63786.45570.0118
3.519.65410.49050.96530.92471.910717.94200.0069
415.86090.08430.42510.74880.77356.85300.0112
4.512.06770.18370.00440.36302.45670.81190.0022
58.27440.06041.12220.64950.227611.90090.0109
5.54.48120.19570.34700.52980.045312.49690.0273
60.68800.24230.00840.34790.31831.40780.0087
6.53.10530.20891.28660.50070.59136.25710.0013
76.89850.26960.27930.40130.864217.34610.0036
7.510.69170.30160.01560.33631.13727.05170.0027
814.48500.29621.46370.40751.41020.61320.0093
8.518.27820.33480.22100.31571.274111.70230.0044
922.07140.36130.02660.32651.956112.69550.0084
9.520.34210.36721.65500.34350.72821.60650.0054
1016.54890.39710.17140.25390.45526.05840.0024
Table 5. Experimental results of computed distances comparison between seven algorithms. There are WGN signals with S N R of 10   dB in s B ( n ) .
Table 5. Experimental results of computed distances comparison between seven algorithms. There are WGN signals with S N R of 10   dB in s B ( n ) .
Test Distance
(Real Distance) (m)
Computed Distance
(m)
DFTPDECCRatioCZTZFFTDFBDWC
22.24062.03782.02882.01422.01122.04051.9947
2.52.24282.38212.39952.43052.43242.39722.5130
33.07783.27253.15693.24133.25703.23443.0044
3.53.45253.27253.44283.30913.30913.27573.4788
44.79374.19133.93454.14484.14204.38194.0037
4.54.53974.93404.52864.89354.89724.71894.4921
54.79354.93405.07974.89354.90864.77025.0177
5.54.79514.93405.47754.94864.91705.01595.5240
66.08486.04236.08346.15546.16616.31215.9786
6.56.79016.72416.26676.80376.79906.80236.5111
76.78566.72417.07006.77726.79126.67986.9593
7.57.56087.63507.67547.52137.53077.65267.5062
87.82347.98077.96557.86417.79927.64658.0192
8.58.38998.71368.34748.36708.39598.57498.5271
99.17528.71369.04209.17089.20799.01468.9694
9.59.72209.75429.75079.72759.71929.74599.4933
109.73249.75429.75079.73749.71639.736710.0143
Table 6. Maximum and minimum A E of each algorithm.
Table 6. Maximum and minimum A E of each algorithm.
DFTPDECCRatioCZTZFFTDFBDWC
Maximum A E (m)0.79370.56600.25070.55140.58300.48410.0407
Minimum A E (m)0.03970.01930.02250.01420.01120.01460.0037
Table 7. The comparison of average computation time consuming between seven algorithms.
Table 7. The comparison of average computation time consuming between seven algorithms.
Sample Time
(μs)
Computation Time Consuming (s)
DFTPDECCRatioCZTZFFTDFBDWC
1000.04230.04490.04580.04520.04410.02432.3891
2000.04440.05010.04680.04640.04600.02595.6616
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hao, Y.; Song, P.; Wang, X.; Pan, Z. A Spectrum Correction Algorithm Based on Beat Signal of FMCW Laser Ranging System. Sensors 2021, 21, 5057. https://doi.org/10.3390/s21155057

AMA Style

Hao Y, Song P, Wang X, Pan Z. A Spectrum Correction Algorithm Based on Beat Signal of FMCW Laser Ranging System. Sensors. 2021; 21(15):5057. https://doi.org/10.3390/s21155057

Chicago/Turabian Style

Hao, Yi, Ping Song, Xuanquan Wang, and Zhikang Pan. 2021. "A Spectrum Correction Algorithm Based on Beat Signal of FMCW Laser Ranging System" Sensors 21, no. 15: 5057. https://doi.org/10.3390/s21155057

APA Style

Hao, Y., Song, P., Wang, X., & Pan, Z. (2021). A Spectrum Correction Algorithm Based on Beat Signal of FMCW Laser Ranging System. Sensors, 21(15), 5057. https://doi.org/10.3390/s21155057

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop