Journal of Applied Mathematics

Volume 2014 (2014), Article ID 469308, 10 pages

http://dx.doi.org/10.1155/2014/469308

## Modified Block Pulse Functions for Numerical Solution of Stochastic Volterra Integral Equations

Department of Mathematics, Islamic Azad University, Karaj Branch, Karaj, Iran

Received 13 December 2013; Accepted 22 January 2014; Published 13 March 2014

Academic Editor: Kai Diethelm

Copyright © 2014 K. Maleknejad et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

We present a new technique for solving numerically stochastic Volterra integral equation based on modified block pulse functions. It declares that the rate of convergence of the presented method is faster than the method based on block pulse functions. Efficiency of this method and good degree of accuracy are confirmed by a numerical example.

#### 1. Introduction

The numerical study and simulation of stochastic Volterra integral equations (SVIEs) have been an active field of research for the past years [1–7]. Most SVIEs do not have analytic solutions and hence it is of great importance to provide numerical schemes. Numerical schemes to stochastic differential equations (SDEs) have been well developed [8–12]. However, there are still few papers discussing the numerical solutions for stochastic Volterra integral equations.

Study in economics, sociology, and various biological and medical models leads to the stochastic Volterra integral equations. These systems are dependent on a noise source, on a Gaussian white noise, so that modeling such phenomena naturally requires the use of various stochastic Volterra integral equations.

In this paper, we consider the linear stochastic Volterra integral equation: where , , , and , for , are the stochastic processes defined on the same probability space with a filtration that is increasing and right continuous and contains all -null sets. is unknown random function and is a standard Brownian motion defined on the probability space and is the Itô integral. Numerous papers have been focusing on existence solution of (1) [13–15].

The paper [3] solves stochastic Volterra integral equations by block pulse functions (BPFs) and [4] applies this method for solving -dimensional stochastic Itô Volterra integral equations. However, BPFs are very common in use; it seems that their convergence is weak. Maleknejad and Rahimi apply in [16] Modified Block Pulse Functions (MBPFs) to solve Volterra integral equation of the first kind numerically. Here, we use this method for solving SVIEs.

This paper is organized as follows. In the rest of this section we describe some general concepts concerning the block pulse functions and epsilon modified block pulse functions and some concepts related to stochastic and Itô integral. Section 2 is devoted to stochastic integration operational matrix. In Section 3, the method is employed to solve stochastic integral equations. Section 4 discusses error analysis of this method. Section 5 gives numerical example. Finally, Section 6 provides the conclusion of this work.

##### 1.1. Block Pulse Functions

BPFs have been variously studied [16–18] and applied for solving different problems. The goal of this section is to recall notations and definition of the BPFs that are used in the next sections.

The block pulse functions are defined on the time interval by where and for convenience we put .

The block pulse functions on have the following properties:(1)disjointness: for where is Kronecker delta;(2)orthogonality: (3)completeness: for every when approaches infinity, Parseval’s identity holds: where Also the Fourier coefficients and the block pulse functions depend on . The set of block pulse functions may be written as a vector of dimension : From the above representation and disjointness property, it follows that where is an -dimensional vector and . Let be an matrix so that where is a vector with elements equal to the diagonal entries of .

The expansion of a function over with respect to , , is given by where and is defined by (6).

Let . It is expanded with respect to BPFs as where and are - and -dimensional BPFs vectors, respectively, and is the block pulse coefficient matrix with the below : For convenience, we put .

Now, integration operational matrix is considered and computed: Since is equal to at midpoint of , we can approximate , for by . Therefore where the th component is . As a result where is operational matrix of integration that is given by So,

##### 1.2. Epsilon Modified Block Pulse Functions (EMBPFs)

A set of epsilon modified block pulse functions , , on the interval are defined as for , and

Similar to BPFs, the most important properties of EMBPFs are(1)disjointness: where ;(2)orthogonality: if we put , (3)completeness: where and is length of interval .

With defining , we have Similar to BPFs, where the operational matrix of EMBPFs is given by and we have the following approximation:

##### 1.3. Stochastic Concepts of Itô Integral

*Definition 1 (Brownian motion process). *A real-valued stochastic process is called Brownian motion, if it satisfies the following properties:(i)independence of increments: , , is independent of the past, that is, of , or of , the -field generated by ;(ii)normal increments: has normal distribution with mean 0 and variance ;(iii)continuity of paths: , are continuous functions of .

*Definition 2. *Let be an increasing family of -algebras of subsets of . A process from to is called -adapted if for each the function is -measurable [19].

