Fundamental Limitations of Discrete-Time Adaptive Nonlinear Control

Size: px
Start display at page:

Download "Fundamental Limitations of Discrete-Time Adaptive Nonlinear Control"

Transcription

1 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 44, NO. 9, SEPTEMBER In hs case, neher of he roos of (3) les on he agnary axs. For k > 0 and > 0, because of he second exresson n (3) a leas one of wo egenvalues les n he lef half-lane. Fro he frs exresson n (3), he su of he arguens of and has o be. Thus oher egenvalues us le n he lef half-lane. Hence he syse s sable. Thus an addon of dssave force o he syse n (9) does no desablze he syse. For k = 0(sac nsably), one of he wo egenvalues wll be a he orgn. For k < 0 (wh > 0), he su of he arguens of and s. Hence oher egenvalues us le n he rgh-hand lane n order o sasfy he exressons n (3). Thus (7) sablzed by gyroscoc forces (9) for k < 0 can be agan desablzed by he addon of dssave forces (3). For exale, he o s sablzed by a gyroscoc oen as long as he sn s suffcenly large. Evenually he sn s decreased by frcon (dssave load) o ake he o unsable. V. CONCLUSIONS In hs aer, sably ssues of arx second-order dynacal syses are dscussed. The necessary and suffcen condons of asyoc sably for e-nvaran syses n arx second-order for under varous dynac loadngs (conservave/nonconservave) are derved and a hyscal nerreaon s resened. The sably condons n he sense of Lyaunov are also derved and analyzed. As he condons are drec n ers of hyscal araeers of he syse, he effec of dfferen loadngs on he syse sably s lucd n he arx second-order for aroach. The condons are shown o be useful n he desgnng conrollers for arx second-order syses. REFERENCES [] R. E. Skelon, Adave orhogonal flers for coensaon of odel errors n arx second-order syses, J. Gudance, Conr. Dyna.,. 4, Mar [] J.-N. Juang and P. G. Magha, Robus egensyse assgnen for second order dynac syses, n Proc. AIAA Srucural Dynacs Secals Conf., 99, [3] A. M. Dwekar and R. K. Yedavall, Sar srucure conrol n arx second order for, n Proc. Norh Aercan Conf. Sar Srucures and Maerals, 995, [4] L. S. Sheh, C. D. Shh, and R. E. Yaes, Soe suffcen and soe necessary condons for he sably of ulvarable syses, Trans. ASME, vol. 00,. 4 8, Se [5] S. K. Shrvasava and S. Pradee, Sably of uldensonal lnear e-varyng syses, J. Gudance, Conr., Dyna., , Oc [6] P. Hsu and J. Wu, Sably of second-order uldensonal lnear e-varyng syses, J. Gudance, Conr. Dyna., vol. 4, , Se. 99. [7] L. S. Sheh, M. M. Meho, and H. M. Db, Sably of he second order arx olynoals, IEEE Trans. Auoa. Conr., vol. 3,. 3 33, Mar [8] F. R. Ganacher, Alcaons of Theory of Marces. New York: Wley, 959. [9] P. Lancaser, Labda Marces and Vbrang Syses. New York: Pergaon, 964. [0] A. M. Dwekar and R. K. Yedavall, Aeroelasc analyss hrough arx second order for, n Proc. AIAA Gudance, Navgaon and Conrol Conf., Aug. 995, [] H. Zegler, Prncles of Srucural Sably. London, U.K.: Blasdell, 966. [] R. Bellan, Inroducon o Marx Analyss. New York: McGraw-Hll, 970. [3] R. A. Horn and C. R. Johnson, Marx Analyss. New York: Cabrdge Unv. Press, 990. [4] E. L. Dowell and H. C. Curss, Jr., A Modern Course n Aeroelascy. London, U.K.: Kluwer, 989. [5] Z. P. Bazan and L. Cedoln, Sably of Srucures. New York: Oxford Unv., 99. Fundaenal Laons of Dscree-Te Adave Nonlnear Conrol Lang-Lang Xe and Le Guo Absrac A arcular olynoal s nroduced n hs aer whch can be used o deerne under wha condons a ycal class of dscreee nonlnear syses wh unceranes n boh araeers and noses s no sablzable by feedback, hus deonsrang he fundaenal laons of dscree-e adave nonlnear conrol. As a consequence, s shown ha for nonlnear syses wh unknown araeers and noses, he syses ay ndeed be nonsablzable, n general, whenever he usual lnear growh condon s relaxed and he nuber of unknown araeers s large, even hough he corresondng nose-free syses are globally sablzable. Index Ters Adave conrol, dscree-e, global sablzably, nonlnear dynacs, sochasc syses. I. INTRODUCTION Adave lnear conrol and he relaed ssues have been he an focus of adave conrol over he as several decades (see, e.g., [] [3], [6], [8], [9], and []). In recen years, aes have been ade oward a heory of adave nonlnear conrol. If he nonlneary s only nvolved n he nu ar, or f he ouu ar of a syse s nonlnear bu has a lnear growh rae, hen s farly well known ha he exsng adave conrol ehods can sll be aled as long as he unknown araeers ener he syse lnearly, wheher he syse s descrbed n connuous-e or dscree-e (see, e.g., [3], [5], and [6]). However, he suaon changes draacally when one aes o deal wh syses wh ouu nonlneares havng growh raes faser han lnear. Neher of he exsng ehods are useful, nor do he slares beween adave conrol of connuous- and dscree-e syses rean. For a large class of connuous-e nonlnear syses, nonlnear-dang and/or back-seng aroaches can be successfully used n adave conrol desgn regardless of he growh rae of he nonlneares (cf., e.g., []). Ths s so even n he case where exernal dsurbances exs (c.f. [4] and [4]). For exale, consder he followng connuouse sochasc conrol odel: dy =(y b + u )d + dw; b > 0 () where s an unknown araeer and w s a sandard Brownan oon, and y and u are he syse ouu and nu sgnals, resecvely. Then can be shown easly by usng he Io forula ha () can be a.s. globally sablzed by he nonlnear dang conrol u = 0y 0 y jy b j for any b > 0. Manuscr receved May 6, 998. Recoended by Assocae Edor, M. Krsc. Ths work was suored by he Naonal Naural Scence Foundaon of Chna and he Naonal Key Projec of Chna. The auhors are wh he Insue of Syses Scence, Chnese Acadey of Scences, Bejng, 00080, Chna (e-al: Lguo@ss03.ss.ac.cn). Publsher Ie Idenfer S (99)0656-X /99$ IEEE

