Abstract

Overlapping signals separation is a difficult problem, where time windowing is unable to separate signals overlapping in time and frequency domain filtering is unable to separate signals with overlapping spectra. In this work, a simulation under MATLAB is implemented to illustrate the concept of overlapping signals. We propose an approach for resolving overlapping signals based on Fourier transform and inverse Fourier transform. The proposed approach is tested under MATLAB, and the simulation results validate the effectiveness and the accuracy of the proposed approach. The approach is developed using Gerchberg superresolution technique to cope with signals with low signal-to-noise ratio. For practical work, an echo shape determination is required to apply the proposed technique. The experimental results show accurate localization of multiple targets.

1. Introduction

In measurement systems based on pulse echo techniques, some problems may occur when received signal contains many overlapping reflected echoes due to many factors: multiple targets, structure of the propagation media that may be consisting of several layers; and the complex physical properties of the propagation path. Therefore, it is important to find an approach for overlapping echoes separation encountered in a wide range of applications such as military and biomedical applications, [15].

Overlapping problem occurs in both time and frequency domains. Time overlapping signals are encountered when time delay between echoes is shorter than the emitted signal’s duration. Overlapping is due to a combination of layer thickness or distance between reflectors [1, 5].

A difficulty arises in separating the composite signal, when two or more of the individual component signals have spectra that are the same or similar. This is also a well-known difficulty in musical sound separation when the harmonics of two or more pitched instruments have similar frequencies [6].

The frequency-phase relationship of the individual signals is incoherent causing constructive and destructive interference, resulting in irregular characteristics of the spectrum of the composite signal rather than the individual signals [1]. So the relative phase of overlapping signals is found to influence the precision of frequency determination and plays a critical role in the observed amplitude of the combined signal; thus the amplitude and the phase of individual signals become unobservable [6, 7].

So it is a challenging problem in signal processing to obtain the spectrum and the time localization of the individual component signals, as they overlap in both time and frequency domains, because it is difficult for traditional methods like spectrum method and time-frequency method to separate overlapping signals. Time domain windowing is unable to separate signals overlapping in time. Frequency domain filtering is unable to separate spectrally overlapped signals due to the similarity between the spectra of the signals.

Previous studies related to signal separation have emerged in recent years. Several different approaches such as short-time Fourier transform, Wigner-Ville distribution, discrete wavelet transform, discrete cosine transform, and chirplet transform, and Fractional Fourier transform have been proposed to cope with the overlapping problem [1, 8].

In this paper, we propose a method for resolving overlapping signals based on Fourier transform and inverse Fourier transform. Proposed method performance for noisy signals is improved using Gerchberg superresolution technique.

To apply the proposed method in practical application (directed acoustic transmission-reception system in our case), it is necessary to obtain the echo shape derived from experiments, so Hilbert transform is used to extract the received signal envelop, and the extracted envelop is used to modulate the transmitted signal.

This paper is organized as follows: Section 2 introduces the concept of overlapping signals and illustrates the problem encountered by a simulation under MATLAB. In Section 3, the proposed technique is introduced and illustrated by simulation results. Section 4 presents experimental echo shape determination by extracting the received signal envelop using Hilbert transform, where the transmitted signal is modulated by the extracted envelop. Section 5 shows experimental examples results of applying the proposed technique to received signal, using acoustic transmission-reception system, and multiple targets. Finally, conclusion is drawn in Section 6.

2. Concept of Overlapping Signals

The concept of overlapping signals is illustrated by a simulation implemented in MATLAB. The composite signal is a summation of a number of individual component signals (three overlapping echoes); the individual component signals have the same frequency (2 kHz), different magnitudes, and different time delays as shown in Figure 1. The summation of overlapping signals is shown in Figure 2. Figure 3 shows the spectra of the three individual component signals after Fourier transformation, and Figure 4 shows the spectrum of the composite signal after Fourier transformation.

In all simulations, the following parameters will be used: the sampling frequency  kHz and the acquisition time  sec, so the number of samples , discrete Fourier transform (DFT) computed with a 100000-point FFT, and the frequency increment  Hz.

As illustrated in Figure 4, frequency domain analysis of time overlapping signals using Fourier transformation does not produce the spectrum of the individual component signals. The individual component signals have regular spectra and main lobe centred at (2 kHz), but the spectrum of the composite signal has distortions and irregular characteristics.

3. Proposed Method

