Next Article in Journal
Parkinson’s Disease EMG Data Augmentation and Simulation with DCGANs and Style Transfer
Next Article in Special Issue
Laboratory Testing of FBGs for Pipeline Monitoring
Previous Article in Journal
Compensation of Rotary Encoders Using Fourier Expansion-Back Propagation Neural Network Optimized by Genetic Algorithm
Previous Article in Special Issue
Three-Dimensional Imaging by Frequency-Comb Spectral Interferometry
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

FMCW Laser Fuze Multiple Scattering Model and Accurate Fixed-Distance Algorithm in a Smoke Environment

School of Mechatronical Engineering, Beijing Institute of Technology, Beijing 100081, China
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(9), 2604; https://doi.org/10.3390/s20092604
Submission received: 7 April 2020 / Revised: 26 April 2020 / Accepted: 28 April 2020 / Published: 3 May 2020
(This article belongs to the Special Issue Development and Applications of Optical Sensing)

Abstract

:
In a smoke environment, suspended particles can scatter and absorb laser photons, making target echo signals extremely weak and difficult to extract and identify, which causes obvious difficulty in fixed-distance of laser fuze. In this paper, the multiple scattering model of frequency-modulated-continuous-wave (FMCW) laser fuze in a smoke environment was established. This model simulates multi-path propagation and multiple scattering of photons. At the same time, we use the correntropy spectral density (CSD) algorithm for accurate fixed-distance of FMCW laser fuze. The absolute error of distance does not exceed 0.15 m in smoke interference environment.

1. Introduction

Fixed-distance fuze is a device that measures the distance to the target by radio, laser, ultra-wideband, millimeter wave, and detonates the warhead at the best explosive point. Pulse laser fuze has been widely used for its characteristics of excellent anti-electromagnetic interference capability, simple structure, and low cost. However, pulse laser fuze is susceptible to interference from suspended particles such as smoke and clouds, which makes it difficult to extract and identify target signals [1,2]. Compared with pulse laser, the frequency modulated continuous wave (FMCW) laser ranging system depends on the beat frequency instead of the magnitude and time of the echo signal so that it has a natural advantage in target detection under dense fog conditions [3,4]. Therefore, it is very meaningful to improve the anti-interference ability of FMCW laser fuze in smoke environments.
Research has been done on the target echo characteristics of laser in the interference environment of suspended particles. Grabner et al. studied the effect of extinction coefficients of different concentrations of fog and haze on laser transmission [5]. Lundin et al. studied the transmission characteristics of FMCW laser signals in solid scattering media [6]. Zhang et al. studied the characteristics of smoke backscatter interference signals of FMCW laser fuze, and proposed a method for detecting intermediate frequency signals based on a normalized threshold [7,8]. In our earlier work, a Monte Carlo simulation model of FMCW laser transmission in a smoke interference environment was established [9]. The parameters of the aerosol scattering coefficient and extinction coefficient change randomly during the dynamic approach. In addition, some traditional denoising methods, including empirical mode decomposition (EMD) direct wavelet thresholding, EMD interval thresholding, correlation-based EMD partial reconstruction, discrete Fourier transform (DFT), and wavelet transformation, are investigated to provide a comparison with the cross-power-spectral-density (CPSD) algorithm [10]. Both the simulation and experiment results showed that the CPSD algorithm has better performance than other algorithms for the anti-interference of FMCW laser detection. However, the CPSD system is based on four photodetectors, and the complex and expensive disadvantages limit its application in FMCW laser fuze. This paper provides a new CSD algorithm for fixed-distance. In addition, the accuracy achieved by CPSD using four detectors can be carried out with only one detector.
Corretropy is used to compare the degree of similarity between the two sets of signals. It was first proposed by Jose C. Principe, a professor of electrical and computer engineering at the University of Florida, and analyzed in detail [11,12]. Garde et al. conducted systematic theoretical analysis, applied research on correntropy, and defined the Fourier transform of correntropy as the Corretropy Spectral Density (CSD) function [13,14,15,16,17]. The correntropy performs well in signal processing, image processing, artificial intelligence, etc. Chen et al. designed an adaptive filter based on the generalized maximum correntropy criterion to maximize the correlation between the input signal and the output signal, thereby effectively suppressing noise [18]. Yu et al. proposed a time delay estimation algorithm based on the weighted correntropy spectral density and this method was accurate in the presence of Gaussian or impulsive noise with a low generalized signal-to-noise ratio [19]. In view of the advantages of CSD algorithm in determining the similarity of two sets of signals, this paper chooses the CSD algorithm to process an intermediate frequency signal.
In this paper, we established a laser multipath propagation model in smoke interference environment in Section 2. The process of using the CSD algorithm to determine the distance was described in Section 3. In Section 4, the feasibility and performance of the CSD algorithm were evaluated and compared with FFT algorithm by MATLAB simulation (Matlab R2019a, MathWorks, Natick, MA, USA) in various environments. In Section 5, we established the system of fixed-distance of FMCW laser fuze and carried out the experiment. The experimental results were consistent with the simulation results. Finally, the conclusions are brought in Section 6.

2. FMCW Laser Multiple Scattering Model in the Smoke Interference Environment