2 778 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 44, NO. 9, SEPTEMBER 999 However, he connuous-e aroaches do no work n he dscree-e case when he nonlnear funcon has a growh rae faser han lnear, as observed by several auhors (see, e.g., [0] and [6]). Thus a queson naurally arses: Can we fnd a sablzng adave conroller n hs case? A generally negave answer was recenly gven n [7], where a crcal sablzably henoenon was found for a class of nonlnear conrol syses. To be recse, for he followng ycal conrol odel: y + = y b + u + w + ; b > 0 () where s an unknown araeer and fw g s a Gaussan whe nose sequence, has been shown n [7] ha () s no a.s. globally sablzable f and only f b 4, or he followng nequaly has a soluon: z 0 bz + b 0; z (; b): (3) The above resul clearly deonsraes he laons of adave conrol n he dscree-e case and shows ha he dscree-e robles are uch ore colcaed. Ths aer s anly concerned wh dscree-e syses. We shall sudy nonlnear odels wh ul-unknown araeers, whch are exensons of he scalar araeer case () as suded n [7]. Corresondng o (3), we shall nroduce a generalzed olynoal P (z) n he ularaeer case, whch wll be used o deerne when a nonlnear syse s no sablzable by feedback, hus deonsrang he fundaenal laons of adave nonlnear conrol n he dscree-e case. The boo lne of our an resul s soewha unexeced whch shows ha for dscree-e nonlnear syses wh rando unceranes n boh araeers and dsurbances, global adave sablzaon s ossble, n general, whou osng he lnear growh condon on he nonlneares of he syses. II. MAIN RESULTS Consder he followng dscree-e olynoal nonlnear regresson odel: y + = y b + y b + + n y b + u + w + ; 0 (4) where y and u are he syse ouu and nu sgnals, resecvely, ( n) are unknown araeers, and w s he nose sgnal. Assue he followng. A) b ( n) are nonnegave real nubers akng (4) eanngful and sasfyng b >b > >b n > 0: A) fw g s a Gaussan whe nose sequence wh dsrbuon N (0; ). A3) The unknown araeer vecor =[ ; ; n] s ndeenden of fw g and has a Gaussan dsrbuon N (; I n). Our objecve s o sudy he global sablzably of (4) under he above condons. Frs, we gve a recse defnon of sablzably. Defnon : Le fy ; 0 g be he feld generaed by he observaons y ; 0. Syse (4) s sad o be a.s. globally sablzable, f here exss a feedback conrol u F y = fy ; 0 g; =0; ; ; such ha for any nal condon y 0 R, l su T! =T T = y < ; a.s. Reark : We reark ha he global sablzaon of (4) s a rval ask n eher he case where s known or he case where he nose s free (.e., w 0). To be recse, f were known, we can u u 0( y b + y b + + ny b ); whch obvously globally sablzes he syse. In he case where s unknown bu he nose s free (w 0), we can oban he rue value of he araeer by solvng n ndeenden lnear equaons. For exale, f n he frs (n +)ses we choose fu ; 0 ng o be ndeendenly dencally dsrbued rando varables wh robably densy funcon (x), hen he rue value of he araeer can be obaned easly by solvng he followng lnear equaon: A =[y 0 u ;y 3 0 u ; ;y n+ 0 u n ] where A =(y b ) nn s a nonsngular arx (cf. [5]). Hence, agan, we can ake he conrol as u = 0( y b + y b + + n y b ) for > n, whch globally sablzes he nose-free syse. For ore general araerc-src-feedback odels wh no nose, a relaed bu ore colcaed wo-hase aroach can also be aled o desgn a globally sablzng adave conroller regardless of he growh rae of he nonlneares (cf., [7]). Unforunaely, he an drawback of he above wo-hase aroach s ha he resulng adave conroller s no robus wh resec o nose. In fac, he resence of nose wll even change he sablzably of dscree-e nonlnear syses draacally f he growh rae of he nonlneares s faser han lnear, as wll be shown by he followng heore ogeher wh s corollares. Theore : Under Assuons A) A3), syse (4) s no a.s. globally sablzable whenever he followng nequaly: P (z) < 0; z (; b ) (5) has a soluon, where P (z) s a olynoal defned by P (z) =z n+ 0b z n +(b 0b )z n0 ++(b n00b n )z+b n : (6) The roof s gven n he nex secon. Reark : Obvously, for n =, P (z) concdes wh he quadrac olynoal n (3). Noe ha a rval necessary condon for (5) o have a soluon s b >, and when b (4) s always a.s. globally sablzable (see [5]). To undersand he lcaons of Theore, we now gve soe dealed dscussons on (5). Corollary : If b ( n) sasfes b > and 0 < b 0b + b =( b 0) ; n0; hen (5) has a soluon whenever n log(( b +)=( b 0 ))= log b. Consequenly, whenever b > and he nuber of unknown araeers n s suably large, here always exs 0 <b n <b n0 < <b such ha (4) s no a.s. globally sablzable. The roof of hs corollary s gven n he nex secon. Reark 3: By Corollary we know ha he usual lnear growh condon osed on he nonlnear funcon f () of he general conrol odel y + = f (y ; ;y 0) +u + w + ; R n (7) canno be essenally relaxed n general for global adave sablzaon, unless addonal condons on he nuber n and he srucure of f () are osed. Corollary : Le b >, hen for n > + log(=(b 0 ))= log(b =), (5) has a soluon for any fb g sasfyng b n < b n0 < <b <b. On he oher hand, f b, hen for any n, here always exs b n <b n0 < <b <b such ha (5) has no soluon. Corollary 3: For any n and any b >b > >b n > 0, we have he followng. n ) A necessary condon for (5) o have a soluon s = b > 4: ) A suffcen condon for (5) o have a soluon s eher b > 4, n or = b > (n + )( + (=n)) n : The roofs of Corollares and 3 are n [5], due o sace laons. III. PROOF OF THE MAIN RESULTS We frs resen he roof of Theore, whch s refaced wh wo leas.

