A Constrained Optimization Approach for LCP

In this paper, LCP is converted to an equivalent nonsmooth nonlinear equation system H(x,y) = 0 by using the famous NCP function-Fischer-Burmeister function. Note that some equations in H(x, y) = 0 are nonsmooth and nonlinear hence difficult to solve while the others are linear hence easy to solve. Then we further convert the nonlinear equation system H(x, y) = 0 to an optimization problem with linear equality constraints. After that we study the conditions under which the K-T points of the optimization problem are the solutions of the original LCP and propose a method to solve the optimization problem. In this algorithm, the search direction is obtained by solving a strict convex programming at each iterative point, However, our algorithm is essentially different from traditional SQP method. The global convergence of the method is proved under mild conditions. In addition, we can prove that the algorithm is convergent superlinearly under the conditions: M is P0 matrix and the limit point is a strict complementarity solution of LCP. Preliminary numerical experiments are reported with this method.     

A derivative-free algorithm for systems of nonlinear inequalities

Recently a new derivative-free algorithm has been proposed for the solution of linearly constrained finite minimax problems. This derivative-free algorithm is based on a smoothing technique that allows one to take into account the non-smoothness of the max function. In this paper, we investigate, both from a theoretical and computational point of view, the behavior of the minmax algorithm when used to solve systems of nonlinear inequalities when derivatives are unavailable. In particular, we show an interesting property of the algorithm, namely, under some mild conditions regarding the regularity of the functions defining the system, it is possible to prove that the algorithm locates a solution of the problem after a finite number of iterations. Furthermore, under a weaker regularity condition, it is possible to show that an accumulation point of the sequence generated by the algorithm exists which is a solution of the system. Moreover, we carried out numerical experimentation and comparison of the method against a standard pattern search minimization method. The obtained results confirm that the good theoretical properties of the method correspond to interesting numerical performance. Moreover, the algorithm compares favorably with a standard derivative-free method, and this seems to indicate that extending the smoothing technique to pattern search algorithms can be beneficial.     

A boundary element method to price time-dependent double barrier options

In this paper we propose a new method for pricing double-barrier options with moving barriers under the Black-Scholes and the CEV models. First of all, by applying a variational technique typical of the boundary element method, we derive an integral representation of the double-barrier option price in which two of the integrand functions are not given explicitly but must be obtained solving a system of Volterra integral equations of the first kind. Second, we develop an ad hoc numerical method to regularize and solve the system of integral equations obtained. Several numerical experiments are carried out showing that the overall algorithm is extraordinarily fast and accurate, even if the barriers are not differentiable functions. Moreover the numerical method presented in this paper performs significantly better than the finite difference approach.     

Modified Newton methods for solving a semismooth reformulation of monotone complementarity problems

In this paper, we propose a Newton-type method for solving a semismooth reformulation of monotone complementarity problems. In this method, a direction-finding subproblem, which is a system of linear equations, is uniquely solvable at each iteration. Moreover, the obtained search direction always affords a direction of sufficient decrease for the merit function defined as the squared residual for the semismooth equation equivalent to the complementarity problem. We show that the algorithm is globally convergent under some mild assumptions. Next, by slightly modifying the direction-finding problem, we propose another Newton-type method, which may be considered a restricted version of the first algorithm. We show that this algorithm has a superlinear, or possibly quadratic, rate of convergence under suitable assumptions. Finally, some numerical results are presented. Supported by Research Fellowships of the Japan Society for the Promotion of Science for Young Scientists. Supported in part by the Scientific Research Grant-in-Aid from the Ministry of Education, Science and Culture, Japan.     

一个新的SQP方法及其超线性收敛性

由Ｗｉｌｓｏｎ,Ｈａｎ,Ｐｏｗｅｌｌ发展的ＳＱＰ技术是解非线性规划的最有效的方法之一,但是,如果其中的二次子规划问题无可行解或者其搜索方向向量无界,该方法ａｎ和Ｂｕｒｋｅ「３」,周广路「２」分别对二次规划问题作了修正,克服了上述矛盾,本文在「２」的基础上,进上步修正,证明在Ａｒｍｉｊｏ搜索下算法具有全局收敛性,并通过解一辅助线性方程组,利用弧式搜索,得出该方法具有超线性收敛性。     

