Chebyshev operational matrix method for solving multi-order fractional ordinary differential equations

The aim of this article is to present an analytical approximation solution for linear and nonlinear multi-order fractional differential equations (FDEs) by extending the application of the shifted Chebyshev operational matrix. For this purpose, we convert FDE into a counterpart system and then using proposed method to solve the resultant system. Our results in solving four different linear and nonlinear FDE, confirm the accuracy of proposed method.     

On a boundary value problem for delay differential equations of population dynamics and Chebyshev approximation

In this paper, a boundary value problem for delay differential equations of population dynamics is considered. We obtain approximate solutions by using Chebyshev polynomial series and Newton–Raphson's procedure and give the error estimation. The method of the error estimation has been obtained in an existence theorem proved by a part of the authors. We carry out some numerical experiments by a computer language MATLAB.     

The use of Chebyshev cardinal functions for the solution of a partial differential equation with an unknown time-dependent coefficient subject to an extra measurement

A numerical technique is presented for the solution of a parabolic partial differential equation with a time-dependent coefficient subject to an extra measurement. The method is derived by expanding the required approximate solution as the elements of Chebyshev cardinal functions. Using the operational matrix of derivative, the problem can be reduced to a set of algebraic equations. From the computational point of view, the solution obtained by this method is in excellent agreement with those obtained by previous works and also it is efficient to use.     

Fractional-order Legendre functions for solving fractional-order differential equations

In this article, a general formulation for the fractional-order Legendre functions (FLFs) is constructed to obtain the solution of the fractional-order differential equations. Fractional calculus has been used to model physical and engineering processes that are found to be best described by fractional differential equations. Therefore, an efficient and reliable technique for the solution of them is too important. For the concept of fractional derivative we will adopt Caputo’s definition by using Riemann–Liouville fractional integral operator. Our main aim is to generalize the new orthogonal functions based on Legendre polynomials to the fractional calculus. Also a general formulation for FLFs fractional derivatives and product operational matrices is driven. These matrices together with the Tau method are then utilized to reduce the solution of this problem to the solution of a system of algebraic equations. The method is applied to solve linear and nonlinear fractional differential equations. Illustrative examples are included to demonstrate the validity and applicability of the presented technique.     

The use of Chebyshev cardinal functions for solution of the second‐order one‐dimensional telegraph equation

A numerical technique is presented for the solution of the second order one‐dimensional linear hyperbolic equation. This method uses the Chebyshev cardinal functions. The method consists of expanding the required approximate solution as the elements of Chebyshev cardinal functions. Using the operational matrix of derivative, the problem is reduced to a set of algebraic equations. Some numerical examples are included to demonstrate the validity and applicability of the technique. The method is easy to implement and produces very accurate results. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

Numerical solution of the system of nonlinear Volterra integro-differential equations with nonlinear differential part by the operational Tau method and error estimation

The operational Tau method, a well-known method for solving functional equations, is employed to solve a system of nonlinear Volterra integro-differential equations with nonlinear differential part. In addition, an error estimation of the method is presented. Some cases of the mentioned equations are solved as examples to illustrate the ability and reliability of the method. The results reveal that the method is very effective and convenient.     

A new approach of the Chebyshev wavelets method for partial differential equations with boundary conditions of the telegraph type

In this paper, we develop an accurate and efficient Chebyshev wavelets method for solution of partial differential equations with boundary conditions of the telegraph type. In the proposed method we have employed mutually the operational matrices of integration and differentiation to get numerical solutions of such equations. The power of this manageable method is confirmed. Moreover the use of Chebyshev wavelet is found to be accurate, simple and fast.     

Wavelets method for solving systems of nonlinear singular fractional Volterra integro-differential equations

This paper presents a computational method for solving a class of system of nonlinear singular fractional Volterra integro-differential equations. First, existences of a unique solution for under studying problem is proved. Then, shifted Chebyshev polynomials and their properties are employed to derive a general procedure for forming the operational matrix of fractional derivative for Chebyshev wavelets. The application of this operational matrix for solving mentioned problem is explained. In the next step, the error analysis of the proposed method is investigated. Finally, some examples are included for demonstrating the efficiency of the proposed method.     

Application of fuzzy differential transform method for solving fuzzy Volterra integral equations

This paper deals with the solutions of fuzzy Volterra integral equations with separable kernel by using fuzzy differential transform method (FDTM). If the equation considered has a solution in terms of the series expansion of known functions, this powerful method catches the exact solution. To this end, we have obtained several new results to solve mentioned problem when FDTM has been applied. In order to show this capability and robustness, some fuzzy Volterra integral equations are solved in detail as numerical examples.     

Spectral methods for pantograph-type differential and integral equations with multiple delays

We analyze the convergence properties of the spectral method when used to approximate smooth solutions of delay differential or integral equations with two or more vanishing delays. It is shown that for the pantograph-type functional equations the spectral methods yield the familiar exponential order of convergence. Various numerical examples are used to illustrate these results.      

New predictor-corrector approach for nonlinear fractional differential equations: error analysis and stability

