Lifeng Wang¹, Yunpeng Ma^1,2, Yongqiang Yang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.101, No.2, pp. 97-111, 2014, DOI:10.3970/cmes.2014.101.097

Abstract In this paper, a numerical method based on the Legendre polynomials is presented for a class of variable order fractional differential equation. We adopt the Coimbra variable order fractional operator, which can be viewed as a Caputo-type definition. Three different kinds of operational matrixes with Legendre polynomials are derived. A truncated the Legendre polynomials series together with the products of several dependent matrixes are utilized to reduce the variable order fractional differential equation to a system of algebraic equations. The solution of this system gives the approximation solution for the truncated limited n. An error analysis technique is also given.… More >

Baofeng Li¹

CMES-Computer Modeling in Engineering & Sciences, Vol.93, No.3, pp. 187-202, 2013, DOI:10.3970/cmes.2013.093.187

Abstract In this article, the second kind Chebyshev wavelet method is presented for solving a class of multi-order fractional differential equations (FDEs) with variable coefficients. We first construct the second kind Chebyshev wavelet, prove its convergence and then derive the operational matrix of fractional integration of the second kind Chebyshev wavelet. The operational matrix of fractional integration is utilized to reduce the fractional differential equations to a system of algebraic equations. In addition, illustrative examples are presented to demonstrate the efficiency and accuracy of the proposed method. More >

Zhou YH^1,2, Wang XM², Wang JZ^1,2 , Liu XJ²

CMES-Computer Modeling in Engineering & Sciences, Vol.77, No.2, pp. 137-160, 2011, DOI:10.3970/cmes.2011.077.137

Abstract In this paper, we present an efficient wavelet-based algorithm for solving a class of fractional vibration, diffusion and wave equations with strong nonlinearities. For this purpose, we first suggest a wavelet approximation for a function defined on a bounded interval, in which expansion coefficients are just the function samplings at each nodal point. As the fractional differential equations containing strong nonlinear terms and singular integral kernels, we then use Laplace transform to convert them into the second type Voltera integral equations with non-singular kernels. Certain property of the integral kernel and the ability of explicit wavelet approximation to the nonlinear… More >

Caibin Zeng¹, Qigui Yang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.59, No.1, pp. 65-78, 2010, DOI:10.3970/cmes.2010.059.065

Abstract In this paper, a fractional order model is established to describe HIV internal viral dynamics involving HAART effect. First, the model is proved to possess non-negative solutions as desired in any population dynamics. Then, a detailed analysis is carried out to study the stability of equilibrium points. Numerical simulations are presented to illustrate the stability analysis. More >

Y. T. Gu¹, P. Zhuang², F. Liu³

CMES-Computer Modeling in Engineering & Sciences, Vol.56, No.3, pp. 303-334, 2010, DOI:10.3970/cmes.2010.056.303

Abstract Recently, the numerical modelling and simulation for anomalous subdiffusion equation (ASDE), which is a type of fractional partial differential equation(FPDE) and has been found with widely applications in modern engineering and sciences, are attracting more and more attentions. The current dominant numerical method for modelling ASDE is Finite Difference Method (FDM), which is based on a pre-defined grid leading to inherited issues or shortcomings. This paper aims to develop an implicit meshless approach based on the radial basis functions (RBF) for numerical simulation of the non-linear ASDE. The discrete system of equations is obtained by using the meshless shape functions… More >

Chunhai Kou¹, Ye Yan², Jian Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.39, No.3, pp. 301-318, 2009, DOI:10.3970/cmes.2009.039.301

Abstract In the paper, stability for fractional order differential equations is studied. Then the result obtained is applied to analyse the stability of equilibrium for the model of HIV. More >