TRACKING CONTROL OF MULTIPLE MOBILE ROBOTS: A CASE STUDY OF INTER-ROBOT COLLISION-FREE PROBLEM

Size: px
Start display at page:

Download "TRACKING CONTROL OF MULTIPLE MOBILE ROBOTS: A CASE STUDY OF INTER-ROBOT COLLISION-FREE PROBLEM"

Transcription

1 265 Asian Journal of Control, Vol. 4, No. 3, pp , September 22 TRACKING CONTROL OF MULTIPLE MOBILE ROBOTS: A CASE STUDY OF INTER-ROBOT COLLISION-FREE PROBLEM Jurachart Jongusuk an Tsutomu Mita ABSTRACT Tracking control problem of multiple mobile robots is consiere. Our system is compose of a reference an two follower robots of unicycle type, which have their own priority numbers. The purpose is to control the two followers so that the reference is tracke with arbitrary esire clearance an also to avoi an inter-robot collision, which can occur ranomly uring the control process subjecte to communication range limitation. We introuce two switching controllers an appropriately using their convergence properties leaing to the collision-free movement. Simulation results prove efficiency of our control techniques. KeyWors: Multiple mobile robots, feeback, collision, tracking control, terminal attractor. I. INTRODUCTION Control of multiple mobile robots has been subject of a consierable research effort over the last few years. The rapi progress in this fiel has been the interplay of systems, theories an problems [1]. Three examples are introuce as the recent research trens, namely foraging, box-pushing an traffic control, foraging aresses the problem of learning an istribute algorithm; box-pushing focuses particularly on cooperative manipulation; an traffic control aresses the problem occurring when multiple robots move in a common environment that typically attempt to avoi collisions. Despite the vast amount of papers publishe on the collision avoiance motivate by the case of traffic control, the majority has concentrate on centralize path an motion plannings, which mostly utilize fee-forwar control in orer to accomplish a goal, leaing also to a computational time consumption, see for examples Arkin [2] an Lavalle [3], while less attention has been pai to control of the multiple mobile robots by a concept of feeback control. Desai [5] has stuie the problem of control an formation where a notation of l-l control is originally presente, however, the issue of sensor capability is absent. l-l control gives a motivation of using relative ynamics in the system, which will be further extene in our paper. In this paper, we focus on tasks in which the multiple mobile robots are require to follow a given trajectory in Manuscript receive January 3, 21; accepte July 18, 21. The authors are with Department of Control an Systems Engineering, Tokyo Institute of Technology, Japan. leaer-follower formation while avoiing collision among themselves. In aition, because of the high-banwith of inter-robot communication, consieration of communication capability becomes necessary. Here, a limite range of communication is consiere. The problem of control an coorination for multiple mobile robots is formulate using ecentralize controllers. We assume that there are no obstacle in the work space, which simplifies the problem to avoiing inter-robot collision. We also propose a simple technique that helps in ecision process of collision avoiance base on geometric analysis. This proposition correspons to the issue of Latombe [4] state that collision avoiance is an intrinsically geometric problem in configuration space. The paper is organize as follows: Section 2 escribes moel of the wheele mobile robot of a unicycle type an formulates the tracking problem for a group of three mobile robots without the presence of obstacle. Section 3 introuces a simple collision etection moel, virtual robot tracking control (VR control) as a main controller an collision avoiance control (l-l control), which is the extene version of the original work of Desai [5]. Analysis of each controllers properties an combination is also provie, which leas to the main issue. Section 4 contains simulation results that illustrate the performance of the propose control strategies. The paper conclues with the summary of results an recommenations for further research. II. PROBLEM FORMULATION A unicycle-type mobile robot, of which kinematic moel is efine by

2 Asian Journal of Control, Vol. 4, No. 3, September q i = B i u i = cos θ i sin θ i 1 u i (1) where q i =[x i, y i, θ i ] T an u i =[υ i, ω i ] T is consiere. This moel captures the characteristic of a class of restricte mobile robot which is what we nee to investigate. The robot are assume to satisfy pure rolling an non-slipping conitions, which lea to nonholonomic velocity constraints, x i sin θ i y i cos θ i =, (2) x i cos θ i + y i sin θ i = υ i (3) III. CONTROLLER DESIGN AND ANALYSIS We arrive here with assumptions efine in the last section. We will first apply a VR tracking control approach 1 to the system to complete the requirement in item (1) of the problem while aing an ability of collision checking. The check is one in a sufficiently small sampling time. When a robot etects possibility to collie with another, it checks if its priority number is higher (small value) or not. If not, it has to change the control to l-l control 2 in orer to avoi the collision an complete the requirement of item (2) of the problem. It is easy to unerstan the algorithm above by the flowchart shown in Fig Assumptions We summarize here all the assumptions use in this paper. (i) System is homogeneous, i.e., all robots are of the same moel escribe in (1) an satisfy the velocity constraints, (2) an (3). (ii) The workspace (xy-plane) is flat an contains no obstacle. This comforts us in fining solution focusing on only inter-robot collision problem. (iii) The reference robot follows a smooth trajectory, an maintains positive velocity. We nee this assumption in orer to guarantee that our problem is tracking not stabilizing control, see [7]. (iv) Each follower robot is inexe by ifferent priority number. Then, system becomes hierachy an we can manage the behavior of each robot uring switching control as will be seen in algorithm in Section 3. (v) Each robot can extract any necessary information from its communication equipment with negligible elay. In practical case, high-spee ceiling camera shoul be use to track the absolute position of each robot. 2.2 Problem statement To keep our investigation on the way, the problem is preliminarily state as follows. Statement 1. Giving initial positions an orientations, q i () for the follower robot i, an the motion of reference robot, it is require that control law u i gives a motion satisfying assumptions (i)-(v) such that, as t, 1 formation is establishe, 2. no collision occurs among robot i an any robot j, an 4. a whole motion satisfies the limitation of communication range. 1 See etails in Section See etails in Section 3.3. Fig. 1. Algorithm flowchart. Referring to this algorithm, we nee to esign a simple collision etecting technique an introuce two feeback control approaches, VR tracking control an l-l control, that also fulfill the requirement in item (3) of the problem, which are the keys of the main issue. 3.1 Detecting collision A simple technique use to etect the collision between any two robots can be realize by moeling both

