Impacts of frequency increment errors on frequency diverse array beampattern

Abstract

Different from conventional phased array, which provides only angle-dependent beampattern, frequency diverse array (FDA) employs a small frequency increment across the antenna elements and thus results in a range angle-dependent beampattern. However, due to imperfect electronic devices, it is difficult to ensure accurate frequency increments, and consequently, the array performance will be degraded by unavoidable frequency increment errors. In this paper, we investigate the impacts of frequency increment errors on FDA beampattern. We derive the beampattern errors caused by deterministic frequency increment errors. For stochastic frequency increment errors, the corresponding upper and lower bounds of FDA beampattern error are derived. They are verified by numerical results. Furthermore, the statistical characteristics of FDA beampattern with random frequency increment errors, which obey Gaussian distribution and uniform distribution, are also investigated.

Gao et al. EURASIP Journal on Advances in Signal Processing (2015) 2015:34 DOI 10.1186/s13634-015-0217-y

RESEARCH Open Access

Impacts of frequency increment errors on frequency diverse array beampattern

Kuandong Gao*, Hui Chen, Huaizong Shao, Jingye Cai and Wen-Qin Wang

Abstract

Different from conventional phased array, which provides only angle-dependent beampattern, frequency diverse array (FDA) employs a small frequency increment across the antenna elements and thus results in a range angle-dependent beampattern. However, due to imperfect electronic devices, it is difficult to ensure accurate frequency increments, and consequently, the array performance will be degraded by unavoidable frequency increment errors. In this paper, we investigate the impacts of frequency increment errors on FDA beampattern.

We derive the beampattern errors caused by deterministic frequency increment errors. For stochastic frequency increment errors, the corresponding upper and lower bounds of FDA beampattern error are derived. They are verified by numerical results. Furthermore, the statistical characteristics of FDA beampattern with random frequency increment errors, which obey Gaussian distribution and uniform distribution, are also investigated.

Keywords: Frequency diverse array (FDA); Frequency increment errors; Beampattern error; FDA beampattern; FDA radar

1 Introduction

Beampattern is widely used to assess the performance of phased arrays [1]. However, a limitation of phased array is that the beam steering is fixed at one angle for all the ranges [2]. Recently, a flexible array called frequency diverse array (FDA) has been proposed [3-5]. Different from phased array, a small frequency increment, as compared to the carrier frequency, is applied across the array elements [6]. This small frequency increment results in a range angle-dependent beampattern [7,8]. Several investigations have been carried on FDA radars. The time and angle periodicity of FDA beampattern was analyzed in [9].A linear FDA was proposed in [10] for forward-looking radar ground moving target indication. The multipath characteristics of FDA radar over a ground plane were investigated and compared with phased array in [11].FDA radar full-wave simulation and implementation with linear frequency-modulated continuous waveform were presented in [12,13]. In [14], we have investigated the FDA

*Correspondence: mailto: [email protected]

Web End [email protected]

Equal contributorsSchool of Communication and Information Engineering, University of Electronic Science and Technology of China, 2008, Road XiYuan, Chengdu, 611731 Sichuan, China

Cramr-Rao lower bounds for estimating direction, range, and velocity. More recent work about the applications of FDA in MIMO radars can be found in [15-19].

Perfect frequency increments are often assumed in existing literatures [20]. However, in an actual array system, there will have imperfect errors including element position errors, mutual coupling, phase errors, and frequency increment errors [21-24]. Some results have been reported about the impacts of element position error, mutual coupling and phase error on beampattern, and direction-of-arrival (DOA) estimation performances. Since FDA beampattern has similar properties with conventional phased array for the impacts of element position errors, mutual coupling, and phase errors, this paper considers only the impacts of frequency increment errors on FDA beampattern. Since FDA beampattern is dependent on the angle and range, it has a potential for target localization , which is different from traditional time-of-arrival (TOA) and angle-of-arrival (AOA)-based localization [25-28]. The contributions can be summarized as follows.(i) More tighter bounds of FDA beampattern deviation are derived. (ii) Statistical analysis of FDA beampattern in terms of expectation value, variance, and probability density function (PDF) are provided.

2015 Gao et al.; licensee Springer. This is an Open Access article distributed under the terms of the Creative Commons

Attribution License http://creativecommons.org/licenses/by/4.0

