arxiv: v1 [math.oc] 8 Sep 2015

Size: px
Start display at page:

Download "arxiv: v1 [math.oc] 8 Sep 2015"

Transcription

1 Non-lnear Graden Algorhm for Parameer Esmaon: Exended verson Juan G Rueda-Escobedo and Jame A oreno arxv: v [mahoc] 8 Sep 05 Absrac Graden algorhms are classcal n adapve conrol and parameer esmaon For nsananeous quadrac cos funcons hey lead o a lnear me-varyng dynamc sysem ha converges exponenally under perssence of excaon condons In hs paper we consder nsananeous non-quadrac cos funcons, for whch he graden algorhm leads o non-lnear and non Lpschz me-varyng dynamcs, whch are homogeneous n he sae We show ha under perssence of excaon condons hey also converge globally, unformly and asympocally Compared o he lnear counerpar, hey accelerae he convergence and can provde for fne-me or fxed-me sably I INTRODUCTION Ths work s an exended verson of [] In hs paper he proof of he clams n [] are presened n he appendces Ths proof are omed n he frs work because of s lengh, nsead, he space s used o dscuss and clarfy he resuls A classcal lnear paramerc model s gven by y := u T θ 0, where θ 0 R n are he parameers, u R n s he regressor and y R s he measured sgnal Oupu y along wh an esmae of he parameers ˆθ s used o buld he oupu esmaon error e := u T ˆθ y, whch can be rewren as e = u T x, where x := ˆθ θ 0 s he parameer esmaon error The am s o use e o drve x o zero Snce Perssence of Excaon PE of u s equvalen o he Unform and Complee Observably of he assocaed lnear dynamcal sysem θ 0 = 0, y = u T θ 0 [], [3], s a necessary condon o assure unform and robus convergence of any algorhm In parcular, for he Lnear Graden Descen and he Recursve Leas Square ehods PE s a necessary and suffcen condon for exponenal convergence [], [3] In fac, he use of correcon erms lnear n e canno provde for convergence faser han exponenal So, f acceleraed convergence s desred, algorhms usng nonlnear correcng erms are requred In recen years, J G Rueda-Escobedo and J A oreno are wh Elécrca y Compuacón, Insuo de Ingenería, Unversdad Naconal Auónoma de éxco, 0450 éxco DF, exco, JRuedaE@ngenunammx; JorenoP@unammx he use of homogeneous sysems and homogeneous correcon erms for conrol and observaon purposes have been very successful n provdng fne me and fxed me convergence [4], [5] Furhermore, homogeneous hgher order sldng modes HOS provde for dsconnuous correcon erms, whch provde no only fne me convergence bu also nsensvy o mached, bounded perurbaons [6], [7], [8] Homogeney has been mporan for hese resuls, snce provdes useful properes: eg local asympoc sably s equvalen o global fne me sably for sysems of negave homogeney degree [9], [6], [0] ovaed by hese resuls we propose n hs work non-lnear esmaon algorhms, wh nonlnear correcon erms n e, ha lead o me-varyng dynamcal sysems The proposed schemes can be obaned as he negave graden of an nsananeous convex non smooh funcon of he oupu esmaon error The resulng error sysems are homogeneous or homogeneous n he b-lm n he esmaon error x [9], [], bu me-varyng Unlke me nvaran homogeneous sysems, me-varyng sysems do no possess such srong properes For example, neher local asympoc sably mples global asympoc sably, nor negave homogeney degree mples fne me convergence [] Some local asympoc sably resuls for me-varyng sysems homogeneous n he sae have been obaned, when he homogeney degree s zero [], and for posve homogeney degree [], usng averagng echnques We show n hs paper, n sraghforward a drec manner, ha he proposed algorhms converge globally and asympocally under he PE condon of he regressor When n =, he algorhm s able o converge n fne-me and s also able o esmae a me-varyng parameer oreover, addng homogeneous erms of posve degree he esmaon error converges n fxedme, e he convergence me s upper bounded by a consan ndependen of he nal esmaon error For n > global, unform and asympoc sably can be assured and acceleraon for large nal condons can be obaned

2 II OTIVATION AND PROBLE STATEENT The esmaon problem may be regarded as a mnmzaon problem In hs conex a common cornersone s o se a convex cos funcon of he oupu esmaon error In hs work we choose he followng srucure for he cos funcon Jˆθ p := u T ˆθ y p+ p + Here he exponen p > 0 s a parameer o be chosen The erm ˆθ Jˆθ p shows he drecon n whch he parameer esmaed needs o change Wh hs dea n mnd we propose he followng algorhm ˆθ := ˆθ Jˆθ p = u T ˆθ y p sgn u T ˆθ y u For he sake of readably le us defne w p := w p sgnw Wh hs convenon he algorhm s rewren as ˆθ = u T ˆθ y p u We denoe as compose algorhm an algorhm ha resuls from addng vecor felds of he form n o aval of he dynamcs ras assocaed o each p Ths leads o ˆθ = u T ˆθ y p u = To analyse he convergence of he algorhms o θ 0 he dynamcs of he esmaon error x = ˆθ θ 0 s needed Hence s necessary o compue he me dervave of x In he case of he sngle algorhm n hs error dynamcs s as follow ẋ = u T x p u 3 Consderng f, x = u T x p u and replacng he argumen x by ɛx, wh ɛ > 0, we have he followng relaon f, ɛx = ɛ p ɛf, x Takng hs no consderaon and followng Defnon n [], we can sad ha he sysem 3 s homogeneous wh homogeney degree p When 0 < p < he homogeney degree s negave; for p = s zero; and for p > he homogeney degree s posve Now, he convergence analyss for hs class of algorhms reduces o esablsh he sably and he aracvy of he orgn of a me-varan homogeneous sysem Repeang he same analyss for he compose algorhm he error dynamcs s ẋ = u T x p u 4 = Please, noce ha hs sysem s no longer homogeneous A quck sably check wh a quadrac Lyapunov funcon V x := xt x yelds o global unform sably GUS of he orgn of 3 and 4 The me dervave of V s presened for boh cases n order o show s negave sem-defnveness: along 3 V = u T x p+ 0, 5a along 4 V = u T x p+ 0 5b = Noce ha he erm u T x can vansh ousde he se { u = 0} { x = 0} To asser he unform asympoc sably UAS of he orgn of 3 and 4 u needs o be PE In [3] a convenen descrpon of he Perssen Excaon s gven Ths descrpon s expressed as a lower bound of he negral of he regressor and s found convenen because fs well n he sudy of he esmaon-error convergence For convenence he defnon s reproduced below Defnon : Le u : R + R n be a pecewse connuous funcon I s sad ha u s of PE f here exs T > 0 and ɛ > 0 such ha T +T u T sw ds ɛ, for all w R n wh w = To conclude hs secon, we presen some defnons of me-varyng sysems Defnon : Le a me-varyng sysem be represened by ẋ = fx, where f0, = 0 for all Le Ω be a conneced open subse of R n, such ha 0 Ω The pon x = 0 s Unformly fne-me sable UFTS f s unformly sable and for any x 0 Ω here exs 0 Tx 0 < + such ha x, 0, x 0 = 0 for all 0 + Tx 0 Also, f Ω = R n hen x = 0 s sad o be globally unformly fne-me sable GUFTS Unformly fxed-me sable UFxTS f s GUFTS and exs T < + such ha T Tx 0 for every x 0 R n III AIN RESULT In hs secon he sably of 3 and 4 s presened Two cases are recognzed: when only one parameer needs o be esmaed scalar case and when are more han one vecor case Ths dvson s done because he resuls n he vecor case do no reflec ceran phenomena ha occurred n he scalar one The proof of he heorems can be found n he Appendces

