Smartphone sensor based orientation determination for indoor navigation

Size: px
Start display at page:

Download "Smartphone sensor based orientation determination for indoor navigation"

Transcription

1 Smartphone sensor based orientation determination for indoor naviation LBS Conference Andreas Ettliner Research Group Enineerin Geodesy Contact:

2 Outline Motivation and Methodoloical approach Orientation and Smartphone Measurement Data Kalman Filter for Orientation Determination Test Measurements and Results Conclusion and Outloo 2

3 Motivation and Methodoloical approach Indoor naviation Applications are aumented by 2D-PDR ( Pedestrian Dead Reconin ) - Distance d i and yaw ψ i à Position P i Critical part à calculation of users movin direction respectively yaw Methodoloical approach Calculate yaw with Kalman filter Detect turns of the user Smooth trajectory when user is walin straiht Reduce systematics of the smartphone sensors ψ 2 ψ 3 ψ 1 d 1 d 2 P 2 P 1 d 3 P 4 ψ 4 d4 P 0 P 3 3

4 Orientation and Smartphone Measurement Data Orientation Rotation of Sensor-Coordinate System x y z accordin to a superordinate Coordinate System X Y Z Different representations of orientation Elements of the Direction-Cosine-Matrix Euler anles Quaternions Euler anles pitch θ roll φ and yaw ψ Disadvantae: imbal loc 4

5 Orientation and Smartphone Measurement Data Pitch and Roll Part of the orientation Used to level measurement data Necessary for calculation of yaw Accelerometer: a x a y a z Calculation of roll and pitch directly ϕ = ( a ) ( ) y z θ = f a f a a a x y z Gravity à low-pass filter tan( ϕ ) = tan ( θ ) = a a a y z y sin a x ( ϕ) + a cos( ϕ) z Gyroscope: ω x ω y ω z Calculation of the anular chanes of roll and pitch ( 1 1) ( 1) ϕ = f ω ω ω ϕ θ x y z θ = f ωy ωz ϕ Initial values are necessary ( sin( 1) cos( 1) ) tan( 1) ϕ = ω ϕ + ω ϕ θ + ω y z x ( ) sin( ) θ = ω cos ϕ ω ϕ y 1 z 1 ϕ = ϕ +Δt ϕ 1 θ = θ +Δt θ 1 5

6 Orientation and Smartphone Measurement Data Yaw Manetometer: m x m y m z Calculation of yaw directly ψ Gyroscope: ω x ω y ω z Calculation of anular chane of yaw = ( x y z ϕ θ) f m m m ( ) ψ = f ωy ωz ϕ 1 θ 1 Initial value necessary mx cosθ sinθsinφ sinθcosφ mx m Y 0 cosφ sinφ m = y m Z sinθ cosθsinφ cosθcosφ m z tan ( ψ ) m = m X Y ( sin( ) ( )) ( ) 1 1 cos 1 cos 1 ψ = ωy ϕ + ωz ϕ θ ψ = ψ +Δt ψ 1 Systematic effects Manetic disturbance Gyro drift 6

7 Kalman Filter for Orientation Determination System equations State vector contains φ θ and ψ Control inputs: yro measurements ω y and ω z Disturbance: w φ w θ and w ψ Random wal model f :ϕ =ϕ + w ϕ f2 :θ =θ 1+ w 1 θ f3:ψ =ψ 1+ w 1 ψ User is turnin User wals straiht Δt f :ψ =ψ + sin ϕ ω + cosϕ ω + w ψ y 1 z 1 cosθ 1 Propaation of covariance = T T + B B + C C T T T xx xx $$ 1 uu ww T B C = δ f3 δ f 3 1 δϕ δθ = δ f3 δ f3 0 δωy δω z =Δt T Transition Control Disturbance 7

8 Kalman Filter for Orientation Determination Indicators for chanin system equation Statistical test on innovations d Yaw-rate H : ( ) 0 0 E d = HA: E( d) 0 d 2 P{ χ1;1 α H0} = 1 α σ 2 ψ 2 d ψ Innovation test responds > 0.3s Yaw-rate > 30 /s 8

