Abstract

Unlike phased array bringing angle-dependent and range-independent beamforming, a new kind of radar system called frequency diverse array (FDA) has been widely used for its ability to generate angle-range-dependent beampattern, which provides a broad prospect of radar application. Nevertheless, FDA beampattern is coupling in the angle and range responses, which results in the difficulty for FDA to differentiate the target in two domains simultaneously. In this paper, a stepped frequency pulse FDA (SFP-FDA) radar is proposed to wipe off the angle-range coupling in FDA and improve the target’s angle and range localization accuracy. Compared to the double-pulse FDA, the proposed system can achieve a better performance in range resolution and multitarget identification. The Cramér-Rao lower bound (CRLB) is derived and used to analyze the performance of the proposed system and double-pulse FDA radar. Numerical results are implemented to verify the validity of the proposed scheme.

1. Introduction

Target localization is to find a way to calculate the position parameters of the target, and the angle and range parameters do draw much attention in localization [1]. The phased array radar [2, 3] has been widely employed for its capability to orientate the direction of the target [4, 5]. Due to its range ambiguity [6], phased array radar cannot distinguish targets from different ranges in a particular angle. To address the current needs for controlling range-dependent transmit energy distribution, a more agile array radar called frequency diverse array (FDA) radar is put forward by Antonik et al. [79]. By using a tiny frequency offset among antenna elements, FDA can achieve an angle-range-dependent beamforming [10]. The degree of freedom is increased, and this bidimensional property can be widely utilized, such as interference rejection [11], moving target indication [12], range ambiguity resolution, and target localization [13, 14].

Because of its wide application prospects, FDA has attracted extensive attention in recent years [1518]. A technology for FDA antenna systems that could generate angle-range-dependent beamforming is put forward in [19]. The multipath features of FDA on the ground propagation are studied in [20]. The Cramér-Rao lower bounds (CRLBs) of FDA for estimating range, angle, and velocity are derived in [21]. Some researchers are trying to find a transmit beampattern, which can last for a period of time in focusing signal power on the expected location [22]. A FDA structure based on periodic triangular frequency-modulated continuous waveform is put forward in [23]. To realize low probability of intercept, the cognitive FDA MIMO radar is presented in [24]. Global optimization algorithm like genetic algorithm is discussed to generate a thumbtack-shaped beampattern [25]. A multicarrier nonlinear frequency modulation FDA system based on pseudorandom frequency offset is presented in [26]. The system is simplified by omitting the optimization algorithm, and signal power can focus on different targets during desired time period. To attain a required radiation performance, such as the element placements and frequency offsets, the sparse FDA technology based on artificial bee colony optimizer is proposed in [27]. A new beampattern synthesis method for FDA radar is presented in [28], which incorporating an arc-tangent function-modulated frequency scheme.

The target angles and ranges are jointly estimated by the MUSIC method in [29]. A multitarget localization algorithm, which incorporating the coprime arrays and the coprime frequency increments for the sparse FDA radar, is proposed in [14]. To realize unambiguous frequency pattern identification, a simple criterion derived by an eigenvector-based method for target localization is presented in [30]. Unlike [14, 29, 30], a concept of double-pulse FDA radar with frequency increments being zero and nonzero is presented in [31] to achieve the target localization. Two steps are needed in this method to estimate the angle and range, respectively. In the first step, a zero frequency increment FDA, i.e., phased array, is utilized to estimate the target in the angle dimension. Then, the angle information is used as the prior knowledge and the range of the target is localized by another pulse with nonzero frequency increment.

In this paper, we propose a FDA transmit scheme called stepped frequency pulse FDA (SFP-FDA). In the proposed scheme, the SFP-FDA radar, like FDA radar, radiating a signal with a bit of frequency increment between the array elements in the first pulse, the difference is that there is another frequency interval between adjacent pulses. Due to the angle and range coupling of the beampattern, both the two dimensions of target cannot be estimated by FDA radar directly. Hence, the double-pulse FDA [31] which aims at angle-range localization of targets is proposed. The proposed system in this paper can be regarded as a kind of transformation or upgrading of [31], taking part of the amount of computation as a cost. It performs better in range resolution and multitarget localization. The data model and beamforming principle have been derived. The CRLB estimation of the proposed radar and the double-pulse radar are derived and simulated. And their performances are compared and analyzed.

The remaining sections are organized as follows. Section 2 gives some important details about FDA radar and phased array radar. Section 3 presents the mathematical principle and design technique of SFP-FDA radar. The SFP-FDA radar’s data model and the details of the targets’ localization in angle and range domains are proposed in Section 4. Section 5 derives the CRLB estimation performance. Next, Section 6 provides the numerical results and discussions. Finally, the conclusion is drawn in Section 7.

2. Frequency Diverse Array and Phased Array Radar

Consider the case where the frequency fed to th element of an element uniform linear array is [710]. where and are the signal operating frequency and the frequency interval between adjacent elements, respectively. The element spacing denoted as which is taken as where is the velocity of light and is the wavelength with respect to the operating frequency . The geometry of FDA is illustrated in Figure 1.