3.1. Theoretical Background

Assuming that the transmitted signal is ,where is the function shape of the transmitted pulse and is the transmission frequency. The received signal could, for instance, be several () scaled and time delayed copies of a known signal :where is the amplitude of the shape of the th echo pulse and is its time delay. By taking the Fourier Transform of (2), we getwhere is the Fourier Transform of (1). We define the transform function in the frequency domain by

By taking the Inverse Fourier Transform of (4), we get

Equation (5) shows that the resulting transform function in the time domain is several () pulses with amplitude and time delay .

By applying this method on signal illustrated in Figure 2, where is the composite signal, we get transform functions in frequency domain and time domain as shown in Figures 5 and 6.

Note that time delays and amplitudes of echoes, obtained in Figure 6, are in agreement with those given in Figure 1.

3.2. Effect of Signal-to-Noise Ratio (SNR)

In practical application, signals could not be as the simulated transmitted and received signals due to the noise, frequency shift, and system response.

Unfortunately, the proposed method is very sensitive to noise, which may cause a catastrophic effect on results; to illustrate this point, a white noise with signal-to-noise ratio () will be added to the composite signal shown in Figure 2, so the noisy composite signal becomes as shown in Figure 7. Transform functions and in frequency and time domain are shown in Figures 8 and 9, respectively.

Narrow pulses with high amplitude shown in Figure 8 are caused by dividing by zeros (of ), resulting in periodic sinusoidal signals in all the time domain, as shown in Figure 9; these signals cover pulses at time delays of our interest.

To cope with problem, (4) could be developed as the following form:where is a small real value related to SNR, representing the ratio of the power spectral density of the noise and that of the signal [9], is the complex conjugate of , , and are Fourier transforms of cross-correlation and autocorrelation functions, respectively [10]. This technique used in (6) is similar to the Wiener Filter used in image processing [9], and it is powerful in measurement conditions with low ratio, where the autocorrelation function contains the same information about the signal, and, by using the cross-correlation function, it is possible to obtain an increase in the ratio (the Wiener-Khintchine theorem) [11, 12].

Figure 10 shows the transform function of previous example according to (6) with .

As shown in Figure 10, the technique in used in (6) has reduced the noise effect (the effect of dividing by zero), but this is not enough to recognize pulses of interest in transform function , so an improvement based on super resolution technique is required.

3.3. Superresolution Technique

The proposed method performance is improved using Gerchberg superresolution technique [13], as illustrated in the following steps:(1)The received signal is windowed by a time window , with zero padding outside the window, where and are limits of time window of interest; this is done by detection algorithm or by dividing time domain of the received signal into predefined time-length frames.(2)Taking the Fourier transform of the signal obtained in step 1, (6) is used with a suitable predefined value for ; then we get .(3)Transform function obtained in step 2 is filtered by time band filter , where and are limits of the same time window in step 1, is the filtered transform function of , and is defined as the following:where and are frequency band limits of interest (around the transmission frequency ).(4)Taking the inverse Fourier transform of , obtained transform function is windowed by the same time window , but, to compensate the filter delay used in step 3, time window is shifted by the same filter delay and then is obtained.(5)Taking Fourier transform of to get and , (7) is used to get a new .(6) is filtered by time band filter , and is the inverse Fourier transform of .(7)A predefined number of iterations are used to get the best resolution for by successive steps: taking the Fourier transform of the windowed with filter delay compensation to get and , (7) is used to get a new , filtering by time band filter to get new , taking inverse Fourier transform of to get new .

Modified transform functions , (red curve), and its absolute value (blue curve) of the previous example, in frequency and time domain, of the third-order iteration are shown in Figures 11 and 12, respectively.

Time delays obtained in Figure 12 are in agreement with those obtained in Figure 6 and given in Figure 1.

In practical application, the problem encountered is to get the reference transmitted signal that really propagates in the medium, where the received signal is considered to be several scaled and time delayed copies of it. So we have to construct heuristically the reference transmitted signal that has the same transmission frequency but modulated by shape function (pulse envelope); it is supposed that echoes have the same shape but with scaled amplitude.

In the following paragraph, we propose an approach to extracting the pulse envelope from the received signal using Hilbert transform.

4. Experimental Envelope Detection

The Hilbert transform facilitates the formation of the analytic signal . The analytic signal is useful in the area of communications, particularly in pass band signal processing.

