一类变系数分数阶微分方程的数值解法（英文）

本文研究了一类变系数分数阶微分方程的数值解法问题.利用Cheyshev小波推导出的分数阶微分方程的算子矩阵把分数阶微分方程转换为代数方程组.同时给出了Cheyshev小波基的收敛性和误差估计表达式,并给出数值算例说明所提方法的精确性和有效性.     

一个新的分数阶比较定理

研究了Caputo和Riemann-Liouville两型分数阶微分方程的比较定理.首先,讨论了一类线性分数阶微分不等式解得非负性.其次,引入单边Lipschitz条件,将微分方程解的比较问题化为线性微分不等式非负解问题,通过线性分数阶微分方程的求解,得到分数阶比较定理.最后,为进一步说明结论,给出了两个数值仿真例子.     

分数阶尺度函数的分解与重构算法

引入分数阶多分辨分析与分数阶尺度函数的概念.运用时频分析方法与分数阶小波变换,研究了分数阶正交小波的构造方法,得到分数阶正交小波存在的充要条件.给出分数阶尺度函数与小波的分解与重构算法,算法比经典的尺度函数与小波的分解与重构算法更具有一般性.     

分数阶常微分方程迭代方法的解

通过采用分数阶积分与导数的复合,把分数阶常微分方程转化为积分方程.构造出迭代格式,证明它的收敛性,进一步给出近似解的误差估计.并给出数值例子.     

分数阶积分微分方程的近似解

本文给出了分数阶积分微分方程的一种新的解法.利用未知函数的泰功多项式展开将分数阶积分微分方程近拟转化为一个涉及未知函数及其n阶导数的线性方程组.数值例子表明该方法的有效性.     

分数阶延迟微分方程的样条配置方法

首次利用三次样条配置方法采用直接法求解了一类非线性分数阶延迟微分方程初值问题,并给出了方法的局部截断误差和若干数值算例.数值结果表明方法求解分数阶延迟微分方程初值问题是非常有效的,结果对于未来研究分数阶延迟微分方程的数值方法具有重要的意义.     

随机分数阶微分方程初值问题基于模拟方程法的数值求解

基于模拟方程法,提出了一种求解随机分数阶微分方程初值问题的数值方法.考虑含两个分数阶导数项的微分方程,引入两个线性的、非耦合的随机模拟方程,利用它们解构原方程,借助Laplace变换及逆变换,得到方程解的积分表达式,同时建立起两个模拟方程之间的联系,结合初始状态,得到求解随机微分方程初值问题的数值迭代算法.作为特例,对于含两个分数阶导数项线性常微分方程的初值问题,给出了基于模拟方程法的数值解法的显式结果.该方法是稳定的,它的误差仅存在于积分近似时的截断误差和计算软件的舍入误差.应用实例说明了数值方法在确定和随机情形的有效性和准确性.     

线性随机分数阶微分方程Euler方法的弱收敛性与弱稳定性

本文主要研究了线性随机分数阶微分方程Euler方法的弱收敛性与弱稳定性.首先构造了数值求解线性随机分数阶微分方程的Euler方法,然后证明该方法是弱稳定的和α阶弱收敛的,文末给出的数值算例验证了所获得的理论结果的正确性.     

几类分数阶多项比例延迟微分方程的Jacobi配置方法

本文利用Jacobi配置方法数值求解几类分数阶多项比例延迟微分方程初值问题,给出相应的误差分析,并利用若干数值算例验证了相应的理论结果,表明Jacobi配置方法求解这几类分数阶比例延迟方程是高效的.同时,也为分数阶泛函微分方程的数值算法提供新的研究思路.     

分数阶R-L微分方程初值问题解的存在唯一性

研究了R-L导数定义下的分数维微分方程初值问题解的存在性及其唯一性,给出了方程的Peano存在定理和不等式定理,基于逐次逼近的方法,利用对分数阶R-L微夯方程构造的Tonelli序列和Ascoli引理证明分数阶R-L微分方程解的存在性,根据分数阶不等式定理证明了分数阶R-L微分方程解的唯一性.     

The operational matrix of fractional integration for shifted Chebyshev polynomials

A new shifted Chebyshev operational matrix (SCOM) of fractional integration of arbitrary order is introduced and applied together with spectral tau method for solving linear fractional differential equations (FDEs). The fractional integration is described in the Riemann–Liouville sense. The numerical approach is based on the shifted Chebyshev tau method. The main characteristic behind the approach using this technique is that only a small number of shifted Chebyshev polynomials is needed to obtain a satisfactory result. Illustrative examples reveal that the present method is very effective and convenient for linear multi-term FDEs.     

Solving a nonlinear fractional differential equation using Chebyshev wavelets

Chebyshev wavelet operational matrix of the fractional integration is derived and used to solve a nonlinear fractional differential equations. Some examples are included to demonstrate the validity and applicability of the technique.     

On the applicability of Genocchi wavelet method for different kinds of fractional‐order differential equations with delay

