Tactical Ballistic Missile Tracking using the Interacting Multiple Model Algorithm
|
|
- Melvyn Perry
- 6 years ago
- Views:
Transcription
1 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 tracking a tactical ballistic missile is complicated by the varying target dynamics in the boost, eo-atmospheric and endo-atmospheric phases of flight A single Kalman filter tuned for constant velocity or constant acceleration is not based upon the correct underlying physical model of these dynamics By including models with all of the correct dynamics, an Interacting Multiple Model (IMM) simultaneously weights all of them and adapts to the one most closely matching the data based upon measurement residuals A multiple-sensor application of this algorithm requires either a single IMM driven by measurements from all sensors (measurement fusion) or an IMM for each sensor driven by its own measurements, followed by fusion across sensors (track fusion) This paper develops a tactical ballistic missile tracker within an IMM framework and gives an eample of the measurement fusion approach for a simulated trajectory and simulated sensor geometry Keywords: Tracking, filtering, Interacting Multiple Model, Tactical Ballistic Missile, Kalman filter 1 Introduction The Interacting Multiple Model (IMM) approach to target tracking has been in use for over a decade, mainly in the area of air defense, in which the goal is to reduce the lags that develop while tracking highly maneuvering manned aircraft These lags develop when the underlying motion model for the target is constant velocity (CV) and the motion deviates substantially from this model, as in a maneuver The simplest IMM for this application is a bank of tracking filters (usually Kalman or Etended Kalman filters (EKFs)), in which each model is tuned to a different acceleration, by means of the Kalman filter process noise A large value of process noise is used for a large acceleration and a small value for a small acceleration The track outputs of the multiple models included in the IMM are combined linearly, with weights that depend upon the likelihood that a measurement fits the assumption of each of the models The number of models to use in the IMM is largely a matter of eperiment, but most implementations use two, or at most three For the tactical ballistic missile (TBM) application in this paper, there are three models, corresponding to the three regions of a TBM trajectory: boost, eo-atmospheric (ballistic) and endo-atmospheric (re-entry) The boost model is a 9-state EKF, in Cartesian coordinates centered on the sensor declared to be local, for the purpose of composite tracking The state elements are position, velocity and acceleration The ballistic model is a 6-state EKF with gravity and Coriolis terms The state elements are position and velocity State propagation, however, includes gravity and Coriolis forces, even though the state does not contain acceleration The re-entry model is a 7-state EKF, identical to the ballistic state but has a 7 th element, the (inverse) ballistic coefficient In this paper, the issue of contact-to-track association is not addressed That is, all contacts from all sensors are assumed to be correctly associated to a track Also, the practical issue of track initiation is not addressed It is assumed that the first two contacts arriving in chronological order (from either one or two different sensors) initiate the track and all further contacts are used for track update Simulated truth trajectories and radar contact data were generated by numerically solving the nd order differential equations of motion for a unitary missile with the Runge-Kutta method and adding simulated measurement noise The code provided in [1] was used to compute the truth trajectories A multiple boost missile was simulated by allowing multiple burn periods, with given weight-thrust characteristics and eperiencing atmospheric drag The missile simulated for this paper was not intended to be representative of any real TBM, but serves for illustrative purposes [1] The single IMM in this study fused contacts from multiple sensors by injecting a chronological stream of associated measurement reports (AMRs); ie, measurements which have been assumed to be associated to the same track The study was limited to two sensors, one defined as local and one remote The local sensor is the sensor in which the track state is defined Prediction and filtering in each IMM EKF involved both the local and 84
2 remote measurement matri H, defined in the EKF formulation [] followed in this study specifying the acceleration terms which affect the transition step of the EKF IMM Formulation Figure 1 shows the architecture of the IMM algorithm used in this study At each time step k, a linear combination of the previous outputs (states and covariances) is input into each model Also, the current measurement is input into each model and residuals are computed, along with corresponding likelihood functions Normalized likelihood functions are used as weights in a linear combination of current model outputs to form the desired blended track state and covariance outputs, while each model outputs are stored for the net iteration In addition, model probabilities (µ) are updated for the net iteration µ BOOST (k-1), µ EXO (k-1), µ ENDO (k-1) X BOOST k-1 k-1 XEXO k-1 k-1 XENDO k-1 k-1 3 TBM Equations The TBM in this study is assumed to be a point target with three degrees of freedom, its position in space This is an obvious simplification of the real TBM problem, ignoring the additional three rotational degrees of freedom The acceleration in an Earth-Centered Earth-Fied (ECEF) is given by * 1 a = ρv vc ˆ D A g mg R * r rˆ * * * * * ( ω v) ω ( ω r) (1) M k Interaction (Miing) X BOOST k-1 k-1 X EXO k-1 k-1 X ENDO k-1 k-1 Boost Model Λ BOOST (k) Eo Model Endo Model Λ EXO (k) Λ ENDO (k) Model Probability Update State Estimate Combination where the first term is atmospheric drag, acting opposite in direction to velocity, the second term is gravity, the third term is Coriolis acceleration and the last term is centrifugal acceleration The missile mass is m, the gravitational acceleration is g (98 m/s ), the velocity is v, shown with the unit vector in its direction, the missile effective surface area is A and the (zero-lift) drag is C D It is further assumed that there is no acceleration due to lift The Earth s equatorial radius is R ( m) [4], the missile position is r and the Earth s rotational velocity is [ -5 rad/s) [4] Although Equation (1) is written for simplicity for a spherical Earth, the tracker developed in this study incorporated the WGS-84 Earth Model [4] The tracker was developed in the local East-North-Up (ENU) coordinate system In this coordinate system, the acceleration components (for a spherical Earth) are given by X BOOST k k X EXO k k X ENDO k k Figure 1 TBM IMM Architecture X k k µ BOOST (k), µ EXO (k), µ ENDO (k) In this figure, the state vectors (and corresponding covariance matri, not shown for simplicity) for all models are combined linearly with the model probabilities as coefficients The mathematical details of the IMM and the EKF can be found in many references [,3] and are not repeated here The mathematics related to the TBM problem involves ρg gr = R * β 3 n z u y * * R + R + ( ω γ ω γ ) u n 3 ( ω γ ω γ ) 3 1 ( ω ω ) ρg gr y y = R * y β 3 u e z * * R + R + e u ( ω ω ) () (3) 85
3 ρg z = R * z β gr ( R + z) * * 3 R + R + ( ω ω ) ( ω γ ω γ ) n 1 e y e where the ENU components of ω at geodetic latitude µ are given by ωe ωn = = ω cos µ n (4), (5) F () = f = y y (1) ωu = ω sin µ the components of vector γ are given by γ 1 = ωn γ = ωu ( R + z) ωu y ω ( R + z) e γ 3 = ωe y ωn the ballistic coefficient is given by, (6) β = mg (7) C D A and R * is the vector to the missile from the origin of the local ENU coordinate system The endo-atmospheric track state includes a drag term α which is defined as 1/β These acceleration components are used in the prediction step of the EKF for the 6-state eo-atmospheric and 7-state endo-atmospheric models The predicted state is given by 1 ˆ = ~ + f () ~ t + F()()( ~ f ~ t) (8) where f () ~ is the time derivative of the state vector, given by f () [ y y z z α ] T = (9) for the endo-atmospheric model and the same vector for the eo-atmospheric model without the drag term The function F () is a Jacobian matri of partial derivatives given by (for the endo-atmospheric model) (The α terms are dropped for the eo-atmospheric model) This matri is a very complicated combination of position (target and sensor), velocity and acceleration terms, which will not be presented here for the sake of brevity As mentioned above, a WGS-84 Earth model is included in these terms, allowing for a non-spherical Earth in the tracking algorithm [5] 4 Results This study developed a multi-sensor TBM tracker within an IMM framework and applied it to a simulated trajectory The results of the study measured the ability of a multiple-sensor TBM track to provide better track accuracy than a single-sensor TBM track for a hypothetical missile, sensors and sensor geometry For this purpose, one sensor was chosen to be close to the launch point of the TBM and one sensor was chosen to be close to the impact point of the TBM The sensors were chosen to have long enough detection range to follow the missile along all phases of its trajectory In reality, both sensors would not necessarily track the entire trajectory and further benefit would be derived from the viewpoint of track continuity, with one sensor filling in when the other sensor drops track However, this study only considered the track accuracy benefit from overlapping coverage Future work will evaluate the benefits derived from multiple geographically dispersed sensors Figure shows the modeled TBM trajectory in ENU coordinates with origin at Sensor 1 86
4 This was done in order to compare results with the single-sensor case The first comparison provided the same total amount of data, while the second comparison provided twice the amount of data, because each sensor in the network was assumed to have the same capability as if it acted alone Both sensors were assumed perfect, ecept for measurement noise, with no false alarms or missed detections and no bias (either sensor or gridlock) The 1 σ measurement noise was assumed to be 1 m in range and 5 mrad in bearing and elevation It is also noted that the sensors are assumed to measure position only, but not velocity Figure Simulated TBM Trajectory The missile was launched at t= from a range of approimately 4 km from Sensor 1, placed near the impact point Sensor was located near the launch point Figure 3 shows the speed and velocity components of the missile in Sensor 1 coordinates The speed and z- component clearly show the two boost stages, one lasting from t= to t=3 seconds and the net from t=3 to t=138 seconds This is followed by the ballistic motion with z- component decreasing to zero at apogee and then by negative z-component as the missile falls to Earth Atmospheric drag acts to slow it down as it continues, so that all velocity components approach zero Figure 4 shows the corresponding acceleration and components The two boost stages appear as two impulses, followed by constant ballistic acceleration of approimately 1g in the negative z-directon (It is actually slightly less than 1g because the acceleration due to gravity decreases with altitude) The decceleration due to drag appears as a positive acceleration because the negative speed is getting less negative, so that the speed profile (Figure 3) has a positive slope, equal to the acceleration AMRs were processed by the multiple EKFs in the TBM IMM filter in chronological order received The approach taken was to quantify the tracking performance using Sensor 1 alone and to quantify the benefit of an additional sensor Three cases were considered : Case 1 had Sensor 1 alone, with an update interval of 1 second Case had both sensors, but each had an update interval of seconds, staggered by 1 second, giving a network update interval of 1 second ; ie, the same total amount of data as for Case 1, but from geometrically diverse sources Case 3 had both sensors, but each was allowed to have an update interval of 1 second (staggered by 5 second), as in Case 1 The TBM IMM algorithm was applied to this set of data For brevity, the plots below compare Case 1 and Case, while the error plots compare all three cases Figures 5 and 6 show the track and truth along the trajectory for Cases 1 and, respectively Ecept for initial fluctuations in the boost phase of Case 1, the track appears reasonably close to the truth Figure 6 shows that measurement fusion improves the track, especially in boost phase Figure 3 TBM Velocity Figure 4 TBM Acceleration 87
5 Figure 5 TBM Truth and Single Sensor IMM Track (Case 1) Figure 7 IMM Model Probabilities transitions occur at the proper times and are very sharp The only noisy region in these curves is in the endo-atmospheric motion, where the algorithm is slightly confused in selecting between the eo- and endo- models, which differ only in the presence of a small non-zero element in the 7th component of the track state in the endo-atmospheric model Figures 8 and 9 show the IMM-derived velocity components and speed corresponding to the boost stage truth data in Figure 3, for Case 1 and, respectively The z- component and the speed are positive, while the - and y- components are negative Figure 6 TBM Truth and Two-Sensor IMM Track (Case ) The eo-atmospheric two-sensor track is smoother than the track with Sensor 1 updates alone, but appears not to quite reach apogee The endo-atmospheric two-sensor track appears approimately the same as for Sensor 1 alone Figure 7 shows the IMM model probabilities for the three models vs time At the bottom is shown the target acceleration (magnitude curve from Figure 4) for comparison As seen from this acceleration plot, the model Figure 8 IMM and True Boost Stage Velocity (Case 1) 88
6 Acceleration estimates from the IMM are shown in Figures 1 and 13, only for the boost stage The boost stage filter is the only filter with acceleration included in the state vector The smooth velocity behavior shown in Figures 9 through 11 for the other stages will suffice to demonstrate good performance for those stages, rather than to compute acceleration with a noisy divided difference approimation Figure 9 IMM and True Boost Stage Velocity (Case ) A comparison shows that the additional sensor smoothes the fluctuations The corresponding curves for the eo- and endo- atmospheric regions differ less from each other on this scale, compared to the boost region and so only the curves for Case are presented in Figures 1 and 11 Figure 1 IMM and True Boost Acceleration (Case 1) Figure 1 IMM and True Ballistic Velocity (Case ) Figure 13 IMM and True Boost Acceleration (Case ) An improvement in Sensor 1 s track in the boost phase due to the contribution of Sensor is seen by comparing these figures The peak acceleration in the first boost stage is detected and estimated with only a few seconds delay in Figure 13 compared to a delay of approimately 15 seconds in Figure 1 In addition, the acceleration estimate itself is closer to the truth for both boost stages Also, the fluctuations in Figure 1 are reduced in Figure 13 Figure 11 IMM and True Re-entry Velocity (Case ) Figure 14 shows the estimate of the ballistic coefficient during the endo-atmospheric phase Also plotted is the value of α corresponding to the value of β = 1/α built into the simulated truth trajectory (5 lb/ft ) It was computed 89
7 for both cases and the results were essentially the same A non-zero estimate begins to be computed as the endoatmospheric model increases in importance, as seen in Figure 7 Figure 14 Ballistic Coefficient Estimation To complete the study, attention is turned to Figures These figures show position and speed errors (root-sumsquared error between truth and track) over the entire trajectory for the three cases, as well as corresponding estimation error (square root of filtered covariance trace) Statistics were computed over time, for each of the three regions of target motion The covariance error is very Figure 16 Position and Speed Errors (Case ) Figure 17 Position and Speed Errors (Case 3) smooth and its close correspondence to the track error shows the near consistency of the IMM Figure 15 Position and Speed Errors (Case 1) Case 1 shows the performance of Sensor 1 alone, with a 1-second update interval Sensor 1 was near the impact point, far from the launch point The large error in its track during the boost phase resulted from the large magnification 83
8 of the angular error over the long range between Sensor 1 and the launch point This large relative boost-phase error also appeared in velocity Cases and 3 demonstrate the improvement in tracking error from measurement fusion in this eample The improvement was demonstrated in position and speed errors, as well A comparison of the composite -second update track (Case ) with the single-sensor track (Case 1) showed that boost and eo-atmospheric position errors decreased for the composite compared to Sensor 1 alone, with a slight increase in the endo-atmospheric region The 1-second composite track (Case 3) showed improvement over Case 1 for all three regions Speed errors decreased in all three regions for both Case and 3 [3] Bar Shalom, Y and KC Chang, Tracking a Maneuvering Target Using Input Estimation Versus the Interacting Multiple Model Algorithm, IEEE Transactions on Aerospace and Electronic Systems, AES-5, (March 1989), 96-3 [4] Ballistic Missile Defense Organization (BMDO), Ballistic Missile Defense (BMD) Navigation Standard, BMD-P-SD-9--A, June 3, 1993 [5] Ballistic Missile Defense Organization (BMDO) and Office of Naval Research (ONR), Functional Description of the Ballistic Missile Defense (BMD) Benchmark, Version 4, December 1,1 These results are summarized in Tables 1 and, below Table 1 Mean Position Error (m) Boost Eo Endo Case Case Case Table Mean Speed Error (m/s) Boost Eo Endo Case Case Case Conclusions A TBM tracking algorithm was developed within the framework of the IMM and applied to simulated data from multiple sensors A tracking accuracy improvement compared to a single sensor was observed when a composite track was computed for all three phases of TBM motion 6 References [1] Zarchan, P, Tactical and Strategic Missile Guidance,Vol 176, Progress in Astronautics and Aeronautics, AIAA, Washington, DC, 1997 [] Castella, FR, Multisensor, Multisite Tracking Filter, IEE Proc Radar,Sonar,Navig,Vol 141, No,April
Target 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 informationv are uncorrelated, zero-mean, white
6.0 EXENDED KALMAN FILER 6.1 Introduction One of the underlying assumptions of the Kalman filter is that it is designed to estimate the states of a linear system based on measurements that are a linear
More informationTracking and Identification of Multiple targets
Tracking and Identification of Multiple targets Samir Hachour, François Delmotte, Eric Lefèvre, David Mercier Laboratoire de Génie Informatique et d'automatique de l'artois, EA 3926 LGI2A first name.last
More informationA New Nonlinear Filtering Method for Ballistic Target Tracking
th International Conference on Information Fusion Seattle, WA, USA, July 6-9, 9 A New Nonlinear Filtering Method for Ballistic arget racing Chunling Wu Institute of Electronic & Information Engineering
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 informationNonlinear Estimation Techniques for Impact Point Prediction of Ballistic Targets
Nonlinear Estimation Techniques for Impact Point Prediction of Ballistic Targets J. Clayton Kerce a, George C. Brown a, and David F. Hardiman b a Georgia Tech Research Institute, Georgia Institute of Technology,
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 informationMULTI-MODEL FILTERING FOR ORBIT DETERMINATION DURING MANOEUVRE
MULTI-MODEL FILTERING FOR ORBIT DETERMINATION DURING MANOEUVRE Bruno CHRISTOPHE Vanessa FONDBERTASSE Office National d'etudes et de Recherches Aérospatiales (http://www.onera.fr) B.P. 72 - F-92322 CHATILLON
More informationHeterogeneous Track-to-Track Fusion
Heterogeneous Track-to-Track Fusion Ting Yuan, Yaakov Bar-Shalom and Xin Tian University of Connecticut, ECE Dept. Storrs, CT 06269 E-mail: {tiy, ybs, xin.tian}@ee.uconn.edu T. Yuan, Y. Bar-Shalom and
More informationFuzzy Inference-Based Dynamic Determination of IMM Mode Transition Probability for Multi-Radar Tracking
Fuzzy Inference-Based Dynamic Determination of IMM Mode Transition Probability for Multi-Radar Tracking Yeonju Eun, Daekeun Jeon CNS/ATM Team, CNS/ATM and Satellite Navigation Research Center Korea Aerospace
More informationSeparable warhead mathematical model of Supersonic & Hypersonic Re-entry Vehicles
16 th International Conference on AEROSPACE SCIENCES & AVIATION TECHNOLOGY, ASAT - 16 May 26-28, 2015, E-Mail: asat@mtc.edu.eg Military Technical College, Kobry Elkobbah, Cairo, Egypt Tel : +(202) 24025292
More informationA Sufficient Comparison of Trackers
A Sufficient Comparison of Trackers David Bizup University of Virginia Department of Systems and Information Engineering P.O. Box 400747 151 Engineer's Way Charlottesville, VA 22904 Donald E. Brown University
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 informationBoost phase tracking with an unscented filter
Boost phase tracking with an unscented filter James R. Van Zandt a a MITRE Corporation, MS-M210, 202 Burlington Road, Bedford MA 01730, USA ABSTRACT Boost phase missile tracking is formulated as a nonlinear
More informationRocket Science 102 : Energy Analysis, Available vs Required
Rocket Science 102 : Energy Analysis, Available vs Required ΔV Not in Taylor 1 Available Ignoring Aerodynamic Drag. The available Delta V for a Given rocket burn/propellant load is ( ) V = g I ln 1+ P
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 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 informationA MULTI-SENSOR FUSION TRACK SOLUTION TO ADDRESS THE MULTI-TARGET PROBLEM
Approved for public release; distribution is unlimited. A MULTI-SENSOR FUSION TRACK SOLUTION TO ADDRESS THE MULTI-TARGET PROBLEM By Dr. Buddy H. Jeun, Jay Jayaraman Lockheed Martin Aeronautical Systems,
More information1 The satellite altimeter measurement
1 The satellite altimeter measurement In the ideal case, a satellite altimeter measurement is equal to the instantaneous distance between the satellite s geocenter and the ocean surface. However, an altimeter
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 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 informationV Requirements for a Gun Assisted Launch to Circular Orbit
V Requirements for a Gun Assisted Launch to Circular Orbit Gerry Flanagan The Alna Space Program May 12, 2011 Introduction and Assumptions An earth-based gun can be used to send a projectile into space,
More informationTrack-to-Track Fusion Architectures A Review
Itzhack Y. Bar-Itzhack Memorial Symposium on Estimation, Navigation, and Spacecraft Control, Haifa, Israel, October 14 17, 2012 Track-to-Track Architectures A Review Xin Tian and Yaakov Bar-Shalom This
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 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 informationDead Reckoning navigation (DR navigation)
Dead Reckoning navigation (DR navigation) Prepared by A.Kaviyarasu Assistant Professor Department of Aerospace Engineering Madras Institute Of Technology Chromepet, Chennai A Navigation which uses a Inertial
More informationDesign of Nearly Constant Velocity Track Filters for Brief Maneuvers
4th International Conference on Information Fusion Chicago, Illinois, USA, July 5-8, 20 Design of Nearly Constant Velocity rack Filters for Brief Maneuvers W. Dale Blair Georgia ech Research Institute
More informationNavigation Mathematics: Kinematics (Earth Surface & Gravity Models) EE 570: Location and Navigation
Lecture Navigation Mathematics: Kinematics (Earth Surface & ) EE 570: Location and Navigation Lecture Notes Update on March 10, 2016 Aly El-Osery and Kevin Wedeward, Electrical Engineering Dept., New Mexico
More informationA Novel Maneuvering Target Tracking Algorithm for Radar/Infrared Sensors
Chinese Journal of Electronics Vol.19 No.4 Oct. 21 A Novel Maneuvering Target Tracking Algorithm for Radar/Infrared Sensors YIN Jihao 1 CUIBingzhe 2 and WANG Yifei 1 (1.School of Astronautics Beihang University
More informationMinimum Time Ascent Phase Trajectory Optimization using Steepest Descent Method
IJCTA, 9(39), 2016, pp. 71-76 International Science Press Closed Loop Control of Soft Switched Forward Converter Using Intelligent Controller 71 Minimum Time Ascent Phase Trajectory Optimization using
More informationFEATURE AIDED SWITCHING MODEL SET AP- PROACH FOR MANEUVERING TARGET TRACKING
Progress In Electromagnetics Research B, Vol. 45, 251 268, 2012 FEATURE AIDED SWITCHING MODEL SET AP- PROACH FOR MANEUVERING TARGET TRACKING J.-P. Fan, Y.-L. Zhu, S.-J. Fan *, H.-Q. Fan, and Q. Fu ATR
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 informationAdvanced ESM Angle Tracker: Volume I Theoretical Foundations
Naval Research Laboratory Washington, DC 075-50 NRL/FR/574--06-0,4 Advanced ESM Angle Tracker: Volume I Theoretical Foundations Edward N. Khoury Surface Electronic Warfare Systems Branch Tactical Electronic
More informationQuaternion based Extended Kalman Filter
Quaternion based Extended Kalman Filter, Sergio Montenegro About this lecture General introduction to rotations and quaternions. Introduction to Kalman Filter for Attitude Estimation How to implement and
More informationInitial Trajectory and Atmospheric Effects
Initial Trajectory and Atmospheric Effects G. Flanagan Alna Space Program July 13, 2011 Introduction A major consideration for an earth-based accelerator is atmospheric drag. Drag loses mean that the gun
More informationDATA FUSION III: Estimation Theory
DATA FUSION III: Estimation Theory Date: March 30, 2006 Time: 5:00 7:30 PM Location: B-300-2-3 (AAR-400) (Main Building, 2nd floor, near freight elevators) Instructor: Dr. James K Beard Credits: 1 Course
More informationTARGET TRACKING AND DATA FUSION: How to Get the Most Out of Your Sensors and make a living out of it FUSION 2017
TARGET TRACKING AND DATA FUSION: How to Get the Most Out of Your Sensors and make a living out of it AN OVERVIEW OF TRACKING ALGORITHMS FOR CLUTTERED AND MULTITARGET-MULTISENSOR ENVIRONMENTS Yaakov Bar-Shalom,
More informationIntegrated Estimator/Guidance Law Design for Improved Ballistic Missile Defense
Integrated Estimator/Guidance Law Design for Improved Ballistic Missile Defense Josef Shinar 2, Yaakov Oshman 3 and Vladimir Turetsky 4 Faculty of Aerospace Engineering Technion, Israel Institute of Technology,
More informationUsing North Korean Missile Development to Enhance Student Interest in Kinematics and Newtonian Mechanics
Using North Korean Missile Development to Enhance Student Interest in Kinematics and Newtonian Mechanics Kendall Mallory Aims Community College University of Northern Colorado Recent Events Missile Development
More informationLecture 2. Lecture 1. Forces on a rotating planet. We will describe the atmosphere and ocean in terms of their:
Lecture 2 Lecture 1 Forces on a rotating planet We will describe the atmosphere and ocean in terms of their: velocity u = (u,v,w) pressure P density ρ temperature T salinity S up For convenience, we will
More informationFIBER OPTIC GYRO-BASED ATTITUDE DETERMINATION FOR HIGH- PERFORMANCE TARGET TRACKING
FIBER OPTIC GYRO-BASED ATTITUDE DETERMINATION FOR HIGH- PERFORMANCE TARGET TRACKING Elias F. Solorzano University of Toronto (Space Flight Laboratory) Toronto, ON (Canada) August 10 th, 2016 30 th AIAA/USU
More informationSliding Window Test vs. Single Time Test for Track-to-Track Association
Sliding Window Test vs. Single Time Test for Track-to-Track Association Xin Tian Dept. of Electrical and Computer Engineering University of Connecticut Storrs, CT 06269-257, U.S.A. Email: xin.tian@engr.uconn.edu
More information5.12 The Aerodynamic Assist Trajectories of Vehicles Propelled by Solar Radiation Pressure References...
1 The Two-Body Problem... 1 1.1 Position of the Problem... 1 1.2 The Conic Sections and Their Geometrical Properties... 12 1.3 The Elliptic Orbits... 20 1.4 The Hyperbolic and Parabolic Trajectories...
More informationAS3010: Introduction to Space Technology
AS3010: Introduction to Space Technology L E C T U R E S 8-9 Part B, Lectures 8-9 23 March, 2017 C O N T E N T S In this lecture, we will look at factors that cause an orbit to change over time orbital
More informationINNOVATIVE STRATEGY FOR Z9 REENTRY
INNOVATIVE STRATEGY FOR Z9 REENTRY Gregor Martens*, Elena Vellutini**, Irene Cruciani* *ELV, Corso Garibaldi, 34 Colleferro (Italy) **Aizoon, Viale Città d Europa 681, 144, Roma (Italy) Abstract Large
More informationLecture Outline. Target Tracking: Lecture 3 Maneuvering Target Tracking Issues. Maneuver Illustration. Maneuver Illustration. Maneuver Detection
REGLERTEKNIK Lecture Outline AUTOMATIC CONTROL Target Tracking: Lecture 3 Maneuvering Target Tracking Issues Maneuver Detection Emre Özkan emre@isy.liu.se Division of Automatic Control Department of Electrical
More informationPrediction of ESTSP Competition Time Series by Unscented Kalman Filter and RTS Smoother
Prediction of ESTSP Competition Time Series by Unscented Kalman Filter and RTS Smoother Simo Särkkä, Aki Vehtari and Jouko Lampinen Helsinki University of Technology Department of Electrical and Communications
More informationModeling Realistic Projectile Motion
Modeling Realistic Projectile Motion Matthew Geske 1. Introduction The calculation of simple projectile motion is one of the most basic exercises in physics. Finding, for example, the range of a thrown
More informationThrust acceleration estimation using an on-line non-linear recursive least squares algorithm
7 Thrust acceleration estimation using an on-line non-linear recursive least squares algorithm N Ghahramani, A Naghash 2, and F Towhidkhah 2 Department of Aerospace Engineering, Amirkabir University of
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 informationImproved Space Target Tracking Through Bias Estimation From In-situ Celestial Observations
Improved Space Target Tracking Through Bias Estimation From In-situ Celestial Observations A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy at
More informationSimple Harmonic Motion. Harmonic motion due to a net restoring force directly proportional to the displacement
Simple Harmonic Motion Harmonic motion due to a net restoring force directly proportional to the displacement Eample: Spring motion: F = -k Net Force: d! k = m dt " F = ma d dt + k m = 0 Equation of motion
More informationTwo Levels Data Fusion Filtering Algorithms of Multimode Compound Seeker
TELKOMNKA, Vol.11, No.11, November 13, pp. 674~68 e-ssn: 87-78X 674 Two Levels Data Fusion Filtering Algorithms of Multimode Compound Seeker Guodong Zhang*, Zhong Liu, Ke Li College of Electronic Engineering,
More information= o + t = ot + ½ t 2 = o + 2
Chapters 8-9 Rotational Kinematics and Dynamics Rotational motion Rotational motion refers to the motion of an object or system that spins about an axis. The axis of rotation is the line about which the
More informationCHAPTER 3 PERFORMANCE
PERFORMANCE 3.1 Introduction The LM-3A performance figures given in this chapter are based on the following assumptions: Launching from XSLC (Xichang Satellite Launch Center, Sichuan Province, China),
More informationChapter 8. Dynamics II: Motion in a Plane
Chapter 8. Dynamics II: Motion in a Plane Chapter Goal: To learn how to solve problems about motion in a plane. Slide 8-2 Chapter 8 Preview Slide 8-3 Chapter 8 Preview Slide 8-4 Chapter 8 Preview Slide
More informationCombined Particle and Smooth Variable Structure Filtering for Nonlinear Estimation Problems
14th International Conference on Information Fusion Chicago, Illinois, USA, July 5-8, 2011 Combined Particle and Smooth Variable Structure Filtering for Nonlinear Estimation Problems S. Andrew Gadsden
More informationWhy do we care? Measurements. Handling uncertainty over time: predicting, estimating, recognizing, learning. Dealing with time
Handling uncertainty over time: predicting, estimating, recognizing, learning Chris Atkeson 2004 Why do we care? Speech recognition makes use of dependence of words and phonemes across time. Knowing where
More informationTechnology of Rocket
Technology of Rocket Parts of Rocket There are four major parts of rocket Structural system Propulsion system Guidance system Payload system Structural system The structural system of a rocket includes
More informationLow-Cost Integrated Navigation System
Institute of Flight System Dynamics Technische Universität München Prof. Dr.-Ing. Florian Holzapfel Eecuted at Robotics and Automation Lab Centre for Intelligent Information Processing Systems University
More informationSatellite Orbital Maneuvers and Transfers. Dr Ugur GUVEN
Satellite Orbital Maneuvers and Transfers Dr Ugur GUVEN Orbit Maneuvers At some point during the lifetime of most space vehicles or satellites, we must change one or more of the orbital elements. For example,
More informationInteractive Multiple Model Hazard States Prediction for Unmanned Aircraft Systems (UAS) Detect and Avoid (DAA)
Interactive Multiple Model Hazard States Prediction for Unmanned Aircraft Systems (UAS) Detect and Avoid (DAA) Adriano Canolla Michael B. Jamoom, and Boris Pervan Illinois Institute of Technology This
More informationA comparison of estimation accuracy by the use of KF, EKF & UKF filters
Computational Methods and Eperimental Measurements XIII 779 A comparison of estimation accurac b the use of KF EKF & UKF filters S. Konatowski & A. T. Pieniężn Department of Electronics Militar Universit
More informationTHE METEOROLOGICAL ROCKET SYSTEM FOR ATMOSPHERIC RESEARCH
THE METEOROLOGICAL ROCKET SYSTEM FOR ATMOSPHERIC RESEARCH Komissarenko Alexander I. (1) Kuznetsov Vladimir M. (1) Filippov Valerii V. (1) Ryndina Elena C. (1) (1) State Unitary Enterprise KBP Instrument
More informationDevelopment and Flight Testing of Energy Management Algorithms for Small-Scale Sounding Rockets
Development and Flight Testing of Energy Management Algorithms for Small-Scale Sounding Rockets and Shane Robinson The development, implementation, and ight results for a navigation algorithm and an energy
More informationState Estimation for Autopilot Control of Small Unmanned Aerial Vehicles in Windy Conditions
University of Colorado, Boulder CU Scholar Aerospace Engineering Sciences Graduate Theses & Dissertations Aerospace Engineering Sciences Summer 7-23-2014 State Estimation for Autopilot Control of Small
More informationEE 570: Location and Navigation
EE 570: Location and Navigation Navigation Mathematics: Kinematics (Earth Surface & Gravity Models) Aly El-Osery Kevin Wedeward Electrical Engineering Department, New Mexico Tech Socorro, New Mexico, USA
More informationOptimized Trajectory Shaping Guidance for an Air-to-Ground Missile Launched from a Gunship. Craig Phillips Ernie Ohlmeyer Shane Sorenson
Optimized Trajectory Shaping Guidance for an Air-to-Ground Missile Launched from a Gunship Craig Phillips Ernie Ohlmeyer Shane Sorenson Overview Mission Scenario Notional Munition Concept Guidance Laws
More informationUsing the Kalman Filter to Estimate the State of a Maneuvering Aircraft
1 Using the Kalman Filter to Estimate the State of a Maneuvering Aircraft K. Meier and A. Desai Abstract Using sensors that only measure the bearing angle and range of an aircraft, a Kalman filter is implemented
More informationLecture 7: Optimal Smoothing
Department of Biomedical Engineering and Computational Science Aalto University March 17, 2011 Contents 1 What is Optimal Smoothing? 2 Bayesian Optimal Smoothing Equations 3 Rauch-Tung-Striebel Smoother
More informationORBIT DETERMINATION THROUGH GLOBAL POSITIONING SYSTEMS: A LITERATURE SURVEY
ORBIT DETERMINATION THROUGH GLOBAL POSITIONING SYSTEMS: A LITERATURE SURVEY ABDUL MANARVI, MASTER S STUDENT DEPARTMENT OF AEROSPACE ENGINEERING EMBRY-RIDDLE AERONAUTICAL UNIVERSITY DAYTONA BEACH, FL, U.S.A
More informationEKF/UKF Maneuvering Target Tracking using Coordinated Turn Models with Polar/Cartesian Velocity
EKF/UKF Maneuvering Target Tracking using Coordinated Turn Models with Polar/Cartesian Velocity Michael Roth, Gustaf Hendeby, and Fredrik Gustafsson Dept. Electrical Engineering, Linköping University,
More informationEvaluation of different wind estimation methods in flight tests with a fixed-wing UAV
Evaluation of different wind estimation methods in flight tests with a fixed-wing UAV Julian Sören Lorenz February 5, 2018 Contents 1 Glossary 2 2 Introduction 3 3 Tested algorithms 3 3.1 Unfiltered Method
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 informationA Multiple Target Range and Range-Rate. Tracker Using an Extended Kalman Filter and. a Multilayered Association Scheme
A Multiple Target Range and Range-Rate Tracker Using an Extended Kalman Filter and a Multilayered Association Scheme A thesis submitted by Leah Uftring In partial fufillment of the degree requirements
More informationHigh-Power Rocketry. Calculating the motion of a rocket for purely vertical flight.
High-Power Rocketry Calculating the motion of a rocket for purely vertical flight. Phase I Boost phase: motor firing (rocket losing mass), going upwards faster and faster (accelerating upwards) Phase II
More informationCRB-based optimal radar placement for target positioning
CRB-based optimal radar placement for target positioning Kyungwoo Yoo and Joohwan Chun KAIST School of Electrical Engineering Yuseong-gu, Daejeon, Republic of Korea Email: babooovv@kaist.ac.kr Chungho
More informationPHYSICS 221, FALL 2009 EXAM #1 SOLUTIONS WEDNESDAY, SEPTEMBER 30, 2009
PHYSICS 221, FALL 2009 EXAM #1 SOLUTIONS WEDNESDAY, SEPTEMBER 30, 2009 Note: The unit vectors in the +x, +y, and +z directions of a right-handed Cartesian coordinate system are î, ĵ, and ˆk, respectively.
More informationCHAPTER 3 PERFORMANCE
PERFORMANCE 3.1 Introduction The LM-3B performance figures given in this chapter are based on the following assumptions: Launching from XSLC (Xichang Satellite Launch Center, Sichuan Province, China),
More informationOPTIMAL INTEGRATION OF ESTIMATION AND GUIDANCE FOR INTERCEPTORS
TECHNION RESEARCH AND DEVELOPMENT FOUNDATION LTD. OPTIMAL INTEGRATION OF ESTIMATION AND GUIDANCE FOR INTERCEPTORS AFOSR Contract No. F61775-01-C-0007 Final Technical Report covering the period of October
More informationArtificial Intelligence
Artificial Intelligence Roman Barták Department of Theoretical Computer Science and Mathematical Logic Summary of last lecture We know how to do probabilistic reasoning over time transition model P(X t
More informationOnboard Engine FDI in Autonomous Aircraft Using Stochastic Nonlinear Modelling of Flight Signal Dependencies
Onboard Engine FDI in Autonomous Aircraft Using Stochastic Nonlinear Modelling of Flight Signal Dependencies Dimitrios G. Dimogianopoulos, John D. Hios and Spilios D. Fassois Stochastic Mechanical Systems
More informationAdaptive 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 informationPerformance Prediction of Multisensor Tracking Systems for Single Maneuvering Targets
Performance Prediction of Multisensor Tracing Systems for Single Maneuvering Targets WILLIAM D. BLAIR PAUL A. MICELI Studying the performance of multisensor tracing systems against maneuvering targets
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 informationPHYSICS. Chapter 8 Lecture FOR SCIENTISTS AND ENGINEERS A STRATEGIC APPROACH 4/E RANDALL D. KNIGHT Pearson Education, Inc.
PHYSICS FOR SCIENTISTS AND ENGINEERS A STRATEGIC APPROACH 4/E Chapter 8 Lecture RANDALL D. KNIGHT Chapter 8. Dynamics II: Motion in a Plane IN THIS CHAPTER, you will learn to solve problems about motion
More informationAn Empirical Study into the Robustness of Split Covariance Addition (SCA) for Human Motion Tracking
An Empirical Study into the Robustness of Split Covariance Addition (SCA) for Human Motion Tracking Simon J. Julier and Joseph J. LaViola Jr. Abstract Accurate human body tracking is extremely important
More informationUlrich Walter. Astronautics. The Physics of Space Flight. 2nd, Enlarged and Improved Edition
Ulrich Walter Astronautics The Physics of Space Flight 2nd, Enlarged and Improved Edition Preface to Second Edition Preface XVII Acknowledgments XIX List of Symbols XXI XV 1 Rocket Fundamentals 1 1.1 Rocket
More informationAgile Missile Controller Based on Adaptive Nonlinear Backstepping Control
Agile Missile Controller Based on Adaptive Nonlinear Backstepping Control Chang-Hun Lee, Tae-Hun Kim and Min-Jea Tahk 3 Korea Advanced Institute of Science and Technology(KAIST), Daejeon, 305-70, Korea
More informationTime Series Prediction by Kalman Smoother with Cross-Validated Noise Density
Time Series Prediction by Kalman Smoother with Cross-Validated Noise Density Simo Särkkä E-mail: simo.sarkka@hut.fi Aki Vehtari E-mail: aki.vehtari@hut.fi Jouko Lampinen E-mail: jouko.lampinen@hut.fi Abstract
More informationOPTIMAL TIME TRANSFER
OPTIMAL TIME TRANSFER James R. Wright and James Woodburn Analytical Graphics, Inc. 220 Valley Creek Blvd, Exton, PA 19341, USA Abstract We have designed an optimal time transfer algorithm using GPS carrier-phase
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 information6. Wiener and Kalman Filters
6. Wiener and Kalman Filters 8 3 6.1. The Wiener Filter. The theory of filtering of stationary time series for a variety of purposes was constructed by Norbert Wiener in the 1940s for continuous time processes
More informationWIND INSTABILITY WHAT BARROWMAN LEFT OUT
WIND INSTABILITY WHAT BARROWMAN LEFT OUT by Robert Galejs galejs@ll.mit.edu In Centuri TIR-33 (reprinted in the March 98 issue of High Power Rocketry), Jim Barrowman outlined a method for the determination
More information3 Space curvilinear motion, motion in non-inertial frames
3 Space curvilinear motion, motion in non-inertial frames 3.1 In-class problem A rocket of initial mass m i is fired vertically up from earth and accelerates until its fuel is exhausted. The residual mass
More informationFundamentals of attitude Estimation
Fundamentals of attitude Estimation Prepared by A.Kaviyarasu Assistant Professor Department of Aerospace Engineering Madras Institute Of Technology Chromepet, Chennai Basically an IMU can used for two
More informationTwo dimensional rate gyro bias estimation for precise pitch and roll attitude determination utilizing a dual arc accelerometer array
Rochester Institute of Technology RIT Scholar Works Theses Thesis/Dissertation Collections -- Two dimensional rate gyro bias estimation for precise pitch and roll attitude determination utilizing a dual
More informationOrbital Anomaly Detection and Application of Space Object
Orbital Anomaly Detection and Application of Space Object Lei CHEN *, Xianzong BAI, Taotao ZHANG College of Aerospace and Materials Engineering, National University of Defense Technology, Changsha 473,China
More informationThe Shifted Rayleigh Filter for 3D Bearings-only Measurements with Clutter
The Shifted Rayleigh Filter for 3D Bearings-only Measurements with Clutter Attila Can Özelçi EEE Department Imperial College London, SW7 BT attila.ozelci@imperial.ac.uk Richard Vinter EEE Department Imperial
More informationUAV Navigation: Airborne Inertial SLAM
Introduction UAV Navigation: Airborne Inertial SLAM Jonghyuk Kim Faculty of Engineering and Information Technology Australian National University, Australia Salah Sukkarieh ARC Centre of Excellence in
More information