Extraction of Musical Pitches from Recorded Music. Mark Palenik
|
|
- Jemima Hicks
- 5 years ago
- Views:
Transcription
1 Extraction of Musical Pitches from Recorded Music Mark Palenik ABSTRACT Methods of determining the musical pitches heard by the human ear hears when recorded music is played were investigated. The ultimate goal was to find a good method for accurately converting recordings into midi, thereby gaining a complete set of notation information. For this project, only methods for determining pitch information, as opposed to rhythm or key signature information, were investigated. Various methods were investigated with the hopes of getting the best response from complex polyphony. I. Introduction When objects vibrate, they produce a mechanical disturbance in the surrounding medium called sound. This disturbance propagates through the medium as a wave whose properties are determined by the motion of the object creating the disturbance and properties of the medium itself. The human ear is the product of millions of years of evolution in an environment filled with a medium, air, that contains many sound waves produced by various sources. As such, it became very advantageous to be able to discriminate between different sound sources and interpret their meanings. Although it isn't all that difficult to a computer with a proper microphone to produce an accurate recording of a sound wave, software, at present, isn't even close to being capable of replicating the type of interpretation of sound that the human ear and brain are capable of together. There has been a lot of research over the years into voice and recognition, and despite the fact that the technology has its limitations, there have been some impressive results. There has been considerably less research, however, into musical recognition, most likely because the technology is considered less important and useful. First, various methods of decomposing waves into frequency components and extracting phase information were explored. This stage of the program was meant to emulate the function of the cochlea in the ear, which performs a similar decomposition of sound waves. From there, the frequency data was sent through a peak search, to find the location of the dominant frequencies. These frequencies were then fed through an algorithm which fit them to a tunable musical scale. Algorithms for analyzing phase information and spectral comparisons were written to help determine audible pitches from overtones and noise, but have not as of yet been implemented. II. Background Looking at a simplified model of human hearing, it could be said that sound recognition occurs in two stages. First, there's the input of sound into the ear, which decomposes sound waves into their constituent harmonics. Second of all, there is the interpretation of this sound by the brain, which filters out extraneous sounds, and matches the relevant ones to some kind of meaning. This project emulates both stages, relying on computational, rather than
2 mechanical methods for many of the aspects in the first step. The cochlea in the ear is a curled piece of material whose thickness varies along its length, giving it a different resonant frequency at each point. Thus, the location of resonance on the cochlea can be used by the brain to determine the frequency of incoming sound waves. Attempts have been made to emulate the function of the cochlea by using a series of overlapping bandpass filters, while other, more common methods of frequency decomposition involve using the Fourier transform, or one of its variants. The first method I explored for frequency decomposition is similar to the bandpass filter method. However, this method was not based on previous research, but rather, something derived independently, from researching cochlear function. It makes use of a Green's function to simulate harmonic oscillator response to a sound wave. Part III: Program Methods The program written for this project decomposed waves into frequency components, did some further signal analysis, found the peak frequencies, and matched them to corresponding notes. The method for doing each of these things is described below. Frequency Decomposition - Green's Function Green's functions are used for solving differential equations that include a term for a time dependent force. The exact solutions for an arbitrary force function are found by taking the limiting case of finite step impulse functions as they approach infinitesimal delta functions. However, if instead of taking the limit as the width of the step size goes to zero, the solution is left as a discrete sum over small, but finite steps, the resulting function is perfect for use with a recorded digital signal. The Green's function used finds the response of a damped harmonic oscillator with a natural frequency w 0. The resulting function is of the form: T B [ t tn] a[] t e sin w[ t tn] w where T is the duration of each impulse, a[t] is the acceleration as a function of time, t is the current time, w is the frequency of the oscillator after dampening, and B is a dampening constant. The function is summed over tn. A set of linearly spaced oscillators were created and their output amplitude measured after responding to a portion of the recorded sound. The amplitudes were then plotted to create a spectrum. Amplitudes were calculated by the phasor addition of each term in the sum. This required doing a second sum over cos(w[t-tn]) and calculating the vector length of the two sums. Only around 250 oscillators are necessary to get well resolved peaks, which is good since large numbers of oscillators cause significant slowdown calculation time. Since the acceleration, rather than force, appears in the function, an appropriate method of determining the mass of the oscillator must be determined, since the amplitude of the recorded sound is proportional to its force, and thus will have a different effect on objects of different masses.
3 The boundary conditions for an actual cochlea are somewhat difficult to emulate. The resonant frequency changes smoothly as a function of position, making it difficult to really treat the object as a collection of discrete oscillators. Most likely, some form of Cauchy boundary conditions on a collection of connected strings would be the best way to model its behavior. For the purposes of this project, however, the boundary conditions were simplified to assume that each oscillator was a string of the same tension, obeying Dirichlet boundary conditions with the ends fixed at one wavelength. These conditions were plugged into the wave equation to produce effective masses for each oscillator, which ended up being proportional to 1/w 0. Testing the equation with a constant mass for each oscillator produced a very large response bias at low frequencies. Inserting a mass term proportional to 1/w 0 corrected the bias and produced similar results the Fourier transform. The function was also tested with and without a dampening term. Inserting a dampening term, B, effectively creates a windowing function on the input signal, causing nearby frequencies to blend together, while removing most of the low frequency noise. The best results occurred with a low dampening constant, around 10 or 15, which removed a small amount of the more difficult to compensate for noise, while retaining good peak separation for nearby frequencies. The latter part is especially important, since for pitches near the bottom of the audible range, the frequencies of neighboring notes can be very close together. Fast Fourier Transform (FFT) Another method explored was the fast Fourier transform, or FFT. The main advantage of the FFT is that it performs very quickly. A regular Fourier transform is calculated by: 2pi ikx F k = f x e [ ] [ ] dx A fast Fourier transform speeds up this calculation for discrete functions so that it can be done in time proportional to 2Nln(N), rather than 2N 2, where N is the number of data points to be analyzed. It also makes use of properties of signals, like the Nyquist frequency, which determines the highest frequency that needs to be used in the calculation. The FFT, however, produces too many extraneous peaks, and does not give good enough peak resolution to be used for this project. Since a musical piece will in theory, change with time, the recorded sample must be broken up into small pieces when analyzed. The FFT does not perform well enough on sample sizes on the order of 1024 samples to be used. Even on larger sample sizes, the Green's function method produced better results, while taking considerably longer to calculate. To improve the frequency response of a FFT, a windowing function can be used. I experimented with the Hanning window, which is of the form: 1 pi x 1 + cos 2 a In many cases, this does improve results somewhat, but it has the same
4 negative effects as the dampening term in the Green's function. It causes nearby peaks to be blended together, which can mask notes that should be distinguishable. The Continuous Wavelet Transform (CWT) To help resolve some of the problems with the FFT, a continuous wavelet transform was used. Although created mainly to retrieve good time dependent information from a sample, the continuous wavelet transform also provided better frequency information at a fixed time than the FFT. This was due to the fact that the frequencies and window sizes tested were not chosen to have orthogonality, but rather, for the best possible frequency response. While this removed any possibility of signal reconstruction, that was not important for the analysis being done. The continuous wavelet transform decomposes a sample into wavelets, rather than sine and cosine functions. The Fourier transform can be thought of as the projection of a function onto a set of orthogonal sine and cosine functions. The continuous wavelet transform uses the same projection operation, but uses a series of wavelets as its basis. A wavelet is a function with some sort of frequency, like a sine or cosine function, but the bulk of the function is confined to a narrow window. The frequency of the wavelet is proportional to its size, so larger wavelets are used to measure lower frequency components of a signal, and smaller wavelets measure higher frequency components. The wavelet can be slid forward by a small increment to change the point in time where the signal is being analyzed. The uncertainty in the time at which the various frequency components appear is equal to the wavelet size of each frequency. Because of this, when processing a digital signal, it is only necessary to slide the wavelet forward in increments proportional to its size. For this project, I used wavelets made up of simple sine and cosine functions cut off by a rectangular window. However, I found that making the wavelength of the sine and cosine functions equal to the size of the wavelet did not produce the best frequency response. Rather, I found that, for a wavelet size L, a wavelength varying linearly from L to 2L as L increased from its minimum to its maximum size, provided very good peak resolution, while still allowing relatively small wavelet sizes for low frequencies. I experimented a with a Hanning window again, using a CWT. The Hanning window created the same problems with frequency resolution as before, except magnified to a greater extent. This is most likely because the wavelets I used did not go to zero at both ends of the window like they do in a Fourier series. Multiplying by a Hanning window, which forces both ends of the window to zero would have introduced additional distortion. Low frequency response with the Hanning window did not even resemble the actual frequency components of the signal. All in all, the CWT gave the best frequency resolution, and better time resolution than any other method. I ended up using a maximum window size of around 2000 sample, but because of the nature of the transform, most frequencies are analyzed at a smaller window size. The CWT gave marginally better frequency response than the Green's function method on windows of a similar size. The CWT, however, performed significantly slower.
5 Extracting Phase Information Green's Function Since the code for calculating harmonic oscillator response was already written, it was a simple matter to modify it for use in determining phase information. As cited in Joe Yasi's paper on extracting phase harmonics from a digital signal, a Fourier series does not provide accurate phase information over a small sample size. I confirmed this fact experimentally, and decided to go with a method similar to his bandpass filter method for determining phase. Instead of a bandpass filter, I used a harmonic oscillator with a natural frequency at the frequency being analyzed. Like before, a summation was done over sine and cosine functions. An additional step was then used, which took the arc tangent of the sine summation divided by the cosine summation, to give a result for the phase of the oscillator response in terms of radians. The output varied from -pi to pi. With a few additional calculations, the phase between 0 and pi was found. To determine the phase offset of the signal, rather than the phase at a particular time, the following code was used: If XofTOut(n) < 2 * * (Freq * n * Tau - Int(Freq * n * Tau)) Then XofTOut(n) = 2 * XofTOut(n) - 2 * * (Freq * n * Tau - Int(Freq * n * Tau)) Else XofTOut(n) = XofTOut(n) - 2 * * (Freq * n * Tau - Int(Freq * n * Tau)) End If where XofTOut is an array that stores the phase information at a point in time, n. The result is the phase offset of the frequency component from the start of the window. There was one additional hurdle to overcome in this section of the program. At every point in time at which the phase was calculated, the Green's functions would need to be re-summed, which was impossible to do in any realistic time frame. However, looking at the sum once more: T B [ t tn] a[] t e sin w[ t tn] w It can be seen that each term in the sum varies by a constant phase factor, T. tn increases by T with each consecutive term in the sum. The current time, t, also increases by T whenever the phase calculation is moved forward in time. This means that given a phase at a particular point in time, P(t), P(t + T) is just the same sum with a phase offset, plug one additional term. If the phase of every term in the sum at t + T can be increased by T, it would only take one additional calculation to determine the phase at time t + T. This type of phase promotion is impossible to do with only a sine function. However, given the exponential function e ikx = cos(x) + isin(x), increasing the phase is a simple matter of multiplying by e it. This corresponds to a raising operator in phase space, or a rotation operator in coordinate space. The phase, then is just the arc tangent of the imaginary part of the function over the real part, with the same additional calculations as described before. Peak Search The next step after determining frequency and phase information of the signal was to run a peak search over the frequency spectrum. I used the method
6 described on the Thermogalactic website (a company that makes spectral analysis software and spectroscopy equipment) with a few minor modifications. The average difference between data points is computed, and stored as a value called noise. Then, the program searches for any slopes that pop up above the noise value, and calls it the start of a peak. The end of the peak is marked at the point after the slope goes down by an amount greater than the noise value, then levels off to a smaller value. The actual peak is calculated by taking the center of mass of the three highest points between the start and endpoints. From there, baselines, and other pieces of information are calculated. Spectrum Comparisons I also wrote a piece of code that does a spectral peak search to match similar spectra in order of decreasing similarity and return a percent match value. Since the human ear is bad at distinguishing between frequencies close to each other, I thought the brain might use a method similar to this in music recognition to recognize chords when the intervals between the notes were small. The idea was to test the spectrum calculated against a list of possible spectra, computed from the CWT of a combination of pure tones, however this process was never implemented, due to time constraints. Matching Peaks to Notes After the peak search takes place, the peaks need to be matched to notes on the musical scale, and not just frequencies. To do this, the program calculates a musical scale, based on an adjustable input for middle A (440 hz). The notes in the scale calculated are A, A#, B, C, C#, D, D#, E, F, F#, G, and G# for as many octaves above and below middle A as specified by an input parameter. No allowances are made in the program for spelling a note with a flat instead of a sharp, or for using double sharps, double flats, etc. The determination of the spelling of a piece of music is a complex task for a computer that would be an entire project in itself. The peak frequencies are matched to the nearest note, and for especially close matches, a true value is returned to indicate that the peak is more likely a real note, and not noise, although it may still be a harmonic. IV. More about the program Also included in the program is a function for inverting matrices, a function for doing a polynomial fit of n degrees, based on Euler's method (as described in the textbook Linear Algera by Leon), Direct Sound code for loading and playing.wav files, and a network which I had written earlier and thought I might use in this project, but was unable to get to perform in a useful manner. Below are some screenshots from the program:
7 Frequency information from a trumpet playing a 440 hz A phase information of a trumpet playing a 440 hz A at the 1760 Hz peak Frequency information of an inverted C major chord played on the piano. (the note information obtained by the program is accurate in this image) V. Conclusion So far, the program doesn't perform a whole lot of useful functions. It has a decent algorithm for determining frequency composition of a wave in discrete time intervals, and a good base of functions for a more advanced program. It still needs a lot of work to put it into a more useable form, which I was not able to do, due to time constraints. Ideally, I would like the program to be able to convert a full recorded song into a midi file, with minimal errors. Implementation of a sliding wavelet in time, to gain time as well as frequency information, and a method for determining rhythms need to be added. I will most likely work on adding these features on my own at a later date. Acknowledgements
8 I would like to thank Robert Bristow-Johnson for providing a simulation filter for the Fletcher-Munson curve at 0 db, despite the fact that I never used it. I would also like to thank Professor Errede for running this class and answering a few questions I had while doing my project. References Classical Dynamics of Particles and Systems: Fifth Edition, Thorton&Marion. Linear Algebra with Applications, Leon The Engineers Guide to Wavelet Analysis, Robi Polikar Cochlear Fluids Laborotory, Alec N. Salt, PhD An Algorithm for Extracting the Relative Phase Harmonics from a Periodic Digital Signal, Joe Yasi.
Sound is the human ear s perceived effect of pressure changes in the ambient air. Sound can be modeled as a function of time.
2. Physical sound 2.1 What is sound? Sound is the human ear s perceived effect of pressure changes in the ambient air. Sound can be modeled as a function of time. Figure 2.1: A 0.56-second audio clip of
More informationTRANSFORMS / WAVELETS
RANSFORMS / WAVELES ransform Analysis Signal processing using a transform analysis for calculations is a technique used to simplify or accelerate problem solution. For example, instead of dividing two
More informationMusical Acoustics, C. Bertulani. Musical Acoustics. Lecture 14 Timbre / Tone quality II
1 Musical Acoustics Lecture 14 Timbre / Tone quality II Odd vs Even Harmonics and Symmetry Sines are Anti-symmetric about mid-point If you mirror around the middle you get the same shape but upside down
More informationSignals A Preliminary Discussion EE442 Analog & Digital Communication Systems Lecture 2
Signals A Preliminary Discussion EE442 Analog & Digital Communication Systems Lecture 2 The Fourier transform of single pulse is the sinc function. EE 442 Signal Preliminaries 1 Communication Systems and
More informationCOMP 546, Winter 2017 lecture 20 - sound 2
Today we will examine two types of sounds that are of great interest: music and speech. We will see how a frequency domain analysis is fundamental to both. Musical sounds Let s begin by briefly considering
More informationIntroduction to Wavelets Michael Phipps Vallary Bhopatkar
Introduction to Wavelets Michael Phipps Vallary Bhopatkar *Amended from The Wavelet Tutorial by Robi Polikar, http://users.rowan.edu/~polikar/wavelets/wttutoria Who can tell me what this means? NR3, pg
More informationOrthonormal bases and tilings of the time-frequency plane for music processing Juan M. Vuletich *
Orthonormal bases and tilings of the time-frequency plane for music processing Juan M. Vuletich * Dept. of Computer Science, University of Buenos Aires, Argentina ABSTRACT Conventional techniques for signal
More informationPHYSICS LAB. Sound. Date: GRADE: PHYSICS DEPARTMENT JAMES MADISON UNIVERSITY
PHYSICS LAB Sound Printed Names: Signatures: Date: Lab Section: Instructor: GRADE: PHYSICS DEPARTMENT JAMES MADISON UNIVERSITY Revision August 2003 Sound Investigations Sound Investigations 78 Part I -
More informationFFT 1 /n octave analysis wavelet
06/16 For most acoustic examinations, a simple sound level analysis is insufficient, as not only the overall sound pressure level, but also the frequency-dependent distribution of the level has a significant
More informationLecture 3 Complex Exponential Signals
Lecture 3 Complex Exponential Signals Fundamentals of Digital Signal Processing Spring, 2012 Wei-Ta Chu 2012/3/1 1 Review of Complex Numbers Using Euler s famous formula for the complex exponential The
More informationME scope Application Note 01 The FFT, Leakage, and Windowing
INTRODUCTION ME scope Application Note 01 The FFT, Leakage, and Windowing NOTE: The steps in this Application Note can be duplicated using any Package that includes the VES-3600 Advanced Signal Processing
More informationFFT analysis in practice
FFT analysis in practice Perception & Multimedia Computing Lecture 13 Rebecca Fiebrink Lecturer, Department of Computing Goldsmiths, University of London 1 Last Week Review of complex numbers: rectangular
More informationThe Fast Fourier Transform
The Fast Fourier Transform Basic FFT Stuff That s s Good to Know Dave Typinski, Radio Jove Meeting, July 2, 2014, NRAO Green Bank Ever wonder how an SDR-14 or Dongle produces the spectra that it does?
More informationChapter 18. Superposition and Standing Waves
Chapter 18 Superposition and Standing Waves Particles & Waves Spread Out in Space: NONLOCAL Superposition: Waves add in space and show interference. Do not have mass or Momentum Waves transmit energy.
More informationStructure of Speech. Physical acoustics Time-domain representation Frequency domain representation Sound shaping
Structure of Speech Physical acoustics Time-domain representation Frequency domain representation Sound shaping Speech acoustics Source-Filter Theory Speech Source characteristics Speech Filter characteristics
More informationBiomedical Signals. Signals and Images in Medicine Dr Nabeel Anwar
Biomedical Signals Signals and Images in Medicine Dr Nabeel Anwar Noise Removal: Time Domain Techniques 1. Synchronized Averaging (covered in lecture 1) 2. Moving Average Filters (today s topic) 3. Derivative
More informationUNIVERSITY OF TORONTO Faculty of Arts and Science MOCK EXAMINATION PHY207H1S. Duration 3 hours NO AIDS ALLOWED
UNIVERSITY OF TORONTO Faculty of Arts and Science MOCK EXAMINATION PHY207H1S Duration 3 hours NO AIDS ALLOWED Instructions: Please answer all questions in the examination booklet(s) provided. Completely
More informationCSC475 Music Information Retrieval
CSC475 Music Information Retrieval Sinusoids and DSP notation George Tzanetakis University of Victoria 2014 G. Tzanetakis 1 / 38 Table of Contents I 1 Time and Frequency 2 Sinusoids and Phasors G. Tzanetakis
More informationNotes on Fourier transforms
Fourier Transforms 1 Notes on Fourier transforms The Fourier transform is something we all toss around like we understand it, but it is often discussed in an offhand way that leads to confusion for those
More informationFourier and Wavelets
Fourier and Wavelets Why do we need a Transform? Fourier Transform and the short term Fourier (STFT) Heisenberg Uncertainty Principle The continues Wavelet Transform Discrete Wavelet Transform Wavelets
More informationA mechanical wave is a disturbance which propagates through a medium with little or no net displacement of the particles of the medium.
Waves and Sound Mechanical Wave A mechanical wave is a disturbance which propagates through a medium with little or no net displacement of the particles of the medium. Water Waves Wave Pulse People Wave
More informationStandard Octaves and Sound Pressure. The superposition of several independent sound sources produces multifrequency noise: i=1
Appendix C Standard Octaves and Sound Pressure C.1 Time History and Overall Sound Pressure The superposition of several independent sound sources produces multifrequency noise: p(t) = N N p i (t) = P i
More informationACOUSTICS. Sounds are vibrations in the air, extremely small and fast fluctuations of airpressure.
ACOUSTICS 1. VIBRATIONS Sounds are vibrations in the air, extremely small and fast fluctuations of airpressure. These vibrations are generated from sounds sources and travel like waves in the water; sound
More informationComparison of a Pleasant and Unpleasant Sound
Comparison of a Pleasant and Unpleasant Sound B. Nisha 1, Dr. S. Mercy Soruparani 2 1. Department of Mathematics, Stella Maris College, Chennai, India. 2. U.G Head and Associate Professor, Department of
More informationProf. Feng Liu. Fall /04/2018
Prof. Feng Liu Fall 2018 http://www.cs.pdx.edu/~fliu/courses/cs447/ 10/04/2018 1 Last Time Image file formats Color quantization 2 Today Dithering Signal Processing Homework 1 due today in class Homework
More informationAcoustics and Fourier Transform Physics Advanced Physics Lab - Summer 2018 Don Heiman, Northeastern University, 1/12/2018
1 Acoustics and Fourier Transform Physics 3600 - Advanced Physics Lab - Summer 2018 Don Heiman, Northeastern University, 1/12/2018 I. INTRODUCTION Time is fundamental in our everyday life in the 4-dimensional
More informationLaboratory Assignment 4. Fourier Sound Synthesis
Laboratory Assignment 4 Fourier Sound Synthesis PURPOSE This lab investigates how to use a computer to evaluate the Fourier series for periodic signals and to synthesize audio signals from Fourier series
More informationChapter 2. Meeting 2, Measures and Visualizations of Sounds and Signals
Chapter 2. Meeting 2, Measures and Visualizations of Sounds and Signals 2.1. Announcements Be sure to completely read the syllabus Recording opportunities for small ensembles Due Wednesday, 15 February:
More informationChapter 12. Preview. Objectives The Production of Sound Waves Frequency of Sound Waves The Doppler Effect. Section 1 Sound Waves
Section 1 Sound Waves Preview Objectives The Production of Sound Waves Frequency of Sound Waves The Doppler Effect Section 1 Sound Waves Objectives Explain how sound waves are produced. Relate frequency
More informationVariations in Waveforms and Energy Spectra between Musical Instruments
Mahalia Lotz Dr. Grant Gustafson MATH 2270 5/6/2016 Variations in Waveforms and Energy Spectra between Musical Instruments Sound occurs when particles are displaced by some initial motion to create a wave-like
More informationResponse spectrum Time history Power Spectral Density, PSD
A description is given of one way to implement an earthquake test where the test severities are specified by time histories. The test is done by using a biaxial computer aided servohydraulic test rig.
More informationFundamentals of Digital Audio *
Digital Media The material in this handout is excerpted from Digital Media Curriculum Primer a work written by Dr. Yue-Ling Wong (ylwong@wfu.edu), Department of Computer Science and Department of Art,
More informationWhat is Sound? Part II
What is Sound? Part II Timbre & Noise 1 Prayouandi (2010) - OneOhtrix Point Never PSYCHOACOUSTICS ACOUSTICS LOUDNESS AMPLITUDE PITCH FREQUENCY QUALITY TIMBRE 2 Timbre / Quality everything that is not frequency
More informationSignals and Systems Lecture 9 Communication Systems Frequency-Division Multiplexing and Frequency Modulation (FM)
Signals and Systems Lecture 9 Communication Systems Frequency-Division Multiplexing and Frequency Modulation (FM) April 11, 2008 Today s Topics 1. Frequency-division multiplexing 2. Frequency modulation
More informationPerception of pitch. Definitions. Why is pitch important? BSc Audiology/MSc SHS Psychoacoustics wk 4: 7 Feb A. Faulkner.
Perception of pitch BSc Audiology/MSc SHS Psychoacoustics wk 4: 7 Feb 2008. A. Faulkner. See Moore, BCJ Introduction to the Psychology of Hearing, Chapter 5. Or Plack CJ The Sense of Hearing Lawrence Erlbaum,
More informationThe Discrete Fourier Transform. Claudia Feregrino-Uribe, Alicia Morales-Reyes Original material: Dr. René Cumplido
The Discrete Fourier Transform Claudia Feregrino-Uribe, Alicia Morales-Reyes Original material: Dr. René Cumplido CCC-INAOE Autumn 2015 The Discrete Fourier Transform Fourier analysis is a family of mathematical
More informationTrigonometric functions and sound
Trigonometric functions and sound The sounds we hear are caused by vibrations that send pressure waves through the air. Our ears respond to these pressure waves and signal the brain about their amplitude
More information6 Sampling. Sampling. The principles of sampling, especially the benefits of coherent sampling
Note: Printed Manuals 6 are not in Color Objectives This chapter explains the following: The principles of sampling, especially the benefits of coherent sampling How to apply sampling principles in a test
More informationIntroduction to Wavelets. For sensor data processing
Introduction to Wavelets For sensor data processing List of topics Why transform? Why wavelets? Wavelets like basis components. Wavelets examples. Fast wavelet transform. Wavelets like filter. Wavelets
More informationApplication of Fourier Transform in Signal Processing
1 Application of Fourier Transform in Signal Processing Lina Sun,Derong You,Daoyun Qi Information Engineering College, Yantai University of Technology, Shandong, China Abstract: Fourier transform is a
More informationProperties and Applications
Properties and Applications What is a Wave? How is it Created? Waves are created by vibrations! Atoms vibrate, strings vibrate, water vibrates A wave is the moving oscillation Waves are the propagation
More informationMichael F. Toner, et. al.. "Distortion Measurement." Copyright 2000 CRC Press LLC. <
Michael F. Toner, et. al.. "Distortion Measurement." Copyright CRC Press LLC. . Distortion Measurement Michael F. Toner Nortel Networks Gordon W. Roberts McGill University 53.1
More informationIntroduction to Wavelet Transform. Chapter 7 Instructor: Hossein Pourghassem
Introduction to Wavelet Transform Chapter 7 Instructor: Hossein Pourghassem Introduction Most of the signals in practice, are TIME-DOMAIN signals in their raw format. It means that measured signal is a
More informationIntroduction to Telecommunications and Computer Engineering Unit 3: Communications Systems & Signals
Introduction to Telecommunications and Computer Engineering Unit 3: Communications Systems & Signals Syedur Rahman Lecturer, CSE Department North South University syedur.rahman@wolfson.oxon.org Acknowledgements
More informationPerception of pitch. Importance of pitch: 2. mother hemp horse. scold. Definitions. Why is pitch important? AUDL4007: 11 Feb A. Faulkner.
Perception of pitch AUDL4007: 11 Feb 2010. A. Faulkner. See Moore, BCJ Introduction to the Psychology of Hearing, Chapter 5. Or Plack CJ The Sense of Hearing Lawrence Erlbaum, 2005 Chapter 7 1 Definitions
More informationAP Physics B (Princeton 15 & Giancoli 11 & 12) Waves and Sound
AP Physics B (Princeton 15 & Giancoli 11 & 12) Waves and Sound Preview What are the two categories of waves with regard to mode of travel? Mechanical Electromagnetic Which type of wave requires a medium?
More informationFrequency analysis put into practice
Technically, guitar strings, audio amplifiers, filters or rotating shafts are one and the same, namely signal sources. These contain substantial information. The content is decoded during the oscilloscopic
More informationTopic 2. Signal Processing Review. (Some slides are adapted from Bryan Pardo s course slides on Machine Perception of Music)
Topic 2 Signal Processing Review (Some slides are adapted from Bryan Pardo s course slides on Machine Perception of Music) Recording Sound Mechanical Vibration Pressure Waves Motion->Voltage Transducer
More informationPerception of pitch. Definitions. Why is pitch important? BSc Audiology/MSc SHS Psychoacoustics wk 5: 12 Feb A. Faulkner.
Perception of pitch BSc Audiology/MSc SHS Psychoacoustics wk 5: 12 Feb 2009. A. Faulkner. See Moore, BCJ Introduction to the Psychology of Hearing, Chapter 5. Or Plack CJ The Sense of Hearing Lawrence
More informationMusic 171: Amplitude Modulation
Music 7: Amplitude Modulation Tamara Smyth, trsmyth@ucsd.edu Department of Music, University of California, San Diego (UCSD) February 7, 9 Adding Sinusoids Recall that adding sinusoids of the same frequency
More informationAPPENDIX MATHEMATICS OF DISTORTION PRODUCT OTOACOUSTIC EMISSION GENERATION: A TUTORIAL
In: Otoacoustic Emissions. Basic Science and Clinical Applications, Ed. Charles I. Berlin, Singular Publishing Group, San Diego CA, pp. 149-159. APPENDIX MATHEMATICS OF DISTORTION PRODUCT OTOACOUSTIC EMISSION
More informationSubtractive Synthesis without Filters
Subtractive Synthesis without Filters John Lazzaro and John Wawrzynek Computer Science Division UC Berkeley lazzaro@cs.berkeley.edu, johnw@cs.berkeley.edu 1. Introduction The earliest commercially successful
More informationDrum Transcription Based on Independent Subspace Analysis
Report for EE 391 Special Studies and Reports for Electrical Engineering Drum Transcription Based on Independent Subspace Analysis Yinyi Guo Center for Computer Research in Music and Acoustics, Stanford,
More informationSound, acoustics Slides based on: Rossing, The science of sound, 1990.
Sound, acoustics Slides based on: Rossing, The science of sound, 1990. Acoustics 1 1 Introduction Acoustics 2! The word acoustics refers to the science of sound and is a subcategory of physics! Room acoustics
More informationReading: Johnson Ch , Ch.5.5 (today); Liljencrants & Lindblom; Stevens (Tues) reminder: no class on Thursday.
L105/205 Phonetics Scarborough Handout 7 10/18/05 Reading: Johnson Ch.2.3.3-2.3.6, Ch.5.5 (today); Liljencrants & Lindblom; Stevens (Tues) reminder: no class on Thursday Spectral Analysis 1. There are
More informationCHAPTER 6 INTRODUCTION TO SYSTEM IDENTIFICATION
CHAPTER 6 INTRODUCTION TO SYSTEM IDENTIFICATION Broadly speaking, system identification is the art and science of using measurements obtained from a system to characterize the system. The characterization
More informationLecture 7: Superposition and Fourier Theorem
Lecture 7: Superposition and Fourier Theorem Sound is linear. What that means is, if several things are producing sounds at once, then the pressure of the air, due to the several things, will be and the
More information21/01/2014. Fundamentals of the analysis of neuronal oscillations. Separating sources
21/1/214 Separating sources Fundamentals of the analysis of neuronal oscillations Robert Oostenveld Donders Institute for Brain, Cognition and Behaviour Radboud University Nijmegen, The Netherlands Use
More informationModulation. Digital Data Transmission. COMP476 Networked Computer Systems. Analog and Digital Signals. Analog and Digital Examples.
Digital Data Transmission Modulation Digital data is usually considered a series of binary digits. RS-232-C transmits data as square waves. COMP476 Networked Computer Systems Analog and Digital Signals
More informationDiscrete Fourier Transform
6 The Discrete Fourier Transform Lab Objective: The analysis of periodic functions has many applications in pure and applied mathematics, especially in settings dealing with sound waves. The Fourier transform
More informationCHAPTER 5. Digitized Audio Telemetry Standard. Table of Contents
CHAPTER 5 Digitized Audio Telemetry Standard Table of Contents Chapter 5. Digitized Audio Telemetry Standard... 5-1 5.1 General... 5-1 5.2 Definitions... 5-1 5.3 Signal Source... 5-1 5.4 Encoding/Decoding
More informationMusic. Sound Part II
Music Sound Part II What is the study of sound called? Acoustics What is the difference between music and noise? Music: Sound that follows a regular pattern; a mixture of frequencies which have a clear
More informationPhotone Sound Design Tutorial
Photone Sound Design Tutorial An Introduction At first glance, Photone s control elements appear dauntingly complex but this impression is deceiving: Anyone who has listened to all the instrument s presets
More informationFrequency Division Multiplexing Spring 2011 Lecture #14. Sinusoids and LTI Systems. Periodic Sequences. x[n] = x[n + N]
Frequency Division Multiplexing 6.02 Spring 20 Lecture #4 complex exponentials discrete-time Fourier series spectral coefficients band-limited signals To engineer the sharing of a channel through frequency
More informationSHAKER TABLE SEISMIC TESTING OF EQUIPMENT USING HISTORICAL STRONG MOTION DATA SCALED TO SATISFY A SHOCK RESPONSE SPECTRUM
SHAKER TABLE SEISMIC TESTING OF EQUIPMENT USING HISTORICAL STRONG MOTION DATA SCALED TO SATISFY A SHOCK RESPONSE SPECTRUM By Tom Irvine Email: tomirvine@aol.com May 6, 29. The purpose of this paper is
More informationAcoustics, signals & systems for audiology. Week 4. Signals through Systems
Acoustics, signals & systems for audiology Week 4 Signals through Systems Crucial ideas Any signal can be constructed as a sum of sine waves In a linear time-invariant (LTI) system, the response to a sinusoid
More informationDefinition of Sound. Sound. Vibration. Period - Frequency. Waveform. Parameters. SPA Lundeen
Definition of Sound Sound Psychologist's = that which is heard Physicist's = a propagated disturbance in the density of an elastic medium Vibrator serves as the sound source Medium = air 2 Vibration Periodic
More informationSpectrum Analysis: The FFT Display
Spectrum Analysis: The FFT Display Equipment: Capstone, voltage sensor 1 Introduction It is often useful to represent a function by a series expansion, such as a Taylor series. There are other series representations
More informationPattern Recognition. Part 6: Bandwidth Extension. Gerhard Schmidt
Pattern Recognition Part 6: Gerhard Schmidt Christian-Albrechts-Universität zu Kiel Faculty of Engineering Institute of Electrical and Information Engineering Digital Signal Processing and System Theory
More informationVocal Command Recognition Using Parallel Processing of Multiple Confidence-Weighted Algorithms in an FPGA
Vocal Command Recognition Using Parallel Processing of Multiple Confidence-Weighted Algorithms in an FPGA ECE-492/3 Senior Design Project Spring 2015 Electrical and Computer Engineering Department Volgenau
More informationG(f ) = g(t) dt. e i2πft. = cos(2πf t) + i sin(2πf t)
Fourier Transforms Fourier s idea that periodic functions can be represented by an infinite series of sines and cosines with discrete frequencies which are integer multiples of a fundamental frequency
More informationTHE BENEFITS OF DSP LOCK-IN AMPLIFIERS
THE BENEFITS OF DSP LOCK-IN AMPLIFIERS If you never heard of or don t understand the term lock-in amplifier, you re in good company. With the exception of the optics industry where virtually every major
More informationDepartment of Electronic Engineering NED University of Engineering & Technology. LABORATORY WORKBOOK For the Course SIGNALS & SYSTEMS (TC-202)
Department of Electronic Engineering NED University of Engineering & Technology LABORATORY WORKBOOK For the Course SIGNALS & SYSTEMS (TC-202) Instructor Name: Student Name: Roll Number: Semester: Batch:
More informationME 365 EXPERIMENT 8 FREQUENCY ANALYSIS
ME 365 EXPERIMENT 8 FREQUENCY ANALYSIS Objectives: There are two goals in this laboratory exercise. The first is to reinforce the Fourier series analysis you have done in the lecture portion of this course.
More informationSAMPLING THEORY. Representing continuous signals with discrete numbers
SAMPLING THEORY Representing continuous signals with discrete numbers Roger B. Dannenberg Professor of Computer Science, Art, and Music Carnegie Mellon University ICM Week 3 Copyright 2002-2013 by Roger
More informationSignals, Sound, and Sensation
Signals, Sound, and Sensation William M. Hartmann Department of Physics and Astronomy Michigan State University East Lansing, Michigan Л1Р Contents Preface xv Chapter 1: Pure Tones 1 Mathematics of the
More information14 fasttest. Multitone Audio Analyzer. Multitone and Synchronous FFT Concepts
Multitone Audio Analyzer The Multitone Audio Analyzer (FASTTEST.AZ2) is an FFT-based analysis program furnished with System Two for use with both analog and digital audio signals. Multitone and Synchronous
More informationTHE CITADEL THE MILITARY COLLEGE OF SOUTH CAROLINA. Department of Electrical and Computer Engineering. ELEC 423 Digital Signal Processing
THE CITADEL THE MILITARY COLLEGE OF SOUTH CAROLINA Department of Electrical and Computer Engineering ELEC 423 Digital Signal Processing Project 2 Due date: November 12 th, 2013 I) Introduction In ELEC
More information3/23/2015. Chapter 11 Oscillations and Waves. Contents of Chapter 11. Contents of Chapter Simple Harmonic Motion Spring Oscillations
Lecture PowerPoints Chapter 11 Physics: Principles with Applications, 7 th edition Giancoli Chapter 11 and Waves This work is protected by United States copyright laws and is provided solely for the use
More informationObjectives. Abstract. This PRO Lesson will examine the Fast Fourier Transformation (FFT) as follows:
: FFT Fast Fourier Transform This PRO Lesson details hardware and software setup of the BSL PRO software to examine the Fast Fourier Transform. All data collection and analysis is done via the BIOPAC MP35
More informationCreating Digital Music
Chapter 2 Creating Digital Music Chapter 2 exposes students to some of the most important engineering ideas associated with the creation of digital music. Students learn how basic ideas drawn from the
More informationDigital Image Processing
In the Name of Allah Digital Image Processing Introduction to Wavelets Hamid R. Rabiee Fall 2015 Outline 2 Why transform? Why wavelets? Wavelets like basis components. Wavelets examples. Fast wavelet transform.
More informationFourier Signal Analysis
Part 1B Experimental Engineering Integrated Coursework Location: Baker Building South Wing Mechanics Lab Experiment A4 Signal Processing Fourier Signal Analysis Please bring the lab sheet from 1A experiment
More informationBasic Electronics Learning by doing Prof. T.S. Natarajan Department of Physics Indian Institute of Technology, Madras
Basic Electronics Learning by doing Prof. T.S. Natarajan Department of Physics Indian Institute of Technology, Madras Lecture 26 Mathematical operations Hello everybody! In our series of lectures on basic
More informationMath and Music: Understanding Pitch
Math and Music: Understanding Pitch Gareth E. Roberts Department of Mathematics and Computer Science College of the Holy Cross Worcester, MA Topics in Mathematics: Math and Music MATH 110 Spring 2018 March
More informationMassachusetts Institute of Technology Dept. of Electrical Engineering and Computer Science Fall Semester, Introduction to EECS 2
Massachusetts Institute of Technology Dept. of Electrical Engineering and Computer Science Fall Semester, 2006 6.082 Introduction to EECS 2 Lab #2: Time-Frequency Analysis Goal:... 3 Instructions:... 3
More informationChapter 3 Data and Signals 3.1
Chapter 3 Data and Signals 3.1 Copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display. Note To be transmitted, data must be transformed to electromagnetic signals. 3.2
More informationLecture 13. Introduction to OFDM
Lecture 13 Introduction to OFDM Ref: About-OFDM.pdf Orthogonal frequency division multiplexing (OFDM) is well-known to be effective against multipath distortion. It is a multicarrier communication scheme,
More information8.3 Basic Parameters for Audio
8.3 Basic Parameters for Audio Analysis Physical audio signal: simple one-dimensional amplitude = loudness frequency = pitch Psycho-acoustic features: complex A real-life tone arises from a complex superposition
More informationAnalog Synthesizer: Functional Description
Analog Synthesizer: Functional Description Documentation and Technical Information Nolan Lem (2013) Abstract This analog audio synthesizer consists of a keyboard controller paired with several modules
More informationInterpolation Error in Waveform Table Lookup
Carnegie Mellon University Research Showcase @ CMU Computer Science Department School of Computer Science 1998 Interpolation Error in Waveform Table Lookup Roger B. Dannenberg Carnegie Mellon University
More informationMusical Acoustics, C. Bertulani. Musical Acoustics. Lecture 13 Timbre / Tone quality I
1 Musical Acoustics Lecture 13 Timbre / Tone quality I Waves: review 2 distance x (m) At a given time t: y = A sin(2πx/λ) A -A time t (s) At a given position x: y = A sin(2πt/t) Perfect Tuning Fork: Pure
More informationDemonstrate understanding of wave systems. Demonstrate understanding of wave systems. Achievement Achievement with Merit Achievement with Excellence
Demonstrate understanding of wave systems Subject Reference Physics 3.3 Title Demonstrate understanding of wave systems Level 3 Credits 4 Assessment External This achievement standard involves demonstrating
More information2. When is an overtone harmonic? a. never c. when it is an integer multiple of the fundamental frequency b. always d.
PHYSICS LAPP RESONANCE, MUSIC, AND MUSICAL INSTRUMENTS REVIEW I will not be providing equations or any other information, but you can prepare a 3 x 5 card with equations and constants to be used on the
More informationFrom Last Time Wave Properties. Description of a Wave. Water waves? Water waves occur on the surface. They are a kind of transverse wave.
From Last Time Wave Properties Amplitude is the maximum displacement from the equilibrium position Wavelength,, is the distance between two successive points that behave identically Period: time required
More informationKey Vocabulary: Wave Interference Standing Wave Node Antinode Harmonic Destructive Interference Constructive Interference
Key Vocabulary: Wave Interference Standing Wave Node Antinode Harmonic Destructive Interference Constructive Interference 1. Work with two partners. Two will operate the Slinky and one will record the
More informationDefinitions. Spectrum Analyzer
SIGNAL ANALYZERS Spectrum Analyzer Definitions A spectrum analyzer measures the magnitude of an input signal versus frequency within the full frequency range of the instrument. The primary use is to measure
More informationAnalysis and design of filters for differentiation
Differential filters Analysis and design of filters for differentiation John C. Bancroft and Hugh D. Geiger SUMMARY Differential equations are an integral part of seismic processing. In the discrete computer
More informationWaves & Energy Transfer. Introduction to Waves. Waves are all about Periodic Motion. Physics 11. Chapter 11 ( 11-1, 11-7, 11-8)
Waves & Energy Transfer Physics 11 Introduction to Waves Chapter 11 ( 11-1, 11-7, 11-8) Waves are all about Periodic Motion. Periodic motion is motion that repeats after a certain period of time. This
More informationTerminology (1) Chapter 3. Terminology (3) Terminology (2) Transmitter Receiver Medium. Data Transmission. Direct link. Point-to-point.
Terminology (1) Chapter 3 Data Transmission Transmitter Receiver Medium Guided medium e.g. twisted pair, optical fiber Unguided medium e.g. air, water, vacuum Spring 2012 03-1 Spring 2012 03-2 Terminology
More information