Robust Static Attitude Determination via Robust Optimization

Size: px
Start display at page:

Download "Robust Static Attitude Determination via Robust Optimization"

Transcription

1 Milano (Italy) August 8 - September, Robust Static Attitude Deteration via Robust Optimization Shakil Ahmed Eric C. Kerrigan Department of Electrical and Electronic Engineering, Imperial College London, UK ( shakil.ahmed8@imperial.ac.uk). Department of Aeronautics and Department of Electrical and Electronic Engineering, Imperial College London, UK ( e.kerrigan@imperial.ac.uk). Abstract: We address the problem of robust attitude deteration using a static approach. In contrast to the dynamic approach, a static approach does not depend on the system dynamics. This approach only requires measurements of some vectors, such as the earth magnetic field, sun vector, etc in two different coordinate frames. These vectors, obtained from some sensor or a mathematical model, may not be accurate due to sensor errors and modeling inaccuracies. We consider all such errors as infinity-norm bounded uncertainties and our main focus is to obtain an attitude estimate, which is least sensitive to such uncertainties. We formulate a robust optimization (RO) problem with a quadratic cost and nonlinear constraints and propose a solution using a quaternion approach with an affine uncertainty parameterization. We transform the RO problem into a suboptimal imization problem, which is non-convex butcan be solved with a good initial guess using a nonlinear optimization solver. The results show a significant advantage of the proposed approach for the worst case uncertainties. Keywords: Transformation matrices, attitude algorithms, modeling errors, least square estimation, robust estimation. INTRODUCTION The main focus of this work is on robust attitude deteration in the presence of uncertainties using a static approach, whichis basedonvector information. Onemain advantage of this approach is that it always gives an attitude estimate independent of the dynamic system nonlinearities, making it especially useful for highly nonlinear systems, such as a tumbling satellite, where dynamic approaches based on linear or nonlinear filtering suffer from divergence issues due to lack of good a priori state estimate (Crassidis, Markley, andcheng, 7). Anadded benefit of static estimation in such cases is that it can provide a good state initialization for dynamic filters, reducing the likelihood of divergence. To compute the attitude of an object i.e. its orientation in space with respect to some known reference, two coordinate frames are needed. One, which is fixed to the body of the object, is called the bodyframe, while the secondis called the reference frame. Selection of the reference frame normally depends upon the control system requirements. Formally, attitude of an object is defined as a coordinate transformation that transforms reference coordinates into the body coordinates (Sidi, ). This transformation is obtained through a proper orthogonal transformation matrix C R 3 3, also known as a direction cosine or attitude matrix. Orthogonality of the transformation matrix dictates the constraint C T C = I, while the proper transformation matrix is defined by having the deterant equal to + to preserve the orientation in a rotation, thus imposing the constraint det(c) = +. The vector quantities normally used in static attitude deteration are the earth magnetic field (EMF), sun vector, star vector, etc. Information of these vectors is required both in the body and reference frames in order to detere the attitude or transformation matrix. Normally the vectors in the body frame are measured by some sensor installed on the body of the object, while the same vector information in the reference frame is obtained from some mathematical model. Withonly onepairof information, i.e. a vector quantity in the bodyandreference frame, deteration of the attitude matrix remains an underdetered problem with infinite solutions. Addingasecond pair makes this problem overdetered (Shuster and Oh, 98). A commonly used approach for such problems is to findanoptimalsolutionthatimizes theweightedleast square cost, first proposed by Wahba (965) for satellite applications, given as: C w i b i Cr i subject to C T C = I, det(c) =, where b i R 3 represents i th measurement in the body frame and i =,..,n, n being the total number of sensors, r i R 3 is the corresponding vector in the reference frame obtained from some model, w i R are non-negative weights while represents vector -norm. Many efficient solutions of this problem for satellite applications can be found in the literature, such as Keat (977), Shuster and Oh (98), Markley (988), Mortari (997) and Mortari (). () Copyright by the International Federation of Automatic Control (IFAC) 587

2 Milano (Italy) August 8 - September, A general approach used in solving the imization problem () is to convert it into an equivalent imization problem. For this, consider the cost function used in () for imization and expand it as w i b i Cr i = w i {b T i b i +ri T r i } w i b T i Cr i. Neglecting the constant term as it will not change argument of the optimization problem, an equivalent imization problem can be written as C w i b T i Cr i subject to C T C = I, det(c) =. To solve this imization problem the most popular approach is based on Davenport s q-method (Keat, 977; Shuster and Oh, 98). Two important steps of the q- method are given now, which will be used for solving the problem addressed in this work. Step : Find an equivalent formulation of () in terms of a quaternion. This new formulation, firstreported in Keat (977), states that the imization of () is equivalent to the following problem (see Appendix A for derivation): q q T Kq subject to q T q =, where q R 4 represents a 4-parameter quaternion, while K R 4 4 is a symmetric, indefinite and traceless matrix given as [ ] B K := T + B tr(b)i z z T, (4) tr(b) where B R 3 3 := n w ib i r T i and z R 3 := n w i(b i r i ), while represents cross product and tr( ) represents trace operator. Step : Using the new formulation, the non-concave imization problem can be converted into an eigenvalue problem. For this we first add the constraint q T q = using a Lagrange multiplier λ in (3) that gives () (3) f(q,λ) = q T Kq λ(q T q ), (5) To obtain a stationary point, we solve f/ q = and f/ λ = and obtain the following expression, that has the same formulation as the standard eigenvalue problem Kq = λq. (6) Four eigenvectors of the symmetric matrix K are possible solutions of (6); however, the eigenvector corresponding to the imum eigenvalue is the solution that will solve (3) (Shuster and Oh, 98),(Markley and Mortari, ), i.e. Kq opt = λ (K)q opt, (7) where q opt is the solution to (3) and λ (K) is the imum eigenvalue of the matrix K. Most of the work in static satellite attitude deteration is based on this result and many efficient algorithms have been proposed, such as QUEST (Shuster and Oh, 98), ESOQ (Mortari, 997), ESOQ (Mortari, ). A survey paper by Markley andmortari () provides ageneraldescriptionofmany of these algorithms. The use of a quaternion in the new formulation imparts many benefits. It not only reduces the number of optimization variables from nine to four, but also avoids the constraint det(c) = + of () being inherent in the definition of quaternion. However, the main benefit obtained is the transformation of the optimization problem into an eigenvalue problem. This classical approach, based on () does not directly address the issue of robustness in the measured and model vectors. Mostly, sensitivity analysis have been presented for different algorithms with analytical expressions of the imum error covariance under stochastic variation in the measurement vector, such as Shuster and Oh (98), Shuster (6), Markley (988), Markley (8). However, modeling errors in the reference vectors are generally not considered. These errors could be significant. In the case of the earth magnetic field, which is a commonly used sensor especially in satellite applications, errors between sophisticated models and the actual field can be significant (McLean et al., December, 4; Roithmayr, 4). The use of simple models, such as low order IGRF (Roithmayr, 4), which are normally preferred due to less computational cost, results in a less accurate earth magnetic field vector in the reference frame, leading to errors in the attitude estimate. Attitude inaccuracy is further increased due to sensor errors, which are mainly due to noise and installation issues. The magnetic field sensing in the post launch tumbling phase is also a big source of uncertainty. We consider all such errors as - norm bounded uncertainties in the input vectors. The main contribution of this paper is consideration of the measurement and modeling uncertainties in the problem formulation for static attitude deteration. We formulate a robust optimization (RO) problem for the uncertain inputs. The formulated - RO problem is relaxed and transformed to a suboptimal imization problem using a suitable affine uncertainty parameterization. The resultingnon-convex imization problem with nonlinear cost and constraints is solved using a nonlinear optimization solver with a good initial guess.. ROBUST PROBLEM DESCRIPTION To formulate a robust attitude deteration problem, we will represent an uncertain measurement vector in the body frame with b i B i R 3 and an uncertain reference vector with r i R i R 3, i =,...,n, where B i and R i are bounded uncertainty sets. To formulate the problem of detering the best uncertainty immunized transformation matrix for attitude deteration, we will use the weighted least square approach as used in () for the noal problem. We define the robust problem as C bi B i, r i R i, i =,..., n w i bi C r i subject to C T C = I,det(C) =. In order to take advantage of using a quaternion to simplify the optimization problem, as a first step, we (8) 588

