Tim Habigt. May 30, 2014

Size: px
Start display at page:

Download "Tim Habigt. May 30, 2014"

Transcription

1 Slide 1/26 Sound-Source Localization Tim Habigt May 2014 Sound-Source Localization Tim Habigt Lehrstuhl für Datenverarbeitung May 30, 2014

2 Slide 2/26 Sound-Source Localization Tim Habigt May 2014 Table of contents Sound source localization Microphone arrays / binaural robots Time delay estimation Beamforming Binaural localization

3 Slide 3/26 Sound-Source Localization Tim Habigt May 2014 Microphone arrays? (a) Array? Figure: Acoustic locators (around 1920) (b) Vertical localization

4 Slide 4/26 Sound-Source Localization Tim Habigt May 2014 Microphone arrays (a) Microsoft Kinect Figure: Commercial arrays (b) Playstation Eye

5 Microphone arrays (a) CCRL array (b) MIT array with 1020 microphones Figure: Larger arrays Slide 5/26 Sound-Source Localization Tim Habigt May 2014

6 Slide 6/26 Sound-Source Localization Tim Habigt May 2014 Binaural localization Figure: Binaural localization

7 Videos Slide 7/26 Sound-Source Localization Tim Habigt May 2014

8 Slide 8/26 Sound-Source Localization Tim Habigt May 2014 Signal propagation Figure: Coordinate systems 1 1 courtesy of Marko Durkovic

9 Slide 9/26 Sound-Source Localization Tim Habigt May 2014 Signal propagation Figure: Time delay 2 2 courtesy of Marko Durkovic

10 Slide 10/26 Sound-Source Localization Tim Habigt May 2014 Signal propagation x 0 (t) = x 1 (t) = s(t) s(t t) t = d c sin(α) Speed of sound c = 343 m s

11 Slide 10/26 Sound-Source Localization Tim Habigt May 2014 Signal propagation x 0 (t) = x 1 (t) = s(t) s(t t) t = d c sin(α) Speed of sound c = 343 m s

12 Slide 10/26 Sound-Source Localization Tim Habigt May 2014 Signal propagation x 0 (t) = x 1 (t) = s(t) s(t t) t = d c sin(α) Speed of sound c = 343 m s

13 Slide 11/26 Sound-Source Localization Tim Habigt May 2014 Localization techniques Time difference of arrival (TDOA) Steered response power (beamforming) HRTF-based binaural localization

14 Slide 12/26 Sound-Source Localization Tim Habigt May 2014 TDOA-based localization Cross-correlation Auto-correlation R xy [m] = R xx [m] = n= n= m denotes the shift denotes complex conjugation x[n]y [n + m] (1) x[n]x [n + m] (2)

15 Slide 13/26 Sound-Source Localization Tim Habigt May 2014 TDOA-based localization Cross-correlation demo

16 Slide 14/26 Sound-Source Localization Tim Habigt May 2014 TDOA-based localization Figure: Time delay 3 Drawbacks: Binaural localization via time delays is ambiguous ( cone of confusion ) 3 courtesy of Marko Durkovic

17 Slide 15/26 Sound-Source Localization Tim Habigt May 2014 Beamforming Figure: Beam pattern 4 4 courtesy of Dr. Andrew Greensted

18 Slide 16/26 Sound-Source Localization Tim Habigt May 2014 Beamforming Figure: Signal summation 5 5 courtesy of Dr. Andrew Greensted

19 Slide 17/26 Sound-Source Localization Tim Habigt May 2014 Beamforming Signal summation demo

20 Slide 18/26 Sound-Source Localization Tim Habigt May 2014 Beamforming Figure: Beamforming 6 6 courtesy of Dr. Andrew Greensted

21 Slide 19/26 Sound-Source Localization Tim Habigt May 2014 Time shift x(t t) = x(ω)e jω t j2πf t = x(ω)e Frequency domain processing demo

22 Slide 20/26 Sound-Source Localization Tim Habigt May 2014 Multiple microphones x 0 (t) = x 1 (t) = x 2 (t) = s(t) s(t t) s(t 2 t) x 0 (ω) = 1 s(ω) x 1 (ω) = e jω t s(ω) x 2 (ω) = e jω2 t s(ω) 1 x(ω) = e jω t s(ω) e jω2 t x(ω) = a(α d )s(ω) a is called steering vector

23 Slide 20/26 Sound-Source Localization Tim Habigt May 2014 Multiple microphones x 0 (t) = x 1 (t) = x 2 (t) = s(t) s(t t) s(t 2 t) x 0 (ω) = 1 s(ω) x 1 (ω) = e jω t s(ω) x 2 (ω) = e jω2 t s(ω) 1 x(ω) = e jω t s(ω) e jω2 t x(ω) = a(α d )s(ω) a is called steering vector

24 Slide 20/26 Sound-Source Localization Tim Habigt May 2014 Multiple microphones x 0 (t) = x 1 (t) = x 2 (t) = s(t) s(t t) s(t 2 t) x 0 (ω) = 1 s(ω) x 1 (ω) = e jω t s(ω) x 2 (ω) = e jω2 t s(ω) 1 x(ω) = e jω t s(ω) e jω2 t x(ω) = a(α d )s(ω) a is called steering vector