9 Kalman Filter for Orientation Determination $ $ 1 ay a z ϕ $ $ $ l = θ 2 ( x y z = a a a ) ψ $ $ $ 3 ( mx my mz ax ay az ) Separate Kalman filter for acceleration measurements à a$ $ $ x ay az Observations Pitch roll and yaw are introduced as directly observed aa $$ ( ) Adaptive covariance matrix of the observations aa $$ 33 T ll F F σ m * I 33 = 0 F 1 1 K $ a x mz 2 3 = K a $ x mz 3 3 K a$ x m z ( m ) 2 2 m = m +Δ 1 σ σ uv uv Δ = m 1 m m 1 9

10 Test Measurements and Results 3 scenarios One smartphone à Samsun Galaxy S4 Ground truth from indoors mobile toolit Evaluation Inner and outer accuracy Comparison to results from orientation sensor of android Scenario 2 Scenario 1 Scenario 3 10

11 Test Measurements and Results Scenario 1 Roll Random wal model Smooth filter result Time shift ~2s No reference values Smartphone was held in textin mode Results are reasonable Pitch 11

12 Test Measurements and Results Scenario 1 Yaw filter result No drift Minimized influence of manetic disturbance Constant offset Outer accuracy: 16.2 Mean deviation from round truth Inner accuracy: 5.4 Mean standard deviation of deviations Error in round truth 12

13 Test Measurements and Results Scenario 2 Reduced influence of manetic disturbance Outer accuracy 8.9 Inner accuracy

14 Test Measurements and Results Scenario 3 Outer accuracy 12.2 Inner accuracy 3.6 Properties of deviations Varyin most of the time 14

15 Results Comparison of Scenario 1 and Scenario 2 Offset and residuals differ 15