3 Milano (Italy) August 8 - September, reformulate (8) introducing the quaternion q using the same approach used to derive (3). We get { } w i ( b T i q b i + r i T r i ) q T Kq bi B i, r i R i, i =,..., n subject to q T q =, where (9) [ B + K := BT tr( B)I ] z z T tr( B), B := w i bi r i T, z := w i ( b i r i ). () In this formulation, similar to (4), K R 4 4 is also a symmetric indefinite matrix, while B R 3 3 and z R 3 depend on uncertain input vectors. In contrast to (), the first term in (9) is a function of optimization variables and cannot be neglected in the optimization problem formulation. 3. UNCERTAINTY MODELING To develop a tractable method of solving (9), we define an affine uncertainty parameterization of the uncertainty sets B i and R i. A general description of such an affine parameterization is given now (Ben-Tal, Ghaoui, and Nemirovski, 9). Let β i := [β i,β i,β i3 ] T and ρ i := [ρ i,ρ i,ρ i3 ] T are vectors of perturbation variables used for uncertainty parameterization, then B i := { b i = b i + R i := { r i = r i + β ilˆbl : β i γ bi }, l= ρ ilˆr l : ρ i γ ri }, () l= where γ bi and γ ri represent uncertainty bounds for the input vectors in the body frame and reference frame, respectively. This type of uncertainty is called an interval uncertainty and correspondingperturbation set represents a box (Ben-Tal et al., 9). The interval uncertainty can be considered a suitable representation for a mix of different type of errors, such as sensor noise, having a stochastic interpretation, biases andmisalignments, which are constant and modeling inaccuracies, etc. However, bounds for each measurement or model vector should be carefully chosen, as unnecessary large values may result in large residual between the non-robust and the robust solution for noal cases. As we considered all vectors in the body and reference frame to be unit vectors, we choose ˆb l and ˆr l as ˆb = ˆr = [ ] T,ˆb = ˆr = [ ] T,ˆb 3 = ˆr 3 = [ ] T. With this choice, γ bi and γ ri represent the imum possible change in each value of b i and r i as a percentage of the input vector magnitude. Although, each perturbation variable in the body and reference frame has different bound, we can use normalization to obtain a single perturbation variable for uncertainty parameterization. For this, we normalize β il and ρ il with γ bi and γ ri, respectively. The new normalized perturbation vector is δ := [δ,δ,δ 3 ], such that δ. Using this normalization, we can represent the above mentioned sets B i and R i as B i = { b i = b i + R i = { r i := r i + δ l bil : δ }, l= δ l r il : δ }, () l= where b il := γ biˆbl and r il := γ riˆr l. 4. SOLUTION OF THE ROBUST PROBLEM To solve (9), we transform the robust - problem into aimization problem. Inthis regard, firstly, we will derive an expression for the uncertain matrix K, similar to the K matrix of (3), using the presented uncertainty model. This will provide a basis for the subsequent results that describe the simplified but a suboptimal formulation obtained for the robust problem (9). Using the uncertainty model described in Section 3, matrix K defined in (), can be written as K = K + (δ l Kl + δl K l s + δ l δ k Kc jk ), (3) l= where the superscripts s and c represent matrices with square and cross terms in δ,k = l+ if k 3 or k = l+ n if k > 3. The matrices K l, K s l, K c lk R4 4 are given as: [( Bl + K l := B l T tr( B ) ] l )I z l z l T tr( B, l ) [( K l s st Bs := l + B l tr( B ) ] l s )I z l s z l st tr( B, l s ) [( K lk c ct Bc := lk + B lk tr( B ) ] lk)i c z lk c tr( B. (4) lk) c z ct lk The other terms used in these expressions are givens as: B l := B c lk := z l := z lk := w i (b i r il T + b il ri T ), B l s := w i ( b il r ik T + b ik r il), T w i bil r il, T w i {(b i r il ) + ( b il r i )}, z l s := w i {( b il r ik ) + w i ( b il r il ), w i {( b ik r il ). (5) It can be observed that Kl depends on measurements, but Ks l and K lk c depend on uncertainty modeling and remain fixed for different measurement vectors, provided the uncertainty description remains same. We now present two results using (3), that will be used to reformulate (9) in terms of the proposed uncertainty model. 589

4 Milano (Italy) August 8 - September, Lemma. The formulation given in (9) is equivalent to q subject to q T q =, { q T Kq+ δ (pt δ + δ T Qδ)} (6) where p R 3 and Q R 3 3 depend on q and are defined as c q T K q p := c q T K q, c 3 q T K3 q c s q T Ks q qt Kc q qt Kc 3 q Q := qt Kc q c s q T Ks q qt Kc 3 q qt Kc 3 q. qt Kc 3 q c s 3 q T Ks 3 q In this equation, c l, c s l are constants that depend on the uncertain vectors and are given as: c l := c s l := w i (b T b i il + b T ilb i + ri T r il + r ilr T i ), w i ( b T il b il + r il T r il ). (7) Proof. We start with the cost given in (9). Expanding the first term, we get w i ( b T b i i + r i T r i ) = w i (b T i b i + ri T r i ) + (δ l c l + δl c s l ) (8) Using (8) and (3) and joining similar terms in δ l we can write the required expression. The robust problem (6) is always feasible because the solution of the noal problem (3) always exists (Shuster and Oh, 98). However, finding the optimal imum value of p T δ + δ T Qδ in (6) is difficult as the term is nonconcave in δ (Q is positive semidefinite). Based on this fact, we detere an upper bound on the imum of (p T δ + δ T Qδ) over δ. The result is given in the following lemma. Lemma. An upper bound on the imization term appearing in (6) is l= δ (pt δ + δ T Qδ) p + 3λ (Q) (9) Proof. The left hand side in (9) can be written as δ (pt δ + δ T Qδ) δ pt δ + δ δt Qδ () Using the Hölder dual norm (Boyd and Vandenberghe, 4), the first term on the right hand side in () is given as δ pt δ = p. () For the second term appearing in (), since Q is a symmetric matrix, we can write the imum eigenvalue of Q as (Boyd and Vandenberghe, 4) λ (Q) = sup δ T Qδ δ. () Hence, we first replace the -norm in the second term on the right hand side of () with the -norm using the inequality δ 3 δ for δ R 3. We can write δ δt Qδ δ δ T Qδ 3 3λ (Q), Using () and (3), we can write (9). (3) Lemma gives an upper bound on the imum value of the quadratic term p T δ+δ T Qδ and can be used to simplify (6). Theorem. The - problem (6) can be approximated with the following imization problem q,u,u,u 3 q T Kq + u + u + u 3 + 3λ (Q) subject to q T q =, u j p j u j,j =,,3, where p j,j =,,3 are elements of vector p. (4) Proof. In formulation (6) given in Lemma, replacing the term with the upper bound given in Lemma and using the fact that a set of n+ linear inequalities u j x j u j, j u j,j =,...,3 represent the nonlinear inequality j x j (Ben-Tal, Ghaoui, and Nemirovski, 9, definition.3.), we can write (4). The argument q of the approximate robust problem (4) approaches the argument of the noal problem (3) when the bound on the uncertainty approaches zero. This becomes evident if, instead of taking δ, we assume a bound α on each element of the normalized vector δ. In that case we can write (9) as δ α (pt δ + δ T Qδ) α p + 3α λ (Q). (5) If α tends to zero, the right hand side of (5) will also tend to zero, making the robust problem equivalent to (3). Since the cost and constraints in (4) are nonlinear, we solve it using a nonlinear optimization solver with a good initial guess. This initial guess for q is obtained from the non-robust form of the problem using the q-method or QUEST (Shuster and Oh, 98), etc. 5. SIMULATION RESULTS We consider the problem of attitude estimation for a low cost CubeSat (Heidt et al., ), a pico class of satellite having three units each of dimension cm 3 and weight less than kg. CubeSats are mainly used for scientific missions and currently most of the satellites are beingdeveloped inuniversities. Beinglow cost, interestin suchsatellites is increasingfor awiderclass ofapplications; however, there are many challenges from an estimation and control point of view to successfully achieve a given objective. The CubeSat under consideration is assumed to be moving in a circular orbit having a radius of 65 km from the center of the earth. We are considering that only two measurements are available in the satellite, namely the earth magnetic field and the sun vector. For the earth magnetic field, two magnetometers are installed, one inside the satellite, which is mainly used in the post- 58

5 Milano (Italy) August 8 - September, launch phase when the satellite is recovering from launch disturbances, while the second is installed on an extended boom, which is deployed once the satellite has de-tumbled and achieved an equilibrium. The sun vector is sensed by a pair of sun sensors installed on the satellite body. Both of these measurements are in the body frame. We used the orbit frame as the reference frame. For the reference earth magnetic field vector the first order IGRF model (Roithmayr, 4) is used, while the reference sun vector is obtained using a simplified sun model based on sun ephemeris (Wertz, 978). Sensor measurements are not accurate due to many factors, such as sensor noise, misalignments due to installation errors, biases, etc. Especially in the post-launch tumbling phase measurement errors further increase due to the use of an internal magnetometer installed on-board the satellite, which suffers from interaction with the magnetic field generated by surrounding electronics. Similarly, the reference vectors for both measurements are also not exact since the mathematical models used to obtain measurement vectors in the reference frame are based on loworder approximations. For simulation purpose, we consider that the bound on uncertainty in the measurements and reference vectors i.e. γ bi and γ ri, respectively to be % of the input magnitude. We take equal weights, such that n w i =. We need to solve the optimization problem formulated in (4). Since this formulation includes nonlinear objective function and constraints, we used the nonlinear optimization solver (fcon) of MATLAB. Two plots are presented to see the effectiveness of the proposed methodology. Firstly, the problem was solved for a given uncertain data set for one instant. To obtain the uncertain data set, we added a uniformally distributed random error in the range of ±γ bi and ±γ ri in each measured and model vectors obtained from a simulation without error. The vectors with added uncertainty are given as: b = [ ] T, r = [ ] T, b = [ ] T, r = [ ] T. Figures presents a comparison between the non-robust and the robust approach for one time instant. For this comparison, anon-robustsolution was obtained satisfying (3), while a robust solution was found satisfying (4). Then, two components ofthenormalized uncertainty vectorδ, i.e. δ and δ are varied, while keeping δ 3 constant for simulation purpose, and a variation in the cost J := n b i C r i is calculated for both solutions. This cost is plotted along the z-axis for both solutions, while x- and y-axis represent δ and δ. It can be observed from both plots that the cost J obtained using the noal solution has large variation when δ and δ are close to set bounds, while the J obtained using the robust solution shows the best immunization against uncertainty. However, the cost obtained for the robust solution has an offset when compared with the cost value obtained using the noal solution without error. Figure presents a performance comparison of the robust and non-robust approaches in the presence of uncertain- J δ.5.5 δ non-robust robust Fig.. Robust and non-robust design comparison for given set of uncertain vectors at one instant. φ ψ θ Time (sec).5 non-robust robust actual Fig.. Comparison of attitude angles obtained using nonrobust and robust algorithms. The dotted line shows the data without errors while the other two cases include errors within the uncertainty bound ±γ bi or ±γ ri. ties, using in-orbit data obtained from nonlinear closedloop simulation for the satellite. The error free data obtained from simulation was corrupted by adding uniformly distributed random errors in the range of ±γ bi and ±γ ri in corresponding vectors. The simulation was run with a sample time of second with initial roll, pitch and yaw body rates of.5,.5 and deg/s and roll, pitch and yaw angles of 3,, deg, respectively. It can be observed that due to the uncertainties in the input information, the non-robust approach gives large errors in the estimated attitude angles, while the robust approach gives much better performance, limiting the imum attitude error to a smaller band. 6. CONCLUSIONS We presented a formulation for robust attitude deteration with -norm bounded uncertainty in the input vectors. The robust optimization problem, using a weighted 58

6 Milano (Italy) August 8 - September, least square approach, has a nonlinear cost and constraints. The formulated - optimization problem was transformed into a suboptimal imization problem due to the use of an upper bound and solved using a nonlinear optimization solver with a good initial guess obtained from the noal solution. The results showed that the robust approach has significant benefit over the noal approach in the situations when inputs are uncertain with boundeduncertainty. Thebenefitis generallyimum in the worst case scenarios. However, the robust solution for the case without error has a larger residual than the nonrobust solution, which is expected in the robust optimal solutions. However, the main issue is increased computational burden in solving the optimization problem. Further work will focus on finding a computationally efficient solution of this problem. ACKNOWLEDGEMENTS The first author is thankful to the Higher Education Commission of Pakistan for supporting his PhD studies. REFERENCES Ben-Tal, A., Ghaoui, L.E., and Nemirovski, A. (9). Robust Optimization. Princeton University Press, USA. Boyd, S. and Vandenberghe, L. (4). Convex Optimization. Cambridge University Press. Crassidis, J.L., Markley, F.L., and Cheng, Y. (7). Survey of nonlinear attitude estimation methods. Journal of Guidance, Control and Dynamics, 3(), 8. Heidt, H., Suari, J.P., Moore, A.S., Nakasuka, S., and Twiggs, R.J. (). Cubesat: A new generation of pico-satellite for education and industry low cost space experimentation. In 5th Annual/USU Conference on Small Satellite. SSC-VIIIb-5. Keat, J.E. (977). Analysis of least-squares attitude deteration routine DOAOP. Technical Report NASA/CSC/TM-77/634. Markley, F.L. (988). Attitude deteration using vector observations and the singular value decomposition. Journal of Astronatutical Sciences, 36(3), Markley, F.L. (8). Optimal attitude matrix from two vector measurements. Journal of Guidance, Control and Dynamics, 3(3), Markley, F.L. and Mortari, D. (). New developments in quaternion estimation from vector observation. In Advances in the Astronautical Sciences, volume 6, p McLean, S., Macmillan, S., Maus, S., Lesur, V., Thomson, A., and Dater, D. (December, 4). The US/UK world magnetic model for 5-. Technical Report NOAA, NESDIS/NGDC-. Mortari, D. (997). ESOQ: A closed form solution to the Wahba problem. Journal of the Astronautical Sciences, 45(), Mortari, D. (). Second estimator of the optimal quaternion. Journal of Guidance, Control and Dynamics, 3(5), Roithmayr, C.M. (4). Contributions of spherical harmonics to magnetic and gravitational fields. Technical Report NASA/TM Shuster, M.D. and Oh, S. (98). Three-axis attitude deteration from vector observation. Journal of guidance and Control, 4(), Shuster, M.D. (6). The generalized wahba problem. Journal of Astronautical Sciences, 54(), Sidi, M.J. (). Spacecraft Dynamics and Control. Cambridge University Press, UK, nd edition. Wahba, G. (965). A least square estimate of spacecraft attitude. SIAM Review, 7(3), 49. Wertz, J.R. (978). Spacecraft AttitudeDeteration and Control. D. Reidel Publishing Company, Holland. Appendix A. DAVENPORT TRANSFORMATION To derive the Davenport transformation, consider the cost function given in () w i b T i Cr i = tr(w i b T i Cr i ). (A.) Now we use two properties of the trace. Firstly trace is invariant under cyclic permutations, and secondly, tr ( i A i) = i tr(a i) A R n n. Using these properties we can write w i b T i Cr i = tr(cb T ), (A.) where B := n w ib i ri T. Now we represent C using the quaternion q R 4 := [ ] T, q T q 4 where q := [q q q 3 ] T, written as (Shuster and Oh, 98) C = (q4 q T q)i + qq T + q 4 Q. (A.3) In the above equation Q := q, where represents vector cross product, is given as [ q3 ] q Q = q 3 q. (A.4) q q Substituting (A.3) in (A.), we get w i b T i Cr i = (q4 q T q)tr(b T ) + tr(qq T B T ) + q 4 tr(qb T ). (A.5) The second term on the right hand side of the above equation can be written as tr(qq T B T ) = q T (B T + B)q. The last term can be written as (A.6) q 4 tr(qb T ) = q 4 (q T z + z T q), (A.7) where z := n w i(b i r i ). Substituting (A.6) and (A.7) in (A.5), we get w i b T i Cr i = [ ] [ ] q T B T + B tr(b)i z q q 4 z tr(b)][ T q 4 = q T Kq, (A.8) which is the required form of (A.), where K R 4 4 is a symmetric, traceless matrix given as [ ] B K := T + B tr(b)i z z T. (A.9) tr(b) 58

HOW TO ESTIMATE ATTITUDE FROM VECTOR OBSERVATIONS

HOW TO ESTIMATE ATTITUDE FROM VECTOR OBSERVATIONS AAS 99-427 WAHBA S PROBLEM HOW TO ESTIMATE ATTITUDE FROM VECTOR OBSERVATIONS F. Landis Markley* and Daniele Mortari The most robust estimators minimizing Wahba s loss function are Davenport s q method

More information

Generalized Attitude Determination with One Dominant Vector Observation

Generalized Attitude Determination with One Dominant Vector Observation Generalized Attitude Determination with One Dominant Vector Observation John L. Crassidis University at Buffalo, State University of New Yor, Amherst, New Yor, 460-4400 Yang Cheng Mississippi State University,

More information

Design Architecture of Attitude Determination and Control System of ICUBE

Design Architecture of Attitude Determination and Control System of ICUBE Design Architecture of Attitude Determination and Control System of ICUBE 9th Annual Spring CubeSat Developers' Workshop, USA Author : Co-Author: Affiliation: Naqvi Najam Abbas Dr. Li YanJun Space Academy,

More information

Attitude Estimation Based on Solution of System of Polynomials via Homotopy Continuation

Attitude Estimation Based on Solution of System of Polynomials via Homotopy Continuation Attitude Estimation Based on Solution of System of Polynomials via Homotopy Continuation Yang Cheng Mississippi State University, Mississippi State, MS 3976 John L. Crassidis University at Buffalo, State

More information

Adaptive Unscented Kalman Filter with Multiple Fading Factors for Pico Satellite Attitude Estimation

Adaptive Unscented Kalman Filter with Multiple Fading Factors for Pico Satellite Attitude Estimation Adaptive Unscented Kalman Filter with Multiple Fading Factors for Pico Satellite Attitude Estimation Halil Ersin Söken and Chingiz Hajiyev Aeronautics and Astronautics Faculty Istanbul Technical University

More information

Attitude Determination using Infrared Earth Horizon Sensors

Attitude Determination using Infrared Earth Horizon Sensors SSC14-VIII-3 Attitude Determination using Infrared Earth Horizon Sensors Tam Nguyen Department of Aeronautics and Astronautics, Massachusetts Institute of Technology 77 Massachusetts Avenue, Cambridge,

More information

Lecture: Examples of LP, SOCP and SDP

Lecture: Examples of LP, SOCP and SDP 1/34 Lecture: Examples of LP, SOCP and SDP Zaiwen Wen Beijing International Center For Mathematical Research Peking University http://bicmr.pku.edu.cn/~wenzw/bigdata2018.html wenzw@pku.edu.cn Acknowledgement:

More information

QUATERNION ATTITUDE ESTIMATION USING VECTOR OBSERVATIONS

QUATERNION ATTITUDE ESTIMATION USING VECTOR OBSERVATIONS QUATERNION ATTITUDE ESTIMATION USING VECTOR OBSERVATIONS F. Landis Markley 1 and Daniele Mortari 2 ABSTRACT This paper contains a critical comparison of estimators minimizing Wahba s loss function. Some

More information

DECENTRALIZED ATTITUDE ESTIMATION USING A QUATERNION COVARIANCE INTERSECTION APPROACH

DECENTRALIZED ATTITUDE ESTIMATION USING A QUATERNION COVARIANCE INTERSECTION APPROACH DECENTRALIZED ATTITUDE ESTIMATION USING A QUATERNION COVARIANCE INTERSECTION APPROACH John L. Crassidis, Yang Cheng, Christopher K. Nebelecky, and Adam M. Fosbury ABSTRACT This paper derives an approach

More information

Pointing Control for Low Altitude Triple Cubesat Space Darts

Pointing Control for Low Altitude Triple Cubesat Space Darts Pointing Control for Low Altitude Triple Cubesat Space Darts August 12 th, 2009 U.S. Naval Research Laboratory Washington, D.C. Code 8231-Attitude Control System James Armstrong, Craig Casey, Glenn Creamer,

More information

Relaxations and Randomized Methods for Nonconvex QCQPs

Relaxations and Randomized Methods for Nonconvex QCQPs Relaxations and Randomized Methods for Nonconvex QCQPs Alexandre d Aspremont, Stephen Boyd EE392o, Stanford University Autumn, 2003 Introduction While some special classes of nonconvex problems can be

More information

An Efficient and Robust Singular Value Method for Star Pattern Recognition and Attitude Determination

An Efficient and Robust Singular Value Method for Star Pattern Recognition and Attitude Determination The Journal of the Astronautical Sciences, Vol. 52, No. 1 2, January June 2004, p. 000 An Efficient and Robust Singular Value Method for Star Pattern Recognition and Attitude Determination Jer-Nan Juang

More information

Semidefinite and Second Order Cone Programming Seminar Fall 2012 Project: Robust Optimization and its Application of Robust Portfolio Optimization

Semidefinite and Second Order Cone Programming Seminar Fall 2012 Project: Robust Optimization and its Application of Robust Portfolio Optimization Semidefinite and Second Order Cone Programming Seminar Fall 2012 Project: Robust Optimization and its Application of Robust Portfolio Optimization Instructor: Farid Alizadeh Author: Ai Kagawa 12/12/2012

More information

Applications of Robust Optimization in Signal Processing: Beamforming and Power Control Fall 2012

Applications of Robust Optimization in Signal Processing: Beamforming and Power Control Fall 2012 Applications of Robust Optimization in Signal Processing: Beamforg and Power Control Fall 2012 Instructor: Farid Alizadeh Scribe: Shunqiao Sun 12/09/2012 1 Overview In this presentation, we study the applications

More information

Space Surveillance with Star Trackers. Part II: Orbit Estimation

Space Surveillance with Star Trackers. Part II: Orbit Estimation AAS -3 Space Surveillance with Star Trackers. Part II: Orbit Estimation Ossama Abdelkhalik, Daniele Mortari, and John L. Junkins Texas A&M University, College Station, Texas 7783-3 Abstract The problem

More information

Quaternion Data Fusion

Quaternion Data Fusion Quaternion Data Fusion Yang Cheng Mississippi State University, Mississippi State, MS 39762-5501, USA William D. Banas and John L. Crassidis University at Buffalo, State University of New York, Buffalo,

More information

Lecture 7: Weak Duality

Lecture 7: Weak Duality EE 227A: Conve Optimization and Applications February 7, 2012 Lecture 7: Weak Duality Lecturer: Laurent El Ghaoui 7.1 Lagrange Dual problem 7.1.1 Primal problem In this section, we consider a possibly

More information

Distributionally Robust Convex Optimization

Distributionally Robust Convex Optimization Submitted to Operations Research manuscript OPRE-2013-02-060 Authors are encouraged to submit new papers to INFORMS journals by means of a style file template, which includes the journal title. However,

More information

Lecture Note 5: Semidefinite Programming for Stability Analysis

Lecture Note 5: Semidefinite Programming for Stability Analysis ECE7850: Hybrid Systems:Theory and Applications Lecture Note 5: Semidefinite Programming for Stability Analysis Wei Zhang Assistant Professor Department of Electrical and Computer Engineering Ohio State

More information

arxiv: v1 [math.oc] 10 Oct 2014

arxiv: v1 [math.oc] 10 Oct 2014 A convex solution to Psiaki s first joint attitude and spin-rate estimation problem James Saunderson Pablo A. Parrilo Alan S. Willsky arxiv:1410.2841v1 [math.oc] 10 Oct 2014 October 13, 2014 Abstract We

More information

ELE539A: Optimization of Communication Systems Lecture 15: Semidefinite Programming, Detection and Estimation Applications

ELE539A: Optimization of Communication Systems Lecture 15: Semidefinite Programming, Detection and Estimation Applications ELE539A: Optimization of Communication Systems Lecture 15: Semidefinite Programming, Detection and Estimation Applications Professor M. Chiang Electrical Engineering Department, Princeton University March

More information

(3.1) a 2nd-order vector differential equation, as the two 1st-order vector differential equations (3.3)

(3.1) a 2nd-order vector differential equation, as the two 1st-order vector differential equations (3.3) Chapter 3 Kinematics As noted in the Introduction, the study of dynamics can be decomposed into the study of kinematics and kinetics. For the translational motion of a particle of mass m, this decomposition

More information

Selected Examples of CONIC DUALITY AT WORK Robust Linear Optimization Synthesis of Linear Controllers Matrix Cube Theorem A.

Selected Examples of CONIC DUALITY AT WORK Robust Linear Optimization Synthesis of Linear Controllers Matrix Cube Theorem A. . Selected Examples of CONIC DUALITY AT WORK Robust Linear Optimization Synthesis of Linear Controllers Matrix Cube Theorem A. Nemirovski Arkadi.Nemirovski@isye.gatech.edu Linear Optimization Problem,

More information

Handout 8: Dealing with Data Uncertainty

Handout 8: Dealing with Data Uncertainty MFE 5100: Optimization 2015 16 First Term Handout 8: Dealing with Data Uncertainty Instructor: Anthony Man Cho So December 1, 2015 1 Introduction Conic linear programming CLP, and in particular, semidefinite

More information

Strong Duality in Robust Semi-Definite Linear Programming under Data Uncertainty

Strong Duality in Robust Semi-Definite Linear Programming under Data Uncertainty Strong Duality in Robust Semi-Definite Linear Programming under Data Uncertainty V. Jeyakumar and G. Y. Li March 1, 2012 Abstract This paper develops the deterministic approach to duality for semi-definite

More information

IAC-11-C1.5.9 INERTIA-FREE ATTITUDE CONTROL OF SPACECRAFT WITH UNKNOWN TIME-VARYING MASS DISTRIBUTION

IAC-11-C1.5.9 INERTIA-FREE ATTITUDE CONTROL OF SPACECRAFT WITH UNKNOWN TIME-VARYING MASS DISTRIBUTION 6nd International Astronautical Congress, Cape Town, SA. Copyright by the International Astronautical Federation. All rights reserved IAC--C.5.9 INERTIA-FREE ATTITUDE CONTROL OF SPACECRAFT WITH UNKNOWN

More information

Filter Design for Linear Time Delay Systems

Filter Design for Linear Time Delay Systems IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 49, NO. 11, NOVEMBER 2001 2839 ANewH Filter Design for Linear Time Delay Systems E. Fridman Uri Shaked, Fellow, IEEE Abstract A new delay-dependent filtering

More information

Attitude Determination using Infrared Earth Horizon Sensors

Attitude Determination using Infrared Earth Horizon Sensors Attitude Determination using Infrared Earth Horizon Sensors Tam N. T. Nguyen Department of Aeronautics and Astronautics Massachusetts Institute of Technology 28 th Annual AIAA/USU Conference on Small Satellites

More information

Angular Velocity Determination Directly from Star Tracker Measurements

Angular Velocity Determination Directly from Star Tracker Measurements Angular Velocity Determination Directly from Star Tracker Measurements John L. Crassidis Introduction Star trackers are increasingly used on modern day spacecraft. With the rapid advancement of imaging

More information

On the Power of Robust Solutions in Two-Stage Stochastic and Adaptive Optimization Problems

On the Power of Robust Solutions in Two-Stage Stochastic and Adaptive Optimization Problems MATHEMATICS OF OPERATIONS RESEARCH Vol. 35, No., May 010, pp. 84 305 issn 0364-765X eissn 156-5471 10 350 084 informs doi 10.187/moor.1090.0440 010 INFORMS On the Power of Robust Solutions in Two-Stage

More information

Multiplicative vs. Additive Filtering for Spacecraft Attitude Determination

Multiplicative vs. Additive Filtering for Spacecraft Attitude Determination Multiplicative vs. Additive Filtering for Spacecraft Attitude Determination F. Landis Markley, NASA s Goddard Space Flight Center, Greenbelt, MD, USA Abstract The absence of a globally nonsingular three-parameter

More information

In-orbit magnetometer bias and scale factor calibration

In-orbit magnetometer bias and scale factor calibration Int. J. Metrol. Qual. Eng. 7, 104 (2016) c EDP Sciences 2016 DOI: 10.1051/ijmqe/2016003 In-orbit magnetometer bias and scale factor calibration Chingiz Hajiyev Faculty of Aeronautics and Astronautics,

More information

Global Geomagnetic Field Models from DMSP Satellite Magnetic Measurements

Global Geomagnetic Field Models from DMSP Satellite Magnetic Measurements Global Geomagnetic Field Models from DMSP Satellite Magnetic Measurements Patrick Alken Stefan Maus Arnaud Chulliat Manoj Nair Adam Woods National Geophysical Data Center, NOAA, Boulder, CO, USA 9 May

More information

Eigenstructure Assignment for Helicopter Hover Control

Eigenstructure Assignment for Helicopter Hover Control Proceedings of the 17th World Congress The International Federation of Automatic Control Eigenstructure Assignment for Helicopter Hover Control Andrew Pomfret Stuart Griffin Tim Clarke Department of Electronics,

More information

Simultaneous Adaptation of the Process and Measurement Noise Covariances for the UKF Applied to Nanosatellite Attitude Estimation

Simultaneous Adaptation of the Process and Measurement Noise Covariances for the UKF Applied to Nanosatellite Attitude Estimation Preprints of the 9th World Congress The International Federation of Automatic Control Simultaneous Adaptation of the Process and Measurement Noise Covariances for the UKF Applied to Nanosatellite Attitude

More information

A Hierarchy of Suboptimal Policies for the Multi-period, Multi-echelon, Robust Inventory Problem

A Hierarchy of Suboptimal Policies for the Multi-period, Multi-echelon, Robust Inventory Problem A Hierarchy of Suboptimal Policies for the Multi-period, Multi-echelon, Robust Inventory Problem Dimitris J. Bertsimas Dan A. Iancu Pablo A. Parrilo Sloan School of Management and Operations Research Center,

More information

On the Power of Robust Solutions in Two-Stage Stochastic and Adaptive Optimization Problems

On the Power of Robust Solutions in Two-Stage Stochastic and Adaptive Optimization Problems MATHEMATICS OF OPERATIONS RESEARCH Vol. xx, No. x, Xxxxxxx 00x, pp. xxx xxx ISSN 0364-765X EISSN 156-5471 0x xx0x 0xxx informs DOI 10.187/moor.xxxx.xxxx c 00x INFORMS On the Power of Robust Solutions in

More information

Predictive Attitude Estimation Using Global Positioning System Signals

Predictive Attitude Estimation Using Global Positioning System Signals Predictive Attitude Estimation Using Global Positioning System Signals John L. Crassidis Department of Mechanical Engineering he Catholic University of America Washington, DC 64 F. Landis Markley, E. Glenn

More information

Convex Optimization M2

Convex Optimization M2 Convex Optimization M2 Lecture 8 A. d Aspremont. Convex Optimization M2. 1/57 Applications A. d Aspremont. Convex Optimization M2. 2/57 Outline Geometrical problems Approximation problems Combinatorial

More information

OPTIMAL CONTROL AND ESTIMATION

OPTIMAL CONTROL AND ESTIMATION OPTIMAL CONTROL AND ESTIMATION Robert F. Stengel Department of Mechanical and Aerospace Engineering Princeton University, Princeton, New Jersey DOVER PUBLICATIONS, INC. New York CONTENTS 1. INTRODUCTION

More information

Attitude Estimation Version 1.0

Attitude Estimation Version 1.0 Attitude Estimation Version 1. Francesco Farina May 23, 216 Contents 1 Introduction 2 2 Mathematical background 2 2.1 Reference frames and coordinate systems............. 2 2.2 Euler angles..............................

More information

OPTIMAL ESTIMATION of DYNAMIC SYSTEMS

OPTIMAL ESTIMATION of DYNAMIC SYSTEMS CHAPMAN & HALL/CRC APPLIED MATHEMATICS -. AND NONLINEAR SCIENCE SERIES OPTIMAL ESTIMATION of DYNAMIC SYSTEMS John L Crassidis and John L. Junkins CHAPMAN & HALL/CRC A CRC Press Company Boca Raton London

More information

1 Robust optimization

1 Robust optimization ORF 523 Lecture 16 Princeton University Instructor: A.A. Ahmadi Scribe: G. Hall Any typos should be emailed to a a a@princeton.edu. In this lecture, we give a brief introduction to robust optimization

More information

Attitude Control Simulator for the Small Satellite and Its Validation by On-orbit Data of QSAT-EOS

Attitude Control Simulator for the Small Satellite and Its Validation by On-orbit Data of QSAT-EOS SSC17-P1-17 Attitude Control Simulator for the Small Satellite and Its Validation by On-orbit Data of QSAT-EOS Masayuki Katayama, Yuta Suzaki Mitsubishi Precision Company Limited 345 Kamikmachiya, Kamakura

More information

Satellite Attitude Determination with Attitude Sensors and Gyros using Steady-state Kalman Filter

Satellite Attitude Determination with Attitude Sensors and Gyros using Steady-state Kalman Filter Satellite Attitude Determination with Attitude Sensors and Gyros using Steady-state Kalman Filter Vaibhav V. Unhelkar, Hari B. Hablani Student, email: v.unhelkar@iitb.ac.in. Professor, email: hbhablani@aero.iitb.ac.in

More information

Riccati difference equations to non linear extended Kalman filter constraints

Riccati difference equations to non linear extended Kalman filter constraints International Journal of Scientific & Engineering Research Volume 3, Issue 12, December-2012 1 Riccati difference equations to non linear extended Kalman filter constraints Abstract Elizabeth.S 1 & Jothilakshmi.R

More information

Static Output Feedback Stabilisation with H Performance for a Class of Plants

Static Output Feedback Stabilisation with H Performance for a Class of Plants Static Output Feedback Stabilisation with H Performance for a Class of Plants E. Prempain and I. Postlethwaite Control and Instrumentation Research, Department of Engineering, University of Leicester,

More information

Distributionally robust optimization techniques in batch bayesian optimisation

Distributionally robust optimization techniques in batch bayesian optimisation Distributionally robust optimization techniques in batch bayesian optimisation Nikitas Rontsis June 13, 2016 1 Introduction This report is concerned with performing batch bayesian optimization of an unknown

More information

Optimization of Orbital Transfer of Electrodynamic Tether Satellite by Nonlinear Programming

Optimization of Orbital Transfer of Electrodynamic Tether Satellite by Nonlinear Programming Optimization of Orbital Transfer of Electrodynamic Tether Satellite by Nonlinear Programming IEPC-2015-299 /ISTS-2015-b-299 Presented at Joint Conference of 30th International Symposium on Space Technology

More information

Simplified Filtering Estimator for Spacecraft Attitude Determination from Phase Information of GPS Signals

Simplified Filtering Estimator for Spacecraft Attitude Determination from Phase Information of GPS Signals WCE 7, July - 4, 7, London, U.K. Simplified Filtering Estimator for Spacecraft Attitude Determination from Phase Information of GPS Signals S. Purivigraipong, Y. Hashida, and M. Unwin Abstract his paper

More information

Influence Analysis of Star Sensors Sampling Frequency on Attitude Determination Accuracy

Influence Analysis of Star Sensors Sampling Frequency on Attitude Determination Accuracy Sensors & ransducers Vol. Special Issue June pp. -8 Sensors & ransducers by IFSA http://www.sensorsportal.com Influence Analysis of Star Sensors Sampling Frequency on Attitude Determination Accuracy Yuanyuan

More information

Analytical Mechanics. of Space Systems. tfa AA. Hanspeter Schaub. College Station, Texas. University of Colorado Boulder, Colorado.

Analytical Mechanics. of Space Systems. tfa AA. Hanspeter Schaub. College Station, Texas. University of Colorado Boulder, Colorado. Analytical Mechanics of Space Systems Third Edition Hanspeter Schaub University of Colorado Boulder, Colorado John L. Junkins Texas A&M University College Station, Texas AIM EDUCATION SERIES Joseph A.

More information

Robust Dual-Response Optimization

Robust Dual-Response Optimization Yanıkoğlu, den Hertog, and Kleijnen Robust Dual-Response Optimization 29 May 1 June 1 / 24 Robust Dual-Response Optimization İhsan Yanıkoğlu, Dick den Hertog, Jack P.C. Kleijnen Özyeğin University, İstanbul,

More information

Multiobjective H 2 /H /impulse-to-peak synthesis: Application to the control of an aerospace launcher

Multiobjective H 2 /H /impulse-to-peak synthesis: Application to the control of an aerospace launcher Multiobjective H /H /impulse-to-peak synthesis: Application to the control of an aerospace launcher D. Arzelier, D. Peaucelle LAAS-CNRS, 7 Avenue du Colonel Roche, 3 77 Toulouse, Cedex 4, France emails:

More information

Handout 6: Some Applications of Conic Linear Programming

Handout 6: Some Applications of Conic Linear Programming ENGG 550: Foundations of Optimization 08 9 First Term Handout 6: Some Applications of Conic Linear Programming Instructor: Anthony Man Cho So November, 08 Introduction Conic linear programming CLP, and

More information

arxiv:math/ v1 [math.oc] 17 Sep 2006

arxiv:math/ v1 [math.oc] 17 Sep 2006 Global Attitude Estimation using Single Direction Measurements Taeyoung Lee, Melvin Leok, N. Harris McClamroch, and Amit Sanyal arxiv:math/0609481v1 [math.oc] 17 Sep 2006 Abstract A deterministic attitude

More information

Efficient robust optimization for robust control with constraints Paul Goulart, Eric Kerrigan and Danny Ralph

Efficient robust optimization for robust control with constraints Paul Goulart, Eric Kerrigan and Danny Ralph Efficient robust optimization for robust control with constraints p. 1 Efficient robust optimization for robust control with constraints Paul Goulart, Eric Kerrigan and Danny Ralph Efficient robust optimization

More information

Distributionally Robust Optimization with ROME (part 1)

Distributionally Robust Optimization with ROME (part 1) Distributionally Robust Optimization with ROME (part 1) Joel Goh Melvyn Sim Department of Decision Sciences NUS Business School, Singapore 18 Jun 2009 NUS Business School Guest Lecture J. Goh, M. Sim (NUS)

More information

A Hierarchy of Polyhedral Approximations of Robust Semidefinite Programs

A Hierarchy of Polyhedral Approximations of Robust Semidefinite Programs A Hierarchy of Polyhedral Approximations of Robust Semidefinite Programs Raphael Louca Eilyan Bitar Abstract Robust semidefinite programs are NP-hard in general In contrast, robust linear programs admit

More information

The Generalized Wahba Problem

The Generalized Wahba Problem The Journal of the Astronautical Sciences, Vol. 54, No., April June 006, pp. 45 59 The Generalized Wahba Problem Malcolm D. Shuster To generalize is to be an idiot. William Blake (757 87) People who like

More information

Extension of Farrenkopf Steady-State Solutions with Estimated Angular Rate

Extension of Farrenkopf Steady-State Solutions with Estimated Angular Rate Extension of Farrenopf Steady-State Solutions with Estimated Angular Rate Andrew D. Dianetti and John L. Crassidis University at Buffalo, State University of New Yor, Amherst, NY 46-44 Steady-state solutions

More information

Robust Fisher Discriminant Analysis

Robust Fisher Discriminant Analysis Robust Fisher Discriminant Analysis Seung-Jean Kim Alessandro Magnani Stephen P. Boyd Information Systems Laboratory Electrical Engineering Department, Stanford University Stanford, CA 94305-9510 sjkim@stanford.edu

More information

Automated Tuning of the Nonlinear Complementary Filter for an Attitude Heading Reference Observer

Automated Tuning of the Nonlinear Complementary Filter for an Attitude Heading Reference Observer Automated Tuning of the Nonlinear Complementary Filter for an Attitude Heading Reference Observer Oscar De Silva, George K.I. Mann and Raymond G. Gosine Faculty of Engineering and Applied Sciences, Memorial

More information

Here represents the impulse (or delta) function. is an diagonal matrix of intensities, and is an diagonal matrix of intensities.

Here represents the impulse (or delta) function. is an diagonal matrix of intensities, and is an diagonal matrix of intensities. 19 KALMAN FILTER 19.1 Introduction In the previous section, we derived the linear quadratic regulator as an optimal solution for the fullstate feedback control problem. The inherent assumption was that

More information

Robust linear optimization under general norms

Robust linear optimization under general norms Operations Research Letters 3 (004) 50 56 Operations Research Letters www.elsevier.com/locate/dsw Robust linear optimization under general norms Dimitris Bertsimas a; ;, Dessislava Pachamanova b, Melvyn

More information

2nd Symposium on System, Structure and Control, Oaxaca, 2004

2nd Symposium on System, Structure and Control, Oaxaca, 2004 263 2nd Symposium on System, Structure and Control, Oaxaca, 2004 A PROJECTIVE ALGORITHM FOR STATIC OUTPUT FEEDBACK STABILIZATION Kaiyang Yang, Robert Orsi and John B. Moore Department of Systems Engineering,

More information

3D 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 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 information

Canonical Problem Forms. Ryan Tibshirani Convex Optimization

Canonical Problem Forms. Ryan Tibshirani Convex Optimization Canonical Problem Forms Ryan Tibshirani Convex Optimization 10-725 Last time: optimization basics Optimization terology (e.g., criterion, constraints, feasible points, solutions) Properties and first-order

More information

Linear Feedback Control Using Quasi Velocities

Linear 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 information

REACTION WHEEL CONFIGURATIONS FOR HIGH AND MIDDLE INCLINATION ORBITS

REACTION WHEEL CONFIGURATIONS FOR HIGH AND MIDDLE INCLINATION ORBITS REACTION WHEEL CONFIGURATIONS FOR HIGH AND MIDDLE INCLINATION ORBITS Zuliana Ismail and Renuganth Varatharajoo Department of Aerospace Engineering, Universiti Putra Malaysia, Malaysia E-Mail: zuliana.ismail@gmail.com

More information

Stability of linear time-varying systems through quadratically parameter-dependent Lyapunov functions

Stability of linear time-varying systems through quadratically parameter-dependent Lyapunov functions Stability of linear time-varying systems through quadratically parameter-dependent Lyapunov functions Vinícius F. Montagner Department of Telematics Pedro L. D. Peres School of Electrical and Computer

More information

SYNTHESIS OF ROBUST DISCRETE-TIME SYSTEMS BASED ON COMPARISON WITH STOCHASTIC MODEL 1. P. V. Pakshin, S. G. Soloviev

SYNTHESIS OF ROBUST DISCRETE-TIME SYSTEMS BASED ON COMPARISON WITH STOCHASTIC MODEL 1. P. V. Pakshin, S. G. Soloviev SYNTHESIS OF ROBUST DISCRETE-TIME SYSTEMS BASED ON COMPARISON WITH STOCHASTIC MODEL 1 P. V. Pakshin, S. G. Soloviev Nizhny Novgorod State Technical University at Arzamas, 19, Kalinina ul., Arzamas, 607227,

More information

Graph and Controller Design for Disturbance Attenuation in Consensus Networks

Graph and Controller Design for Disturbance Attenuation in Consensus Networks 203 3th International Conference on Control, Automation and Systems (ICCAS 203) Oct. 20-23, 203 in Kimdaejung Convention Center, Gwangju, Korea Graph and Controller Design for Disturbance Attenuation in

More information

Least Squares with Examples in Signal Processing 1. 2 Overdetermined equations. 1 Notation. The sum of squares of x is denoted by x 2 2, i.e.

Least Squares with Examples in Signal Processing 1. 2 Overdetermined equations. 1 Notation. The sum of squares of x is denoted by x 2 2, i.e. Least Squares with Eamples in Signal Processing Ivan Selesnick March 7, 3 NYU-Poly These notes address (approimate) solutions to linear equations by least squares We deal with the easy case wherein the

More information

ATTITUDE CONTROL MECHANIZATION TO DE-ORBIT SATELLITES USING SOLAR SAILS

ATTITUDE CONTROL MECHANIZATION TO DE-ORBIT SATELLITES USING SOLAR SAILS IAA-AAS-DyCoSS2-14-07-02 ATTITUDE CONTROL MECHANIZATION TO DE-ORBIT SATELLITES USING SOLAR SAILS Ozan Tekinalp, * Omer Atas INTRODUCTION Utilization of solar sails for the de-orbiting of satellites is

More information

Research Article An Equivalent LMI Representation of Bounded Real Lemma for Continuous-Time Systems

Research Article An Equivalent LMI Representation of Bounded Real Lemma for Continuous-Time Systems Hindawi Publishing Corporation Journal of Inequalities and Applications Volume 28, Article ID 67295, 8 pages doi:1.1155/28/67295 Research Article An Equivalent LMI Representation of Bounded Real Lemma

More information

Investigation of the Attitude Error Vector Reference Frame in the INS EKF

Investigation of the Attitude Error Vector Reference Frame in the INS EKF Investigation of the Attitude Error Vector Reference Frame in the INS EKF Stephen Steffes, Jan Philipp Steinbach, and Stephan Theil Abstract The Extended Kalman Filter is used extensively for inertial

More information

Robust portfolio selection under norm uncertainty

Robust portfolio selection under norm uncertainty Wang and Cheng Journal of Inequalities and Applications (2016) 2016:164 DOI 10.1186/s13660-016-1102-4 R E S E A R C H Open Access Robust portfolio selection under norm uncertainty Lei Wang 1 and Xi Cheng

More information

H Estimation. Speaker : R.Lakshminarayanan Guide : Prof. K.Giridhar. H Estimation p.1/34

H Estimation. Speaker : R.Lakshminarayanan Guide : Prof. K.Giridhar. H Estimation p.1/34 H Estimation Speaker : R.Lakshminarayanan Guide : Prof. K.Giridhar H Estimation p.1/34 H Motivation The Kalman and Wiener Filters minimize the mean squared error between the true value and estimated values

More information

Largest dual ellipsoids inscribed in dual cones

Largest dual ellipsoids inscribed in dual cones Largest dual ellipsoids inscribed in dual cones M. J. Todd June 23, 2005 Abstract Suppose x and s lie in the interiors of a cone K and its dual K respectively. We seek dual ellipsoidal norms such that

More information

Feedback Control of Spacecraft Rendezvous Maneuvers using Differential Drag

Feedback Control of Spacecraft Rendezvous Maneuvers using Differential Drag Feedback Control of Spacecraft Rendezvous Maneuvers using Differential Drag D. Pérez 1 and R. Bevilacqua Rensselaer Polytechnic Institute, Troy, New York, 1180 This work presents a feedback control strategy

More information

SATELLITE ORBIT ESTIMATION USING EARTH MAGNETIC FIELD MEASUREMENTS

SATELLITE ORBIT ESTIMATION USING EARTH MAGNETIC FIELD MEASUREMENTS International Journal of Engineering and echnology, Vol. 3, No., 6, pp. 63-71 63 SAELLIE ORBI ESIMAION USING EARH MAGNEIC FIELD MEASUREMENS Mohammad Nizam Filipsi, Renuganth Varatharajoo Department of

More information

Deterministic Relative Attitude Determination of Formation Flying Spacecraft

Deterministic Relative Attitude Determination of Formation Flying Spacecraft Deterministic Relative Attitude Determination of Formation Flying Spacecraft Michael S. Andrle, Baro Hyun, John L. Crassidis, and Richard Linares University at Buffalo, State University of New York, Amherst,

More information

Convex Optimization. (EE227A: UC Berkeley) Lecture 6. Suvrit Sra. (Conic optimization) 07 Feb, 2013

Convex Optimization. (EE227A: UC Berkeley) Lecture 6. Suvrit Sra. (Conic optimization) 07 Feb, 2013 Convex Optimization (EE227A: UC Berkeley) Lecture 6 (Conic optimization) 07 Feb, 2013 Suvrit Sra Organizational Info Quiz coming up on 19th Feb. Project teams by 19th Feb Good if you can mix your research

More information

15 x 104 Q (5,6) (t) and Q (6,6) (t)

15 x 104 Q (5,6) (t) and Q (6,6) (t) PERIODIC H SYNTHESIS FOR SPACECRAFT ATTITUDE DETERMINATION AND CONTROL WITH A VECTOR MAGNETOMETER AND MAGNETORQUERS Rafaψl Wiśniewski Λ Jakob Stoustrup Λ Λ Department of Control Engineering, Aalborg University,

More information

NEW EUMETSAT POLAR SYSTEM ATTITUDE MONITORING SOFTWARE

NEW EUMETSAT POLAR SYSTEM ATTITUDE MONITORING SOFTWARE NEW EUMETSAT POLAR SYSTEM ATTITUDE MONITORING SOFTWARE Pablo García Sánchez (1), Antonio Pérez Cambriles (2), Jorge Eufrásio (3), Pier Luigi Righetti (4) (1) GMV Aerospace and Defence, S.A.U., Email: pgarcia@gmv.com,

More information

Operations Research Letters

Operations Research Letters Operations Research Letters 37 (2009) 1 6 Contents lists available at ScienceDirect Operations Research Letters journal homepage: www.elsevier.com/locate/orl Duality in robust optimization: Primal worst

More information

Attitude Determination for NPS Three-Axis Spacecraft Simulator

Attitude Determination for NPS Three-Axis Spacecraft Simulator AIAA/AAS Astrodynamics Specialist Conference and Exhibit 6-9 August 4, Providence, Rhode Island AIAA 4-5386 Attitude Determination for NPS Three-Axis Spacecraft Simulator Jong-Woo Kim, Roberto Cristi and

More information

Interval solutions for interval algebraic equations

Interval solutions for interval algebraic equations Mathematics and Computers in Simulation 66 (2004) 207 217 Interval solutions for interval algebraic equations B.T. Polyak, S.A. Nazin Institute of Control Sciences, Russian Academy of Sciences, 65 Profsoyuznaya

More information

Infrared Earth Horizon Sensors for CubeSat Attitude Determination

Infrared Earth Horizon Sensors for CubeSat Attitude Determination Infrared Earth Horizon Sensors for CubeSat Attitude Determination Tam Nguyen Department of Aeronautics and Astronautics Massachusetts Institute of Technology Outline Background and objectives Nadir vector

More information

Further Results on Model Structure Validation for Closed Loop System Identification

Further Results on Model Structure Validation for Closed Loop System Identification Advances in Wireless Communications and etworks 7; 3(5: 57-66 http://www.sciencepublishinggroup.com/j/awcn doi:.648/j.awcn.735. Further esults on Model Structure Validation for Closed Loop System Identification

More information

Advances in Convex Optimization: Theory, Algorithms, and Applications

Advances in Convex Optimization: Theory, Algorithms, and Applications Advances in Convex Optimization: Theory, Algorithms, and Applications Stephen Boyd Electrical Engineering Department Stanford University (joint work with Lieven Vandenberghe, UCLA) ISIT 02 ISIT 02 Lausanne

More information

On construction of constrained optimum designs

On construction of constrained optimum designs On construction of constrained optimum designs Institute of Control and Computation Engineering University of Zielona Góra, Poland DEMA2008, Cambridge, 15 August 2008 Numerical algorithms to construct

More information

Auxiliary signal design for failure detection in uncertain systems

Auxiliary signal design for failure detection in uncertain systems Auxiliary signal design for failure detection in uncertain systems R. Nikoukhah, S. L. Campbell and F. Delebecque Abstract An auxiliary signal is an input signal that enhances the identifiability of a

More information

Quaternion-Based Tracking Control Law Design For Tracking Mode

Quaternion-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 information

EE5138R: Problem Set 5 Assigned: 16/02/15 Due: 06/03/15

EE5138R: Problem Set 5 Assigned: 16/02/15 Due: 06/03/15 EE538R: Problem Set 5 Assigned: 6/0/5 Due: 06/03/5. BV Problem 4. The feasible set is the convex hull of (0, ), (0, ), (/5, /5), (, 0) and (, 0). (a) x = (/5, /5) (b) Unbounded below (c) X opt = {(0, x

More information

A NONLINEARITY MEASURE FOR ESTIMATION SYSTEMS

A NONLINEARITY MEASURE FOR ESTIMATION SYSTEMS AAS 6-135 A NONLINEARITY MEASURE FOR ESTIMATION SYSTEMS Andrew J. Sinclair,JohnE.Hurtado, and John L. Junkins The concept of nonlinearity measures for dynamical systems is extended to estimation systems,

More information

H 2 Adaptive Control. Tansel Yucelen, Anthony J. Calise, and Rajeev Chandramohan. WeA03.4

H 2 Adaptive Control. Tansel Yucelen, Anthony J. Calise, and Rajeev Chandramohan. WeA03.4 1 American Control Conference Marriott Waterfront, Baltimore, MD, USA June 3-July, 1 WeA3. H Adaptive Control Tansel Yucelen, Anthony J. Calise, and Rajeev Chandramohan Abstract Model reference adaptive

More information

ADAPTIVE FILTER THEORY

ADAPTIVE FILTER THEORY ADAPTIVE FILTER THEORY Fifth Edition Simon Haykin Communications Research Laboratory McMaster University Hamilton, Ontario, Canada International Edition contributions by Telagarapu Prabhakar Department

More information

Robust Optimization for Risk Control in Enterprise-wide Optimization

Robust Optimization for Risk Control in Enterprise-wide Optimization Robust Optimization for Risk Control in Enterprise-wide Optimization Juan Pablo Vielma Department of Industrial Engineering University of Pittsburgh EWO Seminar, 011 Pittsburgh, PA Uncertainty in Optimization

More information