3 267 J. Jongusuk an T. Mita: Tracking Control of Multiple Mobile Robots: A Case Stuy of Inter-Robot robots to be covere with ouble circles centering at the control points, referring to Fig. 2. The soli ones cover the entire robots with raius D while the outer otte circles with raius D + is esigne for the require clearances between robots. For convenience, we efine the istance between two robots in the following efinition. Definition 1. Let (x i, y i ) an (x j, y j ) enote the control points of robot i an robot j, respectively. Distance between robots is efine by ρ ij = (x i x j ) 2 +(y i y j ) 2 (4) Then, we efine a function that can be use to check the collision between robots i an j as follows. f ij = ρ ij2 4(D + ) 2 f ij > safety collision (5) Here, f ij mathematically represents a conition for item (2) of the problem statement. Note that the choice of shoul be esigne with appropriate value. A large increases safety level but ecreases the smoothness of the robot because it etects the possibility of collision more often, which yiels more switching of control inputs. 3.2 Virtual robot tracking control This approach offers a tracking control with arbitrary clearance between robots, which is, in fact, necessary for navigating the multiple mobile robots while maintaining the formation as we nee. However, the collision that can occur uring transient state of the control is not consiere. We give the etail of this approach in this section. To easily separate many robot configuration symbols, we will use the following efinitions thereafter. q r = [x r, y r, θ r ] T refers to a reference robot, q i = [x i, y i, θ i ] T, i = 1, 2,, refers to follower robot i an q υi =[x υi, y υi, θ υi ] T is for VR of follower robot i. In orer to separate follower from reference robot avoiing collision, a concept of VR is efine with aitional objective to force error parameters to zero when system approaches final time, as follows. Definition 2. Virtual robot (VR): a hypothetical robot whose orientation is ientical to that of its corresponing follower robot, but position is place apart from it by the preefine r l clearances. The symbols l, r enote longituinal clearance an clearance along rear wheel axis, respectively. The relation between VR an the follower robot is written as follows. x υi = x i r sin θ i + l cos θ i y υi = y i + r cos θ i + l sin θ i θ υi = θ i (6) Note that the above equation is satisfie when the follower is on VR s right-han sie. From (6) an assumptions (2), (3), we can erive the kinematic moel of VR as follows. q υi = cos θ i r cos θ i l sin θ i sin θ i r sin θ i + l cos θ i u i = 1 B vi 1 u i (7) See proof in appenix A. Figure 3 escribes the parameters of VR on xy-plane. A stanar technique of I/O linearization is use to generate a control law for each follower [9]. A property inherite in this controller is escribe by the following statement. Property 1. VR tracking control gives solution that exponentially converges in an internal shape x e an y e. D Fig. 2. Collision avoiance moel. Fig. 3. VR tracking moel.

4 Asian Journal of Control, Vol. 4, No. 3, September Proof. An error moel is formulate as z ei = z υi z r = B vi u i b r υ r (8) where B vi is the matrix efine in (7), υ r enotes reference robot s translational velocity an u i enotes follower robot s velocity vector, an z υi = x υi y υi, z r = x r y r, b r = cos θ r sin θ r After applying I/O linearization, we obtain u i = B vi 1 (b r υ r λz ei ) (9) where λ enotes positive-constant iagonal matrix. The solution of this controlle system, is z ei = e λt z ei () (1) which itself tells that the solution exponentially converges in an internal shape x e an y e. Moreover, solution is wellefine because parameter l is esigne not to be zero, an hence, B vi 1 is non-singular. Curve in Fig. 3 illustrates this property. In aition, zero ynamics given in term of θ i must be prove to be stable in orer to guarantee the use of control law in (9). The ifferential equation for θ i is given by θ i =( λ 1x ei υ r cos θ r l or can be written in the form ) sin θ i +( λ 2y ei υ r sin θ r ) cos θ l i (11) θ i = k 1 sin (θ i + k 2 ) (12) Corollary 1. Property 1 ensures that applying (9) to follower robot will automatically avoi collision with reference robot if the conition below is satisfie, ρ ri (t ) > 2 ( + D) (15) where ρ ri enotes istance between reference an follower robot i an t is initial time. Proof. Recall that (1) ensures exponential convergence of ρ ri. Also, the preetermine r l clearance implies that ρ ri r 2 + l 2. Obviously, if robots o not with collision, i.e. (15) is satisfie, the motions of follower robots are free of collision with reference robot. Remark. Corollary 1 only guarantees collision-free motion between reference an one of follower robot. Hence, there is no guarantee of the motion between follower robots themselves. 3.3 l-l control The aim of this control is to maintain the esire lengths, l 13 an l 23 of consiere robot (Robot 3 in Fig. 4) from its two leaers. Referring to the main algorithm, at any collision occurrence, the lowest priority robot will be switche from VR tracking control to this control. l-l control is originally presente in Desai [5]; however, the convergence of states is not establishe within finite time an is aime for maintaining formation. Here, we utilize l-l control as collision avoiance controller, which requires that formation must be establishe within finite time. The kinematic equations for the Robot 3 is given as follows, where k 1 = 1 l (λ 1 x e υ r cos θ r ) 2 +( λ 2 y e + υ r sin θ r ) 2 k 2 = arctan λ 2 y e + υ r sin θ r λ 1 x e υ r cos θ r The trajectory, θ e w i = k 2 arcsin r is shown to be k 1 exponentially stable by linearization [1] of (12) as follows. θ i = k 1 cos (θ i + k 2 ) θ i = θ i e θ i (13) = k 1 2 w r 2 θ i (14) 1 a 2 an enotes a evia- Note that cos(arcsin (a)) = tion operator. Fig. 4. Notation for l-l control moel.

5 269 J. Jongusuk an T. Mita: Tracking Control of Multiple Mobile Robots: A Case Stuy of Inter-Robot an i 13 i 23 = cos γ 1 D sin γ 1 cos γ 2 D sin γ 2 υ 3 ω 3 υ 1 cos ψ 13 υ 2 cos ψ 23 = B ll u 3 υ ll (16) Again, zero ynamics in term of θ 3 can be prove stable, clearly when two leaers follows parallel straight lines with the same velocities or υ 1 = υ 2 = υ >, ω 1 = ω 2 =, θ 1 () = θ 2 () = θ. From (16), (19), we have θ 3 = υ cos ψ 13 cos (ψ 23 + θ θ 3 ) D sin (ψ 13 ψ 23 ) θ 3 = ω 3 where γ i = θ i + ψ i3 θ 3, (i = 1, 2). Remark. In (16), B ll is singular only if γ 2 γ 1 = nπ, (n =, 1,, 1, ) which etermines the three robots lie on the same line connecting them. We will refer it to singular line later in the context. For the purpose of collision avoiance, a property below is realize by I/O linearization. Property 2. l-l control gives solution that monotonically converges in an internal shape l 13 an l 23 within finite time. Proof. Require that the controlle variables l 13 an l 23 monotonically converge to l 13 an l 23, respectively, within finite time T r, we set i 13 = α 1 (l 13 i 23 α 2 (l 23 l 13 ) p q l 23 ) p q (17) = α l e (18) which is ientical to a terminal attractor moel where (p, q) can be, for example, set to (2, 3). See Appenix B for etails. By arbitrary setting reaching time T r, we obtain a formula for fining α 1, α 2 as follows. α 1 = sign (l 13 l 13 ()) l 1 13 l 13 () 3 /( T r 3 ) α 2 = sign (l 23 l 23 ()) l 1 23 l 23 () 3 /( T r 3 ) The sign( ) an are use to avoi complex solutions. We can generate a feeback control law as follows. u 3 = B ll 1 (υ ll + α l e ), t T r (19) The solution of this controlle system is l e = (l i3 l i3 ()) 1 3 α 3 i 3 t, i =1,2, (2) which proves the properties. Note that the term in parenthesis is linear. υ cos ψ 23cos (ψ 13 + θ θ 3 ) D sin (ψ 13 ψ 23 ) = υ D sin (θ θ 3 ) (21) Lyapunov caniate, V = 1 cos(θ θ 3 ) guarantees the stability of (21). In case of more complicate motions of two leaers, the stability is guarantee by observing that ρ 12 exponentially ecreases accoring to Property 1. This means l 13, l 23 are always etermine or the only singularity γ 1 γ 2 = is not passe. Remark. Note that (2) also gives an interesting convergent motion in an internal shape l 13 an l 23. It matches with our require motion when collision is etecte, i.e., when Robot 3 gets close to the two leaers, with the control (19), it tens to monotonically separate from the leaers avoiing collision. 3.4 Combination of VR an l-l control We can formulate a control of the whole system using the combination of three techniques as state in Section The flow of control process is formerly presente in this section; however, conitions in which the control process can be performe properly are not analyze. We summarize all limitations occurs from the combination as follows. 1) r 2 + l 2 > 4(D + ) 2 (22) 2) l i3 >2(D + ), i =1,2 (23) 3) ρ 12 < l 13 + l 23 (24) Conition (22) is extracte from the first requirement in item (1) of the problem statement an the esign of collision etection in Section 3.1. Conition (23) comes from the restriction of the robots not to get closer than safety regions when l-l control is performe. Again, in l- l control phase, l 13 an l 23 must be esigne in such a way that three robots can form a triangle formation at the en of control, which yiels (24). As state in Section 3.3, when three robots lie on the same line connecting them, the control law u 3 in (19) is not etermine. In aition, the requirement in item (3) of the problem states that communication range is limite, namely,