The signal radiated by the th element can be expressed as where is the complex weight related to th transmitted signal. Total electric field observed by a point target from far field can be written as where is the distance between target and th element; it can be given as , where and are the angle and range of the target, respectively.

The phase of th signal is

In comparison with the reference and th signal, the phase difference can be given as

The array factor is not just affected by the angle , the range , and frequency increment but also plays an important role. Since and , the 4th term of (6) is small and ignorable. So, (6) becomes

After rearranging terms, (4) can be expressed as

The steering vector can be given as where is the transpose operator. Define the element weighting vector as . In beamspace design sense, the common factor can also be ignored in the discussion; (8) simplifies to where is the conjugate transpose operator. Consider the case that , the transmit beampattern is given by which is the beampattern of FDA radar. Note that if , the beampattern becomes which is only angle dependent just like conventional phased array radar.

3. Mathematical Principle and Design Technique of SFP-FDA Radar

Consider an element FDA, the frequency transmitted on the th element in the first pulse of SFP-FDA radar is the same as that of FDA radar, and the frequency of the waveform radiated from each element is incremented by a small amount from pulse to pulse. As shown in Figure 2, the frequency in the th pulse of the th element can be given as where is the frequency interval between adjacent pulses, and is the number of the transmit pulses.

Adopting the signal of the th element, the electric field observed by can be taken as where is the pulse width, and is the complex weight of th pulse associated with th transmitted signal. Substituting and (13) into (14) yields

Assuming and , then let satisfied with . Due to the fact that is a periodic function with a period of , so can be approximately obtained. Furthermore, in the sense of amplitude, we can utilize as approximation. Then, (15) can be equivalently reformulated as where . Similarly, and ; also, we assumed that , and then the terms and are small and can be ignored. The common factor can also be ignored in the sense of beamspace design. And the corresponding transmit beampattern can be written as

Provided that the weights of signals’ pulses with the same order number in each element are equal. And and are defined as the signal weighting vector and the element weighting vector, respectively. Taking as the Kronecker product operator, (17) can be rewritten as where and

This implies that by optimizing and , the transmit beampattern can be schemed out.

Consider the case that , the beampattern (17) is then shown to be where

Note that when , we can get , that is, the expression of the conventional FDA. And if , (20) is equivalent to (12). Therefore, the phased array radar and FDA radar can be regarded as two special forms of SFP-FDA. Assume that , , and , where is the greatest common divisor (GCD) of and . For a given range and time, namely, and , we then have ; (20) will be equivalent to conventional FDA. Under these circumstances, both the FDA radar and the SFP-FDA radar have the same expression; hence, they possess the same angle resolution. According to [32, 33], the function can be defined as

Its 3 dB resolution can be calculated as

To calculate the range resolution, (11) can be simplified to

Then, the 3 dB range resolution of FDA can be given as

Similarly, (20) can be simplified to

Assuming that and , (26) can be equivalently reformulated as where . Taking and as an example, assuming , we can plot Figure 3.

Then, the 3 dB range resolution of SFP-FDA can be given as

The range resolution depends on , , , and . The second factor of (26) is satisfied with

Then, we will have , where and are the values of when (Figure 3) in SFP-FDA and FDA, respectively; we can also obtain identical conclusion by comparing (25) and (28). Therefore, SFP-FDA outperforms FDA in range resolution.

4. Target Localization in Angle-Range Dimensions

Consider the SFP-FDA with elements and pulses, the transmitted signal can be expressed as

The reflected signal is formed as a result of delay and attenuation of the transmitted signal. where is the complex reflection amplitudes from the target in ; and are the signal weighting vector and transmit steering vector of th array element, respectively. is the propagation delay in th array element. is the baseband waveform. We assumed

The demodulated baseband received signal can be written as where and are the received steering vector and normalized Gaussian noise vector, respectively. The steering vector of transmit-receive can be written as

The traditional nonadaptive receive beamforming is used, i.e., the weighting vector , where and denote the far field target’s angle and range, respectively. The received beampattern can be given as [34]

Note that if there are multiple targets such as three targets with the coordinates being , , and , the weighting vector will change to

According to analysis of Section 3 and Section 4, the data model of received signal can be written as where , and is the complex amplitudes.

For a nonadaptive beamformer, the weighting vector and the angle and range of the target can be estimated as

Owing to the decoupling transmit beampattern of the SFP-FDA radar, the angle and range of the target can be simultaneously estimated by utilizing only one step.

5. Performance Analysis

An analytical comparison of the performance of the proposed system and the double-pulse FDA radar in terms of CRLB is conducted; we rewrite (37) as where and are the equivalent receiver steering vector and the signal to noise ratio, respectively. And the mean and variance of are zero and . and are defined under the circumstances of SFP-FDA radar and the double-pulse FDA radar, respectively. A closed-form CRLB can be deducted approximately by arranging the steering vectors.

