Abstract

An investigation on the improvement of the resolution of a radar target identification system is presented in this paper. Degradation of resolution is mainly due to influence factors associated with antennas, including the strong coupling between transmitting and receiving antennas and the variation in the antenna response. A filtering technique was therefore introduced to mitigate the underlying problem. In the technique, the antenna effects were filtered out of the total response backscattered from the objects in the radar target identification system. The short-time matrix pencil method (STMPM) was then employed to extract the poles from the backscattered response in order to identify the object. Simulation and experimentation examples are illustrated to confirm the improvement of the resolution by filtering the antenna effects. The simulation and experimentation were divided into several categories, that is, different antennas and differently shaped objects, in order to validate the advantage of filtering the antenna effects. They were setup in order to demonstrate that the poles obtained from performing the STMPM without the filtering technique were mainly because of the antenna rather than the object’s characteristic. The results showed that the resolution of the identification was significantly increased when performing pole extraction and filtering the antenna effects.

1. Introduction

Since over half a century, radar target identification systems have been intensely investigated in research areas of electromagnetic (EM) wave propagations and signal processing. There are many potential applications of the system, such as air defense, ground penetrating radar (GPR), and chipless RFID identification systems [15]. The challenging problem in the system is identifying unknown objects with different shapes and constituents, buried or not underground. The singularity expansion method (SEM) is one of the most employed methods to characterize the electromagnetic response of a radar target. The SEM was first introduced by observing the time-domain signal backscattered from the target [6]. By illuminating the target with the wideband EM pulse, the scattered response in the time domain is composed of two successive parts, namely, early-time and late-time responses. The first portion comes from the direct reflection of the incident EM pulse on the target. The second is due to the resonance phenomena of the target. A key to the SEM is to model the late-time portion as a sum of the complex exponentials, called the poles. Therefore, there is a need for a scheme to be employed in order to extract the poles from the late-time response.

The most popular method that has been widely employed for pole extraction is the matrix pencil method (MPM) because of its performance in dealing with noise [7]. However, the main drawback of applying MPM for identification is that it requires correct isolation between the early-time and late-time responses. In 2008, the time frequency domain was used to observe the evolution of partial resonances in the early-time responses and the commencement of the late-time response [8]. Thus, this approach can overcome the weakness of requiring the correct isolation between the early-time and late-time responses. Recently, the short-time matrix pencil method (STMPM) has been also proposed to overcome this weakness by sliding a proper time window of the conventional MPM’s (5). It has been applied for detection in chipless RFID systems. An alternative method proposed in [9], based on a sliding-window analysis, applied the Kronecker’s theorem to determine the characteristic changes in the Hankel matrix, including the matrix rank and the distribution of eigenvalues. The change is an intrinsic indicator of the stage of the target response that is currently windowed. This novel method can determine the commencement of the late-time response without a priori knowledge of the target geometry or orientation. More recently, we proposed a simple technique to estimate the commencement of the late-time response by using a cross correlation between the backscattering and transmitted waves [10]. Using the STMPM along with the proposed late-time estimation, the extracted poles will converge to a stable constant value at the position of the superimposed pattern between the transmitted waves and the late-time response. Thus, the correct extracted poles can be achieved without human skill or experience. This proposed technique was then applied to the GPR system in order to achieve the automatic detection and classification [11].

Besides to the precise pole extraction, suitable RF components are required for the radar target identification system. A popular solution is to use ultra-wideband components, including antennas and RF chain [12]. Unfortunately, the use of ultra-wideband components is a high-cost solution and designing microwave components that possess the characteristics of no variation in the ultra-wideband frequency range becomes a more serious problem. An alternative candidate has been proposed to mitigate the antenna effects by using the deconvolution technique [13, 14]. In addition, antenna calibration was also employed to resolve the problems of antenna dispersion and multiple reflections in radar systems [15, 16]. However, the antennas used in these literatures were designed for a particular purpose, especially radar systems and their bandwidths were ultra-wideband. For example, the horn antenna can operate in the frequency range of 800 MHz to 10 GHz [15]. Moreover, the effect of the use of different antenna types has not been evaluated, and nobody has seriously taken the pole resolution into account. A solution to resolve this problem due to the antenna effects is introduced here along with the investigation conducted for the pole resolution improvement.

