Optimal control of connected vehicle systems

Size: px
Start display at page:

Download "Optimal control of connected vehicle systems"

Transcription

1 53rd IEEE Conference on Decision and Conrol December 15-17, 214. Los Angeles, California, USA Opimal conrol of conneced vehicle sysems Jin I. Ge and Gábor Orosz Absrac In his paper, linear quadraic racking LQ) is used o opimize he conrol gains for conneced cruise conrol CCC). We consider a vehicle sring where he CCC vehicle a he ail receives posiion and velociy signals hrough wireless vehicle-o-vehicle V2V) communicaion from oher vehicles ahead ha are no equipped wih CCC). An opimal feedback law is obained by minimizing a cos funcion defined by headway and velociy errors and he acceleraion of he CCC vehicle on an infinie horizon. We show ha he feedback gains can be obained recursively as signals from vehicles farher ahead become available, and ha he gains decay exponenially wih he number of cars beween he source of he signal and he CCC vehicle. he effecs of he cos funcion on he head-oail sring sabiliy are invesigaed and he robusness agains variaions in human parameers is esed. he analyical resuls are verified by numerical simulaions. I. INRODUCION Conneced cruise conrol CCC) has been proposed o mainain smooh raffic flow in heerogeneous conneced vehicle sysems by exploiing vehicle-o-vehicle V2V) communicaion [1]. he CCC conroller receives informaion abou he moion of muliple vehicles ahead, and acuaes he vehicle or assiss he driver based on hese signals. he influence of conneciviy srucures, signal ypes, packe drops, and communicaion delays on he longiudinal moion of vehicular chains ha include CCC vehicles has been invesigaed [2] [5]. Our goal here is o opimize he feedback gains in order o maximize he benefi of conneciviy and reduce he complexiy of uning gains individually in large sysems; see [6], [7] for iniial aemps using simple configuraions. Moreover, he design parameers should be chosen so ha addiional performance requiremenssuch as sring sabiliy) are saisfied. In his paper we opimize he gains of a CCC vehicle ha receives posiion and velociy informaion from muliple human-driven vehicles ahead. he goal of opimizaion is o obain a CCC conroller ha ensures he sabiliy of uniform raffic flow i.e. he aenuaion of perurbaions along he vehicular chain), while minimizing velociy and headway error and acceleraion of he CCC vehicle. his problem is solved by using linear quadraic racking LQ) wih design parameers being he weighs on he error erms and he acceleraion erm in he cos funcion. We show ha he gains of he opimized conroller follow he spaial causaliy of raffic sysems: informaion from vehicles farher downsream have less influence on he CCC vehicle and does *his work was suppored by he Naional Science Foundaion Award Number ) Jin I. Ge and Gábor Orosz are wih he Deparmen of Mechanical Engineering, Universiy of Michigan, Ann Arbor, Michigan 4819, USA. Corresponding gejin@umich.edu, orosz@umich.edu no change he feedback laws on signals from closer vehicles. he opimal gains are deermined by he weighs used in he opimizaion design parameers) and he driver parameers of oher vehicles. he range of design parameers ensuring head-o-ail sring sabiliy, and heir robusness agains variaions of driver parameers are also demonsraed. Finally, simulaions are performed o demonsrae he effeciveness of he opimal design. II. CONNECED CAR-FOLLOWING MODELS We consider a chain of n+1 vehicles raveling on a single lane as shown in Fig. 1a). he ail vehicle he las vehicle of he chain) implemens a CCC algorihm using posiion and velociy signals received hrough V2V communicaion from n preceding vehicles, while oher vehicles are human driven and only ransmi informaion abou heir moion. he dynamics of he CCC vehicle is modeled by ḣ 1 ) = v 2 ) v 1 ), v 1 ) = u), 1) where he do sands for differeniaion wih respec o ime, h 1 is he headway i.e., he bumper-o-bumper disance beween he CCC vehicle and he vehicle immediaely ahead), and v 1 is he velociy of he CCC vehicle; see Fig. 1a). Finally, u) is he conrol inpu ha will be designed using LQ based on he velociy and headway of oher vehicles he laer obained from posiion informaion). For simpliciy, we consider ha vehicles i = 2,..., n are idenical and are described by he car-following model ḣ i ) = v i+1 ) v i ), v i ) = α V h i )) v i ) ) + β v i+1 ) v i ) ), 2) ha can be obained as a simplificaion of he physics-based model presened in [1]. Here h i and v i denoe he headway and velociy of vehicle i; see Fig. 1a). he firs erm in he second equaion represens he driver s inenion o drive a a disance-dependen velociy given by V h)), while he second erm represens he driver s aim o mach he velociy o ha of he vehicle immediaely ahead. he corresponding gains are denoed α and β. We remark ha he proposed algorihm can also be applied in case of non-idenical drivers as well. he desired velociy in 2) is deermined by he range policy if h h s, v V h) = max 2 1 cos ) ) π h hs h go h s if h s < h < h go, v max if h h go, 3) /14/$ IEEE 417

