On higher-order semi-explicit symplectic partitioned Runge-Kutta methods for constrained Hamiltonian systems

Summary. In this paper we generalize the class of explicit partitioned Runge-Kutta (PRK) methods for separable Hamiltonian systems to systems with holonomic constraints. For a convenient analysis of such schemes, we first generalize the backward error analysis for systems in to systems on manifolds embedded in . By applying this analysis to constrained PRK methods, we prove that such methods will, in general, suffer from order reduction as well-known for higher-index differential-algebraic equations. However, this order reduction can be avoided by a proper modification of the standard PRK methods. This modification increases the number of projection steps onto the constraint manifold but leaves the number of force evaluations constant. We also give a numerical comparison of several second, fourth, and sixth order methods. Received May 5, 1995 / Revised version received February 7, 1996     

Hamilton系统的连续有限元法

利用常微分方程的连续有限元法,对非线性Hamilton系统证明了连续一次、二次有限元法分别是2阶和3阶的拟辛格式,且保持能量守恒;连续有限元法是辛算法对线性Hamilton系统,且保持能量守恒．在数值计算上探讨了辛性质和能量守恒性,与已有的辛算法进行对比,结果与理论相吻合．     

Fault detection based on fractional order models: Application to diagnosis of thermal systems

The aim of this paper is to propose diagnosis methods based on fractional order models and to validate their efficiency to detect faults occurring in thermal systems. Indeed, it is first shown that fractional operator allows to derive in a straightforward way fractional models for thermal phenomena. In order to apply classical diagnosis methods, such models could be approximated by integer order models, but at the expense of much higher involved parameters and reduced precision. Thus, two diagnosis methods initially developed for integer order models are here extended to handle fractional order models. The first one is the generalized dynamic parity space method and the second one is the Luenberger diagnosis observer. Proposed methods are then applied to a single-input multi-output thermal testing bench and demonstrate the methods efficiency for detecting faults affecting thermal systems.     

Half-explicit Runge-Kutta methods with explicit stages for differential-algebraic systems of index 2

Usually the straightforward generalization of explicit Runge-Kutta methods for ordinary differential equations to half-explicit methods for differential-algebraic systems of index 2 results in methods of orderq≤2. The construction of higher order methods is simplified substantially by a slight modification of the method combined with an improved strategy for the computation of the algebraic solution components. We give order conditions up to orderq=5 and study the convergence of these methods. Based on the fifth order method of Dormand and Prince the fifth order half-explicit Runge-Kutta method HEDOP5 is constructed that requires the solution of 6 systems of nonlinear equations per step of integration.     

An indirect Lyapunov approach to robust stabilization for a class of linear fractional-order system with positive real uncertainty

The paper is concerned with the problem of the robust stabilization for a class of fractional order linear systems with positive real uncertainty under Riemann–Liouville (RL) derivatives. Firstly, by utilizing the continuous frequency distributed model of the fractional integrator, the fractional order system is expressed as an infinite dimensional integral order system. And via using indirect Lyapunov approach and linear matrix inequality techniques, sufficient condition for robust asymptotic stability of the fractional order systems and design methods of the state feedback controller are presented. Secondly, by using matrixs singular value decomposition technique the static output feedback controller and observer-based controller for asymptotically stabilizing the fractional order systems are derived. Finally, the validity of the proposed methods are demonstrated by numerical examples.     

High order finite volume schemes based on reconstruction of states for solving hyperbolic systems with nonconservative products. Applications to shallow-water systems

This paper is concerned with the development of high order methods for the numerical approximation of one-dimensional nonconservative hyperbolic systems. In particular, we are interested in high order extensions of the generalized Roe methods introduced by I. Toumi in 1992, based on WENO reconstruction of states. We also investigate the well-balanced properties of the resulting schemes. Finally, we will focus on applications to shallow-water systems.

     

A technique to construct symmetric variable-stepsize linear multistep methods for second-order systems

