Shape-Based Similarity Query for Trajectory of Mobile Objects

Size: px
Start display at page:

Download "Shape-Based Similarity Query for Trajectory of Mobile Objects"

Transcription

1 Shape-Based Similariy Query for Trajecory of Mobile Objecs uaka anagisawa, Jun-ichi Akahani, and Tesuji Saoh NTT Communicaion Science Laboraories, NTT Corporaion hp:// Absrac. In his paper, we describe an efficien indexing mehod for a shape-based similariy search of he rajecory of dynamically changing locaions of people and mobile objecs.in order o manage rajecories in daabase sysems, we define a daa model of rajecories as direced lines in a space, and he similariy beween rajecories is defined as he Euclidean disance beween direced discree lines.our proposed similariy query can be used o find ineresed paerns embedded ino he rajecories, for example, he rajecories of mobile cars in a ciy may include paerns for expecing raffic jams.furhermore, we propose an efficien indexing mehod o rerieve similar rajecories for a query by combining a spaial indexing echnique (R + -Tree) and a dimension reducion echnique, which is called PAA (Piecewise Approximae Aggregae).The indexing mehod can efficienly rerieve rajecories whose shape in a space is similar o he shape of a candidae rajecory from he daabase. 1 Inroducion Recenly, many locaion sensors such as GPS have been developed, and we can obain he rajecory of users and moving objecs using hese sensors [12]. Trajecory daa are widely used in locaion-aware sysems [1], car navigaion sysems, and oher locaion-based informaion sysems, ha can provide services according o a user s curren locaion. These applicaions have sored in hem a lo of rajecories, and hese rajecories may include ineresing individual paerns of each user. For example, by analyzing rajecories of users who work in a building, we can find imporan passages, rooms, sairs, and oher faciliies ha are used frequenly. The resul of he analysis can be used for he managemen and mainenance of he buildings. In he case of a navigaion sysem, a driver can check he roue o a ciy by referring o he rajecories of oher users who have driven o he ciy before. In anoher case, we can sudy characerisics o improve performance in a spor by analyzing he moion daa measured by he sensors aached o he bodies of opspors players. There have been many sudies on managing mobile objecs daa (MOD) [2] [8] [11] [14] [16]. One of he mos ineresing of hese is he developmen of efficien M.-S. Chen e al. (Eds.): MDM 2003, LNCS 2574, pp , c Springer-Verlag Berlin Heidelberg 2003

2 64 uaka anagisawa, Jun-ichi Akahani, and Tesuji Saoh mehod o rerieve objecs, which is indicaed by eiher a spaioemporal range query or a spaio-emporal neares neighbor query. Boh queries are defined as he disance beween he rajecory of a mobile objec and an indicaed poin in a space. For example, he range query is generally defined as he query for rerieving all objecs which passed wihin a given disance of an indicaed poin, such as rerieve all of he people who walked wihin one mile of he buildings a he ime. The range query can also be defined as he query o rerieve all objecs ha passed wihin an indicaed polygon. In boh cases, he query is defined using he disance beween figures in a space. These disance-based queries are useful in locaion managemen of mobile objecs [15], however, hese queries do no have enough power o analyze he paern of he objecs moion. As menioned above, because we are ineresed in he exracion of he individual moving paerns of each objec from he rajecories, i is necessary o developmore powerful ools o analyze he rajecories. Hence, we propose shape-based queries of rajecories in space for he analysis, for insance, rerieve all objecs ha have a similar shape o he rajecory where a user walked in a shop. Using his query, we may classify he cusomers in he shopbased on heir shape paerns of he rajecories. In oher words, our approach is based on he shape similariy beween lines, while he exising approaches adop he disance beween poins as he key o rerieve required objecs. I is difficul o define he similariy beween lines in a space. However, we found his useful idea hrough research of ime series daabases [6] [7] [10]. The ime series daabase sysems can sore ime series daa such as emperaure, economic indicaors, populaion, wave signals, and so on, in addiion o supporing queries for exracing paerns from he ime series daa. Mos of he ime series daabase sysems adop he Euclidean disance beween wo ime daa sequences [7] for analysis; if wo sequences, c, c,aregivenas w 1,w 2,...,w n and w 1,w 2,...,w n, he similariy can be defined as D(c, c )= (w 1 w 1 ) (w n w n )2. (In Secion 3, we describe he similariy in deail) Because rajecory is a ype of ime series daa, he ime series daabase is able o deal wih rajecory. However, rajecory no only has a ime series daa feaure, bu also has a direced line in space feaure. For example, i is difficul for he ime series daabase o find daa for a geographic and spaial query. Therefore, in his paper, we presen a daa model for rajecories of mobile daa, and a query based on he disance beween wo rajecories by exending he similariy used in he ime series daabase sysems. Moreover, we propose a new indexing mehod for rerieving required rajecories by queries based on our defined disance beween rajecories. In Secion 2, we describe our proposed daa model for he rajecory. Secion 3 describes he disance beween wo discree direced lines for calculaing similariies beween wo rajecories. In Secion 4, we presen boh he processing mehod for our proposed query and an indexing

