Robust Control Design for Nonlinear Uncertain Systems with an Unknown Time-Varying Control Direction

Size: px
Start display at page:

Download "Robust Control Design for Nonlinear Uncertain Systems with an Unknown Time-Varying Control Direction"

Transcription

1 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 4, NO. 3, MARCH [7] L. Hsu, Variable srucure model reference adapive conrol using only I/O measuremen: General case, IEEE Trans. Auoma. Conr., vol. 35, no. 11, pp , [8] L. Hsu, A. D. Araújo, and R. R. Cosa, Analysis and design of I/O based variable srucure adapive conrol, IEEE Trans. Auoma. Conr., vol. 39, no. 1, pp. 4 1, [9] L. Hsu, A. D. Araújo, and F. Lizarralde, New resuls on VS-MRAC: Design and sabiliy analysis, in Proc. Amer. Conr. Conf., San Francisco, 1993, pp [10] L. Hsu and R. R. Cosa, Variable srucure model reference adapive conrol using only inpu and oupu measuremen: Par I, In. J. Conr., vol. 49, no., pp , [11] L. Hsu and F. Lizarralde, Redesign and sabiliy analysis of VS-MRAC sysems, in Proc. Amer. Conr. Conf., Chicago, 199, pp [1], Experimenal resuls on variable srucure adapive robo conrol wihou velociy measuremen, in Proc. Amer. Conr. Conf., Seale, 1995, pp [13] P. Ioannou and K. Tsakalis, A robus direc adapive conroller, IEEE Trans. Auoma. Conr., vol. 31, no. 11, pp , [14] S. M. Naik, P. R. Kumar, and B. E. Ydsie, Robus coninuous-ime adapive conrol by parameer projecion, IEEE Trans. Auoma. Conr., vol. 37, no., pp , 199. [15] K. Narendra and A. Annaswamy, Sable Adapive Sysems. Englewood Cliffs, NJ: Prenice-Hall, [16] S. S. Sasry and M. Bodson, Adapive Conrol: Sabiliy, Convergence and Robusness. Englewood Cliffs, NJ: Prenice-Hall, [17] V. I. Ukin, Sliding Modes and Their Applicaion in Variable Srucure Sysems. MIR, [18] A. C. Wu, L. C. Fu, and C. F Hsu, Robus MRAC for plan wih arbirary relaive degree using variable srucure design, in Proc. Amer. Conr. Conf., Chicago, 199, pp Robus Conrol Design for Nonlinear Uncerain Sysems wih an Unknown Time-Varying Conrol Direcion Joseph Kalous and Z. Qu Absrac In his paper, a coninuous robus conrol design approach is proposed for firs-order nonlinear sysems whose dynamics conain boh nonlinear uncerainy and an unknown ime-varying conrol direcion. The so-called conrol direcion is he muliplier of he conrol erm in he dynamic equaion, and i effecively represens he direcion of moion under any given conrol. A nonlinear robus conrol is designed o on-line and coninuously idenify sign changes of he unknown conrol direcion and o guaranee sabiliy of uniform ulimae boundedness. The proposed robus conrol design requires only hree condiions: he nominal sysem is sable, he conrol direcion is smooh, and he uncerainy in he sysem is bounded by a known funcion. The necessiy of hese condiions is esablished in his paper. Coninuiy of he proposed robus conrol is achieved by using a so-called shifing law ha changes smoohly he sign of robus conrol and racks he change of he unknown conrol direcion. Analysis and design is shown by using he Lyapunov s direc mehod. Index Terms Lyapunov sabiliy, robus conrol, uncerain sysems. and sabilizing robus conrollers have been developed, for example, he resuls in [1], [3], [10], [16] [19], and [1]. The uncerainies admissible in hese resuls can be of eiher mached ype or unmached ype, and hey are generally bounded by known nonlinear funcions of he sae. However, inpu-relaed uncerainies analyzed in he exising resuls are assumed o be of fixed and known signs; ha is, he direcion of moion of he sysem is fixed and known under any choice of conrol. The assumpion of priori knowledge of he conrol direcion being known makes robus conrol design sraighforward. In [5], he problem of designing robus conrol was sudied for he class of sysems wih an unknown bu consan and boundedaway-from-zero conrol direcion. The objecive of his paper is o generalize he resuls in [5] so ha sysems wih an unknown and ime-varying conrol direcion can be reaed. This exension is nonrivial since he ime-varying conrol direcion is allowed o cross zero an infinie number of imes and since idenificaion mus now be done fas and coninuously o keep racking he changing conrol direcion. As before, he Lyapunov s direc mehod is uilized o proceed wih he design and o ensure he sabiliy of uniform ulimae boundedness. Similar o he resuls inroduced in [5], coninuiy of he proposed robus conrol is achieved by designing he so-called shifing law, which coninuously and smoohly seers he sign of our robus conrol. In he area of adapive conrol, many resuls dealing wih sysems conaining an unknown consan conrol direcion, referred o herein as he unknown high-frequency gain, have been obained. For example, Nussbaum [15] proposed for linear ime-invarian sysems an adapive conrol algorihm whose sign changes an infinie number of imes as is argumen ends o infiniy. In [1], Mudge and his coworkers developed an adapive conrol design scheme for general linear sysems wihou requiring he knowledge of high-frequency gain by following he same principle as ha in [15]. In he effor of exending hese resuls o nonlinear sysems having unknown conrol direcions, wo adapive conrol schemes have been proposed for simple firs-order nonlinear sysems [], [6]. In [] and [6], a hyseresis mehod was used o consruc an adapive conrol ha jumps a finie number of imes over he infinie ime horizon. Along he line of nonlinear robus conrol, a coninuous robus conrol design scheme was developed in [5] for uncerain nonlinear sysems wih an unknown bu consan conrol direcion. Compared wih he exising resul, he proposed resul is he firs aemp in addressing he conrol design problem in he presence of an unknown bu smooh conrol direcion, whose sign change is of unlimied number and may occur anyime. The paper is broken ino he following secions. Sysem descripion, necessary assumpions, and he problem of designing robus conrol are described in Secion II. Robus conrol design is proceeded in Secion III by he Lyapunov s direc mehod. A simulaion example is presened in Secion IV. I. INTRODUCTION For more han one decade, robus conrol design of nonlinear uncerain sysems has been an acive area of research. Several imporan classes of sabilizable uncerain sysems have been idenified, Manuscrip received March 14, 1995; revised Sepember 8, 1995, April 0, 1996, and July 1, J. Kalous is wih Loral Vough Sysems, Dallas, TX USA. Z. Qu is wih he Deparmen of Elecrical and Compuer Engineering, Universiy of Cenral Florida, Orlando, FL 3816 USA. Publisher Iem Idenifier S (97) II. PROBLEM FORMULATION We shall consider he nonlinear uncerain sysem described by _x = f (x; ) +1f 0 (x; ) +a()u(x; ) (1) where x < is he sae, u < is he conrol inpu, f (x; ) and 1f 0 (x; ) are known and unknown dynamics of he sysem, respecively. Symbol a() characerizes he ime-varying direcion of moion wih respec o a given conrol. The sysem is chosen o be scalar in order o highligh he developmen of robus conrol design /97$ IEEE

