Echo cancellation by deforming sound waves through inverse convolution R. Ay 1 ward DeywrfmzMf o/ D/g 0001, Gauteng, South Africa
|
|
- Shon Russell
- 6 years ago
- Views:
Transcription
1 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 acoustical hearing discomforts in medium to large environments (e.g. auditoriums and concert halls). Discomforts arise from both phase distortion (echo) and perceived frequency transformations (Comb filters) Borwickfl], Ando[3]. These distortions are caused by the physical constraints imposed on sound transmission within large environments. Deconvolution was applied to the measured signals to recover the original sound from the environmentally distorted sound, thus separating it from the transformation system function (Room impulse response). Several samples were used for different positions within the hall. The different methods investigated for the deconvolution process were cross-correlation, cepstmm and adaptive filter techniques, Oppcnlicim[2]. Upon separation the inverse of the system transfer function was determined and used in the latter part of the study to pre-deform the original sound. Only a part of the system response was used. The first ins constructively contribute to the intelligibility of speech and masked towards the end of the impulse response by the background noise at 40 dba, Borwickfl]. The critical part between and 2 sec was then used as the basis of the filter algorithm design. Because of the length of this part of the impulse response as well as the real time processing constraints a FIR filler could not be implemented, Ando[3, Tohyama[4]. Instead thefilterwas designed in the frequency domain. During transmission the environment now acting on the pre-deformcd (filtered) sound renders better quality speech that is easier to understand, effectively removing part of the deformation. Finally the psychoacouslics of the total system were evaluated by a panel of listeners in the auditorium (concert hall) and the system was found to be effective for some targeted areas within the particular hall. The impact of these results for addressing acoustical problems will be discussed. 1 Introduction Looking at the traditional methods of improving the acoustic properties of a room one finds a wide variety of solutions. Some designs are simple while others are rather elaborate but most often these solutions are all expensive,
2 122 Computational Acoustics and Its Environmental Applications Ando[3]. The need for this special treatment of listening rooms, lecture halls or auditoria, lies with the ease of understanding a speaker at such a venue. The first phenomenon is frequency related. Early reflections act as comb filters. The impulse response for the early reflections is given by Ando[3] as With g < I the reflections decrease exponentially. Taking the Fourier transform of equation (1) gives (0 (2) obtained by using the geometric series properties. The effect of (2) can be seen in Fig. 1. Fig.1 Comb Filter Frequency Response The second phenomenon is that of masking. Strong echoes effectively mask the direct sound. Consonants are especially prone to be masked by echoed vowels. Optimum signal to noise ratio for speech intelligibility is 25dB, which is also true for the ratio between the direct sound and the echoes or reverberation, Borwick[l].
3 Computational Acoustics and Its Environmental Applications 123 Taking the above into consideration the objective of the study was to find an inverse model for part of the impulse and pre-deform the signal before transmission. 2 Room impulse response The room impulse response was measured with a Sennheiser microphone on a standard analog to digital converter card for a PC and afterwards imported into MATLAB. The data was measured at 22 khz with 16 bit resolution. The impulse was generated with an electric spark generator from Grozier Technical Systems. The room response x,[n] of a local class room was obtained by taking the average over M=1024 measured responses, Berwick [1]. This gives signal to noise ratio improvement of 30dB (3) as each doubling in the number of measured responses used, give rise to 6 db improvement of signal and only 3 db 'improvement' of the noise. -improvement - 10* logm (3) Similarly the excitation pulse Xj[n] was measured in an anechoic room and averaged over 1024 sample signals and extended to the same length as Xr[n]. The measured responses were extended to double their original length and filled with zeros to ensure that the discrete signals are causal, Oppenheim [2], Using the properties of the cepstrum the measured room response x,[n] was deconvoluted with the aid of the measured impulse x,[n] (4) (5) (6) (7). (5) (6) = log[ A\ (w)] - log[a\ (7) x,[nj Xj[nj DFT[.] TT / \ LOG[ DFT[ ] 'V / LOG[.] log H,(co) ] ~ EXP[.] xv^v 1DFT[.] Cx[n] Fig. 2 Block diagram for cepstrum-deconvolution.
4 124 Computational Acoustics and Its Environmental Applications Using an M-file based on (Fig 2) to implement equations (5), (6) and (7) the new impulse response x%[n] was obtained. Any artifacts caused by the excitation pulse are now effectively removed from the measured room response. The newly obtained room impulse response (Fig 3) was used in the design of the inverse filter. Fig. 3 Room Impulse Response 2 3 Samples (1=1/22050) x10 3 Inverse response and modelling Enclosed spaces are beneficial to speech intelligibility. The reverberation of small rooms contribute to the total energy reaching the ears as the ears and brain integrate sounds for the P* msec. There after distinct echoes can be heard. In large rooms the early echoes and parts reverberation work destructively on the direct sound waves, Borwick[l]. As mentioned in the introduction some of the reflections are much larger than parts of the direct sound waves. The non sounding consonants are masked and speech intelligibility suffers. The critical section is therefore defined as the period after 25 msec until the response is masked at 40 dba by the background noise (Fig 3). The preferred response is obtained by cutting away the room response %a[n] after 25msec and padding the signal with zeros. The preferred response (Fig 4) is used in the design of the inverse filter to remove the effects of the critical section (Fig 3).
5 Computational Acoustics and Its Environmental Applications Fig. 4 Desired Room Impulse Response CL E Samples (1=1/22050) x10 The inverse response z[n] is obtained using the procedure outlined in Fig. 2. Using equation (8) the result can be verified. with (8) Z(co) is the inverse frequency response of the modified room frequency response XRmod(co) and is used as the required inverse filter. The inverse filter obtained is not directly usable. The zeros of the room response become the poles of the inverse response during the above mentioned process, Oppenheim[2] & Tohyama[4]. As some of the zeros fall outside the unit circle the converted poles are also outside the unit circle and cause the inverse response to be unstable and non causal. To be usable the time response has to be shifted samples. This was done by changing the slope of the phase by multiplying the phase with 2%(21028), Phillips [5] and equation (9). >G(co)e (9) This results in a time response with t = 0 in the middle of the samples instead of at the first sample. This response is a non causal time response. The new time response z[n] can be seen in Fig. 5
6 126 Computational Acoustics and Its Environmental Applications Fig. 5 Inverse Filter Impulse Response o.: Samples (1=1/22050) Implementation and results Implementing z[n] as a convolution filter is a bit unpractical as it consists of samples. This amounts to large size memory for any digital signal processor (DSP) system. As floating point arithmetic is needed for the calculations almost 1 Mbytes of ram is needed, 2 x samples x (80 bits / 8), to implement a circular buffer convolution filter. Sampling at 22 khz the filter will have a 4 second delay before the first results appear. As a real time solution this system can not be implemented. As a non causal system it is further not of practical value. Turning to the frequency domain for answers the memory and real time constraints are once more impairing the practical implementation of the filter. It was decided in the end to implement it as an off-line system. Samples of speech were recorded and manipulated off-line. The spoken words used were shorter than 2 sec. This limitation is due to the causality constraint. The speech samples Sa(co) were filtered in the frequency domain using equation (10). (10)
7 Computational Acoustics and Its Environmental Applications 127 SAmod(co) was transformed back to the time domain using the inverse discrete Fourier transform. These 'filtered' sound samples were then played back in the room where the original samples of the room impulse response were measured. The results were not as effective as originally envisaged. The reason for this is the unstable poles of the inverse filter outside the unit circle. Working with the minimum phase alone does not provide an adequate answer either. 5 Conclusion It has been shown that the echoes can be eliminated alhough the effect is localised to a single position at this time. More robust systems would involve exponetial windowing and elimination of the poles outside the unit circle. Although only limited answers could be found, the author is confident that the results obtained are both encouraging and useful. Traditional methods are still better but future work might bring about less expensive and practical solutions. Current research is directed at the further investigation of the all-pass and minimum phase responses of the inverse filter as well as windowing techniques. Another part of the study is concerned with the modelling of the inverse response by finding a simplified pole-zero descriptions for the inverse response. References 1. Borwick, J. (ed). Loudspeaker and Headphone Handbook, Butterworth, London, Boston, Singapore, Sydney, Toronto and Wellington, Oppenheim, A.V. & Schafer R W Discrete-time Signal Processing, Prentise-Hall, Engelwood Cliffs, New Jersey, Ando, Y. Concert Hall Acoustics, Springer-Verlag, Berlin, Heidelberg New York and Tokyo, Tohyama, M, Suzuki, H. & Ando, Y. The Nature and Technology of Acoustic Space, Academic Press, London, San Diego, NewYork, Boston, Sydney, Tokyo and Toronto, Phillips, C.L. & Parr J.M. Signals, Systems and Transforms, Prentice Hall, Englewood Cliffs, New Jersey, 1995.
BLIND DEREVERBERATION USING SHORT TIME CEPSTRUM FRAME SUBTRACTION. J.S. van Eeghem (l), T. Koike (2) and M.Tohyama (1)
FIFTH INTERNATIONAL CONGRESS ON SOUND DECEMBER 15-18, 1997 ADELAIDE, SOUTH AUSTRALIA AND VIBRATION BLIND DEREVERBERATION USING SHORT TIME CEPSTRUM FRAME SUBTRACTION J.S. van Eeghem (l), T. Koike (2) and
More informationDISCRETE-TIME SIGNAL PROCESSING
THIRD EDITION DISCRETE-TIME SIGNAL PROCESSING ALAN V. OPPENHEIM MASSACHUSETTS INSTITUTE OF TECHNOLOGY RONALD W. SCHÄFER HEWLETT-PACKARD LABORATORIES Upper Saddle River Boston Columbus San Francisco New
More informationInvertibility of a room impulse response
Invertibility of a room impulse response Stephen T. Neely a) Jont B. Allen Acoustics Research Department, Bell Laboratories, Murray Hill, New Jersey 07974 (Received 6 July 1979) t When a conversation takes
More informationABSTRACT 1. INTRODUCTION
THE EXCESS PHASE IN LOUDSPEAKER/ROOM TRANSFER FUNCTIONS: Can it be ignored in Equalization Tasks? Lars G. Johansen & Per Rubak DSP-group, Dept. of Communication Technology Institute of Electronic Systems,
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 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 informationSound Listener s perception
Inversion of Loudspeaker Dynamics by Polynomial LQ Feedforward Control Mikael Sternad, Mathias Johansson and Jonas Rutstrom Abstract- Loudspeakers always introduce linear and nonlinear distortions in a
More informationDHANALAKSHMI COLLEGE OF ENGINEERING DEPARTMENT OF ELECTRICAL AND ELECTRONICS ENGINEERING EC2314- DIGITAL SIGNAL PROCESSING UNIT I INTRODUCTION PART A
DHANALAKSHMI COLLEGE OF ENGINEERING DEPARTMENT OF ELECTRICAL AND ELECTRONICS ENGINEERING EC2314- DIGITAL SIGNAL PROCESSING UNIT I INTRODUCTION PART A Classification of systems : Continuous and Discrete
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 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 informationZ - 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 informationE : Lecture 1 Introduction
E85.2607: Lecture 1 Introduction 1 Administrivia 2 DSP review 3 Fun with Matlab E85.2607: Lecture 1 Introduction 2010-01-21 1 / 24 Course overview Advanced Digital Signal Theory Design, analysis, and implementation
More informationAN INVERTIBLE DISCRETE AUDITORY TRANSFORM
COMM. MATH. SCI. Vol. 3, No. 1, pp. 47 56 c 25 International Press AN INVERTIBLE DISCRETE AUDITORY TRANSFORM JACK XIN AND YINGYONG QI Abstract. A discrete auditory transform (DAT) from sound signal to
More informationDFT & Fast Fourier Transform PART-A. 7. Calculate the number of multiplications needed in the calculation of DFT and FFT with 64 point sequence.
SHRI ANGALAMMAN COLLEGE OF ENGINEERING & TECHNOLOGY (An ISO 9001:2008 Certified Institution) SIRUGANOOR,TRICHY-621105. DEPARTMENT OF ELECTRONICS AND COMMUNICATION ENGINEERING UNIT I DFT & Fast Fourier
More informationDIGITAL SIGNAL PROCESSING LECTURE 1
DIGITAL SIGNAL PROCESSING LECTURE 1 Fall 2010 2K8-5 th Semester Tahir Muhammad tmuhammad_07@yahoo.com Content and Figures are from Discrete-Time Signal Processing, 2e by Oppenheim, Shafer, and Buck, 1999-2000
More informationAcoustic Research Institute ARI
Austrian Academy of Sciences Acoustic Research Institute ARI System Identification in Audio Engineering P. Majdak piotr@majdak.com Institut für Schallforschung, Österreichische Akademie der Wissenschaften;
More information/ (2π) X(e jω ) dω. 4. An 8 point sequence is given by x(n) = {2,2,2,2,1,1,1,1}. Compute 8 point DFT of x(n) by
Code No: RR320402 Set No. 1 III B.Tech II Semester Regular Examinations, Apr/May 2006 DIGITAL SIGNAL PROCESSING ( Common to Electronics & Communication Engineering, Electronics & Instrumentation Engineering,
More informationVibration Testing. an excitation source a device to measure the response a digital signal processor to analyze the system response
Vibration Testing For vibration testing, you need an excitation source a device to measure the response a digital signal processor to analyze the system response i) Excitation sources Typically either
More informationLast 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 informationSpeech Signal Representations
Speech Signal Representations Berlin Chen 2003 References: 1. X. Huang et. al., Spoken Language Processing, Chapters 5, 6 2. J. R. Deller et. al., Discrete-Time Processing of Speech Signals, Chapters 4-6
More informationRadar Systems Engineering Lecture 3 Review of Signals, Systems and Digital Signal Processing
Radar Systems Engineering Lecture Review of Signals, Systems and Digital Signal Processing Dr. Robert M. O Donnell Guest Lecturer Radar Systems Course Review Signals, Systems & DSP // Block Diagram of
More informationSignal 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 informationDSP Design Lecture 2. Fredrik Edman.
DSP Design Lecture Number representation, scaling, quantization and round-off Noise Fredrik Edman fredrik.edman@eit.lth.se Representation of Numbers Numbers is a way to use symbols to describe and model
More informationEIGENFILTERS FOR SIGNAL CANCELLATION. Sunil Bharitkar and Chris Kyriakakis
EIGENFILTERS FOR SIGNAL CANCELLATION Sunil Bharitkar and Chris Kyriakakis Immersive Audio Laboratory University of Southern California Los Angeles. CA 9. USA Phone:+1-13-7- Fax:+1-13-7-51, Email:ckyriak@imsc.edu.edu,bharitka@sipi.usc.edu
More informationFrequency Domain Speech Analysis
Frequency Domain Speech Analysis Short Time Fourier Analysis Cepstral Analysis Windowed (short time) Fourier Transform Spectrogram of speech signals Filter bank implementation* (Real) cepstrum and complex
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 informationGeorge Mason University ECE 201: Introduction to Signal Analysis Spring 2017
Assigned: March 20, 2017 Due Date: Week of April 03, 2017 George Mason University ECE 201: Introduction to Signal Analysis Spring 2017 Laboratory Project #6 Due Date Your lab report must be submitted on
More informationDiscrete-time signals and systems
Discrete-time signals and systems 1 DISCRETE-TIME DYNAMICAL SYSTEMS x(t) G y(t) Linear system: Output y(n) is a linear function of the inputs sequence: y(n) = k= h(k)x(n k) h(k): impulse response of the
More informationChapter 9. Linear Predictive Analysis of Speech Signals 语音信号的线性预测分析
Chapter 9 Linear Predictive Analysis of Speech Signals 语音信号的线性预测分析 1 LPC Methods LPC methods are the most widely used in speech coding, speech synthesis, speech recognition, speaker recognition and verification
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 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 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 Architectural Acoustics Session 1pAAa: Advanced Analysis of Room Acoustics:
More informationFeature 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 informationVibration Testing. Typically either instrumented hammers or shakers are used.
Vibration Testing Vibration Testing Equipment For vibration testing, you need an excitation source a device to measure the response a digital signal processor to analyze the system response Excitation
More informationResponses of Digital Filters Chapter Intended Learning Outcomes:
Responses of Digital Filters Chapter Intended Learning Outcomes: (i) Understanding the relationships between impulse response, frequency response, difference equation and transfer function in characterizing
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 informationVALLIAMMAI ENGINEERING COLLEGE. SRM Nagar, Kattankulathur DEPARTMENT OF INFORMATION TECHNOLOGY. Academic Year
VALLIAMMAI ENGINEERING COLLEGE SRM Nagar, Kattankulathur- 603 203 DEPARTMENT OF INFORMATION TECHNOLOGY Academic Year 2016-2017 QUESTION BANK-ODD SEMESTER NAME OF THE SUBJECT SUBJECT CODE SEMESTER YEAR
More informationFESI DOCUMENT A5 Acoustics in rooms
FESI DOCUMENT A5 Acoustics in rooms FileName: A5 Acoustics in rooms - English version.pdf o o Abstract: The first theory developed was Sabine\'s formula (1902) and it is the basis of the so-called \"classic
More informationConverting Infinite Impulse Response Filters to Parallel Form
PUBLISHED IN THE IEEE SIGNAL PROCESSING MAGAZINE, VOL. 35, NO. 3, PP. 124 13, MAY 218 1 Converting Infinite Impulse Response Filters to Parallel Form Balázs Bank, Member, IEEE Abstract Discrete-time rational
More informationODEON APPLICATION NOTE Calibration of Impulse Response Measurements
ODEON APPLICATION NOTE Calibration of Impulse Response Measurements Part 2 Free Field Method GK, CLC - May 2015 Scope In this application note we explain how to use the Free-field calibration tool in ODEON
More informationMITIGATING UNCORRELATED PERIODIC DISTURBANCE IN NARROWBAND ACTIVE NOISE CONTROL SYSTEMS
17th European Signal Processing Conference (EUSIPCO 29) Glasgow, Scotland, August 24-28, 29 MITIGATING UNCORRELATED PERIODIC DISTURBANCE IN NARROWBAND ACTIVE NOISE CONTROL SYSTEMS Muhammad Tahir AKHTAR
More informationA new Audacity feature: room objective acustical parameters calculation module
A new Audacity feature: room objective acustical parameters calculation module Simone CAMPANINI and Angelo FARINA Departement of Industrial Engineering - University of Parma viale G.P.Usberti, 181/A 431
More informationMethods for Synthesizing Very High Q Parametrically Well Behaved Two Pole Filters
Methods for Synthesizing Very High Q Parametrically Well Behaved Two Pole Filters Max Mathews Julius O. Smith III Center for Computer Research in Music and Acoustics (CCRMA) Department of Music, Stanford
More informationDSP Configurations. responded with: thus the system function for this filter would be
DSP Configurations In this lecture we discuss the different physical (or software) configurations that can be used to actually realize or implement DSP functions. Recall that the general form of a DSP
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 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 informationHow to manipulate Frequencies in Discrete-time Domain? Two Main Approaches
How to manipulate Frequencies in Discrete-time Domain? Two Main Approaches Difference Equations (an LTI system) x[n]: input, y[n]: output That is, building a system that maes use of the current and previous
More informationA SPARSENESS CONTROLLED PROPORTIONATE ALGORITHM FOR ACOUSTIC ECHO CANCELLATION
6th European Signal Processing Conference (EUSIPCO 28), Lausanne, Switzerland, August 25-29, 28, copyright by EURASIP A SPARSENESS CONTROLLED PROPORTIONATE ALGORITHM FOR ACOUSTIC ECHO CANCELLATION Pradeep
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 informationAn Fir-Filter Example: Hanning Filter
An Fir-Filter Example: Hanning Filter Josef Goette Bern University of Applied Sciences, Biel Institute of Human Centered Engineering - microlab Josef.Goette@bfh.ch February 7, 2018 Contents 1 Mathematical
More informationDigital Signal Processing
COMP ENG 4TL4: Digital Signal Processing Notes for Lecture #21 Friday, October 24, 2003 Types of causal FIR (generalized) linear-phase filters: Type I: Symmetric impulse response: with order M an even
More informationDepartment of Electrical and Computer Engineering Digital Speech Processing Homework No. 6 Solutions
Problem 1 Department of Electrical and Computer Engineering Digital Speech Processing Homework No. 6 Solutions The complex cepstrum, ˆx[n], of a sequence x[n] is the inverse Fourier transform of the complex
More informationLet H(z) = P(z)/Q(z) be the system function of a rational form. Let us represent both P(z) and Q(z) as polynomials of z (not z -1 )
Review: Poles and Zeros of Fractional Form Let H() = P()/Q() be the system function of a rational form. Let us represent both P() and Q() as polynomials of (not - ) Then Poles: the roots of Q()=0 Zeros:
More informationEE482: Digital Signal Processing Applications
Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu EE482: Digital Signal Processing Applications Spring 2014 TTh 14:30-15:45 CBC C222 Lecture 02 DSP Fundamentals 14/01/21 http://www.ee.unlv.edu/~b1morris/ee482/
More informationSource localization and separation for binaural hearing aids
Source localization and separation for binaural hearing aids Mehdi Zohourian, Gerald Enzner, Rainer Martin Listen Workshop, July 218 Institute of Communication Acoustics Outline 1 Introduction 2 Binaural
More informationThe practical application of G and C 50 in classrooms
The practical application of G and C 50 in classrooms Jack Harvie-Clark 1, Nicholas Dobinson 2 1 Apex Acoustics Ltd 2 Apex Acoustics Ltd ABSTRACT Reverberation time remains the primary indicator of room
More informationVel Tech High Tech Dr.Ranagarajan Dr.Sakunthala Engineering College Department of ECE
Subject Code: EC6502 Course Code:C302 Course Name: PRINCIPLES OF DIGITAL SIGNAL PROCESSING L-3 : T-1 : P-0 : Credits 4 COURSE OBJECTIVES: 1. To learn discrete Fourier transform and its properties 2. To
More informationModeling Measurement Uncertainty in Room Acoustics P. Dietrich
Modeling Measurement Uncertainty in Room Acoustics P. Dietrich This paper investigates a way of determining and modeling uncertainty contributions in measurements of room acoustic parameters, which are
More informationACOUSTIC CLARITY AND AUDITORY ROOM SIZE PERCEPTION. Densil Cabrera 1. Sydney, NSW 2006, Australia
ICSV14 Cairns Australia 9-12 July, 2007 ACOUSTIC CLARITY AND AUDITORY ROOM SIZE PERCEPTION Densil Cabrera 1 1 Faculty of Architecture, Design and Planning, University of Sydney Sydney, NSW 2006, Australia
More informationCorrelator I. Basics. Chapter Introduction. 8.2 Digitization Sampling. D. Anish Roshi
Chapter 8 Correlator I. Basics D. Anish Roshi 8.1 Introduction A radio interferometer measures the mutual coherence function of the electric field due to a given source brightness distribution in the sky.
More informationEE 225D LECTURE ON DIGITAL FILTERS. University of California Berkeley
University of California Berkeley College of Engineering Department of Electrical Engineering and Computer Sciences Professors : N.Morgan / B.Gold EE225D Digital Filters Spring,1999 Lecture 7 N.MORGAN
More informationEnhancement of Acoustically Reverberant Speech Using Cepstral Methods
Enhancement of Acoustically Reverberant Speech Using Cepstral Methods Duncan Charles Bees A thesis submitted to the Faculty of Graduate Studies and Research in partial fulfillment of the requirements for
More informationPredicting 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 informationTransform analysis of LTI systems Oppenheim and Schafer, Second edition pp For LTI systems we can write
Transform analysis of LTI systems Oppenheim and Schafer, Second edition pp. 4 9. For LTI systems we can write yœn D xœn hœn D X kd xœkhœn Alternatively, this relationship can be expressed in the z-transform
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 informationLecture 19 IIR Filters
Lecture 19 IIR Filters Fundamentals of Digital Signal Processing Spring, 2012 Wei-Ta Chu 2012/5/10 1 General IIR Difference Equation IIR system: infinite-impulse response system The most general class
More informationStatistical and Adaptive Signal Processing
r Statistical and Adaptive Signal Processing Spectral Estimation, Signal Modeling, Adaptive Filtering and Array Processing Dimitris G. Manolakis Massachusetts Institute of Technology Lincoln Laboratory
More 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 informationChirp Transform for FFT
Chirp Transform for FFT Since the FFT is an implementation of the DFT, it provides a frequency resolution of 2π/N, where N is the length of the input sequence. If this resolution is not sufficient in a
More informationAccuracy of speech transmission index predictions based on the reverberation time and signal-to-noise ratio Galbrun, Laurent; Kitapci, Kivanc
Heriot-Watt University Heriot-Watt University Research Gateway Accuracy of speech transmission index predictions based on the reverberation time and signal-to-noise ratio Galbrun, Laurent; Kitapci, Kivanc
More informationDESIGN OF QUANTIZED FIR FILTER USING COMPENSATING ZEROS
DESIGN OF QUANTIZED FIR FILTER USING COMPENSATING ZEROS Nivedita Yadav, O.P. Singh, Ashish Dixit Department of Electronics and Communication Engineering, Amity University, Lucknow Campus, Lucknow, (India)
More informationFOR many applications in room acoustics, when a microphone
IEEE TRANSACTIONS ON SPEECH AND AUDIO PROCESSING, VOL. 8, NO. 3, MAY 2000 311 Equalization in an Acoustic Reverberant Environment: Robustness Results Biljana D. Radlović, Student Member, IEEE, Robert C.
More informationOn the Frequency-Domain Properties of Savitzky-Golay Filters
On the Frequency-Domain Properties of Savitzky-Golay Filters Ronald W Schafer HP Laboratories HPL-2-9 Keyword(s): Savitzky-Golay filter, least-squares polynomial approximation, smoothing Abstract: This
More informationPS403 - Digital Signal processing
PS403 - Digital Signal processing 6. DSP - Recursive (IIR) Digital Filters Key Text: Digital Signal Processing with Computer Applications (2 nd Ed.) Paul A Lynn and Wolfgang Fuerst, (Publisher: John Wiley
More informationReverberation Impulse Response Analysis
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
More informationINFLUENCE OF TIME-VARIANCE IN AUDITORIUM ON IMPULSE RESPONSE MEASUREMENT
INFLUENCE OF TIME-VARIANCE IN AUDITORIUM ON IMPULSE RESPONSE MEASUREMENT PACS: 43.55.MC Satoh Fumiaki* 1 ; Nagayama Mitsuru* 1 ; Tachibana Hideki* 2 *1 Chiba Institution of Technology Tsudanuma 2-17-1,
More informationAN SVD-BASED MIMO EQUALIZER APPLIED TO THE AURALIZATION OF AIRCRAFT NOISE IN A CABIN SI- MULATOR
AN SVD-BASED MIMO EQUALIZER APPLIED TO THE AURALIZATION OF AIRCRAFT NOISE IN A CABIN SI- MULATOR Luiz F. O. Chamon, Giuliano S. Quiqueto, Sylvio R. Bistafa Noise and Vibration Group Mechanical Engineering
More informationLecture - 1 Introduction to Adaptive Filters
Adaptive Signal Processing Prof. M. Chakraborty Department of Electronics & Electrical Communication Engineering Indian Institute of Technology, Kharagpur Lecture - 1 Introduction to Adaptive Filters So,
More informationDiscrete-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 informationA.1 THE SAMPLED TIME DOMAIN AND THE Z TRANSFORM. 0 δ(t)dt = 1, (A.1) δ(t)dt =
APPENDIX A THE Z TRANSFORM One of the most useful techniques in engineering or scientific analysis is transforming a problem from the time domain to the frequency domain ( 3). Using a Fourier or Laplace
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 informationDEPARTMENT OF EI DIGITAL SIGNAL PROCESSING ASSIGNMENT 1
This PDF is Created by Simpo PDF Creator unregistered version - http://wwwsimpopdfcom Study hard, for the well is deep, and our brains are shallow DEPARTMENT OF EI DIGITAL SIGNAL PROCESSING ASSIGNMENT
More informationRoll No. :... Invigilator s Signature :.. CS/B.Tech (EE-N)/SEM-6/EC-611/ DIGITAL SIGNAL PROCESSING. Time Allotted : 3 Hours Full Marks : 70
Name : Roll No. :.... Invigilator s Signature :.. CS/B.Tech (EE-N)/SEM-6/EC-611/2011 2011 DIGITAL SIGNAL PROCESSING Time Allotted : 3 Hours Full Marks : 70 The figures in the margin indicate full marks.
More informationR13 SET - 1
R13 SET - 1 III B. Tech II Semester Regular Examinations, April - 2016 DIGITAL SIGNAL PROCESSING (Electronics and Communication Engineering) Time: 3 hours Maximum Marks: 70 Note: 1. Question Paper consists
More informationOn the Use of A Priori Knowledge in Adaptive Inverse Control
54 IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS PART I: FUNDAMENTAL THEORY AND APPLICATIONS, VOL 47, NO 1, JANUARY 2000 On the Use of A Priori Knowledge in Adaptive Inverse Control August Kaelin, Member,
More informationNon-linear decays in simple spaces and their possible exploitation
akutek www.akutek.info PRESENTS Non-linear decays in simple spaces and their possible exploitation by Mike Barron ABSTRACT Standard room acoustic theory is based on the assumption that the sound field
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 informationL8: Source estimation
L8: Source estimation Glottal and lip radiation models Closed-phase residual analysis Voicing/unvoicing detection Pitch detection Epoch detection This lecture is based on [Taylor, 2009, ch. 11-12] Introduction
More informationChapter 3 Data Acquisition and Manipulation
1 Chapter 3 Data Acquisition and Manipulation In this chapter we introduce z transf orm, or the discrete Laplace Transform, to solve linear recursions. Section 3.1 z-transform Given a data stream x {x
More informationDigital Filter Implementation 1
Digital Filter Implementation 1 V 5, November 13, 2017 Christian Feldbauer, Bernhard Geiger, Josef Kulmer, {kulmer}@tugraz.at Signal Processing and Speech Communication Laboratory, www.spsc.tugraz.at Inffeldgasse
More informationMeasuring 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 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 informationInternational Journal of Advanced Research in Computer Science and Software Engineering
Volume 2, Issue 11, November 2012 ISSN: 2277 128X International Journal of Advanced Research in Computer Science and Software Engineering Research Paper Available online at: www.ijarcsse.com Acoustic Source
More informationLAB 6: FIR Filter Design Summer 2011
University of Illinois at Urbana-Champaign Department of Electrical and Computer Engineering ECE 311: Digital Signal Processing Lab Chandra Radhakrishnan Peter Kairouz LAB 6: FIR Filter Design Summer 011
More informationFast Deconvolution of Multi-Channel Systems using Regularisation
Fast Deconvolution of Multi-Channel Systems using Regularisation by Ole Kirkeby, Philip A. Nelson, Hareo Hamada*, and Felipe Orduna-Bustamante # Institute of Sound & Vibration Research, Southampton University,
More informationFINITE PRECISION EFFECTS 1. FLOATING POINT VERSUS FIXED POINT 3. TYPES OF FINITE PRECISION EFFECTS 4. FACTORS INFLUENCING FINITE PRECISION EFFECTS
FINITE PRECISION EFFECTS 1. FLOATING POINT VERSUS FIXED POINT 2. WHEN IS FIXED POINT NEEDED? 3. TYPES OF FINITE PRECISION EFFECTS 4. FACTORS INFLUENCING FINITE PRECISION EFFECTS 5. FINITE PRECISION EFFECTS:
More informationA REVERBERATOR BASED ON ABSORBENT ALL-PASS FILTERS. Luke Dahl, Jean-Marc Jot
Proceedings of the COST G-6 Conference on Digital Audio Effects (DAFX-00), Verona, Italy, December 7-9, 000 A REVERBERATOR BASED ON ABSORBENT ALL-PASS FILTERS Lue Dahl, Jean-Marc Jot Creative Advanced
More informationAnalog vs. discrete signals
Analog vs. discrete signals Continuous-time signals are also known as analog signals because their amplitude is analogous (i.e., proportional) to the physical quantity they represent. Discrete-time signals
More informationDIGITAL SIGNAL PROCESSING UNIT III INFINITE IMPULSE RESPONSE DIGITAL FILTERS. 3.6 Design of Digital Filter using Digital to Digital
DIGITAL SIGNAL PROCESSING UNIT III INFINITE IMPULSE RESPONSE DIGITAL FILTERS Contents: 3.1 Introduction IIR Filters 3.2 Transformation Function Derivation 3.3 Review of Analog IIR Filters 3.3.1 Butterworth
More informationOn Information Maximization and Blind Signal Deconvolution
On Information Maximization and Blind Signal Deconvolution A Röbel Technical University of Berlin, Institute of Communication Sciences email: roebel@kgwtu-berlinde Abstract: In the following paper we investigate
More information