Noise reduction of continuous wave radar and pulse radar using matched filter and wavelets

This paper analyzes noise reduction using matched filter and wavelet transform in the signals of continuous wave radar and pulse radar. The denoising application of wavelets has been used in spectrum cleaning of atmospheric radar signals. Matched filter has a strong anti-noise ability; it can also achieve accurate pulse compression in a very noisy environment. This paper analyzes the algorithms of matched filter and wavelets that are used in radar signal processing to reduce the noise. The simulation results indicate that matched filter has a strong anti-noise ability for pulse radar and wavelet for continuous wave radar.


Introduction
Chaos is the very complicated behavior of a low-order dynamical system, because it is both nonlinear and deterministic [1,2]. It demonstrates a strong notion, permitting the use of a simple deterministic system to illustrate highly irregular fluctuations exhibited by physical phenomena encountered in nature. Recently, some engineering applications of chaos have been reported in literature [3][4][5][6]. These can be grouped under two broadly defined categories [3,4]. One group is synthesis of chaotic signals, which includes signal masking and spread-spectrum communications. Another is analysis of chaotic signals. It exploits the fact that some physical phenomena allow the use of a chaotic model.
For more than 20 years, many methods for smoothing denoising filters have been proposed. Various standard spatial domain filters have been proposed to reduce noise [7]. The 2-D Gaussian filter [8] and the median filter [9] have also been widely applied. The most commonly used is the Lee filter [10], but it requires knowledge of interferometric coherence. The adaptive contoured-window filter method [11] and the two-stage filtering method [12] can also reduce noise effectively.
Wavelet analysis [5] and matched filter [6] are two of the most significant tools in the field of signal processing in the last few decades. We analyze both techniques and show how such an individual can improve the quality of radar-received signals in a noisy environment for both types of radar, i.e., continuous wave radar (radars continuously transmit a high-frequency signal and the reflected energy is also received and processed continuously) and pulse radar (transmits high power, high-frequency pulses toward the target, and it waits for the echo of the transmitted signal before it transmits a new pulse). The problem addressed here concerns the denoising of the radar-received signal immersed in noise. Several simulations have been performed to verify the algorithm for both types of radar. All of the simulations give the same time delay, and even the received signal is attenuated more than 90%. The simulation results indicate that the algorithm is effective and robust even when the receiver receives a very weak signal.
The remainder of the paper is organized as follows: In Section 2, we briefly describe the literature of matched filter and the characteristics of matched filter related to radar noise reduction. Section 3 presents wavelets denoising technique, while the simulations and data analysis are described in Section 4. Finally, we make some conclusions about our comparison related to noise reduction.

Theory of matched filter
Matched filter is not a specific type of filter, but a theoretical frame work. It is an ideal filter that processes a received signal to minimize the effect of noise. Therefore, it optimizes the signal-to-noise ratio (SNR) of the filtered signal [13].
The original matched filter was proposed by D. O. North [14]. Later, many authors [15][16][17] tried to improve the performance of matched filter. However, the conventional matched-filter detector has some drawbacks [18]. To reduce the drawbacks, the exhaustivesearch matched filter (EMF) detector and the near-optimal MF (NMF) detector are proposed in [19]. Recently, [20] proposed a developed NMF (D-NMF) to improve the performance of NMF.
The radar-received signal r(t) contains two components s i (t) and n i (t) that represent the certain signal (e.g., targets) and noise, respectively, i.e., r(t) = s i (t) + n i (t). The matched filter h(t) yielding the output y(t) = s o (t) + n o (t) is to generate the maximal ratio of s o (T) and n o (T) in the sampling values at time T. Where s o (t) and n o (t) are the outputs of s i (t) and n i (t), respectively, after the matched filter shown in Figure 1.
In Figure 1, s i (t) represents the target signal we are attempting to detect and n i (t) represents the additive white Gaussian noise in the system. s o (t) and n o (t) represent the outputs of the matched filter by s i (t) and n i (t). Matched filter maximizes the SNR, the ratio of the power of s o (t), and the power of n o (t) according to the Schwarz inequality.
Here, we suppose that the noise n i (t) additive white Gaussian noise, whose power spectrum is N/2, and the spectrum of the target's signal s i (t) is [13] The above equation is the Fourier transform of s i (t). The output of the matched filter y(t) also contains two components representing the target's signal and noise, The average power of the noise equals the value of autocorrelation function, which is Now, according to the definition of SNR, time T is Using the Schwarz inequality, we get Hence, The numerator of the above equation denotes power of the signal according to Parseval's theorem. From Equation 5, the matched filter maximizes the SNR of the filtered signal and has an impulse response that is a reverse time-shifted version of the input signal. So, to obtain the maximum SNR, we need the time delay, D.
With the use of this time delay, D, we obtain the output of the cross-correlation between transmitted signals and received signals.