Web End =(http://creativecommons.org/licenses/by/4.0) , which permits unrestricted use, distribution, and reproduction

in any medium, provided the original work is properly credited.

Gao et al. EURASIP Journal on Advances in Signal Processing (2015) 2015:34 Page 2 of 12

The rest of this paper is organized as follows. Section 2 formulates the data models of FDA radar without and with frequency increment errors, respectively. Section 3 analyzes the impacts of deterministic frequency increment errors on FDA beampattern. Thereafter, the impacts of random frequency increment errors are investigated in Section 4. Finally, simulation results are provided in Section 5, and conclusions are drawn in Section 6.

2 FDA beampattern without frequency increment errors

Suppose an N-element uniform linear FDA with inter-element spacing denoted as d. The radiated frequency from the nth element is as follows:

fn = f0 + n f , n = 0, 1, 2, . . . , N 1 (1) where f0 and f are the carrier frequency and frequency increment, respectively. Taking the first element as the reference for the array, under far-field condition, one might express the direct wave component of the electric field emitted from the FDA at the observation point (, r) as [17]:

A(, r, t) =

N1

section as a special case. By definition, the elements are excited with uniform amplitude, but they are allowed to have a phase progression across the array. These specifications translate to the following expressions for dn and an:

dn = nd and an = ejna (5)

where a stands for the phase progression. Submitting Equations 1 and 5 into Equation 4 yields:

A(, r, t) = ( |w0 )

N1

t

r+nd sin

c0

n=0ejna ej2

(f0+n f )

r

N1

=

( |w0 )r ej

0

n=0ejna ej2n f t

r+nd sin

c0

j2f0 nd sin

c0

(6)

where 0 = 2f0t 2f0 rc0 . For notation convenience, we define

= ( |w0 ) ej2f0tj2f0

rc0 , Equation 6 can then

be rewritten as:

A(, r, t) =

N1

n=0ejnaej2n t f r+nd sinc0 f f0 d sinc0

(7)

Since ndsin r, Equation 7 can be reformulated as:

A(, r, t)

r

n=0ann ( |wn ) ejwn

t r+dn sinc0

r (2)

where N is the number of FDA elements, an represents the complex excitation coefficient for the nth element, n ( |wn ) stands for the far-field vector radiation pattern for the nth element at range r and angular frequency wn =

2fn, c0 is the light speed, dn is the element position of the nth element reference relative to the first element, and t is the time parameter. In accordance with the far-field

assumption, ejwn

t r+dn sinc0 /r corresponds to the delayed carrier with free space loss.

To interpret the effect of frequency diversity within the scope of an array factor, we should factor the vector element pattern out of Equation 2. This can indeed be done under certain conditions. Assuming all elements in the FDA are identical, we can eliminate the frequency dependence in the element factor. Since r dn, we have:

n(|wn) (|w0). (3) where w0 is the carrier angular frequency. So Equation 2 can be rewritten as:

A(, r, t) = ( |w0 )

N1

N1

r

n=0ejnaej2n t f rc0 f f0 d sinc0

(8)

=

(, r, t)

where

w =

1 eja . . . ej(N1)a

T (9)

and

v(, r, t) =

1 ej2

t f rc0 f f0 d sinc0 . . .

ej2(N1)

t f rc0 f f0 d sinc0 T

.

n=0an ejwn

t

r . (4)

Further simplification becomes possible by considering particular FDA arrangements that are simple to handle and yet able to provide valuable insight. A uniform linear FDA utilizing discrete, linear frequency increments is such a practical configuration, and it is examined in this

(10)

with [ ]T and [ ]H being the transpose operator and hermitian transpose operator, respectively. For simplicity,

= 1 is assumed in the following discussions. Figure 1 shows the ideal linear uniform FDA beampattern with t = 1/ f .

2.1 FDA beampattern with frequency increment errors

Suppose the frequency increment errors being = [ 0 1 . . . N1 ], Equation 8 can be rewritten as:

A(, r) =

r+dn sin

c0

wH v(, r, t)r (11)

Gao et al. EURASIP Journal on Advances in Signal Processing (2015) 2015:34 Page 3 of 12

Figure 1 Ideal two-dimension FDA beampattern with N = 16,

f = 30 KHz, t = 1/(3 104), a = 0, and w = 1N1.

where

v(, r, t) =

According to Equation 13, it can be known that:

C+(r, t) = diag

j2

(N )rc0 + t(N )

ej2 (N )rc0 +t(N ) + O()

diag

j2

(N )rc0 + t(N )

ej2 (N )rc0 +t(N ) .

(15)

where N =[ 0, 1, 2, . . . , N 1] denotes element number of FDA radar, and O() represents for the high order terms about , which can be ignored for Equation 15. Using Equation 15, we can rewrite the beampattern (13) as:

A(, r, t) =

wHC(r, t)v(, r, t)

r

1 ej2 t( f +1)rc0 ( f +1)f0d sin c0

. . .

wHv(, r, t) + wHC+(r, t)v(, r, t) r

= A(, r, t) +

=

ej2(N1)

t( f +N1)

wHC+(r, t)v(, r, t) r

= A(, r, t) + A(r, t)

(16)

rc0 ( f +N1)f0

d sin c0

T

H C+(r,t)v(,r,t)

r denotes the FDA beam-pattern deviation. Figure 2 shows the FDA beampattern with frequency increment errors. The array parameters are the same as that for Figure 1. Due to the frequency increment errors, the FDA beampattern sidelobes are different from Figure 1. Stronger sidelobe peaks will make the FDA energy scattering and consequently degrade the array performance.

3 Impacts of deterministic frequency increment errors

Firstly, we consider uniform frequency increment errors,i.e., = [0 1 . . . N1] = [0 . . . ]. In this case, Equation 13 can be rewritten as:

A(, r)=

wHdiag ( v(r, t)) v(, r, t) r

=

N

=

1 ej2 t f rc0 f f0d sin c0

ej2 t1rc0 1 . . .

ej2(N1)

where A(r, t) = w

t f

rc0 f f0

d sin c0

ej2

N1)(tN1rc0 N1 T.

(12)

Then Equation 11 can be reformulated as:

A(, r) =

wH v(, r, t) r

=

wH(v(, r, t) v(r, t)) r

(13)

=

wHdiag( v(r, t))v(, r, t) r

=

wHC(r, t)v(, r, t)

r

where

v(r, t) =

1 ej2 t rc0 . . . ej2(N1) t rc0 T

and

C = diag ( v(r, t))

and denotes Hadamard product. Note that diag( v(r, t)) denotes a diagonal matrix with v(r, t) being its diagonal elements. For small frequency increment errors, Taylor series expansion about C is performed as follows:

C(r, t) = I + C+(r, t). (14)

n=1 ej2(n1)

t( f +) rc0 ( f +)f0 d sinc0

r

sin

N 2

t

f + rc0 ( f + ) f0 d sinc0

sin

=

k r

t( f +) rc0 ( f +)f0 d sinc0 2

sin

N 2

f

t+ f t

fc0

r+ f r

f0 d sinc0

=

k r

sin

f

t+ f t fc0 r+ f r f0 d sinc0

2

(17)

Gao et al. EURASIP Journal on Advances in Signal Processing (2015) 2015:34 Page 4 of 12

Figure 2 Two-dimension FDA beampattern with frequency increment error n and n N(0, 2) with = 500.

where k = ejN

t( f +) rc0 ( f +)f0 d sinc0 and w = 1N1.

Equation 17 arrives the maximum value when:

f t + f t f c0

and the angle error is as follows:

= arcsin

c0 2 f f0d t

r + f r f0 d sinc0 =2m

m = 0, 1, 2, . . . .

(18)

If no frequency increment error exists, the FDA main-lobe will pass the location

0, c0 f

. (22)

If the range error and angle error are required to be smaller than 50 m and 0.5, respectively, the frequency increment error should be smaller than 150 Hz.

Figure 3a,b shows the FDA beampattern with uniform frequency increment errors , frequency increment f =

30 KHz, and t = 1/ f . It can be noticed that the phase error for = 0.5 is / f = 5 103. Note that since the range attenuation exists, the range error has a small deviation from the calculation of Equation 21.

4 Impacts of stochastic frequency increment errors

In this section, we investigate the statistical properties of FDA beampattern with random frequency increment errors.

4.1 Error boundary of FDA beampattern

Suppose the frequency increment errors are random. In this case, it is not possible to derive a deterministic beam-pattern expression like Equation 17. Here, we are interested in the maximum beampattern deviation, namely:

max

t,r, | A(t, r, )| . (23) According to Equation 16 and using the Cauchy-

Schwarz inequality, we have:

max

t,r, | A(, r, t)|

1r wmaxC+ 2 v

f . However, when there are frequency increment errors, the location

0, c0 f

may be not at the beampattern maximum point. According to Equation 18, the corresponding range error is as follows:

r =

f r. (19)

Since the phase error caused by time error can be equivalently regarded as range error and angle error, Equation 18 can be reformulated as:

f t + f t f c0

at t =

1

r + f r f0 d sin c0

= ft +

f t

f c0

r + f r f0 d sin c0

= ft

f c0

r + f r c02 f 2 t f0 d c0

sin c0 2 f f0d t

= 2m, m = 0, 1, 2, . . .

(20)

Therefore, the corresponding range error is as follows:

r =

f r

c02 f 2 t (21)

(, r) 2

Gao et al. EURASIP Journal on Advances in Signal Processing (2015) 2015:34 Page 5 of 12

50 0 50

The maximum deviation, in a linear scale, is the same over the whole beampattern. When only C+ (and thus C) is known, or when the influence of unknown factors on C+, e.g., aging or temperature, comes into play, it makes sense to consider as random matrix, with some statistical model for the entries of C+. We then have to calculate:

max C+ 2 (26)

for the random matrix C+ to find an upper bound.

max C+ 2 = max 2

rc0 + t diag(N )

2

a

0

ideal / f=5e3 / f=2e3 / f=1e3

5

10

Beampattern (dB)

15

20

1.5

1.55

1.6

1.65

1.7

0.5 0

0.5 1

25

30

35

rc0 + t

N 2

40

= max 2

Angle (Degree)

0

= 2

rc0 + t

b

N 2 |max| .

(27)

When applying this result to an array with uniform weighting wmax = 1, Equation 24 leads to:

Amax =

ideal /

f=5e3 /

f=2e3 /

f=1e3

5

Beampattern (dB)

10

15

N 2 |max| N. (28)

For a constant time t = 1/ f , the time error caused by frequency increment error will become amplitude error of

FDA beampattern, as shown follows:

Amax =

1.45

1.5

1.55

1.6

1.65

9940 9960

9980 10000

10020

20

2 r

rc0 + t

25

30

N 2 |max| N. (29)

It can be noticed from Equation 29 that the errors caused by time and by range counterbalance the total beampattern error. Therefore, we delete 1/ f . Then, Equation 29 can be reformulated as:

Amax

35 0.7 0.8 0.9 1 1.1 1.2 1.3 x 104

Range(m)

2 r

rc0 +

1 f

Figure 3 FDA beampattern profiles with frequency increment errors.(a) In angle dimension. (b) In range dimension.

=

1r wmaxC+ 2

N

N

2 r

r c0

1r |wmax| C+ 2

N 2 |max| N

(24)

where || and 2 denote the absolute value and Euclidean norm, respectively, and wmax denotes the maximum element of w. Note that compared with [29], |wmax|

|w|2, the bound expressed in Equation 24 is much tighter. According to Equation 16, C+ can be rewritten as:

C+ 2 =

diag

j2

(30)

Equation 30 gives FDA beampattern error bound which indicates its worst case. Since it is caused by the frequency increment errors, the maximum device errors which produce the FDA frequency increment are regulated by the bound. This gives guideline about device selections and predicts the possible FDA beampattern derivation bound.

4.2 Statistical properties of FDA beampattern

In most cases, we know the statistical properties of frequency increment errors and would like to derive the respective properties of impacted beampattern. In this subsection, we will give the formulas of expectation value, variance, and PDF about FDA beampattern based on different PDFs of the frequency increment errors.

= Amax

=

2c0 N 2 |

max| N.

(N )rc0 + t(N )

ej2((N )rc0 +t(N ))

2

(25)

=

diag

j2

rc0 + t (N )

2

= 2

rc0 + t

diag(N

)

2.

Gao et al. EURASIP Journal on Advances in Signal Processing (2015) 2015:34 Page 6 of 12

4.2.1 Expected value

Using Equations 8 and 13, it is straightforward to show that:

E { A(, r, t)} = E

and maximum values max and max. We can get the expected value of nth element in diag(C) as follows:

E {Cnn(r, t)} = E

ej2(n1) t rc0

wHv(, r, t)

= wHE {C(r, t)} v(, r, t) wHv(, r, t)(31)

where E{} denotes the expected value.

Assume that the frequency increment error n of the nth element satisfies the Gaussian statistical model n

N(0, 2) and the PDF is as follows:

fn() =

1 2 e

wHC(r, t)v(, r, t)

max

1

=

2max d

(36)

By using Equation 36 onto 31, the expected value of FDA beampattern deviation with uniform distribution frequency increment error can be reformulated as:

E { A(, r, t)} = wHE {C(r, t)} v(, r, t) wHv(, r, t)

= wH

max

ej2(n1)

t rc0

sin(2(n 1)

t r c0

max)

.

= 2(n 1)

t r c0

max

2

22 . (32)

Therefore, the expected value of the nth element in diag(C) can be calculated as [30]:

E {Cnn(r, t)} = E

ej2(n1) t rc0

diag (bu) I

v(, r, t)

(37)

where

1,

sin

2

t r c0

max

bu =

, . . . ,

2

+ ej2(n1)

t rc0 e

t r c0

1 2

max

=

2

22 d

(38)

Similarly, the expected value of the FDA beampattern deviation is dependent on the maximum value of frequency increment error and not 0 for the uniform distribution model. From Equations 34 and 37, the expected values of FDA beampattern have range offset and angle offset, which is caused by time dependent error on the two distributions.

4.2.2 Variance

The beampattern deviation variance var { A(, r, t)} equals to the beampattern variance:

var { A(, r, t)} = var {A(, r, t)}

= var

sin

2(N 1)

t r c0

max

2(n1)2

t rc0 2

.

2

.

= e2

2(N 1)

t r c0

max

(33)

where Cnn denotes the nth element of diag(C). When applying Equation 33 onto 31, we can get the expected value of FDA beampattern deviation with Gaussian distribution frequency increment errors as follows:

E { A(, r, t)} = wHE {C(r, t)} v(, r, t) wHv(, r, t)

= wH

diag

bg I v(, r, t)(34)

where

wHC(r, t)v(, r, t)

.

(35)

and I denotes N N unit matrix. In this case, the expected value of the FDA beampattern deviation is dependent on the variance of frequency increment errors and is not 0 for the Gaussian distribution model.

For another case, we assume that n is uniform distribution, i.e., n u(max, max) with minimum

bg =

1, e22

t rc0 22, . . . , e22(N1)2 t rc0 2 2

. (39)

As we deal with vectors and matrices, we utilize the covariance matrix cov{}, and var{} = cov{} for any scalar. We thus have:

var { A(, r, t)} = cov

wHC(r, t)v(, r, t)

(40)

Using the Kronecker product , the vectorization transformation vec() and the identity:vec {ABC} = CT Avec {B} . (41)

Since

wHC(r, t)v(, r, t) = vT(, r, t) wHvec {C(r, t)} (42)

Gao et al. EURASIP Journal on Advances in Signal Processing (2015) 2015:34 Page 7 of 12

we use c = vec {C(r, t)} and the vector t(, r, t) = vT(, r, t) wH, which utilizing:

cov {AX} = Acov {X} AH (43)

to yield:

var { A(, r, t)} = t(, r, t)cov {c(r, t)} tH(, r, t). (44)

Since C is a diagonal matrix and its entries are independent random variables, cov(c) is a diagonal matrix and has non-zero value with the nNth entry. We then have:

var { A(, r, t)} = ktk(, r, t)tHk(, r, t)cov{c(r, t)}kk

(45)

where k = (n 1)N + 1 denotes the kth entry of vector. So calculating the var { A(, r)} can be equivalently obtained by calculating cov{c(r)}kk. This again is a very general result that we can evaluate and review for different statistical models for cov{c(r)}kk.

cov{c(r, t)}kk =

+

+ E(C(r, t))nnE(C(r, t))nnf ()d

=

(47)

where [ ] denotes the conjugate operator. Here, the

E (C(r, t)nn) = e2

2(n1)2

+ C(r, t)nnC(r, t)nnf ()d E(C(r, t))2nn

t rc0 2

2

is uti

lized. Applying Equation 32, the first term of Equation 47 can be rewritten as:

= E

C(r, t)nn

+ C(r, t)nnC(r, t)nnf ()d =

1 2

+ ej2(n1)

t rc0 ej2(n1)

t rc0 e

222 d = 1

(48)

Applying it to Equation 47 yields:

cov{c(r, t)}kk = 1 E(C(r, t))2nn

= 1 e4

2(n1)2

t rc0 2

2 (49)

where k is the kth element of c(r, t), which is used to distinguish from n. Consequently, we can get the FDA beampattern deviation variance:

var { A(, r, t)} =

ktk(, r, t)tHk(, r, t)

1 e42(n1)2 t rc0 2 2

.(50)

Similarly, we can get the FDA beampattern variance with the random frequency increment errors, which are uniformly distributed as aforementioned:

var { A(, r, t)} =

ktk(, r, t)tHk(, r, t)

1

+ c(r, t)kk E(c(r, t)kk)

2f ()d

(46)

Assume that all the random frequency increment errors have the same distribution. For the first case, frequency increment error n of the nth element satisfies the Gaussian statistical model n N(0, 2). According to the

Equations 13 and 36, 46 can be rewritten as:

cov{c(r, t)}kk =

=

+ C(r, t)nn E(C(r, t)nn)

2f ()d

+ C(r, t)nn E(C(, r, t)nn)

sin(2(n 1)

t rc0 max)

2(n 1)

.

(51)

Compared with Equations 50 and 51, we can see that the FDA beampattern variances at both kinds of random frequency increment error distributions are dependent on the statistics properties. Given the same frequency increment errors, the variances will decrease when the number of elements is increased.

4.2.3 PDF

We have shown how to calculate the expectation value and the variance of the beampattern deviation and the beampattern itself. A remained question is: What form will their respective PDF have? In this subsection, we will

C(r, t)nn E(C(r, t)nn)

f ()d

t rc0 max

2

=

+ C(r, t)nnC(r, t)nn E(C(r, t)nn)C(r, t)nn

C(r, t)nnE(C(r, t)nn) +E(C(r, t)nn)E(C(r, t))nn

f ()d

+ C(r, t)nnC(r, t)nnf ()d

= + E(C(r, t))nnC((r, t)nn)f ()d

+ C(r, t)nnE(C(r, t))nnf ()d

Gao et al. EURASIP Journal on Advances in Signal Processing (2015) 2015:34 Page 8 of 12

derive the PDF of beampattern based on the PDF of the frequency increment errors.

For the first case, we investigate the PDF of beampattern when the frequency increment errors obey the Gaussian distribution as aforementioned. According to Equation 13, we derive the PDF of C. It is known that:

Cnn = ej2(n1)

t rc0

where vn denotes the nth element of v. Defining vn =

vn(,r,t)r for notation convenience, we can find

%1 g22 + jg2 = A(, r, t)

N

n=3

%

1 g2n + jgn

vnv1

/v2.

(57)

Let

a2 = imag

= cos

2(n 1)

t r c0

A

(, r)

N

&n=3

$

1 g2n + jgn

vn v1

(52)

j sin

2(n 1)

t r c0

= imag( A(, r)) imag

N

&n=3

$

1 g2n + jgn

vn v1

= real(Cnn) j

imag(Cnn)

where real(Cnn) = cos

2(n 1)

t r c0

and

imag(Cnn) = sin

2(n 1)

t r c0

denote the real part and imaginary part of Cnn, respectively. Consequently, we have:

=

arcsin(imag(Cnn))

2(n 1) t r c0

and v2 = cos +jsin, where is the phase of v2, utilizing Equation 57, we can get:

g2 cos +

%1 g22 sin = a2. (58)

Solving Equation 58, it yields:

g2 = a2 cos + sin $1 a22

= 2( Ai)

(59)

arcsin(gn)

where Ai = imag

A(, r)

denotes the imaginary part of

A. And we define:g3 =g3 = 3(g3)

. . .

gN =gN = N(gN).

(53)

where gn = imag(Cnn) denotes the imaginary part of Cnn. As we known that if a variable x has PDF fx(x), the PDF of y = hg(x) is as follows:fy(y) = fx[ h(y)]

h (y)

= 2(n 1)

t r c0

= h(gn)

(60)

The joint PDF of the variables above is as follows:

fAig3...gN ( Aig3 . . . gN)

= fg1g3...gN

(54)

where x = h(y) is the inverse function of hg(x), fy() denotes the PDF responding to y, and [ ] denotes the derivation operation. We have:

fgn(gn) = f[ h(gn)]

h (g

2 Ai(, r, t)

, 3(g3), . . . , N(gN)

|J|

(61)

n)

h (g

, 1 < gn < 1

(55)

According to Equations 59 and 60, we can get that:

J =

2 Ai

2

g3 . . .

2

gN

=

1 2 e

h2(gn) 22

n)

2 Ai

2

g3 . . .

2

gN

3 Ai

3

g3 . . .

3

gN

=

where

. . . . . . . . . . . .

N

Ai

N

g3 . . .

0 1 . . . 0 . . . . . . . . . . . .

0 0 . . . 1

1

N

gN

(62)

h (gn) =

2(n 1)

$1 g2n t r c0

= cos

sin

%1 a22

denotes the derivation of h(gn). Utilizing gn, Equation 13 can be rewritten as:

A(, r, t) =

wHC(, r, t)v(, r, t) r

=

where 1Ai = cos

sin

%1a22

N

n=2

%

1 g2n + jgn

vn(, r, t)r +

v1(, r, t) r

. Since the each element of C is

responding to only one frequency increment error which is random and independent, gn is the independent with other gm with m = 1, 2, . . . , N 1 and m = n. Therefore,

Equation 61 can be rewritten as:

fg2g3...gN

2( Ai(, r, t)), 3(g3), . . . , N(gN)

= fg2

=

N

n=2

%

1 g2n + jgn

vn + v1

(63)

(56)

2( Ai(, r, t))

fg3

g3 . . . fgN

gN

Gao et al. EURASIP Journal on Advances in Signal Processing (2015) 2015:34 Page 9 of 12

Moreover, we can get the PDF of FDA beampattern as:

fAi ( Ai) =

1

a

5

1

. . .

1

1

fg2g3...gN

2( Ai(, r, t)), 3(g3), . . . , N(gN)

0

Bideal

BRandom

Bbound

cos

sin

%1 a22

5

Beampattern (dB)

dg3 . . . dgN .

10

(64)

Utilizing the same approach, we can get the PDF of real part of FDA beampattern. Since the relationship between real part and imaginary part is complex, so we cannot give the PDF of the whole FDA beampattern.

Similarly, we can derive the FDA beampattern variance with the random frequency increment errors, which are uniformly distributed as aforementioned. The PDF of gn is

as follows:

fgn(gn) = f[ h(gn)]

h (g

n)

15

20

25

Angle (Degree)

b

5

Bideal

B

Random

Bbound

0

1

2max

1

(65)

and the PDF of FDA beampattern imaginary part is as follows:

fAi ( Ai) =

1

$1

=

,

2(n 1)

t r c0

g2n

1 < gn < 1

Beampattern (dB)

5

10

15

20

1

. . .

1

1

fg2g3...gN

2( Ai(, r, t)), 3(g3), . . . , N(gN)

25 0.8 0.9 1 1.1 1.2x 104

cos

Range (m)

sin

%1 a22

dg3 . . . dgN .

Figure 4 Profiles for the ideal, the bounds, and the errors FDA beampattern. (a) In angle dimension. (b) In range dimension.

(66)

5 Simulation results5.1 Example 1: FDA beampattern bound

Consider a 16-element uniform linear FDA with half of wavelength spacing between neighbor elements. The center frequency f0 is 10 GHz, and the increment frequency f is 30 KHz. The target is located at the range 10 km, angle 0, and the time is on 1/ f . Figure 4a,b shows the comparisons of ideal FDA beampattern Bideal,

FDA beampattern upper and lower bound Bbound, and

FDA beampattern with random frequency increment errors BRandom. In the FDA beampattern with random fre

quency increment errors, the frequency increment errors are Gaussian distribution, i.e., n N(0, 902) and the random FDA beampatterns are based on 50 independent Monte Carlo simulation runs. It can be shown that the present bounds hold for all the simulated beampattern realizations for both Figure 4a,b. Simulations results show that the bounds are tighter in range dimension than

50 0 50

that in the angle dimension which is caused by range counterbalance of the total error in beampattern.

Furthermore, we give two figures about expectation and variance about FDA beampattern errors caused by frequency increment errors, shown as in Figures 5 and 6, respectively. Figure 5 compares the expectation value of theoretical result and empirical result. The theoretical result is based on the method of Section 4.2.1. The empirical result is based on the 10,000 independent Monte Carlo simulation runs. The comparisons show that when the FDA beampattern errors have smaller value, two results are close to each other. The variances of theoretical result and empirical result are shown in Figure 6. Similar with Figure 5, the theoretical result is based on the method of Section 4.2.2. The empirical result is based on the 1,000 independent Monte Carlo simulation runs. Different from the expectation, the empirical result fluctuates around the theoretical result along with increase of .

Gao et al. EURASIP Journal on Advances in Signal Processing (2015) 2015:34 Page 10 of 12

50

Theoretical

Empirical

100

Beampattern error (dB)

150

200

250

101 100 101 102

Figure 5 The comparison of Gaussian random FDA beampattern error expectation value.

5.2 Example 2: FDA beampattern PDF

Consider a uniform linear FDA with four elements. In this example, FDA beampattern amplitude error does not divide into the factor r0, and its value is very large compared with that of Figures 5 and 6. Other array parameters are the same to that of example 1. Figure 7 shows PDFs of theoretical and empirical results, and Figure 8 shows the PDFs of imaginary part of random FDA beampattern with different s. It can be known that lower

enjoys the smaller beampattern errors, and PDF curves are centrosymmetric with the center about Ai = 0. One might use this figure and Equation 64 to specify tolerance requirements of frequency increment errors to fulfill a given beampttern requrement with certain probability. For instance, if probability of the beampttern errors at the domain (0.5 to 0.5) is not smaller than 0.95, the standard deviations of for every elements of FDA array will be required to be no larger than 60.

40

50

Beampattern error (dB)

60

70

80

Theoretical

Empirical

90

100

110

101 100 101 102

Figure 6 The comparison of Gaussian random FDA beampattern error variance value.

Gao et al. EURASIP Journal on Advances in Signal Processing (2015) 2015:34 Page 11 of 12

1.4

Theoretical Empirical

1.2

1

0.8

PDF

0.6

0.4

0.2

0

1 0.5 0 0.5 1

Ai

Figure 7 The PDF of Gaussian random FDA beampattern error.

6 Conclusions

In this paper, we have investigated the impacts of frequency increment errors on FDA beampattern based on deterministic errors and random errors. For uniform and linear deterministic frequency increment errors, the specific beampattern error formulations are provided, which gives guideline for device selection. For the stochastic frequency increment errors, we have derived a very tight

worst-case boundary of the FDA beampattern. Simulation results show that all the random beampatterns are held for the derived bounds, and the worst-case boundary is helpful to the FDA system design. At last, we derived the statistical properties of the expectation, variance, and PDF. They can be used to analyze the probability of FDA beam-pattern fluctuations for the given distribution frequency increment errors.

2

=60 =90 =150

1.5

PDF

1

0.5

0

1 0.5 0 0.5 1

Ai

Figure 8 The PDF of imaginary part of random FDA beampattern.

Gao et al. EURASIP Journal on Advances in Signal Processing (2015) 2015:34 Page 12 of 12

Competing interests

The authors declare that they have no competing interests.

Authors contributions

All authors formulated and discussed the idea together. Additionally, KG wrote the paper. All authors read and approved the final manuscript.

Acknowledgements

This work was supported in part by the Program for New Century Excellent Talents in University under grant no. NCET-12-0095.

Received: 3 November 2014 Accepted: 12 March 2015

References

1. BDV Veen, KM Buckley, Beamforming: a versatile approach to spatial filtering. IEEE ASSP Mag. 5(2), 424 (1988)

2. J Xie, H Li, Z He, C Li, A robust adaptive beamforming method based on the matrix reconstruction against a large doa mismatch. EURASIP J. Adv. Signal Process. 2014(91), 110 (2014)

3. P Antonik, MC Wicks, HD Griffiths, CJ Baker, in Proceedings of the IEEE Radar Conference. Frequency diverse array radars (IEEE, Verona, NY, 2006),pp. 215217

4. MC Wicks, P Antonik, Frequency diverse array with independent modulation of frequency, amplitude, and phase. in January 15, 2008. (U.S.A Patent 7,319,427, USA)

5. P Antonik, MC Wicks, Frequency diverse array with independent modulation of frequency, amplitude, and phase. in June 5, 2008. U.S.A Patent 7, Application 20080129584, USA

6. P Antonik, MC Wicks, Method and apparatus for a frequency diverse array. in March 31, 2009. U.S.A Patent 7.511, 665B2, USA

7. P Antonik, An investigation of a frequency diverse array. PhD thesis, University College London (2009)

8. S Brady, Frequency diverse array radar: signal characterization and measurement accuacy. PhD thesis, Air Force Institute of Technology (2010)

9. M Secmen, S Demir, A Hizal, T Eker, in Proceedings of the IEEE Radar Conference. Frequency diverse array antenna with periodic time modulated pattern in range and angle (IEEE, Boston, MA, 2007), pp. 427430

10. P Baizert, TB Hale, MA Temple, MC Wicks, Forward-looking radar GMTI benefits using a linear frequency diverse array. Electron. Lett. 42(22), 13111312 (2006)

11. C Cetintepe, S Demir, Multipath characteristics of frequency diverse arrays over a ground plane. IEEE Trans. Antennas Propagation. 62(7), 35673574 (2014)

12. J Shin, J-H Choi, J Kim, J Yang, W Lee, J So, C Cheon, in Microwave Conference Proceedings (APMC). Full-wave simulation of frequency diverse array antenna using the FDTD method (IEEE, Seoul, Korea, 2013),pp. 10701072

13. T Eker, S Demir, A Hizal, Exploitation of linear frequency modulated continuous waveform (LFMCW) for frequency diverse arrays. IEEE Trans. Antennas Propagation. 61(7), 35463553 (2013)

14. YB Wang, WQ Wang, HZ Shao, Frequency diverse array Cramr-Rao lower bounds for estimating direction, range and velocity. Int J Antennas Propagation. 2014, 110 (2014)

15. V Ravenni, in Proceedings of the 4th European Radar Conference. Performance evaluations of frequency diversity radar system (IEEE, Munich, Germany, 2009), pp. 436439

16. L Zhuang, XZ Liu, Application of frequency diversity to suppress grating lobes in coherent MIMO radar with separated subapertures. EURASIP J. Adv. Signal Process. 2009, 110 (2009)

17. WQ Wang, Phased-MIMO radar with frequency diversity for range-dependent beamforming. EURASIP J. Adv. Signal Process. 13(4), 13201328 (2013)

18. W-Q Wang, HZ Shao, Range-angle localization of targets by a double-pulse frequency diverse array radar. IEEE J. Sel. Top. Signal Process. 8(1), 106114 (2014)

19. W-Q Wang, HC So, Transmit subaperturing for range and angle estimation in frequency diverse array radar. IEEE Trans. Signal Process. 62(8), 20002011 (2014)

20. E Yazdian, S Gazor, MH Bastani, Limiting spectral distribution of the sample covariance matrix of the windowed array data. EURASIP J. Adv. Signal Process. 2013(42), 115 (2013)

21. J Xie, Z He, H Li, J Li, 2D DOA estimation with sparse uniform circular arrays in the presence of mutual coupling. EURASIP J. Adv. Signal Process. 2011(127), 118 (2011)

22. M Khodja, A Belouchrani, K Abed-Meraim, Performance analysis for time-frequency music algorithm in presence of both additive noise and array calibration errors. EURASIP J. Adv. Signal Process. 2012(94), 111 (2012)

23. Y Han, D Zhang, A recursive Bayesian beamforming for steering vector uncertainties. EURASIP J. Adv. Signal Process. 1(108), 110 (2014)

24. S Henault, SK Podilchak, SM Mikki, YMM Antar, A methodology for mutual coupling estimation and compensation in antennas. IEEE Trans. Antennas Propagation. 61(3), 11191132 (2013)

25. H Chen, B Liu, P Huang, J Liang, Y Gu, Mobility-assisted node localization based on TOA measurements without time synchronization in wireless sensor networks. MONET. 17(1), 9099 (2012)

26. H Chen, G Wang, Z Wang, H-C So, HV Poor, Non-line-of-sight node localization based on semi-definite programming in wireless sensor networks. IEEE Trans. Wireless Commun. 11(1), 108116 (2012)

27. H Chen, Q Shi, R Tan, HV Poor, K Sezaki, Mobile element assisted cooperative localization for wireless sensor networks with obstacles. IEEE Trans. Wireless Commun. 9(3), 956963 (2010)

28. W Zhang, Q Yin, H Chen, F Gao, N Ansari, Distributed angle estimation for localization in wireless sensor networks. IEEE Trans. Wireless Commun. 12(2), 527537 (2013)

29. CM Schmid, S Schuster, R Feger, A Stelzer, On the effects of calibration errors and mutual coupling on the beam pattern of an antenna array. IEEE Trans. Antennas Propagation. 61(8), 40634071 (2013)

30. H Chen, F Gao, MH Martins, JL P Huang, Accurate and efficient node localization for mobile sensor networks. MONET. 18(1), 141147 (2013)

Submit your manuscript to a journal and benet from:

7 Convenient online submission7 Rigorous peer review7 Immediate publication on acceptance7 Open access: articles freely available online 7 High visibility within the eld7 Retaining the copyright to your article

Submit your next manuscript at 7 springeropen.com