VIRTUAL STRUCTURE BASED SPACECRAFT FORMATION CONTROL WITH FORMATION FEEDBACK

Size: px
Start display at page:

Download "VIRTUAL STRUCTURE BASED SPACECRAFT FORMATION CONTROL WITH FORMATION FEEDBACK"

Transcription

1 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 of Electrical an Computer Engineering Brigham Young University Provo, UT 8 fweiren, bearg@ee.byu.eu Abstract ormation control for multiple vehicles has become an active research area in recent years. Generally there are three approaches to this problem, namely leaerfollowing, behavioral, an virtual structure approaches. In this paper, formation control ieas for multiple spacecraft using virtual structure approach are presente. If there is no formation feeback from spacecraft to the virtual structure, the spacecraft will get out of formation when the virtual structure moves too fast for the spacecraft to track or the total system must sacrifice convergence spee in orer to keep the spacecraft in formation. The spacecraft may also get out of formation when the system is affecte by internal or external isturbances. To overcome these rawbacks, a novel way of introucing formation feeback from spacecraft to the virtual structure is illustrate in etail. An application of these ieas to multiple spacecraft interferometers in eep space is given. Introuction ormation control for multiple vehicles has become an active research area in recent years. Applications in this area inclue the coorination of multiple robots, UAVs, satellites, aircraft, an spacecraft. 5 Many papers have been publishe to eal with ifferent control strategies, schemes, an applications of multiple vehicle control. While the applications are ifferent, the funamental approaches to formation control are similar: the common theme being the coorination of multiple vehicles to accomplish an objective. Generally there are three approaches to multi-vehicle coorination reporte in the literature, namely leaerfollowing, behavioral, an virtual structure approaches. In the leaer following approach,,, one of the agents is esignate as the leaer, with the rest of the agents esignate as followers. The leaer tracks a pre-efine trajectory, an the followers track a transforme version of the leaer s states. The avantage of leaer following is that group behavior is irecte by specifying the behavior of a single quantity: the leaer. The isavantage is that there is no explicit feeback to the formation. Another isavantage is that the leaer is a single point of failure for the formation. In the remainer of the paper, we use the following efinition: 7 The group feeback from vehicles to the formation is referre to as formation feeback. In the behavioral approach, 8 several esire behaviors are prescribe for each agent. The basic iea is to make the control action of each agent a weighte average of the control for each behavior. Possible behaviors inclue collision avoiance, obstacle avoiance, goal seeking, an formation keeping. The avantage of the behavioral approach is that it is natural to erive control strategies when agents have multiple competing objectives. In aition, there is explicit feeback to the formation since each agent reacts accoring to the position of its neighbors. Another avantage is that the behavioral approach lens itself naturally to a ecentralize implementation. The isavantage is that the group behavior cannot be explicitly efine, rather the group behavior is sai to emerge. In aition, it is ifficult to analyze the behavioral approach mathematically an guarantee its group stability. In the virtual structure approach,, 7, the entire formation is treate as a single structure. The virtual structure can evolve as a rigi boy in a given irection with some given orientation an maintain a rigi geometric relationship among multiple vehicles. The avantage of the virtual structure approach is that it is fairly easy to prescribe a coorinate behavior of the group. The isavantage is that requiring the formation to act as a virtual structure limits the class of potential applications of this approach. Another isavantage is that its current evelopment lens itself to a centralize control implementation. In the case of the application of synthesizing multiple Copyright by the, Inc. All rights reserve.

2 spacecraft interferometers in eep space, it is esirable to have a constellation of spacecraft act as a single rigi boy in orer to image stars in eep space. As a result, it is suitable to choose the virtual structure approach to accomplish formation maneuvers. In the remainer of the paper we use the term formation an virtual structure interchangeably. In general, there is a ilemma when there is no feeback applie from spacecraft to the virtual structure. On the one han, if the virtual structure evolves too fast, the spacecraft cannot track their esire trajectories accurately an they will get out of formation. On the other han, the virtual structure might be slowe own sufficiently to allow the spacecraft to track their trajectories accurately. However, this results in unreasonably slow formation ynamics. Also, when performing formation maneuvers, the total system is often isturbe by internal or external factors. or example, some spacecraft may fail for a perio of time ue to mechanical or electrical malfunctions or may isintegrate from the formation ue to external isturbances in eep space. If there is no formation feeback from spacecraft to the formation, the faile or isintegrate spacecraft will be left behin while the other spacecraft still keep moving towars their final goals, an the entire system cannot ajust to maintain formation. ormation feeback from spacecraft to the virtual structure provies a goo compromise between formation keeping an convergence spee as well as improve group stability an robustness. In Ref. the authors introuce a coorination architecture for spacecraft formation control which subsumes leaer-following, behavioral, an virtual structure approaches to the multi-agent coorination problem. In Ref. formation maneuvers are easily prescribe an group stability is guarantee, but formation feeback is not inclue. In Ref. 7 formation feeback is use for the coorinate control problem for multiple robots. This paper is aime as a further evelopment for Ref. an Ref. 7. The main contribution of this paper is to propose a novel iea of introucing formation feeback from spacecraft to the virtual structure an apply this iea to the spacecraft interferometry problem so that formation keeping is guarantee throughout the maneuver an the total system robustness is improve. The ecentralize control implementation of the virtual structure approach nees to be explore in the future. The outline of this paper is as follows. In section we introuce spacecraft ynamics. In section we escribe virtual structure equations of motion for spacecraft. In section we present formation control strategies with formation feeback. In section 5 we illustrate simulation results for spacecraft formation control. By comparing the results with an without formation feeback, we emonstrate the superiority of the system with formation feeback over the one without formation feeback. Spacecraft Dynamics In this paper each spacecraft is moele as a rigi boy, with r i, v i, q i an! i representing the position, velocity, unit quaternion, an angular velocity of the ith spacecraft, where r i, v i, an! i are vectors an q i is a unit quaternion use to represent the attitue of a rigi boy. We will represent r i, v i, an! i in terms of their components in the inertial frame C O. or simplicity, we use the same symbol to enote a vector an its corresponing coorinate frame representation in the remainer of the paper. Euler s theorem for rigi boy rotations states that the general isplacement of a rigi boy with one point fixe is a rotation about some axis. Let z represent a unit vector in the irection of rotation, calle the eigenaxis, an let represent the angle of rotation about z, calle the Euler angle. The unit quaternion representing this rotation is given by q = [z T sin(=); cos(=)] T = [~q T ; q] T, where ~q is a vector with its components represente in the given coorinate frame an q is a scalar. It is easy to see that q an q represent the same attitue. To simplify our iscussion in the remainer of the paper, we assume q. Given a vector p, the corresponing cross-prouct operation p is efine as p = p p p p p p where p =[p ; p ; p ] T in terms of its components in the given coorinate frame. If we let C O O be a rotation matrix that represents the orientation of the frame C O with respect to C O, then r O = C O Or O, where r O an r O represent vector r in terms of its components in the frame C O an C O separately. The relationship between the unit quaternion q represente in the frame C O an the rotation matrix C O O is efine as 5 ; C O O =(q )I +~q ~q T q~q : The relationship between two rotation matrices C O O an C OO is given as C OO = C T O O =(q )I +~q ~q T +q~q : The multiplication of two quaternions is given by the formula q a q b = Q(q b )q a, where Q(q b ) = q b I ~q b ~q b ~q T : Let q b q be the inverse of the quaternion q given by q ~q ~q b = =. Suppose that q q the unit quaternions q an q represent the actual attitue

