EE 570: Location and Navigation

Size: px
Start display at page:

Download "EE 570: Location and Navigation"

Transcription

1 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 Egieerig Departmet Embry-Riddle Aeroautical Uivesity Prescott, Arizoa, USA April 4, 2016 Attitude NAV Velocity Error Nav Positio Error Summary Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

2 NAV Attitude Error C b = C b (Ωb ib Ωb i) = d dt [ (I + [δψ ] b ])Ĉ b = (I + [δ ψ b ])Ĉ b Ωb b = [δ ψ b ]Ĉ b + (I + [δ ψ b ]) Ĉ b = (I + [δ ψ b ])Ĉ b (ˆΩ b b + δωb ib δωb i) (I + [δ ψ b ])Ĉ b ˆΩ b b + Ĉ b (δωb ib δωb i) Attitude NAV Velocity Error Nav Positio Error Summary Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

3 NAV Attitude Error C b = C b (Ωb ib Ωb i) = d dt [ (I + [δψ ] b ])Ĉ b = (I + [δ ψ b ])Ĉ b Ωb b = [δ ψ b ]Ĉ b + (I + [δ ψ b ]) Ĉ b = (I + [δ ψ b ])Ĉ b (ˆΩ b b + δωb ib δωb i) (I + [δ ψ b ])Ĉ b ˆΩ b b + Ĉ b (δωb ib δωb i) Attitude NAV Velocity Error Nav Positio Error Summary Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

4 NAV Attitude Error C b = C b (Ωb ib Ωb i) = d dt [ (I + [δψ ] b ])Ĉ b = (I + [δ ψ b ])Ĉ b Ωb b = [δ ψ b ]Ĉ b + (I + [δ ψ b ]) Ĉ b = (I + [δ ψ b ])Ĉ b (ˆΩ b b + δωb ib δωb i) (I + [δ ψ b ])Ĉ b ˆΩ b b + Ĉ b (δωb ib δωb i) Attitude NAV Velocity Error Nav Positio Error Summary Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

5 NAV Attitude Error C b = C b (Ωb ib Ωb i) = d dt [ (I + [δψ ] b ])Ĉ b = (I + [δψ b ])Ĉ b Ωb b = [δ ψ b ]Ĉb + (I + [δ ψ b ]) Ĉ b = (I + [δψ b ])Ĉ b (ˆΩ b b + δωb ib δωb i) (I + [δψ b ])Ĉ b ˆΩ b b + Ĉ b (δωb ib δωb i) [δψ b ]δωb b 0 Attitude NAV Velocity Error Nav Positio Error Summary Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

6 NAV Attitude Error C b = C b (Ωb ib Ωb i) = d dt [ (I + [δψ ] b ])Ĉ b = (I + [δψ b ])Ĉ b Ωb b = [δ ψ b ]Ĉb + (I + [δ ψ b ]) Ĉ b = (I + [δψ b ])Ĉ b (ˆΩ b b + δωb ib δωb i) = (I + [δψ b ])Ĉ b ˆΩ b b + Ĉ b (δωb ib δωb i) Attitude NAV Velocity Error Nav Positio Error Summary Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

7 NAV Attitude Error C b = C b (Ωb ib Ωb i) = d dt [ (I + [δψ ] b ])Ĉ b = (I + [δ ψ b ])Ĉ b Ωb b = [δ ψ b ]Ĉ b + (I + [δ ψ b ]) Ĉ b = (I + [δ ψ b ])Ĉ b (ˆΩ b b + δωb ib δωb i) (I + [δ ψ b ])Ĉ b ˆΩ b b + Ĉ b (δωb ib δωb i) [δ ψ b ] = Ĉ b (δωb b δωb i)ĉ b = [Ĉ b (δ ω b ib δ ω b i) ] (1) Attitude NAV Velocity Error Nav Positio Error Summary Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

8 NAV Attitude Error C b = C b (Ωb ib Ωb i) = d dt [ (I + [δψ ] b ])Ĉ b = (I + [δ ψ b ])Ĉ b Ωb b = [δ ψ b ]Ĉ b + (I + [δ ψ b ]) Ĉ b = (I + [δ ψ b ])Ĉ b (ˆΩ b b + δωb ib δωb i) (I + [δ ψ b ])Ĉ b ˆΩ b b + Ĉ b (δωb ib δωb i) [δ ψ b ] = Ĉ b (δωb b δωb i)ĉ b = [Ĉ b (δ ω b ib δ ω b i) ] (1) δ ψ b = Ĉ b (δωb ib δ ω b i) (2) Attitude NAV Velocity Error Nav Positio Error Summary Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

9 NAV Attitude Error (cot.) Computig δ ω b i Recall that ω b i = C b ω i Expressig the above equatio i terms of estimates we get ˆ ω i b + δ ω b i = Ĉ b (I [δψ b ])(ˆ ω i + δ ω i) δ ω b i Ĉ b (δ ω i [δψ b ]ˆ ω i) Substitutig this result i Equatio 2 δ ψ b = ˆ Ω i δψ b + Ĉ b δ ω b ib δ ω i (3) Attitude NAV Velocity Error Nav Positio Error Summary Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

10 NAV Attitude Error (cot.) Computig δ ω i Usig Taylor series ad retaiig the first order terms where r eb = [L b, λ b, h b ] T ad δ ω i = ˆ ω i δ r ˆ r eb + ˆ ω i δ v eb ˆ v eb eb ω ie cos ˆL b + ˆ v eb,e R ˆ ω i = ˆ ω ie + ˆ ω E (ˆL b )+ĥ b e = ˆ v eb,n R N (ˆL b )+ĥb ω ie si ˆL b ˆ v eb,e ta ˆL b R E (ˆL b )+ĥ b Attitude NAV Velocity Error Nav Positio Error Summary Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

11 NAV Attitude Error (cot.) Defie terms for simplicity ˆ ω i ˆ r eb ˆ ω i ˆ v eb ˆ Ω i = [(ˆ ω ie + ˆ ω e) ] = F ψψ (4) RE (ˆLb)+ĥb = = F ψv (5) RN(ˆLb)+ĥb ta 0 ˆLb 0 RE (ˆLb)+ĥb ˆ v eb,e ω ie si ˆL b 0 (RE (ˆLb)+ĥb)2 = ˆ v 0 0 eb,n = F ψr (6) (RN(ˆLb)+ĥb)2 ω ie cos ˆL ˆ v b + eb,e 0 ˆ v eb,e ta ˆLb (RE (ˆLb)+ĥb) cos 2 ˆLb (RE (ˆLb)+ĥb) 2 Attitude NAV Velocity Error Nav Positio Error Summary Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

12 NAV Attitude Error (cot.) Fial Expressio δ ψ b = F ψψ δ ψ b + F ψv δ v eb + F ψr δ r eb + Ĉ b δ ω b ib (7) where F ψψ, F ψv ad F ψr are defied by Equatios 4, 5 ad 6, respectively. Attitude NAV Velocity Error Nav Positio Error Summary Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

