{"id":201,"date":"2012-08-08T21:25:31","date_gmt":"2012-08-08T21:25:31","guid":{"rendered":"http:\/\/www-sop.inria.fr\/members\/Fabien.Campillo\/?page_id=201"},"modified":"2024-10-06T17:47:58","modified_gmt":"2024-10-06T17:47:58","slug":"smc-demos","status":"publish","type":"page","link":"http:\/\/localhost:8888\/wordpress\/software\/smc-demos\/","title":{"rendered":"Sequential Monte Carlo (particle filtering) for tracking"},"content":{"rendered":"\n\n\n<p>A set of <code>Matlab<\/code> software demonstrations of particle filtering for tracking applications, developed in the context of lectures given from 2005 to 2013:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>SAGEM, Argenteuil \u2014 12\u201313 juin, 3\u20134 juillet, 19 septembre. 2003 [<a href=\"http:\/\/localhost:8888\/wordpress\/wp-content\/uploads\/pdf\/campillo2003a.pdf\" data-type=\"attachment\" data-id=\"2707\">slides<\/a>],<\/li>\n\n\n\n<li>Centre de Comp\u00e9tence Technique (CCT) et Centre National d&#8217;\\&#8217;Etudes Spatiales (CNES) Toulouse, Novembre 24, 2003,<\/li>\n\n\n\n<li>\u00c9cole Chercheurs en \u201ctraitement du signal\u201d, Inria Rennes, 18 au 20 octobre 2004,<\/li>\n\n\n\n<li>IAP DYSCO Study Day, Dynamical systems, control and optimization \u2013 University of Louvain, Mons. Friday May 24, 2013 [<a href=\"http:\/\/localhost:8888\/wordpress\/wp-content\/uploads\/pdf\/campillo2013e.pdf\" data-type=\"attachment\" data-id=\"1943\">slides<\/a>],<\/li>\n<\/ul>\n\n\n\n<p>and in many seminars and in doctoral training programs. I also used these <code>Matlab<\/code> software demonstrations in a course I taught in Toulon from 2004 to 2008 as part of the Master of Science and Technology at the University of the South Toulon-Var [<a href=\"http:\/\/localhost:8888\/wordpress\/wp-content\/uploads\/pdf\/premiers_pas_en_filtrage_particulaire.pdf\">pdf<\/a>, the 2006 version if in the CEL website <a href=\"http:\/\/cel.archives-ouvertes.fr\/inria-00103752\">here<\/a>].<\/p>\n\n\n\n<p>See\u00a0<a href=\"https:\/\/gitlab.inria.fr\/campillo\/smc-demos\">gitlab repository<\/a>.<\/p>\n\n\n\n<p>The code has been registered in 2009 at the Agence pour la Protection des Programmes (APP) [reference IDDN.FR.001.28003.000.S.P.2009.000.31235].<\/p>\n\n\n\n<p>SMC are a class of algorithms for approximate inference in dynamic models. They are used to estimate the state of a system over time, given a sequence of noisy measurements. The basic idea of SMC methods is to represent the probability distribution over the state of the system at each time step using a set of weighted samples, or particles. The particles are propagated through time using the dynamics of the system, and their weights are updated based on the measurements that are made. The weights are used to approximate the true distribution, and they can be used to compute various statistics of interest, such as the mean and covariance of the state. SMC methods are useful because they allow us to make inferences about the state of a system in a computationally efficient manner, without requiring us to compute the true distribution explicitly. They have been applied in many fields.<\/p>\n\n\n\n<h1 class=\"wp-block-heading\">Context<\/h1>\n\n\n\n<p>In all these examples a mobile is describe by its state <img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/localhost:8888\/wordpress\/wp-content\/ql-cache\/quicklatex.com-351de704673400097afff0f28b07b11d_l3.png\" class=\"ql-img-inline-formula \" alt=\"&#88;&#95;&#116;\" title=\"Rendered by QuickLaTeX.com\" height=\"15\" width=\"20\" style=\"vertical-align: -3px;\"\/> in <img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/localhost:8888\/wordpress\/wp-content\/ql-cache\/quicklatex.com-770f2472b162241a2dd7eec50baa137a_l3.png\" class=\"ql-img-inline-formula \" alt=\"&#92;&#109;&#97;&#116;&#104;&#98;&#98;&#123;&#82;&#125;&#94;&#110;\" title=\"Rendered by QuickLaTeX.com\" height=\"12\" width=\"21\" style=\"vertical-align: 0px;\"\/> at time <img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/localhost:8888\/wordpress\/wp-content\/ql-cache\/quicklatex.com-6b97bb0f65c75b6cc0fba1868749478d_l3.png\" class=\"ql-img-inline-formula \" alt=\"&#116;\" title=\"Rendered by QuickLaTeX.com\" height=\"12\" width=\"6\" style=\"vertical-align: 0px;\"\/> which is a Markov process described by a Markov transition kernel <img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/localhost:8888\/wordpress\/wp-content\/ql-cache\/quicklatex.com-b9b58c39b8edee516087228c48f0f7a4_l3.png\" class=\"ql-img-inline-formula \" alt=\"&#81;&#95;&#116;&#40;&#120;&#44;&#100;&#120;&#39;&#41;\" title=\"Rendered by QuickLaTeX.com\" height=\"19\" width=\"75\" style=\"vertical-align: -5px;\"\/> with:<\/p>\n<p><p class=\"ql-center-displayed-equation\" style=\"line-height: 19px;\"><span class=\"ql-right-eqno\"> &nbsp; <\/span><span class=\"ql-left-eqno\"> &nbsp; <\/span><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/localhost:8888\/wordpress\/wp-content\/ql-cache\/quicklatex.com-9317f06167a117ddd3fe3361d1584ecb_l3.png\" height=\"19\" width=\"229\" class=\"ql-img-displayed-equation \" alt=\"&#92;&#91; &#81;&#95;&#116;&#40;&#120;&#44;&#65;&#41;&#32;&#61;&#32;&#92;&#109;&#97;&#116;&#104;&#98;&#98;&#123;&#80;&#125;&#40;&#88;&#95;&#116;&#92;&#105;&#110;&#32;&#65;&#124;&#88;&#95;&#48;&#61;&#120;&#41; &#92;&#93;\" title=\"Rendered by QuickLaTeX.com\"\/><\/p><\/p>\n<p>for example it could be given by an SDE<\/p>\n<p><p class=\"ql-center-displayed-equation\" style=\"line-height: 19px;\"><span class=\"ql-right-eqno\"> &nbsp; <\/span><span class=\"ql-left-eqno\"> &nbsp; <\/span><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/localhost:8888\/wordpress\/wp-content\/ql-cache\/quicklatex.com-5a124d36e8fa2a3f3a6a265d52967994_l3.png\" height=\"19\" width=\"212\" class=\"ql-img-displayed-equation \" alt=\"&#92;&#91; &#100;&#32;&#88;&#95;&#116;&#32;&#61;&#32;&#102;&#40;&#88;&#95;&#116;&#41;&#92;&#44;&#100;&#116;&#32;&#43;&#32;&#103;&#40;&#88;&#95;&#116;&#41;&#92;&#44;&#100;&#66;&#95;&#116; &#92;&#93;\" title=\"Rendered by QuickLaTeX.com\"\/><\/p><\/p>\n<p>The observation process is <img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/localhost:8888\/wordpress\/wp-content\/ql-cache\/quicklatex.com-89ebcdbd21e41e70166e16676e71dffd_l3.png\" class=\"ql-img-inline-formula \" alt=\"&#89;&#95;&#107;\" title=\"Rendered by QuickLaTeX.com\" height=\"15\" width=\"17\" style=\"vertical-align: -3px;\"\/> which is the observation of <img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/localhost:8888\/wordpress\/wp-content\/ql-cache\/quicklatex.com-351de704673400097afff0f28b07b11d_l3.png\" class=\"ql-img-inline-formula \" alt=\"&#88;&#95;&#116;\" title=\"Rendered by QuickLaTeX.com\" height=\"15\" width=\"20\" style=\"vertical-align: -3px;\"\/> at time <img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/localhost:8888\/wordpress\/wp-content\/ql-cache\/quicklatex.com-4211092fe5a03aff3e8810a9c7f20d60_l3.png\" class=\"ql-img-inline-formula \" alt=\"&#116;&#95;&#107;\" title=\"Rendered by QuickLaTeX.com\" height=\"15\" width=\"13\" style=\"vertical-align: -3px;\"\/> given by:<\/p>\n<p><p class=\"ql-center-displayed-equation\" style=\"line-height: 19px;\"><span class=\"ql-right-eqno\"> &nbsp; <\/span><span class=\"ql-left-eqno\"> &nbsp; <\/span><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/localhost:8888\/wordpress\/wp-content\/ql-cache\/quicklatex.com-f1736859f16353bc4060f33f6de29d6d_l3.png\" height=\"19\" width=\"249\" class=\"ql-img-displayed-equation \" alt=\"&#92;&#91; &#32;&#32;&#32;&#89;&#95;&#107;&#32;&#61;&#32;&#104;&#40;&#88;&#95;&#123;&#116;&#95;&#107;&#125;&#41;&#43;&#92;&#116;&#101;&#120;&#116;&#114;&#109;&#123;&#111;&#98;&#115;&#101;&#114;&#118;&#97;&#116;&#105;&#111;&#110;&#32;&#110;&#111;&#105;&#115;&#101;&#125; &#92;&#93;\" title=\"Rendered by QuickLaTeX.com\"\/><\/p><\/p>\n<p>with <img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/localhost:8888\/wordpress\/wp-content\/ql-cache\/quicklatex.com-90d55f50f199c3ecf9af22628c1b7785_l3.png\" class=\"ql-img-inline-formula \" alt=\"&#48;&#61;&#116;&#95;&#48;&#60;&#116;&#95;&#49;&#60;&#116;&#95;&#50;&#60;&#92;&#99;&#100;&#111;&#116;&#115;\" title=\"Rendered by QuickLaTeX.com\" height=\"15\" width=\"165\" style=\"vertical-align: -3px;\"\/>.<\/p>\n\n\n\n<h1 class=\"wp-block-heading\">Target tracking with angle only and obstacles<\/h1>\n\n\n\n<div class=\"page\" title=\"Page 40\">\n<div class=\"layoutArea\">\n<div class=\"column\">\n<p><em><span style=\"font-size: revert;\"><img loading=\"lazy\" decoding=\"async\" class=\"size-medium wp-image-2450 alignleft\" src=\"http:\/\/localhost:8888\/wordpress\/wp-content\/uploads\/pictures\/demo3-geometry_eng-281x300.png\" alt=\"\" width=\"281\" height=\"300\" srcset=\"http:\/\/localhost:8888\/wordpress\/wp-content\/uploads\/pictures\/demo3-geometry_eng-281x300.png 281w, http:\/\/localhost:8888\/wordpress\/wp-content\/uploads\/pictures\/demo3-geometry_eng.png 554w\" sizes=\"(max-width: 281px) 100vw, 281px\" \/><\/span><\/em><\/p>\n<p><em><span style=\"font-size: revert;\">The original idea for this very nice example came from<\/span><span style=\"font-size: revert;\"> <a href=\"http:\/\/www.simonmaskell.com\">Simon Maskell<\/a>.<\/span><\/em><\/p>\n<\/div>\n<div class=\"column\">\n<p>The idea is to follow a mobile on the plane. We have <img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/localhost:8888\/wordpress\/wp-content\/ql-cache\/quicklatex.com-6046e01a81406b0611d267d68b760dfb_l3.png\" class=\"ql-img-inline-formula \" alt=\"&#83;\" title=\"Rendered by QuickLaTeX.com\" height=\"12\" width=\"12\" style=\"vertical-align: 0px;\"\/> observation stations (squares) and <img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/localhost:8888\/wordpress\/wp-content\/ql-cache\/quicklatex.com-649dd83fe5ed56224cf7f675eed3a2b9_l3.png\" class=\"ql-img-inline-formula \" alt=\"&#76;\" title=\"Rendered by QuickLaTeX.com\" height=\"12\" width=\"12\" style=\"vertical-align: 0px;\"\/> obstacles (black segments) that hide the movement of the mobile.<\/p>\n<p>The measures are the angle of the line of sight plus noise.&nbsp;<\/p>\n<p>In the video, we place 1)&nbsp;the observations stations, 2) the obstacles, then 3) we plot the trajectory of the mobile, and 4) we simulate the observations, finally 5) we plot the evolution of the particles of the approximation of the nonlinear filter.<\/p>\n<p>We notice the ability of the particle filter to handle non Gaussian situations.&nbsp;<\/p>\n<\/div>\n<\/div>\n<\/div>\n\n\n\n<div style=\"height:22px\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n\n<div class=\"wp-block-columns is-layout-flex wp-container-core-columns-is-layout-1 wp-block-columns-is-layout-flex\">\n<div class=\"wp-block-column is-layout-flow wp-block-column-is-layout-flow\" style=\"flex-basis:25%\"><\/div>\n\n\n\n<div class=\"wp-block-column is-layout-flow wp-block-column-is-layout-flow\" style=\"flex-basis:50%\">\n<figure class=\"wp-block-video\"><video controls src=\"http:\/\/localhost:8888\/wordpress\/wp-content\/uploads\/movies\/smc-demo-4-tracking-with-obstacles.mov\"><\/video><\/figure>\n<\/div>\n\n\n\n<div class=\"wp-block-column is-layout-flow wp-block-column-is-layout-flow\" style=\"flex-basis:25%\"><\/div>\n<\/div>\n\n\n\n<h1 class=\"wp-block-heading\">A mobile in Manhattan<\/h1>\n\n\n\n<div class=\"wp-block-columns is-layout-flex wp-container-core-columns-is-layout-3 wp-block-columns-is-layout-flex\">\n<div class=\"wp-block-column is-layout-flow wp-block-column-is-layout-flow\" style=\"flex-basis:100%\">\n\n\n<div class=\"wp-block-columns are-vertically-aligned-bottom is-layout-flex wp-container-core-columns-is-layout-2 wp-block-columns-is-layout-flex\">\n<div class=\"wp-block-column is-vertically-aligned-bottom is-layout-flow wp-block-column-is-layout-flow\">\n<figure class=\"wp-block-image size-full is-resized\"><a href=\"http:\/\/localhost:8888\/wordpress\/wp-content\/uploads\/pictures\/smc-demo5-zone.png\"><img loading=\"lazy\" decoding=\"async\" width=\"996\" height=\"950\" src=\"http:\/\/localhost:8888\/wordpress\/wp-content\/uploads\/pictures\/smc-demo5-zone.png\" alt=\"\" class=\"wp-image-2466\" style=\"width:312px;height:297px\" srcset=\"http:\/\/localhost:8888\/wordpress\/wp-content\/uploads\/pictures\/smc-demo5-zone.png 996w, http:\/\/localhost:8888\/wordpress\/wp-content\/uploads\/pictures\/smc-demo5-zone-300x286.png 300w, http:\/\/localhost:8888\/wordpress\/wp-content\/uploads\/pictures\/smc-demo5-zone-768x733.png 768w, http:\/\/localhost:8888\/wordpress\/wp-content\/uploads\/pictures\/smc-demo5-zone-624x595.png 624w\" sizes=\"(max-width: 996px) 100vw, 996px\" \/><\/a><figcaption class=\"wp-element-caption\">G\u00e9olocalisation of a mobile<\/figcaption><\/figure>\n<\/div>\n\n\n\n<div class=\"wp-block-column is-vertically-aligned-bottom is-layout-flow wp-block-column-is-layout-flow\">\n<figure class=\"wp-block-image size-large is-resized\"><a href=\"http:\/\/localhost:8888\/wordpress\/wp-content\/uploads\/pictures\/smc-demo5-map.png\"><img loading=\"lazy\" decoding=\"async\" width=\"1024\" height=\"863\" src=\"http:\/\/localhost:8888\/wordpress\/wp-content\/uploads\/pictures\/smc-demo5-map-1024x863.png\" alt=\"\" class=\"wp-image-2467\" style=\"width:284px;height:239px\" srcset=\"http:\/\/localhost:8888\/wordpress\/wp-content\/uploads\/pictures\/smc-demo5-map-1024x863.png 1024w, http:\/\/localhost:8888\/wordpress\/wp-content\/uploads\/pictures\/smc-demo5-map-300x253.png 300w, http:\/\/localhost:8888\/wordpress\/wp-content\/uploads\/pictures\/smc-demo5-map-768x647.png 768w, http:\/\/localhost:8888\/wordpress\/wp-content\/uploads\/pictures\/smc-demo5-map-1536x1294.png 1536w, http:\/\/localhost:8888\/wordpress\/wp-content\/uploads\/pictures\/smc-demo5-map-624x526.png 624w, http:\/\/localhost:8888\/wordpress\/wp-content\/uploads\/pictures\/smc-demo5-map.png 1550w\" sizes=\"(max-width: 1024px) 100vw, 1024px\" \/><\/a><figcaption class=\"wp-element-caption\">Attenuation map.<\/figcaption><\/figure>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n\n\n\n<p>It is about geolocating a mobile, a person using a cell phone, in a simplified Manhattan. The mobile can move in the street (white zone) or in the buildings (grey zones).<\/p>\n<p>The area is covered by several antennas represented by green dots. To each of these antennas corresponds an attenuation map (bottom figure). This map describes in each point of the area the attenuation of the signal transmitted by the mobile and received by the antenna. In this simplified version it is assumed that this signal is zero <br \/>when the mobile is in a building, and that it decreases when the mobile moves away from the antenna.<\/p>\n<p>In the following video we position the antennas (green dots) and then we define the trajectory of the mobile: it can move in the street or in buildings.<\/p>\n\n\n\n<div style=\"height:45px\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n\n<div class=\"wp-block-columns is-layout-flex wp-container-core-columns-is-layout-4 wp-block-columns-is-layout-flex\">\n<div class=\"wp-block-column is-layout-flow wp-block-column-is-layout-flow\" style=\"flex-basis:25%\"><\/div>\n\n\n\n<div class=\"wp-block-column is-layout-flow wp-block-column-is-layout-flow\" style=\"flex-basis:50%\">\n<figure class=\"wp-block-video\"><video controls src=\"http:\/\/localhost:8888\/wordpress\/wp-content\/uploads\/movies\/smc-demo-7-geolocation.mov\"><\/video><\/figure>\n<\/div>\n\n\n\n<div class=\"wp-block-column is-layout-flow wp-block-column-is-layout-flow\" style=\"flex-basis:25%\"><\/div>\n<\/div>\n","protected":false},"excerpt":{"rendered":"<p>A set of Matlab software demonstrations of particle filtering for tracking applications, developed in the context of lectures given from 2005 to 2013: and in many seminars and in doctoral training programs. I also used these Matlab software demonstrations in a course I taught in Toulon from 2004 to 2008 as part of the Master [&hellip;]<\/p>\n","protected":false},"author":1,"featured_media":0,"parent":11,"menu_order":0,"comment_status":"closed","ping_status":"closed","template":"page-templates\/full-width.php","meta":{"footnotes":""},"class_list":["post-201","page","type-page","status-publish","hentry"],"_links":{"self":[{"href":"http:\/\/localhost:8888\/wordpress\/wp-json\/wp\/v2\/pages\/201"}],"collection":[{"href":"http:\/\/localhost:8888\/wordpress\/wp-json\/wp\/v2\/pages"}],"about":[{"href":"http:\/\/localhost:8888\/wordpress\/wp-json\/wp\/v2\/types\/page"}],"author":[{"embeddable":true,"href":"http:\/\/localhost:8888\/wordpress\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"http:\/\/localhost:8888\/wordpress\/wp-json\/wp\/v2\/comments?post=201"}],"version-history":[{"count":49,"href":"http:\/\/localhost:8888\/wordpress\/wp-json\/wp\/v2\/pages\/201\/revisions"}],"predecessor-version":[{"id":3213,"href":"http:\/\/localhost:8888\/wordpress\/wp-json\/wp\/v2\/pages\/201\/revisions\/3213"}],"up":[{"embeddable":true,"href":"http:\/\/localhost:8888\/wordpress\/wp-json\/wp\/v2\/pages\/11"}],"wp:attachment":[{"href":"http:\/\/localhost:8888\/wordpress\/wp-json\/wp\/v2\/media?parent=201"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}