ANALYSIS OF REAL-TIME DATA FROM INSTRUMENTED STRUCTURES
|
|
- Mervyn Terence Beasley
- 5 years ago
- Views:
Transcription
1 3 h World Conference on Earhquake Engineering Vancouver, B.C., Canada Augus -6, 004 Paper No. 93 ANAYSIS OF REA-IME DAA FROM INSRUMENED SRUCURES Erdal SAFAK SUMMARY Analyses of real-ie daa fro insruened srucures require ools and echniques ha are differen han hose used for riggered daa. Real-ie onioring involves analysis of vibraion daa ha are low in apliudes and high in noise. Moreover, real-ie analysis requires ehods ha can be applied in real ie and recursively. his paper presens soe of he basics in processing and analyzing real-ie vibraion daa. he opics discussed include uilizaion of running ie windows, racking ean and ean-square values, derending, filering, and syse idenificaion. he ehods presened are based on he concep of opial adapive filering and recursive leas-squares esiaion. INRODUCION In ajoriy of srucures insruened for vibraion onioring, he recorders are se o rigger only during large-apliude oions, such as hose generaed by oderae o large earhquakes. Recenly, a few srucures have been insruened o provide coninuous daa in real ie, recording no only large-apliude oions bu also sall-apliude oions generaed by abien loads. he ain objecive in coninuous onioring is o rack any changes in srucural characerisics in order o deec daage. his ype of onioring is coonly known as Srucural Healh Monioring. Fourier-based specral analysis ehods have been he priary ool o analyze daa fro insruened srucures. Mehods ha are used o analyze large-apliude vibraion daa are no always appropriae o analyze abien vibraion daa. In ers of daa analysis, large-apliude records (e.g., earhquake records) have uch beer signal-o-noise raios han abien records. However, hey are ypically ransien, exhibi ie varying (i.e., nonsaionary) eporal and frequency characerisics, and can have nonlinear properies. In coparison, abien records have he advanages of being infiniely long duraion, saionary, and in os cases linear. However, hey conain high level of noise. he long duraion and saionariy of abien daa allow us o uilize saisical signal processing ools, which can copensae for he adverse effecs of low signal-o-noise raios. his paper presens soe of he basics in real-ie daa analysis, including he uilizaion of running ie windows, racking ean and ean-square values, derending, filering, and syse idenificaion. he ephasis is on ehods ha can be applied recursively, which are ore appropriae for real-ie analysis. Research Srucural Engineer, U.S. Geological Survey, Pasadena, California, USA. E-ail: safak@usgs.gov
2 RACKING IME VARIAIONS IN SIGNA PROPERIES Coninuous onioring requires coninuous and auoaed daa processing and analysis. he ehods ha are used for analysis should be able o adap and accoun for any changes in signal characerisics. he siples and os sraighforward approach o analyze coninuous daa is he block-daa approach. In his ehod, he records are handled in blocks of specified lengh. Each block is processed and analyzed as soon as i is full, and while he daa for he nex block are being acquired. More efficien ways o analyze coninuous daa can be developed by uilizing running ie windows. Running windows are in essence weighing funcions ha ephasize recen daa, while gradually deephasizing pas daa. he windows ensure ha any propery calculaed fro daa conains easureens ha are relevan o he curren sae of he srucure. he wo widely used weighing funcions are exponenially decaying windows and sliding recangular windows. he exponenially decaying window is defined as λ i wi (, ) wih i 0,,,, and wi (, ) λ = = = () λ where λ is known as he forgeing facor wih 0.0 < λ <.0. he window applies exponenially decaying weighs o pas daa poins. he eory ie consan, N 0, of he window is defined as he nuber of sapling poins over which he characerisics of he srucure can be assued o reain consan. I can be approxiaed as: N 0 λ () his equaion provides a siple crierion for he selecion of λ.. ypically, λ=0.900 ~ Sliding recangular windows consider only a liied nuber of pas daa poins wih equal weighs. I is defined as: wi (, ) = for i wih 0,,,, and wi (, ) = (3) where is he window lengh and he eory ie consan. RACKING OF MEAN VAUE he ean value in a vibraion record represens he saic coponen of he srucure s response. Norally, he ean value of vibraion records should be zero. In real-ie onioring, he ean value ay flucuae around zero due o various facors, such as iperfecions in sensors, environenal facors (e.g., wind, rain, ec.), naural changes in srucure (e.g., srucural odificaions, change in loads, aging, ec.), or peranen daage afer an exree even (e.g., inelasic deforaions during an earhquake). herefore, i is iporan ha he changes in ean value are racked accuraely. We can consider he calculaion of ean value of a signal as a proble of a weighed leassquares fi of a consan o he record by iniizing he following error funcion:
3 [ ] (4) ε() = wi (,) xi () () where x(i) is he signal, w(,i) is he weighing funcion, and () is he ean value. For he weighing funcions described above, he iniizaion gives he following expressions for he ean value a ie : λ i For exponenially decaying window: ( ) = λ x() i λ For sliding regangular window: ( ) = xi ( ) + (5) Above expressions can be pu ino a recursive for such ha he ean value a ie is calculaed fro he ean value a ie -. he recursive fors are (Safak, 004) For exponenially decaying window: ( ) = λ ( ) + ( λ) x( i) For sliding regangular window: ( ) = ( ) + y() y( ) [ ] (6) Copuaionally, he recursive fors are ore appropriae for real-ie daa. RACKING OF MEAN-SQUARE VAUE Mean-square (MS) value is anoher iporan paraeer ha needs o be racked in real-ie onioring. MS value provides inforaion on average vibraion apliudes and is one of he key paraeers o deec changes in srucural response and exciaion. Siilar o he ean value, he MS value, s(), a ie can be calculaed by iniizing he following equaion: ε ( ) = wi (, ) x ( i) s ( ) (7) which resuls in he following recursive expressions for s() (Safak, 004): For exponenially decaying window: s( ) = λ s ( ) + ( λ) x( i) For sliding regangular window: s( ) = s( ) + y ( ) y ( ) (8) DERENDING Derending, also known as baseline correcion, reoves a linear rend fro records. Again, for real-ie daa i should be applied by uilizing running ie windows, so ha he rend ha is being reoved is he curren one. he equaions for derending are developed by deerining he bes sraigh-line fi, i.e., y()=a x() + b, o daa and incorporaing he ie windows of Eqs. and 3. he coefficiens a and b of he sraigh line are deerined by iniizing he following error funcion wih respec o a and b:
4 ε() = wi (,) xi () a xi () b [ ] (9) he explici expressions for a and b can be found in Safak (004). Noe ha because of he er b in he equaion, derending also reoves he ean value fro he daa. FIERING Band-pass filering Filering is required o reduce noise and eliinae unwaned frequency coponens in records. he sandard approach o noise reducion has been o use band-pass filers. Band-pass filers reove he frequency coponens ha are believed o be doinaed by noise (ypically he very low and very high ends of he frequency band). Alhough he noise is sill presen in he reaining frequency band, i is assued ha he signal-o-noise raio is high enough so ha he noise can be ignored. Band-pass filering, as well as low- or high-pass filering, can be done by using recursive iedoain filers of he following for (e.g., reer, 976): n k l (0) k= l= 0 y() = a y ( k) + b x ( l) where y() is he filered signal and x() is he original signal. he filer paraeers a k and b l and he filer orders and n are deerined by he corner frequencies, and he rae of decay of filer apliudes near he corner frequencies. A unique se of paraeers can be calculaed for any specified filer. For exaple, for a second order band-pass Buerworh filer wih corner frequencies of 0. Hz and 0 Hz and a sapling inerval of 0.0 second, he filer paraeers, calculaed by using MAAB (MahWorks, 003), are: a= , a= , b0=0.05, b= 0, b = () Opial filering Mos of real-ie vibraion daa consis of abien vibraions, which are characerized by low apliudes and low signal-o-noise raios. Analyses of such noisy signals require he applicaions of sophisicaed filers so ha he aoun of inforaion exraced fro he daa is axiized. A general class of such filers is known as opial filers. he key difference beween band-pass filers and opial filers is ha opial filers reduce he noise over he enire frequency band. We will presen he general concep of opial filering based on he principle of Recursive eas Squares (RS) approxiaion. Siilar filers can be developed by using he eas Mean Square (MS) approxiaion, as well as by uilizing sae-space forulaions, leading o Kalan Filers (KF) and Exended Kalan Filers (EKF). Due o space liiaions, hese alernaive fors will no be discussed; hey can be found in a large nuber of references ha are available in he lieraure (e.g., Widrow and Searns, 985; Braer and Siffling, 989). he basic philosophy in opial filering is ha he recorded response is coposed of a signal (i.e., he acual response) plus noise. Since he signal represens he response of a dynaic syse o soe exciaion, i is reasonable o expec ha response values a discree ies will follow a paern (e.g., a daped sinusoid), and herefore will be correlaed wih each oher. For linear syses, his paern can be represened by wriing he response a ie as a linear cobinaion
5 of pas responses; ha is x () = ak x ( k) + n () () k= he firs er on he righ-hand side represens he correlaed porion of he record, i.e., he signal, and he second er he uncorrelaed porion of he signal, i.e., he noise. Hence, he proble becoes deerining he opial values of a k such ha he prediced value of x() is as close o is recorded value as possible. Assuing ha n() is a zero-ean rando process and all he pas values of x() are known, he bes esiae, E[x()], of x() can be wrien as Ex [ ( )] = ak x ( k) (3) k= he error in he esiae, ε(), is defined as he difference beween he recorded value and he esiaed value of x(); and he oal error, V, is he su of squared errors over he record lengh, N: N k ε (4) k= = ε() = x() + a x( k) and V = () where N is he nuber of sapling poins. By aking V / ak = 0 for k =,, we obain he following arix equaion for a k : xx K xx xx a r (0) r ( ) r () M = M O M M a rxx( ) rxx(0) r xx( ) (5) where rxx ( j) = xkxk ( ) ( j) (6) k= is he auocorrelaion funcion of x(). By defining he following vecors and he arix, we can wrie i in a ore copac for as: a rxx () rxx (0) K rxx ( ) θ = R f where θ =, f, R M = M = M O M a rxx( ) r xx( ) rxx(0) (7)
6 A key paraeer ha needs o be seleced o apply hese equaions is, he filer order. A filer wih oo sall does no accuraely represen he signal, whereas a filer wih oo large ay ry o represen noise as well as he signal. here several crieria available in he lieraure o selec (see, Sodersro, 987). A sipler and ore sraighforward selecion can be ade by ploing he variaion of V = ε () wih (Safak, 004). his su ypically shows a fas drop wih increasing, and hen level off. he value where he su sars o level off can be aken as he opial filer order. Eq. 4 is based on he iniizaion of forward squared predicion errors because we used he pas values of x() o predic is curren value. Siilar equaions can be developed by iniizing he backward predicion errors. We can also iniize he su of forward and backward predicion errors, which resuls in he so-called Burg s ehod (Burg, 968). he relaionships beween he forward and backward iniizaion forulaions lead o very fas recursive algorihs for he soluion of filer coefficiens, such as evinson-durbin algorih (evinson, 947; Durbin, 959) and laice filers (Griffihs, 977; Makhoul, 977). An iporan assupion ade in he derivaion of Eq. 4 is ha x() is a saionary signal. In oher words, he eporal and frequency characerisics of x() does no change significanly wih ie, and herefore he auocorrelaion funcion r xx (Eq. 5) is he funcion of he ie lag only beween he wo coponens. he assupion of saionariy is appropriae for vibraions under abien forces and wind loads, bu no for vibraions under ransien loads such as earhquakes or blas loads. Adapive filering Adapive filering is required when he eporal and frequency characerisics of he signal change wih ie, i.e., he assupion of saionariy is no longer valid. Equaions for adapive filering can be developed by siple incorporaing a weighing funcion (e.g., Eq. or 3) in he error funcion for he leas-squares esiaion. (Eq. 4). ha leads o he following error funcion: k (8) k= V() = w(,) i x() i + a x( i k) Noe ha V is now a funcion of ie. By denoing φ () = [ x( ), x( ),, x( )] (9) and using he exponenially decaying for for w(,i), he paraeer vecor θ = [ a, a,, a ] can be calculaed recursively by he following se of equaions (jung, 999):
7 θ() = θ( ) + () x() φ () θ( ) P ( ) φ ( ) () = λ + φ () P ( ) φ( P ( ) φ( ) φ ( P ) ( ) P () = P ( ) λ λ + φ () P ( ) φ( (0) As indicaed, he filer paraeers now are ie varying. SYSEM IDENIFICAION Syse idenificaion refers o he exracion of srucural characerisics fro recorded signals. he basic approach o syse idenificaion is very siilar o ha used for opial filering. We search for an opial filer ha convers he recorded exciaion signal ino he recorded response signals. he filer idenified represens he dynaic characerisics of he srucure. As will be shown laer, here is a one-o-one correspondence beween he filer paraeers and he odal frequencies and daping raios of he srucure. If we do no have access o he exciaion signal (e.g., ground acceleraions in he case of an earhquake-induced vibraions), we assue ha he exciaion is a zero-ean wind-band rando process. his case is known as oupu-only idenificaion. We presen he basics of syse idenificaion for boh oupu-only and inpuoupu cases below. Oupu-only syse idenificaion: he forulaion for oupu-only syse idenificaion is idenical o ha for opial filering. We search for a filer ha, when applied o daa, will give a residual ha is as close o a zero-ean whie-noise sequence as possible. Maheaically, his leads o he sae error funcion and he iniizaion proble defined by Eq. 4 (for saionary signals) or by Eq. 8 (for nonsaionary signals). he filer paraeers are calculaed fro he resuling equaions, Eq. 5 or Eq. 0. In order o calculae he dynaic properies of he srucure fro he filer coefficiens, we firs calculae he roos of he following equaion: + az + + a z = 0 () he roos of he equaion are called he poles, p k, of he ransfer funcion for he srucure. For srucures wih posiive daping, all he roos are in coplex conjugae pairs, resuling in / disinc odal frequencies, f k, and daping raios, ξ k. he f k, and ξ k are calculaed fro he following equaions (Safak, 99): ( p ) + ( p ) / ( p ) ln / od( ) ln / od( ) ξk =, fk = in Hz, k =,, / () arg( p ) ln / od( ) πξk k where od( ) and arg( ) denoe he odulus and arguen of he coplex poles, and is he sapling inerval in seconds. Inpu-oupu syse idenificaion: If he inpu (i.e., he exciaion) signal is also recorded, we deerine filer paraeers by aking
8 he exciaion signal as he inpu o he filer and aching he oupu-signal in he leas-squares sense. he filer equaion for he inpu-oupu case is: n k l (3) k= l= 0 x () = a x ( k) + b u ( l) + n () where u() denoes he inpu signal. he filer paraeers, a k and b l can be deerined recursively fro Eq. 0, by siply changing he definiions of he vecors θ() and φ() as follows: θ ( ) = [ a, a,, a, b0, b,, bn] φ( ) = [ x ( ), x ( ),, x ( ), u ( ), u ( ),, u ( n)] he naural frequencies and daping raios are calculaed fro a k coefficiens as discussed above; b l coefficiens are relaed o he odal paricipaion facors (for deails, see Safak, 988; 99). Since θ() is ie-varying, his algorih can rack any changes in he odal frequencies and daping raios. SUMMARY AND CONCUSIONS Analyses of real-ie daa fro insruened srucures require fas and efficien algorihs. Coninuous records differ fro riggered records in ha hey are ypically very low in apliudes and high in noise. Mehods used o analyze real-ie daa should be able o exrac he signal buried in noise, and iniize he effecs of noise over he enire frequency band. he ehods should also be in a fora ha can be applied recursively in real ie, and be able o adap o rack any ie variaions in signal properies. Opial adapive filers offer one of he bes ools o analyze real-ie daa. hese filers exrac he signal fro record by searching hidden correlaion paerns in he daa. Such filers are used no only for noise reducion over he enire frequency band, bu also for syse idenificaion, because he forulaion for boh leads o he iniizaion of sae error funcions. Opial adapive filers can be expressed in a large nuber of differen fors, and known under differen naes depending on he forulaion and error iniizaion algorihs used in he developen, such as RS filers, MS filers, and Kalan Filers. REFERENCES. Braer, K. and Siffling, G. Kalan-Bucy Filers, Arech House, Norwood, MA, Burg, J.P. A new analysis echnique for ie-series daa, NAO Advanced Sudy Insiue on Signal Processing, Enschede, he Neherlands, Durbin, J. Efficien esiaors of paraeers in oving average odels, Bioerica, 46:306-36, Griffihs,.J. A coninuously adapive filer ipleened as a laice srucure, Proc. 977 IEEE In. Conf. on Acousics, Speech and Signal Processing, , evinson, N. he Wiener rs error crierion filer design and predicion, J. Mah. (4)
9 Phys., 5:6-78, jung,. Syse Idenificaion: heory for he User, Prenice Hall PR, Upper Saddle River, NJ, Makhoul, J. Sable and efficien laice ehods for linear predicion, IEEE rans. Acousics, Speech and Signal Processing, ASSP-5:43-48, he MahWorks, Inc. MAAB he anguage of echnical Copuing, Version 6, Naick, MA, Safak, E. (988). Analysis of recordings in srucural engineering: Adapive filering, predicion, and conrol, Open-File Repor , U. S. Geological Survey, Menlo Park, California. 0. Safak, E. (99). Idenificaion of linear srucures using discree-ie filers, Journal of Srucural Engineering, ASCE, Vol.7, No.0, Ocober 99, pp Safak, E. (004). Analysis of real-ie vibraion daa (in preparaion).. Sodersro,. Model srucure deerinaion, in Encyclopedia of Syses and Conrol (M. Singh, ed.), Pergaon Press, Elsford, NY, reer, S.A. Inroducion o Discree-ie Signal Processing, John Wiley and Sons, New York, NY, Widrow, B. and Searns, S.D. Adapive Signal Processing, Prenice-Hall, Inc., Englewood Cliffs, NJ, 985.
Introduction to Numerical Analysis. In this lesson you will be taken through a pair of techniques that will be used to solve the equations of.
Inroducion o Nuerical Analysis oion In his lesson you will be aen hrough a pair of echniques ha will be used o solve he equaions of and v dx d a F d for siuaions in which F is well nown, and he iniial
More informationTIME DELAY BASEDUNKNOWN INPUT OBSERVER DESIGN FOR NETWORK CONTROL SYSTEM
TIME DELAY ASEDUNKNOWN INPUT OSERVER DESIGN FOR NETWORK CONTROL SYSTEM Siddhan Chopra J.S. Laher Elecrical Engineering Deparen NIT Kurukshera (India Elecrical Engineering Deparen NIT Kurukshera (India
More information1. Calibration factor
Annex_C_MUBDandP_eng_.doc, p. of pages Annex C: Measureen uncerainy of he oal heigh of profile of a deph-seing sandard ih he sandard deviaion of he groove deph as opography er In his exaple, he uncerainy
More informationOn the approximation of particular solution of nonhomogeneous linear differential equation with Legendre series
The Journal of Applied Science Vol. 5 No. : -9 [6] วารสารว ทยาศาสตร ประย กต doi:.446/j.appsci.6.8. ISSN 53-785 Prined in Thailand Research Aricle On he approxiaion of paricular soluion of nonhoogeneous
More informationApplication of Homotopy Analysis Method for Solving various types of Problems of Partial Differential Equations
Applicaion of Hooopy Analysis Mehod for olving various ypes of Probles of Parial Differenial Equaions V.P.Gohil, Dr. G. A. anabha,assisan Professor, Deparen of Maheaics, Governen Engineering College, Bhavnagar,
More informationThus the force is proportional but opposite to the displacement away from equilibrium.
Chaper 3 : Siple Haronic Moion Hooe s law saes ha he force (F) eered by an ideal spring is proporional o is elongaion l F= l where is he spring consan. Consider a ass hanging on a he spring. In equilibriu
More informationFourier Series & The Fourier Transform. Joseph Fourier, our hero. Lord Kelvin on Fourier s theorem. What do we want from the Fourier Transform?
ourier Series & The ourier Transfor Wha is he ourier Transfor? Wha do we wan fro he ourier Transfor? We desire a easure of he frequencies presen in a wave. This will lead o a definiion of he er, he specru.
More informationVehicle Arrival Models : Headway
Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where
More informationConnectionist Classifier System Based on Accuracy in Autonomous Agent Control
Connecionis Classifier Syse Based on Accuracy in Auonoous Agen Conrol A S Vasilyev Decision Suppor Syses Group Riga Technical Universiy /4 Meza sree Riga LV-48 Lavia E-ail: serven@apollolv Absrac In his
More informationMultivariate Auto-Regressive Model for Groundwater Flow Around Dam Site
Mulivariae uo-regressive Model for Groundwaer Flow round Da Sie Yoshiada Mio ), Shinya Yaaoo ), akashi Kodaa ) and oshifui Masuoka ) ) Dep. of Earh Resources Engineering, Kyoo Universiy, Kyoo, 66-85, Japan.
More informationChapter 15. Time Series: Descriptive Analyses, Models, and Forecasting
Chaper 15 Time Series: Descripive Analyses, Models, and Forecasing Descripive Analysis: Index Numbers Index Number a number ha measures he change in a variable over ime relaive o he value of he variable
More informationKriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 Kriging Models Predicing Arazine Concenraions in Surface Waer Draining Agriculural Waersheds Paul L. Mosquin, Jeremy Aldworh, Wenlin Chen Supplemenal Maerial Number
More informationLecture 23 Damped Motion
Differenial Equaions (MTH40) Lecure Daped Moion In he previous lecure, we discussed he free haronic oion ha assues no rearding forces acing on he oving ass. However No rearding forces acing on he oving
More informationDiebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles
Diebold, Chaper 7 Francis X. Diebold, Elemens of Forecasing, 4h Ediion (Mason, Ohio: Cengage Learning, 006). Chaper 7. Characerizing Cycles Afer compleing his reading you should be able o: Define covariance
More informationSub Module 2.6. Measurement of transient temperature
Mechanical Measuremens Prof. S.P.Venkaeshan Sub Module 2.6 Measuremen of ransien emperaure Many processes of engineering relevance involve variaions wih respec o ime. The sysem properies like emperaure,
More information1 Widrow-Hoff Algorithm
COS 511: heoreical Machine Learning Lecurer: Rob Schapire Lecure # 18 Scribe: Shaoqing Yang April 10, 014 1 Widrow-Hoff Algorih Firs le s review he Widrow-Hoff algorih ha was covered fro las lecure: Algorih
More informationChapter 9 Sinusoidal Steady State Analysis
Chaper 9 Sinusoidal Seady Sae Analysis 9.-9. The Sinusoidal Source and Response 9.3 The Phasor 9.4 pedances of Passive Eleens 9.5-9.9 Circui Analysis Techniques in he Frequency Doain 9.0-9. The Transforer
More informationOBJECTIVES OF TIME SERIES ANALYSIS
OBJECTIVES OF TIME SERIES ANALYSIS Undersanding he dynamic or imedependen srucure of he observaions of a single series (univariae analysis) Forecasing of fuure observaions Asceraining he leading, lagging
More informationPhasor Estimation Algorithm Based on the Least Square Technique during CT Saturation
Journal of Elecrical Engineering & Technology Vol. 6, No. 4, pp. 459~465, 11 459 DOI: 1.537/JEET.11.6.4. 459 Phasor Esiaion Algorih Based on he Leas Square Technique during CT Sauraion Dong-Gyu Lee*, Sang-Hee
More informationCHAPTER 2 Signals And Spectra
CHAPER Signals And Specra Properies of Signals and Noise In communicaion sysems he received waveform is usually caegorized ino he desired par conaining he informaion, and he undesired par. he desired par
More informationb denotes trend at time point t and it is sum of two
Inernaional Conference on Innovaive Applicaions in Engineering and Inforaion echnology(iciaei207) Inernaional Journal of Advanced Scienific echnologies,engineering and Manageen Sciences (IJASEMSISSN: 2454356X)
More informationNature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time.
Supplemenary Figure 1 Spike-coun auocorrelaions in ime. Normalized auocorrelaion marices are shown for each area in a daase. The marix shows he mean correlaion of he spike coun in each ime bin wih he spike
More informationTP A.14 The effects of cut angle, speed, and spin on object ball throw
echnical proof echnical proof TP A.14 The effecs of cu angle, speed, and spin on objec ball hrow supporing: The Illusraed Principles of Pool and illiards hp://billiards.colosae.edu by Daid G. Alciaore,
More informationSTATE-SPACE MODELLING. A mass balance across the tank gives:
B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing
More informationStable block Toeplitz matrix for the processing of multichannel seismic data
Indian Journal of Marine Sciences Vol. 33(3), Sepember 2004, pp. 215-219 Sable block Toepliz marix for he processing of mulichannel seismic daa Kiri Srivasava* & V P Dimri Naional Geophysical Research
More informationModal identification of structures from roving input data by means of maximum likelihood estimation of the state space model
Modal idenificaion of srucures from roving inpu daa by means of maximum likelihood esimaion of he sae space model J. Cara, J. Juan, E. Alarcón Absrac The usual way o perform a forced vibraion es is o fix
More informationDecision Tree Learning. Decision Tree Learning. Decision Trees. Decision Trees: Operation. Blue slides: Mitchell. Orange slides: Alpaydin Humidity
Decision Tree Learning Decision Tree Learning Blue slides: Michell Oulook Orange slides: Alpaydin Huidiy Sunny Overcas Rain ral Srong Learn o approxiae discree-valued arge funcions. Sep-by-sep decision
More informationCHAPTER 10 VALIDATION OF TEST WITH ARTIFICAL NEURAL NETWORK
175 CHAPTER 10 VALIDATION OF TEST WITH ARTIFICAL NEURAL NETWORK 10.1 INTRODUCTION Amongs he research work performed, he bes resuls of experimenal work are validaed wih Arificial Neural Nework. From he
More information( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is
UNIT IMPULSE RESPONSE, UNIT STEP RESPONSE, STABILITY. Uni impulse funcion (Dirac dela funcion, dela funcion) rigorously defined is no sricly a funcion, bu disribuion (or measure), precise reamen requires
More informationTime series Decomposition method
Time series Decomposiion mehod A ime series is described using a mulifacor model such as = f (rend, cyclical, seasonal, error) = f (T, C, S, e) Long- Iner-mediaed Seasonal Irregular erm erm effec, effec,
More informationNavneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi
Creep in Viscoelasic Subsances Numerical mehods o calculae he coefficiens of he Prony equaion using creep es daa and Herediary Inegrals Mehod Navnee Saini, Mayank Goyal, Vishal Bansal (23); Term Projec
More informationMapping in Dynamic Environments
Mapping in Dynaic Environens Wolfra Burgard Universiy of Freiburg, Gerany Mapping is a Key Technology for Mobile Robos Robos can robusly navigae when hey have a ap. Robos have been shown o being able o
More informationOn Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature
On Measuring Pro-Poor Growh 1. On Various Ways of Measuring Pro-Poor Growh: A Shor eview of he Lieraure During he pas en years or so here have been various suggesions concerning he way one should check
More informationSIGNALS AND SYSTEMS LABORATORY 8: State Variable Feedback Control Systems
SIGNALS AND SYSTEMS LABORATORY 8: Sae Variable Feedback Conrol Syses INTRODUCTION Sae variable descripions for dynaical syses describe he evoluion of he sae vecor, as a funcion of he sae and he inpu. There
More informationGeorey E. Hinton. University oftoronto. Technical Report CRG-TR February 22, Abstract
Parameer Esimaion for Linear Dynamical Sysems Zoubin Ghahramani Georey E. Hinon Deparmen of Compuer Science Universiy oftorono 6 King's College Road Torono, Canada M5S A4 Email: zoubin@cs.orono.edu Technical
More informationRC, RL and RLC circuits
Name Dae Time o Complee h m Parner Course/ Secion / Grade RC, RL and RLC circuis Inroducion In his experimen we will invesigae he behavior of circuis conaining combinaions of resisors, capaciors, and inducors.
More informationSome Basic Information about M-S-D Systems
Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,
More informationPractice Problems - Week #4 Higher-Order DEs, Applications Solutions
Pracice Probles - Wee #4 Higher-Orer DEs, Applicaions Soluions 1. Solve he iniial value proble where y y = 0, y0 = 0, y 0 = 1, an y 0 =. r r = rr 1 = rr 1r + 1, so he general soluion is C 1 + C e x + C
More informationKINEMATICS IN ONE DIMENSION
KINEMATICS IN ONE DIMENSION PREVIEW Kinemaics is he sudy of how hings move how far (disance and displacemen), how fas (speed and velociy), and how fas ha how fas changes (acceleraion). We say ha an objec
More information3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon
3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of
More informationMechanical Fatigue and Load-Induced Aging of Loudspeaker Suspension. Wolfgang Klippel,
Mechanical Faigue and Load-Induced Aging of Loudspeaker Suspension Wolfgang Klippel, Insiue of Acousics and Speech Communicaion Dresden Universiy of Technology presened a he ALMA Symposium 2012, Las Vegas
More informationM x t = K x F t x t = A x M 1 F t. M x t = K x cos t G 0. x t = A x cos t F 0
Forced oscillaions (sill undaped): If he forcing is sinusoidal, M = K F = A M F M = K cos G wih F = M G = A cos F Fro he fundaenal heore for linear ransforaions we now ha he general soluion o his inhoogeneous
More information23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes
Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals
More informationChapter 2. First Order Scalar Equations
Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.
More informationLINEAR TIME SERIES ANALYSIS
LINEAR TIME ERIE ANALYI ivaraane, N. Indian Agriculural aisics Research Insiue, New Delhi-00 Tie series daa The variable conaining observaions over ie is called a ie series variable and he daase as ie
More information2.1 Harmonic excitation of undamped systems
2.1 Haronic exciaion of undaped syses Sanaoja 2_1.1 2.1 Haronic exciaion of undaped syses (Vaienaaoan syseein haroninen heräe) The following syse is sudied: y x F() Free-body diagra f x g x() N F() In
More informationHomework 2 Solutions
Mah 308 Differenial Equaions Fall 2002 & 2. See he las page. Hoework 2 Soluions 3a). Newon s secon law of oion says ha a = F, an we know a =, so we have = F. One par of he force is graviy, g. However,
More information5. Stochastic processes (1)
Lec05.pp S-38.45 - Inroducion o Teleraffic Theory Spring 2005 Conens Basic conceps Poisson process 2 Sochasic processes () Consider some quaniy in a eleraffic (or any) sysem I ypically evolves in ime randomly
More information- The whole joint distribution is independent of the date at which it is measured and depends only on the lag.
Saionary Processes Sricly saionary - The whole join disribuion is indeenden of he dae a which i is measured and deends only on he lag. - E y ) is a finie consan. ( - V y ) is a finie consan. ( ( y, y s
More informationLicenciatura de ADE y Licenciatura conjunta Derecho y ADE. Hoja de ejercicios 2 PARTE A
Licenciaura de ADE y Licenciaura conjuna Derecho y ADE Hoja de ejercicios PARTE A 1. Consider he following models Δy = 0.8 + ε (1 + 0.8L) Δ 1 y = ε where ε and ε are independen whie noise processes. In
More information10. State Space Methods
. Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he
More informationUnderwater Target Tracking Based on Gaussian Particle Filter in Looking Forward Sonar Images
Journal of Copuaional Inforaion Syses 6:4 (00) 480-4809 Available a hp://www.jofcis.co Underwaer Targe Tracing Based on Gaussian Paricle Filer in Looing Forward Sonar Iages Tiedong ZHANG, Wenjing ZENG,
More information2.160 System Identification, Estimation, and Learning. Lecture Notes No. 8. March 6, 2006
2.160 Sysem Idenificaion, Esimaion, and Learning Lecure Noes No. 8 March 6, 2006 4.9 Eended Kalman Filer In many pracical problems, he process dynamics are nonlinear. w Process Dynamics v y u Model (Linearized)
More informationLab 10: RC, RL, and RLC Circuits
Lab 10: RC, RL, and RLC Circuis In his experimen, we will invesigae he behavior of circuis conaining combinaions of resisors, capaciors, and inducors. We will sudy he way volages and currens change in
More informationNotes on Kalman Filtering
Noes on Kalman Filering Brian Borchers and Rick Aser November 7, Inroducion Daa Assimilaion is he problem of merging model predicions wih acual measuremens of a sysem o produce an opimal esimae of he curren
More informationIntroduction D P. r = constant discount rate, g = Gordon Model (1962): constant dividend growth rate.
Inroducion Gordon Model (1962): D P = r g r = consan discoun rae, g = consan dividend growh rae. If raional expecaions of fuure discoun raes and dividend growh vary over ime, so should he D/P raio. Since
More informationApplication of Speed Transform to the diagnosis of a roller bearing in variable speed
Applicaion of Speed Transform o he diagnosis of a roller bearing in variable speed Julien Roussel 1, Michel Hariopoulos 1, Edgard Sekko 1, Cécile Capdessus 1 and Jérôme Anoni 1 PRISME laboraory 1 rue de
More informationReferences are appeared in the last slide. Last update: (1393/08/19)
SYSEM IDEIFICAIO Ali Karimpour Associae Professor Ferdowsi Universi of Mashhad References are appeared in he las slide. Las updae: 0..204 393/08/9 Lecure 5 lecure 5 Parameer Esimaion Mehods opics o be
More informationLecture 3: Exponential Smoothing
NATCOR: Forecasing & Predicive Analyics Lecure 3: Exponenial Smoohing John Boylan Lancaser Cenre for Forecasing Deparmen of Managemen Science Mehods and Models Forecasing Mehod A (numerical) procedure
More informationLAPLACE TRANSFORM AND TRANSFER FUNCTION
CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION Professor Dae Ryook Yang Spring 2018 Dep. of Chemical and Biological Engineering 5-1 Road Map of he Lecure V Laplace Transform and Transfer funcions
More informationAn Extension to the Tactical Planning Model for a Job Shop: Continuous-Time Control
An Exension o he Tacical Planning Model for a Job Shop: Coninuous-Tie Conrol Chee Chong. Teo, Rohi Bhanagar, and Sephen C. Graves Singapore-MIT Alliance, Nanyang Technological Univ., and Massachuses Insiue
More informationJoint Spectral Distribution Modeling Using Restricted Boltzmann Machines for Voice Conversion
INTERSPEECH 2013 Join Specral Disribuion Modeling Using Resriced Bolzann Machines for Voice Conversion Ling-Hui Chen, Zhen-Hua Ling, Yan Song, Li-Rong Dai Naional Engineering Laboraory of Speech and Language
More informationIntroduction to Mechanical Vibrations and Structural Dynamics
Inroducion o Mechanical Viraions and Srucural Dynaics The one seeser schedule :. Viraion - classificaion. ree undaped single DO iraion, equaion of oion, soluion, inegraional consans, iniial condiions..
More informationModule 4: Time Response of discrete time systems Lecture Note 2
Module 4: Time Response of discree ime sysems Lecure Noe 2 1 Prooype second order sysem The sudy of a second order sysem is imporan because many higher order sysem can be approimaed by a second order model
More informationFinancial Econometrics Jeffrey R. Russell Midterm Winter 2009 SOLUTIONS
Name SOLUTIONS Financial Economerics Jeffrey R. Russell Miderm Winer 009 SOLUTIONS You have 80 minues o complee he exam. Use can use a calculaor and noes. Try o fi all your work in he space provided. If
More informationINTRODUCTION TO MACHINE LEARNING 3RD EDITION
ETHEM ALPAYDIN The MIT Press, 2014 Lecure Slides for INTRODUCTION TO MACHINE LEARNING 3RD EDITION alpaydin@boun.edu.r hp://www.cmpe.boun.edu.r/~ehem/i2ml3e CHAPTER 2: SUPERVISED LEARNING Learning a Class
More informationLecture 20: Riccati Equations and Least Squares Feedback Control
34-5 LINEAR SYSTEMS Lecure : Riccai Equaions and Leas Squares Feedback Conrol 5.6.4 Sae Feedback via Riccai Equaions A recursive approach in generaing he marix-valued funcion W ( ) equaion for i for he
More information8. Basic RL and RC Circuits
8. Basic L and C Circuis This chaper deals wih he soluions of he responses of L and C circuis The analysis of C and L circuis leads o a linear differenial equaion This chaper covers he following opics
More informationWall. x(t) f(t) x(t = 0) = x 0, t=0. which describes the motion of the mass in absence of any external forcing.
MECHANICS APPLICATIONS OF SECOND-ORDER ODES 7 Mechanics applicaions of second-order ODEs Second-order linear ODEs wih consan coefficiens arise in many physical applicaions. One physical sysems whose behaviour
More informationLinear Response Theory: The connection between QFT and experiments
Phys540.nb 39 3 Linear Response Theory: The connecion beween QFT and experimens 3.1. Basic conceps and ideas Q: How do we measure he conduciviy of a meal? A: we firs inroduce a weak elecric field E, and
More informationViscous Damping Summary Sheet No Damping Case: Damped behaviour depends on the relative size of ω o and b/2m 3 Cases: 1.
Viscous Daping: && + & + ω Viscous Daping Suary Shee No Daping Case: & + ω solve A ( ω + α ) Daped ehaviour depends on he relaive size of ω o and / 3 Cases:. Criical Daping Wee 5 Lecure solve sae BC s
More information23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes
Half-Range Series 2.5 Inroducion In his Secion we address he following problem: Can we find a Fourier series expansion of a funcion defined over a finie inerval? Of course we recognise ha such a funcion
More informationLecture 28: Single Stage Frequency response. Context
Lecure 28: Single Sage Frequency response Prof J. S. Sih Conex In oday s lecure, we will coninue o look a he frequency response of single sage aplifiers, saring wih a ore coplee discussion of he CS aplifier,
More informationL07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS. NA568 Mobile Robotics: Methods & Algorithms
L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS NA568 Mobile Roboics: Mehods & Algorihms Today s Topic Quick review on (Linear) Kalman Filer Kalman Filering for Non-Linear Sysems Exended Kalman Filer (EKF)
More information14 Autoregressive Moving Average Models
14 Auoregressive Moving Average Models In his chaper an imporan parameric family of saionary ime series is inroduced, he family of he auoregressive moving average, or ARMA, processes. For a large class
More informationStationary Time Series
3-Jul-3 Time Series Analysis Assoc. Prof. Dr. Sevap Kesel July 03 Saionary Time Series Sricly saionary process: If he oin dis. of is he same as he oin dis. of ( X,... X n) ( X h,... X nh) Weakly Saionary
More informationSTRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN
Inernaional Journal of Applied Economerics and Quaniaive Sudies. Vol.1-3(004) STRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN 001-004 OBARA, Takashi * Absrac The
More information1. VELOCITY AND ACCELERATION
1. VELOCITY AND ACCELERATION 1.1 Kinemaics Equaions s = u + 1 a and s = v 1 a s = 1 (u + v) v = u + as 1. Displacemen-Time Graph Gradien = speed 1.3 Velociy-Time Graph Gradien = acceleraion Area under
More informationInstitute for Mathematical Methods in Economics. University of Technology Vienna. Singapore, May Manfred Deistler
MULTIVARIATE TIME SERIES ANALYSIS AND FORECASTING Manfred Deisler E O S Economerics and Sysems Theory Insiue for Mahemaical Mehods in Economics Universiy of Technology Vienna Singapore, May 2004 Inroducion
More informationA generalization of the Burg s algorithm to periodically correlated time series
A generalizaion of he Burg s algorihm o periodically correlaed ime series Georgi N. Boshnakov Insiue of Mahemaics, Bulgarian Academy of Sciences ABSTRACT In his paper periodically correlaed processes are
More informationPhysics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle
Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,
More informationV AK (t) I T (t) I TRM. V AK( full area) (t) t t 1 Axial turn-on. Switching losses for Phase Control and Bi- Directionally Controlled Thyristors
Applicaion Noe Swiching losses for Phase Conrol and Bi- Direcionally Conrolled Thyrisors V AK () I T () Causing W on I TRM V AK( full area) () 1 Axial urn-on Plasma spread 2 Swiching losses for Phase Conrol
More informationA Generalization of Student s t-distribution from the Viewpoint of Special Functions
A Generalizaion of Suden s -disribuion fro he Viewpoin of Special Funcions WOLFRAM KOEPF and MOHAMMAD MASJED-JAMEI Deparen of Maheaics, Universiy of Kassel, Heinrich-Ple-Sr. 4, D-343 Kassel, Gerany Deparen
More information12: AUTOREGRESSIVE AND MOVING AVERAGE PROCESSES IN DISCRETE TIME. Σ j =
1: AUTOREGRESSIVE AND MOVING AVERAGE PROCESSES IN DISCRETE TIME Moving Averages Recall ha a whie noise process is a series { } = having variance σ. The whie noise process has specral densiy f (λ) = of
More informationPosition, Velocity, and Acceleration
rev 06/2017 Posiion, Velociy, and Acceleraion Equipmen Qy Equipmen Par Number 1 Dynamic Track ME-9493 1 Car ME-9454 1 Fan Accessory ME-9491 1 Moion Sensor II CI-6742A 1 Track Barrier Purpose The purpose
More informationON THE BEAT PHENOMENON IN COUPLED SYSTEMS
8 h ASCE Specialy Conference on Probabilisic Mechanics and Srucural Reliabiliy PMC-38 ON THE BEAT PHENOMENON IN COUPLED SYSTEMS S. K. Yalla, Suden Member ASCE and A. Kareem, M. ASCE NaHaz Modeling Laboraory,
More informationChapter 2. Models, Censoring, and Likelihood for Failure-Time Data
Chaper 2 Models, Censoring, and Likelihood for Failure-Time Daa William Q. Meeker and Luis A. Escobar Iowa Sae Universiy and Louisiana Sae Universiy Copyrigh 1998-2008 W. Q. Meeker and L. A. Escobar. Based
More informationI-Optimal designs for third degree kronecker model mixture experiments
Inernaional Journal of Saisics and Applied Maheaics 207 2(2): 5-40 ISSN: 2456-452 Mahs 207 2(2): 5-40 207 Sas & Mahs www.ahsjournal.co Received: 9-0-207 Acceped: 20-02-207 Cheruiyo Kipkoech Deparen of
More informationTypes of Exponential Smoothing Methods. Simple Exponential Smoothing. Simple Exponential Smoothing
M Business Forecasing Mehods Exponenial moohing Mehods ecurer : Dr Iris Yeung Room No : P79 Tel No : 788 8 Types of Exponenial moohing Mehods imple Exponenial moohing Double Exponenial moohing Brown s
More informationProblem set 2 for the course on. Markov chains and mixing times
J. Seif T. Hirscher Soluions o Proble se for he course on Markov chains and ixing ies February 7, 04 Exercise 7 (Reversible chains). (i) Assue ha we have a Markov chain wih ransiion arix P, such ha here
More informationChapter 15: Phenomena. Chapter 15 Chemical Kinetics. Reaction Rates. Reaction Rates R P. Reaction Rates. Rate Laws
Chaper 5: Phenomena Phenomena: The reacion (aq) + B(aq) C(aq) was sudied a wo differen emperaures (98 K and 35 K). For each emperaure he reacion was sared by puing differen concenraions of he 3 species
More informationMath Final Exam Solutions
Mah 246 - Final Exam Soluions Friday, July h, 204 () Find explici soluions and give he inerval of definiion o he following iniial value problems (a) ( + 2 )y + 2y = e, y(0) = 0 Soluion: In normal form,
More informationApplication of a Stochastic-Fuzzy Approach to Modeling Optimal Discrete Time Dynamical Systems by Using Large Scale Data Processing
Applicaion of a Sochasic-Fuzzy Approach o Modeling Opimal Discree Time Dynamical Sysems by Using Large Scale Daa Processing AA WALASZE-BABISZEWSA Deparmen of Compuer Engineering Opole Universiy of Technology
More informationSolutions to Odd Number Exercises in Chapter 6
1 Soluions o Odd Number Exercises in 6.1 R y eˆ 1.7151 y 6.3 From eˆ ( T K) ˆ R 1 1 SST SST SST (1 R ) 55.36(1.7911) we have, ˆ 6.414 T K ( ) 6.5 y ye ye y e 1 1 Consider he erms e and xe b b x e y e b
More informationIntroduction to AC Power, RMS RMS. ECE 2210 AC Power p1. Use RMS in power calculations. AC Power P =? DC Power P =. V I = R =. I 2 R. V p.
ECE MS I DC Power P I = Inroducion o AC Power, MS I AC Power P =? A Solp //9, // // correced p4 '4 v( ) = p cos( ω ) v( ) p( ) Couldn' we define an "effecive" volage ha would allow us o use he same relaionships
More informationMethodology. -ratios are biased and that the appropriate critical values have to be increased by an amount. that depends on the sample size.
Mehodology. Uni Roo Tess A ime series is inegraed when i has a mean revering propery and a finie variance. I is only emporarily ou of equilibrium and is called saionary in I(0). However a ime series ha
More informationThe equation to any straight line can be expressed in the form:
Sring Graphs Par 1 Answers 1 TI-Nspire Invesigaion Suden min Aims Deermine a series of equaions of sraigh lines o form a paern similar o ha formed by he cables on he Jerusalem Chords Bridge. Deermine he
More informationLecture 18 GMM:IV, Nonlinear Models
Lecure 8 :IV, Nonlinear Models Le Z, be an rx funcion of a kx paraeer vecor, r > k, and a rando vecor Z, such ha he r populaion oen condiions also called esiain equaions EZ, hold for all, where is he rue
More information(1) (2) Differentiation of (1) and then substitution of (3) leads to. Therefore, we will simply consider the second-order linear system given by (4)
Phase Plane Analysis of Linear Sysems Adaped from Applied Nonlinear Conrol by Sloine and Li The general form of a linear second-order sysem is a c b d From and b bc d a Differeniaion of and hen subsiuion
More informationLecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.
Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in
More information