13 Velocity v eb = C b f b ib + g b (Ω e + 2Ω ie) v e eb (8) ˆ v eb = Ĉ b ˆ f b ib + ˆ g b (ˆΩ e + 2ˆΩ ie)ˆ v eb = (I [δ ψ b ])C b ( f b ib δ f b ib ) + ˆ g b (ˆΩ e + 2ˆΩ ie)ˆ v eb (9) Computig the δ v eb we obtai δ v eb = v eb ˆ v eb = [δ ψ b ]C b f b ib + Ĉ b δ f b ib + δ g b + (ˆΩ e + 2ˆΩ ie)δ v eb (δω e + 2δΩ ie)ˆ v eb Attitude NAV Velocity Error Nav Positio Error Summary Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

14 Velocity Cot. (ˆΩ e + 2ˆΩ ie) The term (ˆΩ e + 2ˆΩ ie ) is dervied as 0 2ω ie si ˆL b ˆ v eb,e ta ˆLb RE (ˆLb)+ĥb RN(ˆLb)+ĥb [(ˆ ω e + 2ˆ ω ie) ] = 2ω ie si ˆL b + ˆ v eb,e ta ˆLb 0 2ω ie cos ˆL b + ˆ v eb,e = F ΩΩ (10) RE (ˆLb)+ĥb RE (ˆLb)+ĥb ˆ v eb,n 2ω ie cos ˆL b ˆ v eb,e 0 RN(ˆLb)+ĥb RE (ˆLb)+ĥb ˆ v eb,n Attitude NAV Velocity Error Nav Positio Error Summary Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

15 Velocity Cot. (δ ˆΩ e + 2δ ˆΩ ie)ˆ v eb The term (δ ˆΩ e + 2δ ˆΩ ie )ˆ v eb is dervied as (δ ˆΩ e + 2δ ˆΩ ie)ˆ v eb = F ΩΩ ˆ r eb δ r ebˆ v eb + F ΩΩ δ v ebˆ v ˆ v eb (11) eb Attitude NAV Velocity Error Nav Positio Error Summary Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

16 Velocity Cot. Gravity Error δ g b 2g 0(ˆL b ) r e es (ˆL b ) δh b (12) Attitude NAV Velocity Error Nav Positio Error Summary Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

17 Velocity Error Cot. After tos of algebra δ v eb = F vψδ ψ b + F vvδ v eb + F vr δ r eb + Ĉ b δ f b ib (13) where F vψ, F vv ad F vr are defied by Equatios 14, 15 ad 16, respectively. Attitude NAV Velocity Error Nav Positio Error Summary Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

18 Velocity Cot. F vψ, F vv ad F vr F vr = ˆv eb,d (RN(ˆLb)+ĥb) F vv = ˆv eb,e ta ˆLb + 2ω (RE (ˆLb)+ĥb) ie si ˆL b 2ˆv eb,n RN(ˆLb)+ĥb [ ] F vψ = (Ĉ b ˆ f b ib ) (ˆv eb,e )2 sec 2 ˆLb 2ˆv ω ie cos ˆL b 0 RE (ˆLb)+ĥb eb,e 2ˆv eb,e ta ˆLb 2ω (RE (ˆLb)+ĥb) ie si ˆL ˆv eb,n b RN(ˆLb)+ĥb ˆv eb,n ta ˆLb+ˆv eb,d ˆv eb,e + 2ω RN(ˆLb)+ĥb (RE (ˆLb)+ĥb) ie cos ˆL b 2ˆv eb,e 2ω (RE (ˆLb)+ĥb) ie cos ˆL b 0 ˆv eb,n ˆv eb,e sec2 ˆLb + 2ˆv RE (ˆLb)+ĥb eb,n ω ie cos ˆL b 2ˆv eb,d ω ie si ˆL b 0 ˆv eb,n ˆv eb,e ta ˆLb+ˆv eb,e ˆv eb,d (RE (ˆLb)+ĥb)2 2ˆv eb,e ω ie si ˆL b 0 (ˆv eb,e )2 ta ˆLb ˆv eb,n ˆv eb,d (RE (ˆLb)+ĥb)2 (RN(ˆLb)+ĥb)2 ( ) (ˆv eb,e )2 + (ˆv eb,n )2 2g0(ˆLb) (RE (ˆLb)+ĥb) 2 (RN(ˆLb)+ĥb) 2 res e (ˆLb) (14) (15) (16) Attitude NAV Velocity Error Nav Positio Error Summary Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

19 Positio ˆ r eb = ˆ v eb,n R N (ˆL b )+ĥb ˆ v eb,e cos ˆL b (R E ( ˆL ) b )+ĥ b ) ˆ v eb,d (17) Computig δ r eb usig taylor series expasio ad retaiig oly the first order terms δ r eb = F rψδ ψ b + F rv δ v eb + F rr δ r eb (18) where F rψ, F rv ad F rr are defied by Equatios 19, 20 ad 21, respectively. Attitude NAV Velocity Error Nav Positio Error Summary Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

20 Positio Cot. F rψ, F rv ad F rr F rψ = (19) R N(ˆL b)+ĥb F rv = (R E (ˆL b)+ĥ b) cos ˆL b (20) 0 0 F rr = ˆv eb,e ta ˆL b 0 (R N(ˆL b)+ĥ b) cos ˆL b ˆv eb,n (R N(ˆL b)+ĥ b) 2 ˆv eb,e (R N(ˆL b)+ĥ b) 2 cos ˆL b (21) Attitude NAV Velocity Error Nav Positio Error Summary Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

21 Summary - i terms of δ f b ib, δ ω b ib δ ψ b F δ v ψψ F ψv F ψr δ ψb 0 Ĉ b b = F vψ F vv F vr δ v b + Ĉ b 0 δ f b ib (22) δ ω b δ r b F rv F rr δ r b ib 0 0 Attitude NAV Velocity Error Nav Positio Error Summary Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

22 Notatio Used Truth value Measured value x x Estimated or computed value ˆ x Error δ x = x ˆ x Basic Defiitios Liearizatio Iertial Measuremets Attitude Error Estimates of Sesor Measuremets Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

23 Notatio Used Truth value Measured value Estimated or computed value x x Nothig above ˆ x Error δ x = x ˆ x Basic Defiitios Liearizatio Iertial Measuremets Attitude Error Estimates of Sesor Measuremets Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

24 Notatio Used Truth value Measured value Estimated or computed value x x ˆ x Use tilde Error δ x = x ˆ x Basic Defiitios Liearizatio Iertial Measuremets Attitude Error Estimates of Sesor Measuremets Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

25 Notatio Used Truth value Measured value Estimated or computed value x x ˆ x Error δ x = x ˆ x Use hat Basic Defiitios Liearizatio Iertial Measuremets Attitude Error Estimates of Sesor Measuremets Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

26 Notatio Used Truth value Measured value x x Estimated or computed value ˆ x Error δ x δ x = x ˆ x ˆ x Basic Defiitios Liearizatio Iertial Measuremets Attitude Error Estimates of Sesor Measuremets Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

27 Liearizatio usig Taylor Series Expasio Give a o-liear system x = f ( x, t) Basic Defiitios Liearizatio Iertial Measuremets Attitude Error Estimates of Sesor Measuremets Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

28 Liearizatio usig Taylor Series Expasio Give a o-liear system x = f ( x, t) Let s assume we have a estimate of x, i.e., ˆ x such that x = ˆ x + δ x x = ˆ x + δ x = f (ˆ x + δ x, t) (23) Basic Defiitios Liearizatio Iertial Measuremets Attitude Error Estimates of Sesor Measuremets Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