2 i = 2,..., n, and linearize 2) abou he equilibrium 4): h i ) = ṽ i+1 ) ṽ i ), ṽ i ) = α f hi ) ṽ i ) ) + β hi ), 5) Fig. 1. a): A chain of n + 1 vehicles wih a CCC vehicle a he ail receiving signals from oher vehicles via V2V communicaion. b): he nonlinear range policy 3) used in his paper. which is shown in Fig. 1b). he desired velociy is zero for small headways h h s ) and equal o he maximum speed v max for large headways h h go ). Beween hese, i increases wih he headway monoonically. o ensure smooh longiudinal dynamics, he funcion 3) and is derivaive are chosen o be coninuous a h s and h go. Here we consider v max = 3 [m/s], h s = 5 [m], h go = 35 [m] ha corresponds o realisic raffic daa [8]. Many oher range policies may be chosen, bu he qualiaive dynamics remain similar if he above characerisics are kep [2], [9], [1]. We remark ha model 1,2,3) may no adequaely describe longiudinal dynamics when vehicles are driven near physical limis, e.g., ire force sauraion due o emergency braking or severe maneuvering. In his paper we focus on he poenials of wireless conneciviy, while undersanding such limiaions is lef for fuure research. III. LINEAR QUADRAIC RACKING OF UNIFORM RAFFIC FLOW In his secion, he opimal conrol problem under disurbance is formulaed, where he CCC vehicle is racking he uniform flow sae. he cos funcion is consruced in order o minimize he headway and velociy errors and he acceleraion of he CCC vehicle. he soluion gives he gains for he CCC vehicle wih respec o he headways and velociies of oher vehicles. he dynamics of he conneced vehicle sysem 1,2) is invesigaed in he viciniy of an equilibrium where all vehicles ravel wih he same consan velociy and mainain consan headways. While he equilibrium velociy v is deermined by he head vehicle he firs vehicle in he chain), he equilibrium headway h i is obained for each non-ccc vehicle using a range policy v = V i h i ), i = 2,..., n. When considering idenical range policies cf. 3)), he vehicles are equidisan and we obain he uniform flow equilibrium h i ) h, v i ) v = V h ), 4) for i = 2,..., n. We define headway perurbaions h i ) = h i ) h and velociy perurbaions ṽ i ) = v i ) v, for i = 2,..., n. Here f = V h ) is he derivaive of he range policy a he equilibrium and he corresponding ime headway is h = 1/f. In his paper, we use h, v ) = 2 [m], 15 [m/s]), which resuls in he maximum slope f = π/2 [1/s] corresponding o he minimum ime headway h = 2/π.64 [s]; cf. 3) wih v max = 3 [m/s], h s = 5 [m], and h go = 35 [m]. For he CCC vehicle, we define headway perurbaion h 1 ) = h 1 ) h 1 and velociy perurbaion ṽ 1 ) = v 1 ) v. hen 1) yields he linearized dynamics h 1 ) = ṽ 2 ) ṽ 1 ), ṽ 1 ) = u). 6) Le s define he sae x = [ h 1, ṽ 1,..., h n, ṽ n ] R 2n, and wrie dynamics 5,6) in he form ẋ) = Ax) + Bu) + Dṽ n+1 ), 7) where u) is he inpu, ṽ n+1 ) is he disurbance, and he coefficien marices ake he form A 1 A 2 B 1 A 3 A 4 A =......, B =., D =. A 3 A 4 A 3 D 1 8) where he block marices are given by [ ] [ ] [ ] A 1 =, A 2 =, A 4 =, β [ ] [ 1 A 3 = αf, B α β 1 =, D 1] 1 = [ 1 β]. 9) Since our goal is o rack he uniform flow equilibrium x cf. 4)) under velociy disurbance ṽ n+1 ) from he head vehicle, we minimize he cos funcion J τ u, x) = τ x )Qx) + r u 2 ) ) d. 1) he firs erm corresponds o he variaion of he headways and velociies which we call racking errors, he second erm corresponds o he magniude of he CCC vehicle s acceleraion, and τ is he ime horizon we use τ laer). he weigh marix Q is chosen o be diagonal, ha is, Q = diag [q 1, q 2,..., q 2n 1, q 2n ]), 11) where q 2i 1 and q 2i are he weighs on he headway and velociy errors for vehicle i, respecively. Since only vehicle 1 has he CCC conroller, hi, ṽ i, i = 2,..., n are no conrollable and he choice of q 2i 1, q 2i, i = 2,..., n does no influence he opimal conrol inpu. hus we se q 2i 1 = q 2i =, i = 2,..., n. 418

3 3 P a) w b) Fig. 2. a): he soluion P) of 13) wih A = 1, B = 1, Q = 1, r = 1, τ = 3 [s]. b): he corresponding soluion w) of 14) wih perurbaion ṽ n+1 ) = sin2), using P) shown in panel a) red dashed curve) and using he approximaion P = P) blue solid curve). Based on he linear quadraic racking heory [11], he soluion of he opimal conrol problem 7,8,1,11) is given by u) = 1 r B P)x) + w) ), 12) where P) R 2n 2n is a symmeric, posiive definie marix ha saisfies he Riccai differenial equaion Ṗ) = 1 r P)BB P) A P) P)A Q, 13) wih end boundary condiion Pτ) =, while w) R 2n is he soluion of ẇ) = A 1 r BB P) ) w) P)Dṽn+1 ), 14) wih end boundary condiion wτ) =. Noe ha wihou disurbance, i.e., ṽ n+1 ), he LQ problem 7,1) simplifies o an LQR problem wih inpu u) = 1 r B P)x), where P) is he soluion of 13). According o [11], when he sysem 7,8) is sabilizable, hen 13) has uniformly bounded soluion. Moreover, in he infinie-ime horizon τ, he soluion P) can be approximaed by a consan marix given by he algebraic Riccai equaion A P + PA + Q 1 r PBB P =. 15) From now on we use P wihou o denoe his imeindependen marix. Subsiuing P ino 14), and considering ha A 1 r BB P) is Hurwiz, he boundary value problem can be viewed as an iniial value problem wih a special iniial condiion ha eliminaes he ransiens. A simple example is given in Fig. 2 for he simplified problem ẋ) = x) + u) + ṽ n+1 ), x, u R using he weighs Q = 1, r = 1. Panel a) shows ha P) is approximaely consan when τ. Panel b) shows ha using consan P = P) insead of P) in 14) only influences w) near ime τ. hus for large τ, P) can be approximaed by P = P) and he laer can be used o calculae w). Le us define he noaion P = Π Π 1n..... Π n1... Π nn, 16) where Π ij = Π ji R 2 2, i, j = 1,..., n. Considering infinie ime horizon τ, he feedback law 12) gives he acceleraion of he CCC vehicle in 6) as u) = n αi hi ) + β i ṽ i ) ) 1 r w 2), 17) i=1 where w 2 ) denoes he second elemen of w) and he gains α i and β i are given by α i = 1 r Π 1i[2, 1], β i = 1 r Π 1i[2, 2], 18) for i = 1,..., n, where [k, l] represens he elemen in he k h row and l h column. Due o he paricular form of he coefficien marices 8), for i = j = 1, 15) yields Π 11 B 1 Π 11 Π 11 A 1 A 1 Π 11 diag[q 1, q 2 ]) =, 19) where B 1 = 1 r B 1B 1 ; c.f. 9). Moreover, using 15) i can be shown ha he firs row of block marices in P cf. 16)) saisfy he recursive equaions A 1 Π 11 B 1 ) Π12 + Π 12 A 3 = Π 11 A 2, A 1 Π 11 B 1 ) Π1j + Π 1j A 3 = Π 1j 1) A 4, 2) where j = 3,..., n. Also, he second row of block marices saisfy he recursive equaions A 3 Π 22 + Π 22 A 3 = Π 21 B 1 Π 12 Π 12A 2 A 2 Π 12, A 3 Π 2j + Π 2j A 3 = Π 21 B 1 Π 1j Π 2j 1) A 4 A 2 Π 1j, 21) where j = 3,..., n and Π 21 = Π 12. For he remaining n 2 rows of block marices, we obain he recursive equaions A 3 Π ij +Π ij A 3 = Π i1 B 1 Π 1j Π ij 1) A 4 A 4 Π i 1)j, 22) where i = 3,..., n, j = i,..., n and Π i1 = Π 1i. hus, he soluion of he Riccai equaion 15) can be obained by solving 19,2,21,22) consecuively. In he physically realisic case q 1, q 2, r >, he only feasible soluion of 19) is given by q Π 11 = 1 q q1 3r q 1 r q 1 r q 2 r + 2, 23) q 1 r 3 and hus 18) yields α 1 = q 1 /r, β 1 = q 2 /r + 2 q 1 /r. 24) Moreover, according o 18), he gains α i, β i obained from 2) can be rewrien as vecπ 12 ) = M vecπ 11 ), vecπ 1i ) = MvecΠ 1i 1) ), 25) 419