3 an the esire attitue of a rigi boy respectively, then the attitue error is given by q e = q q ~qe = : q e The translational ynamics for the spacecraft are C C C _r i =v i M i _v i =f i ; () where M i is the mass of the ith spacecraft, an f i is the control force. The rotational ynamics for the spacecraft are _q i =! i ~q i + q i! i _q i =! i ~q i () C O given by igure : Coorinate frame geometry. C J i _! i =! i J i! i + i ; where J i is the moment of inertia of the ith spacecraft, an i is the control torque on the ith spacecraft. Virtual Structure Equations Of Motion or Spacecraft In the virtual structure approach, we treat the entire formation as a rigi boy with place-holers fixe in the formation to represent the esire position an attitue for each spacecraft. As the virtual structure evolves in time, the place-holers trace out trajectories for each spacecraft to track. The relative orientation of each place-holer within the formation is fixe with respect to each other. 7 The control is erive in four steps: first, the esire ynamics of the virtual structure are efine, secon, the motion of the virtual structure is translate into the esire motion for each spacecraft, thir, tracking controls for each spacecraft are erive, an finally, formation feeback is introuce from each spacecraft to the virtual structure. The coorinate frame geometry is shown in igure. rame C O is an inertial frame. Since the formation can be thought of as a rigi boy with inertial position r, velocity v, attitue q, an angular velocity!, we efine the formation reference frame C locate at r with an orientation given by q with respect to the inertial frame C O. We also have one reference frame C i imbee in each spacecraft to represent the configuration of each spacecraft. Each spacecraft can be represente either by position r i, velocity v i, unit quaternion attitue q i, an angular velocity! i with respect to the inertial frame C O or by r i, v i, q i, an! i with respect to the formation reference frame C. Virtual structure equations of motion for each placeholer, that is, the esire motion for each spacecraft are r i (t) =r (t) +C O (t)r i (t) v i (t) =v (t) +C O (t)v i (t) +! (t) (C O (t)r i (t)) () q i (t) =q (t)q i (t)! i (t) =! (t) +C O (t)! i (t); where C O (t) is the rotation matrix of the frame C O with respect to C, an superscript above a vector means a esire state for each spacecraft. C O is given as C O = C T O =(q )I +~q ~q T +q ~q : Group maneuvers that preserve formation shape can be achieve as a succession of elementary formation maneuvers. Therefore, we will introuce virtual structure equations of motion for spacecraft via elementary formation maneuvers. The elementary formation maneuvers inclue translations, rotations, an expansions/contractions. Let (t) = [ (t); (t); (t)] T with its components represent the expansion/contraction rates along each formation reference frame axis. An expansion/contraction matrix is efine as (t) =iag((t)), which is a iagonal matrix. Generally all parameters in () can vary with time. However, if we are concerne with formation maneuvers that preserve the overall formation shape, ri, v, i qi, an! i are constant. To realize elementary formation maneuvers, we can vary r an v to translate the formation, vary q (q can be transforme to the rotation matrix C O.) an! to rotate the formation, an replace ri in the first an secon equations in () by (t)r i an replace v i in the secon equation in () by _(t)r i to expan or contract the formation. Arbitrary formation maneuvers can be realize by varying r (t), v (t), q (t),! (t), (t), an (t) _ simultaneously. In the case of preserving the overall formation shape uring