In this paper, we present the investigation of the resolution improvement in terms of poles for a radar target identification system. A filtering technique will be introduced to mitigate the ground bounce and antenna effects due to coupling between transmitting and receiving antennas and the variation in the antenna response. The STMPM was employed to extract poles from the scattering response. A few of the antennas under consideration was the off-the-shelf products which were not designed for the radar applications. The performances of the radar identification system using off-the-shelf and not off-the-shelf antennas were compared. The simulation and experimentation were setup to demonstrate if the resolution degradation was mainly because of the antenna effects rather than objects. The results show that the resolution of the identification is then significantly increased when performing pole extraction along with filtering the antenna effects. Following this, filtering the antenna effects for a bistatic radar system is then introduced along with derivation of its equations in Section 2, in which a target identification technique based on the STMPM is discussed. Simulations of the target identification with and without filtering the antenna effects are presented in Section 3. The experimentation results are also shown in order to confirm the superiority of the proposed approach in Section 4. Finally, conclusions are drawn in Section 5.

2. System Model

2.1. Filtering the Antenna Effects for the Bistatic Radar

Antennas are essential components in a radar target identification system. Typically, there is a variation in the antenna response throughout the range of operating frequency. The influence factor associated with antennas results in not only waveform distortion but also resolution degradation in target identification. In this section, a mathematical description of the proposed technique being employed to resolve the underlying problem is given. Figure 1 shows a basic model of a bistatic radar system. The receiving (Rx) antenna is placed with the transmitting (Tx) antenna as a bistatic radar configuration. The electromagnetic (EM) pulse s(t) is fed to the transmitting antenna, whose impulse response is . When the output signal of the transmitting antenna strikes an object which is modeled by the impulse response of hs(t), the signal scattered from the object will be received by the receiving antenna whose impulse response is . The scattering signal r(t) received by the receiving antenna contains the desired information of the object characteristic and the undesired information, such as the free space and antenna effect, as expressed as where pt(t) and pr(t) are the transfer functions of free space for the incident and reflected signals, respectively. The hc is the loss between the transmitting and receiving antennas, and “” denotes and the convolution. Reconsidering Figure 1, the generated signal s(t) was transmitted by using the transmitting antenna and then propagated through the free space channel pt(t) to strike an object hs(t). The signal backscattered from the object propagates through the free space channel pr(t) and was then received by the receiving antenna . This behavior corresponds to the first term in (1). On the other hand, the transmitted signal corresponding to s(t) directly propagates to the receiving antenna. This direct path from the transmitting to the receiving antenna corresponds to the second term in (1).

Based on the SEM principle, the extracted complex frequency is a tool for aspect-independent target identification. Thus, the free space loss does not affect the extracted complex frequency. Their transfer functions can be omitted from the equation, which is rewritten as

Here, the term includes the mutual coupling between the transmitting and receiving antennas and can be substituted with which is the antenna effects. Practically, can be obtained by measuring an empty room without objects. The received signal is therefore given as

Another choice to filter the antenna effects in the stepped frequency continuous wave (SFCW) radar system is to perform the signal processing technique in the frequency domain. By taking a Fourier transform of (3), the frequency response of the target is derived according to the following equation: where Hs(ω), R(ω), S(ω), and are the Fourier transforms of hs(t), r(t), s(t), and , respectively. So far, one can find the object characteristic by filtering the antenna effects given by (3) or (4).

2.2. Target Identifications

