Optimal control of connected vehicle systems
|
|
- Kathlyn Short
- 5 years ago
- Views:
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
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 information10. 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 informationSTATE-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
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 informationLecture 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 informationTraveling 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 informationStability 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 informationSome 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 informationFor 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 informationSliding 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 informationA 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 informationClass 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 informationPhysics 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 information4.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 information4. 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 informationMath 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 informationChapter 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 information15. 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 informationt 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 information1. 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 informationVehicle 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 information4.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 informationMean-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 informationTwo 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 informationEstimation 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 informationdi 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 informationINDEX. 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 informationSignal 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 informationWEEK-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 informationMatlab 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 informationLAPLACE 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 informationSection 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 informationGuest 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 informationSliding 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)
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 informationWall. 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 informationRobotics 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 informationDifferential 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 informationContinuous 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 informationProblem 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 informationKinematics 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 informationLecture 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 informationModule 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 information3.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 informationChapter 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 information2. 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 informationEXERCISES 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 informationGround 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 informationODEs 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 informationSOLUTIONS 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 informationChapter 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 informationLinear 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 informationVanishing 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 informationInventory 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)
"#"$%&#'()"#&'(*%)+,&',-)./)1-*) #$%&'()*+,&',-.%,/)*+,-&1*#$)()5*6$+$%*,7&*-'-&1*(,-&*6&,7.$%$+*&%'(*8$&',-,%'-&1*(,-&*6&,79*(&,%: ;..,*&1$&$.$%&'()*1$$.,'&',-9*(&,%)?%*,('&5
More informationOptimal 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 information6.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 informationChapter 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 informationEE363 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 informationh[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 informationIB 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 informationThe 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 informationSection 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 informationMath 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 informationA 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 informationLaplace 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 informationdt = 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 informationON 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 informationNotes 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 informationOn-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 informationNumerical 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 informationLet 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 informationModule 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 informationKINEMATICS 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 informationSimulation-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 information23.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 information2.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 informationMulti-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 informationTime 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 informationIn 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 information5. 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 informationEECE 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 informationAn 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 informationPade 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 informationSystem 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 informationCHAPTER 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 informationChapter #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 information8. 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 informationA 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 informationMacroeconomic 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 informationKEY. 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 information6.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 informationElectrical 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 informationTheory 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 informationAfter 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 informationChapter 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)
.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 informationBifurcation 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 informationOnline 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