25 Slide 20/26 Sound-Source Localization Tim Habigt May 2014 Multiple microphones x 0 (t) = x 1 (t) = x 2 (t) = s(t) s(t t) s(t 2 t) x 0 (ω) = 1 s(ω) x 1 (ω) = e jω t s(ω) x 2 (ω) = e jω2 t s(ω) 1 x(ω) = e jω t s(ω) e jω2 t x(ω) = a(α d )s(ω) a is called steering vector

26 Slide 21/26 Sound-Source Localization Tim Habigt May 2014 Beamforming Figure: Beamformer

27 Slide 22/26 Sound-Source Localization Tim Habigt May 2014 Delay and Sum Beamforming y[n] = w H x[n] To align all the microphone signals choose w = a(θ d )

28 Slide 23/26 Sound-Source Localization Tim Habigt May 2014 Delay and Sum Beamforming Delay and Sum Beamforming demo

29 Slide 24/26 Sound-Source Localization Tim Habigt May 2014 Delay and Sum Beamforming Figure: Spatial aliasing No aliasing if d c f

30 Binaural localization Head-Related Transfer Function (HRTF) HRTF describes influence of head and pinna on the incoming sound wave Spatial dependency Slide 25/26 Sound-Source Localization Tim Habigt May 2014

31 Slide 26/26 Sound-Source Localization Tim Habigt May 2014 Cross-convolution localization Convolve microphone signal with all pairs of HRTFs (left and right ear switched) Correct location η maximizes similarity

COMP 546. Lecture 20. Head and Ear. Thurs. March 29, 2018

COMP 546. Lecture 20. Head and Ear. Thurs. March 29, 2018 COMP 546 Lecture 20 Head and Ear Thurs. March 29, 2018 1 Impulse function at t = 0. I X, Y, Z, t = δ(x X 0, Y Y 0, Z Z 0, t) To define an impulse function properly in a continuous space requires more math.

More information

Conventional beamforming

Conventional beamforming INF5410 2012. Conventional beamforming p.1 Conventional beamforming Slide 2: Beamforming Sven Peter Näsholm Department of Informatics, University of Oslo Spring semester 2012 svenpn@ifi.uio.no Office telephone

More information

COMP 546. Lecture 21. Cochlea to brain, Source Localization. Tues. April 3, 2018

COMP 546. Lecture 21. Cochlea to brain, Source Localization. Tues. April 3, 2018 COMP 546 Lecture 21 Cochlea to brain, Source Localization Tues. April 3, 2018 1 Ear pinna auditory canal cochlea outer middle inner 2 Eye Ear Lens? Retina? Photoreceptors (light -> chemical) Ganglion cells

More information

Spatial sound. Lecture 8: EE E6820: Speech & Audio Processing & Recognition. Columbia University Dept. of Electrical Engineering

Spatial sound. Lecture 8: EE E6820: Speech & Audio Processing & Recognition. Columbia University Dept. of Electrical Engineering EE E6820: Speech & Audio Processing & Recognition Lecture 8: Spatial sound 1 Spatial acoustics 2 Binaural perception 3 Synthesizing spatial audio 4 Extracting spatial sounds Dan Ellis

More information

Delay-and-Sum Beamforming for Plane Waves

Delay-and-Sum Beamforming for Plane Waves Delay-and-Sum Beamforming for Plane Waves ECE 6279: Spatial Array Processing Spring 2011 Lecture 6 Prof. Aaron D. Lanterman School of Electrical & Computer Engineering Georgia Institute of Technology AL:

More information

Sound impulse. r = (X X 0 ) 2 + (Y Y 0 ) 2 + (Z Z 0 ) 2. and

Sound impulse. r = (X X 0 ) 2 + (Y Y 0 ) 2 + (Z Z 0 ) 2. and Sound impulse Consider an isolated perturbation of air pressure at 3D point (X o, Y o, Z o ) and at time t = t 0, for example, due to some impact. Or, you can imagine a digital sound generation system

More information

APPLICATION OF MVDR BEAMFORMING TO SPHERICAL ARRAYS

APPLICATION OF MVDR BEAMFORMING TO SPHERICAL ARRAYS AMBISONICS SYMPOSIUM 29 June 2-27, Graz APPLICATION OF MVDR BEAMFORMING TO SPHERICAL ARRAYS Anton Schlesinger 1, Marinus M. Boone 2 1 University of Technology Delft, The Netherlands (a.schlesinger@tudelft.nl)

More information

6.003 (Fall 2011) Quiz #3 November 16, 2011

6.003 (Fall 2011) Quiz #3 November 16, 2011 6.003 (Fall 2011) Quiz #3 November 16, 2011 Name: Kerberos Username: Please circle your section number: Section Time 2 11 am 3 1 pm 4 2 pm Grades will be determined by the correctness of your answers (explanations

More information

Review of Discrete-Time System

Review of Discrete-Time System Review of Discrete-Time System Electrical & Computer Engineering University of Maryland, College Park Acknowledgment: ENEE630 slides were based on class notes developed by Profs. K.J. Ray Liu and Min Wu.

More information

A LOCALIZATION METHOD FOR MULTIPLE SOUND SOURCES BY USING COHERENCE FUNCTION