The SEM was introduced by observing the transient response of an object [6]. When an object was illuminated by an ultra-wideband EM wave, the scattering transient response given in (3) can be divided into two successive portions as expressed by where rET(t) and rLT(t) denote the early-time and late-time responses of the received signal, respectively. η(t) denotes the noise in the system, uncorrelated with either responses. The first response comes from the direct reflection of the incident wave on the object surface. The second one is due to the object’s resonance phenomena which can be separated into internal and external modes. The internal resonances are caused by the internal waves experiencing multiple internal reflections while the surface creeping waves cause the external resonance. The signal model of the late-time response from an object can be formulated as a sum of damped exponentials of the following form: where is an ith pole including the damping factor σi and natural frequency ωi. Ri and M are the residues (complex amplitudes) and the number of poles, respectively. The pole is a tool for aspect-independent object identification, depending upon the geometry and physical properties of the object. The commencement of the late-time period coincides with the end of the forced response, due to scattering from discrete object surfaces, and is therefore dependent on the object’s geometry and its orientation with reference to the excitation wavefront. On the other hand, the residues are aspect dependent since they represent the coupling between incident electric wave and induced currents on the object. Following the pole extraction using the STMPM in [5], the windowed late-time response can be written as in which T is the shifting time and

Taking the logarithmic scale, (8) can be expressed as

According to (9), the logarithmic magnitude of the residues linearly decreases versus T with slope αi. To obtain the damping factor, one can find the slope of the logarithm of the residues in the time frequency scale as well.

3. Performance Validation via Electromagnetic Simulations

Simulation examples are illustrated to confirm that the proposed technique, which filters the antenna effects out of the total received backscattering signal, can improve the pole resolution for the radar target identification. The simulations were conducted by using the Higher Order Basic Integral Equation Solver (HOBBIES) EM software. This was done in the frequency domain with the step frequency (Δf) of 1 MHz. The frequency was swept from 1 MHz to 1200 MHz, and the number of the sample frequency (Nf) was 1201. In simulations, two same antennas were placed together as transmitting and receiving antennas with the spacing of 100 cm. The object was located in the free space at 200 cm away from the center of the transmitting and receiving antennas. Simulations were divided into four categories regarding to the antenna and object’s shape. The antennas which were bow-tie and cone antennas were chosen from the example list of the EM simulation software. The return loss of the bow-tie and cone antennas was below −10 dB from 416 MHz to 487 MHz and from 120 MHz to 160 MHz, respectively. The bandwidths of these antennas were somewhat narrow, compared with those of the typical radar antenna. In the simulations, the object was constructed with two different shapes, that is, PEC box and plate. The dimensions of the PEC box and plate in cm were 20 × 20 × 30 and 50 × 60, respectively. In this section, simulation examples are distinguished according to the target identification with and without filtering the antenna effects as follows.

3.1. Target Identification without Filtering the Antenna Effects

Figure 2 shows the frequency response obtained from simulations of the target identification without filtering the antenna effects. It is clearly seen that the responses obtained from the use of the bow-tie and cone antennas were significantly different although the same object was under our consideration. On the other hand, the responses of the PEC box and plate were similar when the same antenna type was employed in the simulation. These results were not as expected for target identification. It implies that the antenna was a more important influence than the object’s shape.

To confirm the target identification, poles extracted from the scattering responses should be determined. The inverse fast Fourier transform (iFFT) was used to transform the frequency responses to time-domain responses. Poles were extracted from these time-domain responses by using the STMPM. Figures 3(a) and 3(b) show the imaginary (frequency) and real (damping factor) parts of the extracted poles, respectively. In the figure, the frequency shows better stability in the late-time than the damping factor. The damping factor was therefore achieved by finding the slope of the normal logarithm of the residues rather than directly finding the damping factor [5]. Figure 3(c) shows the magnitude of the residues in the logarithm scale. In general, the poles of objects with different shapes and/or constituents should be different. Nevertheless, the figures illustrate that the extracted poles from the PEC box and plate were almost identical for both cases of the bow-tie and cone antennas. The poles of the same object’s shape in the simulation with different antennas must be the same; however, the figures showed that the simulation with different antennas achieved different poles. Thus, these poles cannot be directly employed for target identification. This confirms the need to filter the antenna effects leading to the resolution degradation of target identification.

3.2. Target Identification with Filtering of the Antenna Effects

