Correlation-based Algorithm for Multi-Dimensional Single-Tone Frequency Estimation
|
|
- MargaretMargaret Bishop
- 5 years ago
- Views:
Transcription
1 Correlation-based Algorithm for Multi-Dimensional Single-Tone Frequency Estimation Weize Sun, H.C. So and Lanxin Lin Department of Electronic Engineering, City University of Hong Kong Abstract In this paper, parameter estimation for a R-dimensional (R-D) single cisoid with R 2 in additive white Gaussian noise is addressed. By exploiting the correlation of the data samples, we construct R single-tone sequences which contain the R-D frequency parameters. Based on linear prediction and weighted linear squares techniques, two proposals are developed for fast and accurate frequency estimation from each constructed sequence. The two devised estimators are proved to be asymptotically unbiased while their variances achieve Cramér-Rao lower bound when the signal-to-noise ratio and/or data length tend to infinity. Computer simulations are also included to compare the proposed approach with conventional R-D harmonic retrieval schemes in terms of mean square error performance and computational complexity. Keywords: frequency estimation, multi-dimensional spectral analysis, fast algorithm, correlation 1. Introduction Estimating the parameters of sinusoidal components from a finite number of noisy discrete-time measurements has been an important research topic be- Preprint submitted to Signal Processing September 25, 2012
2 cause it corresponds to many science and engineering problems [1] [3]. Apart from the standard one-dimensional signal model [4] [6], multi-dimensional spectral estimation [7] in fact has many applications such as array processing [8] [9], nuclear magnetic resonance (NMR) spectroscopy [10], wireless communication channel estimation [11] [12] as well as detection and localization of multiple targets using multiple-input multiple-output (MIMO) radar [13]. The crucial step in R-dimensional (R-D) sinusoidal parameter estimation, withr 2, istofindthefrequency parameterswhicharenonlinearintheobserved data. For 2-D frequency estimation, maximum-likelihood (ML) [14] [15] and subspace [8] [13], [16] [17] approaches are commonly used choices. For higher dimensional signals, the latter seems to be the only practical solution although its computational requirement is very high. State-of-the-art subspace methods include MUSIC [9], decoupled root-music [17], unitary ESPRIT (UE) [8], rank reduction estimator (RARE) [12], multi-dimensional folding (MDF) [11] and improved MDF (IMDF) [16], and they usually need eigenvalue decomposition (EVD) or singular value decomposition (SVD) in the algorithm implementation. By effectively utilizing the correlation of the observed data, we contribute to devising an accurate frequency estimator for a R-D single-tone with less computational load in this work. The rest of the paper is organized as follows. To facilitate the understanding of our underlying idea, we first develop the correlation-based approach for a 2-D cisoid in the presence of white Gaussian noise in Section 2. By exploiting correlation on the 2-D raw data, a novel single-tone sequence whose frequency corresponds to one of the two dimensions is constructed. The generalized weighted linear predictor (GWLP) [18], which can provide optimum 2
3 performance for 1-D single-tone in white noise, is then applied for estimation of the two frequency parameters in a separable manner. Note that linear prediction (LP) and weighted least squares (WLS) are two key ingredients in the GWLP. Two estimators are devised, and their mean and variance expressions are produced. It is worthy to point out that the derivation of the weighting matrix and performance analysis are different from those of [18]. Section 3 generalizes the proposed solutions to higher dimensional signals with the use of tensor algebra [19]. Simulation results are included in Section 4 to corroborate the theoretical development and to compare the correlation-based approach with the approximate iterative quadratic ML (AIQML) [15], IMDF [16] and UE [8] algorithms as well as Cramér-Rao lower bound (CRLB) [20]. It is demonstrated that the proposed methods are more computationally efficient and their performance is close to CRLB particularly when the signalto-noise ratio (SNR) is sufficiently high. Finally, conclusions are drawn in Section Two-dimensional Frequency Estimation We first define the notation as follows. Scalars, vectors, matrices and tensors are denoted by italic, bold lower-case, bold upper-case and bold calligraphic symbols, respectively. The magnitude, angle, real and imaginary parts of a complex scalar a are a, (a), R{a} and I{a}, while the variable, noise-free value and estimate of a vector a are represented by ã, ā and â, respectively. The mean and variance of â are denoted by E{â} and var(â). The I i is the i i identity matrix and 0 i j is the i j zero matrix. The operators for transpose, complex conjugate, conjugate transpose, inverse and 3
4 Khatri-Rao product are T,, H, 1, and, respectively. The mth element of a M 1 and (m,n) entry of A M N are denoted by [a] m and [A] m,n, respectively, whileweusea m1,m 2,...,m R torepresent the(m 1,m 2,...,m R )entry of a R-D tensor A M 1 M 2 M R where R 3. where In matrix form, the 2-D signal model with size M 1 M 2 is expressed as: Y = X+Q (1) [X] m1,m 2 = γe j(m 1ω 1 +m 2 ω 2 ), m r = 1,2,,M r, r = 1,2 (2) is the noise-free tone sample, and [Q] m1,m 2 is a zero-mean complex white Gaussian process with unknown variance σq 2. Moreover, γ is the complex amplitude while ω 1 ( π,π) and ω 2 ( π,π) are the frequencies, and they are all unknown constants. Here, the task is to determine ω 1 and ω 2 from the observed matrix Y with M = M 1 M 2 samples. The main idea of the algorithm development is to exploit correlation to convert Y into a vector which contains the information of ω 1 or ω 2. That is, ω 1 and ω 2 are estimated in a separable manner. We notice that X is of rank 1 and can be factorized as X = γg 1 g2 T (3) [ ] T where g r = e jωr e j2ωr e jmrωr, r = 1,2 (4) ] Defining X = [x 1 x 2 x M2 and writing Y and Q in the same way, we easily obtain the LP property: x 1 = e jω 2 γg 1 (5) x m = e jω 2 x m 1 (6) 4
5 As a result, pre-multiplying X by x H 1 gives γ 2 M 1 e jω 2 g2 T, which is a singletone vector with frequency ω 2. Since X is not available, we use its noisy version which is given by the zero-lag correlation between all columns of Y and its first column vector to achieve frequency estimation: z T (2) = yh 1 Y (7) with elements: [z (2) ] m2 = y H 1 y m2 = y H 1 (x m2 +q m2 ) (8) Note that the subscript of (2) in z refers to ω 2. According to the LP property, we obtain y H 1 x m = e jω 2 y H 1 x m 1. Hence (8) satisfies: [z (2) ] m2 e jω 2 [z (2) ] m2 1 = y H 1 ( qm2 e jω 2 ) q m2 1, m2 = 2,3,,M 2 (9) Denotingz (2)1 = [ T ] T [z (2) ] 1 [z (2) ] 2 [z (2) ] M2 1] andz(2)2 = [[z (2) ] 2 [z (2) ] 3 [z (2) ] M2 and assigning ρ 2 = e jω 2, we construct the LP error vector based on (9): e (2) = z (2)2 ρ 2 z (2)1 (10) The ω 2 is then estimated by minimizing the following WLS cost function: J( ρ 2 ) = ( ) H ( ) z (2)2 ρ 2 z (2)1 W M2 1 z(2)2 ρ 2 z (2)1 (11) where W M2 1 C (M 2 1) (M 2 1) is a weighting matrix. Differentiating J( ρ 2 ) with respect to ρ 2 and then setting the resultant expression to zero, we get the estimate of ρ 2 : ˆρ 2 = zh (2)1 W M 2 1z (2)2 z H (2)1 W M 2 1z (2)1 (12) 5
6 In Appendix A, we have derived the optimal weighting matrix, denoted by W M2 1(ω 2 ), which is a function of the unknown frequency ω 2 : W M2 1(ω 2 ) = [ E{e (2) e H (2)} ] 1 = P 1 (ω 2 ) (13) where 2+(M/M 2 1)/(SNR+1) e jω e jω 2 2 e jω P(ω 2 ) = e jω 2 2 e jω e jω 2 2 (14) with SNR = γ 2 /σq 2. In the following, two iterative algorithms are proposed to approximate the conceptual solution of (12). As SNR is unknown, our first proposal is to set [P(ω 2 )] 1,1 = 2 which is valid particularly when SNR is sufficiently large, namely, SNR >> M/M 2. Underthisapproximation, P(ω 2 )hasaclosed-form inverse and the phase of z H (2)1 W M 2 1(ω 2 )z (2)1 in (12) is zero [18]. Employing the substitution of ˆρ 2 = e jˆω 2, (12) is simplified to ˆω 2 = (z H (2)1 W M 2 1(ω 2 )z (2)2 ) (15) where the elements of W M2 1(ω 2 ) are: [W M2 1(ω 2 )] m,n = M 2min(m,n) mn M 2 e j(m n)ω 2 (16) m = 1,2,,M 2 1, n = 1,2,,M 2 1 6
7 In this study, we apply the GWLP [18] algorithm which updates ˆω 2 and W M2 1(ω 2 ) in an iterative manner. We refer this estimator to as correlation- 1 (C-1) method. The steps of finding ˆω 2 are summarized in Table 1. The estimation of ˆω 1 is performed using the same procedure on Y T. (i) Compute the z (2) using (7) (ii) Obtain an initial value of ˆω 2 using (15) with [W M2 1(ω 2 )] m,n = 0 for m n in (16), which is in fact the weighted linear predictor estimate [21] (iii) Construct W M2 1(ω 2 ) according to (16) with ω 2 = ˆω 2 (iv) Compute an updated ˆω 2 using (15) (v) Repeat Steps (iii) (iv) until a stopping criterion is reached. In this study, we stop for τ iterations. Table 1: Estimation algorithm for ω 2 In Appendix B, we have shown that when the frequency estimation error is sufficiently small, the mean and variance of ˆω r, r = 1,2, are: E{ˆω r } ω r (17) var(ˆω r ) 6 ( γ 2 +σ 2 Q) σ 2 Q ( M/Mr γ 2 +σ 2 Q) Mr (M 2 r 1) γ 2 (18) Equation (17) means that the algorithm is an asymptotically unbiased estimator. The CRLB for ω r, denoted by CRLB(ω r ), is [20]: CRLB(ω r ) = 6σ 2 Q M(M 2 r 1) γ 2 (19) 7
8 To study the limit of the C-1 estimator, we define the efficiency, which is the ratio of CRLB to the variance, denoted by: eff 1 (ω r ) = CRLB(ω ( r) γ 2 +M r /Mσ 2 ( ) Q) = (SNR+M r/m) var(ˆω r ) γ 2 +σq 2 (SNR+1) (20) whosevalueisbetween 0and1. Itiseasy toseethatthec-1methodachieves optimum performance when SNR, with eff 1 (ω r ) = 1. On the other hand, the above development may not be appropriate when M/M r, r = 1,2, is large particularly when SNR >> M/M r does not hold. To tackle this scenario, our second approximation on (12) is to discard [z (2) ] 1 of (7) in the computation. Let z (r)s be z (r) without the first element and follow (9) (16), we obtain ˆω r = (z H (r)s1 W M r 2(ω r )z (r)s2 ) (21) [ T where z (r)s1 = [z (r)s ] 1 [z (r)s ] 2 [z (r)s ] Mr 2] (22) [ T z (r)s2 = [z (r)s ] 2 [z (r)s ] 3 [z (r)s ] Mr 1] (23) We refer this estimator to as correlation-2 (C-2) method. As a result, the estimate of (21) is obtained in a similar manner as in Table 1. Following Appendix A and replacing M r by M r 1, we see that the frequency estimate is also asymptotically unbiased and its variance is: var(ˆω r ) 6 ( γ 2 +σ 2 Q) σ 2 Q ( M/Mr γ 2 +σ 2 Q) (Mr 1)[(M r 1) 2 1] γ 2 (24) Analogous to (20), the efficiency of this estimator is eff 2 (ω r ) = (M r 1)[(M r 1) 2 1](SNR+M r /M) M r (M 2 r 1)(SNR+1) (0,1) (25) and we see that eff 2 (ω r ) = 1 when both SNR and M r tend to infinity. 8
9 3. Extension to Higher Dimension In this section, we generalize our development to higher dimensional signals with R 3, where tensor representation is employed. The observed R-D sinusoidal signal is modeled as: Y = X +Q (26) where x m1,m 2,...,m R = γ R e jωrmr, m r = 1,2,...,M r, r = 1,2,...,R, (27) r=1 The Y M 1 M 2 M R is the tensorial structured data set with length M r along the rth dimension. The tensor X is the signal component where γ and ω r ( π,π) represent the unknown complex amplitude and frequency in the rth dimension. On the other hand, the entries in Q are zero-mean complex white Gaussian noises with unknown variances σq 2. Now the task is to find {ω r } from the M = R r=1 M r samples of Y. [ ] T, By writing g r = e jωr e jωr2 e jωrmr to align with the presentation in Section 2 we define the rth unfolding of X as the transpose version of [19]: [X] (r) = γ(g r+1 g R g 1 g 2 g r 1 )gr T ] = [x (r)1 x (r)2 x (r)mr (28) and multiplying [X] (r) by x H (r)1 : c (r) = x H (r)1 [X] (r) = M M r γ 2 e jωr g T r (29) 9
10 ] Analogousto(7) (9),wedefine[Y] (r) = [y (r)1 y (r)2 y (r)mr and[q] (r) = ] [q (r)1 q (r)2 q (r)mr and construct where z T (r) = y H (r)1[y] (r) (30) [z (r) ] mr = e jω 2 [z (r) ] mr 1 +y H (r)1( q(r)mr e jω 2 q (r)mr 1), mr = 2,3,,M r Generalizing (11), the WLS cost function is: J( ρ r ) = ( ) H ( ) z (r)2 ρ r z (r)1 W Mr 1 z(r)2 ρ r z (r)1 (31) (32) [ T where ρ r = e j ωr while z (r)1 = [z (r) ] 1 [z (r) ] 2 [z (r) ] Mr 1] and z(r)2 = [ ] T. [z (r) ] 2 [z (r) ] 3 [z (r) ] Mr Following (10) (15) and Appendix A with the assumption of SNR >> M/M r, the C-1 estimate of ω r, r = 1,2,,R, is: ˆω r = (z H (r)1w Mr 1(ω r )z (r)2 ) (33) Similarly, by discarding the first element in z (r), r = 1,2,,R, the C-2 frequency estimate has the same form as in (14). In Appendix B, we have proved that the higher dimensional C-1 and C-2 estimates are asymptotically unbiased and their variances equal (18) and (24), respectively. In summary, according to the development, the C-1 estimator is more preferable when SNR >> M/M r while the C-2 method is a better choice if M r >> 1, r = 1,2,,R. Finally, the complexity orders of the proposed methods are investigated as follows. In our study, the iterative procedure stops when the number of 10
11 iterations reaches the value of τ. According to Table 1, the major computations in the rth dimension has a complexity of O(τM r ), implying that the orders of complexity for the C-1 and C-2 algorithms are O( R r=1 τm r) and O( R r=1 τ(m r 1)), respectively. As a comparison, the complexity order of the AIQML [15] method is O( R r=1 τm r) while the IMDF [16] and UE [8] schemes have complexities of O(4k t (L E L E /L 1 )K E ) and O(2k t L E K E ), respectively, wherek t isaconstantdependsonthedesignofthesvdalgorithm, L E r +KE r 1 = M r, L E r > 1, KE r > 1, L E = R r=1 LE r and KE = R r=1 KE r. It can be roughly seen that when τ is small, the correlation-based methods are more computationally efficient. 4. Numerical Examples Computer simulations have been conducted to evaluate the R-D singetone frequency estimation performance of the correlation-based approach in terms of mean square frequency error (MSFE) and complexity performance. We also include comparisons with the AIQML [15], IMDF [16] and UE [8] algorithms as well as CRLB [20]. Three iterations, that is, τ = 3, are employed as the stopping criterion in the proposed and AIQML methods because no significant improvement is observed for more iterations. The tone magnitude is γ = 1 and its phase varies from π to π in each independent trial. We properly scale the zero-mean white Gaussian noise data to produce different SNR conditions, where SNR = γ 2 /σq 2 = 1/σ2, and all results provided are averages of 1000 independent. In the first test, we solely compare the performance of the C-1 and C-2 methods. Figures 1 and 2 show the MSFEs versus SNR for the 2-D and 3-D 11
12 ] [ ] ] [ ] cases with [ω 1 ω 2 = π and [ω 1 ω 2 ω 3 = π, respectively. The derived variances and CRLB are also included to verify the theoretical findings and investigate the algorithm optimality. Three choices of N, namely, N = 5, N = 10 and N = 20, are studied where M r = N for all dimensions. We see from both figures that the performance of the C-1 and C-2 estimators agrees with the analytic expressions of (18) and (24), respectively, when SNR 4dB. Although the C-1 method is generally superior to the C-2 algorithm, it is observed that the latter can give smaller MSFE at SNR = 10dB when N = 20 at R = 3. Furthermore, comparing to the situation when the data length is smaller, namely, N=5, the performance of the C-2 method is close to the CRLB for a larger N, say N=20, which also align with the development in Section 2. Inthesecond test, comparison ismadewiththeaiqml andimdfmethods in 2-D estimation with M 1 = M 2 = 10. The frequency parameters are ] [ ] [ω 1 ω 2 = π and the results are shown in Figure 3. We observe that the IMDF method has the best threshold performance at SNR = 6dB, while the AIQML and proposed methods have threshold SNRs at 2dB and 0dB, respectively. Moreover, the MSFEs of the AIQML and C-1 schemes attain the CRLB at SNR 10dB. On the other hand, the average computational times of the AIQML, IMDF, C-1 and C-2 algorithms for a single run are s, s, s and s, respectively, indicating the attractiveness of the proposed approach in terms of complexity. It is seen that the computational times also agree with the complexity order analysis in Section 3. In the third test, a comparative study in 3-D estimation is made with 12
13 the UE and IMDF algorithms. The results for M 1 = M 2 = M 3 = 10 and ] [ ] [ω 1 ω 2 ω 3 = π are plotted in Figure 4. At SNR 2dB, the UE and IMDF methods outperform the proposed estimators. However, when SNR > 5dB, the MSFE of the C-1 scheme achieves the CRLB. The average computational times for the UE, IMDF, C-1 and C-2 algorithms for a single run are s, s, s and s, which also agree with the complexity study in Section 3. Finally, we perform comparison for different R [3,5] with M r = 4, r = 1,2,,R, and the results are shown in Figure 5. The SNR is set to 20dB and all the frequencies varies from 0.99π to 0.99π in each independent run. It is seen that the C-1 algorithm performs the best and the theoretical calculations of (18) and (24) are again validated. Their average single run computational times are tabulated in Table 2, which agree with Section 3 and indicate the computational advantage of the proposed methodology particularly when the data dimension is large. R UE IMDF C-1 C Table 2: Average computational times versus R 5. Conclusion A fast and accurate multi-dimensional single frequency estimation approach based on correlation has been devised. Two estimators, C-1 and C-2, 13
14 which are ideal for high SNR and/or large sample conditions, are proposed. The main idea in the algorithm development is to construct single-tone sequences which contain the frequency parameters by exploiting correlation of the observed data. Frequency estimation for each dimension is performed separately by using LP and WLS techniques. It is proved that the variances of the C-1 and C-2 estimators attain CRLB in the presence of white Gaussian noise, when the SNR and/or data length tend to infinity. Furthermore, the proposed correlation approach is more computationally attractive than the AIQML, UE and IMDF schemes, and provides superior estimation accuracy at higher SNRs. As a future work, we will extend the correlation-based methodology to estimation of multiple tones. Appendix A We now derive the weighting matrix for the general signal model of (26). According to (31), the LP residue error vector for the rth dimension is e (r) = z (r)1 ρ r z (r)2 (A.1) where[e (r) ] nr = y(r)1 H ξ (r)n r,n r = 1,2,,M r 1,ξ (r)nr = ( ) q (r)(nr+1) ρ r q (r)nr andρ r = e jωr. Applying thegauss-markov theorem [22], the optimal weighting matrix which is a function of ρ r, denoted by W Mr 1(ρ r ), is: W Mr 1(ρ r ) = [ E { e (r) e(r)}] H 1 = B 1 (ρ r ) (A.2) where [B(ρ r )] m,n =E{[e (r) ] m [e (r) ] n} = E{y H (r)1 ξ (r)m ξ H (r)n y (r)1 } =[K 1 ] m,n +[K 2 ] m,n +[K 3 ] m,n +[K 4 ] m,n, m,n = 1,2,,M r 1 (A.3) 14
15 with [K 1 ] m,n =E{x H (r)1 ξ (r)mξ H (r)n x (r)1} (A.4) [K 2 ] m,n =E{q H (r)1 ξ (r)mξ H (r)n q (r)1} (A.5) [K 3 ] m,n =E{x H (r)1 ξ (r)mξ H (r)n q (r)1} (A.6) [K 4 ] m,n =E{q H (r)1 ξ (r)mξ H (r)n x (r)1} (A.7) On the other hand, since Q contains zero-mean independent and identically distributed Gaussian samples, we get: E{x H (r)1 q (r)p} = E{q H (r)p q (r)q} = 0 (A.8) and E{q H (r)pq (r)p } = M M r σ 2 Q (A.9) for any p q, p,q = 1,2,,M r, which result in K 3 = K 4 = 0 (Mr 1) (M r 1) and E{ξ (r)m ξ H (r)n} = σq 2(1+ ρ r 2 )I M, m = n Mr σq 2ρ r I, M Mr m = n+1 σq 2ρ ri M, Mr m = n 1 0, otherwise (A.10) Therefore (1+ ρ r 2 ) ρ r ρ K 1 = M r (1+ ρ r 2 ) ρ r 0 0 γ 2 σq M r 0 0 ρ r (1+ ρ r 2 ) ρ ρ r (1+ ρ r 2 ) (A.11) 15
16 TocalculateK 2, wenoticethatwhenp 1, q (r)1 andξ (r)p areindependent to each other, meaning that E{q H (r)1 ξ (r)p} = 0. Hence [K 2 ] m,n = [K 1 ] m,n σq 2/ γ 2 for any m 1 or n 1, and [K 2 ] 1,1 = E{q H (r)1 ξ (r)1ξ H (r)1 q (r)1} = M/M r (1 + M/M r ρ r 2 )σq 4. Substituting K 1, K 2 and ρ r = e jωr into (A.3) yields B(ω r ) = M M r ( γ 2 +σ 2 Q)σ 2 QP(ω r ) (A.12) where 2+(M/M r 1)σQ 2/( γ 2 +σq 2) ejωr e jωr 2 e jωr 0 0 P(ω r ) = e 2 e jωr e jωr 2 (A.13) As M/M r ( γ 2 + σq 2)σ2 Q is just a scalar, we can simply write W M r 1(ρ r ) = W Mr 1(ω r ) = P 1 (ω r ). Note that (7) (8) correspond to the special case of 2-D estimation. Appendix B The bias and variance of ˆω r, r = 1,2,,R, for the general model of (26) are now analyzed. Let the bias of ˆω r be ω r = ˆω r ω r. Following [18], for sufficiently high SNR conditions, we have ω r var(ˆω r ) κ r z H (r)1 W M r 1(ω r ) z (r)1 E{κ 2 r } (E{ z H (r)1 W M r 1(ω r ) z (r)1 }) 2 (B.1) (B.2) 16
17 where κ r I { e jωr ( z H (r)1w Mr 1(ω r ) z (r)2 + z H (r)1w Mr 1(ω r ) z (r)2 )} (B.3) with z (r)i = z (r)i + z (r)i, i = 1,2. Notice that z T (r) = (x (r)1 +q (r)1 ) H [x (r)1 x (r)2 x (r)mr ]and z T(r) = (x (r)1+q (r)1 ) H [q (r)1 q (r)2 q (r)mr ]. According to (A.8) (A.9), we obtain E{x H (r)(m) (x (r)1 +q (r)1 )(x (r)1 +q (r)1 ) H x (r)(n) } =x H (r)(m)x (r)1 x H (r)1x (r)(n) +x H (r)(m)e{q (r)1 )(q (r)1 ) H }x (r)(n) =e j(m n) M M r ( M M r γ 2 +σ 2 Q ) γ 2 (B.4) and [ ] T E{ z (r) } = M/M r σq (B.5) Therefore the expected value of the denominator of (B.1) is: E{ z H (r)1 W M r 1(ω r ) z (r)1 } M r 1 m=1 M r 1 n=1 E{x H (r)m (x (r)1 +q (r)1 )(M r min(m,n) mn)e j(m n)ω (x (r)1 +q (r)1 ) H x (r)n } ( = M ( M Mr 1 γ 2 +σ 2 Mr 2 Q M ) γ 2 (M r m m 2 )+ r m=1 ) M r 1 n 1 (M r m mn) n=1 m=1 M r 1 m=1 m 1 (M r n mn)+ = M(M2 r 1) ( M γ 2 +σ 12 M Q) γ 2 2 (B.6) r On the other hand, κ r can be written as: n=1 κ r =t r1 +t r2 (B.7) 17
18 where tr1 = e jωr z H (r)1 W M r 1(ω r ) z (r)2 e jωr z H (r)2 W M r 1(ω r ) z (r)1 2j = zh (r)2 W M r 1(ω r )( z (r)2 ρ r z (r)1 ) 2j (B.8) t r2 = e jωr z T (r)2 W M r 1 (ω r) z (r)1 ejωr z T (r)1 W M r 1 (ω r) z (r)2 2j = ( z (r)2 ρ r z (r)1 ) H W Mr 1(ω r ) z (r)2 2j (B.9) because z (r)2 = ρ r z (r)1, ρ r = e jωr and W Mr 1(ω r ) = W H M r 1 (ω r). As E{ z (r)2 } = 0 (Mr 1) 1, we have: E{t r1 } = 1 2j E{ρ r z H (r)2 W M r 1(ω r ) z (r)1 }+0 = 1 2j E{ zh (r)1 W M r 1(ω r ) z (r)1 } = 1 2j M r 1 m=1 M r 1 n=1 (x (r)1 +q (r)1 ) H q (r)n } = 1 2j = 1 2j = 1 2j M r 1 m=1 M r 1 m=1 M r 1 n=1 E{x H (r)m (x (r)1 +q (r)1 )(M r min(m,n) mn)e j(m n)ω E{x H (r)m(x (r)1 +q (r)1 )(M r m)e j(1 m)ω (x (r)1 +q (r)1 ) H q (r)1 }+0 (M r m) E{x H (r)1 (x (r)1 +q (r)1 )(x (r)1 +q (r)1 ) H q (r)1 } (M r n) E{x H (r)1 x (r)1q H (r)1 q (r)1}+0 =j M2 (M r 1) 4M r γ 2 σ 2 Q (B.10) and E{t r2 } = E{t r1 }. This means that E{ ω r } E{κ r } = 0 when SNR is sufficiently high, indicating the asymptotically unbiasedness of ˆω r. To calculate (B.2), we notice that E{κ 2 r } =E{t r1t r2 }+E{t r2t r1 }+E{t r1t r1 }+E{t r2t r2 } =2E{t r1 t r2 }+E{t 2 r1 }+E{t2 r2 } (B.11) 18
19 where E{t 2 r1} = E{t 2 r2} = M4 (M r 1) 2 γ 4 σ 4 16Mr 2 Q (B.12) and E{t r1 t r2 } = 1 4 E{ zh (r)2 W M r 1(ω r )( z (r)2 ρ r z (r)1 )( z (r)2 ρ r z (r)1 ) H W Mr 1(ω r ) z (r)2 } 1 4 E{ zh (r)2 W M r 1(ω r )B(ω r )W Mr 1(ω r ) z (r)2 } = 1 4 E{ zh (r)2 W M r 1(ω r ) M M r ( γ 2 +σ 2 Q )σ2 Q W 1 M r 1 (ω r)w Mr 1(ω r ) z (r)2 } = M 4M r ( γ 2 +σ 2 Q )σ2 Q E{ zh (r)1 W M r 1(ω r ) z (r)1 } (B.13) For sufficiently high SNRs, E{t 2 r1} E{t r1 t r2 }, and substituting (B.6) and (B.11) (B.13) into (B.2), we obtain (18). Acknowledgement The work described in this paper was fully supported by a grant from the Research Grants Council of the Hong Kong Special Administrative Region, China [Project No. CityU ]. References [1] S.L. Marple, Digital Spectral Analysis with Applications, Englewood Cliffs, NJ: Prentice-Hall, 1987 [2] S.M. Kay, Modern Spectral Estimation: Theory and Application, Englewood Cliffs, NJ: Prentice-Hall, 1988 [3] P. Stoica and R. Moses, Spectral Analysis of Signals, Upper Saddle River, NJ: Prentice-Hall,
20 [4] F.K.W. Chan, H.C. So, W.H. Lau and C.F. Chan, Structured total least squares approach for efficient frequency estimation, Signal Processing, vol.91, no.4, pp , Apr [5] K.W.K. Lui and H.C. So, Two-stage autocorrelation approach for accurate single sinusoidal frequency estimation, Signal Processing, vol.88, no.7, pp , Jul [6] Z. Zhou and H.C. So, Linear prediction approach to oversampling parameter estimation for multiple complex sinusoids, Signal Processing, vol.92, no.6, pp , Jun [7] J.H. McClellan, Multidimensional spectral estimation, Proceedings of the IEEE, vol.70, no.9, pp , Sep [8] M. Haardt and J.A. Nossek, Simultaneous Schur decomposition of several nonsymmetric matrices to achieve automatic pairing in multidimensional harmonic retrieval problems, IEEE Transactions on Signal Processing, vol.40, no.1, pp , Jan [9] H.L. van Trees, Optimum Array Processing: Detection, Estimation, and Modulation Theory. New York: Wiley, 2002, pt. IV [10] Y. Li, J. Razavilar and K.J.R. Liu, A high-resolution technique for multidimensional NMR spectroscopy, IEEE Transactions on Biomedical Engineering, vol.45, pp.78-86, Jan [11] K.N. Mokios, N.D. Sidiropoulos, M. Pesavento and C.F. Mecklenbrauker, On 3-D harmonic retrieval for wireless channel sounding, 20
21 Proc., IEEE International Conference on Acoustics, Speech and Signal Processing, vol. 2, pp , Montreal, Quebec, Canada, May 2004 [12] M. Pesavento, C.F. Mecklenbraker and J.F. Bome, Multi-dimensional rank reduction estimator for parametric MIMO channel models, EURASIP J. Appl. Signal Process., vol. 2004, pp , Sep [13] D. Nion and N.D. Sidiropoulos, Tensor algebra and multidimensional harmonic retrieval in signal processing for MIMO radar IEEE Transactions on Signal Processing, vol.58, no.11, pp , Nov [14] M.P. Clark and L.L. Scharf, Two-dimensional modal analysis based on maximum likelihood, IEEE Transactions on Signal Processing, vol.42, no.6, pp , Jun [15] H.C. So and F.K.W. Chan, Approximate maximum-likelihood algorithms for two-dimensional frequency estimation of a complex sinusoid, IEEE Transactions on Signal Processing, vol.54, no.8, pp , Aug [16] J. Liu and X. Liu, An eigenvector-based approach for multidimensional frequency estimation with improved identifiability, IEEE Transactions on Signal Processing, vol.54, no.12, pp , Dec [17] R. Boyer, Decoupled root-music algorithm for multidimensional harmonic retrieval, IEEE Workshop on Signal Processing Advances in Wireless Communications, pp , Recife, Brazil, Jul
22 [18] H.C. So and F.K.W. Chan, A generalized weighted linear predictor frequency estimation approach for a complex sinusoid, IEEE Transactions on Signal Processing, vol.54, no.4, pp , Apr [19] L. de Lathauwer, B. de Moor, and J. Vanderwalle, A multilinear singular value decomposition, SIAM J. Matrix Anal. Appl., vol. 21, no. 4, pp , 2000 [20] R. Boyer, Deterministic asymptotic Cramer-Rao lower bound for the multidimensional harmonic model, Signal Processing, vol. 88, no. 12, pp , Dec [21] S. Kay, A fast and accurate single frequency estimator, IEEE Trans. Acoust., Speech, Signal Processing, vol.37, pp , Dec [22] S.M. Kay, Fundamentals of Statistical Signal Processing: Estimation Theory, Englewood Cliffs, NJ: Prentice Hall,
23 0 10 C 1 C 2 (13) (17) CRLB Mean Square Error (db) N=10 N=5 60 N= SNR (db) Figure 1: Average mean square frequency error versus SNR at M 1 = M 2 = N 23
24 Mean Square Error (db) N=5 N=10 N=20 C 1 C 2 (13) (17) CRLB SNR (db) Figure 2: Average mean square frequency error versus SNR at M 1 = M 2 = M 3 = N 24
25 Mean Square Error (db) C 1 C 2 AIQML [15] IMDF [16] (13) (17) CRLB SNR (db) Figure 3: Average mean square frequency error versus SNR for 2-D case 25
26 Mean Square Error (db) C 1 C 2 UE [8] IMDF [16] (13) (17) CRLB SNR (db) Figure 4: Average mean square frequency error versus SNR for 3-D case 26
27 Mean Square Error (db) C 1 C 2 UE [8] IMDF [16] (13) (17) CRLB Number of Dimensions Figure 5: Average mean square frequency error versus R at M r = 4 and SNR= 20dB 27
Multidimensional Sinusoidal Frequency Estimation Using Subspace and Projection Separation Approaches
Multidimensional Sinusoidal Frequency Estimation Using Subspace and Projection Separation Approaches 1 Longting Huang, Yuntao Wu, H. C. So, Yanduo Zhang and Lei Huang Department of Electronic Engineering,
More informationFAST AND ACCURATE DIRECTION-OF-ARRIVAL ESTIMATION FOR A SINGLE SOURCE
Progress In Electromagnetics Research C, Vol. 6, 13 20, 2009 FAST AND ACCURATE DIRECTION-OF-ARRIVAL ESTIMATION FOR A SINGLE SOURCE Y. Wu School of Computer Science and Engineering Wuhan Institute of Technology
More informationTensor-based Subspace Approach for Multi-Dimensional Harmonic Retrieval
Accurate and Computationally Efficient 1 Tensor-based Subspace Approach for Multi-Dimensional Harmonic etrieval Weize Sun and H.C. So Abstract In this paper, parameter estimation for -dimensional (-D)
More informationProgress In Electromagnetics Research, PIER 102, 31 48, 2010
Progress In Electromagnetics Research, PIER 102, 31 48, 2010 ACCURATE PARAMETER ESTIMATION FOR WAVE EQUATION F. K. W. Chan, H. C. So, S.-C. Chan, W. H. Lau and C. F. Chan Department of Electronic Engineering
More information2290 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 53, NO. 7, JULY 2005
2290 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL 53, NO 7, JULY 2005 Linear Prediction Approach for Efficient Frequency Estimation of Multiple Real Sinusoids: Algorithms and Analyses H C So, Member, IEEE,
More informationSubspace Approach for Two-Dimensional Parameter Estimation of Multiple Damped Sinusoids
Subspace Approach for Two-Dimensional Parameter Estimation of Multiple Damped Sinusoids 1 Frankie K.W. Chan, H.C. So and Weize Sun ckwf@hkexperts.com, hcso@ee.cityu.edu.hk, weizesun2@student.cityu.edu.hk
More informationCOMPLEX CONSTRAINED CRB AND ITS APPLICATION TO SEMI-BLIND MIMO AND OFDM CHANNEL ESTIMATION. Aditya K. Jagannatham and Bhaskar D.
COMPLEX CONSTRAINED CRB AND ITS APPLICATION TO SEMI-BLIND MIMO AND OFDM CHANNEL ESTIMATION Aditya K Jagannatham and Bhaskar D Rao University of California, SanDiego 9500 Gilman Drive, La Jolla, CA 92093-0407
More information/16/$ IEEE 1728
Extension of the Semi-Algebraic Framework for Approximate CP Decompositions via Simultaneous Matrix Diagonalization to the Efficient Calculation of Coupled CP Decompositions Kristina Naskovska and Martin
More informationOn Moving Average Parameter Estimation
On Moving Average Parameter Estimation Niclas Sandgren and Petre Stoica Contact information: niclas.sandgren@it.uu.se, tel: +46 8 473392 Abstract Estimation of the autoregressive moving average (ARMA)
More informationFrequency estimation by DFT interpolation: A comparison of methods
Frequency estimation by DFT interpolation: A comparison of methods Bernd Bischl, Uwe Ligges, Claus Weihs March 5, 009 Abstract This article comments on a frequency estimator which was proposed by [6] and
More informationRobust Implementation of the MUSIC algorithm Zhang, Johan Xi; Christensen, Mads Græsbøll; Dahl, Joachim; Jensen, Søren Holdt; Moonen, Marc
Aalborg Universitet Robust Implementation of the MUSIC algorithm Zhang, Johan Xi; Christensen, Mads Græsbøll; Dahl, Joachim; Jensen, Søren Holdt; Moonen, Marc Published in: I E E E International Conference
More informationLinear Prediction Approach to Oversampling Parameter Estimation for Multiple Complex Sinusoids
Linear Prediction Approach to Oversampling Parameter Estimation for Multiple Complex Sinusoids Zhenhua ZHOU and H C SO Department of Electronic Engineering City University of Hong Kong, Tat Chee Avenue,
More informationImproved Unitary Root-MUSIC for DOA Estimation Based on Pseudo-Noise Resampling
140 IEEE SIGNAL PROCESSING LETTERS, VOL. 21, NO. 2, FEBRUARY 2014 Improved Unitary Root-MUSIC for DOA Estimation Based on Pseudo-Noise Resampling Cheng Qian, Lei Huang, and H. C. So Abstract A novel pseudo-noise
More informationTensor MUSIC in Multidimensional Sparse Arrays
Tensor MUSIC in Multidimensional Sparse Arrays Chun-Lin Liu 1 and P. P. Vaidyanathan 2 Dept. of Electrical Engineering, MC 136-93 California Institute of Technology, Pasadena, CA 91125, USA cl.liu@caltech.edu
More informationthe tensor rank is equal tor. The vectorsf (r)
EXTENSION OF THE SEMI-ALGEBRAIC FRAMEWORK FOR APPROXIMATE CP DECOMPOSITIONS VIA NON-SYMMETRIC SIMULTANEOUS MATRIX DIAGONALIZATION Kristina Naskovska, Martin Haardt Ilmenau University of Technology Communications
More informationsine wave fit algorithm
TECHNICAL REPORT IR-S3-SB-9 1 Properties of the IEEE-STD-57 four parameter sine wave fit algorithm Peter Händel, Senior Member, IEEE Abstract The IEEE Standard 57 (IEEE-STD-57) provides algorithms for
More informationJ. Liang School of Automation & Information Engineering Xi an University of Technology, China
Progress In Electromagnetics Research C, Vol. 18, 245 255, 211 A NOVEL DIAGONAL LOADING METHOD FOR ROBUST ADAPTIVE BEAMFORMING W. Wang and R. Wu Tianjin Key Lab for Advanced Signal Processing Civil Aviation
More informationDesign Criteria for the Quadratically Interpolated FFT Method (I): Bias due to Interpolation
CENTER FOR COMPUTER RESEARCH IN MUSIC AND ACOUSTICS DEPARTMENT OF MUSIC, STANFORD UNIVERSITY REPORT NO. STAN-M-4 Design Criteria for the Quadratically Interpolated FFT Method (I): Bias due to Interpolation
More informationCP DECOMPOSITION AND ITS APPLICATION IN NOISE REDUCTION AND MULTIPLE SOURCES IDENTIFICATION
International Conference on Computer Science and Intelligent Communication (CSIC ) CP DECOMPOSITION AND ITS APPLICATION IN NOISE REDUCTION AND MULTIPLE SOURCES IDENTIFICATION Xuefeng LIU, Yuping FENG,
More informationEfficient and Accurate Rectangular Window Subspace Tracking
Efficient and Accurate Rectangular Window Subspace Tracking Timothy M. Toolan and Donald W. Tufts Dept. of Electrical Engineering, University of Rhode Island, Kingston, RI 88 USA toolan@ele.uri.edu, tufts@ele.uri.edu
More informationEstimation of the Optimum Rotational Parameter for the Fractional Fourier Transform Using Domain Decomposition
Estimation of the Optimum Rotational Parameter for the Fractional Fourier Transform Using Domain Decomposition Seema Sud 1 1 The Aerospace Corporation, 4851 Stonecroft Blvd. Chantilly, VA 20151 Abstract
More informationMULTIPLE-CHANNEL DETECTION IN ACTIVE SENSING. Kaitlyn Beaudet and Douglas Cochran
MULTIPLE-CHANNEL DETECTION IN ACTIVE SENSING Kaitlyn Beaudet and Douglas Cochran School of Electrical, Computer and Energy Engineering Arizona State University, Tempe AZ 85287-576 USA ABSTRACT The problem
More informationOptimal Time Division Multiplexing Schemes for DOA Estimation of a Moving Target Using a Colocated MIMO Radar
Optimal Division Multiplexing Schemes for DOA Estimation of a Moving Target Using a Colocated MIMO Radar Kilian Rambach, Markus Vogel and Bin Yang Institute of Signal Processing and System Theory University
More informationPASSIVE NEAR-FIELD SOURCE LOCALIZATION BASED ON SPATIAL-TEMPORAL STRUCTURE
Progress In Electromagnetics Research C, Vol. 8, 27 41, 29 PASSIVE NEAR-FIELD SOURCE LOCALIZATION BASED ON SPATIAL-TEMPORAL STRUCTURE Y. Wu Wuhan Institute of Technology Wuhan 4373, China H. C. So Department
More informationA New Subspace Identification Method for Open and Closed Loop Data
A New Subspace Identification Method for Open and Closed Loop Data Magnus Jansson July 2005 IR S3 SB 0524 IFAC World Congress 2005 ROYAL INSTITUTE OF TECHNOLOGY Department of Signals, Sensors & Systems
More informationA Subspace Approach to Estimation of. Measurements 1. Carlos E. Davila. Electrical Engineering Department, Southern Methodist University
EDICS category SP 1 A Subspace Approach to Estimation of Autoregressive Parameters From Noisy Measurements 1 Carlos E Davila Electrical Engineering Department, Southern Methodist University Dallas, Texas
More informationA New High-Resolution and Stable MV-SVD Algorithm for Coherent Signals Detection
Progress In Electromagnetics Research M, Vol. 35, 163 171, 2014 A New High-Resolution and Stable MV-SVD Algorithm for Coherent Signals Detection Basma Eldosouky, Amr H. Hussein *, and Salah Khamis Abstract
More informationA CLOSED FORM SOLUTION TO SEMI-BLIND JOINT SYMBOL AND CHANNEL ESTIMATION IN MIMO-OFDM SYSTEMS
A CLOSED FORM SOLUTION TO SEMI-BLIND JOINT SYMBOL AND CHANNEL ESTIMATION IN MIMO-OFDM SYSTEMS Kefei Liu 1, João Paulo C. L. da Costa 2, André L. F. de Almeida 3, and H.C. So 1 1 Department of Electronic
More informationIEEE copyright notice
IEEE copyright notice Personal use of this material is permitted. However, permission to reprint/republish this material for advertising or promotional purposes or for creating new collective works for
More informationLOW COMPLEXITY COVARIANCE-BASED DOA ESTIMATION ALGORITHM
LOW COMPLEXITY COVARIANCE-BASED DOA ESTIMATION ALGORITHM Tadeu N. Ferreira, Sergio L. Netto, and Paulo S. R. Diniz Electrical Engineering Program COPPE/DEL-Poli/Federal University of Rio de Janeiro P.O.
More informationLocal Strong Convexity of Maximum-Likelihood TDOA-Based Source Localization and Its Algorithmic Implications
Local Strong Convexity of Maximum-Likelihood TDOA-Based Source Localization and Its Algorithmic Implications Huikang Liu, Yuen-Man Pun, and Anthony Man-Cho So Dept of Syst Eng & Eng Manag, The Chinese
More informationEIE6207: Estimation Theory
EIE6207: Estimation Theory Man-Wai MAK Dept. of Electronic and Information Engineering, The Hong Kong Polytechnic University enmwmak@polyu.edu.hk http://www.eie.polyu.edu.hk/ mwmak References: Steven M.
More informationDOA Estimation of Quasi-Stationary Signals Using a Partly-Calibrated Uniform Linear Array with Fewer Sensors than Sources
Progress In Electromagnetics Research M, Vol. 63, 185 193, 218 DOA Estimation of Quasi-Stationary Signals Using a Partly-Calibrated Uniform Linear Array with Fewer Sensors than Sources Kai-Chieh Hsu and
More informationCORE CONSISTENCY DIAGNOSTIC AIDED BY RECONSTRUCTION ERROR FOR ACCURATE ENUMERATION OF THE NUMBER OF COMPONENTS IN PARAFAC MODELS
CORE CONSISTENCY DIAGNOSTIC AIDED BY RECONSTRUCTION ERROR FOR ACCURATE ENUMERATION OF THE NUMBER OF COMPONENTS IN PARAFAC MODELS Kefei Liu 1, H.C. So 1, João Paulo C. L. da Costa and Lei Huang 3 1 Department
More informationDOA Estimation using MUSIC and Root MUSIC Methods
DOA Estimation using MUSIC and Root MUSIC Methods EE602 Statistical signal Processing 4/13/2009 Presented By: Chhavipreet Singh(Y515) Siddharth Sahoo(Y5827447) 2 Table of Contents 1 Introduction... 3 2
More informationMean Likelihood Frequency Estimation
IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 48, NO. 7, JULY 2000 1937 Mean Likelihood Frequency Estimation Steven Kay, Fellow, IEEE, and Supratim Saha Abstract Estimation of signals with nonlinear as
More informationX. Zhang, G. Feng, and D. Xu Department of Electronic Engineering Nanjing University of Aeronautics & Astronautics Nanjing , China
Progress In Electromagnetics Research Letters, Vol. 13, 11 20, 2010 BLIND DIRECTION OF ANGLE AND TIME DELAY ESTIMATION ALGORITHM FOR UNIFORM LINEAR ARRAY EMPLOYING MULTI-INVARIANCE MUSIC X. Zhang, G. Feng,
More informationPERFORMANCE ANALYSIS OF COARRAY-BASED MUSIC AND THE CRAMÉR-RAO BOUND. Mianzhi Wang, Zhen Zhang, and Arye Nehorai
PERFORANCE ANALYSIS OF COARRAY-BASED USIC AND THE CRAÉR-RAO BOUND ianzhi Wang, Zhen Zhang, and Arye Nehorai Preston. Green Department of Electrical and Systems Engineering, Washington University in St.
More informationMANY digital speech communication applications, e.g.,
406 IEEE TRANSACTIONS ON AUDIO, SPEECH, AND LANGUAGE PROCESSING, VOL. 15, NO. 2, FEBRUARY 2007 An MMSE Estimator for Speech Enhancement Under a Combined Stochastic Deterministic Speech Model Richard C.
More informationFast adaptive ESPRIT algorithm
Fast adaptive ESPRIT algorithm Roland Badeau, Gaël Richard, Bertrand David To cite this version: Roland Badeau, Gaël Richard, Bertrand David. Fast adaptive ESPRIT algorithm. Proc. of IEEE Workshop on Statistical
More informationDirection of Arrival Estimation: Subspace Methods. Bhaskar D Rao University of California, San Diego
Direction of Arrival Estimation: Subspace Methods Bhaskar D Rao University of California, San Diego Email: brao@ucsdedu Reference Books and Papers 1 Optimum Array Processing, H L Van Trees 2 Stoica, P,
More informationUpper Bounds on MIMO Channel Capacity with Channel Frobenius Norm Constraints
Upper Bounds on IO Channel Capacity with Channel Frobenius Norm Constraints Zukang Shen, Jeffrey G. Andrews, Brian L. Evans Wireless Networking Communications Group Department of Electrical Computer Engineering
More information2-D SENSOR POSITION PERTURBATION ANALYSIS: EQUIVALENCE TO AWGN ON ARRAY OUTPUTS. Volkan Cevher, James H. McClellan
2-D SENSOR POSITION PERTURBATION ANALYSIS: EQUIVALENCE TO AWGN ON ARRAY OUTPUTS Volkan Cevher, James H McClellan Georgia Institute of Technology Atlanta, GA 30332-0250 cevher@ieeeorg, jimmcclellan@ecegatechedu
More informationA METHOD OF ADAPTATION BETWEEN STEEPEST- DESCENT AND NEWTON S ALGORITHM FOR MULTI- CHANNEL ACTIVE CONTROL OF TONAL NOISE AND VIBRATION
A METHOD OF ADAPTATION BETWEEN STEEPEST- DESCENT AND NEWTON S ALGORITHM FOR MULTI- CHANNEL ACTIVE CONTROL OF TONAL NOISE AND VIBRATION Jordan Cheer and Stephen Daley Institute of Sound and Vibration Research,
More informationA Fast Algorithm for. Nonstationary Delay Estimation. H. C. So. Department of Electronic Engineering, City University of Hong Kong
A Fast Algorithm for Nonstationary Delay Estimation H. C. So Department of Electronic Engineering, City University of Hong Kong Tat Chee Avenue, Kowloon, Hong Kong Email : hcso@ee.cityu.edu.hk June 19,
More information-Dimensional ESPRIT-Type Algorithms for Strictly Second-Order Non-Circular Sources and Their Performance Analysis
4824 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 62, NO. 18, SEPTEMBER 15, 2014 -Dimensional ESPRIT-Type Algorithms for Strictly Second-Order Non-Circular Sources Their Performance Analysis Jens Steinwt,
More informationSIGNALS with sparse representations can be recovered
IEEE SIGNAL PROCESSING LETTERS, VOL. 22, NO. 9, SEPTEMBER 2015 1497 Cramér Rao Bound for Sparse Signals Fitting the Low-Rank Model with Small Number of Parameters Mahdi Shaghaghi, Student Member, IEEE,
More informationThis is the published version of a paper presented at IEEE International Instrumentation and Measurement Technology Conference (I2MTC), 2013.
http://www.diva-portal.org This is the published version of a paper presented at IEEE International Instrumentation and Measurement Technology Conference (IMTC), 3. Citation for the original published
More information5682 IEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 55, NO. 12, DECEMBER /$ IEEE
5682 IEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 55, NO. 12, DECEMBER 2009 Hyperplane-Based Vector Quantization for Distributed Estimation in Wireless Sensor Networks Jun Fang, Member, IEEE, and Hongbin
More informationThreshold Effects in Parameter Estimation from Compressed Data
PAKROOH et al.: THRESHOLD EFFECTS IN PARAMETER ESTIMATION FROM COMPRESSED DATA 1 Threshold Effects in Parameter Estimation from Compressed Data Pooria Pakrooh, Student Member, IEEE, Louis L. Scharf, Life
More informationFast multilinear Singular Values Decomposition for higher-order Hankel tensors
Fast multilinear Singular Values Decomposition for higher-order Hanel tensors Maxime Boizard, Remy Boyer, Gérard Favier, Pascal Larzabal To cite this version: Maxime Boizard, Remy Boyer, Gérard Favier,
More informationGeneralization Propagator Method for DOA Estimation
Progress In Electromagnetics Research M, Vol. 37, 119 125, 2014 Generalization Propagator Method for DOA Estimation Sheng Liu, Li Sheng Yang, Jian ua uang, and Qing Ping Jiang * Abstract A generalization
More informationEXTENSION OF NESTED ARRAYS WITH THE FOURTH-ORDER DIFFERENCE CO-ARRAY ENHANCEMENT
EXTENSION OF NESTED ARRAYS WITH THE FOURTH-ORDER DIFFERENCE CO-ARRAY ENHANCEMENT Qing Shen,2, Wei Liu 2, Wei Cui, Siliang Wu School of Information and Electronics, Beijing Institute of Technology Beijing,
More informationADAPTIVE FILTER THEORY
ADAPTIVE FILTER THEORY Fourth Edition Simon Haykin Communications Research Laboratory McMaster University Hamilton, Ontario, Canada Front ice Hall PRENTICE HALL Upper Saddle River, New Jersey 07458 Preface
More informationPitch Estimation and Tracking with Harmonic Emphasis On The Acoustic Spectrum
Downloaded from vbn.aau.dk on: marts 31, 2019 Aalborg Universitet Pitch Estimation and Tracking with Harmonic Emphasis On The Acoustic Spectrum Karimian-Azari, Sam; Mohammadiha, Nasser; Jensen, Jesper
More informationTwo-Dimensional Sparse Arrays with Hole-Free Coarray and Reduced Mutual Coupling
Two-Dimensional Sparse Arrays with Hole-Free Coarray and Reduced Mutual Coupling Chun-Lin Liu and P. P. Vaidyanathan Dept. of Electrical Engineering, 136-93 California Institute of Technology, Pasadena,
More informationRobust Subspace DOA Estimation for Wireless Communications
Robust Subspace DOA Estimation for Wireless Communications Samuli Visuri Hannu Oja ¾ Visa Koivunen Laboratory of Signal Processing Computer Technology Helsinki Univ. of Technology P.O. Box 3, FIN-25 HUT
More informationLocal detectors for high-resolution spectral analysis: Algorithms and performance
Digital Signal Processing 15 (2005 305 316 www.elsevier.com/locate/dsp Local detectors for high-resolution spectral analysis: Algorithms and performance Morteza Shahram, Peyman Milanfar Electrical Engineering
More informationReal-Valued Khatri-Rao Subspace Approaches on the ULA and a New Nested Array
Real-Valued Khatri-Rao Subspace Approaches on the ULA and a New Nested Array Huiping Duan, Tiantian Tuo, Jun Fang and Bing Zeng arxiv:1511.06828v1 [cs.it] 21 Nov 2015 Abstract In underdetermined direction-of-arrival
More information2542 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 62, NO. 10, MAY 15, Keyong Han, Student Member, IEEE, and Arye Nehorai, Fellow, IEEE
2542 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 62, NO. 10, MAY 15, 2014 Nested Vector-Sensor Array Processing via Tensor Modeling Keyong Han, Student Member, IEEE, and Arye Nehorai, Fellow, IEEE Abstract
More informationCOMPARISON OF MODEL ORDER SELECTION TECHNIQUES FOR HIGH-RESOLUTION PARAMETER ESTIMATION ALGORITHMS
COMPARISON OF MODEL ORDER SELECTION TECHNIQUES FOR HIGH-RESOLUTION PARAMETER ESTIMATION ALGORITHMS João Paulo C. L. da Costa, Arpita Thare 2, Florian Roemer, and Martin Haardt Ilmenau University of Technology
More information4674 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 58, NO. 9, SEPTEMBER 2010
4674 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 58, NO. 9, SEPTEMBER 2010 A Low Complexity Blind Estimator of Narrowband Polynomial Phase Signals Alon Amar, Member, IEEE, Amir Leshem, Senior Member,
More informationarxiv: v1 [math.st] 13 Dec 2016
BEST WIDELY LINEAR UNBIASED ESTIMATOR FOR REAL VALUED PARAMETER VECTORS O. Lang and M. Huemer arxiv:1612.04060v1 [math.st 13 Dec 2016 ABSTRACT For classical estimation with an underlying linear model the
More informationEFFICIENT CONSTRUCTION OF POLYPHASE SEQUENCES WITH OPTIMAL PEAK SIDELOBE LEVEL GROWTH. Arindam Bose and Mojtaba Soltanalian
EFFICIENT CONSTRUCTION OF POLYPHASE SEQUENCES WITH OPTIAL PEAK SIDELOBE LEVEL GROWTH Arindam Bose and ojtaba Soltanalian Department of Electrical and Computer Engineering University of Illinois at Chicago
More informationJoint High-Resolution Fundamental Frequency and Order Estimation. Christensen, Mads Græsbøll; Jakobsson, Andreas; Jensen, Søren Holdt
Downloaded from vbnaaudk on: januar 18, 219 Aalborg Universitet Joint High-Resolution Fundamental Frequency and Order Estimation Christensen, Mads Græsbøll; Jakobsson, Andreas; Jensen, Søren Holdt Published
More informationIMPROVED BLIND 2D-DIRECTION OF ARRIVAL ESTI- MATION WITH L-SHAPED ARRAY USING SHIFT IN- VARIANCE PROPERTY
J. of Electromagn. Waves and Appl., Vol. 23, 593 606, 2009 IMPROVED BLIND 2D-DIRECTION OF ARRIVAL ESTI- MATION WITH L-SHAPED ARRAY USING SHIFT IN- VARIANCE PROPERTY X. Zhang, X. Gao, and W. Chen Department
More informationTHE PROBLEM of estimating the parameters of a twodimensional
IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL 46, NO 8, AUGUST 1998 2157 Parameter Estimation of Two-Dimensional Moving Average Rom Fields Joseph M Francos, Senior Member, IEEE, Benjamin Friedler, Fellow,
More informationPerformance Analysis of an Adaptive Algorithm for DOA Estimation
Performance Analysis of an Adaptive Algorithm for DOA Estimation Assimakis K. Leros and Vassilios C. Moussas Abstract This paper presents an adaptive approach to the problem of estimating the direction
More informationPerformance Comparison of Two Implementations of the Leaky. LMS Adaptive Filter. Scott C. Douglas. University of Utah. Salt Lake City, Utah 84112
Performance Comparison of Two Implementations of the Leaky LMS Adaptive Filter Scott C. Douglas Department of Electrical Engineering University of Utah Salt Lake City, Utah 8411 Abstract{ The leaky LMS
More informationStatistical and Adaptive Signal Processing
r Statistical and Adaptive Signal Processing Spectral Estimation, Signal Modeling, Adaptive Filtering and Array Processing Dimitris G. Manolakis Massachusetts Institute of Technology Lincoln Laboratory
More informationGeneralized Triangular Transform Coding
Generalized Triangular Transform Coding Ching-Chih Weng, Chun-Yang Chen, and P. P. Vaidyanathan Dept. of Electrical Engineering, MC 136-93 California Institute of Technology, Pasadena, CA 91125, USA E-mail:
More informationDetection & Estimation Lecture 1
Detection & Estimation Lecture 1 Intro, MVUE, CRLB Xiliang Luo General Course Information Textbooks & References Fundamentals of Statistical Signal Processing: Estimation Theory/Detection Theory, Steven
More informationNew Phase Difference Measurement Method for Non-Integer Number of Signal Periods Based on Multiple Cross-Correlations
Send Orders for Reprints to reprints@benthamscience.ae The Open Automation and Control Systems Journal, 15, 7, 1537-154537 Open Access ew Phase Difference Measurement Method for on-integer umber of Signal
More informationVirtual Array Processing for Active Radar and Sonar Sensing
SCHARF AND PEZESHKI: VIRTUAL ARRAY PROCESSING FOR ACTIVE SENSING Virtual Array Processing for Active Radar and Sonar Sensing Louis L. Scharf and Ali Pezeshki Abstract In this paper, we describe how an
More informationDiffusion LMS Algorithms for Sensor Networks over Non-ideal Inter-sensor Wireless Channels
Diffusion LMS Algorithms for Sensor Networs over Non-ideal Inter-sensor Wireless Channels Reza Abdolee and Benoit Champagne Electrical and Computer Engineering McGill University 3480 University Street
More informationJoint Direction-of-Arrival and Order Estimation in Compressed Sensing using Angles between Subspaces
Aalborg Universitet Joint Direction-of-Arrival and Order Estimation in Compressed Sensing using Angles between Subspaces Christensen, Mads Græsbøll; Nielsen, Jesper Kjær Published in: I E E E / S P Workshop
More informationOptimum Sampling Vectors for Wiener Filter Noise Reduction
58 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 50, NO. 1, JANUARY 2002 Optimum Sampling Vectors for Wiener Filter Noise Reduction Yukihiko Yamashita, Member, IEEE Absact Sampling is a very important and
More informationOn the Convergence of the Inverses of Toeplitz Matrices and Its Applications
180 IEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 49, NO. 1, JANUARY 2003 On the Convergence of the Inverses of Toeplitz Matrices Its Applications Feng-Wen Sun, Member, IEEE, Yimin Jiang, Member, IEEE,
More informationTechnique for Numerical Computation of Cramér-Rao Bound using MATLAB
Technique for Numerical Computation of Cramér-Rao Bound using MATLAB Hing Cheung So http://www.ee.cityu.edu.hk/~hcso Department of Electronic Engineering City University of Hong Kong H. C. So Page 1 The
More informationA NEW BAYESIAN LOWER BOUND ON THE MEAN SQUARE ERROR OF ESTIMATORS. Koby Todros and Joseph Tabrikian
16th European Signal Processing Conference EUSIPCO 008) Lausanne Switzerland August 5-9 008 copyright by EURASIP A NEW BAYESIAN LOWER BOUND ON THE MEAN SQUARE ERROR OF ESTIMATORS Koby Todros and Joseph
More informationSUBSPACE based direction-of-arrival (DOA) estimation
IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL 64, NO 16, AUGUST 15, 2016 4127 Enhanced PUMA for Direction-of-Arrival Estimation and Its Performance Analysis Cheng Qian, Student Member, IEEE, LeiHuang, Senior
More informationAdaptive beamforming for uniform linear arrays with unknown mutual coupling. IEEE Antennas and Wireless Propagation Letters.
Title Adaptive beamforming for uniform linear arrays with unknown mutual coupling Author(s) Liao, B; Chan, SC Citation IEEE Antennas And Wireless Propagation Letters, 2012, v. 11, p. 464-467 Issued Date
More informationTruncation Strategy of Tensor Compressive Sensing for Noisy Video Sequences
Journal of Information Hiding and Multimedia Signal Processing c 2016 ISSN 207-4212 Ubiquitous International Volume 7, Number 5, September 2016 Truncation Strategy of Tensor Compressive Sensing for Noisy
More informationMAXIMUM A POSTERIORI ESTIMATION OF SIGNAL RANK. PO Box 1500, Edinburgh 5111, Australia. Arizona State University, Tempe AZ USA
MAXIMUM A POSTERIORI ESTIMATION OF SIGNAL RANK Songsri Sirianunpiboon Stephen D. Howard, Douglas Cochran 2 Defence Science Technology Organisation PO Box 500, Edinburgh 5, Australia 2 School of Mathematical
More informationFree Probability, Sample Covariance Matrices and Stochastic Eigen-Inference
Free Probability, Sample Covariance Matrices and Stochastic Eigen-Inference Alan Edelman Department of Mathematics, Computer Science and AI Laboratories. E-mail: edelman@math.mit.edu N. Raj Rao Deparment
More informationPerformance Analysis of Coarray-Based MUSIC and the Cramér-Rao Bound
Performance Analysis of Coarray-Based MUSIC and the Cramér-Rao Bound Mianzhi Wang, Zhen Zhang, and Arye Nehorai Preston M. Green Department of Electrical & Systems Engineering Washington University in
More informationSparse Sensing in Colocated MIMO Radar: A Matrix Completion Approach
Sparse Sensing in Colocated MIMO Radar: A Matrix Completion Approach Athina P. Petropulu Department of Electrical and Computer Engineering Rutgers, the State University of New Jersey Acknowledgments Shunqiao
More informationTitle without the persistently exciting c. works must be obtained from the IEE
Title Exact convergence analysis of adapt without the persistently exciting c Author(s) Sakai, H; Yang, JM; Oka, T Citation IEEE TRANSACTIONS ON SIGNAL 55(5): 2077-2083 PROCESS Issue Date 2007-05 URL http://hdl.handle.net/2433/50544
More informationSpatial Smoothing and Broadband Beamforming. Bhaskar D Rao University of California, San Diego
Spatial Smoothing and Broadband Beamforming Bhaskar D Rao University of California, San Diego Email: brao@ucsd.edu Reference Books and Papers 1. Optimum Array Processing, H. L. Van Trees 2. Stoica, P.,
More informationA Block-Jacobi Algorithm for Non-Symmetric Joint Diagonalization of Matrices
A Block-Jacobi Algorithm for Non-Symmetric Joint Diagonalization of Matrices ao Shen and Martin Kleinsteuber Department of Electrical and Computer Engineering Technische Universität München, Germany {hao.shen,kleinsteuber}@tum.de
More informationEfficient Line Search Methods for Riemannian Optimization Under Unitary Matrix Constraint
Efficient Line Search Methods for Riemannian Optimization Under Unitary Matrix Constraint Traian Abrudan, Jan Eriksson, Visa Koivunen SMARAD CoE, Signal Processing Laboratory, Helsinki University of Technology,
More informationAn Even Order Symmetric B Tensor is Positive Definite
An Even Order Symmetric B Tensor is Positive Definite Liqun Qi, Yisheng Song arxiv:1404.0452v4 [math.sp] 14 May 2014 October 17, 2018 Abstract It is easily checkable if a given tensor is a B tensor, or
More informationAsymptotic Analysis of the Generalized Coherence Estimate
IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 49, NO. 1, JANUARY 2001 45 Asymptotic Analysis of the Generalized Coherence Estimate Axel Clausen, Member, IEEE, and Douglas Cochran, Senior Member, IEEE Abstract
More informationNOISE ROBUST RELATIVE TRANSFER FUNCTION ESTIMATION. M. Schwab, P. Noll, and T. Sikora. Technical University Berlin, Germany Communication System Group
NOISE ROBUST RELATIVE TRANSFER FUNCTION ESTIMATION M. Schwab, P. Noll, and T. Sikora Technical University Berlin, Germany Communication System Group Einsteinufer 17, 1557 Berlin (Germany) {schwab noll
More informationInformation-Preserving Transformations for Signal Parameter Estimation
866 IEEE SIGNAL PROCESSING LETTERS, VOL. 21, NO. 7, JULY 2014 Information-Preserving Transformations for Signal Parameter Estimation Manuel Stein, Mario Castañeda, Amine Mezghani, and Josef A. Nossek Abstract
More informationECE Department, U. of Minnesota, Minneapolis, MN. 3SPICE-ECE-UMN N. Sidiropoulos / SIAM50
( *;:
More informationReceived Signal, Interference and Noise
Optimum Combining Maximum ratio combining (MRC) maximizes the output signal-to-noise ratio (SNR) and is the optimal combining method in a maximum likelihood sense for channels where the additive impairment
More informationThe Solvability Conditions for the Inverse Eigenvalue Problem of Hermitian and Generalized Skew-Hamiltonian Matrices and Its Approximation
The Solvability Conditions for the Inverse Eigenvalue Problem of Hermitian and Generalized Skew-Hamiltonian Matrices and Its Approximation Zheng-jian Bai Abstract In this paper, we first consider the inverse
More informationLecture 7 MIMO Communica2ons
Wireless Communications Lecture 7 MIMO Communica2ons Prof. Chun-Hung Liu Dept. of Electrical and Computer Engineering National Chiao Tung University Fall 2014 1 Outline MIMO Communications (Chapter 10
More informationA GLRT FOR RADAR DETECTION IN THE PRESENCE OF COMPOUND-GAUSSIAN CLUTTER AND ADDITIVE WHITE GAUSSIAN NOISE. James H. Michels. Bin Liu, Biao Chen
A GLRT FOR RADAR DETECTION IN THE PRESENCE OF COMPOUND-GAUSSIAN CLUTTER AND ADDITIVE WHITE GAUSSIAN NOISE Bin Liu, Biao Chen Syracuse University Dept of EECS, Syracuse, NY 3244 email : biliu{bichen}@ecs.syr.edu
More information