The FMCW laser ranging system modulates the light intensity of transmitted laser signal by a signal in which instantaneous frequency varies linearly with time and uses the modulation signal as a local oscillator signal. Due to the time delay of laser transmission in the atmosphere, the instantaneous frequency of the echo signal changes. The echo signal is mixed with the local oscillator signal to obtain the intermediate frequency (IF) signal. When there is no interference, the IF signal should be a cosine signal with a certain frequency. The frequency is related to the target distance and relative speed. Symmetric triangular wave FMCW modulation is widely used because it can directly calculate the distance and speed of the target from the IF signal. The detection principle is shown in Figure 1, where the initial phase of the modulation signal is zero, and the echo signal does not consider the phase shift and amplitude attenuation for the time being, where f 0 is the initial frequency, c is speed of light, A 0 is amplitude of the transmitted signal, τ is time delay between the echo signal and the transmitted signal, and μ is Frequency modulation slope. The calculation formula is μ = 2 B / T . T is modulation period and B is modulation bandwidth.
For the convenience of research, this model only analyzes the up-sweep period. At this time, the up-sweep modulation period T m = 0.5 T and the IF frequency of the up-sweep period satisfy f I F = f I F , u p = f I F , d o w n . Then, there is:
R = f I F , u p + f I F , d o w n T c 8 B = f I F + f I F 2 T m c 8 B = T m c 2 B f I F
The target distance R need to be calculated by f I F .
Existing simulations of FMCW laser fuzes in smoke interference environment assume that both the fuze and target are in a stable smoke environment, assuming that the smoke scattering coefficient and extinction coefficient are fixed values. In fact, as fuze gradually approaches the target covered by smoke, the occurrence of smoke interference is a random dynamic process from scratch. In our earlier work, the aerosol extinction coefficient and scattering coefficient changed dynamically according to the different smoke density variation, aerosol particle distribution, and photon spatial positions. Based on [9], in this paper, we take account of the multiple scattering and multipath propagation of photons, and establish a multiple scattering model of the FMCW laser in the smoke environment. The scattering model and detection process are shown in Figure 2.
According to the transmission process of the laser in the smoke environment, three types of echo signal photons can be received by the fuze receiving system at a certain time, including: (1) photons that are directly transmitted to the target surface after being emitted, and are directly received by the fuse after being reflected by the target; (2) photons that do not reach the target surface scattered by the aerosol and return to the receiving syatem; (3) photons that reach the target surface scattered by the aerosol during transmission and eventually return to the receiving system. At this time, the power P E of the laser echo signal can be expressed as:
P E = P E 1 + P E 2 + P E 3
where P E 1 , P E 2 , and P E 3 correspond to the laser echo signal power generated by the above three different types of photons.
(1)
For the first case, the photon is directly reflected by the target and then received by the receiving optical system. During the multiple scattering process, the optical path length accumulated, and time delay occurred between the echo signal and the transmitted signal. During the approach of the projectile-target encounter, if the forward speed of the fuze is v, and the target is detected for consecutive k times, each time interval is T c , the delay time τ 1 generated by the photon transmitting between fuze and target can be calculated by:
τ 1 = 2 ( R k T c v ) c , k = 0 , 1 , 2
where R is the distance between fuze and target. If the transmission power is P T 1 , the power of the echo signal formed by the first type photons is:
P E 1 ( t ) = A 1 · c o s 2 π f 0 ( t τ 1 ) + π μ ( t τ 1 ) 2
A 1 = T S P T 1 ρ T e σ ( 2 ( R k T c v ) ) cos θ r 1 ( 2 ( R k T c v ) ) 2
T S = D 2 η T T T R
where σ is the extinction coefficient of the atmosphere to the laser, θ r 1 is the direction angle when the photon is received, T S is the coefficient related to the optical system characteristics of the fuze, D is the diameter of aperture of the receiving optical system, η is the emission system efficiency factor, T T is the laser emission solid angle, and T R is the transmission coefficient of the receiving optical system.
(2)
For the second case, received photons don’t reach the target surface and are all scattered by smoke. Assuming that the laser emission power is P T 2 , the photon reaches the surface of the smoke after moving R S , and then returns to the fuze after scattering once. According to the single scattering model [20], if the delay time is not considered, the echo signal power obtained by the scattering is:
P δ V = V 2 T S [ exp ( σ ( 2 R S ) ) c o s θ r ( 2 R S ) ] × γ s c a exp ( γ e x t · l 1 ) l 1 2 P T 2 d V
where γ s c a is the smoke aerosol scattering coefficient, γ e x t is the smoke aerosol extinction coefficient, l 1 is the optical path of photon scattered once, θ r is the direction angle when the photon is received by receiving optical system, and d V represents the unit space volume at which photons are scattered.
We construct multiple scattering models based on single scattering model [20,21], the power expression of a photon after i times scattering is:
P δ V i = P δ V i 1 v i γ s c a · P ( θ i ) exp ( γ e x t · l i ) l i 2 d V i
where θ i is the direction angle at the i-th scattering. According to the Mie scattering theory, the scattering phase function P ( θ ) in a smoke aerosol environment can be considered to be related only to the photon scattering direction angle, using the Henyey–Greenstein phase function [H-G] expressed [22]:
P ( θ ) = 1 g 2 ( 1 + g 2 2 g cos θ ) 3 2
where g is asymmetric factor and its value range is [−1,1].
The delay time τ 2 of the photon through N 2 times of scattering can be calculated by:
τ 2 = 2 R S 2 k T c v c + i = 1 N 2 l i c , k = 0 , 1 , 2
The echo signal power composed of the second type of photons after N 2 times of scattering is:
P E 2 N 2 t = A 2 ( i , N 2 ) · cos { 2 π [ f ( t τ 2 ) + 1 2 μ ( t τ 2 ) 2 ] }
A 2 ( i , N 2 ) = V 1 V N 2 T S P T 1 ρ s γ s c a N 2 cos θ r N 2 p θ 1 p θ 2 · · · p θ N 2 × exp [ σ ( 2 ( R S k T c v ) ) γ e x t i = 1 N 2 l i ] ( 2 R S 2 k T c v ) 2 l 1 2 l 2 2 · · · l N 2 2 d V 1 · · · d V N 2
where ρ s is the reflection coefficient of aerosol material.
In the study, it is assumed that there are M 2 photons scattered by the aerosol into the receiving system; they scattered N 2 times. The signal formed by the i-th scattered photon can be expressed as P E 2 ( i , N 2 ) , the scattering angle is θ i , N 2 , the scattering optical path is l i , N 2 , and the direction angle when photon is received is θ r ( i , N 2 ) . The power of the echo signal formed by the first photon can be expressed as:
P E 2 ( t ) = i = 1 M 2 P E 2 ( i , N 2 ) ( t )
(3)
For the third case, photons that reach the target surface are scattered by the aerosol and eventually return to the fuze receiving system. The scattering process of photons in smoke is consistent with the first case, assuming that the power of the laser-emitting signal is P T 3 . If the photon collides with the target at the i-th scattering, the target is the Lambert scatterer, and the reflection coefficient of the target is ρ T . The direction angle is θ T i , so the laser power reflected by the target is:
P δ V i + 1 = P δ V i × ρ T cos ( θ T i )
The multiple scattering process caused the optical path length to accumulate and generated time delay. The delay time τ 3 of the photon through N 3 times of scattering can be calculated by:
τ 3 = 2 R S 2 k T c v c + i = 1 N 3 l i c , k = 0 , 1 , 2
Therefore, the echo signal power composed of the third type of photons after the N 3 -th scattering can be expressed as:
P E 3 N 3 t = A 3 ( i , N 3 ) · cos { 2 π [ f ( t τ 3 ) + 1 2 μ ( t τ 3 ) 2 ] }
A 3 ( i , N 3 ) = V 1 V N 3 T S P T 3 ρ T cos θ T N 3 ρ s γ s c a N 3 cos θ r N 3 p θ 1 p θ 2 · · · p θ N 3 × exp [ σ ( 2 ( R S k T c v ) ) γ e x t i = 1 N 3 l i ] ( 2 R S 2 k T c v ) 2 l 1 2 l 2 2 · · · l N 3 2 d V 1 · · · d V N 3
If there are M 3 photons of the second type, the echo signal power composed of the second type of photons can be expressed as:
P E 3 ( t ) = i = 1 M 3 P E 3 ( i , N 3 ) ( t )
If the photoelectric conversion efficiency of the photodetector is η e , after the photoelectric converter is converted, the received signal can be expressed as:
s r ( t ) = η e P E 1 ( t ) + P E 2 ( t ) + P E 3 ( t )
In this study, the local oscillator signal expression is:
s L O ( t ) = A 0 cos 2 π ( f 0 t + 1 2 μ t 2 )
According to the principle of FMCW laser detection system, shown in Figure 1, the IF signal of the fuze is obtained by mixing the echo signal with the local oscillator signal. In addition, the random noise of the system is unavoidable. There are various sources of this random noise, including circuit thermal noise, dark current noise, photon shot noise, mixer frequency conversion loss, noise from AD sampling and low pass filtering, and so on. This random noise can be simulated with Gaussian white noise. Therefore, a Gaussian white noise component n ( t ) is added to the IF signal, and the IF signal expression is:
s I F ( t ) = s r ( t ) · s L O ( t ) + n ( t )
If the mixing process of the local oscillator signal and the echo signal does not consider the signal attenuation and filters out the higher harmonic term, the IF signal of the FMCW laser fuze in the approach of the projectile can be expressed as:
s I F ( t ) = η e A 0 2 { A 1 · cos 2 π ( f 0 + μ t ) 2 R 2 k T c v c + A 2 ( i , N 2 ) · cos 2 π ( f 0 + μ t ) 2 R S 2 k T c v c + j = 1 N 2 l i , j c + A 3 ( i , N 3 ) · cos 2 π ( f 0 + μ t ) 2 R S 2 k T c v c + j = 1 N 3 l i , j c } + n m ( t ) = s I F 1 ( t ) + s I F 2 ( t ) + s I F 3 ( t ) + n m ( t )
Under smoke interference conditions, the target echo signal s I F 1 ( t ) , the smoke-scattered target echo signal s I F 2 ( t ) , the IF signal is a mixture of the smoke backscattered echo signal s I F 3 ( t ) , and the random noise. s I F 1 ( t ) is the key signal required for fixed-distance. Ideally, s I F 1 ( t ) is a single frequency cosine signal whose frequency has a linear relationship with the distance between fuze and target.
Visibility is represented by smoke density during simulation. The relationship between smoke density and visibility is as follows [23]:
C = c 0 v x γ
where v x is the visibility and the unit is m; c 0 and γ are constants. Typical values for these constants are c 0 = 37.3 and γ = 1.07 [23,24]. Normally, when the visibility is 10 m, smoke density is 3.2 g/m 3 , and when the visibility is 50 m, the smoke density is 0.57 g/m 3 . In Section 4, we set the smoke density to 3.2 g/m 3 and aerosol partical diameter to 0.1–100 μ m. In Section 5, we used an atmospheric visibility sensor to calibrate the visibility of the experimental environment to 10 m.

