A Stochastic Paradox for Reflected Brownian Motion?
|
|
- Noel Campbell
- 6 years ago
- Views:
Transcription
1 Proceedings of the 9th International Symposium on Mathematical Theory of Networks and Systems MTNS 2 9 July, 2 udapest, Hungary A Stochastic Parado for Reflected rownian Motion? Erik I. Verriest Abstract This paper is pedagogical in nature. It is shown that a stochastic integral with respect to a Wiener process does not always yield a martingale. Consequently, ignoring the dw term in the calculation of epectations, as might be done in a first back-of-the-envelope approach, can lead to false results. We develop a model for rownian motion constrained to >, for which the moments can be computed eactly. There is no parado if one solves the problem correctly, however a blind application of the Itô calculus yields to the parado that the variance may become negative. In order to put the student back on track, we close the paper by stating some sufficient conditions that guarantee that the stochastic integral is a true martingale. I. INTRODUCTION The objective of this paper is to point out to the student that when one uses the Itô calculus, one should be careful about the interpretation of the Itô differentials, and handle the process of echanging the order of taking mathematical epectations and differentiation or integration with care. Typically, one of the first eamples the student may encounter is the scalar linear system with additive noise, in Itô form given by the equation d = a + bdw Applying the Itô differential rule, it is derived that d 2 = 2a 2 + b 2 + 2bdw. 2 Taking epectations gives respectively and Ed = ae E d 2 = 2aE 2 + b 2. one may conclude that de de 2 = ae and = 2aE 2 +b 2 leading to the well known mean and variance differential equations. ṁ = am P = 2aP + b 2. At this point, the student may epect that E df = d E f, and proceed accordingly in similar problems. The problem of reflected rownian motion shows otherwise. First we build a suitable model for reflected rownian motion in Section II and show some of its properties regarding the eit behavior and show that the process is well defined in Sections III and IV. We compute carelessly, I should say the variance in V and find that it gets negative: a parado! Then we resolve the parado in Section VI. E.I. Verriest is with the School of Electrical and Computer Engineering, Georgia Institute of Technology, Atlanta, GA 2-2, USA erik.verriest@ece.gatech.edu II. A MODEL FOR REFLECTED ROWNIAN MOTION In this section we build a model to approimate a reflected rownian motion by a process which is mathematically more tractable. It is assumed that the process starts deterministically at >. If t is large, it should behave as a pure rownian motion, but when t approaches zero, a potential function should push it back towards the positive side. Let us choose for this potential. This process obeys the Itô equation d = + dw. This does not quite satisfy the Lipschitz conditions, but something can be done by taking ǫ, as the domain of validity, and using various bounding arguments. For now, I shall neglect that issue, and proceed determining some properties. Taking epectations of the above equation yields de = E. 6 Net, we derive the Itô differential for the process /. d = 2 d d = 2 + dw + = 2 dw. Now, taking epectations yields d = 7 The initial condition is =. Hence, for all t > we get E =. 8 t Consequently, from which de = 9 E t = + t. Intuitively, this seems plausible, the effect of the nonlocal potential gives a drift towards the right increasing. We ll sharpen our intuition about this process by determining the average eit time and eit probabilities in the net sections. ISN
2 E. I. Verriest A Stochastic Parado in a Model for Reflected rownian Motion? III. EXPECTED EXIT TIME Consider the closed interval [β, ], with < β < <. The epected first passage time is given as the solution to the equation [] u + 2 u = with the boundary conditions uβ = u =. We get u + 2u = 2 u + u = 2 u + u βu β uβ = β 2 2. The left hand side contains With uβ =, we get [u] = u + u. u βuβ [βu β + uβ] β = β β 2 β u = β β 2 β + βu β β. From u =, it follows that Hence, and thus finally βu β = 2 + β + β 2 β 2. u = u β = β β β + + β 2 or, letting τ [β,] be the eit time for the process t, starting at =, from the interval [β, ], of course with β < < Eτ [β,] = β + + β Also this seems intuitively clear. If is either β or, eit is instantaneous. Taking the limit for β, Eτ [,] = 2 2. If, the eit time approaches infinity, no matter how small is. This reinforces the intuition that eit at zero is not possible. IV. EXIT PROAILITIES Consider again the closed interval [β, ], with < β < <. The probability of eit at is given by the solution to the equation [] v + 2 v = with the boundary conditions vβ = and v =. We get v + v = v + v = A. for some constant A. The left hand side is [v] = A v = A + a, where a is another constant. With vβ =, we get a = Aβ, and from v = = A + a = A β, it follows that A = β and thus finally It follows that v = β Pr eit at = β. 6 β β. 7 If β, then Pr eit at. Hence the probability of the process crossing over = is zero. Also this seems intuitively plausible. V. VARIANCE PARADOX In this section the variance for this approimation to the reflected rownian motion is computed. Since the process starts deterministically at, the initial variance is zero. Consider now the process 2 t. With the Itô-rule, Taking epectations d 2 = 2d + = 2 + dw + = + 2dw. 8 de 2 =. Integrating, this yields eplicitly E 2 t = 2 + t. We also have the deterministic equation E t 2 = t
3 Proceedings of the 9th International Symposium on Mathematical Theory of Networks and Systems MTNS 2 9 July, 2 udapest, Hungary The variance is obtained from Var t = E 2 t E t 2 = 2 + t t 2 = t t For t > 2, this variance is negative, a parado! What went wrong? VI. RESOLUTION OF THE PARADOX Consider the d-dimensional rownian motion, where the coordinates of the sample paths are independent standard one-dimensional rownian motions. It is known [2, p.9] that if rt = t, then Pr { rt dr R s } = gt s, rs, rdr, where gt, a, b = t ep a2 + b 2 ab ab d 2 I d 2 b d t and I d/2 is the modified essel function. The probability is conditioned on the filtration R s = {rτ τ s}. This gt, a, b is also the fundamental solution of u t = 2 with boundary condition 2 u b 2 + d b lim b d u b b =. u b, The corresponding Itô equation for the process rt is drt = d + dwt. 2rt Note that for d =, the three-dimensional rownian motion, this gives drt = + dwt, 2 rt which is precisely the equation we were investigating. Equation 2 and its connection with the essel process provides now a different avenue to resolve the parado described in section V. Consider thus the three-dimensional rownian motion M, starting at the point with coordinates,,, where = r. Let rt be the distance of the M from the origin. The stochastics of rt, the distance of the rownian particle from the origin, is identical to the -D dynamics on the halfline for the particle in a / potential. ut this -dimensional problem is easily solved by elementary probability []. The easiest way is to obtain first the density of r 2. From the properties of -D rownian motion, we know that t is normally distributed with mean r and variance t standard case. Likewise, both yt and zt are normally distributed with mean zero and variance t. These three processes are independent. Hence ut = y 2 t+z 2 t has a χ-square distribution with two degrees of freedom. The density is f ut u = ep u Hu, 2 where H is the Heaviside step function. The random variable v = 2 t has density With f vt v = 2 f v + f v Hv. 22 v f t = 2πt ep 2, 2 this gives ep v r 2 f vt v = Hv. 2 v 2πt + erf r Independence of ut and vt results in the convolution formula for the density of their sum r 2 t = ut + vt., f r 2 tρ = = ρ f vt ρ sf ut sds ep ρ r 2 2 /2 r πt ep ρ + r 2 + Hρ. 2 With this density, the epected value of r 2 t is readily computed E r 2 t = ρf r2 tρdρ = r 2 + t. 26 Likewise, we obtain the epected value of rt by E rt = = r π ρfr 2 tρdρ r 2 r ep + + πr 2 + terf r. 27 Finally, the inverse distance, /rt has average value E = f r rt ρ 2 tρdρ = r erf. 28 r Unless t =, this is strictly less than /r. Figure r displays E rt and it can be seen that the naive prediction equality to only holds approimately for small values of t and. It is now easily verified that for the process we considered in Section II, we get indeed, de = E,
4 E. I. Verriest A Stochastic Parado in a Model for Reflected rownian Motion? r~ 2 t~ t~ Fig.. The process r rt Fig. 2. Variance of rt but which is not zero! de = ep 2, 2πt Finally, the correct epression of the variance follows from Varrt = r 2 t rt 2. Figure 2 displays the variance for different values of as function of time. For, the variance has slope 8/π, whereas for we get slope as epected standard rownian motion. The curves in the figure correspond to =., 2,,,, 2 and. The figure below shows the average position as function of time together with the -sigma bounds for the uncertainty σ = t for initialization at =, and 2. VII. FURTHER THOUGHTS The origin of the problem is that in the differential Itô equation we get correctly E d =. However, this does not imply that de =. The integrated form is: E = t E t dws s 2. More detailed arguments, invoking stopping times, show that /mint, τ k approaches /t a.e. as k +, however the convergence is neither monotone nor dominated 2 2 t~ Fig.. Evolution of mean and uncertainty -σ bounds so the only thing one can do to pass to the limit in the integral is to use Fatou s Lemma to get E t = E lim k + mint, τ k lim inf E = /. k + mint, τ k So one gets only inequality instead of the equality E /t = /. Hence interchanging epectation and integration is not implied, as only Fatou s lemma can be 678
5 Proceedings of the 9th International Symposium on Mathematical Theory of Networks and Systems MTNS 2 9 July, 2 udapest, Hungary used to provide the bound E. 2 t Also we note that does not satisfy a Lipshitz condition since y cannot be bounded by C y for any constant C. Hence this process may not be an Itô diffusion. We note that the three dimensional rownian motion is trivially an Itô diffusion. Applying the Itô-rule to the map φ : w = w w 2 w w = w 2 + w2 2 + w2 is dubious, since φ is not C 2 at the origin. Formal application gives i= d w = w + i= w i dw i w. w Now, idw i w has the same finite-dimensional distributions as a one dimensional rownian motion dw [, p. ]. Since [w, w 2, w ] never hits the origin, is an Itô diffusion. It can be shown that the process uw = w is a positive supermartingale [, p.]. If {T n} is an increasing sequence of stopping times, with T n, then the process uwt T n is a martingale, implying that uwt is a local martingale. Since E uwt as t, and E uw =, it follows that uwt cannot be a martingale. VIII. CONCLUSIONS We have illustrated by way of eample, that it is worthwhile to check the intricate details and conditions involved to guarantee that the differential rule and interchanging integrals integration and epectation are applicable. We have shown that carelessness leads to the paradoical situation of resulting in a process with a negative variance. The parado is resolved by developing properties of the three dimensional rownian motion, and its related essel process. For this approimation of the reflected rownian motion by a repulsive potential, eit probabilities and epected eit times from an interval were calculated. Sufficient conditions were presented under which the result of the Itô integral is a true martingale. Acknowledgement The author is indebted to Prof. Andrzej Swiech Math, GT, for valuable discussions on the parado epressed in section V and providing the impetus to unravel the parado. In addition, the reviewers comments, which helped to improve this note, are gratefully acknowledged. REFERENCES [] Z. Schuss, Theory and Application of Stochastic Differential Equations Wiley 98. [2] K. Itô and H.P. McKean, Jr. Diffusion Processes and their Sample Paths. Academic Press, 96. [] A. Papoulis, Probability, Random Variables and Stochastic Processes. McGraw-Hill, 96 []. /Oksendal, Stochastic Differential Equations, Fourth Ed. Springer- Verlag 99. [] P. Protter, Stochatic Integration and Differential Equations, Springer- Verlag 99. We recall the following result []. Let w be a Wiener process rownian motion. Define F t to be the σ-algebra generated by the random variables ws, for s t. Let ft, ω : [, Ω R be such that i t, ω ft, ω is F-measurable, where denotes the σ-algebra on [, ii ft, ω is F t -adapted iii E t s fτ, ω2 dτ < then the integral t fτ, ωdwt, ω in the Itô-sense is a s martingale. In fact the conditions can be relaed [, p.]: Condition ii may be replaced by ii : There eists an increasing family of σ-algebras H t, such that wt is a martingale with respect to H t and ft is H t -adapted. Condition iii can be relaed to: iii { t } Pr fs, ω 2 ds < for all t =. 679
OPTIMAL STOPPING OF A BROWNIAN BRIDGE
OPTIMAL STOPPING OF A BROWNIAN BRIDGE ERIK EKSTRÖM AND HENRIK WANNTORP Abstract. We study several optimal stopping problems in which the gains process is a Brownian bridge or a functional of a Brownian
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 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 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 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 informationON THE PATHWISE UNIQUENESS OF SOLUTIONS OF STOCHASTIC DIFFERENTIAL EQUATIONS
PORTUGALIAE MATHEMATICA Vol. 55 Fasc. 4 1998 ON THE PATHWISE UNIQUENESS OF SOLUTIONS OF STOCHASTIC DIFFERENTIAL EQUATIONS C. Sonoc Abstract: A sufficient condition for uniqueness of solutions of ordinary
More informationUpper Bounds for the Ruin Probability in a Risk Model With Interest WhosePremiumsDependonthe Backward Recurrence Time Process
Λ4flΛ4» ν ff ff χ Vol.4, No.4 211 8fl ADVANCES IN MATHEMATICS Aug., 211 Upper Bounds for the Ruin Probability in a Risk Model With Interest WhosePremiumsDependonthe Backward Recurrence Time Process HE
More informationUNCERTAINTY FUNCTIONAL DIFFERENTIAL EQUATIONS FOR FINANCE
Surveys in Mathematics and its Applications ISSN 1842-6298 (electronic), 1843-7265 (print) Volume 5 (2010), 275 284 UNCERTAINTY FUNCTIONAL DIFFERENTIAL EQUATIONS FOR FINANCE Iuliana Carmen Bărbăcioru Abstract.
More informationSome Terminology and Concepts that We will Use, But Not Emphasize (Section 6.2)
Some Terminology and Concepts that We will Use, But Not Emphasize (Section 6.2) Statistical analysis is based on probability theory. The fundamental object in probability theory is a probability space,
More informationOptimal Sojourn Time Control within an Interval 1
Optimal Sojourn Time Control within an Interval Jianghai Hu and Shankar Sastry Department of Electrical Engineering and Computer Sciences University of California at Berkeley Berkeley, CA 97-77 {jianghai,sastry}@eecs.berkeley.edu
More informationMATH 116, LECTURES 10 & 11: Limits
MATH 6, LECTURES 0 & : Limits Limits In application, we often deal with quantities which are close to other quantities but which cannot be defined eactly. Consider the problem of how a car s speedometer
More informationStochastic integral. Introduction. Ito integral. References. Appendices Stochastic Calculus I. Geneviève Gauthier.
Ito 8-646-8 Calculus I Geneviève Gauthier HEC Montréal Riemann Ito The Ito The theories of stochastic and stochastic di erential equations have initially been developed by Kiyosi Ito around 194 (one of
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 informationAPPLICATION OF THE KALMAN-BUCY FILTER IN THE STOCHASTIC DIFFERENTIAL EQUATIONS FOR THE MODELING OF RL CIRCUIT
Int. J. Nonlinear Anal. Appl. 2 (211) No.1, 35 41 ISSN: 28-6822 (electronic) http://www.ijnaa.com APPICATION OF THE KAMAN-BUCY FITER IN THE STOCHASTIC DIFFERENTIA EQUATIONS FOR THE MODEING OF R CIRCUIT
More informationStochastic Processes
Stochastic Processes A very simple introduction Péter Medvegyev 2009, January Medvegyev (CEU) Stochastic Processes 2009, January 1 / 54 Summary from measure theory De nition (X, A) is a measurable space
More informationSTAT 331. Martingale Central Limit Theorem and Related Results
STAT 331 Martingale Central Limit Theorem and Related Results In this unit we discuss a version of the martingale central limit theorem, which states that under certain conditions, a sum of orthogonal
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 information3. WIENER PROCESS. STOCHASTIC ITÔ INTEGRAL
3. WIENER PROCESS. STOCHASTIC ITÔ INTEGRAL 3.1. Wiener process. Main properties. A Wiener process notation W W t t is named in the honor of Prof. Norbert Wiener; other name is the Brownian motion notation
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 informationStudents were asked what is the slope of the line tangent to the graph of f at. and justify the answer.
18 AB5 Marshall Ransom, Georgia Southern University Problem Overview: Let f be the function defined by f ( ) e cos( ). Part a: Students were asked to find the average rate of change of f on the interval.
More informationONLINE APPENDIX TO: NONPARAMETRIC IDENTIFICATION OF THE MIXED HAZARD MODEL USING MARTINGALE-BASED MOMENTS
ONLINE APPENDIX TO: NONPARAMETRIC IDENTIFICATION OF THE MIXED HAZARD MODEL USING MARTINGALE-BASED MOMENTS JOHANNES RUF AND JAMES LEWIS WOLTER Appendix B. The Proofs of Theorem. and Proposition.3 The proof
More informationPart Two. Diagnostic Test
Part Two Diagnostic Test AP Calculus AB and BC Diagnostic Tests Take a moment to gauge your readiness for the AP Calculus eam by taking either the AB diagnostic test or the BC diagnostic test, depending
More informationThe Uniform Integrability of Martingales. On a Question by Alexander Cherny
The Uniform Integrability of Martingales. On a Question by Alexander Cherny Johannes Ruf Department of Mathematics University College London May 1, 2015 Abstract Let X be a progressively measurable, almost
More informationLimits and the derivative function. Limits and the derivative function
The Velocity Problem A particle is moving in a straight line. t is the time that has passed from the start of motion (which corresponds to t = 0) s(t) is the distance from the particle to the initial position
More informationLecture 17 Brownian motion as a Markov process
Lecture 17: Brownian motion as a Markov process 1 of 14 Course: Theory of Probability II Term: Spring 2015 Instructor: Gordan Zitkovic Lecture 17 Brownian motion as a Markov process Brownian motion is
More informationRichard F. Bass Krzysztof Burdzy University of Washington
ON DOMAIN MONOTONICITY OF THE NEUMANN HEAT KERNEL Richard F. Bass Krzysztof Burdzy University of Washington Abstract. Some examples are given of convex domains for which domain monotonicity of the Neumann
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 informationDaily Lessons and Assessments for AP* Calculus AB, A Complete Course Page 584 Mark Sparks 2012
The Second Fundamental Theorem of Calculus Functions Defined by Integrals Given the functions, f(t), below, use F( ) f ( t) dt to find F() and F () in terms of.. f(t) = 4t t. f(t) = cos t Given the functions,
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 informationEULER MARUYAMA APPROXIMATION FOR SDES WITH JUMPS AND NON-LIPSCHITZ COEFFICIENTS
Qiao, H. Osaka J. Math. 51 (14), 47 66 EULER MARUYAMA APPROXIMATION FOR SDES WITH JUMPS AND NON-LIPSCHITZ COEFFICIENTS HUIJIE QIAO (Received May 6, 11, revised May 1, 1) Abstract In this paper we show
More informationTaylor Series and Asymptotic Expansions
Taylor Series and Asymptotic Epansions The importance of power series as a convenient representation, as an approimation tool, as a tool for solving differential equations and so on, is pretty obvious.
More informationOptimal stopping for non-linear expectations Part I
Stochastic Processes and their Applications 121 (2011) 185 211 www.elsevier.com/locate/spa Optimal stopping for non-linear expectations Part I Erhan Bayraktar, Song Yao Department of Mathematics, University
More informationn E(X t T n = lim X s Tn = X s
Stochastic Calculus Example sheet - Lent 15 Michael Tehranchi Problem 1. Let X be a local martingale. Prove that X is a uniformly integrable martingale if and only X is of class D. Solution 1. If If direction:
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 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 Made Easy
Stochastic Calculus Made Easy Most of us know how standard Calculus works. We know how to differentiate, how to integrate etc. But stochastic calculus is a totally different beast to tackle; we are trying
More informationebay/google short course: Problem set 2
18 Jan 013 ebay/google short course: Problem set 1. (the Echange Parado) You are playing the following game against an opponent, with a referee also taking part. The referee has two envelopes (numbered
More informationTaylor Series and Series Convergence (Online)
7in 0in Felder c02_online.te V3 - February 9, 205 9:5 A.M. Page CHAPTER 2 Taylor Series and Series Convergence (Online) 2.8 Asymptotic Epansions In introductory calculus classes the statement this series
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 informationOptimal Trade Execution with Instantaneous Price Impact and Stochastic Resilience
Optimal Trade Execution with Instantaneous Price Impact and Stochastic Resilience Ulrich Horst 1 Humboldt-Universität zu Berlin Department of Mathematics and School of Business and Economics Vienna, Nov.
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 informationThe Pedestrian s Guide to Local Time
The Pedestrian s Guide to Local Time Tomas Björk, Department of Finance, Stockholm School of Economics, Box 651, SE-113 83 Stockholm, SWEDEN tomas.bjork@hhs.se November 19, 213 Preliminary version Comments
More informationIn economics, the amount of a good x demanded is a function of the price of that good. In other words,
I. UNIVARIATE CALCULUS Given two sets X and Y, a function is a rule that associates each member of X with eactly one member of Y. That is, some goes in, and some y comes out. These notations are used to
More informationSTAT Sample Problem: General Asymptotic Results
STAT331 1-Sample Problem: General Asymptotic Results In this unit we will consider the 1-sample problem and prove the consistency and asymptotic normality of the Nelson-Aalen estimator of the cumulative
More informationRates of Convergence to Self-Similar Solutions of Burgers Equation
Rates of Convergence to Self-Similar Solutions of Burgers Equation by Joel Miller Andrew Bernoff, Advisor Advisor: Committee Member: May 2 Department of Mathematics Abstract Rates of Convergence to Self-Similar
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 informationSolving the Poisson Disorder Problem
Advances in Finance and Stochastics: Essays in Honour of Dieter Sondermann, Springer-Verlag, 22, (295-32) Research Report No. 49, 2, Dept. Theoret. Statist. Aarhus Solving the Poisson Disorder Problem
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 informationA Change of Variable Formula with Local Time-Space for Bounded Variation Lévy Processes with Application to Solving the American Put Option Problem 1
Chapter 3 A Change of Variable Formula with Local Time-Space for Bounded Variation Lévy Processes with Application to Solving the American Put Option Problem 1 Abstract We establish a change of variable
More informationLecture 3: September 10
CS294 Markov Chain Monte Carlo: Foundations & Applications Fall 2009 Lecture 3: September 10 Lecturer: Prof. Alistair Sinclair Scribes: Andrew H. Chan, Piyush Srivastava Disclaimer: These notes have not
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 informationMATH 250 TOPIC 11 LIMITS. A. Basic Idea of a Limit and Limit Laws. Answers to Exercises and Problems
Math 5 T-Limits Page MATH 5 TOPIC LIMITS A. Basic Idea of a Limit and Limit Laws B. Limits of the form,, C. Limits as or as D. Summary for Evaluating Limits Answers to Eercises and Problems Math 5 T-Limits
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 informationSTABILITY. Have looked at modeling dynamic systems using differential equations. and used the Laplace transform to help find step and impulse
SIGNALS AND SYSTEMS: PAPER 3C1 HANDOUT 4. Dr David Corrigan 1. Electronic and Electrical Engineering Dept. corrigad@tcd.ie www.sigmedia.tv STABILITY Have looked at modeling dynamic systems using differential
More informationHW #1 SOLUTIONS. g(x) sin 1 x
HW #1 SOLUTIONS 1. Let f : [a, b] R and let v : I R, where I is an interval containing f([a, b]). Suppose that f is continuous at [a, b]. Suppose that lim s f() v(s) = v(f()). Using the ɛ δ definition
More informationChapter 31. The Laplace Transform The Laplace Transform. The Laplace transform of the function f(t) is defined. e st f(t) dt, L[f(t)] =
Chapter 3 The Laplace Transform 3. The Laplace Transform The Laplace transform of the function f(t) is defined L[f(t)] = e st f(t) dt, for all values of s for which the integral exists. The Laplace transform
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 submartingale / supermartingale property of diffusions in natural scale
On the submartingale / supermartingale property of diffusions in natural scale Alexander Gushchin Mikhail Urusov Mihail Zervos November 13, 214 Abstract Kotani 5 has characterised the martingale property
More informationMFE6516 Stochastic Calculus for Finance
MFE6516 Stochastic Calculus for Finance William C. H. Leon Nanyang Business School December 11, 2017 1 / 51 William C. H. Leon MFE6516 Stochastic Calculus for Finance 1 Symmetric Random Walks Scaled Symmetric
More informationMan Kyu Im*, Un Cig Ji **, and Jae Hee Kim ***
JOURNAL OF THE CHUNGCHEONG MATHEMATICAL SOCIETY Volume 19, No. 4, December 26 GIRSANOV THEOREM FOR GAUSSIAN PROCESS WITH INDEPENDENT INCREMENTS Man Kyu Im*, Un Cig Ji **, and Jae Hee Kim *** Abstract.
More informationNonlinear Control Systems
Nonlinear Control Systems António Pedro Aguiar pedro@isr.ist.utl.pt 3. Fundamental properties IST-DEEC PhD Course http://users.isr.ist.utl.pt/%7epedro/ncs2012/ 2012 1 Example Consider the system ẋ = f
More informationANALYTIC SEMIGROUPS AND APPLICATIONS. 1. Introduction
ANALYTIC SEMIGROUPS AND APPLICATIONS KELLER VANDEBOGERT. Introduction Consider a Banach space X and let f : D X and u : G X, where D and G are real intervals. A is a bounded or unbounded linear operator
More informationAP Calculus AB/IB Math SL2 Unit 1: Limits and Continuity. Name:
AP Calculus AB/IB Math SL Unit : Limits and Continuity Name: Block: Date:. A bungee jumper dives from a tower at time t = 0. Her height h (in feet) at time t (in seconds) is given by the graph below. In
More informationMath 1314 Lesson 4 Limits
Math 1314 Lesson 4 Limits What is calculus? Calculus is the study of change, particularly, how things change over time. It gives us a framework for measuring change using some fairly simple models. In
More informationThe Cameron-Martin-Girsanov (CMG) Theorem
The Cameron-Martin-Girsanov (CMG) Theorem There are many versions of the CMG Theorem. In some sense, there are many CMG Theorems. The first version appeared in ] in 944. Here we present a standard version,
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 informationVerona Course April Lecture 1. Review of probability
Verona Course April 215. Lecture 1. Review of probability Viorel Barbu Al.I. Cuza University of Iaşi and the Romanian Academy A probability space is a triple (Ω, F, P) where Ω is an abstract set, F is
More information3.3.1 Linear functions yet again and dot product In 2D, a homogenous linear scalar function takes the general form:
3.3 Gradient Vector and Jacobian Matri 3 3.3 Gradient Vector and Jacobian Matri Overview: Differentiable functions have a local linear approimation. Near a given point, local changes are determined by
More informationThe Skorokhod problem in a time-dependent interval
The Skorokhod problem in a time-dependent interval Krzysztof Burdzy, Weining Kang and Kavita Ramanan University of Washington and Carnegie Mellon University Abstract: We consider the Skorokhod problem
More informationA NOTE ON STOCHASTIC INTEGRALS AS L 2 -CURVES
A NOTE ON STOCHASTIC INTEGRALS AS L 2 -CURVES STEFAN TAPPE Abstract. In a work of van Gaans (25a) stochastic integrals are regarded as L 2 -curves. In Filipović and Tappe (28) we have shown the connection
More informationMathematical Tripos Part III Michaelmas 2017 Distribution Theory & Applications, Example sheet 1 (answers) Dr A.C.L. Ashton
Mathematical Tripos Part III Michaelmas 7 Distribution Theory & Applications, Eample sheet (answers) Dr A.C.L. Ashton Comments and corrections to acla@damtp.cam.ac.uk.. Construct a non-zero element of
More informationSolution of Partial Differential Equations
Solution of Partial Differential Equations Introduction and Problem Statement Combination of Variables R. Shankar Subramanian We encounter partial differential equations routinely in transport phenomena.
More informationUBC-SFU-UVic-UNBC Calculus Exam Solutions 7 June 2007
This eamination has 15 pages including this cover. UBC-SFU-UVic-UNBC Calculus Eam Solutions 7 June 007 Name: School: Signature: Candidate Number: Rules and Instructions 1. Show all your work! Full marks
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 informationMaximum Process Problems in Optimal Control Theory
J. Appl. Math. Stochastic Anal. Vol. 25, No., 25, (77-88) Research Report No. 423, 2, Dept. Theoret. Statist. Aarhus (2 pp) Maximum Process Problems in Optimal Control Theory GORAN PESKIR 3 Given a standard
More informationSystem Identification & Parameter Estimation
System Identification & Parameter Estimation Wb3: SIPE lecture Correlation functions in time & frequency domain Alfred C. Schouten, Dept. of Biomechanical Engineering (BMechE), Fac. 3mE // Delft University
More informationOPTIMAL SOLUTIONS TO STOCHASTIC DIFFERENTIAL INCLUSIONS
APPLICATIONES MATHEMATICAE 29,4 (22), pp. 387 398 Mariusz Michta (Zielona Góra) OPTIMAL SOLUTIONS TO STOCHASTIC DIFFERENTIAL INCLUSIONS Abstract. A martingale problem approach is used first to analyze
More informationFundamental Inequalities, Convergence and the Optional Stopping Theorem for Continuous-Time Martingales
Fundamental Inequalities, Convergence and the Optional Stopping Theorem for Continuous-Time Martingales Prakash Balachandran Department of Mathematics Duke University April 2, 2008 1 Review of Discrete-Time
More informationLaplace Transform Part 1: Introduction (I&N Chap 13)
Laplace Transform Part 1: Introduction (I&N Chap 13) Definition of the L.T. L.T. of Singularity Functions L.T. Pairs Properties of the L.T. Inverse L.T. Convolution IVT(initial value theorem) & FVT (final
More informationSection 14.1 Vector Functions and Space Curves
Section 14.1 Vector Functions and Space Curves Functions whose range does not consists of numbers A bulk of elementary mathematics involves the study of functions - rules that assign to a given input a
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 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 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 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 informationIt is convenient to introduce some notation for this type of problems. I will write this as. max u (x 1 ; x 2 ) subj. to. p 1 x 1 + p 2 x 2 m ;
4 Calculus Review 4.1 The Utility Maimization Problem As a motivating eample, consider the problem facing a consumer that needs to allocate a given budget over two commodities sold at (linear) prices p
More informationAnswer Key 1973 BC 1969 BC 24. A 14. A 24. C 25. A 26. C 27. C 28. D 29. C 30. D 31. C 13. C 12. D 12. E 3. A 32. B 27. E 34. C 14. D 25. B 26.
Answer Key 969 BC 97 BC. C. E. B. D 5. E 6. B 7. D 8. C 9. D. A. B. E. C. D 5. B 6. B 7. B 8. E 9. C. A. B. E. D. C 5. A 6. C 7. C 8. D 9. C. D. C. B. A. D 5. A 6. B 7. D 8. A 9. D. E. D. B. E. E 5. E.
More informationNon-Zero-Sum Stochastic Differential Games of Controls and St
Non-Zero-Sum Stochastic Differential Games of Controls and Stoppings October 1, 2009 Based on two preprints: to a Non-Zero-Sum Stochastic Differential Game of Controls and Stoppings I. Karatzas, Q. Li,
More informationHarmonic Analysis Homework 5
Harmonic Analysis Homework 5 Bruno Poggi Department of Mathematics, University of Minnesota November 4, 6 Notation Throughout, B, r is the ball of radius r with center in the understood metric space usually
More informationE-Companion to The Evolution of Beliefs over Signed Social Networks
OPERATIONS RESEARCH INFORMS E-Companion to The Evolution of Beliefs over Signed Social Networks Guodong Shi Research School of Engineering, CECS, The Australian National University, Canberra ACT 000, Australia
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 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 informationTopics in fractional Brownian motion
Topics in fractional Brownian motion Esko Valkeila Spring School, Jena 25.3. 2011 We plan to discuss the following items during these lectures: Fractional Brownian motion and its properties. Topics in
More informationSOLVING DIFFERENTIAL EQUATIONS
SOLVING DIFFERENTIAL EQUATIONS ARICK SHAO 1. Solving Ordinary Differential Equations Consider first basic algebraic equations with an unknown, e.g., 3(x + 2) = 18, 3x 2 = cos x. Such equations model some
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 informationAn Introduction to Malliavin calculus and its applications
An Introduction to Malliavin calculus and its applications Lecture 3: Clark-Ocone formula David Nualart Department of Mathematics Kansas University University of Wyoming Summer School 214 David Nualart
More informationIntroduction. So, why did I even bother to write this?
Introduction This review was originally written for my Calculus I class, but it should be accessible to anyone needing a review in some basic algebra and trig topics. The review contains the occasional
More informationPre-Calculus Mathematics Limit Process Calculus
NOTES : LIMITS AND DERIVATIVES Name: Date: Period: Mrs. Nguyen s Initial: LESSON.1 THE TANGENT AND VELOCITY PROBLEMS Pre-Calculus Mathematics Limit Process Calculus The type of it that is used to find
More informationChange detection problems in branching processes
Change detection problems in branching processes Outline of Ph.D. thesis by Tamás T. Szabó Thesis advisor: Professor Gyula Pap Doctoral School of Mathematics and Computer Science Bolyai Institute, University
More informationM151B Practice Problems for Exam 1
M151B Practice Problems for Eam 1 Calculators will not be allowed on the eam. Unjustified answers will not receive credit. 1. Compute each of the following its: 1a. 1b. 1c. 1d. 1e. 1 3 4. 3. sin 7 0. +
More informationLOCATION OF THE PATH SUPREMUM FOR SELF-SIMILAR PROCESSES WITH STATIONARY INCREMENTS. Yi Shen
LOCATION OF THE PATH SUPREMUM FOR SELF-SIMILAR PROCESSES WITH STATIONARY INCREMENTS Yi Shen Department of Statistics and Actuarial Science, University of Waterloo. Waterloo, ON N2L 3G1, Canada. Abstract.
More information