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

Size: px
Start display at page:

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

Transcription

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

2 Ear pinna auditory canal cochlea outer middle inner 2

3 Eye Ear Lens? Retina? Photoreceptors (light -> chemical) Ganglion cells (spikes)?? Optic nerve? 3

4 Eye Ear Lens Outer ear Retina Cochlea Photoreceptors (light -> chemical) Ganglion cells (spikes) hair cells (mechanical -> chemical) Ganglion cells (spikes) Optic nerve VestibuloCochlear nerve 4

5 Basilar Membrane BM fibres have bandpass frequency mechanical responses. 20,000 Hz 20 Hz 5

6 Basilar Membrane: Place code ( tonotopic ) Nerve cells (hair + ganglion) are distributed along the BM. They have similar bandpass frequency response functions. 20,000 Hz 20 Hz 6

7 Bandpass responses (more details next lecture) ,000 7

8 Neural coding of sound in cochlea Basilar membrane responds by vibrating with sound. Hair cells at each BM location release neurotransmitter that signal BM amplitude at that location Ganglion cells respond to neurotransmitter signals by spiking 8

9 Louder sound within frequency band greater amplitude of BM vibration at that location greater release of neurotransmitter by hair cell greater probability of spike at each peak of filtered wave 9

10 Hair cell neurotransmitter release can signal exact timing of BM amplitude peaks for frequencies up to ~2 khz. For higher frequencies, hair cells encode only the envelope of BM vibrations. t 10

11 Timing of ganglion cell spikes: for frequencies up to 2 KHz ( phase locking ) BM vibration Spikes Hair cells release more neurotransmitter at BM amplitude peaks. Ganglion cells respond to neurotransmitter peaks by spiking. This allows exact timing of BM vibrations to be encoded by spikes. 11

12 Ganglion cells cannot spike faster than 500 times per second. So we need many ganglion cells for each hair cell. 3,000 hair cells in each cochlea (left and right) cochlear nerve (to brain) 30,00 ganglion cells in each cochlea 12

13 Volley code Different ganglion cells at same spatial position on BM 13

14 From cochlea to brain stem BRAIN STEM cochlea cochlear nucleus lateral and medial superior olive (LSO, MSO) auditory cortex 14

15 Tonotopic maps lateral superior olive (LSO) high ω low ω medial superior olive (MSO) auditory nerve cochlear nucleus (CN) cochlea 15

16 Binaural Hearing LSO MSO MSO LSO high ω low ω low ω high ω CN CN 16

17 Binaural Hearing MSO combines low frequency signals. LSO MSO MSO LSO high ω low ω low ω high ω CN CN 17

18 Binaural Hearing LSO combines high frequency signals. LSO MSO MSO LSO high ω high ω CN CN 18

19 Levels of Analysis high - what is the task? what problem is being solved? - brain areas and pathways - neural coding - neural mechanisms low 19

20 For high frequency bands, the head casts a shadow the timing of the peaks cannot be accurately coded by the spikes (only the rate of spikes is informative) For low frequency bands, the head casts a weak shadow only the timing of the peaks can be encoded by spikes 20

21 Duplex theory of binaural hearing (Rayleigh, 1907) level differences computed for higher frequencies (ILD -- interaural level differences) timing differences computed for lower frequencies (ITD - interaural timing differences) 21

22 Level differences (high frequencies) Excitatory input comes from the ear on the same side. Inhibitory input comes from ear on the opposite side. high ω LSO MSO MSO + + LSO high ω CN CN 22

23 Timing differences (low frequencies) Sum excitatory input from both sides. Reminiscent of binocular complex cells in V1? LSO MSO MSO LSO high ω CN + low ω + + low ω + CN high ω 23

24 Jeffress Model (1948) for timing differences from left ear E D C B A from right ear C B D A E 24

25 Spike Timing precision required for Jeffress Model? from left ear E D C B A from right ear distance = 1 10 millimetres speed of spike = 10 metres second 1 time = distance speed = millisecond See Exercises 19 Q2c 25

26 Jeffress model remains controversial. It is not known exactly how coincidence detection occurs in MSO. ω A B C D E Coincidence detection for each low frequency band A B C D E A B C D E A B C D E 26

27 Naïve Computational Model of Source Localization (Recall lecture 20) I l (t) = α I r (t τ) + ε(t) shadow delay model error Find the α and τ that minimize T t=1 { I l (t) α I r (t τ) } 2 where τ < 0.5 ms. 27

28 Timing difference: find the τ that maximizes t I l (t) I r (t τ). Level difference: T t=1 I l (t) 2 10 log 10 T I r (t) 2 t=1 28

29 For each low frequency band j, find the τ that maximizes t I left j (t) I right j (t τ). (or use summation model similar to binocular cells or Jeffress model) An estimated value of delay τ in frequency band j is consistent with various possible source directions ( φ, θ ). Similar to cone of confusion, but more general because of frequency dependence 29

30 For each high frequency band j, compute interaural level difference (ILD) : T t=1 I j left (t) 2 ILD j = 10 log 10 T t=1 I j right (t) 2 What does each ILD j tell us? 30

