Charpit System of Partial Differential Equations with Algebraic Constraints

In this paper the Charpit system of partial differential equations with algebraic constraints is considered. So, first the compatibility conditions of a system of algebraic equations and also of the Charpit system of partial differential equations are separately considered. For the combined system of equations of both types sufficient conditions for the existence of a solution are found. They lead to an algorithm for reducing the combined system to a Charpit system of partial differential equations of dimension less than the initial system and without algebraic constraints. Moreover, it is proved that this system identically satisfies the compatibility conditions if so does the initial system.     

New Multivalue Methods for Differential Algebraic Equations

Multivalue methods are slightly different from the general linear methods John Butcher proposed over 30 years ago. Multivalue methods capable of solving differential algebraic equations have not been developed. In this paper, we have constructed three new multivalue methods for solving DAEs of index 1, 2 or 3, which include multistep methods and multistage methods as special cases. The concept of stiff accuracy will be introduced and convergence results will be given based on the stage order of the methods. These new methods have the diagonal implicit property and thus are cheap to implement and will have order 2 or more for both the differential and algebraic components. We have implemented these methods with fixed step size and they are shown to be very successful on a variety of problems. Some numerical experiments with these methods are presented.     

求解延迟微分代数方程的多步Runge-Kutta方法的渐近稳定性

延迟微分代数方程(DDAEs)广泛出现于科学与工程应用领域．本文将多步Runge-Kutta方法应用于求解线性常系数延迟微分代数方程，讨论了该方法的渐近稳定性．数值试验表明该方法对求解DDAEs是有效的．     

一类随机延迟微分代数系统的Euler-Maruyama方法

我们主要构造了数值求解一类1指标随机延迟微分代数系统的Euler-Maruyama方法,并且证明用该方法求解此类问题可达到1/2阶均方收敛.最后的效值试验验证了方法的有效性及所获结论的正确性.     

关于复代数微分方程组解的增长级问题(英文)

This paper is concerned with the order of the solutions of systems of high-order complex algebraic differential equations.By means of Zalcman Lemma,the systems of equations of[1]is extended to more general form.     

Boundary Value Problems for Systems of Functional Differential Equations

Algorithms for finding an approximate solution of boundary value problems for systems of functional ordinary differential equations are studied. Sufficient conditions for consistency and convergence of these methods are given. In the last section, a construction of methods of arbitrary order is presented.     

Optimal Control of One-Dimensional Partial Differential Algebraic Equations with Applications

We present an approach to compute optimal control functions in dynamic models based on one-dimensional partial differential algebraic equations (PDAE). By using the method of lines, the PDAE is transformed into a large system of usually stiff ordinary differential algebraic equations and integrated by standard methods. The resulting nonlinear programming problem is solved by the sequential quadratic programming code NLPQL. Optimal control functions are approximated by piecewise constant, piecewise linear or bang-bang functions. Three different types of cost functions can be formulated. The underlying model structure is quite flexible. We allow break points for model changes, disjoint integration areas with respect to spatial variable, arbitrary boundary and transition conditions, coupled ordinary and algebraic differential equations, algebraic equations in time and space variables, and dynamic constraints for control and state variables. The PDAE is discretized by difference formulae, polynomial approximations with arbitrary degrees, and by special update formulae in case of hyperbolic equations. Two application problems are outlined in detail. We present a model for optimal control of transdermal diffusion of drugs, where the diffusion speed is controlled by an electric field, and a model for the optimal control of the input feed of an acetylene reactor given in form of a distributed parameter system.     

General Linear Methods for Stiff Differential Equations

A new type of general linear method is constructed which combines A-stability or L-stability with ease of implementation. The method is structured in such a manner that its stability region is identical with that of a Runge-Kutta method, using a restriction known as inherent RK stability.This revised version was published online in October 2005 with corrections to the Cover Date.     

解微分方程组的改进尤拉方法的改进

对改进尤拉方法解微分方程组的方法作了改进,改进的算法与原来算法的计算量一样,但精度比较高.     

Convergence of Runge-Kutta Methods for Delay Differential Equations