3. Accurate Fixed-Distance of FMCW Laser Fuze Using a CSD Algorithm

3.1. Correntropy Spectral Density Theory

CSD can compare the similarity of two signals in the frequency domain. Therefore, when using CSD algorithm for fixed-distance, computers are used to calculate IF signal frequency under ideal conditions according to the expected distance, and get the signal waveform, that is, the waveform of the reference IF signal. Ideally, an IF signal is a single frequency cosine signal whose frequency has a linear relationship with the distance between fuze and target. The IF signal under ideal conditions is referred to as the reference IF signal. According to the principle of FMCW laser detection, the reference IF signal s I F e ( t ) is a cosine signal with a certain frequency f I F e . f I F e can be calculated by Equation (24):
f I F e = 2 B T m c R e
where R e is the expected distance.
This paper uses the CSD to compare the similarity in frequency domain between the reference IF signal and the measured IF signal as a criterion for fuze determining distance. If a spectral peak appears at f I F e in correntropy spectral distribution, it means that both signals include this frequency and actual distance from the target is close to the expected distance.
In this paper, the correntropy expression [13,25] of measured IF signal s I F ( t ) and reference IF signal s I F e ( t ) is:
V s I F ( t ) , s I F e ( t ) = E κ σ s I F ( t ) , s I F e ( t ) = E 1 2 π σ exp s I F ( t ) s I F e ( t ) 2 2 σ 2
where σ is the Gaussian kernel bandwidth. In this study, σ is determined by Silverman’s rule of density estimation [26]:
σ = 4 a 5 3 N s 1 5
where a is standard deviation of the sampled signal and N s is the number of samples.
The Mercer kernel function κ σ ( t ) used for correntropy is a normal distribution pattern [25]. When the two input signals s I F ( t ) and s I F e ( t ) are similar in the time domain, the value of s = s I F ( t ) s I F e ( t ) is around 0, and the amplitude information is preserved. If the difference between the two signals is large, the uncorrelated part of the two signals is rapidly attenuated to zero in the correntropy. The area of attenuation and the rate of attenuation can be controlled by the setting of σ . Therefore, only the similar components in the two signals are retained in the correntropy, and the non-correlated parts are discarded, which can greatly reduce the difficulty of information extraction.
In the implementation of laser fuze engineering, s I F ( t ) and s I F e ( t ) are discrete digital signals, which can be expressed as s I F ( n ) n = 1 N s and s I F e ( n ) n = 1 N s , so the correntropy between them can be estimated by the following formula [14]:
V ^ ( m ) = 1 N m + 1 n = m N κ σ s I F ( n ) , s I F e ( n m )
Fourier transform of V ^ ( m ) can obtain the correntropy spectral density of s I F ( n ) n = 1 N s and s I F e ( n ) n = 1 N s to compare the two signals in the frequency domain similarity. The expression is [14]:
C ( f ) = m = + V ^ ( m ) e j 2 π f m 0
Whether the two signals are similar can be discriminated by the presence or absence of a peak at a reference frequency in the correntropy spectral density distribution. When identifying the frequency peaks in the CSD distribution, a rectangular window function such as Equation (29) can be used to intercept data in a specified bandwidth range to determine whether a peak exists:
W ( f ) = 1 , f I F e 1 2 Δ f I F e f f I F e + 1 2 Δ f I F e 0 , other
where Δ f I F e is the effective bandwidth near the peak frequency f I F e . It should be set according to the fixed-distance accuracy requirement of the FMCW laser fuze. For example, when the fixed-distance error is ± Δ R e , Δ f I F e can be calculated by:
Δ f I F e = 2 × 2 B T m c Δ R e
If the effective bandwidth Δ f I F e is too wide, the fixed-distance accuracy will be set at a low level, resulting in the position of laser fuze detonation is not ideal. Due to the high speed of the FMCW laser fuze, the signal processing time is very limited, and a narrow effective bandwidth Δ f I F e may cause the recognition failure of the target signal spectrum peaks.
In the absence of smoke interference, there is only target signal s I F 1 ( n ) and random noise n i ( n ) in s I F ( n ) . At this time, if the frequency of reference IF signal is consistent with the target signal frequency, then s I F 1 ( n ) and s I F e ( t ) have only the difference in amplitude and phase. If the delay between s I F 1 ( n ) and s I F e ( t ) is D n , then there is s I F 1 ( n D n ) = s I F e ( n ) . In the expression of V ^ ( m ) , there is s I F ( n D n ) = s I F e ( n m ) . When m = D n , correntropy has a maximum value. Correspondingly, a distinct peak appears at the frequency f I F e in the correntropy spectral density distribution. At this time, the current distance from target can be determined to satisfy the fixed-distance condition. Under the condition of smoke interference, the target information is hard to extract. In the correntropy, only the similar components in the two signals are retained, and the non-correlated parts are discarded, which can greatly reduce the difficulty of information extraction.

