Xiang Zhe Chen Baixiao

(National Laboratory of Radar Signal Processing, Xidian University, Xi’an 710071, China)

1 Introduction

Frequency Diverse Arrays (FDA) radar has recently drawn much attention among the researchers. FDA differs from the traditional phased array by using a small frequency increment across the array elements, which results in a rangeangle-dependent beam pattern. FDA radar system is first proposed in Ref. [1]. In FDA radar, a uniform interelement frequency offset is applied across the array elements. FDA radar with uniform small and large frequency offset frequency has been investigated in Refs. [1–7]. Small frequency offset has been exploited to generate range-dependent beampattern, while large frequency offset can get independent echoes from the target.

Unlike the phased array, the range-angle dependency of the FDA beampattern allows the radar system to focus the transmit energy in a desired range-angle space. This unique feature of FDA helps to suppress the range-dependent interferences[8]and increases the received SINR consequently. Especially for the mainlobe interference and clutter, the FDA can achieve a significant improvement in SINR against the phased array because the FDA provides the increased Degrees Of Freedom (DOFs) in range domain.However, the FDA beampattern is shown to be periodic in range and time[2], which goes to maximum at multiple time and range values. With this multiple-maximum beampattern, the resulting SINR will be deteriorated when the interferers are located at any of the maxima. To improve SINR, FDA with Time-Dependent Frequency Offset (TDFO-FDA) was proposed to achieve a time-independent beampattern at the target location[9]. Nevertheless, the proposed beampattern is still periodic in range which will result in the loss of SINR. A nonuniformly spaced linear FDA with linear incremental frequency increment has also been studied in Ref. [10], and a nonrepeating beampattern has been obtained for range-angle imaging of targets. A uniformly spaced linear FDA with Logarithmically (Log-FDA) increasing frequency offset is proposed in Ref. [11]. The proposed strategy provides a nonperiodic beampattern with the single-maximum in space. In Ref. [12], the beampattern of FDA who transmits the pulsed signal has been studied.Lately, few more publications have done some work in decoupling the range-angle dependent beampattern of FDA[13–16]. All these papers only address the properties of the FDA beampattern,and they do not study the common rule for the FDA configuration to form a single-maximum transmit beampattern.

With the pioneer work on FDA radar, we aim to decouple the range and angle in the beampattern and provide a nonperiodic beampattern with the single-maximum in the illuminated range-angle space. In this paper, we propose a basic criteria for the FDA configuration, in which the element spacing and frequency increment are configurable, to form a single-maximum beampattern through mathematical analysis. This singlemaximum beampattern, unlike the multiple-maxima beampattern, can help to further suppress range-dependent interferences, causing improved SINR and increased detect ability. The proposed rule for the FDA configuration will be helpful to design the FDA.

The rest of the paper is organized as follows.In Section 2, the basic FDA model has been described and the basic criterion is derived for the FDA configuration to form a single-maximum beampattern through mathematical analysis.Moreover, several specific conditions are introduced. In Section 3, the beampattern has been plotted for the specific conditions discussed in Section 2. Finally, in Section 4 we conclude the paper.

2 Design and Mathematical Analysis of FDA

2.1 System description

Consider an array ofMtransmit elements,we assume that the waveform radiated from each antenna element is identical with a frequency increment, as shown in Fig. 1. The radiated frequency from them-th element is

where Δfmis the frequency increment ofm-th element with reference to the carrier frequencyf0.Specifically, Δf0= 0.

Considering a given far-field point, the phase of the signal transmitted by them-th element can be represented by

where c andrmare the speed of light and the distance between them-th element and the observed point, respectively.

The range difference between individual elements is approximated by

whereθ0is the desired angle,dmis the spacing between them-th element and the first element.Specifically,d0= 0.

So the phase difference between them-th element and the first element is

In Eq. (4), the third term is important because it shows that the FDA radiation pattern depends on both the range and the frequency increment. Taking the first element as the reference for the array, the steering vector can be expressed as

where [·]Tdenotes the transpose operator.

wheredtis the duty cycle, which is usually small.Eq. (7) holds whenis also small.

So under the condition Eq. (7), the timetcan be neglected. So Eq. (5) can be simplified as