3 Shape-Based Similariy Query for Trajecory of Mobile Objecs 65 λ1 λ2 12 λ1 λ4 λ λ4 11 (a) (b) Fig. 1. Trajecory of Mobile Objecs ((a) rajecory in he real world. (b) rajecory sored in a daabase) echnique ha is an exension of Piecewise Aggregae Approximaion (PAA) [7]. Finally he evaluaion of our approach is shown in Secion 5. 2 Trajecory of Mobile Objecs In order o effecively manage mobile objecs, i is necessary o manage he locaion of each objec a each ime. Generally, a locaion managemen sysem can rerieve objecs locaed in he indicaed area a he indicaed ime [2]. However, we are ineresed in he similariy of he rajecory s shape in a space. In order o define he similariy beween rajecories, i is necessary a firs o define he rajecory as a figure drawn in space. Hence, we define he daa model for he rajecory of mobile objecs 1. A real-world rajecory is a direced coninuous line wih a sar and end poin (Figure 1(a)). Given a wo-dimensional space R 2 and a closed ime inerval I λ =[, ]wih<, a rajecory λ is defined as follows, Definiion 1 : Trajecory A rajecory is he image of a coninuous mapping of λ : I λ R 2. This definiion is a emporal exension of he definiion of a simple line described in [3]. Nex, we denoe he lengh of rajecories in R 2 as L S and he inerval of rajecories in emporal space as L T : Definiion 2 : Lengh of Trajecory in Space R 2 The lengh of rajecory λ during a period [ 0, 1 ] is denoed as L S (λ, [ 0, 1 ]) calculaed as follows: 1 (dx ) 2 ( ) 2 dy L S (λ, [ 0, 1 ]) = + d, whereλ() =(x, y) 0 d d The lengh of he whole rajecory is denoed as L S (λ)( = L S (λ, [, ]). 1 For simplificaion of he problem, we jus focus on he rajecory of mobile objecs; in oher words, we do no discuss he daa model of he oher aribues of he objecs, such as shape, name, and so on.

4 66 uaka anagisawa, Jun-ichi Akahani, and Tesuji Saoh Definiion 3 : Temporal Inerval of Trajecory The x =(x, y) is a vecor in space R 2. The emporal inerval of rajecory λ beween x i and x j on λ is defined as follows: L T (λ, [x i, x j ]) = j i, where λ( i )=x i,λ( j )=x j, and i, j I λ [, ] L T (λ) = However, a locaion sensor device such as GPS does no coninuously measure he coordinaes of a mobile objec, bu samples such daa. The measured daa are hus a sequence of coordinaes of posiions shown in Figure 1(b). Hence, we define discree rajecory λ as a discree funcion. Each vecor x i represens a posiion of a mobile objec a each ime T λ = { 0, 1,..., m } in he space. Definiion 4 :Discree Trajecory A discree rajecory is he image of a discree mapping: λ : T λ R 2. A discree rajecory can be represened as a vecor sequence x 1,...,x m, also. If T λ = {1, 2,...,m}, we denoe he discrie rajecory λ as jus a simple vecor sequence x 1,...,x m. Addiionally, where λ( i ) = x i, we inroduce several noaions; T λ(i) = i, λ(i) =x i,and λ is he number of he vecors included in λ ( λ = T λ ). Nex, we define he disance beween wo vecors x, x in R 2. Definiion 5 :Disance of Vecors D(x, x )= (x x ) 2 +(y y ) 2 Tha is, his definiion assumes ha space R 2 is Euclidean. Alhough λ is a discree line, i is necessary o deal wih λ as a coninuous line in a query. In order o saisfy his requiremen, we define a funcion o conver he discree line ino a coninuous line. There are various mehods o calculae an approximae coninuous line from a discree line [4]. In our approach, we adoped he piecewise linear approximaion because of is simpliciy and populariy [15]. Definiion 6 : Piecewise Linear Approximaion of λ λ :[ 0, m ] R 2 is given as λ() if T λ() ={ λ i λ(i i+1 i )+ i+1 λ(i+1 i+1 i )if T λ i mus be seleced under he condiion: i << i+1 In he res of his paper, we mainly discuss he feaures of boh λ and λ. Noe ha in his paper, we only menion he rajecory on R 2, bu our proposed model and echniques can obviously be adaped o he higher dimensional space R n.

5 Shape-Based Similariy Query for Trajecory of Mobile Objecs 67 L1 L2 L1 d L2 i Q k l L3 j L4 (a) (b) Fig. 2. Disance beween rajecories Val L1 Q Q L2 L1 n L3 L4 Time (a) Previous Spaial knn (b) Previous Temporal knn Fig. 3. Exising knn approaches 3 Similariy Query Based on Shapes of Lines 3.1 Shape-Based Approach The similariy query is useful in is own righ as a ool for exploraory daa analysis[7], and i is a significan elemen in many daa mining applicaions. For insance, we may find he opimum arrangemen of iems in a marke by analyzing he rajecories of cusomers walking around in a shop. In addiion o is usefulness in he rajecory daabase, he similariy query is one of he mos ineresing fields in ime series daabases. In ime series daabases, he similariy beween wo ses of ime series daa is ypically measured by he Euclidean disance [6] [7], which can be calculaed efficienly. However, here have been few discussions on he similariy beween wo lines in space because he previous approaches for spaial queries have focused on he disance beween a poin and a line [2] [9] [15]. The ineres of he previous approaches was mainly o find objecs ha pass a poin near he indicaed poin, such as a car passing hrough a sree. On he oher hand, we are ineresed in he shape of he rajecory. In order o calculae shape-based similariies among rajecories, i is necessary o define a new similariy for he rajecories, as shown in Figure 2(b).

6 68 uaka anagisawa, Jun-ichi Akahani, and Tesuji Saoh In general, he similariy query is represened as a k Neares Neighbor Query (knn) [2] [5] [9]. There are wo ypes of exising approaches, one is based on spaial similariies, he oher is based on similariy beween wo ime series daa. The example of he exising spaial knn is illusraed in Figure 3(a). In his case, he answer is L 1,L 2 when K is 2. On he oher hand, he similariy beween wo ime series daa is defined as he Euclidean disance beween wo ime series, where he lengh of each is n. The disance is defined as he Euclidean disance beween wo n-dimensonal vecor daa [7] shown in Figure 3(b). While his disance of he ime series daa is based on shape, he disance is defined only in he case of R 1 T (T =[0, ]), bu no in he case of R n T,shown in Figure 2(b). Since he rajecory has boh spaial and emporal feaures, we consider hree ypes of similariy queries for rajecories as follows: Spaio-Temporal Similariy: based on a spaio-emporal feaure in R 2 T. Spaial Similariy: basedonaspaialonlyfeaureinr 2 wihou emporal feaures. Temporal Similariy: based on a emporal only feaure in R 1 T wihou spaial feaures. In he res of his secion, we define he similariy in he firs wo cases. We do no define he emporal similariy because his similariy is he same as he similariy defined for he ime series daabases. 3.2 Shape-Based Similariy Query As menioned above, he rajecory has a ime series daa feaure. We define he similariy beween wo rajecories in he same manner as for he similariy defined in he ime series query[7]. For he ime series daabase, he similariy of he wo ime series daa, where each has n value, is given by he Euclidean disance beween vecors in R n. In [6] and [7], when here are wo ime series daa, c = w 1,w 2,...,w n, c = w 1,w 2,...,w n, he disance D(c, c ) is defined as follows: D(c, c )= (w 1 w 1 ) (w n w n )2 This definiion can be exended if each vecor x is a vecor in space R 2,when he ime series vecors are = x 1, x 2,...,x n, = x 1, x 2,...,x n,and he disance is D(c, c ). We define he disance beween wo ime series vecors D(, ) by exending he definiion of D(c, c ), as follows: D(, )= D(x 1, x 0 1 ) D(x n, x 0 n )2 Based on his definiion, we consider he shape-based similariy query for rajecories. Here, Λ is he se of discree rjecories sored in he daabase, and each λ i ( λ i Λ) is a discree rajecory, such as λ i = x 1,...,x m.thequery rajecory λ q is given as λ q = x 1,...,x n. The shape-based range query can hen be defined using Λ, λ q, and he previous defined disance beween wo ime series vecors, as follows:

7 Shape-Based Similariy Query for Trajecory of Mobile Objecs 69 Inpu :, λ q and θ (θ is a naural number). Oupu : a, { λ a1,..., λ ak } a. funcion Q range( θ: ineger, λ q, ) : a begin var j : ineger, l := λ q, a := φ; for each λ i in do for j := 1 o λ i l +1do begin λ ij = subsequence( λ i, j, l ); { This funcion will reurn a subsequence of he original sequence λ i, such as x j, x j+1,...,x j+l 1,eachx λ i} if D( λ q, λ ij) <θhen Add λ ij o a; end; reurn end. a; Fig. 4. The process of he shape-based range query of rajecories 9 Query Trajecory 4 8 Subsequence Mach Sored Trajecory A B C (a) (b) Fig. 5. Similariy query for rajecories Definiion 7 : Shape-based Range Query The process for calculaion of he shape-based range query Q range (θ, λ q, Λ) is given in Figure 4. The range query is defined as a subsequence mach of rajecories as shown in Figure 5. In addiion, he neares neighbor query can be defined using our disance beween rajecories. In our definiion, he emporal feaures are no indicaed in he query, however, we consider ha he emporal feaures can be indicaed independenly from he range query. For example, a query Q range ( θ, λ q, Λ) 11:00 < T λai (1) < 12:00 means rerieving subsequences λ ai where he disance beween λ q and λ ai is less han θ. Moreover, he firs vecor in λ ai is measured wihin he inerval [11:00, 12:00].

8 70 uaka anagisawa, Jun-ichi Akahani, and Tesuji Saoh (a) Original Daa (b) Approximae Trajecory d d (c) Temporal Normalizaion (d) Spaial Normalizaion Fig. 6. Normalizaion of rajecories 3.3 Spaio-emporal Disance beween Two Trajecories Our defined disance D(, ) can be used only in he case where each vecor x is measured by he same inerval, ha is = i+1 i (i =1,...,n 1), where i is an inerval from he ime when x i is measured. However, each vecor in he rajecory is no always measured by he same inerval because sensor devices ofen lose he daa. For example, a discree rajecory illusraed in Figure 6(a) has no measured vecors a =5, 7, 9.. Therefore, o calculae he similariy using our definiion, we define a emporal normalized discree rajecory λ for rajecory λ, as follows: Definiion 8 : Temporal Normalized Discree Trajecory Given a rajecory λ defined for ime inerval [ S, E ], and a naural number m, he emporal normalized discree rajecory λ is defined as follows: λ = λ( S ),λ( S + ),...,λ( S + m ),where S + m = E Inuiively, his discree rajecory λ is he re-sampled rajecory per fixed inerval from λ, shown in Figure 6(c). In oher words, λ is generaed by dividing λ ino equal inerval. For discree rajecory λ, we can use he piecewise linear approximaion λ insead of λ. In he case of Figure 6, he emporal normalized discree rajecory (Figure 6(c)) is generaed from he approximae rajecory (Figure 6(b)).

9 Shape-Based Similariy Query for Trajecory of Mobile Objecs 71 Definiion 9 : Spaial-Temporal Similariy beween wo Trajecories Given wo rajecories λ and λ wih he same emporal lengh (i.e. L T (λ) = L T (λ ) ) and a naural number m, he spaio-emporal disance (similariy) D TS (λ, λ ) beween λ and λ is defined as follows: D TS (λ, λ )= 1 m D( m +1 λ (i), λ (i))2, where = L T (λ) m = L T (λ ) m i=0 Noe ha D TS ( λ, λ ) can be defined as D TS ( λ, λ ). In his definiion, he similariy is he Euclidean disance beween rajecories represened as m +1 dimensional vecors, and he inerval of each rajecory is normalized. Using his definiion, i is possible o find rajecories whose shape is more similar o he query rajecory han can be found using previous mehods. 3.4 Spaial Disance beween Two Trajecories In definiion 8, we focused on he shapes of he rajecories in space R 2 T. However, here are cases where he shapes in R 2 (wihou emporal feaures) are jus as imporan, such as for example in he case of finding similar rajecories o hose of a specified user when focusing on he spaial shape. Hence, we also define he spaial similariy beween wo rajecories. Definiion 10 : Spaial Normalized Discree Trajecory Given a rajecory λ and a naural number m, he spaial normalized discree rajecory λ δ is defined as follows; λ δ = λ( 0 ),...,λ( m ),wherel S (λ, [ i 1, i ]) = δ, (i =1,...,m) Similar o λ, λ δ is generaed by dividing λ ino equal spaial lengh δ. In he case of Figure 6, he spaial normalized discree rajecory (Figure 6(d)) is generaed from he approximae rajecory (Figure 6(b)). Definiion 11 : Spaial Similariy beween Trajecories Given wo rajecories λ and λ wih he same spaial lengh (i.e. L S (λ) = L S (λ )) and a naural number m, he spaial disance (similariy) D S (λ, λ ) beween λ and λ is defined as follows: D S (λ, λ )= 1 m D( m +1 λδ (i), λ (i)) 2 where δ = L S(λ) δ m = L S(λ ) m i=0 Using his definiion, i is possible o find he rajecories whose spaial shape is similar o ha of he query rajecory wihou emporal feaures.

10 72 uaka anagisawa, Jun-ichi Akahani, and Tesuji Saoh C1 C1 C1 C1 D(c1, c2) C2 D(c1, c2) C2 C2 C2 (a) (b) (c) Fig. 7. Disance beween wo sequences c 1, c 2, and Disance beween approximae sequences c 1, c 2 4 Indexing Wih our proposed mehod for calculaing similariy beween rajecories, he daabase sysem can find he rajecories ha have similar shapes o he shape of he query rajecory. However, he cos of calculaing our defined similariy is very high, because i is necessary o calculae Euclidean disances beween each poin on he query rajecory and each poin on all rajecories sored in he daabase. In general, daabase sysems sore a lo of rajecories, and he amoun of daa is increasing rapidly. Therefore, i is imporan o reduce he cos of calculaing similariies. In his secion, we presen an indexing mehod o reduce he cos of calculaing similariies, which is based on echniques for reducing he dimensions of vecor daa. 4.1 Piecewise Aggregae Approximaion Piecewise Aggregae Approximaion (PAA) [7] is a echnique for reducing he cos of comparing wo ses of ime series daa. The essenial idea of his echnique is he reducion of he number of compared daa using he lower limi values of ime series daa. Here, we describe only an ouline of his echnique because i was fully presened in [7]. As menioned in Secion 3, he similariy beween wo ime series daa ses can be defined as Euclidean disance beween wo sequences represened as mulidimensional vecors. Even if he query sequence has a shoer lengh m han he candidae sequence, he similariy can be defined as he disance beween he query sequence and each subsequence of he candidae sequence, as illusraed in Figure 3(b). According o he definiion of he similariy beween wo sequences (menioned in Secion 3.3), when he lengh of he query sequence is m and he maximum lengh of a candidae sequence is n, he order of calculaing he similariy is obviously O(mn) for each sored sequence. In order o reduce cos of comparison, i is necessary o reduce he number of compared values in sequence.

11 Shape-Based Similariy Query for Trajecory of Mobile Objecs 73 Original Daa 1. Projecion 2. Approximaion 3. Ploing on Plane 4. Generaing R + -Tree Fig. 8. The process for generaing indexes o rajecories The PAA is a echnique for generaing approximae sequences o efficienly calculae similariy. If he original sequence has n values, he approximae sequence has only k values ( k is a facor of n and k is much less han n ). Each member of he approximae sqeuence c = w,..., w k is given as follows: w i = 1 k ki j=k(i 1)+1 In shor, each w i is calculaed as he average of w k(i 1)+1,...,w ki.moreover,i was proved in [7] ha he approximae sequences c, c have a special relaionship wih he original sequences c, c ( c = c = k): w j D( c, c ) <D(c, c ) This relaionship means he disance beween he approximae sequences is he lower limi of he disance beween he original sequences. For example, in he case of Figure 7(a), he disance beween c 1 and c 2 (Figure 7(b)) is always greaer han he disance beween approximae sequences c 1 and c 2,showninFigure 7(c). Using his resul, he daabase sysem can reduce he number of compared sequences wih he query sequence.

12 74 uaka anagisawa, Jun-ichi Akahani, and Tesuji Saoh 4.2 Exended Indexing Mehod for Shape-Based Similariy Query The PAA is a simple and efficien echnique for reducing he number of compared ime series daa; however, his echnique has only been adaped o he daa in space R 1 T. Hence, we exend his echnique o rajecory daa in space R 2 T. Moreover, we presen an efficien indexing mehod for rajecories by combining wo echniques: PAA and a spaial indexing echnique (R + -Tree [13]) 2 We only describe he case of he spaial similariy query, bu he essenial idea can also be adaped o he spaio-emporal similariy query. Firs, we define he exension of PAA for he 2-dimensional space: Definiion 12 : 2-Dimensional PAA (2D-PAA) Given a normalized spaial rajecory λ δ (L S ( λ δ )=nδ) andk, whichisa facor of n, 2-Dimensional PAA of λ δ is: = x 1, x 2,..., x k such ha x i =( x i, ȳ i ) x i = 1 k ki j=k(i 1)+1 x j, ȳ i = 1 k ki j=k(i 1)+1 y j where x j =(x j,y j ) Inuiively, 2D-PAA can approximae he rajecory s shape by calculaing he cener of he poins conained in he rajecory. Using 2D-PAA and R + -Tree, he indexes o rajecories can be generaed as shown in Figure 8. In order o rerieve rajecories whose disance o he query rajecory is less han θ, he daabase sysem processes he following seps; 1. Calculaing q = x q1,..., x qk for he query rajecory λ q (λ δq ). 2. Searching sequences 1, 2,... (he lengh of each sequence is k) onr + - Tree, such as D S ( i, q ) <θ. 3. Finding answer rajecories λ 1,λ 2,... such as D S (λ i,λ q ) <θ. In our approach, he number of compared daa is reduced in wo seps. The combinaion of boh reducion echniques enables daabases o rerieve answer rajecories efficienly. Therefore he rajecory daabases can suppor he shapebased similariy query wihou a heavy load. Our indexing mehod s performance is evaluaed in Secion 5. 5 Performance Sudy 5.1 Experimenal Seings There are basically hree variables ha could affec he similariy beween rajecories. The firs is he lengh of a query, because as he lengh of query increases, he similariy decreases, bu he cos of calculaing similariies increases. 2 APCA [6], which is an exension of PAA, uses R-Tree based indexing echniques. However, he indexing mehod only uses R-Tree echniques o rerieve ime series daa (in space R 1 T ) efficienly, no o rerieve sequences of vecors such as rajecories in space R 2 T.

13 Shape-Based Similariy Query for Trajecory of Mobile Objecs Time (msec) Lengh Poins Time (msec) Lengh Poins Fig. 9. The comparison of our proposed indexing mehod wih he exising mehod The second variable is he densiy of he poins in a space, since he similariy increases and he number of compared rajecories increases wih he densiy of he poins. The final variable is he complexiy of he rajecory s shape. Generally, he irregulariy of an objec, moion causes greaer complexiy of shapes; in oher words, people walking randomly generae he mos complex shapes. Therefore, we generae sample rajecories by changing hese variables. The number of rajecories is 2-32, he lengh of rajecories is jus 1000 (i.e., he maximum number of poins is 32000), and he sampling inerval is fixed 3. Each rajecory has shapes ha represen people walking freely on a plane while changing speed and direcion, and he frequency of he change in hese values is alered in order o generae various complex shapes. In addiion, he rajecories are embedded ino a fixed area (size is ) wihou any endencies, in oher words, he densiy of poins can be conroled by he number of embedded rajecories. 5.2 Efficiency We have an experimen o assess he performance of our indexing mehod. For he experimen, we implemened wo ypes of engine o find he rajecory which is he neares o he query rajecory. The firs engine checks every poin on all rajecories (wihou any index), while he second engine checks only poins filered by our proposed index. We generaed random rajecories for a query, and measured he calculaion ime required o find he neares rajecory o he generaed query rajecory. Since his is a simple performance evaluaion of our approach, we gave a very simple query for each engine; he lengh of he query is fixed o k. In oher word, he query rajecory λ q has jus k iems, 3 Eveninhecasewhere is no fixed, we can ake similar resuls from he case where is fixed, because fixed rajecories can be mechanically generaed from original rajecories wih he piecewise linear approximaion defined in Secion 2.

14 76 uaka anagisawa, Jun-ichi Akahani, and Tesuji Saoh such as λ q = x 1,...,x k, and he approximae rajecory λ q has only one iem x =(1/k k i=1 x i, 1/k k i=1 y i)wherex i =(x i,y i ). Figure 9 illusraes he experimenal resul in he case where he shape of rajecories is simple. While he lef graph in he figure shows he calculaion ime wih our indexing mehods, he righ graph shows wihou any index. Wih our proposed index, he calculaion ime is 60% 75% less han wihou he index in each siuaion. For example, in he case where he query lengh is 256 and he number of poins is (i is he wors case), he calculaion ime is 3,325 msec wih he index. On he oher hand he calculaion ime wihou he index is 18,206 msec. Alhough he ime aken o generae indexes (overhead) is 3,377 msec when he number of poins is 32,000, i is less han he ime required wihou he index. As a resul of our experimen, he advanage of our indexing mehod is clear. In addiion, we measured he calculaion ime in cases where he shape of rajecories is very complex, such as rajecories where each objec moves almos randomly. In his case, he calculaion ime is 15% 20% greaer han he ime in he simple case; however, he rae of increase is he same in cases wih he index and wihou he index. 6 Conclusion The main conribuion of his paper is he presenaion of an efficien indexing mehod for processing shape-based similariy queries for rajecory daabases. In order o calculae similariy beween rajecories, we defined discree rajecories ha were re-sampled per fixed inerval. Furhermore, we described he performance of our proposed indexing mehod, and we show he advanage of our mehod over he exising mehods. As fuure work, we will implemen he rajecory daabase sysem o evaluae our proposed model, Furhermore, we will develop several real applicaion programs such as a car navigaion sysem, a personal navigaion sysem, and oher locaion-wares. References [1] G.Chen and D.Koz.Caegorizing binary opological relaions beween regions, lines, and poins in geographic daabases.technical Repor TR , A Survey of Conex-Aware Mobile Compuing Research, Dep.of Compuer Science, Darmouh College, [2] H.Chon, D.Agrawal, and A.E.Abbadi. Query processing for moving objecs wih space-ime grid sorage model.in MDM2002 Conference Proceedings, pages , [3] E.Clemenini and P.D.Felice.Topological invarians for lines. IEEE Transacion on Knowledge and Daa Engineering, 10(1):38 54, [4] L.E.Elsgolc. Calculus of Variaions.Pergamon Press LTD, 1961.

15 Shape-Based Similariy Query for Trajecory of Mobile Objecs 77 [5] E.G.Hoel and H.Same. Efficien processing of spaial queries in line segmen daabases. In O.Gunher and H.J.Schek, ediors, SSD 91 Proceedings, volume 525, pages Springer-Verlag, [6] E.Keogh, K.Chakrabari, S.Mehrora, and M.Pazzan. Locally adapive dimensionaliy reducion for indexing large ime series daabases.in SIGMOD2001 Conference Proceedings, pages , [7] E.Keogh, K.Chakrabari, M.Pazzani, and S.Mehrora. Dimensionaliy reducion for fas similariy search in large ime series daabases. Knowledge and Informaion Sysems, 3(3): , [8] G.Kollios, D.Gunopulos, and V.J.Tsoras. On indexing mobile objecs. In SIGMOD 99 Conference Proceedings, pages , [9] G.Kollios, V.J.Tsoras, D.Gunopulos, A.Delis, and M.Hadjielefheriou.Indexing animaed objecs using spaioemporal access mehods. IEEE Transacions on Knowledge and Daa Engineering, 13(5): , [10].-S. Moon, K.-. Whang, and W.-S. Han. General mach: A subsequence maching mehod in ime-series daabases based on generalized windows.in SIGMOD 2002 Conference Proceedings, pages , [11] K.Porkaew, I.Lazaridis, and S.Mehrora. Querying mobile objecs in spaioemporal daabases.in C.S.Jensen, M.Schneider, B.Seeger, and V.J.Tsoras, ediors, SSTD 2001, volume 2121 of Lecure Noes in Compuer Science, pages Springer-Verlag, [12] N.Priyanha, A.Miu, H.Balakrishnan, and S.Teller. The cricke compass for conex-aware mobile applicaions.in MOBICOM2001 Conference Proceedings, pages 1 14, [13] T.Sellis, N.Roussopoulos, and C.Falousos. The R + -ree: A dynamic index for mulidimensional objecs.in VLDB 87 Conference Proceedings, pages 3 11, [14] A.P.Sisla, O.Wolfson, S.Chamberlain, and S.Dao. Modeling and querying moving objecs.in ICDE 97 Proceedings, pages , [15] M.Vazirgiannis and O.Wolfson.A spaioemporal model and language for moving objecs on road neworks. In C.S.Jensen, M.Schneider, B.Seeger, and V.J. Tsoras, ediors, SSTD 2001, volume 2121 of Lecure Noes in Compuer Science, pages Springer-Verlag, [16] O.Wolfson, B.u, S.Chamberlain, and L.Jiang. Moving objecs daabases: Issues and soluions.in Saisical and Scienific Daabase Managemen (SSDM 98) Conference Proceedings, pages , 1998.

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

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

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

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

Linear Response Theory: The connection between QFT and experiments

Linear Response Theory: The connection between QFT and experiments Phys540.nb 39 3 Linear Response Theory: The connecion beween QFT and experimens 3.1. Basic conceps and ideas Q: How do we measure he conduciviy of a meal? A: we firs inroduce a weak elecric field E, and

More information

Retrieval Models. Boolean and Vector Space Retrieval Models. Common Preprocessing Steps. Boolean Model. Boolean Retrieval Model

Retrieval Models. Boolean and Vector Space Retrieval Models. Common Preprocessing Steps. Boolean Model. Boolean Retrieval Model 1 Boolean and Vecor Space Rerieval Models Many slides in his secion are adaped from Prof. Joydeep Ghosh (UT ECE) who in urn adaped hem from Prof. Dik Lee (Univ. of Science and Tech, Hong Kong) Rerieval

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

Overview. COMP14112: Artificial Intelligence Fundamentals. Lecture 0 Very Brief Overview. Structure of this course

Overview. COMP14112: Artificial Intelligence Fundamentals. Lecture 0 Very Brief Overview. Structure of this course OMP: Arificial Inelligence Fundamenals Lecure 0 Very Brief Overview Lecurer: Email: Xiao-Jun Zeng x.zeng@mancheser.ac.uk Overview This course will focus mainly on probabilisic mehods in AI We shall presen

More information

Content-Based Shape Retrieval Using Different Shape Descriptors: A Comparative Study Dengsheng Zhang and Guojun Lu

Content-Based Shape Retrieval Using Different Shape Descriptors: A Comparative Study Dengsheng Zhang and Guojun Lu Conen-Based Shape Rerieval Using Differen Shape Descripors: A Comparaive Sudy Dengsheng Zhang and Guojun Lu Gippsland School of Compuing and Informaion Technology Monash Universiy Churchill, Vicoria 3842

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

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

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

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

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

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

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

Single and Double Pendulum Models

Single and Double Pendulum Models Single and Double Pendulum Models Mah 596 Projec Summary Spring 2016 Jarod Har 1 Overview Differen ypes of pendulums are used o model many phenomena in various disciplines. In paricular, single and double

More information

Christos Papadimitriou & Luca Trevisan November 22, 2016

Christos Papadimitriou & Luca Trevisan November 22, 2016 U.C. Bereley CS170: Algorihms Handou LN-11-22 Chrisos Papadimiriou & Luca Trevisan November 22, 2016 Sreaming algorihms In his lecure and he nex one we sudy memory-efficien algorihms ha process a sream

More information

Analyze patterns and relationships. 3. Generate two numerical patterns using AC

Analyze patterns and relationships. 3. Generate two numerical patterns using AC envision ah 2.0 5h Grade ah Curriculum Quarer 1 Quarer 2 Quarer 3 Quarer 4 andards: =ajor =upporing =Addiional Firs 30 Day 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 andards: Operaions and Algebraic Thinking

More information

Final Spring 2007

Final Spring 2007 .615 Final Spring 7 Overview The purpose of he final exam is o calculae he MHD β limi in a high-bea oroidal okamak agains he dangerous n = 1 exernal ballooning-kink mode. Effecively, his corresponds o

More information

Section 4.4 Logarithmic Properties

Section 4.4 Logarithmic Properties Secion. Logarihmic Properies 59 Secion. Logarihmic Properies In he previous secion, we derived wo imporan properies of arihms, which allowed us o solve some asic eponenial and arihmic equaions. Properies

More information

Logic in computer science

Logic in computer science Logic in compuer science Logic plays an imporan role in compuer science Logic is ofen called he calculus of compuer science Logic plays a similar role in compuer science o ha played by calculus in he physical

More information

Class Meeting # 10: Introduction to the Wave Equation

Class Meeting # 10: Introduction to the Wave Equation MATH 8.5 COURSE NOTES - CLASS MEETING # 0 8.5 Inroducion o PDEs, Fall 0 Professor: Jared Speck Class Meeing # 0: Inroducion o he Wave Equaion. Wha is he wave equaion? The sandard wave equaion for a funcion

More information

Ensemble Confidence Estimates Posterior Probability

Ensemble Confidence Estimates Posterior Probability Ensemble Esimaes Poserior Probabiliy Michael Muhlbaier, Aposolos Topalis, and Robi Polikar Rowan Universiy, Elecrical and Compuer Engineering, Mullica Hill Rd., Glassboro, NJ 88, USA {muhlba6, opali5}@sudens.rowan.edu

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

Section 4.4 Logarithmic Properties

Section 4.4 Logarithmic Properties Secion. Logarihmic Properies 5 Secion. Logarihmic Properies In he previous secion, we derived wo imporan properies of arihms, which allowed us o solve some asic eponenial and arihmic equaions. Properies

More information

5. Stochastic processes (1)

5. Stochastic processes (1) Lec05.pp S-38.45 - Inroducion o Teleraffic Theory Spring 2005 Conens Basic conceps Poisson process 2 Sochasic processes () Consider some quaniy in a eleraffic (or any) sysem I ypically evolves in ime randomly

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

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

ERROR LOCATING CODES AND EXTENDED HAMMING CODE. Pankaj Kumar Das. 1. Introduction and preliminaries

ERROR LOCATING CODES AND EXTENDED HAMMING CODE. Pankaj Kumar Das. 1. Introduction and preliminaries MATEMATIČKI VESNIK MATEMATIQKI VESNIK 70, 1 (2018), 89 94 March 2018 research paper originalni nauqni rad ERROR LOCATING CODES AND EXTENDED HAMMING CODE Pankaj Kumar Das Absrac. Error-locaing codes, firs

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

Speaker Adaptation Techniques For Continuous Speech Using Medium and Small Adaptation Data Sets. Constantinos Boulis

Speaker Adaptation Techniques For Continuous Speech Using Medium and Small Adaptation Data Sets. Constantinos Boulis Speaker Adapaion Techniques For Coninuous Speech Using Medium and Small Adapaion Daa Ses Consaninos Boulis Ouline of he Presenaion Inroducion o he speaker adapaion problem Maximum Likelihood Sochasic Transformaions

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

Unsteady Flow Problems

Unsteady Flow Problems School of Mechanical Aerospace and Civil Engineering Unseady Flow Problems T. J. Craf George Begg Building, C41 TPFE MSc CFD-1 Reading: J. Ferziger, M. Peric, Compuaional Mehods for Fluid Dynamics H.K.

More information

Traveling Waves. Chapter Introduction

Traveling Waves. Chapter Introduction Chaper 4 Traveling Waves 4.1 Inroducion To dae, we have considered oscillaions, i.e., periodic, ofen harmonic, variaions of a physical characerisic of a sysem. The sysem a one ime is indisinguishable from

More information

6.2 Transforms of Derivatives and Integrals.

6.2 Transforms of Derivatives and Integrals. SEC. 6.2 Transforms of Derivaives and Inegrals. ODEs 2 3 33 39 23. Change of scale. If l( f ()) F(s) and c is any 33 45 APPLICATION OF s-shifting posiive consan, show ha l( f (c)) F(s>c)>c (Hin: In Probs.

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

Lecture 4 Notes (Little s Theorem)

Lecture 4 Notes (Little s Theorem) Lecure 4 Noes (Lile s Theorem) This lecure concerns one of he mos imporan (and simples) heorems in Queuing Theory, Lile s Theorem. More informaion can be found in he course book, Bersekas & Gallagher,

More information

The Arcsine Distribution

The Arcsine Distribution The Arcsine Disribuion Chris H. Rycrof Ocober 6, 006 A common heme of he class has been ha he saisics of single walker are ofen very differen from hose of an ensemble of walkers. On he firs homework, we

More information

1.6. Slopes of Tangents and Instantaneous Rate of Change

1.6. Slopes of Tangents and Instantaneous Rate of Change 1.6 Slopes of Tangens and Insananeous Rae of Change When you hi or kick a ball, he heigh, h, in meres, of he ball can be modelled by he equaion h() 4.9 2 v c. In his equaion, is he ime, in seconds; c represens

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

Chapter 3 Boundary Value Problem

Chapter 3 Boundary Value Problem Chaper 3 Boundary Value Problem A boundary value problem (BVP) is a problem, ypically an ODE or a PDE, which has values assigned on he physical boundary of he domain in which he problem is specified. Le

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

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature On Measuring Pro-Poor Growh 1. On Various Ways of Measuring Pro-Poor Growh: A Shor eview of he Lieraure During he pas en years or so here have been various suggesions concerning he way one should check

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

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

Longest Common Prefixes

Longest Common Prefixes Longes Common Prefixes The sandard ordering for srings is he lexicographical order. I is induced by an order over he alphabe. We will use he same symbols (,

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

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

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

THE DISCRETE WAVELET TRANSFORM

THE DISCRETE WAVELET TRANSFORM . 4 THE DISCRETE WAVELET TRANSFORM 4 1 Chaper 4: THE DISCRETE WAVELET TRANSFORM 4 2 4.1 INTRODUCTION TO DISCRETE WAVELET THEORY The bes way o inroduce waveles is hrough heir comparison o Fourier ransforms,

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

10. State Space Methods

10. State Space Methods . Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he

More information

2.1: What is physics? Ch02: Motion along a straight line. 2.2: Motion. 2.3: Position, Displacement, Distance

2.1: What is physics? Ch02: Motion along a straight line. 2.2: Motion. 2.3: Position, Displacement, Distance Ch: Moion along a sraigh line Moion Posiion and Displacemen Average Velociy and Average Speed Insananeous Velociy and Speed Acceleraion Consan Acceleraion: A Special Case Anoher Look a Consan Acceleraion

More information

2. Nonlinear Conservation Law Equations

2. Nonlinear Conservation Law Equations . Nonlinear Conservaion Law Equaions One of he clear lessons learned over recen years in sudying nonlinear parial differenial equaions is ha i is generally no wise o ry o aack a general class of nonlinear

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

This is an example to show you how SMath can calculate the movement of kinematic mechanisms.

This is an example to show you how SMath can calculate the movement of kinematic mechanisms. Dec :5:6 - Kinemaics model of Simple Arm.sm This file is provided for educaional purposes as guidance for he use of he sofware ool. I is no guaraeed o be free from errors or ommissions. The mehods and

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

Two Coupled Oscillators / Normal Modes

Two Coupled Oscillators / Normal Modes Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own

More information

Matrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality

Matrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality Marix Versions of Some Refinemens of he Arihmeic-Geomeric Mean Inequaliy Bao Qi Feng and Andrew Tonge Absrac. We esablish marix versions of refinemens due o Alzer ], Carwrigh and Field 4], and Mercer 5]

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

Isolated-word speech recognition using hidden Markov models

Isolated-word speech recognition using hidden Markov models Isolaed-word speech recogniion using hidden Markov models Håkon Sandsmark December 18, 21 1 Inroducion Speech recogniion is a challenging problem on which much work has been done he las decades. Some of

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

1 Review of Zero-Sum Games

1 Review of Zero-Sum Games COS 5: heoreical Machine Learning Lecurer: Rob Schapire Lecure #23 Scribe: Eugene Brevdo April 30, 2008 Review of Zero-Sum Games Las ime we inroduced a mahemaical model for wo player zero-sum games. Any

More information

12: AUTOREGRESSIVE AND MOVING AVERAGE PROCESSES IN DISCRETE TIME. Σ j =

12: AUTOREGRESSIVE AND MOVING AVERAGE PROCESSES IN DISCRETE TIME. Σ j = 1: AUTOREGRESSIVE AND MOVING AVERAGE PROCESSES IN DISCRETE TIME Moving Averages Recall ha a whie noise process is a series { } = having variance σ. The whie noise process has specral densiy f (λ) = of

More information

) were both constant and we brought them from under the integral.

) were both constant and we brought them from under the integral. YIELD-PER-RECRUIT (coninued The yield-per-recrui model applies o a cohor, bu we saw in he Age Disribuions lecure ha he properies of a cohor do no apply in general o a collecion of cohors, which is wha

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

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

Module 4: Time Response of discrete time systems Lecture Note 2

Module 4: Time Response of discrete time systems Lecture Note 2 Module 4: Time Response of discree ime sysems Lecure Noe 2 1 Prooype second order sysem The sudy of a second order sysem is imporan because many higher order sysem can be approimaed by a second order model

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

Effects of Coordinate Curvature on Integration

Effects of Coordinate Curvature on Integration Effecs of Coordinae Curvaure on Inegraion Chrisopher A. Lafore clafore@gmail.com Absrac In his paper, he inegraion of a funcion over a curved manifold is examined in he case where he curvaure of he manifold

More information

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions Muli-Period Sochasic Models: Opimali of (s, S) Polic for -Convex Objecive Funcions Consider a seing similar o he N-sage newsvendor problem excep ha now here is a fixed re-ordering cos (> 0) for each (re-)order.

More information

Applying Genetic Algorithms for Inventory Lot-Sizing Problem with Supplier Selection under Storage Capacity Constraints

Applying Genetic Algorithms for Inventory Lot-Sizing Problem with Supplier Selection under Storage Capacity Constraints IJCSI Inernaional Journal of Compuer Science Issues, Vol 9, Issue 1, No 1, January 2012 wwwijcsiorg 18 Applying Geneic Algorihms for Invenory Lo-Sizing Problem wih Supplier Selecion under Sorage Capaciy

More information

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals

More information

ADDITIONAL PROBLEMS (a) Find the Fourier transform of the half-cosine pulse shown in Fig. 2.40(a). Additional Problems 91

ADDITIONAL PROBLEMS (a) Find the Fourier transform of the half-cosine pulse shown in Fig. 2.40(a). Additional Problems 91 ddiional Problems 9 n inverse relaionship exiss beween he ime-domain and freuency-domain descripions of a signal. Whenever an operaion is performed on he waveform of a signal in he ime domain, a corresponding

More information

9/9/99 (T.F. Weiss) Signals and systems This subject deals with mathematical methods used to describe signals and to analyze and synthesize systems.

9/9/99 (T.F. Weiss) Signals and systems This subject deals with mathematical methods used to describe signals and to analyze and synthesize systems. 9/9/99 (T.F. Weiss) Lecure #: Inroducion o signals Moivaion: To describe signals, boh man-made and naurally occurring. Ouline: Classificaion ofsignals Building-block signals complex exponenials, impulses

More information

SUPPLEMENTARY INFORMATION

SUPPLEMENTARY INFORMATION SUPPLEMENTARY INFORMATION DOI: 0.038/NCLIMATE893 Temporal resoluion and DICE * Supplemenal Informaion Alex L. Maren and Sephen C. Newbold Naional Cener for Environmenal Economics, US Environmenal Proecion

More information

Numerical Dispersion

Numerical Dispersion eview of Linear Numerical Sabiliy Numerical Dispersion n he previous lecure, we considered he linear numerical sabiliy of boh advecion and diffusion erms when approimaed wih several spaial and emporal

More information

Lecture 20: Riccati Equations and Least Squares Feedback Control

Lecture 20: Riccati Equations and Least Squares Feedback Control 34-5 LINEAR SYSTEMS Lecure : Riccai Equaions and Leas Squares Feedback Conrol 5.6.4 Sae Feedback via Riccai Equaions A recursive approach in generaing he marix-valued funcion W ( ) equaion for i for he

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

Equations of motion for constant acceleration

Equations of motion for constant acceleration Lecure 3 Chaper 2 Physics I 01.29.2014 Equaions of moion for consan acceleraion Course websie: hp://faculy.uml.edu/andriy_danylo/teaching/physicsi Lecure Capure: hp://echo360.uml.edu/danylo2013/physics1spring.hml

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

Efficient Matching for State-Persistent Publish/Subscribe Systems

Efficient Matching for State-Persistent Publish/Subscribe Systems Efficien Maching for Sae-Persisen Publish/Subscribe Sysems Absrac Huber Ka Yau Leung IBM Torono Lab hkyleung@ca.ibm.com Conen-based publish/subscribe sysems allow informaion disseminaion and fine-grained

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

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

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t M ah 5 2 7 Fall 2 0 0 9 L ecure 1 0 O c. 7, 2 0 0 9 Hamilon- J acobi Equaion: Explici Formulas In his lecure we ry o apply he mehod of characerisics o he Hamilon-Jacobi equaion: u + H D u, x = 0 in R n

More information

20. Applications of the Genetic-Drift Model

20. Applications of the Genetic-Drift Model 0. Applicaions of he Geneic-Drif Model 1) Deermining he probabiliy of forming any paricular combinaion of genoypes in he nex generaion: Example: If he parenal allele frequencies are p 0 = 0.35 and q 0