A LOCALIZATION METHOD FOR MULTIPLE SOUND SOURCES BY USING COHERENCE FUNCTION 8th European Signal Processing Conference (EUSIPCO-2) Aalborg, Denmark, August 23-27, 2 A LOCALIZATION METHOD FOR MULTIPLE SOUND SOURCES BY USING COHERENCE FUNCTION Hiromichi NAKASHIMA, Mitsuru KAWAMOTO,

More information

Albenzio Cirillo INFOCOM Dpt. Università degli Studi di Roma, Sapienza

Albenzio Cirillo INFOCOM Dpt. Università degli Studi di Roma, Sapienza Albenzio Cirillo INFOCOM Dpt. Università degli Studi di Roma, Sapienza albenzio.cirillo@uniroma1.it http://ispac.ing.uniroma1.it/albenzio/index.htm ET2010 XXVI Riunione Annuale dei Ricercatori di Elettrotecnica

More information

Module 3. Convolution. Aim

Module 3. Convolution. Aim Module Convolution Digital Signal Processing. Slide 4. Aim How to perform convolution in real-time systems efficiently? Is convolution in time domain equivalent to multiplication of the transformed sequence?

More information

Speaker Tracking and Beamforming

Speaker 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 information

Introduction to Audio and Music Engineering

Introduction to Audio and Music Engineering Introduction to Audio and Music Engineering Lecture 7 Sound waves Sound localization Sound pressure level Range of human hearing Sound intensity and power 3 Waves in Space and Time Period: T Seconds Frequency:

More information

Properties of LTI Systems

Properties of LTI Systems Properties of LTI Systems Properties of Continuous Time LTI Systems Systems with or without memory: A system is memory less if its output at any time depends only on the value of the input at that same

More information

ON THE LIMITATIONS OF BINAURAL REPRODUCTION OF MONAURAL BLIND SOURCE SEPARATION OUTPUT SIGNALS

ON THE LIMITATIONS OF BINAURAL REPRODUCTION OF MONAURAL BLIND SOURCE SEPARATION OUTPUT SIGNALS th European Signal Processing Conference (EUSIPCO 12) Bucharest, Romania, August 27-31, 12 ON THE LIMITATIONS OF BINAURAL REPRODUCTION OF MONAURAL BLIND SOURCE SEPARATION OUTPUT SIGNALS Klaus Reindl, Walter

More information

Signals & Systems Handout #4

Signals & Systems Handout #4 Signals & Systems Handout #4 H-4. Elementary Discrete-Domain Functions (Sequences): Discrete-domain functions are defined for n Z. H-4.. Sequence Notation: We use the following notation to indicate the

More information

A Probability Model for Interaural Phase Difference

A Probability Model for Interaural Phase Difference A Probability Model for Interaural Phase Difference Michael I. Mandel, Daniel P.W. Ellis Department of Electrical Engineering Columbia University, New York, New York {mim,dpwe}@ee.columbia.edu Abstract

More information

FUNDAMENTAL LIMITATION OF FREQUENCY DOMAIN BLIND SOURCE SEPARATION FOR CONVOLVED MIXTURE OF SPEECH

FUNDAMENTAL LIMITATION OF FREQUENCY DOMAIN BLIND SOURCE SEPARATION FOR CONVOLVED MIXTURE OF SPEECH FUNDAMENTAL LIMITATION OF FREQUENCY DOMAIN BLIND SOURCE SEPARATION FOR CONVOLVED MIXTURE OF SPEECH Shoko Araki Shoji Makino Ryo Mukai Tsuyoki Nishikawa Hiroshi Saruwatari NTT Communication Science Laboratories

More information

NOISE 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 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 information

Continuous Fourier transform of a Gaussian Function