2.2 Transmit beampattern analysis

Throughout this paper, we assume a narrowband system where the propagation delays manifest as phase shifts to the transmitted signals and Eq. (7) is satisfied. To steer the maximum at an expected target location (θ0,r0), the complex weights are configured asa(θ0,r0), so the transmit beampattern can be expressed

where [·]Hdenotes the conjugate transpose operator.

It is easy to see that the beam direction will vary as a function of the range and angle, which means the beampattern is range-angle dependent.Since the beampattern is coupled in the range and angle, the target’s range and angle cannot be estimated directly by the FDA beamformer output.Note that the beampattern is also related to Δfmanddm, so the desired single-maximum beampattern can be obtained by setting the proper Δfmanddm.

In Eq. (9), whenm= 0, the exponential term is equal to 1. To obtain the maximum value of the beampattern, the exponential terms should be all equal to 1 form=1，2，···，M-1. So the phase of the exponential term should be the integral multiple ofwhich can be expressed as

whereLmis an integer,e.g.Lm= 0, ±1,···,m=1，2，···，M-1.

Sincedm(sinθ-sinθ0)?r0, the termin Eq. (8) can be neglected. The curves formed by Eq. (11) will be called as range-angle distribution curves throughout the paper. Then Eq. (11) can be approximately expressed as

Rewrite Eq. (12) into matrix form as

To decouple the range and angle, the beampattern should have the unique maximum point in the range-angle distribution diagram, which means the Eq. (13) has the unique solution (θ0,r0). The necessary and sufficient condition of that the Eq. (13) has the unique solution is

where rank (·) is the rank of a matrix,the Eq. (13) has infinite solutions, corresponding to the conventional FDA condition, which will be discussed in detail later.

To satisfy rank (A) = 2,dm/=PΔfm,Pis a constant.

In the following, under the condition ofdm/=PΔfm，Lm=0，±1,···, we make a summary with different parameter configurations:

(1) whenLm= 0, the Eq. (13) has the unique solution (θ0,r0);

(2) whenLm=Qdm/=0,Lm/=SΔfm, the Eq. (10) has the solutions. When|sin(θ0)-kQλ0|≤1, the angle grating lobes will occur at angle arcsin(sin(θ0)-kQλ0) in the beampattern. Otherwise, the Eq. (13) has no solution, resulting in no angle grating lobes;

(3) whenLm=SΔfm/=0，Lm/=Qdm, the Eq.(13) has the (θ0，r0+Skc),k=0，±1，±2，···,which means the range grating lobes will occur atr0+Skcin the beampattern.

Note that in array theory, when the adjacent element spacing is less than half the wavelength, the angle grating lobes will never appear. IfQλ0=±1 andθ0=0°, where the element spacing is0, then the grating lobes will occur at angle ±90°. ButLm=SΔfmcan always be satisfied sinceLmis an integer whose range is[-∞，+∞]. The range grating lobes will always occur at ranger0+Skc in the beampattern. The position of the range grating lobe changes with differentS. For example,S= 0.01, the distance between the grating lobe and the mainlobe is 3000 km. So ifr0±Sc is out of the illuminated range space [Rmin，Rmax], the beampattern has a singlemaximum point (θ0,r0) in the illuminated range space, which means the range and angle are decoupled.

So we can conclude that the beampattern of the FDA is always range-periodic, the grating lobes will always occur at ranger0+Skc,k=0，±1，±2，···. To obtain the single-maximum beampattern in the illuminated range space[Rmin，Rmax], the designing criteria for the FDA isdm/=PΔfm, and 2Sc＞Rmax-Rmin,|sin(θ0)±Qλ0|＞1,k=0，±1，±2，···,Lm=0，±1，···,m=1，2，···,M– 1,Pis a constant,Q,Sare the minimum non-zero constants to satisfy the equationsSΔfm/=0andLm=Qdm/=0.

Once the range and angle is decoupled, the target’s range and angle can be estimated directly by the FDA beamformer output. Also the 2-dementional MUSIC algorithm[16]for estimating the target’s range and angle can be used as well.

2.3 Specific examples