4 formation maneuvers, the equations of motion are given as r i (t) =r (t) +C O (t)(t)r i v i (t) =v (t) +C O (t) _ (t)r i q i (t) =q (t)q i! i (t) =! (t): +! (t) (C O (t)(t)r i ) () Note that! i is zero since q i is constant. The erivatives of the esire states are given by _r i (t) =_r (t) + _ C O (t)(t)r i + C O (t) _ (t)r i _v i (t) =_v (t) + _ C O (t) _ (t)r i + C O (t) (t)r i + _! (t) (C O (t)(t)r i ) (5) +! (t) ( _ C O (t)(t)r i + C O _ (t)r i ) _q i (t) =_q (t)q i _! i (t) =_! (t): rom (), we can see that if the velocity of the formation is zero, that is, v (t) =,! (t) =, an (t) _ =, then the esire velocity of each spacecraft is zero, that is, vi (t) =an! i (t) =. Also, if both the velocity an acceleration of the formation are zero, that is, v (t) =,! (t) =, (t) _ =, _v (t) =, _! (t) =, an (t) =, then both the esire velocity an acceleration of each spacecraft are zero, that is, vi (t) =,! i (t) =, _v i (t) =, an _! i (t) =. ormation Control Strategies With ormation eeback Let X i = [r T i ; vt i ; qt i ;!T i ]T an X i = [ri T ; v T i ; q T i ;! T i ] T represent the states an esire states for the ith spacecraft with respect to the inertial frame C O repectively. Let X i = [ri T ; vt i ; qt i ;!T i ]T T ;v i T ;q i T ;! i T ] T an Xi = [ri represent the states an esire states for the ith spacecraft with respect to the formation frame C respectively. Let X = [r T ;vt ;qt ;!T ;T ; _ T ]T an X = [r T ;vt ;qt ;!T ;T ; _ T ]T represent respectively the states an esire states for the virtual structure with respect to the inertial frame C O. We know that Xi is constant since we want to preserve the formation shape uring the group maneuvers. The aim of the formation maneuver is to evolve X (t) to X (t) while guaranteeing that X i(t) tracks Xi (t). Accoringly, a formation maneuver is efine as follows: A formation maneuver is asymptotically achieve if X (t)! X (t) an X i(t)! Xi (t) as t!. The control law for each spacecraft without formation feeback is given by the following lemma. Lemma. Let X = [r T ; v T ; q T ;! T ] T, X = [r T ; v T ; q T ;! T ] T, an let _r _v _q J _! C A = where (!) = If _v K r (r r ) K v (v v ) (!)q! J! + K q ~q e K! (!! )!!! T an q e = q q.. K r = K T r >, K v = K T v >, K q = K T q >, K! = K T! >,. r L [; ) \ L [; ),. _! +! L [; ) \ L [; ), then X X! as t!. C A ; () Proof: see Ref. an Ref.. In the case of the control law for the virtual structure, we nee to a two equations to () since an _ are use to represent the expansion/contraction rate of the formation, an pairs (K r ; K v ), (k q ; K! ), (K ; K _ ) correspon to translation, rotation, an expansion/contraction gains for the formation respectively. We also assume that q an! satisfy the rotational ynamics. Note that in the simple case when X is constant, the rotational ynamics is satisfie obviously. The control law for the virtual structure is given as follows. Lemma. Let _r _v _q C B _! _ = C B _v _v K r(r r ) K v(v v ) (! )q _! + k q~q e K! (!! ) _ K ( ) K _( _ ; C A _ ) (7) where q e = q q. If K r, K v, K!, K, an K _ are symmetric positive efinite matrices, an k q is a positive scalar, then X X! as t!. Proof: We can rewrite the secon equation in (7) as r r = K r(r r ) K v(_r _r ): Let ~r = r r, then ~r = K r ~r K v _ ~r. Since K r an K v are positive efinite, it is obvious to see that

5 r r! an v v! asymptotically as t!. Rewriting the thir equation in (7) as _~q = (q!! ~q ) _q =!T ~q : (8) Base on the assumption above, this equation is also true for q, so we can get _~q = (q!! ~q ) _q =!T ~q : (9) Let ~q = q q an ~! =!!. rom (8) an (9), we know that _~q = (q!! ~q q! +! ~q ) _~q = (!T ~q! T ~q ): () Let V = ~q T ~q, V = ~!T ~!, an consier the Lyapunov function caniate: V = k q V + V : Differentiating V, we get _V = ~q T _ ~q = ~! T (q ~q ~q ~q q ~q ): After some manipulation, we also know that q e = q q q = ~q + ~q ~q +q ~q ~q T ~q + q q which means that ~q e = q ~q +~q ~q +q ~q. Thus, _V = ~! T ~q e. Rewriting the fourth equation in (7) as ~! _ = k q ~q e K! ~!, an ifferentiating V, we can arrive at V _ = ~! T ~! _ = ~! T (k q~q e K! ~!). Therefore, V _ = k _ q V + _V = ~! T K! ~!. Let =f(~q ; ~! )j V _ =g, an be the largest invariant set containe in.on, V _ =, which implies that ~! since K! is positive efinite. When we plug! =! into the fourth equation in (7), we can show that ~q e = [; ; ] T, which implies that q = q. Therefore, by LaSalle s invariance principle, q q! an!!! asymptotically as t!. Thus X X! as t!. Therefore, X (t)! X (t) an X i(t)! Xi (t) as t!, an formation maneuvers without formation feeback are asymptotically achieve. ; or a secon orer system s + k s + k =,if we efine rise time t r an amping ratio, then natural frequency! n is approximately :8=t r. Therefore, if we let k =! n = (:8=t r) an k =! n = (:8=t r ), the transient specifications for the system are satisfie. We can efine K r, k q, K accoring to k, an efine K v, K!, K _ accoring to k. rom Lemma. an., we can see that the performe maneuver will be achieve an the spacecraft will track their esire states in the en. However, how well the spacecraft will preserve the formation shape uring the maneuver is not guarantee by this control law. or example, the errors r i (t) ri (t) an q i (t) qi (t) for the ith spacecraft may be large uring the maneuver, that is, the spacecraft may get out of the esire formation. If the virtual structure moves too fast, the spacecraft coul fall far behin their esire positions. If the virtual structure moves too slowly, the maneuver cannot be achieve within a short time. Therefore, we introuce formation feeback from the spacecraft to the virtual structure to overcome these rawbacks. We will introuce nonlinear gains in the control law for the virtual structure. Let X = [X T ;XT ; ;XT N ]T an X = [X T ;X T ; ;X T N ]T, where N is the number of spacecraft in the formation. The performance measure is efine as X X. We woul like to esign the nonlinear gains to meet the following requirements. When the spacecraft are out of the esire formation, that is, X X is large, the virtual structure will slow own or stop, allowing the spacecraft to regain formation. When the spacecraft are maintaining formation, that is, X X is small, the virtual structure will keep moving towar its final goal. A caniate for such gains can be efine as = K+K X X, where K = K T > is the gain when there is no formation feeback, an K = K T > is the formation gain which weights the performance measure X X.We can see that X X! )! K X X!)!: We can use nonlinear gains v,!, an _ to replace K v, K!, an K _ in (7), where nonlinear gains are efine as follows. v =K v + K X X! =K! + K X X () _ =K _ + K X X : Of course, we can use ifferent K an rise times for pairs (K r ; v ), (k q ; q ), an (K ; _ ) to change the weights of translation, rotation, an expansion/contraction effects. As a result, nonlinear gains slow 5

