FOURIER-BASED METHODS FOR THE SPECTRAL ANALYSIS OF MUSICAL SOUNDS. Sylvain Marchand
|
|
- Kerrie Cummings
- 6 years ago
- Views:
Transcription
1 FOUIE-BAS METHOS FO THE SPECTAL ANALYSIS OF MUSICAL SOUNS Sylvain Marchand Univerity of Bret, Lab-STICC CNS UM 628, 292 Bret, Brittany, France ABSTACT When dealing with muical ound, the hort-time Fourier tranform prevail and inuoid play a key role, according to both acoutic (vibrating mode) and pychoacoutic (pure tone). The value obtained when decompoing the ignal on the time-frequency atom are uually aigned to their geometrical center, leading to etimation error for the inuoidal parameter. To correct thi, one can exploit the amplitude or phae information, ue the derivative of the analyi window, or thoe of the audio ignal. Thi lead to three method (phae vocoder, pectral reaignment, derivative algorithm) equally efficient: they are in fact different formulation of the bet analyi method baed on the Fourier pectrum. Index Term Sound analyi, inuoidal modeling, Fourier tranform, pectral reaignment 1. INTOUCTION Additive ynthei can be conidered a a pectrum modeling technique. It i originally rooted in Fourier theorem, which tate that any periodic function can be modeled a a um of inuoid at variou amplitude and harmonically related frequencie. Sinuoidal modeling conit in conidering the trajectorie in time of the amplitude and frequency parameter of each inuoid preent in the ound. Thi wa propoed by McAulay and Quatieri [1] for peech ignal and by Smith and Serra [2] for muical ound. Sinuoidal modeling lead to meaningful ound repreentation, uitable e.g. for audio effect (time tretching, pitch hifting, etc.), audio coding, ource eparation, or muic trancription. The inuoidal model being parametric, an important problem i to be able to etimate the model parameter a accurately a poible, to get high quality ound. In thi paper, we focu on etimator baed on the Fourier pectrum and well-uited for muical ound, although other approache exit (e.g. ee [3]) and other ignal could alo be conidered. After a preentation of inuoidal modeling in Section 2, Section 3 decribe three method for the etimation of the inuoidal parameter, and Section 4 how that they are all efficient in practice and in fact equivalent in theory: they are different formulation of the bet Fourier-baed analyi method. Together with Lagrange we tudied in [4] the equivalence of thee etimator in the tationary cae and only for one Thi reearch wa partly upported by the AN agency, eam project (AN-9-CO-6). inuoidal parameter (the frequency). Thi paper can be regarded a an extenion, where pectral reaignment play a central role. 2. SPECTAL SOUN MOELING 2.1. Sinuoidal Modeling Let u conider here the inuoidal model under it mot general expreion, which i a um of complex inuoid / exponential (the partial) with low time-varying amplitude a p and non-harmonically related frequencie ω p (defined a the firt derivative of the phae φ p ). The reulting ignal i thu (t) = P a p (t) exp(jφ p (t)) (1) p=1 where P i the number of partial. Since thi paper focue on the tatitical quality of the parameter etimator rather than their frequency reolution, the ignal model i reduced to only one partial (P = 1). The ubcript notation for the partial i then uele. Let u alo define Π a being the value of a given parameter Π at time, correponding to the center of the analyi frame. The ignal i then (t) = exp (λ + µ t) +j (φ + ω t) (2) }{{}}{{} λ(t)=log(a(t)) φ(t) where µ (the amplitude modulation) i the derivative of λ (the log-amplitude), and ω (the frequency) i the derivative of φ (the phae). Thu, the log-amplitude and the phae are modeled by polynomial of degree 1, which can be viewed either a truncated Taylor expanion of more complicated amplitude and frequency modulation (e.g. tremolo / vibrato), or either a an extenion of the tationary cae where µ = Spectral Analyi The main problem we have to tackle now i the etimation of the model parameter. Thi can be achieved, a in the tationary cae, by uing the hort-time Fourier tranform (STFT): = (τ)w(τ t) exp ( jω(τ t)) dτ (3) where S w i the hort-time pectrum of the ignal.
2 Amplitude Magnitude (db) Frequency (Hz) Fig. 1. Magnitude pectrum of a harmonic ound, reulting from a fat (dicrete) Fourier tranform (FFT). Each partial caue a peak (e.g. ee ) and will be conidered individually Frequency (bin) Fig. 2. Zoom on one peak: The analyzed inuoid (dahed line) i oberved from the pectrum of the analyi window through the dicrete frequencie (olid line) of the Fourier tranform, leading to everal bin with ignificant energy ( ). Note that, a in [], we ue here a lightly modified definition of the STFT. Indeed we let the time reference lide with the window, which i alo the cae in practice when the STFT i implemented uing a liding fat Fourier tranform (FFT). For the ake of implicity, all the mathematical derivation will be done in the continuou domain. However, in practice the ignal are dicrete (with ome ampling frequency F ). The Fourier tranform (FT) will then be replaced by it dicrete verion (FT) of ize N, and the time will be expreed in ample (ample n being at time n/f ) and the frequency in bin (bin m being at frequency mf /N). Many intrumental ound are harmonic, meaning that the frequencie of the partial are multiple of ome fundamental frequency (related to our perception of pitch). The magnitude of the Fourier pectrum exhibit then a erie of peak (ee Fig. 1). Each peak i a local maximum m in the magnitude pectrum and correpond to ome partial. S w involve an analyi window w, uually ymmetric and band-limited in uch a way that for any frequency correponding to one pecific partial, the influence of the other partial can be neglected (in the general cae when P > 1). In the tationary cae (µ = ), the pectrum of the analyi window get imply centered on the frequency ω and multiplied by the complex amplitude = a exp(jφ ) = exp(λ + jφ ), (4) a hown in Fig. 2, which can be regarded a a zoom on one of the peak of the preceding figure. In the non-tationary cae however, conidering Equation (3) at etimation time, we ee that get multiplied by Γ w (ω ω, µ ) where Γ w (ω, µ ) = w(t) exp (µ t + jωt) dt. () In the pecial cae of uing a Gauian window for w, an analytic formula can be derived. Ele, it i alway poible to compute Γ w directly from Equation (). Once the etimated amplitude modulation ˆµ and frequency ˆω are known, the amplitude and phae parameter can eventually be etimated ince ŝ = S w (, ω m ) Γ w (ˆω ω m, ˆµ ), (6) where ω m i the (dicrete) frequency of the local maximum of the magnitude pectrum where the partial i detected. 3. SINUSOIAL ESTIMATION The problem i yet to etimate the amplitude modulation and frequency parameter. In thi ection, we preent 3 method providing etimation function ˆµ and ˆω for thee parameter. In practice, for each detected partial at time (center of the analyi frame) and (dicrete) frequency ω m, the etimate of it parameter are given by: ˆµ = ˆµ(, ω m ), (7) ˆω = ˆω(, ω m ). (8) 3.1. ifference Method (Phae Vocoder) The Fourier-baed approach tarted together with computer muic, about year ago. The phae vocoder introduced by Flanagan and Golden [6] wa already uing the phae of the Fourier pectrum to etimate the frequency of the partial, and more preciely the phae difference of conecutive pectra [7]. Thi imple yet efficient difference approach wa generalized recently to the non-tationary cae [8].
3 Thu i the pectrum of the frame centered at the deired (dicrete) etimation time, and let S w (ω) = S w (t 1/F, ω) be it left (previou, i.e. one ample before) and right (next, i.e. one ample after) neighboring pectra, repectively (F denoting the ampling frequency). Since the log-amplitude and phae difference correpond to the real and imaginary part of the logarithm of pectral ratio, repectively, let u define: λ (S 1, S 2 ) = log S 1 log S 2 = (log (S 1 /S 2 )), (9) φ (S 1, S 2 ) = S 1 S 2 = I (log (S 1 /S 2 )) (1) (S 1 and S 2 denoting two arbitrary complex pectra). Since we can meaure the amplitude of the pectra, we can compute the left and right etimate of the amplitude modulation, and retain their mean a the final etimation: µ = λ (S w, Sw ) F, µ + = λ (S w +, S w ) F, ˆµ = (µ + µ + )/2. (11) Similarly, with the meaured phae of the pectra, we can compute an etimation of the intantaneou frequency: ω = φ (S w, S w ) F, ω + = φ (S w +, S w ) F, ˆω = (ω + ω + )/2. (12) In practice, we are looking for poitive frequencie and ince the phae i meaured modulo 2π, after each call to the φ function we mut apply the phae unwrapping procedure of the phae vocoder, i.e. adding 2π to the reult if lower than Spectral eaignment eaignment wa firt propoed by Kodera et al. [9] and wa generalized by Auger and Flandrin [1] to improve timefrequency repreentation. Uually, the value obtained when decompoing the ignal on the time-frequency atom are aigned to the geometrical center of the cell (center of the analyi window and bin of the Fourier tranform). The reaignment method aign each value to the center of gravity of the cell energy. The method ue the knowledge of the firt derivative w obtained by analytic differentiation of the analyi window w in order to adjut the frequency inide the Fourier tranform bin. Thi approach wa generalized for the amplitude modulation in the non-tationary cae (ee []). The complex hort-time pectrum reulting from the STFT i, in the polar form: = a(t, ω) exp (jφ(t, ω)) (13) where the intantaneou amplitude a and phae φ are realvalued function of time t and frequency ω. By conidering Equation (3), we can eaily derive: t log (S w(t, ω)) = jω S w(t, ω) (14) where w denote the derivative of w. Then, ince the amplitude modulation (rep. frequency) i the derivative of the amplitude (rep. phae), from Equation (13) and (14), we obtain the reaigned parameter: ˆµ(t, ω) = t (log (S w(t, ω))) Sw (t, ω) =, (1) ˆω(t, ω) = t I (log (S w(t, ω))) Sw (t, ω) = ω I. (16) 3.3. erivative Method Together with eainte-catherine in [11], we propoed to ue the ignal derivative to etimate the inuoidal parameter in the tationary cae; and with epalle in [], we generalized thi derivative method to the non-tationary cae. Indeed, conidering Equation (2), ince the derivative of an exponential i an exponential, we have: (t) = (µ + jω ) (t) and thu = µ and I = ω. For thi method to work in the cae of a ignal made of everal partial, we have to witch to the pectral domain and define: S ˆµ = w, (17) S w S ˆω = I w, (18) where S w i the hort-time pectrum of the ignal derivative. A hown in [], in practice thi (dicrete) derivative can be obtained by convolving the dicrete ignal by the following differentiator filter: S w ( 1) n h[n] = F for n, and h() = (19) n of infinite time upport. Thu, in practice, we multiply h by a (finite-length) Hann window. Thi reult in a high-pa filter, and can lead to etimation problem in the high frequencie (above approx. 3/4 of the Nyquit frequency), fortunately above the audible limit (16kHz). So, for each partial m detected in the (dicrete) Fourier pectrum at time t and frequency ω m, together with Equation (7) and (8), we have now 3 way to etimate the amplitude modulation and frequency parameter of the partial: the difference approach with Equation (11)-(12), the reaignment approach with Equation (1)-(16), the derivative approach with Equation (17)-(18). Once thee parameter are known, the other (amplitude and phae) can be etimated in turn uing Equation (6).
4 4. COMPAING THE ESTIMATOS Now the quetion i: Which i the bet approach? 4.1. Experimental eult To quantitatively evaluate the preciion of thee approache for the etimation of all the model parameter, we ran the ame experiment a in [, 8]. We conider dicrete-time ignal, with ampling rate F, each coniting of 1 complex exponential generated according to Equation (2) with an initial amplitude a = 1, and mixed with a Gauian white noie. In our experiment, we et the ampling frequency F = 441Hz, the FFT ize N = 11, and the ignal-to-noie ratio (SN) goe from 2dB to +1dB by tep of db. For each SN and for each analyi method, we tet everal parameter combination: 99 frequencie (ω ) linearly ditributed in the (, 3F /8) interval, 9 phae (φ ) linearly ditributed in ( π, +π), and amplitude modulation (µ ) linearly ditributed in [ 1, +1]. For the analyi window w, we ue the ymmetric Hann window. We focu on the variance of the etimation error over thi tet et (the mean being zero for unbiaed etimator). We then compare the difference method (), the reaignment method (), and two variant of the derivative method: The etimated derivative method (), where the derivative i etimated uing the differentiator filter of Equation (19) of ize 123; and the theoretic derivative method (), where the exact derivative i ued for ince it i analytically known for the tet ignal. The reult of the method can be regarded a the bet performance the method could achieve, at the expene of a longer differentiator filter though. When looking at the reult of thee experiment (ee Fig. 3), we ee that all thee method are very efficient, cloe to the the Cramér-ao bound (), which are the limit to the bet poible performance achievable by an unbiaed etimator given a data et (ee []). In the high SN, the method i biaed becaue of the approximation of the derivative by the finite-length differentiator filter. Applying the pectral reaignment on the dicrete pectrum caue a bia, a noticed by Hainworth [12], degrading the performance of the method. Perhap urpriingly, the implet method i the mot efficient Theoretical Equivalence eaignment and erivative In [4], the reaignment (Section 3.2) and derivative (Section 3.3) method were proven to be theoretically equivalent, at leat a regard the etimation of the frequency in the tationary cae. In [], we generalized the proof of the equivalence to the non-tationary cae, and for the etimation of both the frequency and the amplitude modulation. More preciely, we introduce ρ = τ t which give another (equivalent) expreion for the STFT (ee Equation (3)): = (t + ρ)w(ρ) exp ( jωρ) dρ (2) from which we can derive t log (S w(t, ω)) = S w(t, ω). (21) By conidering Equation (21) intead of Equation (14) in Section 3.2 (reaignment approach), we would have obtained the equation of Section 3.3 (derivative approach). Thu the reaignment approach i equivalent to the derivative approach, at leat in the continuou cae. A different proof, baed on integration by part, can be found in [13] eaignment and ifference It turn out that the reaignment approach i alo equivalent to the difference approach in the dicrete cae. Indeed, in the phae vocoder approach the parameter are etimated by firt-order difference, approximating the differentiation of the pectrum of Equation (1)-(16) in the dicrete-time cae.. CONCLUSIONS Although the three approache we preented in Section 3 are equivalent in theory, the mall difference oberved in practice in Section 4.1 are due to a bia of the reaignment method in the dicrete cae, and of the derivative method when uing a finite-length differentiator filter. The mot efficient approach turn out to be the implet one, baed on firt-order difference, i.e. a rather crude approximation of differentiation... But further invetigation are neceary, ince in theory the reaignment method (rep. the derivative method) require the analyi window (rep. time ignal) to be differentiable, which are a priori different condition. 6. EFEENCES [1]. J. McAulay and T. F. Quatieri, Speech Analyi/Synthei Baed on a Sinuoidal epreentation, IEEE Tran. on Acou., Speech, and Sig. Proc., vol. 34, no. 4, pp , [2] J. O. Smith III and X. Serra, PASHL: An Analyi/Synthei Program for Non-Harmonic Sound baed on a Sinuoidal epreentation, in Proc. Int. Computer Muic Conf., 1987, pp [3]. Badeau, G. ichard, and B. avid, Performance of ESPIT for Etimating Mixture of Complex Exponential Modulated by Polynomial, IEEE Tran. on Sig. Proc., vol. 6, no. 2, pp , 28. [4] S. Marchand and M. Lagrange, On the Equivalence of Phae-Baed Method for the Etimation of Intantaneou Frequency, in Proc. 14th European Conf. on Sig. Proc., 26. [] S. Marchand and Ph. epalle, Generalization of the erivative Analyi Method to Non-Stationary Sinuoidal Modeling, in Proc. Int. Conf. on igital Audio Effect, 28, pp
5 variance of the error (log1 cale) -1 etimation of the amplitude variance of the error (log1 cale) -1 etimation of the amplitude modulation ignal-to-noie ratio (db) (a) ignal-to-noie ratio (db) (b) variance of the error (log1 cale) -1 etimation of the phae variance of the error (log1 cale) -1 etimation of the frequency ignal-to-noie ratio (db) (c) ignal-to-noie ratio (db) (d) Fig. 3. Etimation error a function of the SN for the amplitude (a), amplitude modulation (b), phae (c), and frequency (d), with the reaignment (), difference (), theoretical derivative (), and etimated derivative () method, and with comparion to the Cramér-ao Bound (). [6] J. L. Flanagan and. M. Golden, Phae vocoder, Bell Sytem Tech. Journal, vol. 4, pp , [7] M. B. olon, The Phae Vocoder: A Tutorial, Computer Muic Journal, vol. 1, no. 4, pp , [8] S. Marchand, The Simplet Analyi Method for Non- Stationary Sinuoidal Modeling, in Proc. Int. Conf. on igital Audio Effect, 212, pp [9] K. Kodera,. Gendrin, and C. de Villedary, Analyi of Time-Varying Signal with Small BT Value, IEEE Tran. on Acou., Speech, and Sig. Proc., vol. 26, no. 1, pp , [1] F. Auger and P. Flandrin, Improving the eadibility of Time-Frequency and Time-Scale epreentation by the eaignment Method, IEEE Tran. on Sig. Proc., vol. 43, no., pp , 199. [11] M. eainte-catherine and S. Marchand, High Preciion Fourier Analyi of Sound Uing Signal erivative, Journal of the AES, vol. 48, no. 7/8, pp , 2. [12] S. W. Hainworth, Technique for the Automated Analyi of Muical Audio, Ph.. thei, Univerity of Cambridge, United Kingdom, 23. [13] X. Wen and M. Sandler, Note on Model-Baed Non- Stationary Sinuoid Etimation Method Uing erivative, in Proc. Int. Conf. on igital Audio Effect, 29.
Advanced Digital Signal Processing. Stationary/nonstationary signals. Time-Frequency Analysis... Some nonstationary signals. Time-Frequency Analysis
Advanced Digital ignal Proceing Prof. Nizamettin AYDIN naydin@yildiz.edu.tr Time-Frequency Analyi http://www.yildiz.edu.tr/~naydin 2 tationary/nontationary ignal Time-Frequency Analyi Fourier Tranform
More informationDetermination of the local contrast of interference fringe patterns using continuous wavelet transform
Determination of the local contrat of interference fringe pattern uing continuou wavelet tranform Jong Kwang Hyok, Kim Chol Su Intitute of Optic, Department of Phyic, Kim Il Sung Univerity, Pyongyang,
More informationChapter 2 Sampling and Quantization. In order to investigate sampling and quantization, the difference between analog
Chapter Sampling and Quantization.1 Analog and Digital Signal In order to invetigate ampling and quantization, the difference between analog and digital ignal mut be undertood. Analog ignal conit of continuou
More informationLecture 10 Filtering: Applied Concepts
Lecture Filtering: Applied Concept In the previou two lecture, you have learned about finite-impule-repone (FIR) and infinite-impule-repone (IIR) filter. In thee lecture, we introduced the concept of filtering
More informationDesign of Digital Filters
Deign of Digital Filter Paley-Wiener Theorem [ ] ( ) If h n i a caual energy ignal, then ln H e dω< B where B i a finite upper bound. One implication of the Paley-Wiener theorem i that a tranfer function
More informationSampling and the Discrete Fourier Transform
Sampling and the Dicrete Fourier Tranform Sampling Method Sampling i mot commonly done with two device, the ample-and-hold (S/H) and the analog-to-digital-converter (ADC) The S/H acquire a CT ignal at
More informationGain and Phase Margins Based Delay Dependent Stability Analysis of Two- Area LFC System with Communication Delays
Gain and Phae Margin Baed Delay Dependent Stability Analyi of Two- Area LFC Sytem with Communication Delay Şahin Sönmez and Saffet Ayaun Department of Electrical Engineering, Niğde Ömer Halidemir Univerity,
More informationEE 477 Digital Signal Processing. 4 Sampling; Discrete-Time
EE 477 Digital Signal Proceing 4 Sampling; Dicrete-Time Sampling a Continuou Signal Obtain a equence of ignal ample uing a periodic intantaneou ampler: x [ n] = x( nt ) Often plot dicrete ignal a dot or
More informationThe Influence of the Load Condition upon the Radial Distribution of Electromagnetic Vibration and Noise in a Three-Phase Squirrel-Cage Induction Motor
The Influence of the Load Condition upon the Radial Ditribution of Electromagnetic Vibration and Noie in a Three-Phae Squirrel-Cage Induction Motor Yuta Sato 1, Iao Hirotuka 1, Kazuo Tuboi 1, Maanori Nakamura
More informationDigital Control System
Digital Control Sytem - A D D A Micro ADC DAC Proceor Correction Element Proce Clock Meaurement A: Analog D: Digital Continuou Controller and Digital Control Rt - c Plant yt Continuou Controller Digital
More informationFinding the location of switched capacitor banks in distribution systems based on wavelet transform
UPEC00 3t Aug - 3rd Sept 00 Finding the location of witched capacitor bank in ditribution ytem baed on wavelet tranform Bahram nohad Shahid Chamran Univerity in Ahvaz bahramnohad@yahoo.com Mehrdad keramatzadeh
More informationLecture 7: Testing Distributions
CSE 5: Sublinear (and Streaming) Algorithm Spring 014 Lecture 7: Teting Ditribution April 1, 014 Lecturer: Paul Beame Scribe: Paul Beame 1 Teting Uniformity of Ditribution We return today to property teting
More informationR. W. Erickson. Department of Electrical, Computer, and Energy Engineering University of Colorado, Boulder
R. W. Erickon Department of Electrical, Computer, and Energy Engineering Univerity of Colorado, Boulder Cloed-loop buck converter example: Section 9.5.4 In ECEN 5797, we ued the CCM mall ignal model to
More informationinto a discrete time function. Recall that the table of Laplace/z-transforms is constructed by (i) selecting to get
Lecture 25 Introduction to Some Matlab c2d Code in Relation to Sampled Sytem here are many way to convert a continuou time function, { h( t) ; t [0, )} into a dicrete time function { h ( k) ; k {0,,, }}
More informationSource slideplayer.com/fundamentals of Analytical Chemistry, F.J. Holler, S.R.Crouch. Chapter 6: Random Errors in Chemical Analysis
Source lideplayer.com/fundamental of Analytical Chemitry, F.J. Holler, S.R.Crouch Chapter 6: Random Error in Chemical Analyi Random error are preent in every meaurement no matter how careful the experimenter.
More informationSuggested Answers To Exercises. estimates variability in a sampling distribution of random means. About 68% of means fall
Beyond Significance Teting ( nd Edition), Rex B. Kline Suggeted Anwer To Exercie Chapter. The tatitic meaure variability among core at the cae level. In a normal ditribution, about 68% of the core fall
More informationSpring 2014 EE 445S Real-Time Digital Signal Processing Laboratory. Homework #0 Solutions on Review of Signals and Systems Material
Spring 4 EE 445S Real-Time Digital Signal Proceing Laboratory Prof. Evan Homework # Solution on Review of Signal and Sytem Material Problem.. Continuou-Time Sinuoidal Generation. In practice, we cannot
More informationLecture 8: Period Finding: Simon s Problem over Z N
Quantum Computation (CMU 8-859BB, Fall 205) Lecture 8: Period Finding: Simon Problem over Z October 5, 205 Lecturer: John Wright Scribe: icola Rech Problem A mentioned previouly, period finding i a rephraing
More informationCHAPTER 4 DESIGN OF STATE FEEDBACK CONTROLLERS AND STATE OBSERVERS USING REDUCED ORDER MODEL
98 CHAPTER DESIGN OF STATE FEEDBACK CONTROLLERS AND STATE OBSERVERS USING REDUCED ORDER MODEL INTRODUCTION The deign of ytem uing tate pace model for the deign i called a modern control deign and it i
More informationRECURSIVE VERSIONS OF THE LEVENBERG-MARQUARDT REASSIGNED SPECTROGRAM AND OF THE SYNCHROSQUEEZED STFT
RECURSIVE VERSIONS OF THE LEVENBERG-MARQUARDT REASSIGNED SPECTROGRAM AND OF THE SYNCHROSQUEEZED STFT Dominique Fourer 1, Françoi Auger 1, Patrick Flandrin 2 1 LUNAM Univerity, IREENA, Saint-Nazaire, France
More informationSocial Studies 201 Notes for March 18, 2005
1 Social Studie 201 Note for March 18, 2005 Etimation of a mean, mall ample ize Section 8.4, p. 501. When a reearcher ha only a mall ample ize available, the central limit theorem doe not apply to the
More informationRaneNote BESSEL FILTER CROSSOVER
RaneNote BESSEL FILTER CROSSOVER A Beel Filter Croover, and It Relation to Other Croover Beel Function Phae Shift Group Delay Beel, 3dB Down Introduction One of the way that a croover may be contructed
More informationComparing Means: t-tests for Two Independent Samples
Comparing ean: t-tet for Two Independent Sample Independent-eaure Deign t-tet for Two Independent Sample Allow reearcher to evaluate the mean difference between two population uing data from two eparate
More informationClustering Methods without Given Number of Clusters
Clutering Method without Given Number of Cluter Peng Xu, Fei Liu Introduction A we now, mean method i a very effective algorithm of clutering. It mot powerful feature i the calability and implicity. However,
More informationChapter 4. The Laplace Transform Method
Chapter 4. The Laplace Tranform Method The Laplace Tranform i a tranformation, meaning that it change a function into a new function. Actually, it i a linear tranformation, becaue it convert a linear combination
More informationBogoliubov Transformation in Classical Mechanics
Bogoliubov Tranformation in Claical Mechanic Canonical Tranformation Suppoe we have a et of complex canonical variable, {a j }, and would like to conider another et of variable, {b }, b b ({a j }). How
More informationSolutions. Digital Control Systems ( ) 120 minutes examination time + 15 minutes reading time at the beginning of the exam
BSc - Sample Examination Digital Control Sytem (5-588-) Prof. L. Guzzella Solution Exam Duration: Number of Quetion: Rating: Permitted aid: minute examination time + 5 minute reading time at the beginning
More informationTheoretical Computer Science. Optimal algorithms for online scheduling with bounded rearrangement at the end
Theoretical Computer Science 4 (0) 669 678 Content lit available at SciVere ScienceDirect Theoretical Computer Science journal homepage: www.elevier.com/locate/tc Optimal algorithm for online cheduling
More informationGNSS Solutions: What is the carrier phase measurement? How is it generated in GNSS receivers? Simply put, the carrier phase
GNSS Solution: Carrier phae and it meaurement for GNSS GNSS Solution i a regular column featuring quetion and anwer about technical apect of GNSS. Reader are invited to end their quetion to the columnit,
More informationSection Induction motor drives
Section 5.1 - nduction motor drive Electric Drive Sytem 5.1.1. ntroduction he AC induction motor i by far the mot widely ued motor in the indutry. raditionally, it ha been ued in contant and lowly variable-peed
More informationEfficient Methods of Doppler Processing for Coexisting Land and Weather Clutter
Efficient Method of Doppler Proceing for Coexiting Land and Weather Clutter Ça gatay Candan and A Özgür Yılmaz Middle Eat Technical Univerity METU) Ankara, Turkey ccandan@metuedutr, aoyilmaz@metuedutr
More informationSIMON FRASER UNIVERSITY School of Engineering Science ENSC 320 Electric Circuits II. Solutions to Assignment 3 February 2005.
SIMON FRASER UNIVERSITY School of Engineering Science ENSC 320 Electric Circuit II Solution to Aignment 3 February 2005. Initial Condition Source 0 V battery witch flip at t 0 find i 3 (t) Component value:
More informationAlternate Dispersion Measures in Replicated Factorial Experiments
Alternate Diperion Meaure in Replicated Factorial Experiment Neal A. Mackertich The Raytheon Company, Sudbury MA 02421 Jame C. Benneyan Northeatern Univerity, Boton MA 02115 Peter D. Krau The Raytheon
More informationSocial Studies 201 Notes for November 14, 2003
1 Social Studie 201 Note for November 14, 2003 Etimation of a mean, mall ample ize Section 8.4, p. 501. When a reearcher ha only a mall ample ize available, the central limit theorem doe not apply to the
More information7.2 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 281
72 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 28 and i 2 Show how Euler formula (page 33) can then be ued to deduce the reult a ( a) 2 b 2 {e at co bt} {e at in bt} b ( a) 2 b 2 5 Under what condition
More informationDYNAMIC MODELS FOR CONTROLLER DESIGN
DYNAMIC MODELS FOR CONTROLLER DESIGN M.T. Tham (996,999) Dept. of Chemical and Proce Engineering Newcatle upon Tyne, NE 7RU, UK.. INTRODUCTION The problem of deigning a good control ytem i baically that
More informationThe continuous time random walk (CTRW) was introduced by Montroll and Weiss 1.
1 I. CONTINUOUS TIME RANDOM WALK The continuou time random walk (CTRW) wa introduced by Montroll and Wei 1. Unlike dicrete time random walk treated o far, in the CTRW the number of jump n made by the walker
More informationOne Class of Splitting Iterative Schemes
One Cla of Splitting Iterative Scheme v Ciegi and V. Pakalnytė Vilniu Gedimina Technical Univerity Saulėtekio al. 11, 2054, Vilniu, Lithuania rc@fm.vtu.lt Abtract. Thi paper deal with the tability analyi
More informationImproving the Efficiency of a Digital Filtering Scheme for Diabatic Initialization
1976 MONTHLY WEATHER REVIEW VOLUME 15 Improving the Efficiency of a Digital Filtering Scheme for Diabatic Initialization PETER LYNCH Met Éireann, Dublin, Ireland DOMINIQUE GIARD CNRM/GMAP, Météo-France,
More informationON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION. Xiaoqun Wang
Proceeding of the 2008 Winter Simulation Conference S. J. Maon, R. R. Hill, L. Mönch, O. Roe, T. Jefferon, J. W. Fowler ed. ON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION Xiaoqun Wang
More informationChapter 2: Problem Solutions
Chapter 2: Solution Dicrete Time Proceing of Continuou Time Signal Sampling à 2.. : Conider a inuoidal ignal and let u ample it at a frequency F 2kHz. xt 3co000t 0. a) Determine and expreion for the ampled
More informationDETECTION OF ROOM REFLECTIONS FROM A BINAURAL ROOM IMPULSE RESPONSE
Sampo Vea and Tapio Lokki. 2006. Detection of room reflection from a binaural room impule repone. In: Proceeding of the 9th International Conference on Digital Audio Effect (DAFx 2006). Montreal, Canada.
More informationDesign By Emulation (Indirect Method)
Deign By Emulation (Indirect Method he baic trategy here i, that Given a continuou tranfer function, it i required to find the bet dicrete equivalent uch that the ignal produced by paing an input ignal
More informationSMALL-SIGNAL STABILITY ASSESSMENT OF THE EUROPEAN POWER SYSTEM BASED ON ADVANCED NEURAL NETWORK METHOD
SMALL-SIGNAL STABILITY ASSESSMENT OF THE EUROPEAN POWER SYSTEM BASED ON ADVANCED NEURAL NETWORK METHOD S.P. Teeuwen, I. Erlich U. Bachmann Univerity of Duiburg, Germany Department of Electrical Power Sytem
More informationAn estimation approach for autotuning of event-based PI control systems
Acta de la XXXIX Jornada de Automática, Badajoz, 5-7 de Septiembre de 08 An etimation approach for autotuning of event-baed PI control ytem Joé Sánchez Moreno, María Guinaldo Loada, Sebatián Dormido Departamento
More informationMath Skills. Scientific Notation. Uncertainty in Measurements. Appendix A5 SKILLS HANDBOOK
ppendix 5 Scientific Notation It i difficult to work with very large or very mall number when they are written in common decimal notation. Uually it i poible to accommodate uch number by changing the SI
More informationChapter 12 Simple Linear Regression
Chapter 1 Simple Linear Regreion Introduction Exam Score v. Hour Studied Scenario Regreion Analyi ued to quantify the relation between (or more) variable o you can predict the value of one variable baed
More informationLecture 9: Shor s Algorithm
Quantum Computation (CMU 8-859BB, Fall 05) Lecture 9: Shor Algorithm October 7, 05 Lecturer: Ryan O Donnell Scribe: Sidhanth Mohanty Overview Let u recall the period finding problem that wa et up a a function
More informationChapter 4: Applications of Fourier Representations. Chih-Wei Liu
Chapter 4: Application of Fourier Repreentation Chih-Wei Liu Outline Introduction Fourier ranform of Periodic Signal Convolution/Multiplication with Non-Periodic Signal Fourier ranform of Dicrete-ime Signal
More informationOptimal Coordination of Samples in Business Surveys
Paper preented at the ICES-III, June 8-, 007, Montreal, Quebec, Canada Optimal Coordination of Sample in Buine Survey enka Mach, Ioana Şchiopu-Kratina, Philip T Rei, Jean-Marc Fillion Statitic Canada New
More informationLecture 4 Topic 3: General linear models (GLMs), the fundamentals of the analysis of variance (ANOVA), and completely randomized designs (CRDs)
Lecture 4 Topic 3: General linear model (GLM), the fundamental of the analyi of variance (ANOVA), and completely randomized deign (CRD) The general linear model One population: An obervation i explained
More informationLecture 2: The z-transform
5-59- Control Sytem II FS 28 Lecture 2: The -Tranform From the Laplace Tranform to the tranform The Laplace tranform i an integral tranform which take a function of a real variable t to a function of a
More informationLinearteam tech paper. The analysis of fourth-order state variable filter and it s application to Linkwitz- Riley filters
Linearteam tech paper The analyi of fourth-order tate variable filter and it application to Linkwitz- iley filter Janne honen 5.. TBLE OF CONTENTS. NTOCTON.... FOTH-OE LNWTZ-LEY (L TNSFE FNCTON.... TNSFE
More informationJul 4, 2005 turbo_code_primer Revision 0.0. Turbo Code Primer
Jul 4, 5 turbo_code_primer Reviion. Turbo Code Primer. Introduction Thi document give a quick tutorial on MAP baed turbo coder. Section develop the background theory. Section work through a imple numerical
More informationCOPRIME ARRAYS AND SAMPLERS FOR SPACE-TIME ADAPTIVE PROCESSING
COPRIE ARRAYS AND SAPLERS FOR SPACE-IE ADAPIVE PROCESSING Chun-Lin Liu 1 and P. P. Vaidyanathan 2 Dept. of Electrical Engineering, 136-93 California Intitute of echnology, Paadena, CA 91125, USA E-mail:
More informationCHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS
CHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS 8.1 INTRODUCTION 8.2 REDUCED ORDER MODEL DESIGN FOR LINEAR DISCRETE-TIME CONTROL SYSTEMS 8.3
More informationWavelet Analysis in EPR Spectroscopy
Vol. 108 (2005) ACTA PHYSICA POLONICA A No. 1 Proceeding of the XXI International Meeting on Radio and Microwave Spectrocopy RAMIS 2005, Poznań-Bȩdlewo, Poland, April 24 28, 2005 Wavelet Analyi in EPR
More informationEE 508 Lecture 16. Filter Transformations. Lowpass to Bandpass Lowpass to Highpass Lowpass to Band-reject
EE 508 Lecture 6 Filter Tranformation Lowpa to Bandpa Lowpa to Highpa Lowpa to Band-reject Review from Lat Time Theorem: If the perimeter variation and contact reitance are neglected, the tandard deviation
More informationProperties of Z-transform Transform 1 Linearity a
Midterm 3 (Fall 6 of EEG:. Thi midterm conit of eight ingle-ided page. The firt three page contain variou table followed by FOUR eam quetion and one etra workheet. You can tear out any page but make ure
More informationEE 4443/5329. LAB 3: Control of Industrial Systems. Simulation and Hardware Control (PID Design) The Inverted Pendulum. (ECP Systems-Model: 505)
EE 4443/5329 LAB 3: Control of Indutrial Sytem Simulation and Hardware Control (PID Deign) The Inverted Pendulum (ECP Sytem-Model: 505) Compiled by: Nitin Swamy Email: nwamy@lakehore.uta.edu Email: okuljaca@lakehore.uta.edu
More informationFinite Element Analysis of a Fiber Bragg Grating Accelerometer for Performance Optimization
Finite Element Analyi of a Fiber Bragg Grating Accelerometer for Performance Optimization N. Baumallick*, P. Biwa, K. Dagupta and S. Bandyopadhyay Fiber Optic Laboratory, Central Gla and Ceramic Reearch
More informationTHE EXPERIMENTAL PERFORMANCE OF A NONLINEAR DYNAMIC VIBRATION ABSORBER
Proceeding of IMAC XXXI Conference & Expoition on Structural Dynamic February -4 Garden Grove CA USA THE EXPERIMENTAL PERFORMANCE OF A NONLINEAR DYNAMIC VIBRATION ABSORBER Yung-Sheng Hu Neil S Ferguon
More informationLecture 21. The Lovasz splitting-off lemma Topics in Combinatorial Optimization April 29th, 2004
18.997 Topic in Combinatorial Optimization April 29th, 2004 Lecture 21 Lecturer: Michel X. Goeman Scribe: Mohammad Mahdian 1 The Lovaz plitting-off lemma Lovaz plitting-off lemma tate the following. Theorem
More informationEE/ME/AE324: Dynamical Systems. Chapter 8: Transfer Function Analysis
EE/ME/AE34: Dynamical Sytem Chapter 8: Tranfer Function Analyi The Sytem Tranfer Function Conider the ytem decribed by the nth-order I/O eqn.: ( n) ( n 1) ( m) y + a y + + a y = b u + + bu n 1 0 m 0 Taking
More informationCopyright 1967, by the author(s). All rights reserved.
Copyright 1967, by the author(). All right reerved. Permiion to make digital or hard copie of all or part of thi work for peronal or claroom ue i granted without fee provided that copie are not made or
More informationRoot Locus Diagram. Root loci: The portion of root locus when k assume positive values: that is 0
Objective Root Locu Diagram Upon completion of thi chapter you will be able to: Plot the Root Locu for a given Tranfer Function by varying gain of the ytem, Analye the tability of the ytem from the root
More informationReal Sources (Secondary Sources) Phantom Source (Primary source) LS P. h rl. h rr. h ll. h lr. h pl. h pr
Ecient frequency domain ltered-x realization of phantom ource iet C.W. ommen, Ronald M. Aart, Alexander W.M. Mathijen, John Gara, Haiyan He Abtract A phantom ound ource i a virtual ound image which can
More informationQuestion 1 Equivalent Circuits
MAE 40 inear ircuit Fall 2007 Final Intruction ) Thi exam i open book You may ue whatever written material you chooe, including your cla note and textbook You may ue a hand calculator with no communication
More informationAnnex-A: RTTOV9 Cloud validation
RTTOV-91 Science and Validation Plan Annex-A: RTTOV9 Cloud validation Author O Embury C J Merchant The Univerity of Edinburgh Intitute for Atmo. & Environ. Science Crew Building King Building Edinburgh
More informationBy Xiaoquan Wen and Matthew Stephens University of Michigan and University of Chicago
Submitted to the Annal of Applied Statitic SUPPLEMENTARY APPENDIX TO BAYESIAN METHODS FOR GENETIC ASSOCIATION ANALYSIS WITH HETEROGENEOUS SUBGROUPS: FROM META-ANALYSES TO GENE-ENVIRONMENT INTERACTIONS
More informationRoadmap for Discrete-Time Signal Processing
EE 4G Note: Chapter 8 Continuou-time Signal co(πf Roadmap for Dicrete-ime Signal Proceing.5 -.5 -..4.6.8..4.6.8 Dicrete-time Signal (Section 8.).5 -.5 -..4.6.8..4.6.8 Sampling Period econd (or ampling
More informationRELIABILITY OF REPAIRABLE k out of n: F SYSTEM HAVING DISCRETE REPAIR AND FAILURE TIMES DISTRIBUTIONS
www.arpapre.com/volume/vol29iue1/ijrras_29_1_01.pdf RELIABILITY OF REPAIRABLE k out of n: F SYSTEM HAVING DISCRETE REPAIR AND FAILURE TIMES DISTRIBUTIONS Sevcan Demir Atalay 1,* & Özge Elmataş Gültekin
More informationCodes Correcting Two Deletions
1 Code Correcting Two Deletion Ryan Gabry and Frederic Sala Spawar Sytem Center Univerity of California, Lo Angele ryan.gabry@navy.mil fredala@ucla.edu Abtract In thi work, we invetigate the problem of
More informationμ + = σ = D 4 σ = D 3 σ = σ = All units in parts (a) and (b) are in V. (1) x chart: Center = μ = 0.75 UCL =
Our online Tutor are available 4*7 to provide Help with Proce control ytem Homework/Aignment or a long term Graduate/Undergraduate Proce control ytem Project. Our Tutor being experienced and proficient
More informationNonlinear Single-Particle Dynamics in High Energy Accelerators
Nonlinear Single-Particle Dynamic in High Energy Accelerator Part 6: Canonical Perturbation Theory Nonlinear Single-Particle Dynamic in High Energy Accelerator Thi coure conit of eight lecture: 1. Introduction
More informationLecture #9 Continuous time filter
Lecture #9 Continuou time filter Oliver Faut December 5, 2006 Content Review. Motivation......................................... 2 2 Filter pecification 2 2. Low pa..........................................
More informationFactor Analysis with Poisson Output
Factor Analyi with Poion Output Gopal Santhanam Byron Yu Krihna V. Shenoy, Department of Electrical Engineering, Neurocience Program Stanford Univerity Stanford, CA 94305, USA {gopal,byronyu,henoy}@tanford.edu
More informationSingular perturbation theory
Singular perturbation theory Marc R. Rouel June 21, 2004 1 Introduction When we apply the teady-tate approximation (SSA) in chemical kinetic, we typically argue that ome of the intermediate are highly
More informationA Non-destructive Approach for Noise Reduction in Time Domain
World Applied Science Journal 6 (1): 53-63, 29 ISSN 1818-4952 IDOSI Publication, 29 A Non-detructive Approach for Noie Reduction in ime Domain A. Zehtabian and H. Haanpour School of Information echnology
More informationRECURSIVE LEAST SQUARES HARMONIC IDENTIFICATION IN ACTIVE POWER FILTERS. A. El Zawawi, K. H. Youssef, and O. A. Sebakhy
RECURSIVE LEAST SQUARES HARMONIC IDENTIFICATION IN ACTIVE POWER FILTERS A. El Zawawi, K. H. Youef, and O. A. Sebakhy Department of Electrical Engineering, Alexandria Univerity, Alexandria 21544, Egypt.P.O.
More informationCompact finite-difference approximations for anisotropic image smoothing and painting
CWP-593 Compact finite-difference approximation for aniotropic image moothing and painting Dave Hale Center for Wave Phenomena, Colorado School of Mine, Golden CO 80401, USA ABSTRACT Finite-difference
More informationA Short Note on Hysteresis and Odd Harmonics
1 A Short Note on yterei and Odd armonic Peter R aoput Senior Reearch Scientit ubocope Pipeline Engineering outon, X 7751 USA pmaoput@varco.com June 19, Abtract hi hort note deal with the exitence of odd
More informationA Study on Simulating Convolutional Codes and Turbo Codes
A Study on Simulating Convolutional Code and Turbo Code Final Report By Daniel Chang July 27, 2001 Advior: Dr. P. Kinman Executive Summary Thi project include the deign of imulation of everal convolutional
More informationRefinements to the Model of a Single Woodwind Instrument Tonehole
Proceeding of 20th International Sympoium on Muic Acoutic (Aociated Meeting of the International Congre on Acoutic) 25-31 Augut 2010, Sydney and Katoomba, Autralia Refinement to the Model of a Single Woodwind
More informationJan Purczyński, Kamila Bednarz-Okrzyńska Estimation of the shape parameter of GED distribution for a small sample size
Jan Purczyńki, Kamila Bednarz-Okrzyńka Etimation of the hape parameter of GED ditribution for a mall ample ize Folia Oeconomica Stetinenia 4()/, 35-46 04 Folia Oeconomica Stetinenia DOI: 0.478/foli-04-003
More informationNon-stationary Random Vibration Analysis of Vehicle with Fractional
13 th National Conference on Mechanim and Machine (NaCoMM7), IISc, Bangalore, India, December 1-13, 7 NaCoMM-7-77 Non-tationary Random Vibration Analyi of Vehicle with Fractional Damping Yunqing Zhang
More informationDigital Control System
Digital Control Sytem Summary # he -tranform play an important role in digital control and dicrete ignal proceing. he -tranform i defined a F () f(k) k () A. Example Conider the following equence: f(k)
More informationA Bluffer s Guide to... Sphericity
A Bluffer Guide to Sphericity Andy Field Univerity of Suex The ue of repeated meaure, where the ame ubject are teted under a number of condition, ha numerou practical and tatitical benefit. For one thing
More informationCoupling of Three-Phase Sequence Circuits Due to Line and Load Asymmetries
Coupling of Three-Phae Sequence Circuit Due to Line and Load Aymmetrie DEGO BELLAN Department of Electronic nformation and Bioengineering Politecnico di Milano Piazza Leonardo da inci 01 Milano TALY diego.ellan@polimi.it
More information0 of the same magnitude. If we don t use an OA and ignore any damping, the CTF is
1 4. Image Simulation Influence of C Spherical aberration break the ymmetry that would otherwie exit between overfocu and underfocu. One reult i that the fringe in the FT of the CTF are generally farther
More informationIII.9. THE HYSTERESIS CYCLE OF FERROELECTRIC SUBSTANCES
III.9. THE HYSTERESIS CYCLE OF FERROELECTRIC SBSTANCES. Work purpoe The analyi of the behaviour of a ferroelectric ubtance placed in an eternal electric field; the dependence of the electrical polariation
More informationLinear predictive coding
Linear predictive coding Thi method combine linear proceing with calar quantization. The main idea of the method i to predict the value of the current ample by a linear combination of previou already recontructed
More informationThe type 3 nonuniform FFT and its applications
Journal of Computational Phyic 206 (2005) 1 5 Short Note The type 3 nonuniform FFT and it application June-Yub Lee a, *, Lelie Greengard b a Department of Mathematic, Ewha Woman Univerity, 11-1 Daehyundong,
More informationUnavoidable Cycles in Polynomial-Based Time-Invariant LDPC Convolutional Codes
European Wirele, April 7-9,, Vienna, Autria ISBN 978--87-4-9 VE VERLAG GMBH Unavoidable Cycle in Polynomial-Baed Time-Invariant LPC Convolutional Code Hua Zhou and Norbert Goertz Intitute of Telecommunication
More informationOverflow from last lecture: Ewald construction and Brillouin zones Structure factor
Lecture 5: Overflow from lat lecture: Ewald contruction and Brillouin zone Structure factor Review Conider direct lattice defined by vector R = u 1 a 1 + u 2 a 2 + u 3 a 3 where u 1, u 2, u 3 are integer
More informationAvoiding Forbidden Submatrices by Row Deletions
Avoiding Forbidden Submatrice by Row Deletion Sebatian Wernicke, Jochen Alber, Jen Gramm, Jiong Guo, and Rolf Niedermeier Wilhelm-Schickard-Intitut für Informatik, niverität Tübingen, Sand 13, D-72076
More information5.5 Application of Frequency Response: Signal Filters
44 Dynamic Sytem Second order lowpa filter having tranfer function H()=H ()H () u H () H () y Firt order lowpa filter Figure 5.5: Contruction of a econd order low-pa filter by combining two firt order
More informationS_LOOP: SINGLE-LOOP FEEDBACK CONTROL SYSTEM ANALYSIS
S_LOOP: SINGLE-LOOP FEEDBACK CONTROL SYSTEM ANALYSIS by Michelle Gretzinger, Daniel Zyngier and Thoma Marlin INTRODUCTION One of the challenge to the engineer learning proce control i relating theoretical
More informationLecture 5 Frequency Response of FIR Systems (III)
EE3054 Signal and Sytem Lecture 5 Frequency Repone of FIR Sytem (III Yao Wang Polytechnic Univerity Mot of the lide included are extracted from lecture preentation prepared by McClellan and Schafer Licene
More informationarxiv: v1 [math.ca] 23 Sep 2017
arxiv:709.08048v [math.ca] 3 Sep 07 On the unit ditance problem A. Ioevich Abtract. The Erdő unit ditance conjecture in the plane ay that the number of pair of point from a point et of ize n eparated by
More information