This paper deals with the adaptation of Runge—Kutta methods to the numerical solution of nonstiff initial value problems for delay differential equations. We consider the interpolation procedure that was proposed in In 't Hout [8], and prove the new and positive result that for any given Runge—Kutta method its adaptation to delay differential equations by means of this interpolation procedure has an order of convergence equal to min {p,q}, where p denotes the order of consistency of the Runge—Kutta method and q is the number of support points of the interpolation procedure.This revised version was published online in October 2005 with corrections to the Cover Date.     

On the Global Error of Discretization Methods for Highly-Oscillatory Ordinary Differential Equations

Commencing from a global-error formula, originally due to Henrici, we investigate the accumulation of global error in the numerical solution of linear highly-oscillating systems of the form y+ g(t)y = 0, where g(t) . Using WKB analysis we derive an explicit form of the global-error envelope for Runge-Kutta and Magnus methods. Our results are closely matched by numerical experiments.Motivated by the superior performance of Lie-group methods, we present a modification of the Magnus expansion which displays even better long-term behaviour in the presence of oscillations.     

Multi-Step Maruyama Methods for Stochastic Delay Differential Equations

Abstract

In this article the numerical approximation of solutions of Itô stochastic delay differential equations is considered. We construct stochastic linear multi-step Maruyama methods and develop the fundamental numerical analysis concerning their 𝕃 _p -consistency, numerical 𝕃 _p -stability and 𝕃 _p -convergence. For the special case of two-step Maruyama schemes we derive conditions guaranteeing their mean-square consistency.     

线性方程组的异步松弛迭代法*

本文考虑解线性方程组经典迭代法的异步形式,对系数矩阵为H矩阵,给出了异步迭代过程收敛性的充分条件,这不仅降低了文献[3]对系数矩阵的要求,而且收敛区域比文献[3]的大.     

具非线性扩散系数的中立双曲型偏微分方程系统解的振动定理

考虑一类具非线性扩散系数的中立双曲型偏微分方程系统解的振动性.利用Green公式和Riccati变换,获得了该方程组在两类不同边值条件下振动的若干充分条件.     

Nonlinear Stability of General Linear Methods for Delay Differential Equations

This paper is devoted to a study of nonlinear stability of general linear methods for the numerical solution of delay differential equations in Hilbert spaces. New stability concepts are further introduced. The stability properties of (k,p,q)-algebraically stable general linear methods with piecewise constant or linear interpolation procedure are investigated. We also discuss stability of linear multistep methods viewed as a special subset of the class of general linear methods.     

多体系统动力学微分／代数方程组的一类新的数值分析方法

本文讨论了多体系统动力学微分／代数混合方程组的数值离散问题．首先把参数t并入广义坐标讨论，简化了方程组及其隐含条件的结构，并将其化为指标1的方程组．然后利用方程组的特殊结构，引入一种局部离散技巧并构造了相应的算法．算法结构紧凑，易于编程，具有较高的计算效率和良好的数值性态，且其形式适合于各种数值积分方法的的实施．文末给出了具体算例．     

Nash Equilibria for the Multiobjective Control of Linear Partial Differential Equations

This article is concerned with the numerical solution of multiobjective control problems associated with linear partial differential equations. More precisely, for such problems, we look for the Nash equilibrium, which is the solution to a noncooperative game. First, we study the continuous case. Then, to compute the solution of the problem, we combine finite-difference methods for the time discretization, finite-element methods for the space discretization, and conjugate-gradient algorithms for the iterative solution of the discrete control problems. Finally, we apply the above methodology to the solution of several tests problems.     

复高阶代数微分方程组的亚纯解

利用亚纯函数的Nevanlinna值分布理论,研究了一类高阶代数微分方程组亚纯允许解的存在性问题,得到了一个主要结果.     

多时滞微分方程数值稳定性

考虑了时滞微分方程的初值问题,分析了用线性多步法求解一类滞后型微分系统数值解的稳定性,在一定的Lagrange插值条件下,给出并证明了求解滞后型微分系统的线性多步法数值稳定的充分必要条件.     

亚纯系数的非齐次线性微分方程的振荡解

本文研究关于亚纯系数的非齐次线性微分方程的复振荡,得到方程f(k)+ak-1fk-1+…+a0f=F(a0,a1,…,ak-1和F是亚纯函数)具有一个振荡解空间,其空间中所有解的零点收敛指数为∞,至多除去一个例外值.