For the conventional FDA, Δfm=mΔf,dm=md, Δfanddare configurable parameters to control the frequency increment and the element spacing. WhenLm=nm,n=0，±1，±2，···,m=1，2，···，M-1, we can getUnder this circumstance, the Eq. (13)has infinite solutions. The solutions form the range-angle dependent curves, which have expressions as:

In Eq. (15), the expression is no longer related tom, which means the range-angle curve for different element coincides with each other, as depicted in Fig. 2(a). In the range-angle distribution diagram, the curve is periodic in range,and the range difference between the adjacent curves is c/Δf. The corresponding beampattern is depicted as Fig. 3(a).

For the Expf-FDA, whose frequency increment is exponentially increased,Δfm=(bm-1)Δf,dm=md.bis a configurable constant.The range solution rises whenbgets larger. The range-angle distribution curves and the beampattern are depicted in Fig. 2(b) and Fig. 3(b), respectively.

Likewise, for the Logd-FDA, whose element spacing is logarithmically increased, Δfm=mΔf,dm=log(m+1)d.WhenLm=nm,n=0，±1，±2，···,m=1，2，···，M-1, we can get rankso the range grating lobes will occur in the beampattern as depicted in Fig. 3(c),and the range difference between the grating lobe and the mainlobe is c/Δf. The corresponding range-angle distribution diagram is depicted in Fig. 2(c).

For another kind of FDA, called Expf-Logd-FDA, where the frequency increment is exponentially increased and the element spacing is logarithmically increased, Δfm=(bm-1)Δf,dm=log(m+1)d. The range-angle distribution curves and beampattern are depicted in Fig. 2(d)and Fig. 3(d), respectively.

3 Simulations, Results, and Discussions

Bampattern expressed in Eq. (9) and the range-angle distribution curves expressed in Eq.(11) were simulated and plotted for different kinds of FDA discussed in Section 2. The results are discussed and compared with the different kinds of FDA. The illuminated range space is (0 km,800 km]. To generate these plots, the values of the configurable parameters have been taken as listed in Tab. 1. To avoid angle grating lobes, the parameterdis less than half the wavelength.

Tab. 1 Parameters for simulations

The range-angle distribution curves are depicted in Fig. 2. The curves with different color represent the range-angle distribution for different elements except the reference element in the FDA. It is easy to see that the range-angle distribution of the elements are the same in conventional FDA, and the range difference between the adjacent curves is c/Δf= 300 km, which is consistent with the analysis in Section 2. For the Expf-FDA and Expf-Logd-FDA,Δfm=(1.4m-1)Δf, the range between the first grating lobe and the mainlobe is 3e6 km, which is not in the illuminated space. So the curves of the elements in the range-angle distribution diagram have the single intersection point, which means that the beampattern has the single-maximum point in the illuminated range space. But for the Logd-FDA, whenLm=nm, the distance between the grating lobe and the mainlobe is c/Δf= 300 km according to the analysis in Section 2. So in the range-angle distribution diagram, the curves have 3 intersection points, the range difference between adjacent intersection points is 300 km.

The beampatterns are depicted in Fig. 3.Similar to the range-angle curve in Fig. 2, the beampattern of the conventional FDA is a rangeangle-dependent beampattern. For the Expf-FDA and Expf-Logd-FDA, the beampatterns have a single-maximum point in the illuminated space corresponding to the single intersection point in the range-angle distribution diagram. For the Logd-FDA, the beampattern has 3 range grating lobes corresponding to the 3 intersection points in the range-angle distribution diagram. For the latter 3 FDAs, the contour of the beampattern is an ellipse, which is because the range-angle distribution curves are tightly distributed. The direction of the major axis of the ellipse is corresponding to the most tightness distribution direction of the range-angle distribution curves.

In the meanwhile, we analyze the range and angle solutions of different FDAs in Fig. 4. The Expf-FDA and the Exp-Logd-FDA have the same range solution, and the range solution of conventional FDA and Logd-FDA are the same. That is because the bandwidth across the whole array is the same for each pair. The same situation occurs for the angle solution, if the arrays have the same array aperture, they have the same angle solutions.

We choose the target position at (400 km,20°), and the transmit beampattern is depicted in Fig. 5. Similar to Fig. 3, we can see that the single maximum beam is formed at the target position.

