|
4 | 4 | ============================================================= |
5 | 5 | 11 - Tracking linear Levy transition models with the MPF |
6 | 6 | ============================================================= |
7 | | -In line with the tutorial examples of the Kalman and particle |
8 | | -filters in Stone Soup, a simplified single-target tracking |
9 | | -example without clutter is provided here to demonstrate the |
10 | | -use of linear transition models driven by non-Gaussian Levy |
11 | | -noise, as well as how to perform inference tasks on this class |
| 7 | +In line with the tutorial examples of the Kalman and particle |
| 8 | +filters in Stone Soup, a simplified single-target tracking |
| 9 | +example without clutter is provided here to demonstrate the |
| 10 | +use of linear transition models driven by non-Gaussian Levy |
| 11 | +noise, as well as how to perform inference tasks on this class |
12 | 12 | of models using the Marginalized Particle Filter (MPF). |
13 | 13 | """ |
14 | 14 |
|
|
22 | 22 | from datetime import datetime, timedelta |
23 | 23 | from stonesoup.types.groundtruth import GroundTruthPath, GroundTruthState |
24 | 24 | from stonesoup.models.base_driver import NoiseCase |
25 | | -from stonesoup.models.driver import AlphaStableNSMDriver |
| 25 | +from stonesoup.models.driver import AlphaStableNSMDriver |
26 | 26 | from stonesoup.models.transition.levy_linear import LevyLangevin, CombinedLinearLevyTransitionModel |
27 | 27 |
|
28 | 28 | start_time = datetime.now().replace(microsecond=0) |
|
60 | 60 | # %% |
61 | 61 |
|
62 | 62 | # Random seem for reproducibility |
63 | | -seed = 1 |
| 63 | +seed = 1 |
64 | 64 |
|
65 | 65 | # Driving process parameters |
66 | 66 | mu_W = 0 |
|
103 | 103 | # Note that the animated plotter requires a list of timesteps as an input, and that ``tail_length`` |
104 | 104 | # is set to 0.3. This means that each data point will be on display for 30% of the total |
105 | 105 | # simulation time. The mapping argument is ``[0, 2]`` because those are the :math:`x` and |
106 | | -# :math:`y` position indices from our state vector. |
107 | | -# If a static plotter is preferred, the :class:`Plotterly` class can be used instead. |
| 106 | +# :math:`y` position indices from our state vector/Users/chongzhenyuen/Documents/iib-project/stonesoup/stonesoup/models. |
| 107 | +# If a static plotter is preferred, the :class:`Plotterly`` class can be used instead. |
108 | 108 |
|
109 | 109 | # from stonesoup.plotter import AnimatedPlotterly |
110 | 110 | # plotter = AnimatedPlotterly(timesteps, tail_length=1.0, width=600, height=600) |
|
238 | 238 | # .. [#] Lemke, Tatjana, and Simon J. Godsill, 'Inference for models with asymmetric α -stable noise processes', in Siem Jan Koopman, and Neil Shephard (eds), Unobserved Components and Time Series Econometrics (Oxford, 2015; online edn, Oxford Academic, 21 Jan. 2016) |
239 | 239 | # |
240 | 240 | # .. [#] S. Godsill, M. Riabiz, and I. Kontoyiannis, “The L ́evy state space model,” in 2019 53rd Asilomar Conference on Signals, Systems, and Computers, 2019, pp. 487–494. |
241 | | - |
242 | | - |
243 | | - |
244 | | - |
245 | | - |
0 commit comments