Superlinear/quadratic smoothing Broyden-like method for the generalized nonlinear complementarity problem

In this article, we first reformulate the generalized nonlinear complementarity problem (GNCP) over a polyhedral cone as a smoothing system of equations and then suggest a smoothing Broyden-like method for solving it. The proposed algorithm has to solve only one system of nonhomogeneous linear equations, perform only one line search and update only one matrix per iteration. We show that the iteration sequence generated by the proposed algorithm converges globally and superlinearly under suitable conditions. Furthermore, the algorithm has local quadratic convergence under mild assumptions. Some numerical examples are given to illustrate the performance and efficiency of the presented algorithm.     

Conjugate gradient path method without line search technique for derivative-free unconstrained optimization

In this paper, we propose conjugate gradient path method for solving derivative-free unconstrained optimization. The iterative direction is obtained by constructing and solving quadratic interpolation model of the objective function with conjugate gradient methods. The global convergence and local superlinear convergence rate of the proposed algorithm are established under some reasonable conditions. Finally, the numerical results are reported to show the effectiveness of the proposed algorithm.     

不等式约束优化拟正则条件下的超线性收敛FSSLE算法

对不等式约束优化问题。提出一个可行序列线性方程组（FSSLE）算法。该算法每次迭代只需求解两个具有相同系数矩阵的线性方程组,因而计算量较小。在一定条件下,算法具有全局收敛性。在没有严格互补条件、比强二阶充分条件弱的拟正则条件下,证明了算法具有超线性收敛性并用数值试验表明其有效性。     

Feasible Direction Interior-Point Technique for Nonlinear Optimization

We propose a feasible direction approach for the minimization by interior-point algorithms of a smooth function under smooth equality and inequality constraints. It consists of the iterative solution in the primal and dual variables of the Karush–Kuhn–Tucker first-order optimality conditions. At each iteration, a descent direction is defined by solving a linear system. In a second stage, the linear system is perturbed so as to deflect the descent direction and obtain a feasible descent direction. A line search is then performed to get a new interior point and ensure global convergence. Based on this approach, first-order, Newton, and quasi-Newton algorithms can be obtained. To introduce the method, we consider first the inequality constrained problem and present a globally convergent basic algorithm. Particular first-order and quasi-Newton versions of this algorithm are also stated. Then, equality constraints are included. This method, which is simple to code, does not require the solution of quadratic programs and it is neither a penalty method nor a barrier method. Several practical applications and numerical results show that our method is strong and efficient.     

An improved sequential quadratic programming algorithm for solving general nonlinear programming problems

In this paper, a class of general nonlinear programming problems with inequality and equality constraints is discussed. Firstly, the original problem is transformed into an associated simpler equivalent problem with only inequality constraints. Then, inspired by the ideals of the sequential quadratic programming (SQP) method and the method of system of linear equations (SLE), a new type of SQP algorithm for solving the original problem is proposed. At each iteration, the search direction is generated by the combination of two directions, which are obtained by solving an always feasible quadratic programming (QP) subproblem and a SLE, respectively. Moreover, in order to overcome the Maratos effect, the higher-order correction direction is obtained by solving another SLE. The two SLEs have the same coefficient matrices, and we only need to solve the one of them after a finite number of iterations. By a new line search technique, the proposed algorithm possesses global and superlinear convergence under some suitable assumptions without the strict complementarity. Finally, some comparative numerical results are reported to show that the proposed algorithm is effective and promising.     

Interval optimization based line sampling method for fuzzy and random reliability analysis

For structural system with fuzzy variables as well as random variables, a novel algorithm for obtaining membership function of fuzzy reliability is presented on interval optimization based Line Sampling (LS) method. In the presented algorithm, the value domain of the fuzzy variables under the given membership level is firstly obtained according to their membership functions. Then, in the value domain of the fuzzy variables, bounds of reliability of the structure are obtained by the nesting analysis of the interval optimization, which is performed by modern heuristic methods, and reliability analysis, which is achieved by the LS method in the reduced space of the random variables. In this way the uncertainties of the input variables are propagated to the safety measurement of the structure, and the membership function of the fuzzy reliability is obtained. The presented algorithm not only inherits the advantage of the direct Monte Carlo method in propagating and distinguishing the fuzzy and random uncertainties, but also can improve the computational efficiency tremendously in case of acceptable precision. Several examples are used to illustrate the advantages of the presented algorithm.     

