Mobility-Pattern Based Localization Update Algorithms for Mobile Wireless Sensor Networks

Size: px
Start display at page:

Download "Mobility-Pattern Based Localization Update Algorithms for Mobile Wireless Sensor Networks"

Transcription

1 Mobiliy-Paern Based Localizaion Updae Algorihms for Mobile Wireless Sensor Neworks Mohammad Y. Al-laho, Min Song, Jun Wang Deparmen of Elecrical and Compuer Engineering Old Dominion Universiy 231 Kaufman Hall Norfolk, VA 23529, USA {malla3, msong, Absrac. In mobile wireless sensor neworks, sensors move in he moniored area a any direcion and a any speed. Unlike many oher neworking hoss, sensor nodes do no have global addresses. Very ofen hey are idenified by using a locaion-based addressing scheme. Therefore, i is imporan o have he knowledge of he sensor locaion indicaing where he daa came from. In his paper, we design hree mobiliy-paern based localizaion updae algorihms. Specifically, we divide sensor movemens ino hree saes, Pause, Linear, and Random. Each sae adops differen localizaion updae algorihms. Analyical and simulaion resuls are provided o sudy he localizaion cos and locaion accuracy of he proposed localizaion-updae algorihm in differen mobiliy paerns. The analysis o hese resuls indicaes ha he localizaion cos is minimized and he locaion accuracy is improved. 1 Inroducion In mobile wireless sensor neworks, sensors move in he moniored area a any direcion and a any speed. Unlike many oher neworking hoss, sensor nodes do no have global addresses. Very ofen hey are idenified by using a locaion-based addressing scheme [5]. The rouing process also requires locaion-based naming where he users are more ineresed in querying a locaion of a phenomenon, raher han querying an individual node [5]. Therefore, i is imporan o have he knowledge of he sensor locaion indicaing where he dae came from. Many researches have been conduced o solve he localizaion problem of moving objecs. The mos famous one is he GPS sysem. Unforunaely, GPS is no a good choice for mobile wireless sensor neworks because of is significan power consumpion of he GPS receivers [7]. Inuiively, he localizaion of mobile wireless sensors can be performed by using predicion echniques o anicipae he possible nex sep of a sensor node based on he movemen paern model and he hisory of he movemen paern [4]. Unforunaely, i has been found ha predicion echniques are no always accurae [1]. Especially, when he movemen paern is no predicive (no linear he accuracy decreases significanly. A more reliable mehod is o inegrae he predicion echniques and localizaion updae algorihms. In his case, he localizaion updae algorihms can correc any inaccurae anicipaion. This mehod, however, is raher

2 2 Mohammad Y. Al-laho, Min Song, Jun Wang oo complex and needs many calculaions and filering echniques o refine he possible sampling daa colleced from previous movemens. The energy consumpion and he memory size resricions of a sensor node limi he use of hese predicion mehods. Anoher problem for mobile localizaion schemes in lieraure is ha hey impracically simplify he mobiliy paern and hus hey do no work well when nodes move randomly. The more he mobiliy paern deviaes from he linear paern, he more he locaion accuracy ges inaccurae and he updae cos ges higher. I is necessary o differeniae differen mobiliy paerns in he sysem. For each paern, a differen localizaion-updae scheme is used o ge he mos accurae resuls. I is worh o noe ha one of he criical design issues in mobile sensor localizaion is simpliciy, since wireless sensor nodes have limied power consumpion and compuaional capabiliy due o he small size of a node. In his paper, we design hree mobiliy-paern based localizaion updae algorihms. More specifically, we divide sensor movemens ino hree saes, Pause, Linear, and Random. Each sae adops differen localizaion updae algorihms. The objecives of his design are o 1 manage he localizaion updae period in opimal raes so ha he updae cos is minimized, and 2 keep he locaion error low. The res of he paper is organized as follows. Secion 2 presens he relaed work in localizaion updae algorihm for mobile wireless sensor neworks. In Secion 3, we inroduce he sysem models and our mobiliy-paern based localizaion updae algorihms. The numerical resuls are provided in Secion 4 o sudy he updae cos and locaion accuracy of he proposed localizaion updae algorihms in differen mobiliy paerns. Secion 5 concludes he paper. 2 Relaed Work Locaion updae schemes generally can be classified as 1 ime-based updaes are made a fixed ime inervals, 2 movemen-based updaes are made whenever he number of cell-boundary crossings since he las updae exceeds a specified hreshold, 3 disance-based a mobile updaes is locaion whenever is disance from an expeced locaion exceeds a specified value since is las updae [9], 4 dead reckoning a scheme o predic a locaion based on velociy componens calculaions from previous locaions, and 5 mobiliy-aware-based updaes are made based on he mobiliy paern. As an example of ime-based scheme, saic fixed rae (SFR performs he localizaion periodically wih a fixed ime period T [1]. The energy consumpion of he SFR mehod is independen of mobiliy. However, he accuracy of he locaion varies wih he mobiliy of he sensor. Dynamic velociy monoonic (DVM [1] is an example of dead reckoning. In DVM, a sensor node adaps is localizaion period as a funcion of is mobiliy speed. Based on he speed value a node schedule he nex ime o localize. A parameer α is se o represen he arge maximum error. For each localizaion measuremen, he speed value is compared wih α o esimae how much ime needed o reach he arge maximum error if he node coninues on he same speed. The nex localizaion is scheduled afer ha ime period. Noe ha he assumpion is for consan velociy be-

3 Mobiliy-Paern Based Localizaion Updae Algorihms for Mobile Wireless Sensor Neworks 3 ween wo measuremen poins and his could affec he scheduling ime for localizaion such ha a node reaches values of error higher han he arge maximum error (hreshold error. To accommodae for exremely low/high speeds, DVM assigns an upper and lower limi for localizaion periods. A more complex mehod is called mobiliy aware dead reckoning driven (MADRD [1]. MADRD uses he dead reckoning model (DRM [3, 8] o calculae he Euclidean disance beween he prediced locaion and he acual curren locaion. If he difference exceeds a predeermined hreshold error he process moves ino a low confidence sae where localizaion is performed wih a higher frequency. Oherwise, he process moves ino a high confidence sae where localizaion is performed wih a lower frequency. Since MADRD essenially is based he predicion echnique, is complexiy inherenly limis is wide use. Recenly, Song e al. proposed an algorihm o rack moving objecs in locaion-based services [2]. In heir design, he mobiliy paern is divided ino hree ypes: Pause, Linear, and Random. Differen localizaion updae policies are used for differen mobiliy paerns. However, his design in [2] is more suiable for cellular neworks. For example, he disance is defined in erms of he number of cells. 3 Mobiliy-Paern Based Localizaion Updae Algorihms 3.1 Main Ideas Moivaed by Song s work in [2], we divide he mobiliy paern of mobile sensors ino hree saes: Pause, Linear, and Random. A modified ime-based scheme is used for Pause sae since here is no disance or movemen involved in a pause. The localizaion period is increased as he pause ime increase. To preven a long wai ime here is a maximum period of wai ime ha canno be exceeded. Dead reckoning scheme is mos suiable for linear mobiliy paern since i works for predicive movemen paern where here is no change of velociy or unprediced change of direcion. The localizaion period is incremened as long as he disance beween he prediced locaion and real locaion is wihin a predefined hreshold. The localizaion is se o he iniial localizaion period if he error exceeds he hreshold value. The mos suiable policy for he random movemen paern is disance-based scheme. Our disance-based scheme is differen han he one used in [2]. In [2] he disance is beween he las updae in he las cell and he curren cell. The idea in our design is no o acually measure he disance beween las updae and curren posiion bu raher o predic when he node will cross he limi of his disance based on he acceleraion value measuremen of he node. The acceleraion a is measured periodically. The localizaion period is hen measured as ˆ = d a hresh where d hresh denoes he disance limied o localize. Afer he node localizes, d hresh will be as he radius lengh of a circle where he node is on he origin. The value of (1

4 4 Mohammad Y. Al-laho, Min Song, Jun Wang d hresh will depend on he applicaions. Applicaions ha need more accuracy will have smaller d hresh. For simulaion of discree ime inervals, ˆ should be convered o discree ime inervals by using (2. n= ˆ imeslo During he n ime inervals he locaions are esimaed as follows, ( x = x + v + v imeslo (3 ˆ 1 x ( y y v v imeslo ˆ = 1 + y + where x and y are he new esimaed posiion coordinaes, and x -1 and y -1 are he previous posiion coordinaes, v x and v y are he velociy componens. Since a equals o he incremen or decremen of velociy per second, a is muliplied by he imeslo period o ge he incremen or decremen in velociy vecor for each ime inerval, ˆv = a imeslo. (2 (4 3.2 Mobiliy Models The mobiliy model used in his sudy is Gaussian Markovian Random Waypoin [6]. I is a hybrid of he Random Waypoin and Gaussian Markovian models. In Random Waypoin model he Mobile node chooses a random desinaion in he simulaion area and a speed ha is uniformly disribued. When i reaches he chosen desinaion, i pauses for a specific period of ime. Then he node chooses a random desinaion again. All desinaions are randomly chosen from wihin a predefined area. In Gaussian Markovian model each node iniially is assigned a curren speed and direcion. The speed (s and direcion (d are updaed a fixed inervals as follows: s = αs 1 + (1 α s + 1 α s 2 n n x n 1 (5 d = αd 1 + (1 α d + 1 α d 2 n n x n 1 (6 where α is a value from o 1 and i is he uning parameer used o vary he randomness: α = leads o very random moion, while α = 1 leads o compleely linear moion; s and d are he mean value of speed and direcion as n ; dx n 1 sx n 1 are wo random variables following a Gaussian disribuion. A each ime inerval he nex locaion is calculaed as follows: x = x + s cos d (7 n n 1 n 1 n 1 and

5 Mobiliy-Paern Based Localizaion Updae Algorihms for Mobile Wireless Sensor Neworks 5 y = y + s sin d (8 n n 1 n 1 n 1 The idea of Gaussian Markovian Random Waypoin is o use random waypoin a he macro level and use Gaussian Markovian a he micro level o provide a noisy behavior. Tha is, he mobile node sill chooses a random desinaion locaion and speed, and hen compues he direcion of ravel. However, i follows he model presened in equaions (5 and (6 o ravel in his direcion in small ime seps. The model is iniialized wih s and d being equal o he chosen speed and direcion. The node ravels following his model for an amoun of ime (T equal o wha i would ake for i o reach he curren desinaion using he chosen speed. The mobile node may no reach he desinaion, bu possibly a locaion close o i depending on he amoun of noise. Afer he ime T he node pauses according o a given pause ime, chooses a random desinaion locaion, and repeas he above process. In our sudy, we use he Gaussian Markovian Random Waypoin wih some aleraion. For Linear sae we se α value equal o 1 in equaions (5 and (6 so i behaves as pure linear movemen. However, he dn-1 in (6 should be recalculaed o he new locaion. In Random sae he α value is se for values beween.1 and.9, because a value of zero will make i a compleely random and memoryless model ha does no depend on he previous speed and direcion. Whenever he mobile node goes in pause ime, i auomaically swiches o he pause sae for ha amoun of ime. 3.3 Analyical and Sysem Models Condiions o move from one sae o anoher will be esed in he curren sae so he updae policy is changed wih a sae change. Fig. 1 (a shows he iniial and hree mobiliy paern saes. In he acual algorihm he main es o ransi from sae o anoher is based on velociy informaion. A es for he pause sae is o have same locaion deeced wo imes in a row. A same locaion means he velociy value equals o zero. In Linear sae velociy should be he same for cerain ime, which is consan velociy or equivalenly acceleraion value equals o zero. For Random sae, he es is o have velociy change which indicaes random mobiliy paern. Fig. 1 (b shows he sae-based mobiliy model. We have he seady-sae ransiion probabiliy vecor π for each sae: PP LP RP + PP LP RL + PP LR RP π P = ( PPR PLR( PLP + PLR + PRL + ( PLP + PLR + PPL( PRP + PRL + PLR PPRPLP + PPRPLR + PLRPPL π R = ( PPR PLR( PLP + PLR + PRL + ( PLP + PLR + PPL( PRP + PRL + PLR PPRPRL + PRPPPL + PRLPPL π L = ( PPR PLR( PLP + PLR + PRL + ( PLP + PLR + PPL( PRP + PRL + PLR (9 (1 (11

6 6 Mohammad Y. Al-laho, Min Song, Jun Wang Equaions (9, (1, and (11 are used o calculae he performance measuremen specified in Secion 3.5. Iniial (Deecs curren mobiliy paern Linear (L P LL Localizaion updae saes Velociy change P LR Linear Consan velociy Random P PL P LP P RL Consan velociy Velociy = Velociy change Velociy = P PP Pause (P P PR Random (R P RR P RP Pause Fig. 1. (a Mobiliy paern ransiion diagram, (b sae-based mobiliy model 3.4 Pseudo-codes of he Algorihms Iniial Sae Algorihm 1. define Ini_updaeime // one ime slo 2. define max_ime //maximum updae waiing ime 3. define Lin_hresh 4. define d_hresh 5. main 6. loc1 loc_updae 7. wai for Ini_updaeime // (x 1, y 1 8. loc2 loc_updae // (x 2, y v1 ((x 1 x 2 2 +(y 1 y 2 2 /( 2 1 if (v1= hen 11. sae Pause sae(, v, loc1 12. else 13. wai for Ini_updaeime 14. loc3 loc_updae // (x 3, y v2 (x2 x3²+(y2 y3²/( if (v1=v2 hen 17. sae Linear sae(, loc else 19. sae Random sae(, v -1, loc -1 Pause Sae Algorihm 1. Pause sae(, v -1, loc updae ime Ini_updaeime 3. loop 4. wai for updae ime 5. loc o loc_updae // Localize 6. calculae v o // from loc -1 and loc o 7. if (v o!= for wo consecuive imes hen 8. if v -1!= v o 9. sae Random sae

7 Mobiliy-Paern Based Localizaion Updae Algorihms for Mobile Wireless Sensor Neworks 7 1. else 11. sae Linear sae if (updae_ime < max_ime hen updae_ime++ Linear Sae Algorihm 1. Linear sae(, loc updae_ime Ini_updaeime 3. loc o loc_updae // ( x o, o y Localize 4. calculae v loop // from loc -1 and loc o 6. compue velociy componens vx, vy perform predicion (dead reckoning //(x pred, y pred wai for updae_ime 9. loc o loc_updae // ( x o, y o 1. calculae v o // v -1 is saved from las calculaion 11. if (v o!= v -1 hen 12. if (v o = hen 13. sae Pause sae 14. else 15. sae Random sae 16. compue d beween prediced loc and real loc 17. if (d < Lin_hresh hen updae_ime++ else 2. updae_ime Ini_updaeime Random Sae Algorihm 1. Random sae(, v -1, loc updae_ime Ini_updaeime 3. loop 4. loc o loc_updae // ( x o, y o localize 5. calculae v o // from loc o and loc calculae a // acceleraion from v o and v if (v o = v -1 hen 8. if (v o = hen 9. sae Pause sae 1. else 11. sae Linear sae 12. ˆ = d a hresh 13. n = ˆ ( imeslo 14. ˆv= a imeslo 15. for (updae_ime=n,updae_ime>, updae_ime compue velociy componens v x + ˆv, v y + ˆv 17. perform predicion // (x pred, y pred

8 8 Mohammad Y. Al-laho, Min Song, Jun Wang 3.5 Performance Measuremens The wo performance measuremens are he localizaion updaes cos C u and he locaion imprecision cos C e. The imprecision cos is he average of he Euclidean disance beween he acual locaion and he esimaed locaion. We have C C = π c (12 i u i u i S = π c (13 i e i e i S where c u i and c e i are he normalized updae cos and imprecision cos for each sae, which are defined as follows: where n is he number of ime slos, oherwise. n 1 i cu = c n n 1 = (14 c is 1 if he updae occurs a ime slo and { } i c = d k, i L, R e = (15 n 1 { } i c = d l, i P e = (16 where d is he Euclidean disance beween he esimaed locaion and he real locaion a ime slo. k is he number of predicions occurred for Linear and Random saes policies. For Pause sae policy he imprecision cos equaion is calculaed differenly since Pause sae policy does no have predicions, so in (16 l is he number of updaes occurred during Pause sae duraion ime. 3.6 An Example Consider ransiions from sae o sae for 15 ime slos as following: PPPPPPRRRRRRRRLLLLLLLLLLLRRRRRRPPPPPPPPPLLLLLLLLLLLLLL PPPPLLLLLLRRRRRRRLLLLLRRRRRRRRRPPPPPPRRRRRRRPPPPPP Couning he number of ransiions N ij from sae i o sae j gives P L R P [ Nij ] i, j S = L R

9 Mobiliy-Paern Based Localizaion Updae Algorihms for Mobile Wireless Sensor Neworks 9 The sae ransiion probabiliy Pij can be esimaed as: ˆ Nij Pij = Calculaing each sae ransiion probabiliy we ge he following N k S ik marix: P L R P Pˆ = L R We obain π P =.288, π L =.346, and π R =.366 as esimaes of residence probabiliies of each movemen sae. Given he mobiliy paern a each sae, we can easily have he performance measuremens using equaions (12 hrough (15. 4 Simulaion Resuls Wih he mobiliy-paern based localizaion updae algorihms (MBLUA, we make a rade-off beween updae cos and imprecision cos. We expec ha when updae cos is low (no updaing oo much, imprecision cos will increase indicaing oo much error as a node moves. On he oher hand, if updae cos is high (updae oo ofen, he imprecision cos will decrease indicaing a low error. Our experimen is divided ino wo pars. The firs par concerns he simulaion of he node s acual movemen. Our model is based on GMRW mobiliy model [6]. In order o analyze he algorihms, we modify BonnMoion [1], and use our mobiliy model o generae he workload for differen mobiliy paern scenarios. The nex mobiliy paern (he nex sae is randomly chosen. However, he probabiliy ha he nex mobiliy paern is Linear or Random is higher han he probabiliy of Pause hrough he enire simulaion. The probabiliy of Linear or Random is.4 for each, and he probabiliy of Pause is.2. The second par is localizaion updae algorihms, and hey are run over he workload generaed from he mobiliy model. Localizaion updae algorihms include algorihm for hree saes combined using he Sae-based Mobiliy Model (SMM, and algorihm for each movemen sae. The purpose of each sae algorihm alone is o compare i wih MBLUA. We use a simulaion area of 2 by 2 meers and simulae he mobiliy model over 1 nodes for 9 seconds. The localizaion scheme should no affec our proposed model, since he model akes care of he localizaion period only and i is mobiliy based dependen. Fig. 2 (a and (b show he race resuls for one slow node and one fas node, respecively, including he acual movemen race from modified GMRW model and he MBLUA race from localizaion updae algorihms for wo differen speed ranges. For slow node, he speed range is form.5 m/s o 1 m/s; for fas node, he speed range is from 2 m/s o 4 m/s. The maximum pause ime is 6 seconds. The resuls indicae a perfec mach beween acual movemen and MBLUA race when he movemen is

10 1 Mohammad Y. Al-laho, Min Song, Jun Wang linear and i ends o have some mismach when he node urns abruply. In Fig. 2 (b, he error increases for fas node. Mobiliy Trace Mobiliy Trace 2 Acual MBUA 2 Acual MBLUA Y(m 1 Y(m X(m (a Slow speed (.5-1 m/s (Max pause ime 6 sec. (max_ime = 5 ime slos, Lin_hresh = 4m, and d_hresh = 5m Fig. 2. Mobiliy race for slow and fas speeds X(m (b Fas speed (2-4 m/s (Max pause ime 6 sec. (max_ime = 5 ime slos, Lin_hresh = 4m, and d_hresh = 5m Table 1 shows he updae and imprecision coss for hree updae policies and our model for he fas and he slow movemen wih he same parameers. Pause sae policy has he leas updae cos bu he mos imprecision cos. Depending on he mobiliy moion percenage, if he moion is more linear, hen he lowes imprecision cos is Linear sae policy; and if he moion is more random, hen he lowes imprecision cos is Random sae policy. However, our MBLUA has considerably low updae cos wih an accepable imprecision cos beween Linear and Random policies. Table 1. Policies updae and imprecision coss comparison Policy Fas node (2-4 m/s Slow node (.5-1 m/s Updae cos Imprecision cos Updae cos Imprecision cos MBLUA Pause-sae Linear-sae Random-sae Fig. 3 (a and (b show he updae and imprecision coss for MBLUA model when varying he d_hresh. I was observed ha, when d_hresh is low, he updae cos is high indicaing more frequen updaes, and he imprecision cos is low since locaion knowledge is updaed more frequen. As d_hresh increases, he updae cos decreases and he imprecision cos increases, since a larger d_hresh value resuls in larger wai ime o localize. We have found ha changing Lin_hresh only have lile impac on he resuls because Linear sae policy increase he waiing ime o he maximum ime as long as he predicion error is less han Lin_hresh, which is he case ha he error value is close o zero for linear movemen. From he resuls of Fig. 3, i is observed ha any d_hresh value greaer han 5 does no have effec on he updae cos meanwhile imprecision cos does increase. The reason is ha he maximum ime (max_ime o localize is fixed in he simulaions. As a

11 Mobiliy-Paern Based Localizaion Updae Algorihms for Mobile Wireless Sensor Neworks 11 a resul, when d_hresh is greaer han5, he sensor node will localize afer max_ime period no maer how large d_hresh value is, which sabilizes he updae cos. Consequenly, when he updae cos sabilizes, he main facor of changing he imprecision cos is he acceleraion value. Updae cos Disance hreshold (m Slow node Fas node (a Updae cos vs. d_hresh (max_ime=5 ime slos and Lin_hresh=4m Imprecision cos Disance hreshold (m Slow node Fas node (b Imprecision cos vs. d_hresh (max_ime=5 ime slos and Lin_hresh=4m Fig. 3. Updae and imprecision cos as a funcion of d_hresh Updae cos Slow node Fas node Maximun waiing ime slos (a Updae cos vs. max_ime (Lin_hresh=4m, d_hresh=5m Imprecision cos Maximum waiing ime slos Slow node Fas node (b Imprecision cos vs. max_ime (Lin_hresh=4m, d_hresh=5m Fig. 4. Updae and imprecision coss as a funcion of maximum waiing ime slos (max_ime Fig. 4 shows he updae and imprecision coss for MBLUA model when varying he maximum wai ime o localize. Wih low value of maximum wai ime, he updae cos is high and he imprecision cos is low since we force he algorihm o localize very ofen. As he maximum ime increases he updae cos decreases and he imprecision cos increases because he algorihm will have more possible ime o increase he wai ime o localize and hence will decrease he localizaion updae frequency. Bu his causes more error as he wai ime increases.

12 12 Mohammad Y. Al-laho, Min Song, Jun Wang Our experimens show ha fas movemen has more updae and imprecision coss. When a mobile node moves faser, he esimaion deviaes from he real locaion in higher raes, which causes larger esimaion errors. Moreover, as long as he error value exceeds he predefined error hresholds (d_hresh or Lin_hresh, a node will localize more frequenly. 5 Conclusions We designed hree mobiliy-paern based localizaion updae algorihms for mobile wireless sensor neworks. The main idea is o divide sensor movemens ino hree saes, Pause, Linear, and Random. Based on he naure of each movemen paern, differen localizaion updae algorihm is adoped a each sae. To verify he design, we have performed boh mahemaical analysis and simulaions. The resuls demonsrae how he localizaion cos is minimized and how he locaion accuracy is improved. This research is significan o conserve he power consumpion in sensor nodes and o rack he locaion of mobile sensors in a real-ime manner. References 1. S. Tilak, V. Kolar, N. B. Abu-Ghazaleh, and K. D. Kang, Dynamic Localizaion Conrol for Mobile Sensor Neworks, Proc. of IEEE Inernaional Workshop on Sraegies for Energy Efficiency in Ad Hoc and Sensor Neworks, Apr M. Song, K. Park, J. Ryu, and C. Hwang, Modeling and Tracking Complexly Moving Objecs in Locaion-Based Services, Journal of Informaion Science and Engineering, vol. 2, no. 3, pp , May V. Kumar, and S. R. Das, Performance of dead reckoning-based locaion service for mobile ad hoc neworks, Wireless Communicaions and Mobile Compuing Journal, vol. 4, no. 2, pp , Mar L. Hu, and D. Evans, Localizaion for Mobile Sensor Neworks, Proc. of Tenh Annual Inernaional Conference on Mobile Compuing and Neworking, pp , Sep I. F. Akyildiz, W. Su, Y. Sankarasubramaniam, and E. Cayirci, Wireless sensor neworks: a survey, Compuer Neworks, vol. 38, no. 4, pp , Mar T. Camp, J. Boleng, and V. Davies, A Survey of Mobiliy Models for Ad Hoc Nework Research, Wireless Communicaion & Mobile Compuing, Special issue on Mobile Ad Hoc Neworking: Research, Trends and Applicaions, vol. 2, no. 5, pp , Apr A. Savvides, C. C. Han, and M. B. Srivasava, Dynamic fine-grained localizaion in ad-hoc neworks of sensors, Proc. of he 7h ACM/IEEE Inernaional Conference on Mobile Compuing and Neworking (MobiCom, pp , Jul A. Agarwal, and S. R. Das, Dead reckoning in mobile ad hoc neworks, Proc. of he IEEE Wireless Communicaions and Neworking Conference, Mar V. W. S. Wong, and V. C. M. Leung, An adapive disance-based locaion updae algorihm for nex-generaion PCS neworks, IEEE J. Selec Areas on Communicaions, vol. 19, no. 1, pp , Oc BonnMoion, hp://web.informaik.uni-bonn.de/iv/miarbeier/dewaal/bonnmoion/

Two Popular Bayesian Estimators: Particle and Kalman Filters. McGill COMP 765 Sept 14 th, 2017

Two Popular Bayesian Estimators: Particle and Kalman Filters. McGill COMP 765 Sept 14 th, 2017 Two Popular Bayesian Esimaors: Paricle and Kalman Filers McGill COMP 765 Sep 14 h, 2017 1 1 1, dx x Bel x u x P x z P Recall: Bayes Filers,,,,,,, 1 1 1 1 u z u x P u z u x z P Bayes z = observaion u =

More information

Vehicle Arrival Models : Headway

Vehicle Arrival Models : Headway Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where

More information

CHAPTER 10 VALIDATION OF TEST WITH ARTIFICAL NEURAL NETWORK

CHAPTER 10 VALIDATION OF TEST WITH ARTIFICAL NEURAL NETWORK 175 CHAPTER 10 VALIDATION OF TEST WITH ARTIFICAL NEURAL NETWORK 10.1 INTRODUCTION Amongs he research work performed, he bes resuls of experimenal work are validaed wih Arificial Neural Nework. From he

More information

Modal identification of structures from roving input data by means of maximum likelihood estimation of the state space model

Modal identification of structures from roving input data by means of maximum likelihood estimation of the state space model Modal idenificaion of srucures from roving inpu daa by means of maximum likelihood esimaion of he sae space model J. Cara, J. Juan, E. Alarcón Absrac The usual way o perform a forced vibraion es is o fix

More information

Notes on Kalman Filtering

Notes on Kalman Filtering Noes on Kalman Filering Brian Borchers and Rick Aser November 7, Inroducion Daa Assimilaion is he problem of merging model predicions wih acual measuremens of a sysem o produce an opimal esimae of he curren

More information

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,

More information

15. Vector Valued Functions

15. Vector Valued Functions 1. Vecor Valued Funcions Up o his poin, we have presened vecors wih consan componens, for example, 1, and,,4. However, we can allow he componens of a vecor o be funcions of a common variable. For example,

More information

Zürich. ETH Master Course: L Autonomous Mobile Robots Localization II

Zürich. ETH Master Course: L Autonomous Mobile Robots Localization II Roland Siegwar Margaria Chli Paul Furgale Marco Huer Marin Rufli Davide Scaramuzza ETH Maser Course: 151-0854-00L Auonomous Mobile Robos Localizaion II ACT and SEE For all do, (predicion updae / ACT),

More information

Random Walk with Anti-Correlated Steps

Random Walk with Anti-Correlated Steps Random Walk wih Ani-Correlaed Seps John Noga Dirk Wagner 2 Absrac We conjecure he expeced value of random walks wih ani-correlaed seps o be exacly. We suppor his conjecure wih 2 plausibiliy argumens and

More information

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid

More information

KINEMATICS IN ONE DIMENSION

KINEMATICS IN ONE DIMENSION KINEMATICS IN ONE DIMENSION PREVIEW Kinemaics is he sudy of how hings move how far (disance and displacemen), how fas (speed and velociy), and how fas ha how fas changes (acceleraion). We say ha an objec

More information

Kinematics Vocabulary. Kinematics and One Dimensional Motion. Position. Coordinate System in One Dimension. Kinema means movement 8.

Kinematics Vocabulary. Kinematics and One Dimensional Motion. Position. Coordinate System in One Dimension. Kinema means movement 8. Kinemaics Vocabulary Kinemaics and One Dimensional Moion 8.1 WD1 Kinema means movemen Mahemaical descripion of moion Posiion Time Inerval Displacemen Velociy; absolue value: speed Acceleraion Averages

More information

Some Basic Information about M-S-D Systems

Some Basic Information about M-S-D Systems Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,

More information

Object tracking: Using HMMs to estimate the geographical location of fish

Object tracking: Using HMMs to estimate the geographical location of fish Objec racking: Using HMMs o esimae he geographical locaion of fish 02433 - Hidden Markov Models Marin Wæver Pedersen, Henrik Madsen Course week 13 MWP, compiled June 8, 2011 Objecive: Locae fish from agging

More information

Introduction to Probability and Statistics Slides 4 Chapter 4

Introduction to Probability and Statistics Slides 4 Chapter 4 Inroducion o Probabiliy and Saisics Slides 4 Chaper 4 Ammar M. Sarhan, asarhan@mahsa.dal.ca Deparmen of Mahemaics and Saisics, Dalhousie Universiy Fall Semeser 8 Dr. Ammar Sarhan Chaper 4 Coninuous Random

More information

AP CALCULUS AB 2003 SCORING GUIDELINES (Form B)

AP CALCULUS AB 2003 SCORING GUIDELINES (Form B) SCORING GUIDELINES (Form B) Quesion A blood vessel is 6 millimeers (mm) long Disance wih circular cross secions of varying diameer. x (mm) 6 8 4 6 Diameer The able above gives he measuremens of he B(x)

More information

Estimation of Poses with Particle Filters

Estimation of Poses with Particle Filters Esimaion of Poses wih Paricle Filers Dr.-Ing. Bernd Ludwig Chair for Arificial Inelligence Deparmen of Compuer Science Friedrich-Alexander-Universiä Erlangen-Nürnberg 12/05/2008 Dr.-Ing. Bernd Ludwig (FAU

More information

!!"#"$%&#'()!"#&'(*%)+,&',-)./0)1-*23)

!!#$%&#'()!#&'(*%)+,&',-)./0)1-*23) "#"$%&#'()"#&'(*%)+,&',-)./)1-*) #$%&'()*+,&',-.%,/)*+,-&1*#$)()5*6$+$%*,7&*-'-&1*(,-&*6&,7.$%$+*&%'(*8$&',-,%'-&1*(,-&*6&,79*(&,%: ;..,*&1$&$.$%&'()*1$$.,'&',-9*(&,%)?%*,('&5

More information

Article from. Predictive Analytics and Futurism. July 2016 Issue 13

Article from. Predictive Analytics and Futurism. July 2016 Issue 13 Aricle from Predicive Analyics and Fuurism July 6 Issue An Inroducion o Incremenal Learning By Qiang Wu and Dave Snell Machine learning provides useful ools for predicive analyics The ypical machine learning

More information

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still. Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in

More information

In this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should

In this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should Cambridge Universiy Press 978--36-60033-7 Cambridge Inernaional AS and A Level Mahemaics: Mechanics Coursebook Excerp More Informaion Chaper The moion of projeciles In his chaper he model of free moion

More information

Biol. 356 Lab 8. Mortality, Recruitment, and Migration Rates

Biol. 356 Lab 8. Mortality, Recruitment, and Migration Rates Biol. 356 Lab 8. Moraliy, Recruimen, and Migraion Raes (modified from Cox, 00, General Ecology Lab Manual, McGraw Hill) Las week we esimaed populaion size hrough several mehods. One assumpion of all hese

More information

Tracking. Announcements

Tracking. Announcements Tracking Tuesday, Nov 24 Krisen Grauman UT Ausin Announcemens Pse 5 ou onigh, due 12/4 Shorer assignmen Auo exension il 12/8 I will no hold office hours omorrow 5 6 pm due o Thanksgiving 1 Las ime: Moion

More information

EE 315 Notes. Gürdal Arslan CLASS 1. (Sections ) What is a signal?

EE 315 Notes. Gürdal Arslan CLASS 1. (Sections ) What is a signal? EE 35 Noes Gürdal Arslan CLASS (Secions.-.2) Wha is a signal? In his class, a signal is some funcion of ime and i represens how some physical quaniy changes over some window of ime. Examples: velociy of

More information

Lab 10: RC, RL, and RLC Circuits

Lab 10: RC, RL, and RLC Circuits Lab 10: RC, RL, and RLC Circuis In his experimen, we will invesigae he behavior of circuis conaining combinaions of resisors, capaciors, and inducors. We will sudy he way volages and currens change in

More information

Stability and Bifurcation in a Neural Network Model with Two Delays

Stability and Bifurcation in a Neural Network Model with Two Delays Inernaional Mahemaical Forum, Vol. 6, 11, no. 35, 175-1731 Sabiliy and Bifurcaion in a Neural Nework Model wih Two Delays GuangPing Hu and XiaoLing Li School of Mahemaics and Physics, Nanjing Universiy

More information

STATE-SPACE MODELLING. A mass balance across the tank gives:

STATE-SPACE MODELLING. A mass balance across the tank gives: B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing

More information

Recursive Least-Squares Fixed-Interval Smoother Using Covariance Information based on Innovation Approach in Linear Continuous Stochastic Systems

Recursive Least-Squares Fixed-Interval Smoother Using Covariance Information based on Innovation Approach in Linear Continuous Stochastic Systems 8 Froniers in Signal Processing, Vol. 1, No. 1, July 217 hps://dx.doi.org/1.2266/fsp.217.112 Recursive Leas-Squares Fixed-Inerval Smooher Using Covariance Informaion based on Innovaion Approach in Linear

More information

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon 3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of

More information

SPH3U: Projectiles. Recorder: Manager: Speaker:

SPH3U: Projectiles. Recorder: Manager: Speaker: SPH3U: Projeciles Now i s ime o use our new skills o analyze he moion of a golf ball ha was ossed hrough he air. Le s find ou wha is special abou he moion of a projecile. Recorder: Manager: Speaker: 0

More information

Probabilistic Robotics

Probabilistic Robotics Probabilisic Roboics Bayes Filer Implemenaions Gaussian filers Bayes Filer Reminder Predicion bel p u bel d Correcion bel η p z bel Gaussians : ~ π e p N p - Univariae / / : ~ μ μ μ e p Ν p d π Mulivariae

More information

Augmented Reality II - Kalman Filters - Gudrun Klinker May 25, 2004

Augmented Reality II - Kalman Filters - Gudrun Klinker May 25, 2004 Augmened Realiy II Kalman Filers Gudrun Klinker May 25, 2004 Ouline Moivaion Discree Kalman Filer Modeled Process Compuing Model Parameers Algorihm Exended Kalman Filer Kalman Filer for Sensor Fusion Lieraure

More information

Rapid Termination Evaluation for Recursive Subdivision of Bezier Curves

Rapid Termination Evaluation for Recursive Subdivision of Bezier Curves Rapid Terminaion Evaluaion for Recursive Subdivision of Bezier Curves Thomas F. Hain School of Compuer and Informaion Sciences, Universiy of Souh Alabama, Mobile, AL, U.S.A. Absrac Bézier curve flaening

More information

d = ½(v o + v f) t distance = ½ (initial velocity + final velocity) time

d = ½(v o + v f) t distance = ½ (initial velocity + final velocity) time BULLSEYE Lab Name: ANSWER KEY Dae: Pre-AP Physics Lab Projecile Moion Weigh = 1 DIRECTIONS: Follow he insrucions below, build he ramp, ake your measuremens, and use your measuremens o make he calculaions

More information

Georey E. Hinton. University oftoronto. Technical Report CRG-TR February 22, Abstract

Georey E. Hinton. University oftoronto.   Technical Report CRG-TR February 22, Abstract Parameer Esimaion for Linear Dynamical Sysems Zoubin Ghahramani Georey E. Hinon Deparmen of Compuer Science Universiy oftorono 6 King's College Road Torono, Canada M5S A4 Email: zoubin@cs.orono.edu Technical

More information

RC, RL and RLC circuits

RC, RL and RLC circuits Name Dae Time o Complee h m Parner Course/ Secion / Grade RC, RL and RLC circuis Inroducion In his experimen we will invesigae he behavior of circuis conaining combinaions of resisors, capaciors, and inducors.

More information

Air Traffic Forecast Empirical Research Based on the MCMC Method

Air Traffic Forecast Empirical Research Based on the MCMC Method Compuer and Informaion Science; Vol. 5, No. 5; 0 ISSN 93-8989 E-ISSN 93-8997 Published by Canadian Cener of Science and Educaion Air Traffic Forecas Empirical Research Based on he MCMC Mehod Jian-bo Wang,

More information

Application of a Stochastic-Fuzzy Approach to Modeling Optimal Discrete Time Dynamical Systems by Using Large Scale Data Processing

Application of a Stochastic-Fuzzy Approach to Modeling Optimal Discrete Time Dynamical Systems by Using Large Scale Data Processing Applicaion of a Sochasic-Fuzzy Approach o Modeling Opimal Discree Time Dynamical Sysems by Using Large Scale Daa Processing AA WALASZE-BABISZEWSA Deparmen of Compuer Engineering Opole Universiy of Technology

More information

SEIF, EnKF, EKF SLAM. Pieter Abbeel UC Berkeley EECS

SEIF, EnKF, EKF SLAM. Pieter Abbeel UC Berkeley EECS SEIF, EnKF, EKF SLAM Pieer Abbeel UC Berkeley EECS Informaion Filer From an analyical poin of view == Kalman filer Difference: keep rack of he inverse covariance raher han he covariance marix [maer of

More information

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x WEEK-3 Reciaion PHYS 131 Ch. 3: FOC 1, 3, 4, 6, 14. Problems 9, 37, 41 & 71 and Ch. 4: FOC 1, 3, 5, 8. Problems 3, 5 & 16. Feb 8, 018 Ch. 3: FOC 1, 3, 4, 6, 14. 1. (a) The horizonal componen of he projecile

More information

Dynamic Econometric Models: Y t = + 0 X t + 1 X t X t k X t-k + e t. A. Autoregressive Model:

Dynamic Econometric Models: Y t = + 0 X t + 1 X t X t k X t-k + e t. A. Autoregressive Model: Dynamic Economeric Models: A. Auoregressive Model: Y = + 0 X 1 Y -1 + 2 Y -2 + k Y -k + e (Wih lagged dependen variable(s) on he RHS) B. Disribued-lag Model: Y = + 0 X + 1 X -1 + 2 X -2 + + k X -k + e

More information

Mechanical Fatigue and Load-Induced Aging of Loudspeaker Suspension. Wolfgang Klippel,

Mechanical Fatigue and Load-Induced Aging of Loudspeaker Suspension. Wolfgang Klippel, Mechanical Faigue and Load-Induced Aging of Loudspeaker Suspension Wolfgang Klippel, Insiue of Acousics and Speech Communicaion Dresden Universiy of Technology presened a he ALMA Symposium 2012, Las Vegas

More information

0.1 MAXIMUM LIKELIHOOD ESTIMATION EXPLAINED

0.1 MAXIMUM LIKELIHOOD ESTIMATION EXPLAINED 0.1 MAXIMUM LIKELIHOOD ESTIMATIO EXPLAIED Maximum likelihood esimaion is a bes-fi saisical mehod for he esimaion of he values of he parameers of a sysem, based on a se of observaions of a random variable

More information

The Contradiction within Equations of Motion with Constant Acceleration

The Contradiction within Equations of Motion with Constant Acceleration The Conradicion wihin Equaions of Moion wih Consan Acceleraion Louai Hassan Elzein Basheir (Daed: July 7, 0 This paper is prepared o demonsrae he violaion of rules of mahemaics in he algebraic derivaion

More information

Adaptive Noise Estimation Based on Non-negative Matrix Factorization

Adaptive Noise Estimation Based on Non-negative Matrix Factorization dvanced cience and Technology Leers Vol.3 (ICC 213), pp.159-163 hp://dx.doi.org/1.14257/asl.213 dapive Noise Esimaion ased on Non-negaive Marix Facorizaion Kwang Myung Jeon and Hong Kook Kim chool of Informaion

More information

NEWTON S SECOND LAW OF MOTION

NEWTON S SECOND LAW OF MOTION Course and Secion Dae Names NEWTON S SECOND LAW OF MOTION The acceleraion of an objec is defined as he rae of change of elociy. If he elociy changes by an amoun in a ime, hen he aerage acceleraion during

More information

Announcements. Recap: Filtering. Recap: Reasoning Over Time. Example: State Representations for Robot Localization. Particle Filtering

Announcements. Recap: Filtering. Recap: Reasoning Over Time. Example: State Representations for Robot Localization. Particle Filtering Inroducion o Arificial Inelligence V22.0472-001 Fall 2009 Lecure 18: aricle & Kalman Filering Announcemens Final exam will be a 7pm on Wednesday December 14 h Dae of las class 1.5 hrs long I won ask anyhing

More information

Introduction to Mobile Robotics

Introduction to Mobile Robotics Inroducion o Mobile Roboics Bayes Filer Kalman Filer Wolfram Burgard Cyrill Sachniss Giorgio Grisei Maren Bennewiz Chrisian Plagemann Bayes Filer Reminder Predicion bel p u bel d Correcion bel η p z bel

More information

Lecture 3: Exponential Smoothing

Lecture 3: Exponential Smoothing NATCOR: Forecasing & Predicive Analyics Lecure 3: Exponenial Smoohing John Boylan Lancaser Cenre for Forecasing Deparmen of Managemen Science Mehods and Models Forecasing Mehod A (numerical) procedure

More information

Particle Swarm Optimization

Particle Swarm Optimization Paricle Swarm Opimizaion Speaker: Jeng-Shyang Pan Deparmen of Elecronic Engineering, Kaohsiung Universiy of Applied Science, Taiwan Email: jspan@cc.kuas.edu.w 7/26/2004 ppso 1 Wha is he Paricle Swarm Opimizaion

More information

Inventory Control of Perishable Items in a Two-Echelon Supply Chain

Inventory Control of Perishable Items in a Two-Echelon Supply Chain Journal of Indusrial Engineering, Universiy of ehran, Special Issue,, PP. 69-77 69 Invenory Conrol of Perishable Iems in a wo-echelon Supply Chain Fariborz Jolai *, Elmira Gheisariha and Farnaz Nojavan

More information

Kriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds

Kriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 Kriging Models Predicing Arazine Concenraions in Surface Waer Draining Agriculural Waersheds Paul L. Mosquin, Jeremy Aldworh, Wenlin Chen Supplemenal Maerial Number

More information

Shiva Akhtarian MSc Student, Department of Computer Engineering and Information Technology, Payame Noor University, Iran

Shiva Akhtarian MSc Student, Department of Computer Engineering and Information Technology, Payame Noor University, Iran Curren Trends in Technology and Science ISSN : 79-055 8hSASTech 04 Symposium on Advances in Science & Technology-Commission-IV Mashhad, Iran A New for Sofware Reliabiliy Evaluaion Based on NHPP wih Imperfec

More information

An introduction to the theory of SDDP algorithm

An introduction to the theory of SDDP algorithm An inroducion o he heory of SDDP algorihm V. Leclère (ENPC) Augus 1, 2014 V. Leclère Inroducion o SDDP Augus 1, 2014 1 / 21 Inroducion Large scale sochasic problem are hard o solve. Two ways of aacking

More information

A Hop Constrained Min-Sum Arborescence with Outage Costs

A Hop Constrained Min-Sum Arborescence with Outage Costs A Hop Consrained Min-Sum Arborescence wih Ouage Coss Rakesh Kawara Minnesoa Sae Universiy, Mankao, MN 56001 Email: Kawara@mnsu.edu Absrac The hop consrained min-sum arborescence wih ouage coss problem

More information

m = 41 members n = 27 (nonfounders), f = 14 (founders) 8 markers from chromosome 19

m = 41 members n = 27 (nonfounders), f = 14 (founders) 8 markers from chromosome 19 Sequenial Imporance Sampling (SIS) AKA Paricle Filering, Sequenial Impuaion (Kong, Liu, Wong, 994) For many problems, sampling direcly from he arge disribuion is difficul or impossible. One reason possible

More information

PHYSICS 220 Lecture 02 Motion, Forces, and Newton s Laws Textbook Sections

PHYSICS 220 Lecture 02 Motion, Forces, and Newton s Laws Textbook Sections PHYSICS 220 Lecure 02 Moion, Forces, and Newon s Laws Texbook Secions 2.2-2.4 Lecure 2 Purdue Universiy, Physics 220 1 Overview Las Lecure Unis Scienific Noaion Significan Figures Moion Displacemen: Δx

More information

State-Space Models. Initialization, Estimation and Smoothing of the Kalman Filter

State-Space Models. Initialization, Estimation and Smoothing of the Kalman Filter Sae-Space Models Iniializaion, Esimaion and Smoohing of he Kalman Filer Iniializaion of he Kalman Filer The Kalman filer shows how o updae pas predicors and he corresponding predicion error variances when

More information

Maintenance Models. Prof. Robert C. Leachman IEOR 130, Methods of Manufacturing Improvement Spring, 2011

Maintenance Models. Prof. Robert C. Leachman IEOR 130, Methods of Manufacturing Improvement Spring, 2011 Mainenance Models Prof Rober C Leachman IEOR 3, Mehods of Manufacuring Improvemen Spring, Inroducion The mainenance of complex equipmen ofen accouns for a large porion of he coss associaed wih ha equipmen

More information

Robotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1.

Robotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1. Roboics I April 11, 017 Exercise 1 he kinemaics of a 3R spaial robo is specified by he Denavi-Harenberg parameers in ab 1 i α i d i a i θ i 1 π/ L 1 0 1 0 0 L 3 0 0 L 3 3 able 1: able of DH parameers of

More information

AP Calculus BC Chapter 10 Part 1 AP Exam Problems

AP Calculus BC Chapter 10 Part 1 AP Exam Problems AP Calculus BC Chaper Par AP Eam Problems All problems are NO CALCULATOR unless oherwise indicaed Parameric Curves and Derivaives In he y plane, he graph of he parameric equaions = 5 + and y= for, is a

More information

Ground Rules. PC1221 Fundamentals of Physics I. Kinematics. Position. Lectures 3 and 4 Motion in One Dimension. A/Prof Tay Seng Chuan

Ground Rules. PC1221 Fundamentals of Physics I. Kinematics. Position. Lectures 3 and 4 Motion in One Dimension. A/Prof Tay Seng Chuan Ground Rules PC11 Fundamenals of Physics I Lecures 3 and 4 Moion in One Dimension A/Prof Tay Seng Chuan 1 Swich off your handphone and pager Swich off your lapop compuer and keep i No alking while lecure

More information

Using the Kalman filter Extended Kalman filter

Using the Kalman filter Extended Kalman filter Using he Kalman filer Eended Kalman filer Doz. G. Bleser Prof. Sricker Compuer Vision: Objec and People Tracking SA- Ouline Recap: Kalman filer algorihm Using Kalman filers Eended Kalman filer algorihm

More information

Bias in Conditional and Unconditional Fixed Effects Logit Estimation: a Correction * Tom Coupé

Bias in Conditional and Unconditional Fixed Effects Logit Estimation: a Correction * Tom Coupé Bias in Condiional and Uncondiional Fixed Effecs Logi Esimaion: a Correcion * Tom Coupé Economics Educaion and Research Consorium, Naional Universiy of Kyiv Mohyla Academy Address: Vul Voloska 10, 04070

More information

Lab #2: Kinematics in 1-Dimension

Lab #2: Kinematics in 1-Dimension Reading Assignmen: Chaper 2, Secions 2-1 hrough 2-8 Lab #2: Kinemaics in 1-Dimension Inroducion: The sudy of moion is broken ino wo main areas of sudy kinemaics and dynamics. Kinemaics is he descripion

More information

Chapter 15: Phenomena. Chapter 15 Chemical Kinetics. Reaction Rates. Reaction Rates R P. Reaction Rates. Rate Laws

Chapter 15: Phenomena. Chapter 15 Chemical Kinetics. Reaction Rates. Reaction Rates R P. Reaction Rates. Rate Laws Chaper 5: Phenomena Phenomena: The reacion (aq) + B(aq) C(aq) was sudied a wo differen emperaures (98 K and 35 K). For each emperaure he reacion was sared by puing differen concenraions of he 3 species

More information

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles Diebold, Chaper 7 Francis X. Diebold, Elemens of Forecasing, 4h Ediion (Mason, Ohio: Cengage Learning, 006). Chaper 7. Characerizing Cycles Afer compleing his reading you should be able o: Define covariance

More information

Embedded Systems and Software. A Simple Introduction to Embedded Control Systems (PID Control)

Embedded Systems and Software. A Simple Introduction to Embedded Control Systems (PID Control) Embedded Sysems and Sofware A Simple Inroducion o Embedded Conrol Sysems (PID Conrol) Embedded Sysems and Sofware, ECE:3360. The Universiy of Iowa, 2016 Slide 1 Acknowledgemens The maerial in his lecure

More information

Block Diagram of a DCS in 411

Block Diagram of a DCS in 411 Informaion source Forma A/D From oher sources Pulse modu. Muliplex Bandpass modu. X M h: channel impulse response m i g i s i Digial inpu Digial oupu iming and synchronizaion Digial baseband/ bandpass

More information

A Dynamic Model of Economic Fluctuations

A Dynamic Model of Economic Fluctuations CHAPTER 15 A Dynamic Model of Economic Flucuaions Modified for ECON 2204 by Bob Murphy 2016 Worh Publishers, all righs reserved IN THIS CHAPTER, OU WILL LEARN: how o incorporae dynamics ino he AD-AS model

More information

Sequential Importance Resampling (SIR) Particle Filter

Sequential Importance Resampling (SIR) Particle Filter Paricle Filers++ Pieer Abbeel UC Berkeley EECS Many slides adaped from Thrun, Burgard and Fox, Probabilisic Roboics 1. Algorihm paricle_filer( S -1, u, z ): 2. Sequenial Imporance Resampling (SIR) Paricle

More information

Decimal moved after first digit = 4.6 x Decimal moves five places left SCIENTIFIC > POSITIONAL. a) g) 5.31 x b) 0.

Decimal moved after first digit = 4.6 x Decimal moves five places left SCIENTIFIC > POSITIONAL. a) g) 5.31 x b) 0. PHYSICS 20 UNIT 1 SCIENCE MATH WORKSHEET NAME: A. Sandard Noaion Very large and very small numbers are easily wrien using scienific (or sandard) noaion, raher han decimal (or posiional) noaion. Sandard

More information

Let us start with a two dimensional case. We consider a vector ( x,

Let us start with a two dimensional case. We consider a vector ( x, Roaion marices We consider now roaion marices in wo and hree dimensions. We sar wih wo dimensions since wo dimensions are easier han hree o undersand, and one dimension is a lile oo simple. However, our

More information

Probabilistic Robotics SLAM

Probabilistic Robotics SLAM Probabilisic Roboics SLAM The SLAM Problem SLAM is he process by which a robo builds a map of he environmen and, a he same ime, uses his map o compue is locaion Localizaion: inferring locaion given a map

More information

Solutions for Assignment 2

Solutions for Assignment 2 Faculy of rs and Science Universiy of Torono CSC 358 - Inroducion o Compuer Neworks, Winer 218 Soluions for ssignmen 2 Quesion 1 (2 Poins): Go-ack n RQ In his quesion, we review how Go-ack n RQ can be

More information

Sub Module 2.6. Measurement of transient temperature

Sub Module 2.6. Measurement of transient temperature Mechanical Measuremens Prof. S.P.Venkaeshan Sub Module 2.6 Measuremen of ransien emperaure Many processes of engineering relevance involve variaions wih respec o ime. The sysem properies like emperaure,

More information

Optimal Path Planning for Flexible Redundant Robot Manipulators

Optimal Path Planning for Flexible Redundant Robot Manipulators 25 WSEAS In. Conf. on DYNAMICAL SYSEMS and CONROL, Venice, Ialy, November 2-4, 25 (pp363-368) Opimal Pah Planning for Flexible Redundan Robo Manipulaors H. HOMAEI, M. KESHMIRI Deparmen of Mechanical Engineering

More information

Navneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi

Navneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi Creep in Viscoelasic Subsances Numerical mehods o calculae he coefficiens of he Prony equaion using creep es daa and Herediary Inegrals Mehod Navnee Saini, Mayank Goyal, Vishal Bansal (23); Term Projec

More information

L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS. NA568 Mobile Robotics: Methods & Algorithms

L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS. NA568 Mobile Robotics: Methods & Algorithms L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS NA568 Mobile Roboics: Mehods & Algorihms Today s Topic Quick review on (Linear) Kalman Filer Kalman Filering for Non-Linear Sysems Exended Kalman Filer (EKF)

More information

04. Kinetics of a second order reaction

04. Kinetics of a second order reaction 4. Kineics of a second order reacion Imporan conceps Reacion rae, reacion exen, reacion rae equaion, order of a reacion, firs-order reacions, second-order reacions, differenial and inegraed rae laws, Arrhenius

More information

Probabilistic Robotics SLAM

Probabilistic Robotics SLAM Probabilisic Roboics SLAM The SLAM Problem SLAM is he process by which a robo builds a map of he environmen and, a he same ime, uses his map o compue is locaion Localizaion: inferring locaion given a map

More information

Simulating models with heterogeneous agents

Simulating models with heterogeneous agents Simulaing models wih heerogeneous agens Wouer J. Den Haan London School of Economics c by Wouer J. Den Haan Individual agen Subjec o employmen shocks (ε i, {0, 1}) Incomplee markes only way o save is hrough

More information

Chapter 2. First Order Scalar Equations

Chapter 2. First Order Scalar Equations Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.

More information

From Particles to Rigid Bodies

From Particles to Rigid Bodies Rigid Body Dynamics From Paricles o Rigid Bodies Paricles No roaions Linear velociy v only Rigid bodies Body roaions Linear velociy v Angular velociy ω Rigid Bodies Rigid bodies have boh a posiion and

More information

Applications in Industry (Extended) Kalman Filter. Week Date Lecture Title

Applications in Industry (Extended) Kalman Filter. Week Date Lecture Title hp://elec34.com Applicaions in Indusry (Eended) Kalman Filer 26 School of Informaion echnology and Elecrical Engineering a he Universiy of Queensland Lecure Schedule: Week Dae Lecure ile 29-Feb Inroducion

More information

Computation of the Effect of Space Harmonics on Starting Process of Induction Motors Using TSFEM

Computation of the Effect of Space Harmonics on Starting Process of Induction Motors Using TSFEM Journal of elecrical sysems Special Issue N 01 : November 2009 pp: 48-52 Compuaion of he Effec of Space Harmonics on Saring Process of Inducion Moors Using TSFEM Youcef Ouazir USTHB Laboraoire des sysèmes

More information

Demodulation of Digitally Modulated Signals

Demodulation of Digitally Modulated Signals Addiional maerial for TSKS1 Digial Communicaion and TSKS2 Telecommunicaion Demodulaion of Digially Modulaed Signals Mikael Olofsson Insiuionen för sysemeknik Linköpings universie, 581 83 Linköping November

More information

An recursive analytical technique to estimate time dependent physical parameters in the presence of noise processes

An recursive analytical technique to estimate time dependent physical parameters in the presence of noise processes WHAT IS A KALMAN FILTER An recursive analyical echnique o esimae ime dependen physical parameers in he presence of noise processes Example of a ime and frequency applicaion: Offse beween wo clocks PREDICTORS,

More information

A Shooting Method for A Node Generation Algorithm

A Shooting Method for A Node Generation Algorithm A Shooing Mehod for A Node Generaion Algorihm Hiroaki Nishikawa W.M.Keck Foundaion Laboraory for Compuaional Fluid Dynamics Deparmen of Aerospace Engineering, Universiy of Michigan, Ann Arbor, Michigan

More information

Wednesday, November 7 Handout: Heteroskedasticity

Wednesday, November 7 Handout: Heteroskedasticity Amhers College Deparmen of Economics Economics 360 Fall 202 Wednesday, November 7 Handou: Heeroskedasiciy Preview Review o Regression Model o Sandard Ordinary Leas Squares (OLS) Premises o Esimaion Procedures

More information

Appendix to Creating Work Breaks From Available Idleness

Appendix to Creating Work Breaks From Available Idleness Appendix o Creaing Work Breaks From Available Idleness Xu Sun and Ward Whi Deparmen of Indusrial Engineering and Operaions Research, Columbia Universiy, New York, NY, 127; {xs2235,ww24}@columbia.edu Sepember

More information

Matlab and Python programming: how to get started

Matlab and Python programming: how to get started Malab and Pyhon programming: how o ge sared Equipping readers he skills o wrie programs o explore complex sysems and discover ineresing paerns from big daa is one of he main goals of his book. In his chaper,

More information

Robust estimation based on the first- and third-moment restrictions of the power transformation model

Robust estimation based on the first- and third-moment restrictions of the power transformation model h Inernaional Congress on Modelling and Simulaion, Adelaide, Ausralia, 6 December 3 www.mssanz.org.au/modsim3 Robus esimaion based on he firs- and hird-momen resricions of he power ransformaion Nawaa,

More information

Nature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time.

Nature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time. Supplemenary Figure 1 Spike-coun auocorrelaions in ime. Normalized auocorrelaion marices are shown for each area in a daase. The marix shows he mean correlaion of he spike coun in each ime bin wih he spike

More information

Evaluation of Mean Time to System Failure of a Repairable 3-out-of-4 System with Online Preventive Maintenance

Evaluation of Mean Time to System Failure of a Repairable 3-out-of-4 System with Online Preventive Maintenance American Journal of Applied Mahemaics and Saisics, 0, Vol., No., 9- Available online a hp://pubs.sciepub.com/ajams/// Science and Educaion Publishing DOI:0.69/ajams--- Evaluaion of Mean Time o Sysem Failure

More information

15210 RECORDING TIMER - AC STUDENT NAME:

15210 RECORDING TIMER - AC STUDENT NAME: CONTENTS 15210 RECORDING TIMER - AC STUDENT NAME: REQUIRED ACCESSORIES o C-Clamps o Ruler or Meer Sick o Mass (200 g or larger) OPTIONAL ACCESSORIES o Ticker Tape Dispenser (#15225) o Consan Speed Vehicle

More information

IB Physics Kinematics Worksheet

IB Physics Kinematics Worksheet IB Physics Kinemaics Workshee Wrie full soluions and noes for muliple choice answers. Do no use a calculaor for muliple choice answers. 1. Which of he following is a correc definiion of average acceleraion?

More information

Econ107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8)

Econ107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8) I. Definiions and Problems A. Perfec Mulicollineariy Econ7 Applied Economerics Topic 7: Mulicollineariy (Sudenmund, Chaper 8) Definiion: Perfec mulicollineariy exiss in a following K-variable regression

More information

Particle Swarm Optimization Combining Diversification and Intensification for Nonlinear Integer Programming Problems

Particle Swarm Optimization Combining Diversification and Intensification for Nonlinear Integer Programming Problems Paricle Swarm Opimizaion Combining Diversificaion and Inensificaion for Nonlinear Ineger Programming Problems Takeshi Masui, Masaoshi Sakawa, Kosuke Kao and Koichi Masumoo Hiroshima Universiy 1-4-1, Kagamiyama,

More information