29 Liearizatio usig Taylor Series Expasio Give a o-liear system x = f ( x, t) Let s assume we have a estimate of x, i.e., ˆ x such that x = ˆ x + δ x x = ˆ x + δ x = f (ˆ x + δ x, t) (23) Usig Taylor series expasio f (ˆ x + δ x, t) = ˆ x f ( x, t) + δ x = f (ˆ x, t) + x δ x + H.O.T x=ˆ x ˆ x f ( x, t) + x δ x x=ˆ x Basic Defiitios Liearizatio Iertial Measuremets Attitude Error Estimates of Sesor Measuremets Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

30 Liearizatio usig Taylor Series Expasio Give a o-liear system x = f ( x, t) Let s assume we have a estimate of x, i.e., ˆ x such that x = ˆ x + δ x x = ˆ x + δ x = f (ˆ x + δ x, t) (23) Usig Taylor series expasio f (ˆ x + δ x, t) = ˆ x f ( x, t) + δ x = f (ˆ x, t) + x δ x + H.O.T = x=ˆ x ˆ x f ( x, t) + x δ x x=ˆ x Basic Defiitios Liearizatio Iertial Measuremets Attitude Error Estimates of Sesor Measuremets Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

31 Liearizatio usig Taylor Series Expasio Give a o-liear system x = f ( x, t) Let s assume we have a estimate of x, i.e., ˆ x such that x = ˆ x + δ x x = ˆ x + δ x = f (ˆ x + δ x, t) (23) Usig Taylor series expasio f (ˆ x + δ x, t) = ˆ x f ( x, t) + δ x = f (ˆ x, t) + x δ x + H.O.T x=ˆ x ˆ x f ( x, t) + x δ x x=ˆ x δ x f ( x, t) x δ x (24) x=ˆ x Basic Defiitios Liearizatio Iertial Measuremets Attitude Error Estimates of Sesor Measuremets Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

32 Actual Measuremets Iitially the accelerometer ad gyroscope measuremets, f b ib ad ω b ib, respectively, will be modeled as f b ib = f b ib + f b ib (25) ω b ib = ω b ib + ω b ib (26) where f b ib ad ω b ib are the specific force ad agular rates, respectively; ad f b ib ad ω b ib represets the errors. I later lectures we will discuss more detailed descriptio of these errors. Basic Defiitios Liearizatio Iertial Measuremets Attitude Error Estimates of Sesor Measuremets Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

33 Actual Measuremets Iitially the accelerometer ad gyroscope measuremets, f b ib ad ω b ib, respectively, will be modeled as f b ib = f b ib + f b ib (25) these terms may ω b ib = ω b ib + ω b be expaded further ib (26) where f b ib ad ω b ib are the specific force ad agular rates, respectively; ad f b ib ad ω b ib represets the errors. I later lectures we will discuss more detailed descriptio of these errors. Basic Defiitios Liearizatio Iertial Measuremets Attitude Error Estimates of Sesor Measuremets Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

34 Error Modelig Example Accelerometers f b ib = b a + (I + M a ) f b ib + l a + w a (27) Gyroscopes ω b ib = b g + (I + M g ) ω b ib + G g f b ib + w g (28) Basic Defiitios Liearizatio Iertial Measuremets Attitude Error Estimates of Sesor Measuremets Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

35 Error Modelig Example Accelerometers f b ib = b a + (I + M a ) f b ib + l a + w a (27) Biases Gyroscopes ω b ib = b g + (I + M g ) ω b ib + G g f b ib + w g (28) Basic Defiitios Liearizatio Iertial Measuremets Attitude Error Estimates of Sesor Measuremets Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

36 Error Modelig Example Accelerometers f b ib = b a + (I + M a ) f b ib + l a + w a (27) Misaligmet ad SF Errors Gyroscopes ω b ib = b g + (I + M g ) ω b ib + G g f b ib + w g (28) Basic Defiitios Liearizatio Iertial Measuremets Attitude Error Estimates of Sesor Measuremets Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

37 Error Modelig Example Accelerometers f b ib = b a + (I + M a ) f b ib + l a + w a (27) No-liearity Gyroscopes ω b ib = b g + (I + M g ) ω b ib + G g f b ib + w g (28) Basic Defiitios Liearizatio Iertial Measuremets Attitude Error Estimates of Sesor Measuremets Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

38 Error Modelig Example Accelerometers f b ib = b a + (I + M a ) f b ib + l a + w a (27) Gyroscopes G-Sesitivity ω b ib = b g + (I + M g ) ω b ib + G g f b ib + w g (28) Basic Defiitios Liearizatio Iertial Measuremets Attitude Error Estimates of Sesor Measuremets Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

39 Error Modelig Example Accelerometers f b ib = b a + (I + M a ) f b ib + l a + w a (27) Noise Gyroscopes ω b ib = b g + (I + M g ) ω b ib + G g f b ib + w g (28) Basic Defiitios Liearizatio Iertial Measuremets Attitude Error Estimates of Sesor Measuremets Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

40 Pos, Vel, Force ad Agular Rate Errors Positio error Velocity error Specific force errors δ r γ βb = r γ βb ˆ r γ βb (29) δ v γ βb = v γ βb ˆ v γ βb (30) δ f b ib = f b ib ˆ f b ib (31) ef b ib = f b ib ˆ f b ib = δf b ib (32) Agular rate errors δ ω b ib = ω b ib ˆ ω b ib (33) e ω b ib = ω b ib ˆ ω b ib = δ ω b ib (34) Basic Defiitios Liearizatio Iertial Measuremets Attitude Error Estimates of Sesor Measuremets Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

41 Attitude Error Defiitio Defie δc γ b = C γ b Ĉ γ b = e [δ ψ γ γb ] I + [δψ γ γb ] (35) This is the error i attitude resultig from errors i estimatig the agular rates. Basic Defiitios Liearizatio Iertial Measuremets Attitude Error Estimates of Sesor Measuremets Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

42 Attitude Error Properties The attitude error is a multiplicative small agle trasformatio from the actual frame to the computed frame Ĉ γ b = (I [δ ψ γ γb ])C γ b (36) Similarly, C γ b = (I + [δ ψ γ γb ])Ĉ γ b (37) Basic Defiitios Liearizatio Iertial Measuremets Attitude Error Estimates of Sesor Measuremets Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

43 Specific Force ad Agular Rates Similarly we ca attempt to estimate the specific force ad agular rate by applyig correctio based o our estimate of the error. ˆ f b ib = f b ib ˆ f b ib (38) ˆ ω b ib = ω b ib ˆ ω b ib (39) where ˆ f b ib ad ˆ ω b ib are the accelerometer ad gyroscope estimated calibratio values, respectively. Basic Defiitios Liearizatio Iertial Measuremets Attitude Error Estimates of Sesor Measuremets Aly El-Osery, Kevi Wedeward (NMT) EE 570: Locatio ad Navigatio April 4, / 23

EE 570: Location and Navigation

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

EE 570: Location and Navigation

EE 570: Location and Navigation 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

More information

EE 570: Location and Navigation

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

EE 570: Location and Navigation

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

EE 565: Position, Navigation and Timing

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

EE 570: Location and Navigation

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