3 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 44, NO. 9, SEPTEMBER Lea : Le fc ; g and fd ; qg be wo sequences of osve nubers sasfyng c > c > > c ; d > d > > d q > 0; q and c c + + ; 0 ; for soe >0: Also, le f ; ; g and fj ; ;j g be wo sequences of negers arbrarly aken fro f; ; ; qg, wh k 6= l and j k 6= j l, for k 6= l. If here exss an neger 0( 0 ) such ha eher 6= 0 or j 6= 0, hen c d +d c d +d c d +d c c d c d n(;) n (d 0d ) c d where by defnon d q+ = 0. Proof: By he assuons d >d > >d q > 0 and k 6= l ;j k 6= j l for k 6= l, s easy o see ha k= d k 0 (d + d j ) 0; 8 [; ]: (8) k= Moreover, whou loss of generaly assue ha 0 =nff: 6= or j 6= ; g: Then k= d k 0 (d + d j ) k= =d 0 (d + d j ) d 0 d + > 0: (9) Now, by he assuon c c + + ;c and (8) and (9), we have c d c d c d c d +d c d +d c d +d d 0(d +d ) d 0(d +d ) d 0(d = c c c +d ) (d 0d 0d ) (d 0d 0d ) d 0d c c c 0d 0 (d 0d 0d ) c + c = 0 c (d 0d 0d ) c = n(;) c (d 0d 0d ) = n(;) (d 0d 0d ) c n(;)(d 0d ) c n(;) n (d 0d ) c : (d 0d 0d ) (d 0d 0d ) Hence Lea s rue. The followng lea lays a key role n he roof of Theore. Lea : Assue ha for soe > 0 and, jy j jy 0j + ;= ; ; ;and ha he nal condon jy 0j s suffcenly large, hen he deernan of he arx P 0 + sasfes 0 0n+ jp+j n+ (0) where by defnon y =for <0, and P+ s defned recursvely by wh P + = P 0 P'' P +' P ' ; P 0 = I () ' =[y b ;y b ; ;y b ] ; b >b > >b n > 0: () Snce he roof of hs lea s raher nvolved, we gve n he Aendx. Proof of Theore : We only need o rove ha f (5) has a soluon, hen for any feedback conrol u F y, here always exss an nal condon y 0 and a se D 0 wh osve robably such ha he ouu sgnal y of he closed-loo conrol syse ends o nfny a a rae faser han exonenal on D 0. Now consder he followng sae sace equaon ( 0): + = ; 0 = ; y + = ' + u + w + (3) where s he unknown araeer vecor defned n Assuon A3) and ' s defned by (). By our Assuons A) and A3) and he fac ha u F y, we know ha (3) s a condonal Gaussan odel, and hence he condonal execaon ^ = E[jF y ] can be generaed by he Kalan fler and he condonal covarance arx of he esaon error ( 0 ^ ) can be generaed by he Rcca equaon () or where ~ = 0 ^ (see, e.g., [3, Sec. 3.]). Nex, by (3) we know ha E[ ~ ~ jf y ]=P (4) y + = ' ~ +(' ^ + u )+w + : (5) Consequenly, by he fac ha E[ ~ jf y ]=0and E[w +jf y ]=0 follows fro (4) and (5) ha for any u F y E[y +jf y ]=' P ' +(' ^ + u ) + ' P ' +: (6) Furherore, by he arx nverson forula follows fro () ha P 0 + = P 0 + ' ' : Then jp+j 0 = jp 0 + ' ' j = jp 0 (I + P ' ' )j = jp 0 j( + ' P ' ): Hence, follows fro hs and (6) ha Le us defne D 0 = E[y+jF y 0 jp+ ] j jp 0 ; 8 0: (7) j =0 f!: E[y +jf y ] ( +)5= y +g: Then by he condonal Gaussan roery of he sequence fy g,a coleely slar arguen as ha used n [7, Aendx B] shows ha Prob(D 0 ) > 0. Hence, by (7) we have y+ jp +j 0 ( +) 5= jp 0 j ; 0; on D 0: (8) Now, le z 0 (; b ) be a soluon of (5). We roceed o rove ha on D 0 jy jjy 0j z ; =; ; : (9)

