Generalized Sidelobe Canceller and MVDR Power Spectrum Estimation. Bhaskar D Rao University of California, San Diego
|
|
- Brittney Bates
- 5 years ago
- Views:
Transcription
1 Generalized Sidelobe Canceller and MVDR Power Spectrum Estimation Bhaskar D Rao University of California, San Diego brao@ucsd.edu
2 Reference Books 1. Optimum Array Processing, H. L. Van Trees 2. Stoica, P., & Moses, R. L. (2005). Spectral analysis of signals (Vol. 1). Upper Saddle River, NJ: Pearson Prentice Hall.
3 Narrow-band Signals Assumptions D 1 x(ω c, n) = V(ω c, k s )F s [n] + V(ω c, k l )F l [n] + Z[n] l=1 F s [n], F l [n], l = 1,.., D 1, and Z[n] are zero mean E( F s [n] 2 ) = p s, and E( F l [n] 2 = p l, l = 1,..., D 1, and E(Z[n]Z H [n]) = σ 2 z I All the signals/sources are uncorrelated with each other and over time: E(F l [n]f m[p]) = p l δ[l m]δ[n p] and E(F l [n]f s [p]) = 0 The sources are uncorrelated with the noise: E(Z[n]F l [m]) = 0
4 Relationship between MVDR, MPDR, and LMMSE beamformers W mvdr = 1 D 1 x(n) = V(k s )F s [n] + V(k l )F l [n] + Z[n] S 1 Vs H S 1 n V s V s W lmmse = p s S 1 x l=1 n V s, W mpdr = 1 S 1 Vs H S 1 x V s x V s and For Uncorrelated sources W mvdr = W mpdr, W lmmse = cw mpdr with c = p s (V H s S 1 x V s ) = c > 0. Output of MPDR beamformer is of the form r[n] = W H mpdr x[n] = F s [n] + q[n], where q[n] = W H mpdr I[n] Now we find a LMMSE estimate of F s [n] given r[n]. with c = p s (V H s S 1 x V s ) ˆF s [n] = c r[n]
5
6 Special Cases Case 1: σ 2 z dominates the interference W mvdr V s This is often referred to as coherent combining. Case 2: Z[n] dominates and E(Z[n]Z H [n]) = S z = diag(σ 2 0, σ2 1,..., σ2 N 1 ) W mvdr S 1 z V s Sensors weighted by the quality of the measurements they provide Case 3: (D 1) Interferers dominate Interferers are cancelled. r[n] = W H mvdr x[n] = F s [n]
7 Generalized Sidelobe Cancelers (GSC) Efficient Structure to implement a beamformer when there are constraints. Generalize the MPDR to include M constraints, where M N and g C M 1. C H W = g If C = V s and g = 1, we get the MPDR/MVDR distortionless constraint. Beamformer Design: Minimize beamformer output power subject to multiple constraints. min W W H S x W subject to C H W = g. Using Lagrange multipliers, we can show W o = S 1 x C(C H S 1 C) 1 g x
8 Singular Value Decomposition (SVD) Assumptions: C C N M, N M (constraints M fewer than sensors N) and rank(c) = M (columns of C are linearly independent) [ ] C = UΣV H S = [U 1, U 2 ] V H = U 0 1 SV H 1. U C N N and V C M M are orthonormal matrices, i.e. UU H = U H U = I N N and VV H = V H V = I M M 2. R(C) = R(U 1 ), R(C H ) = R(V ), and N (C H ) = R(U 2 ). 3. S = diag(σ) C M M and σ 1 σ 2.. σ M > 0. Since the beamformer W C N 1, W = U 1 α + U 2 β = U 1 α BW a = W q + W q where B is any orthonormal set with R(B) = R(U 2 ) and W a C (N M) 1 vector. Will show α C M 1 or W q can be pre-determined by the constraints and only β C (N M) 1 or Wq be determined from data. has to
9 GSC structure Computing W q C H W = VSU H 1 W = VSU H 1 (U 1 α + U 2 β) = VSα = g So α = S 1 V H g leading to W q = U 1 α = U 1 S 1 V H g = (C H ) + g Looking at the beamformer output r[n] = W H x[n] = W H q x[n] + (W q ) H x[n] = r q [n] W H a B H x[n] = r q [n] W H a r B [n] r q [n] and r B [n] can be computed from x[n] knowing the constraints, and hence can be assumed to be readily available. Problem is one of predicting r q [n] using r B [n]. Choose W a to minimize the mean squared error W a = Σ 1 r B r B Σ rb r q
10 GSC block diagram Computing W a LMMSE estimate W a = Σ 1 r B r B Σ rb r q. Can use an adaptive algorithm like the LMS or RLS
11 Receive versus Transmit Beamforming So far we have considered receive beamforming: Source transmits F s [n] and e jωc τ l F s [n] is received by sensor l Transmit Beamforming: How to transmit a signal s[n] from the array to a single sensor receiver? Imagine we transmit s[n] from sensor l. Signal received by receiver is e jωc τ l s[n] If all sensors transmit the same signal s[n], the transmitted signal is 1 1 x[n] =. s[n]. 1 The transmit beamformer is W = [1, 1,..., 1] T The signal received is by superposition principle N 1 l=0 τ l s[n] + z[n] = [1, 1,.., 1] V e jωc s s[n] + z[n] = }{{} xt [n]v s + z l [n] W T
12 Transmit BF If sensor l transmits w l s[n], the signal transmitted by the array is w 0 w 1 x[n] =. s[n]. w N 1 The transmit beamformer is W = [w 0, w 1,..., w N 1 ] T The signal received is, by superposition principle, x T [n]v s + z l [n] = W T V s s[n] + z[n] = N 1 l=0 w le jωc τ l s[n] + z[n] To maximize the signal received, we try to insure the components add up coherently, i.e. w l e jωc τ l Due to power constraints, we may have to divide the power among the sensors Summary: Transmmit signal is V s s[n] and transmit BF is V s
13 SNR gain due to Transmit BF Single sensor (at origin) transmitting: Received signal is s[n] + z[n]. SNR = σ2 s σ 2 z Assuming transmit BF is 1 N V s, the received signal is given by N 1 l=0 1 N e +jωc τ l e jωc τ l s[n] + z[n] = Ns[n] + z[n] SNR gain due to transmit BF is N. SNR BF = Nσ2 s σ 2 z How about transmitting to multiple users?
14 Transmitting to D users Challenge: Minimize interference, i.e. signal s l [n] meant for receiver l is not received by receiver k, k l. We need their location k l. Let us assume we know that. Transmitted signal D 1 x[n] = W l s l [n] = W s[n] l=0 where W = [W 0, W 1,.., W D 1 ] and s[n] = [s 0 [n], s 1 [n],..., s D 1 ] T. W l chosen to insure receiver l receives signal s l [n] and not by the others. Zero-Forcing Precoder W V (V T V ) 1 where V = [V 0, V 1,..., V D 1 ] Signal received by user l is r l [n] = x T [n]v l + z l [n] = s T [n]w T V l + z l [n] = s T [n](v H V) 1 V H V l + z l [n] = s T [n](v H V) 1 V H Ve l + z l [n] where e l = [0, }{{} 1, 0,..., 0] T position l = s l [n] + z l [n]
15 Power Spectrum Estimation of a Time Series Quality of Power Spectrum estimate given a time series of length N depends on Quality of bandpass filter (filter length L, assume FIR) Quality of power estimate, i.e. number of samples at the output of the bandpass filter.
16 MVDR Power Spectrum FIR filter of order L: H(z) = h 0 + h 1 z 1... h L 1 z (L 1) Impulse response vector h = [h 0, h 1,..., h L 1 ] T. Bandpass Filter defined by a distortionless constraint, i.e. H(e jω0 ) = 1. Here ω 0 is the frequency where the PSD is of interest H(e jωo ) = 1 or h H V(e jω0 ) = 1 where V(e jω0 ) = [1, e jω0,..., e j(l 1)ω0 ] T. Bandpass filter design: Minimize the filter output power E( y[n] 2 ), where L 1 y[n] = hl x[n l] = h H x[n], where x[n] = [x[n], x[n 1],..., x[n L+1]] T. l=1 Note that E( y[n] 2 ) = h H R xx h, where R xx = E(x[n]x H [n]) is a Toeplitz correlation matrix MVDR filter: h = 1 V H (e jω 0 )R 1 xx V(e jω 0 ) R 1 xx V(e jω0 ) MVDR Power spectrum: P(e jω0 ) = 1 V H (e jω 0 )R 1 xx V(e jω 0 )
17 Direction of Arrival (DOA) Estimation D 1 x[n] = V l F l [n] + Z[n], n = 0, 1,..., L 1 l=0 Here we have chosen to denote V s by V 0 for simplicity of notation. Note V l = V(k l ). We will refer to the L vector measurements as L snapshots. Goal: Estimate the direction of arrival k l or equivalently (θ l, φ l ), l = 0,..., D 1. Can also consider estimating the source powers p l, l = 0, 1,..., D 1, and noise variance σ 2 z. ULA : N 1 j V(ψ) = [e 2 ψ,..., e j N 1 2 ψ ] T N 1 j = e 2 ψ [1, e jψ, e j2ψ,..., e j(n 1)ψ ] T Note for a ULA, V(ψ) = JV (ψ), where J is all zeros except for ones along the anti-diagonal. Also note J = J 2.
18 ULA and Time Series x 0 [n] x 1 [n]. x N 1 [n] D 1 = l=0 1 e jψ l. e j(n 1)ψ l e j N 1 2 ψ l F l [n]+z[n], n = 0, 1,..., L 1 Considerable similarity with a time series problems composed of complex exponentials x[n] = D 1 l=0 c le jω l n + z[n], n = 0, 1,..., (M 1). In vector form x[0] 1 z[0] x[1] e jω l z[1]. x[m 1] D 1 = l=0. e j(m 1)ω l c l +. z[m 1] The time series, sum of exponentials, problem with M samples is equivalent to a ULA with M sensors and one (vector) measurement or one snapshot
19 Estimation of the Frequencies of the Exponential In array processing we have a limited number of sensors N but can compensate potentially with L snapshots. The time series problem is equivalent to one snapshot but a potentially large array if the number of samples M is large. One complex exponential (Noise free): x[n] = c 0 e jω0n, n = 0, 1,..., M 1. Interested in the frequency ω 0. Potential solution: Take the FFT and look at the peak. Another option is to note that x[n] = c 0 e jω0(n 1) e jω0 = e jω0 x[n 1]. In the noise free case, 2 samples is enough. Can solve a least squares problem min a M 1 n=1 x[n] ax[n 1] 2 and find ω 0 from the estimated a.
20 More than one exponentials Suppose we have two exponentials: x[n] = c 0 e jω0n + c 1 e jω1n, n = 0, 1,.., M 1. A FFT based approach may fail to resolve the frequencies if they are close and the number of samples M is small. Thumb rule ω 0 ω 1 > 2π M. Can show x[n] = a 1 x[n 1] + a 2 x[n 2] and so one can solve a 1 and a 2 using a least squares technique. In the absence of noise we just need 4 samples to set a system of equations to solve for the parameters. [ ] [ ] [ ] x[n] x[n 1] x[n 2] a1 = x[n 1] x[n 2] x[n 3] a 2 How to get the frequencies from the parameters a 1 and a 2? Can show that the zeros of the filter H(z) = 1 a 1 z 1 a 2 z 2 are located at e jω0 and e jω1.
21 D exponentials Suppose we have a series composed of D exponentials: x[n] = D 1 l=0 c le jω l n, n = 0, 1,.., M 1. A FFT based approach can be used with drawbacks as before: resolution and large number of samples may be needed. Can show x[n] = D k=1 a kx[n k] and so one can solve a k, k = 1,.., D using a least squares technique. In the absence of noise we just need 2D samples to set up a system of equations to solve for the parameters. x[n] x[n 1]... x[n D] a 1 x[n 1] x[n 2]... x[n D 1] a 2. x[n D + 1] =..... x[n D]... x[n 2D 1] How to get the frequencies from the parameters a k, k = 1,..D? Can show that the zeros of the filter H(z) = 1 D k=1 a kz k are located at e jω l, l = 0,.., D 1. a D
22 Relationship of Zeros of Filters to e jω l Claim: A signal with D exponentials, i.e. x[n] = D 1 l=0 c le jω l n, is exactly predictable from its D past samples, i.e. x[n] = D k=1 a kx[n k] or x[n] D k=1 a kx[n k] = 0 Consider a FIR filter H(z) = 1 P k=1 b kx[n k], P D, with input x[n] and denote its output by e[n] e[n] = x[n] P b k x[n k] = k=1 If we minimize n e[n] 2, what is the result? D H(e ω k )c k e jω k n The output will be zero. This is achieved by making H(e ω k ) = 0 at the D frequencies. Since the filter is of order P D, we can select the filter zeros to meet this criteria and minimize the output. k=1 The excess (P D) zeros can be placed arbitrarily.
Direction 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 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 informationMVDR and MPDR. Bhaskar D Rao University of California, San Diego
MVDR ad MPDR Bhaskar D Rao Uiversity of Califoria, Sa Diego Email: brao@ucsd.edu Referece Books. Optimum Array Processig, H. L. Va Trees 2. Stoica, P., & Moses, R. L. (2005). Spectral aalysis of sigals
More informationLinear Optimum Filtering: Statement
Ch2: Wiener Filters Optimal filters for stationary stochastic models are reviewed and derived in this presentation. Contents: Linear optimal filtering Principle of orthogonality Minimum mean squared error
More informationMicrophone-Array Signal Processing
Microphone-Array Signal Processing, c Apolinárioi & Campos p. 1/30 Microphone-Array Signal Processing José A. Apolinário Jr. and Marcello L. R. de Campos {apolin},{mcampos}@ieee.org IME Lab. Processamento
More informationBeamforming. A brief introduction. Brian D. Jeffs Associate Professor Dept. of Electrical and Computer Engineering Brigham Young University
Beamforming A brief introduction Brian D. Jeffs Associate Professor Dept. of Electrical and Computer Engineering Brigham Young University March 2008 References Barry D. Van Veen and Kevin Buckley, Beamforming:
More informationADAPTIVE ANTENNAS. SPATIAL BF
ADAPTIVE ANTENNAS SPATIAL BF 1 1-Spatial reference BF -Spatial reference beamforming may not use of embedded training sequences. Instead, the directions of arrival (DoA) of the impinging waves are used
More informationOverview of Beamforming
Overview of Beamforming Arye Nehorai Preston M. Green Department of Electrical and Systems Engineering Washington University in St. Louis March 14, 2012 CSSIP Lab 1 Outline Introduction Spatial and temporal
More informationCHAPTER 3 ROBUST ADAPTIVE BEAMFORMING
50 CHAPTER 3 ROBUST ADAPTIVE BEAMFORMING 3.1 INTRODUCTION Adaptive beamforming is used for enhancing a desired signal while suppressing noise and interference at the output of an array of sensors. It is
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 informationCo-Prime Arrays and Difference Set Analysis
7 5th European Signal Processing Conference (EUSIPCO Co-Prime Arrays and Difference Set Analysis Usham V. Dias and Seshan Srirangarajan Department of Electrical Engineering Bharti School of Telecommunication
More informationSpeaker Tracking and Beamforming
Speaker Tracking and Beamforming Dr. John McDonough Spoken Language Systems Saarland University January 13, 2010 Introduction Many problems in science and engineering can be formulated in terms of estimating
More informationIndependent Component Analysis (ICA) Bhaskar D Rao University of California, San Diego
Independent Component Analysis (ICA) Bhaskar D Rao University of California, San Diego Email: brao@ucsdedu References 1 Hyvarinen, A, Karhunen, J, & Oja, E (2004) Independent component analysis (Vol 46)
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 informationIII.C - Linear Transformations: Optimal Filtering
1 III.C - Linear Transformations: Optimal Filtering FIR Wiener Filter [p. 3] Mean square signal estimation principles [p. 4] Orthogonality principle [p. 7] FIR Wiener filtering concepts [p. 8] Filter coefficients
More informationMassachusetts Institute of Technology
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.011: Introduction to Communication, Control and Signal Processing QUIZ, April 1, 010 QUESTION BOOKLET Your
More informationA ROBUST BEAMFORMER BASED ON WEIGHTED SPARSE CONSTRAINT
Progress In Electromagnetics Research Letters, Vol. 16, 53 60, 2010 A ROBUST BEAMFORMER BASED ON WEIGHTED SPARSE CONSTRAINT Y. P. Liu and Q. Wan School of Electronic Engineering University of Electronic
More informationMicrophone-Array Signal Processing
Microphone-Array Signal Processing, c Apolinárioi & Campos p. 1/115 Microphone-Array Signal Processing José A. Apolinário Jr. and Marcello L. R. de Campos {apolin},{mcampos}@ieee.org IME Lab. Processamento
More information1. Determine if each of the following are valid autocorrelation matrices of WSS processes. (Correlation Matrix),R c =
ENEE630 ADSP Part II w/ solution. Determine if each of the following are valid autocorrelation matrices of WSS processes. (Correlation Matrix) R a = 4 4 4,R b = 0 0,R c = j 0 j 0 j 0 j 0 j,r d = 0 0 0
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 informationSample ECE275A Midterm Exam Questions
Sample ECE275A Midterm Exam Questions The questions given below are actual problems taken from exams given in in the past few years. Solutions to these problems will NOT be provided. These problems and
More informationProbability Space. J. McNames Portland State University ECE 538/638 Stochastic Signals Ver
Stochastic Signals Overview Definitions Second order statistics Stationarity and ergodicity Random signal variability Power spectral density Linear systems with stationary inputs Random signal memory Correlation
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 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 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 informationPolynomial Root-MUSIC Algorithm for Efficient Broadband Direction Of Arrival Estimation
Polynomial Root-MUSIC Algorithm for Efficient Broadband Direction Of Arrival Estimation William Coventry, Carmine Clemente, and John Soraghan University of Strathclyde, CESIP, EEE, 204, George Street,
More informationSingle-User MIMO systems: Introduction, capacity results, and MIMO beamforming
Single-User MIMO systems: Introduction, capacity results, and MIMO beamforming Master Universitario en Ingeniería de Telecomunicación I. Santamaría Universidad de Cantabria Contents Introduction Multiplexing,
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 informationDOA Estimation of Uncorrelated and Coherent Signals in Multipath Environment Using ULA Antennas
DOA Estimation of Uncorrelated and Coherent Signals in Multipath Environment Using ULA Antennas U.Somalatha 1 T.V.S.Gowtham Prasad 2 T. Ravi Kumar Naidu PG Student, Dept. of ECE, SVEC, Tirupati, Andhra
More informationAdaptive beamforming. Slide 2: Chapter 7: Adaptive array processing. Slide 3: Delay-and-sum. Slide 4: Delay-and-sum, continued
INF540 202 Adaptive beamforming p Adaptive beamforming Sven Peter Näsholm Department of Informatics, University of Oslo Spring semester 202 svenpn@ifiuiono Office phone number: +47 22840068 Slide 2: Chapter
More informationELEC E7210: Communication Theory. Lecture 10: MIMO systems
ELEC E7210: Communication Theory Lecture 10: MIMO systems Matrix Definitions, Operations, and Properties (1) NxM matrix a rectangular array of elements a A. an 11 1....... a a 1M. NM B D C E ermitian transpose
More informationEE538 Final Exam Fall 2007 Mon, Dec 10, 8-10 am RHPH 127 Dec. 10, Cover Sheet
EE538 Final Exam Fall 2007 Mon, Dec 10, 8-10 am RHPH 127 Dec. 10, 2007 Cover Sheet Test Duration: 120 minutes. Open Book but Closed Notes. Calculators allowed!! This test contains five problems. Each of
More informationPolynomial Matrix Formulation-Based Capon Beamformer
Alzin, Ahmed and Coutts, Fraser K. and Corr, Jamie and Weiss, Stephan and Proudler, Ian K. and Chambers, Jonathon A. (26) Polynomial matrix formulation-based Capon beamformer. In: th IMA International
More informationEEM 409. Random Signals. Problem Set-2: (Power Spectral Density, LTI Systems with Random Inputs) Problem 1: Problem 2:
EEM 409 Random Signals Problem Set-2: (Power Spectral Density, LTI Systems with Random Inputs) Problem 1: Consider a random process of the form = + Problem 2: X(t) = b cos(2π t + ), where b is a constant,
More informationTHE estimation of covariance matrices is a crucial component
1 A Subspace Method for Array Covariance Matrix Estimation Mostafa Rahmani and George K. Atia, Member, IEEE, arxiv:1411.0622v1 [cs.na] 20 Oct 2014 Abstract This paper introduces a subspace method for the
More informationSECTION FOR DIGITAL SIGNAL PROCESSING DEPARTMENT OF MATHEMATICAL MODELLING TECHNICAL UNIVERSITY OF DENMARK Course 04362 Digital Signal Processing: Solutions to Problems in Proakis and Manolakis, Digital
More informationConvolutive Transfer Function Generalized Sidelobe Canceler
IEEE TRANSACTIONS ON AUDIO, SPEECH AND LANGUAGE PROCESSING, VOL. XX, NO. Y, MONTH 8 Convolutive Transfer Function Generalized Sidelobe Canceler Ronen Talmon, Israel Cohen, Senior Member, IEEE, and Sharon
More informationarxiv: v1 [cs.ce] 3 Apr 2016
Two Dimensional Angle of Arrival Estimation arxiv:160400594v1 [csce] 3 Apr 2016 1 Introduction Santhosh Kumar, Pradip Sircar Department of Electrical Engineering Indian Institute of Technology Kanpur Kanpur
More informationBeamspace Adaptive Beamforming and the GSC
April 27, 2011 Overview The MVDR beamformer: performance and behavior. Generalized Sidelobe Canceller reformulation. Implementation of the GSC. Beamspace interpretation of the GSC. Reduced complexity of
More informationASYMPTOTIC PERFORMANCE ANALYSIS OF DOA ESTIMATION METHOD FOR AN INCOHERENTLY DISTRIBUTED SOURCE. 2πd λ. E[ϱ(θ, t)ϱ (θ,τ)] = γ(θ; µ)δ(θ θ )δ t,τ, (2)
ASYMPTOTIC PERFORMANCE ANALYSIS OF DOA ESTIMATION METHOD FOR AN INCOHERENTLY DISTRIBUTED SOURCE Jooshik Lee and Doo Whan Sang LG Electronics, Inc. Seoul, Korea Jingon Joung School of EECS, KAIST Daejeon,
More informationELEN E4810: Digital Signal Processing Topic 11: Continuous Signals. 1. Sampling and Reconstruction 2. Quantization
ELEN E4810: Digital Signal Processing Topic 11: Continuous Signals 1. Sampling and Reconstruction 2. Quantization 1 1. Sampling & Reconstruction DSP must interact with an analog world: A to D D to A x(t)
More informationEs e j4φ +4N n. 16 KE s /N 0. σ 2ˆφ4 1 γ s. p(φ e )= exp 1 ( 2πσ φ b cos N 2 φ e 0
Problem 6.15 : he received signal-plus-noise vector at the output of the matched filter may be represented as (see (5-2-63) for example) : r n = E s e j(θn φ) + N n where θ n =0,π/2,π,3π/2 for QPSK, and
More informationDETECTION theory deals primarily with techniques for
ADVANCED SIGNAL PROCESSING SE Optimum Detection of Deterministic and Random Signals Stefan Tertinek Graz University of Technology turtle@sbox.tugraz.at Abstract This paper introduces various methods for
More informationA SEMI-BLIND TECHNIQUE FOR MIMO CHANNEL MATRIX ESTIMATION. AdityaKiran Jagannatham and Bhaskar D. Rao
A SEMI-BLIND TECHNIQUE FOR MIMO CHANNEL MATRIX ESTIMATION AdityaKiran Jagannatham and Bhaskar D. Rao Department of Electrical and Computer Engineering University of California, San Diego La Jolla, CA 9093-0407
More informationEE538 Final Exam Fall :20 pm -5:20 pm PHYS 223 Dec. 17, Cover Sheet
EE538 Final Exam Fall 005 3:0 pm -5:0 pm PHYS 3 Dec. 17, 005 Cover Sheet Test Duration: 10 minutes. Open Book but Closed Notes. Calculators ARE allowed!! This test contains five problems. Each of the five
More informationAn Adaptive Sensor Array Using an Affine Combination of Two Filters
An Adaptive Sensor Array Using an Affine Combination of Two Filters Tõnu Trump Tallinn University of Technology Department of Radio and Telecommunication Engineering Ehitajate tee 5, 19086 Tallinn Estonia
More informationSolutions. Number of Problems: 10
Final Exam February 4th, 01 Signals & Systems (151-0575-01) Prof. R. D Andrea Solutions Exam Duration: 150 minutes Number of Problems: 10 Permitted aids: One double-sided A4 sheet. Questions can be answered
More information1.1.3 The narrowband Uniform Linear Array (ULA) with d = λ/2:
Seminar 1: Signal Processing Antennas 4ED024, Sven Nordebo 1.1.3 The narrowband Uniform Linear Array (ULA) with d = λ/2: d Array response vector: a() = e e 1 jπ sin. j(π sin )(M 1) = 1 e jω. e jω(m 1)
More informationarxiv: v1 [cs.it] 6 Nov 2016
UNIT CIRCLE MVDR BEAMFORMER Saurav R. Tuladhar, John R. Buck University of Massachusetts Dartmouth Electrical and Computer Engineering Department North Dartmouth, Massachusetts, USA arxiv:6.272v [cs.it]
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 informationThe Probability Distribution of the MVDR Beamformer Outputs under Diagonal Loading. N. Raj Rao (Dept. of Electrical Engineering and Computer Science)
The Probability Distribution of the MVDR Beamformer Outputs under Diagonal Loading N. Raj Rao (Dept. of Electrical Engineering and Computer Science) & Alan Edelman (Dept. of Mathematics) Work supported
More informationADAPTIVE FILTER THEORY
ADAPTIVE FILTER THEORY Fifth Edition Simon Haykin Communications Research Laboratory McMaster University Hamilton, Ontario, Canada International Edition contributions by Telagarapu Prabhakar Department
More informationReview of Fundamentals of Digital Signal Processing
Chapter 2 Review of Fundamentals of Digital Signal Processing 2.1 (a) This system is not linear (the constant term makes it non linear) but is shift-invariant (b) This system is linear but not shift-invariant
More informationNotes on Linear Minimum Mean Square Error Estimators
Notes on Linear Minimum Mean Square Error Estimators Ça gatay Candan January, 0 Abstract Some connections between linear minimum mean square error estimators, maximum output SNR filters and the least square
More informationQuantization and Compensation in Sampled Interleaved Multi-Channel Systems
Quantization and Compensation in Sampled Interleaved Multi-Channel Systems The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters. Citation
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 informationx[n] = x a (nt ) x a (t)e jωt dt while the discrete time signal x[n] has the discrete-time Fourier transform x[n]e jωn
Sampling Let x a (t) be a continuous time signal. The signal is sampled by taking the signal value at intervals of time T to get The signal x(t) has a Fourier transform x[n] = x a (nt ) X a (Ω) = x a (t)e
More informationThe Selection of Weighting Functions For Linear Arrays Using Different Techniques
The Selection of Weighting Functions For Linear Arrays Using Different Techniques Dr. S. Ravishankar H. V. Kumaraswamy B. D. Satish 17 Apr 2007 Rajarshi Shahu College of Engineering National Conference
More informationLeast Square Es?ma?on, Filtering, and Predic?on: ECE 5/639 Sta?s?cal Signal Processing II: Linear Es?ma?on
Least Square Es?ma?on, Filtering, and Predic?on: Sta?s?cal Signal Processing II: Linear Es?ma?on Eric Wan, Ph.D. Fall 2015 1 Mo?va?ons If the second-order sta?s?cs are known, the op?mum es?mator is given
More informationUSING STATISTICAL ROOM ACOUSTICS FOR ANALYSING THE OUTPUT SNR OF THE MWF IN ACOUSTIC SENSOR NETWORKS. Toby Christian Lawin-Ore, Simon Doclo
th European Signal Processing Conference (EUSIPCO 1 Bucharest, Romania, August 7-31, 1 USING STATISTICAL ROOM ACOUSTICS FOR ANALYSING THE OUTPUT SNR OF THE MWF IN ACOUSTIC SENSOR NETWORKS Toby Christian
More informationThe goal of the Wiener filter is to filter out noise that has corrupted a signal. It is based on a statistical approach.
Wiener filter From Wikipedia, the free encyclopedia In signal processing, the Wiener filter is a filter proposed by Norbert Wiener during the 1940s and published in 1949. [1] Its purpose is to reduce the
More informationEE 5407 Part II: Spatial Based Wireless Communications
EE 5407 Part II: Spatial Based Wireless Communications Instructor: Prof. Rui Zhang E-mail: rzhang@i2r.a-star.edu.sg Website: http://www.ece.nus.edu.sg/stfpage/elezhang/ Lecture II: Receive Beamforming
More informationAcoustic MIMO Signal Processing
Yiteng Huang Jacob Benesty Jingdong Chen Acoustic MIMO Signal Processing With 71 Figures Ö Springer Contents 1 Introduction 1 1.1 Acoustic MIMO Signal Processing 1 1.2 Organization of the Book 4 Part I
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 informationProf. Dr.-Ing. Armin Dekorsy Department of Communications Engineering. Stochastic Processes and Linear Algebra Recap Slides
Prof. Dr.-Ing. Armin Dekorsy Department of Communications Engineering Stochastic Processes and Linear Algebra Recap Slides Stochastic processes and variables XX tt 0 = XX xx nn (tt) xx 2 (tt) XX tt XX
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 informationCoprime Coarray Interpolation for DOA Estimation via Nuclear Norm Minimization
Coprime Coarray Interpolation for DOA Estimation via Nuclear Norm Minimization Chun-Lin Liu 1 P. P. Vaidyanathan 2 Piya Pal 3 1,2 Dept. of Electrical Engineering, MC 136-93 California Institute of Technology,
More informationReduced-Rank Multi-Antenna Cyclic Wiener Filtering for Interference Cancellation
Reduced-Rank Multi-Antenna Cyclic Wiener Filtering for Interference Cancellation Hong Zhang, Ali Abdi and Alexander Haimovich Center for Wireless Communications and Signal Processing Research Department
More information1 Linearity and Linear Systems
Mathematical Tools for Neuroscience (NEU 34) Princeton University, Spring 26 Jonathan Pillow Lecture 7-8 notes: Linear systems & SVD Linearity and Linear Systems Linear system is a kind of mapping f( x)
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 informationChapter 2 Wiener Filtering
Chapter 2 Wiener Filtering Abstract Before moving to the actual adaptive filtering problem, we need to solve the optimum linear filtering problem (particularly, in the mean-square-error sense). We start
More informationNear Optimal Adaptive Robust Beamforming
Near Optimal Adaptive Robust Beamforming The performance degradation in traditional adaptive beamformers can be attributed to the imprecise knowledge of the array steering vector and inaccurate estimation
More informationEE482: Digital Signal Processing Applications
Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu EE482: Digital Signal Processing Applications Spring 2014 TTh 14:30-15:45 CBC C222 Lecture 11 Adaptive Filtering 14/03/04 http://www.ee.unlv.edu/~b1morris/ee482/
More informationROBUST ADAPTIVE BEAMFORMING BASED ON CO- VARIANCE MATRIX RECONSTRUCTION FOR LOOK DIRECTION MISMATCH
Progress In Electromagnetics Research Letters, Vol. 25, 37 46, 2011 ROBUST ADAPTIVE BEAMFORMING BASED ON CO- VARIANCE MATRIX RECONSTRUCTION FOR LOOK DIRECTION MISMATCH R. Mallipeddi 1, J. P. Lie 2, S.
More informationVery useful for designing and analyzing signal processing systems
z-transform z-transform The z-transform generalizes the Discrete-Time Fourier Transform (DTFT) for analyzing infinite-length signals and systems Very useful for designing and analyzing signal processing
More informationEE 225D LECTURE ON DIGITAL FILTERS. University of California Berkeley
University of California Berkeley College of Engineering Department of Electrical Engineering and Computer Sciences Professors : N.Morgan / B.Gold EE225D Digital Filters Spring,1999 Lecture 7 N.MORGAN
More informationArray Signal Processing Algorithms for Beamforming and Direction Finding
Array Signal Processing Algorithms for Beamforming and Direction Finding This thesis is submitted in partial fulfilment of the requirements for Doctor of Philosophy (Ph.D.) Lei Wang Communications Research
More informationCS281 Section 4: Factor Analysis and PCA
CS81 Section 4: Factor Analysis and PCA Scott Linderman At this point we have seen a variety of machine learning models, with a particular emphasis on models for supervised learning. In particular, we
More informationNew Cramér-Rao Bound Expressions for Coprime and Other Sparse Arrays
New Cramér-Rao Bound Expressions for Coprime and Other Sparse Arrays Chun-Lin Liu 1 P. P. Vaidyanathan 2 1,2 Dept. of Electrical Engineering, MC 136-93 California Institute of Technology, Pasadena, CA
More informationNotes on singular value decomposition for Math 54. Recall that if A is a symmetric n n matrix, then A has real eigenvalues A = P DP 1 A = P DP T.
Notes on singular value decomposition for Math 54 Recall that if A is a symmetric n n matrix, then A has real eigenvalues λ 1,, λ n (possibly repeated), and R n has an orthonormal basis v 1,, v n, where
More informationEE482: Digital Signal Processing Applications
Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu EE482: Digital Signal Processing Applications Spring 2014 TTh 14:30-15:45 CBC C222 Lecture 11 Adaptive Filtering 14/03/04 http://www.ee.unlv.edu/~b1morris/ee482/
More informationAdvanced Digital Signal Processing -Introduction
Advanced Digital Signal Processing -Introduction LECTURE-2 1 AP9211- ADVANCED DIGITAL SIGNAL PROCESSING UNIT I DISCRETE RANDOM SIGNAL PROCESSING Discrete Random Processes- Ensemble Averages, Stationary
More informationUnderdetermined DOA Estimation Using MVDR-Weighted LASSO
Article Underdetermined DOA Estimation Using MVDR-Weighted LASSO Amgad A. Salama, M. Omair Ahmad and M. N. S. Swamy Department of Electrical and Computer Engineering, Concordia University, Montreal, PQ
More informationParametric Signal Modeling and Linear Prediction Theory 1. Discrete-time Stochastic Processes (cont d)
Parametric Signal Modeling and Linear Prediction Theory 1. Discrete-time Stochastic Processes (cont d) Electrical & Computer Engineering North Carolina State University Acknowledgment: ECE792-41 slides
More information12.4 Known Channel (Water-Filling Solution)
ECEn 665: Antennas and Propagation for Wireless Communications 54 2.4 Known Channel (Water-Filling Solution) The channel scenarios we have looed at above represent special cases for which the capacity
More informationCast of Characters. Some Symbols, Functions, and Variables Used in the Book
Page 1 of 6 Cast of Characters Some s, Functions, and Variables Used in the Book Digital Signal Processing and the Microcontroller by Dale Grover and John R. Deller ISBN 0-13-081348-6 Prentice Hall, 1998
More informationZ subarray. (d,0) (Nd-d,0) (Nd,0) X subarray Y subarray
A Fast Algorithm for 2-D Direction-of-Arrival Estimation Yuntao Wu 1,Guisheng Liao 1 and H. C. So 2 1 Laboratory for Radar Signal Processing, Xidian University, Xian, China 2 Department of Computer Engineering
More informationPerformance Analysis for Strong Interference Remove of Fast Moving Target in Linear Array Antenna
Performance Analysis for Strong Interference Remove of Fast Moving Target in Linear Array Antenna Kwan Hyeong Lee Dept. Electriacal Electronic & Communicaton, Daejin University, 1007 Ho Guk ro, Pochen,Gyeonggi,
More informationMITIGATING UNCORRELATED PERIODIC DISTURBANCE IN NARROWBAND ACTIVE NOISE CONTROL SYSTEMS
17th European Signal Processing Conference (EUSIPCO 29) Glasgow, Scotland, August 24-28, 29 MITIGATING UNCORRELATED PERIODIC DISTURBANCE IN NARROWBAND ACTIVE NOISE CONTROL SYSTEMS Muhammad Tahir AKHTAR
More informationRadar Systems Engineering Lecture 3 Review of Signals, Systems and Digital Signal Processing
Radar Systems Engineering Lecture Review of Signals, Systems and Digital Signal Processing Dr. Robert M. O Donnell Guest Lecturer Radar Systems Course Review Signals, Systems & DSP // Block Diagram of
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 informationCCNY. BME I5100: Biomedical Signal Processing. Stochastic Processes. Lucas C. Parra Biomedical Engineering Department City College of New York
BME I5100: Biomedical Signal Processing Stochastic Processes Lucas C. Parra Biomedical Engineering Department CCNY 1 Schedule Week 1: Introduction Linear, stationary, normal - the stuff biology is not
More informationRobust Adaptive Beamforming Based on Low-Complexity Shrinkage-Based Mismatch Estimation
1 Robust Adaptive Beamforming Based on Low-Complexity Shrinkage-Based Mismatch Estimation Hang Ruan and Rodrigo C. de Lamare arxiv:1311.2331v1 [cs.it] 11 Nov 213 Abstract In this work, we propose a low-complexity
More informationHIGH RESOLUTION DOA ESTIMATION IN FULLY COHERENT ENVIRONMENTS. S. N. Shahi, M. Emadi, and K. Sadeghi Sharif University of Technology Iran
Progress In Electromagnetics Research C, Vol. 5, 35 48, 28 HIGH RESOLUTION DOA ESTIMATION IN FULLY COHERENT ENVIRONMENTS S. N. Shahi, M. Emadi, and K. Sadeghi Sharif University of Technology Iran Abstract
More informationX. Cross Spectral Analysis
X. Cross Spectral Analysis Cross Spectrum We have already dealt with the crosscovariance function (ccvf): ˆ 1 C (k) = X(i)Y(i + k) N N k i= 0 The Fourier transform of the ccvf is called the cross spectrum.
More informationExamination with solution suggestions SSY130 Applied Signal Processing
Examination with solution suggestions SSY3 Applied Signal Processing Jan 8, 28 Rules Allowed aids at exam: L. Råde and B. Westergren, Mathematics Handbook (any edition, including the old editions called
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 informationDigital Image Processing Lectures 25 & 26
Lectures 25 & 26, Professor Department of Electrical and Computer Engineering Colorado State University Spring 2015 Area 4: Image Encoding and Compression Goal: To exploit the redundancies in the image
More information26. Filtering. ECE 830, Spring 2014
26. Filtering ECE 830, Spring 2014 1 / 26 Wiener Filtering Wiener filtering is the application of LMMSE estimation to recovery of a signal in additive noise under wide sense sationarity assumptions. Problem
More informationSignal Design for Band-Limited Channels
Wireless Information Transmission System Lab. Signal Design for Band-Limited Channels Institute of Communications Engineering National Sun Yat-sen University Introduction We consider the problem of signal
More information