16 Conclusion and Outloo Conclusions Systematic effects of the sensors are minimized Turns are detected correctly Residuals are rarely reproducible and not normally distributed Filter result ~20% better than android solution Inner accuracy ~60% better than outer accuracy Scenario Outer accuracy Inner accuracy Orienta(on filter Android Orienta(on filter 1 (straiht) (straiht with turns) (polyonal repeatedly) Outloo Tests on different smartphones Improvin round truth Additional orientation information to overcome remainin systematic effects Adaption of test statistics 16

17 Andreas Ettliner Research Group Enineerin Geodesy Contact: LBS Conference

Application of state observers in attitude estimation using low-cost sensors

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

Inertial Odometry using AR Drone s IMU and calculating measurement s covariance

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

Attitude Estimation Version 1.0

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

More information

Quaternion based Extended Kalman Filter

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

Attitude Estimation for Augmented Reality with Smartphones

Attitude Estimation for Augmented Reality with Smartphones Attitude Estimation for Augmented Reality with Smartphones Thibaud Michel Pierre Genevès Hassen Fourati Nabil Layaïda Université Grenoble Alpes, INRIA LIG, GIPSA-Lab, CNRS June 13 th, 2017 http://tyrex.inria.fr/mobile/benchmarks-attitude

More information

Two dimensional rate gyro bias estimation for precise pitch and roll attitude determination utilizing a dual arc accelerometer array

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

EE565:Mobile Robotics Lecture 6

EE565: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 information

Linear Motion. Miroslav Mihaylov. February 13, 2014

Linear Motion. Miroslav Mihaylov. February 13, 2014 Linear Motion Miroslav Mihaylov February 13, 2014 1 Vector components Vector A has manitude A and direction θ with respect to the horizontal. On Fiure 1 we chose the eastbound as a positive x direction

More information

XI PHYSICS M. AFFAN KHAN LECTURER PHYSICS, AKHSS, K. https://promotephysics.wordpress.com

XI PHYSICS M. AFFAN KHAN LECTURER PHYSICS, AKHSS, K. https://promotephysics.wordpress.com XI PHYSICS M. AFFAN KHAN LECTURER PHYSICS, AKHSS, K affan_414@live.com https://promotephysics.wordpress.com [MOTION IN TWO DIMENSIONS] CHAPTER NO. 4 In this chapter we are oin to discuss motion in projectile

More information

Lecture 9: Modeling and motion models

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

Winter 2017 Ma 1b Analytical Problem Set 2 Solutions

Winter 2017 Ma 1b Analytical Problem Set 2 Solutions 1. (5 pts) From Ch. 1.10 in Apostol: Problems 1,3,5,7,9. Also, when appropriate exhibit a basis for S. Solution. (1.10.1) Yes, S is a subspace of V 3 with basis {(0, 0, 1), (0, 1, 0)} and dimension 2.

More information

Applications Linear Control Design Techniques in Aircraft Control I

Applications Linear Control Design Techniques in Aircraft Control I Lecture 29 Applications Linear Control Design Techniques in Aircraft Control I Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Topics Brief Review

More information

UAVBook Supplement Full State Direct and Indirect EKF

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

State Estimation and Motion Tracking for Spatially Diverse VLC Networks

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

(a) 1m s -2 (b) 2 m s -2 (c) zero (d) -1 m s -2

(a) 1m s -2 (b) 2 m s -2 (c) zero (d) -1 m s -2 11 th Physics - Unit 2 Kinematics Solutions for the Textbook Problems One Marks 1. Which one of the followin Cartesian coordinate system is not followed in physics? 5. If a particle has neative velocity

More information

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

Methodology of modelling the flexural-torsional vibrations in transient states of the rotating power transmission systems

Methodology of modelling the flexural-torsional vibrations in transient states of the rotating power transmission systems Methodology of modelling the flexural-torsional vibrations in transient states of the rotating power transmission systems Tomasz Matyja Faculty of Transport Department of Logistics and Aviation Technologies

More information

Full-dimension Attitude Determination Based on Two-antenna GPS/SINS Integrated Navigation System

Full-dimension Attitude Determination Based on Two-antenna GPS/SINS Integrated Navigation System Full-dimension Attitude Determination Based on Two-antenna GPS/SINS Integrated Navigation System Lifan Zhang Institute of Integrated Automation School of Electronic and Information Engineering Xi an Jiaotong

More information

State Estimation for Autopilot Control of Small Unmanned Aerial Vehicles in Windy Conditions

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

CE 530 Molecular Simulation

CE 530 Molecular Simulation CE 530 Molecular Simulation Lecture 7 Beyond Atoms: Simulating Molecules David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu Review Fundamentals units, properties, statistical

More information

UAV Coordinate Frames and Rigid Body Dynamics

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

Road bank angle estimation for three wheel tilting vehicle using multi model estimator

Road bank angle estimation for three wheel tilting vehicle using multi model estimator Road ban angle estimation for three wheel tilting vehicle using multi model estimator Kicheol Jeong, Seibum Choi,a KAIST, Department of Mechanical Engineering, 3 Deajeon, Republic of Korea Abstract. In

More information

KINEMATICS PREVIOUS EAMCET BITS ENGINEERING PAPER

KINEMATICS PREVIOUS EAMCET BITS ENGINEERING PAPER KINEMATICS PREVIOUS EAMCET BITS ENGINEERING PAPER. A body is projected vertically upwards at time t = 0 and is seen at a heiht at time t and t seconds durin its fliht. The maximum heiht attained is [ =

More information

KALMAN FILTER IMPLEMENTATION FOR UNMANNED AERIAL VEHICLES NAVIGATION DEVELOPED WITHIN A GRADUATE COURSE

KALMAN FILTER IMPLEMENTATION FOR UNMANNED AERIAL VEHICLES NAVIGATION DEVELOPED WITHIN A GRADUATE COURSE ALMAN FILER IMPLEMENAION FOR UNMANNED AERIAL VEHICLES NAVIGAION DEVELOPED WIHIN A GRADUAE COURSE Vasko Sazdovski, atjana olemishevska-guulovska, Mile Stankovski Institute of ASE at Faculty of EE,, St.

More information

NAWCWPNS TM 8128 CONTENTS. Introduction Two-Dimensinal Motion Three-Dimensional Motion Nonrotating Spherical Earth...

NAWCWPNS TM 8128 CONTENTS. Introduction Two-Dimensinal Motion Three-Dimensional Motion Nonrotating Spherical Earth... CONTENTS Introduction... 3 Two-Dimensinal Motion... 3 Three-Dimensional Motion... 5 Nonrotating Spherical Earth...10 Rotating Spherical Earth...12 WGS84...14 Conclusion...14 Appendixes: A. Kalman Filter...15

More information

Motion Locus Analysis to Detect Rotation

Motion Locus Analysis to Detect Rotation International Journal of Information and Electronics Engineering, Vol. 7, No. 6, November 07 Motion Locus Analysis to Detect Rotation Toshiki Iso Abstract For mobile applications such as command interfaces

More information

NAVAL POSTGRADUATE SCHOOL THESIS

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

An Adaptive Filter for a Small Attitude and Heading Reference System Using Low Cost Sensors

An Adaptive Filter for a Small Attitude and Heading Reference System Using Low Cost Sensors An Adaptive Filter for a Small Attitude and eading Reference System Using Low Cost Sensors Tongyue Gao *, Chuntao Shen, Zhenbang Gong, Jinjun Rao, and Jun Luo Department of Precision Mechanical Engineering

More information

Kinetics of Spatial Mechanisms: Kinetics Equation

Kinetics of Spatial Mechanisms: Kinetics Equation Lecture Note (6): Kinetics of Spatial Mechanisms: Kinetics Equation Objectives - Kinetics equations - Problem solving procedure - Example analysis Kinetics Equation With respect to the center of gravity

More information

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

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

More information

Fundamentals of attitude Estimation

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

Problem Set 5 Math 213, Fall 2016

Problem Set 5 Math 213, Fall 2016 Problem Set 5 Math 213, Fall 216 Directions: Name: Show all your work. You are welcome and encouraged to use Mathematica, or similar software, to check your answers and aid in your understanding of the

More information

Inertial Odometry on Handheld Smartphones

Inertial Odometry on Handheld Smartphones Inertial Odometry on Handheld Smartphones Arno Solin 1 Santiago Cortés 1 Esa Rahtu 2 Juho Kannala 1 1 Aalto University 2 Tampere University of Technology 21st International Conference on Information Fusion

More information

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

( ) = x( u, v) i + y( u, v) j + z( u, v) k

( ) = x( u, v) i + y( u, v) j + z( u, v) k Math 8 ection 16.6 urface Integrals The relationship between surface integrals and surface area is much the same as the relationship between line integrals and arc length. uppose f is a function of three

More information

Presenter: Siu Ho (4 th year, Doctor of Engineering) Other authors: Dr Andy Kerr, Dr Avril Thomson

Presenter: Siu Ho (4 th year, Doctor of Engineering) Other authors: Dr Andy Kerr, Dr Avril Thomson The development and evaluation of a sensor-fusion and adaptive algorithm for detecting real-time upper-trunk kinematics, phases and timing of the sit-to-stand movements in stroke survivors Presenter: Siu

More information

Refinements to the General Methodology Behind Strapdown Airborne Gravimetry

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

Miscellaneous. Regarding reading materials. Again, ask questions (if you have) and ask them earlier

Miscellaneous. 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 information

Dynamic Modeling of Fixed-Wing UAVs

Dynamic Modeling of Fixed-Wing UAVs Autonomous Systems Laboratory Dynamic Modeling of Fixed-Wing UAVs (Fixed-Wing Unmanned Aerial Vehicles) A. Noth, S. Bouabdallah and R. Siegwart Version.0 1/006 1 Introduction Dynamic modeling is an important

More information

Adaptive Kalman Filter for MEMS-IMU based Attitude Estimation under External Acceleration and Parsimonious use of Gyroscopes

Adaptive Kalman Filter for MEMS-IMU based Attitude Estimation under External Acceleration and Parsimonious use of Gyroscopes Author manuscript, published in "European Control Conference ECC (214" Adaptive Kalman Filter for MEMS-IMU based Attitude Estimation under External Acceleration and Parsimonious use of Gyroscopes Aida

More information

Inertial Navigation System Aiding Using Vision

Inertial Navigation System Aiding Using Vision Air Force Institute of Technology AFIT Scholar Theses and Dissertations 3-21-2013 Inertial Navigation System Aiding Using Vision James O. Quarmyne Follow this and additional works at: https://scholar.afit.edu/etd

More information

Unscented Kalman filter and Magnetic Angular Rate Update (MARU) for an improved Pedestrian Dead-Reckoning

Unscented Kalman filter and Magnetic Angular Rate Update (MARU) for an improved Pedestrian Dead-Reckoning Unscented Kalman filter and Magnetic Angular Rate Update (MARU) for an improved Pedestrian Dead-Reckoning Francisco Zampella, Mohammed Khider, Patrick Robertson and Antonio Jiménez Centre for Automation

More information

Lecture 10 Reprise and generalized forces

Lecture 10 Reprise and generalized forces Lecture 10 Reprise and generalized forces The Lagrangian Holonomic constraints Generalized coordinates Nonholonomic constraints Generalized forces we haven t done this, so let s start with it Euler Lagrange

More information

Unscented Kalman filter and Magnetic Angular Rate Update (MARU) for an improved Pedestrian Dead-Reckoning

Unscented Kalman filter and Magnetic Angular Rate Update (MARU) for an improved Pedestrian Dead-Reckoning Submitted to the 212 IEEE/ION Position, Location and Navigation Symposium Unscented Kalman filter and Magnetic Angular Rate Update (MARU) for an improved Pedestrian Dead-Reckoning Francisco Zampella, Mohammed

More information

Engineering Notebook

Engineering Notebook Engineering Notebook miscellaneous problems and solutions compiled by C Bond Vol1 no9 Actuator with zero bend angle Rotary Actuator Equations Given two points at which the track radius and gap skew angle

More information

VN-100 Velocity Compensation

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

Research Article Low-Cost MEMS-Based Pedestrian Navigation Technique for GPS-Denied Areas

Research Article Low-Cost MEMS-Based Pedestrian Navigation Technique for GPS-Denied Areas Journal of Sensors Volume 213, Article ID 1979, 1 pages http://dx.doi.org/1.1155/213/1979 Research Article Low-Cost MEMS-Based Pedestrian Navigation Technique for GPS-Denied Areas Abdelrahman Ali and Naser

More information

Attitude determination method using single-antenna GPS, Gyro and Magnetometer

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

ECE Spring Prof. David R. Jackson ECE Dept. Notes 10

ECE Spring Prof. David R. Jackson ECE Dept. Notes 10 ECE 6345 Spring 215 Prof. David R. Jackson ECE Dept. Notes 1 1 Overview In this set of notes we derive the far-field pattern of a circular patch operating in the dominant TM 11 mode. We use the magnetic

More information

PART ONE DYNAMICS OF A SINGLE PARTICLE

PART ONE DYNAMICS OF A SINGLE PARTICLE PART ONE DYNAMICS OF A SINGLE PARTICLE 1 Kinematics of a Particle 1.1 Introduction One of the main goals of this book is to enable the reader to take a physical system, model it by using particles or rigid

More information

with any accuracy. The parameter a of cubic substances is directly proportional to the spacing d of any particular set of lattice planes:

with any accuracy. The parameter a of cubic substances is directly proportional to the spacing d of any particular set of lattice planes: The change in solute concentration or temperature produce only a small change in lattice parameter precise parameter measurements is needed to measure these quantities with any accuracy. The parameter

More information

Motion in Three Dimensions

Motion in Three Dimensions Motion in Three Dimensions We ve learned about the relationship between position, velocity and acceleration in one dimension Now we need to extend those ideas to the three-dimensional world In the 1-D

More information

Lecture 37: Principal Axes, Translations, and Eulerian Angles

Lecture 37: Principal Axes, Translations, and Eulerian Angles Lecture 37: Principal Axes, Translations, and Eulerian Angles When Can We Find Principal Axes? We can always write down the cubic equation that one must solve to determine the principal moments But if

More information

THE USE OF KALMAN FILTRATION TO ESTIMATE CHANGES OF TRUNK INCLINATION ANGLE DURING WEIGHTLIFTING 1. INTRODUCTION

THE USE OF KALMAN FILTRATION TO ESTIMATE CHANGES OF TRUNK INCLINATION ANGLE DURING WEIGHTLIFTING 1. INTRODUCTION JOURNAL OF MEDICAL INFORMATICS & TECHNOLOGIES Vol. 15/2010, ISSN 1642-6037 Kalman filtration, filter algorithm, accelerometric sensor Grzegorz SAPOTA 1, Anna SAPOTA 1, Zygmunt WRÓBEL 1 THE USE OF KALMAN

More information

Lesson Rigid Body Dynamics

Lesson Rigid Body Dynamics Lesson 8 Rigid Body Dynamics Lesson 8 Outline Problem definition and motivations Dynamics of rigid bodies The equation of unconstrained motion (ODE) User and time control Demos / tools / libs Rigid Body

More information

THE COMPOUND ANGLE IDENTITIES

THE COMPOUND ANGLE IDENTITIES TRIGONOMETRY THE COMPOUND ANGLE IDENTITIES Question 1 Prove the validity of each of the following trigonometric identities. a) sin x + cos x 4 4 b) cos x + + 3 sin x + 2cos x 3 3 c) cos 2x + + cos 2x cos