6 own or spee up the virtual structure base on how far out of the esire formation the spacecraft are. The control law for the virtual structure with formation feeback is given as follows. Lemma. Let _r _v _q _! _ = C B _v _v K r(r r ) v(v v ) (! )q _! + k q~q e! (!! ) _ K ( ) _ ( _ _ ) ; C A () where q e = q q. If v,!, an _ are given by (), then X X! as t!. Proof: We can follow the same proceure as Lemma. except that we use nonlinear gains v,!, an _ to replace the linear gains K v, K!, an K respectively _ everywhere in the proof. Since v,!, an _ are positive efinite, we can show that X X! as t!. Combine with the control law for each spacecraft, we know that X (t)! X (t) an X i(t)! Xi (t) as t!. ormation maneuvers with formation feeback are asymptotically achieve. When X (t) is specifie for the virtual structure, X (t) will track X (t) accoring to the control law for the virtual structure with formation feeback. If the formation moves too fast, X X will increase. As a result of the formation feeback, the virtual structure will slow own for the spacecraft to track their esire states, that is, to keep the formation. Thus X X will ecrease corresponingly, an the formation can keep moving towar its goal with a reasonable spee. As this couple proceure procees with time, the formation maneuver will be asymptotically achieve. 5 Simulation Results In this section we will consier a group of three spacecraft each with mass given 5 Kg. The esire original positions of the three spacecraft are given by r = [8; ; ]T, r = [; 8; ]T, r = [; ; 8]T meters an the esire original attitues are given by q = q = q = [; ; ; ]T with respect to the formation frame C. We suppose that the three spacecraft start from rest with some initial errors. The three spacecraft will perform a formation maneuver of a combination of translation, rotation, an expansion. The formation will start from rest with inertial position r () = [; ; ] T an inertial attitue q () = [; ; ; ] T to the esire position r = [; ; ]T an esire attitue q = [ut sin(=); cos(=)] T, where u = [= p ; = p ; = p ] T, an expan.5 times the original size. r i r i (m) abs( r i r i+ r i+ r i+ ) (m).5.5 (a) spacecraft # spacecraft # 8 x (c) i= i= 8 q i q i+ (b) spacecraft # spacecraft # 8 8 x () i= i= 8 igure : Position an attitue errors without formation feeback (convergence time: 9 sec). In simulation, we will plot absolute position an attitue errors as well as relative position an attitue errors for each spacecraft to the time when the system converges. When X X + X X < :, we say that the system has converge. Absolute position error is represente by absolute ifference between actual position an esire position for each spacecraft. Absolute attitue error is represente by absolute ifference between actual attitue an esire attitue for each spacecraft. Since the formation shape is an equilateral triangle an the three spacecraft keep the same attitue in the formation, we use absolute ifference between lengths of the sies in the equilateral triangle to represent the relative position error an absolute ifference between the attitue of each spacecraft to represent the relative attitue error. If the formation is preserve exactly, the relative position an attitue errors shoul be zero. In this section, we use a subscript i ( i ) which is efine moulo to represent the states for the ith spacecraft. or each figure in this section, in part (a), we plot absolute position errors represente by r i ri. In part (b), we plot absolute attitue errors represente by q i qi. In part (c), we plot relative position errors represente byjkr i r i+ k kr i+ r i+ kj. In part (), we plot relative attitue errors represente by kq i q i+ k. Note that sometimes some curves may coincie with each other. igure shows the formation maneuver without formation feeback. igure shows the formation maneu-