Cross-correlation to find time delay D
Correlation is the process to determine degree of 'fit' between two waveforms and to determine the time at which the maximum correlation coefficient or 'best fit' occurs [21][22][23]. For the radar system, if we correlate between the transmitted signal and the received signal, then we get the time difference between the transmitted and received signals. We consider the transmitted signal to be x(n), and then the returned signal r(n) may be modeled as: where w(n) is assumed to be the additive noise during the transmission, α is the attenuation factor, and D is the delay which is the time taken for the signal to travel from the transmitter to the target and back to the receiver.
The cross-correlation between the transmitted signal, x(n), and the received signal, r(n), is [21] From Equation 7, the maximum value of the crosscorrelation will occur at l = D, which is our interest in cross-correlation from which we can get the time delay, D. For the multiple targets, we get the multiple number of D from Equation 7. For example, if there are n targets then we can get n number of delays such as D 1 , D 2 , D 3 ,…, D n .

Matched filter for radar
For pulse radars, consider pulse width τ ρ and τ κ is the time that a target is illuminated by the radars. Thus, we can write r(t) as From Equation 5, and for Equation 8, we can write the pulse radar as [24] where P t is the peak transmitted power of radar, G is the antenna gain, σ is the radar cross section (RCS), R is the range which electromagnetic wave transmits, λ is the wavelength, K = 1.38 × 10 -23 J/K is Boltzmann's constant, B is bandwidth, T 0 = 290 K is the operating temperature of antenna, F is the noise figure of receiver, and L denotes as radar losses.
For continuous wave radars, radar equation can be written as [24] where P CW is the continuous wave average transmitted power and T DWELL is the dwell interval.

Wavelet denoising
For removing noise and extracting signal from any data, wavelet analysis is one of the most important methods. The wavelet denoising application has been used in spectrum cleaning of the atmospheric signals. There are different types of wavelets available like Morlet, Coiflet, Mexican hat, Symlet, Biorthogonal, and Haar, which have their own specifications such as filter coefficients and reconstruction filter coefficients. In this paper, to eliminate noise embedded in the radar signal 'sym8, ' wavelets have been used. The goal of this study is to denoise the radar signal. One often encounters the term 'denoising' in recent wavelet literature, described in an informal way with various schemes that attempt to reject noise by damping or thresholding in the wavelet domain [25,26]. The threshold of wavelet coefficient has near optimal noise reduction for different kinds of signals. Wavelets have many advantages over fast Fourier transform. Fourier analysis has a major drawback, which is that time information is lost, when transforming to the frequency domain. Thus, it is impossible to tell when a particular event took place under Fourier analysis. Wavelet analysis is capable of revealing aspects of data that other signal analysis techniques, aspects such trends, breakdown points, discontinuities in higher derivatives, and self-similarity, are unable to reveal. Wavelet analysis can often denoise a signal without appreciable degradation. Wavelet transform performs a correlation analysis. Therefore, the output is expected to be maximal when the input signal most resembles the mother wavelet.

Wavelet transform
According to the definition of wavelet transform [27], for function f(t), wavelet transform coefficient W f (a,τ) Here, f(t), ψ (a,τ) (t) is the wavelet basis function, ψ Ã t−τ a À Á is a conjugate of wavelet basis function, τ is the amount of shift, and α is scale.

Wavelet denoising
The wavelet denoising procedure proceeds in three steps: Step 1 Signal decomposing: Choose the wavelet basis function, and to determine the decomposition level N, get the coarse and detail coefficients by DWT.  Step 2 Threshold detail coefficients: For each level from 1 to N, compare the detail coefficients to threshold values.
Step 3 Reconstructing the signal: Reconstruct the denoised signal based on the original approximation coefficients of level N and the modified detail coefficients of levels from 1 to N.