2 394 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 4, NO. 3, MARCH 1997 in he presence of an unknown, ime-varying conrol direcion. I is obvious ha (1) can be made of a vecor and ha, based on he resuls in [5], i can also be exended o be of higher order. The conrol direcion, sign[a()] of funcion a(); is unknown, and i may swich over ime among 01, 0, and 1 arbirarily. For idenifiabiliy, an upper bound will be imposed on he rae of change of funcion a(); and he necessiy of such an assumpion can easily be shown by conradicion. Since a() is allowed o change is sign, he case ha a() may be zero all he ime mus be considered. Noe ha a() 0 implies ha here is no conrol acion, no maer wha u is. In his case, moion and sabiliy of he sysem are governed solely by known and unknown dynamics, bu hese dynamics, if unsable, mus be compensaed. For his reason, wo necessary assumpions mus be made. Firs, known dynamics f (x; ) have o be sable (for he case ha a() 0 and 1f (x; ) 0). Second, he uncerainy in he sysem has o have he propery ha for some well defined funcion 1f (x; ); 1f 0 (x; ) = a()1f(x; ): Wihou a() being he common facor beween he erms of conrol and uncerainy, sabilizabiliy canno be guaraneed if he bounding funcion of he uncerainy is arbirary. Tha is, if funcion a() is muliplied only a he fron of conrol u; he sysem is sable only if he uncerainy is small in he sense ha i is dominaed by he robus sabiliy margin. This laer case is no considered in he paper, since sabiliy is a hand simply by leing u =0; herefore, he conrol design is rivial. Wih hese facs in mind, we can define he following assumpions for sysems in he form of (1). Assumpion 1: Under a coninuous conrol u o be devised, he dynamics of (1) are Caraheodory. Tha is, (1) has a classical soluion under any coninuous conrol. Assumpion : Funcion f (x; ) characerizing he dynamics of he nominal sysem is exponenially sable. Tha is, by he Lyapunov converse heorem, here exiss a differeniable funcion V 1(x; ) such ha, for some posiive consans 1; ; 3; and 4 ; 1 jxj V 1 (x; ) jxj ; j =@j 3 jxj; and =@ +( =)f (x; ) 0 4 jxj : Assumpion 3: Lumped uncerainy 1f 0 (x; ) can be rewrien as 1f 0 (x; ) = a()1f(x; ), and 1f (x; ) is bounded as, for all (x; ); j1f (x; )j (x); where (1) is known, coninuous, and locally uniformly bounded wih respec o x: Assumpion 4: Funcion a() is bounded and normalized in he sense ha ja()j 1:In addiion, a() is smooh in he sense ha is firs-order derivaive is bounded from above as j _a()j for a known consan >0: Assumpion 5: For all possible choices of he uncerainy 1f 0 (x; ) bounded by ( + )(x) for some consan >0; (1) wih u 0 does no have a finie escape ime ha is infiniesimal. Mos sysems in he form of (1) saisfy Assumpion 5. To see his conclusion, consider he choice of bounding funcion (x) ha (x) =jxj p for some p>0:then, i follows from Assumpion and from u =0ha, along every rajecory of (1), he ime derivaive of Lyapunov funcion V 1 (x) saisfies _V jxj +(+)jxj(x) 0 4 V p+1 V p+1= 1 : I is obvious from he upper bound on he growh rae ( _ V 1 ) ha Lyapunov funcion V 1, and herefore he sae x, will no become infinie insananeously. This resul holds for any choice of posiive real number p: Noe ha by he Taylor series expansion, polynomial funcions form a local basis for all differeniable funcions. Neverheless, o avoid he complicaion of proving Assumpion 5 for all possible expressions of bounding funcion (1); we choose o make he assumpion. Once (x) is given, one can simply proceed wih he proof of Assumpion 5 using he above argumen. Assumpion 5 is needed for on-line idenificaion of he unknown and ime-varying direcion. If Assumpion 5 does no hold, no idenificaion algorihm can idenify insananeously and keep perfecly racking an unknown changing conrol direcion so ha a robus conrol of he correc sign can be implemened. The reason ha he bounding funcion used in he assumpion is ( + )(x) raher han (x) is ha robus conrol mus have is magniude no smaller han (x) for compensaing for uncerainy of size (x) and ha he wors case of robus conrol iniially having a wrong sign and consequenly being desabilizing mus be considered. In he case ha he finie escape ime is infiniesimal, here is no conrol of finie magniude ha can pull he sysem back, even if he conrol direcion can be idenified insananeously. Based on hese assumpions, we can proceed wih robus conrol design in he nex secion. III. ROBUST CONTROL DESIGN The main approach for designing robus conrol for nonlinear uncerain sysems has been he Lyapunov s direc mehod. In is applicaion, deerminisic funcions are used o bound uncerainies, and several classes of uncerain sysems have been shown o be globally sabilizable [1], [3], [10], [16], [19], [1]. More recen developmen along his line of work is he recursive-inerlacing design procedure [17], [18] which is generic enough o overcome he common srucural condiions imposed on uncerain sysems. There are wo common feaures among hese resuls. Firs, global sabilizing robus conrollers are devised if bounding funcions of uncerainies are available. The selecion of a robus conroller is done hrough a Lyapunov argumen such ha robus conrol dominaes uncerainies in boh magniude and sign. Second, in hese resuls he conrol direcions are known in he sense ha afer choosing a robus conroller, he sign of he robus conrol deermines he signs of he ime derivaives of boh he sae and he Lyapunov funcion. I was shown in [5] ha he dominaing effec of robus conrol can be used o idenify a consan unknown conrol direcion and ha a coninuous robus conrol can be devised whose sign is changed smoohly. In his paper, robus conrol design is sudied for sysems wih an unknown and ime-varying conrol direcion, which can be viewed as an exension of [5]. Unlike he problem of dealing wih a fixed unknown conrol direcion, he ime-varying conrol direcion considered here is no bounded away from zero, he number of is sign changes is unlimied over ime, and a successful idenificaion algorihm mus rack closely is sign changes and be capable of handling he poenial problem of singulariy a zero crossings. Thus, he proposed exension is no rivial. Like he exising resuls, he design of a sabilizing robus conrol in he presence of a ime-varying and unknown conrol direcion involves a Lyapunov argumen and he concep of having he conrol dominae he unknown dynamics. As in [5], an iniial guess is made on he sign of he unknown conrol direcion, idenificaion of he unknown conrol direcion is hen carried ou on-line coninuously, and a robus conrol is designed whose sign is modified according o he idenified conrol direcion. The key resuls in his sudy are ha despie he fac ha he unknown conrol direcion may change is sign an infinie number of imes, idenificaion and racking can be done using he concep of size dominaion, and since a coninuous conrol usually produces a beer ransien response han disconinuous conrols, a smooh ransiion of he sign of he robus conrol can be achieved. Inuiively, he basic idea emanaes as follows: begin wih some guess ^a() of a(); design a robus conrol

3 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 4, NO. 3, MARCH in erms of ^a(); idenify a() on line; formulae shifing law _^a() so ha he robus conrol is coninuous; and le all hree designs be coordinaed and inegraed in such a way ha robus sabiliy is ensured. Robus conrol design in he presence of a ime-varying conrol direcion consiss of hree inegraed pars: 1) selecion of robus conrol using he Lyapunov s direc mehod in a similar fashion as ha wih a known conrol direcion; ) a scheme ha on-line idenifies he unknown ime-varying conrol direcion; 3) formulae a coninuous shifing law ha changes he sign of robus conrol smoohly in racking he idenified sign of he imevarying conrol direcion. These hree pars are deailed as follows. Par I: In his par, a candidae of robus conrol is seleced. As shown in [1] and [19], robus conrol in he absence of an unknown conrol direcion can be seleced o be 1 j(x)j + u r = 0 j(x)j (x)(x) () + 3 where >0is a design parameer, and (x) =( =)(x): Noe ha ju r j = j(x)j + j(x)j j(x)j + 3 (x) (x): (3) The above inequaliy holds since ( )(j(x)j + 3 ) j(x)j j(x)j j(x)j + j(x)j (4) in which he Cauchy Scharwz inequaliy a + b ab is used. In he presence of he ime-varying, unknown conrol direcion, we propose he following robus conrol: u(x; ^a()) = ^a() () 1 u r (5) where u r is defined in (), ^a() represens on-line esimae of funcion a(); and () >>0 is a design ime funcion such ha for all 0 j^a()j (): (6) Funcion (); o be defined shorly, is o normalize robus conrol (5) so ha is magniude is compaible wih conrol () and herefore wih he class of uncerainies under consideraion. I follows from (3) and (6) ha ja()[1f (x; ) +u]j(x)+ju rj( )(x): (7) Therefore, by choosing small such ha 4; we can always ensure ha he condiion in Assumpion 5 holds, which is why (6) was imposed on conrol (5). Noe ha besides he requiremen of 4; in (5) or () can be chosen freely by he designer, and he smaller is, he smaller he ulimae bound on he sae. As inroduced in [5], he concep of size dominaion means ha in a Lyapunov argumen, he erm associaed wih he conrol dominaes in magniude, and hus in sign, he erm(s) associaed wih uncerainies. If size dominaion is possible, he sign of he ime derivaive of he Lyapunov funcion is dicaed by ha of he erm associaed wih he robus conrol, and herefore sabiliy of he conrolled sysem can be ensured. Since ^a() may become zero and since () in robus conrol (5) mus be bounded from zero, we choose funcion () o be j^a()j if j^a()j () = (8) if j^a()j <: The exac role of funcion ; a design consan, will become clear shorly. Obviously, choice (8) of () makes (6) hold. In he case ha he conrol direcion is ime invarian, seing () === j^a( 0 )jreduces conrol (5) o ha in [5] (in [5], ^a( 0 ) is se o be eiher 1 or 01 so ha () =1): As will be shown in he upcoming sabiliy analysis, size dominaion of robus conrol over uncerainy is achieved in a Lyapunov argumen if u > 1f(x; ) (9) or equivalenly j^a()j () 1 j(x)j + j(x)j j(x)j j(x)j + 3 j(x)j: Therefore, we know ha he size dominaion is accomplished under he following condiions: j(x)j (10) and j^a()j : (11) Noe ha esimae ^a() has o be able o cross zero in order o rack he idenified version of a() and ha whenever j^a()j <; funcion ^a() may soon have zero crossing. The size dominaion in erms of (9) is no longer possible, and he sign of a() may no be idenified correcly. Since insabiliy will occur if he robus conrol provides a wrong conrol sign for a period of ime, ^a() mus rack a() closely. However, i is impossible for ^a() o correcly rack a() a all ime. To eliminae he possibiliy of insabiliy, we shall devise a sign idenificaion scheme which successfully idenifies he sign of a() unless eiher a() is abou o change is sign or x is very small. I is by limiing he occurence of idenificaion failure o be only wihin he described regions in he spaces of x and a() ha he sabiliy of he overall sysem can be guaraneed. Par II: In his par, a scheme o on-line idenify he conrol direcion is synhesized using he concep of size dominaion. Since he conrol direcion is unknown, any conrol will fail o ensure sabiliy wihou a correc idenificaion of he conrol direcion. If a() is consan, idenificaion of he conrol direcion needs o be performed only for a very small inerval of ime as was proposed in [5], and robus conrol can hen be operaed forever based on he one-ime idenificaion. Since he conrol direcion considered in his paper is a smooh ime funcion, is sign changes, possibly an infinie number of imes, mus be accouned for coninuously. In order o idenify sign[a()]; muliplying = on boh sides of (1) and hen inegraing he boh sides yield, for all [ 1; ] [ 0 ;1) dx 0 f(x();) d = a( ) 1f(x();)+ u(x()) d: (1) Noe ha he inegrand funcions are all coninuous and ha he sign of erm (=)u is he negaive of he sign of ^a(): Recall ha under condiions defined in (10) and (11), u(x) will dominae uncerainy 1f (x; ) in magniude as specified by (9). Therefore, i follows from he mean value heorem ha under condiions (10) and (11), he sign of a() can be idenified by he following equaion if boh a() and ^a() do no change heir signs in he inerval [ 1 ; ]: sign[a()] = 0sign[^a()]sign dx 0 f(x();) d : (13) Wih size dominaion of robus conrol over uncerainy, (13) can be used o idenify he correc conrol direcion, so long as a() and