More information

Open loop vs Closed Loop. Example: Open Loop. Example: Feedforward Control. Advanced Control I

Open loop vs Closed Loop. Example: Open Loop. Example: Feedforward Control. Advanced Control I Open loop vs Closed Loop Advanced I Moor Command Movemen Overview Open Loop vs Closed Loop Some examples Useful Open Loop lers Dynamical sysems CPG (biologically inspired ), Force Fields Feedback conrol

More information

Interpretation of special relativity as applied to earth-centered locally inertial

Interpretation of special relativity as applied to earth-centered locally inertial Inerpreaion of special relaiviy as applied o earh-cenered locally inerial coordinae sysems in lobal osiioning Sysem saellie experimens Masanori Sao Honda Elecronics Co., Ld., Oyamazuka, Oiwa-cho, Toyohashi,

More information

A DELAY-DEPENDENT STABILITY CRITERIA FOR T-S FUZZY SYSTEM WITH TIME-DELAYS

A DELAY-DEPENDENT STABILITY CRITERIA FOR T-S FUZZY SYSTEM WITH TIME-DELAYS A DELAY-DEPENDENT STABILITY CRITERIA FOR T-S FUZZY SYSTEM WITH TIME-DELAYS Xinping Guan ;1 Fenglei Li Cailian Chen Insiue of Elecrical Engineering, Yanshan Universiy, Qinhuangdao, 066004, China. Deparmen