EE 570: Location and Navigation

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

EE 565: Position, Navigation, and Timing

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

EE 570: Location and Navigation

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

Navigation Mathematics: Kinematics (Earth Surface & Gravity Models) EE 570: Location and Navigation

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

Navigation Mathematics: Kinematics (Coordinate Frame Transformation) EE 570: Location and Navigation

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

EE 570: Location and Navigation

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

EE 570: Location and Navigation

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

The Mathematical Model and the Simulation Modelling Algoritm of the Multitiered Mechanical System

The Mathematical Model and the Simulation Modelling Algoritm of the Multitiered Mechanical System The Mathematical Model ad the Simulatio Modellig Algoritm of the Multitiered Mechaical System Demi Aatoliy, Kovalev Iva Dept. of Optical Digital Systems ad Techologies, The St. Petersburg Natioal Research

More information

TERRAIN REFERENCED UAV NAVIGATION WITH LIDAR A COMPARISON OF SEQUENTAIL PROCESSING AND BATCH PROCESSING ALGORITHMS

TERRAIN REFERENCED UAV NAVIGATION WITH LIDAR A COMPARISON OF SEQUENTAIL PROCESSING AND BATCH PROCESSING ALGORITHMS 8 H INERNAIONAL CONGRESS OF HE AERONAUICAL SCIENCES ERRAIN REFERENCED UAV NAVIGAION WIH LIDAR A COPARISON OF SEQUENAIL PROCESSING AND BACH PROCESSING ALGORIHS Yeoha Hwag, i-jea ah KAIS (Korea Advaced Istitute

More information

Presentation of complex number in Cartesian and polar coordinate system

Presentation of complex number in Cartesian and polar coordinate system a + bi, aεr, bεr i = z = a + bi a = Re(z), b = Im(z) give z = a + bi & w = c + di, a + bi = c + di a = c & b = d The complex cojugate of z = a + bi is z = a bi The sum of complex cojugates is real: z +

More information

1. Linearization of a nonlinear system given in the form of a system of ordinary differential equations

1. Linearization of a nonlinear system given in the form of a system of ordinary differential equations . Liearizatio of a oliear system give i the form of a system of ordiary differetial equatios We ow show how to determie a liear model which approximates the behavior of a time-ivariat oliear system i a

More information

FAILURE CRITERIA: MOHR S CIRCLE AND PRINCIPAL STRESSES

FAILURE CRITERIA: MOHR S CIRCLE AND PRINCIPAL STRESSES LECTURE Third Editio FAILURE CRITERIA: MOHR S CIRCLE AND PRINCIPAL STRESSES A. J. Clark School of Egieerig Departmet of Civil ad Evirometal Egieerig Chapter 7.4 b Dr. Ibrahim A. Assakkaf SPRING 3 ENES

More information

Learning Theory: Lecture Notes

Learning Theory: Lecture Notes Learig Theory: Lecture Notes Kamalika Chaudhuri October 4, 0 Cocetratio of Averages Cocetratio of measure is very useful i showig bouds o the errors of machie-learig algorithms. We will begi with a basic

More information

Physics 324, Fall Dirac Notation. These notes were produced by David Kaplan for Phys. 324 in Autumn 2001.

Physics 324, Fall Dirac Notation. These notes were produced by David Kaplan for Phys. 324 in Autumn 2001. Physics 324, Fall 2002 Dirac Notatio These otes were produced by David Kapla for Phys. 324 i Autum 2001. 1 Vectors 1.1 Ier product Recall from liear algebra: we ca represet a vector V as a colum vector;

More information

Notes The Incremental Motion Model:

Notes The Incremental Motion Model: The Icremetal Motio Model: The Jacobia Matrix I the forward kiematics model, we saw that it was possible to relate joit agles θ, to the cofiguratio of the robot ed effector T I this sectio, we will see

More information

De Moivre s Theorem - ALL

De Moivre s Theorem - ALL De Moivre s Theorem - ALL. Let x ad y be real umbers, ad be oe of the complex solutios of the equatio =. Evaluate: (a) + + ; (b) ( x + y)( x + y). [6]. (a) Sice is a complex umber which satisfies = 0,.

More information

Chapter 4 : Laplace Transform

Chapter 4 : Laplace Transform 4. Itroductio Laplace trasform is a alterative to solve the differetial equatio by the complex frequecy domai ( s = σ + jω), istead of the usual time domai. The DE ca be easily trasformed ito a algebraic

More information

MAT 271 Project: Partial Fractions for certain rational functions

MAT 271 Project: Partial Fractions for certain rational functions MAT 7 Project: Partial Fractios for certai ratioal fuctios Prerequisite kowledge: partial fractios from MAT 7, a very good commad of factorig ad complex umbers from Precalculus. To complete this project,

More information

EE 570: Location and Navigation: Theory & Practice

EE 570: Location and Navigation: Theory & Practice EE 570: Locatio ad Naigatio: Thory & Practic Naigatio Ssors ad INS Mchaizatio NMT EE 570: Locatio ad Naigatio: Thory & Practic Slid 1 of 13 Naigatio Ssors ad INS Mchaizatio Naigatio Equatios Cas 3: Na

More information

MCT242: Electronic Instrumentation Lecture 2: Instrumentation Definitions

MCT242: Electronic Instrumentation Lecture 2: Instrumentation Definitions Faculty of Egieerig MCT242: Electroic Istrumetatio Lecture 2: Istrumetatio Defiitios Overview Measuremet Error Accuracy Precisio ad Mea Resolutio Mea Variace ad Stadard deviatio Fiesse Sesitivity Rage

More information

Physics 7440, Solutions to Problem Set # 8

Physics 7440, Solutions to Problem Set # 8 Physics 7440, Solutios to Problem Set # 8. Ashcroft & Mermi. For both parts of this problem, the costat offset of the eergy, ad also the locatio of the miimum at k 0, have o effect. Therefore we work with

More information

Signal Processing in Mechatronics. Lecture 3, Convolution, Fourier Series and Fourier Transform

Signal Processing in Mechatronics. Lecture 3, Convolution, Fourier Series and Fourier Transform Sigal Processig i Mechatroics Summer semester, 1 Lecture 3, Covolutio, Fourier Series ad Fourier rasform Dr. Zhu K.P. AIS, UM 1 1. Covolutio Covolutio Descriptio of LI Systems he mai premise is that the

More information

Solutions to Final Exam Review Problems

Solutions to Final Exam Review Problems . Let f(x) 4+x. Solutios to Fial Exam Review Problems Math 5C, Witer 2007 (a) Fid the Maclauri series for f(x), ad compute its radius of covergece. Solutio. f(x) 4( ( x/4)) ( x/4) ( ) 4 4 + x. Sice the

More information

Generalizing the DTFT. The z Transform. Complex Exponential Excitation. The Transfer Function. Systems Described by Difference Equations

Generalizing the DTFT. The z Transform. Complex Exponential Excitation. The Transfer Function. Systems Described by Difference Equations Geeraliig the DTFT The Trasform M. J. Roberts - All Rights Reserved. Edited by Dr. Robert Akl 1 The forward DTFT is defied by X e jω = x e jω i which = Ω is discrete-time radia frequecy, a real variable.

More information

f(w) w z =R z a 0 a n a nz n Liouville s theorem, we see that Q is constant, which implies that P is constant, which is a contradiction.

f(w) w z =R z a 0 a n a nz n Liouville s theorem, we see that Q is constant, which implies that P is constant, which is a contradiction. Theorem 3.6.4. [Liouville s Theorem] Every bouded etire fuctio is costat. Proof. Let f be a etire fuctio. Suppose that there is M R such that M for ay z C. The for ay z C ad R > 0 f (z) f(w) 2πi (w z)

More information

5.74 TIME-DEPENDENT QUANTUM MECHANICS

5.74 TIME-DEPENDENT QUANTUM MECHANICS p. 1 5.74 TIME-DEPENDENT QUANTUM MECHANICS The time evolutio of the state of a system is described by the time-depedet Schrödiger equatio (TDSE): i t ψ( r, t)= H ˆ ψ( r, t) Most of what you have previously

More information

Classical Mechanics Qualifying Exam Solutions Problem 1.

Classical Mechanics Qualifying Exam Solutions Problem 1. Jauary 4, Uiversity of Illiois at Chicago Departmet of Physics Classical Mechaics Qualifyig Exam Solutios Prolem. A cylider of a o-uiform radial desity with mass M, legth l ad radius R rolls without slippig

More information

6.003 Homework #3 Solutions

6.003 Homework #3 Solutions 6.00 Homework # Solutios Problems. Complex umbers a. Evaluate the real ad imagiary parts of j j. π/ Real part = Imagiary part = 0 e Euler s formula says that j = e jπ/, so jπ/ j π/ j j = e = e. Thus the

More information

HE ATOM & APPROXIMATION METHODS MORE GENERAL VARIATIONAL TREATMENT. Examples:

HE ATOM & APPROXIMATION METHODS MORE GENERAL VARIATIONAL TREATMENT. Examples: 5.6 4 Lecture #3-4 page HE ATOM & APPROXIMATION METHODS MORE GENERAL VARIATIONAL TREATMENT Do t restrict the wavefuctio to a sigle term! Could be a liear combiatio of several wavefuctios e.g. two terms:

More information

ECE-S352 Introduction to Digital Signal Processing Lecture 3A Direct Solution of Difference Equations

ECE-S352 Introduction to Digital Signal Processing Lecture 3A Direct Solution of Difference Equations ECE-S352 Itroductio to Digital Sigal Processig Lecture 3A Direct Solutio of Differece Equatios Discrete Time Systems Described by Differece Equatios Uit impulse (sample) respose h() of a DT system allows

More information

Matrix Algebra 2.2 THE INVERSE OF A MATRIX Pearson Education, Inc.

Matrix Algebra 2.2 THE INVERSE OF A MATRIX Pearson Education, Inc. 2 Matrix Algebra 2.2 THE INVERSE OF A MATRIX MATRIX OPERATIONS A matrix A is said to be ivertible if there is a matrix C such that CA = I ad AC = I where, the idetity matrix. I = I I this case, C is a

More information

x !1! + 1!2!

x !1! + 1!2! 4 Euler-Maclauri Suatio Forula 4. Beroulli Nuber & Beroulli Polyoial 4.. Defiitio of Beroulli Nuber Beroulli ubers B (,,3,) are defied as coefficiets of the followig equatio. x e x - B x! 4.. Expreesio

More information

Modified Decomposition Method by Adomian and. Rach for Solving Nonlinear Volterra Integro- Differential Equations

Modified Decomposition Method by Adomian and. Rach for Solving Nonlinear Volterra Integro- Differential Equations Noliear Aalysis ad Differetial Equatios, Vol. 5, 27, o. 4, 57-7 HIKARI Ltd, www.m-hikari.com https://doi.org/.2988/ade.27.62 Modified Decompositio Method by Adomia ad Rach for Solvig Noliear Volterra Itegro-

More information

Chapter 10 Advanced Topics in Random Processes

Chapter 10 Advanced Topics in Random Processes ery Stark ad Joh W. Woods, Probability, Statistics, ad Radom Variables for Egieers, 4th ed., Pearso Educatio Ic.,. ISBN 978--3-33-6 Chapter Advaced opics i Radom Processes Sectios. Mea-Square (m.s.) Calculus

More information

Signal Processing in Mechatronics

Signal Processing in Mechatronics Sigal Processig i Mechatroics Zhu K.P. AIS, UM. Lecture, Brief itroductio to Sigals ad Systems, Review of Liear Algebra ad Sigal Processig Related Mathematics . Brief Itroductio to Sigals What is sigal

More information

System Diagnostics using Kalman Filter Estimation Error

System Diagnostics using Kalman Filter Estimation Error System Diagostics usig Kalma Filter Estimatio Error Prof. Seugchul Lee, Seugtae Park, Heechag Kim, Hyusuk Huh ICMR2015 (11/25/2015) Machie Health Diagostics Desig Desig DNA + + Blue Prit phagocytes lymphocytes

More information

Guiding-center transformation 1. δf = q ( c v f 0, (1) mc (v B 0) v f 0. (3)

Guiding-center transformation 1. δf = q ( c v f 0, (1) mc (v B 0) v f 0. (3) Guidig-ceter trasformatio 1 This is my otes whe readig Liu Che s book[1]. 1 Vlasov equatio The liearized Vlasov equatio is [ t + v x + q m E + v B c ] δf = q v m δe + v δb c v f, 1 where f ad δf are the

More information

5.80 Small-Molecule Spectroscopy and Dynamics

5.80 Small-Molecule Spectroscopy and Dynamics MIT OpeCourseWare http://ocw.mit.edu 5.80 Small-Molecule Spectroscopy ad Dyamics Fall 2008 For iformatio about citig these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Lecture # 33 Supplemet

More information

Quantum Theory Assignment 3

Quantum Theory Assignment 3 Quatum Theory Assigmet 3 Assigmet 3.1 1. Cosider a spi-1/ system i a magetic field i the z-directio. The Hamiltoia is give by: ) eb H = S z = ωs z. mc a) Fid the Heiseberg operators S x t), S y t), ad

