Unknown Input High Gain Observer for Fault Detection and Isolation of Uncertain Systems

Size: px
Start display at page:

Download "Unknown Input High Gain Observer for Fault Detection and Isolation of Uncertain Systems"

Transcription

1 Engneerng Leer, 7:, EL_7 08 Unknown Inp Hgh Gan Observer or Fal Deecon and Isolaon o Unceran Sysems Sharddn Mondal*, G Chakrabory and K Bhaacharyya Absrac An nknown np hgh gan observer (UIHGO) based componen al deecon and solaon (FDI) echnqe s presened Frs, a redced order UIHGO s derved or a lnear sysem whose parameers are nceran o some exen he observer gan s deermned by solvng he well-known algebrac Rcca eqaon (ARE) hen, sng a bank o sch observers, a FDI algorhm s devsed o deec and solae he componen al (e, paramerc al) o an nceran sysem he FDI algorhm consss o wo seps In he rs sep, he deecon o al and he solaon o aly regon are accomplshed and n he nex sep, he aly parameer s solaed rom he aly regon Eecveness o he proposed observer as well as he FDI echnqe s shown wh he help o a nmercal example Index erms Unknown np hgh gan observer; componen al; al deecon and solaon; nceran sysem; parameer esmaon I INRODUCION Wh he rsng demands o hgh relably and saey o advanced processes lke avoncs, nclear power saons, aomobles ec have led o ncreasng reqremens o developng new mehods o spervson and monorng as a par o overall process conrol scheme Deren al deecon and solaon (FDI) schemes have been developed or avodng alre o he plans Model based al deecon echnqes (lke Kalman ler or observer based) have receved ncreasng aenon ollowng he poneerng work o Beard [] he FDI concep sng observers or Kalman lers s devsed based on he assmpon ha he mahemacal model o a sysem s perecly known In realy, however hs assmpon does no hold becase he parameers o a process are n general nceran or me varyng Agan he characerscs o dsrbances or nose are no compleely known; hence hey canno be perecly modeled here s always a msmach beween he acal process model and s mahemacal model (even here s no al n he process), whch somemes prodces alse alarms corrpng he Manscrp receved Aprl 5, 009 S Mondal s wh he Robocs and Bo-mecharoncs Laboraory, Deparmen o Mechancal Engneerng, Pohang Unversy o Scence and echnology (Posech), Pohang, , Soh Korea, el No: /84, Fax: , e-mal: sharm@posechackr G Chakrabory and K Bhaacharyya are wh he Deparmen o Mechancal Engneerng, Indan Inse o echnology, Kharagpr, Pn- 70, Inda, e-mals: goam@mechkgpernen, kng@mechkgpernen * Correspondng ahor perormance o he FDI echnqe o avod alse alarms, he FDI mehod shold be robs e, nsensve o modelng nceranes B he algorhm shold no be oo robs o gnore he al e, a sgncanly large varaon o he parameer vales Over he years, varos knds o robs al deecon and solaon echnqes have been developed o dagnose deren ypes o als lke sensors, acaors or componens [-4, 7, 9, 0,, 6-0] Frank [8], n a srvey paper, descrbed deren ypes o observer based robs al dagnoss echnqes Paon and Chen [5] dscssed varos robsness sses relaed observer based al dagnoss echnqes Lnear marx neqaly (LMI) based robs al deecon echnqes or nceran sysems have been developed n [9, 0] he dencaon based FDI echnqes have been sed by many researchers [,, ] o deec parameer als Daley and Wang [5] sed a hgh gan observer, whch was developed by Peersen and Hollo [4], as a ool or sensor al deecon In he presen work, an nknown np hgh gan observer (UIHGO) based componen (e, parameer) al deecon and solaon echnqe s derved Frs, an nknown np hgh gan observer s developed or a lnear nceran sysem Sch ype o nknown np observers has wde applcaons n modern conrol sysems he nceranes (modelng or paramerc or boh) are navodable Nex, sng a bank o sch observers, a parameer al deecon and solaon echnqe s devsed or a paramercally nceran sysem on he assmpons ha sensors and acaors are al ree Snce he hgh gan observer [5, 4] s robs agans parameer nceranes o some exen, he FDI echnqe s also robs agans he nceranes he FDI algorhm works n wo seps In sep-, he deecon o al and solaon o aly zone s accomplshed sng a bank o UIHGOs In he nex sep, aly parameer s solaed by parameer solaon mehod In he presen work, a par o he sysem parameers (e, he parameers o he aly sbsysem) s esmaed n sep- and only when a al occrs n he sysem In hs respec he complexy o al solaon s drascally redced n comparson wh sandard parameer dencaon echnqe [,, ] all he parameers o a sysem are esmaed a every me sep rrespecve o he occrrence o any al and he esmaed vales are compared wh her nomnal vales A nmercal example s presened o demonsrae he eecveness o he proposed observer as well as he FDI echnqe he basc mehodology o desgnng he nknown np hgh gan observer or an nceran sysem s dscssed n (Advance onlne pblcaon: May 009)