More information

Lecture. Aided INS EE 570: Location and Navigation. 1 Overview. 1.1 ECEF as and Example. 1.2 Inertial Measurements

Lecture. 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 information

OUTCOME 2 KINEMATICS AND DYNAMICS

OUTCOME 2 KINEMATICS AND DYNAMICS Unit 60: Dynamics of Machines Unit code: H/601/1411 QCF Level:4 Credit value:15 OUTCOME 2 KINEMATICS AND DYNAMICS TUTORIAL 3 GYROSCOPES 2 Be able to determine the kinetic and dynamic parameters of mechanical

More information

Poncelet s porism and periodic triangles in ellipse 1

Poncelet s porism and periodic triangles in ellipse 1 Poncelet s porism and periodic triangles in ellipse 1 Vladimir Georgiev, Veneta Nedyalkova 1 Small historical introduction One of the most important and beautiful theorems in projective geometry is that

More information

Robot Control Basics CS 685

Robot Control Basics CS 685 Robot Control Basics CS 685 Control basics Use some concepts from control theory to understand and learn how to control robots Control Theory general field studies control and understanding of behavior

More information

Quadrotor Modeling and Control for DLO Transportation

Quadrotor Modeling and Control for DLO Transportation Quadrotor Modeling and Control for DLO Transportation Thesis dissertation Advisor: Prof. Manuel Graña Computational Intelligence Group University of the Basque Country (UPV/EHU) Donostia Jun 24, 2016 Abstract