4 .8 α i i.8 β i Fig. 3. he opimized headway and velociy gains α i, β i, i = 2,..., n of he CCC vehicle in a n + 1)-car plaoon for n = 5 red circles) and for n = 1 blue crosses). where vecπ 1i ), i = 3,..., n, is a vecor obained by wriing columns of Π 1i ino a vecor and we have M = I A ) 1 Π 11 B 1 + A 1A 3 I) 4 I), M = I A ) 1 Π 11 B 1 + A 1A 3 I) 2 I). herefore, i 26) vecπ 1i ) = M i 2 M vecπ 11 ) 27) is a map beween Π 11 and Π 1i, for i = 2,..., n, and hus α i, β i, i = 2,..., n can be obained as funcions of α 1, β 1, α, β, f ; cf. 9,18,26). Moreover, 27) indicaes ha he feedback gains of vehicle i are deermined by he carfollowing dynamics of vehicles 1, 2,..., i 1 and are no influenced by he dynamics of vehicles i+1,..., n+1. his propery means ha our CCC design is scalable, since he values of gains can be kep consan regardless how many vehicles ahead are moniored. Finally, we noe ha 21,22) are only needed o obain w); cf. 14). One may show ha he eigenvalues of M in 27) are inside he uni circle for human gains α >, β >. herefore, 27) is a conracing map in realisic scenarios and α i and β i converge o zero following geomeric series as i increases. his indicaes ha he CCC vehicle relies more on signals obained from closer vehicles, and his characerisic behavior is no influenced by he choice of weighs q 1 and q 2 in he cos funcion 1,11). As an example, we consider q 1 = 2 [1/s 2 ], q 2 = 4, r = 1 [s 2 ] and assume non-ccc vehicles wih gains α =.6 [1/s] and β =.9 [1/s]. In his case 24) gives he gains α [1/s 2 ], β [1/s]. he exponenial decay of he gains α i, β i wih he vehicle index i = 2,..., n is demonsraed graphically in Fig. 3 for a 5 + 1) vehicle chain red circles) and for a 1 + 1) vehicle chain blue crosses). his is suppored by he fac ha he eigenvalues of M in 27) are λ 1 =.61, λ 2 =.37 and λ 3,4 =, which are locaed inside he uni circle. We remark ha he gains on signals for vehicles 7 1 are small, indicaing ha close-o-opimal design can be obained when only observing approximaely 5 6 vehicles ahead. In his sense, he benefis of increasing he number of cars ahead saurae, and having very long connecions may no be favorable as hey only make he nework srucure more complicaed. IV. HEAD-O-AIL SRING SABILIY In his secion, we discuss how he opimal CCC design influences he longiudinal sabiliy of he conneced vehicle sysem. In paricular, we analyze plan sabiliy and sring sabiliy. he plan sabiliy of a CCC vehicle is given as follows: suppose ha he vehicles whose signals are used by a CCC vehicle are driven a he same consan velociy, hen he velociy of he CCC vehicle approaches his consan velociy. Sring sabiliy characerizes he aenuaion of velociy flucuaions as hey propagae upsream [12]. However, sring sabiliy can only be ensured for he CCC vehicle as we have no conrol over he non-ccc vehicles). herefore, we evaluae he head-o-ail sring sabiliy, i.e., compare he velociy flucuaions of he head vehicle and he CCC vehicle a he ail. Noice ha his definiion allows ha vehicles in he middle may amplify he velociy flucuaions of vehicles ahead. Despie he presence of such inra-plaoon sring insabiliy, a CCC vehicle can be used o ensure heado-ail sring sabiliy. We consider he velociy perurbaion ṽ n+1 of he head vehicle as he inpu and he velociy perurbaion ṽ 1 of he ail vehicle as he oupu. Since perurbaion signals can be represened using Fourier componens and superposiion holds for linear sysems, he head-o-ail sring sabiliy is ensured when sinusoidal signals are aenuaed beween he head and he ail vehicles for all exciaion frequencies. hus, we consider he periodic exciaion ṽ n+1 ) = v amp sinω) and seady-sae soluion of 7,12,14,15), see Fig. 2b) for illusraion. hen aking he Laplace ransform of 14) wih zero iniial condiion and reversed ime while using P) P leads o W s) = si A + 1 r BB P) PDṼn+1s), 28) where W s) is he Laplace ransform of w), Ṽ n+1 s) is he Laplace ransform of ṽ n+1 ). aking he Laplace ransform of he sysem 7,12) wih zero iniial condiions and eliminaing he velociies of he oher vehicles and he headways, we obain he head-o-ail ransfer funcion Γs) = Ṽ1s) Ṽ n+1 s) = 1 α 1 Γ n 1 s) + sf n+1s) G 1 s) G 1 s) n + αi + β i s α i )Γ s) ) ) Γ n i s). 29) i=2 Here Ṽ1s) and Ṽn+1s) denoe he Laplace ransform of ṽ 1 ) and ṽ n+1 ), respecively, and Γ s) = F s) G s), G s) = s 2 + α + β)s + αf F s) = βs + αf, G 1 s) = s 2 β 1 s + α 1, F n+1 s) = α n + ββ n )s + Π 1n [1, 1] + βπ 1n [1, 2]. 3) Plan sabiliy a he linear level is deermined by he denominaor of he ransfer funcion in 29). he sysem 411