4 396 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 4, NO. 3, MARCH 1997 ^a() do no change heir signs in he inerval [ 1 ; ]: If here is a sign change in eiher a() or ^a(); he resul from (13) mus be considered o be wrong. Since ^a() is known, is sign change will be known as well so ha proper precauion can be buil ino he idenificaion and robus conrol schemes o overcome any possible wrong idenificaion. In fac, as will be shown in he upcoming sabiliy analysis, he soluion is simply o make he convergence of ^a() o an esimae of sign[a()] fas enough, regardless of wheher he esimae of sign[a()] is righ or wrong. On he oher hand, since a() is unknown, is sign change canno be prediced. Therefore, (13) may poenially produce a wrong sign a any ime, which is a major sabiliy problem. Our soluion o his sabiliy problem is o selec judiciously he inegraion inerval [ 1 ; ] and o use he fac ha he nominal sysem is exponenially sable. Specifically, we choose he inegraion inerval as follows: firs le = and 1 = 0 T so ha he inegraion inerval is a sliding window of he curren ime insan. Second, selec he inerval such ha is lengh T saisfies he following inequaliy: T< (14) where is defined in Assumpion 4, and 0 < 1 is he design parameer inroduced in (8). The role of consan allows us o analyze he effec of zero crossing of a() and o ake proper precauions in our robus conrol scheme. Specifically, if ja()j (namely ja( )j) is greaer han or equal o ; here will be no sign change in he inerval s [ 1; ]because of he choice of T in (14). Tha is, as far as a() is concerned, he conrol direcion will be idenified correcly if ja()j >: (15) I is when ja( )j < ha here may be a sign change occurring in he inerval. Thus, assuming ha a correc idenificaion of conrol direcion implies sabiliy (which will be guaraneed by designs of he robus conrol and he shifing law), robus sabiliy can be esablished by sudying wheher he conrolled sysem is sable if ja()j <: Recalling ha he nominal sysem obained by seing a() = 0is exponenial, we know ha robus sabiliy can be esablished for he sysem (wih he robus conrol of a wrong sign) by simply leing be sufficienly small. Deails on selecing will be given in he sabiliy proof. In fac, such a reamen of zoning is inuiive since, if ja()j < 1; sign idenificaion is no possible anyway as a() is ime varying. Summarizing he above discussions, we know ha he conrol direcion found from (13) is correc if (10), (11), and (15) hold. Par III: As soon as an esimae of sign[a()] becomes available, a shifing law can be consruced so ha he sign of robus conrol is changed smoohly from a possibly wrongly guessed conrol direcion o he idenified conrol direcion and is hen updaed all he ime. The shifing law given in (16), shown a he boom of he page, is chosen such ha he overall sysem is robusly sable during he ime period in which smooh ransiions of he sign of robus conrol are accomplished, where he iniial guess of he conrol direcion a(); denoed by ^a( 0 ); is chosen o be eiher 1 or 0 or 01. k s is a consan gain saisfying k s p 4 (17) and funcion w() is he idenificaion resul of sign[a()] defined by w() =0sign[^a()] sign 0T dx 0 f(x();) d : 0T Implemenaion of shifing law (16) requires w(): No maer wheher he idenificaion resul w() is he correc value of sign[a()] or no, w() has only hree discree values: 01, 1, and 0. As soon as w() picks up is new value, w() will remain o be a consan unil one of inequaliies (10), (11), and (15) fails. During he period ha w() is consan, we have according o (16) ha d[^a() 0 w()] =[^a()0w()] d^a d d 0 k s p j^a()0w()j: Solving he above differenial inequaliy, we know ha he esimae ^a() will converge o w() in a coninuous fashion, in a finie ime, and wihou any overshoo, ha is j^a()j 1: (18) If w() = sign[a()]; he esimae ^a() will converge quickly o he correc conrol direcion. In fac, he gain k s defined in (17) for shifing law (16) is chosen such ha he ime for ^a() o converge o sign[a()] is less han T defined in (14). This choice ensures ha ^a() is capable of racking sign[a()] excep for he case ha a() is in he zone of ja()j < and is on he verge of changing is sign. In he case ha he idenificaion resul is wrong, he esimae may converge o a wrong direcion. Sabiliy analysis mus be done o show ha any wrong idenificaion is emporary and ha he sysem sae will no escape from he sabiliy region anyime. Since (10) has been incorporaed ino he shifing law, we have w() = sign[a()] by implemening he shifing law if (11) and (15) hold. Tha is, sabiliy mus be concluded for he cases ha hese wo inequaliies fail. Global sabiliy of (1) under robus conrol (5) and shifing law (16) is analyzed by he Lyapunov s direc mehod using he Lyapunov funcion candidae V (x; ^a()) = V 1(x; ) +V (^a();) (19) where V 1(1) is ha given in Assumpion, and 1 V (^a();)= [a()^a() 0ja()j] ja()j = 1 [ja()j3=^a ()+ja()j 3= 0 sign[a()]ja()j 3=^a()]: Noe ha for any differenial funcion q(); jq()j p and sign[q ()]jq()j p wih p>1 are also differeniable. The following heorem shows he sabiliy of uniform and ulimae boundedness for (1). Theorem: Assume ha (1) saisfies Assumpions 1 5. Then, under robus conrol (5) and shifing law (16), he sysem sae can be made uniformly and ulimaely bounded by leing be small enough. if j(x)j () 0 if j(x)j < _^a() = 0sign[^a() 0 w()] 1 1 p k s j(x)j (16)