More information

with Application to Autonomous Vehicles

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

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

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

More information

Practice Differentiation Math 120 Calculus I Fall 2015

Practice Differentiation Math 120 Calculus I Fall 2015 . x. Hint.. (4x 9) 4x + 9. Hint. Practice Differentiation Math 0 Calculus I Fall 0 The rules of differentiation are straightforward, but knowing when to use them and in what order takes practice. Although

More information

Rotations about the coordinate axes are easy to define and work with. Rotation about the x axis by angle θ is. R x (θ) = 0 cos θ sin θ 0 sin θ cos θ

Rotations about the coordinate axes are easy to define and work with. Rotation about the x axis by angle θ is. R x (θ) = 0 cos θ sin θ 0 sin θ cos θ Euler Angle Formulas David Eberly Magic Software, Inc. http://www.magic-software.com Created: December 1, 1999 1 Introduction Rotations about the coordinate axes are easy to define and work with. Rotation

More information

Understanding trigonometric functions and complex numbers

Understanding trigonometric functions and complex numbers Understanding trigonometric functions and complex numbers Qian Kemao School of Computer Engineering, Nanyang Technological University, Singapore 639798 mkmqian@ntu.edu.sg ABSTRACT Two topics, trigonometric

More information

Inertial navigation for divers

Inertial navigation for divers Inertial navigation for divers Problem presented by Nick Bushell VR Technology Executive Summary A SCUBA diver would like to know his/her position underwater, relative to the dive start position. GPS is