5 α =.6 [1/s] q 2 a) q 2 α =.9 [1/s] b) β =.6 [1/s] q 1 q 1 β =.9 [1/s] q 2 c) q 2 A B C d) q 1 q 1 Fig. 4. Sring sabiliy chars of a 5+1)-car plaoon in he q 1, q 2 )-plane for differen human parameers α, β as indicaed. he sring sable domains are shaded. is linearly plan sable, if and only if all soluions of he characerisic equaion G n s)g 1 s) = are locaed in he lef half complex plane. Noice ha plan sabiliy is only influenced by he human parameers α, β and he CCC gains α 1, β 1. Using Rouh-Hurwiz crieria, we obain he condiions for plan sabiliy Fig. 5. Magniude of ransfer funcion as a funcion of he exciaion frequency. Panels a c) correspond o poins marked A C in Fig. 4 c). β a) q 2 = 1 β b) q 2 = 5 β c) q 2 = 1 q1 = 1 α α α β d) q 1 = 2 β e) q 1 = 5 β f) q 1 = 1 q2 = 1 α >, α + β >, α 1 >, β 1 <. 31) In he following analysis, we only consider plan sable human parameers α, β. Soluion 23) provides he gains α 1, β 1 ha have o saisfy 31). hese can be used o obain all oher gains α i and β i, i = 2,..., n; see 18,27). A he linear level he necessary and sufficien condiion of head-o-ail sring sabiliy is given by Γiω) 2 1 = ω 2 fω) <, ω >, 32) where Γiω) is defined by 29,3); see [4], [13]. he order of fω) increases wih he number of vehicles n. Sring sabiliy is violaed when he maximum of fω) is larger han, and hus, he sring sabiliy boundary is given by he equaions fω cr ) =, fω cr ) ω =, 33) subjec o 2 fω cr ) ω 2 <, where ω cr indicaes he locaion of he maximum of fω). o obain sring sabiliy chars, we solve 33) numerically and plo he sring sabiliy boundary in he q 1, q 2 )-plane and in he α, β)-plane. Noe ha in pracical ranges of human parameers α and β he sring insabiliy only occurs a zero frequency i.e., ω cr = ). When he human parameers α and β change, he range of weighs q 1, q 2, r ha resuls in a sring sabiliy changes. Wihou loss of generaliy, we fix r = 1 [s 2 ] and only consider he change of weighs q 1, q 2. As observed in he previous secion, gains on vehicles i, i > 6, are small, and herefore we consider n = 5. α α α Fig. 6. Sring sabiliy chars of a 5 + 1)-car plaoon in he α, β)-plane. he noaion is he same as in Fig. 4. In Fig. 4 we fix he human parameers α =.6,.9 [1/s] and β =.6,.9 [1/s] and shade he sring sable domains in he q 1, q 2 )-plane. I has been shown [3] ha wihou CCC, he sysem is sring unsable when α + 2β 2f >, 34) and hus, we have a sring unsable sysem for q 1 = q 2 =. Fig. 4a) shows ha increasing he weighs q 1, q 2 on racking errors is beneficial for sring sabiliy, and ha no sring sabiliy exiss for q 1 1, q 2 4. However, q 1 and q 2 canno be chosen independenly. Similar resuls are observed in panels b d). Comparing he four panels, one can noice ha increasing eiher α or β increases he size of he sring sable domain, while he minimum q 2 ensuring sring sabiliy decreases o zero. However, for q 2 = he range of sring sable q 1 is very small. We remark ha sring sabiliy loss in he conneced vehicle sysem analyzed here only happens a zero frequency. Fig. 5 demonsraes such sabiliy loss for he poins A C marked in Fig. 4c). As q 2 decreases, he magniude of ransfer funcion 29) increases and i exceeds 1 in he lowfrequency domain in case C. he robusness of opimized CCC designs is invesigaed 4111