Some previous works show that symmetric fixed- and variable-stepsize linear multistep methods for second-order systems which do not have any parasitic root in their first characteristic polynomial give rise to a slow error growth with time when integrating reversible systems. In this paper, we give a technique to construct variable-stepsize symmetric methods from their fixed-stepsize counterparts, in such a way that the former have the same order as the latter. The order and symmetry of the integrators obtained is proved independently of the order of the underlying fixed-stepsize integrators. As this technique looks for efficiency, we concentrate on explicit linear multistep methods, which just make one function evaluation per step, and we offer some numerical comparisons with other one-step adaptive methods which also show a good long-term behaviour.

     

Arbitrary-order functionally fitted energy-diminishing methods for gradient systems

It is well known that for gradient systems in Euclidean space or on a Riemannian manifold, the energy decreases monotonically along solutions. In this letter we derive and analyse functionally fitted energy-diminishing methods to preserve this key property of gradient systems. It is proved that the novel methods are energy-diminishing and can achieve damping for very stiff gradient systems. We also show that the methods can be of arbitrarily high order and discuss their implementations. A numerical test is reported to illustrate the efficiency of the new methods in comparison with three existing numerical methods in the literature.     

Two-step numerical methods for parabolic differential equations

A family of two-stepA-stable methods of maximal order for the numerical solution of ordinary differential systems is developed. If these methods are applied to the stiff, large systems which originate from linear parabolic differential equations they yield a large, sparse set of linear algebraic equations of special form. This set is considerably easier to solve than the algebraic equations which are obtained when using diagonal Obrechkoff methods, which are one-step,A-stable and of maximal order     

A family of Steffensen type methods with seventh-order convergence

In this paper, based on some known fourth-order Steffensen type methods, we present a family of three-step seventh-order Steffensen type iterative methods for solving nonlinear equations and nonlinear systems. For nonlinear systems, a development of the inverse first-order divided difference operator for multivariable function is applied to prove the order of convergence of the new methods. Numerical experiments with comparison to some existing methods are provided to support the underlying theory.     

Switching adaptive controllers to control fractional‐order complex systems with unknown structure and input nonlinearities

This article investigates the chaos control problem for the fractional‐order chaotic systems containing unknown structure and input nonlinearities. Two types of nonlinearity in the control input are considered. In the first case, a general continuous nonlinearity input is supposed in the controller, and in the second case, the unknown dead‐zone input is included. In each case, a proper switching adaptive controller is introduced to stabilize the fractional‐order chaotic system in the presence of unknown parameters and uncertainties. The control methods are designed based on the boundedness property of the chaotic system's states, where, in the proposed methods the nonlinear/linear dynamic terms of the fractional‐order chaotic systems are assumed to be fully unknown. The analytical results of the mentioned techniques are proved by the stability analysis theorem of fractional‐order systems and the adaptive control method. In addition, as an application of the proposed methods, single input adaptive controllers are adopted for control of a class of three‐dimensional nonlinear fractional‐order chaotic systems. And finally, some numerical examples illustrate the correctness of the analytical results. © 2014 Wiley Periodicals, Inc. Complexity 21: 211–223, 2015     

波动方程的一类显式辛格式

1．引言和预备知识本文主要考虑如下波动方程初边值问题的数值方法．一般的哈密顿系统可写成那么，系统（1．幻称为可分的哈密顿系统．众所周知，方程（1．la）在引进新变量。t＝v后，它变成一类可分的线性哈密顿系统：它的哈密顿函数为H一三矿（V“＋。乏）咖，并且在离散＊。＝。。。后，（1．4）是一个方程组，u，V是向量．人们早就知道，初边值问题（1．1）在应用隐式中点公式（辛格式）进行数值解时，能保持真解许多重要性质，并且格式是无条件稳定的，但美中不足的是，在每积分步，它要求解一个线性方程组，当考虑大的离散系统时，…     

Exploiting structure in the construction of DIMSIMs

