A stochastic differential equation, or SDE for short, is a differential equation driven by one or more stochastic processes. For example, in physics, a Langevin equation describing the motion of a point in n-dimensional phase space is of the form
The dynamics are described by the functions , and the problem is to find a solution for X, given its value at an initial time. What distinguishes this from an ordinary differential equation are random noise terms and, consequently, solutions to the Langevin equation are stochastic processes. It is difficult to say exactly how should be defined directly, but we can suppose that their integrals are continuous with independent and identically distributed increments. A candidate for such a process is standard Brownian motion and, up to constant scaling factor and drift term, it can be shown that this is the only possibility. However, Brownian motion is nowhere differentiable, so the original noise terms do not have well defined values. Instead, we can rewrite equation (1) is terms of the Brownian motions. This gives the following SDE for an n-dimensional process
where are independent Brownian motions. This is to be understood in terms of the differential notation for stochastic integration. It is known that if the functions are Lipschitz continuous then, given any starting value for X, equation (2) has a unique solution. In this post, I give a proof of this using the basic properties of stochastic integration as introduced over the past few posts.
First, in keeping with these notes, equation (2) can be generalized by replacing the Brownian motions and time t by arbitrary semimartingales. As always, we work with respect to a complete filtered probability space . In integral form, the general SDE for a cadlag adapted process is as follows,
Here, are semimartingales. For example, as well as Brownian motion, Lévy processes are often also used. If are -measurable random variables, they simply specify the starting value . More generally, in (3), we can allow N to be any cadlag and adapted process, which acts as a `source term’ in the SDE. Furthermore, rather than just being functions of X at time t, we allow to be a function of the process X at all times up to t. For example, it could depend on it maximum so far, or a running average. In that case, we instead impose a functional Lipschitz condition as follows. The notation is used for the maximum of a process X, and denotes the set of all cadlag and adapted n-dimensional processes.
- (P1) is a map from to the set of predictable and –integrable processes.
- (P2) There is a constant K such that
for all times and .
The special case where is a Lipschitz continuous function of clearly satisfies both (P1) and (P2). The uniqueness theorem for SDEs with Lipschitz continuous coefficients is as follows.
Theorem 1 Suppose that satisfy properties (P1), (P2). Then, there is a unique satisfying SDE (3).
Recall that, throughout these notes, we are identifying any processes which agree on a set of probability one so that, here, uniqueness refers to uniqueness up to evanescence.
After existence and uniqueness, an important property of solutions to SDEs is stability. That is, small changes to the coefficients only have a small effect on the solution. The notion of a `small change’ first needs to be made precise by specifying a topology on these terms. For the coefficients , uniform convergence on bounded sets will be used. Given maps for r=1,2,…, I will say that uniformly on bounded sets if, for each
tends to zero as r goes to infinity. For the solutions X and source term N, which are cadlag processes, the appropriate topology is that of convergence uniformly on compacts in probability (ucp convergence). Recall that if in probability for each t. The stability of solutions to SDEs under small changes to the coefficients can then be stated precisely.
Theorem 2 Suppose that satisfy properties (P1), (P2) and let X be the unique solution to the SDE (3) given by Theorem 1. Suppose furthermore that and are sequences such that and uniformly on bounded sets. Then, any sequence of processes satisfying
converges ucp to X as r goes to infinity.
Equivalently, the map from to the solution X satisfying SDE (3) is continuous under the respective topologies.
Proof of Existence and Uniqueness
Throughout this subsection, assume that N is a cadlag adapted n-dimensional process, are semimartingales and that satisfy properties (P1), (P2) above. To simplify the notation a bit, for any , write to denote the following n-dimensional process,
So, and equation (3) just says that X=F(X). Existence and uniqueness of solutions to this SDE is equivalent to F having a unique fixed point.
The `standard’ proof for Lipschitz continuous coefficients, at least when the driving semimartingales are Brownian motions, makes use of the Ito isometry to construct a norm on the cadlag processes under which F is a contraction. The result then follows from the contraction mapping theorem. Although this method can also be applied for arbitrary semimartingales, it is more difficult and some rather advanced results on semimartingale decompositions are required. However, the aim here is to show how existence and uniqueness follows for all semimartingales as a consequence of the basic properties of the stochastic integration, so a different approach is taken. Furthermore, the method used of constructing discrete approximations making the error term X-F(X) small is closer to methods which might be employed in practice when explicitly simulating solutions to (3).
First, the following inequality bounds the values of a positive semimartingale X in terms of the maximum value of the stochastic integral .
Proof: For any continuous decreasing function , is a cadlag decreasing process and integration by parts gives the following,
It is easily seen that is constant over any interval on which and, therefore, the integral above with respect to is zero over such intervals. So, we can replace the integrand by and, again applying integration by parts,
As the process is cadlag, there must exist an such that either or . In the first case, evaluating (4) at time s gives
The final equality follows from the fact that must be constant on the interval [s,t]. Similarly, in the case where , the same inequality results from evaluating (4) at time s-.
We are only interested in the special case . The next step is to apply the generalized Ito formula to .
The inequality follows from concavity of the logarithm and, putting shows that the summation above is non-positive,
The general idea behind the proofs given here is that if X,Y are processes such that the error terms X-F(X) and Y-F(Y) to the SDE (3) are small then
The Lipschitz property implies that the integrand is bounded by , and the following result will be used to show that .
for predictable processes and cadlag adapted processes converging ucp to zero as r goes to infinity. Then, .
Proof: By definition of ucp convergence, for any the limit holds as r goes to infinity. Therefore, there exists a sequence such that . So, by stopping the processes as soon as , we may assume that for all r.
Now, writing ,
for predictable processes , which are uniformly bounded by K. Applying Lemma 3 to the nonnegative semimartingale ,
However, applying integration by parts to ,
The inequality shows that the integrands on the right hand side of the above expression are all bounded by in the first integral and by in the second integral. Hence, by dominated convergence, if is a sequence of real numbers then
in probability as . This shows that the sequence is bounded in probability as and, by exponentiating, is bounded in probability. However, and therefore, . Finally, as required.
Now, with the aid of the previous result, it is possible to construct approximate solutions to the SDE (3) making the error X-F(X) as small as we like.
Proof: The idea is to define X inductively as a piecewise constant function across finite time intervals, while adjusting the time step sizes to force the error to remain within the tolerance . To do this, define stopping times and processes by and, for ,
Then, increases to some limit as r goes to infinity, and we can define the process X on the time interval by whenever . By definition, on (see Figure 1). To complete the proof, it just needs to be shown that almost surely.
First, we show that X cannot explode at any finite time. The process is bounded by , and choosing a sequence of real numbers , the stopped processes can be written as,
The final term on the right hand side tends ucp to zero as r goes to infinity, and the Lipschitz property implies that is bounded by . Lemma 4 then gives . So, for any time , the sequence is bounded in probability and, therefore, is almost surely bounded on the interval whenever . In particular, is a locally bounded and hence is an X-integrable process. So, whenever , the limit
exists with probability one.
However, from the definition of X,
contradicting the convergence of . So almost surely, as required.
Next, if we construct a sequence of processes such that the error terms go to zero, then they necessarily converge to the solution to the SDE.
Proof: Setting and then
and as r and s go to infinity. Furthermore, by Lipschitz continuity, . So, Lemma 4 says that as r,s go to infinity.
By completeness under ucp convergence, for some . Passing to a subsequence if necessary, we may suppose that tends to X uniformly on bounded intervals. So, dominated convergence gives
Existence and uniqueness of solutions, as stated by Theorem 1, now follows from the previous lemmas.
Proof: First, for uniqueness, suppose that and are two such solutions. Forming the infinite sequence by alternating these, for r even and for r odd, Lemma 6 says that this is convergent, so X=Y.
Finally, the proof of Theorem 2 also follows from the results above.
Theorem 8 Suppose that satisfy properties (P1), (P2) and let X be the unique solution to the SDE (3) given by Theorem 1. Suppose furthermore that and are sequences such that and uniformly on bounded sets. Then, any sequence of processes satisfying
converges ucp to X as r goes to infinity.
Proof: For any , let be the first time at which ,
As is uniformly bounded by L over the interval , dominated convergence implies that the first term inside the final parenthesis converges ucp to zero as r goes to infinity. Also, by Lipschitz continuity, is bounded by . So, applying Lemma 4 to the above expression gives . It remains to show that the non-stopped processes also converge ucp to zero.
Fix any . Note that if and then .
By choosing L large, this can be made as small as we like. Finally,
By ucp convergence of the stopped processes the first term on the right hand side vanishes as r goes to infinity, and the second term can be made as small as we like by choosing L large. Therefore, in probability, as required.