3.2. FMCW Laser Fuze Fixed-Distance Processing Using a CSD Algorithm

The basic steps of the CSD-based fixed-distance algorithm are as follows:
(1)
Obtain the currently measured IF signal.
(2)
Calculate the reference frequency f I F e according to the distance R e , and generate a reference IF signal with a frequency f I F e .
(3)
Calculating V ^ ( m ) of the currently measured intermediate frequency signal and the reference intermediate frequency signal.
(4)
Do Fourier transform for the V ^ ( m ) to get C f .
(5)
Determine whether the peak appears in the CSD processing result. If the peak appears at the reference intermediate frequency f I F e , step 6 ) is entered; if there is no peak, the first step is returned to start the detection again.
(6)
Recalculate the distance according to the motion speed of the fuze and the detection time interval, and begin the second detection at step 1 ) until the peak of the continuous measurement is stable, which proves that the target is detected.
Figure 2 shows the process of continuous detection. For the reference IF signal, if the initial value of R e is R 1 , the complete expression of f I F e is:
f I F e = 2 B T m c R e = 2 B T m c R 1 k 1 T c v k = 1 , 2
where k is the number of detections, T c is the detection time of the fuze, and v is the motion speed of the fuze.

4. Simulation

This section studies the CSD algorithm to identify target signals in the case of target without smoke interference, smoke interference but no target, and the target is covered by smoke. The CSD algorithm processing results are compared with the FFT algorithm processing results to further verify the advantages of the CSD algorithm in terms of under smoke interference conditions. The specific parameters used in simulation are shown in Table 1. According to the requirements, the fixed-distance accuracy of the fuze is ±30 cm, so Δ f I F e is set to 4 KHz calculated by Equation (30). Other parameters setting is consistent with [9]. Table 2 shows theoretical frequency corresponding to 3 m, 5 m, 9 m, and 12 m calculated according to Equation (24).
In this section, when expected distance is 3 m, the simulation results are compared when the distance is 3 m and 5 m between fuze and target. When expected distance is 9 m, the simulation results are compared when the distance is 9 m and 12 m between fuze and target. Each simulation result includes two consecutive detetion, the red curve represents the first detection, and the black curve represents the second detection.

