Reverberation Impulse Response Analysis
|
|
- Annabella Joseph
- 5 years ago
- Views:
Transcription
1 MUS424: Signal Processing Techniques for Digital Audio Effects Handout #18 Jonathan Abel, David Berners April 27, 24 Lecture #9: April 27, 24 Lecture Notes 7 Reverberation Impulse Response Analysis 1 CCRMA Lobby Impulse Response.8.6 direct path.4 early reflections late-field reverberation time - milliseconds Figure 1: CCRMA Lobby Impulse Response showing direct path, early reflections, and latefield reverberation components.
2 2 MUS424: Handout #18 Reverberation Impulse Response Components Direct Path. The direct path is the arrival of source signals along a straight-line path from the source. The arrival time and amplitudeare fixed by the sourcelistener distance; filtering dependent on the source direction and orientation are also present. Early Reflections. Typically a small number of specular reflections from environment surfaces will arrive at the listener well separated in time or amplitude from other reflected energy. These so-called early reflections convey a sense of the environment geometry and size. Late-Field Reverberation. After a period of time, source signals have interacted sufficiently with room objects and surfaces that the density of arrivals is great enough to be indistinguishible from noise. This late-field reverberation is further characterized by a frequency-dependent exponential decay determined by the environment size and materials present.
3 MUS424: Handout #18 3 Sabine Theory Space with volume V : perfectly reflecting surfaces non-absorbing air Window or area A. Consider what happens after the energy is thoroughly mixed, and the energy density w(t, x), initially a function of time t and position x, is only a function of time. What is w(t)?
4 4 MUS424: Handout #18 Room Mixing Time To generate m reflections, roughly speaking, source signals must propagate long enough that a sphere with radius equal to the propagation distance contains roughly m times the room volume, 4 3 π r3 mix = mv. If m reflections are needed to thoroughly mix a room, then the mixing time τ mix is roughly τ mix = r mix /c, r mix = ( 3 4π mv )1/3, where c is the speed of sound in air. For m = 1, typical dimension V 1/3 = 1 meters, the mixing time is roughly τ mix 2 milliseconds.
5 MUS424: Handout #18 5 Energy Change during t Compute the change in total energy during a small interval of time, t, assuming the energy density is constant throughout the volume V. The total energy at time t+ t is the energy density scaled by the volume, V w(t + t). The total energy at time t is similarly V w(t). The difference is what leaves through the aperture, V w(t + t) = V w(t) γ(w(t),a, t).
6 6 MUS424: Handout #18 Computing γ: The Two-Dimensional Case. Consider the 2D case. A small volume of energy a distance d from the aperture, and traveling at angle ϕ relative to the aperture normal, will escape during t if its direction of travel takes it through the aperture and d c t cosϕ, where c is the speed of sound. The total energy leaving along ϕ is then w(t) A c t cosϕ
7 MUS424: Handout #18 7 Computing γ: The Two-Dimensional Case. Under the assumption that the room energy is thoroughly mixed, the probability that a given bit of energy is traveling along ϕ is (ϕ) = 1 2π. The energy lost during t is γ(w(t),a, t) = π/2 π/2 = 1 w(t) A c. π w(t) A c t cosϕ (ϕ)dϕ
8 8 MUS424: Handout #18 Computing γ: Circular Aperture. The energy lost during t traveling along direction (ϕ, θ) is w(t) A c t cosϕ. As a result, the total energy lost during t is γ(w(t), A, t) = π/2 2π w(t) A c t cosϕ (ϕ,θ)dθdϕ, where (ϕ, θ) is the probability that any bit of energy is traveling along direction (ϕ, θ).
9 MUS424: Handout #18 9 Computing γ: Circular Aperture. Assuming all directions on the sphere are equally likely, (ϕ,θ) = sinϕ 4π. The energy lost during t is therefore γ(w(t), A, t) w(t) A c t = 4π w(t) A c t = 2 w(t) A c t =. 4 π/2 2π π/2 cosϕ sinϕdθ dϕ, cosϕ sin ϕdϕ,
10 1 MUS424: Handout # w_ exp(-t / \tau) time / decay time Energy Density Behavior. By conservation of energy, V w(t + t) = V w(t) gcaw(t) t, where g is a geometric factor. Rearranging terms, w(t + t) w(t) = gca t V w(t), and taking the limit t, we have dw dt = 1 τ w(t), τ = V gca, where τ is the so-called characteristic decay time. For a well mixed room with w(t = ) = w, w(t) = w e t/τ.
11 MUS424: Handout # log_{1}(w_ exp(2 ln(.1) t / T_{6}) time / T_{6} Energy Density Behavior. The energy density follows an exponential decay to zero, w(t) = w e t/τ, τ = V gca. The time constant is proportional to room volume, and inversely proportional to window area. The reverberation time or T 6 is the time it takes the energy density to decay to 6 db below its initial value, We have w(t) = w e t/τ = w e 2ln(.1)t/T 6. T 6 = 2ln(.1) V gca (.161 S/m)V A.
12 12 MUS424: Handout #18 Multiple Windows. In the presence of multiple noninteracting windows, the energy losses add, V w(t + t) = V w(t) γ(w(t), t,a 1,...,A N ), and γ = ( A i )g c w(t). i The reverberation time is accordingly reduced, V T 6 (A i,i = 1,...,N) = 2ln(.1) gc i A. i
13 MUS424: Handout #18 13 Materials Patches vs. Windows. Absorption Coefficients, S material frequency marble brick concrete block plywood cork glass window drapery carpet hardwood grass Room surfaces typically absorb a portion of impinging acoustic energy. The portion depends on the surface material and is a function of frequency. To compute the reverberation time associated with a patch of material, its area a is scaled by its absorbing power S, A effective = a S. Note that S = 1 for an open window, and S = for a perfectly reflecting wall. The reverberation time is then V T 6 = 2 ln(.1) gc i a. is i
14 14 MUS424: Handout # T_{6} - 6-dB decay time, various low-frequency absorption settings frequency - khz Frequency-Dependent T 6. The EMT14 Plate Reverberator provides control over low-frequency reverberation time.
15 MUS424: Handout # ANSI standard air absorption at 25 degrees C and 5% relative humidity frequency - khz 1 1 Air Absorption. Air absorbs sound at a frequency-dependent rate proportional to the distance traveled. During a time interval t, every unit volume of energy density experiences the same attenuation, exp{ α(ω) t}.
16 16 MUS424: Handout #18 Air Absorption. In a room of volume V with uniform energy density w(t), the total energy absorbed via air absorption during the time interval t is γ air = V w(t) [1 exp{ α(ω) t}], V w(t) α(ω) t, assuming t 1. Again, by conservation of energy, V w(t + t) = V w(t) γ surfaces γ air, and, substituting for the absorbed energy terms, ( ) V w(t+ t) = V w(t) gc i A i + V α(ω) Rearranging, and taking the interval t to zero, ( dw gc dt = i A ) i + α(ω) w(t). V Again, we have an exponential decay, this time with characteristic decay time τ = V gc i A i + V α(ω). w(t) t.
17 MUS424: Handout #18 17 Reverberation Time as a Function of Room Size. Consider a room with volume V, absorbing materials with effective areas A i, and characteristic dimension l = V 1/3. The reverberation time is 1 T 6 (ω) = 2 ln(.1) gc i A i/v + α(ω), 1 = 2 ln(.1) gcσ(ω)/l + α(ω), where the σ(ω) is the sum of effective areas, normalized by l 2. Note that the reverberation time is the harmonic mean of the materials and air reverberation times, [ ] 1 1 T 6 (ω) = T 6 surfaces (ω) + 1 T 6 air (ω). This means that in any frequency band, the shorter reverberation times dominate. Also, air absorption provides an upper limit on reverberation time as a function of frequency. As the room gets larger, the contribution of the air becomes more important: Larger rooms will have longer reverberation times with darker tails.
18 18 MUS424: Handout #18.5 impulse response response energy profile, 25-msec. frames time - milliseconds Impulse Response Energy Envelope. Reverberation impulse responses are nonstationary noise processes with exponentially decaying variances. How do you estimate the variance of a noise process? You average sample variances: For an averaging width β, the smoothed energy envelope of an impulse response h(t) is P(t;β) = 1 β t+β/2 n=t β/2 h 2 (n).
19 MUS424: Handout #18 19 Impulse Response Energy Envelope. In computing the energy envelope, P(t; β) = 1 β t+β/2 n=t β/2 h 2 (n), it is important to pick the averaging window to be wide enough to suppress noise, but not so wide that the estimate is biased. For a late-field decay, E { h 2 (t) } = w exp( t/τ). In a neighborhood of t, t [t β/2,t + β/2] we have exp( t/τ) exp( t /τ) (1 t/τ). provided β τ. Accordingly, when β τ. E{P(t;β)} = w exp( t/τ)
20 2 MUS424: Handout #18 response power spectrum, 46-msec. interval between frames frequency - khz 1 1 Energy Decay Relief frequency - khz Impulse Response Time-Frequency Analysis. STFT, P(ω,t). Energy Decay Relief (EDR), EDR(ω,t) = Loudness Spectrogram. t P(ω,τ)dτ.
21 MUS424: Handout #18 21 Impulse Response Time-Frequency Analysis. response spectra, 25-msec. frames. power - db frequency - Bark
22 22 MUS424: Handout #18 Impulse Response Spectrogram Model. response power spectrum, Bark-spaced frequencies time - milliseconds The smoothed measured late-field spectrogram is modeled as a nonstationary noise process with mean E { PH (ω,t) } = q(ω) 2 exp( 2t/τ(ω)) + P N (ω), where q(ω) is the magnitude of the initial late-field equalization, τ(ω) is the late-field decay rate, and P N (ω) is the measurement noise spectrum.
23 MUS424: Handout #18 23 Equalization and Reverberation Time Estimation. measured, modeled response energy profile time - milliseconds Estimate q(ω), τ(ω) separately for each frequency. (These estimates form sufficient statistics if they can be accurately estimated.) Note that the db initial equalization and reverberation time are linear in the db smoothed power spectrum. P H (ω,t) = q(ω) 2 exp( 2t/τ(ω)) + ν(ω, t), Denotingby P H (ω,t) the model q(ω) 2 exp( 2t/τ(ω)), η H (ω,t) = 1log 1 {P H (ω,t) [1 + ν(ω,t)/p H (ω,t)]}. Assuming small measurement noise, η H (ω,t) 1log 1 P H (ω,t) + ν(ω, t)/p H (ω, t).
24 24 MUS424: Handout #18 Equalization and Reverberation Time Estimation. Denote by θ the unknown equalization and decay rate, [ ] ln q(ω) θ =, 1/τ(ω) and by η θ the stack of hypothesized band powers for a given equalization and decay rate, ln q(ω) t /τ(ω) η θ =.. ln q(ω) t N /τ(ω) Given noisy measurements, we have the following equation η(ω,t i ) = η θ (ω,t i ) + ɛ(ω,t i ) at a set of times t i, choose the equalization magnitude and decay rate at the ones minimizing the sum of square equation errors, ˆθ = argmin θ J(θ), J(θ) = [η θ η] [η θ η].
25 MUS424: Handout #18 25 Equalization and Reverberation Time Estimation. measured, modeled response energy profile time - milliseconds Denoting by B the basis B = [ 1 t ], the equalization and decay rate estimate is ˆθ = (B B) 1 B η, with the estimated db energy as a function of time being the projection onto B of the measured db band power as a function of time, ˆη = Bˆθ = B(B B) 1 B η.
26 26 MUS424: Handout #18 Estimate Statistics. measured, modeled response energy profile time - milliseconds Assuming the model is valid, and the equation error zero mean, the estimate is unbiased with variance } { Var {ˆθ = (B B) 1 B E [ η η ˆθ ][ η η ˆθ ] } B(B B), = σ 2 (B B) 1, in the case of an i.i.d. equation error column, { E [ η η ˆθ ][ η η ˆθ ] } = σ 2 I.
27 MUS424: Handout #18 27 Impulse Response Histrograms. 25 early impulse response histogram response amplitude late impulse response histogram response amplitude Any segment of the late field is approximately colored Gaussian noise.
28 28 MUS424: Handout #18 The Central Limit Theorem. The sum of a sufficient number of i.i.d. random variables is approximately Gaussian. Rationale: The sum of two independent random variables is the convolution of their probability densities. Defining the moment generating function µ(t) of a probability density (x) as the Fourier Transform of the probability density, the sum of n i.i.d. random variables has its moment generating function µ n (t). Because the probability density is real, positive and integrates to one, the moment generating function is approximately µ(t) (1 2σ 2 t 2 ) in a small neighborhood about t =. Raising µ(t) to the Nth power, we see that it and therefore its generating density is approximately Gaussian, µ(t) n (1 (2nσ 2 t 2 )/n) n exp{ 2nσ 2 t 2 }.
29 MUS424: Handout #18 29 Echo Density Profile. 1 impulse response time - milliseconds echo density profile, 2-msec. frames time - msec. A simple measure of impulse response echo density as a function of time is the percentage of samples in a frame which are outside a standard deviation from the mean. In the late field, this measure is expected to be that of a Gaussian, about 3%. In the presence of early reflections, the measure is expected to be lower, as the energy (and therefore standard deviation) is concentrated in the reflections.
Last time: small acoustics
Last time: small acoustics Voice, many instruments, modeled by tubes Traveling waves in both directions yield standing waves Standing waves correspond to resonances Variations from the idealization give
More informationPeaking and Shelving Filter Properties
MUS424: Signal Processing Techniques for Digital Audio Effects Handout #22 Jonathan Abel, David Berners May 18, 24 Lecture #14: May 18, 24 Lecture Notes 12 Peaking and Shelving Filter Properties Shelf
More informationConvention Paper Presented at the 121st Convention 2006 October 5 8 San Francisco, CA, USA
Audio Engineering Society Convention Paper Presented at the st Convention 6 October 5 8 San Francisco, CA, USA This convention paper has been reproduced from the author s advance manuscript, without editing,
More informationNoise in enclosed spaces. Phil Joseph
Noise in enclosed spaces Phil Joseph MODES OF A CLOSED PIPE A 1 A x = 0 x = L Consider a pipe with a rigid termination at x = 0 and x = L. The particle velocity must be zero at both ends. Acoustic resonances
More informationEE1.el3 (EEE1023): Electronics III. Acoustics lecture 18 Room acoustics. Dr Philip Jackson.
EE1.el3 EEE1023): Electronics III Acoustics lecture 18 Room acoustics Dr Philip Jackson www.ee.surrey.ac.uk/teaching/courses/ee1.el3 Room acoustics Objectives: relate sound energy density to rate of absorption
More informationECE 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 informationAnalysis and synthesis of room reverberation based on a statistical time-frequency model
Analysis and synthesis of room reverberation based on a statistical time-frequency model Jean-Marc Jot, Laurent Cerveau, Olivier Warusfel IRCAM. 1 place Igor-Stravinsky. F-75004 Paris, France. Tel: (+33)
More informationPhysics 115 Lecture 20. Reflection and Reverberation March 9, 2018
Physics 115 Lecture 20 Reflection and Reverberation March 9, 2018 1 Attenuation by absorption The wave intensity decreases due to internal friction among the molecules of the medium Attenuation is described
More informationCepstral Deconvolution Method for Measurement of Absorption and Scattering Coefficients of Materials
Cepstral Deconvolution Method for Measurement of Absorption and Scattering Coefficients of Materials Mehmet ÇALIŞKAN a) Middle East Technical University, Department of Mechanical Engineering, Ankara, 06800,
More informationIntroduction 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 informationMobile Radio Communications
Course 3: Radio wave propagation Session 3, page 1 Propagation mechanisms free space propagation reflection diffraction scattering LARGE SCALE: average attenuation SMALL SCALE: short-term variations in
More informationCMPT 889: Lecture 2 Sinusoids, Complex Exponentials, Spectrum Representation
CMPT 889: Lecture 2 Sinusoids, Complex Exponentials, Spectrum Representation Tamara Smyth, tamaras@cs.sfu.ca School of Computing Science, Simon Fraser University September 26, 2005 1 Sinusoids Sinusoids
More informationSound, acoustics Slides based on: Rossing, The science of sound, 1990, and Pulkki, Karjalainen, Communication acoutics, 2015
Acoustics 1 Sound, acoustics Slides based on: Rossing, The science of sound, 1990, and Pulkki, Karjalainen, Communication acoutics, 2015 Contents: 1. Introduction 2. Vibrating systems 3. Waves 4. Resonance
More informationThe Hilbert Transform
The Hilbert Transform David Hilbert 1 ABSTRACT: In this presentation, the basic theoretical background of the Hilbert Transform is introduced. Using this transform, normal real-valued time domain functions
More informationSinusoids. Amplitude and Magnitude. Phase and Period. CMPT 889: Lecture 2 Sinusoids, Complex Exponentials, Spectrum Representation
Sinusoids CMPT 889: Lecture Sinusoids, Complex Exponentials, Spectrum Representation Tamara Smyth, tamaras@cs.sfu.ca School of Computing Science, Simon Fraser University September 6, 005 Sinusoids are
More informationunit 4 acoustics & ultrasonics
unit 4 acoustics & ultrasonics acoustics ACOUSTICS Deals with the production, propagation and detection of sound waves Classification of sound: (i) Infrasonic 20 Hz (Inaudible) (ii) Audible 20 to 20,000Hz
More informationSound radiation and sound insulation
11.1 Sound radiation and sound insulation We actually do not need this chapter You have learned everything you need to know: When waves propagating from one medium to the next it is the change of impedance
More informationTheoretical Seismology. Astrophysics and Cosmology and Earth and Environmental Physics. FTAN Analysis. Fabio ROMANELLI
Theoretical Seismology Astrophysics and Cosmology and Earth and Environmental Physics FTAN Analysis Fabio ROMANELLI Department of Mathematics & Geosciences University of Trieste romanel@units.it 1 FTAN
More informationRecent topics in acoustic scattering coefficient determination for wall surfaces
Toronto, Canada International Symposium on Room Acoustics 2013 June 9-11 Recent topics in acoustic scattering coefficient determination for wall surfaces Tetsuya Sakuma (sakuma@k.u-tokyo.ac.jp) Hyojin
More informationCMPT 318: Lecture 5 Complex Exponentials, Spectrum Representation
CMPT 318: Lecture 5 Complex Exponentials, Spectrum Representation Tamara Smyth, tamaras@cs.sfu.ca School of Computing Science, Simon Fraser University January 23, 2006 1 Exponentials The exponential is
More informationBasic Study on a Laboratory Measurement Method of the Normal-Incidence Scattering Coefficient
Basic Study on a Laboratory Measurement Method of the Normal-Incidence Scattering Coefficient Hyojin LEE, Tetsuya SAKUMA Graduate School of Frontier Sciences, University of Tokyo, Kashiwa, Japan. 1. Introduction
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 informationOptimal Speech Enhancement Under Signal Presence Uncertainty Using Log-Spectral Amplitude Estimator
1 Optimal Speech Enhancement Under Signal Presence Uncertainty Using Log-Spectral Amplitude Estimator Israel Cohen Lamar Signal Processing Ltd. P.O.Box 573, Yokneam Ilit 20692, Israel E-mail: icohen@lamar.co.il
More informationMultiple Random Variables
Multiple Random Variables Joint Probability Density Let X and Y be two random variables. Their joint distribution function is F ( XY x, y) P X x Y y. F XY ( ) 1, < x
More informationDigital Band-pass Modulation PROF. MICHAEL TSAI 2011/11/10
Digital Band-pass Modulation PROF. MICHAEL TSAI 211/11/1 Band-pass Signal Representation a t g t General form: 2πf c t + φ t g t = a t cos 2πf c t + φ t Envelope Phase Envelope is always non-negative,
More informationL6: Short-time Fourier analysis and synthesis
L6: Short-time Fourier analysis and synthesis Overview Analysis: Fourier-transform view Analysis: filtering view Synthesis: filter bank summation (FBS) method Synthesis: overlap-add (OLA) method STFT magnitude
More informationDIRECTION ESTIMATION BASED ON SOUND INTENSITY VECTORS. Sakari Tervo
7th European Signal Processing Conference (EUSIPCO 9) Glasgow, Scotland, August 4-8, 9 DIRECTION ESTIMATION BASED ON SOUND INTENSITY VECTORS Sakari Tervo Helsinki University of Technology Department of
More information13.1 Ion Acoustic Soliton and Shock Wave
13 Nonlinear Waves In linear theory, the wave amplitude is assumed to be sufficiently small to ignore contributions of terms of second order and higher (ie, nonlinear terms) in wave amplitude In such a
More informationCEPSTRAL ANALYSIS SYNTHESIS ON THE MEL FREQUENCY SCALE, AND AN ADAPTATIVE ALGORITHM FOR IT.
CEPSTRAL ANALYSIS SYNTHESIS ON THE EL FREQUENCY SCALE, AND AN ADAPTATIVE ALGORITH FOR IT. Summarized overview of the IEEE-publicated papers Cepstral analysis synthesis on the mel frequency scale by Satochi
More informationRevision of Lecture 4
Revision of Lecture 4 We have completed studying digital sources from information theory viewpoint We have learnt all fundamental principles for source coding, provided by information theory Practical
More informationChapter 17: Waves II. Sound waves are one example of Longitudinal Waves. Sound waves are pressure waves: Oscillations in air pressure and air density
Sound waves are one example of Longitudinal Waves Sound waves are pressure waves: Oscillations in air pressure and air density Before we can understand pressure waves in detail, we need to understand what
More informationJ10M.1 - Rod on a Rail (M93M.2)
Part I - Mechanics J10M.1 - Rod on a Rail (M93M.2) J10M.1 - Rod on a Rail (M93M.2) s α l θ g z x A uniform rod of length l and mass m moves in the x-z plane. One end of the rod is suspended from a straight
More informationChapter 5 Frequency Domain Analysis of Systems
Chapter 5 Frequency Domain Analysis of Systems CT, LTI Systems Consider the following CT LTI system: xt () ht () yt () Assumption: the impulse response h(t) is absolutely integrable, i.e., ht ( ) dt< (this
More informationSIMPLE APPLICATION OF STI-METHOD IN PREDICTING SPEECH TRANSMISSION IN CLASSROOMS. Jukka Keränen, Petra Larm, Valtteri Hongisto
SIMPLE APPLICATION OF STI-METHOD IN PREDICTING SPEECH TRANSMISSION IN CLASSROOMS Jukka Keränen, Petra Larm, Valtteri Hongisto Finnish Institute of Occupational Health Laboratory of Ventilation and Acoustics
More informationECE6604 PERSONAL & MOBILE COMMUNICATIONS. Week 3. Flat Fading Channels Envelope Distribution Autocorrelation of a Random Process
1 ECE6604 PERSONAL & MOBILE COMMUNICATIONS Week 3 Flat Fading Channels Envelope Distribution Autocorrelation of a Random Process 2 Multipath-Fading Mechanism local scatterers mobile subscriber base station
More informationNoise Robust Isolated Words Recognition Problem Solving Based on Simultaneous Perturbation Stochastic Approximation Algorithm
EngOpt 2008 - International Conference on Engineering Optimization Rio de Janeiro, Brazil, 0-05 June 2008. Noise Robust Isolated Words Recognition Problem Solving Based on Simultaneous Perturbation Stochastic
More informationStochastic Structural Dynamics Prof. Dr. C. S. Manohar Department of Civil Engineering Indian Institute of Science, Bangalore
Stochastic Structural Dynamics Prof. Dr. C. S. Manohar Department of Civil Engineering Indian Institute of Science, Bangalore Lecture No. # 33 Probabilistic methods in earthquake engineering-2 So, we have
More informationFrequency- and Time-Domain Spectroscopy
Frequency- and Time-Domain Spectroscopy We just showed that you could characterize a system by taking an absorption spectrum. We select a frequency component using a grating or prism, irradiate the sample,
More informationSinusoidal Modeling. Yannis Stylianou SPCC University of Crete, Computer Science Dept., Greece,
Sinusoidal Modeling Yannis Stylianou University of Crete, Computer Science Dept., Greece, yannis@csd.uoc.gr SPCC 2016 1 Speech Production 2 Modulators 3 Sinusoidal Modeling Sinusoidal Models Voiced Speech
More informationChapter 3 Room acoustics
Chapter 3 Room acoustics Acoustic phenomena in a closed room are the main focal problems in architectural acoustics. In this chapter, fundamental theories are described and their practical applications
More informationBME 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 informationBasic principles of the seismic method
Chapter 2 Basic principles of the seismic method In this chapter we introduce the basic notion of seismic waves. In the earth, seismic waves can propagate as longitudinal (P) or as shear (S) waves. For
More informationHandout 10: Inductance. Self-Inductance and inductors
1 Handout 10: Inductance Self-Inductance and inductors In Fig. 1, electric current is present in an isolate circuit, setting up magnetic field that causes a magnetic flux through the circuit itself. This
More information6.003: Signals and Systems Lecture 18 April 13, 2010
6.003: Signals and Systems Lecture 8 April 3, 200 6.003: Signals and Systems Discrete-Time Frequency Representations Signals and/or Systems Two perspectives: feedback and control (focus on systems) X +
More informationROOM ACOUSTICS THREE APPROACHES 1. GEOMETRIC RAY TRACING SOUND DISTRIBUTION
ROOM ACOUSTICS THREE APPROACHES 1. GEOMETRIC RAY TRACING. RESONANCE (STANDING WAVES) 3. GROWTH AND DECAY OF SOUND 1. GEOMETRIC RAY TRACING SIMPLE CONSTRUCTION OF SOUND RAYS WHICH OBEY THE LAWS OF REFLECTION
More informationEcho cancellation by deforming sound waves through inverse convolution R. Ay 1 ward DeywrfmzMf o/ D/g 0001, Gauteng, South Africa
Echo cancellation by deforming sound waves through inverse convolution R. Ay 1 ward DeywrfmzMf o/ D/g 0001, Gauteng, South Africa Abstract This study concerns the mathematical modelling of speech related
More informationOSE801 Engineering System Identification. Lecture 09: Computing Impulse and Frequency Response Functions
OSE801 Engineering System Identification Lecture 09: Computing Impulse and Frequency Response Functions 1 Extracting Impulse and Frequency Response Functions In the preceding sections, signal processing
More informationWireless Communications
NETW701 Wireless Communications Dr. Wassim Alexan Winter 2018 Lecture 2 NETW705 Mobile Communication Networks Dr. Wassim Alexan Winter 2018 Lecture 2 Wassim Alexan 2 Reflection When a radio wave propagating
More informationBayesian room-acoustic modal analysis
Bayesian room-acoustic modal analysis Wesley Henderson a) Jonathan Botts b) Ning Xiang c) Graduate Program in Architectural Acoustics, School of Architecture, Rensselaer Polytechnic Institute, Troy, New
More informationPart 7. Nonlinearity
Part 7 Nonlinearity Linear System Superposition, Convolution re ( ) re ( ) = r 1 1 = r re ( 1 + e) = r1 + r e excitation r = r() e response In the time domain: t rt () = et () ht () = e( τ) ht ( τ) dτ
More informationElec4621 Advanced Digital Signal Processing Chapter 11: Time-Frequency Analysis
Elec461 Advanced Digital Signal Processing Chapter 11: Time-Frequency Analysis Dr. D. S. Taubman May 3, 011 In this last chapter of your notes, we are interested in the problem of nding the instantaneous
More informationProblem Sheet 1 Examples of Random Processes
RANDOM'PROCESSES'AND'TIME'SERIES'ANALYSIS.'PART'II:'RANDOM'PROCESSES' '''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''Problem'Sheets' Problem Sheet 1 Examples of Random Processes 1. Give
More informationω 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 informationTest Report. RI Acoustic Lab. Measurement of Sound Absorption Coefficient for RockDelta NoiStop Noise Barrier. 20 Feb. 07
Test Report RI Acoustic Lab Measurement of Sound Absorption Coefficient for RockDelta NoiStop Noise Barrier 20 Feb. 07 Title Measurement of Sound Absorption Coefficient for RockDelta NoiStop Noise Barrier
More informationFeature extraction 2
Centre for Vision Speech & Signal Processing University of Surrey, Guildford GU2 7XH. Feature extraction 2 Dr Philip Jackson Linear prediction Perceptual linear prediction Comparison of feature methods
More informationL29: Fourier analysis
L29: Fourier analysis Introduction The discrete Fourier Transform (DFT) The DFT matrix The Fast Fourier Transform (FFT) The Short-time Fourier Transform (STFT) Fourier Descriptors CSCE 666 Pattern Analysis
More informationVer Chap Lecture 15- ECE 240a. Q-Switching. Mode Locking. ECE 240a Lasers - Fall 2017 Lecture Q-Switch Discussion
ing Ver Chap. 9.3 Lasers - Fall 2017 Lecture 15 1 ing ing (Cavity Dumping) 1 Turn-off cavity - prevent lasing 2 Pump lots of energy into upper state - use pulsed pump 3 Turn cavity back on - all the energy
More informationProceedings of Meetings on Acoustics
Proceedings of Meetings on Acoustics Volume 19, 2013 http://acousticalsociety.org/ ICA 2013 Montreal Montreal, Canada 2-7 June 2013 Signal Processing in Acoustics Session 4aSP: Sensor Array Beamforming
More informationHeriot-Watt University
Heriot-Watt University Distinctly Global www.hw.ac.uk Thermodynamics By Peter Cumber Prerequisites Interest in thermodynamics Some ability in calculus (multiple integrals) Good understanding of conduction
More informationThe Dyson equation is exact. ( Though we will need approximations for M )
Lecture 25 Monday Nov 30 A few final words on Multiple scattering Our Dyson equation is or, written out, < G(r r ') >= G o (r r ') + dx dy G o (r x)m (x y) < G(y r ') > The Dyson equation is exact. ( Though
More informationAPPENDIX B. Noise Primer
APPENDIX B Noise Primer NOISE PRIMER TABLE OF CONTENTS 1. INTRODUCTION...1 2. BASIC SOUND PRINCIPLES...1 2.1. SOUND AS A WAVE...1 2.2. SOUND PRESSURE LEVEL...2 2.2.1. Decibel Math...4 2.2.2. Descriptive
More informationIntroduction to Time-Frequency Distributions
Introduction to Time-Frequency Distributions Selin Aviyente Department of Electrical and Computer Engineering Michigan State University January 19, 2010 Motivation for time-frequency analysis When you
More informationPROPAGATION PHASE REPRESENTATION IN 3D SPACE USING POLES AND ZEROS IN COMPLEX FREQUENCY PLANE. ( address of lead author)
ICSV14 Cairns Australia 9-12 July, 2007 PROPAGATION PHASE REPRESENTATION IN 3D SPACE USING POLES AND ZEROS IN COMPLEX FREQUENCY PLANE Yoshinori Takahashi 1, Mikio Tohyama 2, and Kazunori Miyoshi 1 1 Kogakuin
More informationSignals, 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 informationInvestigating the Relationship Between Cavendish Temperature Fluctuation and Torsional Oscillation
Investigating the Relationship Between Cavendish Temperature Fluctuation and Torsional Oscillation John Grasel 4 March 2010 Abstract The Cavendish apparatus measures the gravitational attraction between
More informationLecture Notes 5: Multiresolution Analysis
Optimization-based data analysis Fall 2017 Lecture Notes 5: Multiresolution Analysis 1 Frames A frame is a generalization of an orthonormal basis. The inner products between the vectors in a frame and
More informationResearch Letter An Algorithm to Generate Representations of System Identification Errors
Research Letters in Signal Processing Volume 008, Article ID 5991, 4 pages doi:10.1155/008/5991 Research Letter An Algorithm to Generate Representations of System Identification Errors Wancheng Zhang and
More informationContinuous Time Signal Analysis: the Fourier Transform. Lathi Chapter 4
Continuous Time Signal Analysis: the Fourier Transform Lathi Chapter 4 Topics Aperiodic signal representation by the Fourier integral (CTFT) Continuous-time Fourier transform Transforms of some useful
More informationLecture 19 Optical MEMS (1)
EEL6935 Advanced MEMS (Spring 5) Instructor: Dr. Huikai Xie Lecture 19 Optical MEMS (1) Agenda: Optics Review EEL6935 Advanced MEMS 5 H. Xie 3/8/5 1 Optics Review Nature of Light Reflection and Refraction
More informationFDTD analysis on the sound insulation performance of wall system with narrow gaps
FDTD analysis on the sound insulation performance of wall system with narrow gaps Takumi Asakura a Shinichi Sakamoto b Institute of Industrial Science, The University of Tokyo. Komaba 4-6-, Meguro-ku,
More informationComputer Vision & Digital Image Processing
Computer Vision & Digital Image Processing Image Restoration and Reconstruction I Dr. D. J. Jackson Lecture 11-1 Image restoration Restoration is an objective process that attempts to recover an image
More informationPHASED-ARRAY FOCUSING WITH LONGITUDINAL GUIDED WAVES IN A VISCOELASTIC COATED HOLLOW CYLINDER
PHASED-ARRAY FOCUSING WITH LONGITUDINAL GUIDED WAVES IN A VISCOELASTIC COATED HOLLOW CYLINDER W. Luo, J. L. Rose, J. K. Van Velsor, J. Mu Department of Engineering Science & Mechanics, The Pennsylvania
More information9231 FURTHER MATHEMATICS
CAMBRIDGE INTERNATIONAL EXAMINATIONS GCE Advanced Level MARK SCHEME for the May/June 201 series 921 FURTHER MATHEMATICS 921/21 Paper 2, maximum raw mark 100 This mark scheme is published as an aid to teachers
More informationCITY UNIVERSITY LONDON. MSc in Information Engineering DIGITAL SIGNAL PROCESSING EPM746
No: CITY UNIVERSITY LONDON MSc in Information Engineering DIGITAL SIGNAL PROCESSING EPM746 Date: 19 May 2004 Time: 09:00-11:00 Attempt Three out of FIVE questions, at least One question from PART B PART
More informationLecture 3 - Design of Digital Filters
Lecture 3 - Design of Digital Filters 3.1 Simple filters In the previous lecture we considered the polynomial fit as a case example of designing a smoothing filter. The approximation to an ideal LPF can
More informationRandomly Modulated Periodic Signals
Randomly Modulated Periodic Signals Melvin J. Hinich Applied Research Laboratories University of Texas at Austin hinich@mail.la.utexas.edu www.la.utexas.edu/~hinich Rotating Cylinder Data Fluid Nonlinearity
More informationImaging with Ambient Noise II
Imaging with Ambient Noise II George Papanicolaou Stanford University Michigan State University, Department of Mathematics Richard E. Phillips Lecture Series April 21, 2009 With J. Garnier, University
More informationChapter 2. Room acoustics
Chapter 2. Room acoustics Acoustic phenomena in a closed room are the main focal problems in architectural acoustics. In this chapter, fundamental theories are described and their practical applications
More informationCompensating for attenuation by inverse Q filtering. Carlos A. Montaña Dr. Gary F. Margrave
Compensating for attenuation by inverse Q filtering Carlos A. Montaña Dr. Gary F. Margrave Motivation Assess and compare the different methods of applying inverse Q filter Use Q filter as a reference to
More informationLoudspeaker Choice and Placement. D. G. Meyer School of Electrical & Computer Engineering
Loudspeaker Choice and Placement D. G. Meyer School of Electrical & Computer Engineering Outline Sound System Design Goals Review Acoustic Environment Outdoors Acoustic Environment Indoors Loudspeaker
More informationMassachusetts Institute of Technology Department of Electrical Engineering and Computer Science : Discrete-Time Signal Processing
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.34: Discrete-Time Signal Processing OpenCourseWare 006 ecture 8 Periodogram Reading: Sections 0.6 and 0.7
More information4.2 MODELLING OF SOUND FIELDS IN ROOMS. OVERVIEW
CHAPTER 4 Room acoustics 4.1 INTRODUCTION In talking about the concept of room acoustics we shall include all aspects of the behaviour of sound in a room, covering both the physical aspects as well as
More information1) The K.E and P.E of a particle executing SHM with amplitude A will be equal to when its displacement is:
1) The K.E and P.E of a particle executing SHM with amplitude A will be equal to when its displacement is: 2) The bob of simple Pendulum is a spherical hallow ball filled with water. A plugged hole near
More informationBasic System. Basic System. Sonosite 180. Acuson Sequoia. Echo occurs at t=2z/c where c is approximately 1500 m/s or 1.5 mm/µs
Bioengineering 280A Principles of Biomedical Imaging Fall Quarter 2007 Ultrasound Lecture Sonosite 80 From Suetens 2002 Acuson Sequoia Basic System Basic System Echo occurs at t=2/c where c is approximately
More informationEFFECTS OF PERMEABILITY ON SOUND ABSORPTION AND SOUND INSULATION PERFORMANCE OF ACOUSTIC CEILING PANELS
EFFECTS OF PERMEABILITY ON SOUND ABSORPTION AND SOUND INSULATION PERFORMANCE OF ACOUSTIC CEILING PANELS Kento Hashitsume and Daiji Takahashi Graduate School of Engineering, Kyoto University email: kento.hashitsume.ku@gmail.com
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 informationA system that is both linear and time-invariant is called linear time-invariant (LTI).
The Cooper Union Department of Electrical Engineering ECE111 Signal Processing & Systems Analysis Lecture Notes: Time, Frequency & Transform Domains February 28, 2012 Signals & Systems Signals are mapped
More informationWavelet Transform. Figure 1: Non stationary signal f(t) = sin(100 t 2 ).
Wavelet Transform Andreas Wichert Department of Informatics INESC-ID / IST - University of Lisboa Portugal andreas.wichert@tecnico.ulisboa.pt September 3, 0 Short Term Fourier Transform Signals whose frequency
More informationFiltering and Edge Detection
Filtering and Edge Detection Local Neighborhoods Hard to tell anything from a single pixel Example: you see a reddish pixel. Is this the object s color? Illumination? Noise? The next step in order of complexity
More informationREAL-TIME TIME-FREQUENCY BASED BLIND SOURCE SEPARATION. Scott Rickard, Radu Balan, Justinian Rosca. Siemens Corporate Research Princeton, NJ 08540
REAL-TIME TIME-FREQUENCY BASED BLIND SOURCE SEPARATION Scott Rickard, Radu Balan, Justinian Rosca Siemens Corporate Research Princeton, NJ 84 fscott.rickard,radu.balan,justinian.roscag@scr.siemens.com
More informationElectromagnetic Waves. Chapter 33 (Halliday/Resnick/Walker, Fundamentals of Physics 8 th edition)
PH 222-3A Spring 2007 Electromagnetic Waves Lecture 22 Chapter 33 (Halliday/Resnick/Walker, Fundamentals of Physics 8 th edition) 1 Chapter 33 Electromagnetic Waves Today s information age is based almost
More informationI. Rayleigh Scattering. EE Lecture 4. II. Dipole interpretation
I. Rayleigh Scattering 1. Rayleigh scattering 2. Dipole interpretation 3. Cross sections 4. Other approximations EE 816 - Lecture 4 Rayleigh scattering is an approximation used to predict scattering from
More informationHeadphone 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 informationContents. Digital Signal Processing, Part II: Power Spectrum Estimation
Contents Digital Signal Processing, Part II: Power Spectrum Estimation 5. Application of the FFT for 7. Parametric Spectrum Est. Filtering and Spectrum Estimation 7.1 ARMA-Models 5.1 Fast Convolution 7.2
More informationMedian Filter Based Realizations of the Robust Time-Frequency Distributions
TIME-FREQUENCY SIGNAL ANALYSIS 547 Median Filter Based Realizations of the Robust Time-Frequency Distributions Igor Djurović, Vladimir Katkovnik, LJubiša Stanković Abstract Recently, somenewefficient tools
More informationACOUSTICAL MEASUREMENTS BY ADAPTIVE SYSTEM MODELING
ACOUSTICAL MEASUREMENTS BY ADAPTIVE SYSTEM MODELING PACS REFERENCE: 43.60.Qv Somek, Branko; Dadic, Martin; Fajt, Sinisa Faculty of Electrical Engineering and Computing University of Zagreb Unska 3, 10000
More information3.2 Complex Sinusoids and Frequency Response of LTI Systems
3. Introduction. A signal can be represented as a weighted superposition of complex sinusoids. x(t) or x[n]. LTI system: LTI System Output = A weighted superposition of the system response to each complex
More informationLecture 1: Pragmatic Introduction to Stochastic Differential Equations
Lecture 1: Pragmatic Introduction to Stochastic Differential Equations Simo Särkkä Aalto University, Finland (visiting at Oxford University, UK) November 13, 2013 Simo Särkkä (Aalto) Lecture 1: Pragmatic
More informationUtility of Correlation Functions
Utility of Correlation Functions 1. As a means for estimating power spectra (e.g. a correlator + WK theorem). 2. For establishing characteristic time scales in time series (width of the ACF or ACV). 3.
More information