An efficient approach based on radial basis functions for solving stochastic fractional differential equations

In this paper, we present a collocation method based on Gaussian Radial Basis Functions (RBFs) for approximating the solution of stochastic fractional differential equations (SFDEs). In this equation the fractional derivative is considered in the Caputo sense. Also we prove the existence and uniqueness of the presented method. Numerical examples confirm the proficiency of the method.


Introduction
Fractional calculus introduced because it can fill the existing gap for describing a large amount of work in engineering [1,2], and different phenomena in nature such as biology, physics [3,4]. Mathematicians and physicists have been created numerous articles about fractional differential equations(FDEs) for finding analytical and numerical methods, including Adomian Decomposition Method [5], Variational Iteration Method [6,7], Homotopy perturbation Method [8] and homotopy analysis method [9]. H. Rezazadeh et al. have generalized the Floquet system to the fractional Floquet system in 2016 [10]. Stochastic Differential Equations (SDEs) models play a great role in various sciences such as physics, economics, biology, chemistry and finance [11][12][13][14][15]. The reader has at least knowledge about independence, expected values and variances and also basic definitions of stochastic, that is necessary to read articles in this field [16]. M.Khodabin et al. approximate solution of stochastic Volterra integral equations in 2014 [17], also R.Ezzati et al. work on a stochastic operational matrix based on block pulse functions in 2014 [18]. We introduce SFDEs [19]: for 0 a 1 and t 2 0; T ½ , where D a is the Caputo fractional derivative of order a which will be defined later. r is Max amplitude of noise also R t t 0 gðt; sÞdwðsÞ is the stochastic term, that produce some noise in our result, throughout the paper we putting r ¼ 1. SFDEs play a remarkable role for physical applications in nature [20][21][22].
The motivation of this paper is to extend the application of the RBF to solve SFDEs. The layout of the paper is the following. In Sect. 2 some essential definitions of fractional calculus is proposed. In Sect. 3 we explain using RBFs method for SFDEs and prove the existence and uniqueness of the presented method. In Sect. 4 various examples are solved to illustrate the effectiveness of the proposed method. Also a conclusion is given in the last section.

Preliminaries and notations
In this section, we give some basic definitions and properties of fractional calculus which are defined as follow [4].
The Caputo fractional derivative of order m is defined as x [ 0 where D p is the classical differential operator of order p.

Remark 2.2 For the Caputo derivative we have
Definition 2.4 Let ðX; F; qÞ be a probability space with a normal filteration ðF t Þ t ! 0 and w ¼ fwðtÞ : t ! 0g be a Brownian motion defined over this filtered probability space. Consider the following SFDE D a uðt; mÞ ¼ f ðt; uðt; mÞÞ þ Z t 0 gðt; s; mÞdwðsÞ; uð0; mÞ ¼ u 0 for t 2 ½0; T, and m 2 X. For simplicity of notation we drop the variable m so we have the following equation gðt; sÞdwðsÞ; therefore, we have Also we admit the following assumptions.
Assumption 2.5 Suppose f and g are L 2 measurable functions satisfying f ðm; xÞ À f ðm; yÞ j j K 1 x À y j j; gðt; mÞ À gðt; nÞ j j K 2 ð m À n j jÞ; ð3Þ for some constants K 1 ; K 2 ; K 3 ; K 4 and for every x; y 2 R and 0 m; n t T ¼ 1: Stochastic integral 2.6 Now we should explain the approximation of the stochastic term. White noise is known as the derivative of the brownian motion W(s) [30], so we approximating the term dw s dt .
where c i $ Nð0; 1Þ is introducing by

Collocation method based on RBFs for solving SFDEs
The Radial basis functions method has been known as a powerful tool for solving ordinary, partial and fractional differential equations and also integral equations and etc. So in this section we use this method for solving (2). Before that we consider some preliminaries.

Interpolation by RBFs
Let ft 1 ; . . .; t N g be a given set of distinct points in 0; T ½ R. Then the approximation of a function u(t) using RBFs uðtÞ ¼ uð t k kÞ, can be written in the following form [32,33]  Using RBFs for solving SFDEs (2) Let D ¼ ½0; T R and u : Cð½0; TÞ ! R and also suppose that u N is the approximation of u based on these functions so we can write N is the number of nodal points within the domain D and c i denotes the shape parameter. Also we know there are different kinds of RBFs, but in this research we need only one of them with titled Gaussian that we represented as follow: The collocation method based on RBF basis for solving (2) can be written in the following form: where ðt À mÞ aÀ1 dm K 1 aCðaÞ ju 1 À u 2 j: On the other side using the assumption of the lemma we know that K 1 aCðaÞ \1, therefore the proof is completed. h

Illustrative example
In this section, we solve SFDEs using RBFs and Galerkin method [19], these equations don't have exact solution so we use numerical approximation for sufficiently small partition on t. While we have this relation je exact À e RBFs j ¼ je exact À e Galerkin þ e Galerkin À e RBFs j je exact À e Galerkin j þ je Galerkin À e RBFs j then we approximate the solution in the form of RMSError [34] as follow In this paper, we introduce RMSError after 50 and 60 times run the program for different points and r.
Example 4.1 Consider the following SFDE: uð0Þ ¼ 1: We have Tables 1, 2 and 3 after 60 times run the program with n ¼ 17.
Example 4.2 Consider the following SFDE: D a uðtÞ þ uðtÞ ¼ 2t 2Àa Cð3 À aÞ À t 1Àa Cð2 À aÞ þ t 2 À t þ r Z t 0 sdwðsÞ; uð0Þ ¼ 0; 0\a 1: For a ¼ 1 2 and various r we have Table 4 after 50 times run the program with nodal points n ¼ 12 and for a ¼ 3 2 and different value for r we have Table 5 after 50 times run the program with n ¼ 11.    In this work, the accuracy of approximate solution, when taking larger n and smaller r, is expected that more accurate the approximate results.

Conclusion
The main goal of this work was to purpose an efficient algorithm for the stochastic fractional differential equations. In this paper, while we don't have exact solution for SFDEs we used RBFs to approximate the solution of these kind of equations. In addition, we discussed about existence and uniqeness of the presented method. The present RMS Error in the tables shows that the results are highly accurate in comparison with another method using by Galerkin algorithm (Tables 6, 7).
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://crea tivecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.