31 Recall head related impulse response function (HRIR) from last lecture.. If the source direction is (q, f), and g j t is the filter for band j. then j I left t; φ, θ = g j t h left (t; φ, θ) I src t; φ, θ j I right t; φ, θ = g j t h right (t; φ, θ) I src t; φ, θ 31

32 Take the Fourier transform and apply convolution theorem : j I left ω; φ, θ = g j ω h left (ω; φ, θ) I src ω; φ, θ j I right ω; φ, θ = g j ω h right (ω; φ, θ) I src ω; φ, θ 32

33 Take the Fourier transform and apply convolution theorem : j I left ω; φ, θ = g j ω h left (ω; φ, θ) I src ω; φ, θ j I right ω; φ, θ = g j ω h right (ω; φ, θ) I src ω; φ, θ If there is just one source direction (φ, θ), then for each frequency ω within band j j I left j I right ω h left (ω; φ, θ) ω h right (ω; φ, θ) 33

34 One can show using Parseval s theorem of Fourier transforms that if h left (ω; φ, θ) and h right (ω; φ, θ) are approximately constant within band j, then: T t=1 I j left (t) 2 j h left ( φ, θ) 2 T t=1 I j right (t) 2 j h right ( φ, θ) 2 34

35 One can show using Parseval s theorem of Fourier transforms that if h left (ω; φ, θ) and h right (ω; φ, θ) are approximately constant within band j, then: T t=1 I j left (t) 2 j h left ( φ, θ) 2 T t=1 I j right (t) 2 j h right ( φ, θ) 2 The ear can measure this and can infer source directions ( φ, θ) that are consistent with it. 35

36 Interaural Level Difference (db) as a function of (φ, θ) for two fixed ω. 700 Hz 11,000 Hz Each iso-contour in each frequency band is consistent with a measured level difference (db). 36

37 Monaural spectral cues (Spatial localization with one ear?) I j t; φ, θ = g j t h(t; φ, θ) I src t; φ, θ I j ω; φ, θ = g j ω h j (ω; φ, θ) I src ω; φ, θ Pattern of peaks and notches across bands will be due to HRTF, not to the source. If the source is noise, then all frequencies make the same contribution on average. 37

38 Pinnal notch frequency varies with elevation of source e.g. in the medial plane. I j ω; φ, θ = g j ω h j (ω; φ, θ) I src ω; φ, θ 0 0 Azimuth θ = 0 HRTF from last lecture e.g. medial plane 38

39 Levels of Analysis high - what is the task? what problem is being solved? Source localization using level and timing differences within frequency channels. - brain areas and pathways (cochlea to CN to MSO and LSO in the brainstem) - neural coding (gave sketch only) - neural mechanisms low 39

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

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

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

Lecture 17 Source Separation

Lecture 17 Source Separation EECS 225D Audio Signal Processing in Humans and Machines Lecture 17 Source Separation 2012-3-19 Professor Nelson Morgan today s lecture by John Lazzaro www.icsi.berkeley.edu/eecs225d/spr12/ EECS 225D L17:

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

Time-Series Analysis for Ear-Related and Psychoacoustic Metrics

Time-Series Analysis for Ear-Related and Psychoacoustic Metrics Time-Series Analysis for Ear-Related and Psychoacoustic Metrics V. Mellert, H. Remmers, R. Weber, B. Schulte-Fortkamp how to analyse p(t) to obtain an earrelated parameter? general remarks on acoustical

More information

Computational Perception. Sound Localization 1

Computational Perception. Sound Localization 1 Computational Perception 15-485/785 January 17, 2008 Sound Localization 1 Orienting sound localization visual pop-out eye/body movements attentional shift 2 The Problem of Sound Localization What are the

More information

The Spectral Differences in Sounds and Harmonics. Through chemical and biological processes, we, as humans, can detect the movement and

The Spectral Differences in Sounds and Harmonics. Through chemical and biological processes, we, as humans, can detect the movement and John Chambers Semester Project 5/2/16 The Spectral Differences in Sounds and Harmonics Through chemical and biological processes, we, as humans, can detect the movement and vibrations of molecules through

More information

Logic. Intro to Neuroscence: Neuromorphic Engineering 12/9/2013. (c) S. Liu and T. Delbruck, Inst. of Neuroinformatics, UZH-ETH Zurich 1

Logic. Intro to Neuroscence: Neuromorphic Engineering 12/9/2013. (c) S. Liu and T. Delbruck, Inst. of Neuroinformatics, UZH-ETH Zurich 1 Introductory Course in Neuroscience Neuromorphic Shih-Chii Liu Inst. of Neuroinformatics http://www.ini.uzh.ch/~shih/wiki/doku.php?id=introneuro What is neuromorphic engineering? It consists of embodying

More information

B.S. (distinction), Cornell University (1993) at the. June 1995

B.S. (distinction), Cornell University (1993) at the. June 1995 A Computational Model of Spatial Hearing by Keith Dana Martin B.S. (distinction), Cornell University (1993) Submitted to the Department of Electrical Engineering and Computer Science in partial fulllment

More information

SENSORY PROCESSES PROVIDE INFORMATION ON ANIMALS EXTERNAL ENVIRONMENT AND INTERNAL STATUS 34.4