On the other hand, where

The mean and covariance of can be respectively expressed as

The Fisher information matrix (FIM) can be obtained by [35]. where

Then, the FIM of SFP-FDA can be written as where

By calculating the inverse matrix of FIM, the CRLB for the target’s estimation in angle and range domains can be obtained as where represents the element in the th row and th column of the matrix. The performance comparison of the two radars is presented in simulations.

6. Numerical Results

Consider a linear array of 16 elements with interelement spacing . The carrier frequency is , whereas the frequency interval between adjacent elements and pulses is . The number of pulses is , and pulse width is . The additive noise is modeled as complex Gaussian zero-mean spatially and temporally white random sequences with identical variance at each antenna element.

6.1. One Target

In the first part of this section, suppose that the target of interest is located in the angle and the slant range . The comparative target responses of the FDA radar and SFP-FDA are shown in Figure 4. In the FDA radar, owing to the angle-range coupling, the target cannot be localized effectively in both domains. Using the double-pulse FDA approach proposed in [31], the target angle can be estimated to be by using the first pulse which can be seen as phased array, and then the target range can be estimated to be by taking advantage of the second pulse with FDA, as shown in Figure 5. Similarly, using the method mentioned in Section 4, as shown in Figure 6, we can localize the target in angle-range dimensions as . Note that the pulse number , the element number , and the frequency increments and all affect the range resolution.

While compared with the double-pulse FDA radar, our proposed scheme has lower sidelobe and more concentrated energy in the mainlobe area. Besides, different from FDA radar, SFP-FDA radar has a range-angle decoupling beampattern.

SFP-FDA radar can estimate the coordinates of target by only one step, whereas the double-pulse FDA radar demands two steps.

By comparing Figure 5 with Figure 6, it can be observed that, despite the equal angle response of both the radars, the SFP-FDA outperforms the FDA in the range profile, i.e., the former has a narrower mainlobe and much lower sidelobe.

Then, we plotted Figure 7, which compared the 3 dB resolution performance of FDA and SFP-FDA. Clearly, the advantage of SFP-FDA in range dimension can be observed, which is consistent with the theoretical analysis in Section 3.

6.2. Two Targets with Different Angle and Different Range

Assuming that the two targets are located in coordinates and . Figure 8 gives the target response of FDA radar and SFP-FDA radar. Since the two targets have different angle, it can be noticed that due to the angle-range coupling, it is very difficult to separate the peaks in both angle and range domains from the FDA beamformer output.

By utilizing the first step of the double-pulse FDA radar, the angle of the targets can be received, namely, and . But when the second step is carried out, there are two peaks in each angle, as shown in Figure 9; thus, four targets have been obtained, and we cannot distinguish which ones are the real ones. The four targets can then be localized in angle-range dimensions as , , , and . Furthermore, by this means, we cannot accurately identify whether there are 2 targets, 3 targets, or even 4 targets.

Figures 10 and 11 give the target response of SFP-FDA radar, the two targets can be directly identified both in the angle and range dimensions, and then they can be localized in and .

It is evident that the two targets can be easily distinguished by the proposed scheme, which demonstrates the effectiveness of the SFP-FDA radar in resolving multiple targets. This simulation example shows that the SFP-FDA radar achieves better performance than the double-pulse FDA radar.

6.3. Estimation Performance

In this subsection, as compared to the double-pulse FDA radar, we simulated the performance of the proposed angle-range localization approach in terms of CRLB. Figure 12 shows the comparative performance of localization as a function of the SNR between the proposed scheme and the double-pulse FDA radar.

Figures 12(a) and 12(b) show the CRLB in the angle and range domain, respectively. It can be seen that the SFP-FDA radar can achieve a better estimation performance in both angle and range dimension, even when the number of its pulses is the same as that of the double-pulse FDA radar, that is, when . Moreover, it is evident that the performance of both in the angle and range domains is improved as the SNR and the number of pulses increase. As shown in Figure 12, the SFP-FDA with performs better than the SFP-FDA with . It is clear that the SFP-FDA radar achieves better estimation performance at the cost of much higher computation complexity.

7. Conclusion

Since the beampattern produced by conventional FDA is angle-range coupling, a novel scheme called stepped frequency pulse FDA radar is proposed in this paper, which can decouple the angle and range responses in the transmit beamforming. Only one step procedure is needed to localize the targets in angles and ranges simultaneously, and the range resolution of targets outperforms that of double-pulse FDA radar. The proposed system has a stronger ability in identifying the multiple targets with different angles and different ranges. It has also been shown that the proposed radar achieves a better CRLB performance than the double-pulse FDA radar, though at the expense of a certain amount of computational complexity. Certainly, further investigations will consider the Doppler of targets. With new application potentials, combining with netted radar networks and MIMO radar technology can effectively improve the radar performance, and it can also be investigated in a future work.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

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

Acknowledgments

This work was supported by the National Natural Science Foundation of China (grant no. 61601502).