arxiv: v1 [cs.ce] 3 Apr 2016

Transcription

1 Two Dimensional Angle of Arrival Estimation arxiv: v1 [csce] 3 Apr Introduction Santhosh Kumar, Pradip Sircar Department of Electrical Engineering Indian Institute of Technology Kanpur Kanpur , Uttar Pradesh, India We present a new method for the estimation of two dimensional (2D) angles of arrival (AOAs), namely, azimuth and incidence angles of multiple narrowband signals of same frequency in the far field of antenna array In this report, we propose an algorithm which uses an L-shape array configuration [1] We extend the polynomial based approach presented in [2] to two dimensional AOA estimation We Consider an L-shaped antenna array as shown in Figure 1 in the XZ plane [1] Figure 1: The L-shaped array configuration The gain of an element at a position (x, y, z) for a signal with reference to 1

2 origin from the direction (θ, φ) is given by a(θ, φ) = e j(2π/λ)(x sin θ cos φ+y sin θ sin φ+z cos θ) (1) Suppose that there are q narrow band sources, with same wavelength λ impinging on each sensor An lth source has an incidence angle θ l and an azimuth angle φ l, l = 1, 2,, q The signal received at the ith element of the Z-subarray is given by [1] z i (t) = a z,i (θ l, φ l )s l (t) + n z,i (t) (2) a z,i (θ l, φ l ) = e jβd cos θ l = e j2π(i 1)(d/λ) cos θ l = e j(i 1)ψ l ψ l = 2π(d/λ) cos θ l (3) s l (t) is the complex envelope of the lth signal arriving from the direction (θ l, φ l ) and n z,i (t) is the zero-mean additive white Gaussian noise (AWGN) whose properties are mentioned in Section 22 The signal received at the ith element of the X-subarray is given by x i (t) = a x,i (θ l, φ l )s l (t) + n x,i (t) (4) a x,i (θ l, φ l ) = e j2π(i 1)(d/λ) sin θ l cos φ l = e j(i 1)ξ l, ξ l = 2π(d/λ) sin θ l cos φ l (5) and n x,i (t) is the zero-mean AWGN If we collect outputs of all the sensors on the Z-subarray into an m 1 vector, we have z(t) = a z (θ l, φ l )s l (t) + n z (t) (6) z(t) = [z 1 (t) z 2 (t) z m (t)] T a z (θ l, φ l ) is the m 1 steering vector of the Z-subarray towards the direction (θ l, φ l ) given by a z (θ l, φ l ) = [a z,1 (θ l, φ l a z,2 (θ l, φ l ) a z,m (θ l, φ l )] T 2

3 and We can rewrite (6) as n z (t) = [ n z,1 (t) n z,2 (t) n z,m (t) ] T a z,1(θ 1, φ 1) a z,1(θ 2, φ 2) a z,1(θ q, φ q) a z,2(θ 1, φ 1) a z,2(θ 2, φ 2) a z,2(θ q, φ q) z(t) = a z,m(θ 1, φ 1) a z,m(θ 2, φ 2) a z,m(θ q, φ q) Similarly for the X-subarray, we have s 1(t) s 2(t) s q(t) + n z,1(t) n z,2(t) n z,m(t) (7) x(t) = a x (θ l, φ l )s l (t) + n x (t) (8) x(t) = [x 1 (t) x 2 (t) x m (t)] T and a x (θ l, φ l ) = [a x,1 (θ l, φ l a x,2 (θ l, φ l ) a x,m (θ l, φ l )] T n x (t) = [n x,1 (t) n x,2 (t) n x,m (t)] T Note that (8) can be rewritten as a x,1(θ 1, φ 1) a x,1(θ 2, φ 2) a x,1(θ q, φ q) a x,2(θ 1, φ 1) a x,2(θ 2, φ 2) a x,2(θ q, φ q) x(t) = a x,m(θ 1, φ 1) a x,m(θ 2, φ 2) a x,m(θ q, φ q) s 1(t) s 2(t) s q(t) + n x,1(t) n x,2(t) n x,m(t) (9) 2 Proposed Method Consider the Z-subarray: On substituting the values of a z,i (θ l, φ l ) in (7) for i = 1, 2,, m and l = 1, 2,, q, we have z(t) = e jψ1 e jψ2 e jψq e j(m 1)ψ1 e j(m 1)ψ2 e j(m 1)ψq {ψ l } are defined in (3) s 1 (t) s 2 (t) s q (t) + n z,1 (t) n z,2 (t) n z,m (t) (10) 3

4 By collecting M equispaced snapshots from each sensor, a data matrix Z is formed as follows s 1 (t 1 ) s 1 (t 2 ) s 1 (t M ) e jψ1 e jψ2 e jψq s 2 (t 1 ) s 2 (t 2 ) s 2 (t M ) Z = e j(m 1)ψ1 e j(m 1)ψ2 e j(m 1)ψq s q (t 1 ) s q (t 2 ) s q (t M ) or Z = + n z,1 (t 1 ) n z,1 (t 2 ) n z,1 (t M ) n z,2 (t 1 ) n z,2 (t 2 ) n z,2 (t M ) n z,m (t 1 ) n z,m (t 2 ) n z,m (t M ) (11) Z = A z S + N z (12) z 1 (t 1 ) z 1 (t 2 ) z 1 (t M ) z 2 (t 1 ) z 2 (t 2 ) z 2 (t M ) z m (t 1 ) z m (t 2 ) z m (t M ) The matrices A z and S are unknown, and are not rank-deficient by assumption Now, let y 1 = e jψ1, y 2 = e jψ2,, y q = e jψq and assume that y 1, y 2,, y q are the roots of a polynomial equation P(y) = 1 + c 1 y + c 2 y c m 1 y m 1 = 0 (13) Now our aim is to calculate the coefficients which satisfy this assumption and then to solve for the roots 21 Method without Noise From (13) and (11), it can easily be shown that the following equations hold for no noise case z 1 (t 1 ) + c 1 z 2 (t 1 ) + c 2 z 3 (t 1 ) + + c m 1 z m (t 1 ) = 0 z 1 (t 2 ) + c 1 z 2 (t 2 ) + c 2 z 3 (t 2 ) + + c m 1 z m (t 2 ) = 0 z 1 (t M ) + c 1 z 2 (t M ) + c 2 z 3 (t M ) + + c m 1 z m (t M ) = 0 4

5 After simple manipulations, c 1 z 2 (t 1 ) + c 2 z 3 (t 1 ) + + c m 1 z m (t 1 ) = z 1 (t 1 ) c 1 z 2 (t 2 ) + c 2 z 3 (t 2 ) + + c m 1 z m (t 2 ) = z 1 (t 2 ) c 1 z 2 (t M ) + c 2 z 3 (t M ) + + c m 1 z m (t M ) = z 1 (t M ) The above set of equations can be written in matrix form as z 2 (t 1 ) z 3 (t 1 ) z m (t 1 ) c 1 z 2 (t 2 ) z 3 (t 2 ) z m (t 2 ) or z 2 (t M ) z 3 (t M ) z m (t M ) c 2 c m 1 = z 1 (t 1 ) z 1 (t 2 ) z 1 (t M ) (14) P C = P 1 (15) This is an overdetermined system of equations The best solution in the least square sense for the coefficient vector C is given by P # is the Moore-Penrose pseudoinverse of P, C = P # P 1 (16) P # = (P H P ) 1 P H Substitute c 1, c 2,, c m 1 obtained in (16) into (13) and solve for the roots of P(y) There are (m 1) roots, and among those only q solutions are of our interest From the (m 1) roots choose q roots say {γ l }, whose magnitudes are nearer to unity: Then, the estimate of ψ l is γ l = e j2π(d/λ) cos θ l l = 1, 2,, q (17) ˆψ l = angle{γ l } (18) Follow the similar approach for the X-subarray to calculate {β l } defined by β l = e j2π(d/λ) sin θ l cos φ l l = 1, 2,, q (19) and hence ˆξ l = angle{β l } (20) Now the estimates of θ l and φ l are given by ( ) ˆψ ˆθ l = sin 1 l 2π(d/λ) ( ) ˆξ ˆφ l = cos 1 l (2π(d/λ) sin ˆθ l ) l = 1, 2,, q (21) l = 1, 2,, q (22) 5

6 22 Method with Noise Note that in (6) and (8), the vectors n z (t) and n x (t) are random noise components on Z and X subarrays respectively The noise is assumed to be additive, temporally white Gaussian with zero mean This assumption is a valid assumption in real-world applications Mathematically, we can express the above stated properties of the noise as follows, E[n z (t)] = E[n x (t)] = 0 (23) and E[n z (t)n H z (t)] = E[n x (t)n H x (t)] = σ 2 I (24) E[n z (t)n T z (t)] = E[n x (t)n T x (t)] = 0 (25) E[n z (t i )n H z (t j )] = E[n z (t i )n T z (t j )] = 0 for i j (26) E[n x (t i )n H x (t j )] = E[n x (t i )n T x (t j )] = 0 for i j (27) E denotes the expectation operation, H denotes the Hermitian operation, T denotes the transposition, and t i and t j are the ith and jth time instants respectively The noise has an negative impact on the accuracy of estimates of the polynomial coefficients In order to reduce the effect of the noise, we follow the same procedure except for the calculation of inverse of the matrix P in (15) we use the singular value decomposition(svd) as discussed below Let the SVD of P be P = UΣV H (28) U and V are the unitary matrices containing the singular vectors of P and Σ is a diagonal matrix of same dimensions as P, and it is given by Σ = U H P V (29) The diagonal elements in Σ are non-negative and are sorted in decreasing order, ie, Σ = diag{τ 1, τ 2,, τ q,, τ m 1 } (30) τ 1 τ 2 τ q τ m 1 are the singular values of P To reduce the effect of noise, replace the last m q 1 singular values with zero in calculating the inverse of P, ie, Σ 1 = diag{τ 1, τ 2,, τ q, 0, } (31) and Σ 1 1 = diag{1/τ 1, 1/τ 2,, 1/τ q, 0, } (32) 6

7 The dimensions of Σ 1 1 are same as that of P H The lease squares approximation of P 1 is given by [3] hence, the estimate of coefficient vector in (15) is P 1 = V Σ 1 1 U H (33) C = V Σ 1 1 U H P 1 (34) Follow the same approach for X-subarray to obtain the coefficient vector and the roots The computation of AOAs from the roots is same as that for the noiseless case References [1] N Tayem and H M Kwon, L-shape 2-dimensional arrival angle estimation with propagator method, IEEE Trans Antennas and Propagation, vol 53, no 5, pp , May 2005 [2] P Singh and P Sircar, Time delays and angles of arrival estimation using known signals, Signal, Image and Video Processing, vol 6, no 2, pp , Feb 2011 [3] P Dewilde and E F Deprettere, Singular value decomposition: An introduction, In E F Deprettere (ed), SVD and Signal Processing: Algorithms, Applications, and Architectures, pp 3 41 Elsevier Science, North Holland,