A smoothing method for second order cone complementarity problem

In this paper, the second order cone complementarity problem is studied. Based on a perturbed symmetrically smoothing function, which has coerciveness under proper conditions, we present a smoothing Newton method for this problem. The boundedness of the level set can be obtained from the coerciveness, which plays an important role in the convergence analysis. Furthermore, the proposed algorithm for the reformulation has no restrictions on the starting point and solves only one system of equations. Preliminary numerical results indicate that the algorithm is effective.     

The interior proximal extragradient method for solving equilibrium problems

In this article we present a new and efficient method for solving equilibrium problems on polyhedra. The method is based on an interior-quadratic proximal term which replaces the usual quadratic proximal term. This leads to an interior proximal type algorithm. Each iteration consists in a prediction step followed by a correction step as in the extragradient method. In a first algorithm each of these steps is obtained by solving an unconstrained minimization problem, while in a second algorithm the correction step is replaced by an Armijo-backtracking linesearch followed by an hyperplane projection step. We prove that our algorithms are convergent under mild assumptions: pseudomonotonicity for the two algorithms and a Lipschitz property for the first one. Finally we present some numerical experiments to illustrate the behavior of the proposed algorithms.     

Nonmonotone Inexact Newton Method for the Extended Linear Complementarity Problem

This paper presents a nonmonotone inexact Newton-type method for the extended linear complementarity problem (ELCP). We first reformulate the optimization system of the ELCP problem into a system of smoothed equations. Then we solve this system by a nonmonotone inexact Newton-type algorithm. The global convergence is obtained and numerical tests for some classes of ELCP include linear complementarity, horizontal linear complementarity, and generalized linear complementarity problems are also given to show the e?ciency of the proposed algorithm.     

Hybrid event-triggered and impulsive control for time-delay systems

In this paper, we study the problem of hybrid event-triggered control for a class of nonlinear time-delay systems. Using a Razumikhin-type input-to-state stability result for time-delay systems, we design an event-triggered control algorithm to stabilize the given time-delay system. In order to exclude Zeno behavior, we combine the impulsive control mechanism with our event-triggered strategy. In this sense, the proposed algorithm is a hybrid impulsive and event-triggered strategy. Sufficient conditions for the stabilization of the nonlinear systems with time delay are obtained by using Lyapunov method and Razumikhin technique. Numerical simulations are provided to show the effectiveness of our theoretical results.     

Finite-difference method for singular nonlinear systems

This paper presents a method for solving nonlinear system with singular Jacobian at the solution. The convergence rate in the case of singularity deteriorates and one way to accelerate convergence is to form bordered system. A local algorithm, with finite-difference approximations, for forming and solving such system is proposed in this paper. To overcome the need that initial approximation has to be very close to the solution, we also propose a method which is a combination of descent method with finite-differences and local algorithm. Some numerical results obtained on relevant examples are presented.     

多约束条件下三状态复杂系统的优化分配

研究了三状态复杂系统在多约束条件下可靠度的问题.以串-并联系统为研究对象,利用选取重要度来提高系统可靠度的方法,获得了重要度对系统可靠度有着重要影响的结果,并通过算法、例题对实例进行了验证.     

A self-adaptive trust region method for the extended linear complementarity problems

By using some NCP functions, we reformulate the extended linear complementarity problem as a nonsmooth equation. Then we propose a self-adaptive trust region algorithm for solving this nonsmooth equation. The novelty of this method is that the trust region radius is controlled by the objective function value which can be adjusted automatically according to the algorithm. The global convergence is obtained under mild conditions and the local superlinear convergence rate is also established under strict complementarity conditions. This work is supported by National Natural Science Foundation of China (No. 10671126) and Shanghai Leading Academic Discipline Project (S30501).     

一个等式约束问题的拟Newton—信赖域型方法及其收敛性

在[1]中,Vardi提出一个信赖域方法,而收敛性证明却是在精确λ-搜索下给出的,本文在[1]的基础上提出一个新的算法-拟Newton-信赖域型算法,并证明该算法是全局收敛的,通过利用二阶修正技术去修正该算法,我们证明了该算法是局部超线性收敛的。