4 780 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 44, NO. 9, SEPTEMBER 999 We ado he nducon arguen. Frs, we consder he case where =. Snce P 0 0 = I, we have by () jp 0 j = jp ' 0 ' 0j =+k' 0 k > 0 : Therefore, by (8) we have for jy 0 j, jy j jp 0 j=jp 0 0 j > jy 0j b > jy 0j z,ond 0. Hence (9) s rue for =. Now le us assue ha for soe, jy jjy 0j z ; = ; ; ;; on D 0 ; hen by Lea, follows ha and 0n+ jp 0 +j 3 0 y0n b jp 0 Consequenly, by (8) we have on D 0 jy +j = = j 3 0 jp +j ( +) 5= jp 0 j 3( +) 5= 0n+ 0 y0n: b = 0n+ 0 0n = 3( jyjb jy 0n+j b +) 5=4 jy 0j b jy 0nj b 3( +) 5=4 jy j b jy 0j b 0b jy 0n+j b 0b jy 0nj 0b : (0) However, by he nducon assuon we have jy 0j jy j z ; ; and so by b + 0 b < 0( n) jy 0j b 0b jy j (b 0b )z ; n: () Noe ha hs nequaly also holds for >, snce y j =for j<0 by defnon. Hence, follows fro (0) and () ha on D 0 jy +j 3( +) 5=4 b +(b 0b )z ++(b 0b )z 0b z jy j 3( +) 5=4 jyj0z P (z ) jy j z [jy 0j 0z P (z ) ] z jy j z jy j z 3( +) 5=4 where he las nequaly holds for suffcenly large jy 0 j because 0z 0n 0 P (z 0 ) > 0. Hence, by nducon, (9) s rue. Thus for all large nal condons jy 0 j, he ouu rocess jy j dverges o nfny a a rae faser han exonenal on D 0. Ths colees he roof of Theore. Proof of Corollary : Take z 0 = b (; b ). We need only verfy ha P (z 0 ) < 0. By he assuon for n, we have z n 0 = ( b ) n ( b + )=( b 0 ): Hence by (6) and he condons on b, we have P (z 0 ) <z n 0 (z 0 0 b )+ b b 0 (z n )+b = 0 z n 0 b b 0 + b z0 n 0 b 0 z b = 0 z n b 0 b b 0 + b 0 z n 0 0 b b 0 + b b = 0 b 0 z n b 0 + b + b b + 0 b 0 b 0 0 b + =0: Hence, Corollary holds. IV. CONCLUDING REMARKS I s farly well known ha for nonlnear sochasc syses descrbed by nonlnear regresson odels wh lnear unknown araeers, a globally sablzng adave conroller can be desgned whenever he nonlnear funcon [say f (x)] nvolved has a lnear growh rae,.e., jf (x)j = O(jxj); as jxj!. However, n conras o he connuous-e case, essenal dffcules eerge n he dscree-e case when he nonlnear funcon has a growh rae faser han lnear. In fac, he nonlnear growh rae has been he crux n dscree-e adave nonlnear conrol for years. Naurally, one would ask he followng quesons: ) Can we reove he usual lnear growh condon n he dscree-e case? ) How far can we go fro lnear growh o nonlnear growh for global sablzaon? A frs se n hs drecon was recenly ade n [7], where was shown ha n he unknown scalar araeer case (n =), he nonlnear conrol syse n queson s globally sablzable f and only f jf (x)j = O(jxj b ) wh b<4. In he resen aer, we have deal wh he general ularaeer case (n ) by consderng he olynoal regresson odel descrbed by (4). By nroducng a new and ore general olynoal (6), we have found a creron abou suaons where (4) s no globally sablzable (Theore ). Based on ha, varous exlc cases are dscussed n Corollares 3. Perhas he os rearkable consequence of our an resul s he followng lcaon for general nonlnear regresson odels. I s ossble n general o essenally relax he usual lnear growh condon for global sablzaon, unless addonal condons are osed (see Reark 3). APPENDIX Proof of Lea : By he arx nverson forula follows fro () ha P 0 + = P 0 + ' ' ; hence, we have (a) shown a he boo of he nex age. b +b Now, le us denoe () = [y ; ; ; ], n; 0, and le (0) = e,.e., he h colun of he deny arx I, hen jp 0 +j = de =0 (); =0 (); ; =0 By he eleenary roeres of deernans, we have jp 0 +j = ; ; ; =0 n() : de( ( ); ( ); ; n(n)): ()

5 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 44, NO. 9, SEPTEMBER I s clear ha f n he grou ( ; ; ;n) here are a leas wo negers havng he sae value (bu dfferen fro 0), hen de( ( ); ;n(n)) = 0. So n he dscussons below we wll exclude hs knd of zero-valued deernan. We roceed o rove (0) by consderng wo cases searaely. Case I) <n0: In hs case, n order ha de( ( ); ( ); ;n(n)) 6= 0, he nuber of (0) s n ( ; ; ;n) us a leas be (n 0 0 ) and he oher negers us be dsnc. Then each er n he exanson of he nonzero deernan de( ( ); ( ); ;n(n)) conans a os [n0(n00)] = ( +)dfferen facors, whose general for s ; + (3) where 6= l ;j 6= jl, and k 6= kl for 6= l. Noe ha one such er s y b 0 0 (fro he roducs of he an dagonal eleens of he arx [ (); ; +(0); +(0); ; n(0)] ), and s dfferen fro oher ers. Now, we roceed o rove ha he absolue value of any oher er s no greaer han =[jy 0 j n(;)n (b 0b ) ] 0 0 : We dvde hs roof no wo subcases. Subcase ): If n he general for (3) = +, hen we can rewre (3) as 0 0. By our assuons, s easy o see ha Lea s alcable, hence we have 0 0 n(;)n (b 0b ) jy 0j 0 0 : Subcase ): If n he general for (3) < +, hen we can add soe ( +; ; +) o ( ; ; ) so ha ( ; ; ; + ; ; + ) T (; 0 ; ; 0): Here and hereafer, we use T (; 0 ; ; 0) o denoe he class of all eruaons of he neger sequence (; 0 ; ; 0). Then by + <n, we can choose suable j + ; ;j + ;k + ; ;k + f; ; ;ngwh j 6= jl, k 6= kl for 6= l; ; l + such ha he er s dfferen fro y b Subcase ), we have 0 0. Thus by he concluson of n(;)n (b 0b ) jy 0j Hence, he desred nequaly s rue. 0 : Now, rewre () as jp 0 + j = R + y b 0 0 ; where R denoes he suaon of all he ers dfferen fro 0 0. I s obvous ha R has a os [( +) n n! 0 ] (n n n! 0 ) ers. Hence, by he resuls roved above, we oban n n n! 0 jr j n(;) n (b 0b ) jy 0 j 0 : Therefore, by choosng he nal value jy 0j large enough, we can ake jr j less han b (=)y 0 0. Consequenly (0) follows. Case II) n 0 : Frs of all, any nonzero deernan de( ( ); ( ); ; n(n)); 0 ; ; n can be exanded as he suaon of n! ers whose general for s ; where (j ; ;j n) T (; ;n); k f0; 0; ;g; k n, and as noed before, any wo k s canno have he sae value dfferen fro (0). Obvously one such er s y b 0 0n+ [one er n de( (); ; n ( 0 n + ))]. We now show ha for any oher ers ( k 6= k or j k 6= k; for soe k n), he followng nequaly holds: j j (=(+)) n (b 0b ) jy 0n+j 0n+: (4) We roceed o rove hs by consderng hree subcases. Subcase ): If ( ; ; ; n ) T (; 0; ;0n +), hen we can drecly aly Lea o ge j j n(;) n (b 0b ) jy 0n+ j =(+) n (b 0b ) jy 0n+ j 0n+ 0n+: Subcase ): If n he grou ( ; ; ; n) here are a leas wo negers less han ( 0 n +), hen we ay fnd ( 0 ; 0 0 ; ; n ) 0 T (; 0 ; ;0n +) such ha l l ( l n) and s dfferen fro y b Consequenly, by he concluson of Subcase ), we have (=(+)) n (b 0b ) jy 0n+ j 0 0n+ : 0n+: jp 0 + j = P = ' ' = + =0 =0 =0 =0 +. =0 =0. =0 = =0 (a)