Figure 13 shows the schematic representation of experimental setup used for envelope detection.

This experimental setup is used in a closed and empty Basketball Hall that is considered a nearly anechoic environment, where microphone and loudspeaker are located about 2 m above the ground, and the distance from the microphone or the loudspeaker to the closest object is greater than 5 m.

First, we generate 10 periods of mathematical sine wave at ( Hz) frequency with a rectangular envelope as a transmission pulse (pulse width = 5 msec). Using a laptop, this pulse is transmitted through a sound card, an acoustic power amplifier, and a loudspeaker. Using a microphone, the pulse is received directly and acquired through the sound card; the received signal is processed to obtain its analytic signal , that is,where is the instantaneous amplitude and is the Hilbert Transform of , which is given by [14, 15]

Figure 14 shows the mathematical generated signal. Figure 15 shows the experimental reference signal and its extracted envelope using Hilbert transform, implemented under MATLAB. This envelope derived from experiments (received signal shifted to time 0), supposed to be similar to echo shape, is used to modulate mathematical sine wave at frequency to be the reference transmitted signal. This model of transmitted signal and reflected echoes takes into account the response of the transmitter and receiver sensors, the number of cycles of the excitation pulse, and the propagation medium characteristics.

5. Experimental Examples

The proposed method is experimentally tested using single, double, and triples targets and a directed acoustic transmission-reception system with a parabolic dish with a focal length of 35 cm. The parabolic dish has an elliptical aperture with a major axis of 65 cm and a minor axis of 60 cm. Figure 16 shows a schematic representation of experimental setup used for proposed method verification.

In all experiments, signals are sampled at the sampling frequency  kHz over an acquisition time  sec, so the number of samples acquired is samples, discrete Fourier transform (DFT) is computed with a 100000-point FFT, and the frequency increment is  Hz.

5.1. Single Target

At first, a single board of wood with dimensions 90 × 60 cm is used as a target, which is located at about 330 cm from the directed acoustic system, as shown in Figure 16. The modified transform function (red curve) and its absolute value (blue curve) are shown in Figure 17.

From Figure 17, target distance is = ((19.27 msec × 340 m))/2 = 327.6 cm; the result is satisfactorily close to the real value.

5.2. Double Target

Two boards of wood with dimensions 60 × 40 cm and 90 × 60 cm are used as a target (the smaller is in front, and the distance between them is about 30 cm); the target is located at about 8 m from the directed acoustic system, as shown in Figure 16. Figures 18, 19, and 20 show the received signal, the transform function with , and the modified transform function of third-order iteration, respectively.

From Figure 20, target distance is = ((46.69 msec × 340 m))/2 = 7.94 m, and distance between two targets is = (((48.5 − 46.69) m·sec × 340 m))/2 = 30.77 cm; results are satisfactorily close to real values.

5.3. Triple Target

Three boards of wood with dimensions 60 × 40 cm, 90 × 60 cm and 160 × 140 cm are used as a target (the smallest board is the closest to the system and the biggest one is the farthest); the target is located at about 9.25 m from the directed acoustic system, as shown in Figure 16, the distance between the first board and the second one is about 45 cm, and the distance between the second board and the third one is about 80 cm; Figure 21 shows the received signal, Figure 22 shows the transform with , and the modified transform function of third-order iteration is shown in Figure 23.

From Figure 23, target distance is = ((54.44 m·sec × 340 m))/2 = 9.255 m, the distance between first and second target is = (((57.19 − 54.44) m·sec × 340 m))/2 = 46.75 cm, and the distance between second and third target is = (((62.13 − 57.19) m·sec × 340 m))/2 = 83.98 cm; results are satisfactorily close to real values.

6. Conclusion

This paper presents a method to separate overlapping echoes using a superresolution algorithm based on experimental echo shape modelling.

The simulation results implemented under MATLAB show a good performance of the method at low SNR.

The method is experimentally tested using acoustic transmission-reception system and multiple targets; the method provides accurate localization of multiple targets due to a good accuracy in determining the time delays of individual echoes reflected from each target. Although amplitudes of individual echoes, obtained by this method, are different from the real amplitudes and smaller than them, each echo’s amplitude relative to others is respected.

This study adds a new possibility of overlapping echoes separation and can be widely used in a range of applications such as military and biomedical applications.

Developments can be added to improve SNR and determine the amplitude of individual echoes that allow reconstructing the composite signal.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.