7 r i r i (m) (a) spacecraft # spacecraft # q i (b) spacecraft # spacecraft # r i r i (m) 8 (a) spacecraft # spacecraft # q i (b) spacecraft # spacecraft # abs( r i r i+ r i+ r i+ ) (m) x (c) i= i= 5 5 q i+ 8 x () i= i= 5 5 abs( r i r i+ r i+ r i+ ) (m) 7 5 (c) i= i= 8 q i () i= i= 8 igure : Position an attitue errors with formation feeback (convergence time: 5 sec). ver with formation feeback. By comparing each part of igure an igure, we can see that the maximum absolute an relative errors of the system without formation feeback is larger than that of the one with formation feeback. Also the system without formation feeback converges faster than the one with formation feeback when we choose the same rise time for both of them. When we ecrease the rise time in igure to let the system converge faster, the corresponing errors will increase significantly. Similarly, we can also increase the rise time to ecrease the errors, but the system will converge more slowly. In igure, since the system has formation feeback, we can choose smaller rise time than that in igure to let the system converge faster while the errors are still maintaine within a reasonable range. In igure an 5, we simulate the formation maneuver results when spacecraft # fails from 5th to th secon with an without formation feeback respectively. Since there is no formation feeback in igure, the virtual structure keeps moving towar its final goal even if one of the spacecraft fails for some time. As a result, spacecraft # cannot track its esire states satisfactorily, an the system has very large absolute an relative errors uring the perio when spacecraft # fails. In fact, in this case the spacecraft are out of formation for some time. However, in igure 5, since there is formation feeback, the virtual structure slows own to preserve the formation when one of the spacecraft fails for a perio of time. As a result, the system in igure 5 has much smaller absolute an relative errors than the one in igure. The formation is preserve much better than that in igure even if spacecraft # fails for 5 secons. Within the same range of error, the system with formation feeback can choose smaller rise time, an thus igure : Position an attitue errors without formation feeback when spacecraft # fails for 5 secons (convergence time: 9 sec). converge faster than the one without formation feeback. Within the same range of convergence spee, the system with formation feeback will have smaller errors than the one without formation feeback. We know that absolute an relative errors will ecrease when the formation gain K increases, but the convergence spee will ecrease corresponingly. At the same time, when the formation gain K ecreases, the system will converge faster, but the absolute an relative errors will increase corresponingly. We also know by simulation that it is har to choose goo rise time beforehan to achieve a goo performance in the system without formation feeback. However, a wie range of rise times work well in the system with formation feeback. Conclusion In this paper we have investigate a novel iea of introucing formation feeback uner the scheme of virtual structures through a etaile application of this iea to the problem of synthesizing multiple spacecraft in eep space. Introucing formation feeback from spacecraft to the formation has several avantages. irst, the system can achieve a goo performance in improving convergence spee an ecreasing maneuver errors. Secon, formation feeback as a sense of group stability an robustness to the whole system. Thir, formation feeback improves the robustness with respect to choosing gains for ifferent spacecraft. inally, formation feeback makes formation keeping more robust to synchronization issues an the variability on each spacecraft. 7

8 r i r i (m) abs( r i r i+ r i+ r i+ ) (m) (a) spacecraft # spacecraft # (c) i= i= 5 5 q i q i (b) spacecraft # spacecraft # () i= i= 5 5 igure 5: Position an attitue errors with formation feeback when spacecraft # fails for 5 secons (convergence time: 5 sec). Acknowlegements This work was fune by Jet Propulsion Laboratory, California Institute of Technology uner grant #9-5. References [] P. K. C. Wang, Navigation strategies for multiple autonomous mobile robots moving in formation, Journal of Robotic Systems, vol. 8, no., pp , 99. Control an Dynamics, vol., pp. 8 5, January 999. [7] B. Young, R. Bear, an J. Kelsey, A control scheme for improving multi-vehicle formation maneuvers, in Proceeings of the American Control Conference, (Arlington, VA), June. [8] J. Lawton, R. Bear, an. Haaegh, Elementary attitue formation maneuvers via behavior-base control, in Proceeings of the AIAA Guiance, Navigation an Control Conference, (Denver, CO), August. AIAA Paper No. AIAA--. [9] T. Balch an R. C. Arkin, Behavior-base formation control for multirobot teams, IEEE Transactions on Robotics an Automation, vol., pp. 9 99, December 998. [] J. Lawton, B. Young, an R. Bear, A ecentralize approach to elementary formation maneuvers, in IEEE International Conference on Robotics an Automation, (San ransico, CA), April. [] R. W. Bear, J. Lawton, an. Y. Haaegh, A feeback architecture for formation control, IEEE Transactions on Control Systems Technology, vol. 9, pp , November. [] P. C. Hughes, Spacecraft Attitue Dynamics. John Wiley & Sons, 98. [] J. T.-Y. Wen an K. Kreutz-Delgao, The attitue control problem, IEEE Transactions on Automatic Control, vol., pp. 8, October 99. [] M. A. Lewis an K.-H. Tan, High precision formation control of mobile robots using virtual structures, Autonomous Robots, vol., pp. 87, 997. [] B. Wie, H. Weiss, an A. Apapostathis, Quaternion feeback regulator for spacecraft eigenaxis rotations, AIAA Journal of Guiance an Control, vol., pp. 75 8, May 989. [] P. K. C. Wang an. Y. Haaegh, Coorination an control of multiple microspacecraft moving in formation, The Journal of the Astronautical Sciences, vol., no., pp. 5 55, 99. [5]. Y. Haaegh, W.-M. Lu, an P. K. C. Wang, Aaptive control of formation flying spacecraft for interferometry, in IAC, IAC, 998. [] P. Wang,. Haaegh, an K. Lau, Synchronize formation rotation an attitue control of multiple free-flying spacecraft, AIAA Journal of Guiance, 8

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

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

Attitude Regulation About a Fixed Rotation Axis

Attitude Regulation About a Fixed Rotation Axis AIAA Journal of Guidance, Control, & Dynamics Revised Submission, December, 22 Attitude Regulation About a Fixed Rotation Axis Jonathan Lawton Raytheon Systems Inc. Tucson, Arizona 85734 Randal W. Beard

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

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

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

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

Calculus in the AP Physics C Course The Derivative

Calculus in the AP Physics C Course The Derivative Limits an Derivatives Calculus in the AP Physics C Course The Derivative In physics, the ieas of the rate change of a quantity (along with the slope of a tangent line) an the area uner a curve are essential.

More information

Short Intro to Coordinate Transformation

Short Intro to Coordinate Transformation Short Intro to Coorinate Transformation 1 A Vector A vector can basically be seen as an arrow in space pointing in a specific irection with a specific length. The following problem arises: How o we represent

More information

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

