Numerical methods for Hamilton Jacobi functional differential equations

Initial and initial boundary value problems for first order partial functional differential equations are considered. Explicit difference schemes of the Euler type and implicit difference methods are investigated. The following theoretical aspects of the methods are presented. Sufficient conditions for the convergence of approximate solutions are given and comparisons of the methods are presented. It is proved that assumptions on the regularity of given functions are the same for both the methods. It is shown that conditions on the mesh for explicit difference schemes are more restrictive than suitable assumptions for implicit methods. There are implicit difference schemes which are convergent and corresponding explicit difference methods are not convergent. Error estimates for both the methods are construted.     

一般线性方法的线性与非线性稳定性间的关系

１．引言１９６３年,Ｄａｈｌｑｕｉｓｔ以一类线性问题为模型提出了Ａ－稳定性概念,此后有关如何判断方法是否Ａ－稳定或确定其稳定域的研究十分活跃,袁兆鼎等~［１］中对此有详细的讨论．Ｂｕｒｒａｇｅ与Ｂｕｔｃｈｅｒ［２］以一类非线性问题为模型,就一般线性方法引入了代数稳定性概念．Ｂｕｔｃｈｅｒ［４］探讨了代数稳定性与Ａ－稳定性间的内在联系．为确保代数稳定性蕴涵Ａ－稳定性,Ｂｕｔｃｈｅｒ［５］进一步要求代数稳定性定义中涉及的矩阵Ｇ是正定的。然而这样一来,正如李寿佛［６］中指出的那样,许多ＡＮ稳定且按［４…     

Expanding stability regions of explicit advective-diffusive finite difference methods by Jacobi preconditioning

Explicit time differencing methods for solving differential equations are advantageous in that they are easy to implement on a computer and are intrinsically very parallel. The disadvantage of explicit methods is the severe restrictions that are placed on stable time-step intervals. Stability bounds for explicit time differencing methods on advective–diffusive problems are generally determined by the diffusive part of the problem. These bounds are very small and implicit methods are used instead. The linear systems arising from these implicit methods are generally solved by iterative methods. In this article we develop a methodology for increasing the stability bounds of standard explicit finite differencing methods by combining explicit methods, implicit methods, and iterative methods in a novel way to generate new time-difference schemes, called preconditioned time-difference methods. A Jacobi preconditioned time differencing method is defined and analyzed for both diffusion and advection–diffusion equations. Several computational examples of both linear and nonlinear advective-diffusive problems are solved to demonstrate the accuracy and improved stability limits. © 1995 John Wiley & Sons, Inc.     

时滞积分微分方程的Rosenbrock方法

本文研究时滞积分微分方程的数值方法.通过改造现有常及离散型延迟微分方程的数值方法,并匹配以适当数值求积公式,构造了求解时滞积分微分方程的Rosenbrock方法,导出了其稳定性准则.数值例子阐明了所获方法的计算有效性.     

Fourth-Order Splitting Methods for Time-Dependant Differential Equations

This study was suggested by previous work on the simulation of evolution equations with scale-dependent processes,e.g.,wave-propagation or heat-transfer,that are modeled by wave equations or heat equations.Here,we study both parabolic and hyperbolic equations.We focus on ADI (alternating direction implicit) methods and LOD (locally one-dimensional) methods,which are standard splitting methods of lower order,e.g.second-order.Our aim is to develop higher-order ADI methods,which are performed by Richardson extrapolation,Crank-Nicolson methods and higher-order LOD methods,based on locally higher-order methods.We discuss the new theoretical results of the stability and consistency of the ADI methods.The main idea is to apply a higher- order time discretization and combine it with the ADI methods.We also discuss the dis- cretization and splitting methods for first-order and second-order evolution equations. The stability analysis is given for the ADI method for first-order time derivatives and for the LOD (locally one-dimensional) methods for second-order time derivatives.The higher-order methods are unconditionally stable.Some numerical experiments verify our results.     

多组变量典型相关分析SSQCOR准则的数值求解

本文研究多组变量相关分析SSQCOR.准则的数值方法.从KKT条件出发引入了Gauss-Seidel型方法,从SSQCOR出发引入了交替变量法.证明了前者是后者的非精确形式,都具有单调上升性.为了提高得到全局解的可能性,引入了初始点策略.用实际数据和模拟数据进行了数值试验以说明算法的有效性.     

Using error-bounds for hyperpower methods to calculate inclusions for the inverse of a matrix

By means of error-bounds for the well-known hyperpower methods for approximating the inverse of a matrix we define inclusion methods for the inverse matrix. These methods are using machine interval operations and are giving guaranteed inclusions for the inverse matrix whenever the convergence of the applied hyperpower methods can be shown. In comparison with the very efficient interval Schulz's method in the literature, our methods are more efficient in terms of the efficiency index. Some numerical examples are given.     

A comparison of statistical methods for prenatal screening for Down syndrome

Currently, prenatal screening for Down Syndrome (DS) uses the mother's age as well as three biochemical markers for risk prediction. Risk calculations for the biochemical markers use a quadratic discriminant function. In this paper we compare several classification procedures to quadratic discrimination methods for biochemical-based DS risk prediction, based on data from a prospective multicentre prenatal screening study. We investigate alternative methods including linear discriminant methods, logistic regression methods, neural network methods, and classification and regression-tree methods. Several experiments are performed, and in each experiment resampling methods are used to create training and testing data sets. The procedures on the test data set are summarized by the area under their receiver operating characteristic curves. In each experiment this process is repeated 500 times and then the classification procedures are compared. We find that several methods are superior to the currently used quadratic discriminant method for risk estimation for these data. The implications of these results for prenatal screening programs are discussed.     