SENSORY PROCESSES PROVIDE INFORMATION ON ANIMALS EXTERNAL ENVIRONMENT AND INTERNAL STATUS 34.4 SENSORY PROCESSES PROVIDE INFORMATION ON ANIMALS EXTERNAL ENVIRONMENT AND INTERNAL STATUS 34.4 INTRODUCTION Animals need information about their external environments to move, locate food, find mates,

More information

Response-Field Dynamics in the Auditory Pathway

Response-Field Dynamics in the Auditory Pathway Response-Field Dynamics in the Auditory Pathway Didier Depireux Powen Ru Shihab Shamma Jonathan Simon Work supported by grants from the Office of Naval Research, a training grant from the National Institute

More information

ABSTRACT A MODEL OF THE AVIAN SUPERIOR OLIVARY. Department of Electrical Engineering

ABSTRACT A MODEL OF THE AVIAN SUPERIOR OLIVARY. Department of Electrical Engineering ABSTRACT Title of Thesis: A MODEL OF THE AVIAN SUPERIOR OLIVARY NUCLEUS Degree candidate: Raul Rodriguez Esteban Degree and year: Master of Science, 2002 Thesis directed by: Professor Jonathan Z. Simon

More information

15 Grossberg Network 1

15 Grossberg Network 1 Grossberg Network Biological Motivation: Vision Bipolar Cell Amacrine Cell Ganglion Cell Optic Nerve Cone Light Lens Rod Horizontal Cell Retina Optic Nerve Fiber Eyeball and Retina Layers of Retina The

More information

Neural coding Ecological approach to sensory coding: efficient adaptation to the natural environment