4.1. Simulation without Smoke Interference

Figure 3 and Figure 4 show the simulation results under the condition of no smoke interference.
It can be seen from the simulation results that, when there is no smoke interference, if the distance between the fuze and the target is equal to the expected distance, both the CSD algorithm and the FFT algorithm processing results have obvious peaks at the reference intermediate frequency, both satisfying the fixed-distance condition. However, the amplitude of the spectral peak of the CDS algorithm is much larger than the peak of the FFT algorithm. When the distance between the fuze and the target is not equal to the expected distance, the processing results of the two methods have no peak near the reference intermediate frequency, indicating that the fixed-distance condition is not satisfied. Since the correntropy is always greater than zero, when the correntropy is Fourier transformed to obtain the correntropy spectral density, zero-frequency interference occurs in the result. Zero-frequency interference can be eliminated by the window function without affecting the identification of the peak of the target signal.

4.2. Simulation with Smoke Interference but No Target

Figure 5 and Figure 6 show the simulation results under the condition of smoke interference but no target.
When there is no target, it can be found that there is no obvious peak in the result of the CSD algorithm from the results of simulations. However, affected by smoke interference, the result of FFT algorithm has multiple peaks, which can’t achieve fixed-distance and may cause false alarms.

4.3. Simulation of Target Being Covered by Smoke

Figure 7 and Figure 8 show the simulation results under the condition of target being covered by smoke.
From simulation results Figure 7 and Figure 8, it can be found that the CSD algorithm results in Figure 7b and Figure 8b showing obvious peaks at the reference intermediate frequency, satisfying fixed-distance condition. Figure 7e and Figure 8e show smaller or no peaks. It can be considered to be caused by smoke interference and does not satisfy the fixed-distance condition. The FFT algorithm is greatly affected by smoke interference, and multiple peaks appear. It is impossible to judge whether the distance determination condition satisfies.
In summary, the CSD algorithm can effectively suppress smoke interference. When performing correntropy calculations, the Gaussian kernel function rapidly attenuates the components whose frequency is not equal to f I F e and the target signal with the frequency f I F e is retained, thereby suppressing the smoke interference. The FFT algorithm cannot effectively suppress the interference of smoke.

5. Experiment

This section builds the FMCW laser fuze test platform, and conducts experiments in different environments. A typical semiconductor laser ( λ = 940 nm, 50 mW) was used in the FMCW laser fuze. After being collimated by the optical system, the laser spot size is about 2.5 mm × 1.0 mm, and beam-divergence angle is less than 3 degrees. An avalanche photo-diode detector (APD, AD500-10) was used as the laser detector. In addition, the target in the experiment is a white aluminum plate (1.5 m × 1 m), and its reflectivity is about 80 % . The collected IF signal data are transmitted to the computer for processing to verify the anti-interference effect and feasibility of the CSD algorithm. The prototype of the FMCW laser fuze is shown in Figure 9a. Figure 9b show the data processing platform.
In experiments, the fuze is facing the target, and the target surface is perpendicular to the emitted laser beam. The fuze adopts a triangular wave frequency modulation system, and the DDS generates a chirp signal with a bandwidth of 50 MHz, and the modulation period is 100 μ s, wherein the up-and-down frequency sweep period is 50 μ s that is, T m = 50 μ s. For the convenience of research, the AD sampling frequency is set to 20.48 MHz, the single sampling period is 400 μ s, the sampling point is 8196 points, and the cutoff frequency of the IF signal spectrum is 150 kHz.

5.1. Performance Test of Fuze Prototype under Conditions of No Smoke Interference

Figure 10 and Figure 11 show the experiment results under the condition of no smoke interference, when fuze is 3 m and 9 m away from the target. The time domain waveforms, frequency domain waveforms of IF signal and the processing results using the CSD algorithm are included in the figures.
The spectrum of IF signal is obtained by fast Fourier transform. It can be found that, in the correntropy spectral density distribution, as shown in Figure 10c and Figure 11c, the spectral peak amplitudes at frequencies of 20.48 kHz and 60.15 kHz are significantly higher than those of other frequency peaks, indicating that the correntropy of the measured IF signal and the reference IF signal reaches maximum value. The similarity is the largest. Table 2 shows that the theoretical frequencies corresponding to 3 m and 9 m are 20.01 kHz and 60.04 kHz which are used as the conventional true values. Based on the proportional relationship between the difference value and the conventional true value, the relative error of the frequency does not exceed 5 % . The absolute error of distance can also be calculated by Equation (24). The absolute error does not exceed 0.15 m.

5.2. Performance Test of Fuze Prototype under Smoke Interference