5 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 4, NO. 3, MARCH Proof: Differeniaing (19) along he rajecory of (1) yields f(x; ) + a() 1f(x; ) +a() u [1 + ^a ()] ja()j sign[a()] 0 3 ^a() ja()j _a() + a() [a()^a() 0ja()j]_^a ja()j 0 4 jxj +ja()j1j(x)j 0a()^a() j(x)j + j(x)j () j(x)j + 3 j(x)j+ 3 4 [1 + ^a ()] ja()j sign[a()] 0 3 ^a() ja()j _a() + a() ja()j [a()^a() 0ja()j]_^a: (0) Noe ha ^a() is bounded by one as shown by (18). Thus, sabiliy may also be sudied by Lyapunov funcion candidae V 1 (x; ): I follows ha _V 1 + f(x; ) +a() 1f(x; ) + u 0 4 jxj +ja()j1j(x)j 0a()^a() j(x)j + j(x)j () j(x)j + 3 j(x)j: (1) We shall sudy he sabiliy resul by invesigaing he following four cases. Case I j(x)j <: In his case, (10) does no hold, and herefore robus conrol does no dominae he uncerainy in size. In his case, w() may no be sign[a()]: I follows from Assumpion 4, shifing law (16), (6), and (18) ha he ime derivaive (0) of Lyapunov funcion V is bounded from above as _V 0 4 jxj +ja()j1j(x)j+ja()j 1 j^a()j () j(x)j + j(x)j j(x)j + 3 j(x)j [1 + ^a ()] ja()j + 3 j^a()j ja()j j _a()j + a() [a()^a() 0ja()j]_^a ja()j 0 4 jxj ++3: Noing ha ^a() is frozen because of _^a() = 0;here is no need o have a negaive definie erm associaed wih he variables in V (^a();):in fac, i follows from he same inequaliies ha he ime derivaive (1) of Lyapunov funcion V 1 is bounded from above as _V jxj +ja()j1j(x)j+ja()j 1 j^a()j () 0 4 jxj +: j(x)j + j(x)j j(x)j + 3 j(x)j From eiher of he above resuls, one can conclude using he sabiliy heorem in [1] ha he sysem is uniformly and ulimaely bounded. Case II ja()j : Noe ha alhough he esimae ^a() of he conrol direcion may converge o a wrong value, ^a() is bounded by one as indicaed by (18). I follows from Assumpion, (6), (4), (8), and (18) ha he ime derivaive (1) of Lyapunov funcion is bounded from above as _V 1 0 4jxj +1j(x)j j(x)j 0 4 jxj jxj(x): I follows from he resuls in [0] ha he sysem is semiglobal sable in he sense ha he sysem is uniformly ulimaely bounded and ha he sabiliy region can be made arbirarily large by leing be small enough. This is rue wheher he on-line idenificaion is done correcly or no. Case III j^a()j : In his case, he esimae ^a() is confined o be in a smaller se han is ulimae bound defined by (18). Since ^a() is small, he conribuion of he robus conrol may no be large enough o compensae for he uncerainy. Thus, his is he case where insabiliy can poenially occur. However, by he choice of gain k s in (17), his case will las a very small period of ime during which Assumpion 5 is saisfied [as ensured by (7)]. Therefore, by choosing small enough, no only can he ime inerval during which his case occurs be small enough, bu also he sae of he sysem will say in a finie region, specifically, he semiglobal sabiliy region defined in Case II. Case IV: j(x)j ; j^a()j >; and ja()j >: This case is he complemen of he union of he previous hree cases. Since condiions (10), (15), and (11) hold, idenificaion is done correcly. Tha is _^a() =0sign[^a() 0 sign[a()]] 1 1 p k s j(x)j : () () I follows from Assumpion 4, shifing law (), (4), and (18) ha he ime derivaive (0) of Lyapunov funcion V is bounded from above as _V 0 4 jxj +ja()j1j(x)j 0a()^a() j(x)j + j(x)j j(x)j j(x)j + 3 () +3+ Noe ha <() 1 and ha since j(x)j >: Hence a() [a()^a() 0ja()j]_^a: ja()j j(x)j + j(x)j j(x)j _V 0 4 jxj +[ja()j0a()^a()] 1 j(x)j + j(x)j j(x)j j(x)j + 3 () +3 + a() [a()^a() 0ja()j]_^a ja()j 0 4jxj 0k sja()^a() 0ja()jj +3 from which sabiliy of uniform and ulimae boundedness can be concluded. In summary, one can conclude he sabiliy of uniform ulimae boundedness for (1) by simply combining he above four cases. In he nex secion, a simulaion example is presened o illusrae our proposed robus conrol design.

6 398 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 4, NO. 3, MARCH 1997 Fig. 1. Sysem performance under a() = sin() and x(0) = : Fig.. Shifing law ^a (solid sign[a()]; dashed ^a(); doed u): Fig. 3. a() =0sin() and x(0) = (doed x; dashed ^a; solid sign[a()]): IV. SIMULATION To illusrae he proposed robus conrol scheme under an unknown ime-varying conrol direcion, simulaions were performed for (1) wih he following choices: known dynamics f(x; ) =0kx wih k =0:5and V 1 = 1 x ; uncerainy 1f =10sin()+[5+5cos()]x +10x 3 which is chosen o be unsable; bounding funcion (x) is chosen o be (x) = 10(1+x +jxj 3 ); unknown conrol direcion is given by ime-varying funcion a() = 6sin() (boh cases were simulaed). The rae of change is bounded by =; design parameers in shifing law (16) and in robus conrol (5) are =0:01;k s = 10; and =0:1; iniial guess of he conrol direcion is se o be ^a(0) = 0 (for all he cases simulaed); iniial condiions of he sae were chosen o be x(0) = 61 (boh cases were simulaed). The simulaions were carried ou using SIMNON. Sysem performance under all possible combinaions of iniial condiions of sae x and unknown conrol direcion a() are shown in Figs. 1, 3, 4, and 5. In Fig. 1, rajecories of he sae, he esimae of he conrol direcion, and he robus conrol are shown separaely. In Figs. 3 5, rajecories of x; ^a; and sign[a()] are combined ino one plo for each case. Conrol u is omied from hese figures for briefness since conrol rajecories are similar o ha in Fig. 1 and since he magniude of he conrol, if u is incorporaed ino he figures, will make he oher rajecories obscure. I should be noed from Fig. 4 ha he speed of ^a racking sign[a()] becomes slower as he sae x converges o he origin. This phenomenon is expeced from (10). To beer illusrae he behavior of he shifing law and is corresponding sabiliy propery, a porion of Fig. 1 (around =)is enlarged o be Fig. in which sign[a()] firs swiches is sign, ^a() soon follows, and in beween, he conrol emporarily has a wrong sign and drives he sysem in a wrong direcion (see he firs plo

7 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 4, NO. 3, MARCH Fig. 4. a() =0sin() and x(0) = 0 (doed x; dashed ^a; solid sign[a()]): Fig. 5. a() = sin() and x(0) = 0 (doed x; dashed ^a; solid sign[a()]): in Fig. 1). Obviously, he simulaion shows ha his increase in jxj does no yield any insabiliy. In fac, he emporary increase of he magniude of he sae makes i possible o idenify he curren sign of a(), idenificaion in urn will provide he correc conrol direcion, and he sae is hen forced o converge back o he origin. Similar resuls can be seen around all poins a which he sign of a() is swiched. V. CONCLUSION Robus conrol of a class of nonlinear uncerain sysems is sudied. A sysem under consideraion conains no only nonlinear uncerainy bounded by a known nonlinear funcion of he sae bu also an unknown ime-varying conrol direcion. Wihou any priori knowledge of he conrol direcion excep for smoohness, a hreepar design is proposed for consrucing robus conrol via he Lyapunov s direc mehod. A resuling robus conrol is coninuous and consiss of an on-line algorihm o idenify coninuously he conrol direcion, a shifing law ha changes he sign of robus conrol smoohly, and a robus conrol modified from he sandard robus conroller. Essenially, he design only requires he bounding funcion of unknown dynamics and he maximum rae of change of he conrol direcion. Sabiliy of uniform ulimae boundedness is guaraneed. REFERENCES [1] M. J. Corless and G. Leimann, Coninuous sae feedback guaraneeing uniform ulimae boundedness for uncerain dynamical sysems, IEEE Trans. Auoma. Conr., vol. 6, pp , [] F. Giri, P. A. Ioannou, and F. Ahmed-Zaid, A sable adapive conrol scheme for firs order plans wih no priori knowledge on he parameers, IEEE Trans. Auoma. Conr., vol. 38, no. 5, pp , [3] S. Guman, Uncerain dynamical sysems A Lyapunov min max approach, IEEE Trans. Auoma. Conr., vol. 4, no. 3, [4] U. Helmke and D. Prazel-Wolers, Sabiliy and robusness properies of universal adapive conrollers for firs order linear sysems, Inern. J. Conr., vol. 48, pp , [5] J. Kalous and Z. Qu, Coninuous robus conrol design for nonlinear uncerain sysems wihou a priori knowledge of he conrol direcion, IEEE Trans. Auoma. Conr., vol. 40, no., pp. 76 8, [6] R. Lozano and B. Brogliao, Adapive conrol of a simple nonlinear sysem wihou a priori informaion on he plan parameer, IEEE Trans. Auoma. Conr., vol. 37, no. 1, pp , 199. [7] R. Lozano, J. Collado, and S. Mondie, Model reference adapive conrol wihou a priori knowledge of he high frequency gain, IEEE Trans. Auoma. Conr., vol. 35, no. 1, pp , [8] R. T. M Closkey and R. M. Murray, Exponenial convergence of nonholonomic sysems: Some analysis ools, in Proc Conf. Decision Conr., 1993, pp [9], Experimens in exponenial sabilizaion of a mobile robo owing a railer, in Proc Amer. Conr. Conf., 1994, pp [10] R. Marino and P. Tomei, Robus sabilizaion of feedback linearizable ime-varying uncerain nonlinear sysems, Auomaica, vol. 9, pp , [11] B. Marensson, Remarks on adapive sabilizaion of firs order nonlinear sysems, Sys. Conr. Le., vol. 14, pp. 1 7, [1] D. R. Mudge and A. S. Morse, Adapive sabilizaion of linear sysems wih unknown high frequency gain, IEEE Trans. Auoma. Conr., vol. 30, pp , [13] R. M. Murray, Z. Li, and S. S. Sasry, A Mahemaical Inroducion o Roboic Manipulaion. New York: CRC, [14] R. M. Murray and S. S. Sasry, Nonholonomic moion planning: Seering using sinusoids, IEEE Trans. Auoma. Conr., vol. 38, no. 5, pp , [15] R. D. Nussbaum, Some remarks on he conjecure in parameer adapive conrol, Sys. Conr. Le., vol. 3, no. 5, pp , [16] Z. Qu, Global sabilizaion of nonlinear sysems wih a class of unmached uncerainies, Sys. Conr. Le., vol. 18, no. 3, pp , 199. [17], Robus conrol of nonlinear uncerain sysems wihou generalized maching condiions, IEEE Trans. Auoma. Conr., vol. 40, no. 8, pp , [18], A new generic procedure of designing robus conrol for nonlinear uncerain sysems: Beyond he back-sepping design, in Proc Wkshp. Robus Conrol via Variable Srucure Lyapunov Techniques, Beneveno, Ialy, Sep. 1994, pp [19], Robus conrol of nonlinear uncerain sysems under generalized maching condiions, Auomaica, vol. 9, pp , [0] Z. Qu and J. Dorsey, Robus conrol of generalized dynamic sysems wihou maching condiions, Trans. ASME, J. Dyn. Sys., Meas., Conr., vol. 113, pp , [1] J. J. E. Sloine and K. Hedrick, Robus inpu oupu feedback linearizaion, Inern. J. Conr., vol. 57, pp , 1993.