Neural coding Ecological approach to sensory coding: efficient adaptation to the natural environment Neural coding Ecological approach to sensory coding: efficient adaptation to the natural environment Jean-Pierre Nadal CNRS & EHESS Laboratoire de Physique Statistique (LPS, UMR 8550 CNRS - ENS UPMC Univ.

More information

Figure by MIT OCW. After Melcher and Kiang, HST MODELING ISSUES IN HEARING AND SPEECH SPRING Vertex - Ipsilateral Earbar

Figure by MIT OCW. After Melcher and Kiang, HST MODELING ISSUES IN HEARING AND SPEECH SPRING Vertex - Ipsilateral Earbar Harvard-MIT Division of Health Sciences and Technology HST.750: Modeling Issues in Speech and Hearing, Spring 2006 Course Directors: Dr. Christopher Shera and Dr. Jennifer Melcher HST 750 - MODELIG ISSUES

More information

Identify the structure labelled 1.

Identify the structure labelled 1. Identify the structure labelled 1. Identify the structure labelled 1. Cornea Identify the structure labelled 2. Identify the structure labelled 2. Sclera 1 2 Identify the structure labelled 1. Pupil Identify

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

Feature extraction 1

Feature extraction 1 Centre for Vision Speech & Signal Processing University of Surrey, Guildford GU2 7XH. Feature extraction 1 Dr Philip Jackson Cepstral analysis - Real & complex cepstra - Homomorphic decomposition Filter

More information

Tim Habigt. May 30, 2014

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

More information

Which of the six boxes below cannot be made from this unfolded box? (There may be more than ffi. ffi

Which of the six boxes below cannot be made from this unfolded box? (There may be more than ffi. ffi Which of the six boxes below cannot be made from this unfolded box? (There may be more than one.) @ ffi w @ @ ffi c Biophysics @ YORK Peripheral sensory transduction Christopher Bergevin (York University,

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

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

Adaptation in the Neural Code of the Retina

Adaptation in the Neural Code of the Retina Adaptation in the Neural Code of the Retina Lens Retina Fovea Optic Nerve Optic Nerve Bottleneck Neurons Information Receptors: 108 95% Optic Nerve 106 5% After Polyak 1941 Visual Cortex ~1010 Mean Intensity

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

Chapter 27 Interaural Time Difference Thresholds as a Function of Frequency

Chapter 27 Interaural Time Difference Thresholds as a Function of Frequency Chapter 27 Interaural ime Difference hresholds as a Function of Frequency William M. Hartmann, Larisa Dunai, and ianshu Qu Abstract Different models of the binaural system make different predictions for

More information

Tuning tuning curves. So far: Receptive fields Representation of stimuli Population vectors. Today: Contrast enhancment, cortical processing

Tuning tuning curves. So far: Receptive fields Representation of stimuli Population vectors. Today: Contrast enhancment, cortical processing Tuning tuning curves So far: Receptive fields Representation of stimuli Population vectors Today: Contrast enhancment, cortical processing Firing frequency N 3 s max (N 1 ) = 40 o N4 N 1 N N 5 2 s max

More information

INFORMATION-THEORETIC PERFORMANCE ANALYSIS OF AZIMUTHAL LOCALIZATION FOR SPATIAL-AUDITORY DISPLAY OF BEAMFORMED SONAR DATA. Jason E.

INFORMATION-THEORETIC PERFORMANCE ANALYSIS OF AZIMUTHAL LOCALIZATION FOR SPATIAL-AUDITORY DISPLAY OF BEAMFORMED SONAR DATA. Jason E. INFORMATION-THEORETIC PERFORMANCE ANALYSIS OF AZIMUTHAL LOCALIZATION FOR SPATIAL-AUDITORY DISPLAY OF BEAMFORMED SONAR DATA Jason E. Summers Applied Research in Acoustics (ARiA) LLC Washington, D.C. jason.e.summers@ariacoustics.com

More information

Sensors. Sensory Physiology. Transduction. Types of Environmental Stimuli. Chemoreception. Taste Buds (Contact Chemoreceptors)

Sensors. Sensory Physiology. Transduction. Types of Environmental Stimuli. Chemoreception. Taste Buds (Contact Chemoreceptors) Sensors Sensory Physiology Chapter 13 Detect changes in environmental conditions Primary Sensors neurons modified to undergo action potentials in response to specific stimuli (e.g. chemical, mechanical)

More information

Time Varying Loudness as a Means of Quantifying Loudness in Nonlinearly Propagated Acoustical Signals. S. Hales Swift. Kent L. Gee, faculty advisor

Time Varying Loudness as a Means of Quantifying Loudness in Nonlinearly Propagated Acoustical Signals. S. Hales Swift. Kent L. Gee, faculty advisor Time Varying Loudness as a Means of Quantifying Loudness in Nonlinearly Propagated Acoustical Signals by S. Hales Swift Kent L. Gee, faculty advisor A capstone project report submitted to the faculty of

More information

5. Algorithms for the Evaluation of Interaural Phase Differences ("Phase-Difference-Cocktail-Party-Processor")

5. Algorithms for the Evaluation of Interaural Phase Differences (Phase-Difference-Cocktail-Party-Processor) - 46-5. Algorithms for the Evaluation of Interaural Phase Differences ("Phase-Difference-Cocktail-Party-Processor") 5.. Requirements for a Cocktail-Party-Processor A system for the analysis of ear signals

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

Comparison between the equalization and cancellation model and state of the art beamforming techniques

Comparison between the equalization and cancellation model and state of the art beamforming techniques Comparison between the equalization and cancellation model and state of the art beamforming techniques FREDRIK GRAN 1,*,JESPER UDESEN 1, and Andrew B. Dittberner 2 Fredrik Gran 1,*, Jesper Udesen 1,*,

More information

ENSC327 Communications Systems 2: Fourier Representations. Jie Liang School of Engineering Science Simon Fraser University

ENSC327 Communications Systems 2: Fourier Representations. Jie Liang School of Engineering Science Simon Fraser University ENSC327 Communications Systems 2: Fourier Representations Jie Liang School of Engineering Science Simon Fraser University 1 Outline Chap 2.1 2.5: Signal Classifications Fourier Transform Dirac Delta Function

More information

Identification and separation of noises with spectro-temporal patterns

Identification and separation of noises with spectro-temporal patterns PROCEEDINGS of the 22 nd International Congress on Acoustics Soundscape, Psychoacoustics and Urban Environment: Paper ICA2016-532 Identification and separation of noises with spectro-temporal patterns

More information

c 2011 Jeremy A. Tidemann

c 2011 Jeremy A. Tidemann c 2011 Jeremy A. Tidemann CHARACTERIZATION OF THE HEAD-RELATED TRANSFER FUNCTION USING CHIRP AND MAXIMUM LENGTH SEQUENCE EXCITATION SIGNALS BY JEREMY A. TIDEMANN THESIS Submitted in partial fulfillment

More information

Basilar membrane models

Basilar membrane models Technical University of Denmark Basilar membrane models 2 nd semester project DTU Electrical Engineering Acoustic Technology Spring semester 2008 Group 5 Troels Schmidt Lindgreen s073081 David Pelegrin

More information

Chapter 1 Fundamentals of Sound Waves -1

Chapter 1 Fundamentals of Sound Waves -1 Chapter 1 Fundamentals of Sound Waves -1 Sound Sections What is Sound? Producing a Sound Wave Characteristics of Sound Waves The Speed of Sound Intensity of Sound Waves What is Sound? Sound is such a common

More information

REAL TIME CALCULATION OF THE HEAD RELATED TRANSFER FUNCTION BASED ON THE BOUNDARY ELEMENT METHOD

REAL TIME CALCULATION OF THE HEAD RELATED TRANSFER FUNCTION BASED ON THE BOUNDARY ELEMENT METHOD REAL TIME CALCULATION OF THE HEAD RELATED TRANSFER FUNCTION BASED ON THE BOUNDARY ELEMENT METHOD Shiro Ise Makoto Otani Department of Architecture and Architectural Systems Faculty of Engineering, Kyoto

More information

Consider the following spike trains from two different neurons N1 and N2:

Consider the following spike trains from two different neurons N1 and N2: About synchrony and oscillations So far, our discussions have assumed that we are either observing a single neuron at a, or that neurons fire independent of each other. This assumption may be correct in

More information

Introduction Principles of Signaling and Organization p. 3 Signaling in Simple Neuronal Circuits p. 4 Organization of the Retina p.

Introduction Principles of Signaling and Organization p. 3 Signaling in Simple Neuronal Circuits p. 4 Organization of the Retina p. Introduction Principles of Signaling and Organization p. 3 Signaling in Simple Neuronal Circuits p. 4 Organization of the Retina p. 5 Signaling in Nerve Cells p. 9 Cellular and Molecular Biology of Neurons

More information

Quantitative Understanding in Biology Module IV: ODEs Lecture III: Fourier Analysis

Quantitative Understanding in Biology Module IV: ODEs Lecture III: Fourier Analysis Quantitative Understanding in Biology Module IV: ODEs Lecture III: Fourier Analysis Introduction Throughout this course we have seen examples of complex mathematical phenomena being represented as linear

More information

Quantitative Understanding in Biology Fourier Analysis and Signal Processing

Quantitative Understanding in Biology Fourier Analysis and Signal Processing Quantitative Understanding in Biology Fourier Analysis and Signal Processing Representing Mathematical Functions as Linear Combinations of Basis Functions Throughout this course we have seen examples of

More information

Binaural Interaction in the Nucleus Laminaris of the Barn Owl: A Quantitative Model

Binaural Interaction in the Nucleus Laminaris of the Barn Owl: A Quantitative Model Research Report BrainWorks 1991-01 ISSN 1615-2042 Binaural Interaction in the Nucleus Laminaris of the Barn Owl: A Quantitative Model Sonja Grün, Hermann Wagner, Catherine Carr, Ad Aertsen Neurobiology

More information

Cell division takes place next to the RPE. Neuroblastic cells have the capacity to differentiate into any of the cell types found in the mature retina

Cell division takes place next to the RPE. Neuroblastic cells have the capacity to differentiate into any of the cell types found in the mature retina RPE is a monolayer of hexagonal shaped neural epithelial cells that have the same embryological origin as the neural retina. They mature before the neural retina and play a key role in metabolic support

More information

The Bayesian Brain. Robert Jacobs Department of Brain & Cognitive Sciences University of Rochester. May 11, 2017

The Bayesian Brain. Robert Jacobs Department of Brain & Cognitive Sciences University of Rochester. May 11, 2017 The Bayesian Brain Robert Jacobs Department of Brain & Cognitive Sciences University of Rochester May 11, 2017 Bayesian Brain How do neurons represent the states of the world? How do neurons represent

More information

Transformation of stimulus correlations by the retina

Transformation of stimulus correlations by the retina Transformation of stimulus correlations by the retina Kristina Simmons (University of Pennsylvania) and Jason Prentice, (now Princeton University) with Gasper Tkacik (IST Austria) Jan Homann (now Princeton

More information

Fourier Series and the Discrete Fourier Expansion

Fourier Series and the Discrete Fourier Expansion 2 2.5.5 Fourier Series and the Discrete Fourier Expansion Matthew Lincoln Adrienne Carter sillyajc@yahoo.com December 5, 2 Abstract This article is intended to introduce the Fourier series and the Discrete

More information

A COMPUTATIONAL SOFTWARE FOR NOISE MEASUREMENT AND TOWARD ITS IDENTIFICATION

A COMPUTATIONAL SOFTWARE FOR NOISE MEASUREMENT AND TOWARD ITS IDENTIFICATION A COMPUTATIONAL SOFTWARE FOR NOISE MEASUREMENT AND TOWARD ITS IDENTIFICATION M. SAKURAI Yoshimasa Electronic Inc., Daiichi-Nishiwaki Bldg., -58-0 Yoyogi, Shibuya, Tokyo, 5-0053 Japan E-mail: saku@ymec.co.jp

More information

Limulus. The Neural Code. Response of Visual Neurons 9/21/2011

Limulus. The Neural Code. Response of Visual Neurons 9/21/2011 Crab cam (Barlow et al., 2001) self inhibition recurrent inhibition lateral inhibition - L16. Neural processing in Linear Systems: Temporal and Spatial Filtering C. D. Hopkins Sept. 21, 2011 The Neural

More information

Some sensory receptors are specialized neurons while others are specialized cells that regulate neurons Figure 50.4

Some sensory receptors are specialized neurons while others are specialized cells that regulate neurons Figure 50.4 1 2 3 4 5 6 7 8 9 10 Sensory and Motor Mechanisms Chapter 50 Sensory receptors transduce stimulus energy and transmit signals to the central nervous system Sensory Pathways Sensory pathways have four basic

More information

Linear and Non-Linear Responses to Dynamic Broad-Band Spectra in Primary Auditory Cortex

Linear and Non-Linear Responses to Dynamic Broad-Band Spectra in Primary Auditory Cortex Linear and Non-Linear Responses to Dynamic Broad-Band Spectra in Primary Auditory Cortex D. J. Klein S. A. Shamma J. Z. Simon D. A. Depireux,2,2 2 Department of Electrical Engineering Supported in part

More information

2A1H Time-Frequency Analysis II

2A1H Time-Frequency Analysis II 2AH Time-Frequency Analysis II Bugs/queries to david.murray@eng.ox.ac.uk HT 209 For any corrections see the course page DW Murray at www.robots.ox.ac.uk/ dwm/courses/2tf. (a) A signal g(t) with period

More information

A Comparison of Computational Precedence Models for Source Separation in Reverberant Environments

A Comparison of Computational Precedence Models for Source Separation in Reverberant Environments A Comparison of Computational Precedence Models for Source Separation in Reverberant Environments CHRISTOPHER HUMMERSONE, 1 AES Member, RUSSELL MASON, 1 AES Member, AND c.hummersone@surrey.ac.uk r.mason@surrey.ac.uk

More information

The effect of impedance on interaural azimuth cues derived from a spherical head model a)

The effect of impedance on interaural azimuth cues derived from a spherical head model a) The effect of impedance on interaural azimuth cues derived from a spherical head model a) Bradley E. Treeby, b Roshun M. Paurobally, and Jie Pan Centre for Acoustics, Dynamics and Vibration, School of

More information

An EM Algorithm for Localizing Multiple Sound Sources in Reverberant Environments

An EM Algorithm for Localizing Multiple Sound Sources in Reverberant Environments An EM Algorithm for Localizing Multiple Sound Sources in Reverberant Environments Michael I. Mandel, Daniel P. W. Ellis LabROSA, Dept. of Electrical Engineering Columbia University New York, NY {mim,dpwe}@ee.columbia.edu

More information

Transaural Audio - The reproduction of binaural signals over loudspeakers. Fabio Kaiser

Transaural Audio - The reproduction of binaural signals over loudspeakers. Fabio Kaiser Transaural Audio - The reproduction of binaural signals over loudspeakers Fabio Kaiser Outline 1 Introduction 2 Inversion of non-minimum phase filters Inversion techniques 3 Implementation of CTC 4 Objective

More information

Spatial Vision: Primary Visual Cortex (Chapter 3, part 1)

Spatial Vision: Primary Visual Cortex (Chapter 3, part 1) Spatial Vision: Primary Visual Cortex (Chapter 3, part 1) Lecture 6 Jonathan Pillow Sensation & Perception (PSY 345 / NEU 325) Princeton University, Spring 2015 1 Chapter 2 remnants 2 Receptive field:

More information

Combination Tones As a Phenomenon of the Central Level in Auditory System

Combination Tones As a Phenomenon of the Central Level in Auditory System P a g e 1 Combination Tones As a Phenomenon of the Central Level in Auditory System Tadeusz ZIĘBAKOWSKI West Pomeranian University of Technology Institute of Manufacturing Engineering Al. Piastów 19, 70-310

More information

Obtaining objective, content-specific room acoustical parameters using auditory modeling

Obtaining objective, content-specific room acoustical parameters using auditory modeling Obtaining objective, content-specific room acoustical parameters using auditory modeling Jasper van Dorp Schuitman Philips Research, Digital Signal Processing Group, Eindhoven, The Netherlands Diemer de

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

Exploring Nonlinear Oscillator Models for the Auditory Periphery

Exploring Nonlinear Oscillator Models for the Auditory Periphery Exploring Nonlinear Oscillator Models for the Auditory Periphery Andrew Binder Dr. Christopher Bergevin, Supervisor October 31, 2008 1 Introduction Sound-induced vibrations are converted into electrical

More information

Physical Acoustics. Hearing is the result of a complex interaction of physics, physiology, perception and cognition.

Physical Acoustics. Hearing is the result of a complex interaction of physics, physiology, perception and cognition. Physical Acoustics Hearing, auditory perception, or audition is the ability to perceive sound by detecting vibrations, changes in the pressure of the surrounding medium through time, through an organ such

More information

ω 0 = 2π/T 0 is called the fundamental angular frequency and ω 2 = 2ω 0 is called the

ω 0 = 2π/T 0 is called the fundamental angular frequency and ω 2 = 2ω 0 is called the he ime-frequency Concept []. Review of Fourier Series Consider the following set of time functions {3A sin t, A sin t}. We can represent these functions in different ways by plotting the amplitude versus

More information

3D Sound Synthesis using the Head Related Transfer Function

3D Sound Synthesis using the Head Related Transfer Function University of Tennessee, Knoxville Trace: Tennessee Research and Creative Exchange Masters Theses Graduate School 12-2004 3D Sound Synthesis using the Head Related Transfer Function Dayu Yang University

More information

Seminar 5. Biophysics of the senses

Seminar 5. Biophysics of the senses Seminar 5 Biophysics of the senses Vision. The eye as a compound lens. Limitations of visual acuity. Colour vision. Imperfect human vision. Correction of vision. Audition. Physics of sound waves. Parameters

More information

An Analytic Model of the Cochlea and Functional Interpretations. Samiya A Alkhairy

An Analytic Model of the Cochlea and Functional Interpretations. Samiya A Alkhairy An Analytic Model of the Cochlea and Functional Interpretations by Samiya A Alkhairy Submitted to the Harvard-MIT Program in Health Sciences and Technology in partial fulfillment of the requirements for

More information

Signals, Instruments, and Systems W5. Introduction to Signal Processing Sampling, Reconstruction, and Filters

Signals, Instruments, and Systems W5. Introduction to Signal Processing Sampling, Reconstruction, and Filters Signals, Instruments, and Systems W5 Introduction to Signal Processing Sampling, Reconstruction, and Filters Acknowledgments Recapitulation of Key Concepts from the Last Lecture Dirac delta function (

More information

Biosciences in the 21st century

Biosciences in the 21st century Biosciences in the 21st century Lecture 1: Neurons, Synapses, and Signaling Dr. Michael Burger Outline: 1. Why neuroscience? 2. The neuron 3. Action potentials 4. Synapses 5. Organization of the nervous

More information

Sensing and acting Bats use sonar to detect their prey. Both bats and moths have complex sensory systems that facilitate their survival

Sensing and acting Bats use sonar to detect their prey. Both bats and moths have complex sensory systems that facilitate their survival Sensing and acting Bats use sonar to detect their prey Moths, a common prey for bats can detect the bat s sonar and attempt to flee Both bats and moths have complex sensory systems that facilitate their

More information

TEL AVIV UNIVERSITY. THE IBY AND ALADAR FLEISCHMAN FACULTY OF ENGINEERING The Zandman-Slaner Graduate School of Engineering STOCHASTIC PROPERTIES OF

TEL AVIV UNIVERSITY. THE IBY AND ALADAR FLEISCHMAN FACULTY OF ENGINEERING The Zandman-Slaner Graduate School of Engineering STOCHASTIC PROPERTIES OF TEL AVIV UNIVERSITY THE IBY AND ALADAR FLEISCHMAN FACULTY OF ENGINEERING The Zandman-Slaner Graduate School of Engineering STOCHASTIC PROPERTIES OF COINCIDENCE-DETECTOR NEURAL CELLS AND THEIR IMPLICATIONS

More information

Copyright Warning & Restrictions

Copyright Warning & Restrictions Copyright Warning & Restrictions The copyright law of the United States (Title 17, United States Code) governs the making of photocopies or other reproductions of copyrighted material. Under certain conditions

More information

Experimental Guided Spherical. Harmonics based Head-Related. Transfer Function Modeling

Experimental Guided Spherical. Harmonics based Head-Related. Transfer Function Modeling Experimental Guided Spherical Harmonics based Head-Related Transfer Function Modeling Mengqiu Zhang M.E. (Xidian University, Xi an, Shanxi, China) B.E. (China Jiliang University, Hangzhou, Zhejiang, China)

More information

Discrete Fourier Transform

Discrete Fourier Transform Last lecture I introduced the idea that any function defined on x 0,..., N 1 could be written a sum of sines and cosines. There are two different reasons why this is useful. The first is a general one,

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

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

2A1H Time-Frequency Analysis II Bugs/queries to HT 2011 For hints and answers visit dwm/courses/2tf

2A1H Time-Frequency Analysis II Bugs/queries to HT 2011 For hints and answers visit   dwm/courses/2tf Time-Frequency Analysis II (HT 20) 2AH 2AH Time-Frequency Analysis II Bugs/queries to david.murray@eng.ox.ac.uk HT 20 For hints and answers visit www.robots.ox.ac.uk/ dwm/courses/2tf David Murray. A periodic

More information

Cutting Edge in Spatial Audio

Cutting Edge in Spatial Audio EAA Euroregio - EAA Winter School 18-21 March 2013 in Merano AIA-DAGA 2013 Conference on Acoustics including the 40th Italian (AIA) Annual Conference on Acoustics and the 39th German Annual Conference

More information

Encoding and Enhancing Acoustic Information at the First Stages of the Auditory System

Encoding and Enhancing Acoustic Information at the First Stages of the Auditory System UNIVERSITY OF PENNSYLVANIA Encoding and Enhancing Acoustic Information at the First Stages of the Auditory System A thesis proposal submitted in partial satisfaction of the requirement for the degree Doctor

More information

Being a Physicist Unit 5. Summary Sheets. Gleniffer High School

Being a Physicist Unit 5. Summary Sheets. Gleniffer High School Being a Physicist Unit 5 Summary Sheets Gleniffer High School 0 Experiences & Outcomes I can explain how sound vibrations are carried by waves through air, water and other materials SCN 2-11a By exploring

More information

Hair Cells: The Sensory Transducers of the Inner Ear

Hair Cells: The Sensory Transducers of the Inner Ear Chapter 1 Hair Cells: The Sensory Transducers of the Inner Ear Hair cells are specialized cells that transform a mechanical motion into changes in membrane potential. Such changes, whereby one form of

More information

Year 7 Physics Learning Cycle 4 Overview Learning Cycle Overview: How has our understanding of waves changed the life of the average English person?

Year 7 Physics Learning Cycle 4 Overview Learning Cycle Overview: How has our understanding of waves changed the life of the average English person? Year 7 Physics Learning Cycle 4 Overview Learning Cycle Overview: How has our understanding of waves changed the life of the average English person? Line of enquiry 1: Hypothesis 1: Hypothesis 2: Hypothesis

More information

Sound 2: frequency analysis

Sound 2: frequency analysis COMP 546 Lecture 19 Sound 2: frequency analysis Tues. March 27, 2018 1 Speed of Sound Sound travels at about 340 m/s, or 34 cm/ ms. (This depends on temperature and other factors) 2 Wave equation Pressure

More information

A Source Localization/Separation/Respatialization System Based on Unsupervised Classification of Interaural Cues

A Source Localization/Separation/Respatialization System Based on Unsupervised Classification of Interaural Cues A Source Localization/Separation/Respatialization System Based on Unsupervised Classification of Interaural Cues Joan Mouba and Sylvain Marchand SCRIME LaBRI, University of Bordeaux 1 firstname.name@labri.fr

More information

TECHNICAL RESEARCH REPORT

TECHNICAL RESEARCH REPORT TECHNICAL RESEARCH REPORT Dynamics of Neural Responses in Ferret Primary Auditory Cortex: I. Spectro-Temporal Response Field Characterization by Dynamic Ripple Spectra by D.A. Depireux, J.Z. Simon, D.J.

More information

Waves Junior Science. Easy to read Version

Waves Junior Science. Easy to read Version Waves Junior Science Easy to read Version 1a Waves transfer energy Waves are a means of transferring energy from one place to another without also transferring matter. Some waves need a medium (matter)

More information

Performance analysis of the Ormia ochracea s coupled ears

Performance analysis of the Ormia ochracea s coupled ears Performance analysis of the Ormia ochracea s coupled ears Murat Akcakaya a and Arye Nehorai b Department of Electrical and Systems Engineering, Washington University in St. Louis, St. Louis, Missouri 6313

More information

Information Processing in the Interaural Time Difference Pathway of the Barn Owl

Information Processing in the Interaural Time Difference Pathway of the Barn Owl Information Processing in the Interaural Time Difference Pathway of the Barn Owl Thesis by Gestur Björn Christianson In Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy California

More information

HEARING DISTANCE: A LOW-COST MODEL FOR NEAR-FIELD BINAURAL EFFECTS

HEARING DISTANCE: A LOW-COST MODEL FOR NEAR-FIELD BINAURAL EFFECTS th European Signal Processing Conference (EUSIPCO 12) Bucharest, Romania, August 27-31, 12 HEARING DISTANCE: A LOW-COST MODEL FOR NEAR-FIELD BINAURAL EFFECTS Simone Spagnol IUAV - University of Venice

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

The development of a whole-head human finiteelement model for simulation of the transmission of bone-conducted sound

The development of a whole-head human finiteelement model for simulation of the transmission of bone-conducted sound The development of a whole-head human finiteelement model for simulation of the transmission of bone-conducted sound You Chang, Namkeun Kim and Stefan Stenfelt Journal Article Original Publication: N.B.:

More information

Tracking of Spread Spectrum Signals

Tracking of Spread Spectrum Signals Chapter 7 Tracking of Spread Spectrum Signals 7. Introduction As discussed in the last chapter, there are two parts to the synchronization process. The first stage is often termed acquisition and typically

More information

the ability to do work or cause change (work is force exerted on an object causing it to move a distance)

the ability to do work or cause change (work is force exerted on an object causing it to move a distance) Vocabulary Terms - Energy energy the ability to do work or cause change (work is force exerted on an object causing it to move a distance) heat Heat is a form of energy that flows between two substances

More information

DEA 350: HUMAN FACTORS: THE AMBIENT ENVIRONMENT (Spring 2008)

DEA 350: HUMAN FACTORS: THE AMBIENT ENVIRONMENT (Spring 2008) DEA 350: HUMAN FACTORS: THE AMBIENT ENVIRONMENT (Spring 2008) NAME: HOMEWORK II Due Date: 24 th April in class (NOTE: You will need to consult the readings as well as your class notes to complete the homework.

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

Exercises. Chapter 1. of τ approx that produces the most accurate estimate for this firing pattern.

Exercises. Chapter 1. of τ approx that produces the most accurate estimate for this firing pattern. 1 Exercises Chapter 1 1. Generate spike sequences with a constant firing rate r 0 using a Poisson spike generator. Then, add a refractory period to the model by allowing the firing rate r(t) to depend

More information

Micro-Flow in a bundle of micro-pillars. A. Keißner, Ch. Brücker

Micro-Flow in a bundle of micro-pillars. A. Keißner, Ch. Brücker Micro-Flow in a bundle of micro-pillars A. Keißner, Ch. Brücker Institute of Mechanics and Fluid Dynamics, University of Freiberg, TU Freiberg, Germany, Email: armin.keissner@imfd.tu-freiberg.de Abstract

More information

Will Penny. 21st April The Macroscopic Brain. Will Penny. Cortical Unit. Spectral Responses. Macroscopic Models. Steady-State Responses

Will Penny. 21st April The Macroscopic Brain. Will Penny. Cortical Unit. Spectral Responses. Macroscopic Models. Steady-State Responses The The 21st April 2011 Jansen and Rit (1995), building on the work of Lopes Da Sliva and others, developed a biologically inspired model of EEG activity. It was originally developed to explain alpha activity

More information

Physics Worksheet Sound and Light Section: Name:

Physics Worksheet Sound and Light Section: Name: Do Now: What is common between sound and light? What are the differences between sound and light? Sound Waves 1. Define Sound Waves from three different perspectives 6. Speed of sound =. Formula: v =.

More information

is a What you Hear The Pressure Wave sets the Ear Drum into Vibration.

is a What you Hear The Pressure Wave sets the Ear Drum into Vibration. is a What you Hear The ear converts sound energy to mechanical energy to a nerve impulse which is transmitted to the brain. The Pressure Wave sets the Ear Drum into Vibration. electroencephalogram v S

More information