带有弱奇性核的多项分数阶非线性随机微分方程的改进Euler-Maruyama格式

针对一类带有弱奇性核的多项分数阶非线性随机微分方程构造了改进Euler-Maruyama (EM)格式,并证明了该格式的强收敛性.具体地,利用随机积分解的充分条件,将此多项分数阶随机微分方程等价地转化为随机Volterra 积分方程的形式,详细推导出对应的改进EM格式,并对该格式进行了强收敛性分析,其强收敛阶为α_m-α_m-1,其中α_i为分数阶导数的指标,且满足0<α₁<…<α_m-1<α_m<1.最后,通过数值实验验证了理论分析结果的正确性.     

一类带乘性噪声随机分数阶微分方程Euler方法的弱收敛性与弱稳定性

本文主要在带加性噪声随机分数阶微分方程的基础上,研究了一类更为困难的带乘性噪声随机分数阶微分方程Euler方法的弱收敛性与弱稳定性,并得到了类似的结论.首先构造了数值求解带乘性噪声随机分数阶微分方程的Euler方法,然后证明当分数阶α满足0α1/2时,该方法是1/2-α阶弱收敛的和弱稳定的,文末数值试验的结果验证了理论结果的正确性.     

时间分数阶扩散方程的数值解法

分数阶微分方程在许多应用科学上比整数阶微分方程更能准确地模拟自然现象.考虑时间分数阶扩散方程,将一阶的时间导数用分数阶导数α(0<α<1)替换,给出了一种计算有效的隐式差分格式,并证明了这个隐式差分格式是无条件稳定和无条件收敛的,最后用数值例子说明差分格式是有效的.     

非线性随机分数阶延迟积分微分方程Euler-Maruyama方法的强收敛性

本文研究了一类新的模型问题:非线性随机分数阶延迟积分微分方程.当方程中的漂移项和扩散项满足全局Lipschitz条件和线性增长条件时,基于压缩映射原理给出了该方程解存在唯一的充分条件.由于理论求解的困难,构造了一种数值方法(Euler-Maruyama方法),并证得强收敛阶为α-1/2,α∈(1/2,1].最后通过数值试验,验证了这一理论结果.     

分数阶波方程的数值解法

首先,把分数阶波方程转换成等价的积分-微分方程;然后,利用带权的分数阶矩形公式和紧差分算子分别对时间和空间方向进行离散.证明了当权重为1/2时,时间方向的收敛阶为α,其中α(1α2)为Caputo导数的阶数.利用Gronwall不等式,证明了数值格式的收敛性和稳定性.数值例子进一步表明了数值格式的有效性.     

求解带有非线性边界条件的涡流方程的A-φ解耦有限元格式

本文研究边界条件符合幂指数型非线性关系H × n = n × (|E × n|^α-1E × n)(0 < α ≤ 1)的涡流方程.使用A-φ耦合有限元格式数值求解这类问题具有较高精度,但计算开销大. A-φ解耦有限元计算格式能够在每个时间步上分别求解矢量A和标量φ,以此降低计算规模,提高计算效率.我们证明了解耦格式中解的存在唯一性,并且给出了它的误差估计.最后给出的数值实验证明了本文所提供的解耦算法是稳定和有效的.     

带有弱奇性核的多项分数阶非线性随机微分方程解的适定性

本文主要研究一类带有多项分数阶Caputo导数的非线性随机微分方程初值问题的解的适定性.具体地,首先把多项分数阶随机微分方程等价地转化为随机Volterra积分方程;然后,给出了该随机积分方程的Euler-Maruyama (EM)格式;最后,借助于该EM格式,证明了多项分数阶随机微分方程的解的适定性.     

非线性中立型随机延迟微分方程随机θ方法的稳定性

本文研究非线性中立型随机延迟微分方程随机θ方法的均方稳定性.在方程解析解均方稳定的条件下,证明了如下结论:当θ∈[0,1/2)时,随机θ方法对于适当小的时间步长是均方稳定的;当θ∈[1/2,1]时,随机θ方法对于任意步长都是均方稳定的.数值结果验证了所获结论的正确性.     

随机微分方程改进的分裂步单支θ方法的强收敛性

本文提出了一个改进的分裂步单支θ方法,在漂移项系数满足单边Lipschitz条件下,证明了当数值方法的参数θ满足1/2 ≤ θ ≤ 1时,该数值方法对于这类随机微分方程是强收敛的,并在现有文献的基础上将方法的收敛阶从1/2阶提高到1阶;当0 ≤ θ ≤ 1/2时,若漂移项系数进一步满足线性增长条件,该数值方法也是强收敛的,收敛阶为1阶.文末的数值试验验证了理论结果的正确性.     

一类分数阶Langevin方程block-by-block算法的数值分析

分数阶Langevin方程有重要的科学意义和工程应用价值，基于经典block-by-block算法，求解了一类含有Caputo导数的分数阶Langevin方程的数值解.Block-by-block算法通过引入二次Lagrange基函数插值，构造出逐块收敛的非线性方程组，通过在每一块耦合求得分数阶Langevin方程的数值解.在0<α<1条件下，应用随机Taylor展开证明block-by-block算法是3+α阶收敛的，数值试验表明在不同α和时间步长h取值下，block-by-block算法具有稳定性和收敛性，克服了现有方法求解分数阶Langevin方程速度慢精度低的缺点，表明block-by-block算法求解分数阶Langevin方程是高效的.     

Higher order numerical methods for solving fractional differential equations

