Multiple Target Tracking With a 2-D Radar Using the JPDAF Algorithm and Combined Motion Model

Size: px
Start display at page:

Download "Multiple Target Tracking With a 2-D Radar Using the JPDAF Algorithm and Combined Motion Model"

Transcription

1 Mulple Targe Trang Wh a -D Radar Ung he JPDAF Algorhm and Combned Moon Model Alabar Gor Daronolae, Mohammad Bagher Menha, Al Doomohammad ABSTRACT Mulple arge rang (MTT aen no aoun a one of he mo mporan op n rang arge wh radar. In h paper, he MTT problem ued for emang he poon of mulple arge when a -D radar employed o gaher meauremen. To do o, he Jon Probabl Daa Aoaon Fler (JPDAF approah appled o rang he poon of mulple arge. To haraerze he moon of eah arge, wo model are ued. Fr, a mple near onan veloy model ondered and hen o enhane he rang performane, peally, when arge mae maneuverng movemen a varable veloy model propoed. In addon, a ombned model alo propoed o mgae he maneuverng movemen beer. Th new model gve an advanage o explore he movemen of he maneuverng obe whh ommon n many rang problem. Smulaon reul how he uperory of he new moon model and effe n he rang performane of mulple arge. KEYWORDS Mulple arge rang, JPDAF algorhm, daa aoaon, maneuverng movemen. ITRODUCTIO The problem of emang he poon of a arge n phyal envronmen an mporan problem n he mlary feld. Indeed, beaue of he nrn error embedded n he daa olleed by meauremen enor, here range and bearng fnder, ome operaon hould be ondued on he rude daa o reah more aurae nowledge abou he poon of a arge. owaday, varou algorhm are propoed o modfy he auray of loalzaon algorhm. In he feld of arge rang, he Kalman fler approah [] play a gnfan role. Th approah ue a uggeed model derbng he arge movemen and he meauremen olleed by a -D radar and, hen, provde an aurae emaon of he arge poon. ow, he Kalman fler grealy ued n many rang problem []. Alhough h fler ha provded a general raegy o olve loalzaon problem, uffer from ome weanee. The mo gnfan problem of he Kalman mehod weane o he ae emaon of many real yem n he preene of evere nonlneare. Aually, n general nonlnear yem, he poeror probably deny funon of ae gven meauremen anno be preened by a Gauan deny funon and, herefore, he Kalman mehod may reul n unafaory repone. Moreover, he Kalman fler very enve o he nal ondon of ae. Aordngly, Kalman baed mehod do no provde aurae reul n many ommon global loalzaon problem. More reenly, parle fler have been nrodued o emae non-gauan, non-lnear dynam proee [5].The ey dea of parle fler o repreen he ae by e of ample (or parle. The maor advanage of h ehnque ha an repreen mul-modal ae dene, a propery whh ha been hown o nreae he robune of he underlyng ae emaon proe. Moreover, parle fler are le enve o he nal value of ae han he Kalman mehod. Th problem very onderable n applaon n whh he nal value of ae are no nown uh a global loalzaon problem. The aforemenoned power of parle fler n he ae emaon of nonlnear yem ha made hem be appled wh grea ue o dfferen ae emaon problem nludng vual rang [6], moble robo loalzaon [7], onrol and hedulng [8], arge rang [9] and dynam probabl newor []. Mulple arge rang anoher arave ue n he feld of aeropae engneerng. Th problem an be Correpondng Auhor, Al A.Gor wh he Deparmen of Eleral And Compuer Engneerng, MMaer Unvery, Hamlon, Canada (e-mal: gorda@mmaer.a. M. B. Menha wh he Deparmen of Eleral Engneerng, Amrabr Unvery of Tehnology, Tehran, Iran (e-mal:menha@au.a.r. Al Doomohammad wh he Deparmen of Eleral Engneerng, Amrabr Unvery of Tehnology, Tehran, Iran. Amrabr / MISC / Vol. 4 / o. / Fall 9 43

2 nrodued a rang arge poon by daa reeved from a -D radar. Alhough h problem appear mlar o he ommon loalzaon algorhm n whh he dynamal model of he arge nown, he radonal approahe anno be ued beaue he radar doe no have ae o he nernal enor daa of eah arge ued n loalzaon algorhm o pred he fuure poon of he arge. Moreover, ung he parle fler o emae he poon of all arge, nown a ombned ae pae, no raable beaue he omplexy of h approah grow exponenally wh he number of obe beng raed. To reolve he problem of he la of nernal enor daa, varou moon model have been propoed. The mple model propoed he near onan veloy model. Th model ha been grealy ued n many rang applaon uh a arraf rang [], people rang [] and ome oher arge rang applaon [3]. Bede h model, ome oher ruure have been propoed uable n ome oher peal applaon. ear onan aeleraon model a modfed veron of he prevou model whh an be appled o rang maneuverng obe [4]. Alo, nerave mulple model (IMM fler repreen anoher ool o ra he movemen of hghly maneuverng obe where he former propoed approahe are no relable [5]. Reenly, many raege have been alo propoed n he leraure o addre he dffule aoaed wh mularge rang. The ey dea of he leraure o emae eah obe poon eparaely. To do o, he meauremen hould be agned o eah arge. Therefore, daa aoaon ha been nrodued o relae eah arge o aoaed meauremen. By ombnng daa aoaon onep wh flerng and ae emaon algorhm, he mulple arge rang ue an be deal wh eaer han applyng ae emaon proedure o he ombned ae pae model. Among he varou mehod, Jon Probabl Daa Aoaon Fler (JPDAF algorhm aen no aoun a one of he mo arave approahe. The radonal JPDAF [6] ue Exended Kalman Flerng (EKF on wh he daa aoaon o emae ae of mulple obe. However, he performane of he algorhm degrade a he nonlneare beome more evere. To enhane he performane of he JPDAF agan he evere nonlneary of dynamal and enor model, reenly, raege have been propoed o ombne he JPDAF wh parle ehnque nown a Mone Carlo JPDAF, o aommodae general nonlnear and non-gauan model [7]. The newer veron of he above-menoned algorhm ue an approah named a of gang o redue he ompuaonal o whh one of he mo erou problem n onlne rang algorhm [7]. Bede he Mone Carlo JPDAF algorhm, ome oher mehod have been propoed n he leraure. In [7], a omprehenve urvey ha been preened on he mo well-nown algorhm applable n he feld of mulple arge rang. In h paper, he Mone Carlo JPDAF algorhm ued for he mulple arge rang problem. A a modfaon o he reen leraure, dfferen moon model are eed o evaluae he performane of rang. Bede he radonal near onan veloy model, we wll e he effeny of he oher wo model a well n he paper. Fr, a varable veloy model ued o haraerze he movemen of arge. Then, we propoe a ombned approah n whh he moon model a lnear ombnaon of onan and varable veloy model. In addon, a mple algorhm repreened o deermne he wegh agned o eah moon model. I hown ha ome nformaon of he JPDAF algorhm an be ued o deermne he aforemenoned faor. The goal, here, o ufy he effeny of he JPDAF algorhm ung ombned model ompared wh eah ndvdual model. The re of h paper organzed a follow. Seon II derbe he Mone Carlo JPDAF algorhm for mulple arge rang. The JPDAF for mulple arge rang he op of eon III. In h eon, we wll preen how he moon of arge n a real envronmen an be derbed. To do o, near onan veloy and aeleraon model are nrodued o derbe he movemen of moble arge. Then, wll be hown how daa reeved from a -D radar an be ued for rang movng obe. In eon IV, varou mulaon reul are preened o uppor he performane of he algorhm. In h eon, he maneuverng movemen of he arge alo ondered. The reul ealy approve he uable performane of he JPDAF algorhm n rang he moon of moble arge n he mulaed envronmen. Fnally, eon V onlude he paper.. THE JPDAF ALGORITHM FOR MULTIPLE TARGET TRACKIG Th eon preen he JPDAF algorhm ondered a he mo ueful and wdely appled raegy for mul-arge rang under daa aoaon unerany. The Mone Carlo veron of he JPDAF algorhm ue ommon parle fler approah o emae he poeror deny funon of ae gven meauremen p ( x y :. ow, onder he problem of rang obe. x denoe he ae of hee obe a me where,,.., he he arge number. Furhermore, he movemen of eah arge an be ondered by a general nonlnear ae pae model a follow: x f ( x + v + ` y g( x + w ( 44 Amrabr / MISC / Vol. 4 / o. / Fall 9

3 where v and w are whe noe wh he ovarane mare Q and R, repevely. Alo, n he above equaon, f and g are he general nonlnear funon repreenng he dynam behavor of he arge and he enor model. The JPDAF algorhm reurvely updae he margnal flerng drbuon for eah arge p ( x y :,,,.., nead of ompung he on flerng drbuon p ( X y :, X x,..., x. To ompue he above drbuon funon, ome pon have o be aen no onderaon a follow: - I very mporan o now how o agn eah arge ae o a meauremen. Indeed, he enor provde a e of meauremen a eah me ep. The oure of hee meauremen an be he arge or he durbane alo nown a luer. Therefore, a peal proedure needed o agn eah arge o aoaed meauremen. Th proedure degnaed a Daa Aoaon ondered a he ey age of he JPDAF algorhm whh derbed n he nex eon. - Beaue he JPDAF algorhm reurvely updae he emaed ae, a reurve oluon hould be appled o updae ae a eah ample me. Tradonally, Kalman flerng ha been a rong ool for reurvely emang ae of he arge n mul-arge rang enaro. Reenly, parle fler on wh he daa aoaon raegy have provded beer emaon, peally, when he enor model nonlnear. Wh regard o he above pon, he followng eon derbe how parle fler paralleled wh he daa aoaon onep an deal wh he mul-arge rang problem. A. The Parle Fler For Onlne Sae Emaon Conder he problem of onlne ae emaon a ompung he poeror probably deny funon p ( x y :. To provde a reurve formulaon for ompung he above deny funon, he followng age are preened: - Predon age: he predon ep done ndependenly for eah arge a: x y x x x y dx ( : x - Updae age: he updae ep an be derbed a follow: p ( x y : y x x y: (3 The parle fler emae he probably deny funon x y : by amplng from a pef drbuon funon a follow: (4 x y w~ δ ( x x : Amrabr / MISC / Vol. 4 / o. / Fall 9 where,,..., he ample number, w ~ he normalzed mporane wegh and δ he dela dra funon. In he above equaon, he ae x ampled from he propoal deny funon q( x x, y :. By ubung he above equaon n ( and he fa ha he ae are drawn from he propoal funon q, he reurve equaon for he predon ep an be wren a follow:, x x (5 α α, q( x x, y: where x, he -h ample of x.ow, by ung (3 he updae age an be expreed a he reurve adumen of mporane wegh a follow: w α y x (6 By repeang he above proedure a eah me ep he equenal mporane amplng (SIS for onlne ae emaon preened below: The SIS Algorhm. For : nalze he ae x, predon weghα and mporane wegh w.. A eah me ep proeed he followng age: 3. Sample he ae from he propoal deny funon a follow: x ~ q( x x, y : (7 4. Updae predon wegh ung (5. 5. Updae mporane wegh ung (6. 6. ormalze mporane wegh a follow: w (8 w~ w 7. Se + and go o ep 3. In he above, for mply, he ndex,, ha been elmnaed. The man falure of he SIS algorhm he degeneray problem. Tha, afer a few eraon one of he normalzed mporane rao end o, whle he remanng rao go o zero. Th problem aue he varane of he mporane wegh o nreae ohaally over me [9]. To avod he degeneray of he SIS algorhm, a eleon (re-amplng age may be ued o elmnae ample wh low mporane wegh and mulply ample wh hgh mporane wegh. There are many approahe o mplemen re-amplng age [9]. Among hem, he redual re-amplng provde a raghforward raegy o remedy he degeneray problem n he SIS algorhm. The redual reamplng on wh he SIS algorhm an be preened a 45

4 he SIR algorhm a follow: The SIR Algorhm. For : nalze he ae x, predon weghα and mporane wegh w.. A eah me ep, do he followng: 3. perform he SIS algorhm o ample ae x and ompue normalzed mporane wegh w ~. 4. Che he re-amplng reron: a. If eff > hreh follow SIS age Ele b. Implemen redual-reamplng approah o mulply/uppre ae wh hgh/low wegh.. Se new normalzed mporane wegh a w~. 5. Se + and go o ep 3. In above, eff a reron heng he degeneray of algorhm whh an be wren a: (9 eff ( w 46 In [9], a omprehenve duon ha been made on how o mplemen he redual re-amplng age. Bede he SIR algorhm, ome oher approahe have been propoed n he leraure o enhane he qualy of he SIR algorhm uh a Marov Chan Mone Carlo parle fler, Hybrd SIR and auxlary parle fler []. Alhough hee mehod are more aurae han he ommon SIR algorhm, ome oher problem uh a he ompuaonal o are he mo gnfan reaon ha, n many real me applaon uh a onlne rang, he radonal SIR algorhm appled o reurve ae emaon. B. Daa Aoaon In he la eon, he SIR algorhm wa preened for onlne ae emaon. However, he maor problem of he propoed algorhm how o ompue he lelhood funon p ( y x. To do o, an aoaon hould be made beween he meauremen and arge. Brefly, he aoaon an be aegorzed no arge o meauremen and meauremen o arge: Defnon : A arge o meauremen aoaon denoed a ( T M defend by ~ λ { ~ r, M, M T } where ~ r ~ r,..., ~ r K wh: If he -h arge undeeed r If he -h arge ha generaed he -h meauremen where,,...,m and m he number of meauremen a eah me ep and,,..,k whh K he number of arge. Defnon : n a mlar fahon, he meauremen o arge aoaon ( M T defned by λ { r, M, M T } where r r,..., rm and r gven below. If he -h meauremen no aoaed wh any arge r If he -h arge In he boh above equaon, aoaed wh he -h arge meauremen due o arge and m T he number of m he number of luer. I very eay o how ha boh defnon are equvalen, however he dmenon of he arge o meauremen aoaon le han ha of he meauremen o arge aoaon. Therefore, n h paper, he arge o meauremen aoaon ued. ow, he lelhood funon for eah arge an be wren a: m ( y x β + β y x In he above equaon, β defned a he probably ha he h meauremen agned o he h arge. Hene, β wren a: β p ( ~ r y : ( Before elaborang more on he above equaon, he followng defnon needed. Defnon 3: We defne he e Λ ~ a all poble agnmen whh an be made beween meauremen and arge. For example, onder a 3-arge rang problem. Aume ha followng aoaon are reognzed beween arge and meauremen: r, r 3, r 3 ow, he e Λ ~ an be led a gven n able I. TABLE ALL POSSIBLE ASSOCIATIOS BETWEE TARGETS AD MEASUREMETS FOR EXAMPLE TARGET TARGET TARGET Amrabr / MISC / Vol. 4 / o. / Fall 9

5 Here, mean ha he arge ha no been deeed. By ung he above defnon, Eq.( an be rewren a: ~ p ( ~ r y : ~ ~ Λ ~ r p y, ( λ ( λ : The rgh de of he above equaon wren a follow: ~ ~ ~ y y:, λ y:, λ λ y: y: (3 ~ ~ y y:, λ λ where we have ued he fa ha he aoaon veor ~ λ no dependen on he hory of meauremen. Eah of he deny funon n he rgh de of he above equaon aed a follow: - ~ λ The above deny funon beome ~ p ( λ ( ~ ( ~ p r, M, M T (4 p r M, M M M T The ompuaon of eah deny funon n he rgh de of he above equaon raghforward a fully dued n [7]. ~ - y y :, λ Beaue he arge are uppoed o be ndependen, he above deny funon an be wren a follow: ~ M r (5 y y, λ V p y y T ( max ( : : Γ where Vmax he maxmum volume whh n he lne of gh of he enor, Γ he e of all vald r meauremen o arge aoaon and p ( y y : he predve lelhood for he r h arge. ow, onder r. The predve lelhood for he -h arge an be formulaed a: p ( y y y x x y dx (6 : : Boh deny funon n he rgh hand de of he above equaon are emaed ung he ample drawn from he propoed deny funon. However, he man problem how o deermne he aoaon beween meauremen ( y and arge. To do o, a of gang mehod propoed n he followng algorhm:. Conder x, Sof Gang,,,,, a ample drawn from he propoed deny funon.. For :m, do he followng ep for he -h meauremen. 3. For :K, do he followng ep: 4. Compue µ a follow: µ ~α g( x,, where g he enor model and wegh a preened n he SIR algorhm. 5. Compue σ a: σ R +, ~ α he normalzed (, (, g( x µ g( x ~ α µ T (7 (8 6. Compue he dane o he -h arge by: T (9 d ( y µ ( σ ( y µ 7. If d < ε, agn he -h meauremen o he -h arge. I eay o how ha he predve lelhood gven n (6 an be approxmaed by: p ( y y : ( µ, σ ( µ, σ a normal drbuon wh mean where ( µ and he ovarane marx σ. By ompung he predve lelhood, he lelhood deny funon ealy emaed. In he nex ubeon, he JPDAF algorhm preened for mul-arge rang problem. C. The JPDAF Algorhm The mahemaal foundaon of he JPDAF algorhm were dued n he la eon. ow, we are ready o propoe he full JPDAF algorhm for he problem of mulple arge rang. To do o, eah arge haraerzed by he dynamal model nrodued n (. The JPDAF algorhm preened below. The JPDAF Algorhm x,. Inalzaon: nalze ae for eah arge a where,,.., and,,...,k, he predve mporane wegh α and mporane wegh w,,.. A eah me ep, proeed a: 3. For :, do he followng ep. 4. For :K, ondu he followng proedure for eah arge. 5. Sample he new ae from he propoed deny funon a follow: (, x x, y, x ~ q : ( Amrabr / MISC / Vol. 4 / o. / Fall 9 47

6 6. Updae he predve mporane wegh a: (,,,, p x x α α,, q x x, y ( : ( Then, normalze he predve mporane wegh. 7. Ue he ampled ae and new obervaon, y, o onue he aoaon veor for eah arge a R { m, y } where ( refer o he aoaon beween he -h arge and he -h meauremen. To do o, ue he of gang proedure derbed n he prevou ubeon. 8. Conue all poble aoaon for he arge and eablh he e Γ ~ a derbed n defnon Ue (3 and ompue β for eah meauremen where l,,..,m and m he number of meauremen.. Ue ( o ompue he lelhood rao for eah arge a (, p y x.. Compue he mporane wegh for eah arge a:,, (, w α p y x (3 Then, normalze he mporane wegh a follow:, w (4 w~,, w. Proeed wh he re-amplng age. To do o, mplemen he mlar proedure derbed n he SIR algorhm. Aferward, for eah arge he reampled ae an be preened a follow: {,, m (5 w~ (,, x }, x 3. Se he me a + and go o ep 3. In he nex eon, he JPDAF algorhm for mulple ae an be ued for mulple arge rang. In addon, ome well-nown moon model are hen preened o ra he moon of a moble arge. 3. THE JPDAF ALGORITHM FOR MULTIPLE TARGET TRACKIG WITH A -D RADAR ow, we wan o ue he JPDAF algorhm for he problem of mulple arge rang wh a -D radar. To do o, Radar yem ue modulaed waveform and dreve anenna o ranm Eleromagne energy no a pef volume n pae o earh for arge. Obe (arge whn a earh volume wll refle poron of h energy (radar reurn or ehoe ba o he radar. Thee ehoe are hen proeed by he radar reever o exra arge nformaon uh a range, veloy, angular poon and oher arge denfyng haraer. To emae he poon of eah arge by daa olleed va a -D radar, we need a moon model o derbe he behavor of he arge. The mple mehod near onan veloy model, preened a: X + AX + Bu (6 X [ x x y y ] where he yem marxe are defned a: T A T T (7 T B T T where T refer o he ample me. In he above equaon, u whe noe wh zero mean and an arbrary ovarane marx. Beaue he model ondered a a near onan veloy, he ovarane of he addve noe anno be o large. Indeed, h model uable for he movemen of arge wh a onan veloy, whh ommon n many applaon uh a ommeral arraf pah rang [8] and people rang []. The movemen of a arge an be modeled by he above model n many ondon. However, n ome peal ae, he arge movemen anno be haraerzed by a mple near onan veloy model. For example, n many problem, arge may ondu a maneuverng movemen. Therefore, he moon raeory of he arge o omplex ha he elemenary onan veloy model wll no reul n afaory repone. To overome h problem, he varable veloy model propoed. The ey dea behnd h model o ue he arge aeleraon a anoher ae varable whh hould be emaed a well a ompung he veloy and poon of he arge. In oher word, he new ae veor defned a, X x y [ x x a y y a ], where a and are he arge aeleraon along wh he x and y ax, repevely. ow, he near onan aeleraon model an be wren a wha menoned n (7, exep for he yem marxe are defned a follow [4]: x a y 48 Amrabr / MISC / Vol. 4 / o. / Fall 9

7 T a T a a A a ex T (8 ex T b b b B b b3 b3 where parameer of he above equaon are defned a follow: b3 ( ex T a b3, b ( T a, a b (9 T b a In he above equaon, a onan value. The above model an be ued o ra he moon raeory of a maneuverng obe. Moreover, ombnng he reul of he propoed model may enhane he performane of rang. Th dea an be preened a follow: ~ ~ a ~ v X αx + ( α X (3 where X ~ a and X ~ are he emaon of he near onan v aeleraon and near onan veloy model, repevely, and α a number beween and. Bede he above approahe, ome oher mehod have been propoed o deal wh rang of maneuverng obe uh a IMM fler [5] whh very applable n mlary heme where arraf hange her aeleraon uddenly. 4. SIMULATIO RESULTS To evaluae he effeny of he JPDAF algorhm, a 3- arge rang enaro ondered. To prepare he enaro and mplemen mulaon, he parameer are preened for arge ruure and mulaon envronmen a TABLE II. ow, he followng expermen are ondued o evaluae he performane of he JPDAF algorhm n varou uaon. TABLE II SIMULATIO PARAMETERS Parameer Derpon n umber of enor R The maxmum range of max mulaon envronmen Q The ovarane marx of he enor addve noe Sample me for runnng mulaon Smulaon maxmum max runnng me A. Loal Trang for on-maneuverng Movemen The nal poon of he radar e a 5,5, π. 8,4,π, Alo, he poon of arge are ondered a [ ],3, π π 3,,9,, 4 repevely. To run he mulaon, ample me,, and mulaon maxmum runnng me, max, are e a and, repevely. Fg. how he raeore generaed for eah arge. The radar provde he dane from he a obe, wall or obale, and he dynam arge. The radar daa an be repreened a a par of [ r,φ] where φ he angle beween he enor and he referene oordnaon whoe relaonhp wh he arge poon o o [ x x + ( y y ] ( r o (3 y y φ aran o x x o o where [ x, y ] denoe he poon of he enor. Alo, he ovarane of he enor noe uppoed o be:. 5 4 (3 ow, he JPDAF framewor ued o ra he raeory of eah arge. In h mulaon, he luer are aumed o be unformly drbued over he oal volume of he envronmen. To apply he near onan veloy model a derbed n (7, he nal value of he veloy hoen unformly beween and - a follow: x u(, y u(, (33 where u(. he unform drbuon funon. Fg. preen he mulaon reul ung he JPDAF algorhm Amrabr / MISC / Vol. 4 / o. / Fall 9 49

8 wh 5 parle for eah arge rang. Smulaon reul ufy he ably of he JPDAF algorhm n modfyng he nrn error embedded n he daa olleed from he meauremen enor. Fg.. The generaed raeory for arge 5. COCLUSIOS In h paper, he JPDAF algorhm wa fully developed for mulple arge rang. Afer nrodung he mahemaal ba of he JPDAF algorhm, we howed how o apply he JPDAF algorhm o he problem of mulple arge rang; h led o hree dfferen model propoed n he paper. ear onan veloy model, near onan aeleraon model were fr derbed and, hen, a ombnaon of he wo former approahe wa preened o enhane he rang auray. Aferward, mulaon reul were arred ou n wo dfferen enaro. Fr, rang he nonmaneuverng moon of arge wa explored. Laer, he JPDAF algorhm wa eed on rang he moon of arge wh maneuverng movemen. Smulaon reul howed ha he ombned approah provde beer performane han he oher wo propoed raege. Fg.. The emaed raeory ung he JPDAF algorhm B. Loal Trang for Maneuverng Movemen ow, onder a maneuverng movemen for eah arge. Therefore, a varable veloy model ondered for eah arge movemen. Fg. 3 how he generaed raeory for eah arge. ow, we apply he JPDAF algorhm for rang eah arge poon. To do o, 3 approahe are mplemened. ear onan veloy model, near onan aeleraon model and ombned model derbed n he la eon are ued o emae he ae of he arge. Fg. 4-6 how he reul for eah arge eparaely. To ompare he effeny of eah raegy, he obaned error preened for eah approah n able III, IV where he followng reron ued for emang he error: max z ~ z ( e z max, z, { x y } (3 The reul how beer performane of he ombned mehod ompared wh he oher approahe. Fg. 3. The generaed raeory for maneuverng arge Fg. 4. The emaed raeory for arge TABLE III THE ERROR FOR ESTIMATIG x Robo um Conan Conan Combned Aeleraon Veloy Amrabr / MISC / Vol. 4 / o. / Fall 9

9 TABLE IV THE ERROR FOR ESTIMATIG y Robo um Conan Conan Combned Aeleraon Veloy Fg. 6. The emaed raeory for arge 3 Fg. 5. The emaed raeory for arge 6. REFERECES [] R. E. Kalman, and R. S. Buy, ew Reul n Lnear Flerng and Predon, Tran. Ameran Soey of Mehanal Engneer, Sere D, Journal of Ba Engneerng, Vol. 83D, pp. 95 8, 96. [] Brano R, Saneev Arulampalam, and el Gordon, Beyond he Kalman Fler: Parle Fler for Trang Applaon, Areh Houe, Boon, London, 4. [3] R. Segwar and I. R. ourbahh, Inoduon o Auonomou Moble Robo, MIT Pre, 4. [4] A. Howard, Mul- robo Smulaneou Loalzaon and Mappng ung Parle Fler, Robo: Sene and Syem I, pp. 8, 5. [5].J. Gordon, D.J. Salmond, and A.F.M. Smh, A ovel Approah o onlnear/non-gauan Bayean Sae Emaon, IEE Proeedng F, 4(:7 3, 993. [6] M. Iard, and A. Blae, Conour Trang by Soha Propagaon of Condonal Deny, In Pro. of he European Conferene of Compuer Von, 996. [7] D. Fox, S. Thrun, F. Dellaer, andw. Burgard, Parle Fler for Moble Robo Loalzaon, Sequenal Mone Carlo Mehod n Prae. Sprnger Verlag, ew Yor,. [8] C. Hue, J. P. L. Cadre, and P. Perez, Sequenal Mone Carlo Mehod for Mulple Targe Trang and Daa Fuon, IEEE Tranaon on Sgnal Proeng, Vol. 5, O., February. [9] B. R, S. Arulampalam, and. Gordon, Beyond he Kalman Fler, Areh Houe, 4. [] K. Kanazawa, D. Koller, and S.J. Ruell, Soha Smulaon Algorhm for Dynam Probabl ewor, In Pro. of he h Annual Conferene on Unerany n AI (UAI, Monreal, Canada, 995. [] F. Guafon, F. Gunnaron,. Bergman, U. Forell, J. Janon, R. Karlon, and P. J. ordlund, Parle Fler for Poonng, avgaon, and Trang, IEEE Tranaon on Sgnal Proeng, Vol. 5, o., February. [] D. Shulz, W. Burgard, D. Fox, and A. B. Cremer, People Trang wh a Moble Robo Ung Sample-baed Jon Probabl Daa Aoaon Fler, Inernaonal Journal of Robo Reearh (IJRR, (, 3. [3] M. S. Arulampalam, S. Maell,. Gordon, and T. Clapp, A Tuoral on Parle Fler for Onlne onlnear/on-gauan Bayean Trang, IEEE Tranaon on Sgnal Proeng, Vol. 5, o., February. [4]. Ioma,. Ihmura, T. Hguh, and H. Maeda, Parle Fler Baed Mehod for Maneuverng Targe Trang, IEEE Inernaonal Worhop on Inellgen Sgnal Proeng, Budape, Hungary, May 4-5, pp.3 8,. [5] R. R. Pre, V. P. Jlov, and X. R. L, A omparave udy of mulple model algorhm for maneuverng arge rang, Pro. 5 SPIE Conf. Sgnal Proeng, Senor Fuon, and Targe Reognon XIV, Orlando, FL, Marh 5. [6] T. E. Formann, Y. Bar-Shalom, and M. Sheffe, Sonar Trang of Mulple Targe Ung Jon Probabl Daa Aoaon, IEEE Journal of Oean Engneerng, Vol. 8, pp , 983. [7] J. Vermaa, S. J. Godll, and P. Perez, Mone Carlo Flerng for Mul-Targe Trang and Daa Aoaon, IEEE Tranaon on Aeropae and Eleron Syem, Vol. 4, o., pp , January 5. [8] O. Fran, J. eo, J. Guvan, and S. Shedng, Mulple Targe Trang Ung Sequenal Mone Carlo Mehod and Saal Daa Aoaon, Proeedng of he 3 IEEE/IRSJ, Inernaonal Conferene on Inellgen Robo and Syem, La Vega, evada, Oober 3. [9] Jao F. G. de Frea, Bayean Mehod for eural ewor, PhD he, Trny College, Unvery of Cambrdge, 999. [] P. D. Moral, A. Doue, and A. Jara, Sequenal Mone Carlo Sampler, J. R. Sa. So. B, 68, Par 3, pp , 6. Amrabr / MISC / Vol. 4 / o. / Fall 9 5

(,,, ) (,,, ). In addition, there are three other consumers, -2, -1, and 0. Consumer -2 has the utility function

(,,, ) (,,, ). In addition, there are three other consumers, -2, -1, and 0. Consumer -2 has the utility function MACROECONOMIC THEORY T J KEHOE ECON 87 SPRING 5 PROBLEM SET # Conder an overlappng generaon economy le ha n queon 5 on problem e n whch conumer lve for perod The uly funcon of he conumer born n perod,

More information

Pendulum Dynamics. = Ft tangential direction (2) radial direction (1)

Pendulum Dynamics. = Ft tangential direction (2) radial direction (1) Pendulum Dynams Consder a smple pendulum wh a massless arm of lengh L and a pon mass, m, a he end of he arm. Assumng ha he fron n he sysem s proporonal o he negave of he angenal veloy, Newon s seond law

More information

V.Abramov - FURTHER ANALYSIS OF CONFIDENCE INTERVALS FOR LARGE CLIENT/SERVER COMPUTER NETWORKS

V.Abramov - FURTHER ANALYSIS OF CONFIDENCE INTERVALS FOR LARGE CLIENT/SERVER COMPUTER NETWORKS R&RATA # Vol.) 8, March FURTHER AALYSIS OF COFIDECE ITERVALS FOR LARGE CLIET/SERVER COMPUTER ETWORKS Vyacheslav Abramov School of Mahemacal Scences, Monash Unversy, Buldng 8, Level 4, Clayon Campus, Wellngon

More information

Fall 2010 Graduate Course on Dynamic Learning

Fall 2010 Graduate Course on Dynamic Learning Fall 200 Graduae Course on Dynamc Learnng Chaper 4: Parcle Flers Sepember 27, 200 Byoung-Tak Zhang School of Compuer Scence and Engneerng & Cognve Scence and Bran Scence Programs Seoul aonal Unversy hp://b.snu.ac.kr/~bzhang/

More information

WiH Wei He

WiH Wei He Sysem Idenfcaon of onlnear Sae-Space Space Baery odels WH We He wehe@calce.umd.edu Advsor: Dr. Chaochao Chen Deparmen of echancal Engneerng Unversy of aryland, College Par 1 Unversy of aryland Bacground

More information

COMPUTER SCIENCE 349A SAMPLE EXAM QUESTIONS WITH SOLUTIONS PARTS 1, 2

COMPUTER SCIENCE 349A SAMPLE EXAM QUESTIONS WITH SOLUTIONS PARTS 1, 2 COMPUTE SCIENCE 49A SAMPLE EXAM QUESTIONS WITH SOLUTIONS PATS, PAT.. a Dene he erm ll-ondoned problem. b Gve an eample o a polynomal ha has ll-ondoned zeros.. Consder evaluaon o anh, where e e anh. e e

More information

Lecture Notes 4: Consumption 1

Lecture Notes 4: Consumption 1 Leure Noes 4: Consumpon Zhwe Xu (xuzhwe@sju.edu.n) hs noe dsusses households onsumpon hoe. In he nex leure, we wll dsuss rm s nvesmen deson. I s safe o say ha any propagaon mehansm of maroeonom model s

More information

A Demand System for Input Factors when there are Technological Changes in Production

A Demand System for Input Factors when there are Technological Changes in Production A Demand Syem for Inpu Facor when here are Technologcal Change n Producon Movaon Due o (e.g.) echnologcal change here mgh no be a aonary relaonhp for he co hare of each npu facor. When emang demand yem

More information

Electromagnetic waves in vacuum.

Electromagnetic waves in vacuum. leromagne waves n vauum. The dsovery of dsplaemen urrens enals a peular lass of soluons of Maxwell equaons: ravellng waves of eler and magne felds n vauum. In he absene of urrens and harges, he equaons

More information

Cooling of a hot metal forging. , dt dt

Cooling of a hot metal forging. , dt dt Tranen Conducon Uneady Analy - Lumped Thermal Capacy Model Performed when; Hea ranfer whn a yem produced a unform emperaure drbuon n he yem (mall emperaure graden). The emperaure change whn he yem condered

More information

In the complete model, these slopes are ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL. (! i+1 -! i ) + [(!") i+1,q - [(!

In the complete model, these slopes are ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL. (! i+1 -! i ) + [(!) i+1,q - [(! ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL The frs hng o es n wo-way ANOVA: Is here neracon? "No neracon" means: The man effecs model would f. Ths n urn means: In he neracon plo (wh A on he horzonal

More information

BIOINFORMATICS. Finding Exact Optimal Motif in Matrix Representation by Partitioning. Henry C.M. Leung, Francis Y.L. Chin and Bethany M.Y.

BIOINFORMATICS. Finding Exact Optimal Motif in Matrix Representation by Partitioning. Henry C.M. Leung, Francis Y.L. Chin and Bethany M.Y. BIOINFORMIS Vol. no. 25, page 5 do:.93/bonforma/b283 Fndng Exa Opmal Mof n Marx Repreenaon by Paronng Henry.M. Leung, Fran Y.L. hn and Behany M.Y. han Deparmen of ompuer Sene, he Unvery of Hong Kong, Pokfulam

More information

Cubic Bezier Homotopy Function for Solving Exponential Equations

Cubic Bezier Homotopy Function for Solving Exponential Equations Penerb Journal of Advanced Research n Compung and Applcaons ISSN (onlne: 46-97 Vol. 4, No.. Pages -8, 6 omoopy Funcon for Solvng Eponenal Equaons S. S. Raml *,,. Mohamad Nor,a, N. S. Saharzan,b and M.

More information

Online Appendix for. Strategic safety stocks in supply chains with evolving forecasts

Online Appendix for. Strategic safety stocks in supply chains with evolving forecasts Onlne Appendx for Sraegc safey socs n supply chans wh evolvng forecass Tor Schoenmeyr Sephen C. Graves Opsolar, Inc. 332 Hunwood Avenue Hayward, CA 94544 A. P. Sloan School of Managemen Massachuses Insue

More information

CS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 4

CS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 4 CS434a/54a: Paern Recognon Prof. Olga Veksler Lecure 4 Oulne Normal Random Varable Properes Dscrmnan funcons Why Normal Random Varables? Analycally racable Works well when observaon comes form a corruped

More information

HEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD

HEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD Journal of Appled Mahemacs and Compuaonal Mechancs 3, (), 45-5 HEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD Sansław Kukla, Urszula Sedlecka Insue of Mahemacs,

More information

)-interval valued fuzzy ideals in BF-algebras. Some properties of (, ) -interval valued fuzzy ideals in BF-algebra, where

)-interval valued fuzzy ideals in BF-algebras. Some properties of (, ) -interval valued fuzzy ideals in BF-algebra, where Inernaonal Journal of Engneerng Advaned Researh Tehnology (IJEART) ISSN: 454-990, Volume-, Issue-4, Oober 05 Some properes of (, )-nerval valued fuzzy deals n BF-algebras M. Idrees, A. Rehman, M. Zulfqar,

More information

THE PREDICTION OF COMPETITIVE ENVIRONMENT IN BUSINESS

THE PREDICTION OF COMPETITIVE ENVIRONMENT IN BUSINESS THE PREICTION OF COMPETITIVE ENVIRONMENT IN BUSINESS INTROUCTION The wo dmensonal paral dfferenal equaons of second order can be used for he smulaon of compeve envronmen n busness The arcle presens he

More information

Simultaneous Localization of Mobile Robot and Multiple Sound Sources Using Microphone Array

Simultaneous Localization of Mobile Robot and Multiple Sound Sources Using Microphone Array Smulaneou Localzaon o Moble Robo and Mulple Sound Source Ung Mcrophone Array Jwu-Sheng Hu, Chen-Yu Chan, Cheng-Kang Wang Deparmen o Elecrcal and Conrol Engneerng aonal Chao-ung Unvery Hnchu, awan, R.O.C.

More information

Fundamentals of PLLs (I)

Fundamentals of PLLs (I) Phae-Locked Loop Fundamenal of PLL (I) Chng-Yuan Yang Naonal Chung-Hng Unvery Deparmen of Elecrcal Engneerng Why phae-lock? - Jer Supreon - Frequency Synhe T T + 1 - Skew Reducon T + 2 T + 3 PLL fou =

More information

WALKING PLANNING AND CONTROL FOR A BIPED ROBOT UPSTAIRS

WALKING PLANNING AND CONTROL FOR A BIPED ROBOT UPSTAIRS WALKING PLANNING AND CONTROL FOR A BIPED ROBOT UPSTAIRS Cenbo Yn, Dongua Zeng and Le Xao 2 Sool o Meanal and Power Engneerng, Nanjng Unvery o Tnology, Nanjng, Cna ynenbo@nju.edu.n, engdongua98@63.om 2

More information

Foundations of State Estimation Part II

Foundations of State Estimation Part II Foundaons of Sae Esmaon Par II Tocs: Hdden Markov Models Parcle Flers Addonal readng: L.R. Rabner, A uoral on hdden Markov models," Proceedngs of he IEEE, vol. 77,. 57-86, 989. Sequenal Mone Carlo Mehods

More information

John Geweke a and Gianni Amisano b a Departments of Economics and Statistics, University of Iowa, USA b European Central Bank, Frankfurt, Germany

John Geweke a and Gianni Amisano b a Departments of Economics and Statistics, University of Iowa, USA b European Central Bank, Frankfurt, Germany Herarchcal Markov Normal Mxure models wh Applcaons o Fnancal Asse Reurns Appendx: Proofs of Theorems and Condonal Poseror Dsrbuons John Geweke a and Gann Amsano b a Deparmens of Economcs and Sascs, Unversy

More information

Variants of Pegasos. December 11, 2009

Variants of Pegasos. December 11, 2009 Inroducon Varans of Pegasos SooWoong Ryu bshboy@sanford.edu December, 009 Youngsoo Cho yc344@sanford.edu Developng a new SVM algorhm s ongong research opc. Among many exng SVM algorhms, we wll focus on

More information

( ) () we define the interaction representation by the unitary transformation () = ()

( ) () we define the interaction representation by the unitary transformation () = () Hgher Order Perurbaon Theory Mchael Fowler 3/7/6 The neracon Represenaon Recall ha n he frs par of hs course sequence, we dscussed he chrödnger and Hesenberg represenaons of quanum mechancs here n he chrödnger

More information

SSRG International Journal of Thermal Engineering (SSRG-IJTE) Volume 4 Issue 1 January to April 2018

SSRG International Journal of Thermal Engineering (SSRG-IJTE) Volume 4 Issue 1 January to April 2018 SSRG Inernaonal Journal of Thermal Engneerng (SSRG-IJTE) Volume 4 Iue 1 January o Aprl 18 Opmal Conrol for a Drbued Parameer Syem wh Tme-Delay, Non-Lnear Ung he Numercal Mehod. Applcaon o One- Sded Hea

More information

FTCS Solution to the Heat Equation

FTCS Solution to the Heat Equation FTCS Soluon o he Hea Equaon ME 448/548 Noes Gerald Reckenwald Porland Sae Unversy Deparmen of Mechancal Engneerng gerry@pdxedu ME 448/548: FTCS Soluon o he Hea Equaon Overvew Use he forward fne d erence

More information

[ ] 2. [ ]3 + (Δx i + Δx i 1 ) / 2. Δx i-1 Δx i Δx i+1. TPG4160 Reservoir Simulation 2018 Lecture note 3. page 1 of 5

[ ] 2. [ ]3 + (Δx i + Δx i 1 ) / 2. Δx i-1 Δx i Δx i+1. TPG4160 Reservoir Simulation 2018 Lecture note 3. page 1 of 5 TPG460 Reservor Smulaon 08 page of 5 DISCRETIZATIO OF THE FOW EQUATIOS As we already have seen, fne dfference appromaons of he paral dervaves appearng n he flow equaons may be obaned from Taylor seres

More information

A Dynamic Algorithm for Longest Common Subsequence Problem using Ant Colony Optimization Technique

A Dynamic Algorithm for Longest Common Subsequence Problem using Ant Colony Optimization Technique A Dynam Algorhm or Longe Common Subequene roblem ung An Colony Opmzaon ehnque Arndam Chaudhur Leurer Mahema and Compuer Sene Meghnad Saha Inue o ehnology Nazrabad Uhhepoa Kolaa Inda Emal: arndam_hau@yahoo.o.n

More information

CTLS 4 SNR. Multi Reference CTLS Method for Passive Localization of Radar Targets

CTLS 4 SNR. Multi Reference CTLS Method for Passive Localization of Radar Targets دا ند رعا ل» ی و ناوری ج ه ع ی و ی «ع وم 79-85 9 C 4 * Donloaded from ad.r a 9:06 +040 on Frda arch nd 09-4 - - - - (9/06/4 : 90/05/7 : ) DOA. DOA. C DOA.. C.. C SR.. C.C DOA : ul Reference C ehod for

More information

Robustness Experiments with Two Variance Components

Robustness Experiments with Two Variance Components Naonal Insue of Sandards and Technology (NIST) Informaon Technology Laboraory (ITL) Sascal Engneerng Dvson (SED) Robusness Expermens wh Two Varance Componens by Ana Ivelsse Avlés avles@ns.gov Conference

More information

Filtrage particulaire et suivi multi-pistes Carine Hue Jean-Pierre Le Cadre and Patrick Pérez

Filtrage particulaire et suivi multi-pistes Carine Hue Jean-Pierre Le Cadre and Patrick Pérez Chaînes de Markov cachées e flrage parculare 2-22 anver 2002 Flrage parculare e suv mul-pses Carne Hue Jean-Perre Le Cadre and Parck Pérez Conex Applcaons: Sgnal processng: arge rackng bearngs-onl rackng

More information

EECE 301 Signals & Systems Prof. Mark Fowler

EECE 301 Signals & Systems Prof. Mark Fowler EECE 31 Signal & Syem Prof. Mark Fowler Noe Se #27 C-T Syem: Laplace Tranform Power Tool for yem analyi Reading Aignmen: Secion 6.1 6.3 of Kamen and Heck 1/18 Coure Flow Diagram The arrow here how concepual

More information

Ordinary Differential Equations in Neuroscience with Matlab examples. Aim 1- Gain understanding of how to set up and solve ODE s

Ordinary Differential Equations in Neuroscience with Matlab examples. Aim 1- Gain understanding of how to set up and solve ODE s Ordnary Dfferenal Equaons n Neuroscence wh Malab eamples. Am - Gan undersandng of how o se up and solve ODE s Am Undersand how o se up an solve a smple eample of he Hebb rule n D Our goal a end of class

More information

Solution in semi infinite diffusion couples (error function analysis)

Solution in semi infinite diffusion couples (error function analysis) Soluon n sem nfne dffuson couples (error funcon analyss) Le us consder now he sem nfne dffuson couple of wo blocks wh concenraon of and I means ha, n a A- bnary sysem, s bondng beween wo blocks made of

More information

Computational results on new staff scheduling benchmark instances

Computational results on new staff scheduling benchmark instances TECHNICAL REPORT Compuaonal resuls on new saff shedulng enhmark nsanes Tm Curos Rong Qu ASAP Researh Group Shool of Compuer Sene Unersy of Nongham NG8 1BB Nongham UK Frs pulshed onlne: 19-Sep-2014 las

More information

Matrix reconstruction with the local max norm

Matrix reconstruction with the local max norm Marx reconrucon wh he local max norm Rna oygel Deparmen of Sac Sanford Unvery rnafb@anfordedu Nahan Srebro Toyoa Technologcal Inue a Chcago na@cedu Rulan Salakhudnov Dep of Sac and Dep of Compuer Scence

More information

8.1. a) For step response, M input is u ( t) Taking inverse Laplace transform. as α 0. Ideal response, K c. = Kc Mτ D + For ramp response, 8-1

8.1. a) For step response, M input is u ( t) Taking inverse Laplace transform. as α 0. Ideal response, K c. = Kc Mτ D + For ramp response, 8-1 8. a For ep repone, inpu i u, U Y a U α α Y a α α Taking invere Laplae ranform a α e e / α / α A α 0 a δ 0 e / α a δ deal repone, α d Y i Gi U i δ Hene a α 0 a i For ramp repone, inpu i u, U Soluion anual

More information

Computer Robot Vision Conference 2010

Computer Robot Vision Conference 2010 School of Compuer Scence McGll Unversy Compuer Robo Vson Conference 2010 Ioanns Rekles Fundamenal Problems In Robocs How o Go From A o B? (Pah Plannng) Wha does he world looks lke? (mappng) sense from

More information

Inferring Human Upper Body Motion

Inferring Human Upper Body Motion Inferrng Human Upper Body Moon Jang Gao, Janbo h Roboc Inue C arnege Mellon Unvery Abrac We preen a new algorhm for auomac nference of human upper body moon n a naural cene. The nal moon cue are fr deeced

More information

Multiple Failures. Diverse Routing for Maximizing Survivability. Maximum Survivability Models. Minimum-Color (SRLG) Diverse Routing

Multiple Failures. Diverse Routing for Maximizing Survivability. Maximum Survivability Models. Minimum-Color (SRLG) Diverse Routing Mulple Falure Dvere Roung for Maxmzng Survvably One-falure aumpon n prevou work Mulple falure Hard o provde 100% proecon Maxmum urvvably Maxmum Survvably Model Mnmum-Color (SRLG) Dvere Roung Each lnk ha

More information

A. Inventory model. Why are we interested in it? What do we really study in such cases.

A. Inventory model. Why are we interested in it? What do we really study in such cases. Some general yem model.. Inenory model. Why are we nereed n? Wha do we really udy n uch cae. General raegy of machng wo dmlar procee, ay, machng a fa proce wh a low one. We need an nenory or a buffer or

More information

Multiple Regressions and Correlation Analysis

Multiple Regressions and Correlation Analysis Mulple Regreon and Correlaon Analy Chaper 4 McGraw-Hll/Irwn Copyrgh 2 y The McGraw-Hll Compane, Inc. All rgh reerved. GOALS. Decre he relaonhp eween everal ndependen varale and a dependen varale ung mulple

More information

Normal Random Variable and its discriminant functions

Normal Random Variable and its discriminant functions Noral Rando Varable and s dscrnan funcons Oulne Noral Rando Varable Properes Dscrnan funcons Why Noral Rando Varables? Analycally racable Works well when observaon coes for a corruped snle prooype 3 The

More information

CS286.2 Lecture 14: Quantum de Finetti Theorems II

CS286.2 Lecture 14: Quantum de Finetti Theorems II CS286.2 Lecure 14: Quanum de Fne Theorems II Scrbe: Mara Okounkova 1 Saemen of he heorem Recall he las saemen of he quanum de Fne heorem from he prevous lecure. Theorem 1 Quanum de Fne). Le ρ Dens C 2

More information

u(t) Figure 1. Open loop control system

u(t) Figure 1. Open loop control system Open loop conrol v cloed loop feedbac conrol The nex wo figure preen he rucure of open loop and feedbac conrol yem Figure how an open loop conrol yem whoe funcion i o caue he oupu y o follow he reference

More information

ANALYSIS AND MODELING OF HYDROLOGIC TIME SERIES. Wasserhaushalt Time Series Analysis and Stochastic Modelling Spring Semester

ANALYSIS AND MODELING OF HYDROLOGIC TIME SERIES. Wasserhaushalt Time Series Analysis and Stochastic Modelling Spring Semester ANALYSIS AND MODELING OF HYDROLOGIC TIME SERIES Waerhauhal Tme Sere Analy and Sochac Modellng Sprng Semeer 8 ANALYSIS AND MODELING OF HYDROLOGIC TIME SERIES Defnon Wha a me ere? Leraure: Sala, J.D. 99,

More information

Lecture 18: The Laplace Transform (See Sections and 14.7 in Boas)

Lecture 18: The Laplace Transform (See Sections and 14.7 in Boas) Lecure 8: The Lalace Transform (See Secons 88- and 47 n Boas) Recall ha our bg-cure goal s he analyss of he dfferenal equaon, ax bx cx F, where we emloy varous exansons for he drvng funcon F deendng on

More information

Lecture 11 SVM cont

Lecture 11 SVM cont Lecure SVM con. 0 008 Wha we have done so far We have esalshed ha we wan o fnd a lnear decson oundary whose margn s he larges We know how o measure he margn of a lnear decson oundary Tha s: he mnmum geomerc

More information

Mechanics Physics 151

Mechanics Physics 151 Mechancs Physcs 5 Lecure 9 Hamlonan Equaons of Moon (Chaper 8) Wha We Dd Las Tme Consruced Hamlonan formalsm H ( q, p, ) = q p L( q, q, ) H p = q H q = p H = L Equvalen o Lagrangan formalsm Smpler, bu

More information

On One Analytic Method of. Constructing Program Controls

On One Analytic Method of. Constructing Program Controls Appled Mahemacal Scences, Vol. 9, 05, no. 8, 409-407 HIKARI Ld, www.m-hkar.com hp://dx.do.org/0.988/ams.05.54349 On One Analyc Mehod of Consrucng Program Conrols A. N. Kvko, S. V. Chsyakov and Yu. E. Balyna

More information

Mechanics Physics 151

Mechanics Physics 151 Mechancs Physcs 5 Lecure 9 Hamlonan Equaons of Moon (Chaper 8) Wha We Dd Las Tme Consruced Hamlonan formalsm Hqp (,,) = qp Lqq (,,) H p = q H q = p H L = Equvalen o Lagrangan formalsm Smpler, bu wce as

More information

GENERATING CERTAIN QUINTIC IRREDUCIBLE POLYNOMIALS OVER FINITE FIELDS. Youngwoo Ahn and Kitae Kim

GENERATING CERTAIN QUINTIC IRREDUCIBLE POLYNOMIALS OVER FINITE FIELDS. Youngwoo Ahn and Kitae Kim Korean J. Mah. 19 (2011), No. 3, pp. 263 272 GENERATING CERTAIN QUINTIC IRREDUCIBLE POLYNOMIALS OVER FINITE FIELDS Youngwoo Ahn and Kae Km Absrac. In he paper [1], an explc correspondence beween ceran

More information

Volatility Interpolation

Volatility Interpolation Volaly Inerpolaon Prelmnary Verson March 00 Jesper Andreasen and Bran Huge Danse Mares, Copenhagen wan.daddy@danseban.com brno@danseban.com Elecronc copy avalable a: hp://ssrn.com/absrac=69497 Inro Local

More information

ELIMINATION OF DOMINATED STRATEGIES AND INESSENTIAL PLAYERS

ELIMINATION OF DOMINATED STRATEGIES AND INESSENTIAL PLAYERS OPERATIONS RESEARCH AND DECISIONS No. 1 215 DOI: 1.5277/ord1513 Mamoru KANEKO 1 Shuge LIU 1 ELIMINATION OF DOMINATED STRATEGIES AND INESSENTIAL PLAYERS We udy he proce, called he IEDI proce, of eraed elmnaon

More information

Chapter 6: AC Circuits

Chapter 6: AC Circuits Chaper 6: AC Crcus Chaper 6: Oulne Phasors and he AC Seady Sae AC Crcus A sable, lnear crcu operang n he seady sae wh snusodal excaon (.e., snusodal seady sae. Complee response forced response naural response.

More information

J i-1 i. J i i+1. Numerical integration of the diffusion equation (I) Finite difference method. Spatial Discretization. Internal nodes.

J i-1 i. J i i+1. Numerical integration of the diffusion equation (I) Finite difference method. Spatial Discretization. Internal nodes. umercal negraon of he dffuson equaon (I) Fne dfference mehod. Spaal screaon. Inernal nodes. R L V For hermal conducon le s dscree he spaal doman no small fne spans, =,,: Balance of parcles for an nernal

More information

H = d d q 1 d d q N d d p 1 d d p N exp

H = d d q 1 d d q N d d p 1 d d p N exp 8333: Sacal Mechanc I roblem Se # 7 Soluon Fall 3 Canoncal Enemble Non-harmonc Ga: The Hamlonan for a ga of N non neracng parcle n a d dmenonal box ha he form H A p a The paron funcon gven by ZN T d d

More information

Robust and Accurate Cancer Classification with Gene Expression Profiling

Robust and Accurate Cancer Classification with Gene Expression Profiling Robus and Accurae Cancer Classfcaon wh Gene Expresson Proflng (Compuaonal ysems Bology, 2005) Auhor: Hafeng L, Keshu Zhang, ao Jang Oulne Background LDA (lnear dscrmnan analyss) and small sample sze problem

More information

Single-loop System Reliability-Based Design & Topology Optimization (SRBDO/SRBTO): A Matrix-based System Reliability (MSR) Method

Single-loop System Reliability-Based Design & Topology Optimization (SRBDO/SRBTO): A Matrix-based System Reliability (MSR) Method 10 h US Naonal Congress on Compuaonal Mechancs Columbus, Oho 16-19, 2009 Sngle-loop Sysem Relably-Based Desgn & Topology Opmzaon (SRBDO/SRBTO): A Marx-based Sysem Relably (MSR) Mehod Tam Nguyen, Junho

More information

Department of Economics University of Toronto

Department of Economics University of Toronto Deparmen of Economcs Unversy of Torono ECO408F M.A. Economercs Lecure Noes on Heeroskedascy Heeroskedascy o Ths lecure nvolves lookng a modfcaons we need o make o deal wh he regresson model when some of

More information

Appendix H: Rarefaction and extrapolation of Hill numbers for incidence data

Appendix H: Rarefaction and extrapolation of Hill numbers for incidence data Anne Chao Ncholas J Goell C seh lzabeh L ander K Ma Rober K Colwell and Aaron M llson 03 Rarefacon and erapolaon wh ll numbers: a framewor for samplng and esmaon n speces dversy sudes cology Monographs

More information

Performance Comparison of Bivariate Copulas on the CUSUM and EWMA Control Charts

Performance Comparison of Bivariate Copulas on the CUSUM and EWMA Control Charts Proeedngs of he World Congress on Engneerng and Compuer Sene 5 Vol II WCECS 5, Oober -3, 5, San Franso, USA Performane Comparson of Bvarae Copulas on he CUSUM and EWMA Conrol Chars Sasgarn Kuvaana, Saowan

More information

Approximate Analytic Solution of (2+1) - Dimensional Zakharov-Kuznetsov(Zk) Equations Using Homotopy

Approximate Analytic Solution of (2+1) - Dimensional Zakharov-Kuznetsov(Zk) Equations Using Homotopy Arcle Inernaonal Journal of Modern Mahemacal Scences, 4, (): - Inernaonal Journal of Modern Mahemacal Scences Journal homepage: www.modernscenfcpress.com/journals/jmms.aspx ISSN: 66-86X Florda, USA Approxmae

More information

Advanced time-series analysis (University of Lund, Economic History Department)

Advanced time-series analysis (University of Lund, Economic History Department) Advanced me-seres analss (Unvers of Lund, Economc Hsor Dearmen) 3 Jan-3 Februar and 6-3 March Lecure 4 Economerc echnues for saonar seres : Unvarae sochasc models wh Box- Jenns mehodolog, smle forecasng

More information

Solution. The straightforward approach is surprisingly difficult because one has to be careful about the limits.

Solution. The straightforward approach is surprisingly difficult because one has to be careful about the limits. ose ad Varably Homewor # (8), aswers Q: Power spera of some smple oses A Posso ose A Posso ose () s a sequee of dela-fuo pulses, eah ourrg depedely, a some rae r (More formally, s a sum of pulses of wdh

More information

Existence and Uniqueness Results for Random Impulsive Integro-Differential Equation

Existence and Uniqueness Results for Random Impulsive Integro-Differential Equation Global Journal of Pure and Appled Mahemacs. ISSN 973-768 Volume 4, Number 6 (8), pp. 89-87 Research Inda Publcaons hp://www.rpublcaon.com Exsence and Unqueness Resuls for Random Impulsve Inegro-Dfferenal

More information

Development of a New Optimal Multilevel Thresholding Using Improved Particle Swarm Optimization Algorithm for Image Segmentation

Development of a New Optimal Multilevel Thresholding Using Improved Particle Swarm Optimization Algorithm for Image Segmentation Inernaonal Journal of Elerons Engneerng, (1), 010, pp. 63-67 Developmen of a New Opmal Mullevel Thresholdng Usng Improved Parle Swarm Opmzaon Algorhm for Image Segmenaon P.D. Sahya 1 & R. Kayalvzh 1 Deparmen

More information

Control Systems. Mathematical Modeling of Control Systems.

Control Systems. Mathematical Modeling of Control Systems. Conrol Syem Mahemacal Modelng of Conrol Syem chbum@eoulech.ac.kr Oulne Mahemacal model and model ype. Tranfer funcon model Syem pole and zero Chbum Lee -Seoulech Conrol Syem Mahemacal Model Model are key

More information

The Maxwell equations as a Bäcklund transformation

The Maxwell equations as a Bäcklund transformation ADVANCED ELECTROMAGNETICS, VOL. 4, NO. 1, JULY 15 The Mawell equaons as a Bäklund ransformaon C. J. Papahrsou Deparmen of Physal Senes, Naval Aademy of Greee, Praeus, Greee papahrsou@snd.edu.gr Absra Bäklund

More information

Performance Analysis for a Network having Standby Redundant Unit with Waiting in Repair

Performance Analysis for a Network having Standby Redundant Unit with Waiting in Repair TECHNI Inernaonal Journal of Compung Scence Communcaon Technologes VOL.5 NO. July 22 (ISSN 974-3375 erformance nalyss for a Nework havng Sby edundan Un wh ang n epar Jendra Sngh 2 abns orwal 2 Deparmen

More information

Lecture 11: Stereo and Surface Estimation

Lecture 11: Stereo and Surface Estimation Lecure : Sereo and Surface Emaon When camera poon have been deermned, ung rucure from moon, we would lke o compue a dene urface model of he cene. In h lecure we wll udy he o called Sereo Problem, where

More information

THERMODYNAMICS 1. The First Law and Other Basic Concepts (part 2)

THERMODYNAMICS 1. The First Law and Other Basic Concepts (part 2) Company LOGO THERMODYNAMICS The Frs Law and Oher Basc Conceps (par ) Deparmen of Chemcal Engneerng, Semarang Sae Unversy Dhon Harano S.T., M.T., M.Sc. Have you ever cooked? Equlbrum Equlbrum (con.) Equlbrum

More information

XMAP: Track-to-Track Association with Metric, Feature, and Target-type Data

XMAP: Track-to-Track Association with Metric, Feature, and Target-type Data XMAP: Track-o-Track Aocaon wh Merc, Feaure, Targe-ype Daa J. Ferry Meron, Inc. Reon, VA, U.S.A. ferry@mec.com Abrac - The Exended Maxmum A Poeror Probably XMAP mehod for rack-o-rack aocaon baed on a formal,

More information

Randomized Perfect Bipartite Matching

Randomized Perfect Bipartite Matching Inenive Algorihm Lecure 24 Randomized Perfec Biparie Maching Lecurer: Daniel A. Spielman April 9, 208 24. Inroducion We explain a randomized algorihm by Ahih Goel, Michael Kapralov and Sanjeev Khanna for

More information

Chapter 5 Signal-Space Analysis

Chapter 5 Signal-Space Analysis Chaper 5 Sgnal-Space Analy Sgnal pace analy provde a mahemacally elegan and hghly nghful ool for he udy of daa ranmon. 5. Inroducon o Sacal model for a genec dgal communcaon yem n eage ource: A pror probable

More information

On computing differential transform of nonlinear non-autonomous functions and its applications

On computing differential transform of nonlinear non-autonomous functions and its applications On compung dfferenal ransform of nonlnear non-auonomous funcons and s applcaons Essam. R. El-Zahar, and Abdelhalm Ebad Deparmen of Mahemacs, Faculy of Scences and Humanes, Prnce Saam Bn Abdulazz Unversy,

More information

Advances in Applied Acoustics (AIAA) Volume 2 Issue 1, February 2013

Advances in Applied Acoustics (AIAA) Volume 2 Issue 1, February 2013 Avane n Apple Aou AIAA) Volume Iue, Februar 0 wwwaaa-ournalorg Moellng Daa Dperon Degree n Auoma Robu Emaon for Mulvarae Gauan Mxure Moel wh an Applaon o No Speeh Proeng Dale Wu *, Haqng Wu Deparmen of

More information

e-journal Reliability: Theory& Applications No 2 (Vol.2) Vyacheslav Abramov

e-journal Reliability: Theory& Applications No 2 (Vol.2) Vyacheslav Abramov June 7 e-ournal Relably: Theory& Applcaons No (Vol. CONFIDENCE INTERVALS ASSOCIATED WITH PERFORMANCE ANALYSIS OF SYMMETRIC LARGE CLOSED CLIENT/SERVER COMPUTER NETWORKS Absrac Vyacheslav Abramov School

More information

NON-HOMOGENEOUS SEMI-MARKOV REWARD PROCESS FOR THE MANAGEMENT OF HEALTH INSURANCE MODELS.

NON-HOMOGENEOUS SEMI-MARKOV REWARD PROCESS FOR THE MANAGEMENT OF HEALTH INSURANCE MODELS. NON-HOOGENEOU EI-AKO EWA POCE FO THE ANAGEENT OF HEATH INUANCE OE. Jacque Janen CEIAF ld Paul Janon 84 e 9 6 Charlero EGIU Fax: 32735877 E-mal: ceaf@elgacom.ne and amondo anca Unverà a apenza parmeno d

More information

5.2 GRAPHICAL VELOCITY ANALYSIS Polygon Method

5.2 GRAPHICAL VELOCITY ANALYSIS Polygon Method ME 352 GRHICL VELCITY NLYSIS 52 GRHICL VELCITY NLYSIS olygon Mehod Velociy analyi form he hear of kinemaic and dynamic of mechanical yem Velociy analyi i uually performed following a poiion analyi; ie,

More information

Methods of Improving Constitutive Equations

Methods of Improving Constitutive Equations Mehods o mprovng Consuve Equaons Maxell Model e an mprove h ne me dervaves or ne sran measures. ³ ª º «e, d» ¼ e an also hange he bas equaon lnear modaons non-lnear modaons her Consuve Approahes Smple

More information

Method of Characteristics for Pure Advection By Gilberto E. Urroz, September 2004

Method of Characteristics for Pure Advection By Gilberto E. Urroz, September 2004 Mehod of Charaerss for Pre Adveon By Glbero E Urroz Sepember 004 Noe: The followng noes are based on lass noes for he lass COMPUTATIONAL HYDAULICS as agh by Dr Forres Holly n he Sprng Semeser 985 a he

More information

Discussion Session 2 Constant Acceleration/Relative Motion Week 03

Discussion Session 2 Constant Acceleration/Relative Motion Week 03 PHYS 100 Dicuion Seion Conan Acceleraion/Relaive Moion Week 03 The Plan Today you will work wih your group explore he idea of reference frame (i.e. relaive moion) and moion wih conan acceleraion. You ll

More information

Dynamic Team Decision Theory. EECS 558 Project Shrutivandana Sharma and David Shuman December 10, 2005

Dynamic Team Decision Theory. EECS 558 Project Shrutivandana Sharma and David Shuman December 10, 2005 Dynamc Team Decson Theory EECS 558 Proec Shruvandana Sharma and Davd Shuman December 0, 005 Oulne Inroducon o Team Decson Theory Decomposon of he Dynamc Team Decson Problem Equvalence of Sac and Dynamc

More information

Graduate Macroeconomics 2 Problem set 5. - Solutions

Graduate Macroeconomics 2 Problem set 5. - Solutions Graduae Macroeconomcs 2 Problem se. - Soluons Queson 1 To answer hs queson we need he frms frs order condons and he equaon ha deermnes he number of frms n equlbrum. The frms frs order condons are: F K

More information

Gradient Flow Independent Component Analysis

Gradient Flow Independent Component Analysis Graden Fow Independen Componen Anay Mun Sanaćevć and Ger Cauwenbergh Adapve Mcroyem ab John Hopkn Unvery Oune Bnd Sgna Separaon and ocazaon Prncpe of Graden Fow : from deay o empora dervave Equvaen ac

More information

Introduction ( Week 1-2) Course introduction A brief introduction to molecular biology A brief introduction to sequence comparison Part I: Algorithms

Introduction ( Week 1-2) Course introduction A brief introduction to molecular biology A brief introduction to sequence comparison Part I: Algorithms Course organzaon Inroducon Wee -2) Course nroducon A bref nroducon o molecular bology A bref nroducon o sequence comparson Par I: Algorhms for Sequence Analyss Wee 3-8) Chaper -3, Models and heores» Probably

More information

Pen Tip Position Estimation Using Least Square Sphere Fitting for Customized Attachments of Haptic Device

Pen Tip Position Estimation Using Least Square Sphere Fitting for Customized Attachments of Haptic Device for Cuomed Ahmen of Hp Deve Mno KOEDA nd Mhko KAO Deprmen of Compuer Sene Ful of Informon Sene nd Ar Ok Elero-Communon Unver Kok 30-70, Shjonwe, Ok, 575-0063, JAPA {koed, 0809@oeu.jp} Ar In h pper, mehod

More information

Online EM Algorithm for Background Subtraction

Online EM Algorithm for Background Subtraction Avalable onlne a www.cencedrec.com Proceda Engneerng 9 (0) 64 69 0 Inernaonal Workhop on Informaon and Elecronc Engneerng (IWIEE) Onlne E Algorhm for Background Subracon Peng Chen a*, Xang Chen b,bebe

More information

Discrete Markov Process. Introduction. Example: Balls and Urns. Stochastic Automaton. INTRODUCTION TO Machine Learning 3rd Edition

Discrete Markov Process. Introduction. Example: Balls and Urns. Stochastic Automaton. INTRODUCTION TO Machine Learning 3rd Edition EHEM ALPAYDI he MI Press, 04 Lecure Sldes for IRODUCIO O Machne Learnng 3rd Edon alpaydn@boun.edu.r hp://www.cmpe.boun.edu.r/~ehem/ml3e Sldes from exboo resource page. Slghly eded and wh addonal examples

More information

TSS = SST + SSE An orthogonal partition of the total SS

TSS = SST + SSE An orthogonal partition of the total SS ANOVA: Topc 4. Orhogonal conrass [ST&D p. 183] H 0 : µ 1 = µ =... = µ H 1 : The mean of a leas one reamen group s dfferen To es hs hypohess, a basc ANOVA allocaes he varaon among reamen means (SST) equally

More information

Optimal Filtering for Linear Discrete-Time Systems with Single Delayed Measurement

Optimal Filtering for Linear Discrete-Time Systems with Single Delayed Measurement 378 Hong-Guo Inernaonal Zhao, Journal Huan-Shu of Conrol, Zhang, Auomaon, Cheng-Hu an Zhang, Syem, an vol. Xn-Mn 6, no. Song 3, pp. 378-385, June 28 Opmal Flerng for Lnear Dcree-me Syem h Sngle Delaye

More information

Outline. Probabilistic Model Learning. Probabilistic Model Learning. Probabilistic Model for Time-series Data: Hidden Markov Model

Outline. Probabilistic Model Learning. Probabilistic Model Learning. Probabilistic Model for Time-series Data: Hidden Markov Model Probablsc Model for Tme-seres Daa: Hdden Markov Model Hrosh Mamsuka Bonformacs Cener Kyoo Unversy Oulne Three Problems for probablsc models n machne learnng. Compung lkelhood 2. Learnng 3. Parsng (predcon

More information

Lecture 2 M/G/1 queues. M/G/1-queue

Lecture 2 M/G/1 queues. M/G/1-queue Lecure M/G/ queues M/G/-queue Posson arrval process Arbrary servce me dsrbuon Sngle server To deermne he sae of he sysem a me, we mus now The number of cusomers n he sysems N() Tme ha he cusomer currenly

More information

Should Exact Index Numbers have Standard Errors? Theory and Application to Asian Growth

Should Exact Index Numbers have Standard Errors? Theory and Application to Asian Growth Should Exac Index umbers have Sandard Errors? Theory and Applcaon o Asan Growh Rober C. Feensra Marshall B. Rensdorf ovember 003 Proof of Proposon APPEDIX () Frs, we wll derve he convenonal Sao-Vara prce

More information

Math 128b Project. Jude Yuen

Math 128b Project. Jude Yuen Mah 8b Proec Jude Yuen . Inroducon Le { Z } be a sequence of observed ndependen vecor varables. If he elemens of Z have a on normal dsrbuon hen { Z } has a mean vecor Z and a varancecovarance marx z. Geomercally

More information

Testing a new idea to solve the P = NP problem with mathematical induction

Testing a new idea to solve the P = NP problem with mathematical induction Tesng a new dea o solve he P = NP problem wh mahemacal nducon Bacground P and NP are wo classes (ses) of languages n Compuer Scence An open problem s wheher P = NP Ths paper ess a new dea o compare he

More information

Laplace Transformation of Linear Time-Varying Systems

Laplace Transformation of Linear Time-Varying Systems Laplace Tranformaon of Lnear Tme-Varyng Syem Shervn Erfan Reearch Cenre for Inegraed Mcroelecronc Elecrcal and Compuer Engneerng Deparmen Unvery of Wndor Wndor, Onaro N9B 3P4, Canada Aug. 4, 9 Oulne of

More information

Laplace Transform. Inverse Laplace Transform. e st f(t)dt. (2)

Laplace Transform. Inverse Laplace Transform. e st f(t)dt. (2) Laplace Tranform Maoud Malek The Laplace ranform i an inegral ranform named in honor of mahemaician and aronomer Pierre-Simon Laplace, who ued he ranform in hi work on probabiliy heory. I i a powerful

More information