6 by varying he human parameers α and β and he resuls are summarized in Fig. 6. We fix q 1 = 1 [1/s 2 ] in panels a c) and q 2 = 1 in panels d f), and shade he sring sable domains. Increasing α and β improves sring sabiliy in each case. For fixed q 1, increasing q 2 enlarges he sring sabile area cf. panels a) and b)), bu oo large q 2 resuls in smaller sring sabiliy region cf. panels b) and c)). For fixed q 2, sring sabiliy increases wih larger q 1 ; see d f). We remark ha sring sabiliy may be los when increasing q 1 even furher, alhough q 1 < 1 is considered o be physically realisic. hus weighing heavily on eiher h 1 or v 1 is derimenal for sring sabiliy. Finally, o evaluae he performance of our CCC algorihm, we consider a 5+1)-car sysem wih sring unsable human parameers and invesigae he evoluion of headway and velociy errors by numerical simulaions. Fig. 7 shows he simulaion resuls of he 5+1)-car sysem wih gains generaed using differen design parameers q 1 and q 2. he simulaion resuls are presened for he parameers corresponding o poins A and C in Fig. 4c), while using he disurbance signal ṽ n+1 ) = v amp sinω) wih ampliude v amp = 1 [m/s] and frequency ω =.3 [rad/s]; see Fig. 5a,c) for he amplificaion plos. he simulaion resuls demonsrae ha case A is head-o-ail sring sable, as he CCC vehicle s velociy flucuaion hick blue curve) has smaller ampliude han he velociy inpu dashed curve) in Fig. 7a). On he oher hand, he CCC vehicle s velociy flucuaion in Fig. 7c) has larger ampliude han he velociy inpu, indicaing sring insabiliy. Noe ha in boh cases he ampliude of velociy flucuaions are amplified by non-ccc vehicles, because he human parameers α =.6 [1/s], β =.9 [1/s] are sring unsable; see 34). Sill, he CCC vehicle is able o mainain sring sabiliy when q 1, q 2 are chosen appropriaely. On he oher hand, he comparison of panels b) and d) shows a rade-off. While an increased weigh on velociy error q 2 ensures sring sabiliy, he relaive weigh on q 1 decreases, and he headway error increases hough sill aenuaed compared o head vehicle). V. CONCLUSION In his paper, we proposed a conneced cruise conrol design based on linear quadraic racking and analyzed he head-o-ail sring sabiliy. I was shown ha he gains depend on he human parameers and he design parameers in he cos funcion. We found ha he opimal gains on preceding vehicles are no influenced by dynamics of vehicles farher downsream, and ha he gains decrease wih he number of cars beween he CCC vehicle and he signaling vehicle. he opimized CCC is shown o be able o sabilize an oherwise sring unsable sysems when he weighs on he headway and velociy errors are chosen appropriaely. he design was robus agains variaions of human parameers, and he resuls were verified using numerical simulaions. Fuure research includes opimizing nonlinear CCC algorihms while considering more complicaed conneciviy srucures and imperfec communicaion v a) v c) 21 h b) h d) Fig. 7. Velociy and headway responses of a 5 + 1)-car vehicle sring wih human parameers α =.6 [1/s], β =.9 [1/s]. a,b): using design parameers q 1 = 2 [1/s 2 ], q 2 = 4. c,d): using design parameers q 1 = 2 [1/s 2 ], q 2 = 1. hick curves denoe he CCC vehicle, hin curves are used for non-ccc vehicles, while dashed curves indicae he velociy disurbance of he head vehicle. REFERENCES [1] G. Orosz, Conneced cruise conrol: modeling, delay effecs, and nonlinear behavior, Vehicle Sysem Dynamics, p. submied, 214. [2] G. Orosz, R. E. Wilson, and G. Sépán, raffic jams: dynamics and conrol, Philosophical ransacions of he Royal Sociey A, vol. 368, no. 1928, pp , 21. [3] L. Zhang and G. Orosz, Designing nework moifs in conneced vehicle sysems: delay effecs and sabiliy, in Proceedings of he ASME Dynamical Sysems and Conrol Conference. ASME, 213. [4] J. I. Ge and G. Orosz, Dynamics of conneced vehicle sysems wih delayed acceleraion feedback, ransporaion Research Par C: Emerging echnologies, vol. 46, pp , 214. [5] W. B. Qin, M. M. Gomez, and G. Orosz, Sabiliy analysis of auonomous cruise conrol wih sochasic delay, in Proceedings of American Conrol Conference, 214. [6] J. Ploeg, D. Shukla, N. van de Wouw, and H. Nijmeijer, Conroller synhesis for sring sabiliy of vehicle plaoons, IEEE ransacions on Inelligen ransporaion Sysems, vol. 15, no. 2, pp , 214. [7] M. Wang, W. Daamen, S. P. Hoogendoorn, and B. van Arem, Rolling horizon conrol framework for driver assisance sysems. par II: Cooperaive sensing and cooperaive conrol, ransporaion Research Par C: Emerging echnologies, vol. 4, pp , 214. [8] B. S. Kerner, he Physics of raffic. Springer, 24. [9] D. C. Gazis, R. Herman, and R. B. Pos, Car-following heory of seadysae raffic flow, Operaions Research, vol. 7, no. 4, pp , [1] A. D. May, raffic flow fundamenals. Prenice Hall, 199. [11] E. D. Sonag, Mahemaical Conrol heory: Deerminisic Finie Dimensional Sysems. Springer, [12] P. Seiler, A. Pan, and K. Hedrick, Disurbance propagaion in vehicle srings, Auomaic Conrol, IEEE ransacions on, vol. 49, no. 1, pp , 24. [13] L. Zhang and G. Orosz, Moif-based analysis of conneced vehicle sysems: delay effecs and sabiliy, Auomaica, p. submied,

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

10. State Space Methods

10. State Space Methods . Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he

More 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

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

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

Traveling Waves. Chapter Introduction

Traveling Waves. Chapter Introduction Chaper 4 Traveling Waves 4.1 Inroducion To dae, we have considered oscillaions, i.e., periodic, ofen harmonic, variaions of a physical characerisic of a sysem. The sysem a one ime is indisinguishable from

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

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

For example, the comb filter generated from. ( ) has a transfer function. e ) has L notches at ω = (2k+1)π/L and L peaks at ω = 2π k/l,

For example, the comb filter generated from. ( ) has a transfer function. e ) has L notches at ω = (2k+1)π/L and L peaks at ω = 2π k/l, Comb Filers The simple filers discussed so far are characeried eiher by a single passband and/or a single sopband There are applicaions where filers wih muliple passbands and sopbands are required The

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

A car following model for traffic flow simulation

A car following model for traffic flow simulation Inernaional Journal of Applied Mahemaical Sciences ISSN 0973-076 Volume 9, Number (206), pp. -9 Research India Publicaions hp://www.ripublicaion.com A car following model for raffic flow simulaion Doudou

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

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

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

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

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems. Mah 2250-004 Week 4 April 6-20 secions 7.-7.3 firs order sysems of linear differenial equaions; 7.4 mass-spring sysems. Mon Apr 6 7.-7.2 Sysems of differenial equaions (7.), and he vecor Calculus we need

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

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

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

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

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

