UNBALANCED MACHINE FAULT DETECTION USING INSTANTANEOUS FREQUENCY
|
|
- Claribel O’Brien’
- 5 years ago
- Views:
Transcription
1 158 UNBALANCED MACHINE FAULT DETECTION USING INSTANTANEOUS FREQUENCY Dhay Arifiato ad Adi Rahmadiasyah Departmet of Egieerig Physics, Faculty of Idustrial Techology Sepuluh Nopember Istitute of Techology Campus ITS Sukolilo Surabaya dhay@ep.its.ac.id Abstract This paper presets a ovel techique to aalyze acoustical sigal emitted from rotatig machiery by usig istataeous frequecy (IF). IF is defied as the derivative of phase with respect to time from the aalytic sigal to track a timevaryig parameter. This techique essetially is derived from the spectrum of short-time Fourier trasform (STFT) ad the remapped that spectrum to derivatio of widow used i STFT aalysis. As a result, the time ad frequecy ca be located simultaeously with better represetatio tha spectrogram. I this work, we focused o the ubalaced rotatig machie due to ueve mass at the ceter of rotatio. The soud patter of the vibratio machie simulator was recorded i Pusdiklat Migas Cepu with kow ubalaced mass to simulate damage coditios. The more mass was added to the rotor, the louder soud of ubalaced would be ad this also implied the higher value of the amplitudes of the fudametal frequecy. For ubalaced coditio with 7.8gr of mass, the domiat fudametal frequecy at 600 rpm, was at Hz (15.41 db SPL). Keywords : rotatig machiery, istataeous frequecy, fudametal frequecy, soud patter 1. INTRODUCTION Rotatig machiery is commoly foud i daily life from home appliaces (e.g., bleder) to large-scale idustry (e.g., electric turbie). Essetially, this type of machiery is drive by electrical iput ad trasformed ito mechaical eergy (rotatio). The rotatio ad its iteractio betwee movig parts durig cyclic motio cause vibratio ad emit audible soud. Due to mechaical imperfectio such as wear-time, frictio, the each rotatig machiery has a uique soud characteristic. This research focuses o aalyzig the emitted vibratio soud of its particular frequecy, i which, the damage of oe particular part of the rotatig machie will emit a specific frequecy. Recetly, LPA method was used i previous research [3]. This techique is based o the assumptio that the slow chages of sigal spectrum i a short rage of widow (< 30 msec). However, the vibratio of rotatig machie is cosiderably lower tha huma laryx. Therefore, the probability of frequecy estimatio error is larger due to the differece from of the ma s soud spectrum ad machie i a short rage of widow. Istataeous frequecy, which is defied as the derivative of phase with respect of time, is a ovel tool to estimate time-varyig parameter with high accuracy. This work also exteds previous ivestigatio by usig a simulator of a real plat [1]. The results were obtaied by determiig the value of fudametal frequecy of ormal machie ad compared it with the other damage part. We used vibratio machie simulator i Pusdiklat Migas Cepu after kowig the kid of machie fault for validatio.. MODEL, ANALYSIS, DESIGN, AND IMPLEMENTATION Rotatig machiery is used to trasform the electrical eergy ito mechaical eergy i form of rotatio. I this research we focused o the damaged typically foud o a rotatig machiery called ubalace. This class of damage is caused by a ueve mass at rotor (the cetre of rotatio)..1 Short Time Fourier Trasform (STFT) STFT is the Fourier Trasform covoluted with a widow fuctio. Fourier Trasformatio from the sigal x(t) chaged ito discrete from is expressed as, L 0 X (e ) = x(m)ω( m)e m (1) m= 0 The eq. (1) has two differet iterpretatios : usig the is costat, X (e ) is same as discrete Fourier Trasformatio which multiplied by the widow fuctio, expressed as,
2 08 Ubalaced Machie Fault Detectio usig Istataeous Frequecy- Dhay Arifiato 159 π 1 θ X (e ) = θ π X(e )* W(e )e d () π I this iterpretatio X ( e ) is Fourier Trasform from x(m)ω(-m), meas that, * < m <. It X ( e ) = X ( e ) X ( e ) (3) eq. (3) ca be expressed ito real ad imagiary part, X ( e ) = a ( ω) + jb ( ω). Usig certai value of ω, STFT ca be viewed as a bak of badpass filters with i is filter idex of ith-filter, i X (e ) = e i x i () (4) ad the impulse respose of sigle filter is expressed i h( ) = w( ) e with : 1. w( ) 0 for 0 L 1. Badwidth from h() is ω 0, 0, ω < ω where, 0 W ( e ) = 0, other ad ω 0 is the ceter frequecy of filter h() []. There are two kids of spectrogram, arrowbad ad widebad. Narrowbad is produced by wide widow, so it will produce the good resolutio at frequecy ad decreasig resolutio i time dimesio. O the other had, the widebad spectrogram, the usig of short widow will produce good resolutio i time dimesio while decreasig resolutio i frequecy dimesio.. Widowig Practically sigal is just moitored at certai period of time. Thus, the sigal spectrum must be approached usig some fiite data. If there is a sigal x 0 (t) with T 0 duratio, where T 0 T ad T = 1/Fs, so the ability to differ the frequecy will be restricted at F = 1/Fs. It some sample L is the result from Fs samplig, if the value of duratio x 0 (t) is loger tha T 0, so widow ω() with L duratio eed to be used, expressed as : ˆ = ω (5) x( ) x( ) ( ).3 Istataeous Frequecy The istataeous frequecy (IF) is a frequecy that is fuctio of time associated with the sigal which describes the frequecy value i istat time. The istataeous frequecy amplitude spectrum (IFAS) provides better the harmoic structure represetatio of speech sigal tha the STFT amplitude spectrum. As ca be see i Fig. 1, the IF techique ca elimiate Heiseberg s ucertaity which is the exact locatio of frequecy ca be determiate at its exact locatio i time simultaeously (represeted by spikes). IFAS is obtaied by rearragig the sigal x(t) covolved with widow fuctio w(t), τ X ( ω, t) = w' ( τ t) e x( τ ) dτ (6) t w (t) fuctio is the derivatio of widow w(t) with respect to time. Figure 1. Defiitio of istataeous frequecy amplitude spectrum [].4 Research Desig The experimetal uit of this research is vibratio machie simulator foud i the vibratio machie simulator room at Pusdiklat Migas Cepu. The machie specificatio is: Dimesio : legth = 89 cm, width = 0 cm, height = 31 cm Motor type is sigle phase AC electric Miimum rotatioal speed is 964 rpm ad maximum is 309 rpm at ormal coditio. We set the damaged as follows, 1. Ubalace coditio without mass. Ubalace coditio with.4 gr of mass 3. Ubalace coditio with 4.8 gr of mass The damaged level of ubalace ca be determied ca be determied by observig the value of amplitude at frequecy 1x rpm (machie rotatio) of vibratio spectra at each damaged coditio. Thus, the higher value of amplitude, the higher level of damaged due to the ubalace mass. The schematic of addig ubalaced mass was coducted, show i Figure. We set three differet variatios of rotatioal speed, 1160, 1840 ad 600 rpm. I order to kow the level damaged of each
3 160 4 th Iteratioal Coferece Iformatio & Commuicatio Techology ad System ubalace coditio coducted ubalace phase measuremet (for ubalace coditio without mass oly) ad the vibratio level (vibratio spectrum) usig VIBXPERT tool i Pusdiklat Migas Cepu. The level damaged at vertical misaligmet coditio was measured by usig ROTALIGN as baselie. Figure. Machiery scheme ad locatio of mass The soud pressure level measuremet was coducted usig the stadard of Oil Compaies Materials Associatio (OCMA) method for small source. Measuremet was take at 30 cm far from machie usig SLM (Rio Type NL-31) coected with computer. The steps are : Figure 3. Measuremet surface ad locatio of measuremet poit [5] I1 I I 3 0 D = + + (7) Do = commo dimesio, I 1 = legth (cm), I = width (cm), I 3 = high (cm) ad, r = Do (8) Usig equatio (7) ad (8), the value of Do was obtaied 55,147 cm, ad the, r = 110,194 cm. 1 The locatio of measuremet poits (show i Fig. 3) is determied by calculatig b ad h, where b is distace betwee poit ad source = 0.9r ad h is poit high = 0.4r. Result of the calculatio is b = 99.6 cm ad h =44.1 cm. Soud pressure level measured at each poit with average value, 0,1L1 0,1L 0,1L L = AV 10log (10) where L av is average SPL, L 1 L is SPL each poit (db), ad is amout of measuremet poit. The maximum frequecy of soud sigal was observed first to determie sample frequecy to fulfill the Nyquist criterio. The soud sigal machie rage is Hz ad weakeig amplitude occurs at frequecy over 5500 Hz. The samplig frequecy was set to 1105 Hz to avoid possible aliasig. The first step is STFT calculatio by dividig the sigal x(t) ito segmets usig shifted half of the legth of the widow alog samples. The spectrogram of the sigal shows the value of frequecy for each measuremet (ordiate) ad time axis i secod (abscissa). The shape ad width of the aalysis widow are importat features. It is ecessary to use the widow whose value of side lobe roll off is high to prevet loss due to frequecy leak. We chose Haig widow. Usig the fact that the istataeous phase agle of the iput sequece is the (uwrapped) agle of the aalytic sigal, the istataeous frequecy is the rate of chage of the istataeous phase agle with respect to time. First, the derivatio will depart from widowed Fourier trasform poit of view to estimate the istataeous frequecy. By meas of the otatio of istataeous frequecy amplitude spectrum, the the harmoic structure of sigal ca be revealed to determie the fudametal frequecy. First we calculated the STFT spectrum of the sigal x(t), the remapped the spectrum ito its respective IF bad. This bad was derivatio of widow fuctio with respect to time. I other words, the coefficiets i the movig widow method are calculated from the movig widow trasform (eq. ). The time-frequecy represetatio is costructed from the squared magitude of these coefficiets. The IF was computed for each STFT bi from partial derivatives of the phase spectrum as expressed by eq. (6). The upper ad lower boud of the IF value iside each frequecy bi were the frequecy cut-off of the badpass filters. If the IF value is outside the filter, the it will be excluded from further calculatio.
4 08 Ubalaced Machie Fault Detectio usig Istataeous Frequecy- Dhay Arifiato RESULTS For SPL (soud pressure level, i db) measuremets, the samples are the value of L eq for 10 secods for each poit ad damaged coditio for three variatio of rotatioal speed (1160, 1840, ad 600 rpm) preseted below : figure is at 1160 rpm, 1840 rpm for middle ad at the lower pael is at 600 rpm, respectively. As it ca be see i Fig. 5, for lower speed, the spectrum of the motor ubalaced was udetected. As the speed icreased, the up-ad-dow spectrum of ubalaced appeared obviously. Table 1. Average SPL of each coditio Average SPL (db) Speed (RPM) Without mass,4 gr 4,8 gr 7,8 Gr Nrml Miss Miss vert , ,4 69,4 68, 75, ,4 7,9 73,3 73,3 71,9 77, ,8 78,9 81,7 81,7 73,9 8,7 The SPL sample measuremet is used to fid out machie SPL distributio at which poit will be used for recordig the ext machie soud sigal. As the ubalaced mass grew, the soud itesity emitted was also louder. We preset oly oe figure of the spectrogram. Figure 4 shows a arrowbad spectrogram of the ubalaced with 4.8 gram additioal mass to the rotor at 600 rpm. As it ca be see, the spectrogram idicates the spectrum of the sigal below 500 Hz. However, it fails to locate precisely the characteristic frequecy of the motor due to smeared frequecy ad poor temporal resolutio. Figure 4. Spectrogram of ubalace with 4.8 gr of mass at 600 rpm I first coditio, the motor was set i ubalace without additioal mass for baselie coditio. Fudametal frequecy was observed at about Hz at 1160 rpm ad at about Hz at 600 rpm. I the ext step,.4, 4.8, ad 7.8 gram of ubalaced mass was added to the machie for each coditio. Figure 5, 6 ad 7 show the patter of frequecy locatio of the.4, 4.8, ad 7.8 gram of ubalaced mass with three differet rotatio speeds. The rotatio speed at upper pael of each Figure 5. Soud patter of ubalace with.4 gram of mass at 1160rpm(upper pael), 1840 rpm (middle) ad 600 rpm.
5 16 4 th Iteratioal Coferece Iformatio & Commuicatio Techology ad System It was apparet that the domiat frequecy was below 100 Hz. The spectrum was also oscillatig with respect to time. I the lower pael of Fig. 5, with the same amout of ubalaced mass ad the rotatio speed was icreased to the highest allowable speed of the motor, the amplitude of lower spectrum became higher tha the middle pael. It is also obvious that the side-spectrum (about 450 Hz) appears. This implies that the more eergy was required to drive the motor with additioal ubalaced mass as the speed icreased. Fudametal frequecy was observed at about 4.7 Hz ad 49.9 Hz i 600 rpm. Whe more ubalaced mass was added to the rotor, the damaged spectrum appeared earlier i the 1840 rpm (see middle pael of Fig. 6). Its amplitude was also comparable to the case of.4 gram ubalaced mass with 600 rpm (lower pael of Fig. 5). However, the siusoid-like oscillatio of the spectrum became deser which idicated higher oscillatio frequecy. It is also easily observed the side spectrum at about 300 Hz, which is ot easily observed with spectrogram. The other sidefrequecy is appeared at about 450 Hz while the mai ad the first side frequecy at 300 Hz have higher amplitude (i.e., itesity). Figure 6. Soud patter of ubalace with 4.8 gram of mass at 1160rpm(upper pael), 1840 rpm (middle) ad 600 rpm. Figure 7. Soud patter of ubalace with 7.8 gram of mass at 1160rpm(upper pael), 1840 rpm (middle) ad 600 rpm.
6 08 Ubalaced Machie Fault Detectio usig Istataeous Frequecy- Dhay Arifiato 163 Fudametal frequecy appeared at about.47 Hz at 1840 rpm ad Hz at 600 rpm for 4.8 gram ubalaced mass case. I the third case, 7.8 gram ubalaced mass was added ad the results are show i Fig. 7. Both the domiat spectrum ad side-spectrum were observed at 1840 rpm speed. However, the domiat spectrum was shifted upward to higher frequecy. I the lower pael of Fig. 7, the domiat spectrum at 600 rpm shifted eve higher tha that of at 1840 rpm. Two side-spectrums appear clearly with siusoid patter at about 50 Hz ad 400 Hz. The frequecy of oscillatio patter at 400 Hz is the same as the domiat spectrum. This suggests that the additioal mass may have iduced a ew vibratig place at the rotor. Compared to the spectrogram, the gradual additio of the ubalaced mass, the patters were also observable which was blurred dark lie patter i rage Hz of frequecy, as if it was oly oe domiat spectrum. The oly obvious differece was, as the ubalaced mass icreased the domiat spectrum appeared darker. Usig IF-based timefrequecy represetatio, the cyclic movemet of the rotatig part of the machie ca be viewed clearly ad the exact umeric value of frequecy ca be obtaied precisely compared to that of the spectrogram. This techique shows a prospective practical usage i the idustry for maiteace ad early failure detectio. [] Arifiato, Dhay. Fudametal Frequecy Estimatio ad Voiced / Uvoiced Determiatio Usig Istataeous Frequecy. Doctoral Dissertatio, Departemet of Iformatio Processig Tokyo Istitute of Techology, Tokyo, 005. [3] Hasa, Fuad. 00. Study of Patter Recogitio i Rotatig Machie Fault Detectio usig Liear Predictio (i Idoesia), udergraduate thesis, Dept. Egieerig Physics ITS, Surabaya, 00. [4] Uited Natios Eviromet Programme, Eergy Efficiecy Guide for Idustry i Asia- Electrical Motor (i Idoesia) USA, 006. [5] Yag, S.J. ad A. J. Elliso Machiery Noise Measuremet. New York, USA : Claredo Press Oxford. 4. CONCLUDING REMARKS I this study, we applied istataeous frequecy based time-frequecy represetatio to show the spectrum cotet of gradual additioal ubalaced mass to a rotor of a motor. We also compared the results with the spectrogram. The ubalaced mass was added to the rotor.4, 4.8, ad 7.8 gram, as well as the rotatio speed was varied to 1600, 1840, ad 600 rpm, respectively. Usig istataeous frequecy techique, the domiat spectrum was clearly observed. As the rotatio speed was icreased ad more ubalaced mass was added, the spectrum amplitude was gettig higher. It was also clearly show the side-spectrum also appeared. This suggests that the ubalaced mass may have iduced additioal vibratio elsewhere i the machie simultaeously with the mai vibratio at the rotor. REFERENCES [1] Ada, Ahmad Rifqi Applicatio of Istataeous Frequecy Techique for Fault Detectio of Electrical Motor. (i Idoesia), udergraduate thesis, Dept. Egieerig Physics ITS, Surabaya, 006.
Direction of Arrival Estimation Method in Underdetermined Condition Zhang Youzhi a, Li Weibo b, Wang Hanli c
4th Iteratioal Coferece o Advaced Materials ad Iformatio Techology Processig (AMITP 06) Directio of Arrival Estimatio Method i Uderdetermied Coditio Zhag Youzhi a, Li eibo b, ag Hali c Naval Aeroautical
More informationSpectral Analysis. This week in lab. Next classes: 3/26 and 3/28. Your next experiment Homework is to prepare
Spectral Aalysis This week i lab Your ext experimet Homework is to prepare Next classes: 3/26 ad 3/28 Aero Testig, Fracture Toughess Testig Read the Experimets 5 ad 7 sectios of the course maual Spectral
More information577. Estimation of surface roughness using high frequency vibrations
577. Estimatio of surface roughess usig high frequecy vibratios V. Augutis, M. Sauoris, Kauas Uiversity of Techology Electroics ad Measuremets Systems Departmet Studetu str. 5-443, LT-5368 Kauas, Lithuaia
More informationMechatronics. Time Response & Frequency Response 2 nd -Order Dynamic System 2-Pole, Low-Pass, Active Filter
Time Respose & Frequecy Respose d -Order Dyamic System -Pole, Low-Pass, Active Filter R 4 R 7 C 5 e i R 1 C R 3 - + R 6 - + e out Assigmet: Perform a Complete Dyamic System Ivestigatio of the Two-Pole,
More informationEE / EEE SAMPLE STUDY MATERIAL. GATE, IES & PSUs Signal System. Electrical Engineering. Postal Correspondence Course
Sigal-EE Postal Correspodece Course 1 SAMPLE STUDY MATERIAL Electrical Egieerig EE / EEE Postal Correspodece Course GATE, IES & PSUs Sigal System Sigal-EE Postal Correspodece Course CONTENTS 1. SIGNAL
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science. BACKGROUND EXAM September 30, 2004.
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Departmet of Electrical Egieerig ad Computer Sciece 6.34 Discrete Time Sigal Processig Fall 24 BACKGROUND EXAM September 3, 24. Full Name: Note: This exam is closed
More informationSignal Processing in Mechatronics. Lecture 3, Convolution, Fourier Series and Fourier Transform
Sigal Processig i Mechatroics Summer semester, 1 Lecture 3, Covolutio, Fourier Series ad Fourier rasform Dr. Zhu K.P. AIS, UM 1 1. Covolutio Covolutio Descriptio of LI Systems he mai premise is that the
More informationSignals & Systems Chapter3
Sigals & Systems Chapter3 1.2 Discrete-Time (D-T) Sigals Electroic systems do most of the processig of a sigal usig a computer. A computer ca t directly process a C-T sigal but istead eeds a stream of
More informationCUMULATIVE DAMAGE ESTIMATION USING WAVELET TRANSFORM OF STRUCTURAL RESPONSE
CUMULATIVE DAMAGE ESTIMATION USING WAVELET TRANSFORM OF STRUCTURAL RESPONSE Ryutaro SEGAWA 1, Shizuo YAMAMOTO, Akira SONE 3 Ad Arata MASUDA 4 SUMMARY Durig a strog earthquake, the respose of a structure
More informationx[0] x[1] x[2] Figure 2.1 Graphical representation of a discrete-time signal.
x[ ] x[ ] x[] x[] x[] x[] 9 8 7 6 5 4 3 3 4 5 6 7 8 9 Figure. Graphical represetatio of a discrete-time sigal. From Discrete-Time Sigal Processig, e by Oppeheim, Schafer, ad Buck 999- Pretice Hall, Ic.
More informationVibratory Motion. Prof. Zheng-yi Feng NCHU SWC. National CHung Hsing University, Department of Soil and Water Conservation
Vibratory Motio Prof. Zheg-yi Feg NCHU SWC 1 Types of vibratory motio Periodic motio Noperiodic motio See Fig. A1, p.58 Harmoic motio Periodic motio Trasiet motio impact Trasiet motio earthquake A powerful
More informationIntroduction to Signals and Systems, Part V: Lecture Summary
EEL33: Discrete-Time Sigals ad Systems Itroductio to Sigals ad Systems, Part V: Lecture Summary Itroductio to Sigals ad Systems, Part V: Lecture Summary So far we have oly looked at examples of o-recursive
More informationPractical Spectral Anaysis (continue) (from Boaz Porat s book) Frequency Measurement
Practical Spectral Aaysis (cotiue) (from Boaz Porat s book) Frequecy Measuremet Oe of the most importat applicatios of the DFT is the measuremet of frequecies of periodic sigals (eg., siusoidal sigals),
More informationOrthogonal Gaussian Filters for Signal Processing
Orthogoal Gaussia Filters for Sigal Processig Mark Mackezie ad Kiet Tieu Mechaical Egieerig Uiversity of Wollogog.S.W. Australia Abstract A Gaussia filter usig the Hermite orthoormal series of fuctios
More informationOBJECTIVES. Chapter 1 INTRODUCTION TO INSTRUMENTATION FUNCTION AND ADVANTAGES INTRODUCTION. At the end of this chapter, students should be able to:
OBJECTIVES Chapter 1 INTRODUCTION TO INSTRUMENTATION At the ed of this chapter, studets should be able to: 1. Explai the static ad dyamic characteristics of a istrumet. 2. Calculate ad aalyze the measuremet
More informationThe z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j
The -Trasform 7. Itroductio Geeralie the complex siusoidal represetatio offered by DTFT to a represetatio of complex expoetial sigals. Obtai more geeral characteristics for discrete-time LTI systems. 7.
More informationDiscrete-Time Signals and Systems. Signals and Systems. Digital Signals. Discrete-Time Signals. Operations on Sequences: Basic Operations
-6.3 Digital Sigal Processig ad Filterig..8 Discrete-ime Sigals ad Systems ime-domai Represetatios of Discrete-ime Sigals ad Systems ime-domai represetatio of a discrete-time sigal as a sequece of umbers
More informationOlli Simula T / Chapter 1 3. Olli Simula T / Chapter 1 5
Sigals ad Systems Sigals ad Systems Sigals are variables that carry iformatio Systemstake sigals as iputs ad produce sigals as outputs The course deals with the passage of sigals through systems T-6.4
More informationEXPERIMENT OF SIMPLE VIBRATION
EXPERIMENT OF SIMPLE VIBRATION. PURPOSE The purpose of the experimet is to show free vibratio ad damped vibratio o a system havig oe degree of freedom ad to ivestigate the relatioship betwee the basic
More information[ ] sin ( ) ( ) = 2 2 ( ) ( ) ( ) ˆ Mechanical Spectroscopy II
Solid State Pheomea Vol. 89 (003) pp 343-348 (003) Tras Tech Publicatios, Switzerlad doi:0.408/www.scietific.et/ssp.89.343 A New Impulse Mechaical Spectrometer to Study the Dyamic Mechaical Properties
More informationFIR Filters. Lecture #7 Chapter 5. BME 310 Biomedical Computing - J.Schesser
FIR Filters Lecture #7 Chapter 5 8 What Is this Course All About? To Gai a Appreciatio of the Various Types of Sigals ad Systems To Aalyze The Various Types of Systems To Lear the Skills ad Tools eeded
More informationA. Basics of Discrete Fourier Transform
A. Basics of Discrete Fourier Trasform A.1. Defiitio of Discrete Fourier Trasform (8.5) A.2. Properties of Discrete Fourier Trasform (8.6) A.3. Spectral Aalysis of Cotiuous-Time Sigals Usig Discrete Fourier
More informationis also known as the general term of the sequence
Lesso : Sequeces ad Series Outlie Objectives: I ca determie whether a sequece has a patter. I ca determie whether a sequece ca be geeralized to fid a formula for the geeral term i the sequece. I ca determie
More informationSignal Processing in Mechatronics
Sigal Processig i Mechatroics Zhu K.P. AIS, UM. Lecture, Brief itroductio to Sigals ad Systems, Review of Liear Algebra ad Sigal Processig Related Mathematics . Brief Itroductio to Sigals What is sigal
More informationModule 11: Applications : Linear prediction, Speech Analysis and Speech Enhancement Prof. Eliathamby Ambikairajah Dr. Tharmarajah Thiruvaran School
Module : Applicatios : Liear predictio, Speech Aalysis ad Speech Ehacemet Prof. Eliathamby Ambiairajah Dr. Tharmarajah Thiruvara School of Electrical Egieerig & Telecommuicatios The Uiversity of New South
More informationUNDERWATER OBJECT CLASSIFICATION BY MEANS OF AN ACOUSTIC METHOD EUGENIUSZ KOZACZKA
UNDERWATER OBJECT CLASSIFICATION BY MEANS OF AN ACOUSTIC METHOD EUGENIUSZ KOZACZKA Naval Uiversity of Gdyia 81-13 Gdyia, Śmidowicza 69, Polad Gdańsk Uiversity of Techology 8-95 Gdańsk, Narutowicza 11/1,
More informationSolutions. Number of Problems: 4. None. Use only the prepared sheets for your solutions. Additional paper is available from the supervisors.
Quiz November 4th, 23 Sigals & Systems (5-575-) P. Reist & Prof. R. D Adrea Solutios Exam Duratio: 4 miutes Number of Problems: 4 Permitted aids: Noe. Use oly the prepared sheets for your solutios. Additioal
More informationANALYSIS OF EXPERIMENTAL ERRORS
ANALYSIS OF EXPERIMENTAL ERRORS All physical measuremets ecoutered i the verificatio of physics theories ad cocepts are subject to ucertaities that deped o the measurig istrumets used ad the coditios uder
More informationStatistical Fundamentals and Control Charts
Statistical Fudametals ad Cotrol Charts 1. Statistical Process Cotrol Basics Chace causes of variatio uavoidable causes of variatios Assigable causes of variatio large variatios related to machies, materials,
More informationFIR Filter Design: Part II
EEL335: Discrete-Time Sigals ad Systems. Itroductio I this set of otes, we cosider how we might go about desigig FIR filters with arbitrary frequecy resposes, through compositio of multiple sigle-peak
More informationFrequency Response of FIR Filters
EEL335: Discrete-Time Sigals ad Systems. Itroductio I this set of otes, we itroduce the idea of the frequecy respose of LTI systems, ad focus specifically o the frequecy respose of FIR filters.. Steady-state
More informationMeasurement uncertainty of the sound absorption
Measuremet ucertaity of the soud absorptio coefficiet Aa Izewska Buildig Research Istitute, Filtrowa Str., 00-6 Warsaw, Polad a.izewska@itb.pl 6887 The stadard ISO/IEC 705:005 o the competece of testig
More informationProbability, Expectation Value and Uncertainty
Chapter 1 Probability, Expectatio Value ad Ucertaity We have see that the physically observable properties of a quatum system are represeted by Hermitea operators (also referred to as observables ) such
More informationSpring 2014, EE123 Digital Signal Processing
Aoucemets EE3 Digital Sigal Processig Last time: FF oday: Frequecy aalysis with DF Widowig Effect of zero-paddig Lecture 9 based o slides by J.M. Kah Spectral Aalysis with the DF Spectral Aalysis with
More informationAnalysis of the No-Load Characteristic of the Moving Coil Linear Compressor
Purdue Uiversity Purdue e-pubs Iteratioal Compressor Egieerig Coferece School of Mechaical Egieerig 008 Aalysis of the No-Load Characteristic of the Movig Coil Liear Compressor Yigbai Xie North Chia Electric
More informationResearch Article Health Monitoring for a Structure Using Its Nonstationary Vibration
Advaces i Acoustics ad Vibratio Volume 2, Article ID 69652, 5 pages doi:.55/2/69652 Research Article Health Moitorig for a Structure Usig Its Nostatioary Vibratio Yoshimutsu Hirata, Mikio Tohyama, Mitsuo
More informationExam. Notes: A single A4 sheet of paper (double sided; hand-written or computer typed)
Exam February 8th, 8 Sigals & Systems (5-575-) Prof. R. D Adrea Exam Exam Duratio: 5 Mi Number of Problems: 5 Number of Poits: 5 Permitted aids: Importat: Notes: A sigle A sheet of paper (double sided;
More informationLinear time invariant systems
Liear time ivariat systems Alejadro Ribeiro Dept. of Electrical ad Systems Egieerig Uiversity of Pesylvaia aribeiro@seas.upe.edu http://www.seas.upe.edu/users/~aribeiro/ February 25, 2016 Sigal ad Iformatio
More informationCHAPTER 10 INFINITE SEQUENCES AND SERIES
CHAPTER 10 INFINITE SEQUENCES AND SERIES 10.1 Sequeces 10.2 Ifiite Series 10.3 The Itegral Tests 10.4 Compariso Tests 10.5 The Ratio ad Root Tests 10.6 Alteratig Series: Absolute ad Coditioal Covergece
More informationA statistical method to determine sample size to estimate characteristic value of soil parameters
A statistical method to determie sample size to estimate characteristic value of soil parameters Y. Hojo, B. Setiawa 2 ad M. Suzuki 3 Abstract Sample size is a importat factor to be cosidered i determiig
More informationComplex Algorithms for Lattice Adaptive IIR Notch Filter
4th Iteratioal Coferece o Sigal Processig Systems (ICSPS ) IPCSIT vol. 58 () () IACSIT Press, Sigapore DOI:.7763/IPCSIT..V58. Complex Algorithms for Lattice Adaptive IIR Notch Filter Hog Liag +, Nig Jia
More information2C09 Design for seismic and climate changes
2C09 Desig for seismic ad climate chages Lecture 02: Dyamic respose of sigle-degree-of-freedom systems I Daiel Grecea, Politehica Uiversity of Timisoara 10/03/2014 Europea Erasmus Mudus Master Course Sustaiable
More informationNumerical Methods in Fourier Series Applications
Numerical Methods i Fourier Series Applicatios Recall that the basic relatios i usig the Trigoometric Fourier Series represetatio were give by f ( x) a o ( a x cos b x si ) () where the Fourier coefficiets
More informationDynamic System Response
Solutio of Liear, Costat-Coefficiet, Ordiary Differetial Equatios Classical Operator Method Laplace Trasform Method Laplace Trasform Properties 1 st -Order Dyamic System Time ad Frequecy Respose d -Order
More informationChapter 2 Feedback Control Theory Continued
Chapter Feedback Cotrol Theor Cotiued. Itroductio I the previous chapter, the respose characteristic of simple first ad secod order trasfer fuctios were studied. It was show that first order trasfer fuctio,
More informationWarped, Chirp Z-Transform: Radar Signal Processing
arped, Chirp Z-Trasform: Radar Sigal Processig by Garimella Ramamurthy Report o: IIIT/TR// Cetre for Commuicatios Iteratioal Istitute of Iformatio Techology Hyderabad - 5 3, IDIA Jauary ARPED, CHIRP Z
More informationProvläsningsexemplar / Preview TECHNICAL REPORT INTERNATIONAL SPECIAL COMMITTEE ON RADIO INTERFERENCE
TECHNICAL REPORT CISPR 16-4-3 2004 AMENDMENT 1 2006-10 INTERNATIONAL SPECIAL COMMITTEE ON RADIO INTERFERENCE Amedmet 1 Specificatio for radio disturbace ad immuity measurig apparatus ad methods Part 4-3:
More informationZ - Transform. It offers the techniques for digital filter design and frequency analysis of digital signals.
Z - Trasform The -trasform is a very importat tool i describig ad aalyig digital systems. It offers the techiques for digital filter desig ad frequecy aalysis of digital sigals. Defiitio of -trasform:
More informationUNIT-I. 2. A real valued sequence x(n) is anti symmetric if a) X(n)=x(-n) b) X(n)=-x(-n) c) A) and b) d) None Ans: b)
DIGITAL SIGNAL PROCESSING UNIT-I 1. The uit ramp sequece is Eergy sigal b) Power sigal c) Either Eergy or Power sigal d) Neither a Power sigal or a eergy sigal As: d) 2. A real valued sequece x() is ati
More informationChapter 7 z-transform
Chapter 7 -Trasform Itroductio Trasform Uilateral Trasform Properties Uilateral Trasform Iversio of Uilateral Trasform Determiig the Frequecy Respose from Poles ad Zeros Itroductio Role i Discrete-Time
More informationCh3 Discrete Time Fourier Transform
Ch3 Discrete Time Fourier Trasform 3. Show that the DTFT of [] is give by ( k). e k 3. Determie the DTFT of the two sided sigal y [ ],. 3.3 Determie the DTFT of the causal sequece x[ ] A cos( 0 ) [ ],
More informationChapter 10: Power Series
Chapter : Power Series 57 Chapter Overview: Power Series The reaso series are part of a Calculus course is that there are fuctios which caot be itegrated. All power series, though, ca be itegrated because
More information6.003 Homework #12 Solutions
6.003 Homework # Solutios Problems. Which are rue? For each of the D sigals x [] through x 4 [] below), determie whether the coditios listed i the followig table are satisfied, ad aswer for true or F for
More informationJitter Transfer Functions For The Reference Clock Jitter In A Serial Link: Theory And Applications
Jitter Trasfer Fuctios For The Referece Clock Jitter I A Serial Lik: Theory Ad Applicatios Mike Li, Wavecrest Ady Martwick, Itel Gerry Talbot, AMD Ja Wilstrup, Teradye Purposes Uderstad various jitter
More informationRun-length & Entropy Coding. Redundancy Removal. Sampling. Quantization. Perform inverse operations at the receiver EEE
Geeral e Image Coder Structure Motio Video (s 1,s 2,t) or (s 1,s 2 ) Natural Image Samplig A form of data compressio; usually lossless, but ca be lossy Redudacy Removal Lossless compressio: predictive
More informationChapter 7: The z-transform. Chih-Wei Liu
Chapter 7: The -Trasform Chih-Wei Liu Outlie Itroductio The -Trasform Properties of the Regio of Covergece Properties of the -Trasform Iversio of the -Trasform The Trasfer Fuctio Causality ad Stability
More informationStopping oscillations of a simple harmonic oscillator using an impulse force
It. J. Adv. Appl. Math. ad Mech. 5() (207) 6 (ISSN: 2347-2529) IJAAMM Joural homepage: www.ijaamm.com Iteratioal Joural of Advaces i Applied Mathematics ad Mechaics Stoppig oscillatios of a simple harmoic
More informationChapter 12 - Quality Cotrol Example: The process of llig 12 ouce cas of Dr. Pepper is beig moitored. The compay does ot wat to uderll the cas. Hece, a target llig rate of 12.1-12.5 ouces was established.
More informationResearch on real time compensation of thermal errors of CNC lathe based on linear regression theory Qiu Yongliang
d Iteratioal Coferece o Machiery, Materials Egieerig, Chemical Egieerig ad Biotechology (MMECEB 015) Research o real time compesatio of thermal errors of CNC lathe based o liear regressio theory Qiu Yogliag
More informationResearch Article Health Monitoring for a Structure Using Its Nonstationary Vibration
Hidawi Publishig Corporatio Advaces i Acoustics ad Vibratio Volume 2, Article ID 69652, 5 pages doi:.55/2/69652 Research Article Health Moitorig for a Structure Usig Its Nostatioary Vibratio Yoshimutsu
More informationMONITORING THE STABILITY OF SLOPES BY GPS
MONITORING THE STABILITY OF SLOPES BY GPS Prof. S. Sakurai Costructio Egieerig Research Istitute Foudatio, Japa Prof. N. Shimizu Dept. of Civil Egieerig, Yamaguchi Uiversity, Japa ABSTRACT The stability
More information6.003 Homework #12 Solutions
6.003 Homework # Solutios Problems. Which are rue? For each of the D sigals x [] through x 4 [] (below), determie whether the coditios listed i the followig table are satisfied, ad aswer for true or F
More informationA new iterative algorithm for reconstructing a signal from its dyadic wavelet transform modulus maxima
ol 46 No 6 SCIENCE IN CHINA (Series F) December 3 A ew iterative algorithm for recostructig a sigal from its dyadic wavelet trasform modulus maxima ZHANG Zhuosheg ( u ), LIU Guizhog ( q) & LIU Feg ( )
More informationCourse 4: Preparation for Calculus Unit 1: Families of Functions
Course 4: Preparatio for Calculus Uit 1: Families of Fuctios Review ad exted properties of basic fuctio families ad their uses i mathematical modelig Develop strategies for fidig rules of fuctios whose
More informationECE 8527: Introduction to Machine Learning and Pattern Recognition Midterm # 1. Vaishali Amin Fall, 2015
ECE 8527: Itroductio to Machie Learig ad Patter Recogitio Midterm # 1 Vaishali Ami Fall, 2015 tue39624@temple.edu Problem No. 1: Cosider a two-class discrete distributio problem: ω 1 :{[0,0], [2,0], [2,2],
More informationStatistics 511 Additional Materials
Cofidece Itervals o mu Statistics 511 Additioal Materials This topic officially moves us from probability to statistics. We begi to discuss makig ifereces about the populatio. Oe way to differetiate probability
More informationA PROCEDURE TO MODIFY THE FREQUENCY AND ENVELOPE CHARACTERISTICS OF EMPIRICAL GREEN'S FUNCTION. Lin LU 1 SUMMARY
A POCEDUE TO MODIFY THE FEQUENCY AND ENVELOPE CHAACTEISTICS OF EMPIICAL GEEN'S FUNCTION Li LU SUMMAY Semi-empirical method, which divides the fault plae of large earthquake ito mets ad uses small groud
More informationDiscrete Orthogonal Moment Features Using Chebyshev Polynomials
Discrete Orthogoal Momet Features Usig Chebyshev Polyomials R. Mukuda, 1 S.H.Og ad P.A. Lee 3 1 Faculty of Iformatio Sciece ad Techology, Multimedia Uiversity 75450 Malacca, Malaysia. Istitute of Mathematical
More informationMath 113 Exam 3 Practice
Math Exam Practice Exam 4 will cover.-., 0. ad 0.. Note that eve though. was tested i exam, questios from that sectios may also be o this exam. For practice problems o., refer to the last review. This
More informationChapter 2 Systems and Signals
Chapter 2 Systems ad Sigals 1 Itroductio Discrete-Time Sigals: Sequeces Discrete-Time Systems Properties of Liear Time-Ivariat Systems Liear Costat-Coefficiet Differece Equatios Frequecy-Domai Represetatio
More informationSpring 2014, EE123 Digital Signal Processing
Aoucemets EE3 Digital Sigal Processig Lecture 9 Lab part I ad II posted will post III today or tomorrow Lab-bash uesday -3pm Cory hree shorter Midterms: / i class / i class /3 (or BD) i class / or / (BD)
More informationDIGITAL MEASUREMENT OF POWER SYSTEM HARMONIC MAGNITUDE AND PHASE ANGLE
DIGIL MESUREMEN OF POWER SYSEM HRMONIC MGNIUDE ND PHSE NGLE R Micheletti (, R Pieri ( ( Departmet of Electrical Systems ad utomatio, Uiversity of Pisa, Via Diotisalvi, I-566 Pisa, Italy Phoe +39 5 565,
More informationNCSS Statistical Software. Tolerance Intervals
Chapter 585 Itroductio This procedure calculates oe-, ad two-, sided tolerace itervals based o either a distributio-free (oparametric) method or a method based o a ormality assumptio (parametric). A two-sided
More information3.2 Properties of Division 3.3 Zeros of Polynomials 3.4 Complex and Rational Zeros of Polynomials
Math 60 www.timetodare.com 3. Properties of Divisio 3.3 Zeros of Polyomials 3.4 Complex ad Ratioal Zeros of Polyomials I these sectios we will study polyomials algebraically. Most of our work will be cocered
More information6.3 Testing Series With Positive Terms
6.3. TESTING SERIES WITH POSITIVE TERMS 307 6.3 Testig Series With Positive Terms 6.3. Review of what is kow up to ow I theory, testig a series a i for covergece amouts to fidig the i= sequece of partial
More information7-1. Chapter 4. Part I. Sampling Distributions and Confidence Intervals
7-1 Chapter 4 Part I. Samplig Distributios ad Cofidece Itervals 1 7- Sectio 1. Samplig Distributio 7-3 Usig Statistics Statistical Iferece: Predict ad forecast values of populatio parameters... Test hypotheses
More informationMETHOD OF FUNDAMENTAL SOLUTIONS FOR HELMHOLTZ EIGENVALUE PROBLEMS IN ELLIPTICAL DOMAINS
Please cite this article as: Staisław Kula, Method of fudametal solutios for Helmholtz eigevalue problems i elliptical domais, Scietific Research of the Istitute of Mathematics ad Computer Sciece, 009,
More information8. СОВЕТУВАЊЕ. Охрид, септември ANALYSIS OF NO LOAD APPARENT POWER AND FREQUENCY SPECTRUM OF MAGNETIZING CURRENT FOR DIFFERENT CORE TYPES
8. СОВЕТУВАЊЕ Охрид, 22 24 септември Leoardo Štrac Frajo Keleme Kočar Power Trasformers Ltd. ANALYSS OF NO LOAD APPARENT POWER AND FREQENCY SPECTRM OF MAGNETZNG CRRENT FOR DFFERENT CORE TYPES ABSTRACT
More informationFilter banks. Separately, the lowpass and highpass filters are not invertible. removes the highest frequency 1/ 2and
Filter bas Separately, the lowpass ad highpass filters are ot ivertible T removes the highest frequecy / ad removes the lowest frequecy Together these filters separate the sigal ito low-frequecy ad high-frequecy
More informationWave Phenomena Physics 15c
Wave Pheomea Physics 5c Lecture Fourier Aalysis (H&L Sectios 3. 4) (Georgi Chapter ) Admiistravia! Midterm average 68! You did well i geeral! May got the easy parts wrog, e.g. Problem (a) ad 3(a)! erm
More informationStatistical Intervals for a Single Sample
3/5/06 Applied Statistics ad Probability for Egieers Sixth Editio Douglas C. Motgomery George C. Ruger Chapter 8 Statistical Itervals for a Sigle Sample 8 CHAPTER OUTLINE 8- Cofidece Iterval o the Mea
More informationActivity 3: Length Measurements with the Four-Sided Meter Stick
Activity 3: Legth Measuremets with the Four-Sided Meter Stick OBJECTIVE: The purpose of this experimet is to study errors ad the propagatio of errors whe experimetal data derived usig a four-sided meter
More informationACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS 1 MATH00030 SEMESTER / Statistics
ACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS 1 MATH00030 SEMESTER 1 018/019 DR. ANTHONY BROWN 8. Statistics 8.1. Measures of Cetre: Mea, Media ad Mode. If we have a series of umbers the
More informationThe Discrete Fourier Transform
The iscrete Fourier Trasform The discrete-time Fourier trasform (TFT) of a sequece is a cotiuous fuctio of!, ad repeats with period. I practice we usually wat to obtai the Fourier compoets usig digital
More informationFrequency Domain Filtering
Frequecy Domai Filterig Raga Rodrigo October 19, 2010 Outlie Cotets 1 Itroductio 1 2 Fourier Represetatio of Fiite-Duratio Sequeces: The Discrete Fourier Trasform 1 3 The 2-D Discrete Fourier Trasform
More informationRandom Variables, Sampling and Estimation
Chapter 1 Radom Variables, Samplig ad Estimatio 1.1 Itroductio This chapter will cover the most importat basic statistical theory you eed i order to uderstad the ecoometric material that will be comig
More informationa b c d e f g h Supplementary Information
Supplemetary Iformatio a b c d e f g h Supplemetary Figure S STM images show that Dark patters are frequetly preset ad ted to accumulate. (a) mv, pa, m ; (b) mv, pa, m ; (c) mv, pa, m ; (d) mv, pa, m ;
More informationA STUDY OF VIBRATION MEASURING AND FATIGUE ANALYSIS FOR CANTILEVER BEAMS
Joural of Techology, Vol. 3, No., pp. 47-56 (07) 47, * LabView ANSYS A STUDY OF VIBRATION MEASURING AND FATIGUE ANALYSIS FOR CANTILEVER BEAMS Yuo-Ter Tsai, * Hsie-Yag Li Departmet of Mechaical Egieerig
More informationMathematical Description of Discrete-Time Signals. 9/10/16 M. J. Roberts - All Rights Reserved 1
Mathematical Descriptio of Discrete-Time Sigals 9/10/16 M. J. Roberts - All Rights Reserved 1 Samplig ad Discrete Time Samplig is the acquisitio of the values of a cotiuous-time sigal at discrete poits
More informationModule 18 Discrete Time Signals and Z-Transforms Objective: Introduction : Description: Discrete Time Signal representation
Module 8 Discrete Time Sigals ad Z-Trasforms Objective:To uderstad represetig discrete time sigals, apply z trasform for aalyzigdiscrete time sigals ad to uderstad the relatio to Fourier trasform Itroductio
More informationADVANCED DIGITAL SIGNAL PROCESSING
ADVANCED DIGITAL SIGNAL PROCESSING PROF. S. C. CHAN (email : sccha@eee.hku.hk, Rm. CYC-702) DISCRETE-TIME SIGNALS AND SYSTEMS MULTI-DIMENSIONAL SIGNALS AND SYSTEMS RANDOM PROCESSES AND APPLICATIONS ADAPTIVE
More informationMicroscopic Theory of Transport (Fall 2003) Lecture 6 (9/19/03) Static and Short Time Properties of Time Correlation Functions
.03 Microscopic Theory of Trasport (Fall 003) Lecture 6 (9/9/03) Static ad Short Time Properties of Time Correlatio Fuctios Refereces -- Boo ad Yip, Chap There are a umber of properties of time correlatio
More informationSequences A sequence of numbers is a function whose domain is the positive integers. We can see that the sequence
Sequeces A sequece of umbers is a fuctio whose domai is the positive itegers. We ca see that the sequece 1, 1, 2, 2, 3, 3,... is a fuctio from the positive itegers whe we write the first sequece elemet
More informationComparison Study of Series Approximation. and Convergence between Chebyshev. and Legendre Series
Applied Mathematical Scieces, Vol. 7, 03, o. 6, 3-337 HIKARI Ltd, www.m-hikari.com http://d.doi.org/0.988/ams.03.3430 Compariso Study of Series Approimatio ad Covergece betwee Chebyshev ad Legedre Series
More informationEE422G Homework #13 (12 points)
EE422G Homework #1 (12 poits) 1. (5 poits) I this problem, you are asked to explore a importat applicatio of FFT: efficiet computatio of covolutio. The impulse respose of a system is give by h(t) (.9),1,2,,1
More informationEVALUATION OF GLASS FIBER/EPOXY INTERFACIAL STRENGTH BY THE CRUCIFORM SPECIMEN METHOD
EVALUATION OF GLASS FIBER/EPOX INTERFACIAL STRENGTH B THE CRUCIFORM SPECIMEN METHOD Ju KOANAGI, Hajime KATO, Akihiro KASHIMA, uichi IGARASHI, Keichi WATANABE 3, Ichiro UENO 4 ad Shiji OGIHARA 4 Istitute
More informationENGI 4421 Confidence Intervals (Two Samples) Page 12-01
ENGI 44 Cofidece Itervals (Two Samples) Page -0 Two Sample Cofidece Iterval for a Differece i Populatio Meas [Navidi sectios 5.4-5.7; Devore chapter 9] From the cetral limit theorem, we kow that, for sufficietly
More information11 Correlation and Regression
11 Correlatio Regressio 11.1 Multivariate Data Ofte we look at data where several variables are recorded for the same idividuals or samplig uits. For example, at a coastal weather statio, we might record
More informationSeptember 2012 C1 Note. C1 Notes (Edexcel) Copyright - For AS, A2 notes and IGCSE / GCSE worksheets 1
September 0 s (Edecel) Copyright www.pgmaths.co.uk - For AS, A otes ad IGCSE / GCSE worksheets September 0 Copyright www.pgmaths.co.uk - For AS, A otes ad IGCSE / GCSE worksheets September 0 Copyright
More informationProblem 1. Problem Engineering Dynamics Problem Set 9--Solution. Find the equation of motion for the system shown with respect to:
2.003 Egieerig Dyamics Problem Set 9--Solutio Problem 1 Fid the equatio of motio for the system show with respect to: a) Zero sprig force positio. Draw the appropriate free body diagram. b) Static equilibrium
More information