6 Asian Journal of Control, Vol. 4, No. 3, September R as maximum range. These two facts enable us to efine the accessible area, shae ones in Fig. 5 an Fig. 6, where the low priority robot can move on uring collision avoiance. Accessible area comes from the intersection between the circle R an the area evie by singular line in which the lowest priority robot exists. Hence, in l-l control phase, the esign of l 13 an l 23 (equivalent to the esign of P x ) is simply one by setting the values within accessible area. Physically, P x is the point where follower is suppose to be move to at the en of l-l control. In fact, the esign of P x can be classifie into two cases: esire target is outsie accessible area, Fig. 5, an insie the area, Fig. 6. The two cases can be istinguishe by the conition below. (For the sake of simplicity, we assume that reference robot moves along the positive y-irection.) if θ a < θ b Target is outsie accessible area else Target is insie accessible area where θ a, θ b [, 2π) are measure in clockwise fashion. θ a is an angle measure from follower 2 to follower 1 centere at reference robot while θ b is measure from follower 2 to its formation target TG 2. Remark. In case like Fig. 5, it is strongly recommene to efine P x behin follower 1 to ensure that collision will harly or not repeat after switching to VR tracking control. 3.5 Main issue Finally, we summarize all results in this paper as a main issue in theorem below. Theorem 1. A system of three mobile robots whose moels are efine by (1), operating in an environment escribe by assumptions (i)-(v) an controlle by VR tracking control law: u i in (9) as a formation controller, an l-l control law: u 3 in (19) as a collision controller, gives a motion of the form, Reference Follower 2 ρ ri R, i =1,2 (25) while avoiing inter-robot collision an maintaining a preetermine formation in the steay state. R is positive scalar corresponing to maximum communication length of reference robot. TG 2 θ b TG 1 θ a P x Proof. It can be simply checke that the istribution of the three robots is guarantee to be insie a circle of raius R centere at reference robot by Property 1, Property 2 an esign process in Section 3.4 along the motion if the two follower robots are initially in this circle. Follower 1 R Fig. 5. Target TG 2 is outsie accessible area : θ a < θ b. Reference Follower 1 IV. SIMULATION RESULTS The aim of these simulations is to valiate our tracking control metho for a multiple mobile robot (threerobot case). Our goal is to make the three robots form a leaer-wing motion. Two cases when target is initially locate outsie an insie accessible area are shown. We from the former case where initial parameters are set as follows. θ b θ a Common: Tr =1.5s, D = 2.25, = 1, R = 25 Reference: q r () = [,, π 2 ]T, u r = [5, ] T TG 2 TG 1 P x Follower 1: q 1 () = [6, 6, π 2 ]T, (r, l) = (6, 6) Follower 2: q 2 () = [1, 2, 3π 4 ]T, (r, l) = ( 6, 6) Follower 2 Fig.6. Target TG 2 is insie accessible area: θ a > θ b. R Assume also that reference robot has the highest priority an follower robot 2 has the lowest one. From Fig. 8, we see that motion of follower robot 2 is kept within the maximum communication range R with

7 271 J. Jongusuk an T. Mita: Tracking Control of Multiple Mobile Robots: A Case Stuy of Inter-Robot a esire formation establishe in the steay state. See snapshot in Fig. 7 for an entire motion of system. Accoring to Table 1, four of collision occurrences between the two followers are foun. The first two occurences fall to the case when target is outsie accessible area, while the successive ones fall to the case when target is insie the area. The number of collision can be ecrease by moifying l 13, l 23 an T r in each occurrence, an here we simply use constant values. On the other han, the case when target is initially locate insie accessible area is one by setting, Follower 1: q1() = [8, 3, 5π 4 ]T 8 Reference Follower 2: q2() = [8, 9, 5π 8 ]T Follower 2 Follower 1 where other parameters are ientical to the last example. Again, the motion of follower robot 2 (Figs. 9, 1) is guarantee to be insie communication range. The esign of l 13, l 23 for this case gives only one collision, referring to Table 2. y Reference 1 6 Follower 2 Follower x Fig. 7. Snapshot of motion (target is outsie accessible area) 3 y Relative istance r-l clearance Max. Range R x Fig. 9. Snapshot of motion (target is insie accessible area) r-irection l-irection time Fig. 8. Distance from reference to follower 2 in r l irection when target is outsie. Table 1. Collision information (Target outsie). Relative istance r-l clearance Max. Range R 6 No. Info. t l 13 l 23 case outsie outsie insie insie -2 r-irection l-irection time Fig. 1. Distance from reference to follower 2 in r l irection when target is insie.