In this section, the simulation results are shown to confirm the need to filter the antenna effects. Equation (4) was employed to filter the antenna effects out of the total received response. Poles were extracted again from the time-domain responses by using the STMPM and shown in Figure 4. In the figure, the natural frequency of the PEC box converged to an average of the stable values of 264 MHz and 261 MHz at late-time tLT = 22 ns when employing the bow-tie and cone antennas, respectively. The averaged values of the extracted frequency of the PEC plate after the late-time response were 173 MHz and 169 MHz, respectively. The damping factor seen in Figure 4(b) begins to stick together when tLT = 22 ns with the averaged values of −6.3 × 108 and −6.1 × 108 Hz for the PEC box and −3.3 × 108 and −3.5 × 108 Hz for the PEC plate when employing the bow-tie and cone antennas, respectively. Note that the extracted poles including frequencies and damping factors were different when determining the differently shaped objects. Figure 4(c) shows the magnitude of the residues. In the figure, the solid line denotes the range of the constant slope of the line fitted to the natural logarithm of the magnitudes of residues by the linear regression. The slope is constant after tLT = 22 ns corresponding to the position of the stability for the extracted natural frequency and damping factor. This is used to confirm the commencement of the late-time response. Reconsidering Figure 4, the shading area denotes the region of the poles extracted from the late-time portions. Poles of the differently shaped objects in different regions showed a significant improvement in the pole resolution.

Table 1 summarizes the averaged pole values obtained from the STMPM with and without filtering the antenna effects. The averaged pole values were calculated from the average of the eighteen samples of the successive poles which were extracted in the late time. In the case of the STMPM without filtering the antenna effects, the averaged poles obtained from the use of the different antennas were significantly different although the same object was under consideration. The averaged poles of different objects were similar when the same antenna type was employed. After using filtering the antenna effects, the averaged poles of the same objects were similar and vice versa although either the same antennas or different antennas were under consideration. Moreover, these poles could be compared with the poles of the objects under plane wave illumination instead of the use of antennas and can be referred to as the reference poles. The table shows that poles obtained from performing the STMPM along with filtering the antenna effects were almost identical to the reference ones. These guarantee that the resolution of the identification was significantly increased when performing the pole extraction along with filtering the antenna effects.

At this end, the resolution of the target identification system must be determined through the Euclidean distance between the poles of the PEC box and plate in the simulations with and without filtering antenna effects. The eighteen samples of the successive poles obtained from the late-time responses of each object were selected for target classification, as shown in Figure 5. The Euclidean distance was calculated from the distance between each averaged pole. The shading symbols denote the averaged value of the poles of each simulation case. The circle line stands for the region of target classification, which was calculated by a half of the Euclidean distance. The Euclidean distances of the pole extraction with and without filtering the antenna were 290 MHz and 51 MHz when employing the bow-tie antennas. The Euclidean distances of 290 MHz and 10 MHz were achieved when employing the cone antennas in the simulations with and without filtering the antenna effects, respectively. The results clearly show that the Euclidean distance for target classification with filtering the antenna effects is longer than that of the standard method. This indicates that the resolution of target classification was improved by filtering the antenna effects.

4. Experimental Validation

The experimental setup with a bistatic configuration was conducted in order to validate the pole resolution improvement of the radar target identification, as depicted in Figure 6. This was done in an anechoic chamber in order to avoid reflection and other EM waves. For the radar target identification, the results obtained from the simulations and experimentation are not directly compared. Antennas used in the simulations demonstrate the proposed method conceptually. However, in our experiment, log-periodic, bow-tie, and Vivaldi antennas were selected based on the availability at the time of research and their dimensions in cm were 28.5 × 23, 17 × 27 × 13.5, and 55 × 70, respectively. The return loss of the log-periodic [17], bow-tie [18], and Vivaldi [19] antennas was below −10 dB in the range of 400–1000, 600–4200, and 240–17,000 MHz, respectively. Two of the same antennas were used as transmitting and receiving antennas. The distance between the transmitting and receiving antennas was 60 cm. The object was placed in such a direction that the electric field was propagated at one meter away from the center of the transmitting and receiving antennas. Here, two different shapes of objects, that is, brass plate and box, were determined in order to illustrate the identification. The dimensions of the brass plate and box were 50 × 60 cm2 and 20 × 20 × 30 cm3, respectively. A Rohde and Schwarz ZVB20 vector network analyzer was used to perform the experiment. The scattering coefficient S21 was measured during the frequency sweep in the range of 200 MHz to 1200 MHz for a total of 101 samples. Each measurement was repeated ten times to achieve accurate results.