This section verifies the anti-interference ability of FMCW laser fuze in the smoke environment. The test scenario schematic is shown in Figure 12.
The experiment room size is 15 m × 4 m wide, test window size is 0.3 m × 0.3 m and confined diffusion space which is 2 m × 4 m wide. Gradually release smoke from smoke source in confined space to reduce visibility. We used an atmospheric visibility sensor to calibrate the visibility of the experimental environment to 10 m. The laser beam emitted by the prototype enters the smoke diffusion space through the test window. Because high speed movement of fuze cannot be simulated indoors, this paper lists the results of two detections of the target under the same conditions. Figure 13 and Figure 14 show the experiment results under the condition of smoke interference, when fuze is 3 m and 9 m away from the target.
From the experimental results, it can be found that the smoke interference increases the noise in the IF signal, and the time domain waveform is obviously distorted. The noise in the spectrum increases, and a high-amplitude interference peak appears on the left of the target signal peak. The smoke interference signal may be used as the target signal to cause false alarm. Under smoke interference conditions, significant peaks can still appear at the reference frequency in the CSD algorithm processing results. The peak frequencies at 3 m and 9 m from the target are 19.6 kHz and 59.5 kHz, respectively. Table 2 shows that the theoretical frequencies corresponding to 3 m and 9 m are 20.01 kHz and 60.04 kHz. Based on the same error analysis in Section 5.1, the relative error between measured peak frequency and theoretical frequency does not exceed 5%, and the absolute error of distance calculated by Equation (24) does not exceed 0.15 m.

6. Conclusions

In this paper, a multiple scattering model of photons in smoke environment is established. This model takes account of the multiple scattering of photons and multi-path propagation. The CSD algorithm is used to process the data, and an accurate fixed-distance algorithm of FMCW laser fuze is proposed. The experimental results show that the CSD algorithm has good anti-interference ability in smoke environments, and our experimental results are consistent with the simulation results. The relative error between measured peak frequency and theoretical value does not exceed 5%, and the absolute error of distance does not exceed 0.15 m. Our detection system uses single-quadrant and single-detector achieving high accuracy, so our system is simpler and easier to implement in engineering.

Author Contributions

C.S. managed the design and development of the model and algorithm, and contributed to the creation of the main ideas in the paper. Y.C. handled the writing of the manuscript. B.L. focused on the experimental task (i.e., setting up the system, as well as the design and realization of the experiments). All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grant No. 61973038).

Acknowledgments

Yabo Duan and Zhihua Pang conducted the data analysis and interpretation tasks together.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Arora, V. Proximity Fuzes: Theory and Techniques; Defence Research and Development Organisation, Ministry of Defence: New Delhi, India, 2010.
  2. Pfennigbauer, M.; Wolf, C.; Weinkopf, J.; Ullrich, A. Online Waveform Processing for Demanding Target Situations; International Society for Optics and Photonics, RIEGL Laser Measurement Systems GmbH: Horn, Austria, 2014; Volume 9080, p. 90800J. [Google Scholar]
  3. Von der Lippe, C.M.; Liu, J.J. Advanced optical fuzing technology. In Optical Technologies for Arming, Safing, Fuzing, and Firing; Optics and Photonics: San Diego, CA, USA, 2005; Volume 5871, pp. 123–130. [Google Scholar]
  4. Liu, J.J.; Christian, M. Advanced high-bandwidth optical fuzing technology. In Nanophotonics for Communication: Materials and Devices II; International Society for Optics and Photonics: Boston, MA, USA, 2005; Volume 6017, p. 60170F. [Google Scholar]
  5. Grabner, M.; Kvicera, V. The wavelength dependent model of extinction in fog and haze for free space optical communication. Opt. Express 2011, 19, 3379–3386. [Google Scholar] [CrossRef] [PubMed]
  6. Mei, L.; Lundin, P.; Andersson-Engels, S.; Svanberg, S.; Somesfalean, G. Characterization and validation of the frequency-modulated continuous-wave technique for assessment of photon migration in solid scattering media. Appl. Phys. B 2012, 109, 467–475. [Google Scholar] [CrossRef]
  7. Zhang, W.; Li, Y.; Huang, Z. Research on the characteristics of fog backscattering signals for frequency modulated continuous wave laser fuze. Optik 2016, 127, 9046–9055. [Google Scholar] [CrossRef]
  8. Zhang, W.; Li, Y.; Huang, Z.; Ma, C. Fog backscattering interference suppression algorithm for FMCW laser fuze based on normalized frequency spectrum threshold. Optik 2017, 131, 188–193. [Google Scholar] [CrossRef]
  9. Liu, B.; Song, C.; Duan, Y. The characteristics simulation of FMCW laser backscattering signals. Opt. Rev. 2018, 25, 197–204. [Google Scholar] [CrossRef]
  10. Liu, B.; Song, C.; Li, Y.; Duan, Y. A frequency-modulated-continuous-wave laser detection system based on the four-quadrant photodetector. Appl. Phys. B 2018, 124, 186. [Google Scholar] [CrossRef]
  11. Santamaría, I.; Pokharel, P.P.; Principe, J.C. Generalized correlation function: Definition, properties, and application to blind equalization. IEEE Trans. Signal Process. 2006, 54, 2187–2197. [Google Scholar] [CrossRef] [Green Version]
  12. Principe, J.C. Information Theoretic Learning: Renyi’s Entropy and Kernel Perspectives; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
  13. Liu, W.; Pokharel, P.P.; Príncipe, J.C. Correntropy: Properties and Applications in Non-Gaussian Signal Processing. IEEE Trans. Signal Process. 2007, 55, 5286–5298. [Google Scholar] [CrossRef]
  14. Garde, A.; Karlen, W.; Ansermino, J.M.; Dumont, G.A. Estimating respiratory and heart rates from the correntropy spectral density of the photoplethysmogram. PLoS ONE 2014, 9, e86427. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Özbek, M.E.; Savacı, F.A. Correntropy function for fundamental frequency determination of musical instrument samples. Expert Syst. Appl. 2011, 38, 10025–10030. [Google Scholar] [CrossRef]
  16. Hassan, M.; Terrien, J.; Marque, C.; Karlsson, B. Comparison between approximate entropy, correntropy and time reversibility: Application to uterine electromyogram signals. Med. Eng. Phys. 2011, 33, 980–986. [Google Scholar] [CrossRef] [PubMed]
  17. Garde, A.; Sörnmo, L.; Jané, R.; Giraldo, B.F. Correntropy-based spectral characterization of respiratory patterns in patients with chronic heart failure. IEEE Trans. Biomed. Eng. 2010, 57, 1964–1972. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Chen, B.; Xing, L.; Zhao, H.; Zheng, N.; Prı, J.C. Generalized correntropy for robust adaptive filtering. IEEE Trans. Signal Process. 2016, 64, 3376–3387. [Google Scholar] [CrossRef] [Green Version]
  19. Yu, L.; Qiu, T.S.; Song, A.M. A Time delay estimation algorithm based on the weighted correntropy spectral density. Circuits Syst. Signal Process. 2017, 36, 1115–1128. [Google Scholar] [CrossRef]
  20. Reilly, D.M.; Warde, C. Temporal characteristics of single-scatter radiation. JOSA 1979, 69, 464–470. [Google Scholar] [CrossRef]
  21. Elshimy, M.A.; Hranilovic, S. Nonline-of-sight single-scatter propagation model for noncoplanar geometries. JOSA A 2011, 28, 420–428. [Google Scholar] [CrossRef]
  22. Toublanc, D. Henyey–Greenstein and Mie phase functions in Monte Carlo radiative transfer computations. Appl. Opt. 1996, 35, 3270–3274. [Google Scholar] [CrossRef]
  23. Ryde, J.; Hillier, N. Performance of laser and radar ranging devices in adverse environmental conditions. J. Field Robot. 2009, 26, 712–727. [Google Scholar] [CrossRef]
  24. Widzyk-Capehart, E. Seeing through Dust and Water Vapor: Millimeter Wave Radar Sensors for Mining Applications. J. Field Robot. 2007, 24, 527–557. [Google Scholar]
  25. Xu, J. Nonlinear Signal Processing Based on Reproducing Kernel Hilbert Space. Ph.D. Thesis, University of Florida, Gainesville, FL, USA, 2007. [Google Scholar]
  26. Dehnad, K. Density Estimation for Statistics and Data Analysis; CRC Press: Boca Raton, FL, USA, 1987. [Google Scholar]