More information

Chapter 7: The z-transform. Chih-Wei Liu

Chapter 7: The z-transform. Chih-Wei Liu Chapter 7: The -Trasform Chih-Wei Liu Outlie Itroductio The -Trasform Properties of the Regio of Covergece Properties of the -Trasform Iversio of the -Trasform The Trasfer Fuctio Causality ad Stability

More information

Affine Structure from Motion

Affine Structure from Motion Affie Structure from Motio EECS 598-8 Fall 24! Foudatios of Computer Visio!! Istructor: Jaso Corso (jjcorso)! web.eecs.umich.edu/~jjcorso/t/598f4!! Readigs: FP 8.2! Date: /5/4!! Materials o these slides

More information

Building an NMR Quantum Computer

Building an NMR Quantum Computer Buildig a NMR Quatum Computer Spi, the Ster-Gerlach Experimet, ad the Bloch Sphere Kevi Youg Berkeley Ceter for Quatum Iformatio ad Computatio, Uiversity of Califoria, Berkeley, CA 9470 Scalable ad Secure

More information

PHY4905: Nearly-Free Electron Model (NFE)

PHY4905: Nearly-Free Electron Model (NFE) PHY4905: Nearly-Free Electro Model (NFE) D. L. Maslov Departmet of Physics, Uiversity of Florida (Dated: Jauary 12, 2011) 1 I. REMINDER: QUANTUM MECHANICAL PERTURBATION THEORY A. No-degeerate eigestates