Sliding Mode Controller for Unstable Systems

Sliding Mode Controller for Unstable Systems S. SIVARAMAKRISHNAN e al., Sliding Mode Conroller for Unsable Sysems, Chem. Biochem. Eng. Q. 22 (1) 41 47 (28) 41 Sliding Mode Conroller for Unsable Sysems S. Sivaramakrishnan, A. K. Tangirala, and M.

More information

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,

More information

di Bernardo, M. (1995). A purely adaptive controller to synchronize and control chaotic systems.

di Bernardo, M. (1995). A purely adaptive controller to synchronize and control chaotic systems. di ernardo, M. (995). A purely adapive conroller o synchronize and conrol chaoic sysems. hps://doi.org/.6/375-96(96)8-x Early version, also known as pre-prin Link o published version (if available):.6/375-96(96)8-x

More information

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals

More information

Vehicle Arrival Models : Headway

Vehicle Arrival Models : Headway Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where

More information

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3 and d = c b - b c c d = c b - b c c This process is coninued unil he nh row has been compleed. The complee array of coefficiens is riangular. Noe ha in developing he array an enire row may be divided or

More information

2. Nonlinear Conservation Law Equations

2. Nonlinear Conservation Law Equations . Nonlinear Conservaion Law Equaions One of he clear lessons learned over recen years in sudying nonlinear parial differenial equaions is ha i is generally no wise o ry o aack a general class of nonlinear

More information

(1) (2) Differentiation of (1) and then substitution of (3) leads to. Therefore, we will simply consider the second-order linear system given by (4)

(1) (2) Differentiation of (1) and then substitution of (3) leads to. Therefore, we will simply consider the second-order linear system given by (4) Phase Plane Analysis of Linear Sysems Adaped from Applied Nonlinear Conrol by Sloine and Li The general form of a linear second-order sysem is a c b d From and b bc d a Differeniaion of and hen subsiuion

More information

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid

More information

15. Vector Valued Functions

15. Vector Valued Functions 1. Vecor Valued Funcions Up o his poin, we have presened vecors wih consan componens, for example, 1, and,,4. However, we can allow he componens of a vecor o be funcions of a common variable. For example,

More information

A DELAY-DEPENDENT STABILITY CRITERIA FOR T-S FUZZY SYSTEM WITH TIME-DELAYS

A DELAY-DEPENDENT STABILITY CRITERIA FOR T-S FUZZY SYSTEM WITH TIME-DELAYS A DELAY-DEPENDENT STABILITY CRITERIA FOR T-S FUZZY SYSTEM WITH TIME-DELAYS Xinping Guan ;1 Fenglei Li Cailian Chen Insiue of Elecrical Engineering, Yanshan Universiy, Qinhuangdao, 066004, China. Deparmen

More information

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon 3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of

More information

Some Basic Information about M-S-D Systems

Some Basic Information about M-S-D Systems Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,

More information

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients Secion 3.5 Nonhomogeneous Equaions; Mehod of Undeermined Coefficiens Key Terms/Ideas: Linear Differenial operaor Nonlinear operaor Second order homogeneous DE Second order nonhomogeneous DE Soluion o homogeneous

More information

In this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should

In this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should Cambridge Universiy Press 978--36-60033-7 Cambridge Inernaional AS and A Level Mahemaics: Mechanics Coursebook Excerp More Informaion Chaper The moion of projeciles In his chaper he model of free moion

More information

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions Muli-Period Sochasic Models: Opimali of (s, S) Polic for -Convex Objecive Funcions Consider a seing similar o he N-sage newsvendor problem excep ha now here is a fixed re-ordering cos (> 0) for each (re-)order.

More information

STATE-SPACE MODELLING. A mass balance across the tank gives:

STATE-SPACE MODELLING. A mass balance across the tank gives: B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing

More information

Hamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation:

Hamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation: M ah 5 7 Fall 9 L ecure O c. 4, 9 ) Hamilon- J acobi Equaion: Weak S oluion We coninue he sudy of he Hamilon-Jacobi equaion: We have shown ha u + H D u) = R n, ) ; u = g R n { = }. ). In general we canno

More information

Stability and Bifurcation in a Neural Network Model with Two Delays

Stability and Bifurcation in a Neural Network Model with Two Delays Inernaional Mahemaical Forum, Vol. 6, 11, no. 35, 175-1731 Sabiliy and Bifurcaion in a Neural Nework Model wih Two Delays GuangPing Hu and XiaoLing Li School of Mahemaics and Physics, Nanjing Universiy

More information

Expert Advice for Amateurs

Expert Advice for Amateurs Exper Advice for Amaeurs Ernes K. Lai Online Appendix - Exisence of Equilibria The analysis in his secion is performed under more general payoff funcions. Wihou aking an explici form, he payoffs of he

More information

4.6 One Dimensional Kinematics and Integration

4.6 One Dimensional Kinematics and Integration 4.6 One Dimensional Kinemaics and Inegraion When he acceleraion a( of an objec is a non-consan funcion of ime, we would like o deermine he ime dependence of he posiion funcion x( and he x -componen of

More information

Properties Of Solutions To A Generalized Liénard Equation With Forcing Term

Properties Of Solutions To A Generalized Liénard Equation With Forcing Term Applied Mahemaics E-Noes, 8(28), 4-44 c ISSN 67-25 Available free a mirror sies of hp://www.mah.nhu.edu.w/ amen/ Properies Of Soluions To A Generalized Liénard Equaion Wih Forcing Term Allan Kroopnick

More information

u(x) = e x 2 y + 2 ) Integrate and solve for x (1 + x)y + y = cos x Answer: Divide both sides by 1 + x and solve for y. y = x y + cos x

u(x) = e x 2 y + 2 ) Integrate and solve for x (1 + x)y + y = cos x Answer: Divide both sides by 1 + x and solve for y. y = x y + cos x . 1 Mah 211 Homework #3 February 2, 2001 2.4.3. y + (2/x)y = (cos x)/x 2 Answer: Compare y + (2/x) y = (cos x)/x 2 wih y = a(x)x + f(x)and noe ha a(x) = 2/x. Consequenly, an inegraing facor is found wih

More information

Guest Lectures for Dr. MacFarlane s EE3350 Part Deux

Guest Lectures for Dr. MacFarlane s EE3350 Part Deux Gues Lecures for Dr. MacFarlane s EE3350 Par Deux Michael Plane Mon., 08-30-2010 Wrie name in corner. Poin ou his is a review, so I will go faser. Remind hem o go lisen o online lecure abou geing an A

More information

4. Advanced Stability Theory

4. Advanced Stability Theory Applied Nonlinear Conrol Nguyen an ien - 4 4 Advanced Sabiliy heory he objecive of his chaper is o presen sabiliy analysis for non-auonomous sysems 41 Conceps of Sabiliy for Non-Auonomous Sysems Equilibrium