Figure 1. Frequency-modulated-continuous-wave (FMCW) laser detection principle. IF signal of FMCW laser fuze is obtained by mixing the echo signal with the local oscillator signal.
Figure 1. Frequency-modulated-continuous-wave (FMCW) laser detection principle. IF signal of FMCW laser fuze is obtained by mixing the echo signal with the local oscillator signal.
Sensors 20 02604 g001
Figure 2. Multiple scattering model and detection process, simulating the movement of fuze with speed of v and a time interval of T c between two detections.
Figure 2. Multiple scattering model and detection process, simulating the movement of fuze with speed of v and a time interval of T c between two detections.
Sensors 20 02604 g002
Figure 3. Simulation results when expected distance is 3 m without smoke interference. When fuze is 3 m away from the target, (a) is time domain waveform of the IF signal; the processing results of CSD algorithm and FFT algorithm are (b,c), respectively. When fuze is 5 m away from the target, (d) is time domain waveform of IF signal; the processing results of CSD algorithm and FFT algorithm are (e,f), respectively.
Figure 3. Simulation results when expected distance is 3 m without smoke interference. When fuze is 3 m away from the target, (a) is time domain waveform of the IF signal; the processing results of CSD algorithm and FFT algorithm are (b,c), respectively. When fuze is 5 m away from the target, (d) is time domain waveform of IF signal; the processing results of CSD algorithm and FFT algorithm are (e,f), respectively.
Sensors 20 02604 g003
Figure 4. Simulation results when expected distance is 9 m without smoke interference. When fuze is 9 m away from the target, (a) is time domain waveform of IF signal; the processing result of CSD algorithm and FFT algorithm are (b,c), respectively. When fuze is 12 m away from the target, (d) is time domain waveform of IF signal; the processing result of CSD algorithm and FFT algorithm are (e,f), respectively.
Figure 4. Simulation results when expected distance is 9 m without smoke interference. When fuze is 9 m away from the target, (a) is time domain waveform of IF signal; the processing result of CSD algorithm and FFT algorithm are (b,c), respectively. When fuze is 12 m away from the target, (d) is time domain waveform of IF signal; the processing result of CSD algorithm and FFT algorithm are (e,f), respectively.
Sensors 20 02604 g004
Figure 5. Simulation results when expected distance is 3 m with smoke interference but no target. When fuze is 3 m away from smoke, (a) is time domain waveform of the IF signal; (b) is the processing result of the CSD algorithm; (c) is the processing result of the FFT algorithm.
Figure 5. Simulation results when expected distance is 3 m with smoke interference but no target. When fuze is 3 m away from smoke, (a) is time domain waveform of the IF signal; (b) is the processing result of the CSD algorithm; (c) is the processing result of the FFT algorithm.
Sensors 20 02604 g005
Figure 6. Simulation results when expected distance is 9 m with Smoke interference but no target. When fuze is 9 m away from smoke, (a) is time domain waveform of the IF signal; (b) is the processing result of CSD algorithm; (c) is the processing result of the FFT algorithm.
Figure 6. Simulation results when expected distance is 9 m with Smoke interference but no target. When fuze is 9 m away from smoke, (a) is time domain waveform of the IF signal; (b) is the processing result of CSD algorithm; (c) is the processing result of the FFT algorithm.
Sensors 20 02604 g006
Figure 7. Simulation results when the target is covered by smoke with expected distance 3 m. When fuze is 3 m away from the target, (a) is time domain waveform of IF signal; the processing result of CSD algorithm and FFT algorithm are (b,c), respectively. When fuze is 5 m away from the target, (d) is time domain waveform of IF signal; the processing results of CSD algorithm and FFT algorithm are (e,f), respectively.
Figure 7. Simulation results when the target is covered by smoke with expected distance 3 m. When fuze is 3 m away from the target, (a) is time domain waveform of IF signal; the processing result of CSD algorithm and FFT algorithm are (b,c), respectively. When fuze is 5 m away from the target, (d) is time domain waveform of IF signal; the processing results of CSD algorithm and FFT algorithm are (e,f), respectively.
Sensors 20 02604 g007
Figure 8. Simulation results when the target is covered by smoke with expected distance 9 m. When fuze is 9 m away from the target, (a) is time domain waveform of IF signal; the processing result of CSD algorithm and FFT algorithm are (b,c), respectively. When fuze is 12 m away from the target, (d) is time domain waveform of IF signal; the processing result of CSD algorithm and FFT algorithm are (e,f), respectively.
Figure 8. Simulation results when the target is covered by smoke with expected distance 9 m. When fuze is 9 m away from the target, (a) is time domain waveform of IF signal; the processing result of CSD algorithm and FFT algorithm are (b,c), respectively. When fuze is 12 m away from the target, (d) is time domain waveform of IF signal; the processing result of CSD algorithm and FFT algorithm are (e,f), respectively.
Sensors 20 02604 g008
Figure 9. (a) shows the composition of FMCW laser fuze prototype; (b) the collected IF signal data are transmitted to the computer for processing.
Figure 9. (a) shows the composition of FMCW laser fuze prototype; (b) the collected IF signal data are transmitted to the computer for processing.
Sensors 20 02604 g009
Figure 10. Experiment result when expected distance is 3 m with no smoke interference. (a) is time domain waveform of IF signal; (b) is spectrum of IF signal; (c) is processing result of the CSD algorithm.
Figure 10. Experiment result when expected distance is 3 m with no smoke interference. (a) is time domain waveform of IF signal; (b) is spectrum of IF signal; (c) is processing result of the CSD algorithm.
Sensors 20 02604 g010
Figure 11. Experiment result when expected distance is 9 m with no smoke interference; (a) is time domain waveform of IF signal; (b) is spectrum of IF signal; (c) is processing results of the CSD algorithm.
Figure 11. Experiment result when expected distance is 9 m with no smoke interference; (a) is time domain waveform of IF signal; (b) is spectrum of IF signal; (c) is processing results of the CSD algorithm.
Sensors 20 02604 g011
Figure 12. Schematic test scenario of fuze prototype under smoke interference. The smoke interference environment test was conducted in a closed room, and confined diffusion space was isolated by partitions.
Figure 12. Schematic test scenario of fuze prototype under smoke interference. The smoke interference environment test was conducted in a closed room, and confined diffusion space was isolated by partitions.
Sensors 20 02604 g012
Figure 13. Experiment result when expected distance is 3 m under smoke interference. (a) is time domain waveform of IF signal; (b) is a spectrum of the IF signal; (c) is the processing result of the CSD algorithm. Each simulation result includes two detetions; the red curve represents the first detection, and the black curve represents the second detection.
Figure 13. Experiment result when expected distance is 3 m under smoke interference. (a) is time domain waveform of IF signal; (b) is a spectrum of the IF signal; (c) is the processing result of the CSD algorithm. Each simulation result includes two detetions; the red curve represents the first detection, and the black curve represents the second detection.
Sensors 20 02604 g013
Figure 14. Experiment results when expected distance is 9 m under smoke interference. Figure layout is the same as Figure 13. Each simulation result includes two detections, the red curve represents the first detection, and the black curve represents the second detection.
Figure 14. Experiment results when expected distance is 9 m under smoke interference. Figure layout is the same as Figure 13. Each simulation result includes two detections, the red curve represents the first detection, and the black curve represents the second detection.
Sensors 20 02604 g014
Table 1. Simulation parameters.
Table 1. Simulation parameters.
ParametersValue
modulation period T m 50 μ s
modulation bandwidth B50 MHz
detection time of fuze T c 1 ms
the speed of light c299,792,458 m/s
speed of fuze v200 m/s
Table 2. The distance and frequency.
Table 2. The distance and frequency.
DistanceFrequency
3 m20.01 kHz
5 m33.36 kHz
9 m60.04 kHz
12 m80.06 kHz

Share and Cite

MDPI and ACS Style

Song, C.; Cui, Y.; Liu, B. FMCW Laser Fuze Multiple Scattering Model and Accurate Fixed-Distance Algorithm in a Smoke Environment. Sensors 2020, 20, 2604. https://doi.org/10.3390/s20092604

AMA Style

Song C, Cui Y, Liu B. FMCW Laser Fuze Multiple Scattering Model and Accurate Fixed-Distance Algorithm in a Smoke Environment. Sensors. 2020; 20(9):2604. https://doi.org/10.3390/s20092604

Chicago/Turabian Style

Song, Chengtian, Ying Cui, and Bohu Liu. 2020. "FMCW Laser Fuze Multiple Scattering Model and Accurate Fixed-Distance Algorithm in a Smoke Environment" Sensors 20, no. 9: 2604. https://doi.org/10.3390/s20092604

APA Style

Song, C., Cui, Y., & Liu, B. (2020). FMCW Laser Fuze Multiple Scattering Model and Accurate Fixed-Distance Algorithm in a Smoke Environment. Sensors, 20(9), 2604. https://doi.org/10.3390/s20092604

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