4 Conclusions

In this paper, we propose a basic criteria for the FDA configuration to provide a single-maximum beampattern in the illuminated space. The single-maximum beampattern can be generated by configuring the element spacing and frequency increment of the FDA. Through the analysis, we can find out that the beampattern of FDA is always range-periodic. Results show that a singlemaximum beampattern can be generated with the corresponding criteria by choosing the proper frequency increment. As this paper describes the transmitter only, designing an appropriate receiver is our future work.

[1]Antonik P, Wicks M C, and Griffiths H D. Range dependent beamforming using element level waveform diversity[C]. International Waveform Diversity Design Conference, Las Vegas, NV, USA, Jan. 2006: 22–27.

[2]Secmen M, Demir S, and Hizal A. Frequency diverse array antenna with periodic time modulated pattern in range and angle[C]. IEEE Radar Conference, Boston, MA, USA, Apr.2007: 427–430.

[3]Antonik P, Wicks M C, and Griffiths H D. Multi-mission,multi-mode waveform diversity[C]. IEEE Radar Conference,Verona, NY, USA, Apr. 2006: 580–582.

[4]Zhuang L and Liu X Z. Precisely beam steering for frequency diverse arrays based on frequency offset selection[C]. International Radar Conference, Bordeaux,France, Oct. 2009: 1–4.

[5]Chen Y G, Li Y T, and Wu Y H. Research on the linear frequency diverse array performance[C]. IEEE International Conference on Signal Processing, Beijing, China, Oct. 2010:2324–2327.

[6]Sammartino P F, Backer C J, Griffiths H D. Frequency diverse MIMO techniques for radar[J].IEEE Transactions on Aerospace and Electronic Systems, 2013, 49(1): 201–222.DOI: 10.1109/TAES.2013.6404099.

[7]Xu J W, Liao G S, and Zhu S Q. Receive beamforming of frequency diverse array radar systems[C]. 31th URSI General Assembly and Scientific Symposium (URSIGASS),Beijing, China, Aug. 2014: 1–5.

[8]Wang W Q. Range-angle dependent transmit beampattern synthesis for linear frequency diverse arrays[J].IEEE Transactions on Antennas Propagation, 2014, 61(8):4073–4081.

[9]Khan W, Qureshi I M. Frequency diverse array radar with time dependent frequency offset[J].IEEE Antennas and Wireless Propagation Letters, 2014, 13: 758–761. DOI:10.1109/LAWP.2014.2315215.

[10]Wang W Q, So H C, Shao H. Nonuniform frequency diverse array for range-angle imaging of targets[J].IEEE Sensors Journal, 2014, 14(8): 2469–2476. DOI: 10.1109/JSEN.2014.2304720.

[11]Khan W, Qureshi I M, Saeed S. Frequency diverse array radarwith logarithmically increasing frequency offset[J].IEEE Antennas and Wireless Propagation Letters, 2015, 14:499–502. DOI: 10.1109/LAWP.2014.2368977.

[12]Xu Y H, Shi X W, Xu J W. Range-angle-dependent beamforming of pulsed frequency diverse array[J].IEEE Transactions on Antennas Propagation, 2015, 63(7):3262–2367. DOI: 10.1109/TAP.2015.2423698.

[13]Gao K, Wang W Q. Decoupled frequency diverse array range-angle-dependent beampattern synthesis using nonlinearly increasing frequency offsets[J].IET Microwaves,Antennas&Propagation, 2016, 10(8): 880–884.

[14]Gao K, Wang W Q. Transmit beamspace design for multicarrier frequency diverse array sensor[J].IEEE Sensors Journal, 2016, 16(14): 5709–5714. DOI: 10.1109/JSEN.2016.2573379.

[15]Shao H, Dai J, Xiong J. Dot-shaped range-angle beampattern synthesis for frequency diverse array[J].IEEE Antennas and Wireless Propagation(published online),2016. DOI: 10.1109/LAWP.2016.2527818.

[16]Schmidt R O. Multiple emitter location and signal parameter estimation[J].IEEE Transactions on Antennas Propagation, 1986, 34(1): 276–280.