《Performance Analysis of AngularSmoothing.doc》由会员分享，可在线阅读，更多相关《Performance Analysis of AngularSmoothing.doc（27页珍藏版）》请在三一办公上搜索。
1、豆丁网精品论文Performance Analysis of Angular-SmoothingBased Root-MUSIC for an L-ShapedAcoustic Vector-Sensor ArrayYougen Xu Zhiwen LiuDepartment of Electronic Engineering, Beijing Institute of Technology, Beijing, PRC 100081yougenxu zwliuAbstractEigenstructure-based direction-of-arrival (DOA) estimation a
2、lgorithms such as Multiple Signal Classification (MUSIC), Root-MUSIC, Estimation of Signal Parameters via Rotational Invariance Techniques (ESPRIT), encounter great difficulty in the presence of perfectly correlated incident signals. For an array composed of a number of translational invariant subar
3、rays such as a uniform linear scalar-sensor array, this problem can be solved by spatial smoothing. An array of identically oriented acoustic vector-sensors can be grouped into four coupled subarrays of identical grid geometry, respectively corresponding to the pressure sensors and differently orien
4、ted velocity-sensors. These four subarrays are angular invariant dependent only on signals direction cosines and an angular smoothing can be exploited for source decorrelation. In this paper, the performance of root-MUSIC incorporated with angular smoothing for correlated source direction finding wi
5、th an L-shaped acoustic vector-sensor array is analyzed in terms of the overall root mean-square errors (RMSE) of DOA estimates. We derive the analytical expression of the RMSE and compare it with simulation results and that of spatial smoothing for a rectangular pressure-sensor array instead.Keywor
6、ds: Antenna arrays, array signal processing, direction-of-arrival estimation, acoustic array1 Introduction1.1 Acoustic Vector-SensorsAcoustic vector-sensors have attracted increasing interest for subspace-based direction-finding8-17. A “complete” acoustic vector-sensor such as a vector-hydrophone an
7、d a microflown consists of a pressure sensor and a collocated triad of three orthogonal velocity sensors 8. These four sensors together measure the scalar acoustic pressure and all three components of the acoustic particle velocity vector of the incident wave-field at a given point. An “incomplete”
8、acoustic vector-sensor consists of a subset of the above four component-sensors, e.g., a three-component hydrophone formed from two orthogonally oriented velocityThis work was supported by the National Natural Science Foundation of China under grant no.60272025, the University Basic Research Foundat
9、ion of Beijing Institute of Technology under grant no.BIT-UBF-200301F07 and no. BIT-UBF-200501F4202, and the Specialized Research Fund for the DoctoralProgram of Higher Education under grant no. 20040007013.hydrophones and a pressure hydrophone 13.It is assumed that the acoustic wave is traveling in
10、 a quiescent, homogeneous, and isotropic fluid, and isfrom a source of azimuth and elevation , where 0 Q ). This proposed angular smoothing scheme represents a counterpart of polarisation smoothing 18 for an array of identically orientedelectromagnetic vector-sensors.Although developed here for a ve
11、ry regular L-shaped acoustic vector-sensor array to apply computationally efficient root-MUSIC, the angular smoothing algorithm is applicable for an array of identically oriented acoustic vector-sensors but with any arbitrary grid geometry. In addition, unlike classical spatial smoothing,the angular
12、 smoothing can be accomplished without spatial-aperture loss. Since it is impossible thatux = 1 ,uy = 1anduz = 1may hold simultaneously, this angular smoothing algorithm can still avoidsignal cancellation in the presence of totally uncorrelated signals. These attractive features are analogous to tho
13、se of the polarisation smoothing algorithm in the context of electromagnetic applications 18.It is important to note that the angular-smoothing scheme generally can only process at most three coherentsignals with distinct DOAs even with a very large number of N xandN y . This shortcoming may beallev
14、iated by incorporating the conventional spatial smoothing. However, such a consideration is beyond the scope of this paper.1.3 The Angular Smoothing-Based Root-MUSICAfter angular smoothing, (8) is equivalent to the covariance matrix with desired rank property of an L-shaped array of conventional sca
15、lar sensors (e.g., the pressure hydrophones), and hence two parallel root-MUSIC can be applied to the two legs of such arrays to obtain the corresponding two directioncosines.x T TthDefineJ = jN ,1, J withJ = O1, , I1 ,jthen column ofI , andNx NyNx M ,n My y xyJ = IN , ON ,N 1 . Then we constructR x
16、 = Jx R (Jx )H ,R y = Jy R (Jy )H, and further eigen-decomposethem to obtain R = E (E )+ E (E )and R = E (E )+ E (E ), whereE ,s ssn n nx x xx H x x x Hy y y y H y y y H xEnns ss n n nscalled the x-leg signal matrix and span a Q -dimensional signal subspace, comprises the Q eigenvectorscorresponding
17、 to the Q largest eigenvalues ofR x , whereas xis composed of the remaining N Qsneigenvectors associated with the smallest eigenvalues and called the x-leg noise matrix, xand xaretwo diagonal matrices whose diagonal entries are respectively the Q largest eigenvalues and the N Qsnssmallest eigenvalue
18、s. The y-leg counterpartsE y ,E y , yand yhave the same definitions. Letn nnn nnP x = E x (E x )H ,P y = E y (E y )H, and bx (, ) = 1,e j x , .,e j (Nx 1)x T(9)by (, ) = 1,e j y , .,e j (Ny 1)y Tx H x xy H y yThen,b (q , q )Pn b (q , q ) = 0 ,b (q , q )Pn b (q , q ) = 0 , and hence the two root-MUSI
19、Cpolynomials are as follows: fx x T x x() = b (1 / ) Pn b ()qy y(10)y y T y y f() = b (1 / ) Pn b ()qx xwhose Q roots of unit modulus are respectivelyv x = expj 2d u/ andvy = expj 2d u/ ,whereq = 1, 2 , .,Q. Then we haveu = cos sin = (2d )1 (v x )and1yx ,qqqx quy,q = sin q sin q = (2dy )Q(vq )( () Q
20、denotes angle operation), and (q , q )can ultimately beobtained via couplingvx ,q q = 1andvy,q q = 1 .zuvector-sensor dyydx xFigure 1. The concerned sensor-source geometry2 Overall Root Mean-Square Errors of Direction Cosine EstimatesPractically, we have no access to the true RyKand hence true R due
21、 to the existence of noise and finitesample support. Commonly, the true Rycan be replaced by the following sample covariance matrix: Ry =1 y(t )yH (t )(11)k kK k = 1where K is the number of independent snapshots, and perturbation inR = 4 J R JTresults inthe error of ultimate DOA estimates.m = 1 m m
22、y mIn this section, we derive the overall root mean-square errors (RMSE) of the direction cosine estimateswhich is defined asQx ,qy,q = E (| u|2 ) + E (| u|2 )(12)q = 1by using root-MUSIC described above.The derivations below are primarily along the lines of 21, and will be focused on the x-leg only
23、. The discussion for the y-leg is similar. First, rewrite the null spectrum along the x-axis as x H x xf (x ) = b (x )Pn b (x ) . By differentiating the perturbedf (x ) , an estimate off (x ) , with respect tox , and using a first-order approximation, it can be derived that 21E(| u|2 ) = 2 1 x ,q 2d 2K (bx ()H Px bx ()2 x x ,q n x ,q 2x ,q 4 (13) m n (x ,q Rm,n x ,q )(x ,q Rn,m x ,q ) HH4m,n =1H H x x+ Rem n (x ,q Rm,n x ,q )(x ,q Rn,m x ,q )wherexx ,q = 2dx ux ,q / m,n =1,b (x q ) = b (x ) / x | = ,andx x H T2,x H x Tx x ,qRm,n = (Rn,m )= Jm Rx ,0 Jn + n (m n)IN,