2 Engneerng Leer, 7:, EL_7 08 secon II he al deecon and solaon algorhm s explaned n secon III In secon IV, a nmercal example s presened o demonsrae he perormances o he proposed mehods he concldng remarks are nclded n secon V II UNKNOWN INPU HIGH GAIN OBSERVER In hs secon, an nknown np hgh gan observer s developed or a lnear nceran sysem he scen condons or he exsence o he observer are provded Consder a lnear me-nvaran sysem wh nknown nps x() = ( A+ A) x() + ( B+ B) () + Ed() () y() = Cx () () n x() R - he sae vecor, m () R - he p measrable np vecor, y() R - he op vecor and q d() R - he nknown np vecor he marces ABC,, and E o sable dmensons are known he marces A and B are he nceranes o he sysem and np marces respecvely hese may be consan or me varyng dependng on he sysem I s assmed ha ( A+ A) s always asympocally sable or all A I hs condon s no sased hen rs a conroller s o be desgned o sablze he sysem I s also assmed ha he sysem sases he rank condon: rank( CE) = rank( E) hs s a basc assmpon or desgnng any nknown np observer Now, sng a sae ransormaon marx, he saes are redened as z() = x () sch ha φ( n r) q E = E E s r q dmensonal marx wh rank( E) = rank( E ) and φ s a nll marx he sysem and op eqaons can be recas as ollows z A A z A A z B = + + z A A z A A z B () B φ( n r) q + + d B E z y() = C C (4) z Now s assmed ha he measremen sgnals are sch ha he ollowng rank condon s sased: rank( C) > rank( E) hs s a necessary condon or desgnng hs observer as he exra measremen sgnals are sed o desgn he redced order observer aer decoplng he nknown nps hs condon allows he rearrangemen o he op eqaon n he ollowng orm wh he help o a ransormaon y = Vy, V s a nonsnglar marx, y as y = C φ z y C C z (5) Now he eqaons () and (5) can be wren n expanded orm as ollows z = Az+ Az+ Az+ Az+ B + B (6) z = Az+ Az + Az+ Az + B+ B+ Ed (7) y C z (8) = y = Cz + Cz (9) Elmnang z rom he eqaon (6) sng he eqaon (9), one ges z = Az + AC ( y Cz) + B + Ed (0) Ed = Az + Az + B wh E - known marx and - nknown sgnal d I can be seen ha C shold be ll rank marx, whch wll be always so as rank( CE) = rank( E) Now, he eqaon (0) can be re-wren n smpled orm as z = Az + B+ Ed, () s s s = A A A C C, Bs = B AC and = y For desgnng an observer, he sysem shold sasy he observably condon: rank( O( AC, )) = n Now one can desgn an observer or he sysems () and (8) o esmae he sae z ˆ ˆ ˆ = s + s + ( () yˆ ˆ = Cz () he observer gan marx K s deermned by solvng he ollowng algebrac Rcca eqaon (ARE) [5, 4] q EE PEE P AP s + PAs + Q+ PC CP + = 0 (4) σ q σ wh K = PC, (5) Q s a pre-chosen posve dene marx and he consans q & σ are speced nmbers I was shown n [4] ha or any σ > 0, here exss q sch ha gan obaned rom he above eqaons wll lead o C( jwi As + KC) E < σ or w R w s he reqency hs condon mples ha he eec o nknown sgnal becomes very small n error dynamcs d or an approprae vale o σ Now he sae ẑ s esmaed rom eqaon (9) sng he esmaed sae ẑ as zˆ = C ( y C z ˆ ) (6) Fnally sng ˆ = { ˆ ˆ } ond o as complees z z z, he esmaed saes ˆx are xˆ = zˆ Wh hs, observer desgn process (Advance onlne pblcaon: May 009)

3 Engneerng Leer, 7:, EL_7 08 III FAUL DEECION AND ISOLAION ALGORIHM In hs secon, a componen al deecon and solaon echnqe or a lnear nceran sysem s derved I consss o wo seps In he rs sep, a se o resdals s generaed wh he help o a bank o nknown np hgh gan observers (UIHGOs) o deec he al and solae he aly regon In he second sep, aly parameer s solaed rom he aly regon Consder a lnear me nvaran sysem as x() = ( A+ A) x() + ( B+ B) () (7) he sgncance o he marces and vecors are same as descrbed n secon II Sppose a al occrs n a componen o he plan he deecon and solaon o he al are carred o n wo seps as ollows Sep-: Deecon and paral solaon o al he aly sysem s wren as x() = ( A+ A + A) x() + ( B+ B + B) (), (8) A and B are he aly pars o he marces A and B respecvely I can be emphaszed ha he magnde o als (e, A and B ) shold be sgncanly larger compared o he magnde o nceranes (e, A and B ) he sae eqaon (8) can now be rearranged as x() = ( A+ A) x() + ( B+ B) () + Ed(), (9) q E s a known marx and d() R s he nknown np sasyng he relaon Ed() = A x() + B () (0) Now he sysem s dvded no N nmbers o sbsysems wh each characerzed by a ew parameers he choce o sbsysems s qe arbrary b here shold no have any common elemens beween hem In a physcal sysem, he sbsysems are chosen based on he physcal proxmy o deren parameers Assme ha he al has been occrred n he -h sbsysem he sysem eqaon consderng he al n he -h sbsysem s wren as x () = ( A+ A) x () + ( B+ B) () + E d () () () () () () he sbscrp () ndcaes ha he al has been consdered n he -h sbsysem he op eqaon or hs sysem s wren as y () = C x (), () () () () he eqaons () and () are smlar o he eqaons () and () Now, ollowng he procedre dscssed n secon II, an nknown np hgh gan observer s desgned o esmae he saes xˆ () Once he saes () () xˆ () calclaed as r () = y () yˆ () = y () C x ˆ () () () () () () () () are esmaed, he resdals are Now, an nknown np observer, properly desgned, can esmae he saes rrespecve o nknown nps So he resdal r, calclaed rom he eqaon (), () () converges whn a bonded vale known as hreshold vale he al occrs n he -h sbsysem or here s no al n he sysem as he eec o possble als n -h sbsysem s consdered as nknown nps In deal case, e, n he absence o nose and parameer nceranes, he resdals shold converge o zero (hogh a small hreshold vale s always se o ake care o errors de o he nmercal lmaons) as n he presen case he convergence akes place whn a hreshold vale, whch agan depends on he amon o nceranes and np sgnal appled o he sysem In hs way, one can deec a al and solae he aly sbsysem sng N nmbers o resdals calclaed wh he help o ha nmbers o UIHGOs However ( N ) sch observers wll be scen o solae a aly sbsysem when N > becase once ( N ) sbsysems are ond al ree, he remanng sbsysem s aomacally dened as he aly one A decson able (as shown n able ) s drawn o solae he aly sbsysem rom observaon o ( N ) resdals able : Decson able or solaon o aly sbsysem Observa on cases For gven resdals ( r() ) and hreshold vales ( ε ): r > ε () () () () () hen se and r ε hen se 0 r r r () () () r( N ) Decson Case No al Case- 0 Fal: SS Case- 0 Fal: SS Case- 0 Fal: SS Case- (N-) 0 Fal: SS (N-) Case-N Fal: SS (N) Sep-: oal solaon o al Once he aly sbsysem s solaed, he aly parameer n he aly sbsysem s dened n hs sep Frs, he eec o he aly sbsysem s smlaed as an nknown np sgnal, say F he relaonshp beween, he () F () parameers o he aly sbsysem, say s, and he saes x() are known and can be wren as F () = ( s, x ()), (4) he ncon s lnear or a lnear sysem he sysem eqaon or hs case becomes x() = ( A+ A) x() + ( B+ B) () + Ed (), (5) (Advance onlne pblcaon: May 009)

4 Engneerng Leer, 7:, EL_7 08 d() = F Wh a measremen marx, an () C observer s hen desgned o esmae he saes Knowng he saes, he nknown np sgnal s esmaed rom he sae eqaons neglecng he nceranes sng he nomnal vales o he parameers o he oher non-aly sbsysems he esmaed sgnal Fˆ () s now sed o esmae he parameers s 's rom he relaon (4), whch s rewren as Fˆ ˆ () = (, s x ()) (6) All he elemens o eqaon (6) excepng he parameers s are known Deren parameer esmaon echnqes can be sed o esmae s rom eqaon (6) However, a very smple logcal approach s appled n he presen work n order o solae he aly elemen Le s consder he k-h parameer as he aly one From he above relaon, s k can be esmaed sng nomnal vales o res o he parameers Mahemacally, sˆ ˆ ˆ k = gs (, s,, sk, sk+, sl, x ( ), F( )) (7) l s he nmber o parameers o he aly sbsysem and g s a nconal I s observed ha n seady sae, he esmaed vales lcae very less he assmpon s correc he movng averages echnqe can be sed o smoohen he lcaon o he esmaed vales de o nceranes I he assmpon s wrong, he esmaed vales vary sgncanly large Now, as he sngle al case s beng consdered, here wll be only one case when he esmaed parameer wll vary less he parclar parameer or whch happens s he aly one In hs way, he aly parameer s solaed In he same way, any paramerc al o any s k sbsysem can be deeced and solaed ollowng he above wo seps he FDI echnqe can be smmarzed n a block dagram as shown gre IV NUMERICAL EXAMPLE Consder a mechancal sysem (as shown n gre ) ha consss o wo mass elemens and hree ses o sprngs and dampers he sae space model o he sysem can be wren as ollows: x() = ( A+ A) x() + ( B+ B) (), , ( K+ K) K ( C+ C) C A= M M M M K ( K + K) C ( C + C) M M M M 000(/ ) = M, x = X X X X, () = F(), B [ ] { } X and mass elemen X are he dsplacemen and velocy o he M respecvely, K - he sness elemen j and C - he dampng coecen (=, and j=,, ) he j marces A and B are he nceranes o he sysem and np marces respecvely Fgre : Mechancal sysem havng wo masses and hree ses o sprng-damper Fgre : Srcre o he FDI algorhm he nmercal vales o he sysem parameers are M = 870 kg, M = 550 kg, K, = N / m K = N / m, K = N / m, C = 500 Ns / m, C = 000 Ns / m and C = 5675 Ns / m he sysem response wh arbrary nal condons s smlaed sng MALAB-SIMULINK oolbox he paramerc nceranes are smlaed n sch a way ha he elemens o he marces A and B der maxmm o ± 5 % rom her nomnal vales A al s now nrodced n he sprng o sness K a he new vale o = 50sec K s se o N/m Now sng he FDI algorhm, (Advance onlne pblcaon: May 009)

5 Engneerng Leer, 7:, EL_7 08 dscssed n secon III, he al s deeced and he aly elemen (here K ) s solaed as ollows Sep-: Deecon and paral solaon o he al Frs, he sysem s dvded no hree sbsysems as ollows: SS : K, & C M ; : SS K & and C SS : K, C & M he nceranes are nrodced as ollows: A= MΣ N and B = MΣ N wh, M = In N = 005 A, N and = 005 B Σ = Σ0sn( w ) Σ 0 = 05 I and w rad/s he marces and = 005 A B are same as A and B excepng he elemens conanng consan erms are replaced wh zeros he snsodal varaon n sysem parameers s nrodced n smlaon he ollowng np sgnal s appled or hs case: = sn( w) wh 0 0 = 00 N and w = rad/s As he sysem s dvded no hree sbsysems, so wo UIHGOs are scen as a par o sep- he observers are desgned or SS and SS he nknown np marces E s and nknown np sgnals d s or hose observers are gven below E [ ], () = d() = A() x() + B() () E () = [ ], d () = A() x() + B() () he op marces are C () = and C() = 0 0 Now applyng he mehod dscssed n secon-ii, wo E s are ond o as E () [ 0 0 ] = and E () [ 0 0 ] = he vales o he nng parameers σ and q are consdered as σ () = 005, σ () = 005, q () = 5 and q he vale o Q s chosen as or boh () = 5 Q = 5I n he observers he gan marces are calclaed by solvng he eqaons (4) and (5) he vales o he observer gans or he above observers are K () = [ ] and K () = [ ] respecvely wo hgh gan observers are hen desgned or he above sysems Fnally he resdals are calclaed and ploed n gre and gre 4 In deal saon (e, al ree and n he absence o parameer nceranes), he resdals shold be zero However n he presen case hese wll no be zero de o presence o parameer nceranes Hence wo small 5 0 ε() = and consan hreshold vales { } 9 { } 0 ns are chosen as he smlaon ε() = 0 0 s carred o applyng xed np sgnals In real saon, adapve hreshold vales [5] shold be chosen as np sgnals vary dependng on operang condons Here he hreshold vales are calclaed n normal operang condon e, when here s no al n he sysem Fgre : Componens o he resdal r () () Fgre 4: Componens o he resdal r () () (Advance onlne pblcaon: May 009)

6 Engneerng Leer, 7:, EL_7 08 As boh he resdals cross he hreshold vales, he exsence o a al s conrmed o solae he aly sbsysem a decson able (able ) s consrced as shown below able : Decson able or solaon o aly sbsysem Observaon Is r () > () r () () ε (se ) or ε (se 0 )? r r () () Decson Case Fal: SS From able-, s seen ha he al s n sbsysem So he nex sep (e, sep-) s carred o o solae he aly parameer Sep-: oal solaon o he al Here he aly sbsysem (SS) s rs replaced wh an nknown orce F () as F () = K( X X) + C( X X ) hen he sysem s remodeled as ollows x() = ( A+ A) x() + ( B+ B) () + Ed() E = 0 0 / M / M and d() = F () I can [ ] be noced ha he parameer nceranes or hs sysem are only n rd and 4 h rows o A and B marces, whch ndcae A and B have non-zero elemens n rd and 4 h rows only Agan he marx E conans non-zero elemens n he same rows For hs smlary here he sysem eqaon s remodeled combnng he nknown nps and nceranes as x() = Ax() + B() + E d () and A, 4x+ B,+ F / M Ec = dc = 0 A4, 4x+ B4, F / M 0 hs s a specal case, whch may no appear or all sysems Now a ll order nknown np observer [6] s desgned 0 0 wh op marx C = o esmae he saes xˆ( ) Usng he esmaed saes xˆ( ) and he nomnal vales o he parameers o sbsysem-, he nknown orce F () s esmaed rom he ollowng relaonshp F ˆ () = M x ˆ + K x ˆ + C x ˆ, ˆ x s calclaed akng he dervave o wh ˆx respec o me Fnally he aly parameers are esmaed sng he ollowng relaon Fˆ () ˆ ˆ ˆ ˆ = K ( x ) ( ) x + C x4 x Frs, he al s assmed o resde n he sness elemen and sng he nomnal vale o C =000 Ns/m, K c c K s esmaed he movng averages are aken o redce he eec o nceranes and nmercal errors o he esmaed vales he nal movng average s aken wh a me wndow o 5 (5-0 seconds o daa) sec Inal daa are no aken o redce he ranson eec ha comes de o nal condons Now wh an ncremen o sec, he movng averages are calclaed po 0 sec o me span and he esmaed vales are ploed n gre 5 he plo shows ha esmaed vales vary very less (maxmm varaon o 5 % rom s mean vale) Now s assmed o be aly and sng he nomnal vale o K =70000 N/m, s esmaed he esmaed vales C aer akng movng averages n smlar manner as n case o are ploed n gre 6 he plo shows ha he esmaed K vales o C C vary wdely (as hgh as 75 % rom s mean vale), whch s becase o he wrong assmpon hs Fgre 5: Esmaed sness Fgre 6: Esmaed dampng coecen (Advance onlne pblcaon: May 009)

7 Engneerng Leer, 7:, EL_7 08 conrms ha aly elemen s K and hereby he al solaon process complees hs s seen ha he FDI scheme works well or he occrrence o a al n sbsysem I can be shown easly ha he mehod works wh eqal ease or he occrrence o any parameer al n any oher sbsysem V CONCLUSIONS An nknown np hgh gan observer (UIHGO) based componen al deecon and solaon (FDI) scheme s presened Frs an UIHGO or a lnear nceran sysem s derved hese ypes o observers have wde applcaons n robs conrol and al dagnoss hen, sng a bank o sch observers, a FDI echnqe s devsed he advanage o he FDI algorhm s ha s capable o esmang als even he parameers are copled n he sysem marx I also redces he complexy o esmang all he parameers a every me nsan nlke exsng dencaon based parameer al dagnoss echnqes he same FDI echnqe can also be sed o deec a al n a nosy sysem or a nonlnear sysem provded oher ypes o nknown np esmaors capable o handlng nose or nonlneary shold be sed Proc 44h IEEE Con on Decson and Conrol, and he Eropean Conrol Con 005, Sevlle, Span, pp 88-84, 005 [4] Peersen, I R and Hollo, C V, Hgh-gan observer approach o dsrbance aenaon sng measremen eedback, In J o Conrol, Vol 48, No 6, pp , 988 [5] Paon, R J and Chen, J, Observer-based al deecon and solaon: robsness and applcaons, Conrol Engneerng Pracce, Vol 5, No 5, pp 67-68, 997 [6] Pg, V, Qevdo, J, Escobe, and Sanc, A, Robs Fal Deecon sng Lnear Inerval Observers, In Proc 5h IFAC Symposm on Fal Deecon, Spervson and Saey o echncal Processes, Washngon DC, pp , 00 [7] Sa M and Gan, Y, A New Approach o Robs Fal Deecon and Idencaon, IEEE rans Aerospace and Elecronc Sysems, Vol 9, No, pp , 99 [8] Shen, L C and Hs, P L, Robs Desgn o Fal Isolaon Observers, Aomaca, Vol 4, No, pp 4-49, 998 [9] Wang, H B, Wang, J L and Lam, J, Robs Fal Deecon Observer Desgn: Ierave LMI Approaches, J o Dynamc Sysems, Measremen, and Conrol, Vol 9, pp 77-8, 007 [0] Zhong, M, Dng, S X, Lam, J and Wang, H, An LMI approach o desgn robs al deecon ler or nceran LI sysems, Aomaca, Vol 9, No, pp , 00 REFERENCES [] Beard, R V, Falre Accommodaon n Lnear Sysems hrogh Sel-Reorganzaon, Ph D dsseraon, Deparmen o Aeronacs and Asronacs, MI, Cambrdge, Massachses, 97 [] Bloch, G, Oladsne, M and homas, P, On-lne al dagnoss o dynamc sysems va robs parameer esmaon, Conrol Eng Pracce, Vol, No, pp , 995 [] Chen, J, Paon, R J and Zhang, H-Y, Desgn o nknown np observers and robs al deecon lers, In J o Conrol, Vol 6, No, pp 85-05, 996 [4] Chow, E Y and Wllsky, A S, Analycal Redndancy and he desgn o robs alre deecon, IEEE rans Aomac Conrol, Vol 9, No 7, pp 60-64, 984 [5] Daley, S and Wang, H, Applcaon o a Hgh Gan Observer o Fal Deecon, In Proc nd IEEE Con on Conrol Applcaons, Vancover, B C, Sepember -6, pp 6-6, 99 [6] Daroach, M, Zasadznsk, M and X, S J, Fll order observers or lnear sysems wh nknown nps, IEEE rans Aomac Conrol, Vol 9, No, pp , 994 [7] Dan, G R and Paon, R J, Robs Fal Deecon n Lnear Sysems sng Lenberger Observers, In Proc UKACC In Con on CONROL 98, Wales Swansea, UK, pp , 998 [8] Frank, P M, Enhancemen o robsness n observer-based al deecon, In J o Conrol, Vol 59, No 4, pp , 994 [9] Ge, W and Fang, C Z, Deecon o aly componens va robs observaon, In J o Conrol, Vol 47, No, pp , 988 [0] Jang, B, Wang, J L and Soh, Y C, Robs Fal Dagnoss or a Class o Lnear Sysems wh Uncerany, In Proc he Amercan Conrol Con, San Dego, Calorna USA, pp , 999 [] Jang,, Khorasan, K and aazol, S, Parameer Esmaon- Based Fal Deecon, Isolaon and Recovery or Nonlnear Saelle Models, IEEE rans Conrol Sysems echnology, Vol 6, No 4, pp , 008 [] Mondal, S, Chakrabory, G and Bhaacharyya, K, An nknown np Kalman ler based componen FDI algorhm and s applcaon n aomoble, In J o Vehcle Aonomos Sysems, Vol 5, No /4, pp 74-87, 007 [] Palma, L B, Coo, F V and da Slva, R N, Dagnoss o Paramerc Fals Based on dencaon and Sascal Mehods, In (Advance onlne pblcaon: May 009)

Unknown Input High Gain Observer for Parametric Fault Detection and Isolation of Dynamical Systems

Unknown Input High Gain Observer for Parametric Fault Detection and Isolation of Dynamical Systems Proceedings o the International MltiConerence o Engineers Compter Scientists 009 Vol II IMECS 009, March 8-0, 009, Hong Kong Unknown Inpt High Gain Observer or Parametric Falt Detection Isolation o Dynamical

More information

Relative controllability of nonlinear systems with delays in control

Relative controllability of nonlinear systems with delays in control Relave conrollably o nonlnear sysems wh delays n conrol Jerzy Klamka Insue o Conrol Engneerng, Slesan Techncal Unversy, 44- Glwce, Poland. phone/ax : 48 32 37227, {jklamka}@a.polsl.glwce.pl Keywor: Conrollably.

More information

CONSISTENT EARTHQUAKE ACCELERATION AND DISPLACEMENT RECORDS

CONSISTENT EARTHQUAKE ACCELERATION AND DISPLACEMENT RECORDS APPENDX J CONSSTENT EARTHQUAKE ACCEERATON AND DSPACEMENT RECORDS Earhqake Acceleraons can be Measred. However, Srcres are Sbjeced o Earhqake Dsplacemens J. NTRODUCTON { XE "Acceleraon Records" }A he presen

More information

Different kind of oscillation

Different kind of oscillation PhO 98 Theorecal Qeson.Elecrcy Problem (8 pons) Deren knd o oscllaon e s consder he elecrc crc n he gre, or whch mh, mh, nf, nf and kω. The swch K beng closed he crc s copled wh a sorce o alernang crren.

More information

Observer Design for Nonlinear Systems using Linear Approximations

Observer Design for Nonlinear Systems using Linear Approximations Observer Desgn for Nonlnear Ssems sng Lnear Appromaons C. Navarro Hernandez, S.P. Banks and M. Aldeen Deparmen of Aomac Conrol and Ssems Engneerng, Unvers of Sheffeld, Mappn Sree, Sheffeld S 3JD. e-mal:

More information

Dynamic Model of the Axially Moving Viscoelastic Belt System with Tensioner Pulley Yanqi Liu1, a, Hongyu Wang2, b, Dongxing Cao3, c, Xiaoling Gai1, d

Dynamic Model of the Axially Moving Viscoelastic Belt System with Tensioner Pulley Yanqi Liu1, a, Hongyu Wang2, b, Dongxing Cao3, c, Xiaoling Gai1, d Inernaonal Indsral Informacs and Comper Engneerng Conference (IIICEC 5) Dynamc Model of he Aally Movng Vscoelasc Bel Sysem wh Tensoner Plley Yanq L, a, Hongy Wang, b, Dongng Cao, c, Xaolng Ga, d Bejng

More information

Including the ordinary differential of distance with time as velocity makes a system of ordinary differential equations.

Including the ordinary differential of distance with time as velocity makes a system of ordinary differential equations. Soluons o Ordnary Derenal Equaons An ordnary derenal equaon has only one ndependen varable. A sysem o ordnary derenal equaons consss o several derenal equaons each wh he same ndependen varable. An eample

More information

Bandlimited channel. Intersymbol interference (ISI) This non-ideal communication channel is also called dispersive channel

Bandlimited channel. Intersymbol interference (ISI) This non-ideal communication channel is also called dispersive channel Inersymol nererence ISI ISI s a sgnal-dependen orm o nererence ha arses ecause o devaons n he requency response o a channel rom he deal channel. Example: Bandlmed channel Tme Doman Bandlmed channel Frequency

More information

DEEP UNFOLDING FOR MULTICHANNEL SOURCE SEPARATION SUPPLEMENTARY MATERIAL

DEEP UNFOLDING FOR MULTICHANNEL SOURCE SEPARATION SUPPLEMENTARY MATERIAL DEEP UNFOLDING FOR MULTICHANNEL SOURCE SEPARATION SUPPLEMENTARY MATERIAL Sco Wsdom, John Hershey 2, Jonahan Le Roux 2, and Shnj Waanabe 2 Deparmen o Elecrcal Engneerng, Unversy o Washngon, Seale, WA, USA

More information

CH.3. COMPATIBILITY EQUATIONS. Continuum Mechanics Course (MMC) - ETSECCPB - UPC

CH.3. COMPATIBILITY EQUATIONS. Continuum Mechanics Course (MMC) - ETSECCPB - UPC CH.3. COMPATIBILITY EQUATIONS Connuum Mechancs Course (MMC) - ETSECCPB - UPC Overvew Compably Condons Compably Equaons of a Poenal Vecor Feld Compably Condons for Infnesmal Srans Inegraon of he Infnesmal

More information

Ordinary Differential Equations in Neuroscience with Matlab examples. Aim 1- Gain understanding of how to set up and solve ODE s

Ordinary Differential Equations in Neuroscience with Matlab examples. Aim 1- Gain understanding of how to set up and solve ODE s Ordnary Dfferenal Equaons n Neuroscence wh Malab eamples. Am - Gan undersandng of how o se up and solve ODE s Am Undersand how o se up an solve a smple eample of he Hebb rule n D Our goal a end of class

More information

by Lauren DeDieu Advisor: George Chen

by Lauren DeDieu Advisor: George Chen b Laren DeDe Advsor: George Chen Are one of he mos powerfl mehods o nmercall solve me dependen paral dfferenal eqaons PDE wh some knd of snglar shock waves & blow-p problems. Fed nmber of mesh pons Moves

More information

Chapter 6 DETECTION AND ESTIMATION: Model of digital communication system. Fundamental issues in digital communications are

Chapter 6 DETECTION AND ESTIMATION: Model of digital communication system. Fundamental issues in digital communications are Chaper 6 DCIO AD IMAIO: Fndaenal sses n dgal concaons are. Deecon and. saon Deecon heory: I deals wh he desgn and evalaon of decson ang processor ha observes he receved sgnal and gesses whch parclar sybol

More information

HEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD

HEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD Journal of Appled Mahemacs and Compuaonal Mechancs 3, (), 45-5 HEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD Sansław Kukla, Urszula Sedlecka Insue of Mahemacs,

More information

, t 1. Transitions - this one was easy, but in general the hardest part is choosing the which variables are state and control variables

, t 1. Transitions - this one was easy, but in general the hardest part is choosing the which variables are state and control variables Opmal Conrol Why Use I - verss calcls of varaons, opmal conrol More generaly More convenen wh consrans (e.g., can p consrans on he dervaves More nsghs no problem (a leas more apparen han hrogh calcls of

More information

OPTIMIZATION OF A NONCONVENTIONAL ENGINE EVAPORATOR

OPTIMIZATION OF A NONCONVENTIONAL ENGINE EVAPORATOR Jornal of KONES Powerran and Transpor, Vol. 17, No. 010 OPTIMIZATION OF A NONCONVENTIONAL ENGINE EVAPORATOR Andre Kovalí, Eml Toporcer Unversy of Žlna, Facly of Mechancal Engneerng Deparmen of Aomove Technology

More information

Chapters 2 Kinematics. Position, Distance, Displacement

Chapters 2 Kinematics. Position, Distance, Displacement Chapers Knemacs Poson, Dsance, Dsplacemen Mechancs: Knemacs and Dynamcs. Knemacs deals wh moon, bu s no concerned wh he cause o moon. Dynamcs deals wh he relaonshp beween orce and moon. The word dsplacemen

More information

On One Analytic Method of. Constructing Program Controls

On One Analytic Method of. Constructing Program Controls Appled Mahemacal Scences, Vol. 9, 05, no. 8, 409-407 HIKARI Ld, www.m-hkar.com hp://dx.do.org/0.988/ams.05.54349 On One Analyc Mehod of Consrucng Program Conrols A. N. Kvko, S. V. Chsyakov and Yu. E. Balyna

More information

Variants of Pegasos. December 11, 2009

Variants of Pegasos. December 11, 2009 Inroducon Varans of Pegasos SooWoong Ryu bshboy@sanford.edu December, 009 Youngsoo Cho yc344@sanford.edu Developng a new SVM algorhm s ongong research opc. Among many exng SVM algorhms, we wll focus on

More information

Hierarchical Sliding Mode Control for Series Double Inverted Pendulums System

Hierarchical Sliding Mode Control for Series Double Inverted Pendulums System Herarchcal Sldng Mode Conrol for Seres Doble Invered Pendlms Sysem Danwe Qan, Janqang Y, Dongbn Zhao, and Ynxng Hao Laboraory of Complex Sysems and Inellgence Scence Inse of Aomaon, Chnese Academy of Scences

More information

Approximate Analytic Solution of (2+1) - Dimensional Zakharov-Kuznetsov(Zk) Equations Using Homotopy

Approximate Analytic Solution of (2+1) - Dimensional Zakharov-Kuznetsov(Zk) Equations Using Homotopy Arcle Inernaonal Journal of Modern Mahemacal Scences, 4, (): - Inernaonal Journal of Modern Mahemacal Scences Journal homepage: www.modernscenfcpress.com/journals/jmms.aspx ISSN: 66-86X Florda, USA Approxmae

More information

ABSTRACT. Keywords: Finite Element, Active Optics, Adaptive Optics, WaveFront Error, System Analysis, Optomechanical 1.

ABSTRACT. Keywords: Finite Element, Active Optics, Adaptive Optics, WaveFront Error, System Analysis, Optomechanical 1. Analyss echne or conrollng Sysem Waveron Error wh Acve/Adapve Opcs Vcor L. Genberg*, Gregory J. Mchels Sgmadyne, 83 Wes Ave, ocheser, NY 14611 *genberg@sgmadyne.com (585)35-746 ABSTACT The lmae goal o

More information

V.Abramov - FURTHER ANALYSIS OF CONFIDENCE INTERVALS FOR LARGE CLIENT/SERVER COMPUTER NETWORKS

V.Abramov - FURTHER ANALYSIS OF CONFIDENCE INTERVALS FOR LARGE CLIENT/SERVER COMPUTER NETWORKS R&RATA # Vol.) 8, March FURTHER AALYSIS OF COFIDECE ITERVALS FOR LARGE CLIET/SERVER COMPUTER ETWORKS Vyacheslav Abramov School of Mahemacal Scences, Monash Unversy, Buldng 8, Level 4, Clayon Campus, Wellngon

More information

GENERATING CERTAIN QUINTIC IRREDUCIBLE POLYNOMIALS OVER FINITE FIELDS. Youngwoo Ahn and Kitae Kim

GENERATING CERTAIN QUINTIC IRREDUCIBLE POLYNOMIALS OVER FINITE FIELDS. Youngwoo Ahn and Kitae Kim Korean J. Mah. 19 (2011), No. 3, pp. 263 272 GENERATING CERTAIN QUINTIC IRREDUCIBLE POLYNOMIALS OVER FINITE FIELDS Youngwoo Ahn and Kae Km Absrac. In he paper [1], an explc correspondence beween ceran

More information

Comb Filters. Comb Filters

Comb Filters. Comb Filters The smple flers dscussed so far are characered eher by a sngle passband and/or a sngle sopband There are applcaons where flers wh mulple passbands and sopbands are requred Thecomb fler s an example of

More information

Real Time Hybrid Simulation using Shaking Tabels

Real Time Hybrid Simulation using Shaking Tabels Real Tme Hybrd Smlaon sng Shakng Tabels Deparmen o Cvl and Envronmenal Engneerng Unversy o Kassel, Germany Olne Inrodcon A ndamenal sbsrcre algorhm wh sb seppng Applcaons o he algorhm Conclsons Inrodcon

More information

In the complete model, these slopes are ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL. (! i+1 -! i ) + [(!") i+1,q - [(!

In the complete model, these slopes are ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL. (! i+1 -! i ) + [(!) i+1,q - [(! ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL The frs hng o es n wo-way ANOVA: Is here neracon? "No neracon" means: The man effecs model would f. Ths n urn means: In he neracon plo (wh A on he horzonal

More information

Online Supplement for Dynamic Multi-Technology. Production-Inventory Problem with Emissions Trading

Online Supplement for Dynamic Multi-Technology. Production-Inventory Problem with Emissions Trading Onlne Supplemen for Dynamc Mul-Technology Producon-Invenory Problem wh Emssons Tradng by We Zhang Zhongsheng Hua Yu Xa and Baofeng Huo Proof of Lemma For any ( qr ) Θ s easy o verfy ha he lnear programmng

More information

Solution in semi infinite diffusion couples (error function analysis)

Solution in semi infinite diffusion couples (error function analysis) Soluon n sem nfne dffuson couples (error funcon analyss) Le us consder now he sem nfne dffuson couple of wo blocks wh concenraon of and I means ha, n a A- bnary sysem, s bondng beween wo blocks made of

More information

Intelligent Motion Control for Electromechanical Servos Using Evolutionary Learning and Adaptation Mechanisms

Intelligent Motion Control for Electromechanical Servos Using Evolutionary Learning and Adaptation Mechanisms Inellgen Moon Conrol or Elecromechancal ervos Usng Evolonary earnng and Adapaon Mechansms ergey Edward yshevsk 1 and Mchael G aonov 1 Deparmen o Elecrcal and Comper Engneerng Prde Unversy a Indanapols,

More information

Real-Time Trajectory Generation and Tracking for Cooperative Control Systems

Real-Time Trajectory Generation and Tracking for Cooperative Control Systems Real-Tme Trajecor Generaon and Trackng for Cooperave Conrol Ssems Rchard Mrra Jason Hcke Calforna Inse of Technolog MURI Kckoff Meeng 14 Ma 2001 Olne I. Revew of prevos work n rajecor generaon and rackng

More information

PHYS 1443 Section 001 Lecture #4

PHYS 1443 Section 001 Lecture #4 PHYS 1443 Secon 001 Lecure #4 Monda, June 5, 006 Moon n Two Dmensons Moon under consan acceleraon Projecle Moon Mamum ranges and heghs Reerence Frames and relae moon Newon s Laws o Moon Force Newon s Law

More information

Cubic Bezier Homotopy Function for Solving Exponential Equations

Cubic Bezier Homotopy Function for Solving Exponential Equations Penerb Journal of Advanced Research n Compung and Applcaons ISSN (onlne: 46-97 Vol. 4, No.. Pages -8, 6 omoopy Funcon for Solvng Eponenal Equaons S. S. Raml *,,. Mohamad Nor,a, N. S. Saharzan,b and M.

More information

John Geweke a and Gianni Amisano b a Departments of Economics and Statistics, University of Iowa, USA b European Central Bank, Frankfurt, Germany

John Geweke a and Gianni Amisano b a Departments of Economics and Statistics, University of Iowa, USA b European Central Bank, Frankfurt, Germany Herarchcal Markov Normal Mxure models wh Applcaons o Fnancal Asse Reurns Appendx: Proofs of Theorems and Condonal Poseror Dsrbuons John Geweke a and Gann Amsano b a Deparmens of Economcs and Sascs, Unversy

More information

Time-interval analysis of β decay. V. Horvat and J. C. Hardy

Time-interval analysis of β decay. V. Horvat and J. C. Hardy Tme-nerval analyss of β decay V. Horva and J. C. Hardy Work on he even analyss of β decay [1] connued and resuled n he developmen of a novel mehod of bea-decay me-nerval analyss ha produces hghly accurae

More information

CS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 4

CS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 4 CS434a/54a: Paern Recognon Prof. Olga Veksler Lecure 4 Oulne Normal Random Varable Properes Dscrmnan funcons Why Normal Random Varables? Analycally racable Works well when observaon comes form a corruped

More information

12d Model. Civil and Surveying Software. Drainage Analysis Module Detention/Retention Basins. Owen Thornton BE (Mech), 12d Model Programmer

12d Model. Civil and Surveying Software. Drainage Analysis Module Detention/Retention Basins. Owen Thornton BE (Mech), 12d Model Programmer d Model Cvl and Surveyng Soware Dranage Analyss Module Deenon/Reenon Basns Owen Thornon BE (Mech), d Model Programmer owen.hornon@d.com 4 January 007 Revsed: 04 Aprl 007 9 February 008 (8Cp) Ths documen

More information

Tools for Analysis of Accelerated Life and Degradation Test Data

Tools for Analysis of Accelerated Life and Degradation Test Data Acceleraed Sress Tesng and Relably Tools for Analyss of Acceleraed Lfe and Degradaon Tes Daa Presened by: Reuel Smh Unversy of Maryland College Park smhrc@umd.edu Sepember-5-6 Sepember 28-30 206, Pensacola

More information

( t) Outline of program: BGC1: Survival and event history analysis Oslo, March-May Recapitulation. The additive regression model

( t) Outline of program: BGC1: Survival and event history analysis Oslo, March-May Recapitulation. The additive regression model BGC1: Survval and even hsory analyss Oslo, March-May 212 Monday May 7h and Tuesday May 8h The addve regresson model Ørnulf Borgan Deparmen of Mahemacs Unversy of Oslo Oulne of program: Recapulaon Counng

More information

Existence of Periodic Solution for a Non-Autonomous Stage-Structured Predator-Prey System with Impulsive Effects

Existence of Periodic Solution for a Non-Autonomous Stage-Structured Predator-Prey System with Impulsive Effects Appled ahemacs 55-6 do:.6/am.. Pblshed Onlne arch (hp://www.scrp.org/jornal/am) Exsence o Perodc Solon or a Non-Aonomos Sage-Srcred Predaor-Prey Sysem wh Implsve Eecs Absrac eng W Zolang Xong Ypng Deng

More information

Decentralised Sliding Mode Load Frequency Control for an Interconnected Power System with Uncertainties and Nonlinearities

Decentralised Sliding Mode Load Frequency Control for an Interconnected Power System with Uncertainties and Nonlinearities Inernaonal Research Journal of Engneerng and echnology IRJE e-iss: 2395-0056 Volume: 03 Issue: 12 Dec -2016 www.re.ne p-iss: 2395-0072 Decenralsed Sldng Mode Load Frequency Conrol for an Inerconneced Power

More information

XIII International PhD Workshop OWD 2011, October Three Phase DC/DC Boost Converter With High Energy Efficiency

XIII International PhD Workshop OWD 2011, October Three Phase DC/DC Boost Converter With High Energy Efficiency X nernaonal Ph Workshop OW, Ocober Three Phase C/C Boos Converer Wh Hgh Energy Effcency Ján Perdľak, Techncal nversy of Košce Absrac Ths paper presens a novel opology of mlphase boos converer wh hgh energy

More information

Dynamic Team Decision Theory. EECS 558 Project Shrutivandana Sharma and David Shuman December 10, 2005

Dynamic Team Decision Theory. EECS 558 Project Shrutivandana Sharma and David Shuman December 10, 2005 Dynamc Team Decson Theory EECS 558 Proec Shruvandana Sharma and Davd Shuman December 0, 005 Oulne Inroducon o Team Decson Theory Decomposon of he Dynamc Team Decson Problem Equvalence of Sac and Dynamc

More information

Stochastic Programming handling CVAR in objective and constraint

Stochastic Programming handling CVAR in objective and constraint Sochasc Programmng handlng CVAR n obecve and consran Leondas Sakalaskas VU Inse of Mahemacs and Informacs Lhana ICSP XIII Jly 8-2 23 Bergamo Ialy Olne Inrodcon Lagrangan & KKT condons Mone-Carlo samplng

More information

2/20/2013. EE 101 Midterm 2 Review

2/20/2013. EE 101 Midterm 2 Review //3 EE Mderm eew //3 Volage-mplfer Model The npu ressance s he equalen ressance see when lookng no he npu ermnals of he amplfer. o s he oupu ressance. I causes he oupu olage o decrease as he load ressance

More information

Research Article Cubic B-spline for the Numerical Solution of Parabolic Integro-differential Equation with a Weakly Singular Kernel

Research Article Cubic B-spline for the Numerical Solution of Parabolic Integro-differential Equation with a Weakly Singular Kernel Researc Jornal of Appled Scences, Engneerng and Tecnology 7(): 65-7, 4 DOI:.96/afs.7.5 ISS: 4-7459; e-iss: 4-7467 4 Mawell Scenfc Pblcaon Corp. Sbmed: Jne 8, Acceped: Jly 9, Pblsed: Marc 5, 4 Researc Arcle

More information

FI 3103 Quantum Physics

FI 3103 Quantum Physics /9/4 FI 33 Quanum Physcs Aleander A. Iskandar Physcs of Magnesm and Phooncs Research Grou Insu Teknolog Bandung Basc Conces n Quanum Physcs Probably and Eecaon Value Hesenberg Uncerany Prncle Wave Funcon

More information

SOME NOISELESS CODING THEOREMS OF INACCURACY MEASURE OF ORDER α AND TYPE β

SOME NOISELESS CODING THEOREMS OF INACCURACY MEASURE OF ORDER α AND TYPE β SARAJEVO JOURNAL OF MATHEMATICS Vol.3 (15) (2007), 137 143 SOME NOISELESS CODING THEOREMS OF INACCURACY MEASURE OF ORDER α AND TYPE β M. A. K. BAIG AND RAYEES AHMAD DAR Absrac. In hs paper, we propose

More information

( ) () we define the interaction representation by the unitary transformation () = ()

( ) () we define the interaction representation by the unitary transformation () = () Hgher Order Perurbaon Theory Mchael Fowler 3/7/6 The neracon Represenaon Recall ha n he frs par of hs course sequence, we dscussed he chrödnger and Hesenberg represenaons of quanum mechancs here n he chrödnger

More information

Robustness Experiments with Two Variance Components

Robustness Experiments with Two Variance Components Naonal Insue of Sandards and Technology (NIST) Informaon Technology Laboraory (ITL) Sascal Engneerng Dvson (SED) Robusness Expermens wh Two Varance Componens by Ana Ivelsse Avlés avles@ns.gov Conference

More information

The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems

The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems Swss Federal Insue of Page 1 The Fne Elemen Mehod for he Analyss of Non-Lnear and Dynamc Sysems Prof. Dr. Mchael Havbro Faber Dr. Nebojsa Mojslovc Swss Federal Insue of ETH Zurch, Swzerland Mehod of Fne

More information

First-order piecewise-linear dynamic circuits

First-order piecewise-linear dynamic circuits Frs-order pecewse-lnear dynamc crcus. Fndng he soluon We wll sudy rs-order dynamc crcus composed o a nonlnear resse one-por, ermnaed eher by a lnear capacor or a lnear nducor (see Fg.. Nonlnear resse one-por

More information

Outline. Probabilistic Model Learning. Probabilistic Model Learning. Probabilistic Model for Time-series Data: Hidden Markov Model

Outline. Probabilistic Model Learning. Probabilistic Model Learning. Probabilistic Model for Time-series Data: Hidden Markov Model Probablsc Model for Tme-seres Daa: Hdden Markov Model Hrosh Mamsuka Bonformacs Cener Kyoo Unversy Oulne Three Problems for probablsc models n machne learnng. Compung lkelhood 2. Learnng 3. Parsng (predcon

More information

3.2 Models for technical systems

3.2 Models for technical systems onrol Laboraory 3. Mahemacal Moelng 3. Moels for echncal sysems 3.. Elecrcal sysems Fg. 3. shows hree basc componens of elecrcal crcs. Varables = me, = volage [V], = crren [A] omponen parameers R = ressance

More information

TSS = SST + SSE An orthogonal partition of the total SS

TSS = SST + SSE An orthogonal partition of the total SS ANOVA: Topc 4. Orhogonal conrass [ST&D p. 183] H 0 : µ 1 = µ =... = µ H 1 : The mean of a leas one reamen group s dfferen To es hs hypohess, a basc ANOVA allocaes he varaon among reamen means (SST) equally

More information

The Elastic Wave Equation. The elastic wave equation

The Elastic Wave Equation. The elastic wave equation The Elasc Wave Eqaon Elasc waves n nfne homogeneos soropc meda Nmercal smlaons for smple sorces Plane wave propagaon n nfne meda Freqency, wavenmber, wavelengh Condons a maeral dsconnes nell s Law Reflecon

More information

. The geometric multiplicity is dim[ker( λi. number of linearly independent eigenvectors associated with this eigenvalue.

. The geometric multiplicity is dim[ker( λi. number of linearly independent eigenvectors associated with this eigenvalue. Lnear Algebra Lecure # Noes We connue wh he dscusson of egenvalues, egenvecors, and dagonalzably of marces We wan o know, n parcular wha condons wll assure ha a marx can be dagonalzed and wha he obsrucons

More information

10. A.C CIRCUITS. Theoretically current grows to maximum value after infinite time. But practically it grows to maximum after 5τ. Decay of current :

10. A.C CIRCUITS. Theoretically current grows to maximum value after infinite time. But practically it grows to maximum after 5τ. Decay of current : . A. IUITS Synopss : GOWTH OF UNT IN IUIT : d. When swch S s closed a =; = d. A me, curren = e 3. The consan / has dmensons of me and s called he nducve me consan ( τ ) of he crcu. 4. = τ; =.63, n one

More information

Notes on the stability of dynamic systems and the use of Eigen Values.

Notes on the stability of dynamic systems and the use of Eigen Values. Noes on he sabl of dnamc ssems and he use of Egen Values. Source: Macro II course noes, Dr. Davd Bessler s Tme Seres course noes, zarads (999) Ineremporal Macroeconomcs chaper 4 & Techncal ppend, and Hamlon

More information

A NEW TECHNIQUE FOR SOLVING THE 1-D BURGERS EQUATION

A NEW TECHNIQUE FOR SOLVING THE 1-D BURGERS EQUATION S19 A NEW TECHNIQUE FOR SOLVING THE 1-D BURGERS EQUATION by Xaojun YANG a,b, Yugu YANG a*, Carlo CATTANI c, and Mngzheng ZHU b a Sae Key Laboraory for Geomechancs and Deep Underground Engneerng, Chna Unversy

More information

NATIONAL UNIVERSITY OF SINGAPORE PC5202 ADVANCED STATISTICAL MECHANICS. (Semester II: AY ) Time Allowed: 2 Hours

NATIONAL UNIVERSITY OF SINGAPORE PC5202 ADVANCED STATISTICAL MECHANICS. (Semester II: AY ) Time Allowed: 2 Hours NATONAL UNVERSTY OF SNGAPORE PC5 ADVANCED STATSTCAL MECHANCS (Semeser : AY 1-13) Tme Allowed: Hours NSTRUCTONS TO CANDDATES 1. Ths examnaon paper conans 5 quesons and comprses 4 prned pages.. Answer all

More information

Model-Based FDI : the control approach

Model-Based FDI : the control approach Model-Baed FDI : he conrol approach M. Saroweck LAIL-CNRS EUDIL, Unver Llle I Olne of he preenaon Par I : model Sem, normal and no normal condon, fal Par II : he decon problem problem eng noe, drbance,

More information

( ) [ ] MAP Decision Rule

( ) [ ] MAP Decision Rule Announcemens Bayes Decson Theory wh Normal Dsrbuons HW0 due oday HW o be assgned soon Proec descrpon posed Bomercs CSE 90 Lecure 4 CSE90, Sprng 04 CSE90, Sprng 04 Key Probables 4 ω class label X feaure

More information

[ ] 2. [ ]3 + (Δx i + Δx i 1 ) / 2. Δx i-1 Δx i Δx i+1. TPG4160 Reservoir Simulation 2018 Lecture note 3. page 1 of 5

[ ] 2. [ ]3 + (Δx i + Δx i 1 ) / 2. Δx i-1 Δx i Δx i+1. TPG4160 Reservoir Simulation 2018 Lecture note 3. page 1 of 5 TPG460 Reservor Smulaon 08 page of 5 DISCRETIZATIO OF THE FOW EQUATIOS As we already have seen, fne dfference appromaons of he paral dervaves appearng n he flow equaons may be obaned from Taylor seres

More information

Implementation of Quantized State Systems in MATLAB/Simulink

Implementation of Quantized State Systems in MATLAB/Simulink SNE T ECHNICAL N OTE Implemenaon of Quanzed Sae Sysems n MATLAB/Smulnk Parck Grabher, Mahas Rößler 2*, Bernhard Henzl 3 Ins. of Analyss and Scenfc Compung, Venna Unversy of Technology, Wedner Haupsraße

More information

Polymerization Technology Laboratory Course

Polymerization Technology Laboratory Course Prakkum Polymer Scence/Polymersaonsechnk Versuch Resdence Tme Dsrbuon Polymerzaon Technology Laboraory Course Resdence Tme Dsrbuon of Chemcal Reacors If molecules or elemens of a flud are akng dfferen

More information

CHAPTER 10: LINEAR DISCRIMINATION

CHAPTER 10: LINEAR DISCRIMINATION CHAPER : LINEAR DISCRIMINAION Dscrmnan-based Classfcaon 3 In classfcaon h K classes (C,C,, C k ) We defned dscrmnan funcon g j (), j=,,,k hen gven an es eample, e chose (predced) s class label as C f g

More information

. The geometric multiplicity is dim[ker( λi. A )], i.e. the number of linearly independent eigenvectors associated with this eigenvalue.

. The geometric multiplicity is dim[ker( λi. A )], i.e. the number of linearly independent eigenvectors associated with this eigenvalue. Mah E-b Lecure #0 Noes We connue wh he dscusson of egenvalues, egenvecors, and dagonalzably of marces We wan o know, n parcular wha condons wll assure ha a marx can be dagonalzed and wha he obsrucons are

More information

M. Y. Adamu Mathematical Sciences Programme, AbubakarTafawaBalewa University, Bauchi, Nigeria

M. Y. Adamu Mathematical Sciences Programme, AbubakarTafawaBalewa University, Bauchi, Nigeria IOSR Journal of Mahemacs (IOSR-JM e-issn: 78-578, p-issn: 9-765X. Volume 0, Issue 4 Ver. IV (Jul-Aug. 04, PP 40-44 Mulple SolonSoluons for a (+-dmensonalhroa-sasuma shallow waer wave equaon UsngPanlevé-Bӓclund

More information

ON THE WEAK LIMITS OF SMOOTH MAPS FOR THE DIRICHLET ENERGY BETWEEN MANIFOLDS

ON THE WEAK LIMITS OF SMOOTH MAPS FOR THE DIRICHLET ENERGY BETWEEN MANIFOLDS ON THE WEA LIMITS OF SMOOTH MAPS FOR THE DIRICHLET ENERGY BETWEEN MANIFOLDS FENGBO HANG Absrac. We denfy all he weak sequenal lms of smooh maps n W (M N). In parcular, hs mples a necessary su cen opologcal

More information

APPROXIMATE ANALYTIC SOLUTIONS OF A NONLINEAR ELASTIC WAVE EQUATIONS WITH THE ANHARMONIC CORRECTION

APPROXIMATE ANALYTIC SOLUTIONS OF A NONLINEAR ELASTIC WAVE EQUATIONS WITH THE ANHARMONIC CORRECTION THE PUBLISHING HOUSE PROCEEDINGS OF THE ROMANIAN ACADEMY Seres A OF THE ROMANIAN ACADEMY Volume 6 Number /5 pp 8 86 APPROXIMATE ANALYTIC SOLUTIONS OF A NONLINEAR ELASTIC WAVE EQUATIONS WITH THE ANHARMONIC

More information

Solving Parabolic Partial Delay Differential. Equations Using The Explicit Method And Higher. Order Differences

Solving Parabolic Partial Delay Differential. Equations Using The Explicit Method And Higher. Order Differences Jornal of Kfa for Maemacs and Compe Vol. No.7 Dec pp 77-5 Solvng Parabolc Paral Delay Dfferenal Eqaons Usng e Eplc Meod And Hger Order Dfferences Asss. Prof. Amal Kalaf Haydar Kfa Unversy College of Edcaon

More information

F-Tests and Analysis of Variance (ANOVA) in the Simple Linear Regression Model. 1. Introduction

F-Tests and Analysis of Variance (ANOVA) in the Simple Linear Regression Model. 1. Introduction ECOOMICS 35* -- OTE 9 ECO 35* -- OTE 9 F-Tess and Analyss of Varance (AOVA n he Smple Lnear Regresson Model Inroducon The smple lnear regresson model s gven by he followng populaon regresson equaon, or

More information

Lecture 6: Learning for Control (Generalised Linear Regression)

Lecture 6: Learning for Control (Generalised Linear Regression) Lecure 6: Learnng for Conrol (Generalsed Lnear Regresson) Conens: Lnear Mehods for Regresson Leas Squares, Gauss Markov heorem Recursve Leas Squares Lecure 6: RLSC - Prof. Sehu Vjayakumar Lnear Regresson

More information

Linear Response Theory: The connection between QFT and experiments

Linear Response Theory: The connection between QFT and experiments Phys540.nb 39 3 Lnear Response Theory: The connecon beween QFT and expermens 3.1. Basc conceps and deas Q: ow do we measure he conducvy of a meal? A: we frs nroduce a weak elecrc feld E, and hen measure

More information

Appendix H: Rarefaction and extrapolation of Hill numbers for incidence data

Appendix H: Rarefaction and extrapolation of Hill numbers for incidence data Anne Chao Ncholas J Goell C seh lzabeh L ander K Ma Rober K Colwell and Aaron M llson 03 Rarefacon and erapolaon wh ll numbers: a framewor for samplng and esmaon n speces dversy sudes cology Monographs

More information

THE PREDICTION OF COMPETITIVE ENVIRONMENT IN BUSINESS

THE PREDICTION OF COMPETITIVE ENVIRONMENT IN BUSINESS THE PREICTION OF COMPETITIVE ENVIRONMENT IN BUSINESS INTROUCTION The wo dmensonal paral dfferenal equaons of second order can be used for he smulaon of compeve envronmen n busness The arcle presens he

More information

A Functional-Link-Based Fuzzy Neural Network for Temperature Control

A Functional-Link-Based Fuzzy Neural Network for Temperature Control Proceedngs of he 7 IEEE Symposm on Fondaons of Compaonal Inellgence (FOCI 7) A Fnconal-Ln-Based Fzzy Neral Neor for emperare Conrol Cheng-Hng Chen *, Chn-eng Ln, Fello, IEEE, and Cheng-Jan Ln, ember, IEEE

More information

The Shapley value for fuzzy games on vague sets

The Shapley value for fuzzy games on vague sets WEA TRANACTIN on INFRMATIN CIENCE APPLICATIN The hapley value or uzzy games on vague ses Fan-Yong Meng* (Correspondng Auhor chool o Managemen Qngdao Technologcal nversy Qngdao 266520 hong Provnce P R Chna

More information

Robust and Accurate Cancer Classification with Gene Expression Profiling

Robust and Accurate Cancer Classification with Gene Expression Profiling Robus and Accurae Cancer Classfcaon wh Gene Expresson Proflng (Compuaonal ysems Bology, 2005) Auhor: Hafeng L, Keshu Zhang, ao Jang Oulne Background LDA (lnear dscrmnan analyss) and small sample sze problem

More information

Department of Economics University of Toronto

Department of Economics University of Toronto Deparmen of Economcs Unversy of Torono ECO408F M.A. Economercs Lecure Noes on Heeroskedascy Heeroskedascy o Ths lecure nvolves lookng a modfcaons we need o make o deal wh he regresson model when some of

More information

FTCS Solution to the Heat Equation

FTCS Solution to the Heat Equation FTCS Soluon o he Hea Equaon ME 448/548 Noes Gerald Reckenwald Porland Sae Unversy Deparmen of Mechancal Engneerng gerry@pdxedu ME 448/548: FTCS Soluon o he Hea Equaon Overvew Use he forward fne d erence

More information

Mechanics Physics 151

Mechanics Physics 151 Mechancs Physcs 5 Lecure 9 Hamlonan Equaons of Moon (Chaper 8) Wha We Dd Las Tme Consruced Hamlonan formalsm H ( q, p, ) = q p L( q, q, ) H p = q H q = p H = L Equvalen o Lagrangan formalsm Smpler, bu

More information

2.1 Constitutive Theory

2.1 Constitutive Theory Secon.. Consuve Theory.. Consuve Equaons Governng Equaons The equaons governng he behavour of maerals are (n he spaal form) dρ v & ρ + ρdv v = + ρ = Conservaon of Mass (..a) d x σ j dv dvσ + b = ρ v& +

More information

Mechanics Physics 151

Mechanics Physics 151 Mechancs Physcs 5 Lecure 9 Hamlonan Equaons of Moon (Chaper 8) Wha We Dd Las Tme Consruced Hamlonan formalsm Hqp (,,) = qp Lqq (,,) H p = q H q = p H L = Equvalen o Lagrangan formalsm Smpler, bu wce as

More information

WiH Wei He

WiH Wei He Sysem Idenfcaon of onlnear Sae-Space Space Baery odels WH We He wehe@calce.umd.edu Advsor: Dr. Chaochao Chen Deparmen of echancal Engneerng Unversy of aryland, College Par 1 Unversy of aryland Bacground

More information

Wronskian Determinant Solutions for the (3 + 1)-Dimensional Boiti-Leon-Manna-Pempinelli Equation

Wronskian Determinant Solutions for the (3 + 1)-Dimensional Boiti-Leon-Manna-Pempinelli Equation Jornal of Appled Mahemacs and Physcs 0 8-4 Pblshed Onlne ovember 0 (hp://www.scrp.org/jornal/jamp) hp://d.do.org/0.46/jamp.0.5004 Wronskan Deermnan Solons for he ( + )-Dmensonal Bo-Leon-Manna-Pempnell

More information

On computing differential transform of nonlinear non-autonomous functions and its applications

On computing differential transform of nonlinear non-autonomous functions and its applications On compung dfferenal ransform of nonlnear non-auonomous funcons and s applcaons Essam. R. El-Zahar, and Abdelhalm Ebad Deparmen of Mahemacs, Faculy of Scences and Humanes, Prnce Saam Bn Abdulazz Unversy,

More information

P R = P 0. The system is shown on the next figure:

P R = P 0. The system is shown on the next figure: TPG460 Reservor Smulaon 08 page of INTRODUCTION TO RESERVOIR SIMULATION Analycal and numercal soluons of smple one-dmensonal, one-phase flow equaons As an nroducon o reservor smulaon, we wll revew he smples

More information

Lecture 9: Dynamic Properties

Lecture 9: Dynamic Properties Shor Course on Molecular Dynamcs Smulaon Lecure 9: Dynamc Properes Professor A. Marn Purdue Unversy Hgh Level Course Oulne 1. MD Bascs. Poenal Energy Funcons 3. Inegraon Algorhms 4. Temperaure Conrol 5.

More information

Graphene nanoplatelets induced heterogeneous bimodal structural magnesium matrix composites with enhanced mechanical properties

Graphene nanoplatelets induced heterogeneous bimodal structural magnesium matrix composites with enhanced mechanical properties raphene nanoplaele nce heerogeneo bmoal rcral magnem marx compoe wh enhance mechancal propere Shln Xang a, b, Xaojn Wang a, *, anoj pa b, Kn W a, Xaoh H a, ngy Zheng a a School of aeral Scence an ngneerng,

More information

Variational method to the second-order impulsive partial differential equations with inconstant coefficients (I)

Variational method to the second-order impulsive partial differential equations with inconstant coefficients (I) Avalable onlne a www.scencedrec.com Proceda Engneerng 6 ( 5 4 Inernaonal Worksho on Aomoble, Power and Energy Engneerng Varaonal mehod o he second-order mlsve aral dfferenal eqaons wh nconsan coeffcens

More information

Fall 2010 Graduate Course on Dynamic Learning

Fall 2010 Graduate Course on Dynamic Learning Fall 200 Graduae Course on Dynamc Learnng Chaper 4: Parcle Flers Sepember 27, 200 Byoung-Tak Zhang School of Compuer Scence and Engneerng & Cognve Scence and Bran Scence Programs Seoul aonal Unversy hp://b.snu.ac.kr/~bzhang/

More information

WebAssign HW Due 11:59PM Tuesday Clicker Information

WebAssign HW Due 11:59PM Tuesday Clicker Information WebAssgn HW Due 11:59PM Tuesday Clcker Inormaon Remnder: 90% aemp, 10% correc answer Clcker answers wll be a end o class sldes (onlne). Some days we wll do a lo o quesons, and ew ohers Each day o clcker

More information

Existence and Uniqueness Results for Random Impulsive Integro-Differential Equation

Existence and Uniqueness Results for Random Impulsive Integro-Differential Equation Global Journal of Pure and Appled Mahemacs. ISSN 973-768 Volume 4, Number 6 (8), pp. 89-87 Research Inda Publcaons hp://www.rpublcaon.com Exsence and Unqueness Resuls for Random Impulsve Inegro-Dfferenal

More information

RELATIONSHIP BETWEEN VOLATILITY AND TRADING VOLUME: THE CASE OF HSI STOCK RETURNS DATA

RELATIONSHIP BETWEEN VOLATILITY AND TRADING VOLUME: THE CASE OF HSI STOCK RETURNS DATA RELATIONSHIP BETWEEN VOLATILITY AND TRADING VOLUME: THE CASE OF HSI STOCK RETURNS DATA Mchaela Chocholaá Unversy of Economcs Braslava, Slovaka Inroducon (1) one of he characersc feaures of sock reurns

More information

Delay Dependent Robust Stability of T-S Fuzzy. Systems with Additive Time Varying Delays

Delay Dependent Robust Stability of T-S Fuzzy. Systems with Additive Time Varying Delays Appled Maemacal Scences, Vol. 6,, no., - Delay Dependen Robus Sably of -S Fuzzy Sysems w Addve me Varyng Delays Idrss Sad LESSI. Deparmen of Pyscs, Faculy of Scences B.P. 796 Fès-Alas Sad_drss9@yaoo.fr

More information

Solution of a diffusion problem in a non-homogeneous flow and diffusion field by the integral representation method (IRM)

Solution of a diffusion problem in a non-homogeneous flow and diffusion field by the integral representation method (IRM) Appled and ompaonal Mahemacs 4; 3: 5-6 Pblshed onlne Febrary 4 hp://www.scencepblshnggrop.com//acm do:.648/.acm.43.3 olon of a dffson problem n a non-homogeneos flow and dffson feld by he negral represenaon

More information

New M-Estimator Objective Function. in Simultaneous Equations Model. (A Comparative Study)

New M-Estimator Objective Function. in Simultaneous Equations Model. (A Comparative Study) Inernaonal Mahemacal Forum, Vol. 8, 3, no., 7 - HIKARI Ld, www.m-hkar.com hp://dx.do.org/.988/mf.3.3488 New M-Esmaor Objecve Funcon n Smulaneous Equaons Model (A Comparave Sudy) Ahmed H. Youssef Professor

More information

EEL 6266 Power System Operation and Control. Chapter 5 Unit Commitment

EEL 6266 Power System Operation and Control. Chapter 5 Unit Commitment EEL 6266 Power Sysem Operaon and Conrol Chaper 5 Un Commmen Dynamc programmng chef advanage over enumeraon schemes s he reducon n he dmensonaly of he problem n a src prory order scheme, here are only N

More information