More information

We 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. International authors and editors We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists 3,800 116,000 120M Open access books available International authors and editors Downloads Our

More information

IMU Filter. Michael Asher Emmanuel Malikides November 5, 2011

IMU Filter. Michael Asher Emmanuel Malikides November 5, 2011 IMU Filter Michael Asher Emmanuel Malikides November 5, 2011 Abstract Despite the ubiquitousness of GPS devices, on board inertial navigation remains important. An IMU like the Sparkfun Ultimate IMU used,

More information

CBE 6333, R. Levicky 1. Orthogonal Curvilinear Coordinates

CBE 6333, R. Levicky 1. Orthogonal Curvilinear Coordinates CBE 6333, R. Levicky 1 Orthogonal Curvilinear Coordinates Introduction. Rectangular Cartesian coordinates are convenient when solving problems in which the geometry of a problem is well described by the

More information

CS491/691: Introduction to Aerial Robotics

CS491/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 information

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

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

More information

Continuous Preintegration Theory for Graph-based Visual-Inertial Navigation

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

Math review. Math review

Math review. Math review Math review 1 Math review 3 1 series approximations 3 Taylor s Theorem 3 Binomial approximation 3 sin(x), for x in radians and x close to zero 4 cos(x), for x in radians and x close to zero 5 2 some geometry