*Definition 3 (see [19]). *Let be the class of functions such that(i) is -measurable, where denotes the Borel -algebra on and is the -algebra on ;(ii) is -adapted, where is the -algebra generated by the random variables , ;(iii).

*Definition 4 (the Itô integral, [19]). *Let ; then the Itô integral of (from to ) is defined by
where is a sequence of elementary functions such that

Theorem 5 (the Itô isometry). *Let ; then
*

*Proof. *See [19].

*Definition 6 (1-dimensional Itô processes, [19]). *Let be 1-dimensional Brownian motion on . A 1-dimensional Itô process (stochastic integral) is a stochastic process on of the form
or
where

*Theorem 7 (the 1-dimensional Itô formula). Let be an Itô process given by (1) and ; then
is again an Itô process, and
where is computed according to the rules
*

*Proof. *See [19].

Moreover, is notation of

*Lemma 8 (the Gronwall inequality). Let be integral with
for , where . Then
*

*For more details see [19, 20].*

*2. Stochastic Integral Operational Matrix for EMBPFs*

*2. Stochastic Integral Operational Matrix for EMBPFs**In this section stochastic integral operational matrix for EMBPFs is considered. For finding vector form of , with EMBPFs, the Itô integral of each single EMBPF can be computed as follows. It is clear that the integrals are stochastic and nondeterministic:
for , and
We approximate(1), by , at midpoint of ;(2) by in at midpoint of ;(3) by in , at midpoint of .*

*
As a result, vector form of , with EMBPFs, is given by
in which the th component is ,
Therefore
where stochastic operational matrix of integration is given bySo, the Itô integral of every function can be approximated as follows:
*

*3. Numerical Solution of SVIEs by EMBPFs*

*3. Numerical Solution of SVIEs by EMBPFs**Here, we modify the method that has been used in [16] by EMBPFs. In the below equation:
we approximate functions , , , and by EMBPFs:
where the vectors , and matrices , are EMBPFs coefficient of , , , and , respectively.*

*Substituting (48) into (47) and using previous relations,
Finally
whereThen
With replacing by , we have a linear system of equations.*

*Now if , there will be numerical answers . Solution is approximated by
*

*4. Error Analysis*

*4. Error Analysis**In this section, error analysis is studied. In the following theorems, for simplicity, we assume and .*

*Theorem 9. If and , then(1) achieves its minimum value;(2) approach pointwise;(3).*

*Proof. *See [16].

*Theorem 10. Assume the following.(1) is continuous and differentiable in , with bounded derivative; that is, .(2), are correspondingly BPFs. MBPFs, MBPFs expansions of base on EMBPFs over interval .(3).*

Then

*Proof. *Trapezoidal rule for integral is
where is error of integration. Suppose and . The representation error when is represented by a series of BPFs over every subinterval , is
where
From (55),

It is obvious that if , then .

So, this error is computed for in interval , .

For this function , so
Then this error with BPFs is .

Similarly, the error when is represented in a series of EMBPFs over every subinterval is
So, the error with EMBPFs is .

For in we have
So, the error is also for .

Now,

We define the representation error between and its 2D-EMBPFs expansion, , over every subregion , is defined as
where

With Taylor’s expansion and similarity to the above discussion,

*Theorem 11. Assume that(1),(2). *

Then

*Proof. *Consider
So,
by the Cauchy-Schwartz inequality, Itô isometry formula, and the linearity of Itô integrals in their integrands.

The first term is satisfied by last theorem:
Now,
Furthermore,
Hence
Then by Gronwall’s inequality, we get

*5. Numerical Example*

*5. Numerical Example**In this section, we present an example for showing the features of the EMBPFs method in this paper. Let denote the EMBP coefficient of exact solution of the given example and let be the EMBP coefficient of computed solution by the presented method. In this example error is defined as
*

*Example 1 (see [3]). *Consider the following linear stochastic Volterra integral equation:

*
with the exact solution , for .*

*The numerical results are shown in Tables 1 and 2. In the tables, is the number of iterations, is error mean, and is standard deviation of error.*

*Table 3 is from [3] for comparison.*

*In some examples by applying BPFs when increases, accuracy decreases, but in EMBPFs we achieve good accuracy by increasing .*

*6. Conclusion*

*6. Conclusion**As some SVIEs cannot be solved analytically, in this paper we present a new technique for solving SVIEs numerically. Here, we consider a modification of the block pulse functions. Some theorems show that if EMBPFs are used for achieving numerical expansions with times more precision, there is no need to increase the number of BPFs, times, which leads to solving a system of equations with times more equations and unknowns. But the results of BPFs solution can be combined with solutions of systems of equations with one more unknown and nearly achieve times more precision. Parallel programming is so useful for this method. Efficiency of this method and good degree of accuracy are confirmed by a numerical example.*

