ORBIT DETERMINATION OF LEO SATELLITES FOR A SINGLE PASS THROUGH A RADAR: COMPARISON OF METHODS
|
|
- Abraham Richards
- 5 years ago
- Views:
Transcription
1 ORBIT DETERMINATION OF LEO SATELLITES FOR A SINGLE PASS THROUGH A RADAR: COMPARISON OF METHODS Zakhary N. Khutorovsky, Sergey Yu. Kamensky and, Nikolay N. Sbytov Vympel Corporation, Moscow, Russia Kyle T. Alfriend Texas A&M University, USA
2 OUTLINE Introduction Least squares algorithm Recurrent algorithm Guarantee approach method
3 Introduction Initial Orbit Determination (IOD) Russia radar sites track each object and report the state vector and covariance matrix as the observation instead of observations. Current IOD approaches were established in the 1960s and 70s. Very little interest until recently. The planned radar fence upgrade and WFOV cameras will grow the catalog from the current objects to an estimated objects. Correlating tracks from these objects will be a major challenge. The accuracy of collisiom probability determination depend on the accuracy state and covariances. Need improved state estimates and covariances from the IOD.
4 The least squares algorithm The least squares method search for the minimum of ( ) Ψ(a) = n 1 (d σd 2 k d k (a)) (α σα 2 k α k (a)) (β k k σβ 2 k β k (a)) 2, k where k=1 d k,α k,β k k-th measurement range, azimuth, elevation (k =1, 2,..., n); t k epochofthek-thmeasurement; σ dk,σ αk,σ βk RMS of the errors d k,α k,β k ; a = a ( t) - state at epoch t. h k (a) =(d k (a), α k (a), β k (a)) values of the parameters of the k-th measurement, calculated using the vector a. ( ) Ψ The solution a min is obtained from a (a min) ( ) 1 2 Ψ and covariance K from K =0.5 a (a min). 2 =0
5 Recurrent algorithm The estimate a min of the least squares method has recurrent structure: the estimate of parameters based on k measurements can be written as a function of the same estimate for k 1 measurements and the k-th measurement. The recurrence feature of the least squares estimate can be extanded to more general case, when the parameters a k =a(t k ) of the system for the time t k (in our case the orbital parameters for the time t k ) in addition to deterministic component have a random one which has Markov s feature (the future for the fixed past depends only on the present ), i.e. p (a k a k 1, a k 2,...) =p (a k a k 1 ) (M) where p (a k...) - conditional probability density a k. The condition (M) is satisfied if a k = f k (a k 1 ) + γ k (M1) where γ k - noise of the system - time independent random perturbations with zero mean and covariation matrices Γ k.
6 For (M), (M1) and independent errors of the measurements for the aposteriori (after acquisition of the measurement x k ) probability density of the parameters a k we have recurrent formulas for density p the first and second moments a, P p (a k x k, a k 1, a k 2,...)=p (a k x k, a k 1 )=p (x k a k ) p (a k, a k 1 )/p (x k ) P 1 k = P 1 + k k 1 H k R 1H k k P 1 k âk = P 1 k k 1âk k 1 + H k R 1x k k, â k k 1 =f k (â k 1 ) F k = f k(â k 1 ) a k 1 H k = h k(â k k 1 ) a k P k k 1 =F k P k 1 F k+γ k If measured parameter is a scalar u w k = w k /(1 + w k h kp k k 1 h k ) P k = P k k 1 w k (P k k 1 h k )(P k k 1 h k ) â k = â k k 1 + w k P k h k (u k h k (â k k 1 ) where w k = 1/σ 2 u k is the weight characteristic of measurement u k.
7 Least Squares and Recurrent Issues Selection of the state a, the dynamic model and the propagation method for state and covariances. Calculation of the initial approximaiton recurrent algorithm, â 0, P 0 for the The technique for calculating the noise matrix of the system Γ k for the recurrent algorithm. Calculation of the initial approximaiton a 0, for minimization of Ψ(a); Selection of technique for reaching the minimum of Ψ(a).
8 The state vector a a =(x, y, z, ẋ, ẏ, ż) for the certain time in the local rectangular coordinate frame (LRCF) related to geographical directions. The basic plane is the plane of local horizon, the x axis is in the plane of local horizon and is directed to the west, the y axis is directed normal to the plane of local horizon upwards ad z axis is in the plane of local horizon and is directed to the north. x = d sin α cos β y = d sin β z = d cos α cos β
9 The dynamic model d(t) d = µ r 3 r 2 ω d ω ( ω r) 3J 2µR 2 e 2r 5 ( r+2( r, ω 0 ) ω 0 5( r, ω ) 0) 2 r r 2 d =(x, y, z), d =(ẋ, ẏ, ż) r =(x b x,y b y,z b z ) Greenwich coordinats of SO r = (x b x ) 2 +(y b y ) 2 +(z b z ) 2 r 1 = X 2 g0 +Y 2 g0 r 0 = X 2 g0 +Y 2 g0 +Z2 g0 z = Z g0 /r 0 r 1 = r 1 /r 0 R = R e (1 α z 2 ) b x =0 b y = r 0 b z =2 α R z r 1 ω = ω e (ω 0x,ω 0y,ω 0z ) ω e = /s ω 0x =0 ω 0y = z + r 1 b z /r 0 ω 0z = r 1 z b z /r 0 α= µ= km 3 /s 2 J 2 = R e = km X g0, Y g0, Z g0 Greenwich coordinats of LRCF origin
10 The propagation method for state and covariances. For solving the system of differential equations d = f( d, d) use the numerical Runge-Cutta methos of the 4th order. a i =(k (i) 1 +2k (i) 2 +2k (i) 3 + k (i) 4 )/6 k (i) 1 =τ i f(a i ) k (i) 2 =τ i f(a i +0.5k (i) 1 ) k (i) 3 =τ i f(a i +0.5k (i) 2 ) k (i) 4 =τ i f(a i +k (i) 3 ) For the operator F k (a k 1 ) we will take the matrix F k with dimensions 6 6 with the following non-zero elements: f i,i =1 for i =1, 2,..., 6, f i,i+3 =τ k for i =1, 2, 3, where τ k = t k t k 1. F n = τ n τ n τ n
11 Calculation of the initial approximaiton â 0, P 0 for the recurrent algorithm, The state and covariance estimate is determined from the 1st two observations d 1 and d 2. The technique for calculating the noise matrix of the system Γ k for the recurrent algorithm. Γ 11,k = Γ k = ( Γ 11,k where 110 σd 2/(τ n 2 ( t ef,d) 4 ) σα 2 /(τ n 2 ( t ef,α) 4 ) σβ 2 /(τ n 2 ( t ef,β) 4 ) ) σ d,σ α,σ β are the RMS values of typical values of the errors of measured parameters, t ef,d, t ef,α, t ef,β efficient memory of the algorithm for parameters d, α, β (parameters of the algorithm selected on experimental basis).
12 Calculation of the initial approximaiton a 0, for minimization of Ψ(a); For the initial approximation a 0 we use the estimate â n, obtained using measurements x 1, x 2,..., x n by the recurrent algorithm. Selection of technique for reaching the minimum of Ψ(a). The major effect of updating by least squares is the better accuracy of the velocity components u=(ẋ, ẏ, ż) of the state vector a. The position parameters u=(x, y, z) of the vector a in fact are not updated. Thus we suggest to make minimization of Ψ(a) =Ψ(u, u) =ψ( u) only by the vector u. For the point of the minimum u min of the function ψ( u) for the case when it is reached by one iteration we suggest the Newton s formula ( ) 1 2 ψ u min = u 0 u 2(a 0) ψ u (a 0) where the first and the second derivatives of the function ψ( u) with respect to parameters ẋ, ẏ, ż are calculated by finite differences method.
13 New (Guarantee) Approach Fundamental problem with the traditional methods (joint and recurrent) is the assumption of uncorrelated measurements and known measurement statistics. Because there is no a priori orbit identifying bad measurements can be difficult. Basis of the Guarantee method is a guaranteed range of the measurement errors. Guarantee means that the algorithm produces not only the orbit estimate but the maximum possible errors in the estimated parameters. Leads to a more accurate estimate of the orbit and uncertainty. Provides an upper bound on the state errors. This feature hin- More computationally intensive. dered early use.
14 Let for the moments t k (t k t k+1 ; k = 1, 2,...n) we acquire measurements u k of certain functions h k (c) of the m-vector of parameters c (m<n), and the errors of the measurements δu k =u k h k (c) are limited from above by constant δ k,max. The limits for errors of the measurements define in the m-dimensional space of parameters c =(c 1, c 2,..., c m ) a domain of possible values of c D n = n {u k δ i,max h k (c) u k δ k,max } k=1 Project this domain to the coordinate of the components of vector c and among the projected points for each axis find the most right c n,r and the most left c n,l. They define the boundaries (maximum and minimum values) for the changes of each component of parameter c. In this case the estimate c n of parameter and maximum errors δ c n,max of this estimate are naturally defined as c n = 1 2 (c n,r + c n,l ) δ c n,max = 1 2 (c n,r c n,l ) When the measured parameters are linearly related to the determined parameter a, finding c n and δ c n,max is a linear programming problem.
15 The Model Example m =1 h(c) =c u k = c + δ k, δ k,max = δ max. D n = n {[ u k δ max,u k +δ max ]} = [ max k=1 c n = 1 2 (max k u k +min k u k ) k u k δ max, min k δ c n,max =δ max 1 2 (max k [ a, b ] - interval with left end a and right end b. u k +δ max ] u k min k u k ) c n does not depend on δ max If we know δ max precisely the value δ c n,max is a correct upper estimate for the error of the estimate c n, for any correlations of the errors of different measurements. For not correlated errors of measurements the estimate c n with regard to accuracy is not inferior and sometimes can be essentially more accurate (!) then the estimate n ĉ n = 1 n u k of the least squares technique. k=1 for uniform distribution: σ cn 1.4 δ max n for triangular distribution: σ cn 0.46 δ max n σĉn 0.58 δ max (1) n σĉn 0.41 δ max (2) n
16 The Real Algorithm Assume there are measurements range d k, azimuth α k and elevation β k at the times t k (k =1, 2,..., n). The error limits are defined be δ d,max,δ α,max,δ β,max. Let epoch be t=0.5 (t 1 +t n ) and orbital parameters are c =(d, α, β, d, α, β). Divide the measurements into n/2 groups. Each group defines two points. The k-th group is (t k, t 0.5n+k). From two points at two different times determine the state or orbit defined by č k =(ď, ˇα, ˇβ, d, ˇ ˇ α, ˇ β)k. The determination of the states č( t) assume 2-body motion. Now computer the measurements that would give this solution including J 2. Integrate the equation for the deviations of r, ṙ from Keplerian orbit with the initial conditions r ( t )= ṙ( t )=0 and r = r k to obtain the errors in the measurements r k (t k ). Then compute the corrected measurements r k =r k r k and transform the r k to the measurement frame. The domain D k is approximated by a six-dimensional parallelepiped with the center at č k and 64 apexes determined by č k ± δ d,max j 1 ± δ α,max j 2 ± δ β,max j 3 ± δ d,max j 4 ± δ α,max j 5 ± δ β,max j 6 where j i (i=1,2,...6) are the lines of the matrix of partial derivatives (c k ) J = (u k, u. The boundary projections c k,l and c k,r of these 0.5n+k) 64 points determine the boundaries of the vector interval (c k,l, c k,r). Then (c l, c r )= n/2 (c k,l, c k,r) c n = 1 2 (c n,r + c n,l ) δ c n,max = 1 2 (c n,r c n,l ) k=1
17 References 1. Bolshakov I.A., Repin V.G. Issues of non-linear filtering, Automatics and telemechanics, v.23, N.4, Shapiro I.I. The Prediction of Ballistic Missile Trajectories From Radar Observation, Massachusetts Institute of Technology Lincoln Laboratory, McGraw-Hill, New York-Toronto-London, Mudrov V.I., Method of least modules, Moscow, Kotov E.O., On the problem of determination of satellite orbit with limited errors of measurements, Space research, v. 19, issue 4, Moscow, 1981 (in Russian). 5. Cramer H., Methods of Mathematical Statistics, Princeton University Press. Princeton, New York, Z.N. Khutorovsky., Efficient memory of the discreet Kalman s filter, USSR Academy of Science, Automatics and Telemechanics, Nauka, Moscow, N.2, 1972, (in Russian). 7. Vallado D.A., Fundamentals of Astrodynamics and Applications, Second Edition, Microcosm Press,
18 2001.
19 SUMMARY Obtaining the best initial orbit is vital for correlating tracks and developing the future space catalog. Philosophy and approach US and Russian methods of initial orbit determination are different. Current Russian methods of initial orbit determination have been reviewed. New method called The Guarantee method has been developed. Provides a more accurate orbit estimate. Provides upper bound on estimated state error. More computationally intensive.
ORBIT DETERMINATION OF LEO SATELLITES FOR A SINGLE PASS THROUGH A RADAR: COMPARISON OF METHODS Z.
ORBIT DETERMINATION OF LEO SATELLITES FOR A SINGLE PASS THROUGH A RADAR: COMPARISON OF METHODS Z. Khutorovsy, S. Kamensy and, N. Sbytov Vympel Corporation, Moscow, Russia K. T. Alfriend Texas A&M University,
More informationA 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 informationSINGULAR VALUE DECOMPOSITION AND LEAST SQUARES ORBIT DETERMINATION
SINGULAR VALUE DECOMPOSITION AND LEAST SQUARES ORBIT DETERMINATION Zakhary N. Khutorovsky & Vladimir Boikov Vympel Corp. Kyle T. Alfriend Texas A&M University Outline Background Batch Least Squares Nonlinear
More informationA.I. Nazarenko 1
Model Study of the Possibilities of Space Debris Cataloging A.I. Nazarenko 1 anazarenko32@mail.ru The objective of this paper is the estimation of the characteristics of the measuring means and the software,
More informationCollision Prediction for LEO Satellites. Analysis of Characteristics
Collision Prediction for LEO Satellites. Analysis of Characteristics Viacheslav F. Fateev Doctor of Science (technical sciences), Professor, Russia, Vympel Corporation, President Sergey A. Sukhanov Doctor
More informationSpace 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 informationDevelopment of an algorithm for the problem of the least-squares method: Preliminary Numerical Experience
Development of an algorithm for the problem of the least-squares method: Preliminary Numerical Experience Sergey Yu. Kamensky 1, Vladimir F. Boykov 2, Zakhary N. Khutorovsky 3, Terry K. Alfriend 4 Abstract
More informationBenefits of a Geosynchronous Orbit (GEO) Observation Point for Orbit Determination
Benefits of a Geosynchronous Orbit (GEO) Observation Point for Orbit Determination Ray Byrne, Michael Griesmeyer, Ron Schmidt, Jeff Shaddix, and Dave Bodette Sandia National Laboratories ABSTRACT Determining
More informationKalman Filters with Uncompensated Biases
Kalman Filters with Uncompensated Biases Renato Zanetti he Charles Stark Draper Laboratory, Houston, exas, 77058 Robert H. Bishop Marquette University, Milwaukee, WI 53201 I. INRODUCION An underlying assumption
More informationPerformance of a Dynamic Algorithm For Processing Uncorrelated Tracks
Performance of a Dynamic Algorithm For Processing Uncorrelated Tracs Kyle T. Alfriend Jong-Il Lim Texas A&M University Tracs of space objects, which do not correlate, to a nown space object are called
More informationWhen Does the Uncertainty Become Non-Gaussian. Kyle T. Alfriend 1 Texas A&M University Inkwan Park 2 Texas A&M University
When Does the Uncertainty Become Non-Gaussian Kyle T. Alfriend Texas A&M University Inkwan Park 2 Texas A&M University ABSTRACT The orbit state covariance is used in the conjunction assessment/probability
More informationDynamic Model of Space Robot Manipulator
Applied Mathematical Sciences, Vol. 9, 215, no. 94, 465-4659 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.12988/ams.215.56429 Dynamic Model of Space Robot Manipulator Polina Efimova Saint-Petersburg
More informationLINEARIZED ORBIT COVARIANCE GENERATION AND PROPAGATION ANALYSIS VIA SIMPLE MONTE CARLO SIMULATIONS
LINEARIZED ORBIT COVARIANCE GENERATION AND PROPAGATION ANALYSIS VIA SIMPLE MONTE CARLO SIMULATIONS Chris Sabol, Paul Schumacher AFRL Thomas Sukut USAFA Terry Alfriend Texas A&M Keric Hill PDS Brendan Wright
More informationMiscellaneous. Regarding reading materials. Again, ask questions (if you have) and ask them earlier
Miscellaneous Regarding reading materials Reading materials will be provided as needed If no assigned reading, it means I think the material from class is sufficient Should be enough for you to do your
More informationRELATIVE NAVIGATION FOR SATELLITES IN CLOSE PROXIMITY USING ANGLES-ONLY OBSERVATIONS
(Preprint) AAS 12-202 RELATIVE NAVIGATION FOR SATELLITES IN CLOSE PROXIMITY USING ANGLES-ONLY OBSERVATIONS Hemanshu Patel 1, T. Alan Lovell 2, Ryan Russell 3, Andrew Sinclair 4 "Relative navigation using
More informationSpacecraft Orbit Anomaly Representation Using Thrust-Fourier-Coefficients with Orbit Determination Toolbox
Spacecraft Orbit Anomaly Representation Using Thrust-Fourier-Coefficients with Orbit Determination Toolbox Hyun Chul Ko and Daniel J. Scheeres University of Colorado - Boulder, Boulder, CO, USA ABSTRACT
More informationCataloging with an Upgraded Space Surveillance Fence
Cataloging with an Upgraded Space Surveillance Fence Felix Hoots, Geoff Pierce, Lester Ford AT&T Hugh Hadley Syracuse Research Integrated solutions from a trusted source Fence Provides Robust, Uncued Surveillance
More informationSTATISTICAL ORBIT DETERMINATION
STATISTICAL ORBIT DETERMINATION Satellite Tracking Example of SNC and DMC ASEN 5070 LECTURE 6 4.08.011 1 We will develop a simple state noise compensation (SNC) algorithm. This algorithm adds process noise
More informationI ml. g l. sin. l 0,,2,...,n A. ( t vs. ), and ( vs. ). d dt. sin l
Advanced Control Theory Homework #3 Student No.597 Name : Jinseong Kim. The simple pendulum dynamics is given by g sin L A. Derive above dynamic equation from the free body diagram. B. Find the equilibrium
More information4 Derivations of the Discrete-Time Kalman Filter
Technion Israel Institute of Technology, Department of Electrical Engineering Estimation and Identification in Dynamical Systems (048825) Lecture Notes, Fall 2009, Prof N Shimkin 4 Derivations of the Discrete-Time
More informationLINEARIZED ORBIT COVARIANCE GENERATION AND PROPAGATION ANALYSIS VIA SIMPLE MONTE CARLO SIMULATIONS
LINEARIZED ORBIT COVARIANCE GENERATION AND PROPAGATION ANALYSIS VIA SIMPLE MONTE CARLO SIMULATIONS Chris Sabol *, Thomas Sukut, Keric Hill, Kyle T. Alfriend, Brendan Wright **, You Li **, and Paul Schumacher
More informationUnderstanding the Differences between LS Algorithms and Sequential Filters
Understanding the Differences between LS Algorithms and Sequential Filters In order to perform meaningful comparisons between outputs from a least squares (LS) orbit determination algorithm and orbit determination
More informationCOVARIANCE DETERMINATION, PROPAGATION AND INTERPOLATION TECHNIQUES FOR SPACE SURVEILLANCE. European Space Surveillance Conference 7-9 June 2011
COVARIANCE DETERMINATION, PROPAGATION AND INTERPOLATION TECHNIQUES FOR SPACE SURVEILLANCE European Space Surveillance Conference 7-9 June 2011 Pablo García (1), Diego Escobar (1), Alberto Águeda (1), Francisco
More informationA Fast Method for Embattling Optimization of Ground-Based Radar Surveillance Network
A Fast Method for Embattling Optimization of Ground-Based Radar Surveillance Network JIANG Hai University of Chinese Academy of Sciences National Astronomical Observatories, Chinese Academy of Sciences
More informationMODELLING OF PERTURBATIONS FOR PRECISE ORBIT DETERMINATION
MODELLING OF PERTURBATIONS FOR PRECISE ORBIT DETERMINATION 1 SHEN YU JUN, 2 TAN YAN QUAN, 3 TAN GUOXIAN 1,2,3 Raffles Science Institute, Raffles Institution, 1 Raffles Institution Lane, Singapore E-mail:
More informationA Stochastic Collocation based. for Data Assimilation
A Stochastic Collocation based Kalman Filter (SCKF) for Data Assimilation Lingzao Zeng and Dongxiao Zhang University of Southern California August 11, 2009 Los Angeles Outline Introduction SCKF Algorithm
More informationNONLINEAR BAYESIAN FILTERING FOR STATE AND PARAMETER ESTIMATION
NONLINEAR BAYESIAN FILTERING FOR STATE AND PARAMETER ESTIMATION Kyle T. Alfriend and Deok-Jin Lee Texas A&M University, College Station, TX, 77843-3141 This paper provides efficient filtering algorithms
More informationMEETING ORBIT DETERMINATION REQUIREMENTS FOR A SMALL SATELLITE MISSION
MEETING ORBIT DETERMINATION REQUIREMENTS FOR A SMALL SATELLITE MISSION Adonis Pimienta-Peñalver, Richard Linares, and John L. Crassidis University at Buffalo, State University of New York, Amherst, NY,
More informationRelevant Applications of Differential Algebra in Astrodynamics
AstroNet-II International Final Conference Tossa de Mar 15-19 June 2015 Relevant Applications of Differential Algebra in Astrodynamics Pierluigi Di Lizia in collaboration with: Roberto Armellin, Alessandro
More informationA Spectral Approach to Linear Bayesian Updating
A Spectral Approach to Linear Bayesian Updating Oliver Pajonk 1,2, Bojana V. Rosic 1, Alexander Litvinenko 1, and Hermann G. Matthies 1 1 Institute of Scientific Computing, TU Braunschweig, Germany 2 SPT
More informationCOMPARISON OF ANGLES ONLY INITIAL ORBIT DETERMINATION ALGORITHMS FOR SPACE DEBRIS CATALOGUING
COMPARISON OF ANGLES ONLY INITIAL ORBIT DETERMINATION ALGORITHMS FOR SPACE DEBRIS CATALOGUING Fran Martinez Fadrique, Alberto Águeda Maté, Joan Jorquera Grau, Jaime Fernández Sánchez, Laura Aivar García
More informationErgodicity in data assimilation methods
Ergodicity in data assimilation methods David Kelly Andy Majda Xin Tong Courant Institute New York University New York NY www.dtbkelly.com April 15, 2016 ETH Zurich David Kelly (CIMS) Data assimilation
More informationOrbital Error Analysis for Surveillance of Space
Orbital Error Analysis for Surveillance of Space N M Harwood Air and Weapons Systems Defence Science and Technology Laboratory, UK nmharwood@dstl.gov.uk M Rutten Intelligence, Surveillance and Reconnaissance
More informationAnalytical Non-Linear Uncertainty Propagation: Theory And Implementation
Analytical Non-Linear Uncertainty Propagation: Theory And Implementation Kohei Fujimoto 1), Daniel J. Scheeres 2), and K. Terry Alfriend 1) May 13, 2013 1) Texas A&M University 2) The University of Colorado
More informationProbabilistic Fundamentals in Robotics. DAUIN Politecnico di Torino July 2010
Probabilistic Fundamentals in Robotics Gaussian Filters Basilio Bona DAUIN Politecnico di Torino July 2010 Course Outline Basic mathematical framework Probabilistic models of mobile robots Mobile robot
More informationV Control Distance with Dynamics Parameter Uncertainty
V Control Distance with Dynamics Parameter Uncertainty Marcus J. Holzinger Increasing quantities of active spacecraft and debris in the space environment pose substantial risk to the continued safety of
More informationORBIT DETERMINATION AND DATA FUSION IN GEO
ORBIT DETERMINATION AND DATA FUSION IN GEO Joshua T. Horwood, Aubrey B. Poore Numerica Corporation, 4850 Hahns Pea Drive, Suite 200, Loveland CO, 80538 Kyle T. Alfriend Department of Aerospace Engineering,
More information8 th US/Russian Space Surveillance Workshop
8 th US/Russian Space Surveillance Workshop Wailea Marriott Resort Wailea, Maui, HI 18-23 April 2010 P. Kenneth Seidelmann General Chair Kyle T. Alfriend US Technical Chair Stanislav Veniaminov Russian
More informationOrbit Determination Using Satellite-to-Satellite Tracking Data
Chin. J. Astron. Astrophys. Vol. 1, No. 3, (2001 281 286 ( http: /www.chjaa.org or http: /chjaa.bao.ac.cn Chinese Journal of Astronomy and Astrophysics Orbit Determination Using Satellite-to-Satellite
More informationData assimilation in high dimensions
Data assimilation in high dimensions David Kelly Courant Institute New York University New York NY www.dtbkelly.com February 12, 2015 Graduate seminar, CIMS David Kelly (CIMS) Data assimilation February
More informationECE276A: Sensing & Estimation in Robotics Lecture 10: Gaussian Mixture and Particle Filtering
ECE276A: Sensing & Estimation in Robotics Lecture 10: Gaussian Mixture and Particle Filtering Lecturer: Nikolay Atanasov: natanasov@ucsd.edu Teaching Assistants: Siwei Guo: s9guo@eng.ucsd.edu Anwesan Pal:
More informationPartial J 2 -Invariance for Spacecraft Formations
Partial J 2 -Invariance for Spacecraft Formations Louis Breger and Jonathan P. How MIT Department of Aeronautics and Astronautics Kyle T. Alfriend Texas A&M University Department of Aerospace Engineering
More informationCENTROID DYNAMICS FOR GROUP OBJECT TRACKING
(Preprint) AAS 15-365 CENTROID DYNAMICS FOR GROUP OBJECT TRACKING Christopher Binz and Liam Healy Immediately following a scenario like a CubeSat deployment, the presence of multiple objects in close proximity
More informationModeling, Dynamics and Control of Spacecraft Relative Motion in a Perturbed Keplerian Orbit
Paper Int l J. of Aeronautical & Space Sci. 16(1), 77 88 (2015) http://dx.doi.org/10.5139/ijass.2015.16.1.77 Modeling, Dynamics and Control of Spacecraft Relative Motion in a Perturbed Keplerian Orbit
More informationCONDENSATION Conditional Density Propagation for Visual Tracking
CONDENSATION Conditional Density Propagation for Visual Tracking Michael Isard and Andrew Blake Presented by Neil Alldrin Department of Computer Science & Engineering University of California, San Diego
More informationTTK4190 Guidance and Control Exam Suggested Solution Spring 2011
TTK4190 Guidance and Control Exam Suggested Solution Spring 011 Problem 1 A) The weight and buoyancy of the vehicle can be found as follows: W = mg = 15 9.81 = 16.3 N (1) B = 106 4 ( ) 0.6 3 3 π 9.81 =
More information2D Image Processing. Bayes filter implementation: Kalman filter
2D Image Processing Bayes filter implementation: Kalman filter Prof. Didier Stricker Kaiserlautern University http://ags.cs.uni-kl.de/ DFKI Deutsches Forschungszentrum für Künstliche Intelligenz http://av.dfki.de
More informationNonlinear Filtering. With Polynomial Chaos. Raktim Bhattacharya. Aerospace Engineering, Texas A&M University uq.tamu.edu
Nonlinear Filtering With Polynomial Chaos Raktim Bhattacharya Aerospace Engineering, Texas A&M University uq.tamu.edu Nonlinear Filtering with PC Problem Setup. Dynamics: ẋ = f(x, ) Sensor Model: ỹ = h(x)
More informationStochastic Models, Estimation and Control Peter S. Maybeck Volumes 1, 2 & 3 Tables of Contents
Navtech Part #s Volume 1 #1277 Volume 2 #1278 Volume 3 #1279 3 Volume Set #1280 Stochastic Models, Estimation and Control Peter S. Maybeck Volumes 1, 2 & 3 Tables of Contents Volume 1 Preface Contents
More informationROBOTICS 01PEEQW. Basilio Bona DAUIN Politecnico di Torino
ROBOTICS 01PEEQW Basilio Bona DAUIN Politecnico di Torino Probabilistic Fundamentals in Robotics Gaussian Filters Course Outline Basic mathematical framework Probabilistic models of mobile robots Mobile
More informationECONOMETRIC METHODS II: TIME SERIES LECTURE NOTES ON THE KALMAN FILTER. The Kalman Filter. We will be concerned with state space systems of the form
ECONOMETRIC METHODS II: TIME SERIES LECTURE NOTES ON THE KALMAN FILTER KRISTOFFER P. NIMARK The Kalman Filter We will be concerned with state space systems of the form X t = A t X t 1 + C t u t 0.1 Z t
More informationRESEARCH OBJECTIVES AND SUMMARY OF RESEARCH
XVIII. DETECTION AND ESTIMATION THEORY Academic Research Staff Prof. H. L. Van Trees Prof. A. B. Baggeroer Graduate Students J. P. Albuquerque J. M. F. Moura L. C. Pusey L. S. Metzger S. Orphanoudakis
More information* is a row matrix * An mxn matrix is a square matrix if and only if m=n * A= is a diagonal matrix if = 0 i
CET MATRICES *A matrix is an order rectangular array of numbers * A matrix having m rows and n columns is called mxn matrix of order * is a column matrix * is a row matrix * An mxn matrix is a square matrix
More informationCOMBINED ADAPTIVE CONTROLLER FOR UAV GUIDANCE
COMBINED ADAPTIVE CONTROLLER FOR UAV GUIDANCE B.R. Andrievsky, A.L. Fradkov Institute for Problems of Mechanical Engineering of Russian Academy of Sciences 61, Bolshoy av., V.O., 199178 Saint Petersburg,
More informationPROBABILITY OF COLLISION WITH SPECIAL PERTURBATIONS DYNAMICS USING THE MONTE CARLO METHOD
AAS 11-435 PROBABILITY OF COLLISION WITH SPECIAL PERTURBATIONS DYNAMICS USING THE MONTE CARLO METHOD Chris Sabol, Christopher Binz and Alan Segerman, Kevin Roe, and Paul W. Schumacher, Jr. INTRODUCTION
More informationA Spherical Coordinate Parametrization for an In-Orbit Bearings-Only Navigation Filter
Proceedings of the EuroGNC 13, nd CEAS Specialist Conference on Guidance, Navigation & Control, Delft University of Technology, Delft, The Netherlands, April 1-1, 13 WeCT1.1 A Spherical Coordinate Parametrization
More informationMinimizing Bearing Bias in Tracking By De-coupled Rotation and Translation Estimates
Minimizing Bearing Bias in Tracking By De-coupled Rotation and Translation Estimates Raman Arora Department of Electrical Engineering University of Washington Seattle, WA 9895 Email: rmnarora@u.washington.edu
More informationSome Results Concerning Uniqueness of Triangle Sequences
Some Results Concerning Uniqueness of Triangle Sequences T. Cheslack-Postava A. Diesl M. Lepinski A. Schuyler August 12 1999 Abstract In this paper we will begin by reviewing the triangle iteration. We
More informationNonlinear State Estimation! Extended Kalman Filters!
Nonlinear State Estimation! Extended Kalman Filters! Robert Stengel! Optimal Control and Estimation, MAE 546! Princeton University, 2017!! Deformation of the probability distribution!! Neighboring-optimal
More informationJimmy Dale Hicks Jr. Auburn, Alabama December 8, Keywords: orbit debris, state estimation, kalman filtering
Performance Comparison of an Extended Kalman Filter and an Iterated Extended Kalman Filter for Orbit Determination of Space Debris with Poor Apriori Information and Intermittent Observations by Jimmy Dale
More informationIterative Learning Control Analysis and Design I
Iterative Learning Control Analysis and Design I Electronics and Computer Science University of Southampton Southampton, SO17 1BJ, UK etar@ecs.soton.ac.uk http://www.ecs.soton.ac.uk/ Contents Basics Representations
More informationThe Use of Short-Arc Angle and Angle Rate Data for Deep-Space Initial Orbit Determination and Track Association
The Use of Short-Arc Angle and Angle Rate Data for Deep-Space Initial Orbit Determination and Track Association Dr. Moriba Jah (AFRL) Mr. Kyle DeMars (UT-Austin) Dr. Paul Schumacher Jr. (AFRL) Background/Motivation
More informationLecture 2c: Satellite Orbits
Lecture 2c: Satellite Orbits Outline 1. Newton s Laws of Mo3on 2. Newton s Law of Universal Gravita3on 3. Kepler s Laws 4. Pu>ng Newton and Kepler s Laws together and applying them to the Earth-satellite
More informationRobotics 2 Target Tracking. Kai Arras, Cyrill Stachniss, Maren Bennewitz, Wolfram Burgard
Robotics 2 Target Tracking Kai Arras, Cyrill Stachniss, Maren Bennewitz, Wolfram Burgard Slides by Kai Arras, Gian Diego Tipaldi, v.1.1, Jan 2012 Chapter Contents Target Tracking Overview Applications
More informationStochastic State Estimation
Chapter 7 Stochastic State Estimation Perhaps the most important part of studying a problem in robotics or vision, as well as in most other sciences, is to determine a good model for the phenomena and
More informationTarget tracking and classification for missile using interacting multiple model (IMM)
Target tracking and classification for missile using interacting multiple model (IMM Kyungwoo Yoo and Joohwan Chun KAIST School of Electrical Engineering Yuseong-gu, Daejeon, Republic of Korea Email: babooovv@kaist.ac.kr
More information1 Kalman Filter Introduction
1 Kalman Filter Introduction You should first read Chapter 1 of Stochastic models, estimation, and control: Volume 1 by Peter S. Maybec (available here). 1.1 Explanation of Equations (1-3) and (1-4) Equation
More information5 Kalman filters. 5.1 Scalar Kalman filter. Unit delay Signal model. System model
5 Kalman filters 5.1 Scalar Kalman filter 5.1.1 Signal model System model {Y (n)} is an unobservable sequence which is described by the following state or system equation: Y (n) = h(n)y (n 1) + Z(n), n
More informationSolutions to exam : 1FA352 Quantum Mechanics 10 hp 1
Solutions to exam 6--6: FA35 Quantum Mechanics hp Problem (4 p): (a) Define the concept of unitary operator and show that the operator e ipa/ is unitary (p is the momentum operator in one dimension) (b)
More information1/30. Rigid Body Rotations. Dave Frank
. 1/3 Rigid Body Rotations Dave Frank A Point Particle and Fundamental Quantities z 2/3 m v ω r y x Angular Velocity v = dr dt = ω r Kinetic Energy K = 1 2 mv2 Momentum p = mv Rigid Bodies We treat a rigid
More informationMathematical Theory of Control Systems Design
Mathematical Theory of Control Systems Design by V. N. Afarias'ev, V. B. Kolmanovskii and V. R. Nosov Moscow University of Electronics and Mathematics, Moscow, Russia KLUWER ACADEMIC PUBLISHERS DORDRECHT
More informationA new unscented Kalman filter with higher order moment-matching
A new unscented Kalman filter with higher order moment-matching KSENIA PONOMAREVA, PARESH DATE AND ZIDONG WANG Department of Mathematical Sciences, Brunel University, Uxbridge, UB8 3PH, UK. Abstract This
More informationDesign of FIR Smoother Using Covariance Information for Estimating Signal at Start Time in Linear Continuous Systems
Systems Science and Applied Mathematics Vol. 1 No. 3 2016 pp. 29-37 http://www.aiscience.org/journal/ssam Design of FIR Smoother Using Covariance Information for Estimating Signal at Start Time in Linear
More informationFundamentals of Astrodynamics and Applications
Fundamentals of Astrodynamics and Applications Third Edition David A. Vallado with technical contributions by Wayne D. McClain Space Technology Library Published Jointly by Microcosm Press Hawthorne, CA
More informationAttitude 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 informationAnalytical Non-Linear Uncertainty Propagation: Theory And Implementation
Analytical Non-Linear Uncertainty Propagation: Theory And Implementation Kohei Fujimoto 1 Texas A&M University, College Station, TX, 77843 Daniel J. Scheeres 2 The University of Colorado-Boulder, Boulder,
More informationA Statistical Study of PITCHf/x Pitched Baseball Trajectories
A Statistical Study of PITCHf/x Pitched Baseball Trajectories Alan M. Nathan Department of Physics, University of Illinois, Urbana, IL 6181 February 8, 28 A statistical study of PITCHf/x trajectories is
More informationON SELECTING THE CORRECT ROOT OF ANGLES-ONLY INITIAL ORBIT DETERMINATION EQUATIONS OF LAGRANGE, LAPLACE, AND GAUSS
AAS 16-344 ON SELECTING THE CORRECT ROOT OF ANGLES-ONLY INITIAL ORBIT DETERMINATION EQUATIONS OF LAGRANGE, LAPLACE, AND GAUSS Bong Wie and Jaemyung Ahn INTRODUCTION This paper is concerned with a classical
More informationState estimation and the Kalman filter
State estimation and the Kalman filter PhD, David Di Ruscio Telemark university college Department of Technology Systems and Control Engineering N-3914 Porsgrunn, Norway Fax: +47 35 57 52 50 Tel: +47 35
More informationSERVO TRACKING OF TARGETS AT SEA
Int. J. Appl. Math. Comput. Sci., 2006, Vol. 16, No. 2, 197 207 SERVO TRACKING OF TARGETS AT SEA ANIS AHMED,,MIECZYSŁAW A. BRDYŚ Department of Electronic, Electrical and Computer Engineering School of
More information2D Image Processing. Bayes filter implementation: Kalman filter
2D Image Processing Bayes filter implementation: Kalman filter Prof. Didier Stricker Dr. Gabriele Bleser Kaiserlautern University http://ags.cs.uni-kl.de/ DFKI Deutsches Forschungszentrum für Künstliche
More informationDIRECTIVITY OPTIMIZATION IN PLANAR SUB-ARRAYED MONOPULSE ANTENNA
Progress In Electromagnetics Research Letters, Vol. 4, 1 7, 2008 DIRECTIVITY OPTIMIZATION IN PLANAR SUB-ARRAYED MONOPULSE ANTENNA P. Rocca, L. Manica, and A. Massa ELEDIA Department of Information Engineering
More informationon space debris objects obtained by the
KIAM space debris data center for processing and analysis of information on space debris objects obtained by the ISON network Vladimir Agapov, Igor Molotov Keldysh Institute of Applied Mathematics RAS
More informationAIR FORCE INSTITUTE OF TECHNOLOGY
Space Based Satellite Tracking and Characterization Utilizing Non-Imaging Passive Sensors THESIS Bradley R. Townsend, Captain, USA AFIT/GA/ENY/08-M06 DEPARTMENT OF THE AIR FORCE AIR UNIVERSITY AIR FORCE
More informationAdaptive State Estimation Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 2018
Adaptive State Estimation Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 218! Nonlinearity of adaptation! Parameter-adaptive filtering! Test for whiteness of the residual!
More informationTactical Ballistic Missile Tracking using the Interacting Multiple Model Algorithm
Tactical Ballistic Missile Tracking using the Interacting Multiple Model Algorithm Robert L Cooperman Raytheon Co C 3 S Division St Petersburg, FL Robert_L_Cooperman@raytheoncom Abstract The problem of
More informationACCURACY ASSESSMENT OF GEOSTATIONARY-EARTH-ORBIT WITH SIMPLIFIED PERTURBATIONS MODELS
ARTIFICIAL SATELLITES, Vol. 51, No. 2 2016 DOI: 10.1515/arsa-2016-0005 ACCURACY ASSESSMENT OF GEOSTATIONARY-EARTH-ORBIT WITH SIMPLIFIED PERTURBATIONS MODELS Lihua Ma, Xiaojun Xu, Feng Pang National Astronomical
More informationRobot Localization and Kalman Filters
Robot Localization and Kalman Filters Rudy Negenborn rudy@negenborn.net August 26, 2003 Outline Robot Localization Probabilistic Localization Kalman Filters Kalman Localization Kalman Localization with
More informationGNSS: Global Navigation Satellite Systems
GNSS: Global Navigation Satellite Systems Global: today the American GPS (Global Positioning Service), http://gps.losangeles.af.mil/index.html the Russian GLONASS, http://www.glonass-center.ru/frame_e.html
More informationECE531 Lecture 11: Dynamic Parameter Estimation: Kalman-Bucy Filter
ECE531 Lecture 11: Dynamic Parameter Estimation: Kalman-Bucy Filter D. Richard Brown III Worcester Polytechnic Institute 09-Apr-2009 Worcester Polytechnic Institute D. Richard Brown III 09-Apr-2009 1 /
More informationSequential Monte Carlo methods for filtering of unobservable components of multidimensional diffusion Markov processes
Sequential Monte Carlo methods for filtering of unobservable components of multidimensional diffusion Markov processes Ellida M. Khazen * 13395 Coppermine Rd. Apartment 410 Herndon VA 20171 USA Abstract
More informationNonlinear State Estimation! Particle, Sigma-Points Filters!
Nonlinear State Estimation! Particle, Sigma-Points Filters! Robert Stengel! Optimal Control and Estimation, MAE 546! Princeton University, 2017!! Particle filter!! Sigma-Points Unscented Kalman ) filter!!
More informationON ENTRY-WISE ORGANIZED FILTERING. E.A. Suzdaleva
ON ENTRY-WISE ORGANIZED FILTERING E.A. Suzdaleva Department of Adaptive Systems Institute of Information Theory and Automation of the Academy of Sciences of the Czech Republic Pod vodárenskou věží 4, 18208
More informationOrbit Determination of Satellite Formations. Terry Alfriend 9 th US Russian Space Surveillance Workshop
Orbit Determination of Satellite Formations Terry Alfriend 9 th US Russian Space Surveillance Workshop Outline What is a Satellite Formation Types of Formations Proposed Approach for Orbit Determination
More informationStochastic Processes. M. Sami Fadali Professor of Electrical Engineering University of Nevada, Reno
Stochastic Processes M. Sami Fadali Professor of Electrical Engineering University of Nevada, Reno 1 Outline Stochastic (random) processes. Autocorrelation. Crosscorrelation. Spectral density function.
More informationRobotics 2 Target Tracking. Giorgio Grisetti, Cyrill Stachniss, Kai Arras, Wolfram Burgard
Robotics 2 Target Tracking Giorgio Grisetti, Cyrill Stachniss, Kai Arras, Wolfram Burgard Linear Dynamical System (LDS) Stochastic process governed by is the state vector is the input vector is the process
More informationInfrared 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 informationPeriodic Dynamic Traffic Assignment Problem
Periodic Dynamic Traffic Assignment Problem Artyom Nahapetyan Center for Applied Optimization, ISE Department University of Florida Joint work with: Siriphong Lawphongpanich and Donald W. Hearn Periodic
More informationResults and Analyses of Debris Tracking from Mt Stromlo
Results and Analyses of Debris Tracking from Mt Stromlo Jizhang Sang 1, Ian Ritchie, Matt Pearson, and Craig Smith EOS Space Systems Pty Ltd, Australia Abstract In the last 2 years, EOS Space Systems has
More informationOutline. 1 Introduction. 2 Problem and solution. 3 Bayesian tracking model of group debris. 4 Simulation results. 5 Conclusions
Outline 1 Introduction 2 Problem and solution 3 Bayesian tracking model of group debris 4 Simulation results 5 Conclusions Problem The limited capability of the radar can t always satisfy the detection
More information