TRACKING CONTROL OF MULTIPLE MOBILE ROBOTS: A CASE STUDY OF INTER-ROBOT COLLISION-FREE PROBLEM 265 Asian Journal of Control, Vol. 4, No. 3, pp. 265-273, September 22 TRACKING CONTROL OF MULTIPLE MOBILE ROBOTS: A CASE STUDY OF INTER-ROBOT COLLISION-FREE PROBLEM Jurachart Jongusuk an Tsutomu Mita

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

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

IN the recent past, the use of vertical take-off and landing

IN the recent past, the use of vertical take-off and landing IEEE TRANSACTIONS ON ROBOTICS, VOL. 27, NO. 1, FEBRUARY 2011 129 Aaptive Position Tracking of VTOL UAVs Anrew Roberts, Stuent Member, IEEE, an Abelhami Tayebi, Senior Member, IEEE Abstract An aaptive position-tracking

More information

Chapter 2 Lagrangian Modeling

Chapter 2 Lagrangian Modeling Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie

More information

Numerical Integrator. Graphics

Numerical Integrator. Graphics 1 Introuction CS229 Dynamics Hanout The question of the week is how owe write a ynamic simulator for particles, rigi boies, or an articulate character such as a human figure?" In their SIGGRPH course notes,

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

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

Euler equations for multiple integrals

Euler equations for multiple integrals Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................

More information

Attitude Control System Design of UAV Guo Li1, a, Xiaoliang Lv2, b, Yongqing Zeng3, c

Attitude Control System Design of UAV Guo Li1, a, Xiaoliang Lv2, b, Yongqing Zeng3, c 4th National Conference on Electrical, Electronics an Computer Engineering (NCEECE 205) Attitue Control ystem Design of UAV Guo Li, a, Xiaoliang Lv2, b, Yongqing eng3, c Automation chool, University of

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

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

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

State-Space Model for a Multi-Machine System

State-Space Model for a Multi-Machine System State-Space Moel for a Multi-Machine System These notes parallel section.4 in the text. We are ealing with classically moele machines (IEEE Type.), constant impeance loas, an a network reuce to its internal

More information

Physics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1

Physics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1 Physics 5153 Classical Mechanics The Virial Theorem an The Poisson Bracket 1 Introuction In this lecture we will consier two applications of the Hamiltonian. The first, the Virial Theorem, applies to systems

More information

Gravitation as the result of the reintegration of migrated electrons and positrons to their atomic nuclei. Osvaldo Domann

Gravitation as the result of the reintegration of migrated electrons and positrons to their atomic nuclei. Osvaldo Domann Gravitation as the result of the reintegration of migrate electrons an positrons to their atomic nuclei. Osvalo Domann oomann@yahoo.com (This paper is an extract of [6] liste in section Bibliography.)

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

Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing

Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing Course Project for CDS 05 - Geometric Mechanics John M. Carson III California Institute of Technology June

More information

Momentum and Energy. Chapter Conservation Principles

Momentum and Energy. Chapter Conservation Principles Chapter 2 Momentum an Energy In this chapter we present some funamental results of continuum mechanics. The formulation is base on the principles of conservation of mass, momentum, angular momentum, an

More information

Lagrangian and Hamiltonian Mechanics

Lagrangian and Hamiltonian Mechanics Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical

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

ECE 422 Power System Operations & Planning 7 Transient Stability

ECE 422 Power System Operations & Planning 7 Transient Stability ECE 4 Power System Operations & Planning 7 Transient Stability Spring 5 Instructor: Kai Sun References Saaat s Chapter.5 ~. EPRI Tutorial s Chapter 7 Kunur s Chapter 3 Transient Stability The ability of

More information

739. Design of adaptive sliding mode control for spherical robot based on MR fluid actuator

739. Design of adaptive sliding mode control for spherical robot based on MR fluid actuator 739. Design of aaptive sliing moe control for spherical robot base on MR flui actuator M. Yue,, B. Y. Liu School of Automotive Engineering, Dalian University of echnology 604, Dalian, Liaoning province,

More information

Math 1B, lecture 8: Integration by parts

Math 1B, lecture 8: Integration by parts Math B, lecture 8: Integration by parts Nathan Pflueger 23 September 2 Introuction Integration by parts, similarly to integration by substitution, reverses a well-known technique of ifferentiation an explores

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

A New Backstepping Sliding Mode Guidance Law Considering Control Loop Dynamics

A New Backstepping Sliding Mode Guidance Law Considering Control Loop Dynamics pp. 9-6 A New Backstepping liing Moe Guiance Law Consiering Control Loop Dynamics V. Behnamgol *, A. Vali an A. Mohammai 3, an 3. Department of Control Engineering, Malek Ashtar University of Technology

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

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

Neural Network Controller for Robotic Manipulator

Neural Network Controller for Robotic Manipulator MMAE54 Robotics- Class Project Paper Neural Network Controller for Robotic Manipulator Kai Qian Department of Biomeical Engineering, Illinois Institute of echnology, Chicago, IL 666 USA. Introuction Artificial

More information

Conservation Laws. Chapter Conservation of Energy

Conservation Laws. Chapter Conservation of Energy 20 Chapter 3 Conservation Laws In orer to check the physical consistency of the above set of equations governing Maxwell-Lorentz electroynamics [(2.10) an (2.12) or (1.65) an (1.68)], we examine the action

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

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

Solutions to Practice Problems Tuesday, October 28, 2008

Solutions to Practice Problems Tuesday, October 28, 2008 Solutions to Practice Problems Tuesay, October 28, 2008 1. The graph of the function f is shown below. Figure 1: The graph of f(x) What is x 1 + f(x)? What is x 1 f(x)? An oes x 1 f(x) exist? If so, what

More information