8 Asian Journal of Control, Vol. 4, No. 3, September Table 2. Collision information (Target insie). No. Info. t l 13 l 23 case insie V. CONCLUSION In this paper, we have stuie strategies for controlling multiple mobile robots (three-robot case) focuse on a collision-free movement. The follower robots track the gait of reference robot using virtual robot tracking control. The extene l-l control technique use to avoi the collision of the three-tuple robots is applie uring finite time interval in such a way that the lowest priority robot can move away from the others. After the collision-free conition is establishe, the controller is switche back to the virtual robot tracking control an maintain the esire formation an so on. The convergences an stabilities of both controllers are throughly stuie, leaing to the main issue in Theorem 1. The problem arises from the nonrobustness of the terminal attractor function use in l-l control, an hence before we move to real experiment, the function shoul be substitute by the more robust one. We are currently working with this issue. We also plan to enlarge our framework to the problems of obstacle avoiance an formation changes, of which we believe the successful strategies are strongly base on the geometric calculation as state in Latombe [4]. Accoring to Theorem 1, the extension of our work to the system of more than three robots is possible particularly when we have a number of groups of three-tuple robots to be controlle. These topics are the subject of our future work. APPENDIX A. Proof of VR kinematic (7) Proof. Using real robot kinematic equation (1) an time erivative of a right-han moel (6), we obtain x υ = x r cos θθ l sin θθ = υ cos θ + ω( r cos θ lsinθ) (A1) y υ = y r sin θθ l cos θθ θ υ = θ = υ sin θ + ω( r sin θ + lcosθ) (A2) = ω (A3) Rewrite the above equations in riftless matrix form yiels (7). B. Terminal attractor To establish the convergence of z ei in arbitrary finite time, T r, consier the following lemma. Lemma 1. The origin of a ifferential equation p z = λz q, λ > is a terminal attractor with a finite reaching time p q (A4) T r = z()1 λ(1 p, p q ) q (, 1) (A5) where p, q I +, an I + = {positive integer} Proof. Direct integration of (A4) yiels z(t) 1 pq = z() 1 pq (1 p q )λt Requiring z(t r ) = implies formula (A5). (A6) Note that T r can be chosen arbitrary from the choice of λ. We shoul take care of the choice p, q, which is suggeste in proposition below. Proposition 1. The choice of p an q in (A4) shoul be selecte such that (p, q) Ω. Ω ={(p, q) p < q, q is o } (A7) Accoring to (A4) an (A5), p q an 1 p q must be esigne to be o-root in orer to avoi 2i-root of negative values. In aition with the assumption that p q <1, Ω is establishe. ACKNOWLEDGEMENTS This work is part of Grant-in-Ai for COE Research Project supporte by the Ministry of Eucation, Science, Sports an Culture, Japan. REFERENCE 1. Cao, Y.U., A.S. Fukunaga an A.B. Kahng, Cooperative Mobile Robotics, Autonomous Robots, Vol. 4, pp (1997). 2. Balch, T. an R. Arkin, Communication in Reactive Multiagent Robotic Systems, Autonomous Robots, Vol. 1, pp (1994). 3. LaValle, S. an S. Hutchinson, Optimal Motion Planning For Multiple Robots Having Inepenent

9 273 J. Jongusuk an T. Mita: Tracking Control of Multiple Mobile Robots: A Case Stuy of Inter-Robot Goals, IEEE Trans. Rob. Autom., Vol. 14, No. 6, pp (1998). 4. Latombe, J.C., Motion Planning: A Journey of Robots, Molecules, Digital Actors, an Other Artifacts, Int. J. Rob. Res., Vol. 18, No. 11, pp (1999). 5. Desai, J.P., J. Ostrowski an V. Kumar, Controlling Formations of Multiple Mobile Robots, IEEE Proc. Int. Conf. Rob. Autom., pp (1998). 6. Desai, J.P., J. Ostrowski an V. Kumar, Control of Changes In Formations For A Team Of Mobile Robots, IEEE Proc. Int. Conf. Rob. Autom., pp (1999). 7. Kanayama, Y., Y. Kimura, F. Miyazaki an T. Noguchi, A Stable Tracking Control Metho for an Autonomous Mobile Robot, IEEE Proc. Int. Conf. Rob. Autom., pp (199). 8. Jiang, Z.P. an H. Nijmeijer, Tracking Control of Mobile Robots: A Case Stuy in Backstepping, Automatica, Vol. 33, No. 7, pp (1997). 9. Slotine, J.E. an W. Li, Applie Nonlinear Control, Prentice-Hall, NJ (1991). 1. Kogo, H. an T. Mita, Introuction to Control System Theory 2 n E., (in Japanese), Jikkyo Press, Japan (1997). 11. Canuas e Wit, C., B. Siciliano an G. Bastin, Theory of Robot Control, Springer Press, Lonon (1997). Jurachart Jongusuk receive B.S. egrees in control engineering at King Mongkut s Institute of Technology, Lakrabang, Thailan in The M.S. egree in the same fiel at Tokyo Institute of Technology (TIT), Japan in 2. He is currently a Ph.D. stuent in Mita laboratory, TIT. His research interests inclue nonlinear control, multitple UAVs formation control. Tsutomu Mita receive the Dr. Eng. egree in electrical engineering from the Tokyo Institute of Technology, Tokyo, Japan, in From 1975 to 1995, he was with the Chiba University, Chiba, Japan, where he was initially appointe as a Lecturer an became an Associate Professor an then a Professor in 1979 an 199, respectively. From December 1991 to September 1992, he was a Visiting Fellow at the Australian National University, Canberra. Since 1995, he has been a Professor in the Department of Control an Systems Engineering, Tokyo Institute of Technology. His research an eucational interests are control theory an its applications to inustrial esign, an robotics. He is the author of several textbooks, incluing An Introuction to H Control (Tokyo, Japan: Shoko-o, 1994).

A Path Planning Method Using Cubic Spiral with Curvature Constraint

A Path Planning Method Using Cubic Spiral with Curvature Constraint A Path Planning Metho Using Cubic Spiral with Curvature Constraint Tzu-Chen Liang an Jing-Sin Liu Institute of Information Science 0, Acaemia Sinica, Nankang, Taipei 5, Taiwan, R.O.C., Email: hartree@iis.sinica.eu.tw

More information

Multi-robot Formation Control Using Reinforcement Learning Method

Multi-robot Formation Control Using Reinforcement Learning Method Multi-robot Formation Control Using Reinforcement Learning Metho Guoyu Zuo, Jiatong Han, an Guansheng Han School of Electronic Information & Control Engineering, Beijing University of Technology, Beijing

More information

VIRTUAL STRUCTURE BASED SPACECRAFT FORMATION CONTROL WITH FORMATION FEEDBACK

VIRTUAL STRUCTURE BASED SPACECRAFT FORMATION CONTROL WITH FORMATION FEEDBACK AIAA Guiance, Navigation, an Control Conference an Exhibit 5-8 August, Monterey, California AIAA -9 VIRTUAL STRUCTURE BASED SPACECRAT ORMATION CONTROL WITH ORMATION EEDBACK Wei Ren Ranal W. Bear Department

More information

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson JUST THE MATHS UNIT NUMBER 10.2 DIFFERENTIATION 2 (Rates of change) by A.J.Hobson 10.2.1 Introuction 10.2.2 Average rates of change 10.2.3 Instantaneous rates of change 10.2.4 Derivatives 10.2.5 Exercises

More information

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments 2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor

More information

