Real Sources (Secondary Sources) Phantom Source (Primary source) LS P. h rl. h rr. h ll. h lr. h pl. h pr
|
|
- Frank McDowell
- 5 years ago
- Views:
Transcription
1 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 be utilized in many application like tereo bae widening, multimedia, and virtual reality engine. Generation of uch a virtual ound image uing adaptive lter that employ the ltered-x NM algorithm i dicued. uch lter are long and complex for the virtual ource to cover the whole audio frequency range. Due to thi complexity, real-time implementation on a ampleby-ample bai i not poible and block proceing technique, uch a block frequency domain adaptive lter, have to be ued. everal optimized verion of the ltered-x algorithm in frequency domain are derived and the dierent implementation are compared. Keyword Active Noie Cancellation, iltered-x, requency Domain Adaptive ilter, oudpeaker, hantom ound ource, ound Reproduction ound ource i a virtual ource which i perceived by the litener a a ource at a dierent poition than that of the real ource. The baic idea behind phantom ource generation i to alter the loudpeaker' input ignal in uch a way that the ound reure evel () at the litener' eardrum generated by the real loudpeaker i the ame a that would be generated by the phantom ource [4]. x w Real ource (econdary ource) w R R hantom ource (rimary ource) I. Introduction h ll h lr h rl h rr l When audio ignal are played by a normal tereo etup, the ound received by a litener contain both directional and ditance information, from which the litener can locate the poition of the ound ource. Experience how that when a tereo ignal i reproduced by two eparate loudpeaker, the perceived ound image would be between the two loudpeaker if the litener i in the center of the litening area. However, in ome application like TV and monitor, the ditance between the loudpeaker ha to be limited in order to obtain a compact arrangement. To give the loudpeaker bae a virtual broadening, we made ue of the phantom ound ource concept. Aphantom. ommen, A. Mathijen, J. Gara and H. He are with the Eindhoven Univerity of Technology, Building Eh 6.33,.O.Box 513, 56 MB Eindhoven, The Netherland, p.c.w.ommen@ele.tue.nl R. Aart i with the hilip Reearch aboratorie WY81, 5656 AA Eindhoven, The Netherland, rmaart@natlab.reearch.philip.com ig. 1. β r et-up for generating a phantom ource. ig. 1 how a ituation where a phantom ound ource p at angle i generated by two real loudpeaker at angle. w and w R are lter which are ued to alter the input ignal to an R. The deign of the lter w and w R i the main tak in generating a phantom ource. Thee lter repreent, in general, long and complex impule repone which are alo poition dependent. Therefore it i dicult to obtain an analytical olution for thee lter. or thi reaon, long adaptive lter are ued which are implemented very eciently in the frequency domain. Thee lter are applied a follow
2 (ee ig. 1): Initially three loudpeaker are ued, one un-ltered at the poition of the phantom ource, the other two ltered by adaptive lter. White noie i played through loudpeaker p while w and w R are adapted to minimize the meaured at the litener' eardrum. The method decribed above can hence be regarded a an active noie cancellation olution. The coecient of w and w R are then frozen and the third loudpeaker i removed. When the two loudpeaker are then driven through the ;w and ;w R, only the illuion of a third ource reult. II. iltered-x algorithm ig. model a implied (ingle point) verion of the ytem hown in ig. 1, with one primary ( )and one econdary (h ) acoutic path. The coecient of the adaptive lter w are updated by an algorithm that i baed on the Normalizeeat Mean quare (NM) algorithm. Thi algorithm update the adaptive weight vector w in the negative direction of the gradient vector 5. It i hown in [1] that an etimate of thi gradient vector i given by x r, in which x i the input ignal vector and r the reidual ignal and the update equation become: w new = w old ; x x r (1) where i the adaptation coecient which i normalized bythevariance of the input ignal x = Efx [k]g. A can be een in ig., the output ignal of the adaptive lter e i rt ltered by the econdary acoutic path h before it i added in the microphone. or thi reaon, we need a lightly dierent algorithm than that given by (1) the o-called ltered-x algorithm [1]. The update part of thi algorithm ue a ltered verion x f intead of the input ignal x itelf. Thi i achieved by ltering x by an etimate of the econdary acoutic path h a hown in ig.. A. Block requency Domain Adaptive iltering Becaue of the ytem complexity, real-time realization of a phantom ource that cover the whole audio frequency range i almot impoible, with the current technology, if carried out on a ample by ample bai. Therefore, block proceing technique have to x h etimation ig.. x f w Update e h ingle point ytem model. r + + MIC be ued. In thi work, the Block requency Domain Adaptive ilter (BDA) [3] ha been utilized. In thi ubection, the BDA will be decribed to introduce the notation and in ubequent ection, the ltered-x part will be dicued in more detail. ig. 3 how the BDA of the ingle point model hown in ig.. Thi gure i almot equivalent to the original BDA implementation. The only dierence i that the input ignal i ltered, in the box `lteredx', before it enter the update part of the algorithm. σ x X [k] f * Δ R[(k-1)] ig. 3. [(k-1)] filtered-x M r[k-+1] overlap M W[k] econdary ource primary ource X[k] E[k] e[k-+1] Block requency Domain Adaptive ilter. Working on a block bai i reected in the fact that the whole implementation i performed on block of data rather than on ample. Each block contain new ample, with 1. The input ignal ample are converted from erial to parallel (box `/') while the output ample are converted again from paral- h +
3 lel to erial (box `/'). A direct conequence of thi block proceing approach i a proceing delay of ample it take ample before the rt block of e can be produced. On the other hand, thi proceing delay ha to be compenated, for ynchronization reaon, in the primary ource path (box ` '). or a gradient update algorithm uch a the NM conidered here, two main operation have to be computed [1]: a (linear) convolution to perform the ltering of the input ignal x with the M adaptive weight a (linear) correlation, between the input ignal x and the reidual ignal r. The rt calculate the lter output, and the econd calculate an etimate of the gradient 5 that i needed for the update of the adaptive weight. or large lter length M thee operation can be carried out very eciently by uing a block proceing approach, uch a the overlap-ave method [], implemented in the frequency domain. In thi cae, at ourier Tranform (T) are ued to perform the ignal tranformation back and forth between time and frequency. In the following the above mentioned convolution and correlation will be decribed when the overlap-ave method i ued (ee ig. 3). The convolution (performed by the right-hand part of ig. 3) of the `innite' input equence (the input ignal ample ) with the `nite' length equence (the M adaptive weight) i broken down into a equence of convolution between two `nite' equence of length. The `innite' input equence i rt divided into overlapping block. The overlap of M ; 1 ample i needed becaue the M adaptive weight have rt to be hifted in the equence of ample. On the other hand, performing the convolution in frequency domain produce a cyclic convolution reult. Therefore the overlap i ued to obtain a linear convolution reult from thi cyclic one. Together with the fact that every block contain new ample, the choen block length equal = M + ; 1. A depicted in ig. 3, the block with input ignal ample are tranformed to the frequency domain with an T of length (` ') reulting in a block that, for notational reaon, i repreented by the tranformed input ignal vector X[k]. The M adaptive weight are rt augmented by ;M zeroe to obtain block of length, then tranformed to frequency domain reulting in W[k]. The cyclic convolution i now performed by element-wie multiplying the two vector X[k] and W[k] E[k] =X[k] W[k]: () where the ymbol i ued to denote element-wie multiplication. Only the lat ample (in time) from thi cyclic convolution repreent the deired linear convolution reult. Therefore, the reult E[k] i tranformed back to time domain by aninvere T ( ;1 ) and the correct ample are ent to the econdary loudpeaker after going through the parallel B1 to erial operation. The correlation operation between the `innite' equence of input ignal ample x and the `nite' equence of reidual ignal ample r i performed by the left-hand part of ig. 3. The only dierence between the correlation and convolution operation i mirroring in time domain. Thi mirror operation can be performed in the ourier domain by imply applying a complex conjugate (`*') operator reulting in the vector X f [(k;1)] 1. At the ame time, conecutive ample of the reidual ignal are collected, padded with zero and tranformed to the frequency domain reulting in the vector R[(k ; 1)]. The, caled, gradient i now calculated by 5[(k ; 1)] = x X f [(k ; 1)] R[(k ; 1)] in which the caling i done by. ince the correlation i performed in the ourier domain, the reult x i alo cyclic and the gradient vector 5[(k ; 1)] ha to be tranformed back to time domain to elect the M correct gradient coecient, a done in the convolution operation. Thee gradient coecient are then augmented by zero and the reulting vector of length i tranformed to frequency domain and ued to update the adaptive weight. 1 Note that the input ignal to the correlation operation, here the ltered-x ignal X f [k], need to be delayed by for ynchronization reaon and the reidual ignal R[(k;1)] contain the ame delay.
4 B. iltered-x in frequency domain The `ltered-x' operation, that i needed for the update algorithm in ig. 3, i depicted in more detail in ig. 4. The function of thi `ltered-x' block i to l- X [k] f Append overlap A-1 B X[k] H B B A h ig. 4. iltered-x implementation type 1. ter (convolve) the input ignal with, an etimate of, the A weight of the econdary path h. Thi convolution i performed in frequency domain uing the overlap-ave method a explained before. Each block contain new input ignal ample while the egment have an overlap of A ; 1 ample. The reulting egment contain B = A + ; 1 input ignal ample. Thee are tranformed to the frequency domain (X[k]) and the convolution i performed by the element-wie vector multiplication with H which i obtained by padding the impule repone h with B ; A +1zero. The correct ample of the cyclic convolution have to be appended in time domain to A-1 previouly calculated ample to form egment of length. Thee are tranformed to frequency domain and reult in X f [k] thatuedinthe correlation operation decribed in the previou ubection. rom ig. 4 it follow that the `ltered-x' operation in itelf cot four T'. In the following ection we will how that thee four T' can be reduced tep by tep to only one. III. Efficient filtered-x The reduction from four T' to one i done in three tep [5]. The rt i reduction to three T' and obtained by combining more carefully the T' in ig. 3 and ig. 4. The next reduction to two T' tem from the fact that the convolution operation, needed for calculating the ltered-x input ignal, and the correlation operation, needed for the gradient etimate, can be interchanged. The lat reduction to one T i obtained by combining again T operation. tep 1: Combining two T' The rt optimization tep i to combine the two eparate ourier tranform that are ued to calculate the tranformed input ignal vector X[k] of length in ig. 3 and length B in ig. 4. The T B in ig. 4 can then be left out. The reult i depicted in ig. 5. X [k] f Append B1-A+1 H ig. 5. iltered-x implementation type. X[k] Here we have aumed that the length A of the econdary path i maller than the length M of the adaptive weight vector W (A <M). By augmenting the impule repone vector h of the econdary path with more zero than before a length vector i obtained that i tranformed to the frequency domain reulting in H. Thi vector i multiplied element-wie by the tranformed input ignal vector X[k] that wa available in ig. 3. The ret of the operation in ig. 5 are exactly the ame a in ig. 4. The reult i one T le in comparion with type 1 implementation hown in ig. 4. tep : Interchanging convolution/correlation The econd optimization tep i to interchange the order of the convolution operation of the `ltered-x' and the correlation operation. In ig. 3, the overlap-ave method ha to be ued in the ltered-x block ince the input ignal i innitely long. Becaue of the cyclic nature of thi convolution, the correct ample have to be elected in time domain coting one invere T and one forward T a hown in ig. 4 and ig. 5. Noting that both the impule repone (A) and the reidual ignal block () are of `nite' length, it i more ecient to rt calculate the convolution of thee two 'nite' equence. Thi will reult in a 'nite' equence of length A + ; 1. Implementing thi operation in frequency domain, with ourier tranform h A
5 of appropriate length, reult in a frequency domain vector that contain no cyclic (incorrect) reult. Thi frequency domain vector can now be ued traight on in the correlation with the 'innite' input ignal a before. The reult of thi approach i depicted in ig. 6. Thi realization can alo be referred to a 'ltered-r' σ x X[k] * H R[(k-1)] B 3 B 3 overlap M+A- M A h - + W[k] filtered-r r[k-+1] overlap M-1 econdary ource h primary ource e[k-+1] X[k] E[k] σ x H * R[(k-1)] B 3 B 3 M A h filtered-r r[k-+1] - + W[k] econdary ource h primary ource overlap M+A- B 3 A-1 X[k] E[k] e[k-+1] ig. 7. iltered-x implementation type 4. ecient implementation of the `ltered-x' (or `lteredr') algorithm in frequency domain that cot one more T compared to the original BDA. ig. 6. iltered-x implementation type 3. becaue the 'ltered-x' operation i in fact calculated by ltering rt the reidual ignal (`ltered-r'). No extra invere ourier tranform i needed, a in ig. 5. We only have to take care of a correct ourier tranform length B 3 = M ++A;. inally we note that only two T' are needed in comparion with ig. 4 and ig. 5. tep3: Combining two T' inally an equivalent tep a decribed in tep can be ued here. The ourier tranform that i needed in ig. 6 to calculate the tranformed input ignal vector X[k] can be combined with the ourier tranform that i needed at the right hand ide of the ame gure for the convolution with the adaptive weight vector. The reult i depicted in ig. 7. rom thi gure it follow that the output of the adaptive lter produce each block A + ; 1 correct reult from which only ample are ued and A; 1 are not ued. However the total reult i that we have introduced an IV. Complexity comparion In the previou ection, the number of T' in the ltered-x BDA implementation i reduced from 9 to 6. However, the T length are not equal and depend on the parameter A, and M. Therefore it i dicult to compare the complexity of the different type jut by looking at the number of T'. ig. 8 how the complexity of the four implementation decribed above. In thi gure we have ued a xed length for the econdary impule repone (h ) A = 18. urthermore we have choen the block length equal to the adaptive lter length (minu one), thu = M ; 1. We have approximated the complexity by the total number of multiplication needed for the computation of one output ample [5]. The total number of multiplication tem from the ued T' of length B, eachhaving a complexity of (B log ( B 8 )+4) (for real input (T) and real output (IT)), the element-wie multiplication (B ; ) and multiplication by a real calar (B). rom thi comparion it follow that, for the choen parameter
6 : iltered x type : iltered x type : background noie : primary loudpeaker 3: B=51 4: B=48 3: iltered x type 3 4 computational complexity : iltered x type power pectrum (db) number of weight M frequency (khz) ig. 8. Computational complexity of the implementation. et, type ha the lowet complexity when M < 76 (curve 1) and type 4 when M > 76 (curve 4). Thu for application where A i relatively mall compared to the adaptive lter length M (e.g. for noie cancellation application with headphone), type 4 i preferred. However when A i larger than ay half the lter length, type i preferred. V. imulation and experimental reult Experiment have been performed [5] on a one point noie cancellation ytem. The et-up for thi ytem conit of a primary loudpeaker, a econdary loudpeaker and a microphone, all placed in an anechoic room, where the primary loudpeaker i placed further away from the microphone than the econdary loudpeaker. The etimated impule repone h i xed during the noie cancellation. ig. 9 how the meaured power pectra of the reidual ignal. Curve 1i the pectrum of the background noie. Curve i the pectrum of the primary loudpeaker. Curve 3 and 4 are the pectra of the reidual ignal after cancellation with repectively M ; 1= = 56 and 14. The average attenuation over a frequency range from 1 to 15 khz i 9 db for M =57and 36 db for M =15. VI. Concluion ig. 9. ower pectrum of the background noie (curve 1), the reidual ignal before cancellation (curve )and after cancellation (curve 3, B=51 and curve 4, B=48). mented in real-time uing block proceing technique in a unitary tranformation domain, uch a ourier domain. everal implementation of the ltered-x block frequency domain adaptive lter are dicued and their complexitie are compared. imulation and experiment have hown that type implementation i the mot ecient when the length of the econdary ource path etimate i longer than half the length of the adaptive lter. When the adaptive lter length i much longer than the econdary path etimate, type 4 ha hown to be the mot ecient. Reference [1] B. Widrow and.d. tearn, Adaptive ignal roceing, rentice-hall, Inc., Englewood Cli, N.J., [] A.V. Oppenheim and R.W. chafer, Digital ignal roceing, rentice-hall, Englewood Cli, N.J., [3].C.W. ommen, Adaptive iltering Method, h.d. thei Eindhoven Univerity of Technology The Netherland, June 199, IBN [4] R. Aart, H. He,. ommen, hantom ource applied to tereo-bae widening uing long frequency domain adaptive lter roceeding of the roric/ieee workhop 1996, pp [5] A.W.M. Mathijen, Generation of phantom ound ource with Block requency Domain Adaptive iltering, Mater thei of Eindhoven Univerity of Technology, eptember Generation of a virtual ound image uing adaptive lter that employ the ltered-x M algorithm ha been dicued. uch complex lter are bet imple-
Design 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 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 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 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 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 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 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 informationAdvanced 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 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 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 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 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: 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 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 information5.5 Sampling. The Connection Between: Continuous Time & Discrete Time
5.5 Sampling he Connection Between: Continuou ime & Dicrete ime Warning: I don t really like how the book cover thi! It i not that it i wrong it jut ail to make the correct connection between the mathematic
More informationLTV System Modelling
Helinki Univerit of Technolog S-72.333 Potgraduate Coure in Radiocommunication Fall 2000 LTV Stem Modelling Heikki Lorentz Sonera Entrum O heikki.lorentz@onera.fi Januar 23 rd 200 Content. Introduction
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 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 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 informationPractice Problems - Week #7 Laplace - Step Functions, DE Solutions Solutions
For Quetion -6, rewrite the piecewie function uing tep function, ketch their graph, and find F () = Lf(t). 0 0 < t < 2. f(t) = (t 2 4) 2 < t In tep-function form, f(t) = u 2 (t 2 4) The graph i the olid
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 informationNONLINEAR CONTROLLER DESIGN FOR A SHELL AND TUBE HEAT EXCHANGER AN EXPERIMENTATION APPROACH
International Journal of Electrical, Electronic and Data Communication, ISSN: 232-284 Volume-3, Iue-8, Aug.-25 NONLINEAR CONTROLLER DESIGN FOR A SHELL AND TUBE HEAT EXCHANGER AN EXPERIMENTATION APPROACH
More informationSINGLE CARRIER BLOCK TRANSMISSION WITHOUT GUARD INTERVAL
SINGLE CARRIER BLOCK TRANSMISSION WITHOUT GUARD INTERVAL Kazunori Hayahi Hideaki Sakai Graduate School of Informatic, Kyoto Univerity Kyoto, JAPAN ABSTRACT Thi paper propoe a imple detection cheme for
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 informationTrellis Shaping Techniques for Satellite Telecommunication Systems
> 1 Trelli Technique for Satellite Telecommunication Sytem M. Álvarez-Díaz 1, M. Neri, C. Moquera 1 and G. E. Corazza 1 Dept. Signal Theory and Communication (TSC), Univerity of Vigo Campu Univeritario
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 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 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 informationSAMPLING. Sampling is the acquisition of a continuous signal at discrete time intervals and is a fundamental concept in real-time signal processing.
SAMPLING Sampling i the acquiition of a continuou ignal at dicrete time interval and i a fundamental concept in real-time ignal proceing. he actual ampling operation can alo be defined by the figure belo
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 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 informationWhite Rose Research Online URL for this paper: Version: Accepted Version
Thi i a repoitory copy of Identification of nonlinear ytem with non-peritent excitation uing an iterative forward orthogonal leat quare regreion algorithm. White Roe Reearch Online URL for thi paper: http://eprint.whiteroe.ac.uk/107314/
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 informationLoss Less Image firmness comparision by DPCM and DPCM with LMS Algorithm
Lo Le Image firmne compariion by DPCM and DPCM with LMS Algorithm Pramod Kumar Rajput 1 and Brijendra Mihra 2 1,2 Department of Electronic & Communication, Nagaji Intitute of Technology and Management
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 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 information84 ZHANG Jing-Shang Vol. 39 of which would emit 5 He rather than 3 He. 5 He i untable and eparated into n + pontaneouly, which can alo be treated a if
Commun. Theor. Phy. (Beijing, China) 39 (003) pp. 83{88 c International Academic Publiher Vol. 39, No. 1, January 15, 003 Theoretical Analyi of Neutron Double-Dierential Cro Section of n+ 11 B at 14. MeV
More informationEE Control Systems LECTURE 6
Copyright FL Lewi 999 All right reerved EE - Control Sytem LECTURE 6 Updated: Sunday, February, 999 BLOCK DIAGRAM AND MASON'S FORMULA A linear time-invariant (LTI) ytem can be repreented in many way, including:
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 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 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 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 informationStudy on the effect of vent on the electroacoustic absorber
INTER-NOISE 016 Study on the effect of vent on the electroacoutic aborber Youngeun Cho 1 ; Semyung Wang 1, Gwangju Intitute of Science and Technology, Republic of Korea ABSTRACT Electroacoutic aborber
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 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 informationRandom vs. Deterministic Deployment of Sensors in the Presence of Failures and Placement Errors
Random v. Determinitic Deployment of Senor in the Preence of Failure and Placement Error Paul Baliter Univerity of Memphi pbalitr@memphi.edu Santoh Kumar Univerity of Memphi antoh.kumar@memphi.edu Abtract
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 informationSuggestions - Problem Set (a) Show the discriminant condition (1) takes the form. ln ln, # # R R
Suggetion - Problem Set 3 4.2 (a) Show the dicriminant condition (1) take the form x D Ð.. Ñ. D.. D. ln ln, a deired. We then replace the quantitie. 3ß D3 by their etimate to get the proper form for thi
More informationECE 3510 Root Locus Design Examples. PI To eliminate steady-state error (for constant inputs) & perfect rejection of constant disturbances
ECE 350 Root Locu Deign Example Recall the imple crude ervo from lab G( ) 0 6.64 53.78 σ = = 3 23.473 PI To eliminate teady-tate error (for contant input) & perfect reection of contant diturbance Note:
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 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 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 informationMoment of Inertia of an Equilateral Triangle with Pivot at one Vertex
oment of nertia of an Equilateral Triangle with Pivot at one Vertex There are two wa (at leat) to derive the expreion f an equilateral triangle that i rotated about one vertex, and ll how ou both here.
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 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 information[Saxena, 2(9): September, 2013] ISSN: Impact Factor: INTERNATIONAL JOURNAL OF ENGINEERING SCIENCES & RESEARCH TECHNOLOGY
[Saena, (9): September, 0] ISSN: 77-9655 Impact Factor:.85 IJESRT INTERNATIONAL JOURNAL OF ENGINEERING SCIENCES & RESEARCH TECHNOLOGY Contant Stre Accelerated Life Teting Uing Rayleigh Geometric Proce
More informationA Simple Approach to Synthesizing Naïve Quantized Control for Reference Tracking
A Simple Approach to Syntheizing Naïve Quantized Control for Reference Tracking SHIANG-HUA YU Department of Electrical Engineering National Sun Yat-Sen Univerity 70 Lien-Hai Road, Kaohiung 804 TAIAN Abtract:
More information1 Basic Equations of the PLLs
1 Baic Equation of the PLL 1.1 INTRODUCTION Phae lock loop (PLL) belong to a larger et of regulation ytem. A an independent reearch and deign field it tarted in the 1950 [1] and gained major practical
More informationModule 4: Time Response of discrete time systems Lecture Note 1
Digital Control Module 4 Lecture Module 4: ime Repone of dicrete time ytem Lecture Note ime Repone of dicrete time ytem Abolute tability i a baic requirement of all control ytem. Apart from that, good
More informationMassachusetts Institute of Technology Dynamics and Control II
I E Maachuett Intitute of Technology Department of Mechanical Engineering 2.004 Dynamic and Control II Laboratory Seion 5: Elimination of Steady-State Error Uing Integral Control Action 1 Laboratory Objective:
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 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 informationCONTROL SYSTEMS, ROBOTICS AND AUTOMATION Vol. VIII Decoupling Control - M. Fikar
DECOUPLING CONTROL M. Fikar Department of Proce Control, Faculty of Chemical and Food Technology, Slovak Univerity of Technology in Bratilava, Radlinkého 9, SK-812 37 Bratilava, Slovakia Keyword: Decoupling:
More informationTrajectory Planning and Feedforward Design for High Performance Motion Systems
Trajectory Planning and Feedforward Deign for High Performance Motion Sytem Paul Lambrecht, Matthij Boerlage, Maarten Steinbuch Faculty of Mechanical Engineering, Control Sytem Technology Group Eindhoven
More informationCorrection for Simple System Example and Notes on Laplace Transforms / Deviation Variables ECHE 550 Fall 2002
Correction for Simple Sytem Example and Note on Laplace Tranform / Deviation Variable ECHE 55 Fall 22 Conider a tank draining from an initial height of h o at time t =. With no flow into the tank (F in
More informationThe Laplace Transform , Haynes Miller and Jeremy Orloff
The Laplace Tranform 8.3, Hayne Miller and Jeremy Orloff Laplace tranform baic: introduction An operator take a function a input and output another function. A tranform doe the ame thing with the added
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 informationNeural Network Linearization of Pressure Force Sensor Transfer Characteristic
Acta Polytechnica Hungarica Vol., No., 006 Neural Network Linearization of Preure Force Senor Tranfer Characteritic Jozef Vojtko, Irena Kováčová, Ladilav Madaráz, Dobrolav Kováč Faculty of Electrical Engineering
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 informationControl of Delayed Integrating Processes Using Two Feedback Controllers R MS Approach
Proceeding of the 7th WSEAS International Conference on SYSTEM SCIENCE and SIMULATION in ENGINEERING (ICOSSSE '8) Control of Delayed Integrating Procee Uing Two Feedback Controller R MS Approach LIBOR
More informationFOURIER-BASED METHODS FOR THE SPECTRAL ANALYSIS OF MUSICAL SOUNDS. Sylvain Marchand
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
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 informationControl Systems Analysis and Design by the Root-Locus Method
6 Control Sytem Analyi and Deign by the Root-Locu Method 6 1 INTRODUCTION The baic characteritic of the tranient repone of a cloed-loop ytem i cloely related to the location of the cloed-loop pole. If
More informationAn Inequality for Nonnegative Matrices and the Inverse Eigenvalue Problem
An Inequality for Nonnegative Matrice and the Invere Eigenvalue Problem Robert Ream Program in Mathematical Science The Univerity of Texa at Dalla Box 83688, Richardon, Texa 7583-688 Abtract We preent
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 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 informationEarth Potential Rise (EPR) Computation for a Fault on Transmission Mains Pole
ACN: 32586675 ABN: 8632586675 NATIONAL ELECTRICAL ENGINEERING CONULTANCY Earth otential Rie (ER) Computation for a Fault on Tranmiion Main ole repared by: M. Naereddine ACN: 32586675 ABN: 8632586675 Abtract
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 informationFast Convolutional Sparse Coding (FCSC)
Fat Convolutional Spare Coding (FCSC) Bailey ong Department of Computer Science Univerity of California, Irvine bhkong@ic.uci.edu Charle C. Fowlke Department of Computer Science Univerity of California,
More informationPIPELINED DIVISION OF SIGNED NUMBERS WITH THE USE OF RESIDUE ARITHMETIC FOR SMALL NUMBER RANGE WITH THE PROGRAMMABLE GATE ARRAY
POZNAN UNIVE RSITY OF TE CHNOLOGY ACADE MIC JOURNALS No 76 Electrical Engineering 03 Robert SMYK* Zenon ULMAN* Maciej CZYŻAK* PIPELINED DIVISION OF SIGNED NUMBERS WITH THE USE OF RESIDUE ARITHMETIC FOR
More informationWhat lies between Δx E, which represents the steam valve, and ΔP M, which is the mechanical power into the synchronous machine?
A 2.0 Introduction In the lat et of note, we developed a model of the peed governing mechanim, which i given below: xˆ K ( Pˆ ˆ) E () In thee note, we want to extend thi model o that it relate the actual
More informationAn Image-encoded Mach-Zehnder Joint Transform Correlator for Polychromatic Pattern Recognition with Multi-level Quantized Reference Functions
Proceeding of the International MultiConference of Engineer and Computer Scientit 008 Vol I IMECS 008, 19-1 March, 008, Hong Kong An Image-encoded Mach-Zehnder Joint Tranform Correlator for Polychromatic
More information11.5 MAP Estimator MAP avoids this Computational Problem!
.5 MAP timator ecall that the hit-or-mi cot function gave the MAP etimator it maimize the a oteriori PDF Q: Given that the MMS etimator i the mot natural one why would we conider the MAP etimator? A: If
More informationInference for Two Stage Cluster Sampling: Equal SSU per PSU. Projections of SSU Random Variables on Each SSU selection.
Inference for Two Stage Cluter Sampling: Equal SSU per PSU Projection of SSU andom Variable on Eac SSU election By Ed Stanek Introduction We review etimating equation for PSU mean in a two tage cluter
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 informationThe Hassenpflug Matrix Tensor Notation
The Haenpflug Matrix Tenor Notation D.N.J. El Dept of Mech Mechatron Eng Univ of Stellenboch, South Africa e-mail: dnjel@un.ac.za 2009/09/01 Abtract Thi i a ample document to illutrate the typeetting of
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 informationSolving Differential Equations by the Laplace Transform and by Numerical Methods
36CH_PHCalter_TechMath_95099 3//007 :8 PM Page Solving Differential Equation by the Laplace Tranform and by Numerical Method OBJECTIVES When you have completed thi chapter, you hould be able to: Find the
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 information1. Intensity of Periodic Sound Waves 2. The Doppler Effect
1. Intenity o Periodic Sound Wae. The Doppler Eect 1-4-018 1 Objectie: The tudent will be able to Deine the intenity o the ound wae. Deine the Doppler Eect. Undertand ome application on ound 1-4-018 3.3
More informationDimensional Analysis A Tool for Guiding Mathematical Calculations
Dimenional Analyi A Tool for Guiding Mathematical Calculation Dougla A. Kerr Iue 1 February 6, 2010 ABSTRACT AND INTRODUCTION In converting quantitie from one unit to another, we may know the applicable
More informationThermal Σ- Modulator: Anemometer Performance Analysis
Intrumentation and Meaurement Technology Conference IMTC 007 Waraw, Poland, May 1-3, 007 Thermal Σ- Modulator: Anemometer Performance Analyi Will R. M. Almeida 1, Georgina M. Freita 1, Lígia S. Palma 3,
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 informationImproved Interference Cancellation Scheme for Two-User Detection of Alamouti Code
Improved Interference Cancellation Scheme for Two-Uer Detection of Alamouti Code Manav R hatnagar and Are Hjørungne Abtract In thi correpondence, we propoe an improved interference cancellation method
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 informationSOURCE CHARACTERIZATION OF INLAND CRUSTAL EARTHQUAKES FOR NEAR-SOURCE GROUND MOTIONS
SOURCE CHARACTERIZATION OF INLAND CRUSTAL EARTHQUAKES FOR NEAR-SOURCE GROUND MOTIONS Hiroe Miyake *, Tomotaka Iwata **, and Kojiro Irikura *** ABSTRACT To obtain the characterized ource model for prediction
More informationAngle Beam Ultrasonic Testing Models and Their Application to Identification and Sizing of Cracks
Angle Beam Ultraonic Teting Model and Their Application to Identification and iing of Crack NDE00 predict. aure. improve. National eminar of INT Chennai,. 7.. 00 www.nde00.org ung-jin ong*, Hee Jun Jung*,
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 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 informationSERIES COMPENSATION: VOLTAGE COMPENSATION USING DVR (Lectures 41-48)
Chapter 5 SERIES COMPENSATION: VOLTAGE COMPENSATION USING DVR (Lecture 41-48) 5.1 Introduction Power ytem hould enure good quality of electric power upply, which mean voltage and current waveform hould
More information