3 A Scalar Case The resuls n he scalar case are sronger han n he vecor one due o he fac ha he produc u T x can only be zero f one or boh of he varables are zero Perssen excaon prevens u o say n zero or o exhb a growng dwellng me n zero Ths mples ha, for sayng n zero, x needs o be zero When he exponen p s chosen n he nerval [0, he algorhm converges n unform fne-me and a bound of he convergence me s gven as an neger mulple of he perssen excaon perod T The nex saemen summarzes hs dscusson Theorem : Le u be a pece-wse connuous funcon of and of PE Le 0 p < and n =, hen he orgn of 3 s globally unformly fne-me sable An upper bound of he convergence me s Tx 0 x0 p pt ɛ p+ T 6 where x 0 s he esmaon error a he nal me; T and ɛ are as n Defnon and denoes he celng funcon When p = he lnear case s obaned and s properes are well know [], [3], for hs s lef ou of he dscusson For p > only global unform asympoc sably GUAS can be assered bu oher neresng propery arse No maer how large he nal error s, he me τc need o reach a smaller level se V x = c s only funcon of c Ths s referred o as escape from nfne n fne me unformly n We gaher hese resuls n he followng heorem Theorem : Le u be a pece-wse connuous funcon of and of PE, p > and n =, hen he orgn of 3 s globally unformly asympocally sable Furhermore, he me needed o escape from nfny o a compac regon V x c s bounded by τc p p T ɛ p+ c p T 7 The sudy of he convergence of he compose algorhm can be done usng he prevous resuls For he scalar case he sably can be assered va a Comparson Lemma for dfferenal nequales [4] In shor, he rajecores of he sysem 4 are below he rajecores correspondng o each of he sngle algorhms wh every solaed exponen p Ths can be seen from 5b whch can be rewren as V = = p + u p+ V p + p + u p+ V p +, =,,, h An mporan case occurs when a leas one exponen s n [0, and anoher one s greaer han one For hs case he me needed for he algorhm o converge o θ 0 s ndependen of he nal esmaon error and he nal me f u s of PE Ths s summarzed n he followng heorem Theorem 3: Le u be a pece-wse connuous funcon of and of PE, also n = Consder he sysem 4 and le P := {p, p,, p h } be he se of he exponens Denoe p m he mnmum elemen n P and p he maxmum Assume ha 0 p m < and p > Le P m be he subse of he exponen smaller han one and P he subse of he exponens greaer han one Defne p P as he exponen whch maxmzes p ɛ p+ and p P m he exponen whch maxmzes p ɛ p+, hen an upper bound of he me needed o converge o zero s T pt ɛ p+ + p p T ɛ p+ T Las bu no leas, a dsconnuous algorhm capable of esmang one varyng parameer s presened Ths algorhm makes use of he regressor sgn and he sgn of he oupu esmaon error ˆθ = L sgn uˆθ y sgn u, 8 wh L > 0 Now assume ha he parameer varaon s bounded, e θ [ γ, γ], γ 0 Also s consdered ha u s of PE and canno say n zero for me nervals, bu can cross The error dynamcs nduced by 8 s ẋ L sgn x + θ 9 whch s a dfferenal ncluson [9] By employng V x = x as Lyapunov Funcon s dervave along he rajecores of 9 s V = L x θx = L + θ sgnx x, whch s negave f L > γ and u 0 If u say n zero he rack of θ s los bu can be recovered laer f u s of PE To guaranee exac rackng, u canno say n zero B Vecor Case When θ 0 s a vecor he se where u T x s zero grows Ths changes he general behavour of our algorhms Only GUAS can be guaraneed n general wh perssen excaon Also, he dsconnuous algorhm whch resul of selecng p = 0 does no converge In he laer case, a sgnal of PE can be consruced for whch he oupu esmaon error e becomes zero n fne

4 me bu x does no reach he orgn The sably properes of he orgn of 3 are summarzed n he followng heorem Theorem 4: Le u be a pece-wse connuous funcon of and of PE, unformly bounded by u, hen he orgn of 3 s globally unformly asympocally sable for any p > 0 Alhough only GUAS can be assered wh PE n general, n he followng secon wo classes of sgnals of PE are presened One of hem guaranees global unform exponenal sably GUES and he oher GUFTS; n boh cases for 0 < p < Ths means ha he sgnals u whch can provde UFTS are n a subse of hose of PE The compose algorhm sll works n he vecor case bu nohng more han GUAS can be clamed In conras o he scalar case he sably of he compose vecor algorhm canno be obaned va he Comparson Lemma The followng heorem synheszes hs dscusson Theorem 5: Le u be a pece-wse connuous funcon of and of PE, unformly bounded by u, hen he orgn of 4 s globally unformly asympocally sable f p p m + where p m s he mnmum exponen n he se {p, p,, p h } and p s he maxmum The exra condon regardng he exponens appears due o he way he proof was done and we hnk s no nrnsc o he sably In he scalar case he Comparson Lemma gves nformaon abou he relaonshp beween he rajecores of he sngle algorhm bu n hs case nohng can be concluded In he nex secon smulaon examples are presened for he compose algorhm In he Fgure 4 and Fgure 5 can be seen ha he Lyapunov funcon s below of he correspondng one for he sngle algorhms bu hs does no need o be rue n general IV EXAPLES In hs secon examples are presened In he scalar case, smulaons of a parameer esmaon process are shown, whereas for he vecor case also he behavour of he error s suded for a specfc class of sgnals A Scalar case Numercal smulaon s performed o llusrae he dfference beween he classc graden algorhm and he famly presened n hs work For hs am we choose: θ 0 = 5, ˆθ0 = 0, u = cos Four algorhms were smulaed: one wh p = 3 4, oher wh p =, a hrd one wh p 3 = 3 and a compose one wh p and p 3 Fgure and Fgure show he behavour of he Lyapunov funcon V = ˆθ θ 0 In he frs fgure he nal value of V and s decayng behavour are presened for all he algorhms As can be seen here s an orderng n he decay: Compose > p = 3 > p = > p = 3 4, and hs s due o he nal value of V 0 >> When he esmaon error becomes smaller he order changes as n Fgure In Fgure he behavour s shown when V < 004 Fg Behavour of V n he Scalar case a In hs regon he decay order s: Compose > p = 3 4 > p = > p = 3 as s expeced from he knd of convergence of every algorhm: u fxed-me > u fne-me > u exponenal > u asympoc Ths order does no change for any fuure me Now a esmaon of a me-varyng parameer s done Fg Behavour of V n he Scalar case b va 8 The smulaon parameers are: θ = cos3 4, ˆθ0 = 0, u = cos + 5, L = 33 Snce he regressor s non-zero for any me, an exac rackng of he parameer s acheved as shown n Fgure 3 B Vecor case Frs smulaon resuls are presened for he esmaon process of hree parameers wh and Four algorhms were smulaed: he frs wh p = 3/4, he second wh p =, he hrd one wh p 3 = 3/ and he las one wh p and p 3 The smulaon resuls are shown n Fgure 4 and Fgure 5 o compare he behavour of he algorhms, when dfferen exponens are chosen Laer, he response of he algorhm for 0 < p < s found for

5 Fg 3 Tme-varyng parameer esmaon process Fg 5 Smulaon resuls vecor case b pecewse consan sgnal and he class of convergence s esablshed for For he smulaon example he nex condons were seleced: u = [ cos cos3 5 cos5 ]T, θ 0 = [ ] T 3 4, ˆθ 0 = [ ]T As n he las secon only he plo of V = ˆθ θ 0 T ˆθ θ 0 s presened In Fgure 4 he decay of V s shown a he begnnng of he process and he decay order s preserved Also n Fgure 5 he change of order ha was found for he scalar case s obaned However we do no see, n general, a change of order for arbrary parameer values Now he response of he sysem 3 when u s Fg 4 Smulaon resuls for vecor case a pecewse consan s analyzed Le us defne a new varable as follows z := u T x, s dynamcs s ż = u T ẋ = u z p snce u = 0 almos everywhere, where denoes he Eucldean norm For smplcy we assume ha he lengh of he nervals where u remans consan has a consan value τ > 0 The soluon of z n he nerval [, +τ] wh u = µ, u µ > 0 and p > 0, p s z = z p p µ p sgnz For 0 < p < he expresson above holds f z p p u and z = 0 oherwse The soluon of z can be used o fnd he soluon of x by nocng ha ẋ = z p u Ths s gven n he forhcomng expresson x = x µ µ µ T x sgnµ T x µ T x p p µ p Resrcng he analyss for 0 < p < s clear ha, f τ s large enough hen x + τ s orhogonal o µ x + τ = x µ µµt x = I n µ µµt x Fx τ, f x p pu p+ m τ hen x + τ becomes orhogonal o µ Ths means ha here exss a ball cenered n zero for whch τ s always large enough o make x orhogonal o any µ If he sequence U = {µ } = s chosen o fulfl he noon of perssen excaon for dscree sysems n [5], hen he orgn s GUAS by Theorem 4 and x can reach any ball cenered n zero n fne me, e always reach he ball n where τ guaranees ha x becomes orhogonal o µ, and hs yelds a dscree sysem whch s descrbed by he followng dfference equaon x k+ = I n µ k µ kµ T k x k Only assumng PE of he sequence U, GUES of he orgn can be concluded makng an analogue analyss o he one shown n he proof for Theorem 4 n Appendx II-D bu akng V k x = x T k x k, V = V k+t V k x nsead of V x and V respecvely Now ake n muual orhogonal vecors {v, v,, v n }, v T v j = 0, j, and exce he sysem wh hem, hen x k+n = 0 x k+n = = n = I n v v v T x k I n n = v v v T x k = 0 A sequence consruced wh hese vecors s of PE and also makes he sysem GUFTS Ths shows ha perssen excaon canno guaranee fne-me convergence n he vecor case bu does no forbd

6 V CONCLUSIONS In hs work a parameer esmaon echnque s presened Wh he proposed algorhms we obaned fneme and fxed-me convergence o he rue parameers However hs properes canno be guaraneed n general A deep sudy of he sgnals ha can asser such mporan properes s sll needed Even hough he algorhms were seleced o make he error dynamcs homogeneous n he sae, hs does no help n he analyss However, he homogeneous nonlneares can enhance he robusness properes of our algorhms wr addve perurbaons n comparson wh classcal approaches ACKNOWLEDGEENT The auhors hank he fnancal suppor from PAPIIT- UNA Programa de Apoyo a Proyecos de Invesgacón e Innovacón Tecnológca, projec IN364; Fondo de Colaboracón II-FI UNA, Projec IISGBAS- -04; CONACyT Consejo Naconal de Cenca y Tecnología, projec 47; and CONACyT CVU: 4970 REFERENCES [] J G Rueda-Escobedo and J A oreno, Non-lnear graden algorhm for parameer esmaon, n Decson and Conrol CDC, 05 IEEE 54h Annual Conference on, 05, p o appear [] B D O Anderson, Exponenal sably of lnear equaons arsng n adapve denfcaon, IEEE Trans Auoma Conr, vol, no, pp 83 88, 977 [3] A P organ and K S Narendra, On he unform asympoc sably of ceran lnear nonauonomous dfferenal equaons, SIA Journal on Conrol and Opmzaon, vol 5, no, pp 5 4, Jan 977 [4] S Bha and D Bernsen, Geomerc homogeney wh applcaons o fne-me sably, ahemacs of Conrol, Sgnals, and Sysems, vol 7, no, pp 0 7, 005 [5] V Andreu, L Praly, and A Asolf, Homogeneous approxmaon, recursve observer desgn and oupu feedback, SIA J Conrol Opm, vol 47, no 4, pp , 008 [6] A Levan, Homogeney approach o hgh-order sldng mode desgn, Auomaca, vol 4, pp , 005 [7] J oreno, Lyapunov approach for analyss and desgn of second order sldng mode algorhms, n Sldng odes afer he frs decade of he s Cenury, ser LNCIS, 4, L Frdman, J oreno, and R Irare, Eds Berln - Hedelberg: Sprnger- Verlag, 0, pp 3 50 [8] E Bernuau, D Efmov, W Perruque, and A Polyakov, On homogeney and s applcaon n sldng mode, Inernaonal Journal of Frankln Inue, vol 35, no 4, pp , 04 [9] A Bacco and L Roser, Lapunov Funcons and Sably n Conrol Theory, ser Communcaon and Conrol Engneerng Berln - Hedelberg: Sprnger-Verlag, 005, ch 5, pp [0] E Bernuau, D Efmov, W Perruque, and A Polyakov, On an exenson of homogeney noon for dfferenal nclusons, n European Conrol Conference, Zurch, Swzerland, Jul 03 [] J Peueman and D Aeyels, Averagng resuls and he sudy of unform asympoc sably of homogeneous dfferenal equaons ha are no fas me-varyng, SIA Journal on Conrol and Opmzaon, vol 37, no 4, pp , 999 [] R T Closkey, An averagng heorem for me-perodc degree zero homogeneous dfferenal equaons, Sysems & Conrol Leers, vol 3, pp 79 8, 997 [3] K S Narendra and A Annaswamy, Perssen excaon n adapve sysems, Inernaonal Journal of Conrol, vol 45, no, pp 7 60, 987 [4] H K Khall, Nonlnear Sysems Upper Saddle Rver, New Jersey 07458: Prence Hall, 00, ch 3, pp 0 04 [5] T-H Lee and K S Narendra, Robus adapve conrol of dscree-me sysems usng perssen excaon, Auomaca, vol 4, no 6, pp , 988 [6] K S Narendra and A Annaswamy, Sable Adapve Sysems Englewood Clffs, New Jersey 0763: Prence Hall, 989, ch, pp 7 74 APPENDIX I PERSISTENT EXCITATION To prove unform asympoc sably raher han unform sably he perssen excaon n he regressor s needed For usng he propery adequaely he nex proposon s developed Proposon : If u s of PE, hen he followng nequaly holds for T, ɛ and w as n Defnon and p 0 +T u T sw p+ ds T ɛ p+ Proof: Applyng he Hölder nequaly o f = u T w and g = on he nerval [, +T ] and usng as Hölder conjugaes we oban p+ and p p+ +T T p u T sw p+ ds usng Defnon leads o +T +T u T sw p+ ds T p+ u T sw ds p+, T p ɛp+ = T ɛ p+ Ths dervaon from he PE s done n he am of easly presen he proof of he heorems A Proof of Theorem APPENDIX II PROOF OF THEORES TO 5 For n = he erm u T x can be rewren as u x = u V replacng hs n 5a yelds V = p+ u p+ V p+ Solvng he dfferenal equaon we have V = V p 0 p p 0 us p+ ds p 0 Ths soluon s vald for p 0 and p For 0 p < he soluon exss f V p 0 p p 0 us p+ ds afer ha V = 0 The PE guaranees ha here exss 0 when he nequaly no longer holds To esmae hs me s suffcen o fnd an neger k such ha he negral of u p+ from

7 0 o 0 + kt s greaer han p / pv p 0 Usng Proposon 0+kT p us p+ ds kt ɛ p+ 0 p V p 0 Solvng for k and akng he leas neger ha fulflls he nequaly yelds p p k = V 0, pt ɛp+ hen he me ha guaranees ha V reaches zero s = k T + 0 Noce ha V 0 can be replaced by x 0 n o oban 6 B Poof of Theorem Consder agan he soluon for V n 0 Snce p >, Equaon 0 can be rewren as V = V p 0 + p p 0 us p+ p ds As u s of PE he denomnaor grows unbounded makng V 0 as Now an esmae of he me needed for V o decreases from V 0 o a value equal or smaller han c s calculaed Subsung V for c n 5a and evaluang he negral from 0 o s clear ha he value of he negral needs o be larger enough o sasfy he nex nequaly 0 us p+ ds p p c p V p 0 Fxng 0 as an neger mulple of T, e 0 = kt, from Proposon we know ha he negral s greaer or equal o kt ɛ p+ Forcng he RHS of he las nequaly o be less han kt ɛ p+ we ge kt ɛ p+ p p c p V p 0 Now solvng for k and akng he smalles neger ha fulfl he nequaly yelds k, p p T ɛ p+ c p V p 0 k = p p T ɛ p+ c p V p 0 Takng he lm when V 0 he bound 7 s found k = p T ɛ p+ p c p C Proof of Theorem 3 From 5b follows ha V u T x p+ = p + u V p +, whch defne a dfferen dfferenal nequaly for each erm n he sum From he Comparson Lemma we known ha he soluon of V s for below of each soluon of V p, where V p = u V p + ; V p ake he form of 0 Le p + P and P m be as n he heorem For each p P we can asser ha V p can escape from nfny o a compac se n fne me and so V Fx he level se as V x = and esmae he me needed o reach τ p p T ɛ p+ T, p P The smalles me ha he algorhm can guaranee s obaned when p ɛ p+ s maxmum Take ha quany as he esmae Now, wh p j P m, we can esmae he me needed for each V pj o converge from he level se V x = o zero Tx 0 p j T ɛ pj+ T, p j P m, x 0 V x0 = Agan, he smalles me ha he algorhm can guaranee s when p j ɛ pj+ s maxmzed Then he me needed by he algorhm o converge s, a mos, he sum of he wo esmaes D Proof of Theorem 4 Followng he dea n [6] a lower bound of he negral of u T x p+ s needed Take he erm u T x and add a zero n he form u T x u T x nsde he absolue value, by means of he rangle nequaly he nex paron s obaned u T x u T x u T x x Rsng boh sdes o p+ and usng he Jensen nequaly afer ha, he nequaly yelds u T x p+ p ut x p+ u T x x p+ Inegrang boh sdes from o +T and usng Proposon o bound he frs erm n he RHS becomes +T u T x p+ d T ɛp+ x p+ p +T u T x x p+ d To bound he magnude of he second erm n he RHS of he las nequaly, he nex procedure may be used

8 assumng u s unformly bounded, e u u : +T u T sx xs p+ ds u p+ +T x xs p+ ds u p+ T sup x xs p+ s [, +T ] u p+ T +T p+ ẋs ds 3 The norm of ẋ n hs case s u T x p u whch s less han u T x p u Hölder nequaly can be used wh f = u T x p and g = n he nerval [, + T ] wh p+ p and p + as Hölder conjugaes, hen +T p T u T sxs p+ ds +T p+ u T sxs p ds Usng hs n 3 and hen he resul n one ges +T u T x p+ d T ɛp+ p x p+ T u p+ +T p u T x p+ d Ths can be rewren as T ɛ p+ p x p+ +T u T x p+ d + T u p+ +T p, u T x p+ d By algebrac manpulaon and defnng z := u T x p+ d he noaon can be smplfed: +T z + T u p+ z p T ɛp+ p x p+ Noce ha he polynomal P z := z + T u p+ z p s a src monooncally ncreasng funcon for z 0 and s nverse exss By denong hs as P z and usng n he nequaly above, one ges +T z = u T x p+ d P T ɛ p+ x p+ p Recallng V from 5a and negrang for he same me nerval s obvous ha +T V d P T ɛ p+ p x p+ < 0 Snce P s also a src monooncally ncreasng funcon and from Theorem 5 n [3] he global unform asympoc sably of he orgn of 3 s sablshed for any p > 0 E Proof of Theorem 5 In 5b several erms of he form u T x p+ appear; an analyss of each erm separaely s done before sudy he full dervave of V By akng he erm u T x and addng a zero n he form u T x u T x nsde he absolue sgn and hen applyng he rangle nequaly he nex paron s found u T x u T x x + u T x, rsng o p + and usng he Jensen Inequaly p ut x p+ u T x x p+ Solvng for u T x p+ [, + T ] +T u T sxs p+ ds p + u T x p+ and negrang beween +T u T sx p+ ds +T u T sx xs p+ ds A lower bound of he frs erm n he RHS of he nequaly can be found usng Proposon The followng nequaly s obaned +T u T sxs p+ ds T ɛp+ p x p + +T u T sx xs p+ ds 4 A upper bound of he magnude of he second erm s que more dffcul Assumng u s unformly bounded by u, e u u, he negral can be bounded as +T u T sx xs p+ ds +T u p+ x xs p+ ds u p+ T sup x xs p+ s [,+T ] u p+ T +T p+ ẋs ds 5 Now an analyss of he negral of ẋ s requred Frs he norm of ẋ s esmaed ẋ = u T x p j u u u T x p j j= j= Inegrang and usng he Jensen nequaly for he funcon p+ +T p + ẋs ds u p + hp +T j= p + u T sxs p j ds 6

9 The resrcon p p m + guaranees ha p+ p j for every, j and hen p+ p p j, + p + p j are Hölder conjugaes Now he Hölder nequaly s used wh f = u T x pj and g = n he nerval [, + T ] and Hölder conjugaes as descrbed before T p p j + +T u T sxs p+ pj ds +T u T sxs p p + j ds Usng 5 n 6, applyng 7 n he sum, one has u p + h p T p +T + T p j j= u p + T +T From 4, 5 and 8 we oban +T u T sxs p+ pj ds u T sxs p+ ds T ɛp+ p x p + +T u p + h p T p + T p j j= 7 p + ẋs ds 8 u T sxs p+ pj ds Defnng z := +T u T sxs p+ ds he las nequaly can be rewren as z + u p + h p T p + T p zp j j T ɛp + x p p+ j= The polynomal n he LHS s a src monooncally ncreasng funcon and zero when z = 0, hen s nverse exss and also s a src monooncally ncreasng funcon By denong he polynomal as P and s nverse as P one ges +T u T sxs p+ ds P T ɛ p + p x p + 9 Now negrang 5b from o + T and usng 9 he nex nequaly for V s obaned +T V sds = = = +T P u T sxs p + ds T ɛ p + p x p + Ths sasfes he assumpons of he Theorem 5 n [3] and guaranees he global unform asympoc sably of he orgn of 4

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

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

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

CS286.2 Lecture 14: Quantum de Finetti Theorems II

CS286.2 Lecture 14: Quantum de Finetti Theorems II CS286.2 Lecure 14: Quanum de Fne Theorems II Scrbe: Mara Okounkova 1 Saemen of he heorem Recall he las saemen of he quanum de Fne heorem from he prevous lecure. Theorem 1 Quanum de Fne). Le ρ Dens C 2

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

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

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

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

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

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

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

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

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

Lecture 18: The Laplace Transform (See Sections and 14.7 in Boas)

Lecture 18: The Laplace Transform (See Sections and 14.7 in Boas) Lecure 8: The Lalace Transform (See Secons 88- and 47 n Boas) Recall ha our bg-cure goal s he analyss of he dfferenal equaon, ax bx cx F, where we emloy varous exansons for he drvng funcon F deendng on

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

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

Comparison of Differences between Power Means 1

Comparison of Differences between Power Means 1 In. Journal of Mah. Analyss, Vol. 7, 203, no., 5-55 Comparson of Dfferences beween Power Means Chang-An Tan, Guanghua Sh and Fe Zuo College of Mahemacs and Informaon Scence Henan Normal Unversy, 453007,

More information

3. OVERVIEW OF NUMERICAL METHODS

3. OVERVIEW OF NUMERICAL METHODS 3 OVERVIEW OF NUMERICAL METHODS 3 Inroducory remarks Ths chaper summarzes hose numercal echnques whose knowledge s ndspensable for he undersandng of he dfferen dscree elemen mehods: he Newon-Raphson-mehod,

More information

Let s treat the problem of the response of a system to an applied external force. Again,

Let s treat the problem of the response of a system to an applied external force. Again, Page 33 QUANTUM LNEAR RESPONSE FUNCTON Le s rea he problem of he response of a sysem o an appled exernal force. Agan, H() H f () A H + V () Exernal agen acng on nernal varable Hamlonan for equlbrum sysem

More information

Chapter Lagrangian Interpolation

Chapter Lagrangian Interpolation Chaper 5.4 agrangan Inerpolaon Afer readng hs chaper you should be able o:. dere agrangan mehod of nerpolaon. sole problems usng agrangan mehod of nerpolaon and. use agrangan nerpolans o fnd deraes and

More information

Chapter 6: AC Circuits

Chapter 6: AC Circuits Chaper 6: AC Crcus Chaper 6: Oulne Phasors and he AC Seady Sae AC Crcus A sable, lnear crcu operang n he seady sae wh snusodal excaon (.e., snusodal seady sae. Complee response forced response naural response.

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

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

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

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

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

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

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

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

. 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

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

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 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

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

Part II CONTINUOUS TIME STOCHASTIC PROCESSES

Part II CONTINUOUS TIME STOCHASTIC PROCESSES Par II CONTINUOUS TIME STOCHASTIC PROCESSES 4 Chaper 4 For an advanced analyss of he properes of he Wener process, see: Revus D and Yor M: Connuous marngales and Brownan Moon Karazas I and Shreve S E:

More information

Should Exact Index Numbers have Standard Errors? Theory and Application to Asian Growth

Should Exact Index Numbers have Standard Errors? Theory and Application to Asian Growth Should Exac Index umbers have Sandard Errors? Theory and Applcaon o Asan Growh Rober C. Feensra Marshall B. Rensdorf ovember 003 Proof of Proposon APPEDIX () Frs, we wll derve he convenonal Sao-Vara prce

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

e-journal Reliability: Theory& Applications No 2 (Vol.2) Vyacheslav Abramov

e-journal Reliability: Theory& Applications No 2 (Vol.2) Vyacheslav Abramov June 7 e-ournal Relably: Theory& Applcaons No (Vol. CONFIDENCE INTERVALS ASSOCIATED WITH PERFORMANCE ANALYSIS OF SYMMETRIC LARGE CLOSED CLIENT/SERVER COMPUTER NETWORKS Absrac Vyacheslav Abramov School

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

Math 128b Project. Jude Yuen

Math 128b Project. Jude Yuen Mah 8b Proec Jude Yuen . Inroducon Le { Z } be a sequence of observed ndependen vecor varables. If he elemens of Z have a on normal dsrbuon hen { Z } has a mean vecor Z and a varancecovarance marx z. Geomercally

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

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

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

Epistemic Game Theory: Online Appendix

Epistemic Game Theory: Online Appendix Epsemc Game Theory: Onlne Appendx Edde Dekel Lucano Pomao Marcano Snscalch July 18, 2014 Prelmnares Fx a fne ype srucure T I, S, T, β I and a probably µ S T. Le T µ I, S, T µ, βµ I be a ype srucure ha

More information

Tight results for Next Fit and Worst Fit with resource augmentation

Tight results for Next Fit and Worst Fit with resource augmentation Tgh resuls for Nex F and Wors F wh resource augmenaon Joan Boyar Leah Epsen Asaf Levn Asrac I s well known ha he wo smple algorhms for he classc n packng prolem, NF and WF oh have an approxmaon rao of

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

How about the more general "linear" scalar functions of scalars (i.e., a 1st degree polynomial of the following form with a constant term )?

How about the more general linear scalar functions of scalars (i.e., a 1st degree polynomial of the following form with a constant term )? lmcd Lnear ransformaon of a vecor he deas presened here are que general hey go beyond he radonal mar-vecor ype seen n lnear algebra Furhermore, hey do no deal wh bass and are equally vald for any se of

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

Sampling Procedure of the Sum of two Binary Markov Process Realizations

Sampling Procedure of the Sum of two Binary Markov Process Realizations Samplng Procedure of he Sum of wo Bnary Markov Process Realzaons YURY GORITSKIY Dep. of Mahemacal Modelng of Moscow Power Insue (Techncal Unversy), Moscow, RUSSIA, E-mal: gorsky@yandex.ru VLADIMIR KAZAKOV

More information

On the numerical treatment ofthenonlinear partial differentialequation of fractional order

On the numerical treatment ofthenonlinear partial differentialequation of fractional order IOSR Journal of Mahemacs (IOSR-JM) e-iss: 2278-5728, p-iss: 239-765X. Volume 2, Issue 6 Ver. I (ov. - Dec.26), PP 28-37 www.osrjournals.org On he numercal reamen ofhenonlnear paral dfferenalequaon of fraconal

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

Performance Analysis for a Network having Standby Redundant Unit with Waiting in Repair

Performance Analysis for a Network having Standby Redundant Unit with Waiting in Repair TECHNI Inernaonal Journal of Compung Scence Communcaon Technologes VOL.5 NO. July 22 (ISSN 974-3375 erformance nalyss for a Nework havng Sby edundan Un wh ang n epar Jendra Sngh 2 abns orwal 2 Deparmen

More information

Mechanics Physics 151

Mechanics Physics 151 Mechancs Physcs 5 Lecure 0 Canoncal Transformaons (Chaper 9) Wha We Dd Las Tme Hamlon s Prncple n he Hamlonan formalsm Dervaon was smple δi δ Addonal end-pon consrans pq H( q, p, ) d 0 δ q ( ) δq ( ) δ

More information

Method of upper lower solutions for nonlinear system of fractional differential equations and applications

Method of upper lower solutions for nonlinear system of fractional differential equations and applications Malaya Journal of Maemak, Vol. 6, No. 3, 467-472, 218 hps://do.org/1.26637/mjm63/1 Mehod of upper lower soluons for nonlnear sysem of fraconal dfferenal equaons and applcaons D.B. Dhagude1 *, N.B. Jadhav2

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

[ ] 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

Lecture 11 SVM cont

Lecture 11 SVM cont Lecure SVM con. 0 008 Wha we have done so far We have esalshed ha we wan o fnd a lnear decson oundary whose margn s he larges We know how o measure he margn of a lnear decson oundary Tha s: he mnmum geomerc

More information

PHYS 705: Classical Mechanics. Canonical Transformation

PHYS 705: Classical Mechanics. Canonical Transformation PHYS 705: Classcal Mechancs Canoncal Transformaon Canoncal Varables and Hamlonan Formalsm As we have seen, n he Hamlonan Formulaon of Mechancs,, are ndeenden varables n hase sace on eual foong The Hamlon

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

Volatility Interpolation

Volatility Interpolation Volaly Inerpolaon Prelmnary Verson March 00 Jesper Andreasen and Bran Huge Danse Mares, Copenhagen wan.daddy@danseban.com brno@danseban.com Elecronc copy avalable a: hp://ssrn.com/absrac=69497 Inro Local

More information

Testing a new idea to solve the P = NP problem with mathematical induction

Testing a new idea to solve the P = NP problem with mathematical induction Tesng a new dea o solve he P = NP problem wh mahemacal nducon Bacground P and NP are wo classes (ses) of languages n Compuer Scence An open problem s wheher P = NP Ths paper ess a new dea o compare he

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

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

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

Motion in Two Dimensions

Motion in Two Dimensions Phys 1 Chaper 4 Moon n Two Dmensons adzyubenko@csub.edu hp://www.csub.edu/~adzyubenko 005, 014 A. Dzyubenko 004 Brooks/Cole 1 Dsplacemen as a Vecor The poson of an objec s descrbed by s poson ecor, r The

More information

Bayes rule for a classification problem INF Discriminant functions for the normal density. Euclidean distance. Mahalanobis distance

Bayes rule for a classification problem INF Discriminant functions for the normal density. Euclidean distance. Mahalanobis distance INF 43 3.. Repeon Anne Solberg (anne@f.uo.no Bayes rule for a classfcaon problem Suppose we have J, =,...J classes. s he class label for a pxel, and x s he observed feaure vecor. We can use Bayes rule

More information

Survival Analysis and Reliability. A Note on the Mean Residual Life Function of a Parallel System

Survival Analysis and Reliability. A Note on the Mean Residual Life Function of a Parallel System Communcaons n Sascs Theory and Mehods, 34: 475 484, 2005 Copyrgh Taylor & Francs, Inc. ISSN: 0361-0926 prn/1532-415x onlne DOI: 10.1081/STA-200047430 Survval Analyss and Relably A Noe on he Mean Resdual

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

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 DEECIO AD EIMAIO: Fundamenal ssues n dgal communcaons are. Deecon and. Esmaon Deecon heory: I deals wh he desgn and evaluaon of decson makng processor ha observes he receved sgnal and guesses

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

II. Light is a Ray (Geometrical Optics)

II. Light is a Ray (Geometrical Optics) II Lgh s a Ray (Geomercal Opcs) IIB Reflecon and Refracon Hero s Prncple of Leas Dsance Law of Reflecon Hero of Aleandra, who lved n he 2 nd cenury BC, posulaed he followng prncple: Prncple of Leas Dsance:

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

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

On elements with index of the form 2 a 3 b in a parametric family of biquadratic elds

On elements with index of the form 2 a 3 b in a parametric family of biquadratic elds On elemens wh ndex of he form a 3 b n a paramerc famly of bquadrac elds Bora JadrevĆ Absrac In hs paper we gve some resuls abou prmve negral elemens p(c p n he famly of bcyclc bquadrac elds L c = Q ) c;

More information

arxiv: v1 [cs.sy] 2 Sep 2014

arxiv: v1 [cs.sy] 2 Sep 2014 Noname manuscrp No. wll be nsered by he edor Sgnalng for Decenralzed Roung n a Queueng Nework Y Ouyang Demoshens Tenekezs Receved: dae / Acceped: dae arxv:409.0887v [cs.sy] Sep 04 Absrac A dscree-me decenralzed

More information

Lecture VI Regression

Lecture VI Regression Lecure VI Regresson (Lnear Mehods for Regresson) Conens: Lnear Mehods for Regresson Leas Squares, Gauss Markov heorem Recursve Leas Squares Lecure VI: MLSC - Dr. Sehu Vjayakumar Lnear Regresson Model M

More information

Online Appendix for. Strategic safety stocks in supply chains with evolving forecasts

Online Appendix for. Strategic safety stocks in supply chains with evolving forecasts Onlne Appendx for Sraegc safey socs n supply chans wh evolvng forecass Tor Schoenmeyr Sephen C. Graves Opsolar, Inc. 332 Hunwood Avenue Hayward, CA 94544 A. P. Sloan School of Managemen Massachuses Insue

More information

An introduction to Support Vector Machine

An introduction to Support Vector Machine An nroducon o Suppor Vecor Machne 報告者 : 黃立德 References: Smon Haykn, "Neural Neworks: a comprehensve foundaon, second edon, 999, Chaper 2,6 Nello Chrsann, John Shawe-Tayer, An Inroducon o Suppor Vecor Machnes,

More information

Density Matrix Description of NMR BCMB/CHEM 8190

Density Matrix Description of NMR BCMB/CHEM 8190 Densy Marx Descrpon of NMR BCMBCHEM 89 Operaors n Marx Noaon Alernae approach o second order specra: ask abou x magnezaon nsead of energes and ranson probables. If we say wh one bass se, properes vary

More information

Existence of Time Periodic Solutions for the Ginzburg-Landau Equations. model of superconductivity

Existence of Time Periodic Solutions for the Ginzburg-Landau Equations. model of superconductivity Journal of Mahemacal Analyss and Applcaons 3, 3944 999 Arcle ID jmaa.999.683, avalable onlne a hp:www.dealbrary.com on Exsence of me Perodc Soluons for he Gnzburg-Landau Equaons of Superconducvy Bxang

More information

Dual Approximate Dynamic Programming for Large Scale Hydro Valleys

Dual Approximate Dynamic Programming for Large Scale Hydro Valleys Dual Approxmae Dynamc Programmng for Large Scale Hydro Valleys Perre Carpener and Jean-Phlppe Chanceler 1 ENSTA ParsTech and ENPC ParsTech CMM Workshop, January 2016 1 Jon work wh J.-C. Alas, suppored

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

Discrete Time Approximation and Monte-Carlo Simulation of Backward Stochastic Differential Equations

Discrete Time Approximation and Monte-Carlo Simulation of Backward Stochastic Differential Equations Dscree Tme Approxmaon and Mone-Carlo Smulaon of Backward Sochasc Dfferenal Equaons Bruno Bouchard Unversé Pars VI, PMA, and CREST Pars, France bouchard@ccrjusseufr Nzar Touz CREST Pars, France ouz@ensaefr

More information

@FMI c Kyung Moon Sa Co.

@FMI c Kyung Moon Sa Co. Annals of Fuzzy Mahemacs and Informacs Volume 8, No. 2, (Augus 2014), pp. 245 257 ISSN: 2093 9310 (prn verson) ISSN: 2287 6235 (elecronc verson) hp://www.afm.or.kr @FMI c Kyung Moon Sa Co. hp://www.kyungmoon.com

More information

This document is downloaded from DR-NTU, Nanyang Technological University Library, Singapore.

This document is downloaded from DR-NTU, Nanyang Technological University Library, Singapore. Ths documen s downloaded from DR-NTU, Nanyang Technologcal Unversy Lbrary, Sngapore. Tle A smplfed verb machng algorhm for word paron n vsual speech processng( Acceped verson ) Auhor(s) Foo, Say We; Yong,

More information

A New Generalized Gronwall-Bellman Type Inequality

A New Generalized Gronwall-Bellman Type Inequality 22 Inernaonal Conference on Image, Vson and Comung (ICIVC 22) IPCSIT vol. 5 (22) (22) IACSIT Press, Sngaore DOI:.7763/IPCSIT.22.V5.46 A New Generalzed Gronwall-Bellman Tye Ineualy Qnghua Feng School of

More information

Genetic Algorithm in Parameter Estimation of Nonlinear Dynamic Systems

Genetic Algorithm in Parameter Estimation of Nonlinear Dynamic Systems Genec Algorhm n Parameer Esmaon of Nonlnear Dynamc Sysems E. Paeraks manos@egnaa.ee.auh.gr V. Perds perds@vergna.eng.auh.gr Ah. ehagas kehagas@egnaa.ee.auh.gr hp://skron.conrol.ee.auh.gr/kehagas/ndex.hm

More information

Handout # 6 (MEEN 617) Numerical Integration to Find Time Response of SDOF mechanical system Y X (2) and write EOM (1) as two first-order Eqs.

Handout # 6 (MEEN 617) Numerical Integration to Find Time Response of SDOF mechanical system Y X (2) and write EOM (1) as two first-order Eqs. Handou # 6 (MEEN 67) Numercal Inegraon o Fnd Tme Response of SDOF mechancal sysem Sae Space Mehod The EOM for a lnear sysem s M X DX K X F() () X X X X V wh nal condons, a 0 0 ; 0 Defne he followng varables,

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

Advanced Machine Learning & Perception

Advanced Machine Learning & Perception Advanced Machne Learnng & Percepon Insrucor: Tony Jebara SVM Feaure & Kernel Selecon SVM Eensons Feaure Selecon (Flerng and Wrappng) SVM Feaure Selecon SVM Kernel Selecon SVM Eensons Classfcaon Feaure/Kernel

More information

Chapter 2 Linear dynamic analysis of a structural system

Chapter 2 Linear dynamic analysis of a structural system Chaper Lnear dynamc analyss of a srucural sysem. Dynamc equlbrum he dynamc equlbrum analyss of a srucure s he mos general case ha can be suded as akes no accoun all he forces acng on. When he exernal loads

More information

Dynamic Team Decision Theory

Dynamic Team Decision Theory Dynamc Team Decson Theory EECS 558 Proec Repor Shruvandana Sharma and Davd Shuman December, 005 I. Inroducon Whle he sochasc conrol problem feaures one decson maker acng over me, many complex conrolled

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

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

( ) [ ] 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

Supporting Information: The integrated Global Temperature change Potential (igtp) and relationships between emission metrics

Supporting Information: The integrated Global Temperature change Potential (igtp) and relationships between emission metrics 2 3 4 5 6 7 8 9 Supporng Informaon: Te negraed Global Temperaure cange Poenal (GTP) and relaonsps beween emsson mercs Glen P. Peers *, Borgar Aamaas, Tere Bernsen,2, Jan S. Fuglesved Cener for Inernaonal

More information

Appendix to Online Clustering with Experts

Appendix to Online Clustering with Experts A Appendx o Onlne Cluserng wh Expers Furher dscusson of expermens. Here we furher dscuss expermenal resuls repored n he paper. Ineresngly, we observe ha OCE (and n parcular Learn- ) racks he bes exper

More information

Lecture 2 M/G/1 queues. M/G/1-queue

Lecture 2 M/G/1 queues. M/G/1-queue Lecure M/G/ queues M/G/-queue Posson arrval process Arbrary servce me dsrbuon Sngle server To deermne he sae of he sysem a me, we mus now The number of cusomers n he sysems N() Tme ha he cusomer currenly

More information

Delay-Range-Dependent Stability Analysis for Continuous Linear System with Interval Delay

Delay-Range-Dependent Stability Analysis for Continuous Linear System with Interval Delay Inernaonal Journal of Emergng Engneerng esearch an echnology Volume 3, Issue 8, Augus 05, PP 70-76 ISSN 349-4395 (Prn) & ISSN 349-4409 (Onlne) Delay-ange-Depenen Sably Analyss for Connuous Lnear Sysem

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