Continuous Fourier transform of a Gaussian Function Continuous Fourier transform of a Gaussian Function Gaussian function: e t2 /(2σ 2 ) The CFT of a Gaussian function is also a Gaussian function (i.e., time domain is Gaussian, then the frequency domain

More information

Headphone Auralization of Acoustic Spaces Recorded with Spherical Microphone Arrays. Carl Andersson

Headphone Auralization of Acoustic Spaces Recorded with Spherical Microphone Arrays. Carl Andersson Headphone Auralization of Acoustic Spaces Recorded with Spherical Microphone Arrays Carl Andersson Department of Civil Engineering Chalmers University of Technology Gothenburg, 2017 Master s Thesis BOMX60-16-03

More information

ECE 598: The Speech Chain. Lecture 5: Room Acoustics; Filters

ECE 598: The Speech Chain. Lecture 5: Room Acoustics; Filters ECE 598: The Speech Chain Lecture 5: Room Acoustics; Filters Today Room = A Source of Echoes Echo = Delayed, Scaled Copy Addition and Subtraction of Scaled Cosines Frequency Response Impulse Response Filter

More information

ANALYSIS OF A PURINA FRACTAL BEAMFORMER. P. Karagiannakis 1, and S. Weiss 1

ANALYSIS OF A PURINA FRACTAL BEAMFORMER. P. Karagiannakis 1, and S. Weiss 1 ANALYSIS OF A PURINA FRACTAL BEAMFORMER P Karagiannakis 1, and S Weiss 1 1 Department of Electronics & Electrical Engineering University of Strathclyde, Glasgow G1 1XW, Scotland {philippkaragiannakis,stephanweiss}@strathacuk

More information

Grades will be determined by the correctness of your answers (explanations are not required).

Grades will be determined by the correctness of your answers (explanations are not required). 6.00 (Fall 20) Final Examination December 9, 20 Name: Kerberos Username: Please circle your section number: Section Time 2 am pm 4 2 pm Grades will be determined by the correctness of your answers (explanations

More information

Beamforming. 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 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 information

Acoustic Source Separation with Microphone Arrays CCNY

Acoustic Source Separation with Microphone Arrays CCNY Acoustic Source Separation with Microphone Arrays Lucas C. Parra Biomedical Engineering Department City College of New York CCNY Craig Fancourt Clay Spence Chris Alvino Montreal Workshop, Nov 6, 2004 Blind

More information

Convention Paper Presented at the 125th Convention 2008 October 2 5 San Francisco, CA, USA

Convention Paper Presented at the 125th Convention 2008 October 2 5 San Francisco, CA, USA Audio Engineering Society Convention Paper Presented at the 125th Convention 2008 October 2 5 San Francisco, CA, USA The papers at this Convention have been selected on the basis of a submitted abstract

More information

Due dates are as mentioned above. Checkoff interviews for PS4 and PS5 will happen together between October 24 and 28, 2012.

Due dates are as mentioned above. Checkoff interviews for PS4 and PS5 will happen together between October 24 and 28, 2012. Problem Set 5 Your answers will be graded by actual human beings (at least that's what we believe!), so don't limit your answers to machine-gradable responses. Some of the questions specifically ask for

More information

ERRATA Discrete-Time Signal Processing, 3e A. V. Oppenheim and R. W. Schafer

ERRATA Discrete-Time Signal Processing, 3e A. V. Oppenheim and R. W. Schafer ERRATA Discrete-Time Signal Processing, 3e A. V. Oppenheim and R. W. Schafer The following were corrected in the second printing; i.e., they are errors found only in the first printing. p.181 In the third

More information

Adaptive beamforming. Slide 2: Chapter 7: Adaptive array processing. Slide 3: Delay-and-sum. Slide 4: Delay-and-sum, continued

Adaptive 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 information

Predicting speech intelligibility in noisy rooms.

Predicting speech intelligibility in noisy rooms. Acknowledgement: Work supported by UK EPSRC Predicting speech intelligibility in noisy rooms. John F. Culling 1, Mathieu Lavandier 2 and Sam Jelfs 3 1 School of Psychology, Cardiff University, Tower Building,

More information

Discrete-Time Signals & Systems

Discrete-Time Signals & Systems Chapter 2 Discrete-Time Signals & Systems 清大電機系林嘉文 cwlin@ee.nthu.edu.tw 03-5731152 Original PowerPoint slides prepared by S. K. Mitra 2-1-1 Discrete-Time Signals: Time-Domain Representation (1/10) Signals

More information

UNIVERSITY OF OSLO. Please make sure that your copy of the problem set is complete before you attempt to answer anything.

UNIVERSITY OF OSLO. Please make sure that your copy of the problem set is complete before you attempt to answer anything. UNIVERSITY OF OSLO Faculty of mathematics and natural sciences Examination in INF3470/4470 Digital signal processing Day of examination: December 9th, 011 Examination hours: 14.30 18.30 This problem set

More information

Microphone-Array Signal Processing

Microphone-Array Signal Processing Microphone-Array Signal Processing, c Apolinárioi & Campos p. 1/27 Microphone-Array Signal Processing José A. Apolinário Jr. and Marcello L. R. de Campos {apolin},{mcampos}@ieee.org IME Lab. Processamento

More information

6.003 Signal Processing

6.003 Signal Processing 6.003 Signal Processing Week 6, Lecture A: The Discrete Fourier Transform (DFT) Adam Hartz hz@mit.edu What is 6.003? What is a signal? Abstractly, a signal is a function that conveys information Signal

More information

Generalized Sidelobe Canceller and MVDR Power Spectrum Estimation. Bhaskar D Rao University of California, San Diego

Generalized Sidelobe Canceller and MVDR Power Spectrum Estimation. Bhaskar D Rao University of California, San Diego Generalized Sidelobe Canceller and MVDR Power Spectrum Estimation Bhaskar D Rao University of California, San Diego Email: brao@ucsd.edu Reference Books 1. Optimum Array Processing, H. L. Van Trees 2.

More information

Fourier Series and Fourier Transforms

Fourier Series and Fourier Transforms Fourier Series and Fourier Transforms EECS2 (6.082), MIT Fall 2006 Lectures 2 and 3 Fourier Series From your differential equations course, 18.03, you know Fourier s expression representing a T -periodic

More information

ENT 315 Medical Signal Processing CHAPTER 2 DISCRETE FOURIER TRANSFORM. Dr. Lim Chee Chin

ENT 315 Medical Signal Processing CHAPTER 2 DISCRETE FOURIER TRANSFORM. Dr. Lim Chee Chin ENT 315 Medical Signal Processing CHAPTER 2 DISCRETE FOURIER TRANSFORM Dr. Lim Chee Chin Outline Introduction Discrete Fourier Series Properties of Discrete Fourier Series Time domain aliasing due to frequency

More information

Polynomial Root-MUSIC Algorithm for Efficient Broadband Direction Of Arrival Estimation

Polynomial 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 information

Towed M-Sequence/ Long HLA Data Analysis

Towed M-Sequence/ Long HLA Data Analysis Towed M-Sequence/ Long HLA Data Analysis Harry DeFerrari University of Miami hdeferrari@rsmas.miamai.edu Last experiment of CALOPS I Five hour tow at 6 Knots (M-sequence 255 digit, 4.08 sec., 250 Hz center

More information

A beamforming system based on the acousto-optic effect

A beamforming system based on the acousto-optic effect A beamforming system based on the acousto-optic effect Antoni Torras-Rosell atr@dfm.dtu.dk Danish Fundamental Metrology A/S, Matematiktorvet 37, 8 Kgs. yngby, Denmark. Finn Jacobsen fja@elektro.dtu.dk

More information

Measuring HRTFs of Brüel & Kjær Type 4128-C, G.R.A.S. KEMAR Type 45BM, and Head Acoustics HMS II.3 Head and Torso Simulators

Measuring HRTFs of Brüel & Kjær Type 4128-C, G.R.A.S. KEMAR Type 45BM, and Head Acoustics HMS II.3 Head and Torso Simulators Downloaded from orbit.dtu.dk on: Jan 11, 219 Measuring HRTFs of Brüel & Kjær Type 4128-C, G.R.A.S. KEMAR Type 4BM, and Head Acoustics HMS II.3 Head and Torso Simulators Snaidero, Thomas; Jacobsen, Finn;

More information

6.003 Signal Processing

6.003 Signal Processing 6.003 Signal Processing Week 6, Lecture A: The Discrete Fourier Transform (DFT) Adam Hartz hz@mit.edu What is 6.003? What is a signal? Abstractly, a signal is a function that conveys information Signal

More information

Discrete-Time David Johns and Ken Martin University of Toronto

Discrete-Time David Johns and Ken Martin University of Toronto Discrete-Time David Johns and Ken Martin University of Toronto (johns@eecg.toronto.edu) (martin@eecg.toronto.edu) University of Toronto 1 of 40 Overview of Some Signal Spectra x c () t st () x s () t xn

More information

MEDE2500 Tutorial Nov-7

MEDE2500 Tutorial Nov-7 (updated 2016-Nov-4,7:40pm) MEDE2500 (2016-2017) Tutorial 3 MEDE2500 Tutorial 3 2016-Nov-7 Content 1. The Dirac Delta Function, singularity functions, even and odd functions 2. The sampling process and

More information

To Separate Speech! A System for Recognizing Simultaneous Speech

To Separate Speech! A System for Recognizing Simultaneous Speech A System for Recognizing Simultaneous Speech John McDonough 1,2,Kenichi Kumatani 2,3,Tobias Gehrig 4, Emilian Stoimenov 4, Uwe Mayer 4, Stefan Schacht 1, Matthias Wölfel 4 and Dietrich Klakow 1 1 Spoken

More information

Bayesian Regularization and Nonnegative Deconvolution for Time Delay Estimation

Bayesian Regularization and Nonnegative Deconvolution for Time Delay Estimation University of Pennsylvania ScholarlyCommons Departmental Papers (ESE) Department of Electrical & Systems Engineering December 24 Bayesian Regularization and Nonnegative Deconvolution for Time Delay Estimation

More information

The Johns Hopkins University Department of Electrical and Computer Engineering Introduction to Linear Systems Fall 2002.

The Johns Hopkins University Department of Electrical and Computer Engineering Introduction to Linear Systems Fall 2002. The Johns Hopkins University Department of Electrical and Computer Engineering 505.460 Introduction to Linear Systems Fall 2002 Final exam Name: You are allowed to use: 1. Table 3.1 (page 206) & Table

More information

ECE 301. Division 2, Fall 2006 Instructor: Mimi Boutin Midterm Examination 3

ECE 301. Division 2, Fall 2006 Instructor: Mimi Boutin Midterm Examination 3 ECE 30 Division 2, Fall 2006 Instructor: Mimi Boutin Midterm Examination 3 Instructions:. Wait for the BEGIN signal before opening this booklet. In the meantime, read the instructions below and fill out

More information

Acoustic Vector Sensor based Speech Source Separation with Mixed Gaussian-Laplacian Distributions

Acoustic Vector Sensor based Speech Source Separation with Mixed Gaussian-Laplacian Distributions Acoustic Vector Sensor based Speech Source Separation with Mixed Gaussian-Laplacian Distributions Xiaoyi Chen, Atiyeh Alinaghi, Xionghu Zhong and Wenwu Wang Department of Acoustic Engineering, School of

More information

Massachusetts Institute of Technology

Massachusetts 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 information

BME 50500: Image and Signal Processing in Biomedicine. Lecture 5: Correlation and Power-Spectrum CCNY

BME 50500: Image and Signal Processing in Biomedicine. Lecture 5: Correlation and Power-Spectrum CCNY 1 BME 50500: Image and Signal Processing in Biomedicine Lecture 5: Correlation and Power-Spectrum Lucas C. Parra Biomedical Engineering Department CCNY http://bme.ccny.cuny.edu/faculty/parra/teaching/signal-and-image/

More information

Signals and Systems Profs. Byron Yu and Pulkit Grover Fall Midterm 2 Solutions

Signals and Systems Profs. Byron Yu and Pulkit Grover Fall Midterm 2 Solutions 8-90 Signals and Systems Profs. Byron Yu and Pulkit Grover Fall 08 Midterm Solutions Name: Andrew ID: Problem Score Max 8 5 3 6 4 7 5 8 6 7 6 8 6 9 0 0 Total 00 Midterm Solutions. (8 points) Indicate whether

More information

ECE 308 SIGNALS AND SYSTEMS SPRING 2013 Examination #2 14 March 2013

ECE 308 SIGNALS AND SYSTEMS SPRING 2013 Examination #2 14 March 2013 ECE 308 SIGNALS AND SYSTEMS SPRING 2013 Examination #2 14 March 2013 Name: Instructions: The examination lasts for 75 minutes and is closed book, closed notes. No electronic devices are permitted, including

More information

Z - Transform. It offers the techniques for digital filter design and frequency analysis of digital signals.

Z - Transform. It offers the techniques for digital filter design and frequency analysis of digital signals. Z - Transform The z-transform is a very important tool in describing and analyzing digital systems. It offers the techniques for digital filter design and frequency analysis of digital signals. Definition

More information

Grades will be determined by the correctness of your answers (explanations are not required).

Grades will be determined by the correctness of your answers (explanations are not required). 6.00 (Fall 2011) Final Examination December 19, 2011 Name: Kerberos Username: Please circle your section number: Section Time 2 11 am 1 pm 4 2 pm Grades will be determined by the correctness of your answers

More information

Problem Value

Problem Value GEORGIA INSTITUTE OF TECHNOLOGY SCHOOL of ELECTRICAL & COMPUTER ENGINEERING FINAL EXAM DATE: 30-Apr-04 COURSE: ECE-2025 NAME: GT #: LAST, FIRST Recitation Section: Circle the date & time when your Recitation

More information

Development of an Efficient Binaural Simulation for the Analysis of Structural Acoustic Data

Development of an Efficient Binaural Simulation for the Analysis of Structural Acoustic Data NASA/CR-22-211753 Development of an Efficient Binaural Simulation for the Analysis of Structural Acoustic Data Aimee L. Lalime and Marty E. Johnson Virginia Polytechnic Institute and State University Blacksburg,

More information

Chapter 5. Fourier Analysis for Discrete-Time Signals and Systems Chapter

Chapter 5. Fourier Analysis for Discrete-Time Signals and Systems Chapter Chapter 5. Fourier Analysis for Discrete-Time Signals and Systems Chapter Objec@ves 1. Learn techniques for represen3ng discrete-)me periodic signals using orthogonal sets of periodic basis func3ons. 2.

More information

EE 4372 Tomography. Carlos E. Davila, Dept. of Electrical Engineering Southern Methodist University

EE 4372 Tomography. Carlos E. Davila, Dept. of Electrical Engineering Southern Methodist University EE 4372 Tomography Carlos E. Davila, Dept. of Electrical Engineering Southern Methodist University EE 4372, SMU Department of Electrical Engineering 86 Tomography: Background 1-D Fourier Transform: F(

More information

Estimation of the Direct-Path Relative Transfer Function for Supervised Sound-Source Localization

Estimation of the Direct-Path Relative Transfer Function for Supervised Sound-Source Localization 1 Estimation of the Direct-Path Relative Transfer Function for Supervised Sound-Source Localization Xiaofei Li, Laurent Girin, Radu Horaud and Sharon Gannot arxiv:1509.03205v3 [cs.sd] 27 Jun 2016 Abstract

More information

PEAT SEISMOLOGY Lecture 12: Earthquake source mechanisms and radiation patterns II

PEAT SEISMOLOGY Lecture 12: Earthquake source mechanisms and radiation patterns II PEAT8002 - SEISMOLOGY Lecture 12: Earthquake source mechanisms and radiation patterns II Nick Rawlinson Research School of Earth Sciences Australian National University Waveform modelling P-wave first-motions

More information

BROADBAND MIMO SONAR SYSTEM: A THEORETICAL AND EXPERIMENTAL APPROACH

BROADBAND MIMO SONAR SYSTEM: A THEORETICAL AND EXPERIMENTAL APPROACH BROADBAND MIMO SONAR SYSTM: A THORTICAL AND XPRIMNTAL APPROACH Yan Pailhas a, Yvan Petillot a, Chris Capus a, Keith Brown a a Oceans Systems Lab., School of PS, Heriot Watt University, dinburgh, Scotland,

More information

ESE 531: Digital Signal Processing

ESE 531: Digital Signal Processing ESE 531: Digital Signal Processing Lec 22: April 10, 2018 Adaptive Filters Penn ESE 531 Spring 2018 Khanna Lecture Outline! Circular convolution as linear convolution with aliasing! Adaptive Filters Penn

More information

x[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

x[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 information

ESE 531: Digital Signal Processing

ESE 531: Digital Signal Processing ESE 531: Digital Signal Processing Lec 8: February 7th, 2017 Sampling and Reconstruction Lecture Outline! Review " Ideal sampling " Frequency response of sampled signal " Reconstruction " Anti-aliasing

More information

! Circular Convolution. " Linear convolution with circular convolution. ! Discrete Fourier Transform. " Linear convolution through circular

! Circular Convolution.  Linear convolution with circular convolution. ! Discrete Fourier Transform.  Linear convolution through circular Previously ESE 531: Digital Signal Processing Lec 22: April 18, 2017 Fast Fourier Transform (con t)! Circular Convolution " Linear convolution with circular convolution! Discrete Fourier Transform " Linear

More information

EE 438 Essential Definitions and Relations

EE 438 Essential Definitions and Relations May 2004 EE 438 Essential Definitions and Relations CT Metrics. Energy E x = x(t) 2 dt 2. Power P x = lim T 2T T / 2 T / 2 x(t) 2 dt 3. root mean squared value x rms = P x 4. Area A x = x(t) dt 5. Average

More information

Lecture 10. Digital Signal Processing. Chapter 7. Discrete Fourier transform DFT. Mikael Swartling Nedelko Grbic Bengt Mandersson. rev.

Lecture 10. Digital Signal Processing. Chapter 7. Discrete Fourier transform DFT. Mikael Swartling Nedelko Grbic Bengt Mandersson. rev. Lecture 10 Digital Signal Processing Chapter 7 Discrete Fourier transform DFT Mikael Swartling Nedelko Grbic Bengt Mandersson rev. 016 Department of Electrical and Information Technology Lund University

More information

An Observer for Phased Microphone Array Signal Processing with Nonlinear Output

An Observer for Phased Microphone Array Signal Processing with Nonlinear Output 2010 Asia-Pacific International Symposium on Aerospace Technology An Observer for Phased Microphone Array Signal Processing with Nonlinear Output Bai Long 1,*, Huang Xun 2 1 Department of Mechanics and

More information

Spatial 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 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 information

EE538 Final Exam Fall :20 pm -5:20 pm PHYS 223 Dec. 17, Cover Sheet

EE538 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 information

Signals & Systems. Chapter 7: Sampling. Adapted from: Lecture notes from MIT, Binghamton University, and Purdue. Dr. Hamid R.

Signals & Systems. Chapter 7: Sampling. Adapted from: Lecture notes from MIT, Binghamton University, and Purdue. Dr. Hamid R. Signals & Systems Chapter 7: Sampling Adapted from: Lecture notes from MIT, Binghamton University, and Purdue Dr. Hamid R. Rabiee Fall 2013 Outline 1. The Concept and Representation of Periodic Sampling

More information

Introduction to Acoustics Exercises

Introduction to Acoustics Exercises . 361-1-3291 Introduction to Acoustics Exercises 1 Fundamentals of acoustics 1. Show the effect of temperature on acoustic pressure. Hint: use the equation of state and the equation of state at equilibrium.

More information

Comparison of RTF Estimation Methods between a Head-Mounted Binaural Hearing Device and an External Microphone

Comparison of RTF Estimation Methods between a Head-Mounted Binaural Hearing Device and an External Microphone Comparison of RTF Estimation Methods between a Head-Mounted Binaural Hearing Device and an External Microphone Nico Gößling, Daniel Marquardt and Simon Doclo Department of Medical Physics and Acoustics

More information

ESE 531: Digital Signal Processing

ESE 531: Digital Signal Processing ESE 531: Digital Signal Processing Lec 8: February 12th, 2019 Sampling and Reconstruction Lecture Outline! Review " Ideal sampling " Frequency response of sampled signal " Reconstruction " Anti-aliasing

More information

Problem Value

Problem Value GEORGIA INSTITUTE OF TECHNOLOGY SCHOOL of ELECTRICAL & COMPUTER ENGINEERING FINAL EXAM DATE: 30-Apr-04 COURSE: ECE-2025 NAME: GT #: LAST, FIRST Recitation Section: Circle the date & time when your Recitation

More information

Array Signal Processing Algorithms for Localization and Equalization in Complex Acoustic Channels

Array Signal Processing Algorithms for Localization and Equalization in Complex Acoustic Channels Array Signal Processing Algorithms for Localization and Equalization in Complex Acoustic Channels Dumidu S. Talagala B.Sc. Eng. (Hons.), University of Moratuwa, Sri Lanka November 2013 A thesis submitted

More information

Data-based Binaural Synthesis Including Rotational and Translatory Head-Movements

Data-based Binaural Synthesis Including Rotational and Translatory Head-Movements Data-based Binaural Synthesis Including Rotational and Translatory Head-Movements Frank Schultz 1 and Sascha Spors 1 1 Institute of Communications Engineering, Universität Rostock, R.-Wagner-Str. 31 (H8),

More information

Final Exam ECE301 Signals and Systems Friday, May 3, Cover Sheet

Final Exam ECE301 Signals and Systems Friday, May 3, Cover Sheet Name: Final Exam ECE3 Signals and Systems Friday, May 3, 3 Cover Sheet Write your name on this page and every page to be safe. Test Duration: minutes. Coverage: Comprehensive Open Book but Closed Notes.

More information

2.161 Signal Processing: Continuous and Discrete Fall 2008

2.161 Signal Processing: Continuous and Discrete Fall 2008 IT OpenCourseWare http://ocw.mit.edu 2.6 Signal Processing: Continuous and Discrete Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. ASSACHUSETTS

More information

Signal representations: Cepstrum

Signal representations: Cepstrum Signal representations: Cepstrum Source-filter separation for sound production For speech, source corresponds to excitation by a pulse train for voiced phonemes and to turbulence (noise) for unvoiced phonemes,

More information

Shift Property of z-transform. Lecture 16. More z-transform (Lathi 5.2, ) More Properties of z-transform. Convolution property of z-transform

Shift Property of z-transform. Lecture 16. More z-transform (Lathi 5.2, ) More Properties of z-transform. Convolution property of z-transform Shift Property of -Transform If Lecture 6 More -Transform (Lathi 5.2,5.4-5.5) then which is delay causal signal by sample period. If we delay x[n] first: Peter Cheung Department of Electrical & Electronic

More information

ADAPTIVE ANTENNAS. SPATIAL BF

ADAPTIVE 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 information

SPACIOUSNESS OF SOUND FIELDS CAPTURED BY SPHERICAL MICROPHONE ARRAYS

SPACIOUSNESS OF SOUND FIELDS CAPTURED BY SPHERICAL MICROPHONE ARRAYS BEN GURION UNIVERSITY OF THE NEGEV FACULTY OF ENGINEERING SCIENCES DEPARTMENT OF ELECTRICAL AND COMPUTER ENGINEERING SPACIOUSNESS OF SOUND FIELDS CAPTURED BY SPHERICAL MICROPHONE ARRAYS THESIS SUBMITTED

More information

IMPROVED MULTI-MICROPHONE NOISE REDUCTION PRESERVING BINAURAL CUES

IMPROVED MULTI-MICROPHONE NOISE REDUCTION PRESERVING BINAURAL CUES IMPROVED MULTI-MICROPHONE NOISE REDUCTION PRESERVING BINAURAL CUES Andreas I. Koutrouvelis Richard C. Hendriks Jesper Jensen Richard Heusdens Circuits and Systems (CAS) Group, Delft University of Technology,

More information

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 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 information

Microphone-Array Signal Processing

Microphone-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 information

Sensor Tasking and Control

Sensor Tasking and Control Sensor Tasking and Control Sensing Networking Leonidas Guibas Stanford University Computation CS428 Sensor systems are about sensing, after all... System State Continuous and Discrete Variables The quantities

More information

ECE 301 Division 1, Fall 2006 Instructor: Mimi Boutin Final Examination

ECE 301 Division 1, Fall 2006 Instructor: Mimi Boutin Final Examination ECE 30 Division, all 2006 Instructor: Mimi Boutin inal Examination Instructions:. Wait for the BEGIN signal before opening this booklet. In the meantime, read the instructions below and fill out the requested

More information

Overview of Beamforming

Overview 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 information

Probabilistic Structure from Sound and Probabilistic Sound Source Localization

Probabilistic Structure from Sound and Probabilistic Sound Source Localization Probabilistic Structure from Sound and Probabilistic Sound Source Localization Chi-Hao Lin and Chieh-Chih Wang Department of Computer Science and Information Engineering Graduate Institute of Networking

More information

EE Homework 13 - Solutions

EE Homework 13 - Solutions EE3054 - Homework 3 - Solutions. (a) The Laplace transform of e t u(t) is s+. The pole of the Laplace transform is at which lies in the left half plane. Hence, the Fourier transform is simply the Laplace

More information

Your solutions for time-domain waveforms should all be expressed as real-valued functions.

Your solutions for time-domain waveforms should all be expressed as real-valued functions. ECE-486 Test 2, Feb 23, 2017 2 Hours; Closed book; Allowed calculator models: (a) Casio fx-115 models (b) HP33s and HP 35s (c) TI-30X and TI-36X models. Calculators not included in this list are not permitted.

More information

Azimuth-elevation direction finding, using one fourcomponent acoustic vector-sensor spread spatially along a straight line

Azimuth-elevation direction finding, using one fourcomponent acoustic vector-sensor spread spatially along a straight line Volume 23 http://acousticalsociety.org/ 169th Meeting of the Acoustical Society of America Pittsburgh, Pennsylvania 18-22 May 2015 Signal Processing in Acoustics: Paper 4aSP4 Azimuth-elevation direction

More information

Down-Sampling (4B) Young Won Lim 11/15/12

Down-Sampling (4B) Young Won Lim 11/15/12 Down-Sampling (B) /5/ Copyright (c) 9,, Young W. Lim. Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version. or any later

More information

New Mexico State University Klipsch School of Electrical Engineering EE312 - Signals and Systems I Fall 2015 Final Exam

New Mexico State University Klipsch School of Electrical Engineering EE312 - Signals and Systems I Fall 2015 Final Exam New Mexico State University Klipsch School of Electrical Engineering EE312 - Signals and Systems I Fall 2015 Name: Solve problems 1 3 and two from problems 4 7. Circle below which two of problems 4 7 you

More information

Estimation of the Direct-Path Relative Transfer Function for Supervised Sound-Source Localization

Estimation of the Direct-Path Relative Transfer Function for Supervised Sound-Source Localization IEEE/ACM TRANSACTIONS ON AUDIO, SPEECH, AND LANGUAGE PROCESSING, VOL. 24, NO. 11, NOVEMBER 2016 2171 Estimation of the Direct-Path Relative Transfer Function for Supervised Sound-Source Localization Xiaofei

More information

ECE 301. Division 3, Fall 2007 Instructor: Mimi Boutin Midterm Examination 3

ECE 301. Division 3, Fall 2007 Instructor: Mimi Boutin Midterm Examination 3 ECE 30 Division 3, all 2007 Instructor: Mimi Boutin Midterm Examination 3 Instructions:. Wait for the BEGIN signal before opening this booklet. In the meantime, read the instructions below and fill out

More information