Traditional threshold function
The major signal information mainly concentrates in the low frequency sub-band of wavelet transform domain.  Noise equally distributes in all wavelet coefficients, so the wavelet transform factor should be larger than the wavelet transform factor of the noises after wavelet decomposition. Therefore, the selection of wavelet threshold is an important step, which also directly impacts on the effect of noise reduction. Different methods have been proposed to choose the threshold. The frequently used thresholding of wavelet coefficients is governed mainly by either soft or hard thresholding function, proposed by Donoho [28]. The soft thresholding is generally referred to as wavelet shrinkage, since it 'shrinks' the coefficients with high amplitude toward zero, whereas the hard thresholding is commonly referred to simply as wavelet thresholding. Given that d jk indicates the value of wavelet coefficient,d jk implies the value of d jk after thresholding function, and T is the threshold value.
The soft thresholding function is defined as T is the threshold and generally can be a function of J and K. The hard thresholding function is defined as Soft thresholding provides smoother results in comparison with the hard thresholding whereas thresholding technique provides better edge preservation in comparison with the soft thresholding technique.
Soft thresholding and hard thresholding have some limitations in denoising of signal [29,30]. The Equation 12 indicates that the reconstructed signal faces oscillation, since the estimated wavelet coefficients are not continuous at position ± T [31]. Although the estimated wavelet coefficients of Equation 13 have good continuity, these coefficients include constant errors [31], which directly influence the accuracy of the reconstructed signal.

An improved threshold function
To overcome the limitations of hard thresholding and soft thresholding denoising methods, an improved thresholding is proposed as follows: Equation 14 improves the reconstruction precision, since it reduces the constant errors. Hence, it enhances the denoising effect. This thresholding function also assures the continuity of estimated wavelet coefficients.

Simulation and performance assessment
The denoising of the received radar signal is simulated in the presence of white Gaussian noise. The effect of  signal parameter changes on the algorithm has been investigated. These parameters include the SNR of the signal. The SNR is defined as the ratio of the signal power to the noise power in the entire period. The following parameters are assumed: sampling frequency = 10 GHz, pulse duration = 8 ns, pulse repetition frequency = 0.24 GHz for the pulse radar ( Figure 2) and transmitted frequency = 10 GHz for the continuous wave radar ( Figure 3) The receiver receives the return from the targets in the present of AWGN.
We recover our transmitted signal from the received signal using the matched filter for both pulse radar and continuous wave radar. For the pulse radar, the recovery signal is almost similar to transmitted signal shown in Figure 2c. But in the case of continuous wave radar, it is almost impossible to understand the shape of the transmitted signal from the recovery signal ( Figure 3). Figures 2 and 3, respectively, except the received signals are attenuated more than 90%. Figure 4 indicates that the matched filter has a strong advantage in anti-noise ability for pulse radar. Unfortunately, the matched filter for continuous wave radar is not superior to recover the transmitted signal (Figure 5c). The reason for this is that the transmitted average power of the continuous wave radar is lower, and the thermal noise relevant to the operating bandwidth and the radar's structure design, so the signal detection threshold of continuous wave radar is lower, and this demands higher receiver sensitivity (Table 1). This is shown in Figure 6 using Equations 9 and 10. Also, applying the matched filter for continuous wave radar limits the performance in three major areas [32]: the effect of nonlinearities in the frequency transmission and the matched filter, variations in the transmitted frequency, and relative velocity of target.

Figures 4 and 5 are similar to
We also apply the wavelet denoising technique to remove the noise from received signal. This technique overcomes the problem of continuous wave radar. Figure 7c shows the noise reduction by traditional wavelet for continuous wave radar. By comparing this with Figure 3c, we can say that more noise is reduced in Figure 7c. Figure 8c shows the improvement in noise reduction by our proposed method. The reason for this is that the SNR is increased (Table 1) when we use the wavelet denoising technique for continuous wave radar. Our proposed wavelet threshold function even further reduced more noise for continuous wave radar. This is shown in Figure 9.

Conclusions
This paper presented a comprehensive comparison between matched filter and wavelet in terms of noise reduction abilities of pulse radar and continuous wave radar. The simulation results show that a significant reduction in noise is achieved for pulse radar by matched filter, but employing matched filter for continuous wave radar increases the difficulty of detection. The use of wavelet denoising technique instead of matched filter for continuous wave radar reduced more noise. Our proposed wavelet threshold function even further reduced more noise for continuous wave radar.  Figure 9 Plots of SNR versus detection range RCS = 1,000 m 2 for continuous wave radar.