The composite Milstein methods for the numerical solution of Stratonovich stochastic differential equations

In this paper, we present two composite Milstein methods for the strong solution of Stratonovich stochastic differential equations driven by d-dimensional Wiener processes. The composite Milstein methods are a combination of semi-implicit and implicit Milstein methods. The criterion for choosing either the implicit or the semi-implicit method at each step of the numerical solution is given. The stability and convergence properties of the proposed methods are analyzed for the linear test equation. It is shown that the proposed methods converge to the exact solution in Stratonovich sense. In addition, the stability properties of our methods are found to be superior to those of the Milstein and the composite Euler methods. The convergence properties for the nonlinear case are shown numerically to be the same as the linear case. Hence, the proposed methods are a good candidate for the solution of stiff SDEs.     

The global Hessenberg and CMRH methods for linear systems with multiple right-hand sides

In this paper, we introduce two new methods for solving large sparse nonsymmetric linear systems with several right-hand sides. These methods are the global Hessenberg and global CMRH methods. Using the global Hessenberg process, these methods are less expensive than the global FOM and global GMRES methods [9]. Theoretical results about the new methods are given, and experimental results that show good performances of these new methods are presented.     

A new two-parameter family of nonlinear conjugate gradient methods

Conjugate gradient methods are an important class of methods for unconstrained optimization, especially for large-scale problems. Recently, they have been much studied. In this paper, we propose a new two-parameter family of conjugate gradient methods for unconstrained optimization. The two-parameter family of methods not only includes the already existing three practical nonlinear conjugate gradient methods, but has other family of conjugate gradient methods as subfamily. The two-parameter family of methods with the Wolfe line search is shown to ensure the descent property of each search direction. Some general convergence results are also established for the two-parameter family of methods. The numerical results show that this method is efficient for the given test problems. In addition, the methods related to this family are uniformly discussed.     

A note on some improvements of the simultaneous methods for determination of polynomial zeros

Applying Gauss-Seidel approach to the improvements of two simultaneous methods for finding polynomial zeros, presented in [9], two iterative methods with faster convergence are obtained. The lower bounds of the R-order of convergence for the accelerated methods are given. The improved methods and their accelerated modifications are discussed in view of the convergence order and the number of numerical operations. The considered methods are illustrated numerically in the example of an algebraic equation.     

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.     

Newton-type methods for constrained optimization with nonregular constraints

The most important classes of Newton-type methods for solving constrained optimization problems are discussed. These are the sequential quadratic programming methods, active set methods, and semismooth Newton methods for Karush-Kuhn-Tucker systems. The emphasis is placed on the behavior of these methods and their special modifications in the case where assumptions concerning constraint qualifications are relaxed or altogether dropped. Applications to optimization problems with complementarity constraints are examined.     

求解界约束优化问题的有效集算法综述

主要介绍了求解界约束优化问题的有效集方法,包括投影共轭梯度法和有效集识别函数法,讨论了各自的优点和不足.最后,指出了有效集法的研究趋势及应用前景.     

Stability properties of numerical methods for solving delay differential equations

Stability properties of numerical methods for delay differential equations are considered. Some suitable definitions for the stability of the numerical methods are included and Runge-Kutta type methods satisfying these properties are tested on a numerical example.     

On the convergence of modified Newton methods for solving equations containing a non-differentiable term

We provide a semilocal convergence analysis for certain modified Newton methods for solving equations containing a non-differentiable term. The sufficient convergence conditions of the corresponding Newton methods are often taken as the sufficient conditions for the modified Newton methods. That is why the latter methods are not usually treated separately from the former. However, here we show that weaker conditions, as well as a finer error analysis than before can be obtained for the convergence of modified Newton methods. Numerical examples are also provided.     

Numerical Methods for the Maxnear Criterion of Multiple-Sets Canonical Analysis

This paper deals with numerical methods for the Maxnear criterion of multiple-sets canonical analysis. Optimality conditions are derived. Upper and lower bounds of the optimal objective function value are presented. Two iterative methods are proposed. One is an alternating variable method, and the other called Gauss-Seidel method is an inexact version of the alternating variable method. Convergence of these methods are analyzed. A starting point strategy is suggested for both methods. Numerical results are presented to demonstrate the efficiency of these methods and the starting point strategy.     

Empirical likelihood confidence intervals for hazard and density functions under right censorship

In this paper, we use smoothed empirical likelihood methods to construct confidence intervals for hazard and density functions under right censorship. Some empirical log-likelihood ratios for the hazard and density functions are obtained and their asymptotic limits are derived. Approximate confidence intervals based on these methods are constructed. Simulation studies are used to compare the empirical likelihood methods and the normal approximation methods in terms of coverage accuracy. It is found that the empirical likelihood methods provide better inference.     

Dual interior point methods for linear semidefinite programming problems

Dual interior point methods for solving linear semidefinite programming problems are proposed. These methods are an extension of dual barrier-projection methods for linear programs. It is shown that the proposed methods converge locally at a linear rate provided that the solutions to the primal and dual problems are nondegenerate.