Dynamical Aspects of Initial/Boundary Value Problems for Ordinary Differential Equations 2014View this Special Issue
Research Article | Open Access
Strong Convergence of the Split-Step -Method for Stochastic Age-Dependent Capital System with Random Jump Magnitudes
We develop a new split-step (SS) method for stochastic age-dependent capital system with random jump magnitudes. The main aim of this paper is to investigate the convergence of the SS method for a class of stochastic age-dependent capital system with random jump magnitudes. It is proved that the proposed method is convergent with strong order 1/2 under given conditions. Finally, an example is simulated to verify the results obtained from theory.
Stochastic partial differential equations are becoming increasingly used to model real-world phenomena in different fields, such as economics, biology, and physics. Recently, the study of the stochastic age-dependent (vintage) capital system has received a great deal of attention. For example, Wang studied stability of solutions for stochastic investment system . Zhang et al. studied the convergence and exponential stability of numerical solutions to the stochastic age-dependent capital system [2, 3].
In the stochastic age-dependent capital system, due to the effects of external environment for capital system, such as innovations in technique, introduction of new products, natural disasters, and changes in laws and government policies, the size of the capital systems increases or decreases drastically. So Poisson jumps with deterministic jump magnitude have been used in stochastic age-dependent population equations. For example, Li et al.  studied the Euler numerical method for stochastic age-dependent population equations with Poisson jumps. L. Wang and X. Wang  analysed the convergence of the semi-implicit Euler method for stochastic age-dependent population equations with Poisson jumps. Rathinasamy et al.  developed the numerical method for stochastic age-dependent population equations with Poisson jump and phase semi-Markovian switching. However, the random jump magnitude is now commonly seen in financial models [7–9]. In this paper, we will consider the following stochastic age-dependent capital system with random jump magnitudes as shown in : where denotes the stock of capital goods of age at time , , , and . is defined as total output produced in year ; also is the age of the capital; the investment in the new capital. The is the appreciation (when ) or depreciation (when ) of the production capacity, and represents the volatility of the capital stock. The value of is the actual jump and is the underlying random variables of the magnitudes, and often it is called “mark" of the jump.
Also is a standard Wiener process. is a scalar Poisson process with intensity . It is assumed that for some there is a constant such that ; that is, the th moment of the jump magnitude is bounded. The maximum physical lifetime of capital , the planning interval of calendar time , the depreciation rate of capital, and the capital density (the initial distribution of capital over age) are given. The denotes the accumulative rate of capital at the moment of , , and is the technical progress at the moment of . This makes that total output produced in year be defined as . In each sector all the firms have an identical neoclassical technology and produce output using labor and capital. The production function is neoclassical, where is the total sum of capital at time and is the labor force.
The integral version of (1) is given by the equation where for fixed .
Since the system (1) does not have closed form solutions, it is necessary to develop numerical methods for (1). Recently, Zhang and Rathinasamy  first derived the numerical solutions for stochastic age-dependent capital system with random jump magnitudes. However, their method belongs to the classic explicit Euler method and has a lower accuracy in  if we do not consider the appropriate step sizes.
Higham and Kloeden  first constructed the split-step backward Euler (SSBE) method for nonlinear stochastic differential equations with Poisson jumps. Tan and Wang  studied the convergence and stability of the SSBE method for linear stochastic delay integrodifferential equations. Ding et al.  developed the split-step method for solving the stochastic differential equations. Rathinasamy  investigated the split-step methods for stochastic age-dependent population equations with Markovian switching. Thus, we can construct the SS method for stochastic age-dependent population equations with random Poisson jumps.
In this paper, we will investigate the convergence of the SS method for system (1). The outline of the paper is as follows. In Section 2, we will introduce some preliminary results which are essential for our analysis. Section 3 will show us the SS method for solving stochastic age-dependent population equations with random Poisson jumps. In Section 4, several lemmas which are useful for our main result are proved. We give the main result that the numerical solutions converge to the true solutions with strong order 1/2 in Section 5. At last, a numerical example is given to verify the results obtained from the theory.
Throughout this paper, it will be denoted by the space of functions that are square-integrable over the domain . Let where is generalized partial derivative with respect to age and is a Sobolev space. such that . is the dual space of . We denote by , and the norms in , , and , respectively, by the scalar product in . , the duality product between and , is defined by
Let be a Wiener process defined on and taking its values in the separable Hilbert space , where is an orthonormal set of eigenvectors of , are mutually independent real Wiener processes with incremental covariance , , and (tr denotes the trace of an operator). For an operator to be the space of all bounded linear operators from into , it is denoted by the Hilbert-Schmidt norm; that is,
Let be the space of all continuous function from into with sup-norm =, , and .
Let be a complete probability space with a filtration , satisfying the usual conditions (i.e., it is increasing and right continuous while contains all -null sets).
Definition 1. Let be the stochastic basis and a Wiener process. Suppose that is a random variable such that . A stochastic process for fixed is said to be a solution on to the stochastic age-structured capital system for if the following conditions are satisfied:(D1) is a -measurable random variable;(D2), , , where denotes the space of all -valued processes (one will write for short) measurable (from into ), and satisfying . Here denotes the space of all continuous functions from to ;(D3)it satisfies the following equation: for all , , a.e. , where the stochastic integrals are understood in the Itô sense.
The parameter is the maximal age of the capital, so , .
As the standing hypotheses we always assume that the following conditions are satisfied.(1), and , .(2)(Lipschitz condition) there exists a positive constant such that and ,(i),(ii).(3) is nonnegative measurable in such that and is nonnegative continuous in such that (4), , , , where is a positive constant.
We note for the following jump process: where and is the indicator function for the set ; that is, .
Lemma 2. There exists a constant for any and for some such that
3. The Split-Step -Method
Let denote the th jump of occurrence time. Suppose, for example, that jumps arrive at distinct, ordered times ; let be a deterministic grid point of . We construct approximate solutions to models of the form (1) at a discrete set of times . This set is the superposition of the random jump times of a Poisson process on [0,] and a deterministic grid and satisfies max. Let , , and denote the increments of the time, Brownian motion, and the Poisson processes, respectively.
For system (1) the split-step approximate solution is defined by the iterative schemewith initial values , , ; is the numerical approximation of with ; the time increment is . Where , if we give , the SS method becomes the SSBE method in . If , the SS method is an explicit method.
To answer the question of the existence of numerical solution, we will give the following lemma.
Lemma 3. Let conditions and hold, and let and ; then the implicit equation (12a) can be solved uniquely for , with probability 1.
When Lemma 3 followed, we find it is convenient to use continuous-time approximation solution in our strong convergence analysis; hence for , we can define the following step functions: where is the largest number such that .
It is easy to verify that ; that is, and coincide with the discrete solutions at the grid points. Hence we refer to as a continuous-time extension of the discrete approximation . So our plan is to prove a strong convergence result for .
4. Several Lemmas
In this section, we will give several lemmas which are useful for the following main result.
Lemma 4. Under the conditions , there are constants and such that
The proof is similar to that in .
Next lemma shows the relationship between and .
Proof. Squaring both sides of (12a) and using the elementary inequality , we find Using the elementary inequality , we obtain Due to and conditions , we can get Taking mathematical expectation for both sides, we can obtain Since , thus ; then by , we have where and . The proof is completed.
The next lemma shows that the continuous numerical solutions have bounded second moments.
Lemma 6. Under conditions , let , ; then there exists a constant such that where , , and .
Proof. From (17), we can get Applying the It formula  to it yields Using conditions and the compensated Poisson process , we have Note where . Consider Taking (29)-(30) into (28), we compute that for some positive constants , , , Now, it follows that for any and by Lemma 5 By Burkholder-Davis-Gundy inequality, we have for some positive constants , , . Substituting (33) into (32) and then by Lemma 5, we can obtain that where Applying the continuous Gronwall inequality, we can easy get with .
The next lemma shows that the continuous-time approximation in (17) remains close to the step functions and in the mean square sense.
Proof. Consider ; we have Squaring both sides and using the element inequality , we have Now, the Cauchy-Schwarz inequality, condition (), and the compensated Poisson process give Taking mathematical expectation, by element inequality , and using the martingale isometry, we have By conditions () and (), we get Since and on , we have Then by Lemmas 5 and 6, we can derive