More information

Chapter 2. First Order Scalar Equations

Chapter 2. First Order Scalar Equations Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.

More information

Designing Information Devices and Systems I Spring 2019 Lecture Notes Note 17

Designing Information Devices and Systems I Spring 2019 Lecture Notes Note 17 EES 16A Designing Informaion Devices and Sysems I Spring 019 Lecure Noes Noe 17 17.1 apaciive ouchscreen In he las noe, we saw ha a capacior consiss of wo pieces on conducive maerial separaed by a nonconducive

More information

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still. Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in

More information

Solutions to Assignment 1

Solutions to Assignment 1 MA 2326 Differenial Equaions Insrucor: Peronela Radu Friday, February 8, 203 Soluions o Assignmen. Find he general soluions of he following ODEs: (a) 2 x = an x Soluion: I is a separable equaion as we

More information

Numerical Dispersion

Numerical Dispersion eview of Linear Numerical Sabiliy Numerical Dispersion n he previous lecure, we considered he linear numerical sabiliy of boh advecion and diffusion erms when approimaed wih several spaial and emporal

More information

Electrical and current self-induction

Electrical and current self-induction Elecrical and curren self-inducion F. F. Mende hp://fmnauka.narod.ru/works.hml mende_fedor@mail.ru Absrac The aricle considers he self-inducance of reacive elemens. Elecrical self-inducion To he laws of

More information