More information

Using Rotational Matrices to Generate Triplex Algebras

Using Rotational Matrices to Generate Triplex Algebras Using Rotational Matrices to Generate Triplex Algebras The real numbers are the dependable breadwinner of the family, the complete ordered field we all rely on. The complex numbers are a slightly flashier

More information

Determining the 3D Spin Axis from TrackMan Data (updated, March 31, 2015)

Determining the 3D Spin Axis from TrackMan Data (updated, March 31, 2015) Determining the 3D Spin Axis from TrackMan Data (updated, March 31, 2015) Alan M. Nathan Department of Physics, University of Illinois, Urbana, IL 61801 A technique is developed to determine the direction

More information

Acceleration and rotation in a pendulum ride, measured using an iphone 4

Acceleration and rotation in a pendulum ride, measured using an iphone 4 F EATURES www.iop.org/journals/physed Acceleration and rotation in a pendulum ride, measured using an iphone 4 Ann-Marie Pendrill 1,2 and Johan Rohlén 1 1 Department of Physics, University of Gothenburg,

More information

Rigid Body Rotation. Speaker: Xiaolei Chen Advisor: Prof. Xiaolin Li. Department of Applied Mathematics and Statistics Stony Brook University (SUNY)

Rigid Body Rotation. Speaker: Xiaolei Chen Advisor: Prof. Xiaolin Li. Department of Applied Mathematics and Statistics Stony Brook University (SUNY) Rigid Body Rotation Speaker: Xiaolei Chen Advisor: Prof. Xiaolin Li Department of Applied Mathematics and Statistics Stony Brook University (SUNY) Content Introduction Angular Velocity Angular Momentum

More information

Adaptive Kalman Filter for MEMS-IMU based Attitude Estimation under External Acceleration and Parsimonious use of Gyroscopes

Adaptive Kalman Filter for MEMS-IMU based Attitude Estimation under External Acceleration and Parsimonious use of Gyroscopes Adaptive Kalman Filter for MEMS-IMU based Attitude Estimation under External Acceleration and Parsimonious use of Gyroscopes Aida Mani, Hassen Fourati, Alain Kibangou To cite this version: Aida Mani, Hassen

More information

Real-time Attitude Estimation Techniques Applied to a Four Rotor Helicopter

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

On the Observability and Self-Calibration of Visual-Inertial Navigation Systems

On the Observability and Self-Calibration of Visual-Inertial Navigation Systems Center for Robotics and Embedded Systems University of Southern California Technical Report CRES-08-005 R B TIC EMBEDDED SYSTEMS LABORATORY On the Observability and Self-Calibration of Visual-Inertial

More information

Question 1: A particle starts at rest and moves along a cycloid whose equation is. 2ay y a

Question 1: A particle starts at rest and moves along a cycloid whose equation is. 2ay y a Stephen Martin PHYS 10 Homework #1 Question 1: A particle starts at rest and moves along a cycloid whose equation is [ ( ) a y x = ± a cos 1 + ] ay y a There is a gravitational field of strength g in the

More information

Physics 2A Chapter 10 - Rotational Motion Fall 2018

Physics 2A Chapter 10 - Rotational Motion Fall 2018 Physics A Chapter 10 - Rotational Motion Fall 018 These notes are five pages. A quick summary: The concepts of rotational motion are a direct mirror image of the same concepts in linear motion. Follow

More information

Mechanics Physics 151

Mechanics Physics 151 Mechanics Phsics 151 Lecture 8 Rigid Bod Motion (Chapter 4) What We Did Last Time! Discussed scattering problem! Foundation for all experimental phsics! Defined and calculated cross sections! Differential

More information

Calibration of a magnetometer in combination with inertial sensors

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

Solution of Math132 Exam1

Solution of Math132 Exam1 Solution of Math3 Exam. ( %) (a) A bacteria has population at time t = an its rate of growth is (t 3 + t + ) bacteria per hour after t hours. What is the population after 4 hours? Suppose p(t) is the bacteria

More information

EE 570: Location and Navigation

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

UNCLASSIFIED. Control of Spinning Symmetric Airframes. Curtis P. Mracek Max Stafford and Mike Unger Raytheon Systems Company Tucson, AZ 85734