Implicit Lyapunov control of closed quantum systems

Implicit Lyapunov control of closed quantum systems Joint 48th IEEE Conference on Decision an Control an 28th Chinese Control Conference Shanghai, P.R. China, December 16-18, 29 ThAIn1.4 Implicit Lyapunov control of close quantum systems Shouwei Zhao, Hai

More information

Linear First-Order Equations

Linear First-Order Equations 5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)

More information

A Second Time Dimension, Hidden in Plain Sight

A Second Time Dimension, Hidden in Plain Sight A Secon Time Dimension, Hien in Plain Sight Brett A Collins. In this paper I postulate the existence of a secon time imension, making five imensions, three space imensions an two time imensions. I will

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

Adaptive Optimal Path Following for High Wind Flights

Adaptive Optimal Path Following for High Wind Flights Milano (Italy) August - September, 11 Aaptive Optimal Path Following for High Win Flights Ashwini Ratnoo P.B. Sujit Mangal Kothari Postoctoral Fellow, Department of Aerospace Engineering, Technion-Israel

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

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

Free rotation of a rigid body 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012

Free rotation of a rigid body 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012 Free rotation of a rigi boy 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012 1 Introuction In this section, we escribe the motion of a rigi boy that is free to rotate

More information

Gyroscopic matrices of the right beams and the discs

Gyroscopic matrices of the right beams and the discs Titre : Matrice gyroscopique es poutres roites et es i[...] Date : 15/07/2014 Page : 1/16 Gyroscopic matrices of the right beams an the iscs Summary: This ocument presents the formulation of the matrices

More information

Lagrangian and Hamiltonian Dynamics

Lagrangian and Hamiltonian Dynamics Lagrangian an Hamiltonian Dynamics Volker Perlick (Lancaster University) Lecture 1 The Passage from Newtonian to Lagrangian Dynamics (Cockcroft Institute, 22 February 2010) Subjects covere Lecture 2: Discussion

More information

Define each term or concept.

Define each term or concept. Chapter Differentiation Course Number Section.1 The Derivative an the Tangent Line Problem Objective: In this lesson you learne how to fin the erivative of a function using the limit efinition an unerstan

More information

Automatica. Global trajectory tracking control of VTOL-UAVs without linear velocity measurements

Automatica. Global trajectory tracking control of VTOL-UAVs without linear velocity measurements Automatica 46 (200) 053 059 Contents lists available at ScienceDirect Automatica journal homepage: www.elsevier.com/locate/automatica Brief paper Global trajectory tracking control of VTOL-UAVs without

More information

Invariant Extended Kalman Filter: Theory and application to a velocity-aided estimation problem

Invariant Extended Kalman Filter: Theory and application to a velocity-aided estimation problem Invariant Extene Kalman Filter: Theory an application to a velocity-aie estimation problem S. Bonnabel (Mines ParisTech) Joint work with P. Martin (Mines ParisTech) E. Salaun (Georgia Institute of Technology)

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

Newton Euler equations in general coordinates

Newton Euler equations in general coordinates Newton Euler equations in general coorinates y ertol ongar an Frank Kirchner Robotics Innoation Center, DFKI GmbH, remen, Germany Abstract For the computation of rigi boy ynamics, the Newton Euler equations

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

Diagonalization of Matrices Dr. E. Jacobs

Diagonalization of Matrices Dr. E. Jacobs Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is

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

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

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

RECENTLY, vertical take-off and landing (VTOL) unmanned

RECENTLY, vertical take-off and landing (VTOL) unmanned IEEE/CAA JOURNAL OF AUTOMATICA SINICA VOL. NO. JANUARY 5 65 Trajectory Tracking of Vertical Take-off an Laning Unmanne Aerial Vehicles Base on Disturbance Rejection Control Lu Wang an Jianbo Su Abstract

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

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

Formation feedback control for multiple spacecraft via virtual structures

Formation feedback control for multiple spacecraft via virtual structures Formation feedback control for multiple spacecraft via virtual structures W. Ren and R.W. Beard Abstract: Formation control ideas for multiple spacecraft using the virtual structure approach are presented.

More information

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

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

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

Reactionless Path Planning Strategies for Capture of Tumbling Objects in Space Using a Dual-Arm Robotic System

Reactionless Path Planning Strategies for Capture of Tumbling Objects in Space Using a Dual-Arm Robotic System Reactionless Path Planning Strategies for Capture of Tumbling Objects in Space Using a Dual-Arm Robotic System Suril V. Shah Robotics Research Center, IIIT Hyeraba, Anhra Praesh, 53, Inia surilshah@iiit.ac.in

More information

Tutorial Test 5 2D welding robot

Tutorial Test 5 2D welding robot Tutorial Test 5 D weling robot Phys 70: Planar rigi boy ynamics The problem statement is appene at the en of the reference solution. June 19, 015 Begin: 10:00 am En: 11:30 am Duration: 90 min Solution.

More information

Left-invariant extended Kalman filter and attitude estimation

Left-invariant extended Kalman filter and attitude estimation Left-invariant extene Kalman filter an attitue estimation Silvere Bonnabel Abstract We consier a left-invariant ynamics on a Lie group. One way to efine riving an observation noises is to make them preserve

More information

Novel Coordinate Transformation and Robust Cooperative Formation Control for Swarms of Spacecraft

Novel Coordinate Transformation and Robust Cooperative Formation Control for Swarms of Spacecraft 4th International Conerence on Spacecrat Formation Flying Missions an echnologies May 18-2 211 St. Hubert Quebec Canaa Novel Coorinate ransormation an Robust Cooperative Formation Control or Swarms o Spacecrat

More information

A Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential

A Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential Avances in Applie Mathematics an Mechanics Av. Appl. Math. Mech. Vol. 1 No. 4 pp. 573-580 DOI: 10.4208/aamm.09-m0946 August 2009 A Note on Exact Solutions to Linear Differential Equations by the Matrix

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

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

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

Differentiation ( , 9.5)

Differentiation ( , 9.5) Chapter 2 Differentiation (8.1 8.3, 9.5) 2.1 Rate of Change (8.2.1 5) Recall that the equation of a straight line can be written as y = mx + c, where m is the slope or graient of the line, an c is the

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

Basic Differentiation Rules and Rates of Change. The Constant Rule

Basic Differentiation Rules and Rates of Change. The Constant Rule 460_00.q //04 4:04 PM Page 07 SECTION. Basic Differentiation Rules an Rates of Change 07 Section. The slope of a horizontal line is 0. Basic Differentiation Rules an Rates of Change Fin the erivative of

More information

Nonlinear Backstepping Control of Permanent Magnet Synchronous Motor (PMSM)

Nonlinear Backstepping Control of Permanent Magnet Synchronous Motor (PMSM) International ournal o Systems Control (Vol.1-010/Iss.1) Merzoug an Benalla / Nonlinear Backstepping Control o ermanent... / pp. 30-34 Nonlinear Backstepping Control o ermanent Magnet Synchronous Motor

More information

The proper definition of the added mass for the water entry problem

The proper definition of the added mass for the water entry problem The proper efinition of the ae mass for the water entry problem Leonaro Casetta lecasetta@ig.com.br Celso P. Pesce ceppesce@usp.br LIE&MO lui-structure Interaction an Offshore Mechanics Laboratory Mechanical

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

Guidance Control of a Small Unmanned Aerial Vehicle with a Delta Wing

Guidance Control of a Small Unmanned Aerial Vehicle with a Delta Wing Proceeings of Australasian Conference on Robotics an Automation, -4 Dec 13, University of New South Wales, Syney Australia Guiance Control of a Small Unmanne Aerial Vehicle with a Delta Wing Shuichi TAJIMA,

More information

Angles-Only Orbit Determination Copyright 2006 Michel Santos Page 1

Angles-Only Orbit Determination Copyright 2006 Michel Santos Page 1 Angles-Only Orbit Determination Copyright 6 Michel Santos Page 1 Abstract This ocument presents a re-erivation of the Gauss an Laplace Angles-Only Methos for Initial Orbit Determination. It keeps close

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

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

Multirate Feedforward Control with State Trajectory Generation based on Time Axis Reversal for Plant with Continuous Time Unstable Zeros

Multirate Feedforward Control with State Trajectory Generation based on Time Axis Reversal for Plant with Continuous Time Unstable Zeros Multirate Feeforwar Control with State Trajectory Generation base on Time Axis Reversal for with Continuous Time Unstable Zeros Wataru Ohnishi, Hiroshi Fujimoto Abstract with unstable zeros is known as

More information

arxiv: v1 [math.oc] 25 Nov 2017

arxiv: v1 [math.oc] 25 Nov 2017 Constraine Geometric Attitue Control on SO(3) Shankar Kulumani* an Taeyoung Lee November 8, 017 arxiv:17199v1 [math.oc] 5 Nov 017 Abstract This paper presents a new geometric aaptive control system with

More information

On Using Unstable Electrohydraulic Valves for Control

On Using Unstable Electrohydraulic Valves for Control Kailash Krishnaswamy Perry Y. Li Department of Mechanical Engineering, University of Minnesota, 111 Church St. SE, Minneapolis, MN 55455 e-mail: kk,pli @me.umn.eu On Using Unstable Electrohyraulic Valves

More information

The Exact Form and General Integrating Factors

The Exact Form and General Integrating Factors 7 The Exact Form an General Integrating Factors In the previous chapters, we ve seen how separable an linear ifferential equations can be solve using methos for converting them to forms that can be easily

More information

Unit quaternion observer based attitude stabilization of a rigid spacecraft without velocity measurement

Unit quaternion observer based attitude stabilization of a rigid spacecraft without velocity measurement Proceedings of the 45th IEEE Conference on Decision & Control Manchester Grand Hyatt Hotel San Diego, CA, USA, December 3-5, 6 Unit quaternion observer based attitude stabilization of a rigid spacecraft

More information

Calculus Class Notes for the Combined Calculus and Physics Course Semester I

Calculus Class Notes for the Combined Calculus and Physics Course Semester I Calculus Class Notes for the Combine Calculus an Physics Course Semester I Kelly Black December 14, 2001 Support provie by the National Science Founation - NSF-DUE-9752485 1 Section 0 2 Contents 1 Average

More information

Further Differentiation and Applications

Further Differentiation and Applications Avance Higher Notes (Unit ) Prerequisites: Inverse function property; prouct, quotient an chain rules; inflexion points. Maths Applications: Concavity; ifferentiability. Real-Worl Applications: Particle

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

Section 2.1 The Derivative and the Tangent Line Problem

Section 2.1 The Derivative and the Tangent Line Problem Chapter 2 Differentiation Course Number Section 2.1 The Derivative an the Tangent Line Problem Objective: In this lesson you learne how to fin the erivative of a function using the limit efinition an unerstan

More information

Systems & Control Letters

Systems & Control Letters Systems & ontrol Letters ( ) ontents lists available at ScienceDirect Systems & ontrol Letters journal homepage: www.elsevier.com/locate/sysconle A converse to the eterministic separation principle Jochen

More information

Angular Velocity Bounds via Light Curve Glint Duration

Angular Velocity Bounds via Light Curve Glint Duration Angular Velocity Bouns via Light Curve Glint Duration Joanna C. Hinks an John L. Crassiis The time uration of specular glint events within a light curve is relate to the angular velocity of the reflecting

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