No Arabic abstract
In this paper, we propose a new approach, based on the so-called modulating functions to estimate the average velocity, the dispersion coefficient and the differentiation order in a space fractional advection dispersion equation. First, the average velocity and the dispersion coefficient are estimated by applying the modulating functions method, where the problem is transferred into solving a system of algebraic equations. Then, the modulating functions method combined with Newtons method is applied to estimate all three parameters simultaneously. Numerical results are presented with noisy measurements to show the effectiveness and the robustness of the proposed method.
In this paper, we give some new thoughts about the classical gradient method (GM) and recall the proposed fractional order gradient method (FOGM). It is proven that the proposed FOGM holds a super convergence capacity and a faster convergence rate around the extreme point than the conventional GM. The property of asymptotic convergence of conventional GM and FOGM is also discussed. To achieve both a super convergence capability and an even faster convergence rate, a novel switching FOGM is proposed. Moreover, we extend the obtained conclusion to a more general case by introducing the concept of p-order Lipschitz continuous gradient and p-order strong convex. Numerous simulation examples are provided to validate the effectiveness of proposed methods.
This contribution deals with identification of fractional-order dynamical systems. System identification, which refers to estimation of process parameters, is a necessity in control theory. Real processes are usually of fractional order as opposed to the ideal integral order models. A simple and elegant scheme of estimating the parameters for such a fractional order process is proposed. This method employs fractional calculus theory to find equations relating the parameters that are to be estimated, and then estimates the process parameters after solving the simultaneous equations. The said simultaneous equations are generated and updated using particle swarm optimization (PSO) technique, the fitness function being the sum of squared deviations from the actual set of observations. The data used for the calculations are intentionally corrupted to simulate real-life conditions. Results show that the proposed scheme offers a very high degree of accuracy even for erroneous data.
Nonlocal operators of fractional type are a popular modeling choice for applications that do not adhere to classical diffusive behavior; however, one major challenge in nonlocal simulations is the selection of model parameters. In this work we propose an optimization-based approach to parameter identification for fractional models with an optional truncation radius. We formulate the inference problem as an optimal control problem where the objective is to minimize the discrepancy between observed data and an approximate solution of the model, and the control variables are the fractional order and the truncation length. For the numerical solution of the minimization problem we propose a gradient-based approach, where we enhance the numerical performance by an approximation of the bilinear form of the state equation and its derivative with respect to the fractional order. Several numerical tests in one and two dimensions illustrate the theoretical results and show the robustness and applicability of our method.
We consider the inverse problem of estimating the spatially varying pulse wave velocity in blood vessels in the brain from dynamic MRI data, as it appears in the recently proposed imaging technique of Magnetic Resonance Advection Imaging (MRAI). The problem is formulated as a linear operator equation with a noisy operator and solved using a conjugate gradient type approach. Numerical examples on experimental data show the usefulness and advantages of the developed algorithm in comparison to previously proposed methods.
This article studies a direct numerical approach for fractional advection-diffusion equations (ADEs). Using a set of cubic trigonometric B-splines as test functions, a differential quadrature (DQ) method is firstly proposed for the 1D and 2D time-fractional ADEs of order $(0,1]$. The weighted coefficients are determined, and with them, the original equation is transformed into a group of general ordinary differential equations (ODEs), which are discretized by an effective difference scheme or Runge-Kutta method. The stability is investigated under a mild theoretical condition. Secondly, based on a set of cubic B-splines, we develop a new Crank-Nicolson type DQ method for the 2D space-fractional ADEs without advection. The DQ approximations to fractional derivatives are introduced and the values of the fractional derivatives of B-splines are computed by deriving explicit formulas. The presented DQ methods are evaluated on five benchmark problems and the concrete simulations of the unsteady propagation of solitons and Gaussian pulse. In comparison with the existing algorithms in the open literature, numerical results finally illustrate the validity and accuracy.