6 78 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 44, NO. 9, SEPTEMBER 999 Subcase 3): If n he grou ( ; ; ;n) here s jus one neger less han ( 0 n +),.e., 0 n; for soe n, hen by jy 0nj jy 0n+j =(+) we have j j 0n +b 0n+ yb (=(+))(b +b ) jy 0n+j 0 (=(+))(b +b ) 0n+ jy 0n+j 0 0n+ jy 0n+j (=(+))b (=(+)) n (b 0b ) jy 0n+j 0n+ : Hence (4) s rue. Thus, slar o he arguens n Case I), we rewre () as jp 0 + j = R + 0 0n+ ; where R denoes he suaon of [( +) n n! 0 ] ers n he deernan exansons, whch are dfferen fro y0n+. b Then by (4) we know ha jr j ( +) n n! 0 (=(+)) n (b 0b ) jy 0n+j ( +) n n! 0 (+) (=(+))n (b 0b ) jy 0 j 0 0n+ 0 0n+ 0 0n+ where he las nequaly holds for suffcenly large jy 0 j. Hence, he roof of Lea s coleed. ACKNOWLEDGMENT The auhors would lke o hank Prof. P. R. Kuar for hs helful suggesons. REFERENCES [] K. J. Åsrö and B. Wenark, Adave Conrol, nd ed. New York: Addson-Wesley, 995. [] P. E. Canes, Lnear Sochasc Syses. New York: Wley, 988. [3] H. F. Chen and L. Guo, Idenfcaon and Sochasc Adave Conrol. Boson, MA: Brkhäuser, 99. [4] H. Deng and M. Krsć, Sochasc nonlnear sablzaon I: A backseng desgn, Sys. Conr. Le., vol. 3, no. 3, , 997. [5] L. L. Xe and L. Guo, Fundaenal laons of adave conrol, Ins. Syses Scence, Chnese Acadey of Scences, Tech. Re., 998. [6] L. Guo, Self-convergence of weghed leas-squares wh alcaons o sochasc adave conrol, IEEE Trans. Auoa. Conr., vol. 4, , Jan [7], On crcal sably of dscree-e adave nonlnear conrol, IEEE Trans. Auoa. Conr., vol. 4, , Nov [8] L. Guo and H. F. Chen, The Åsrö Wenark self-unng regulaor revsed and ELS-based adave rackers, IEEE Trans. Auoa. Conr., vol. 36,. 80 8, July 99. [9] P. A. Ioannou and J. Sun, Robus Adave Conrol. Englewood Clffs, NJ: Prence-Hall, 996. [0] I. Kanellakooulos, A dscree-e adave nonlnear syse, IEEE Trans. Auoa. Conr., vol. 39,. 36 3, Nov [] P. R. Kuar and P. Varaya, Sochasc Syses: Esaon, Idenfcaon and Adave Conrol. Englewood Clffs, NJ: Prence-Hall, 986. [] M. Krsć, I. Kanellakooulos, and P. V. Kokoovć, Nonlnear and Adave Conrol Desgn. New York: Wley, 995. [3] W. Ln and J.-M. Yong, Drec adave conrol of a class of MIMO nonlnear syses, In. J. Conr., vol. 56, no. 5,. 03 0, 99. [4] Z. Pan and T. Basar, Adave conroller desgn for rackng and dsurbance aenuaon n araerc-src-feedback nonlnear syses, IEEE Trans. Auoa. Conr., vol. 43, , Aug [5] C. We and L. Guo, Adave conrol of a class of nonlnear sochasc syses, n Proc. 996 Chnese Conrol Conf., Qngdao, Se. 5 0, 996, [6] P.-C. Yeh and P. V. Kokoovć, Adave conrol of a class of nonlnear dscree-e syses, In. J. Conr., vol. 6, no., , 995. [7] J.-X. Zhao and I. Kanellakooulos, Adave conrol of dscree-e src-feedback nonlnear syses, n Proc. IEEE Conf. Decson and Conrol, San Dego, CA, Dec Coens on he Couaon of Inerval Rouh Aroxans Chy Hwang and Shh-Feng Yang Absrac In recen aers [6], [7], he Rouh aroxaon ehod was exended o derve reduced-order nerval odels for lnear nerval syses. In hs aer, he auhors show ha: ) nerval Rouh aroxans o a hgh-order nerval ransfer funcon deend on he leenaon of nerval Rouh exanson and nverson algorhs; ) nerval Rouh exanson algorhs canno guaranee he success n generang a full nerval Rouh array; 3) soe nerval Rouh aroxans ay no be robusly sable even f he orgnal nerval syse s robusly sable; and 4) an nerval Rouh aroxan s n general no useful for robus conroller desgn because s dynac unceranes (n ers of robus frequency resonses) do no cover hose of he orgnal nerval syse. Index Ters Inerval syses, odel reducon, Rouh aroxaon. I. INTRODUCTION In he las wo decades, he Rouh aroxaon ehod oneered by Huon and Fredland [] and s varans [] [5] have been recevng uch aenon n he feld of odel reducon. The ehod s based on usng he Rouh sably array o derve reduced-order odels for hgh-order lnear syses. The an advanages of he Rouh aroxaon ehod are ha has he ably o yeld sable reduced-order odels for sable orgnal hgh-order syses, o roduce a faly of reduced-order odels of dfferen orders va a sngle se of algebrac couaons, and o oban reduced odels ha rean he frs several e-oens and/or Markov araeers of he orgnal syses. Recenly, he Rouh aroxaon ehod has been exended o derve reduced-order nerval odels for hgh-order nerval ransfer funcons [6], [7]. The exenson s based on usng nerval arhec o erfor Rouh or canoncal connued-fracon exanson and nverson. Surrsngly, s observed fro he leraure [7] Manuscr receved Deceber, 997. Recoended by Assocae Edor, A. Varga. Ths work was suored by he Naonal Scence Councl of he Reublc of Chna under Gran NSC86-4-E C. Hwang s wh he Dearen of Checal Engneerng, Naonal Chung Cheng Unversy, Cha-Y 6, Tawan (e-al: chch@ccunx.ccu.edu.w). S.-F. Yang s wh he Dearen of Checal Engneerng, Naonal Cheng Kung Unversy, Tanan 70, Tawan. Publsher Ie Idenfer S (99) /99$ IEEE

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

