Data Analysis and Machine Learning Lecture 10: Applications of Kalman Filtering
|
|
- Anabel Martin
- 6 years ago
- Views:
Transcription
1 Data Analysis and Machine Learning Lecture 10: Applications of Kalman Filtering Lecturer, Imperial College London Founder, CEO, Thalesians Ltd
2 The Newtonian system I Let r be the position of a particle, v its velocity, and a its (constant) acceleration. If t is the time, then dr/dt dv/dt da/dt = The state of our system x, is the vector (r, v, a) T. r v a Denote the matrix of ones and zeros above by A. Then the evolution of the state is described by the matrix differential equation dx t dt = Ax t. By analogy with the scalar ordinary differential equation (ODE), its solution is x t = e At x 0, where x t = (r t, v t, a t ) T is the state of the system at time t and the matrix exponential is defined by the (matrix) Taylor series expansion e At := j=0 (At) j. j!.
3 The Newtonian system II Note that, for s t T, x t = e A(t s) x s. Let us now discretise this equation using h k := t s as the time interval between the time ticks k 1 and k (for k N ): where x k = F k x k 1, F k := e Ah (Ah k = k ) j j! j=0 = hk h 2 2 k = 1 h k hk 2/2 0 1 h k, since all the powers of A greater than 2 are zero matrices. We observe that this process model fits the Kalman filter paradigm.
4 The variance-scaled Wiener process with drift Consider the stochastic differential equation (SDE) dx t = µ dt + σ dw t. Integrating both sides, we get hence Note that, for s t T, t 0 t t dx s = µ ds+ σ dw s, 0 0 X t = X 0 + µt + σw t. X t = X s + µ(t s)+σ(w t W s ). Let us now discretise this equation using h k := t s as the time interval between the time ticks k 1 and k (for k N ): x k = F k x k 1 + a k + w k, where F k = I, a k = µh k, w k N(0, σ 2 h k ). We observe that this process model also fits the Kalman filter paradigm.
5 Tutorial Consider the one-dimensional Ornstein-Uhlenbeck (OU) process dx t = θ(µ X t ) dt + σ dw t, where X t R, X 0 = x 0, and θ > 0, µ and σ > 0 are constants. Discretise this SDE and formulate the Kalman filter process model. Write a Python function to generate realisations of this process. Implement a Kalman filter in Python. Apply it to your process model and test on generated process realisations (experiment with adding different levels of white noise to the observations).
6 Tutorial: Hints I First, we must solve the OU process SDE. The solution is well-known, so you may simply quote it. However, this is how you would solve it. Recall Itō s lemma: if dx t = µ t dt + σ t dw t, f(x, t) is a function, and Y t := f(x t, t), then ( f dy t = t + f µ t x + (σ ) 2 2 ) f 2 x 2 dt + σ f t x dw t. We apply the lemma for the OU SDE (so that, in our case, µ t = θ(µ x t ) and σ t = σ) and f(x, t) := xe θt. Then f t = θxe θt, f x = eθt, 2 f = x 2 0, and ( dy t = θx t e θt + θ(µ X t )e θt) dt + σe θt dw t On integrating both sides, we get = θµe θt dt + σe θt dw t. t 0 t dy u = θµe θu ds+ 0 X t e θt X 0 = µ(e θt 1)+ t σe θu dw u, 0 t σe θu dw u. 0
7 Tutorial: Hints II Thus the solution of the OU SDE is It is easy to see that, for s t T, t X t = X 0 e θt + µ(1 e θt )+ σe θ(t u) dw u. 0 t X t = X s e θ(t s) + µ(1 e θ(t s) )+ σe θ(t u) dw u. s It is well known that an Itō integral, t s f(u) dw u, of a deterministic integrand, f(u), is a Gaussian random variable with mean 0 and variance t 0 f 2 (u) du. In our case, f(u) = σe θ(t u), and t 0 f 2 (u) du = σ2 2θ ( 1 e 2θ(t s)).
8 Tutorial: Hints III Let us now discretise the equation for X t using h k := t s as the time interval between the time ticks k 1 and k (for k N ): x k = F k x k 1 + a k + w k, ( ( ) ) where F k = e θh k, a k = µ(1 e θh k), w k N 0, σ2 2θ 1 e 2θh k. We observe that this process model fits the Kalman filter paradigm, and we can apply the Kalman filter equations directly.
9 The loglikelihood function for the Kalman filter We can view the Kalman filter as parameterised by some parameter vector θ R d θ, so that F k, a k, Q k, H k, b k, R k are functions of θ: F k (θ), a k (θ), Q k (θ), H k (θ), b k (θ), R k (θ). Then, for all T N, we can show that the following holds: ln(l(θ)) := ln(f(y 1,...,y T θ)) = T i=1 ln(f(y t θ; y 1,...,y t 1 )) = T 2 ln(2π) 1 T ln(det(s t )) 1 T 2 2 i=1 i=1 ỹ T t S 1 t ỹ. We can use this loglikelihood function as our objective function in the estimation of the parameter vector θ. See, for example, tutorials/xfghtmlnode92.html
10 Outlier detection Note that the predicted observation is distributed as ( ) N H k ˆx k k 1, H k P k k 1 Hk T. Just as we can assign a z-score to y k if y is one-dimensional, we can assign a Mahalanobis distance to it (which corresponds to the z-score if y k is one-dimensional). In general, the Mahalanobis norm of a vector y with respect to N (µ, Σ) is given by y N(µ,Σ) = (y µ) T Σ 1 (y µ). It measures the distance of y R m from the centroid (multidimensional mean) of the distribution. y 2 N(µ,Σ) follows the χ2 -distribution with m degrees of freedom. Thus we can set a cut-off for y k, e.g. on the basis of the 0.975th quantile of the χ 2 -distribution.
11 The extended Kalman filter (EKF) Suppose that our n-dimensional state vector evolves according to the process model x k = f k (x k 1 )+w k, and suppose that our m-dimensional observation vector is related to the state by the observation model y k = h k (x k )+v k where f and h are differentiable functions, w k N(0, Q k ) are R n -valued random variables, and v k N(0, R k ) are R m -valued random variables. {x 0, w 1,...,w k, v 1,...,v k } are mutually independent. The extended Kalman filter [McE66, SSM62], [Hay01, Section 1.6] algorithm is as follows: The prediction step: Predicted (prior) signal state estimate: ˆx k k 1 = f k (ˆx k 1 k 1 ). Predicted (prior) error covariance: P k k 1 = F k P k 1 k 1 F k T + Qk. The update (or correction) step: Innovation, a.k.a. observation residual: ỹ k = y k h k (ˆx k k 1 ). Innovation (observation residual) covariance:s k = H k P k k 1 H k T + Rk. (Optimal) Kalman gain: K k = P k k 1 H T k S 1. k Updated (posterior) signal state estimate: ˆx k k = ˆx k k 1 + K k ỹ k. Updated (posterior) error covariance: P k k = (I K k H k )P k k 1. Here F k = f x and H k = ˆxk 1 k 1 h x. So the idea is simple: linearise the model ˆxk k 1 using Jacobians.
12 Credit risk Credit risk is the risk that an obligor does not honour his payment obligations [Sch03, page 1]. This risk is directly or indirectly associated with some events, called credit events. They can be hard (e.g. bankruptcy) or soft (e.g. restructuring). Credit risk can be categorised into credit deterioration and default risk. A bond that has credit risk associated with it is risky. A bond that has no credit risk associated with it is riskless. The credit spread is the difference between the (promised) yield on the risky bond and the yield on riskless bonds [KN15, pages ].
13 Credit spread How can we quantify the credit spread of a risky bond? One way is to solve for z t0 in the following equation: P t0 = n i=1 ( 1+ cf ti (s i t0 + z t0 ) δ where P t0 is the dirty market price of the bond at time, cf ti is the cash flow generated by the bond at time t i, s i is the zero-coupon swap rate of appropriate maturity for this cashflow, δ is the frequency of the cashflows expressed as a fraction of the year. The resulting z t0 is called the zero volatility spread, or Z-spread for short. For a given issuer, we are interested in modelling the term structure of the Z-spreads across a universe of this issuer s bonds. The Z-spread is then viewed as a function of τ : z t0 (τ); τ could be the time of maturity or, for example, the modified duration. ) i,
14 Modelling the Z-spread The function z t0 (τ) can be written as z t0 (τ) = g(τ; θ t0 ), where θ t0 is a d-dimensional vector of parameters. Its dependence on indicates that the Z-spread curve (as a function of τ) evolves as time progresses. One of our tasks, then, is to keep estimating θ t0 as moves on. Suppose that for a particular issuer we have a universe of K bonds with Z-spreads z (1),...,z (K) and maturities (or modified durations, etc.) τ (1),...,τ (K). Each of these Z-spreads may not lie exactly on the Z-spread curve z t0 (τ) due to the idiosyncracies of that particular bond, so we allow there to be an idiosyncratic spread, λ (k), k {1,...,K}: z (k) = z t0 (τ (k) )+λ (k). As indicated by their dependence on, the idiosyncratic spreads also evolve over time. We will have to keep computing their updated estimates as well.
15 Setting up the EKF: the state and observation Our state at is the vector (θ t0 ) 1. (θ t0 ) d x t0 = λ (1)... λ (K) We are observing individual bond prices, thus our observation at time is y t0 = P (k), the dirty market price of bond k for some k {1,...,K}. NB! Chances are that we our market feeds are providing us with clean prices. We should therefore remember to convert these prices to dirty prices before we proceed.
16 Setting up the EKF: the observation model The function h that maps our state to the corresponding observation is given by h (k) (x t0 ) = P (k) (z (k) ) n (k) = i=1 n (k) = i=1 cf (k) t i ( ( ) ) 1+ st i 0 + z (k) i δ (k) cf (k) t i ( ( ) ) 1+ st i 0 + z t0 (τ (k) )+λ (k) i, δ (k) where z t0 (τ (k) ) is, of couse, our curve model function, whose parameters θ t0 appear in our state, evaluated at the maturity of the bond, and λ (k) is the idiosyncratic spread for this particular bond, which is likewise in our state. Thus, given our state vector, we can evaluate h (k) (x t0 ). Note that the equation above also depends on the appropriate zero-coupon swap rates s i. These are fast-moving and can be provided exogenously. The reason why we need an EKF rather than a KF is that h (k) (x t0 ) is a nonlinear function of the state.
17 Setting up the EKF: the observation Jacobian Let the scalar parameter α t0 be a particular element of our parameter vector θ t0, so it is (θ t0 ) j for some j {1,...,d}. Then h (k) = α t0 xt0 h(k) z = z t0 xt0 α t0 xt0 δ n(k) (k) i=1 i cf (k) t i ( ( ) 1+ st i 0 + z t0 (τ (k) )+λ (k) where z α can be computed analytically for many simple curve models. t0 ) i z τ δ (k) α (k t0 Similarly, for j {1,...,K}, h (k) h (k) λ (j) = z t0 xt0 z λ (j) = δ (k) n(k) i cf (k) t i i=1 ( ( t xt0 xt st i +z t0 (τ (k) )+λ (k) ) i, j = k; 0 t δ (k)) 0 0, otherwise.
18 market depth orders liquidity RFQ, RFS, CAT responses other sources of information venues (economic indicators, news feeds, etc.) (D2C, D2D) trading strategy pricing quoting across the board: beta price alpha price spreading skewing (passive hedging) alerts, circuit breakers, recording and audit, visualisation, control hit-lift protection hedging position service optimal execution strategies risk service
19 Simon Haykin, editor. Kalman Filtering and Neural Networks. John Wiley and Sons, Inc., Robert L. Kosowski and Salih N. Neftci. Principles of Financial Engineering. Academic Press, Bruce A. McElhoe. An assessment of the navigation and course corrections for a manned flyby of mars or venus. IEEE Transactions on Aerospace and Electronic Systems, AES-2(4): , Philipp J. Schönbucher. Credit Derivatives Pricing Models: Models, Pricing and Implementation. Wiley Finance, Gerald L. Smith, Stanley F. Schmidt, and Leonard A. McGee. Application of statistical filter theory to the optimal estimation of position and velocity on board a circumlunar vehicle. NASA technical report R-135, National Aeronautics and Space Administration (NASA), 1962.
Brownian Motion. An Undergraduate Introduction to Financial Mathematics. J. Robert Buchanan. J. Robert Buchanan Brownian Motion
Brownian Motion An Undergraduate Introduction to Financial Mathematics J. Robert Buchanan 2010 Background We have already seen that the limiting behavior of a discrete random walk yields a derivation of
More informationGaussian Process Approximations of Stochastic Differential Equations
Gaussian Process Approximations of Stochastic Differential Equations Cédric Archambeau Dan Cawford Manfred Opper John Shawe-Taylor May, 2006 1 Introduction Some of the most complex models routinely run
More information1 Kalman Filter Introduction
1 Kalman Filter Introduction You should first read Chapter 1 of Stochastic models, estimation, and control: Volume 1 by Peter S. Maybec (available here). 1.1 Explanation of Equations (1-3) and (1-4) Equation
More informationStochastic contraction BACS Workshop Chamonix, January 14, 2008
Stochastic contraction BACS Workshop Chamonix, January 14, 2008 Q.-C. Pham N. Tabareau J.-J. Slotine Q.-C. Pham, N. Tabareau, J.-J. Slotine () Stochastic contraction 1 / 19 Why stochastic contraction?
More informationPrediction of ESTSP Competition Time Series by Unscented Kalman Filter and RTS Smoother
Prediction of ESTSP Competition Time Series by Unscented Kalman Filter and RTS Smoother Simo Särkkä, Aki Vehtari and Jouko Lampinen Helsinki University of Technology Department of Electrical and Communications
More informationStochastic Calculus. Kevin Sinclair. August 2, 2016
Stochastic Calculus Kevin Sinclair August, 16 1 Background Suppose we have a Brownian motion W. This is a process, and the value of W at a particular time T (which we write W T ) is a normally distributed
More informationBayes Filter Reminder. Kalman Filter Localization. Properties of Gaussians. Gaussians. Prediction. Correction. σ 2. Univariate. 1 2πσ e.
Kalman Filter Localization Bayes Filter Reminder Prediction Correction Gaussians p(x) ~ N(µ,σ 2 ) : Properties of Gaussians Univariate p(x) = 1 1 2πσ e 2 (x µ) 2 σ 2 µ Univariate -σ σ Multivariate µ Multivariate
More informationAutonomous Navigation for Flying Robots
Computer Vision Group Prof. Daniel Cremers Autonomous Navigation for Flying Robots Lecture 6.2: Kalman Filter Jürgen Sturm Technische Universität München Motivation Bayes filter is a useful tool for state
More informationLecture 1: Pragmatic Introduction to Stochastic Differential Equations
Lecture 1: Pragmatic Introduction to Stochastic Differential Equations Simo Särkkä Aalto University, Finland (visiting at Oxford University, UK) November 13, 2013 Simo Särkkä (Aalto) Lecture 1: Pragmatic
More informationLecture 2: From Linear Regression to Kalman Filter and Beyond
Lecture 2: From Linear Regression to Kalman Filter and Beyond January 18, 2017 Contents 1 Batch and Recursive Estimation 2 Towards Bayesian Filtering 3 Kalman Filter and Bayesian Filtering and Smoothing
More informationFactor Analysis and Kalman Filtering (11/2/04)
CS281A/Stat241A: Statistical Learning Theory Factor Analysis and Kalman Filtering (11/2/04) Lecturer: Michael I. Jordan Scribes: Byung-Gon Chun and Sunghoon Kim 1 Factor Analysis Factor analysis is used
More informationLecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Approximations
Lecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Approximations Simo Särkkä Aalto University, Finland November 18, 2014 Simo Särkkä (Aalto) Lecture 4: Numerical Solution of SDEs November
More informationI forgot to mention last time: in the Ito formula for two standard processes, putting
I forgot to mention last time: in the Ito formula for two standard processes, putting dx t = a t dt + b t db t dy t = α t dt + β t db t, and taking f(x, y = xy, one has f x = y, f y = x, and f xx = f yy
More informationInformation and Credit Risk
Information and Credit Risk M. L. Bedini Université de Bretagne Occidentale, Brest - Friedrich Schiller Universität, Jena Jena, March 2011 M. L. Bedini (Université de Bretagne Occidentale, Brest Information
More informationIntroduction to numerical simulations for Stochastic ODEs
Introduction to numerical simulations for Stochastic ODEs Xingye Kan Illinois Institute of Technology Department of Applied Mathematics Chicago, IL 60616 August 9, 2010 Outline 1 Preliminaries 2 Numerical
More informationNonlinear Parameter Estimation for State-Space ARCH Models with Missing Observations
Nonlinear Parameter Estimation for State-Space ARCH Models with Missing Observations SEBASTIÁN OSSANDÓN Pontificia Universidad Católica de Valparaíso Instituto de Matemáticas Blanco Viel 596, Cerro Barón,
More informationLecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Process Approximations
Lecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Process Approximations Simo Särkkä Aalto University Tampere University of Technology Lappeenranta University of Technology Finland November
More informationFE610 Stochastic Calculus for Financial Engineers. Stevens Institute of Technology
FE610 Stochastic Calculus for Financial Engineers Lecture 3. Calculaus in Deterministic and Stochastic Environments Steve Yang Stevens Institute of Technology 01/31/2012 Outline 1 Modeling Random Behavior
More informationInterest Rate Models:
1/17 Interest Rate Models: from Parametric Statistics to Infinite Dimensional Stochastic Analysis René Carmona Bendheim Center for Finance ORFE & PACM, Princeton University email: rcarmna@princeton.edu
More informationA new unscented Kalman filter with higher order moment-matching
A new unscented Kalman filter with higher order moment-matching KSENIA PONOMAREVA, PARESH DATE AND ZIDONG WANG Department of Mathematical Sciences, Brunel University, Uxbridge, UB8 3PH, UK. Abstract This
More informationA variational radial basis function approximation for diffusion processes
A variational radial basis function approximation for diffusion processes Michail D. Vrettas, Dan Cornford and Yuan Shen Aston University - Neural Computing Research Group Aston Triangle, Birmingham B4
More informationGaussian Process Approximations of Stochastic Differential Equations
Gaussian Process Approximations of Stochastic Differential Equations Cédric Archambeau Centre for Computational Statistics and Machine Learning University College London c.archambeau@cs.ucl.ac.uk CSML
More informationNested Uncertain Differential Equations and Its Application to Multi-factor Term Structure Model
Nested Uncertain Differential Equations and Its Application to Multi-factor Term Structure Model Xiaowei Chen International Business School, Nankai University, Tianjin 371, China School of Finance, Nankai
More informationLinear Filtering of general Gaussian processes
Linear Filtering of general Gaussian processes Vít Kubelka Advisor: prof. RNDr. Bohdan Maslowski, DrSc. Robust 2018 Department of Probability and Statistics Faculty of Mathematics and Physics Charles University
More informationLecture 2: From Linear Regression to Kalman Filter and Beyond
Lecture 2: From Linear Regression to Kalman Filter and Beyond Department of Biomedical Engineering and Computational Science Aalto University January 26, 2012 Contents 1 Batch and Recursive Estimation
More informationStochastic differential equation models in biology Susanne Ditlevsen
Stochastic differential equation models in biology Susanne Ditlevsen Introduction This chapter is concerned with continuous time processes, which are often modeled as a system of ordinary differential
More informationEndogenous Information Choice
Endogenous Information Choice Lecture 7 February 11, 2015 An optimizing trader will process those prices of most importance to his decision problem most frequently and carefully, those of less importance
More informationA new approach for investment performance measurement. 3rd WCMF, Santa Barbara November 2009
A new approach for investment performance measurement 3rd WCMF, Santa Barbara November 2009 Thaleia Zariphopoulou University of Oxford, Oxford-Man Institute and The University of Texas at Austin 1 Performance
More informationStochastic Processes and Advanced Mathematical Finance
Steven R. Dunbar Department of Mathematics 23 Avery Hall University of Nebraska-Lincoln Lincoln, NE 68588-13 http://www.math.unl.edu Voice: 42-472-3731 Fax: 42-472-8466 Stochastic Processes and Advanced
More informationEM-algorithm for Training of State-space Models with Application to Time Series Prediction
EM-algorithm for Training of State-space Models with Application to Time Series Prediction Elia Liitiäinen, Nima Reyhani and Amaury Lendasse Helsinki University of Technology - Neural Networks Research
More informationUsing the Kalman Filter to Estimate the State of a Maneuvering Aircraft
1 Using the Kalman Filter to Estimate the State of a Maneuvering Aircraft K. Meier and A. Desai Abstract Using sensors that only measure the bearing angle and range of an aircraft, a Kalman filter is implemented
More informationNumerical Integration of SDEs: A Short Tutorial
Numerical Integration of SDEs: A Short Tutorial Thomas Schaffter January 19, 010 1 Introduction 1.1 Itô and Stratonovich SDEs 1-dimensional stochastic differentiable equation (SDE) is given by [6, 7] dx
More informationSimulation and Parametric Estimation of SDEs
Simulation and Parametric Estimation of SDEs Jhonny Gonzalez School of Mathematics The University of Manchester Magical books project August 23, 2012 Motivation The problem Simulation of SDEs SDEs driven
More informationNumerical Solutions of ODEs by Gaussian (Kalman) Filtering
Numerical Solutions of ODEs by Gaussian (Kalman) Filtering Hans Kersting joint work with Michael Schober, Philipp Hennig, Tim Sullivan and Han C. Lie SIAM CSE, Atlanta March 1, 2017 Emmy Noether Group
More informationStochastic Volatility and Correction to the Heat Equation
Stochastic Volatility and Correction to the Heat Equation Jean-Pierre Fouque, George Papanicolaou and Ronnie Sircar Abstract. From a probabilist s point of view the Twentieth Century has been a century
More informationAffine Processes. Econometric specifications. Eduardo Rossi. University of Pavia. March 17, 2009
Affine Processes Econometric specifications Eduardo Rossi University of Pavia March 17, 2009 Eduardo Rossi (University of Pavia) Affine Processes March 17, 2009 1 / 40 Outline 1 Affine Processes 2 Affine
More informationDiscretization of SDEs: Euler Methods and Beyond
Discretization of SDEs: Euler Methods and Beyond 09-26-2006 / PRisMa 2006 Workshop Outline Introduction 1 Introduction Motivation Stochastic Differential Equations 2 The Time Discretization of SDEs Monte-Carlo
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 15. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationBayesian inference for stochastic differential mixed effects models - initial steps
Bayesian inference for stochastic differential ixed effects odels - initial steps Gavin Whitaker 2nd May 2012 Supervisors: RJB and AG Outline Mixed Effects Stochastic Differential Equations (SDEs) Bayesian
More informationGeneralized Autoregressive Score Models
Generalized Autoregressive Score Models by: Drew Creal, Siem Jan Koopman, André Lucas To capture the dynamic behavior of univariate and multivariate time series processes, we can allow parameters to be
More informationDUBLIN CITY UNIVERSITY
DUBLIN CITY UNIVERSITY SAMPLE EXAMINATIONS 2017/2018 MODULE: QUALIFICATIONS: Simulation for Finance MS455 B.Sc. Actuarial Mathematics ACM B.Sc. Financial Mathematics FIM YEAR OF STUDY: 4 EXAMINERS: Mr
More informationFiltering bond and credit default swap markets
Filtering bond and credit default swap markets Peter Cotton May 20, 2017 Overview Disclaimer Filtering credit markets A visual introduction A unified state space for bond and CDS markets Kalman filtering
More informationFinancial Econometrics
Financial Econometrics Nonlinear time series analysis Gerald P. Dwyer Trinity College, Dublin January 2016 Outline 1 Nonlinearity Does nonlinearity matter? Nonlinear models Tests for nonlinearity Forecasting
More informationStochastic Integration and Stochastic Differential Equations: a gentle introduction
Stochastic Integration and Stochastic Differential Equations: a gentle introduction Oleg Makhnin New Mexico Tech Dept. of Mathematics October 26, 27 Intro: why Stochastic? Brownian Motion/ Wiener process
More informationIntroduction to Diffusion Processes.
Introduction to Diffusion Processes. Arka P. Ghosh Department of Statistics Iowa State University Ames, IA 511-121 apghosh@iastate.edu (515) 294-7851. February 1, 21 Abstract In this section we describe
More informationIntroduction. Stochastic Processes. Will Penny. Stochastic Differential Equations. Stochastic Chain Rule. Expectations.
19th May 2011 Chain Introduction We will Show the relation between stochastic differential equations, Gaussian processes and methods This gives us a formal way of deriving equations for the activity of
More informationProbabilistic Fundamentals in Robotics. DAUIN Politecnico di Torino July 2010
Probabilistic Fundamentals in Robotics Gaussian Filters Basilio Bona DAUIN Politecnico di Torino July 2010 Course Outline Basic mathematical framework Probabilistic models of mobile robots Mobile robot
More informationB8.3 Mathematical Models for Financial Derivatives. Hilary Term Solution Sheet 2
B8.3 Mathematical Models for Financial Derivatives Hilary Term 18 Solution Sheet In the following W t ) t denotes a standard Brownian motion and t > denotes time. A partition π of the interval, t is a
More informationLeast Squares Estimators for Stochastic Differential Equations Driven by Small Lévy Noises
Least Squares Estimators for Stochastic Differential Equations Driven by Small Lévy Noises Hongwei Long* Department of Mathematical Sciences, Florida Atlantic University, Boca Raton Florida 33431-991,
More informationThe moment-generating function of the log-normal distribution using the star probability measure
Noname manuscript No. (will be inserted by the editor) The moment-generating function of the log-normal distribution using the star probability measure Yuri Heymann Received: date / Accepted: date Abstract
More informationGeometric projection of stochastic differential equations
Geometric projection of stochastic differential equations John Armstrong (King s College London) Damiano Brigo (Imperial) August 9, 2018 Idea: Projection Idea: Projection Projection gives a method of systematically
More informationLecture Note 12: Kalman Filter
ECE 645: Estimation Theory Spring 2015 Instructor: Prof. Stanley H. Chan Lecture Note 12: Kalman Filter LaTeX prepared by Stylianos Chatzidakis) May 4, 2015 This lecture note is based on ECE 645Spring
More informationLecture 4: Introduction to stochastic processes and stochastic calculus
Lecture 4: Introduction to stochastic processes and stochastic calculus Cédric Archambeau Centre for Computational Statistics and Machine Learning Department of Computer Science University College London
More informationCurrent Topics in Credit Risk
Risk Measures and Risk Management EURANDOM, Eindhoven, 10 May 2005 Current Topics in Credit Risk Mark Davis Department of Mathematics Imperial College London London SW7 2AZ www.ma.ic.ac.uk/ mdavis 1 Agenda
More informationThe Kalman Filter ImPr Talk
The Kalman Filter ImPr Talk Ged Ridgway Centre for Medical Image Computing November, 2006 Outline What is the Kalman Filter? State Space Models Kalman Filter Overview Bayesian Updating of Estimates Kalman
More information28 March Sent by to: Consultative Document: Fundamental review of the trading book 1 further response
28 March 203 Norah Barger Alan Adkins Co Chairs, Trading Book Group Basel Committee on Banking Supervision Bank for International Settlements Centralbahnplatz 2, CH 4002 Basel, SWITZERLAND Sent by email
More informationStochastic Processes. M. Sami Fadali Professor of Electrical Engineering University of Nevada, Reno
Stochastic Processes M. Sami Fadali Professor of Electrical Engineering University of Nevada, Reno 1 Outline Stochastic (random) processes. Autocorrelation. Crosscorrelation. Spectral density function.
More informationContinuum Limit of Forward Kolmogorov Equation Friday, March 06, :04 PM
Continuum Limit of Forward Kolmogorov Equation Friday, March 06, 2015 2:04 PM Please note that one of the equations (for ordinary Brownian motion) in Problem 1 was corrected on Wednesday night. And actually
More informationAdvanced Computational Methods in Statistics: Lecture 5 Sequential Monte Carlo/Particle Filtering
Advanced Computational Methods in Statistics: Lecture 5 Sequential Monte Carlo/Particle Filtering Axel Gandy Department of Mathematics Imperial College London http://www2.imperial.ac.uk/~agandy London
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationThe Scaled Unscented Transformation
The Scaled Unscented Transformation Simon J. Julier, IDAK Industries, 91 Missouri Blvd., #179 Jefferson City, MO 6519 E-mail:sjulier@idak.com Abstract This paper describes a generalisation of the unscented
More informationRobotics 2 Target Tracking. Kai Arras, Cyrill Stachniss, Maren Bennewitz, Wolfram Burgard
Robotics 2 Target Tracking Kai Arras, Cyrill Stachniss, Maren Bennewitz, Wolfram Burgard Slides by Kai Arras, Gian Diego Tipaldi, v.1.1, Jan 2012 Chapter Contents Target Tracking Overview Applications
More informationLecture 6: Bayesian Inference in SDE Models
Lecture 6: Bayesian Inference in SDE Models Bayesian Filtering and Smoothing Point of View Simo Särkkä Aalto University Simo Särkkä (Aalto) Lecture 6: Bayesian Inference in SDEs 1 / 45 Contents 1 SDEs
More informationImplicit sampling for particle filters. Alexandre Chorin, Mathias Morzfeld, Xuemin Tu, Ethan Atkins
0/20 Implicit sampling for particle filters Alexandre Chorin, Mathias Morzfeld, Xuemin Tu, Ethan Atkins University of California at Berkeley 2/20 Example: Try to find people in a boat in the middle of
More informationThe 1d Kalman Filter. 1 Understanding the forward model. Richard Turner
The d Kalman Filter Richard Turner This is a Jekyll and Hyde of a document and should really be split up. We start with Jekyll which contains a very short derivation for the d Kalman filter, the purpose
More informationExpectation, variance and moments
Expectation, variance and moments John Appleby Contents Expectation and variance Examples 3 Moments and the moment generating function 4 4 Examples of moment generating functions 5 5 Concluding remarks
More informationContagious default: application of methods of Statistical Mechanics in Finance
Contagious default: application of methods of Statistical Mechanics in Finance Wolfgang J. Runggaldier University of Padova, Italy www.math.unipd.it/runggaldier based on joint work with : Paolo Dai Pra,
More informationOptimization-Based Control
Optimization-Based Control Richard M. Murray Control and Dynamical Systems California Institute of Technology DRAFT v1.7a, 19 February 2008 c California Institute of Technology All rights reserved. This
More informationData assimilation with and without a model
Data assimilation with and without a model Tim Sauer George Mason University Parameter estimation and UQ U. Pittsburgh Mar. 5, 2017 Partially supported by NSF Most of this work is due to: Tyrus Berry,
More informationStochastic Gradient Descent in Continuous Time
Stochastic Gradient Descent in Continuous Time Justin Sirignano University of Illinois at Urbana Champaign with Konstantinos Spiliopoulos (Boston University) 1 / 27 We consider a diffusion X t X = R m
More informationComputational Aspects of Continuous-Discrete Extended Kalman-Filtering
Computational Aspects of Continuous-Discrete Extended Kalman-Filtering Thomas Mazzoni 4-3-7 Abstract This paper elaborates how the time update of the continuous-discrete extended Kalman-Filter EKF can
More informationCLOSE-TO-CLEAN REGULARIZATION RELATES
Worshop trac - ICLR 016 CLOSE-TO-CLEAN REGULARIZATION RELATES VIRTUAL ADVERSARIAL TRAINING, LADDER NETWORKS AND OTHERS Mudassar Abbas, Jyri Kivinen, Tapani Raio Department of Computer Science, School of
More informationROBOTICS 01PEEQW. Basilio Bona DAUIN Politecnico di Torino
ROBOTICS 01PEEQW Basilio Bona DAUIN Politecnico di Torino Probabilistic Fundamentals in Robotics Gaussian Filters Course Outline Basic mathematical framework Probabilistic models of mobile robots Mobile
More informationShort-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility
Short-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility José Enrique Figueroa-López 1 1 Department of Statistics Purdue University Statistics, Jump Processes,
More informationKalman Filter. Predict: Update: x k k 1 = F k x k 1 k 1 + B k u k P k k 1 = F k P k 1 k 1 F T k + Q
Kalman Filter Kalman Filter Predict: x k k 1 = F k x k 1 k 1 + B k u k P k k 1 = F k P k 1 k 1 F T k + Q Update: K = P k k 1 Hk T (H k P k k 1 Hk T + R) 1 x k k = x k k 1 + K(z k H k x k k 1 ) P k k =(I
More informationEKF, UKF. Pieter Abbeel UC Berkeley EECS. Many slides adapted from Thrun, Burgard and Fox, Probabilistic Robotics
EKF, UKF Pieter Abbeel UC Berkeley EECS Many slides adapted from Thrun, Burgard and Fox, Probabilistic Robotics Kalman Filter Kalman Filter = special case of a Bayes filter with dynamics model and sensory
More informationEKF, UKF. Pieter Abbeel UC Berkeley EECS. Many slides adapted from Thrun, Burgard and Fox, Probabilistic Robotics
EKF, UKF Pieter Abbeel UC Berkeley EECS Many slides adapted from Thrun, Burgard and Fox, Probabilistic Robotics Kalman Filter Kalman Filter = special case of a Bayes filter with dynamics model and sensory
More informationState-space Model. Eduardo Rossi University of Pavia. November Rossi State-space Model Fin. Econometrics / 53
State-space Model Eduardo Rossi University of Pavia November 2014 Rossi State-space Model Fin. Econometrics - 2014 1 / 53 Outline 1 Motivation 2 Introduction 3 The Kalman filter 4 Forecast errors 5 State
More informationAmbiguity and Information Processing in a Model of Intermediary Asset Pricing
Ambiguity and Information Processing in a Model of Intermediary Asset Pricing Leyla Jianyu Han 1 Kenneth Kasa 2 Yulei Luo 1 1 The University of Hong Kong 2 Simon Fraser University December 15, 218 1 /
More informationIf we want to analyze experimental or simulated data we might encounter the following tasks:
Chapter 1 Introduction If we want to analyze experimental or simulated data we might encounter the following tasks: Characterization of the source of the signal and diagnosis Studying dependencies Prediction
More informationThe concentration of a drug in blood. Exponential decay. Different realizations. Exponential decay with noise. dc(t) dt.
The concentration of a drug in blood Exponential decay C12 concentration 2 4 6 8 1 C12 concentration 2 4 6 8 1 dc(t) dt = µc(t) C(t) = C()e µt 2 4 6 8 1 12 time in minutes 2 4 6 8 1 12 time in minutes
More informationSensor Tasking and Control
Sensor Tasking and Control Sensing Networking Leonidas Guibas Stanford University Computation CS428 Sensor systems are about sensing, after all... System State Continuous and Discrete Variables The quantities
More informationExercises in stochastic analysis
Exercises in stochastic analysis Franco Flandoli, Mario Maurelli, Dario Trevisan The exercises with a P are those which have been done totally or partially) in the previous lectures; the exercises with
More informationAccurate approximation of stochastic differential equations
Accurate approximation of stochastic differential equations Simon J.A. Malham and Anke Wiese (Heriot Watt University, Edinburgh) Birmingham: 6th February 29 Stochastic differential equations dy t = V (y
More informationTuning of Extended Kalman Filter for nonlinear State Estimation
OSR Journal of Computer Engineering (OSR-JCE) e-ssn: 78-0661,p-SSN: 78-877, Volume 18, ssue 5, Ver. V (Sep. - Oct. 016), PP 14-19 www.iosrjournals.org Tuning of Extended Kalman Filter for nonlinear State
More informationFrom Bayes to Extended Kalman Filter
From Bayes to Extended Kalman Filter Michal Reinštein Czech Technical University in Prague Faculty of Electrical Engineering, Department of Cybernetics Center for Machine Perception http://cmp.felk.cvut.cz/
More informationToward a benchmark GPU platform to simulate XVA
Diallo - Lokman (INRIA) MonteCarlo16 1 / 22 Toward a benchmark GPU platform to simulate XVA Babacar Diallo a joint work with Lokman Abbas-Turki INRIA 6 July 2016 Diallo - Lokman (INRIA) MonteCarlo16 2
More informationState Estimation for Nonlinear Systems using Restricted Genetic Optimization
State Estimation for Nonlinear Systems using Restricted Genetic Optimization Santiago Garrido, Luis Moreno, and Carlos Balaguer Universidad Carlos III de Madrid, Leganés 28911, Madrid (Spain) Abstract.
More informationA NEW NONLINEAR FILTER
COMMUNICATIONS IN INFORMATION AND SYSTEMS c 006 International Press Vol 6, No 3, pp 03-0, 006 004 A NEW NONLINEAR FILTER ROBERT J ELLIOTT AND SIMON HAYKIN Abstract A discrete time filter is constructed
More informationEnsemble Kalman Filter
Ensemble Kalman Filter Geir Evensen and Laurent Bertino Hydro Research Centre, Bergen, Norway, Nansen Environmental and Remote Sensing Center, Bergen, Norway The Ensemble Kalman Filter (EnKF) Represents
More informationTSRT14: Sensor Fusion Lecture 8
TSRT14: Sensor Fusion Lecture 8 Particle filter theory Marginalized particle filter Gustaf Hendeby gustaf.hendeby@liu.se TSRT14 Lecture 8 Gustaf Hendeby Spring 2018 1 / 25 Le 8: particle filter theory,
More informationAutonomous Mobile Robot Design
Autonomous Mobile Robot Design Topic: Extended Kalman Filter Dr. Kostas Alexis (CSE) These slides relied on the lectures from C. Stachniss, J. Sturm and the book Probabilistic Robotics from Thurn et al.
More informationLARGE-SCALE TRAFFIC STATE ESTIMATION
Hans van Lint, Yufei Yuan & Friso Scholten A localized deterministic Ensemble Kalman Filter LARGE-SCALE TRAFFIC STATE ESTIMATION CONTENTS Intro: need for large-scale traffic state estimation Some Kalman
More informationMean-square Stability Analysis of an Extended Euler-Maruyama Method for a System of Stochastic Differential Equations
Mean-square Stability Analysis of an Extended Euler-Maruyama Method for a System of Stochastic Differential Equations Ram Sharan Adhikari Assistant Professor Of Mathematics Rogers State University Mathematical
More informationVast Volatility Matrix Estimation for High Frequency Data
Vast Volatility Matrix Estimation for High Frequency Data Yazhen Wang National Science Foundation Yale Workshop, May 14-17, 2009 Disclaimer: My opinion, not the views of NSF Y. Wang (at NSF) 1 / 36 Outline
More informationA MODEL FOR THE LONG-TERM OPTIMAL CAPACITY LEVEL OF AN INVESTMENT PROJECT
A MODEL FOR HE LONG-ERM OPIMAL CAPACIY LEVEL OF AN INVESMEN PROJEC ARNE LØKKA AND MIHAIL ZERVOS Abstract. We consider an investment project that produces a single commodity. he project s operation yields
More informationStochastic Differential Equations
Chapter 5 Stochastic Differential Equations We would like to introduce stochastic ODE s without going first through the machinery of stochastic integrals. 5.1 Itô Integrals and Itô Differential Equations
More informationStrong uniqueness for stochastic evolution equations with possibly unbounded measurable drift term
1 Strong uniqueness for stochastic evolution equations with possibly unbounded measurable drift term Enrico Priola Torino (Italy) Joint work with G. Da Prato, F. Flandoli and M. Röckner Stochastic Processes
More informationLecture 9. Time series prediction
Lecture 9 Time series prediction Prediction is about function fitting To predict we need to model There are a bewildering number of models for data we look at some of the major approaches in this lecture
More informationOn the Strong Approximation of Jump-Diffusion Processes
QUANTITATIV FINANC RSARCH CNTR QUANTITATIV FINANC RSARCH CNTR Research Paper 157 April 25 On the Strong Approximation of Jump-Diffusion Processes Nicola Bruti-Liberati and ckhard Platen ISSN 1441-81 www.qfrc.uts.edu.au
More information