4.1. Experimental Results without Filtering the Antenna Effects

Figure 7 shows the measured frequency responses without filtering the antenna effects. The figure clearly reveals that the responses obtained from the use of the same antennas were almost identical although the differently shaped objects were determined. On the other hand, the responses of the same objects were different when employing different antennas. The poles extracted by using the STMPM are shown in Figure 8. As seen in the figure, the dispersion in the poles, including the damping factor and frequency of the brass plate and box, was because of the antenna effects and nonidealities in the target identification. The extracted poles should converge to a stable constant value at the commencement of the late-time response [5]. However, in these cases, we could not find the commencement of the late-time responses. Thus, there were no correct poles being employed for the target identification. This guarantees that the problem of the degradation of the identification resolution due to the antenna effects must be resolved.

4.2. Experimental Results with Filtering the Antenna Effects

To improve the pole resolution, the preprocessing technique described by (4) was employed to filter the antenna effects out of the measured frequency response as shown in Figure 7. The resulting frequency response was then transformed to the time-domain response being used in the process of pole extraction. Figure 9 shows the poles extracted from the experimental results by using the STMPM along with filtering the antenna effects. Compared to the STMPM without filtering the antenna effects, Figure 9 clearly shows that the natural frequencies show a more stable behavior along the time axis. In order to classify the objects, ten samples of the successive poles were selected after the late-time of tLT = 20 ns. The selected poles were calculated for their averaged values. Table 2 summarizes the averaged pole values. In fact, we could not calculate the averaged pole values in the case of the STMPM without filtering the antenna effects because the certain commencement of the late-time response was nonexistent. In the table, these poles were not applicable (N/A). With filtering the antenna effects, the averaged poles can be achieved, as shown in the table. Note that the poles shown in Tables 1 and 2 were different. The poles obtained from the simulation and experimentation cannot be directly compared. This is because of the identification sensitivity due to the incomplete consistency between the simulation and experimentation conditions such as the antenna type and object material.

Figure 10 shows the ten samples of the successive poles extracted after tLT = 20 ns by using the STMPM with and without filtering the antenna effects. The shading symbol denotes the averaged pole value of each experimentation case. The Euclidean distances of the averaged poles attained from the STMPM without filtering the antenna effects were 98 MHz, 53 MHz, and 43 MHz when employing the log-periodic, bow-tie, and Vivaldi antennas, respectively. With filtering the antenna effects, the Euclidean distances of 148 MHz, 163 MHz, and 176 MHz were achieved when employing the log-periodic, bow-tie, and Vivaldi antennas, respectively. The significant difference in the extracted poles of different objects indicates that the STMPM with filtering the antenna effects successfully improved the resolution of the target identification. Furthermore, note that the extracted poles of the use of the commercial off-the-shelf log-periodic antenna whose bandwidth was not ultra-wideband were almost identical to those of the bow-tie and Vivaldi antennas which were designed for the particular purpose when performing the STMPM with filtering the antenna effects.

5. Conclusion

In this paper, the investigation of the resolution improvement of radar target identification systems has been presented. The improvement was achieved by filtering the antenna effects due to coupling between the transmitting and receiving antennas and the variation in the antenna response, out of the total backscattering received response. The STMPM was then employed to extract poles from the filtered response. The simulations with different antennas and differently shaped objects were exampled to validate the proposed approach. Furthermore, the experimentations were conducted to confirm the superiority of the proposed approach again. Without filtering the antenna effects, their results showed that one cannot classify objects because extracted poles of differently shaped objects were almost identical. The Euclidean distance of poles of differently shaped objects is larger when performing pole extraction along with filtering the antenna effects. The agreement between simulation and experimental results revealed that filtering the antenna effects can significantly improve the pole resolution.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work has been supported by the Thailand Research Fund through the TRF Senior Research Scholar Program (Grant no. RTA6080008).