arxiv:505.04855v [math.na] 9 May 205 Approximation solution of two-dimensional linear stochastic Volterra integral equation by applying the Haar wavelet M. Fallahpour, M. Khodabin 2, K. Maleknejad 3 Department of Mathematics, Karaj Branch, Islamic Azad University, Karaj, Iran. Abstract. Numerical solution of one-dimensional stochastic integral equations because of the randomness has its own problems, i.e. some of them no have analytical solution or finding their analytic solution is very difficult. This problem for two-dimensional equations is twofold. Thus, finding an efficient way to approximate solutions of these equations is an essential requirement. To begin this important issue in this paper, we will give an efficient method based on Haar wavelet to approximate a solution for the two-dimensional linear stochastic Volterra integral equation. We also give an example to demonstrate the accuracy of the method. Mathematical subject classification: 65C30, 65C20, 60H20, 60H35, 68U20. Keywords: Haar wavelet; Two-dimensional stochastic Volterra integral equation; Brownian motion process; Ito integral.. Introduction As we know, two dimensional ordinary integral equations provide an important tool for modeling a numerous problems in engineering and science [6, 7]. The second kind of two-dimensional integral equations may arise from some problems of nonhomogeneous elasticity and electrostatics [8]. M.Fallahpour@kiau.ac.ir. 2 Corresponding author: M-Khodabin@kiau.ac.ir. 3 Maleknejad@iust.ac.ir.
2 Dobner presented an equivalent formulation of the Dorboux problem as a twodimensional Volterra integral equation [9]. We can also see this kind of equations in contact problems for bodies with complex properties [0, ], and in the theory of radio wave propagation [2], and in the theory of the elastic problem of axial translation of a rigid elliptical disc-inclusion [3], and various physical, mechanical and biological problems. Some numerical schemes have been inspected for resolvent of two-dimensional ordinary integral equations by several probers. Computational complexity of mathematical operations is the most important obstacle for solving ordinary integral equations in higher dimensionas. The Nystrom method [4], collocation method [5, 6, 7], Gauss product quadrature rule method [8], Galerkin method [9], using triangular fuctions [20, 2], Legender polynomial method [22], differential transform method [23], meshless method [24], Bernstein polynomials method [25] and Haar wavelet method [26]. This paper is first focused on proposing a generic framework for numerical solution of two-dimensional ordinary linear Volterra integral equations of second kind. The use of the Haar wavelet for the numerical solution of linear integral equations has previously been discussed in [] and references therein. The paper [] should be considered as a logical continuation of the papers [2 4]. In [2] a new numerical method based on Haar wavelet is introduced for solution of nonlinear one-dimensional Fredholm and Volterra integral equations. In [3] the Haar wavelet method [2] is extended to numerical solution of integro-differential equation. In [4] the Haar wavelet method [2, 3] is improved in terms of efficiency by introducing one-dimentional Haar wavelet approximation of the kernel function. The method [] is fundamentally different from the other numerical methods based on Haar wavelet for the numerical solution of integral equations as it approximates kernel function using Haar wavelet. The general hyperbolic differential equation is defined as [9]. u xy = rx, y, u, u x, u y, x, y B R 2 u ˆB = hx, y, ˆB B, where the domain B and the subset ˆB of the border B are chosen according to the different initial value problems. It s easy to show that the integral form of. is given by two-dimentional Volterra integral equation gx, y = fx, y y x 0 0 K x, y, s, tgs, tdsdt.
Similarly, if we import statistical noise in to., we can obtain two-dimensional linear stochastic Volterra integral equation of the second kind, i.e..2 gx, y = fx, y y x 0 0 y x 0 0 K x, y, s, tgs, tdsdt K 2 x, y, s, tgs, tdbsdbt x, y [0, ] [0, ], s x < t y. 3 where the kernels K x, y, s, t and K 2 x, y, s, t in, 2 are known functions and fx, y is also a known function whereas gx, y is unknown function and is called the solution of two-dimensional stochastic integral equation. The condition s x < t y is necessary for adaptability to the filtration {F t ; 0 t } where F t = σ{bs; 0 s }. Lemma.Put φt, s = Kx, y, s, tgs, t. Let φ be a function in L 2 [0, ] 2. Then there exists a sequence φ n of off-diagonal step functions such that b lim n a b a φt, s φ n t, s 2 dtds = 0. Definition. Let φ L 2 [0, ] 2. Then the double Wiener-Itô integral of φ is defined as b b a a b φt, sdbtdbs = lim n a b a φ n t, sdbtdbs in L 2 Ω. 2. Haar Wavelets A wavelet family ψ j, i y j N,i Z is an orthonormal subfamily of the Hilbert space L 2 R with the property that all function in the wavelet family are generated from a fixed function ψ called mother wavelet through dilations and translations.the wavelet family satisfies the following relation ψ j, i y = 2 j/2 ψ 2 j y i. For Haar wavelet family on the interval [0, we have: {, fory [0, h y = 0, otherwise,
4 fory [α, h i y =, fory [, 0, otherwise, i = 2, 3,..., α = n m n 0.5 n, =, = m m ; m = 2 l, l = 0,,..., n = 0,,..., m. The integer l indicats the level of the wavelet and n is the translation parameter. Any square integrable function fy defined on [0, can be expressed as follows: fy = a i h i y, where a i are real constants. For approximation aim we consider a maximum value L of the integer l, level of the Haar wavelet in the above definition. The integer L is then called maximum level of resolution. We also define integer M = 2 L. Hence for any square integrable function fy we have a finite sum of Haar wavelets as follows: fy The following notation is introduced []: 2. p i, y = y 0 a i h i y. h i udu, where by the definition of Haar wavelet equation 2. reduce to y α, for y [α, p i, y = y, for y [, 0, elsewhere. We can also the following stochastic notation introduce 2.2 q i, y = y 0 h i udbu, where equation 2.2 can be evaluated similarly by the definition of Haar wavelet and is given as follow:
5 By Bα, for y [α, q i, y = 2B Bα By, for y [, 0, elsewhere. 3. Numerical method In this section, proposed numerical method [] will be discussed for two-dimensional linear stochastic Volterra integral equation of the second kind. In the first subsection, we state some results for efficient evalution of two-dimensional Haar wavelet approximations. In the second subsection, we apply these results for finding numerical solutions equation.. For Haar wavelet approximation of a function fx, y of two real variables x and y, we assume that the domain 0 x, y is divided into a grid of size 2M 2N using the following collocation points 3. x m = m 0.5, m =, 2,..., 2M, 2M 3.2 y n = n 0.5, n =, 2,..., 2N. 2N 3. Two-dimensional Haar wavelet system A real-valued function Gx, y of two real variables x and y can be approximated using two-dimensional Haar wavelets basis as [, 27]: 3.3 Gx, y p= q= b p,q h p xh q y. In order to calculate the unknown coefficients b i,j s, the collocation points defined in Eqs. 3. and 3.2 are substituted in Eq. 3.3.Hence, we obtain the following 2M 2N linear system with unknowns b i,j s: 3.4 Gx m, y n = p= q= b p,q h p x m h q y n, m =, 2,..., 2M, n =, 2,..., 2N. The solution of system 3.4 can be calculated from the following theorem. Theorem 2. The solution of the system 3.4 is given below:
6 b i, = b,j = ρ 2N p=α 2M ρ 2 2M p= b, = q= 2M 2N Gx m, y n Gx m, y n p= q= p= q= p= q= 2 Gx m, y n, Gx m, y n Gx m, y n, i = 2, 3,..., 2M,, j = 2, 3,..., 2N, b i,j = where ρ ρ 2 p=α Gx m, y n p= q= 2 Gx m, y n p=α q= 2 Gx m, y n p= Gx m, y n, i = 2, 3,..., 2M, j = 2, 3,..., 2N, 3.5 and similarly, 3.6 Proof. See [2]. α = ρ σ, = ρ σ ρ 2, = ρ σ, ρ = 2M τ, σ = i τ, τ = 2 log 2 i α 2 = ρ 2 σ 2, 2 = ρ 2 σ 2 ρ 2 2, 2 = ρ 2 σ 2, ρ 2 = 2N τ 2, σ 2 = j τ 2, τ 2 = 2 log 2 j..
Consider a function Gx, y, s, t of four variables x, y, s and t. Suppose Gx, y, s, t is approximated using two-dimensional Haar wavelet as follows []: 7 3.7 Gx, y, s, t p= q= b p,q x, yh p sh q t. Substituting the collocation points s i = i 0.5 2M and t j = j 0.5 2N we obtain the linear system 3.8 Gx, y, s i, t j p= q=, i =, 2,..., 2M,, j =, 2,..., 2N, b p,q x, yh p s i h q t j, i =, 2,..., 2M, j =, 2,..., 2N. Corollary. The solution of the system 3.8 for any value of x.y [0, ] is given as follows []: b i, x, y = b,j x, y = b i,j x, y = b, x, y = ρ 2N p=α M 2M ρ 2 ρ ρ 2 p= p=α q= 2M 2N Gx, y, s p, t q Gx, y, s p, t q Gx, y, s p, t q p= p= q= Gx, y, s p, t q, p= q= p= q= 2 Gx, y, s p, t q p=α q= 2 Gx, y, s p, t q Gx, y, s p, t q Gx, y, s p, t q p= q= 2 Gx, y, s p, t q, i = 2, 3,..., 2M,, j = 2, 3,..., 2N, i = 2, 3,..., 2M, j = 2, 3,..., 2N, where α,, and ρ are defined as in Eq. 3.5 and α 2, 2, 2 and ρ 2 are defined as in Eq. 3.6. Corollary 2. Suppose a function Gx, y of two variables x and y is approximated using Haar wavelet approximation given in Eq. 3.3. Suppose further that Gx, y is known at collocation points x m, y m, m =, 2,..., 2M, n =, 2,..., 2N. Then the approximate value of the function Gx, y at any other point of the domain can be,
8 calculated as follows []: Gx, y = 2M 2N j= ρ 2N p= q= 2M ρ 2 Gx m, y m h xh y p=α q= 2M p= j= p= Gx m, y m Gx m, y m ρ ρ 2 p=α Gx m, y m p= q= p= q= 2 Gx m, y m Gx m, y m p= q= 2 Gx m, y m p=α q= 2 Gx m, y m h i xh y h xh j y F x m, y m h i xh j y, where α,, and ρ are defined as in Eq. 3.5 and α 2, 2, 2 and ρ 2 are defined as in Eq. 3.6. 3.2 Two-dimensional linear stochastic Volterra integral equation Consider the two-dimensional linear stochastic Volterra integral equation.2. Assume that the function Kx, y, s, tgs, t is approximated using two-dimensional Haar wavelet as follows: 3.9 K x, y, s, tgs, t 3.0 K 2 x, y, s, tgs, t j= j= b i,j x, yh i sh j t. c i,j x, yh i sh j t. With this approximation Eq..2 can be writen as follows: 3. gx, y = fx, y y x 0 0 y x 0 0 j= b i,j x, yh i sh j tdsdt j= c i,j x, yh i sh j tdbsdbt. Eq. 3. can be written in a more compact form using the notations introduced in equations 2. and 2.2 and is given as follows: gx, y = fx, y j= b i,j x, yp i, xp j, y j= c i,j x, yq i, xq j, y.
Substituting the collocation points given in 3. and 3.2, we obtain the following system of equations: gx m, y n = fx m, y n j= b i,j x m, y n p i, x m p j, y n j= c i,j x m, y n q i, x m q j, y n. Now b i,j, i =, 2,..., 2M, j =, 2,..., 2N and similarly c i,j, i =, 2,..., 2M, j =, 2,..., 2N can be replaced with their expressions given in Corollary and the following system of equations is obtained: 3.2 gx m, y n = fx m, y n p,x m p, y n 2M 2N i=2 j=2 i=2 p i, x m p, y n ρ 2N p, x m p j, y n 2M ρ 2 j=2 p=α q= 2 q, x m q, y n 2M 2N i=2 p i, x m p j, y n ρ ρ 2 p=α q= 2M p= p= q= K x m, y n, s p, t q gs p, t q p= q= K x m, y n, s p, t q gs p, t q K x m, y n, s p, t q gs p, t q p=α 2 K x m, y n, s p, t q gs p, t q p=α K x m, y n, s p, t q gs p, t q p= q= q i, x m q, y n ρ 2N p= q= 2 K x m, y n, s p, t q gs p, t q K x m, y n, s p, t q gs p, t q p= p= q= 2 K 2 x m, y n, s p, t q gs p, t q p=α q= 9 K x m, y n, s p, t q gs p, t q K x m, y n, s p, t q gs p, t q K 2 x m, y n, s p, t q gs p, t q p= q= K 2 x m, y n, s p, t q gs p, t q
0 j=2 i=2 q, x m q j, y n 2M ρ 2 j=2 p=α q= 2 q i, x m q j, y n ρ ρ 2 2M p= p=α 2 K 2 x m, y n, s p, t q gs p, t q p=α K 2 x m, y n, s p, t q gs p, t q p= q= 2 K 2 x m, y n, s p, t q gs p, t q p= p= q= 2 K 2 x m, y n, s p, t q gs p, t q K 2 x m, y n, s p, t q gs p, t q K 2 x m, y n, s p, t q gs p, t q m =, 2,..., 2M, n =, 2,..., 2N. Eq. 3.2 represents 2M 2N system which can be solved using either prevalent methods for solving linear systems. The solution of this system gives values of gx, y at the collocation points. The values of gx, y at points other than collocation points can be calculated using Corollary 2., 4. Numerical Example In this section, the numerical example is given to demonstrate the applicability and accuracy of our method. Consider the following linear 2D stochastic Volterra integral equation of second kind: ux, y = fx, y where y x xyt sus, tdsdt y x 0 0 0 0 xytsus, tdbsdbt fx, y = x y 2 xyx3 4x 2 y 4xy 2 y 3. The solutions mean together confidence interval at the collocation points for the present method for 000 iterative of system 3.2 is shown in Table. In Figs. 4, three-dimensional graphs of the approximate solution for various values of level L are shown. Table :The solutions mean together confidence interval for above example
J M 2M x, y ūx, y Confidence L Interval U 0 2 0.25, 0.75.995.995.995 0.25, 0.375.0677.0675.0679 2 4 0.375, 0.875 2.4348 2.4394 2.43769 0.625, 0.875 2.7259 2.72036 2.7347 0.0625, 0.4375.0498.0498.0498 2 4 8 0.325, 0.6875 2.5292 2.5266 2.538 0.825, 0.9375 2.68057 2.67707 2.68407 0.0325, 0.7875.5492.5488.5495 3 8 6 0.40625, 0.5325 2.25846 2.25532 2.2659 0.7825, 0.96875 2.6747 2.6739 2.6756 0.05625, 0.609375.2632.2629.2635 4 6 32 0.296875, 0.796875 2.25999 2.2589 2.2607 0.859375, 0.984375 2.58876 2.58870 2.5888
2 Figure. Plot of approximate solution of level L = for test example. Figure 2. Plot of approximate solution of level L = 2 for test example.
3 Figure 3. Plot of approximate solution of level L = 3 for test example. Figure 4. Plot of approximate solution of level L = 4 for test example.
4 5. Conclusion As mentioned above, numerical solution of two-dimensional stochastic integral equations because of the randomness is very difficult or sometimes impossible. In this paper, we have successfully developed Haar wavelets numerical method for approximate a solution of two-dimensional linear stochastic Volterra integral equations. The example confirm that the method is considerably fast and highly accurate as sometimes lead to exact solution. Although, theoretically for getting higher accuracy we can set the method with larger values of M and N and also larger of the degree of approximation, p and q, but it leads to solving MN linear systems of size pq pq, that have its difficulties. The method can be improved to be more accurate by using other numerical methods. Mathematica has been used for computations in this paper.
5 References [] I. Aziz, Siraj-ul-Islam, F. Khan, A new method based on Haar wavelet for numerical solution of two-dimensional nonlinear integral equations, J. Comp. Appl. Math. 272 204, 70-80. [2] I. Aziz, Siraj-ul-Islam, New algorithms for numerical solution of nonlinear Fredholm and Volterra integral equations using Haar wavelets, J. Comp. Appl. Math. 239 203 333-345. [3] Siraj-ul-Islam, I. Aziz, M. Fayyaz, A new approach for numerical solution of integro-differential equations via Haar wavelets, Int, J. Comp. Math. 90 203 97-989. [4] Siraj-ul-Islam, I. Aziz, A. Al-Fhaid, An improved method based on Haar wavelets for numerical solution of nonlinear and integro-differential equations of first and higher orders, J. Comp. Appl. Math. 260 204 449-469. [5] Kuo, Hui-Hsiung, Introduction to stochastic integration, Springer ScienceBusiness Media, Inc. 2006. [6] K. E. Atkinson, The numerical solution of integral equations of the second kind, Cambridge University Press 997. [7] A. J. Jerri, Introduction to integral equations with applications, John Wiley and Sons, INC 999. [8] T. S. Sankar, V. I. Fabrikant, Investigations of a two-dimentional integral equation in the theory of elasticity and electrostatics, J. Mec. Theor. Appl. 2 983 285-299. [9] H. J. Dobner, Bounds for the solution of hyperbolic problems, Computing 38 987 209-28. [0] V. M. Aleksandrov, A. V. Manzhirov, Two-dimentional integral equations in applied mechanics of deformable solids, J. Appl. Mech. Tech. Phys. 5 987 46-52. [] A. V. Manzhirov, Contact problems of the interaction between viscoelastic foundations subject to ageing and systems of stamps not applied simultaneously, Prikl. Matem. Mekhan. 4 987 523-535. [2] O. V. Soloviev, Low-frequency radio wave propagation in the earth-ionosphere waveguide disturbed by a large-scale three-dimensional irregularity, Radiophysics and Quantum Electronics 4 998 392-402. [3] M. Rahman, A rigid elliptical disc-inclusion in an elastic solid, subject to a polynomial normal shift, J. Elasticity 66 2002 207-235.
6 [4] H. Guoqiang, W. Jiong, Extrapolation of nystrom solution for two dimentional nonlinear Fredholm integral equations, J. Comp. App. Math. 34 200 259-268. [5] H. Brunner, Collocation methods for Volterra integral and related functional equations, Cambridge University Press, 2004. [6] H. Guoqiang, K. Itayami, K. Sugihara, W. Jiong, Extrapolation method of iterated collocation solution for two-dimentional nonlinear Volterra integral equations, Appl. Math. Comput. 2 2000 49-6. [7] W. Xie, F. R. Lin, A fast numerical solution method for two dimensional Fredholm integral equations of the second kind, App. Num. Math. 59 2009 709-79. [8] S. Bazm, E. Babolian, Numerical solution of nonlinear two-dimensional Fredholm integral equations of the second kind using Gauss product quadrature rules, Commun. Nonlinear Sci. Numer. Simult. 7 202 25-223. [9] G. Han, R. Wang, Richardson extrapolation of iterated discrete Galerkin solution for twodimensional Fredholm integral equations, J. Comp. App. Math. 39 2002 49-63. [20] K. Maleknejad, Z. JafariBehbahani, Application of two-dimensional triangular functions for solving nonlinear class of mixed Volterra-Fredholm integral equations, Math. Comp. Mode. 55 202 833-844. [2] E. Babolian, K. Maleknejad, M. Roodaki, H. Almasieh, Two dimensional triangular functions and their applications to nonlinear 2d Volterra-Fredholm equations, Comp. Math. App. 60 200 7-722. [22] S. Nemati, P. Lima, Y. Ordokhani, Numerical solution of a class of two-dimensional nonlinear Volterra integral equations using legender polynomials, J. Comp. Appl. Math. 242 203 53-69. [23] A. Tari, M. Rahimi, S. Shahmorad, F. Talati, Solving a class of two-dimensional linear and nonlinear Volterra integral equations by the differential transform method, J. Comp. Appl. Math. 228 2009 70-76. [24] P. Assari, H. Adibi, M. Dehghal, A meshless method for solving nonlinear two-dimensional integral equations of the second kind on non-rectangular domains using radial basis functions with error analysis, J. Comp. Appl. Math. 239 203 72-92. [25] M. H. Reihani, Z. Abadi, Rationalized Haar functions method for solving Fredholm and Volterra integral equations, J. Comp. Appl. Math. 200 2007 2-20.
7 [26] F. Hosseini Shekarabi, K. Maleknejad, R. Ezzati, Application of two-dimensional Bernstein polynomials for solving mixed Volterra-Fredholm integral equations, African Mathematical Union and Springer-Verlag Berlin Heidelberg, DOI 0. 007/s 3370-04-0283-6 204. [27] F. Keinert, Wavelets and Multiwavelets, A Crc Press Company Boca Raton London New York Washington, D. C, 2004.