*Conflict of Interests*

*Conflict of Interests**The authors declare that there is no conflict of interests regarding the publication of this paper.*

*References*

*References*

- J. A. D. Appley, S. Devin, and D. W. Reynolds, “Almost sure convergence of solutions of linear stochastic Volterra equations to nonequilibrium limit,”
*Journal of Integral Equations and Applications*, vol. 19, no. 4, pp. 405–437, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. A. Berger and V. J. Mizel, “Volterra equations with Itô integrals. I,”
*Journal of Integral Equations*, vol. 2, no. 3, pp. 187–245, 1980. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - K. Maleknejad, M. Khodabin, and M. Rostami, “Numerical solution of stochastic Volterra integral equations by a stochastic operational matrix based on block pulse functions,”
*Mathematical and Computer Modelling*, vol. 55, no. 3-4, pp. 791–800, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - K. Maleknejad, M. Khodabin, and M. Rostami, “A numerical method for solving
*m*-dimensional stochastic Itô-Volterra integral equations by stochastic operational matrix,”*Computers & Mathematics with Applications*, vol. 63, no. 1, pp. 133–143, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - E. Pardoux and P. Protter, “Stochastic Volterra equations with anticipating coeffcients,”
*The Annals of Probability*, vol. 18, no. 4, pp. 1635–1655, 1990. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Y. Shiota, “A linear stochastic integral equation containing the extended Itô integral,”
*Mathematics Reports, Toyama University*, vol. 9, pp. 43–65, 1986. View at Google Scholar - C. H. Wen and T. S. Zhang, “Improved rectangular method on stochastic Volterra equations,”
*Journal of Computational and Applied Mathematics*, vol. 235, no. 8, pp. 2492–2501, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - D. J. Higham, X. Mao, and A. M. Stuart, “Strong convergence of Euler-type methods for nonlinear stochastic differential equations,”
*SIAM Journal on Numerical Analysis*, vol. 40, no. 3, pp. 1041–1063, 2002. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - M. Khodabin, K. Maleknejad, M. Rostami, and M. Nouri, “Numerical solution of stochastic differential equations by second order Runge-Kutta methods,”
*Mathematical and Computer Modelling*, vol. 53, no. 9-10, pp. 1910–1920, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - P. E. Kloeden and E. Platen,
*Numerical Solution of Stochastic Differential Equations*, Applications of Mathematics, Springer, Berlin, Germany, 1992. View at MathSciNet - X. Mao,
*Stochastic Differential Equations and Applcation*, Horwood Series in Mathematics and Applications, Horwood Publishing Limited, Chichester, UK, 1997. View at MathSciNet - E. Platen, “An introduction to numerical methods for stochastic differential equations,”
*Acta Numerica*, vol. 8, pp. 197–246, 1999. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. Janković and D. Ilić, “One linear analytic approximation for stochastic integrodifferential equations,”
*Acta Mathematica Scientia*, vol. 30, no. 4, pp. 1073–1085, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - W. J. Padgett and C. P. Tsokos, “Existence of a solution of a stochastic integral equation in turbulence theory,”
*Journal of Mathematical Physics*, vol. 12, no. 2, pp. 210–212, 1971. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - R. Subramaniam, K. Balachandran, and J. K. Kim, “Existence of solutions of a stochastic integral equation with an application from the theory of epidemics,”
*Nonlinear Functional Analysis and Applications*, vol. 5, no. 1, pp. 23–29, 2000. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - K. Maleknejad and B. Rahimi, “Modification of block pulse functions and their application to solve numerically Volterra integral equation of the first kind,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 16, no. 6, pp. 2469–2477, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - E. Babolian, K. Maleknejad, M. Mordad, and B. Rahimi, “A numerical method for solving Fredholm-Volterra integral equations in two-dimensional spaces using block pulse functions and an operational matrix,”
*Journal of Computational and Applied Mathematics*, vol. 235, no. 14, pp. 3965–3971, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - K. Maleknejad, M. Shahrezaee, and H. Khatami, “Numerical solution of integral equations system of the second kind by Block-Pulse functions,”
*Applied Mathematics and Computation*, vol. 166, no. 1, pp. 15–24, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - B. K. Øksendal,
*Stochastic Differential Equations: An Introduction with Applications*, Springer, New York, NY, USA, 5th edition, 1998. View at MathSciNet - F. C. Klebaner,
*Intoduction to Stochastic Calculus with Applications*, Monash University, Melbourne, Australia, 2nd edition, 2005. View at MathSciNet

*
*