( ) a system of differential equations with continuous parametrization ( T = R + These look like, respectively:

( ) a system of differential equations with continuous parametrization ( T = R + These look like, respectively: XIII. DIFFERENCE AND DIFFERENTIAL EQUATIONS Ofen funcions, or a sysem of funcion, are paramerized in erms of some variable, usually denoed as and inerpreed as ime. The variable is wrien as a funcion of

More information

Morning Time: 1 hour 30 minutes Additional materials (enclosed):

Morning Time: 1 hour 30 minutes Additional materials (enclosed): ADVANCED GCE 78/0 MATHEMATICS (MEI) Differenial Equaions THURSDAY JANUARY 008 Morning Time: hour 30 minues Addiional maerials (enclosed): None Addiional maerials (required): Answer Bookle (8 pages) Graph

More information

EXERCISES FOR SECTION 1.5

EXERCISES FOR SECTION 1.5 1.5 Exisence and Uniqueness of Soluions 43 20. 1 v c 21. 1 v c 1 2 4 6 8 10 1 2 2 4 6 8 10 Graph of approximae soluion obained using Euler s mehod wih = 0.1. Graph of approximae soluion obained using Euler

More information

Sliding Mode Extremum Seeking Control for Linear Quadratic Dynamic Game

Sliding Mode Extremum Seeking Control for Linear Quadratic Dynamic Game Sliding Mode Exremum Seeking Conrol for Linear Quadraic Dynamic Game Yaodong Pan and Ümi Özgüner ITS Research Group, AIST Tsukuba Eas Namiki --, Tsukuba-shi,Ibaraki-ken 5-856, Japan e-mail: pan.yaodong@ais.go.jp

More information

Predator - Prey Model Trajectories and the nonlinear conservation law

Predator - Prey Model Trajectories and the nonlinear conservation law Predaor - Prey Model Trajecories and he nonlinear conservaion law James K. Peerson Deparmen of Biological Sciences and Deparmen of Mahemaical Sciences Clemson Universiy Ocober 28, 213 Ouline Drawing Trajecories

More information

arxiv: v1 [math.ca] 15 Nov 2016

arxiv: v1 [math.ca] 15 Nov 2016 arxiv:6.599v [mah.ca] 5 Nov 26 Counerexamples on Jumarie s hree basic fracional calculus formulae for non-differeniable coninuous funcions Cheng-shi Liu Deparmen of Mahemaics Norheas Peroleum Universiy

More information

Class Meeting # 10: Introduction to the Wave Equation

Class Meeting # 10: Introduction to the Wave Equation MATH 8.5 COURSE NOTES - CLASS MEETING # 0 8.5 Inroducion o PDEs, Fall 0 Professor: Jared Speck Class Meeing # 0: Inroducion o he Wave Equaion. Wha is he wave equaion? The sandard wave equaion for a funcion

More information

The motions of the celt on a horizontal plane with viscous friction

The motions of the celt on a horizontal plane with viscous friction The h Join Inernaional Conference on Mulibody Sysem Dynamics June 8, 18, Lisboa, Porugal The moions of he cel on a horizonal plane wih viscous fricion Maria A. Munisyna 1 1 Moscow Insiue of Physics and

More information

Ordinary Differential Equations

Ordinary Differential Equations Ordinary Differenial Equaions 5. Examples of linear differenial equaions and heir applicaions We consider some examples of sysems of linear differenial equaions wih consan coefficiens y = a y +... + a

More information

= ( ) ) or a system of differential equations with continuous parametrization (T = R

= ( ) ) or a system of differential equations with continuous parametrization (T = R XIII. DIFFERENCE AND DIFFERENTIAL EQUATIONS Ofen funcions, or a sysem of funcion, are paramerized in erms of some variable, usually denoed as and inerpreed as ime. The variable is wrien as a funcion of

More information

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t... Mah 228- Fri Mar 24 5.6 Marix exponenials and linear sysems: The analogy beween firs order sysems of linear differenial equaions (Chaper 5) and scalar linear differenial equaions (Chaper ) is much sronger

More information

Signal and System (Chapter 3. Continuous-Time Systems)

Signal and System (Chapter 3. Continuous-Time Systems) Signal and Sysem (Chaper 3. Coninuous-Time Sysems) Prof. Kwang-Chun Ho kwangho@hansung.ac.kr Tel: 0-760-453 Fax:0-760-4435 1 Dep. Elecronics and Informaion Eng. 1 Nodes, Branches, Loops A nework wih b

More information

Stability and Convergence in Adaptive Systems *

Stability and Convergence in Adaptive Systems * Sabiliy and Convergence in Adapive Sysems * Margarea Sefanovic, Rengrong Wang, and Michael G. Safonov Universiy of Souhern California, Los Angeles, CA 989-563 Absrac - Sufficien condiions for adapive conrol

More information

The Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales

The Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales Advances in Dynamical Sysems and Applicaions. ISSN 0973-5321 Volume 1 Number 1 (2006, pp. 103 112 c Research India Publicaions hp://www.ripublicaion.com/adsa.hm The Asympoic Behavior of Nonoscillaory Soluions

More information

Math 334 Test 1 KEY Spring 2010 Section: 001. Instructor: Scott Glasgow Dates: May 10 and 11.

Math 334 Test 1 KEY Spring 2010 Section: 001. Instructor: Scott Glasgow Dates: May 10 and 11. 1 Mah 334 Tes 1 KEY Spring 21 Secion: 1 Insrucor: Sco Glasgow Daes: Ma 1 and 11. Do NOT wrie on his problem saemen bookle, excep for our indicaion of following he honor code jus below. No credi will be

More information

Chapter 8 The Complete Response of RL and RC Circuits

Chapter 8 The Complete Response of RL and RC Circuits Chaper 8 The Complee Response of RL and RC Circuis Seoul Naional Universiy Deparmen of Elecrical and Compuer Engineering Wha is Firs Order Circuis? Circuis ha conain only one inducor or only one capacior

More information

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x WEEK-3 Reciaion PHYS 131 Ch. 3: FOC 1, 3, 4, 6, 14. Problems 9, 37, 41 & 71 and Ch. 4: FOC 1, 3, 5, 8. Problems 3, 5 & 16. Feb 8, 018 Ch. 3: FOC 1, 3, 4, 6, 14. 1. (a) The horizonal componen of he projecile

More information

Linear Response Theory: The connection between QFT and experiments

Linear Response Theory: The connection between QFT and experiments Phys540.nb 39 3 Linear Response Theory: The connecion beween QFT and experimens 3.1. Basic conceps and ideas Q: How do we measure he conduciviy of a meal? A: we firs inroduce a weak elecric field E, and

More information

6.2 Transforms of Derivatives and Integrals.

6.2 Transforms of Derivatives and Integrals. SEC. 6.2 Transforms of Derivaives and Inegrals. ODEs 2 3 33 39 23. Change of scale. If l( f ()) F(s) and c is any 33 45 APPLICATION OF s-shifting posiive consan, show ha l( f (c)) F(s>c)>c (Hin: In Probs.

More information

Four Generations of Higher Order Sliding Mode Controllers. L. Fridman Universidad Nacional Autonoma de Mexico Aussois, June, 10th, 2015

Four Generations of Higher Order Sliding Mode Controllers. L. Fridman Universidad Nacional Autonoma de Mexico Aussois, June, 10th, 2015 Four Generaions of Higher Order Sliding Mode Conrollers L. Fridman Universidad Nacional Auonoma de Mexico Aussois, June, 1h, 215 Ouline 1 Generaion 1:Two Main Conceps of Sliding Mode Conrol 2 Generaion

More information

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004 ODEs II, Lecure : Homogeneous Linear Sysems - I Mike Raugh March 8, 4 Inroducion. In he firs lecure we discussed a sysem of linear ODEs for modeling he excreion of lead from he human body, saw how o ransform

More information

Notes on Kalman Filtering

Notes on Kalman Filtering Noes on Kalman Filering Brian Borchers and Rick Aser November 7, Inroducion Daa Assimilaion is he problem of merging model predicions wih acual measuremens of a sysem o produce an opimal esimae of he curren

More information

SOLUTIONS TO ECE 3084

SOLUTIONS TO ECE 3084 SOLUTIONS TO ECE 384 PROBLEM 2.. For each sysem below, specify wheher or no i is: (i) memoryless; (ii) causal; (iii) inverible; (iv) linear; (v) ime invarian; Explain your reasoning. If he propery is no

More information

An introduction to the theory of SDDP algorithm

An introduction to the theory of SDDP algorithm An inroducion o he heory of SDDP algorihm V. Leclère (ENPC) Augus 1, 2014 V. Leclère Inroducion o SDDP Augus 1, 2014 1 / 21 Inroducion Large scale sochasic problem are hard o solve. Two ways of aacking

More information

Non-uniform circular motion *

Non-uniform circular motion * OpenSax-CNX module: m14020 1 Non-uniform circular moion * Sunil Kumar Singh This work is produced by OpenSax-CNX and licensed under he Creaive Commons Aribuion License 2.0 Wha do we mean by non-uniform

More information

KINEMATICS IN ONE DIMENSION

KINEMATICS IN ONE DIMENSION KINEMATICS IN ONE DIMENSION PREVIEW Kinemaics is he sudy of how hings move how far (disance and displacemen), how fas (speed and velociy), and how fas ha how fas changes (acceleraion). We say ha an objec

More information

T L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB

T L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB Elecronic Companion EC.1. Proofs of Technical Lemmas and Theorems LEMMA 1. Le C(RB) be he oal cos incurred by he RB policy. Then we have, T L E[C(RB)] 3 E[Z RB ]. (EC.1) Proof of Lemma 1. Using he marginal

More information

2.160 System Identification, Estimation, and Learning. Lecture Notes No. 8. March 6, 2006

2.160 System Identification, Estimation, and Learning. Lecture Notes No. 8. March 6, 2006 2.160 Sysem Idenificaion, Esimaion, and Learning Lecure Noes No. 8 March 6, 2006 4.9 Eended Kalman Filer In many pracical problems, he process dynamics are nonlinear. w Process Dynamics v y u Model (Linearized)

More information

Random Walk with Anti-Correlated Steps

Random Walk with Anti-Correlated Steps Random Walk wih Ani-Correlaed Seps John Noga Dirk Wagner 2 Absrac We conjecure he expeced value of random walks wih ani-correlaed seps o be exacly. We suppor his conjecure wih 2 plausibiliy argumens and

More information

Final Spring 2007

Final Spring 2007 .615 Final Spring 7 Overview The purpose of he final exam is o calculae he MHD β limi in a high-bea oroidal okamak agains he dangerous n = 1 exernal ballooning-kink mode. Effecively, his corresponds o

More information

Optimality Conditions for Unconstrained Problems

Optimality Conditions for Unconstrained Problems 62 CHAPTER 6 Opimaliy Condiions for Unconsrained Problems 1 Unconsrained Opimizaion 11 Exisence Consider he problem of minimizing he funcion f : R n R where f is coninuous on all of R n : P min f(x) x

More information

Multi-Layer Switching Control

Multi-Layer Switching Control 5 American Conrol Conference June 8-1, 5. Porland, OR, USA FrC7.3 Muli-Layer Swiching Conrol Idin Karuei, Nader Meskin, Amir G. Aghdam Deparmen of Elecrical and Compuer Engineering, Concordia Universiy

More information

Unsteady Flow Problems

Unsteady Flow Problems School of Mechanical Aerospace and Civil Engineering Unseady Flow Problems T. J. Craf George Begg Building, C41 TPFE MSc CFD-1 Reading: J. Ferziger, M. Peric, Compuaional Mehods for Fluid Dynamics H.K.

More information

Matrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality

Matrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality Marix Versions of Some Refinemens of he Arihmeic-Geomeric Mean Inequaliy Bao Qi Feng and Andrew Tonge Absrac. We esablish marix versions of refinemens due o Alzer ], Carwrigh and Field 4], and Mercer 5]

More information

LAPLACE TRANSFORM AND TRANSFER FUNCTION

LAPLACE TRANSFORM AND TRANSFER FUNCTION CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION Professor Dae Ryook Yang Spring 2018 Dep. of Chemical and Biological Engineering 5-1 Road Map of he Lecure V Laplace Transform and Transfer funcions

More information

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is UNIT IMPULSE RESPONSE, UNIT STEP RESPONSE, STABILITY. Uni impulse funcion (Dirac dela funcion, dela funcion) rigorously defined is no sricly a funcion, bu disribuion (or measure), precise reamen requires

More information

Two Popular Bayesian Estimators: Particle and Kalman Filters. McGill COMP 765 Sept 14 th, 2017

Two Popular Bayesian Estimators: Particle and Kalman Filters. McGill COMP 765 Sept 14 th, 2017 Two Popular Bayesian Esimaors: Paricle and Kalman Filers McGill COMP 765 Sep 14 h, 2017 1 1 1, dx x Bel x u x P x z P Recall: Bayes Filers,,,,,,, 1 1 1 1 u z u x P u z u x z P Bayes z = observaion u =

More information

Boundedness and Exponential Asymptotic Stability in Dynamical Systems with Applications to Nonlinear Differential Equations with Unbounded Terms

Boundedness and Exponential Asymptotic Stability in Dynamical Systems with Applications to Nonlinear Differential Equations with Unbounded Terms Advances in Dynamical Sysems and Applicaions. ISSN 0973-531 Volume Number 1 007, pp. 107 11 Research India Publicaions hp://www.ripublicaion.com/adsa.hm Boundedness and Exponenial Asympoic Sabiliy in Dynamical

More information

CHAPTER 2 Signals And Spectra

CHAPTER 2 Signals And Spectra CHAPER Signals And Specra Properies of Signals and Noise In communicaion sysems he received waveform is usually caegorized ino he desired par conaining he informaion, and he undesired par. he desired par

More information

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3. Mah Rahman Exam Review Soluions () Consider he IVP: ( 4)y 3y + 4y = ; y(3) = 0, y (3) =. (a) Please deermine he longes inerval for which he IVP is guaraneed o have a unique soluion. Soluion: The disconinuiies

More information

IB Physics Kinematics Worksheet

IB Physics Kinematics Worksheet IB Physics Kinemaics Workshee Wrie full soluions and noes for muliple choice answers. Do no use a calculaor for muliple choice answers. 1. Which of he following is a correc definiion of average acceleraion?

More information

Summary of shear rate kinematics (part 1)

Summary of shear rate kinematics (part 1) InroToMaFuncions.pdf 4 CM465 To proceed o beer-designed consiuive equaions, we need o know more abou maerial behavior, i.e. we need more maerial funcions o predic, and we need measuremens of hese maerial

More information

Comments on Window-Constrained Scheduling

Comments on Window-Constrained Scheduling Commens on Window-Consrained Scheduling Richard Wes Member, IEEE and Yuing Zhang Absrac This shor repor clarifies he behavior of DWCS wih respec o Theorem 3 in our previously published paper [1], and describes

More information

Vanishing Viscosity Method. There are another instructive and perhaps more natural discontinuous solutions of the conservation law

Vanishing Viscosity Method. There are another instructive and perhaps more natural discontinuous solutions of the conservation law Vanishing Viscosiy Mehod. There are anoher insrucive and perhaps more naural disconinuous soluions of he conservaion law (1 u +(q(u x 0, he so called vanishing viscosiy mehod. This mehod consiss in viewing

More information

STABILITY OF RESET SWITCHING SYSTEMS

STABILITY OF RESET SWITCHING SYSTEMS STABILITY OF RESET SWITCHING SYSTEMS J.P. Paxman,G.Vinnicombe Cambridge Universiy Engineering Deparmen, UK, jpp7@cam.ac.uk Keywords: Sabiliy, swiching sysems, bumpless ransfer, conroller realizaion. Absrac

More information

Georey E. Hinton. University oftoronto. Technical Report CRG-TR February 22, Abstract

Georey E. Hinton. University oftoronto.   Technical Report CRG-TR February 22, Abstract Parameer Esimaion for Linear Dynamical Sysems Zoubin Ghahramani Georey E. Hinon Deparmen of Compuer Science Universiy oftorono 6 King's College Road Torono, Canada M5S A4 Email: zoubin@cs.orono.edu Technical

More information

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature On Measuring Pro-Poor Growh 1. On Various Ways of Measuring Pro-Poor Growh: A Shor eview of he Lieraure During he pas en years or so here have been various suggesions concerning he way one should check

More information

Exponential Weighted Moving Average (EWMA) Chart Under The Assumption of Moderateness And Its 3 Control Limits

Exponential Weighted Moving Average (EWMA) Chart Under The Assumption of Moderateness And Its 3 Control Limits DOI: 0.545/mjis.07.5009 Exponenial Weighed Moving Average (EWMA) Char Under The Assumpion of Moderaeness And Is 3 Conrol Limis KALPESH S TAILOR Assisan Professor, Deparmen of Saisics, M. K. Bhavnagar Universiy,

More information

Math 333 Problem Set #2 Solution 14 February 2003

Math 333 Problem Set #2 Solution 14 February 2003 Mah 333 Problem Se #2 Soluion 14 February 2003 A1. Solve he iniial value problem dy dx = x2 + e 3x ; 2y 4 y(0) = 1. Soluion: This is separable; we wrie 2y 4 dy = x 2 + e x dx and inegrae o ge The iniial

More information

8. Basic RL and RC Circuits

8. Basic RL and RC Circuits 8. Basic L and C Circuis This chaper deals wih he soluions of he responses of L and C circuis The analysis of C and L circuis leads o a linear differenial equaion This chaper covers he following opics

More information

MATH 4330/5330, Fourier Analysis Section 6, Proof of Fourier s Theorem for Pointwise Convergence

MATH 4330/5330, Fourier Analysis Section 6, Proof of Fourier s Theorem for Pointwise Convergence MATH 433/533, Fourier Analysis Secion 6, Proof of Fourier s Theorem for Poinwise Convergence Firs, some commens abou inegraing periodic funcions. If g is a periodic funcion, g(x + ) g(x) for all real x,

More information

Lecture 20: Riccati Equations and Least Squares Feedback Control

Lecture 20: Riccati Equations and Least Squares Feedback Control 34-5 LINEAR SYSTEMS Lecure : Riccai Equaions and Leas Squares Feedback Conrol 5.6.4 Sae Feedback via Riccai Equaions A recursive approach in generaing he marix-valued funcion W ( ) equaion for i for he

More information

Differential Equations

Differential Equations Mah 21 (Fall 29) Differenial Equaions Soluion #3 1. Find he paricular soluion of he following differenial equaion by variaion of parameer (a) y + y = csc (b) 2 y + y y = ln, > Soluion: (a) The corresponding

More information

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow 1 KEY Mah 4 Miderm I Fall 8 secions 1 and Insrucor: Sco Glasgow Please do NOT wrie on his eam. No credi will be given for such work. Raher wrie in a blue book, or on our own paper, preferabl engineering

More information

Chapter 15: Phenomena. Chapter 15 Chemical Kinetics. Reaction Rates. Reaction Rates R P. Reaction Rates. Rate Laws

Chapter 15: Phenomena. Chapter 15 Chemical Kinetics. Reaction Rates. Reaction Rates R P. Reaction Rates. Rate Laws Chaper 5: Phenomena Phenomena: The reacion (aq) + B(aq) C(aq) was sudied a wo differen emperaures (98 K and 35 K). For each emperaure he reacion was sared by puing differen concenraions of he 3 species

More information

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes Some common engineering funcions 2.7 Inroducion This secion provides a caalogue of some common funcions ofen used in Science and Engineering. These include polynomials, raional funcions, he modulus funcion

More information

Two Coupled Oscillators / Normal Modes

Two Coupled Oscillators / Normal Modes Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own

More information

Effects of Coordinate Curvature on Integration

Effects of Coordinate Curvature on Integration Effecs of Coordinae Curvaure on Inegraion Chrisopher A. Lafore clafore@gmail.com Absrac In his paper, he inegraion of a funcion over a curved manifold is examined in he case where he curvaure of he manifold

More information

Oscillation of an Euler Cauchy Dynamic Equation S. Huff, G. Olumolode, N. Pennington, and A. Peterson

Oscillation of an Euler Cauchy Dynamic Equation S. Huff, G. Olumolode, N. Pennington, and A. Peterson PROCEEDINGS OF THE FOURTH INTERNATIONAL CONFERENCE ON DYNAMICAL SYSTEMS AND DIFFERENTIAL EQUATIONS May 4 7, 00, Wilmingon, NC, USA pp 0 Oscillaion of an Euler Cauchy Dynamic Equaion S Huff, G Olumolode,

More information

Problemas das Aulas Práticas

Problemas das Aulas Práticas Mesrado Inegrado em Engenharia Elecroécnica e de Compuadores Conrolo em Espaço de Esados Problemas das Aulas Práicas J. Miranda Lemos Fevereiro de 3 Translaed o English by José Gaspar, 6 J. M. Lemos, IST

More information

arxiv: v1 [math.gm] 4 Nov 2018

arxiv: v1 [math.gm] 4 Nov 2018 Unpredicable Soluions of Linear Differenial Equaions Mara Akhme 1,, Mehme Onur Fen 2, Madina Tleubergenova 3,4, Akylbek Zhamanshin 3,4 1 Deparmen of Mahemaics, Middle Eas Technical Universiy, 06800, Ankara,

More information

Recursive Least-Squares Fixed-Interval Smoother Using Covariance Information based on Innovation Approach in Linear Continuous Stochastic Systems

Recursive Least-Squares Fixed-Interval Smoother Using Covariance Information based on Innovation Approach in Linear Continuous Stochastic Systems 8 Froniers in Signal Processing, Vol. 1, No. 1, July 217 hps://dx.doi.org/1.2266/fsp.217.112 Recursive Leas-Squares Fixed-Inerval Smooher Using Covariance Informaion based on Innovaion Approach in Linear

More information

Problem Set on Differential Equations

Problem Set on Differential Equations Problem Se on Differenial Equaions 1. Solve he following differenial equaions (a) x () = e x (), x () = 3/ 4. (b) x () = e x (), x (1) =. (c) xe () = + (1 x ()) e, x () =.. (An asse marke model). Le p()

More information

Hybrid Control and Switched Systems. Lecture #3 What can go wrong? Trajectories of hybrid systems

Hybrid Control and Switched Systems. Lecture #3 What can go wrong? Trajectories of hybrid systems Hybrid Conrol and Swiched Sysems Lecure #3 Wha can go wrong? Trajecories of hybrid sysems João P. Hespanha Universiy of California a Sana Barbara Summary 1. Trajecories of hybrid sysems: Soluion o a hybrid

More information

1. VELOCITY AND ACCELERATION

1. VELOCITY AND ACCELERATION 1. VELOCITY AND ACCELERATION 1.1 Kinemaics Equaions s = u + 1 a and s = v 1 a s = 1 (u + v) v = u + as 1. Displacemen-Time Graph Gradien = speed 1.3 Velociy-Time Graph Gradien = acceleraion Area under

More information

A New Perturbative Approach in Nonlinear Singularity Analysis

A New Perturbative Approach in Nonlinear Singularity Analysis Journal of Mahemaics and Saisics 7 (: 49-54, ISSN 549-644 Science Publicaions A New Perurbaive Approach in Nonlinear Singulariy Analysis Ta-Leung Yee Deparmen of Mahemaics and Informaion Technology, The

More information

18 Biological models with discrete time

18 Biological models with discrete time 8 Biological models wih discree ime The mos imporan applicaions, however, may be pedagogical. The elegan body of mahemaical heory peraining o linear sysems (Fourier analysis, orhogonal funcions, and so

More information

R.#W.#Erickson# Department#of#Electrical,#Computer,#and#Energy#Engineering# University#of#Colorado,#Boulder#

R.#W.#Erickson# Department#of#Electrical,#Computer,#and#Energy#Engineering# University#of#Colorado,#Boulder# .#W.#Erickson# Deparmen#of#Elecrical,#Compuer,#and#Energy#Engineering# Universiy#of#Colorado,#Boulder# Chaper 2 Principles of Seady-Sae Converer Analysis 2.1. Inroducion 2.2. Inducor vol-second balance,

More information

Inventory Control of Perishable Items in a Two-Echelon Supply Chain

Inventory Control of Perishable Items in a Two-Echelon Supply Chain Journal of Indusrial Engineering, Universiy of ehran, Special Issue,, PP. 69-77 69 Invenory Conrol of Perishable Iems in a wo-echelon Supply Chain Fariborz Jolai *, Elmira Gheisariha and Farnaz Nojavan

More information