Distributed coordination control for multi-robot networks using Lyapunov-like barrier functions

Distributed coordination control for multi-robot networks using Lyapunov-like barrier functions IEEE TRANSACTIONS ON 1 Distribute coorination control for multi-robot networks using Lyapunov-like barrier functions Dimitra Panagou, Dušan M. Stipanović an Petros G. Voulgaris Abstract This paper aresses

More information

Laplacian Cooperative Attitude Control of Multiple Rigid Bodies

Laplacian Cooperative Attitude Control of Multiple Rigid Bodies Laplacian Cooperative Attitue Control of Multiple Rigi Boies Dimos V. Dimarogonas, Panagiotis Tsiotras an Kostas J. Kyriakopoulos Abstract Motivate by the fact that linear controllers can stabilize the

More information

Nonlinear Adaptive Ship Course Tracking Control Based on Backstepping and Nussbaum Gain

Nonlinear Adaptive Ship Course Tracking Control Based on Backstepping and Nussbaum Gain Nonlinear Aaptive Ship Course Tracking Control Base on Backstepping an Nussbaum Gain Jialu Du, Chen Guo Abstract A nonlinear aaptive controller combining aaptive Backstepping algorithm with Nussbaum gain

More information

Exponential Tracking Control of Nonlinear Systems with Actuator Nonlinearity

Exponential Tracking Control of Nonlinear Systems with Actuator Nonlinearity Preprints of the 9th Worl Congress The International Feeration of Automatic Control Cape Town, South Africa. August -9, Exponential Tracking Control of Nonlinear Systems with Actuator Nonlinearity Zhengqiang

More information

6 General properties of an autonomous system of two first order ODE

6 General properties of an autonomous system of two first order ODE 6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x

More information

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:

More information

1 Lecture 20: Implicit differentiation

1 Lecture 20: Implicit differentiation Lecture 20: Implicit ifferentiation. Outline The technique of implicit ifferentiation Tangent lines to a circle Derivatives of inverse functions by implicit ifferentiation Examples.2 Implicit ifferentiation

More information

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control 19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior

More information

Switching Time Optimization in Discretized Hybrid Dynamical Systems

Switching Time Optimization in Discretized Hybrid Dynamical Systems Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set

More information

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+

More information

Experimental Robustness Study of a Second-Order Sliding Mode Controller

Experimental Robustness Study of a Second-Order Sliding Mode Controller Experimental Robustness Stuy of a Secon-Orer Sliing Moe Controller Anré Blom, Bram e Jager Einhoven University of Technology Department of Mechanical Engineering P.O. Box 513, 5600 MB Einhoven, The Netherlans

More information

d dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1

d dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1 Lecture 5 Some ifferentiation rules Trigonometric functions (Relevant section from Stewart, Seventh Eition: Section 3.3) You all know that sin = cos cos = sin. () But have you ever seen a erivation of

More information

FORMATION INPUT-TO-STATE STABILITY. Herbert G. Tanner and George J. Pappas

FORMATION INPUT-TO-STATE STABILITY. Herbert G. Tanner and George J. Pappas Copyright 2002 IFAC 5th Triennial Worl Congress, Barcelona, Spain FORMATION INPUT-TO-STATE STABILITY Herbert G. Tanner an George J. Pappas Department of Electrical Engineering University of Pennsylvania

More information

The Principle of Least Action

The Principle of Least Action Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of

More information

Lecture 2 Lagrangian formulation of classical mechanics Mechanics

Lecture 2 Lagrangian formulation of classical mechanics Mechanics Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,

More information

BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS. Mauro Boccadoro Magnus Egerstedt Paolo Valigi Yorai Wardi

BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS. Mauro Boccadoro Magnus Egerstedt Paolo Valigi Yorai Wardi BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS Mauro Boccaoro Magnus Egerstet Paolo Valigi Yorai Wari {boccaoro,valigi}@iei.unipg.it Dipartimento i Ingegneria Elettronica

More information

Schrödinger s equation.

Schrödinger s equation. Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of

More information

Optimal Variable-Structure Control Tracking of Spacecraft Maneuvers

Optimal Variable-Structure Control Tracking of Spacecraft Maneuvers Optimal Variable-Structure Control racking of Spacecraft Maneuvers John L. Crassiis 1 Srinivas R. Vaali F. Lanis Markley 3 Introuction In recent years, much effort has been evote to the close-loop esign

More information

Total Energy Shaping of a Class of Underactuated Port-Hamiltonian Systems using a New Set of Closed-Loop Potential Shape Variables*

Total Energy Shaping of a Class of Underactuated Port-Hamiltonian Systems using a New Set of Closed-Loop Potential Shape Variables* 51st IEEE Conference on Decision an Control December 1-13 212. Maui Hawaii USA Total Energy Shaping of a Class of Uneractuate Port-Hamiltonian Systems using a New Set of Close-Loop Potential Shape Variables*

More information

Nested Saturation with Guaranteed Real Poles 1

Nested Saturation with Guaranteed Real Poles 1 Neste Saturation with Guarantee Real Poles Eric N Johnson 2 an Suresh K Kannan 3 School of Aerospace Engineering Georgia Institute of Technology, Atlanta, GA 3332 Abstract The global stabilization of asymptotically

More information

Vectors in two dimensions

Vectors in two dimensions Vectors in two imensions Until now, we have been working in one imension only The main reason for this is to become familiar with the main physical ieas like Newton s secon law, without the aitional complication

More information

Introduction to the Vlasov-Poisson system

Introduction to the Vlasov-Poisson system Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its

More information

Distributed Force/Position Consensus Tracking of Networked Robotic Manipulators

Distributed Force/Position Consensus Tracking of Networked Robotic Manipulators 180 IEEE/CAA JOURNAL OF AUTOMATICA SINICA, VOL. 1, NO. 2, APRIL 2014 Distribute Force/Position Consensus Tracking of Networke Robotic Manipulators Lijiao Wang Bin Meng Abstract In this paper, we aress

More information

Dynamic Load Carrying Capacity of Spatial Cable Suspended Robot: Sliding Mode Control Approach

Dynamic Load Carrying Capacity of Spatial Cable Suspended Robot: Sliding Mode Control Approach Int J Avance Design an Manufacturing echnology, Vol. 5/ No. 3/ June - 212 73 Dynamic Loa Carrying Capacity of Spatial Cable Suspene Robot: Sliing Moe Control Approach M. H. Korayem Department of Mechanical

More information

Math 342 Partial Differential Equations «Viktor Grigoryan

Math 342 Partial Differential Equations «Viktor Grigoryan Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite

More information

Summary: Differentiation