More information

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

CHAPTER 2 Signals And Spectra

CHAPTER 2 Signals And Spectra CHAPER Signals And Specra Properies of Signals and Noise In communicaion sysems he received waveform is usually caegorized ino he desired par conaining he informaion, and he undesired par. he desired par

More information

1. VELOCITY AND ACCELERATION

1. VELOCITY AND ACCELERATION 1. VELOCITY AND ACCELERATION 1.1 Kinemaics Equaions s = u + 1 a and s = v 1 a s = 1 (u + v) v = u + as 1. Displacemen-Time Graph Gradien = speed 1.3 Velociy-Time Graph Gradien = acceleraion Area under

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

MODULE - 9 LECTURE NOTES 2 GENETIC ALGORITHMS

MODULE - 9 LECTURE NOTES 2 GENETIC ALGORITHMS 1 MODULE - 9 LECTURE NOTES 2 GENETIC ALGORITHMS INTRODUCTION Mos real world opimizaion problems involve complexiies like discree, coninuous or mixed variables, muliple conflicing objecives, non-lineariy,

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

The average rate of change between two points on a function is d t

The average rate of change between two points on a function is d t SM Dae: Secion: Objecive: The average rae of change beween wo poins on a funcion is d. For example, if he funcion ( ) represens he disance in miles ha a car has raveled afer hours, hen finding he slope

More information

Adaptive Compressive Tracking Based on Perceptual Hash Algorithm Lei ZHANG, Zheng-guang XIE * and Hong-jun LI

Adaptive Compressive Tracking Based on Perceptual Hash Algorithm Lei ZHANG, Zheng-guang XIE * and Hong-jun LI 2017 2nd Inernaional Conference on Informaion Technology and Managemen Engineering (ITME 2017) ISBN: 978-1-60595-415-8 Adapive Compressive Tracking Based on Percepual Hash Algorihm Lei ZHANG, Zheng-guang

More information