UNCLASSIFIED. Control of Spinning Symmetric Airframes. Curtis P. Mracek Max Stafford and Mike Unger Raytheon Systems Company Tucson, AZ 85734 Control of Spinning Symmetric Airframes Curtis P. Mracek Max Stafford and Mike Unger Raytheon Systems Company Tucson, AZ 85734 Several missiles and projectiles have periods during their flight where the

More information

PH 451/551 Quantum Mechanics Capstone Winter 201x

PH 451/551 Quantum Mechanics Capstone Winter 201x These are the questions from the W7 exam presented as practice problems. The equation sheet is PH 45/55 Quantum Mechanics Capstone Winter x TOTAL POINTS: xx Weniger 6, time There are xx questions, for

More information

Multi-layer Flight Control Synthesis and Analysis of a Small-scale UAV Helicopter

Multi-layer Flight Control Synthesis and Analysis of a Small-scale UAV Helicopter Multi-layer Flight Control Synthesis and Analysis of a Small-scale UAV Helicopter Ali Karimoddini, Guowei Cai, Ben M. Chen, Hai Lin and Tong H. Lee Graduate School for Integrative Sciences and Engineering,

More information

Introduction to Mobile Robotics Probabilistic Motion Models

Introduction to Mobile Robotics Probabilistic Motion Models Introduction to Mobile Robotics Probabilistic Motion Models Wolfram Burgard 1 Robot Motion Robot motion is inherently uncertain. How can we model this uncertainty? 2 Dynamic Bayesian Network for Controls,

More information

1 Kalman Filter Introduction

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 information

arxiv: v1 [math.oc] 12 May 2018

arxiv: v1 [math.oc] 12 May 2018 Observability Analysis of Aided NS with Heterogeneous Features of Points, Lines and Planes Yulin Yang - yuyang@udeledu uoquan Huang - ghuang@udeledu Department of Mechanical Engineering University of Delaware,

More information

Laws of gyroscopes / cardanic gyroscope

Laws of gyroscopes / cardanic gyroscope Principle If the axis of rotation of the force-free gyroscope is displaced slightly, a nutation is produced. The relationship between precession frequency or nutation frequency and gyrofrequency is examined

More information

A Machine Learning Approach for Dead-Reckoning Navigation at Sea Using a Single Accelerometer

A Machine Learning Approach for Dead-Reckoning Navigation at Sea Using a Single Accelerometer A Machine Learning Approach for Dead-Reckoning Navigation at Sea Using a Single Accelerometer 1 Roee Diamant, Yunye Jin Electrical and Computer Engineering, The University of British Columbia Institute

More information

A Low-Cost GPS Aided Inertial Navigation System for Vehicular Applications

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

Aerodynamics SYST 460/560. George Mason University Fall 2008 CENTER FOR AIR TRANSPORTATION SYSTEMS RESEARCH. Copyright Lance Sherry (2008)

Aerodynamics SYST 460/560. George Mason University Fall 2008 CENTER FOR AIR TRANSPORTATION SYSTEMS RESEARCH. Copyright Lance Sherry (2008) Aerodynamics SYST 460/560 George Mason University Fall 2008 1 CENTER FOR AIR TRANSPORTATION SYSTEMS RESEARCH Copyright Lance Sherry (2008) Ambient & Static Pressure Ambient Pressure Static Pressure 2 Ambient

More information

NEW EUMETSAT POLAR SYSTEM ATTITUDE MONITORING SOFTWARE

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

More information

Christoffel Symbols. 1 In General Topologies. Joshua Albert. September 28, W. First we say W : λ n = x µ (λ) so that the world

Christoffel Symbols. 1 In General Topologies. Joshua Albert. September 28, W. First we say W : λ n = x µ (λ) so that the world Christoffel Symbols Joshua Albert September 28, 22 In General Topoloies We have a metric tensor nm defined by, Note by some handy theorem that for almost any continuous function F (L), equation 2 still

More information

Physics 312, Winter 2007, Practice Final

Physics 312, Winter 2007, Practice Final Physics 312, Winter 2007, Practice Final Time: Two hours Answer one of Question 1 or Question 2 plus one of Question 3 or Question 4 plus one of Question 5 or Question 6. Each question carries equal weight.

More information