Summary: Differentiation Techniques of Differentiation. Inverse Trigonometric functions The basic formulas (available in MF5 are: Summary: Differentiation ( sin ( cos The basic formula can be generalize as follows: Note: ( sin

More information

Appendix: Proof of Spatial Derivative of Clear Raindrop

Appendix: Proof of Spatial Derivative of Clear Raindrop Appenix: Proof of Spatial erivative of Clear Rainrop Shaoi You Robby T. Tan The University of Tokyo {yous,rei,ki}@cvl.iis.u-tokyo.ac.jp Rei Kawakami Katsushi Ikeuchi Utrecht University R.T.Tan@uu.nl Layout

More information

Some Remarks on the Boundedness and Convergence Properties of Smooth Sliding Mode Controllers

Some Remarks on the Boundedness and Convergence Properties of Smooth Sliding Mode Controllers International Journal of Automation an Computing 6(2, May 2009, 154-158 DOI: 10.1007/s11633-009-0154-z Some Remarks on the Bouneness an Convergence Properties of Smooth Sliing Moe Controllers Wallace Moreira

More information

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent

More information

Separation of Variables

Separation of Variables Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical

More information

Table of Common Derivatives By David Abraham

Table of Common Derivatives By David Abraham Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec

More information

TRAJECTORY TRACKING FOR FULLY ACTUATED MECHANICAL SYSTEMS

TRAJECTORY TRACKING FOR FULLY ACTUATED MECHANICAL SYSTEMS TRAJECTORY TRACKING FOR FULLY ACTUATED MECHANICAL SYSTEMS Francesco Bullo Richar M. Murray Control an Dynamical Systems California Institute of Technology Pasaena, CA 91125 Fax : + 1-818-796-8914 email

More information

Differentiability, Computing Derivatives, Trig Review. Goals:

Differentiability, Computing Derivatives, Trig Review. Goals: Secants vs. Derivatives - Unit #3 : Goals: Differentiability, Computing Derivatives, Trig Review Determine when a function is ifferentiable at a point Relate the erivative graph to the the graph of an

More information

Qubit channels that achieve capacity with two states

Qubit channels that achieve capacity with two states Qubit channels that achieve capacity with two states Dominic W. Berry Department of Physics, The University of Queenslan, Brisbane, Queenslan 4072, Australia Receive 22 December 2004; publishe 22 March

More information

u t v t v t c a u t b a v t u t v t b a

u t v t v t c a u t b a v t u t v t b a Nonlinear Dynamical Systems In orer to iscuss nonlinear ynamical systems, we must first consier linear ynamical systems. Linear ynamical systems are just systems of linear equations like we have been stuying

More information

Implicit Differentiation

Implicit Differentiation Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,

More information

Optimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations

Optimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations Optimize Schwarz Methos with the Yin-Yang Gri for Shallow Water Equations Abessama Qaouri Recherche en prévision numérique, Atmospheric Science an Technology Directorate, Environment Canaa, Dorval, Québec,

More information

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x) Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)

More information

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y Ph195a lecture notes, 1/3/01 Density operators for spin- 1 ensembles So far in our iscussion of spin- 1 systems, we have restricte our attention to the case of pure states an Hamiltonian evolution. Toay

More information

Adaptive Gain-Scheduled H Control of Linear Parameter-Varying Systems with Time-Delayed Elements

Adaptive Gain-Scheduled H Control of Linear Parameter-Varying Systems with Time-Delayed Elements Aaptive Gain-Scheule H Control of Linear Parameter-Varying Systems with ime-delaye Elements Yoshihiko Miyasato he Institute of Statistical Mathematics 4-6-7 Minami-Azabu, Minato-ku, okyo 6-8569, Japan

More information

Sliding mode approach to congestion control in connection-oriented communication networks

Sliding mode approach to congestion control in connection-oriented communication networks JOURNAL OF APPLIED COMPUTER SCIENCE Vol. xx. No xx (200x), pp. xx-xx Sliing moe approach to congestion control in connection-oriente communication networks Anrzej Bartoszewicz, Justyna Żuk Technical University

More information

Interpolated Rigid-Body Motions and Robotics

Interpolated Rigid-Body Motions and Robotics Interpolate Rigi-Boy Motions an Robotics J.M. Selig Faculty of Business, Computing an Info. Management. Lonon South Bank University, Lonon SE AA, U.K. seligjm@lsbu.ac.uk Yaunquing Wu Dept. Mechanical Engineering.

More information

The total derivative. Chapter Lagrangian and Eulerian approaches

The total derivative. Chapter Lagrangian and Eulerian approaches Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function

More information

Lecture 6: Control of Three-Phase Inverters

Lecture 6: Control of Three-Phase Inverters Yoash Levron The Anrew an Erna Viterbi Faculty of Electrical Engineering, Technion Israel Institute of Technology, Haifa 323, Israel yoashl@ee.technion.ac.il Juri Belikov Department of Computer Systems,

More information

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21 Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting

More information

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5

More information

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013 Survey Sampling Kosuke Imai Department of Politics, Princeton University February 19, 2013 Survey sampling is one of the most commonly use ata collection methos for social scientists. We begin by escribing

More information

Designing Information Devices and Systems II Fall 2017 Note Theorem: Existence and Uniqueness of Solutions to Differential Equations

Designing Information Devices and Systems II Fall 2017 Note Theorem: Existence and Uniqueness of Solutions to Differential Equations EECS 6B Designing Information Devices an Systems II Fall 07 Note 3 Secon Orer Differential Equations Secon orer ifferential equations appear everywhere in the real worl. In this note, we will walk through

More information

Sensors & Transducers 2015 by IFSA Publishing, S. L.

Sensors & Transducers 2015 by IFSA Publishing, S. L. Sensors & Transucers, Vol. 184, Issue 1, January 15, pp. 53-59 Sensors & Transucers 15 by IFSA Publishing, S. L. http://www.sensorsportal.com Non-invasive an Locally Resolve Measurement of Soun Velocity

More information

Math 115 Section 018 Course Note

Math 115 Section 018 Course Note Course Note 1 General Functions Definition 1.1. A function is a rule that takes certain numbers as inputs an assigns to each a efinite output number. The set of all input numbers is calle the omain of

More information

Section 7.2. The Calculus of Complex Functions

Section 7.2. The Calculus of Complex Functions Section 7.2 The Calculus of Complex Functions In this section we will iscuss limits, continuity, ifferentiation, Taylor series in the context of functions which take on complex values. Moreover, we will

More information

θ x = f ( x,t) could be written as

θ x = f ( x,t) could be written as 9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)

More information

Potential Fields for Maintaining Connectivity of Mobile Networks

Potential Fields for Maintaining Connectivity of Mobile Networks 82 IEEE TRANSACTIONS ON ROBOTICS, VOL. 2, NO. 4, AUGUST 2007 9] N. Yanai, M. Yamamoto, an Akira Mohri, Anti-sway control for wiresuspene mechanism base on ynamics compensation, in Proc. Int. Conf. Robot.

More information

Multi-agent Systems Reaching Optimal Consensus with Time-varying Communication Graphs

Multi-agent Systems Reaching Optimal Consensus with Time-varying Communication Graphs Preprints of the 8th IFAC Worl Congress Multi-agent Systems Reaching Optimal Consensus with Time-varying Communication Graphs Guoong Shi ACCESS Linnaeus Centre, School of Electrical Engineering, Royal

More information

Lecture 6: Calculus. In Song Kim. September 7, 2011

Lecture 6: Calculus. In Song Kim. September 7, 2011 Lecture 6: Calculus In Song Kim September 7, 20 Introuction to Differential Calculus In our previous lecture we came up with several ways to analyze functions. We saw previously that the slope of a linear

More information

G j dq i + G j. q i. = a jt. and

G j dq i + G j. q i. = a jt. and Lagrange Multipliers Wenesay, 8 September 011 Sometimes it is convenient to use reunant coorinates, an to effect the variation of the action consistent with the constraints via the metho of Lagrange unetermine

More information

NOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy,

NOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy, NOTES ON EULER-BOOLE SUMMATION JONATHAN M BORWEIN, NEIL J CALKIN, AND DANTE MANNA Abstract We stuy a connection between Euler-MacLaurin Summation an Boole Summation suggeste in an AMM note from 196, which

More information

Design A Robust Power System Stabilizer on SMIB Using Lyapunov Theory

Design A Robust Power System Stabilizer on SMIB Using Lyapunov Theory Design A Robust Power System Stabilizer on SMIB Using Lyapunov Theory Yin Li, Stuent Member, IEEE, Lingling Fan, Senior Member, IEEE Abstract This paper proposes a robust power system stabilizer (PSS)

More information

We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors

We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors We are IntechOpen, the worl s leaing publisher of Open Access books Built by scientists, for scientists 4,100 116,000 120M Open access books available International authors an eitors Downloas Our authors

More information

Commun Nonlinear Sci Numer Simulat

Commun Nonlinear Sci Numer Simulat Commun Nonlinear Sci Numer Simulat 14 (2009) 3901 3913 Contents lists available at ScienceDirect Commun Nonlinear Sci Numer Simulat journal homepage: www.elsevier.com/locate/cnsns Dynamics of a ring network

More information

Power Generation and Distribution via Distributed Coordination Control

Power Generation and Distribution via Distributed Coordination Control Power Generation an Distribution via Distribute Coorination Control Byeong-Yeon Kim, Kwang-Kyo Oh, an Hyo-Sung Ahn arxiv:407.4870v [math.oc] 8 Jul 204 Abstract This paper presents power coorination, power

More information

Application of Nonlinear Control to a Collision Avoidance System

Application of Nonlinear Control to a Collision Avoidance System Application of Nonlinear Control to a Collision Avoiance System Pete Seiler Mechanical Eng. Dept. U. of California-Berkeley Berkeley, CA 947-74, USA Phone: -5-64-6933 Fax: -5-64-663 pseiler@euler.me.berkeley.edu

More information

A New Nonlinear H-infinity Feedback Control Approach to the Problem of Autonomous Robot Navigation

A New Nonlinear H-infinity Feedback Control Approach to the Problem of Autonomous Robot Navigation Intell In Syst (15) 1:179 186 DOI 1.17/s493-15-1- ORIGINAL PAPER A New Nonlinear H-infinity Feeback Control Approach to the Problem of Autonomous Robot Navigation G. Rigatos 1 P. Siano Receive: 8 December

More information

Optimal Control of Spatially Distributed Systems

Optimal Control of Spatially Distributed Systems Optimal Control of Spatially Distribute Systems Naer Motee an Ali Jababaie Abstract In this paper, we stuy the structural properties of optimal control of spatially istribute systems. Such systems consist

More information

Control of Generalized Contact Motion and Force in Physical Human-Robot Interaction

Control of Generalized Contact Motion and Force in Physical Human-Robot Interaction 15 IEEE International Conference on Robotics an Automation ICRA Washington State Convention Center Seattle, Washington, May -3, 15 Control of Generalize Contact Motion an Force in Physical Human-Robot

More information

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1 Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition

More information

and from it produce the action integral whose variation we set to zero:

and from it produce the action integral whose variation we set to zero: Lagrange Multipliers Monay, 6 September 01 Sometimes it is convenient to use reunant coorinates, an to effect the variation of the action consistent with the constraints via the metho of Lagrange unetermine

More information

12.5. Differentiation of vectors. Introduction. Prerequisites. Learning Outcomes

12.5. Differentiation of vectors. Introduction. Prerequisites. Learning Outcomes Differentiation of vectors 12.5 Introuction The area known as vector calculus is use to moel mathematically a vast range of engineering phenomena incluing electrostatics, electromagnetic fiels, air flow

More information

Make graph of g by adding c to the y-values. on the graph of f by c. multiplying the y-values. even-degree polynomial. graph goes up on both sides

Make graph of g by adding c to the y-values. on the graph of f by c. multiplying the y-values. even-degree polynomial. graph goes up on both sides Reference 1: Transformations of Graphs an En Behavior of Polynomial Graphs Transformations of graphs aitive constant constant on the outsie g(x) = + c Make graph of g by aing c to the y-values on the graph

More information

VISUAL SERVOING WITH ORIENTATION LIMITS OF A X4-FLYER

VISUAL SERVOING WITH ORIENTATION LIMITS OF A X4-FLYER VISUAL SERVOING WITH ORIENTATION LIMITS OF A X4-FLYER Najib Metni,Tarek Hamel,Isabelle Fantoni Laboratoire Central es Ponts et Chaussées, LCPC-Paris France, najib.metni@lcpc.fr Cemif-Sc FRE-CNRS 2494,

More information

Equilibrium in Queues Under Unknown Service Times and Service Value

Equilibrium in Queues Under Unknown Service Times and Service Value University of Pennsylvania ScholarlyCommons Finance Papers Wharton Faculty Research 1-2014 Equilibrium in Queues Uner Unknown Service Times an Service Value Laurens Debo Senthil K. Veeraraghavan University

More information

OPTIMAL CONTROL PROBLEM FOR PROCESSES REPRESENTED BY STOCHASTIC SEQUENTIAL MACHINE

OPTIMAL CONTROL PROBLEM FOR PROCESSES REPRESENTED BY STOCHASTIC SEQUENTIAL MACHINE OPTIMA CONTRO PROBEM FOR PROCESSES REPRESENTED BY STOCHASTIC SEQUENTIA MACHINE Yaup H. HACI an Muhammet CANDAN Department of Mathematics, Canaale Onseiz Mart University, Canaale, Turey ABSTRACT In this

More information

Predictive control of synchronous generator: a multiciterial optimization approach

Predictive control of synchronous generator: a multiciterial optimization approach Preictive control of synchronous generator: a multiciterial optimization approach Marián Mrosko, Eva Miklovičová, Ján Murgaš Abstract The paper eals with the preictive control esign for nonlinear systems.

More information

An inductance lookup table application for analysis of reluctance stepper motor model

An inductance lookup table application for analysis of reluctance stepper motor model ARCHIVES OF ELECTRICAL ENGINEERING VOL. 60(), pp. 5- (0) DOI 0.478/ v07-0-000-y An inuctance lookup table application for analysis of reluctance stepper motor moel JAKUB BERNAT, JAKUB KOŁOTA, SŁAWOMIR

More information

initial configuration initial configuration end

initial configuration initial configuration end Design of trajectory stabilizing feeback for riftless at systems M. FLIESS y J. L EVINE z P. MARTIN x P. ROUCHON { ECC95 Abstract A esign metho for robust stabilization of riftless at systems aroun trajectories

More information

Situation awareness of power system based on static voltage security region

Situation awareness of power system based on static voltage security region The 6th International Conference on Renewable Power Generation (RPG) 19 20 October 2017 Situation awareness of power system base on static voltage security region Fei Xiao, Zi-Qing Jiang, Qian Ai, Ran

More information

TMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments

TMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments Problem F U L W D g m 3 2 s 2 0 0 0 0 2 kg 0 0 0 0 0 0 Table : Dimension matrix TMA 495 Matematisk moellering Exam Tuesay December 6, 2008 09:00 3:00 Problems an solution with aitional comments The necessary

More information

Calculus of Variations

Calculus of Variations 16.323 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t

More information

Math 1272 Solutions for Spring 2005 Final Exam. asked to find the limit of the sequence. This is equivalent to evaluating lim. lim.

Math 1272 Solutions for Spring 2005 Final Exam. asked to find the limit of the sequence. This is equivalent to evaluating lim. lim. Math 7 Solutions for Spring 5 Final Exam ) We are gien an infinite sequence for which the general term is a n 3 + 5n n + n an are 3 + 5n aske to fin the limit of the sequence. This is equialent to ealuating

More information

Least-Squares Regression on Sparse Spaces

Least-Squares Regression on Sparse Spaces Least-Squares Regression on Sparse Spaces Yuri Grinberg, Mahi Milani Far, Joelle Pineau School of Computer Science McGill University Montreal, Canaa {ygrinb,mmilan1,jpineau}@cs.mcgill.ca 1 Introuction

More information

Differentiability, Computing Derivatives, Trig Review

Differentiability, Computing Derivatives, Trig Review Unit #3 : Differentiability, Computing Derivatives, Trig Review Goals: Determine when a function is ifferentiable at a point Relate the erivative graph to the the graph of an original function Compute

More information

Hyperbolic Moment Equations Using Quadrature-Based Projection Methods

Hyperbolic Moment Equations Using Quadrature-Based Projection Methods Hyperbolic Moment Equations Using Quarature-Base Projection Methos J. Koellermeier an M. Torrilhon Department of Mathematics, RWTH Aachen University, Aachen, Germany Abstract. Kinetic equations like the

More information

12.11 Laplace s Equation in Cylindrical and

12.11 Laplace s Equation in Cylindrical and SEC. 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential 593 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential One of the most important PDEs in physics an engineering

More information

A Comparison between a Conventional Power System Stabilizer (PSS) and Novel PSS Based on Feedback Linearization Technique

A Comparison between a Conventional Power System Stabilizer (PSS) and Novel PSS Based on Feedback Linearization Technique J. Basic. Appl. Sci. Res., ()9-99,, TextRoa Publication ISSN 9-434 Journal of Basic an Applie Scientific Research www.textroa.com A Comparison between a Conventional Power System Stabilizer (PSS) an Novel

More information

Robust Tracking Control of Robot Manipulator Using Dissipativity Theory

Robust Tracking Control of Robot Manipulator Using Dissipativity Theory Moern Applie Science July 008 Robust racking Control of Robot Manipulator Using Dissipativity heory Hongrui Wang Key Lab of Inustrial Computer Control Engineering of Hebei Province Yanshan University Qinhuangao

More information

An Approach for Design of Multi-element USBL Systems

An Approach for Design of Multi-element USBL Systems An Approach for Design of Multi-element USBL Systems MIKHAIL ARKHIPOV Department of Postgrauate Stuies Technological University of the Mixteca Carretera a Acatlima Km. 2.5 Huajuapan e Leon Oaxaca 69000

More information

A check digit system over a group of arbitrary order

A check digit system over a group of arbitrary order 2013 8th International Conference on Communications an Networking in China (CHINACOM) A check igit system over a group of arbitrary orer Yanling Chen Chair of Communication Systems Ruhr University Bochum

More information

Application of the homotopy perturbation method to a magneto-elastico-viscous fluid along a semi-infinite plate

Application of the homotopy perturbation method to a magneto-elastico-viscous fluid along a semi-infinite plate Freun Publishing House Lt., International Journal of Nonlinear Sciences & Numerical Simulation, (9), -, 9 Application of the homotopy perturbation metho to a magneto-elastico-viscous flui along a semi-infinite

More information

6. Friction and viscosity in gasses

6. Friction and viscosity in gasses IR2 6. Friction an viscosity in gasses 6.1 Introuction Similar to fluis, also for laminar flowing gases Newtons s friction law hols true (see experiment IR1). Using Newton s law the viscosity of air uner

More information

Optimal Control of Spatially Distributed Systems

Optimal Control of Spatially Distributed Systems Optimal Control of Spatially Distribute Systems Naer Motee an Ali Jababaie Abstract In this paper, we stuy the structural properties of optimal control of spatially istribute systems. Such systems consist

More information

Diophantine Approximations: Examining the Farey Process and its Method on Producing Best Approximations

Diophantine Approximations: Examining the Farey Process and its Method on Producing Best Approximations Diophantine Approximations: Examining the Farey Process an its Metho on Proucing Best Approximations Kelly Bowen Introuction When a person hears the phrase irrational number, one oes not think of anything

More information

An M/G/1 Retrial Queue with Priority, Balking and Feedback Customers

An M/G/1 Retrial Queue with Priority, Balking and Feedback Customers Journal of Convergence Information Technology Volume 5 Number April 1 An M/G/1 Retrial Queue with Priority Balking an Feeback Customers Peishu Chen * 1 Yiuan Zhu 1 * 1 Faculty of Science Jiangsu University

More information

Mathematics 116 HWK 25a Solutions 8.6 p610

Mathematics 116 HWK 25a Solutions 8.6 p610 Mathematics 6 HWK 5a Solutions 8.6 p6 Problem, 8.6, p6 Fin a power series representation for the function f() = etermine the interval of convergence. an Solution. Begin with the geometric series = + +

More information

Spectral properties of a near-periodic row-stochastic Leslie matrix

Spectral properties of a near-periodic row-stochastic Leslie matrix Linear Algebra an its Applications 409 2005) 66 86 wwwelseviercom/locate/laa Spectral properties of a near-perioic row-stochastic Leslie matrix Mei-Qin Chen a Xiezhang Li b a Department of Mathematics

More information

Harmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method

Harmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method 1 Harmonic Moelling of Thyristor Briges using a Simplifie Time Domain Metho P. W. Lehn, Senior Member IEEE, an G. Ebner Abstract The paper presents time omain methos for harmonic analysis of a 6-pulse

More information