A novel collocation method based on Genocchi wavelet is presented for the numerical solution of fractional differential equations and time‐fractional partial differential equations with delay. In this work, to achieve the approximate solution with height accuracy, we employed the operational matrix of integer derivative and the pseudo‐operational matrix of fractional derivative in Caputo sense. Also, based on Genocchi function properties, we presented delay and pantograph operational matrices of Genocchi wavelet functions (GWFs). Due to operational and pseudo‐operational matrices, the equations under this study can be turned into nonlinear algebraic equations with the unknown GWF coefficients. For illustrating the upper bound of error for the proposed method, we estimate the error in the sense of Sobolev space. In addition, to demonstrate the efficacy of the pseudo‐operational matrix of fractional derivative, we investigate the upper bound of error for the mentioned matrix. Finally, the algorithm based on the proposed approach is implemented for some numerical experiments to confirm accuracy and applicability.     

A fast numerical algorithm based on the Taylor wavelets for solving the fractional integro‐differential equations with weakly singular kernels

In this paper, a fast numerical algorithm based on the Taylor wavelets is proposed for finding the numerical solutions of the fractional integro‐differential equations with weakly singular kernels. The properties of Taylor wavelets are given, and the operational matrix of fractional integration is constructed. These wavelets are utilized to reduce the solution of the given fractional integro‐differential equation to the solution of a linear system of algebraic equations. Also, convergence of the proposed method is studied. Illustrative examples are included to demonstrate the validity and applicability of the technique.     

Numerical solution of two-dimensional time fractional cable equation with Mittag-Leffler kernel

The main motive of this article is to study the recently developed Atangana-Baleanu Caputo (ABC) fractional operator that is obtained by replacing the classical singular kernel by Mittag-Leffler kernel in the definition of the fractional differential operator. We investigate a novel numerical method for the nonlinear two-dimensional cable equation in which time-fractional derivative is of Mittag-Leffler kernel type. First, we derive an approximation formula of the fractional-order ABC derivative of a function t^k using a numerical integration scheme. Using this approximation formula and some properties of shifted Legendre polynomials, we derived the operational matrix of ABC derivative. In the author of knowledge, this operational matrix of ABC derivative is derived the first time. We have shown the efficiency of this newly derived operational matrix by taking one example. Then we solved a new class of fractional partial differential equations (FPDEs) by the implementation of this ABC operational matrix. The two-dimensional model of the time-fractional model of the cable equation is solved and investigated by this method. We have shown the effectiveness and validity of our proposed method by giving the solution of some numerical examples of the two-dimensional fractional cable equation. We compare our obtained numerical results with the analytical results, and we conclude that our proposed numerical method is feasible and the accuracy can be seen by error tables. We see that the accuracy is so good. This method will be very useful to investigate a different type of model that have Mittag-Leffler fractional derivative.     

Numerical Solution for a Variable-Order Fractional Nonlinear Cable Equation via Chebyshev Cardinal Functions

In this paper, a variable-order fractional derivative nonlinear cable equation is considered. It is commonly accepted that fractional differential equations play an important role in the explanation of many physical phenomena. For this reason we need a reliable and efficient technique for the solution of fractional differential equations. This paper deals with the numerical solution of class of fractional partial differential equation with variable coefficient of fractional differential equation in various continues functions of spatial and time orders. Our main aim is to generalize the Chebyshev cardinal operational matrix to the fractional calculus. Finally, illustrative examples are included to demonstrate the validity and applicability of the presented technique.     

Vieta–Fibonacci wavelets: Application in solving fractional pantograph equations

In this paper, the Vieta–Fibonacci wavelets as a new family of orthonormal wavelets are generated. An operational matrix concerning fractional integration of these wavelets is extracted. A numerical scheme is established based on these wavelets and their fractional integral matrix together with the collocation technique to solve fractional pantograph equations. The presented method reduces solving the problem under study into solving a system of algebraic equations. Several examples are provided to show the accuracy of the method.     

Numerical solution of nonlinear fractional-order Volterra integro-differential equations by SCW

Fractional calculus is an extension of derivatives and integrals to non-integer orders and has been widely used to model scientific and engineering problems. In this paper, we describe the fractional derivative in the Caputo sense and give the second kind Chebyshev wavelet (SCW) operational matrix of fractional integration. Then based on above results we propose the SCW operational matrix method to solve a kind of nonlinear fractional-order Volterra integro-differential equations. The main characteristic of this approach is that it reduces the integro-differential equations into a nonlinear system of algebraic equations. Thus, it can simplify the problem of fractional order equation solving. The obtained numerical results indicate that the proposed method is efficient and accurate for this kind equations.     

Haar wavelet operational matrix of fractional order integration and its applications in solving the fractional order differential equations

Haar wavelet operational matrix has been widely applied in system analysis, system identification, optimal control and numerical solution of integral and differential equations. In the present paper we derive the Haar wavelet operational matrix of the fractional order integration, and use it to solve the fractional order differential equations including the Bagley-Torvik, Ricatti and composite fractional oscillation equations. The results obtained are in good agreement with the existing ones in open literatures and it is shown that the technique introduced here is robust and easy to apply.     

Numerical technique based on generalized Laguerre and shifted Chebyshev polynomials

In this study, we present a numerical scheme for solving a class of fractional partial differential equations. First, we introduce psi -Laguerre polynomials like psi-shifted Chebyshev polynomials and employ these newly introduced polynomials for the solution of space-time fractional differential equations. In our approach, we project these polynomials to develop operational matrices of fractional integration. The use of these orthogonal polynomials converts the problem under consideration into a system of algebraic equations. The solution of this system provide us the desired results. The convergence of the proposed method is analyzed. Finally, some illustrative examples are included to observe the validity and applicability of the proposed method.