More information

The time evolution of the state of a quantum system is described by the time-dependent Schrödinger equation (TDSE): ( ) ( ) 2m "2 + V ( r,t) (1.

The time evolution of the state of a quantum system is described by the time-dependent Schrödinger equation (TDSE): ( ) ( ) 2m 2 + V ( r,t) (1. Adrei Tokmakoff, MIT Departmet of Chemistry, 2/13/2007 1-1 574 TIME-DEPENDENT QUANTUM MECHANICS 1 INTRODUCTION 11 Time-evolutio for time-idepedet Hamiltoias The time evolutio of the state of a quatum system

More information

Course 4: Preparation for Calculus Unit 1: Families of Functions

Course 4: Preparation for Calculus Unit 1: Families of Functions Course 4: Preparatio for Calculus Uit 1: Families of Fuctios Review ad exted properties of basic fuctio families ad their uses i mathematical modelig Develop strategies for fidig rules of fuctios whose

More information

MTH112 Trigonometry 2 2 2, 2. 5π 6. cscθ = 1 sinθ = r y. secθ = 1 cosθ = r x. cotθ = 1 tanθ = cosθ. central angle time. = θ t.

MTH112 Trigonometry 2 2 2, 2. 5π 6. cscθ = 1 sinθ = r y. secθ = 1 cosθ = r x. cotθ = 1 tanθ = cosθ. central angle time. = θ t. MTH Trigoometry,, 5, 50 5 0 y 90 0, 5 0,, 80 0 0 0 (, 0) x, 7, 0 5 5 0, 00 5 5 0 7,,, Defiitios: siθ = opp. hyp. = y r cosθ = adj. hyp. = x r taθ = opp. adj. = siθ cosθ = y x cscθ = siθ = r y secθ = cosθ

More information

SOLUTION SET VI FOR FALL [(n + 2)(n + 1)a n+2 a n 1 ]x n = 0,

SOLUTION SET VI FOR FALL [(n + 2)(n + 1)a n+2 a n 1 ]x n = 0, 4. Series Solutios of Differetial Equatios:Special Fuctios 4.. Illustrative examples.. 5. Obtai the geeral solutio of each of the followig differetial equatios i terms of Maclauri series: d y (a dx = xy,

More information

Acoustic Field inside a Rigid Cylinder with a Point Source

Acoustic Field inside a Rigid Cylinder with a Point Source Acoustic Field iside a Rigid Cylider with a Poit Source 1 Itroductio The ai objectives of this Deo Model are to Deostrate the ability of Coustyx to odel a rigid cylider with a poit source usig Coustyx

More information

Ch3 Discrete Time Fourier Transform

Ch3 Discrete Time Fourier Transform Ch3 Discrete Time Fourier Trasform 3. Show that the DTFT of [] is give by ( k). e k 3. Determie the DTFT of the two sided sigal y [ ],. 3.3 Determie the DTFT of the causal sequece x[ ] A cos( 0 ) [ ],

More information

Lecture 25 (Dec. 6, 2017)

Lecture 25 (Dec. 6, 2017) Lecture 5 8.31 Quatum Theory I, Fall 017 106 Lecture 5 (Dec. 6, 017) 5.1 Degeerate Perturbatio Theory Previously, whe discussig perturbatio theory, we restricted ourselves to the case where the uperturbed

More information

On the Derivation and Implementation of a Four Stage Harmonic Explicit Runge-Kutta Method *

On the Derivation and Implementation of a Four Stage Harmonic Explicit Runge-Kutta Method * Applied Mathematics, 05, 6, 694-699 Published Olie April 05 i SciRes. http://www.scirp.org/joural/am http://dx.doi.org/0.46/am.05.64064 O the Derivatio ad Implemetatio of a Four Stage Harmoic Explicit

More information

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

ECE Spring Prof. David R. Jackson ECE Dept. Notes 20 ECE 6341 Sprig 016 Prof. David R. Jackso ECE Dept. Notes 0 1 Spherical Wave Fuctios Cosider solvig ψ + k ψ = 0 i spherical coordiates z φ θ r y x Spherical Wave Fuctios (cot.) I spherical coordiates we

More information

ME NUMERICAL METHODS Fall 2007

ME NUMERICAL METHODS Fall 2007 ME - 310 NUMERICAL METHODS Fall 2007 Group 02 Istructor: Prof. Dr. Eres Söylemez (Rm C205, email:eres@metu.edu.tr ) Class Hours ad Room: Moday 13:40-15:30 Rm: B101 Wedesday 12:40-13:30 Rm: B103 Course

More information

Quiz. Use either the RATIO or ROOT TEST to determine whether the series is convergent or not.

Quiz. Use either the RATIO or ROOT TEST to determine whether the series is convergent or not. Quiz. Use either the RATIO or ROOT TEST to determie whether the series is coverget or ot. e .6 POWER SERIES Defiitio. A power series i about is a series of the form c 0 c a c a... c a... a 0 c a where

More information

Stopping oscillations of a simple harmonic oscillator using an impulse force

Stopping oscillations of a simple harmonic oscillator using an impulse force It. J. Adv. Appl. Math. ad Mech. 5() (207) 6 (ISSN: 2347-2529) IJAAMM Joural homepage: www.ijaamm.com Iteratioal Joural of Advaces i Applied Mathematics ad Mechaics Stoppig oscillatios of a simple harmoic

More information

R is a scalar defined as follows:

R is a scalar defined as follows: Math 8. Notes o Dot Product, Cross Product, Plaes, Area, ad Volumes This lecture focuses primarily o the dot product ad its may applicatios, especially i the measuremet of agles ad scalar projectio ad

More information

Measurement uncertainty of the sound absorption

Measurement uncertainty of the sound absorption Measuremet ucertaity of the soud absorptio coefficiet Aa Izewska Buildig Research Istitute, Filtrowa Str., 00-6 Warsaw, Polad a.izewska@itb.pl 6887 The stadard ISO/IEC 705:005 o the competece of testig

More information

EECS 442 Computer vision. Multiple view geometry Affine structure from Motion

EECS 442 Computer vision. Multiple view geometry Affine structure from Motion EECS 442 Computer visio Multiple view geometry Affie structure from Motio - Affie structure from motio problem - Algebraic methods - Factorizatio methods Readig: [HZ] Chapters: 6,4,8 [FP] Chapter: 2 Some

More information

Signal Processing. Lecture 02: Discrete Time Signals and Systems. Ahmet Taha Koru, Ph. D. Yildiz Technical University.

Signal Processing. Lecture 02: Discrete Time Signals and Systems. Ahmet Taha Koru, Ph. D. Yildiz Technical University. Sigal Processig Lecture 02: Discrete Time Sigals ad Systems Ahmet Taha Koru, Ph. D. Yildiz Techical Uiversity 2017-2018 Fall ATK (YTU) Sigal Processig 2017-2018 Fall 1 / 51 Discrete Time Sigals Discrete

More information

APPENDIX F Complex Numbers

APPENDIX F Complex Numbers APPENDIX F Complex Numbers Operatios with Complex Numbers Complex Solutios of Quadratic Equatios Polar Form of a Complex Number Powers ad Roots of Complex Numbers Operatios with Complex Numbers Some equatios

More information

j=1 dz Res(f, z j ) = 1 d k 1 dz k 1 (z c)k f(z) Res(f, c) = lim z c (k 1)! Res g, c = f(c) g (c)

j=1 dz Res(f, z j ) = 1 d k 1 dz k 1 (z c)k f(z) Res(f, c) = lim z c (k 1)! Res g, c = f(c) g (c) Problem. Compute the itegrals C r d for Z, where C r = ad r >. Recall that C r has the couter-clockwise orietatio. Solutio: We will use the idue Theorem to solve this oe. We could istead use other (perhaps

More information

BASIC STUDY ON THE REAL-TIME THREE-DIMENSIONAL POSITION MEASUREMENT OF GROUND MOVING OBJECTS USING AN ACCELEROMETER AND GYROS

BASIC STUDY ON THE REAL-TIME THREE-DIMENSIONAL POSITION MEASUREMENT OF GROUND MOVING OBJECTS USING AN ACCELEROMETER AND GYROS BASIC STUDY ON THE REAL-TIME THREE-DIMENSIONAL POSITION MEASUREMENT OF GROUND MOVING OBJECTS USING AN ACCELEROMETER AND GYROS Yoshiwo OKAMOTO*, Toshio KOIZUMI**, Yasuyuki SIRAI*** *Dep. Electrical, Electroics

More information

Optimal Estimator for a Sample Set with Response Error. Ed Stanek

Optimal Estimator for a Sample Set with Response Error. Ed Stanek Optial Estiator for a Saple Set wit Respose Error Ed Staek Itroductio We develop a optial estiator siilar to te FP estiator wit respose error tat was cosidered i c08ed63doc Te first 6 pages of tis docuet

More information

Solution of Differential Equation from the Transform Technique

Solution of Differential Equation from the Transform Technique Iteratioal Joural of Computatioal Sciece ad Mathematics ISSN 0974-3189 Volume 3, Number 1 (2011), pp 121-125 Iteratioal Research Publicatio House http://wwwirphousecom Solutio of Differetial Equatio from

More information

Modelling and Simulation of Marine Craft KINEMATICS

Modelling and Simulation of Marine Craft KINEMATICS Modellig ad Simulatio of Marie Craft KINEMATICS Edi Omerdic Seior Research Fellow Outlie Kiematics Referece Frames ECEF Frame NED Frame BODY Frame UNITY Frame UTM Trasformatios BODY NED Euler Agles Quaterios

More information

Σk=1. d n () z = () z Σ dz n. d n. 1 = () z. z (0.1 - ) d f 2 = dz 0 () z = 1 () z, log. f n = n-2 () z = n-1 () z. e log z. = e.

Σk=1. d n () z = () z Σ dz n. d n. 1 = () z. z (0.1 - ) d f 2 = dz 0 () z = 1 () z, log. f n = n-2 () z = n-1 () z. e log z. = e. Series Expasio of Gamma Fuctio & the Reciprocal. Taylor Expasio aroud a Higher Derivative of Gamma Fuctio The formula of the higher derivative of the gamma fuctio & the reciprocal was discovered by Masayuki

More information

Office: JILA A709; Phone ;

Office: JILA A709; Phone ; Office: JILA A709; Phoe 303-49-7841; email: weberjm@jila.colorado.edu Problem Set 5 To be retured before the ed of class o Wedesday, September 3, 015 (give to me i perso or slide uder office door). 1.

More information

TIME-CORRELATION FUNCTIONS

TIME-CORRELATION FUNCTIONS p. 8 TIME-CORRELATION FUNCTIONS Time-correlatio fuctios are a effective way of represetig the dyamics of a system. They provide a statistical descriptio of the time-evolutio of a variable for a esemble

More information

THE KALMAN FILTER RAUL ROJAS

THE KALMAN FILTER RAUL ROJAS THE KALMAN FILTER RAUL ROJAS Abstract. This paper provides a getle itroductio to the Kalma filter, a umerical method that ca be used for sesor fusio or for calculatio of trajectories. First, we cosider

More information

The z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j

The z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j The -Trasform 7. Itroductio Geeralie the complex siusoidal represetatio offered by DTFT to a represetatio of complex expoetial sigals. Obtai more geeral characteristics for discrete-time LTI systems. 7.

More information

Notes 18 Green s Functions

Notes 18 Green s Functions ECE 638 Fall 017 David R. Jackso Notes 18 Gree s Fuctios Notes are from D. R. Wilto, Dept. of ECE 1 Gree s Fuctios The Gree's fuctio method is a powerful ad systematic method for determiig a solutio to

More information

Assignment 2 Solutions SOLUTION. ϕ 1 Â = 3 ϕ 1 4i ϕ 2. The other case can be dealt with in a similar way. { ϕ 2 Â} χ = { 4i ϕ 1 3 ϕ 2 } χ.

Assignment 2 Solutions SOLUTION. ϕ 1  = 3 ϕ 1 4i ϕ 2. The other case can be dealt with in a similar way. { ϕ 2 Â} χ = { 4i ϕ 1 3 ϕ 2 } χ. PHYSICS 34 QUANTUM PHYSICS II (25) Assigmet 2 Solutios 1. With respect to a pair of orthoormal vectors ϕ 1 ad ϕ 2 that spa the Hilbert space H of a certai system, the operator  is defied by its actio

More information

Let A(x) and B(x) be two polynomials of degree n 1:

Let A(x) and B(x) be two polynomials of degree n 1: MI-EVY (2011/2012) J. Holub: 4. DFT, FFT ad Patter Matchig p. 2/42 Operatios o polyomials MI-EVY (2011/2012) J. Holub: 4. DFT, FFT ad Patter Matchig p. 4/42 Efficiet Patter Matchig (MI-EVY) 4. DFT, FFT

More information

Preliminary Examination - Day 1 Thursday, May 12, 2016

Preliminary Examination - Day 1 Thursday, May 12, 2016 UNL - Departmet of Physics ad Astroomy Prelimiary Examiatio - Day Thursday, May, 6 This test covers the topics of Quatum Mechaics (Topic ) ad Electrodyamics (Topic ). Each topic has 4 A questios ad 4 B

More information

Proof of Fermat s Last Theorem by Algebra Identities and Linear Algebra

Proof of Fermat s Last Theorem by Algebra Identities and Linear Algebra Proof of Fermat s Last Theorem by Algebra Idetities ad Liear Algebra Javad Babaee Ragai Youg Researchers ad Elite Club, Qaemshahr Brach, Islamic Azad Uiversity, Qaemshahr, Ira Departmet of Civil Egieerig,

More information

Appendix F: Complex Numbers

Appendix F: Complex Numbers Appedix F Complex Numbers F1 Appedix F: Complex Numbers Use the imagiary uit i to write complex umbers, ad to add, subtract, ad multiply complex umbers. Fid complex solutios of quadratic equatios. Write

More information

CHAPTER 10 INFINITE SEQUENCES AND SERIES

CHAPTER 10 INFINITE SEQUENCES AND SERIES CHAPTER 10 INFINITE SEQUENCES AND SERIES 10.1 Sequeces 10.2 Ifiite Series 10.3 The Itegral Tests 10.4 Compariso Tests 10.5 The Ratio ad Root Tests 10.6 Alteratig Series: Absolute ad Coditioal Covergece

More information

TERMINAL GUIDANCE USING A DOPPLER BEAM SHARPENING RADAR Jeremy A. Hodgson *, BSc(Hons), MSc, MIEE The MathWorks Ltd, Cambridge, UK

TERMINAL GUIDANCE USING A DOPPLER BEAM SHARPENING RADAR Jeremy A. Hodgson *, BSc(Hons), MSc, MIEE The MathWorks Ltd, Cambridge, UK TERMINAL GUIDANCE USING A DOPPLER BEAM SHARPENING RADAR Jeremy A. Hodgso *, BSc(Hos, MSc, MIEE The MathWorks Ltd, Cambridge, UK David W. Lee, BSc(Hos,MSc QietiQ Limited, Farborough, UK Abstract This paper

More information

Q-BINOMIALS AND THE GREATEST COMMON DIVISOR. Keith R. Slavin 8474 SW Chevy Place, Beaverton, Oregon 97008, USA.

Q-BINOMIALS AND THE GREATEST COMMON DIVISOR. Keith R. Slavin 8474 SW Chevy Place, Beaverton, Oregon 97008, USA. INTEGERS: ELECTRONIC JOURNAL OF COMBINATORIAL NUMBER THEORY 8 2008, #A05 Q-BINOMIALS AND THE GREATEST COMMON DIVISOR Keith R. Slavi 8474 SW Chevy Place, Beaverto, Orego 97008, USA slavi@dsl-oly.et Received:

More information

1. Hydrogen Atom: 3p State

1. Hydrogen Atom: 3p State 7633A QUANTUM MECHANICS I - solutio set - autum. Hydroge Atom: 3p State Let us assume that a hydroge atom is i a 3p state. Show that the radial part of its wave fuctio is r u 3(r) = 4 8 6 e r 3 r(6 r).

More information

It should be unbiased, or approximately unbiased. Variance of the variance estimator should be small. That is, the variance estimator is stable.

It should be unbiased, or approximately unbiased. Variance of the variance estimator should be small. That is, the variance estimator is stable. Chapter 10 Variace Estimatio 10.1 Itroductio Variace estimatio is a importat practical problem i survey samplig. Variace estimates are used i two purposes. Oe is the aalytic purpose such as costructig

More information

Created by T. Madas SERIES. Created by T. Madas

Created by T. Madas SERIES. Created by T. Madas SERIES SUMMATIONS BY STANDARD RESULTS Questio (**) Use stadard results o summatios to fid the value of 48 ( r )( 3r ). 36 FP-B, 66638 Questio (**+) Fid, i fully simplified factorized form, a expressio

More information

Sequences, Mathematical Induction, and Recursion. CSE 2353 Discrete Computational Structures Spring 2018

Sequences, Mathematical Induction, and Recursion. CSE 2353 Discrete Computational Structures Spring 2018 CSE 353 Discrete Computatioal Structures Sprig 08 Sequeces, Mathematical Iductio, ad Recursio (Chapter 5, Epp) Note: some course slides adopted from publisher-provided material Overview May mathematical

More information

Microscopic Theory of Transport (Fall 2003) Lecture 6 (9/19/03) Static and Short Time Properties of Time Correlation Functions

Microscopic Theory of Transport (Fall 2003) Lecture 6 (9/19/03) Static and Short Time Properties of Time Correlation Functions .03 Microscopic Theory of Trasport (Fall 003) Lecture 6 (9/9/03) Static ad Short Time Properties of Time Correlatio Fuctios Refereces -- Boo ad Yip, Chap There are a umber of properties of time correlatio

More information

Fall 2011, EE123 Digital Signal Processing

Fall 2011, EE123 Digital Signal Processing Lecture 5 Miki Lustig, UCB September 14, 211 Miki Lustig, UCB Motivatios for Discrete Fourier Trasform Sampled represetatio i time ad frequecy umerical Fourier aalysis requires a Fourier represetatio that

More information

The Gamma function. Marco Bonvini. October 9, dt e t t z 1. (1) Γ(z + 1) = z Γ(z) : (2) = e t t z. + z dt e t t z 1. = z Γ(z).

The Gamma function. Marco Bonvini. October 9, dt e t t z 1. (1) Γ(z + 1) = z Γ(z) : (2) = e t t z. + z dt e t t z 1. = z Γ(z). The Gamma fuctio Marco Bovii October 9, 2 Gamma fuctio The Euler Gamma fuctio is defied as Γ() It is easy to show that Γ() satisfy the recursio relatio ideed, itegratig by parts, dt e t t. () Γ( + ) Γ()

More information

Problem 1. Problem Engineering Dynamics Problem Set 9--Solution. Find the equation of motion for the system shown with respect to:

Problem 1. Problem Engineering Dynamics Problem Set 9--Solution. Find the equation of motion for the system shown with respect to: 2.003 Egieerig Dyamics Problem Set 9--Solutio Problem 1 Fid the equatio of motio for the system show with respect to: a) Zero sprig force positio. Draw the appropriate free body diagram. b) Static equilibrium

More information

MH1101 AY1617 Sem 2. Question 1. NOT TESTED THIS TIME

MH1101 AY1617 Sem 2. Question 1. NOT TESTED THIS TIME MH AY67 Sem Questio. NOT TESTED THIS TIME ( marks Let R be the regio bouded by the curve y 4x x 3 ad the x axis i the first quadrat (see figure below. Usig the cylidrical shell method, fid the volume of

More information

Data Assimilation. Alan O Neill University of Reading, UK

Data Assimilation. Alan O Neill University of Reading, UK Data Assimilatio Ala O Neill Uiversity of Readig, UK he Kalma Filter Kalma Filter (expesive Use model equatios to propagate B forward i time. B B(t Aalysis step as i OI Evolutio of Covariace Matrices (

More information

Appendix K. The three-point correlation function (bispectrum) of density peaks

Appendix K. The three-point correlation function (bispectrum) of density peaks Appedix K The three-poit correlatio fuctio (bispectrum) of desity peaks Cosider the smoothed desity field, ρ (x) ρ [ δ (x)], with a geeral smoothig kerel W (x) δ (x) d yw (x y)δ(y). (K.) We defie the peaks

More information

CALCULUS BASIC SUMMER REVIEW

CALCULUS BASIC SUMMER REVIEW CALCULUS BASIC SUMMER REVIEW NAME rise y y y Slope of a o vertical lie: m ru Poit Slope Equatio: y y m( ) The slope is m ad a poit o your lie is, ). ( y Slope-Itercept Equatio: y m b slope= m y-itercept=

More information

Notes 8 Singularities

Notes 8 Singularities ECE 6382 Fall 27 David R. Jackso Notes 8 Sigularities Notes are from D. R. Wilto, Dept. of ECE Sigularity A poit s is a sigularity of the fuctio f () if the fuctio is ot aalytic at s. (The fuctio does

More information