4.5 Constant Acceleration

4.5 Constant Acceleration 4.5 Consan Acceleraion v() v() = v 0 + a a() a a() = a v 0 Area = a (a) (b) Figure 4.8 Consan acceleraion: (a) velociy, (b) acceleraion When he x -componen of he velociy is a linear funcion (Figure 4.8(a)),

More information

Mean-square Stability Control for Networked Systems with Stochastic Time Delay

Mean-square Stability Control for Networked Systems with Stochastic Time Delay JOURNAL OF SIMULAION VOL. 5 NO. May 7 Mean-square Sabiliy Conrol for Newored Sysems wih Sochasic ime Delay YAO Hejun YUAN Fushun School of Mahemaics and Saisics Anyang Normal Universiy Anyang Henan. 455

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

Estimation of feedback gains and delays in connected vehicle systems

Estimation of feedback gains and delays in connected vehicle systems 26 American Conrol Conference (ACC) Boson Marrio Copley Place July 6-8, 26 Boson, MA, USA Esimaion of feedback gains and delays in conneced vehicle sysems Jin I Ge and Gábor Orosz Absrac In his paper,

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

INDEX. Transient analysis 1 Initial Conditions 1

INDEX. Transient analysis 1 Initial Conditions 1 INDEX Secion Page Transien analysis 1 Iniial Condiions 1 Please inform me of your opinion of he relaive emphasis of he review maerial by simply making commens on his page and sending i o me a: Frank Mera

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

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

Matlab and Python programming: how to get started

Matlab and Python programming: how to get started Malab and Pyhon programming: how o ge sared Equipping readers he skills o wrie programs o explore complex sysems and discover ineresing paerns from big daa is one of he main goals of his book. In his chaper,

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

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

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

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

