EE 570: Location and Navigation
|
|
- Bertha Jacobs
- 5 years ago
- Views:
Transcription
1 EE 570: Location and Navigation Error Mechanization (Tangential) Aly El-Osery 1 Stephen Bruder 2 1 Electrical Engineering Department, New Mexico Tech Socorro, New Mexico, USA 2 Electrical and Computer Engineering Department, Emry-Riddle Aeronautical Univesity Prescott, Arizona, USA March 23, 2014 Attitude Velocity Gravity Position Summary Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
2 Tangential Attitude Error Ċ t = Ct Ω t = Ct (Ω i Ω ie) = d dt [ (I +[δ ψ t t ] ])Ĉt = (I +[δ ψ t t ])Ĉt Ω t = [δ ψ t t ]Ĉt +(I +[δ ψ t t ]) Ĉ t = (I +[δ ψ t t ])Ĉt ( ˆΩ t + δω i δω ie) (I +[δ ψ t t ])Ĉt ˆΩ t + Ĉt (δω i δω ie) Attitude Velocity Gravity Position Summary Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
3 Tangential Attitude Error Ċ t = Ct Ω t = Ct (Ω i Ω ie) = d dt [ (I +[δ ψ t t ] ])Ĉt = (I +[δ ψ t t ])Ĉt Ω t = [δ ψ t t ]Ĉt +(I +[δ ψ t t ]) Ĉ t = (I +[δ ψ t t ])Ĉt ( ˆΩ t + δω i δω ie) (I +[δ ψ t t ])Ĉt ˆΩ t + Ĉt (δω i δω ie) Attitude Velocity Gravity Position Summary Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
4 Tangential Attitude Error Ċ t = Ct Ω t = Ct (Ω i Ω ie) = d dt [ (I +[δ ψ t t ] ])Ĉt = (I +[δ ψ t t ])Ĉt Ω t = [δ ψ t t ]Ĉt +(I +[δ ψ t t ]) Ĉ t = (I +[δ ψ t t ])Ĉt ( ˆΩ t + δω i δω ie) (I +[δ ψ t t ])Ĉt ˆΩ t + Ĉt (δω i δω ie) Attitude Velocity Gravity Position Summary Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
5 Tangential Attitude Error Ċ t = Ct Ω t = Ct (Ω i Ω ie) = d dt [ (I +[δ ψ t t ] ])Ĉt = (I +[δ ψ t t ])Ĉt Ω t = [δ ψ t t ]Ĉt +(I +[δ ψ t t ]) Ĉ t = (I +[δ ψ t t ])Ĉt ( ˆΩ t + δω i δω ie) (I +[δ ψ t t ])Ĉt ˆΩ t + Ĉt (δω i δω ie) [δ ψ t t ]δω t 0 Attitude Velocity Gravity Position Summary Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
6 Tangential Attitude Error Ċ t = Ct Ω t = Ct (Ω i Ω ie) = d dt [ (I +[δ ψ t t ] ])Ĉt = (I +[δ ψ t t ])Ĉt Ω t = [δ ψ t t ]Ĉt +(I +[δ ψ t t ]) Ĉ t = (I +[δ ψ t t ])Ĉt ( ˆΩ t + δω i δω ie) (I +[δ ψ t t ])Ĉt ˆΩ t + Ĉt (δω i δω ie) = Attitude Velocity Gravity Position Summary Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
7 Tangential Attitude Error Ċ t = Ct Ω t = Ct (Ω i Ω ie) = d dt [ (I +[δ ψ t t ] ])Ĉt = (I +[δ ψ t t ])Ĉt Ω t = [δ ψ t t ]Ĉt +(I +[δ ψ t t ]) Ĉ t = (I +[δ ψ t t ])Ĉt ( ˆΩ t + δω i δω ie) (I +[δ ψ t t ])Ĉt ˆΩ t + Ĉt (δω i δω ie) [δ ψ t t ] = Ĉt (δω i δω ie)ĉt = [Ĉ t(δ ω i δ ω ie) ] (1) Attitude Velocity Gravity Position Summary Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
8 Tangential Attitude Error Ċ t = Ct Ω t = Ct (Ω i Ω ie) = d dt [ (I +[δ ψ t t ] ])Ĉt = (I +[δ ψ t t ])Ĉt Ω t = [δ ψ t t ]Ĉt +(I +[δ ψ t t ]) Ĉ t = (I +[δ ψ t t ])Ĉt ( ˆΩ t + δω i δω ie) (I +[δ ψ t t ])Ĉt ˆΩ t + Ĉt (δω i δω ie) [δ ψ t t ] = Ĉt (δω i δω ie)ĉt = [Ĉ t(δ ω i δ ω ie) ] (1) δ ψ t t = Ĉt (δω i δ ω ie) (2) Attitude Velocity Gravity Position Summary Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
9 Tangential Attitude Error (cont.) δ ψ t t = Ĉt (δω i δ ω ie) = Ĉ t δω i Ĉt ( ω ie ˆ ω ie) = Ĉ t δω i (Ĉt C t I) ω t ie = Ĉ t δω i + δ ψ t t ωt ie Attitude Velocity Gravity Position Summary Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
10 Tangential Attitude Error (cont.) δ ψ t t = Ĉt (δω i δ ω ie) = Ĉ t δω i Ĉt ( ω ie ˆ ω ie) = Ĉ t δω i (Ĉt C t I) ω t ie = Ĉ t δω i + δ ψ t t ωt ie δ ψ t t = Ĉt δω i Ω t ie δ ψ t t (3) Attitude Velocity Gravity Position Summary Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
11 Velocity v t t = Ct f i + gt 2Ωt ie v t t (4) ˆ v t t = Ĉt ˆ f i + ˆ g t 2Ωt ˆ v ie t t = (I [δ ψ t t ])Ct ( f i δ f i )+ ˆ g t 2Ωt ˆ v ie t t (5) Attitude Velocity Gravity Position Summary Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
12 Velocity v t t = Ct f i + gt 2Ωt ie v t t (4) ˆ v t t = Ĉt ˆ f i + ˆ g t 2Ωt ˆ v ie t t = (I [δ ψ t t ])Ct ( f i δ f i )+ ˆ g t 2Ωt ˆ v ie t t δ v t t = v t t ˆ v t t = [δ ψ t t ]Ct fi + Ĉt δ f i + δ g t 2Ωt ieδ v t t = [δ ψ t t ]Ĉt ˆ f i + Ĉt δ f i + δ g t 2Ωt ieδ v t t (5) Attitude Velocity Gravity Position Summary Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
13 Velocity v t t = Ct f i + gt 2Ωt ie v t t (4) ˆ v t t = Ĉt ˆ f i + ˆ g t 2Ωt ˆ v ie t t = (I [δ ψ t t ])Ct ( f i δ f i )+ ˆ g t 2Ωt ˆ v ie t t δ v t t = v t t ˆ v t t = [δ ψ t t ]Ct fi + Ĉt δ f i + δ g t 2Ωt ieδ v t t = [δ ψ t t ]Ĉt ˆ f i + Ĉt δ f i + δ g t 2Ωt ieδ v t t (5) δ v t t = [ C ˆ t ˆ f i ]δ ψ t t + Ĉt δ f i + δ g t 2Ωt ˆ v ie t t (6) Attitude Velocity Gravity Position Summary Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
14 Gravity Error δ g t Ĉt e 2g 0 (ˆL ) ˆ r e e res e (ˆL ) ˆ r 2(ˆ r e e e e )T Ĉt e δ r t t (7) Attitude Velocity Gravity Position Summary Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
15 Position r t t = vt t (8) Attitude Velocity Gravity Position Summary Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
16 Position r t t = vt t (8) δ r t t = δ vt t (9) Attitude Velocity Gravity Position Summary Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
17 Summary - in terms of δ f i, δ ω i δ ψ t t Ω t ie δ v t t = [ ˆ C t ˆ f i ] 2Ωt ie Ĉe t 2g 0 ( ˆL ) ˆ r res e e e ( ˆL ) 2(ˆ r e ˆ r e e e )T Ĉt e δ r t t I Ĉ t Ĉ t 0 δ f i δ ω i 0 0 δ ψ t t δ v t t δ r t t + (10) Attitude Velocity Gravity Position Summary Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
18 Notation Used Truth value Measured value x x Estimated or computed value ˆ x Error δ x = x ˆ x Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
19 Notation Used Truth value Measured value Estimated or computed value x x Nothing aove ˆ x Error δ x = x ˆ x Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
20 Notation Used Truth value Measured value Estimated or computed value x x ˆ x Use tilde Error δ x = x ˆ x Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
21 Notation Used Truth value Measured value Estimated or computed value x x ˆ x Error δ x = x ˆ x Use hat Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
22 Notation Used Truth value Measured value x x Estimated or computed value ˆ x Error δ x δ x = x ˆ x ˆ x Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
23 Linearization using Taylor Series Expansion Given a non-linear system x = f( x, t) Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
24 Linearization using Taylor Series Expansion Given a non-linear system x = f( x, t) Let s assume we have an estimate of x, i.e., ˆ x such that x = ˆ x + δ x x = ˆ x + δ x = f(ˆ x + δ x, t) (11) Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
25 Linearization using Taylor Series Expansion Given a non-linear system x = f( x, t) Let s assume we have an estimate of x, i.e., ˆ x such that x = ˆ x + δ x Using Taylor series expansion x = ˆ x + δ x = f(ˆ x + δ x, t) (11) f(ˆ x + δ x, t) = ˆ x + δ x = f(ˆ x, t) + ˆ x + f( x, t) x f( x, t) x δ x x= ˆ x x= ˆ x δ x + H.O.T Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
26 Linearization using Taylor Series Expansion Given a non-linear system x = f( x, t) Let s assume we have an estimate of x, i.e., ˆ x such that x = ˆ x + δ x x = ˆ x + δ x = f(ˆ x + δ x, t) (11) Using Taylor series expansion f(ˆ x + δ x, t) = ˆ x + δ x f( x, t) = f(ˆ x, t) + x = ˆ x f( x, t) + δ x x x= ˆ x x= ˆ x δ x + H.O.T Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
27 Linearization using Taylor Series Expansion Given a non-linear system x = f( x, t) Let s assume we have an estimate of x, i.e., ˆ x such that x = ˆ x + δ x Using Taylor series expansion x = ˆ x + δ x = f(ˆ x + δ x, t) (11) f(ˆ x + δ x, t) = ˆ x + δ x = f(ˆ x, t) + ˆ x + δ x f( x, t) x f( x, t) x f( x, t) x δ x x= ˆ x x= ˆ x x= ˆ x δ x + H.O.T δ x (12) Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
28 Actual Measurements Initially the accelerometer and gyroscope measurements, f i and ω i, respectively, will e modeled as f i = f i + f i (13) ω i = ω i + ω i (14) where f i and ω i are the specific force and angular rates, respectively; and f i and ω i represents the errors. In later lectures we will discuss more detailed description of these errors. Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
29 Actual Measurements Initially the accelerometer and gyroscope measurements, f i and ω i, respectively, will e modeled as f i = f i + f i (13) these terms may ω i = ω i + ω e expanded further i (14) where f i and ω i are the specific force and angular rates, respectively; and f i and ω i represents the errors. In later lectures we will discuss more detailed description of these errors. Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
30 Error Modeling Example Accelerometers f i = a +(I + M a ) f i + nla + w a (15) Gyroscopes ω i = g +(I + M g ) ω i + G g f i + w g (16) Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
31 Error Modeling Example Accelerometers f i = a +(I + M a ) f i + nla + w a (15) Biases Gyroscopes ω i = g +(I + M g ) ω i + G g f i + w g (16) Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
32 Error Modeling Example Accelerometers f i = a +(I + M a ) f i + nla + w a (15) Misalignment and SF Errors Gyroscopes ω i = g +(I + M g ) ω i + G g f i + w g (16) Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
33 Error Modeling Example Accelerometers f i = a +(I + M a ) f i + nla + w a (15) Non-linearity Gyroscopes ω i = g +(I + M g ) ω i + G g f i + w g (16) Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
34 Error Modeling Example Accelerometers f i = a +(I + M a ) f i + nla + w a (15) Gyroscopes G-Sensitivity ω i = g +(I + M g ) ω i + G g f i + w g (16) Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
35 Error Modeling Example Accelerometers f i = a +(I + M a ) f i + nla + w a (15) Noise Gyroscopes ω i = g +(I + M g ) ω i + G g f i + w g (16) Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
36 Pos, Vel, Force and Angular Rate Errors Position error Velocity error δ r γ β = r γ β ˆ r γ β (17) δ v γ β = v γ β ˆ v γ β (18) Specific force errors Angular rate errors δ f i = f i ˆ f i (19) e f i = f i ˆ f i = δ f i (20) δ ω i = ω i ˆ ω i (21) e ω i = ω i ˆ ω i = δ ω i (22) Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
37 Attitude Error Definition Define δc γ = C γ Ĉ γ = e [δ ψ γ γ ] I +[δ ψ γ γ ] (23) This is the error in attitude resulting from errors in estimating the angular rates. Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
38 Attitude Error Properties The attitude error is a multiplicative small angle transformation from the actual frame to the computed frame Ĉ γ = (I [δ ψ γ γ ])C γ (24) Similarly, C γ = (I +[δ ψ γ γ ])Ĉ γ (25) Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
39 Specific Force and Agnular Rates Similarly we can attempt to estimate the specific force and angular rate y applying correction ased on our estimate of the error. ˆ f i = f i ˆ f i (26) ˆ ω i = ω i ˆ ω i (27) where ˆ f i and ˆ ω i are the accelerometer and gyroscope estimated caliration values, respectively. Aly El-Osery, Stephen Bruder (NMT,ERAU) EE 570: Location and Navigation March 23, / 15
EE 570: Location and Navigation
EE 570: Location and Navigation Error Mechanization (ECEF) Aly El-Osery 1 Stephen Bruder 2 1 Electrical Engineering Department, New Mexico Tech Socorro, New Mexico, USA 2 Electrical and Computer Engineering
More informationEE 570: Location and Navigation
EE 570: Location and Navigation Error Mechanization (ECEF) Aly El-Osery Electrical Engineering Department, New Mexico Tech Socorro, New Mexico, USA April 11, 2013 Attitude Velocity Gravity Position Summary
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 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 informationEE 570: Location and Navigation
EE 570: Locatio ad Navigatio Error Mechaizatio (NAV) Aly El-Osery Kevi Wedeward Electrical Egieerig Departmet, New Mexico Tech Socorro, New Mexico, USA I Collaboratio with Stephe Bruder Electrical ad Computer
More informationEE 565: Position, Navigation and Timing
EE 565: Position, Navigation and Timing Navigation Mathematics: Angular and Linear Velocity Kevin Wedeward Aly El-Osery Electrical Engineering Department New Mexico Tech Socorro, New Mexico, USA In Collaboration
More informationEE 570: Location and Navigation
EE 570: Location and Navigation Sensor Technology Stephen Bruder 1 Aly El-Osery 2 1 Electrical and Computer Engineering Department, Embry-Riddle Aeronautical Univesity Prescott, Arizona, USA 2 Electrical
More informationEE 565: Position, Navigation, and Timing
EE 565: Position, Navigation, and Timing Kalman Filtering Example Aly El-Osery Kevin Wedeward Electrical Engineering Department, New Mexico Tech Socorro, New Mexico, USA In Collaboration with Stephen Bruder
More informationEE 570: Location and Navigation
EE 570: Location and Navigation Navigation Mathematics: Coordinate Frames Kevin Wedeward Aly El-Osery Electrical Engineering Department, New Mexico Tech Socorro, New Mexico, USA In Collaboration with Stephen
More informationEE 570: Location and Navigation
EE 570: Location and Navigation Navigation Equations: Nav Mechanization Aly El-Osery Kevin Wedeward Electrical Engineering Department, New Mexico Tech Socorro, New Mexico, USA In Collaboration with Stephen
More informationEE 570: Location and Navigation
EE 570: Location and Navigation Sensor Technology Aly El-Osery Kevin Wedeward Electrical Engineering Department, New Mexico Tech Socorro, New Mexico, USA In Collaboration with Stephen Bruder Electrical
More informationEE 570: Location and Navigation
EE 570: Location and Navigation Sensor Technology Stephen Bruder 1 Aly El-Osery 2 1 Electrical and Computer Engineering Department, Embry-Riddle Aeronautical Univesity Prescott, Arizona, USA 2 Electrical
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 informationNavigation Mathematics: Kinematics (Coordinate Frame Transformation) EE 570: Location and Navigation
Lecture Navigation Mathematics: Kinematics (Coordinate Frame Transformation) EE 570: Location and Navigation Lecture Notes Update on Feruary 16, 2016 Aly El-Osery and Kevin Wedeward, Electrical Engineering
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 informationNMT EE 589 & UNM ME 482/582 ROBOT ENGINEERING. Dr. Stephen Bruder NMT EE 589 & UNM ME 482/582
NM EE 589 & UNM ME 48/58 ROBO ENGINEERING Dr. Stephen Bruder NM EE 589 & UNM ME 48/58 5. Robot Dynamics 5. he Microbot Arm Dynamic Model A Second Dynamic Model Example: he Microbot Robot Arm he Denavit-Hartenberg
More informationEE 521: Instrumentation and Measurements
Aly El-Osery Electrical Engineering Department, New Mexico Tech Socorro, New Mexico, USA August 30, 2009 1 / 19 1 Course Overview 2 Measurement System 3 Noise 4 Error Analysis 2 / 19 Textbooks C.W. De
More informationNavigation Mathematics: Angular and Linear Velocity EE 570: Location and Navigation
Lecture Nvigtion Mthemtics: Angulr n Liner Velocity EE 57: Loction n Nvigtion Lecture Notes Upte on Februry, 26 Kevin Weewr n Aly El-Osery, Electricl Engineering Dept., New Mexico Tech In collbortion with
More informationContinuous Preintegration Theory for Graph-based Visual-Inertial Navigation
Continuous Preintegration Theory for Graph-based Visual-Inertial Navigation Kevin Ecenhoff - ec@udel.edu Patric Geneva - pgeneva@udel.edu Guoquan Huang - ghuang@udel.edu Department of Mechanical Engineering
More informationIntroduction & Laplace Transforms Lectures 1 & 2
Introduction & Lectures 1 & 2, Professor Department of Electrical and Computer Engineering Colorado State University Fall 2016 Control System Definition of a Control System Group of components that collectively
More informationERRATA TO: STRAPDOWN ANALYTICS, FIRST EDITION PAUL G. SAVAGE PARTS 1 AND 2 STRAPDOWN ASSOCIATES, INC. 2000
ERRATA TO: STRAPDOWN ANAYTICS, FIRST EDITION PAU G. SAVAGE PARTS 1 AND 2 STRAPDOWN ASSOCIATES, INC. 2000 Corrections Reported From April 2004 To Present Pg. 10-45 - Section 10.1.3, first paragraph - Frequency
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 informationEE C245 ME C218 Introduction to MEMS Design
EE C245 ME C218 Introduction to MEMS Design Fall 2007 Prof. Clark T.-C. Nguyen Dept. of Electrical Engineering & Computer Sciences University of California at Berkeley Berkeley, CA 94720 Lecture 20: Equivalent
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 informationNMT EE 589 & UNM ME 482/582 ROBOT ENGINEERING. Dr. Stephen Bruder NMT EE 589 & UNM ME 482/582
NMT EE 589 & UNM ME 482/582 ROBOT ENGINEERING Dr. Stephen Bruder NMT EE 589 & UNM ME 482/582 4. Moton Knematcs 4.2 Angular Velocty Knematcs Summary From the last lecture we concluded that: If the jonts
More informationEE 521: Instrumentation and Measurements
Aly El-Osery Electrical Engineering Department, New Mexico Tech Socorro, New Mexico, USA September 23, 2009 1 / 18 1 Sampling 2 Quantization 3 Digital-to-Analog Converter 4 Analog-to-Digital Converter
More informationAN IMPROVED STATIONARY FINE SELF-ALIGNMENT APPROACH FOR SINS USING MEASUREMENT AUGMENTATION
Belo Horizonte, MG, 2 a 24 de Setembro de 214 AN IMPROVED STATIONARY FINE SELF-ALIGNMENT APPROACH FOR SINS USING MEASUREMENT AUGMENTATION Felipe Oliveira e Silva, Elder Moreira Hemerly, Waldemar de Castro
More informationWave Phenomena Physics 15c. Lecture 11 Dispersion
Wave Phenomena Physics 15c Lecture 11 Dispersion What We Did Last Time Defined Fourier transform f (t) = F(ω)e iωt dω F(ω) = 1 2π f(t) and F(w) represent a function in time and frequency domains Analyzed
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 informationLecture 7 Discrete Systems
Lecture 7 Discrete Systems EE 52: Instrumentation and Measurements Lecture Notes Update on November, 29 Aly El-Osery, Electrical Engineering Dept., New Mexico Tech 7. Contents The z-transform 2 Linear
More informationEE 521: Instrumentation and Measurements
Aly El-Osery Electrical Engineering Department, New Mexico Tech Socorro, New Mexico, USA November 1, 2009 1 / 27 1 The z-transform 2 Linear Time-Invariant System 3 Filter Design IIR Filters FIR Filters
More informationLecture Schedule: Week Date Lecture Title
http://elec34.org Sampling and CONVOLUTION 24 School of Information Technology and Electrical Engineering at The University of Queensland Lecture Schedule: Week Date Lecture Title 2-Mar Introduction 3-Mar
More informationMath 266 Midterm Exam 2
Math 266 Midterm Exam 2 March 2st 26 Name: Ground Rules. Calculator is NOT allowed. 2. Show your work for every problem unless otherwise stated (partial credits are available). 3. You may use one 4-by-6
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 informationThe Calculus of Vec- tors
Physics 2460 Electricity and Magnetism I, Fall 2007, Lecture 3 1 The Calculus of Vec- Summary: tors 1. Calculus of Vectors: Limits and Derivatives 2. Parametric representation of Curves r(t) = [x(t), y(t),
More informationPendulum with a vibrating base
Pendulum with a vibrating base Chennai, March 14, 2006 1 1 Fast perturbations- rapidly oscillating perturbations We consider perturbations of a Hamiltonian perturbed by rapid oscillations. Later we apply
More informationClosed-form solution of visual-inertial structure from motion
Closed-form solution of visual-inertial structure from motion Agostino Martinelli To cite this version: Agostino Martinelli. Closed-form solution of visual-inertial structure from motion. International
More informationMathematical models for class-d amplifiers
Mathematical models for class-d amplifiers Stephen Cox School of Mathematical Sciences, University of Nottingham, UK 12 November 2012 Stephen Cox Mathematical models for class-d amplifiers 1/38 Background
More informationTheoretical Physics II B Quantum Mechanics. Lecture 5
Theoretical Physics II B Quantum Mechanics Lecture 5 Frank Krauss February 11, 2014 Solutions to previous control questions 4.1 Determine first the transformation matrix ˆT through ( ˆT ij = ψ i φ j =
More informationNeed for new theory Stern-Gerlach Experiments. Analogy with mathematics of light Feynman s double slit thought experiment
Lecture 3 Need for new theory Stern-Gerlach Experiments Some doubts Analogy with mathematics of light Feynman s double slit thought experiment Three empirical laws To explain one physical phenomenon Radiation
More informationLecture IV: Time Discretization
Lecture IV: Time Discretization Motivation Kinematics: continuous motion in continuous time Computer simulation: Discrete time steps t Discrete Space (mesh particles) Updating Position Force induces acceleration.
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 informationRefinements to the General Methodology Behind Strapdown Airborne Gravimetry
Refinements to the General Methodology Behind Strapdown Airborne Gravimetry AE 8900 MS Special Problems Report Space Systems Design Lab (SSDL) Guggenheim School of Aerospace Engineering Georgia Institute
More informationWe call these coupled differential equations, because as you can see the derivative of θ1 involves θ2 and the derivative of the θ2 involves θ1.
Torsion Waves In Monday s class we saw waves on the above device. By the end of today s class, I hope you will have some idea of how those waves arise from the equations you already know. The above device
More informationMotion in Space. MATH 311, Calculus III. J. Robert Buchanan. Fall Department of Mathematics. J. Robert Buchanan Motion in Space
Motion in Space MATH 311, Calculus III J. Robert Buchanan Department of Mathematics Fall 2011 Background Suppose the position vector of a moving object is given by r(t) = f (t), g(t), h(t), Background
More informationP321(b), Assignement 1
P31(b), Assignement 1 1 Exercise 3.1 (Fetter and Walecka) a) The problem is that of a point mass rotating along a circle of radius a, rotating with a constant angular velocity Ω. Generally, 3 coordinates
More informationPhysics 2514 Lecture 22
Physics 2514 Lecture 22 P. Gutierrez Department of Physics & Astronomy University of Oklahoma Physics 2514 p. 1/15 Information Information needed for the exam Exam will be in the same format as the practice
More informationLecture V: The game-engine loop & Time Integration
Lecture V: The game-engine loop & Time Integration The Basic Game-Engine Loop Previous state: " #, %(#) ( #, )(#) Forces -(#) Integrate velocities and positions Resolve Interpenetrations Per-body change
More informationQuestion 1: Spherical Pendulum
Question 1: Spherical Pendulum Consider a two-dimensional pendulum of length l with mass M at its end. It is easiest to use spherical coordinates centered at the pivot since the magnitude of the position
More informationPhysics Fall Mechanics, Thermodynamics, Waves, Fluids. Lecture 20: Rotational Motion. Slide 20-1
Physics 1501 Fall 2008 Mechanics, Thermodynamics, Waves, Fluids Lecture 20: Rotational Motion Slide 20-1 Recap: center of mass, linear momentum A composite system behaves as though its mass is concentrated
More informationDealing with Rotating Coordinate Systems Physics 321. (Eq.1)
Dealing with Rotating Coordinate Systems Physics 321 The treatment of rotating coordinate frames can be very confusing because there are two different sets of aes, and one set of aes is not constant in
More informationConservation of Angular Momentum
10 March 2017 Conservation of ngular Momentum Lecture 23 In the last class, we discussed about the conservation of angular momentum principle. Using RTT, the angular momentum principle was given as DHo
More informationUAVBook Supplement Full State Direct and Indirect EKF
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
More informationLecture 38: Equations of Rigid-Body Motion
Lecture 38: Equations of Rigid-Body Motion It s going to be easiest to find the equations of motion for the object in the body frame i.e., the frame where the axes are principal axes In general, we can
More informationLecture 3: Kinematics of Particles
Lecture 3: Kinematics of Particles Recap: Numerical integration [Appendix A, 2.1] 2 and 3-Dimensional Motion Resolving vectors in 3-Dimensions Differentiation of vectors leading into a review of Polar
More informationRecap. The bigger the exhaust speed, ve, the higher the gain in velocity of the rocket.
Recap Classical rocket propulsion works because of momentum conservation. Exhaust gas ejected from a rocket pushes the rocket forwards, i.e. accelerates it. The bigger the exhaust speed, ve, the higher
More informationVerification of a Dual-State Extended Kalman Filter with Lidar-Enabled Autonomous Hazard- Detection for Planetary Landers
Marquette University e-publications@marquette Master's Theses (29 -) Dissertations, Theses, and Professional Projects Verification of a Dual-State Extended Kalman Filter with Lidar-Enabled Autonomous Hazard-
More informationPhysics 101 Discussion Week 12 Explanation (2011)
Physics 101 Discussion Week 12 Eplanation (2011) D12-1 Horizontal oscillation Q0. This is obviously about a harmonic oscillator. Can you write down Newton s second law in the (horizontal) direction? Let
More informationAutonomous satellite orbit prediction
Autonomous satellite orbit prediction Mari Seppänen Department of Mathematics Tampere University of Technology P.O. Box 553, FI-33101 Tampere, Finland e-mail: mari.j.seppanen@tut.fi Abstract A method to
More informationLecture 38: Equations of Rigid-Body Motion
Lecture 38: Equations of Rigid-Body Motion It s going to be easiest to find the equations of motion for the object in the body frame i.e., the frame where the axes are principal axes In general, we can
More informationEE 230 Lecture 33. Nonlinear Circuits and Nonlinear Devices. Diode BJT MOSFET
EE 230 Lecture 33 Nonlinear Circuits and Nonlinear Devices Diode BJT MOSFET Review from Last Time: n-channel MOSFET Source Gate L Drain W L EFF Poly Gate oxide n-active p-sub depletion region (electrically
More informationSolutions to Homework 1, Introduction to Differential Equations, 3450: , Dr. Montero, Spring y(x) = ce 2x + e x
Solutions to Homewor 1, Introduction to Differential Equations, 3450:335-003, Dr. Montero, Spring 2009 problem 2. The problem says that the function yx = ce 2x + e x solves the ODE y + 2y = e x, and ass
More informationSpacecraft Attitude Dynamics for Undergraduates
Session 1123 Spacecraft Attitude Dynamics for Undergraduates Dr. Rachel Shinn Embry Riddle Aeronautical University, Prescott, AZ Abstract Teaching spacecraft attitude dynamics to undergraduate students
More informationInertial Navigation and Various Applications of Inertial Data. Yongcai Wang. 9 November 2016
Inertial Navigation and Various Applications of Inertial Data Yongcai Wang 9 November 2016 Types of Gyroscope Mechanical Gyroscope Laser Gyroscope Sagnac Effect Stable Platform IMU and Strapdown IMU In
More informationdt r r r V(x,t) = F(x,t)dx
Quantum Mechanics and Atomic Physics Lecture 3: Schroedinger s Equation: Part I http://www.physics.rutgers.edu/ugrad/361 Prof. Sean Oh Announcement First homework due on Wednesday Sept 14 at the beginning
More informationStabilization of Insect Flight via Sensors of Coriolis Force
1 Problem Stabilization of Insect Flight via Sensors of Coriolis Force Kirk T. McDonald Joseph Henry Laboratories, Princeton University, Princeton, NJ 0544 February 17, 2007) The dynamics of insect flight
More informationRotational Kinematics
Rotational Kinematics 1 Linear Motion Rotational Motion all variables considered positive if motion in counterclockwise direction displacement velocity acceleration angular displacement (Δθ) angular velocity
More informationRigid body simulation. Once we consider an object with spatial extent, particle system simulation is no longer sufficient
Rigid body dynamics Rigid body simulation Once we consider an object with spatial extent, particle system simulation is no longer sufficient Rigid body simulation Unconstrained system no contact Constrained
More informationErrata for Instructor s Solutions Manual for Gravity, An Introduction to Einstein s General Relativity 1st printing
Errata for Instructor s Solutions Manual for Gravity, An Introduction to Einstein s General Relativity st printing Updated 7/7/003 (hanks to Scott Fraser who provided almost all of these.) Statement of
More informationPeriodic functions: simple harmonic oscillator
Periodic functions: simple harmonic oscillator Recall the simple harmonic oscillator (e.g. mass-spring system) d 2 y dt 2 + ω2 0y = 0 Solution can be written in various ways: y(t) = Ae iω 0t y(t) = A cos
More informationPart: Frequency and Time Domain
Numerical Methods Fourier Transform Pair Part: Frequency and Time Domain For more details on this topic Go to Clic on eyword Clic on Fourier Transform Pair You are free to Share to copy, distribute, display
More informationCHEE 319 Tutorial 3 Solutions. 1. Using partial fraction expansions, find the causal function f whose Laplace transform. F (s) F (s) = C 1 s + C 2
CHEE 39 Tutorial 3 Solutions. Using partial fraction expansions, find the causal function f whose Laplace transform is given by: F (s) 0 f(t)e st dt (.) F (s) = s(s+) ; Solution: Note that the polynomial
More informationRobot Dynamics II: Trajectories & Motion
Robot Dynamics II: Trajectories & Motion Are We There Yet? METR 4202: Advanced Control & Robotics Dr Surya Singh Lecture # 5 August 23, 2013 metr4202@itee.uq.edu.au http://itee.uq.edu.au/~metr4202/ 2013
More informationFirst Year Physics: Prelims CP1. Classical Mechanics: Prof. Neville Harnew. Problem Set III : Projectiles, rocket motion and motion in E & B fields
HT017 First Year Physics: Prelims CP1 Classical Mechanics: Prof Neville Harnew Problem Set III : Projectiles, rocket motion and motion in E & B fields Questions 1-10 are standard examples Questions 11-1
More informationAComparativeStudy ofall-accelerometer Strapdownsfor UAVINS
AComparativeStudy ofall-accelerometer Strapdownsfor UAVINS Philippe Cardou(pcardou@cim.mcgill.ca) Jorge Angeles (angeles@cim.mcgill.ca) Department of Mechanical Engineering & Centrefor Intelligent Machines
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 informationGyroscopes and statics
Gyroscopes and statics Announcements: Welcome back from Spring Break! CAPA due Friday at 10pm We will finish Chapter 11 in H+R on angular momentum and start Chapter 12 on stability. Friday we will begin
More informationUsing Two Tri-Axis Accelerometers for Rotational Measurements
Using Two Tri-Axis Accelerometers for Rotational Measurements Introduction In man applications, customers would like to measure rotational motions (angular velocit, angular acceleration) in addition to
More informationApr 14, Calculus with Algebra and Trigonometry II Lecture 20More physics Aprapplications
Calculus with Algebra and Trigonometry II Lecture 20 More physics applications Apr 14, 2015 14, 2015 1 / 14 Motion in two dimensions A particle s motion can be described by specifying how the coordinates
More informationCS491/691: Introduction to Aerial Robotics
CS491/691: Introduction to Aerial Robotics Topic: Midterm Preparation Dr. Kostas Alexis (CSE) Areas of Focus Coordinate system transformations (CST) MAV Dynamics (MAVD) Navigation Sensors (NS) State Estimation
More informationEnergy during a burst of deceleration
Problem 1. Energy during a burst of deceleration A particle of charge e moves at constant velocity, βc, for t < 0. During the short time interval, 0 < t < t its velocity remains in the same direction but
More informationTaylor Series. Math114. March 1, Department of Mathematics, University of Kentucky. Math114 Lecture 18 1/ 13
Taylor Series Math114 Department of Mathematics, University of Kentucky March 1, 2017 Math114 Lecture 18 1/ 13 Given a function, can we find a power series representation? Math114 Lecture 18 2/ 13 Given
More informationCMB beyond a single power spectrum: Non-Gaussianity and frequency dependence. Antony Lewis
CMB beyond a single power spectrum: Non-Gaussianity and frequency dependence Antony Lewis http://cosmologist.info/ Evolution of the universe Opaque Transparent Hu & White, Sci. Am., 290 44 (2004) CMB temperature
More informationPhysics 201, Lecture 18
q q Physics 01, Lecture 18 Rotational Dynamics Torque Exercises and Applications Rolling Motion Today s Topics Review Angular Velocity And Angular Acceleration q Angular Velocity (ω) describes how fast
More informationWe are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors
We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists 3,500 108,000 1.7 M Open access books available International authors and editors Downloads Our
More informationProbing P & T-violation Beyond the Standard Model. Aaron E. Leanhardt
An Electron EDM Search in HfF + : Probing P & T-violation Beyond the Standard Model Aaron E. Leanhardt Experiment: Laura Sinclair, Russell Stutz & Eric Cornell Theory: Ed Meyer & John Bohn JILA, NIST,
More informationGeneralized mathematical model of a linear single-axis accelerometer as an integral part of the inclinometer
Generalized mathematical model of a linear single-axis accelerometer as an integral part of the inclinometer A. N. Krasnov 1, G. Y. Kolovertnov 2, V. E. Lyalin 3 1, 2 Ufa State Petroleum Technological
More informationEE C245 ME C218 Introduction to MEMS Design Fall 2007
EE C45 ME C8 Introduction to MEMS Design Fall 007 Prof. Clark T.C. Nguyen Dept. of Electrical Engineering & Computer Sciences University of California at Berkeley Berkeley, CA 9470 Lecture 3: Input Modeling
More informationNAVAL POSTGRADUATE SCHOOL THESIS
NAVAL POSTGRADUATE SCHOOL MONTEREY, CALIFORNIA THESIS ATTITUDE DETERMINATION FOR THE THREE-AXIS SPACECRAFT SIMULATOR (TASS) BY APPLICATION OF PARTICLE FILTERING TECHNIQUES by Ioannis Kassalias June 2005
More informationE11 Lecture 13: Motors. Professor Lape Fall 2010
E11 Lecture 13: Motors Professor Lape Fall 2010 Overview How do electric motors work? Electric motor types and general principles of operation How well does your motor perform? Torque and power output
More informationIntroductory Physics. Week 2015/05/29
2015/05/29 Part I Summary of week 6 Summary of week 6 We studied the motion of a projectile under uniform gravity, and constrained rectilinear motion, introducing the concept of constraint force. Then
More informationChapter 3: Second Order ODE 3.8 Elements of Particle Dy
Chapter 3: Second Order ODE 3.8 Elements of Particle Dynamics 3 March 2018 Objective The objective of this section is to explain that any second degree linear ODE represents the motion of a particle. This
More informationAS3010: Introduction to Space Technology
AS3010: Introduction to Space Technology L E C T U R E 22 Part B, Lecture 22 19 April, 2017 C O N T E N T S Attitude stabilization passive and active. Actuators for three axis or active stabilization.
More informationResearch Article Relative Status Determination for Spacecraft Relative Motion Based on Dual Quaternion
Mathematical Prolems in Engineering Volume 214, Article ID 62724, 7 pages http://dx.doi.org/1.1155/214/62724 Research Article Relative Status Determination for Spacecraft Relative Motion Based on Dual
More informationOptional Problems on the Harmonic Oscillator
8.04: Quantum Mechanics Professor Allan Adams Massachusetts Institute of Technology Tuesday March 9 Optional Problems on the Harmonic Oscillator. Coherent States Consider a state ϕ α which is an eigenstate
More informationCalibration of a magnetometer in combination with inertial sensors
Calibration of a magnetometer in combination with inertial sensors Manon Kok, Linköping University, Sweden Joint work with: Thomas Schön, Uppsala University, Sweden Jeroen Hol, Xsens Technologies, the
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 informationLecture D21 - Pendulums
J. Peraire 16.07 Dynamics Fall 2004 Version 1.1 ecture D21 - Pendulums A pendulum is a rigid body suspended from a fixed point (hinge) which is offset with respect to the body s center of mass. If all
More informationSearching for gravitational waves from neutron stars
Searching for gravitational waves from neutron stars Ian Jones D.I.Jones@soton.ac.uk General Relativity Group, Southampton University Ian Jones Searching for gravitational waves from neutron stars 1/23
More informationModern Navigation. Thomas Herring
12.215 Modern Navigation Thomas Herring Basic Statistics Summary of last class Statistical description and parameters Probability distributions Descriptions: expectations, variances, moments Covariances
More information