EE 570: Location and Navigation
|
|
- Carmella Casey
- 5 years ago
- Views:
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 Error Mechanization (ECEF) Aly El-Osery Electrical Engineering Department, New Mexico Tech Socorro, New Mexico, USA April 11, 2013 Attitude Velocity Gravity Position Summary
More informationEE 570: Location and Navigation
EE 570: Location and Navigation 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 informationEE 570: Location and Navigation
EE 570: Location and Navigation Error Mechanization (ECEF) Aly El-Osery 1 Stephen Bruder 2 1 Electrical Engineering Department, New Mexico Tech Socorro, New Mexico, USA 2 Electrical and Computer Engineering
More informationEE 570: Location and Navigation
EE 570: Location and Navigation 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 informationEE 570: Location and Navigation
EE 570: Location and Navigation Navigation Equations: Nav Mechanization Aly El-Osery Kevin Wedeward Electrical Engineering Department, New Mexico Tech Socorro, New Mexico, USA In Collaboration with Stephen
More informationEE 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 informationLecture. Aided INS EE 570: Location and Navigation. 1 Overview. 1.1 ECEF as and Example. 1.2 Inertial Measurements
Lecture Aided EE 570: Location and Navigation Lecture Notes Update on April 13, 2016 Aly El-Osery and Kevin Wedeward, Electrical Engineering Dept., New Mexico Tech In collaoration with Stephen Bruder,
More informationEE 570: Location and Navigation
EE 570: 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 informationEE 570: Location and Navigation
EE 570: Location and Navigation Navigation Mathematics: Coordinate Frames Kevin Wedeward Aly El-Osery Electrical Engineering Department, New Mexico Tech Socorro, New Mexico, USA In Collaboration with Stephen
More informationEE 565: Position, Navigation, and Timing
EE 565: Position, Navigation, and Timing Kalman Filtering Example Aly El-Osery Kevin Wedeward Electrical Engineering Department, New Mexico Tech Socorro, New Mexico, USA In Collaboration with Stephen Bruder
More informationEE 570: Location and Navigation
EE 570: Location and Navigation Sensor Technology Aly El-Osery Kevin Wedeward Electrical Engineering Department, New Mexico Tech Socorro, New Mexico, USA In Collaboration with Stephen Bruder Electrical
More informationNavigation Mathematics: Kinematics (Earth Surface & Gravity Models) EE 570: Location and Navigation
Lecture Navigation Mathematics: Kinematics (Earth Surface & ) EE 570: Location and Navigation Lecture Notes Update on March 10, 2016 Aly El-Osery and Kevin Wedeward, Electrical Engineering Dept., New Mexico
More informationNavigation Mathematics: Kinematics (Coordinate Frame Transformation) EE 570: Location and Navigation
Lecture Navigation Mathematics: Kinematics (Coordinate Frame Transformation) EE 570: Location and Navigation Lecture Notes Update on Feruary 16, 2016 Aly El-Osery and Kevin Wedeward, Electrical Engineering
More informationEE 570: Location and Navigation
EE 570: Location and Navigation Sensor Technology Stephen Bruder 1 Aly El-Osery 2 1 Electrical and Computer Engineering Department, Embry-Riddle Aeronautical Univesity Prescott, Arizona, USA 2 Electrical
More informationEE 570: Location and Navigation
EE 570: Location and Navigation 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 informationThe 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 informationTERRAIN 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 informationPresentation 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 information1. 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 informationFAILURE 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 informationLearning 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 informationPhysics 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 informationNotes 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 informationDe 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 informationChapter 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 informationMAT 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 informationEE 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 informationMCT242: 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 informationPhysics 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 informationSignal 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 informationSolutions 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 informationGeneralizing 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 informationf(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 information5.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 informationClassical 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 information6.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 informationHE 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 informationECE-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 informationMatrix 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 informationx !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 informationModified 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 informationChapter 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 informationSignal 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 informationSystem 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 informationGuiding-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 information5.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 informationQuantum 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 informationChapter 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 informationAffine 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 informationBuilding 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 informationPHY4905: 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 informationThe 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 informationCourse 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 informationMTH112 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 informationSOLUTION 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 informationAcoustic 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 informationCh3 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 informationLecture 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 informationOn 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 informationECE 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 informationME 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 informationQuiz. 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 informationStopping 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 informationR 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 informationMeasurement 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 informationEECS 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 informationSignal 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 informationAPPENDIX 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 informationj=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 informationBASIC 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 informationOptimal 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 informationSolution 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 informationModelling 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.
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 informationOffice: 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 informationTIME-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 informationTHE 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 informationThe 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 informationNotes 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 informationAssignment 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 informationLet 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 informationPreliminary 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 informationProof 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 informationAppendix 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 informationCHAPTER 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 informationTERMINAL 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 informationQ-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 information1. 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 informationIt 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 informationCreated 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 informationSequences, 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 informationMicroscopic 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 informationFall 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 informationThe 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 informationProblem 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 informationMH1101 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 informationData 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 informationAppendix 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 informationCALCULUS 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 informationNotes 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