We describe the structure of the nonlinear system for the coefficients of diagonally implicit multistage integration methods for ordinary differential equations. This structure is then utilized in the search for methods of high order and stage order with a given stability function. New methods were obtained by solving the nonlinear systems by state-of-the-art software based on least squares minimization.     

On Stability of Symplectic Algorithms

1.FundamentalDeflnitionsLemma1.Thesolutionofalinearoofinarydtherentialequationwithcon8tantcoeffcientY=AYissta6leifalleigenvalue8ofAhaven0nP6sitivercalpartsandtheeigenvalueswithnullrealpartaresingleroots0ftheminimalp0lynomial.,/P\ThelinearHamiltoniansystemcanbeden0tedasZ=JSZwhereZ=(q),J=(ELs),andtheHamiltonianfuncti0nH(z)=ty.Lemma2.Thesolution80flinearHamiltoniansy8temsarecmticallysta6leifalleigenvaluesofJShavenullrsalpartandaresinglerootsojtheminitnalp0lyno?nial.Definiti0n1.Whenthemo…     

PARALLEL ALGORITHMS OF ONE-LEG METHOD AND ITERATED DEFECT CORRECTION METHODS

The parallel algorithms of iterated defect correction methods (PIDeCM's) are constructed, which are of efficiency and high order B-convergence for general nonlinear stiff systems in ODE'S. As the basis of constructing and discussing PIDeCM's. a class of parallel one-leg methods is also investigated, which are of particular efficiency for linear systems.     

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.     

Efficient high-order methods based on golden ratio for nonlinear systems

We derive new iterative methods with order of convergence four or higher, for solving nonlinear systems, by composing iteratively golden ratio methods with a modified Newton’s method. We use different efficiency indices in order to compare the new methods with other ones and present several numerical tests which confirm the theoretical results.     

An improved tri-coloured rooted-tree theory and order conditions for ERKN methods for general multi-frequency oscillatory systems

This paper develops an improved tri-coloured rooted-tree theory for the order conditions for ERKN methods solving general multi-frequency and multidimensional second-order oscillatory systems. The bottleneck of the original tricoloured rooted-tree theory is the existence of numerous redundant trees. In light of the fact that the sum of the products of the symmetries and the elementary differentials is meaningful, this paper naturally introduces the so-called extended elementary differential mappings. Then, the new improved tri-coloured rooted tree theory is established based on a subset of the original tri-coloured rooted-tree set. This new theory makes all redundant trees disappear, and thus, the order conditions of ERKN methods for general multi-frequency and multidimensional second-order oscillatory systems are reduced greatly. Furthermore, with this new theory, we present some new ERKN methods of order up to four. Numerical experiments are implemented and the results show that ERKN methods can be competitive with other existing methods in the scientific literature, especially when comparatively large stepsizes are used.     

On projected Runge-Kutta methods for differential-algebraic equations

Ascher and Petzold recently introducedprojected Runge-Kutta methods for the numerical solution of semi-explicit differential-algebraic systems of index 2. Here it is shown that such a method can be regarded as the limiting case of a standard application of a Runge-Kutta method with a very small implicit Euler step added to it. This interpretation allows a direct derivation of the order conditions and of superconvergence results for the projected methods from known results for standard Runge-Kutta methods for index-2 differential-algebraic systems, and an extension to linearly implicit differential-algebraic systems.     

Implementation of dynamic programming for chaos control in discrete systems

In this paper the control of discrete chaotic systems by designing linear feedback controllers is presented. The linear feedback control problem for nonlinear systems has been formulated under the viewpoint of dynamic programming. For suppressing chaos with minimum control effort, the system is stabilized on its first order unstable fixed point (UFP). The presented method also could be employed to make any desired nth order fixed point of the system, stable. Two different methods for higher order UFPs stabilization are suggested. Afterwards, these methods are applied to two well-known chaotic discrete systems: the Logistic and the Henon Maps. For each of them, the first and second UFPs in their chaotic regions are stabilized and simulation results are provided for the demonstration of performance.