In this paper we introduce higher order numerical methods for solving fractional differential equations. We use two approaches to this problem. The first approach is based on a direct discretisation of the fractional differential operator: we obtain a numerical method for solving a linear fractional differential equation with order 0<α<1. The order of convergence of the numerical method is O(h ^3?α ). Our second approach is based on discretisation of the integral form of the fractional differential equation and we obtain a fractional Adams-type method for a nonlinear fractional differential equation of any order α>0. The order of convergence of the numerical method is O(h ³) for α≥1 and O(h ^1+2α ) for 0<α≤1 for sufficiently smooth solutions. Numerical examples are given to show that the numerical results are consistent with the theoretical results.     

Stochastic Evolution Equations Driven by a Fractional White Noise

Abstract

In this paper we study stochastic evolution equations driven by a fractional white noise with arbitrary Hurst parameter in infinite dimension. We establish the existence and uniqueness of a mild solution for a nonlinear equation with multiplicative noise under Lipschitz condition by using a fixed point argument in an appropriate inductive limit space. In the linear case with additive noise, a strong solution is obtained. Those results are applied to stochastic parabolic partial differential equations perturbed by a fractional white noise.     

Solving systems of fractional differential equations using differential transform method

This paper presents approximate analytical solutions for systems of fractional differential equations using the differential transform method. The fractional derivatives are described in the Caputo sense. The application of differential transform method, developed for differential equations of integer order, is extended to derive approximate analytical solutions of systems of fractional differential equations. The solutions of our model equations are calculated in the form of convergent series with easily computable components. Some examples are solved as illustrations, using symbolic computation. The numerical results show that the approach is easy to implement and accurate when applied to systems of fractional differential equations. The method introduces a promising tool for solving many linear and nonlinear fractional differential equations.     

Existence of solutions for nonlinear fractional stochastic differential equations

The fractional stochastic differential equations have wide applications in various fields of science and engineering. This paper addresses the issue of existence of mild solutions for a class of fractional stochastic differential equations with impulses in Hilbert spaces. Using fractional calculations, fixed point technique, stochastic analysis theory and methods adopted directly from deterministic fractional equations, new set of sufficient conditions are formulated and proved for the existence of mild solutions for the fractional impulsive stochastic differential equation with infinite delay. Further, we study the existence of solutions for fractional stochastic semilinear differential equations with nonlocal conditions. Examples are provided to illustrate the obtained theory.     

Principal resonance responses of SDOF systems with small fractional derivative damping under narrow-band random parametric excitation

The principal resonance responses of nonlinear single-degree-of-freedom (SDOF) systems with lightly fractional derivative damping of order α (0 < α < 1) subject to the narrow-band random parametric excitation are investigated. The method of multiple scales is developed to derive two first order stochastic differential equation of amplitude and phase, and then to examine the influences of fractional order and intensity of random excitation on the first-order and second-order moment. As an example, the stochastic Duffing oscillator with fractional derivative damping is considered. The effects of detuning frequency parameter, the intensity of random excitation and the fractional order derivative damping on stability are studied through the largest Lyapunov exponent. The corresponding theoretical results are well verified through direct numerical simulations. In addition, the phenomenon of stochastic jump is analyzed for parametric principal resonance responses via finite differential method. The stochastic jump phenomena indicates that the most probable motion is around the larger non-trivial branch of the amplitude response when the intensity of excitation is very small, and the probable motion of amplitude responses will move from the larger non-trivial branch to trivial branch with the increasing of the intensity of excitation. Such stochastic jump can be considered as bifurcation.     

The minimax approach for a class of variable order fractional differential equation

This paper introduces an approximate solution for Liouville‐Caputo variable order fractional differential equations with order 0 < α(t) ≤ 1 . The solution is adapted using a family of fractional‐order Chebyshev functions with unknown coefficients. These coefficients have been obtained by using an optimization approach based on minimax technique and the least pth optimization function. Several linear and nonlinear fractional‐order differential equations are discussed using the proposed technique for fixed and variable order fractional‐order derivatives. Moreover, the response of RC charging circuit with variable order fractional capacitor is studied for different cases. Several comparisons with related published techniques have been added to illustrate the accuracy of the proposed approach.     

Poisson和公式在分数阶热方程中的应用

本文利用Poisson和公式,证明了如下分数阶热方程(D_t~αlu=D_x~2u u(x1 0)=f(x))当f分别为周期函数和f∈S(■)时(速降函数空间),它们的热核满足关系H_t~α(x)=∑n=-∞H_t~α(x+n)进一步,我们把结论推广到更一般的分数阶微分方程和高维情形     

Efficient numerical schemes for the solution of generalized time fractional Burgers type equations

The main aim of this paper is to propose two semi-implicit Fourier pseudospectral schemes for the solution of generalized time fractional Burgers type equations, with an analysis of consistency, stability, and convergence. Under some assumptions, the unconditional stability of the schemes is shown. In implementation of these schemes, the fast Fourier transform (FFT) can be used efficiently to improve the computational cost. Various test problems are included to illustrate the results that we have obtained regarding the proposed schemes. The results of numerical experiments are compared with analytical solutions and other existing methods in the literature to show the efficiency of proposed schemes in both accuracy and CPU time. As numerical solution of fractional stochastic nonlinear partial differential equations driven by Brownian motions are among current related research interests, we report the performance of these schemes on stochastic time fractional Burgers equation as well.     

Numerical approach to differential equations of fractional order

In this paper, the variational iteration method and the Adomian decomposition method are implemented to give approximate solutions for linear and nonlinear systems of differential equations of fractional order. The two methods in applied mathematics can be used as alternative methods for obtaining analytic and approximate solutions for different types of differential equations. In these schemes, the solution takes the form of a convergent series with easily computable components. This paper presents a numerical comparison between the two methods for solving systems of fractional differential equations. Numerical results show that the two approaches are easy to implement and accurate when applied to differential equations of fractional order.