Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio ( )
|
|
- Robert Murphy
- 6 years ago
- Views:
Transcription
1 Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio ( ) Etienne Tanré - Olivier Faugeras INRIA - Team Tosca November 26th, 2014 E. Tanré (INRIA - Team Tosca) Mathematical Methods for Neurosciences November 26th, / 34
2 Outline 1 Motivation 2 Complements 3 Stochastic Differential Equations 4 Link between SDE and PDE 5 Approximation of Solutions 6 Noisy Integrate and Fire Models 7 Complement: Point Poisson Processes E. Tanré (INRIA - Team Tosca) Mathematical Methods for Neurosciences November 26th, / 34
3 Outline 1 Motivation 2 Complements 3 Stochastic Differential Equations 4 Link between SDE and PDE 5 Approximation of Solutions 6 Noisy Integrate and Fire Models 7 Complement: Point Poisson Processes E. Tanré (INRIA - Team Tosca) Mathematical Methods for Neurosciences November 26th, / 34
4 Motivation First models of individual neurons Simple Integrate-and-Fire (SIF) Model Leaky Integrate-and-Fire (LIF) Model E. Tanré (INRIA - Team Tosca) Mathematical Methods for Neurosciences November 26th, / 34
5 A first model of single neuron The integrate-and-fire neuronal model was introduced by Lapicque in The membrane equation where I L (V ): the leak current I syn(v, t): the synaptic current I ext(t): the external current C: membrane capacitance dv (t) C = I L (V ) + I syn (V, t) + I ext (t), dt A spike response is generated whenever the membrane potential reaches a fixed threshold V th. After the spike, V is reset to a fixed value V reset. E. Tanré (INRIA - Team Tosca) Mathematical Methods for Neurosciences November 26th, / 34
6 Two simple examples Simple integrate-and-fire model The simplest form: the IF neuron has no leak current i.e. I L (V ) = 0. Leaky integrate-and-fire model I L (V ) = g L (V (t) E L ), where g L : the leak conductance E L : resting potential E. Tanré (INRIA - Team Tosca) Mathematical Methods for Neurosciences November 26th, / 34
7 Synaptic current First simple model of synaptic transmission An instantaneous rise or fall of the synaptic current at arrival of a presynaptic spike The PSC is described by a delta function with amplitude of efficiency J The total synaptic current stemming from N sym synaptic input channels takes the form N syn I syn (V, t) = I syn (t) = τ m J δ(t ti k ), where τ m = C. g L Asymptotic behaviour i=1 Assume the neuron receives a high barrage of Poissonian distributed and uncorrelated synaptic inputs Assume the amplitude J is small, i.e. k J << V th E L E. Tanré (INRIA - Team Tosca) Mathematical Methods for Neurosciences November 26th, / 34
8 Asymptotic behaviour Current I syn (t)dt µdt + σ τ m dw t, where µ = τ m JN syn ν syn σ 2 = τ m J 2 N syn ν syn where ν syn is the mean activation rate of each synapse. A continuous time limit equation τ m dv (t) = f (V (t))dt + I ext (t)dt + µdt + σ τ m dw t. E. Tanré (INRIA - Team Tosca) Mathematical Methods for Neurosciences November 26th, / 34
9 Simple examples f (V ) = 0, I ext = 0 τ m dv (t) = µdt + σ τ m dw t. The potential V evolves as a Brownian Motion with constant drift. f (V ) = V, I ext = 0 τ m dv (t) = (µ V (t))dt + σ τ m dw t. The potential V evolves as an Ornstein Uhlenbeck process. Spiking times Recall that the considered neuron emits a spike at each time τ its potential hits threshold V th From a mathematical viewpoint, the spiking times are the first hitting time of constant threshold by a stochastic process τ = inf {t > 0, V (t) V th } E. Tanré (INRIA - Team Tosca) Mathematical Methods for Neurosciences November 26th, / 34
10 Outline 1 Motivation 2 Complements 3 Stochastic Differential Equations 4 Link between SDE and PDE 5 Approximation of Solutions 6 Noisy Integrate and Fire Models 7 Complement: Point Poisson Processes E. Tanré (INRIA - Team Tosca) Mathematical Methods for Neurosciences November 26th, / 34
11 Levy s Characterization of Brownian Motion Theorem A stochastic process X = (X t ) t 0 is a standard Brownian Motion if and only if it is a continuous local martingale with [X] t = t. Theorem (Multi-dimensional Version) Let X = (X 1 t,, X n t ) t 0 be continuous local martingales such that [X i, X j ] t = tδ i,j. Then X is a standard n-dimensional Brownian Motion. Remark Condition in previous theorems are obviously characterization of Brownian motions. E. Tanré (INRIA - Team Tosca) Mathematical Methods for Neurosciences November 26th, / 34
12 Itô Isometry Proposition Let H be an adapted process such that Then, ( T E 0 T 0 E ( H 2 s ) ds <. ) 2 T H s dw S = E ( H 2 ) s ds 0 E. Tanré (INRIA - Team Tosca) Mathematical Methods for Neurosciences November 26th, / 34
13 Change of time Theorem (Dubins-Schwartz) Let M be a continuous locale Martingale such that M 0 = 0 a.s., and [M] =. We set τ s = inf{t 0, [M] t s}, then B s = M τs is an F τs -Brownian Motion with M t = B [M]t Remark The Dubins-Schwarz theorem, which shows that continuous martingales with unbounded (as time goes to infinity) quadratic variation ARE Brownian Motion, up to a (stochastic) time change. E. Tanré (INRIA - Team Tosca) Mathematical Methods for Neurosciences November 26th, / 34
14 Outline 1 Motivation 2 Complements 3 Stochastic Differential Equations 4 Link between SDE and PDE 5 Approximation of Solutions 6 Noisy Integrate and Fire Models 7 Complement: Point Poisson Processes E. Tanré (INRIA - Team Tosca) Mathematical Methods for Neurosciences November 26th, / 34
15 Stochastic Differential Equations dx s = b(s, X s )ds + σ(s, X s )dw s, where b and σ are predictable. Solution t t X t = X 0 + b(s, X s )ds + σ(s, X s )dw s, 0 0 E. Tanré (INRIA - Team Tosca) Mathematical Methods for Neurosciences November 26th, / 34
16 Strong Solutions A strong solution of the SDE on the given probability space (Ω, F, F, P) and with respect to the fixed F Brownian motion (W t ) t 0 and initial condition ξ is a process (X t ) t 0 with continuous sample paths and with the following properties X is adapted to the filtration F P(X 0 = ξ) = 1 for every t 0 holds almost surely. ( t ) P b(s, X s ) + σ 2 (s, X s )ds < = 1 0 t t X t = X 0 + b(s, X s )ds + σ(s, X s )dw s 0 0 E. Tanré (INRIA - Team Tosca) Mathematical Methods for Neurosciences November 26th, / 34
17 Weak Solutions A weak solution to SDE is a triple (X, W ), (ω, F, P),F where (Ω, F, P) is a probability space, and F is a filtration of sub-σ-fields of F satisfying the usual conditions. X is a continuous, F-adapted stochastic process W is an F-Brownian motion for every t 0 holds almost surely. ( t ) P b(s, X s ) + σ 2 (s, X s )ds < = 1 0 t t X t = X 0 + b(s, X s )ds + σ(s, X s )dw s 0 0 E. Tanré (INRIA - Team Tosca) Mathematical Methods for Neurosciences November 26th, / 34
18 Strong Existence and Uniqueness Theorem Let T > 0 and b(.,.) : [0, T ] R n R n and σ(.,.) : [0, T ] R n R n m be measurable functions satisfying b(t, x) + σ(t, x) C(1 + x ); x R n, t [0, T ] for some constant C and such that b(t, x) b(t, y) + σ(t, x) σ(t, y) D x y ; x, y R n, t [0, T ]. Then the stochastic differential equation dx t = b(t, X t )dt + σ(t, X t )db t, X 0 = X 0 has a unique solution X, continuous in time, such that X is adapted to the filtration generated by X 0 and the Brownian Motion and [ ] T E X t 2 dt <. 0 E. Tanré (INRIA - Team Tosca) Mathematical Methods for Neurosciences November 26th, / 34
19 Main Ideas of the Proof Uniqueness is obtained thanks to Gronwall Lemma. Existence: Picard iteration scheme. E. Tanré (INRIA - Team Tosca) Mathematical Methods for Neurosciences November 26th, / 34
20 Infinitesimal Generator Associated to a Feller Process Let X be a Feller process; a function f in C 0 is said to belong to the domain D A of the infinitesimal generator of X if the limit Af (x) = lim t 0 E x (f (X t )) f (x) t exists in C 0. The operator A : D A C 0 is called the infinitesimal generator of the process X. Example (Diffusion Processes) The infinitesimal generator associated to the solution of a Stochastic Differential Equation dx t = b(x t )dt + σ(x t )dw t writes with Af (x) = i b i (x) x i f (x) a ij (x) = 2 a ij (x) f (x), x i x j i r σ ik (x)σ kj (x). k=1 E. Tanré (INRIA - Team Tosca) Mathematical Methods for Neurosciences November 26th, / 34 j
21 Outline 1 Motivation 2 Complements 3 Stochastic Differential Equations 4 Link between SDE and PDE 5 Approximation of Solutions 6 Noisy Integrate and Fire Models 7 Complement: Point Poisson Processes E. Tanré (INRIA - Team Tosca) Mathematical Methods for Neurosciences November 26th, / 34
22 Kolmogorov s backward equation Theorem Let f C 2 0 (Rn ). Define then u(t,.) D A for each t and u(t, x) = E x [f (X t )] u t = Au, t > 0, x Rn u(0, x) = f (x), x R n. where the right hand side is to be interpreted as A applied to the function x u(t, x). E. Tanré (INRIA - Team Tosca) Mathematical Methods for Neurosciences November 26th, / 34
23 Feynman-Kac formula Theorem Let f C 2 0 (Rn ) and q C(R n ). Assume that q is lower bounded. Put ( t ) ] v(t, x) = E [exp x q(x s )ds f (X t ). 0 v = Av qv, t > 0, x Rn t v(0, x) = f (x), x R n. E. Tanré (INRIA - Team Tosca) Mathematical Methods for Neurosciences November 26th, / 34
24 Fokker Planck Equation Theorem Let X be an Itô diffusion in R n, solution of the stochastic differential equation dx t = b(x t )dt + σ(x t )dw t. Assume that P x (X t dy) = Γ(t, x, y)dy, for all x R n, t > 0. Assume that y Γ(t, x, y) is smooth for each t and x. Then, Γ satisfies the Kolmogorov forward equation (also known as Fokker Planck eq.) d dt Γ(t, x, y) = A Γ(t, x, y), where A is the adjoint of the operator A A φ(y) = i y i (b i (y)φ(y)) i j 2 y i y j [a ij (y)φ(y)]. Here a = σσ t, that is a ij (x) = r k=1 σ ik(x)σ kj (x). E. Tanré (INRIA - Team Tosca) Mathematical Methods for Neurosciences November 26th, / 34
25 Backward Kolmogorov Equation Theorem Let X be a solution of the Stochastic Differential Equation dx t = b(x t )dt + σ(x t )dw t. Denote Γ(t, x, y) = P(X t = y X 0 = x). Then, Γ Γ (t, x, y) = b(x) t x (t, x, y) a(x) 2 Γ (t, x, y) x 2 = AΓ(t, x, y), where the infinitesimal operator A acts here on the variable x and a = σσ t. E. Tanré (INRIA - Team Tosca) Mathematical Methods for Neurosciences November 26th, / 34
26 Outline 1 Motivation 2 Complements 3 Stochastic Differential Equations 4 Link between SDE and PDE 5 Approximation of Solutions 6 Noisy Integrate and Fire Models 7 Complement: Point Poisson Processes E. Tanré (INRIA - Team Tosca) Mathematical Methods for Neurosciences November 26th, / 34
27 Approximation of Solutions Euler scheme X δ 0 = X 0 X δ (k+1)δ = X δ kδ + b( X δ kδ)δ + σ( X δ kδ)(w (k+1)δ W kδ ) Theorem The numerical scheme is strongly convergent lim E ( X T X T δ ) = 0. δ 0 The numerical scheme is weakly convergent lim Eg(XT ) Eg( X T δ ) = 0. δ 0 E. Tanré (INRIA - Team Tosca) Mathematical Methods for Neurosciences November 26th, / 34
28 Rate of convergence Theorem Assume b and σ are C 4 functions with bounded derivatives, then E ( X T X T δ ) C T δ 1/2. Eg(X T ) Eg( X T δ ) C T δ. Remark: Romberg extrapolation Assume, we have an expansion of the error Eg(X T ) Eg( X T δ ) = CT 1 δ + CT 2 δ 2 + O(δ 3 ). Then a well chosen combination of X δ and X δ/2 gives an order 2 scheme. E. Tanré (INRIA - Team Tosca) Mathematical Methods for Neurosciences November 26th, / 34
29 Outline 1 Motivation 2 Complements 3 Stochastic Differential Equations 4 Link between SDE and PDE 5 Approximation of Solutions 6 Noisy Integrate and Fire Models 7 Complement: Point Poisson Processes E. Tanré (INRIA - Team Tosca) Mathematical Methods for Neurosciences November 26th, / 34
30 Simple Integrate-and-Fire Model Special case I ext = 0 dv t = σdw t. Exercise: Law of the first hitting time Compute the law of ( P (τ a t) = P sup 0 s t τ a := inf {t > 0, W t a}. ) W s a ( ) ( ) = P sup 0 s t W s a, W t a ( + P sup 0 s t W s a, W t < a ) = P (W t a) + P sup W s a, W t W τa < 0 0 s t ( ) = P (W t a) + P sup 0 s t W s a, W t W τa > 0 = 2P (W t a) E. Tanré (INRIA - Team Tosca) Mathematical Methods for Neurosciences November 26th, / 34
31 Ornstein Uhlenbeck Process dv t = λ(v V t )dt + σdw t V 0 = v Solve explicitly the equation Give the law of V t (the conditional law given V 0 ) Make explicit the associated stationary measure. E. Tanré (INRIA - Team Tosca) Mathematical Methods for Neurosciences November 26th, / 34
32 Outline 1 Motivation 2 Complements 3 Stochastic Differential Equations 4 Link between SDE and PDE 5 Approximation of Solutions 6 Noisy Integrate and Fire Models 7 Complement: Point Poisson Processes E. Tanré (INRIA - Team Tosca) Mathematical Methods for Neurosciences November 26th, / 34
33 Point Poisson Process (P.P.P.) Let us consider a set D (e.g. [0, T ] [0, K]). A realisation of a P.P.P. on D with intensity I(t, x) is a set of points of D. For all subset F of D, we denote by N F the number of points of the P.P.P. which are in F. Characterization of a Point Poisson Process for all F D, N F is a random variable (with value in N) with Poisson law of parameter I(t, x)dtdx, F for all subset F and G with empty intersection, N F and N G are independent random variables. A particular case is: the intensity I is equal to 1. For all subset F of D, the number N F of points of the P.P.P. in F is a Poisson random variable with parameter equal to the volume of F. E. Tanré (INRIA - Team Tosca) Mathematical Methods for Neurosciences November 26th, / 34
34 Point Poisson Process Recall that a Poisson random variable X of parameter λ has the following law: P(X = k) = exp( λ) λk k!. In particular, P(X = 0) = exp( λ). Let us consider a process (X t ; 0 t T ) and a non negative function φ, bounded from above by K. In order to simulate an event with probability exp( T 0 φ(x s)ds), we can use Point Poisson Processes. Indeed, the probability that the P.P.P. on [0, T ] [0, K] have no point below the curve t φ(x t ) is precisely exp( T 0 φ(x s)ds). E. Tanré (INRIA - Team Tosca) Mathematical Methods for Neurosciences November 26th, / 34
p 1 ( Y p dp) 1/p ( X p dp) 1 1 p
Doob s inequality Let X(t) be a right continuous submartingale with respect to F(t), t 1 P(sup s t X(s) λ) 1 λ {sup s t X(s) λ} X + (t)dp 2 For 1 < p
More informationMathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio ( )
Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio (2014-2015) Etienne Tanré - Olivier Faugeras INRIA - Team Tosca October 22nd, 2014 E. Tanré (INRIA - Team Tosca) Mathematical
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 informationStochastic Calculus February 11, / 33
Martingale Transform M n martingale with respect to F n, n =, 1, 2,... σ n F n (σ M) n = n 1 i= σ i(m i+1 M i ) is a Martingale E[(σ M) n F n 1 ] n 1 = E[ σ i (M i+1 M i ) F n 1 ] i= n 2 = σ i (M i+1 M
More information1. Stochastic Process
HETERGENEITY IN QUANTITATIVE MACROECONOMICS @ TSE OCTOBER 17, 216 STOCHASTIC CALCULUS BASICS SANG YOON (TIM) LEE Very simple notes (need to add references). It is NOT meant to be a substitute for a real
More informationStochastic Differential Equations.
Chapter 3 Stochastic Differential Equations. 3.1 Existence and Uniqueness. One of the ways of constructing a Diffusion process is to solve the stochastic differential equation dx(t) = σ(t, x(t)) dβ(t)
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 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 informationA Concise Course on Stochastic Partial Differential Equations
A Concise Course on Stochastic Partial Differential Equations Michael Röckner Reference: C. Prevot, M. Röckner: Springer LN in Math. 1905, Berlin (2007) And see the references therein for the original
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 informationKolmogorov Equations and Markov Processes
Kolmogorov Equations and Markov Processes May 3, 013 1 Transition measures and functions Consider a stochastic process {X(t)} t 0 whose state space is a product of intervals contained in R n. We define
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 informationExercises. T 2T. e ita φ(t)dt.
Exercises. Set #. Construct an example of a sequence of probability measures P n on R which converge weakly to a probability measure P but so that the first moments m,n = xdp n do not converge to m = xdp.
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 informationHJB equations. Seminar in Stochastic Modelling in Economics and Finance January 10, 2011
Department of Probability and Mathematical Statistics Faculty of Mathematics and Physics, Charles University in Prague petrasek@karlin.mff.cuni.cz Seminar in Stochastic Modelling in Economics and Finance
More informationTheoretical Tutorial Session 2
1 / 36 Theoretical Tutorial Session 2 Xiaoming Song Department of Mathematics Drexel University July 27, 216 Outline 2 / 36 Itô s formula Martingale representation theorem Stochastic differential equations
More informationPoisson Jumps in Credit Risk Modeling: a Partial Integro-differential Equation Formulation
Poisson Jumps in Credit Risk Modeling: a Partial Integro-differential Equation Formulation Jingyi Zhu Department of Mathematics University of Utah zhu@math.utah.edu Collaborator: Marco Avellaneda (Courant
More information2012 NCTS Workshop on Dynamical Systems
Barbara Gentz gentz@math.uni-bielefeld.de http://www.math.uni-bielefeld.de/ gentz 2012 NCTS Workshop on Dynamical Systems National Center for Theoretical Sciences, National Tsing-Hua University Hsinchu,
More informationIntroduction to Random Diffusions
Introduction to Random Diffusions The main reason to study random diffusions is that this class of processes combines two key features of modern probability theory. On the one hand they are semi-martingales
More informationBrownian Motion. 1 Definition Brownian Motion Wiener measure... 3
Brownian Motion Contents 1 Definition 2 1.1 Brownian Motion................................. 2 1.2 Wiener measure.................................. 3 2 Construction 4 2.1 Gaussian process.................................
More informationStability of Stochastic Differential Equations
Lyapunov stability theory for ODEs s Stability of Stochastic Differential Equations Part 1: Introduction Department of Mathematics and Statistics University of Strathclyde Glasgow, G1 1XH December 2010
More informationThe Smoluchowski-Kramers Approximation: What model describes a Brownian particle?
The Smoluchowski-Kramers Approximation: What model describes a Brownian particle? Scott Hottovy shottovy@math.arizona.edu University of Arizona Applied Mathematics October 7, 2011 Brown observes a particle
More informationLecture 4: Ito s Stochastic Calculus and SDE. Seung Yeal Ha Dept of Mathematical Sciences Seoul National University
Lecture 4: Ito s Stochastic Calculus and SDE Seung Yeal Ha Dept of Mathematical Sciences Seoul National University 1 Preliminaries What is Calculus? Integral, Differentiation. Differentiation 2 Integral
More informationNonlinear representation, backward SDEs, and application to the Principal-Agent problem
Nonlinear representation, backward SDEs, and application to the Principal-Agent problem Ecole Polytechnique, France April 4, 218 Outline The Principal-Agent problem Formulation 1 The Principal-Agent problem
More informationControlled Diffusions and Hamilton-Jacobi Bellman Equations
Controlled Diffusions and Hamilton-Jacobi Bellman Equations Emo Todorov Applied Mathematics and Computer Science & Engineering University of Washington Winter 2014 Emo Todorov (UW) AMATH/CSE 579, Winter
More informationSome SDEs with distributional drift Part I : General calculus. Flandoli, Franco; Russo, Francesco; Wolf, Jochen
Title Author(s) Some SDEs with distributional drift Part I : General calculus Flandoli, Franco; Russo, Francesco; Wolf, Jochen Citation Osaka Journal of Mathematics. 4() P.493-P.54 Issue Date 3-6 Text
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 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 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 informationOn continuous time contract theory
Ecole Polytechnique, France Journée de rentrée du CMAP, 3 octobre, 218 Outline 1 2 Semimartingale measures on the canonical space Random horizon 2nd order backward SDEs (Static) Principal-Agent Problem
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 informationSTATISTICS 385: STOCHASTIC CALCULUS HOMEWORK ASSIGNMENT 4 DUE NOVEMBER 23, = (2n 1)(2n 3) 3 1.
STATISTICS 385: STOCHASTIC CALCULUS HOMEWORK ASSIGNMENT 4 DUE NOVEMBER 23, 26 Problem Normal Moments (A) Use the Itô formula and Brownian scaling to check that the even moments of the normal distribution
More informationA Short Introduction to Diffusion Processes and Ito Calculus
A Short Introduction to Diffusion Processes and Ito Calculus Cédric Archambeau University College, London Center for Computational Statistics and Machine Learning c.archambeau@cs.ucl.ac.uk January 24,
More informationHomogenization with stochastic differential equations
Homogenization with stochastic differential equations Scott Hottovy shottovy@math.arizona.edu University of Arizona Program in Applied Mathematics October 12, 2011 Modeling with SDE Use SDE to model system
More informationPartial Differential Equations with Applications to Finance Seminar 1: Proving and applying Dynkin s formula
Partial Differential Equations with Applications to Finance Seminar 1: Proving and applying Dynkin s formula Group 4: Bertan Yilmaz, Richard Oti-Aboagye and Di Liu May, 15 Chapter 1 Proving Dynkin s formula
More informationKolmogorov equations in Hilbert spaces IV
March 26, 2010 Other types of equations Let us consider the Burgers equation in = L 2 (0, 1) dx(t) = (AX(t) + b(x(t))dt + dw (t) X(0) = x, (19) where A = ξ 2, D(A) = 2 (0, 1) 0 1 (0, 1), b(x) = ξ 2 (x
More informationWalsh Diffusions. Andrey Sarantsev. March 27, University of California, Santa Barbara. Andrey Sarantsev University of Washington, Seattle 1 / 1
Walsh Diffusions Andrey Sarantsev University of California, Santa Barbara March 27, 2017 Andrey Sarantsev University of Washington, Seattle 1 / 1 Walsh Brownian Motion on R d Spinning measure µ: probability
More informationGaussian processes for inference in stochastic differential equations
Gaussian processes for inference in stochastic differential equations Manfred Opper, AI group, TU Berlin November 6, 2017 Manfred Opper, AI group, TU Berlin (TU Berlin) inference in SDE November 6, 2017
More informationBrownian motion. Samy Tindel. Purdue University. Probability Theory 2 - MA 539
Brownian motion Samy Tindel Purdue University Probability Theory 2 - MA 539 Mostly taken from Brownian Motion and Stochastic Calculus by I. Karatzas and S. Shreve Samy T. Brownian motion Probability Theory
More informationConvergence of Particle Filtering Method for Nonlinear Estimation of Vortex Dynamics
Convergence of Particle Filtering Method for Nonlinear Estimation of Vortex Dynamics Meng Xu Department of Mathematics University of Wyoming February 20, 2010 Outline 1 Nonlinear Filtering Stochastic Vortex
More informationElliptic Operators with Unbounded Coefficients
Elliptic Operators with Unbounded Coefficients Federica Gregorio Universitá degli Studi di Salerno 8th June 2018 joint work with S.E. Boutiah, A. Rhandi, C. Tacelli Motivation Consider the Stochastic Differential
More informationDiscretization of Stochastic Differential Systems With Singular Coefficients Part II
Discretization of Stochastic Differential Systems With Singular Coefficients Part II Denis Talay, INRIA Sophia Antipolis joint works with Mireille Bossy, Nicolas Champagnat, Sylvain Maire, Miguel Martinez,
More informationStochastic differential equations in neuroscience
Stochastic differential equations in neuroscience Nils Berglund MAPMO, Orléans (CNRS, UMR 6628) http://www.univ-orleans.fr/mapmo/membres/berglund/ Barbara Gentz, Universität Bielefeld Damien Landon, MAPMO-Orléans
More informationMA8109 Stochastic Processes in Systems Theory Autumn 2013
Norwegian University of Science and Technology Department of Mathematical Sciences MA819 Stochastic Processes in Systems Theory Autumn 213 1 MA819 Exam 23, problem 3b This is a linear equation of the form
More informationDynamical systems with Gaussian and Levy noise: analytical and stochastic approaches
Dynamical systems with Gaussian and Levy noise: analytical and stochastic approaches Noise is often considered as some disturbing component of the system. In particular physical situations, noise becomes
More informationExact Simulation of Diffusions and Jump Diffusions
Exact Simulation of Diffusions and Jump Diffusions A work by: Prof. Gareth O. Roberts Dr. Alexandros Beskos Dr. Omiros Papaspiliopoulos Dr. Bruno Casella 28 th May, 2008 Content 1 Exact Algorithm Construction
More informationOn a class of stochastic differential equations in a financial network model
1 On a class of stochastic differential equations in a financial network model Tomoyuki Ichiba Department of Statistics & Applied Probability, Center for Financial Mathematics and Actuarial Research, University
More informationContinuous dependence estimates for the ergodic problem with an application to homogenization
Continuous dependence estimates for the ergodic problem with an application to homogenization Claudio Marchi Bayreuth, September 12 th, 2013 C. Marchi (Università di Padova) Continuous dependence Bayreuth,
More informationBernardo D Auria Stochastic Processes /12. Notes. March 29 th, 2012
1 Stochastic Calculus Notes March 9 th, 1 In 19, Bachelier proposed for the Paris stock exchange a model for the fluctuations affecting the price X(t) of an asset that was given by the Brownian motion.
More informationlim n C1/n n := ρ. [f(y) f(x)], y x =1 [f(x) f(y)] [g(x) g(y)]. (x,y) E A E(f, f),
1 Part I Exercise 1.1. Let C n denote the number of self-avoiding random walks starting at the origin in Z of length n. 1. Show that (Hint: Use C n+m C n C m.) lim n C1/n n = inf n C1/n n := ρ.. Show that
More informationMATH 56A SPRING 2008 STOCHASTIC PROCESSES 197
MATH 56A SPRING 8 STOCHASTIC PROCESSES 197 9.3. Itô s formula. First I stated the theorem. Then I did a simple example to make sure we understand what it says. Then I proved it. The key point is Lévy s
More informationStochastic Areas and Applications in Risk Theory
Stochastic Areas and Applications in Risk Theory July 16th, 214 Zhenyu Cui Department of Mathematics Brooklyn College, City University of New York Zhenyu Cui 49th Actuarial Research Conference 1 Outline
More information1. Stochastic Processes and filtrations
1. Stochastic Processes and 1. Stoch. pr., A stochastic process (X t ) t T is a collection of random variables on (Ω, F) with values in a measurable space (S, S), i.e., for all t, In our case X t : Ω S
More informationThe Wiener Itô Chaos Expansion
1 The Wiener Itô Chaos Expansion The celebrated Wiener Itô chaos expansion is fundamental in stochastic analysis. In particular, it plays a crucial role in the Malliavin calculus as it is presented in
More informationLecture 12: Diffusion Processes and Stochastic Differential Equations
Lecture 12: Diffusion Processes and Stochastic Differential Equations 1. Diffusion Processes 1.1 Definition of a diffusion process 1.2 Examples 2. Stochastic Differential Equations SDE) 2.1 Stochastic
More informationRegularization by noise in infinite dimensions
Regularization by noise in infinite dimensions Franco Flandoli, University of Pisa King s College 2017 Franco Flandoli, University of Pisa () Regularization by noise King s College 2017 1 / 33 Plan of
More informationSome Tools From Stochastic Analysis
W H I T E Some Tools From Stochastic Analysis J. Potthoff Lehrstuhl für Mathematik V Universität Mannheim email: potthoff@math.uni-mannheim.de url: http://ls5.math.uni-mannheim.de To close the file, click
More informationAlbert N. Shiryaev Steklov Mathematical Institute. On sharp maximal inequalities for stochastic processes
Albert N. Shiryaev Steklov Mathematical Institute On sharp maximal inequalities for stochastic processes joint work with Yaroslav Lyulko, Higher School of Economics email: albertsh@mi.ras.ru 1 TOPIC I:
More informationBackward martingale representation and endogenous completeness in finance
Backward martingale representation and endogenous completeness in finance Dmitry Kramkov (with Silviu Predoiu) Carnegie Mellon University 1 / 19 Bibliography Robert M. Anderson and Roberto C. Raimondo.
More informationSolution for Problem 7.1. We argue by contradiction. If the limit were not infinite, then since τ M (ω) is nondecreasing we would have
362 Problem Hints and Solutions sup g n (ω, t) g(ω, t) sup g(ω, s) g(ω, t) µ n (ω). t T s,t: s t 1/n By the uniform continuity of t g(ω, t) on [, T], one has for each ω that µ n (ω) as n. Two applications
More informationOn the martingales obtained by an extension due to Saisho, Tanemura and Yor of Pitman s theorem
On the martingales obtained by an extension due to Saisho, Tanemura and Yor of Pitman s theorem Koichiro TAKAOKA Dept of Applied Physics, Tokyo Institute of Technology Abstract M Yor constructed a family
More informationBrownian 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 informationLecture 21 Representations of Martingales
Lecture 21: Representations of Martingales 1 of 11 Course: Theory of Probability II Term: Spring 215 Instructor: Gordan Zitkovic Lecture 21 Representations of Martingales Right-continuous inverses Let
More informationWeek 9 Generators, duality, change of measure
Week 9 Generators, duality, change of measure Jonathan Goodman November 18, 013 1 Generators This section describes a common abstract way to describe many of the differential equations related to Markov
More informationBranching Processes II: Convergence of critical branching to Feller s CSB
Chapter 4 Branching Processes II: Convergence of critical branching to Feller s CSB Figure 4.1: Feller 4.1 Birth and Death Processes 4.1.1 Linear birth and death processes Branching processes can be studied
More informationBridging the Gap between Center and Tail for Multiscale Processes
Bridging the Gap between Center and Tail for Multiscale Processes Matthew R. Morse Department of Mathematics and Statistics Boston University BU-Keio 2016, August 16 Matthew R. Morse (BU) Moderate Deviations
More informationSolutions to the Exercises in Stochastic Analysis
Solutions to the Exercises in Stochastic Analysis Lecturer: Xue-Mei Li 1 Problem Sheet 1 In these solution I avoid using conditional expectations. But do try to give alternative proofs once we learnt conditional
More informationFunctional Limit theorems for the quadratic variation of a continuous time random walk and for certain stochastic integrals
Functional Limit theorems for the quadratic variation of a continuous time random walk and for certain stochastic integrals Noèlia Viles Cuadros BCAM- Basque Center of Applied Mathematics with Prof. Enrico
More informationWEYL S LEMMA, ONE OF MANY. Daniel W. Stroock
WEYL S LEMMA, ONE OF MANY Daniel W Stroock Abstract This note is a brief, and somewhat biased, account of the evolution of what people working in PDE s call Weyl s Lemma about the regularity of solutions
More informationSynchrony in Stochastic Pulse-coupled Neuronal Network Models
Synchrony in Stochastic Pulse-coupled Neuronal Network Models Katie Newhall Gregor Kovačič and Peter Kramer Aaditya Rangan and David Cai 2 Rensselaer Polytechnic Institute, Troy, New York 2 Courant Institute,
More informationUniformly Uniformly-ergodic Markov chains and BSDEs
Uniformly Uniformly-ergodic Markov chains and BSDEs Samuel N. Cohen Mathematical Institute, University of Oxford (Based on joint work with Ying Hu, Robert Elliott, Lukas Szpruch) Centre Henri Lebesgue,
More informationStochastic Differential Equations
CHAPTER 1 Stochastic Differential Equations Consider a stochastic process X t satisfying dx t = bt, X t,w t dt + σt, X t,w t dw t. 1.1 Question. 1 Can we obtain the existence and uniqueness theorem for
More informationBernardo D Auria Stochastic Processes /10. Notes. Abril 13 th, 2010
1 Stochastic Calculus Notes Abril 13 th, 1 As we have seen in previous lessons, the stochastic integral with respect to the Brownian motion shows a behavior different from the classical Riemann-Stieltjes
More informationProperties of an infinite dimensional EDS system : the Muller s ratchet
Properties of an infinite dimensional EDS system : the Muller s ratchet LATP June 5, 2011 A ratchet source : wikipedia Plan 1 Introduction : The model of Haigh 2 3 Hypothesis (Biological) : The population
More informationWeak solutions of mean-field stochastic differential equations
Weak solutions of mean-field stochastic differential equations Juan Li School of Mathematics and Statistics, Shandong University (Weihai), Weihai 26429, China. Email: juanli@sdu.edu.cn Based on joint works
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 informationMost Probable Escape Path Method and its Application to Leaky Integrate And Fire Neurons
Most Probable Escape Path Method and its Application to Leaky Integrate And Fire Neurons Simon Fugmann Humboldt University Berlin 13/02/06 Outline The Most Probable Escape Path (MPEP) Motivation General
More informationSTOCHASTIC CALCULUS JASON MILLER AND VITTORIA SILVESTRI
STOCHASTIC CALCULUS JASON MILLER AND VITTORIA SILVESTRI Contents Preface 1 1. Introduction 1 2. Preliminaries 4 3. Local martingales 1 4. The stochastic integral 16 5. Stochastic calculus 36 6. Applications
More informationWeak convergence and large deviation theory
First Prev Next Go To Go Back Full Screen Close Quit 1 Weak convergence and large deviation theory Large deviation principle Convergence in distribution The Bryc-Varadhan theorem Tightness and Prohorov
More informationLower Tail Probabilities and Normal Comparison Inequalities. In Memory of Wenbo V. Li s Contributions
Lower Tail Probabilities and Normal Comparison Inequalities In Memory of Wenbo V. Li s Contributions Qi-Man Shao The Chinese University of Hong Kong Lower Tail Probabilities and Normal Comparison Inequalities
More informationLecture 12. F o s, (1.1) F t := s>t
Lecture 12 1 Brownian motion: the Markov property Let C := C(0, ), R) be the space of continuous functions mapping from 0, ) to R, in which a Brownian motion (B t ) t 0 almost surely takes its value. Let
More informationSimulation methods for stochastic models in chemistry
Simulation methods for stochastic models in chemistry David F. Anderson anderson@math.wisc.edu Department of Mathematics University of Wisconsin - Madison SIAM: Barcelona June 4th, 21 Overview 1. Notation
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 informationExample 4.1 Let X be a random variable and f(t) a given function of time. Then. Y (t) = f(t)x. Y (t) = X sin(ωt + δ)
Chapter 4 Stochastic Processes 4. Definition In the previous chapter we studied random variables as functions on a sample space X(ω), ω Ω, without regard to how these might depend on parameters. We now
More informationReflected Brownian Motion
Chapter 6 Reflected Brownian Motion Often we encounter Diffusions in regions with boundary. If the process can reach the boundary from the interior in finite time with positive probability we need to decide
More information1 Brownian Local Time
1 Brownian Local Time We first begin by defining the space and variables for Brownian local time. Let W t be a standard 1-D Wiener process. We know that for the set, {t : W t = } P (µ{t : W t = } = ) =
More informationON THE POLICY IMPROVEMENT ALGORITHM IN CONTINUOUS TIME
ON THE POLICY IMPROVEMENT ALGORITHM IN CONTINUOUS TIME SAUL D. JACKA AND ALEKSANDAR MIJATOVIĆ Abstract. We develop a general approach to the Policy Improvement Algorithm (PIA) for stochastic control problems
More informationSMSTC (2007/08) Probability.
SMSTC (27/8) Probability www.smstc.ac.uk Contents 12 Markov chains in continuous time 12 1 12.1 Markov property and the Kolmogorov equations.................... 12 2 12.1.1 Finite state space.................................
More informationA Lévy-Fokker-Planck equation: entropies and convergence to equilibrium
1/ 22 A Lévy-Fokker-Planck equation: entropies and convergence to equilibrium I. Gentil CEREMADE, Université Paris-Dauphine International Conference on stochastic Analysis and Applications Hammamet, Tunisia,
More informationExact Simulation of Multivariate Itô Diffusions
Exact Simulation of Multivariate Itô Diffusions Jose Blanchet Joint work with Fan Zhang Columbia and Stanford July 7, 2017 Jose Blanchet (Columbia/Stanford) Exact Simulation of Diffusions July 7, 2017
More informationLAN property for sde s with additive fractional noise and continuous time observation
LAN property for sde s with additive fractional noise and continuous time observation Eulalia Nualart (Universitat Pompeu Fabra, Barcelona) joint work with Samy Tindel (Purdue University) Vlad s 6th birthday,
More informationMultilevel Monte Carlo for Stochastic McKean-Vlasov Equations
Multilevel Monte Carlo for Stochastic McKean-Vlasov Equations Lukasz Szpruch School of Mathemtics University of Edinburgh joint work with Shuren Tan and Alvin Tse (Edinburgh) Lukasz Szpruch (University
More information(2m)-TH MEAN BEHAVIOR OF SOLUTIONS OF STOCHASTIC DIFFERENTIAL EQUATIONS UNDER PARAMETRIC PERTURBATIONS
(2m)-TH MEAN BEHAVIOR OF SOLUTIONS OF STOCHASTIC DIFFERENTIAL EQUATIONS UNDER PARAMETRIC PERTURBATIONS Svetlana Janković and Miljana Jovanović Faculty of Science, Department of Mathematics, University
More informationOptimal Stopping and Maximal Inequalities for Poisson Processes
Optimal Stopping and Maximal Inequalities for Poisson Processes D.O. Kramkov 1 E. Mordecki 2 September 10, 2002 1 Steklov Mathematical Institute, Moscow, Russia 2 Universidad de la República, Montevideo,
More informationFirst passage time for Brownian motion and piecewise linear boundaries
To appear in Methodology and Computing in Applied Probability, (2017) 19: 237-253. doi 10.1007/s11009-015-9475-2 First passage time for Brownian motion and piecewise linear boundaries Zhiyong Jin 1 and
More informationAnnealed Brownian motion in a heavy tailed Poissonian potential
Annealed Brownian motion in a heavy tailed Poissonian potential Ryoki Fukushima Research Institute of Mathematical Sciences Stochastic Analysis and Applications, Okayama University, September 26, 2012
More informationQuantifying Intermittent Transport in Cell Cytoplasm
Quantifying Intermittent Transport in Cell Cytoplasm Ecole Normale Supérieure, Mathematics and Biology Department. Paris, France. May 19 th 2009 Cellular Transport Introduction Cellular Transport Intermittent
More informationMATH 425, HOMEWORK 3 SOLUTIONS
MATH 425, HOMEWORK 3 SOLUTIONS Exercise. (The differentiation property of the heat equation In this exercise, we will use the fact that the derivative of a solution to the heat equation again solves the
More informationIntroduction to multiscale modeling and simulation. Explicit methods for ODEs : forward Euler. y n+1 = y n + tf(y n ) dy dt = f(y), y(0) = y 0
Introduction to multiscale modeling and simulation Lecture 5 Numerical methods for ODEs, SDEs and PDEs The need for multiscale methods Two generic frameworks for multiscale computation Explicit methods
More informationSimulation of conditional diffusions via forward-reverse stochastic representations
Weierstrass Institute for Applied Analysis and Stochastics Simulation of conditional diffusions via forward-reverse stochastic representations Christian Bayer and John Schoenmakers Numerical methods for
More information