EP2200 Queuing theory and teletraffic systems. 3rd lecture Markov chains Birth-death process - Poisson process. Viktoria Fodor KTH EES

EP2200 Queuing theory and teletraffic systems. 3rd lecture Markov chains Birth-death process - Poisson process. Viktoria Fodor KTH EES EP Queung heory and eleraffc sysems 3rd lecure Marov chans Brh-deah rocess - Posson rocess Vora Fodor KTH EES Oulne for oday Marov rocesses Connuous-me Marov-chans Grah and marx reresenaon Transen and

More information

Normal Random Variable and its discriminant functions

Normal Random Variable and its discriminant functions Noral Rando Varable and s dscrnan funcons Oulne Noral Rando Varable Properes Dscrnan funcons Why Noral Rando Varables? Analycally racable Works well when observaon coes for a corruped snle prooype 3 The

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

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

グラフィカルモデルによる推論 確率伝搬法 (2) Kenji Fukumizu The Institute of Statistical Mathematics 計算推論科学概論 II (2010 年度, 後期 )

グラフィカルモデルによる推論 確率伝搬法 (2) Kenji Fukumizu The Institute of Statistical Mathematics 計算推論科学概論 II (2010 年度, 後期 ) グラフィカルモデルによる推論 確率伝搬法 Kenj Fukuzu he Insue of Sascal Maheacs 計算推論科学概論 II 年度 後期 Inference on Hdden Markov Model Inference on Hdden Markov Model Revew: HMM odel : hdden sae fne Inference Coue... for any Naïve

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

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

II The Z Transform. Topics to be covered. 1. Introduction. 2. The Z transform. 3. Z transforms of elementary functions

II The Z Transform. Topics to be covered. 1. Introduction. 2. The Z transform. 3. Z transforms of elementary functions II The Z Trnsfor Tocs o e covered. Inroducon. The Z rnsfor 3. Z rnsfors of eleenry funcons 4. Proeres nd Theory of rnsfor 5. The nverse rnsfor 6. Z rnsfor for solvng dfference equons II. Inroducon The

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

. 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

On Local Existence and Blow-Up of Solutions for Nonlinear Wave Equations of Higher-Order Kirchhoff Type with Strong Dissipation

On Local Existence and Blow-Up of Solutions for Nonlinear Wave Equations of Higher-Order Kirchhoff Type with Strong Dissipation Inernaonal Journal of Modern Nonlnear Theory and Alcaon 7 6-5 h://wwwscrorg/journal/jna ISSN Onlne: 67-987 ISSN Prn: 67-979 On Local Exsence and Blow-U of Soluons for Nonlnear Wave Euaons of Hgher-Order

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

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

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

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

Response of MDOF systems

Response of MDOF systems Response of MDOF syses Degree of freedo DOF: he nu nuber of ndependen coordnaes requred o deerne copleely he posons of all pars of a syse a any nsan of e. wo DOF syses hree DOF syses he noral ode analyss

More information

CONSISTENT ESTIMATION OF THE NUMBER OF DYNAMIC FACTORS IN A LARGE N AND T PANEL. Detailed Appendix

CONSISTENT ESTIMATION OF THE NUMBER OF DYNAMIC FACTORS IN A LARGE N AND T PANEL. Detailed Appendix COSISE ESIMAIO OF HE UMBER OF DYAMIC FACORS I A LARGE AD PAEL Dealed Aendx July 005 hs verson: May 9, 006 Dane Amengual Dearmen of Economcs, Prnceon Unversy and Mar W Wason* Woodrow Wlson School and Dearmen

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

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

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

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

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

THEORETICAL AUTOCORRELATIONS. ) if often denoted by γ. Note that

THEORETICAL AUTOCORRELATIONS. ) if often denoted by γ. Note that THEORETICAL AUTOCORRELATIONS Cov( y, y ) E( y E( y))( y E( y)) ρ = = Var( y) E( y E( y)) =,, L ρ = and Cov( y, y ) s ofen denoed by whle Var( y ) f ofen denoed by γ. Noe ha γ = γ and ρ = ρ and because

More information

MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES. Institute for Mathematical Research, Universiti Putra Malaysia, UPM Serdang, Selangor, Malaysia

MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES. Institute for Mathematical Research, Universiti Putra Malaysia, UPM Serdang, Selangor, Malaysia Malaysan Journal of Mahemacal Scences 9(2): 277-300 (2015) MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES Journal homeage: h://ensemumedumy/journal A Mehod for Deermnng -Adc Orders of Facorals 1* Rafka Zulkal,

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