(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

Wall. x(t) f(t) x(t = 0) = x 0, t=0. which describes the motion of the mass in absence of any external forcing.

Wall. x(t) f(t) x(t = 0) = x 0, t=0. which describes the motion of the mass in absence of any external forcing. MECHANICS APPLICATIONS OF SECOND-ORDER ODES 7 Mechanics applicaions of second-order ODEs Second-order linear ODEs wih consan coefficiens arise in many physical applicaions. One physical sysems whose behaviour

More information

Robotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1.

Robotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1. Roboics I April 11, 017 Exercise 1 he kinemaics of a 3R spaial robo is specified by he Denavi-Harenberg parameers in ab 1 i α i d i a i θ i 1 π/ L 1 0 1 0 0 L 3 0 0 L 3 3 able 1: able of DH parameers of

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

Continuous Time. Time-Domain System Analysis. Impulse Response. Impulse Response. Impulse Response. Impulse Response. ( t) + b 0.

Continuous Time. Time-Domain System Analysis. Impulse Response. Impulse Response. Impulse Response. Impulse Response. ( t) + b 0. Time-Domain Sysem Analysis Coninuous Time. J. Robers - All Righs Reserved. Edied by Dr. Rober Akl 1. J. Robers - All Righs Reserved. Edied by Dr. Rober Akl 2 Le a sysem be described by a 2 y ( ) + a 1

More information

Problem Set #1. i z. the complex propagation constant. For the characteristic impedance:

Problem Set #1. i z. the complex propagation constant. For the characteristic impedance: Problem Se # Problem : a) Using phasor noaion, calculae he volage and curren waves on a ransmission line by solving he wave equaion Assume ha R, L,, G are all non-zero and independen of frequency From

More information

Kinematics Vocabulary. Kinematics and One Dimensional Motion. Position. Coordinate System in One Dimension. Kinema means movement 8.

Kinematics Vocabulary. Kinematics and One Dimensional Motion. Position. Coordinate System in One Dimension. Kinema means movement 8. Kinemaics Vocabulary Kinemaics and One Dimensional Moion 8.1 WD1 Kinema means movemen Mahemaical descripion of moion Posiion Time Inerval Displacemen Velociy; absolue value: speed Acceleraion Averages

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

Module 4: Time Response of discrete time systems Lecture Note 2

Module 4: Time Response of discrete time systems Lecture Note 2 Module 4: Time Response of discree ime sysems Lecure Noe 2 1 Prooype second order sysem The sudy of a second order sysem is imporan because many higher order sysem can be approimaed by a second order model

More 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

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

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

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

Ground Rules. PC1221 Fundamentals of Physics I. Kinematics. Position. Lectures 3 and 4 Motion in One Dimension. A/Prof Tay Seng Chuan

Ground Rules. PC1221 Fundamentals of Physics I. Kinematics. Position. Lectures 3 and 4 Motion in One Dimension. A/Prof Tay Seng Chuan Ground Rules PC11 Fundamenals of Physics I Lecures 3 and 4 Moion in One Dimension A/Prof Tay Seng Chuan 1 Swich off your handphone and pager Swich off your lapop compuer and keep i No alking while lecure

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

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

Chapter 1 Fundamental Concepts

Chapter 1 Fundamental Concepts Chaper 1 Fundamenal Conceps 1 Signals A signal is a paern of variaion of a physical quaniy, ofen as a funcion of ime (bu also space, disance, posiion, ec). These quaniies are usually he independen variables

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

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

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

!!"#"$%&#'()!"#&'(*%)+,&',-)./0)1-*23)

!!#$%&#'()!#&'(*%)+,&',-)./0)1-*23) "#"$%&#'()"#&'(*%)+,&',-)./)1-*) #$%&'()*+,&',-.%,/)*+,-&1*#$)()5*6$+$%*,7&*-'-&1*(,-&*6&,7.$%$+*&%'(*8$&',-,%'-&1*(,-&*6&,79*(&,%: ;..,*&1$&$.$%&'()*1$$.,'&',-9*(&,%)?%*,('&5

More information

Optimal Path Planning for Flexible Redundant Robot Manipulators

Optimal Path Planning for Flexible Redundant Robot Manipulators 25 WSEAS In. Conf. on DYNAMICAL SYSEMS and CONROL, Venice, Ialy, November 2-4, 25 (pp363-368) Opimal Pah Planning for Flexible Redundan Robo Manipulaors H. HOMAEI, M. KESHMIRI Deparmen of Mechanical Engineering

More information

6.003 Homework #9 Solutions

6.003 Homework #9 Solutions 6.00 Homework #9 Soluions Problems. Fourier varieies a. Deermine he Fourier series coefficiens of he following signal, which is periodic in 0. x () 0 0 a 0 5 a k sin πk 5 sin πk 5 πk for k 0 a k 0 πk j

More information

Chapter 3 Boundary Value Problem

Chapter 3 Boundary Value Problem Chaper 3 Boundary Value Problem A boundary value problem (BVP) is a problem, ypically an ODE or a PDE, which has values assigned on he physical boundary of he domain in which he problem is specified. Le

More information

EE363 homework 1 solutions

EE363 homework 1 solutions EE363 Prof. S. Boyd EE363 homework 1 soluions 1. LQR for a riple accumulaor. We consider he sysem x +1 = Ax + Bu, y = Cx, wih 1 1 A = 1 1, B =, C = [ 1 ]. 1 1 This sysem has ransfer funcion H(z) = (z 1)

More information

h[n] is the impulse response of the discrete-time system:

h[n] is the impulse response of the discrete-time system: Definiion Examples Properies Memory Inveribiliy Causaliy Sabiliy Time Invariance Lineariy Sysems Fundamenals Overview Definiion of a Sysem x() h() y() x[n] h[n] Sysem: a process in which inpu signals are

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

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

Section 7.4 Modeling Changing Amplitude and Midline

Section 7.4 Modeling Changing Amplitude and Midline 488 Chaper 7 Secion 7.4 Modeling Changing Ampliude and Midline While sinusoidal funcions can model a variey of behaviors, i is ofen necessary o combine sinusoidal funcions wih linear and exponenial curves

More information

Math Week 15: Section 7.4, mass-spring systems. These are notes for Monday. There will also be course review notes for Tuesday, posted later.

Math Week 15: Section 7.4, mass-spring systems. These are notes for Monday. There will also be course review notes for Tuesday, posted later. Mah 50-004 Week 5: Secion 7.4, mass-spring sysems. These are noes for Monday. There will also be course review noes for Tuesday, posed laer. Mon Apr 3 7.4 mass-spring sysems. Announcemens: Warm up exercise:

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

Laplace Transform and its Relation to Fourier Transform

Laplace Transform and its Relation to Fourier Transform Chaper 6 Laplace Transform and is Relaion o Fourier Transform (A Brief Summary) Gis of he Maer 2 Domains of Represenaion Represenaion of signals and sysems Time Domain Coninuous Discree Time Time () [n]

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

ON QUANTIZATION AND COMMUNICATION TOPOLOGIES IN MULTI-VEHICLE RENDEZVOUS 1. Karl Henrik Johansson Alberto Speranzon,2 Sandro Zampieri

ON QUANTIZATION AND COMMUNICATION TOPOLOGIES IN MULTI-VEHICLE RENDEZVOUS 1. Karl Henrik Johansson Alberto Speranzon,2 Sandro Zampieri ON QUANTIZATION AND COMMUNICATION TOPOLOGIES IN MULTI-VEHICLE RENDEZVOUS 1 Karl Henrik Johansson Albero Speranzon, Sandro Zampieri Deparmen of Signals, Sensors and Sysems Royal Insiue of Technology Osquldas

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

On-line Adaptive Optimal Timing Control of Switched Systems

On-line Adaptive Optimal Timing Control of Switched Systems On-line Adapive Opimal Timing Conrol of Swiched Sysems X.C. Ding, Y. Wardi and M. Egersed Absrac In his paper we consider he problem of opimizing over he swiching imes for a muli-modal dynamic sysem when

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

Let us start with a two dimensional case. We consider a vector ( x,

Let us start with a two dimensional case. We consider a vector ( x, Roaion marices We consider now roaion marices in wo and hree dimensions. We sar wih wo dimensions since wo dimensions are easier han hree o undersand, and one dimension is a lile oo simple. However, our

More information

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n Module Fick s laws of diffusion Fick s laws of diffusion and hin film soluion Adolf Fick (1855) proposed: d J α d d d J (mole/m s) flu (m /s) diffusion coefficien and (mole/m 3 ) concenraion of ions, aoms

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

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

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

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

Multi-scale 2D acoustic full waveform inversion with high frequency impulsive source

Multi-scale 2D acoustic full waveform inversion with high frequency impulsive source Muli-scale D acousic full waveform inversion wih high frequency impulsive source Vladimir N Zubov*, Universiy of Calgary, Calgary AB vzubov@ucalgaryca and Michael P Lamoureux, Universiy of Calgary, Calgary

More information

Time Domain Transfer Function of the Induction Motor

Time Domain Transfer Function of the Induction Motor Sudies in Engineering and Technology Vol., No. ; Augus 0 ISSN 008 EISSN 006 Published by Redfame Publishing URL: hp://se.redfame.com Time Domain Transfer Funcion of he Inducion Moor N N arsoum Correspondence:

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

5. Stochastic processes (1)

5. Stochastic processes (1) Lec05.pp S-38.45 - Inroducion o Teleraffic Theory Spring 2005 Conens Basic conceps Poisson process 2 Sochasic processes () Consider some quaniy in a eleraffic (or any) sysem I ypically evolves in ime randomly

More information

EECE 301 Signals & Systems Prof. Mark Fowler

EECE 301 Signals & Systems Prof. Mark Fowler EECE 3 Signals & Sysems Prof. Mark Fowler Noe Se #2 Wha are Coninuous-Time Signals??? Reading Assignmen: Secion. of Kamen and Heck /22 Course Flow Diagram The arrows here show concepual flow beween ideas.

More information

An Introduction to Malliavin calculus and its applications

An Introduction to Malliavin calculus and its applications An Inroducion o Malliavin calculus and is applicaions Lecure 5: Smoohness of he densiy and Hörmander s heorem David Nualar Deparmen of Mahemaics Kansas Universiy Universiy of Wyoming Summer School 214

More information

Pade and Laguerre Approximations Applied. to the Active Queue Management Model. of Internet Protocol

Pade and Laguerre Approximations Applied. to the Active Queue Management Model. of Internet Protocol Applied Mahemaical Sciences, Vol. 7, 013, no. 16, 663-673 HIKARI Ld, www.m-hikari.com hp://dx.doi.org/10.1988/ams.013.39499 Pade and Laguerre Approximaions Applied o he Acive Queue Managemen Model of Inerne

More information

System of Linear Differential Equations

System of Linear Differential Equations Sysem of Linear Differenial Equaions In "Ordinary Differenial Equaions" we've learned how o solve a differenial equaion for a variable, such as: y'k5$e K2$x =0 solve DE yx = K 5 2 ek2 x C_C1 2$y''C7$y

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

Chapter #1 EEE8013 EEE3001. Linear Controller Design and State Space Analysis

Chapter #1 EEE8013 EEE3001. Linear Controller Design and State Space Analysis Chaper EEE83 EEE3 Chaper # EEE83 EEE3 Linear Conroller Design and Sae Space Analysis Ordinary Differenial Equaions.... Inroducion.... Firs Order ODEs... 3. Second Order ODEs... 7 3. General Maerial...

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

A First Course on Kinetics and Reaction Engineering. Class 19 on Unit 18

A First Course on Kinetics and Reaction Engineering. Class 19 on Unit 18 A Firs ourse on Kineics and Reacion Engineering lass 19 on Uni 18 Par I - hemical Reacions Par II - hemical Reacion Kineics Where We re Going Par III - hemical Reacion Engineering A. Ideal Reacors B. Perfecly

More information

Macroeconomic Theory Ph.D. Qualifying Examination Fall 2005 ANSWER EACH PART IN A SEPARATE BLUE BOOK. PART ONE: ANSWER IN BOOK 1 WEIGHT 1/3

Macroeconomic Theory Ph.D. Qualifying Examination Fall 2005 ANSWER EACH PART IN A SEPARATE BLUE BOOK. PART ONE: ANSWER IN BOOK 1 WEIGHT 1/3 Macroeconomic Theory Ph.D. Qualifying Examinaion Fall 2005 Comprehensive Examinaion UCLA Dep. of Economics You have 4 hours o complee he exam. There are hree pars o he exam. Answer all pars. Each par has

More information

KEY. Math 334 Midterm III Winter 2008 section 002 Instructor: Scott Glasgow

KEY. Math 334 Midterm III Winter 2008 section 002 Instructor: Scott Glasgow KEY Mah 334 Miderm III Winer 008 secion 00 Insrucor: Sco Glasgow Please do NOT wrie on his exam. No credi will be given for such work. Raher wrie in a blue book, or on your own paper, preferably engineering

More information

6.003 Homework #9 Solutions

6.003 Homework #9 Solutions 6.003 Homework #9 Soluions Problems. Fourier varieies a. Deermine he Fourier series coefficiens of he following signal, which is periodic in 0. x () 0 3 0 a 0 5 a k a k 0 πk j3 e 0 e j πk 0 jπk πk e 0

More information

Electrical Circuits. 1. Circuit Laws. Tools Used in Lab 13 Series Circuits Damped Vibrations: Energy Van der Pol Circuit

Electrical Circuits. 1. Circuit Laws. Tools Used in Lab 13 Series Circuits Damped Vibrations: Energy Van der Pol Circuit V() R L C 513 Elecrical Circuis Tools Used in Lab 13 Series Circuis Damped Vibraions: Energy Van der Pol Circui A series circui wih an inducor, resisor, and capacior can be represened by Lq + Rq + 1, a

More information

Theory of! Partial Differential Equations!

Theory of! Partial Differential Equations! hp://www.nd.edu/~gryggva/cfd-course/! Ouline! Theory o! Parial Dierenial Equaions! Gréar Tryggvason! Spring 011! Basic Properies o PDE!! Quasi-linear Firs Order Equaions! - Characerisics! - Linear and

More information

After the completion of this section the student. Theory of Linear Systems of ODEs. Autonomous Systems. Review Questions and Exercises

After the completion of this section the student. Theory of Linear Systems of ODEs. Autonomous Systems. Review Questions and Exercises Chaper V ODE V.5 Sysems of Ordinary Differenial Equaions 45 V.5 SYSTEMS OF FIRST ORDER LINEAR ODEs Objecives: Afer he compleion of his secion he suden - should recall he definiion of a sysem of linear

More information

Chapter 7 Response of First-order RL and RC Circuits

Chapter 7 Response of First-order RL and RC Circuits Chaper 7 Response of Firs-order RL and RC Circuis 7.- The Naural Response of RL and RC Circuis 7.3 The Sep Response of RL and RC Circuis 7.4 A General Soluion for Sep and Naural Responses 7.5 Sequenial

More information

ψ(t) = V x (0)V x (t)

ψ(t) = V x (0)V x (t) .93 Home Work Se No. (Professor Sow-Hsin Chen Spring Term 5. Due March 7, 5. This problem concerns calculaions of analyical expressions for he self-inermediae scaering funcion (ISF of he es paricle in

More information

Bifurcation Analysis of a Stage-Structured Prey-Predator System with Discrete and Continuous Delays

Bifurcation Analysis of a Stage-Structured Prey-Predator System with Discrete and Continuous Delays Applied Mahemaics 4 59-64 hp://dx.doi.org/.46/am..4744 Published Online July (hp://www.scirp.org/ournal/am) Bifurcaion Analysis of a Sage-Srucured Prey-Predaor Sysem wih Discree and Coninuous Delays Shunyi

More information

Online Appendix to Solution Methods for Models with Rare Disasters

Online Appendix to Solution Methods for Models with Rare Disasters Online Appendix o Soluion Mehods for Models wih Rare Disasers Jesús Fernández-Villaverde and Oren Levinal In his Online Appendix, we presen he Euler condiions of he model, we develop he pricing Calvo block,

More information