UAVBook Supplement Full State Direct and Indirect EKF
|
|
- William Holt
- 5 years ago
- Views:
Transcription
1 UAVBook Supplement Full State Direct and Indirect EKF Randal W. Beard March 14, 217 This supplement will explore alternatives to the state estimation scheme presented in the book. In particular, we will derive a full state estimation scheme using both the direct Kalman filter and the indirect Kalman filter. We begin by focusing on the longitudinal states of the aircraft. 1 Direct Kalman Filter for Longitudinal States 1.1 Aircraft and sensor models The model for the dynamics of the aircraft are similar to that presented in the book after removing the lateral states. In particular, Equations become ṗ n = u cos θ + w sin θ ḣ = u sin θ w cos θ u = qw + f x m ẇ = qu + f z m θ = q. We will assume that the sensors that are available are the y-axis gyroscope, the x- and z-axis accelerometers, the static and differential pressure 1
2 sensors, the GPS-north measurement, and the GPS-ground-speed measurement. The models for the sensors, adapted from Equations 7.3, 7.5, 7.9, 7.1, 7.18, and 7.21 are given by y accel,x = f x m + g sin θ + η accel,x y accel,z = f z m g cos θ + η accel,z y gyro,y = q + b y + η gyro,y y abs pres = ρgh AGL + η abs pres y diff pres = ρv a η diff pres y GPS,n = p n + ν n [n] y GPS,Vg = V a + w n + η V. where b y is the unknown bias on the y-axis gyroscope. 1.2 Propagation model for the EKF The longitudinal states will be defined as x = p n, h, u, w, θ, b y, w n. Given the fact that the accelerometers measure total force minus gravity, we can write the equations of motion for x as ṗ n = u cos θ + w sin θ ḣ = u sin θ w cos θ u = y gyro,y b y w g sin θ + y accel,x wη gyro,y η accel,x ẇ = y gyro,y b y u + g cos θ + y accel,z + uη gyro,y η accel,z θ = y gyro,y b y + η gyro,y ḃ y = ẇ n =. 2
3 and If we define u cos θ + w sin θ u sin θ w cos θ fx, y = y gyro,y b y w g sin θ + y accel,x y gyro,y b y u + g cos θ + y accel,z y gyro,y b y then we can write where and G = w 1 u 1, 1 ẋ = fx, y + Gxη s + η d, y gyro,y y = y accel,x y accel,z η gyro,y η s = η accel,x η accel,z is the noise associated with the gyros and accelerometers, and where η i is a noise term associated with model uncertainty. The Jacobian of f can be shown to be cos θ sin θ [ u sin θ + w cos θ] sin θ cos θ [u cos θ + w sin θ] Ax, y = f x = y gyro,y b y g cos θ w y gyro,y b y g sin θ u. 1 3
4 The covariance of the noise term Gη s + η d is given by E[Gη s + η i Gη s + η i ] = GQ s G + Q d, where Q s = diag [σgyro,y, 2 σaccel 2 x, σgyro 2 z ], Q d = diag [σp 2 n, σh, 2 σu, 2 σv, 2 σθ, 2 σb 2 y, σw 2 n ], and where Q d are tuning parameters. 1.3 Sensor Update Equations Static Pressure Sensor The model for the static pressure sensor is Therefore, the Jacobian is given by h static x = ρgh. C static x = ρg Differential Pressure Sensor The model for the differential pressure sensor is h diff x = 1 2 ρv 2 a, where Va 2 = u 2 r + wr, 2 and where u r and w r are the wind relative velocities given at the top of page 57 in the uavbook. For the longitudinal case where w d = we can write the relative velocities as Therefore u r = u w n cos θ w r = w w n sin θ. h diff x = 1 2 ρu w n cos θ 2 + w w n sin θ 2, with an associated Jacobian of C diff x = ρ u r w r u r w n sin θ w r w n cos θ u r cos θ w r sin θ. 4
5 1.3.3 GPS North sensor The model for the GPS north sensor is with associated Jacobian GPS ground-speed sensor h GPS,n x = p n C GPS,n x = 1. The ground-speed is the projection of the inertial velocity vector u, w on to the ground plane. Therefore with associated Jacobian h GPS,Vg x = u cos θ + w sin θ C GPS,Vg x = cos θ sin θ u sin θ + w cos θ. 1.4 Algorithm for Direct Kalman Filter We assume that the gyro, accelerometers, and pressure sensors are sampled at rate T s, which is the same rate that state estimates are required by the controller. We also assume that GPS is sampled at T GPS >> T s. The algorithm for the direct Kalman filter is given as follows.. Initialize Filter. Set ˆx = p n, h, V a,,,,, P = diage 2 p n, e 2 h, e 2 u, e 2 w, e 2 θ, e 2 b y, e 2 w n, where e is the standard deviation of the expected error in. 1. At the fast sample rate T s 1.a. Propagate ˆx and P according to ˆx = fˆx, y P = Aˆx, yp + P A ˆx, y + GˆxQ s Q ˆx + Q d 5
6 1.b. Update ˆx and P with the static pressure sensor according to L = P C staticˆx /σ 2 abs pres + C static ˆx P C staticˆx P + = I LC static ˆx P ˆx + = ˆx + Ly abs pres h static ˆx. 1.c. Update ˆx and P with the differential pressure sensor according to L = P C diffˆx /σ 2 diff pres + C diff ˆx P C diffˆx P + = I LC diff ˆx P ˆx + = ˆx + Ly diff pres h diff ˆx. 2. When GPS measurements are received at T GPS : 2.a. Update ˆx and P with the GPS north sensor according to L = P C GPS,nˆx /σ 2 GPS,n + C GPS,n ˆx P C GPS,nˆx P + = I LC GPS,n ˆx P ˆx + = ˆx + Ly GPS,n h GPS,n ˆx. 2.b. Update ˆx and P with the GPS ground-speed sensor according to L = P C GPS,V g ˆx /σ 2 GPS,V g + C GPS,Vg ˆx P C GPS,V g ˆx P + = I LC GPS,Vg ˆx P ˆx + = ˆx + Ly GPS,Vg h GPS,Vg ˆx. Results for the direct Kalman filter are shown in Figure 1. Red is the commanded signal, blue is truth, and green is the estimate. It can be seen from the figure that, with the possible exception of the ground-speed V g, the filter performs well. The simulation files associated with this example are included on the uavbook website. 6
7 25 2 p n h V g V a α b y w n θ δ t q δ e Figure 1: Estimation results for the direct Kalman filter. 7
8 2 Indirect Kalman Filter for Longitudinal States For the indirect Kalman filter, the idea is to filter the error states, which should satisfy the linear and Gaussian assumptions better than the direct implementation. Let x be the true state, ˆx the estimated state, and let x = x ˆx be the error state. If x satisfies the differential equation ẋ = fx, y + Gxη s + η i where all definitions are identical to the previous section, and if ˆx satisfies ˆx = fˆx, y, then the error state evolves according to x = fx, y + Gxη s + η i fˆx, y. Using the Taylor series expansion up to the linear term fx, y fˆx, y + Aˆx, y x Gxη s Gˆxη s gives x = Aˆx, y x + Gˆxη s + η i. Note that the covariance of the noise term E[Gˆxη s + η i Gˆxη s + η i ] = GˆxQ s Gˆx + Q i, where Q i is a tuning parameter similar to Q d for the direct EKF. The basic idea is to run the Kalman filter on the error state, and then to correct the state estimate at each step by adding in the error. The error state is reset to zero after each correction. The algorithm is outlined below. 2.1 Algorithm for the Indirect Kalman Filter The algorithm for the indirect Kalman filter is given as follows. Terms that differ from the direct Kalman filter are highlighted in blue.. Initialize Filter..1 Initialize the state estimate according to ˆx = p n, h, V a,,,,. 8
9 .2 Initialize the error state estimate according to x=,,,,,,, P = diage 2 p n, e 2 h, e 2 u, e 2 w, e 2 θ, e 2 b y, e 2 w n, where e is the standard deviation of the expected error in. 1. At the fast sample rate T s 1.a. Propagate the estimated state ˆx according to 1.b. Propagate x and P according to ˆx = fˆx, y x= Aˆx, y x P = Aˆx, yp + P A ˆx, y + GˆxQ s Q ˆx + Q d 1.c. Update x and P with the static pressure sensor according to L = P C staticˆx /σ 2 abs pres + C static ˆx P C staticˆx P + = I LC static ˆx P x + = x + Ly abs pres h static ˆx C static ˆx x. 1.d. Update x and P with the differential pressure sensor according to L = P C diffˆx /σ 2 diff pres + C diff ˆx P C diffˆx P + = I LC diff ˆx P x + = x + Ly diff pres h diff ˆx C diff ˆx x. 2. When GPS measurements are received at T GPS : 2.a. Update x and P with the GPS north sensor according to L = P C GPS,nˆx /σ 2 GPS,n + C GPS,n ˆx P C GPS,nˆx P + = I LC GPS,n ˆx P x + = x + Ly GPS,n h GPS,n ˆx C GPS,n ˆx x. 2.b. Update x and P with the GPS ground-speed sensor according to L = P C GPS,V g ˆx /σ 2 GPS,V g + C GPS,Vg ˆx P C GPS,V g ˆx P + = I LC GPS,Vg ˆx P x + = x + Ly GPS,Vg h GPS,Vg ˆx C GPS,Vg ˆx x. 9
10 2.c. Update ˆx according to ˆx + = ˆx + x +. 2.d. Set x =,,,,,,. Results for the indirect Kalman filter are shown in Figure 2. Red is the commanded signal, blue is truth, and green is the estimate. It can be seen from the figure that, with the possible exception of the ground-speed V g, the filter performs almost identical to the direct KF. We should note that the results are obtained when the estimator is in the feedback loop. The simulation files associated with this example are included on the uavbook website. 2 p n h V g θ q V a α b y w n 5 δ t δ e Figure 2: Estimation results for the indirect Kalman filter. 1
11 3 Direct Kalman Filter for Full States In this section we will expand upon the previous discussion to include the full aircraft state. 3.1 Aircraft and sensor models The model for the dynamics of the aircraft are similar to that presented in the book. If we define p = p n, p e, p d, v = u, v, w, Θ = φ, θ, ψ, ω = p, q, r, then we can write the dynamics as ṗ = RΘv v = v ω + a + Rg Ω = SΘω ḃ = 3 1 ẇ = 2 1, where we have used the notation a c b b = c a, c b a and where c θ c ψ c θ s ψ s θ RΩ = s φ s θ c ψ c φ s ψ s φ s θ s ψ + c φ c ψ s φ c θ c φ s θ c ψ + s φ s ψ c φ s θ s ψ s φ c ψ c φ c θ 1 sin φ tan θ cos φ tan θ SΩ = cos φ sin φ, sin φ sec θ cos φ sec θ and where we have used the model for the accelerometer in Equation 7.1 to get 1 m f = a + R Θg, where g =,, g. 11
12 We will assume that the sensors that are available are the gyroscopes, the accelerometers, the static and differential pressure sensors, the GPS north and east measurements, and the GPS ground-speed and course measurement. The models for the sensors are given by y accel = a + η accel y gyro = ω + b + η gyro y abs pres = ρgh AGL + η abs pres y diff pres = ρv a η diff pres y GPS,n = p n + ν n [n] y GPS,e = p e + ν e [n] = V a cos ψ + w n 2 + V a sin ψ + w e 2 + η V y GPS,χ = atan2v a sin ψ + w e, V a cos ψ + w n + η χ. y GPS,Vg 3.2 Propagation model for the EKF Using the gyro and accelerometer models, the equations of motion can be written as Defining ṗ = RΘv v = v y gyro b η gyro + y accel η accel + R Θg Ω = SΘy gyro b η gyro ḃ = 3 1 ẇ = 2 1. x = p, v, Θ, b, w y = y accel, y gyro, we get ẋ = fx, y + G g xη gyro + G a η accel + η, 12
13 where RΘv v y gyro b + y accel + R Θg fx, y = SΘy gyro b , 3 3 v G g x = SΘ G a = I , and where η is the general process noise associated with the model uncertainty and assumed to have covariance Q. The Jacobian of f is [RΘv] 3 3 RΘ Θ Ax, y = f x = 3 3 y gyro b [R Θg] v Θ 3 2 [SΘy 3 3 gyro b] 3 3 SΘ Θ It is straightforward to show that when z R 3 and Az is a 3 3 matrix whose elements are functions of z, then for any w R 3 [Azw] z = Az z 1 w, Az z 2 w, Az z 3 w. 13
14 Therefore [RΘv] Θ [ R Θg ] Θ [SΘy gyro b] Θ = = = RΘ φ v, RΘ θ v, R Θ φ g, R Θ θ g, SΘ y φ gyro b, RΘ ψ v R Θ g ψ SΘ y θ gyro b, cos θ = g cos θ cos φ sin θ sin φ cosθ sin φ sin θ cos φ SΘ ψ gyro b, where RΘ φ is the matrix where each element of RΘ has been differentiated by φ. The covariance of the noise term G g η gyro + G a η accel + η is E[G g η gyro +G a η accel +ηg g η gyro +G a η accel +η ] = G g Q gyro G g +G a Q accel G a +Q, where and where Q gyro = diag[σ 2 gyro,x, σ 2 gyro y, σ 2 gyro z ] Q accel = diag[σ 2 accel,x, σ 2 accel y, σ 2 accel z ], Q = diag[σ 2 p n, σ 2 p e, σ 2 p d, σ 2 u, σ 2 v, σ 2 w, σ 2 φ, σ 2 θ, σ 2 ψ, σ 2 b x, σ 2 b y, σ 2 b z, σ 2 w n, σ 2 w e ] are tuning parameters. 3.3 Sensor Update Equations Static Pressure Sensor The model for the static pressure sensor is Therefore, the Jacobian is given by h static x = ρgh = ρgp d. C static x =,, ρg,,,,,,,,,,,. 14
15 3.3.2 Differential Pressure Sensor The model for the differential pressure sensor is h diff x = 1 2 ρv 2 a, where Va 2 = v R Θw v R Θw, and where w = w n, w e,. Therefore h diff x = 1 2 ρ v R Θw v R Θw. The associated Jacobian is given by C diff x = ρ GPS North sensor [R Θw] Θ The model for the GPS north sensor is 3 1 v R Θw v R Θw 3 1 I 2 2, 2 1 RΘ v R Θw. h GPS,n x = p n with associated Jacobian C GPS,n x = 1,,,,,,,,,,,,, GPS East sensor The model for the GPS east sensor is with associated Jacobian h GPS,e x = p e C GPS,e x =, 1,,,,,,,,,,,,. 15
16 3.3.5 GPS ground-speed sensor The inertial velocity in the world frame, projected onto the north-east plane is given by v g = [ I 2 2, 2 1 ] RΘv. 1 The ground-speed is given by V g v, Θ = v g. The model for the sensor is therefore where P = [ I 2 2, C GPS,Vg x = h GPS,Vg x = P RΘv, ] 2 1, and the associated Jacobian is given by 1 3, v R ΘP V P RΘ, gv,θ Θ, 1 3, 1 2, where Vgv,Θ Θ is computed numerically as V g v, Θ = V gv, Θ + e i V g v, Θ, Θ i where e i is the 3 1 vector with one in the i th location and zeros elsewhere, and is a small number GPS course The course angle is the direction of v g = v g n, v g e given in Equation 1. Therefore, the course angle is given by The model for the sensor is therefore χv, Θ = atan2v g e, v g n. h GPS,χ x = atan2v g e, v g n, and the associated Jacobian is given by C GPS,χ x = 1 3, where χv,v v χv,v v, χv,θ Θ, 1 3, 1 2, and χv,θ Θ are computed numerically as χv, Θ = χv + e i, Θ χv, Θ v i χv, Θ = χv, Θ + e i χv, Θ. Θ i 16
17 3.3.7 Pseudo Measurement for Side Slip Angle With the given sensors, the side-slip angle, or side-to-side velocity is not observable and can therefore drift. To help correct this situation, we can add a pseudo measurement on the side-slip angle by assuming that it is zero. The side slip angle is given by β = v r V a, therefore an equivalent condition is that v r =. The airspeed vector is given by v a = v RΘ w, which implies that v r v, Θ, w = [ 1 ] v RΘ w. Therefore, the model for the pseudo-sensor is given by h β x = v r v, Θ, w, and the associated Jacobian is given by C β x = 1 3, v rv,θ,w v, v rv,θ,w Θ, 1 3, where the partial derivatives are computed numerically as v r v, Θ, w = v rv + e i, Θ, w v r v, Θ, w v i v r v, Θ, w = v rv, Θ + e i, w v r v, Θ, w Θ i v r v, Θ, w = v rv, Θ, w + e i v r v, Θ, w. w i The measurement used in the correction step is y β =. 3.4 Algorithm for Direct Kalman Filter v rv,θ,w w,, We assume that the gyro, accelerometers, and pressure sensors are sampled at rate T s, which is the same rate that state estimates are required by the controller. We also assume that GPS is sampled at T GPS >> T s. 17
18 The algorithm for the direct Kalman filter is given as follows.. Initialize Filter. Set ˆx = p n, p e, p d, V a,,,,, ψ,,,,,, P = diag [e 2 p n, e 2 p e, e 2 p d, e 2 u, e 2 v, e 2 w, e 2 φ, e 2 θ, e 2 ψ, e 2 b x, e 2 b y, e 2 b z, e 2 w n, e 2 w e ], where e is the standard deviation of the expected error in. 1. At the fast sample rate T s 1.a. Propagate ˆx and P according to ˆx = fˆx, y P = Aˆx, yp + P A ˆx, y + G g ˆxQ gyros G g ˆx + G a Q accel G a + Q 1.b. Update ˆx and P with the static pressure sensor according to L = P C staticˆx /σ 2 abs pres + C static ˆx P C staticˆx P + = I LC static ˆx P ˆx + = ˆx + Ly abs pres h static ˆx. 1.c. Update ˆx and P with the differential pressure sensor according to L = P C diffˆx /σ 2 diff pres + C diff ˆx P C diffˆx P + = I LC diff ˆx P ˆx + = ˆx + Ly diff pres h diff ˆx. 1.d. Update ˆx and P with the side-slip pseudo measurement according to L = P C β ˆx /σ 2 β + C β ˆx P C β ˆx P + = I LC β ˆx P ˆx + = ˆx + L h β ˆx. 2. When GPS measurements are received at T GPS : 2.a. Update ˆx and P with the GPS north measurement according to L = P C GPS,nˆx /σ 2 GPS,n + C GPS,n ˆx P C GPS,nˆx P + = I LC GPS,n ˆx P ˆx + = ˆx + Ly GPS,n h GPS,n ˆx. 18
19 2.b. Update ˆx and P with the GPS east measurement according to L = P CGPS,eˆx /σgps,e 2 + C GPS,e ˆx P CGPS,eˆx P + = I LC GPS,e ˆx P ˆx + = ˆx + Ly GPS,e h GPS,e ˆx. 2.c. Update ˆx and P with the GPS groundspeed measurement according to L = P CGPS,V g ˆx /σgps,v 2 g + C GPS,Vg ˆx P CGPS,V g ˆx P + = I LC GPS,Vg ˆx P ˆx + = ˆx + Ly GPS,Vg h GPS,Vg ˆx. 2.d. Update ˆx and P with the GPS course measurement according to L = P CGPS,χˆx /σgps,χ 2 + C GPS,χ ˆx P CGPS,χˆx P + = I LC GPS,χ ˆx P ˆx + = ˆx + Ly GPS,χ h GPS,V χ ˆx. 19
20 4 Indirect Kalman Filter for Full States 4.1 Algorithm for the Indirect Kalman Filter The algorithm for the indirect Kalman filter is given as follows. Terms that differ from the direct Kalman filter are highlighted in blue.. Initialize Filter..1 Initialize the state estimate according to ˆx = p n, p e, p d, V a,,,,, ψ,,,,,..2 Initialize the error state estimate according to x=,,,,,,,,,,,,,, P = diag [e 2 p n, e 2 p e, e 2 p d, e 2 u, e 2 v, e 2 w, e 2 φ, e 2 θ, e 2 ψ, e 2 b x, e 2 b y, e 2 b z, e 2 w n, e 2 w e ], where e is the standard deviation of the expected error in. 1. At the fast sample rate T s 1.a. Propagate the estimated state ˆx according to 1.b. Propagate x and P according to ˆx = fˆx, y x= Aˆx, y x P = Aˆx, yp + P A ˆx, y + G g ˆxQ gyros G g ˆx + G a Q accel G a + Q 1.c. Update x and P with the static pressure sensor according to L = P C staticˆx /σ 2 abs pres + C static ˆx P C staticˆx P + = I LC static ˆx P x + = x + Ly abs pres h static ˆx C static ˆx x. 1.d. Update x and P with the differential pressure sensor according to L = P C diffˆx /σ 2 diff pres + C diff ˆx P C diffˆx P + = I LC diff ˆx P x + = x + Ly diff pres h diff ˆx C diff ˆx x. 2
21 1.e. Update ˆx and P with the side-slip pseudo measurement according to L = P C β ˆx /σ 2 β + C β ˆx P C β ˆx P + = I LC β ˆx P x + = x + L h β ˆx C β ˆx x. 2. When GPS measurements are received at T GPS : 2.a. Update ˆx and P with the GPS north measurement according to L = P C GPS,nˆx /σ 2 GPS,n + C GPS,n ˆx P C GPS,nˆx P + = I LC GPS,n ˆx P x + = x + Ly GPS,n h GPS,n ˆx C GPS,n ˆx x. 2.b. Update ˆx and P with the GPS east measurement according to L = P C GPS,eˆx /σ 2 GPS,e + C GPS,e ˆx P C GPS,eˆx P + = I LC GPS,e ˆx P x + = x + Ly GPS,e h GPS,e ˆx C GPS,e ˆx x. 2.c. Update ˆx and P with the GPS groundspeed measurement according to L = P C GPS,V g ˆx /σ 2 GPS,V g + C GPS,Vg ˆx P C GPS,V g ˆx P + = I LC GPS,Vg ˆx P x + = x + Ly GPS,Vg h GPS,Vg ˆx C GPS,Vg ˆx x. 2.d. Update ˆx and P with the GPS course measurement according to L = P C GPS,χˆx /σ 2 GPS,χ + C GPS,χ ˆx P C GPS,χˆx P + = I LC GPS,χ ˆx P x + = x + Ly GPS,χ h GPS,V χ ˆx C GPS,χ ˆx x. 2.e. Update ˆx according to ˆx + = ˆx + x +. 2.f. Set x =,,,,,,,,,,,,,. References 21
1 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 informationApplication of state observers in attitude estimation using low-cost sensors
Application of state observers in attitude estimation using low-cost sensors Martin Řezáč Czech Technical University in Prague, Czech Republic March 26, 212 Introduction motivation for inertial estimation
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 informationEquations of Motion for Micro Air Vehicles
Equations of Motion for Micro Air Vehicles September, 005 Randal W. Beard, Associate Professor Department of Electrical and Computer Engineering Brigham Young University Provo, Utah 84604 USA voice: (80
More informationLecture. Aided INS EE 570: Location and Navigation. 1 Overview. 1.1 ECEF as and Example. 1.2 Inertial Measurements
Lecture Aided EE 570: Location and Navigation Lecture Notes Update on April 13, 2016 Aly El-Osery and Kevin Wedeward, Electrical Engineering Dept., New Mexico Tech In collaoration with Stephen Bruder,
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 informationRao-Blackwellized Particle Filtering for 6-DOF Estimation of Attitude and Position via GPS and Inertial Sensors
Rao-Blackwellized Particle Filtering for 6-DOF Estimation of Attitude and Position via GPS and Inertial Sensors GRASP Laboratory University of Pennsylvania June 6, 06 Outline Motivation Motivation 3 Problem
More informationEE 570: Location and Navigation
EE 570: Location and Navigation Aided INS Aly El-Osery Kevin Wedeward Electrical Engineering Department, New Mexico Tech Socorro, New Mexico, USA In Collaboration with Stephen Bruder Electrical and Computer
More informationChapter 10. Path Following. Beard & McLain, Small Unmanned Aircraft, Princeton University Press, 2012, Chapter 10, Slide 1
Chapter 10 Path Following Beard & McLain, Small Unmanned Aircraft, Princeton University Press, 2012, Chapter 10, Slide 1 Control Architecture destination, obstacles map path planner waypoints status path
More informationInvariant Extended Kalman Filter: Theory and application to a velocity-aided estimation problem
Invariant Extene Kalman Filter: Theory an application to a velocity-aie estimation problem S. Bonnabel (Mines ParisTech) Joint work with P. Martin (Mines ParisTech) E. Salaun (Georgia Institute of Technology)
More informationHomework: Ball Predictor for Robot Soccer
Homework: Ball Predictor for Robot Soccer ECEn 483 November 17, 2004 1 Introduction The objective of this homework is to prepare you to implement a ball predictor in the lab An extended Kalman filter will
More informationQuadrotor Dynamics and Control Rev 0.1
Brigham Young University BYU ScholarsArchive All Faculty Publications 28-5-5 Quadrotor Dynamics and Control Rev.1 Randal Beard beard@byu.edu Follow this and additional works at: https://scholarsarchive.byu.edu/facpub
More informationUAV Coordinate Frames and Rigid Body Dynamics
Brigham Young University BYU ScholarsArchive All Faculty Publications 24-- UAV oordinate Frames and Rigid Body Dynamics Randal Beard beard@byu.edu Follow this and additional works at: https://scholarsarchive.byu.edu/facpub
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 informationA Study of Covariances within Basic and Extended Kalman Filters
A Study of Covariances within Basic and Extended Kalman Filters David Wheeler Kyle Ingersoll December 2, 2013 Abstract This paper explores the role of covariance in the context of Kalman filters. The underlying
More informationExtension 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 informationNeural Networks Lecture 10: Fault Detection and Isolation (FDI) Using Neural Networks
Neural Networks Lecture 10: Fault Detection and Isolation (FDI) Using Neural Networks H.A. Talebi Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Winter 2011.
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 informationBayes Filter Reminder. Kalman Filter Localization. Properties of Gaussians. Gaussians. Prediction. Correction. σ 2. Univariate. 1 2πσ e.
Kalman Filter Localization Bayes Filter Reminder Prediction Correction Gaussians p(x) ~ N(µ,σ 2 ) : Properties of Gaussians Univariate p(x) = 1 1 2πσ e 2 (x µ) 2 σ 2 µ Univariate -σ σ Multivariate µ Multivariate
More informationTSRT14: Sensor Fusion Lecture 6. Kalman Filter (KF) Le 6: Kalman filter (KF), approximations (EKF, UKF) Lecture 5: summary
TSRT14 Lecture 6 Gustaf Hendeby Spring 217 1 / 42 Le 6: Kalman filter KF approximations EKF UKF TSRT14: Sensor Fusion Lecture 6 Kalman filter KF KF approximations EKF UKF Gustaf Hendeby hendeby@isyliuse
More informationChapter 1. Introduction. 1.1 System Architecture
Chapter 1 Introduction 1.1 System Architecture The objective of this book is to prepare the reader to do research in the exciting and rapidly developing field of autonomous navigation, guidance, and control
More informationAutonomous Mobile Robot Design
Autonomous Mobile Robot Design Topic: Extended Kalman Filter Dr. Kostas Alexis (CSE) These slides relied on the lectures from C. Stachniss, J. Sturm and the book Probabilistic Robotics from Thurn et al.
More informationLecture 9: Modeling and motion models
Sensor Fusion, 2014 Lecture 9: 1 Lecture 9: Modeling and motion models Whiteboard: Principles and some examples. Slides: Sampling formulas. Noise models. Standard motion models. Position as integrated
More informationKinematics for a Three Wheeled Mobile Robot
Kinematics for a Three Wheeled Mobile Robot Randal W. Beard Updated: March 13, 214 1 Reference Frames and 2D Rotations î 1 y î 2 y w 1 y w w 2 y î 2 x w 2 x w 1 x î 1 x Figure 1: The vector w can be expressed
More informationEE565:Mobile Robotics Lecture 6
EE565:Mobile Robotics Lecture 6 Welcome Dr. Ahmad Kamal Nasir Announcement Mid-Term Examination # 1 (25%) Understand basic wheel robot kinematics, common mobile robot sensors and actuators knowledge. Understand
More informationInertial Odometry using AR Drone s IMU and calculating measurement s covariance
Inertial Odometry using AR Drone s IMU and calculating measurement s covariance Welcome Lab 6 Dr. Ahmad Kamal Nasir 25.02.2015 Dr. Ahmad Kamal Nasir 1 Today s Objectives Introduction to AR-Drone On-board
More informationLecture AC-1. Aircraft Dynamics. Copy right 2003 by Jon at h an H ow
Lecture AC-1 Aircraft Dynamics Copy right 23 by Jon at h an H ow 1 Spring 23 16.61 AC 1 2 Aircraft Dynamics First note that it is possible to develop a very good approximation of a key motion of an aircraft
More informationSmartphone sensor based orientation determination for indoor navigation
Smartphone sensor based orientation determination for indoor naviation LBS Conference 15.11.2016 Andreas Ettliner Research Group Enineerin Geodesy Contact: andreas.ettliner@tuwien.ac.at Outline Motivation
More informationAttitude 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 informationChapter 9. Nonlinear Design Models. Beard & McLain, Small Unmanned Aircraft, Princeton University Press, 2012, Chapter 9, Slide 1
Chapter 9 Nonlinear Design Models Beard & McLain, Small Unmanned Aircraft, Princeton University Press, 2012, Chapter 9, Slide 1 Architecture Destination, obstacles Waypoints Path Definition Airspeed, Altitude,
More informationProblem 1: Ship Path-Following Control System (35%)
Problem 1: Ship Path-Following Control System (35%) Consider the kinematic equations: Figure 1: NTNU s research vessel, R/V Gunnerus, and Nomoto model: T ṙ + r = Kδ (1) with T = 22.0 s and K = 0.1 s 1.
More informationMAE 142 Homework #5 Due Friday, March 13, 2009
MAE 142 Homework #5 Due Friday, March 13, 2009 Please read through the entire homework set before beginning. Also, please label clearly your answers and summarize your findings as concisely as possible.
More 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 informationOptimization-Based Control
Optimization-Based Control Richard M. Murray Control and Dynamical Systems California Institute of Technology DRAFT v1.7a, 19 February 2008 c California Institute of Technology All rights reserved. This
More informationComputer Problem 1: SIE Guidance, Navigation, and Control
Computer Problem 1: SIE 39 - Guidance, Navigation, and Control Roger Skjetne March 12, 23 1 Problem 1 (DSRV) We have the model: m Zẇ Z q ẇ Mẇ I y M q q + ẋ U cos θ + w sin θ ż U sin θ + w cos θ θ q Zw
More informationVN-100 Velocity Compensation
VN-100 Velocity Compensation Velocity / Airspeed Aiding for AHRS Applications Application Note Abstract This application note describes how the VN-100 can be used in non-stationary applications which require
More informationForecasting and data assimilation
Supported by the National Science Foundation DMS Forecasting and data assimilation Outline Numerical models Kalman Filter Ensembles Douglas Nychka, Thomas Bengtsson, Chris Snyder Geophysical Statistics
More informationChapter 4 State Estimation
Chapter 4 State Estimation Navigation of an unmanned vehicle, always depends on a good estimation of the vehicle states. Especially if no external sensors or marers are available, more or less complex
More informationInvestigation 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 informationWe provide two sections from the book (in preparation) Intelligent and Autonomous Road Vehicles, by Ozguner, Acarman and Redmill.
We provide two sections from the book (in preparation) Intelligent and Autonomous Road Vehicles, by Ozguner, Acarman and Redmill. 2.3.2. Steering control using point mass model: Open loop commands We consider
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 informationA Stochastic Observability Test for Discrete-Time Kalman Filters
A Stochastic Observability Test for Discrete-Time Kalman Filters Vibhor L. Bageshwar 1, Demoz Gebre-Egziabher 2, William L. Garrard 3, and Tryphon T. Georgiou 4 University of Minnesota Minneapolis, MN
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 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 informationSatellite 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 informationTopography and Force Estimation in Atomic Force Microscopy by State and Parameter Estimation
Topography and Force Estimation in Atomic Force Microscopy by State and Parameter Estimation Michael R. P. Ragazzon, J. Tommy Gravdahl, Kristin Y. Pettersen, Arnfinn A. Eielsen Department of Engineering
More informationEstimation and Control of a Quadrotor Attitude
Estimation and Control of a Quadrotor Attitude Bernardo Sousa Machado Henriques Mechanical Engineering Department, Instituto Superior Técnico, Lisboa, Portugal E-mail: henriquesbernardo@gmail.com Abstract
More informationKalman Filter. Predict: Update: x k k 1 = F k x k 1 k 1 + B k u k P k k 1 = F k P k 1 k 1 F T k + Q
Kalman Filter Kalman Filter Predict: x k k 1 = F k x k 1 k 1 + B k u k P k k 1 = F k P k 1 k 1 F T k + Q Update: K = P k k 1 Hk T (H k P k k 1 Hk T + R) 1 x k k = x k k 1 + K(z k H k x k k 1 ) P k k =(I
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 informationA Low-Cost GPS Aided Inertial Navigation System for Vehicular Applications
A Low-Cost GPS Aided Inertial Navigation System for Vehicular Applications ISAAC SKOG Master of Science Thesis Stockholm, Sweden 2005-03-09 IR-SB-EX-0506 1 Abstract In this report an approach for integration
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 informationCubature Particle filter applied in a tightly-coupled GPS/INS navigation system
Cubature Particle filter applied in a tightly-coupled GPS/INS navigation system Yingwei Zhao & David Becker Physical and Satellite Geodesy Institute of Geodesy TU Darmstadt 1 Yingwei Zhao & David Becker
More informationFlight Test Results for Circular Path Following by Model Predictive Control
Preprints of the 19th World Congress The International Federation of Automatic Control Flight Test Results for Circular Path Following by Model Predictive Control Yoshiro Hamada Taro Tsukamoto Shinji Ishimoto
More informationIntegration of a strapdown gravimeter system in an Autonomous Underwater Vehicle
Integration of a strapdown gravimeter system in an Autonomous Underwater Vehicle Clément ROUSSEL PhD - Student (L2G - Le Mans - FRANCE) April 17, 2015 Clément ROUSSEL ISPRS / CIPA Workshop April 17, 2015
More informationFrom Bayes to Extended Kalman Filter
From Bayes to Extended Kalman Filter Michal Reinštein Czech Technical University in Prague Faculty of Electrical Engineering, Department of Cybernetics Center for Machine Perception http://cmp.felk.cvut.cz/
More informationMEMS Gyroscope Control Systems for Direct Angle Measurements
MEMS Gyroscope Control Systems for Direct Angle Measurements Chien-Yu Chi Mechanical Engineering National Chiao Tung University Hsin-Chu, Taiwan (R.O.C.) 3 Email: chienyu.me93g@nctu.edu.tw Tsung-Lin Chen
More informationwith Application to Autonomous Vehicles
Nonlinear with Application to Autonomous Vehicles (Ph.D. Candidate) C. Silvestre (Supervisor) P. Oliveira (Co-supervisor) Institute for s and Robotics Instituto Superior Técnico Portugal January 2010 Presentation
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 informationOptimization-Based Control
Optimization-Based Control Richard M. Murray Control and Dynamical Systems California Institute of Technology DRAFT v2.1a, January 3, 2010 c California Institute of Technology All rights reserved. This
More informationSIMULTANEOUS STATE AND PARAMETER ESTIMATION USING KALMAN FILTERS
ECE5550: Applied Kalman Filtering 9 1 SIMULTANEOUS STATE AND PARAMETER ESTIMATION USING KALMAN FILTERS 9.1: Parameters versus states Until now, we have assumed that the state-space model of the system
More informationJoint GPS and Vision Estimation Using an Adaptive Filter
1 Joint GPS and Vision Estimation Using an Adaptive Filter Shubhendra Vikram Singh Chauhan and Grace Xingxin Gao, University of Illinois at Urbana-Champaign Shubhendra Vikram Singh Chauhan received his
More informationIsobath following using an altimeter as a unique exteroceptive sensor
Isobath following using an altimeter as a unique exteroceptive sensor Luc Jaulin Lab-STICC, ENSTA Bretagne, Brest, France lucjaulin@gmail.com Abstract. We consider an underwater robot equipped with an
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 informationTopic # Feedback Control
Topic #7 16.31 Feedback Control State-Space Systems What are state-space models? Why should we use them? How are they related to the transfer functions used in classical control design and how do we develop
More informationPH 1A - Fall 2017 Solution 2
PH A - Fall 207 Solution 2 6.24 Let the tension of the cable between be T AB. First, we analyze the total torque τ on the boom with respect to fulcrum C. It should be 0 since the boon is static. τ = T
More informationA Sensor Driven Trade Study for Autonomous Navigation Capabilities
A Sensor Driven Trade Study for Autonomous Navigation Capabilities Sebastián Muñoz and E. Glenn Lightsey The University of Texas at Austin, Austin, TX, 78712 Traditionally, most interplanetary exploration
More informationCombining model-free and model-based Angle of Attack estimation for small fixed-wing UAVs using a standard sensor suite
Combining model-free and model-based Angle of Attack estimation for small fixed-wing UAVs using a standard sensor suite Andreas Wenz, Tor Arne Johansen, Andrea Cristofaro Centre for Autonomous Marine Operations
More informationAttitude determination method using single-antenna GPS, Gyro and Magnetometer
212 Asia-Pacific International Symposium on Aerospace echnology Nov. 13-1, Jeju, Korea Attitude determination method using single-antenna GPS, Gyro and Magnetometer eekwon No 1, Am Cho 2, Youngmin an 3,
More informationChapter 4 The Equations of Motion
Chapter 4 The Equations of Motion Flight Mechanics and Control AEM 4303 Bérénice Mettler University of Minnesota Feb. 20-27, 2013 (v. 2/26/13) Bérénice Mettler (University of Minnesota) Chapter 4 The Equations
More informationState Estimation and Motion Tracking for Spatially Diverse VLC Networks
State Estimation and Motion Tracking for Spatially Diverse VLC Networks GLOBECOM Optical Wireless Communications Workshop December 3, 2012 Anaheim, CA Michael Rahaim mrahaim@bu.edu Gregary Prince gbprince@bu.edu
More informationSpacecraft Angular Rate Estimation Algorithms For Star Tracker-Based Attitude Determination
AAS 3-191 Spacecraft Angular Rate Estimation Algorithms For Star Tracker-Based Attitude Determination Puneet Singla John L. Crassidis and John L. Junkins Texas A&M University, College Station, TX 77843
More informationAdaptive Two-Stage EKF for INS-GPS Loosely Coupled System with Unknown Fault Bias
Journal of Gloal Positioning Systems (26 Vol. 5 No. -2:62-69 Adaptive wo-stage EKF for INS-GPS Loosely Coupled System with Unnown Fault Bias Kwang Hoon Kim Jang Gyu Lee School of Electrical Engineering
More informationAutomated 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 informationReal-time Attitude Estimation Techniques Applied to a Four Rotor Helicopter
43rd IEEE Conference on Decision and Control December 14-17, 4 Atlantis, Paradise Island, Bahamas ThC9.3 Real- Attitude Estimation Techniques Applied to a Four Rotor Helicopter Matthew G. Earl and Raffaello
More informationA Study of the Effects of Stochastic Inertial Sensor Errors. in Dead-Reckoning Navigation
A Study of the Effects of Stochastic Inertial Sensor Errors in Dead-Reckoning Navigation Except where reference is made to the work of others, the work described in this thesis is my own or was done in
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 informationSystem identification and sensor fusion in dynamical systems. Thomas Schön Division of Systems and Control, Uppsala University, Sweden.
System identification and sensor fusion in dynamical systems Thomas Schön Division of Systems and Control, Uppsala University, Sweden. The system identification and sensor fusion problem Inertial sensors
More informationThe Kalman Filter ImPr Talk
The Kalman Filter ImPr Talk Ged Ridgway Centre for Medical Image Computing November, 2006 Outline What is the Kalman Filter? State Space Models Kalman Filter Overview Bayesian Updating of Estimates Kalman
More informationFILTERING SOLUTION TO RELATIVE ATTITUDE DETERMINATION PROBLEM USING MULTIPLE CONSTRAINTS
(Preprint) AAS FILTERING SOLUTION TO RELATIVE ATTITUDE DETERMINATION PROBLEM USING MULTIPLE CONSTRAINTS Richard Linares, John L. Crassidis and Yang Cheng INTRODUCTION In this paper a filtering solution
More informationSelf-Driving Car ND - Sensor Fusion - Extended Kalman Filters
Self-Driving Car ND - Sensor Fusion - Extended Kalman Filters Udacity and Mercedes February 7, 07 Introduction Lesson Ma 3 Estimation Problem Refresh 4 Measurement Udate Quiz 5 Kalman Filter Equations
More informationEvery real system has uncertainties, which include system parametric uncertainties, unmodeled dynamics
Sensitivity Analysis of Disturbance Accommodating Control with Kalman Filter Estimation Jemin George and John L. Crassidis University at Buffalo, State University of New York, Amherst, NY, 14-44 The design
More informationDesign and Flight Performance of the Orion. Pre-Launch Navigation System
Design and Flight Performance of the Orion Pre-Launch Navigation System Renato Zanetti 1, Greg Holt 2, Robert Gay 3, and Christopher D Souza 4 NASA Johnson Space Center, Houston, Texas 77058. Jastesh Sud
More informationTuning of Extended Kalman Filter for nonlinear State Estimation
OSR Journal of Computer Engineering (OSR-JCE) e-ssn: 78-0661,p-SSN: 78-877, Volume 18, ssue 5, Ver. V (Sep. - Oct. 016), PP 14-19 www.iosrjournals.org Tuning of Extended Kalman Filter for nonlinear State
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 informationOptimization-Based Control
Optimization-Based Control Richard M. Murray Control and Dynamical Systems California Institute of Technology DRAFT v1.7b, 23 February 2008 c California Institute of Technology All rights reserved. This
More informationNonlinear Observer Design for Dynamic Positioning
Author s Name, Company Title of the Paper DYNAMIC POSITIONING CONFERENCE November 15-16, 2005 Control Systems I J.G. Snijders, J.W. van der Woude Delft University of Technology (The Netherlands) J. Westhuis
More informationAided Inertial Navigation With Geometric Features: Observability Analysis
Aided Inertial Navigation With Geometric Features: Observability Analysis Yulin Yang - yuyang@udel.edu Guoquan Huang - ghuang@udel.edu Department of Mechanical Engineering University of Delaware, Delaware,
More informationE190Q Lecture 11 Autonomous Robot Navigation
E190Q Lecture 11 Autonomous Robot Navigation Instructor: Chris Clark Semester: Spring 013 1 Figures courtesy of Siegwart & Nourbakhsh Control Structures Planning Based Control Prior Knowledge Operator
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 informationc 2009 John Gregory Warner
c 2009 John Gregory Warner ATTITUDE DETERMINATION AND CONTROL OF NANO-SATELLITES BY JOHN GREGORY WARNER THESIS Submitted in partial fulfillment of the requirements for the degree of Master of Science in
More informationHigh-Gain Observers in Nonlinear Feedback Control. Lecture # 2 Separation Principle
High-Gain Observers in Nonlinear Feedback Control Lecture # 2 Separation Principle High-Gain ObserversinNonlinear Feedback ControlLecture # 2Separation Principle p. 1/4 The Class of Systems ẋ = Ax + Bφ(x,
More informationKalman Filter Computer Vision (Kris Kitani) Carnegie Mellon University
Kalman Filter 16-385 Computer Vision (Kris Kitani) Carnegie Mellon University Examples up to now have been discrete (binary) random variables Kalman filtering can be seen as a special case of a temporal
More informationLie Groups for 2D and 3D Transformations
Lie Groups for 2D and 3D Transformations Ethan Eade Updated May 20, 2017 * 1 Introduction This document derives useful formulae for working with the Lie groups that represent transformations in 2D and
More informationMeasurement Observers for Pose Estimation on SE(3)
Measurement Observers for Pose Estimation on SE(3) By Geoffrey Stacey u4308250 Supervised by Prof. Robert Mahony 24 September 2010 A thesis submitted in part fulfilment of the degree of Bachelor of Engineering
More informationA Crash Course on Kalman Filtering
A Crash Course on Kalman Filtering Dan Simon Cleveland State University Fall 2014 1 / 64 Outline Linear Systems Probability State Means and Covariances Least Squares Estimation The Kalman Filter Unknown
More informationEKF and SLAM. McGill COMP 765 Sept 18 th, 2017
EKF and SLAM McGill COMP 765 Sept 18 th, 2017 Outline News and information Instructions for paper presentations Continue on Kalman filter: EKF and extension to mapping Example of a real mapping system:
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 informationThe SPHERES Navigation System: from Early Development to On-Orbit Testing
The SPHERES Navigation System: from Early Development to On-Orbit Testing Simon Nolet MIT Space Systems Laboratory, Cambridge, MA 2139 The MIT Space Systems Laboratory has developed the Synchronized Position
More informationU RV QW (sin Θ)g + a x. where. D ΦΘ is determined assuming σ Φ, σ Θ are small and using the approximations
1 American Control Conference (ACC) Boston Marriott Copley Place July -8, 1. Boston, MA, USA Aircraft Sensor Fault Detection Using State and Input Estimation Ahmad Ansari and Dennis S. Bernstein Abstract
More informationPose Tracking II! Gordon Wetzstein! Stanford University! EE 267 Virtual Reality! Lecture 12! stanford.edu/class/ee267/!
Pose Tracking II! Gordon Wetzstein! Stanford University! EE 267 Virtual Reality! Lecture 12! stanford.edu/class/ee267/!! WARNING! this class will be dense! will learn how to use nonlinear optimization
More information