( 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

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

An Improved Stabilization Method for Linear Time-Delay Systems

An Improved Stabilization Method for Linear Time-Delay Systems IEEE RASACIOS O AUOMAIC COROL, VOL. 47, O. 11, OVEMBER 191 An Improved Sablzaon Mehod for Lnear me-delay Sysems Emla Frdman and Ur Shaked Absrac In hs noe, we combne a new approach for lnear medelay sysems

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

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

THE PUBLISHING HOUSE PROCEEDINGS OF THE ROMANIAN ACADEMY, Series A, OF THE ROMANIAN ACADEMY Volume 9, Number 1/2008, pp

THE PUBLISHING HOUSE PROCEEDINGS OF THE ROMANIAN ACADEMY, Series A, OF THE ROMANIAN ACADEMY Volume 9, Number 1/2008, pp THE PUBLISHING HOUSE PROCEEDINGS OF THE ROMNIN CDEMY, Seres, OF THE ROMNIN CDEMY Volue 9, Nuber /008, pp. 000 000 ON CIMMINO'S REFLECTION LGORITHM Consann POP Ovdus Unversy of Consana, Roana, E-al: cpopa@unv-ovdus.ro

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

Endogeneity. Is the term given to the situation when one or more of the regressors in the model are correlated with the error term such that

Endogeneity. Is the term given to the situation when one or more of the regressors in the model are correlated with the error term such that s row Endogeney Is he erm gven o he suaon when one or more of he regressors n he model are correlaed wh he error erm such ha E( u 0 The 3 man causes of endogeney are: Measuremen error n he rgh hand sde

More information

Static Output-Feedback Simultaneous Stabilization of Interval Time-Delay Systems

Static Output-Feedback Simultaneous Stabilization of Interval Time-Delay Systems Sac Oupu-Feedback Sulaneous Sablzaon of Inerval e-delay Syses YUAN-CHANG CHANG SONG-SHYONG CHEN Deparen of Elecrcal Engneerng Lee-Mng Insue of echnology No. - Lee-Juan Road a-shan ape Couny 4305 AIWAN

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

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

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

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

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

A DECOMPOSITION METHOD FOR SOLVING DIFFUSION EQUATIONS VIA LOCAL FRACTIONAL TIME DERIVATIVE

A DECOMPOSITION METHOD FOR SOLVING DIFFUSION EQUATIONS VIA LOCAL FRACTIONAL TIME DERIVATIVE S13 A DECOMPOSITION METHOD FOR SOLVING DIFFUSION EQUATIONS VIA LOCAL FRACTIONAL TIME DERIVATIVE by Hossen JAFARI a,b, Haleh TAJADODI c, and Sarah Jane JOHNSTON a a Deparen of Maheacal Scences, Unversy

More information

NPTEL Project. Econometric Modelling. Module23: Granger Causality Test. Lecture35: Granger Causality Test. Vinod Gupta School of Management

NPTEL Project. Econometric Modelling. Module23: Granger Causality Test. Lecture35: Granger Causality Test. Vinod Gupta School of Management P age NPTEL Proec Economerc Modellng Vnod Gua School of Managemen Module23: Granger Causaly Tes Lecure35: Granger Causaly Tes Rudra P. Pradhan Vnod Gua School of Managemen Indan Insue of Technology Kharagur,

More information

Advanced time-series analysis (University of Lund, Economic History Department)

Advanced time-series analysis (University of Lund, Economic History Department) Advanced me-seres analss (Unvers of Lund, Economc Hsor Dearmen) 3 Jan-3 Februar and 6-3 March Lecure 4 Economerc echnues for saonar seres : Unvarae sochasc models wh Box- Jenns mehodolog, smle forecasng

More information

Periodic motions of a class of forced infinite lattices with nearest neighbor interaction

Periodic motions of a class of forced infinite lattices with nearest neighbor interaction J. Mah. Anal. Appl. 34 28 44 52 www.elsever.co/locae/jaa Peroc oons of a class of force nfne laces wh neares neghbor neracon Chao Wang a,b,, Dngban Qan a a School of Maheacal Scence, Suzhou Unversy, Suzhou

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

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

Water Hammer in Pipes

Water Hammer in Pipes Waer Haer Hydraulcs and Hydraulc Machnes Waer Haer n Pes H Pressure wave A B If waer s flowng along a long e and s suddenly brough o res by he closng of a valve, or by any slar cause, here wll be a sudden

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

Chapter 3: Signed-rank charts

Chapter 3: Signed-rank charts Chaer : gned-ran chars.. The hewhar-ye conrol char... Inroducon As menoned n Chaer, samles of fxed sze are aen a regular nervals and he long sasc s hen loed. The queson s: Whch qualy arameer should be

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

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

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

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

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

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

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

2 Aggregate demand in partial equilibrium static framework

2 Aggregate demand in partial equilibrium static framework Unversy of Mnnesoa 8107 Macroeconomc Theory, Sprng 2009, Mn 1 Fabrzo Perr Lecure 1. Aggregaon 1 Inroducon Probably so far n he macro sequence you have deal drecly wh represenave consumers and represenave

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

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

Track Properities of Normal Chain

Track Properities of Normal Chain In. J. Conemp. Mah. Scences, Vol. 8, 213, no. 4, 163-171 HIKARI Ld, www.m-har.com rac Propes of Normal Chan L Chen School of Mahemacs and Sascs, Zhengzhou Normal Unversy Zhengzhou Cy, Hennan Provnce, 4544,

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

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

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

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

Scattering at an Interface: Oblique Incidence

Scattering at an Interface: Oblique Incidence Course Insrucor Dr. Raymond C. Rumpf Offce: A 337 Phone: (915) 747 6958 E Mal: rcrumpf@uep.edu EE 4347 Appled Elecromagnecs Topc 3g Scaerng a an Inerface: Oblque Incdence Scaerng These Oblque noes may

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

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

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

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

Political Economy of Institutions and Development: Problem Set 2 Due Date: Thursday, March 15, 2019.

Political Economy of Institutions and Development: Problem Set 2 Due Date: Thursday, March 15, 2019. Polcal Economy of Insuons and Developmen: 14.773 Problem Se 2 Due Dae: Thursday, March 15, 2019. Please answer Quesons 1, 2 and 3. Queson 1 Consder an nfne-horzon dynamc game beween wo groups, an ele and

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

ON THE ADDITION OF UNITS AND NON-UNITS IN FINITE COMMUTATIVE RINGS

ON THE ADDITION OF UNITS AND NON-UNITS IN FINITE COMMUTATIVE RINGS ROCKY MOUNTAIN JOURNAL OF MATHEMATICS Volume 45, Number 6, 2015 ON THE ADDITION OF UNITS AND NON-UNITS IN FINITE COMMUTATIVE RINGS DARIUSH KIANI AND MOHSEN MOLLAHAJIAGHAEI ABSTRACT. Le R be a fne commuave

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

GMM parameter estimation. Xiaoye Lu CMPS290c Final Project

GMM parameter estimation. Xiaoye Lu CMPS290c Final Project GMM paraeer esaon Xaoye Lu M290c Fnal rojec GMM nroducon Gaussan ure Model obnaon of several gaussan coponens Noaon: For each Gaussan dsrbuon:, s he ean and covarance ar. A GMM h ures(coponens): p ( 2π

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

THERMODYNAMICS 1. The First Law and Other Basic Concepts (part 2)

THERMODYNAMICS 1. The First Law and Other Basic Concepts (part 2) Company LOGO THERMODYNAMICS The Frs Law and Oher Basc Conceps (par ) Deparmen of Chemcal Engneerng, Semarang Sae Unversy Dhon Harano S.T., M.T., M.Sc. Have you ever cooked? Equlbrum Equlbrum (con.) Equlbrum

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

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

THE POLYNOMIAL TENSOR INTERPOLATION

THE POLYNOMIAL TENSOR INTERPOLATION Pease ce hs arce as: Grzegorz Berna, Ana Ceo, The oynoma ensor neroaon, Scenfc Research of he Insue of Mahemacs and Comuer Scence, 28, oume 7, Issue, ages 5-. The webse: h://www.amcm.cz./ Scenfc Research

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

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

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

Pavel Azizurovich Rahman Ufa State Petroleum Technological University, Kosmonavtov St., 1, Ufa, Russian Federation

Pavel Azizurovich Rahman Ufa State Petroleum Technological University, Kosmonavtov St., 1, Ufa, Russian Federation VOL., NO. 5, MARCH 8 ISSN 89-668 ARN Journal of Engneerng and Aled Scences 6-8 Asan Research ublshng Nework ARN. All rghs reserved. www.arnjournals.com A CALCULATION METHOD FOR ESTIMATION OF THE MEAN TIME

More information

Sklar: Sections (4.4.2 is not covered).

Sklar: Sections (4.4.2 is not covered). COSC 44: Dgal Councaons Insrucor: Dr. Ar Asf Deparen of Copuer Scence and Engneerng York Unversy Handou # 6: Bandpass Modulaon opcs:. Phasor Represenaon. Dgal Modulaon Schees: PSK FSK ASK APK ASK/FSK)

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

Homework 8: Rigid Body Dynamics Due Friday April 21, 2017

Homework 8: Rigid Body Dynamics Due Friday April 21, 2017 EN40: Dynacs and Vbraons Hoework 8: gd Body Dynacs Due Frday Aprl 1, 017 School of Engneerng Brown Unversy 1. The earh s roaon rae has been esaed o decrease so as o ncrease he lengh of a day a a rae of

More information

Beyond Balanced Growth : Some Further Results

Beyond Balanced Growth : Some Further Results eyond alanced Growh : Some Furher Resuls by Dens Sec and Helmu Wagner Dscusson Paer o. 49 ay 27 Dskussonsberäge der Fakulä für Wrschafswssenschaf der FernUnversä n Hagen Herausgegeben vom Dekan der Fakulä

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

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

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 If we say wh one bass se, properes vary only because of changes n he coeffcens weghng each bass se funcon x = h< Ix > - hs s how we calculae

More information

Research Article Adaptive Synchronization of Complex Dynamical Networks with State Predictor

Research Article Adaptive Synchronization of Complex Dynamical Networks with State Predictor Appled Mahemacs Volume 3, Arcle ID 39437, 8 pages hp://dxdoorg/55/3/39437 Research Arcle Adapve Synchronzaon of Complex Dynamcal eworks wh Sae Predcor Yunao Sh, Bo Lu, and Xao Han Key Laboraory of Beng

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

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

UNIVERSITAT AUTÒNOMA DE BARCELONA MARCH 2017 EXAMINATION

UNIVERSITAT AUTÒNOMA DE BARCELONA MARCH 2017 EXAMINATION INTERNATIONAL TRADE T. J. KEHOE UNIVERSITAT AUTÒNOMA DE BARCELONA MARCH 27 EXAMINATION Please answer wo of he hree quesons. You can consul class noes, workng papers, and arcles whle you are workng on 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 Hqp (,,) = qp Lqq (,,) H p = q H q = p H L = Equvalen o Lagrangan formalsm Smpler, bu wce as

More information

AT&T Labs Research, Shannon Laboratory, 180 Park Avenue, Room A279, Florham Park, NJ , USA

AT&T Labs Research, Shannon Laboratory, 180 Park Avenue, Room A279, Florham Park, NJ , USA Machne Learnng, 43, 65 91, 001 c 001 Kluwer Acadec Publshers. Manufacured n The Neherlands. Drfng Gaes ROBERT E. SCHAPIRE schapre@research.a.co AT&T Labs Research, Shannon Laboraory, 180 Park Avenue, Roo

More information

January Examinations 2012

January Examinations 2012 Page of 5 EC79 January Examnaons No. of Pages: 5 No. of Quesons: 8 Subjec ECONOMICS (POSTGRADUATE) Tle of Paper EC79 QUANTITATIVE METHODS FOR BUSINESS AND FINANCE Tme Allowed Two Hours ( hours) Insrucons

More information

Graduate Macroeconomics 2 Problem set 5. - Solutions

Graduate Macroeconomics 2 Problem set 5. - Solutions Graduae Macroeconomcs 2 Problem se. - Soluons Queson 1 To answer hs queson we need he frms frs order condons and he equaon ha deermnes he number of frms n equlbrum. The frms frs order condons are: F K

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