Spacecraft Angular Velocity Stabilization Using A Single-Gimbal Variable Speed Control Moment Gyro
|
|
- Jacob Chandler
- 6 years ago
- Views:
Transcription
1 Spacecraft Angular Velocity Stabilization Using A Single-Gimbal Variable Speed Control Moment Gyro Ancil Marsall and Panagiotis Tsiotras Georgia Institute of Tecnology, Atlanta, GA Feedback controllers for te stabilization of te angular velocity vector of a rigid spacecraft using a single-gimbal Variable Speed Control Moment Gyro (VSCMG) are presented. Linearization of te equations of motion sow tat complete attitude stabilization is not possible via linear metods. Noneteless, it is sown tat te linearized angular velocity equations are controllable, and a simple LQR control law is used to locally asymptotically stabilize te angular velocity vector. A Lyapunov-based approac is subsequently used to derive a state feedback control law tat globally asymptotically stabilizes te nonlinear angular velocity system. Introduction Stabilization of te angular velocity equations of a rigid spacecraft wit less tan tree control torques ave been addressed in several papers using various tecniques. In Ref. 1 it was sown, via Lyapunov metods, tat te angular velocity equations can be made locally asymptotically stable about te origin by means of two torques applied along two principal axes. Te control law proposed in Ref. 1 was nonlinear. Reference 2 complemented tese results by sowing asymptotic stability via te construction of a center manifold. A new control law was proposed, and te control law of Ref. 1 was verified. Reference 3 continued tis avenue of researc by sowing tat one external torque, aligned wit a principal axis, could stabilize te angular velocity vector about te origin. Moreover, it was sown tat te controller was robust relative to canges in te parameters defining te In Ref. 4, global asymptotic stability of te angular velocity was proved using a single, linear control law, provided tat te spacecraft as no symmetries. It was also sown tat a single control torque aligned wit a principal axis cannot asymptotically stabilize te system. Reference 5 furter proved tat a body wit an axis of symmetry can be globally asymptotically stabilized using one control torque. Te resulting control law must necessarily be nonlinear. However, no controller was presented in Ref. 5. Reference 6 verified tat te results of Refs. 4 Graduate Student, Scool of Aerospace Engineering. Member AIAA. ancil marsall@ae.gatec.edu Associate Professor, Scool of Aerospace Engineering. Associate FellowAIAA. p.tsiotras@ae.gatec.edu, Tel: (44) , Fax: (44) Corresponding autor. Copyrigt c 23 by P. Tsiotras and A. Marsall. Publised by te American Institute of Aeronautics and Astronautics, Inc. wit permission. and 5 follow easily as an application of te Jurdjevic- Quinn approac. It also included an explicit nonlinear control law wic provided global stability for te axisymmetric case. Reference 7 sowed tat te angular velocity of an axi-symmetric rigid body can be globally asymptotically stabilized by means of a linear feedback wen two control torques act on te body. Oter approaces used to develop globally asymptotically stabilizing controllers for a rigid spacecraft wit two torques include te general metodology of nonlinear zero dynamics in Ref. 8, and te energy tecniques of Ref. 9. On te same token, te autors of Ref. 1 addressed te angular velocity stabilization of an almost axi-symmetric rigid spacecraft wit partial attitude stabilization using two external torques. In te previously mentioned references, te control torques are assumed to be provided by gas jets. Alternatively, internal torques can be generated by momentum (or reaction) weels or control moment gyroscopes (CMGs). Te spin axis of a momentum weel is fixed in te body frame, and te spin rate of te flyweel is varied to produce a torque along te spin axis. In te CMG case te weel speed of te flyweel is kept constant. A gimbal assembly canges te spin axis of te flyweel, tus producing a torque wic is ortogonal to bot te spin and gimbal axes of te CMG. It is well known tat te primary advantage of single-gimbal CMGs over oter momentum excange devices is teir torque amplification property, tat is, te output torque produced from te rotation of te angular momentum vector is muc larger tan te one required for gimbal rotation. Several references discuss te use of CMGs for spacecraft attitude control. See, for example, Refs A complete controllability analysis of te spacecraft equations as been reported in Ref. 15. Tere, it is
2 sown tat te system is not controllable wit less tan tree reaction weels. Krisnan et al 16 provided a control law using two momentum weels for te restricted case of zero angular momentum. Reference 17 developed a control law to stabilize te spin axis of a rigid spacecraft about a specified inertia axis using two reaction weels. Finally, te autors in Ref. 18 applied modern nonlinear control tecniques for detumbling of a spacecraft wit a single momentum weel aligned along one of te spacecraft principal axes. Te use of Variable Speed Control Moment Gyros (VSCMGs) for spacecraft stabilization as received 19, 2 attention recently. A VSCMG can be tougt of as a ybrid device comprised of a momentum (or reaction) weel and a CMG. In particular, te weel speed of a VSCMG is allowed to vary, tus producing an additional torque over a conventional CMG. Tis torque is not fixed in te spacecraft body frame, as in te case of a momentum weel; rater, te direction of te spin axis of te VSCMG is allowed to rotate via a gimbal. An additional torque, perpendicular to te spin and gimbal axes is tus generated, as in te conventional CMG case. Tis additional degree of freedom can be utilized to avoid te gimbal lock singularity tat as plagued traditional CMG clusters. 19 Te use of te additional torque of VSCMGs as also been utilized for attitude control (and energy storage) of spacecraft in Refs In bot te CMG and VSCMG cases presented in te literature to date, a cluster of actuators as been used to provide a sufficient number of torques to acieve complete attitude stabilization (and possibly energy storage). In tis paper, we consider te case of control (stabilization) of a spacecraft via a single VSCMG actuator. Te outline of tis paper is as follows. First, we present te complete equations of motion of a spacecraft wit one VSCMG in an arbitrary orientation. Tese equations are composed of te dynamic and kinematic equations. Next, we linearize te equations of motion about an equilibrium point. Linearization sows tat te spacecraft attitude is uncontrollable wit only one VSCMG. However, te angular velocity equations remain controllable. A simple LQR feedback law is designed to acieve local asymptotic stability at te origin of te linearized angular velocity system. Next, we examine te exact angular velocity equations. We derive a nonlinear control law tat ensures global asymptotic stability of te angular velocity of te spacecraft about te origin using only one VSCMG. Several numerical examples are included to demonstrate te viability of te control algoritms proposed. Equations of Motion Te dynamic equations of motion of a spacecraft wit a cluster of VSCMGs ave been fully derived in te literature Herein, we will use te equations as derived by Ricie et al 21 and Yoon and Tsiotras. 2 In Ref. 2 it is assumed tat te center of mass of eac VSCMG weel coincides wit tat of te gimbal structure; te spacecraft, weels, and gimbal structure are rigid; te flyweels and gimbals are balanced; and te spacecraft rotational motion is decoupled from its translational motion. Figure 1 sows a spacecraft wit a single VSCMG. Te origin of te body frame B, is located at te center of mass of te entire spacecraft. Te gimbal frame G, represented by te ortonormal set of unit vectors, and, is arbitrarily located in te spacecraft platform. Fig. 1 Body Frame O Spacecraft wit a single VSCMG Dynamics in te body frame Specializing te dynamical equations of motion presented in Yoon and Tsiotras 2 to a single VSCMG, we get J ω+ Jω+A g I cg γ +A t Ω γ +A s Ω+ω = (1) were, := Jω + A g I cg γ + A s Ω (2) J(γ) := I B B + A s (γ)i cs A T s (γ)+a t (γ)i ct A T t (γ) +A g I cg A T g (3) Te argument in J( ) is included to denote explicitly te dependence of te spacecraft inertia matrix on te gimbal axis angle, γ. For notational simplicity, in te sequel we will often drop te argument wen it is clear from te context. As a result of (3) it follows tat J = γa t (I cs I ct )A T s + γa s (I cs I ct )A T t (4) Notice tat J is linear in te gimbal rate. In (2) te (column) vector IR 3,istetotal angular momentum of te spacecraft wit respect to te inertial
3 frame, expressed in te body frame. Similarly, ω IR 3 is te angular velocity vector of te spacecraft wit respect to te inertial frame, expressed in te body frame of te spacecraft. Te quantities Ω, and Ω are te magnitudes of te angular velocity and angular acceleration of te VSCMG weel, respectively, about te gimbal spin axis, wit respect to te gimbal frame. Also, γ and γ are te magnitudes of te gimbal rate and gimbal acceleration, respectively, about te gimbal axis, wit respect to te spacecraft platform. In our analysis we will assume, as usual, gimbal rate commands. Tis is also te case for standard CMG s in order to take full advantage of te torque amplification property. Tis implies a velocity steering law for te gimbal. Tus, our control inputs are Ω and γ and we can write γ = u 1, Ω =u2 (5) Often we use u := [u 1, u 2 ] T IR 2 for te combined control vector. Te matrix-valued function J : [, 2π) IR 3 3 provides te inertia matrix of te entire spacecraft, i.e. te spacecraft platform, including te weel and gimbal structure of te VSCMG, given te gimbal angle γ. Notice tat te matrix J(γ) is positive definite for all values of te gimbal angle γ [, 2π). Te constant matrix I B B is te sum of te inertia of te spacecraft platform, and te inertias of te mass centers of te weel, gimbal and spacecraft platform, about te equivalent mass center of te entire spacecraft. I c represents te sum of te inertia scalars of te weel and gimbal structure, i.e. I c = I w + I g, were = {s, t, g}. Te column vectors A IR 3, were = {s, t, g} are te body frame representations of te gimbal frame unit vectors,, and. Observe tat A s and A t are functions of te gimbal angle γ, as follows A s (γ) = A so c γ + A to s γ (6a) A t (γ) = A so s γ + A to c γ (6b) were A so and A to are te values of A s and A t at some initial time, and c γ := cos γ and s γ := sin γ. Finally, for any vector v =[v 1, v 2 v 3 ] T IR 3,te notation v IR 3 3 represents te skew symmetric matrix v 3 v 2 v = v 3 v 1. v 2 v 1 For te details of te derivation of equations (1)-(3), as well as te notation used in tis paper, te interested reader may refer to Ref. 2. Kinematics Witout loss of generality, Euler angles will be used to represent te attitude of te spacecraft. For a Euler angle sequence, te kinematic equations are given by 22 φ = ω 1 + ω 2 sin φ tan θ + ω 3 cos φ tan θ (7a) θ = ω 2 cos φ ω 3 sin φ (7b) ψ = ω 2 sin φ sec θ + ω 3 cos φ sec θ (7c) Linear System Analysis In tis section, we linearize te full nonlinear equations of motion, given by (1), (5) and (7), and examine teir controllability properties. We also present a linear control law wic stabilizes te angular velocity of te linearized system. Te only mild assumption made ere is tat te gimbal acceleration γ is negligible. Linearization Te equilibrium points of Eqs. (1), (5) and (7) are given by ω = φ = θ = ψ = and γ = γ o, Ω = Ω o, were γ o and Ω o are arbitrary constants. From equations (6) we get were A s A sf (γ o )+A tf (γ o ) γ (8a) A t A tf (γ o ) A sf (γ o ) γ, (8b) A sf (γ o ) := A so c γo + A to s γo, (9a) A tf (γ o ) := A so s γo + A to c γo, (9b) and were ( ) represents a small perturbation in te variable from its equilibrium value. Similarly, eac term of Eq. (1) results, to first order, in te following terms J ω J f ω (1a) Jω (1b) A t Ω γ Ω o A tf γ (1c) A s Ω Iws A sf Ω (1d) ω Ω o A sf ω (1e) were J f (γ o ):=I B B +A sf I cs A T sf +A tfi ct A T tf +A gi cg A T g. Te linearization of (1) tus yields ω = A 1 ω + B 1 γ + B 2 Ω (11) were te matrices A 1 IR 3 3, B 1 IR 3 1,andB 2 IR 3 1 are given by A 1 (γ o, Ω o ) := J 1 f Ω o A sf (12a) B 1 (γ o, Ω o ) := J 1 f Ω o A tf (12b) B 2 (γ o, Ω o ) := J 1 f A sf (12c) Note tat tese matrices depend on te equilibrium/reference values γ o and Ω o. Defining te new
4 state variable as x := [ ω T, φ, θ, ψ] T IR 6 and te control as u := [ γ, Ω] T IR 2,wecan express te linearized equations in te familiar form ẋ = A x + B u (13) were te matrices A IR 6 6 and B IR 6 2 are given by [ ] [ ] A1 A := 3 3 B1 B, B := 2 (14) were 1 is te identity matrix. Controllability of Linearized System Here we give two results on te controllability of te linearized complete system in Eqs. (13) and of te linearized angular velocity equations in Eq. (16). Proposition 1 Te linearized system described by Eqs. (13) and (14) is uncontrollable for any γ o [, 2π) and Ω o IR. Te result follows by sowing tat te controllability matrix 23 C o := [B ABA 2 BA 3 BA 4 BA 5 B] (15) as rank five for all γ o [, 2π) and Ω o IR. Since te state dimension is six, te pair (A, B) is uncontrollable. Tis result implies tat it is not possible to use linear tecniques to stabilize te complete attitude of te spacecraft using a single VSCMG. It leaves, owever, open te possibility tat te complete system of equations are controllable in te nonlinear sense. In te sequel we restrict our attention to te angular velocity subsystem. To tis end, define à := A 1 and B := [B 1 B 2 ], and rewrite Eq. (11) as ω = à ω + B u (16) Linearizing te kinematic equations, we get φ = ω 1, θ = ω 2, ψ = ω 3 (17) Proposition 2 Te linearized angular velocity system of Eq. (16) is controllable for all γ o [, 2π), and Ω o. Proof. From Teorem 3.1 of Ref. 23, (Ã, B) iscontrollable if te matrix [à λ1 B] as full row rank for all λ IR. In particular, tis must be true for all te eigenvalues λ i (i =1, 2, 3) of te matrix Ã. Calculating te row-reduced ecelon form of Ã, one obtains tat for λ = λ i (i =1, 2, 3) and Ω o, 1 [à λ1 B] = 1 (18) 1 Te rank of te matrix C o was calculated using te Symbolic Toolbox of MATLAB. 24 were ( ) are algebraic expressions in terms of te components of à tat do not affect te row rank of te matrix. Since te rank of te above matrix is always tree, te linearized angular velocity system of Eq. (11) is controllable for all γ o [, 2π), and Ω o. LQR Controller for te Angular Velocity Subsystem Given te controllable system in Eq. (16), we can find a linear control law via LQR metods. For example, we can determine a static full-state feedback law u = K ω suc tat te performance cost J = [ ω T Q ω + u T R u]dt (19) is minimized subject to te dynamics (16). Te matrix Q must be positive semi-definite wile R must be positive definite. Te control gain matrix K is given by K = R 1 BT P. Te matrix P = P T is positive semi-definite and satisfies te Algebraic Riccati Equation ÃT P + P à P BR 1 BT P + Q =. LQR optimal control designs is by now folklore. Details of can be found, for instance, in Ref. 23. Nonlinear System Analysis Te LQR controller of te previous section ensures asymptotic stability only locally about te equilibrium ω = and for gimbal angles and weel speeds close to teir reference values γ o and Ω o, respectively. Te last restriction is particularly troublesome, since stabilization of te angular velocity vector sould not inge upon γ and Ω being close to γ o and Ω o. As a matter of fact, significant control autority may tend to produce large deviations of te gimbal angle and weel speed from teir reference values; see Eqs. (5). In realistic cases, it is not reasonable to expect tat te states γ and Ω (wose values are of no particular interest, tus are not penalized in (19)) will remain small. For a more compreensive analysis of te stabilization problem, it is terefore necessary to work wit te exact, nonlinear equations of motion. In te sequel we improve on te previous results by finding a control law tat ensures global asymptotic stability for te nonlinear system. We tus also avoid te issue of te restricted (local) validity of te linearized equations due to potentially large deviation of te gimbal angle and weel speed from teir reference values. As in te linear case, in te sequel we assume tat te gimbal acceleration γ is negligible, and te control inputs are γ and Ω. Te dynamic equations of motion are tus given by J ω + Jω + A t Ω γ + A s Ω+ω = (2)
5 were and J are as in (2) and (3). To derive a stabilizing control law for (2), we consider te positive definite, continuously differentiable Lyapunov function V (ω) := 1 2 ωt Jω. Te derivative of V along te trajectories of te system is V (ω) =ω T J ω ωt Jω = ω T J ω + ω T Jω 1 2 ωt Jω = ω T (J ω + Jω) 1 2 ωt Jω = ω T (A t Ω γ + A s Ω+ω ) 1 2 ωt Jω Rewriting J = Φ γ were Φ := A t (I cs I ct )A T s + A s (I cs I ct )A T t and using te fact tat ω T ω = yields V (ω) = ω T A t Ω γ 1 2 ωt Φω γ ω T A s Ω = ω t Ω γ ω s ω t (I cs I ct ) γ ω s Ω were ω s = ω T A s and ω t = ω T A t are te components of te body angular velocity vector ω along te spin and transverse axes of te gimbal frame, respectively, i.e. ω = ω s + ω t + ω g. Proposition 3 Consider te following control law ( ( )) Ics I ct γ = k 1 (Ω) ω t Ω ω s ( ) 2 Ics I ct Ω =k 2 ω s + k 3 ω t k 4 (Ω) + k 1 (Ω) ωt 2 ω s were k 1 : R IR + is any function suc tat k 1 (Ω)Ω 2 is bounded for all Ω IR, k 4 (Ω) := Ωk (Ω), k 2 > and 2 k 2 >k 3. Tis control law globally asymptotically stabilizes te system given by Eq. (2) for all Ω(). Proof. Substituting tis control law in te expression for V (ω) leads to ( ( )) V (ω) = ωt 2 Ics I ct Ω k 1 (Ω) Ω ω s ( ( )) ω s ωt 2 Ics I ct (I cs I ct )k 1 (Ω) Ω ω s ω s (k 2 ω s + k 3 ω t k 4 (Ω) ( ) 2 Ics I ) ct + k 1 (Ω) ωt 2 ω s = ωt 2 Ω 2 k 1 (Ω) k 2 ωs 2 k 3 ω s ω t k 4 (Ω) [ ] [ ] ωt = ωt ω s G(Ω,ω,γ) ω s were te matrix G(Ω,ω,γ) is given by [ G(Ω,ω,γ):= ] 2 sgn(ω t )k 4 (Ω) k 3 Ω 2 k 1 (Ω) k 3 2 sgn(ω t )k 4 (Ω) k 2 It can be easily sown tat G(Ω,ω,γ) for all (Ω,ω,γ) R 2 [, 2π) andg(ω,ω,γ) > forω. It follows tat V. Te last inequality sows tat V, and ence ω is bounded. Terefore, γ and Ω as well as γω are bounded. Moreover, ω is bounded from (2). It follows tat ω, γ and Ω are uniformly continuous and tus V is uniformly continuous as well. From Barbalat s Lemma 25 it follows tat V. Tis implies tat ω s and Ωω t ast. Assume now tat Ω() and tat ω t ω t. Since ω s we ten ave tat after a sufficiently long time, Ω k 3 ω t Ωk (Ω) and te equilibrium Ω = is unstable. Hence, Ωω t, a contradiction. It follows tat, necessarily, ω t ast. Assume now tat ω s =Ωω t. It follows tat γ = Ω = and from Eq. (2) J ω + ω (Jω + A s Ω) = (21) wic, wen expressed in te gimbal frame, becomes J 13 ω g J 23 ωg 2 = (22a) J 23 ω g + J 13 ωg 2 + Ωω g = (22b) J 33 ω g = (22c) From Eq. (22c), we get ω g =. From Eq. (22a) or Eq. (22b), we conclude tat ω g =. Tus, te largest invariant set in {ω : V (ω) = } is te set ω =. Asymptotic stability follows from LaSalle s teorem and global asymptotic stability follows from te radial unboundedness of te function V and te fact tat te previous analysis olds for all initial conditions ω IR 3. Acceleration Steering Law In te actual spacecraft te gimbal control input is a torque (or gimbal acceleration) command, rater tan a gimbal velocity command. Te derived velocity command as to be implemented via an internal servo control loop. A simple implementation of tis idea is to use, say γ = K p ( γ d γ) (23) were K p > and were γ d as in Proposition 3. Tis (proportional) control law will ensure tat te actual gimbal velocity γ approaces te desired command γ d, as t. In practice K p as to be sufficiently large in order for te convergence to take place in a sort interval of time.
6 Numerical Examples In tis section we give some illustrative examples of te control design metods for te angular velocity subsystem using bot te linear and te nonlinear analysis of te previous sections. Bot examples applied te control laws developed earlier to te complete equations of motion in Eqs. (1)-(4) using te acceleration steering law in Eq. (23). Tis was done in order to compare eac control law individually and in relation to eac oter using a realistic evaluation model. Table 1 summarizes te values te moments of inertia and gimbal used in te simulations. Tese values rougly correspond to te spacecraft simulator described in Ref. 26. Te controller gains and te initial conditions are given in Table 2. Te first example corresponds to te LQR control design metod wic was developed from te linearized angular velocity system. Te results are sown in Fig. 2. Te weigting matrices Q and R in (19) were cosen by trial and error to stabilize te system quickly wit suitable damping. Teir values are sown in Table 2. Figures 3 and 4 sow te values of te gimbal angle and weel speed as well as teir rates. Table 1 Moments of inertia values. Symbol Value Units I B B kg m kg m 2 I wt,i wg.24 kg m 2 I gs.93 kg m 2 I gt,i gg.54 kg m 2 A so [ 1,, ] T A to [,.8161,.5779] T A go [,.5779,.8161] T Te second example corresponds to te nonlinear control law of Proposition 3. Te function k 1 (Ω) was cosen as k 1 (Ω) = µ/(1 + Ω 2 ). Te angular velocity istories wit te nonlinear control law are sown in Fig. 5. Te time istory of γ and Ω aresownin Fig. 6. Te time istories of te gimbal angles and te weel speed velocity are sown in Fig. 7. Notice tat te nonlinear controller is more aggressive resulting in larger values for te gimbal angle and weel speed. Since in a pysical system te weel and te gimbal rate commands saturate, it is imperative to modify te nonlinear control law so as to take into account tese saturation effects. Tis is left for future investigation. On te oter and, one may coose to use te nonlinear controller only if te initial conditions be- ω [rad/sec] Time istory of Angular Velocities ω 1 ω 2 ω 3.25 Fig. 2 Numerical simulations wit te LQR control law. γ [deg/sec] Ω [rpm/sec] Time istory of γ and Ω 25 Fig. 3 Time istory of γ and Ω wit te LQR Table 2 Controller gains and initial conditions. Symbol Linear Nonlinear Units ω() [.2,.1,.2] T [.1,.1,.1] T rad/sec γ() γ o = 2 2 deg γ() deg Ω() Ω o = rpm Q diag{1e 4, 1e 4, 1e 4 } R diag{1e 2, 1} µ 8 sec 1 k 2 4 sec 1 k 3 1 sec 1 2 K p 1 1 sec 1 come too large. After te trajectories reac a small neigborood of te origin (and witin te region of attraction of te linear controller), one can ten switc to te LQR controller, wose local performance can be pre-assigned via te optimization criterion (19). In
7 5 Time istory of γ and Ω 15 Time istory of γ and Ω 4 1 γ [deg] γ [deg/sec] Ω[rpm] Fig. 4 Time istory of γ and Ω wit te LQR ω [rad/sec] Time istory of Angular Velocities ω 1 ω 2 ω 3.2 Fig. 5 Numerical simulations wit te nonlinear tis sense, te linear controller acieves (local) performance and stability, wereas te nonlinear controller acts as a safety net to protect te system from large initial conditions. For comparison, in Fig. 8 we sow te results from te numerical simulations of te LQR wit initial conditions ω() = [.1.1.1] T. For tese (large) initial conditions, te LQR does not stabilize te nonlinear system. Finally, Fig. 9 sows a series of snapsots of a spacecraft wit one VSCMG undergoing a detumbling maneuver using te nonlinear control law of Proposition 3. Note tat, as expected, te final orientation of te spacecraft is suc tat te spin axis of te VSCMG is aligned wit te total angular momentum vector, wic remains constant in inertial frame at all times. Conclusions In tis paper, we ave addressed te stabilization problem of a rigid spacecraft wit a single-gimbal variable-speed control moment gyro (VSCMG). Since Ω [rpm/sec] Fig. 6 Time istory of γ and Ω wit te nonlinear γ [deg] Ω[rpm] Time istory of γ and Ω 15 Fig. 7 Time istory of γ and Ω wit te nonlinear no external control torques act on te system, reorientation of te spacecraft is acieved via momentum transfer between te spacecraft platform and te VSCMG. We sowed tat te complete attitude equations are not linearly controllable. Te angular velocity equations are, noneteless controllable. A simple LQR controller was used to locally asymptotically stabilize te angular velocity equations for an arbitrary gimbal frame orientation. Abandoning te restrictive assumptions made in te linear case we developed a control law for te nonlinear system wic ensures global asymptotic stability of te angular velocity equations. Acknowledgement: Support for tis work as been provided by AFOSR troug award no. F References 1 Brockett, R., Asymptotic Stability and Feedback Stabilization, Differential Geometric Control Teory, 1983,
8 ω ω ω ω Time =. sec, γ =. deg Time = 5. sec, γ = 2.4 deg Time = 1. sec, γ = 2. deg Time = 15. sec, γ = 31.5 deg Time = 2. sec, γ = 39.4 deg Time = 3. sec, γ = 49.1 deg Time = 4. sec, γ = 54.2 deg Time = 1. sec, γ = 57.9 deg Fig. 9 Series of snapsots of a detumbling maneuver of a spacecraft wit one VSMG. Te gimbal frame unit vectors,,and, te angular momentum vector, and te angular velocity vector ω are sown (not drawn to scale). Te time of eac snapsot and te gimbal angle γ are also depicted. y Time istory of γ and Ω.2 Fig. 8 Response of LQR controller wit large initial conditions. pp , (R.W. Brockett, R.S. Millman, and H.J. Sussman, eds.). 2 Aeyels, D., Stabilization of a class of nonlinear systems by a smoot feedback control, Systems & Control Letters, Vol. 5, 1985, pp Aeyels, D., Stabilization by smoot feedback of te angular velocity of a rigid body, Systems & Control Letters, Vol. 5, 1985, pp Aeyels, D. and Szarfanski, M., Comments on te stability of te angular velocity of a rigid body, Systems and Control Letters, Vol. 1, 1988, pp Sontag, E. and Sussmann, H., Furter comments on te stabilizability of te angular velocity of a rigid body, Systems & Control Letters, Vol. 12, 1988, pp w1 w2 w3 6 Outbib, R. and Sallet, G., Stabilizability of te angular velocity of a rigid body revisited, Systems & Control Letters, Vol. 18, 1992, pp Andriano, V., Global feedback stabilization of te angular velocity of a symmetric rigid body, Systems & Control Letters, Vol. 2, 1993, pp Brynes, C. and Isidori, A., NewResults and examples in nonlinear feedback stabilization, Systems & Control Letters, Vol. 12, 1989, pp Astolfi, A. and Ortega, R., Energy-Based Stabilization of Angular Velocity of Rigid Body in Failure Configuration, AIAA Journal of Guidance, Control and Dynamics, Vol.25, No. 1, 22, pp Tsiotras, P. and Scleicer, A., Detumbling and Partial Attitude Stabilization of a Rigid Spacecraft Under Actuator Failure, AIAA Guidance, Navigation and Control Conference, Denver, CO, 2, AIAA Paper O, S. and Vadali, S., Feedback Control and Steering Laws for Spacecraft Using Single Gimbal Control Moment Gyros, Journal of te Astronautical Sciences, Vol. 39, No. 2, 1991, pp Sing, S. N. and Bossart, T., Exact Feedback Linearization and Control of Space Station Using CMG, IEEE Transactions on Automatic Control, Vol. 38, No. 1, 1993, pp Ford, K. and Hall, C., Singular Direction Avoidance Steering for Control-Moment Gyros, Journal of Guidance, Control and Dynamics, Vol. 23, No. 4, 2, pp Heiberg, C. and Bailey, D., Precision Spacecraft Pointing Using Single-Gimbal Control Moment Gyroscopes wit Disturbance, Journal of Guidance, Control, and Dynamcis, Vol. 23, No. 1, 2, pp Crouc, P. E., Spacecraft attitude control and stabilization: Applications of geometric control teory to rigid body models, IEEE Transactions on Automatic Control, Vol. 29, No. 4, 1984, pp
9 16 Krisnan, H., McClamroc, H., and Reyanoglu, M., Attitude Stabilization of a Rigid Spacecraft Using Two Momentum Weel Actuators, Journal of Guidance, Control, and Dynamics, Vol. 18, No. 2, 1995, pp Kim, S. and Kim, Y., Spin-Axis Stabilization of a Rigid Spacecraft Using Two Reaction Weels, Journal of Guidance, Control, and Dynamics, Vol. 24, No. 5, 21, pp Bang, H., Myung, H.-S., and Tak, M.-J., Nonlinear Momentum Transfer Control of Spacecraft by Feedback Linearization, Journal of Spacecraft and Rockets, Vol. 39, No. 6, Scaub, H., Vadali, S. R., and Junkins, J. L., Feedback Control Lawfor Variable Speed Control Moment Gyro, Journal of te Astronautical Sciences, Vol. 46, No. 3, 1998, pp Yoon, H. and Tsiotras, P., Spacecraft Adaptive Attitude Control and Power Tracking Variable Speed Control Moment Gyros, Journal of Guidance, Control and Dynamics, Vol. 25, No. 6, 22, pp Ricie, D., Tsiotras, P., and Fausz, J., Simultaneous Attitude Control and Energy Storage using VSCMGs: Teory and Simulation, Proceedings of te American Control Conference, 21, pp , Arlington, VA. 22 Huges, P. C., Spacecraft Attitude Dynamics, Jon Wiley & Sons, NewYork, 1986, pp Zou, K. and Doyle, J., Essentials of Robust Control, Prentice Hall, 1998, pp Matworks, MATLAB: Symbolic Mat Toolbox v 2.1.2, Natwick, MA, Slotine, J. and Li, W., Applied Nonlinear Control, Prentice-Hall Inc., Englewood Cliffs, 1991, pp Jung, D. and Tsiotras, P., A 3-DoF Experimental Test- Bed for Integrated Attitude Dynamics and Control Researc, AIAA Guidance, Navigation and Control Conference, Austin, TX, 23, AIAA Paper
Stabilization of Angular Velocity of Asymmetrical Rigid Body. Using Two Constant Torques
Stabilization of Angular Velocity of Asymmetrical Rigid Body Using Two Constant Torques Hirohisa Kojima Associate Professor Department of Aerospace Engineering Tokyo Metropolitan University 6-6, Asahigaoka,
More informationMixed Control Moment Gyro and Momentum Wheel Attitude Control Strategies
AAS03-558 Mixed Control Moment Gyro and Momentum Wheel Attitude Control Strategies C. Eugene Skelton II and Christopher D. Hall Department of Aerospace & Ocean Engineering Virginia Polytechnic Institute
More informationTracking Rigid Body Motion Using Thrusters and Momentum. Wheels
JAS 199 Tracking Rigid Body Motion Using Thrusters and Momentum Wheels Christopher D. Hall, Panagiotis Tsiotras, and Haijun Shen Abstract Tracking control laws are developed for a rigid spacecraft using
More informationAdaptive Spacecraft Attitude Tracking Control with Actuator Uncertainties
Adaptive Spacecraft Attitude Tracking Control with Actuator Uncertainties Hyungjoo Yoon and Panagiotis Tsiotras Georgia Institute of Technology Atlanta, Georgia 3332-15, USA An adaptive control algorithm
More informationAttitude 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 informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc
More informationSpacecraft Attitude Control with RWs via LPV Control Theory: Comparison of Two Different Methods in One Framework
Trans. JSASS Aerospace Tech. Japan Vol. 4, No. ists3, pp. Pd_5-Pd_, 6 Spacecraft Attitude Control with RWs via LPV Control Theory: Comparison of Two Different Methods in One Framework y Takahiro SASAKI,),
More informationSpacecraft Attitude and Power Control Using Variable Speed Control Moment Gyros. Hyungjoo Yoon
Spacecraft Attitude and Power Control Using Variable Speed Control Moment Gyros A Thesis Presented to The Academic Faculty by Hyungjoo Yoon In Partial Fulfillment of the Requirements for the Degree Doctor
More informationClosed Loop Near Time Optimal Magnetic Attitude Control Using Dynamic Weighted Neural Network
Closed Loop Near ime Optimal Magnetic Attitude Control Using Dynamic Weigted Neural Network Ali Heydari, Seid H. Pourtakdoust Abstract e problem of time optimal magnetic attitude control is treated and
More informationRobust Leader-follower Formation Control of Mobile Robots Based on a Second Order Kinematics Model
Vol. 33, No. 9 ACTA AUTOMATICA SINICA September, 27 Robust Leader-follower Formation Control of Mobile Robots Based on a Second Order Kinematics Model LIU Si-Cai 1, 2 TAN Da-Long 1 LIU Guang-Jun 3 Abstract
More informationQuaternion-Based Tracking Control Law Design For Tracking Mode
A. M. Elbeltagy Egyptian Armed forces Conference on small satellites. 2016 Logan, Utah, USA Paper objectives Introduction Presentation Agenda Spacecraft combined nonlinear model Proposed RW nonlinear attitude
More informationContinuity and Differentiability of the Trigonometric Functions
[Te basis for te following work will be te definition of te trigonometric functions as ratios of te sides of a triangle inscribed in a circle; in particular, te sine of an angle will be defined to be te
More information232 Calculus and Structures
3 Calculus and Structures CHAPTER 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS FOR EVALUATING BEAMS Calculus and Structures 33 Copyrigt Capter 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS 17.1 THE
More informationVolume 29, Issue 3. Existence of competitive equilibrium in economies with multi-member households
Volume 29, Issue 3 Existence of competitive equilibrium in economies wit multi-member ouseolds Noriisa Sato Graduate Scool of Economics, Waseda University Abstract Tis paper focuses on te existence of
More informationIntroduction to Derivatives
Introduction to Derivatives 5-Minute Review: Instantaneous Rates and Tangent Slope Recall te analogy tat we developed earlier First we saw tat te secant slope of te line troug te two points (a, f (a))
More informationLinear Feedback Control Using Quasi Velocities
Linear Feedback Control Using Quasi Velocities Andrew J Sinclair Auburn University, Auburn, Alabama 36849 John E Hurtado and John L Junkins Texas A&M University, College Station, Texas 77843 A novel approach
More information3D Pendulum Experimental Setup for Earth-based Testing of the Attitude Dynamics of an Orbiting Spacecraft
3D Pendulum Experimental Setup for Earth-based Testing of the Attitude Dynamics of an Orbiting Spacecraft Mario A. Santillo, Nalin A. Chaturvedi, N. Harris McClamroch, Dennis S. Bernstein Department of
More informationNumerical Differentiation
Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function
More informationNonlinear Analysis and Control of a Reaction Wheel-based 3D Inverted Pendulum
Nonlinear Analysis and Control of a Reaction Weel-based 3D Inverted Pendulum Micael Muelebac, Gajamoan Moanaraja, and Raffaello D Andrea Abstract Tis paper presents te nonlinear analysis and control design
More informationPre-Calculus Review Preemptive Strike
Pre-Calculus Review Preemptive Strike Attaced are some notes and one assignment wit tree parts. Tese are due on te day tat we start te pre-calculus review. I strongly suggest reading troug te notes torougly
More informationch (for some fixed positive number c) reaching c
GSTF Journal of Matematics Statistics and Operations Researc (JMSOR) Vol. No. September 05 DOI 0.60/s4086-05-000-z Nonlinear Piecewise-defined Difference Equations wit Reciprocal and Cubic Terms Ramadan
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximating a function f(x, wose values at a set of distinct points x, x, x 2,,x n are known, by a polynomial P (x
More informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Matematics, Vol.4, No.3, 6, 4 44. ACCELERATION METHODS OF NONLINEAR ITERATION FOR NONLINEAR PARABOLIC EQUATIONS Guang-wei Yuan Xu-deng Hang Laboratory of Computational Pysics,
More informationThe Verlet Algorithm for Molecular Dynamics Simulations
Cemistry 380.37 Fall 2015 Dr. Jean M. Standard November 9, 2015 Te Verlet Algoritm for Molecular Dynamics Simulations Equations of motion For a many-body system consisting of N particles, Newton's classical
More informationPre-lab Quiz/PHYS 224 Earth s Magnetic Field. Your name Lab section
Pre-lab Quiz/PHYS 4 Eart s Magnetic Field Your name Lab section 1. Wat do you investigate in tis lab?. For a pair of Helmoltz coils described in tis manual and sown in Figure, r=.15 m, N=13, I =.4 A, wat
More informationQuantum Mechanics Chapter 1.5: An illustration using measurements of particle spin.
I Introduction. Quantum Mecanics Capter.5: An illustration using measurements of particle spin. Quantum mecanics is a teory of pysics tat as been very successful in explaining and predicting many pysical
More informationMath 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006
Mat 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006 f(x+) f(x) 10 1. For f(x) = x 2 + 2x 5, find ))))))))) and simplify completely. NOTE: **f(x+) is NOT f(x)+! f(x+) f(x) (x+) 2 + 2(x+) 5 ( x 2
More informationNew Streamfunction Approach for Magnetohydrodynamics
New Streamfunction Approac for Magnetoydrodynamics Kab Seo Kang Brooaven National Laboratory, Computational Science Center, Building 63, Room, Upton NY 973, USA. sang@bnl.gov Summary. We apply te finite
More informationDerivation Of The Schwarzschild Radius Without General Relativity
Derivation Of Te Scwarzscild Radius Witout General Relativity In tis paper I present an alternative metod of deriving te Scwarzscild radius of a black ole. Te metod uses tree of te Planck units formulas:
More informationQuasiperiodic phenomena in the Van der Pol - Mathieu equation
Quasiperiodic penomena in te Van der Pol - Matieu equation F. Veerman and F. Verulst Department of Matematics, Utrect University P.O. Box 80.010, 3508 TA Utrect Te Neterlands April 8, 009 Abstract Te Van
More informationMore on generalized inverses of partitioned matrices with Banachiewicz-Schur forms
More on generalized inverses of partitioned matrices wit anaciewicz-scur forms Yongge Tian a,, Yosio Takane b a Cina Economics and Management cademy, Central University of Finance and Economics, eijing,
More informationPerturbed Feedback Linearization of Attitude Dynamics
008 American Control Conference Westin Seattle Hotel, Seattle, Washington, USA June -3, 008 FrC6.5 Perturbed Feedback Linearization of Attitude Dynamics Abdulrahman H. Bajodah* Abstract The paper introduces
More information1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point
MA00 Capter 6 Calculus and Basic Linear Algebra I Limits, Continuity and Differentiability Te concept of its (p.7 p.9, p.4 p.49, p.55 p.56). Limits Consider te function determined by te formula f Note
More informationGlobal Output Feedback Stabilization of a Class of Upper-Triangular Nonlinear Systems
9 American Control Conference Hyatt Regency Riverfront St Louis MO USA June - 9 FrA4 Global Output Feedbac Stabilization of a Class of Upper-Triangular Nonlinear Systems Cunjiang Qian Abstract Tis paper
More informationWEIGHTING MATRICES DETERMINATION USING POLE PLACEMENT FOR TRACKING MANEUVERS
U.P.B. Sci. Bull., Series D, Vol. 75, Iss. 2, 2013 ISSN 1454-2358 WEIGHTING MATRICES DETERMINATION USING POLE PLACEMENT FOR TRACKING MANEUVERS Raluca M. STEFANESCU 1, Claudiu L. PRIOROC 2, Adrian M. STOICA
More informationChapter 5 FINITE DIFFERENCE METHOD (FDM)
MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential
More informationERROR BOUNDS FOR FINITE-DIFFERENCE METHODS FOR RUDIN OSHER FATEMI IMAGE SMOOTHING
ERROR BOUNDS FOR FINITE-DIFFERENCE METHODS FOR RUDIN OSHER FATEMI IMAGE SMOOTHING JINGYUE WANG AND BRADLEY J. LUCIER Abstract. We bound te difference between te solution to te continuous Rudin Oser Fatemi
More informationA SHORT INTRODUCTION TO BANACH LATTICES AND
CHAPTER A SHORT INTRODUCTION TO BANACH LATTICES AND POSITIVE OPERATORS In tis capter we give a brief introduction to Banac lattices and positive operators. Most results of tis capter can be found, e.g.,
More informationCopyright c 2008 Kevin Long
Lecture 4 Numerical solution of initial value problems Te metods you ve learned so far ave obtained closed-form solutions to initial value problems. A closedform solution is an explicit algebriac formula
More informationCommand shaping for a flexible satellite platform controlled by advanced fly-wheels systems. 1 Introduction
56th International Astronautical Congress, Paper IAC-5-C.1.3 Command shaping for a flexible satellite platform controlled by advanced fly-wheels systems Dario Izzo and Lorenzo Pettazzi ESA, Advanced Concepts
More informationA = h w (1) Error Analysis Physics 141
Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.
More informationParameter Fitted Scheme for Singularly Perturbed Delay Differential Equations
International Journal of Applied Science and Engineering 2013. 11, 4: 361-373 Parameter Fitted Sceme for Singularly Perturbed Delay Differential Equations Awoke Andargiea* and Y. N. Reddyb a b Department
More informationIMU-RGBD Camera Extrinsic Calibration: Observability Analysis and Consistency Improvement
IMU-RGBD Camera Extrinsic Calibration: Observability Analysis and Consistency Improvement Cao X. Guo and Stergios I. Roumeliotis Abstract In tis paper, we address te problem of extrinsically calibrating
More information1 Calculus. 1.1 Gradients and the Derivative. Q f(x+h) f(x)
Calculus. Gradients and te Derivative Q f(x+) δy P T δx R f(x) 0 x x+ Let P (x, f(x)) and Q(x+, f(x+)) denote two points on te curve of te function y = f(x) and let R denote te point of intersection of
More informationDifferentiation in higher dimensions
Capter 2 Differentiation in iger dimensions 2.1 Te Total Derivative Recall tat if f : R R is a 1-variable function, and a R, we say tat f is differentiable at x = a if and only if te ratio f(a+) f(a) tends
More informationOptimal Fault-Tolerant Configurations of Thrusters
Optimal Fault-Tolerant Configurations of Thrusters By Yasuhiro YOSHIMURA ) and Hirohisa KOJIMA, ) ) Aerospace Engineering, Tokyo Metropolitan University, Hino, Japan (Received June st, 7) Fault tolerance
More informationAPPLICATION OF A DIRAC DELTA DIS-INTEGRATION TECHNIQUE TO THE STATISTICS OF ORBITING OBJECTS
APPLICATION OF A DIRAC DELTA DIS-INTEGRATION TECHNIQUE TO THE STATISTICS OF ORBITING OBJECTS Dario Izzo Advanced Concepts Team, ESTEC, AG Noordwijk, Te Neterlands ABSTRACT In many problems related to te
More informationTaylor Series and the Mean Value Theorem of Derivatives
1 - Taylor Series and te Mean Value Teorem o Derivatives Te numerical solution o engineering and scientiic problems described by matematical models oten requires solving dierential equations. Dierential
More informationSolving Continuous Linear Least-Squares Problems by Iterated Projection
Solving Continuous Linear Least-Squares Problems by Iterated Projection by Ral Juengling Department o Computer Science, Portland State University PO Box 75 Portland, OR 977 USA Email: juenglin@cs.pdx.edu
More informationSymmetry Labeling of Molecular Energies
Capter 7. Symmetry Labeling of Molecular Energies Notes: Most of te material presented in tis capter is taken from Bunker and Jensen 1998, Cap. 6, and Bunker and Jensen 2005, Cap. 7. 7.1 Hamiltonian Symmetry
More informationBrazilian Journal of Physics, vol. 29, no. 1, March, Ensemble and their Parameter Dierentiation. A. K. Rajagopal. Naval Research Laboratory,
Brazilian Journal of Pysics, vol. 29, no. 1, Marc, 1999 61 Fractional Powers of Operators of sallis Ensemble and teir Parameter Dierentiation A. K. Rajagopal Naval Researc Laboratory, Wasington D. C. 2375-532,
More informationA MONTE CARLO ANALYSIS OF THE EFFECTS OF COVARIANCE ON PROPAGATED UNCERTAINTIES
A MONTE CARLO ANALYSIS OF THE EFFECTS OF COVARIANCE ON PROPAGATED UNCERTAINTIES Ronald Ainswort Hart Scientific, American Fork UT, USA ABSTRACT Reports of calibration typically provide total combined uncertainties
More informationA general articulation angle stability model for non-slewing articulated mobile cranes on slopes *
tecnical note 3 general articulation angle stability model for non-slewing articulated mobile cranes on slopes * J Wu, L uzzomi and M Hodkiewicz Scool of Mecanical and Cemical Engineering, University of
More informationClick here to see an animation of the derivative
Differentiation Massoud Malek Derivative Te concept of derivative is at te core of Calculus; It is a very powerful tool for understanding te beavior of matematical functions. It allows us to optimize functions,
More informationNumerical analysis of a free piston problem
MATHEMATICAL COMMUNICATIONS 573 Mat. Commun., Vol. 15, No. 2, pp. 573-585 (2010) Numerical analysis of a free piston problem Boris Mua 1 and Zvonimir Tutek 1, 1 Department of Matematics, University of
More information5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems
5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we
More informationarxiv: v1 [math.dg] 4 Feb 2015
CENTROID OF TRIANGLES ASSOCIATED WITH A CURVE arxiv:1502.01205v1 [mat.dg] 4 Feb 2015 Dong-Soo Kim and Dong Seo Kim Abstract. Arcimedes sowed tat te area between a parabola and any cord AB on te parabola
More informationNon-linear Analysis Method of Ground Response Using Equivalent Single-degree-of-freedom Model
Proceedings of te Tent Pacific Conference on Eartquake Engineering Building an Eartquake-Resilient Pacific 6-8 November 25, Sydney, Australia Non-linear Analysis Metod of Ground Response Using Equivalent
More informationLIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION
LIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION LAURA EVANS.. Introduction Not all differential equations can be explicitly solved for y. Tis can be problematic if we need to know te value of y
More informationIMU-RGBD Camera Extrinsic Calibration: Observability Analysis and Consistency Improvement
IMU-RGBD Camera Extrinsic Calibration: Observability Analysis and Consistency Improvement Cao X Guo and Stergios I Roumeliotis Multiple Autonomous Robotic Systems Laboratory Tecnical Report Number -212-1
More informationNONLINEAR SYSTEMS IDENTIFICATION USING THE VOLTERRA MODEL. Georgeta Budura
NONLINEAR SYSTEMS IDENTIFICATION USING THE VOLTERRA MODEL Georgeta Budura Politenica University of Timisoara, Faculty of Electronics and Telecommunications, Comm. Dep., georgeta.budura@etc.utt.ro Abstract:
More informationPoisson Equation in Sobolev Spaces
Poisson Equation in Sobolev Spaces OcMountain Dayligt Time. 6, 011 Today we discuss te Poisson equation in Sobolev spaces. It s existence, uniqueness, and regularity. Weak Solution. u = f in, u = g on
More informationDerivatives of Exponentials
mat 0 more on derivatives: day 0 Derivatives of Eponentials Recall tat DEFINITION... An eponential function as te form f () =a, were te base is a real number a > 0. Te domain of an eponential function
More informationMath 31A Discussion Notes Week 4 October 20 and October 22, 2015
Mat 3A Discussion Notes Week 4 October 20 and October 22, 205 To prepare for te first midterm, we ll spend tis week working eamples resembling te various problems you ve seen so far tis term. In tese notes
More informationChapter 2. Limits and Continuity 16( ) 16( 9) = = 001. Section 2.1 Rates of Change and Limits (pp ) Quick Review 2.1
Capter Limits and Continuity Section. Rates of Cange and Limits (pp. 969) Quick Review..... f ( ) ( ) ( ) 0 ( ) f ( ) f ( ) sin π sin π 0 f ( ). < < < 6. < c c < < c 7. < < < < < 8. 9. 0. c < d d < c
More informationSECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY
(Section 3.2: Derivative Functions and Differentiability) 3.2.1 SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY LEARNING OBJECTIVES Know, understand, and apply te Limit Definition of te Derivative
More informationMath 2921, spring, 2004 Notes, Part 3. April 2 version, changes from March 31 version starting on page 27.. Maps and di erential equations
Mat 9, spring, 4 Notes, Part 3. April version, canges from Marc 3 version starting on page 7.. Maps and di erential equations Horsesoe maps and di erential equations Tere are two main tecniques for detecting
More informationEFFICIENCY OF MODEL-ASSISTED REGRESSION ESTIMATORS IN SAMPLE SURVEYS
Statistica Sinica 24 2014, 395-414 doi:ttp://dx.doi.org/10.5705/ss.2012.064 EFFICIENCY OF MODEL-ASSISTED REGRESSION ESTIMATORS IN SAMPLE SURVEYS Jun Sao 1,2 and Seng Wang 3 1 East Cina Normal University,
More informationlecture 26: Richardson extrapolation
43 lecture 26: Ricardson extrapolation 35 Ricardson extrapolation, Romberg integration Trougout numerical analysis, one encounters procedures tat apply some simple approximation (eg, linear interpolation)
More informationChemical Engineering & Process Techniques
emical Engineering & Process Tecniques eview Article eedback ontrol for Liquid Level in a Gravity-Drained Multi-Tank System Larry K Jang* Department of emical Engineering, alifornia State University, USA
More informationImplicit-explicit variational integration of highly oscillatory problems
Implicit-explicit variational integration of igly oscillatory problems Ari Stern Structured Integrators Worksop April 9, 9 Stern, A., and E. Grinspun. Multiscale Model. Simul., to appear. arxiv:88.39 [mat.na].
More informationExam 1 Review Solutions
Exam Review Solutions Please also review te old quizzes, and be sure tat you understand te omework problems. General notes: () Always give an algebraic reason for your answer (graps are not sufficient),
More informationMath Spring 2013 Solutions to Assignment # 3 Completion Date: Wednesday May 15, (1/z) 2 (1/z 1) 2 = lim
Mat 311 - Spring 013 Solutions to Assignment # 3 Completion Date: Wednesday May 15, 013 Question 1. [p 56, #10 (a)] 4z Use te teorem of Sec. 17 to sow tat z (z 1) = 4. We ave z 4z (z 1) = z 0 4 (1/z) (1/z
More informationChemical Engineering & Process Techniques
emical Engineering & Process Tecniques eview Article eedback ontrol for Liquid Level in a Gravity-Drained Multi-Tank System Larry K. Jang* Department of emical Engineering, alifornia State University,
More information2.8 The Derivative as a Function
.8 Te Derivative as a Function Typically, we can find te derivative of a function f at many points of its domain: Definition. Suppose tat f is a function wic is differentiable at every point of an open
More informationSpacecraft Attitude Control using CMGs: Singularities and Global Controllability
1 / 28 Spacecraft Attitude Control using CMGs: Singularities and Global Controllability Sanjay Bhat TCS Innovation Labs Hyderabad International Workshop on Perspectives in Dynamical Systems and Control
More informationLecture XVII. Abstract We introduce the concept of directional derivative of a scalar function and discuss its relation with the gradient operator.
Lecture XVII Abstract We introduce te concept of directional derivative of a scalar function and discuss its relation wit te gradient operator. Directional derivative and gradient Te directional derivative
More informationMATH1151 Calculus Test S1 v2a
MATH5 Calculus Test 8 S va January 8, 5 Tese solutions were written and typed up by Brendan Trin Please be etical wit tis resource It is for te use of MatSOC members, so do not repost it on oter forums
More informationWYSE Academic Challenge 2004 Sectional Mathematics Solution Set
WYSE Academic Callenge 00 Sectional Matematics Solution Set. Answer: B. Since te equation can be written in te form x + y, we ave a major 5 semi-axis of lengt 5 and minor semi-axis of lengt. Tis means
More informationA finite element approximation for the quasi-static Maxwell Landau Lifshitz Gilbert equations
ANZIAM J. 54 (CTAC2012) pp.c681 C698, 2013 C681 A finite element approximation for te quasi-static Maxwell Landau Lifsitz Gilbert equations Kim-Ngan Le 1 T. Tran 2 (Received 31 October 2012; revised 29
More informationThe entransy dissipation minimization principle under given heat duty and heat transfer area conditions
Article Engineering Termopysics July 2011 Vol.56 No.19: 2071 2076 doi: 10.1007/s11434-010-4189-x SPECIAL TOPICS: Te entransy dissipation minimization principle under given eat duty and eat transfer area
More informationLIMITS AND DERIVATIVES CONDITIONS FOR THE EXISTENCE OF A LIMIT
LIMITS AND DERIVATIVES Te limit of a function is defined as te value of y tat te curve approaces, as x approaces a particular value. Te limit of f (x) as x approaces a is written as f (x) approaces, as
More informationNumerical Analysis MTH603. dy dt = = (0) , y n+1. We obtain yn. Therefore. and. Copyright Virtual University of Pakistan 1
Numerical Analysis MTH60 PREDICTOR CORRECTOR METHOD Te metods presented so far are called single-step metods, were we ave seen tat te computation of y at t n+ tat is y n+ requires te knowledge of y n only.
More informationNonlinear Tracking Control of Underactuated Surface Vessel
American Control Conference June -. Portland OR USA FrB. Nonlinear Tracking Control of Underactuated Surface Vessel Wenjie Dong and Yi Guo Abstract We consider in this paper the tracking control problem
More informationarxiv: v1 [math.oc] 18 May 2018
Derivative-Free Optimization Algoritms based on Non-Commutative Maps * Jan Feiling,, Amelie Zeller, and Cristian Ebenbauer arxiv:805.0748v [mat.oc] 8 May 08 Institute for Systems Teory and Automatic Control,
More information2.11 That s So Derivative
2.11 Tat s So Derivative Introduction to Differential Calculus Just as one defines instantaneous velocity in terms of average velocity, we now define te instantaneous rate of cange of a function at a point
More informationManipulator Dynamics (1) Read Chapter 6
Manipulator Dynamics (1) Read Capter 6 Wat is dynamics? Study te force (torque) required to cause te motion of robots just like engine power required to drive a automobile Most familiar formula: f = ma
More informationMA455 Manifolds Solutions 1 May 2008
MA455 Manifolds Solutions 1 May 2008 1. (i) Given real numbers a < b, find a diffeomorpism (a, b) R. Solution: For example first map (a, b) to (0, π/2) and ten map (0, π/2) diffeomorpically to R using
More informationLecture 15. Interpolation II. 2 Piecewise polynomial interpolation Hermite splines
Lecture 5 Interpolation II Introduction In te previous lecture we focused primarily on polynomial interpolation of a set of n points. A difficulty we observed is tat wen n is large, our polynomial as to
More informationMotion planning for cooperating mobile manipulators
Motion planning for cooperating mobile manipulators Jaydev P. Desai Division of Engineering and Applied Sciences 9 Oxford Street, G8 Pierce Hall Harvard University Cambridge, MA 38, USA Vijay Kumar General
More informationMAE 142 Homework #5 Due Friday, March 13, 2009
MAE 142 Homework #5 Due Friday, March 13, 2009 Please read through the entire homework set before beginning. Also, please label clearly your answers and summarize your findings as concisely as possible.
More information158 Calculus and Structures
58 Calculus and Structures CHAPTER PROPERTIES OF DERIVATIVES AND DIFFERENTIATION BY THE EASY WAY. Calculus and Structures 59 Copyrigt Capter PROPERTIES OF DERIVATIVES. INTRODUCTION In te last capter you
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 34, pp. 14-19, 2008. Copyrigt 2008,. ISSN 1068-9613. ETNA A NOTE ON NUMERICALLY CONSISTENT INITIAL VALUES FOR HIGH INDEX DIFFERENTIAL-ALGEBRAIC EQUATIONS
More informationStabilization of a 3D Rigid Pendulum
25 American Control Conference June 8-, 25. Portland, OR, USA ThC5.6 Stabilization of a 3D Rigid Pendulum Nalin A. Chaturvedi, Fabio Bacconi, Amit K. Sanyal, Dennis Bernstein, N. Harris McClamroch Department
More information*School of Aeronautics and Astronautics, Purdue University, West Lafayette, IN
NEW CONTROL LAWS FOR THE ATTITUDE STABILIZATION OF RIGID BODIES PANAGIOTIS TSIOTRAS *School of Aeronautics and Astronautics, Purdue University, West Lafayette, IN 7907. Abstract. This paper introduces
More informationMAT 145. Type of Calculator Used TI-89 Titanium 100 points Score 100 possible points
MAT 15 Test #2 Name Solution Guide Type of Calculator Used TI-89 Titanium 100 points Score 100 possible points Use te grap of a function sown ere as you respond to questions 1 to 8. 1. lim f (x) 0 2. lim
More informationDepartment of Statistics & Operations Research, Aligarh Muslim University, Aligarh, India
Open Journal of Optimization, 04, 3, 68-78 Publised Online December 04 in SciRes. ttp://www.scirp.org/ournal/oop ttp://dx.doi.org/0.436/oop.04.34007 Compromise Allocation for Combined Ratio Estimates of
More informationUnit 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 informationERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER*
EO BOUNDS FO THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BADLEY J. LUCIE* Abstract. Te expected error in L ) attimet for Glimm s sceme wen applied to a scalar conservation law is bounded by + 2 ) ) /2 T
More information(a) At what number x = a does f have a removable discontinuity? What value f(a) should be assigned to f at x = a in order to make f continuous at a?
Solutions to Test 1 Fall 016 1pt 1. Te grap of a function f(x) is sown at rigt below. Part I. State te value of eac limit. If a limit is infinite, state weter it is or. If a limit does not exist (but is
More information