In this paper, the predictor-corrector approach is used to propose two algorithms for the numerical solution of linear and non-linear fractional differential equations (FDE). The fractional order derivative is taken to be in the sense of Caputo and its properties are used to transform FDE into a Volterra-type integral equation. Simpson''s 3/8 rule is used to develop new numerical schemes to obtain the approximate solution of the integral equation associated with the given FDE. The error and stability analysis for the two methods are presented. The proposed methods are compared with the ones available in the literature. Numerical simulation is performed to demonstrate the validity and applicability of both the proposed techniques. As an application, the problem of dynamics of the new fractional order non-linear chaotic system introduced by Bhalekar and Daftardar-Gejji is investigated by means of the obtained numerical algorithms.     

Numerical solution of Volterra integral and integro-differential equations of convolution type by using operational matrices of piecewise constant orthogonal functions

In this paper, we use operational matrices of piecewise constant orthogonal functions on the interval [0,1)$[0,1)$ to solve Volterra integral and integro-differential equations of convolution type without solving any system. We first obtain Laplace transform of the problem and then we find numerical inversion of Laplace transform by operational matrices. Numerical examples show that the approximate solutions have a good degree of accuracy.     

Method of Lyapunov functions for differential equations with piecewise constant delay

We address differential equations with piecewise constant argument of generalized type [5], [6], [7] and [8] and investigate their stability with the second Lyapunov method. Despite the fact that these equations include delay, stability conditions are merely given in terms of Lyapunov functions; that is, no functionals are used. Several examples, one of which considers the logistic equation, are discussed to illustrate the development of the theory. Some of the results were announced at the 14th International Congress on Computational and Applied Mathematics (ICCAM2009), Antalya, Turkey, in 2009.     

On the convergence of spline collocation methods for solving fractional differential equations

In the first part of this paper we study the regularity properties of solutions of initial value problems of linear multi-term fractional differential equations. We then use these results in the convergence analysis of a polynomial spline collocation method for solving such problems numerically. Using an integral equation reformulation and special non-uniform grids, global convergence estimates are derived. From these estimates it follows that the method has a rapid convergence if we use suitable nonuniform grids and the nodes of the composite Gaussian quadrature formulas as collocation points. Theoretical results are verified by some numerical examples.     

On collocation methods for delay differential and Volterra integral equations with proportional delay

To compute long term integrations for the pantograph differential equation with proportional delay qt, 0 < q ⩽ 1: y′(t) = ay(t) + by(qt) + f(t), y(0) = y ₀, we offer two kinds of numerical methods using special mesh distributions, that is, a rational approximant with ‘quasi-uniform meshes’ (see E. Ishiwata and Y. Muroya [Appl. Math. Comput., 2007, 187: 741-747]) and a Gauss collocation method with ‘quasi-constrained meshes’. If we apply these meshes to rational approximant and Gauss collocation method, respectively, then we obtain useful numerical methods of order p ^* = 2m for computing long term integrations. Numerical investigations for these methods are also presented.      

线性常微分方程的全局误差估计和优化求解方法

线性常微分方程初值问题求解在许多应用中起着重要作用.目前,已存在很多的数值方法和求解器用于计算离散网格点上的近似解,但很少有对全局误差(global error)进行估计和优化的方法.本文首先通过将离散数值解插值成为可微函数用来定义方程的残差;再给出残差与近似解的关系定理并推导出全局误差的上界;然后以最小化残差的二范数为目标将方程求解问题转化为优化求解问题;最后通过分析导出矩阵的结构,提出利用共轭梯度法对其进行求解.之后将该方法应用于滤波电路和汽车悬架系统等实际问题.实验分析表明,本文估计方法对线性常微分方程的初值问题的全局误差具有比较好的估计效果,优化求解方法能够在不增加网格点的情形下求解出线性常微分方程在插值解空间中的全局最优解.     

Approximation, integration and differentiation of time functions using a set of orthogonal hybrid functions (HF) and their application to solution of first order differential equations

Differential equations of different types and orders are of utmost importance for mathematical modeling of control system problems. State variable method uses the concept of expressing n number of first order differential equations in vector matrix form to model and analyze/synthesize control systems.The present work proposes a new set of orthogonal hybrid functions (HF) which evolved from synthesis of sample-and-hold functions (SHF) and triangular functions (TF). This HF set is used to approximate a time function in a piecewise linear manner with the mean integral square error (MISE) much less than block pulse function based approximation which always provides staircase solutions.The operational matrices for integration and differentiation in HF domain are also derived and employed for solving non-homogeneous and homogeneous differential equations of the first order as well as state equations. The results are compared with exact solutions, the 4th order Runge-Kutta method and its further improved versions proposed by Simos [6]. The presented HF domain theory is well supported by a few illustrations.     

B-splines methods with redefined basis functions for solving fourth order parabolic partial differential equations

In this work, we discuss two methods for solving a fourth order parabolic partial differential equation. In Method-I, we decompose the given equation into a system of second order equations and solve them by using cubic B-spline method with redefined basis functions. In Method-II, the equation is solved directly by applying quintic B-spline method with redefined basis functions. Stability of these methods have been discussed. Both methods are unconditionally stable. These methods are tested on four examples. The computed results are compared wherever possible with those already available in literature. We have developed Method-I for fourth order non homogeneous parabolic partial differential equation from which we can obtain displacement and bending moment both simultaneously, while Method-II gives only displacement. The results show that the derived methods are easily implemented and approximate the exact solution very well.     

Oscillation criteria of a certain class of third order nonlinear delay differential equations with damping

In this paper, we are concerned with the oscillation of third order nonlinear delay differential equations of the form

(r₂(t)(r₁(t)y^′)′)′+p(t)y^′+q(t)f(y(g(t)))=0.