A Two-Step Modulus-Based Matrix Splitting Iteration Method Without Auxiliary Variables for Solving Vertical Linear Complementarity Problems

In this paper,a two-step iteration method is established which can be viewed as a generalization of the existing modulus-based methods for vertical linear complementarity problems given by He and Vong(Appl.Math.Lett.134:108344,2022).The convergence analysis of the proposed method is established,which can improve the existing results.Numerical examples show that the proposed method is efficient with the two-step technique.

Parametric variational solution of linear-quadratic optimal control problems with control inequality constraints

A parametric variational principle and the corresponding numerical algo- rithm are proposed to solve a linear-quadratic （LQ） optimal control problem with control inequality constraints. Based on the parametric variational principle, this control prob- lem is transformed into a set of Hamiltonian canonical equations coupled with the linear complementarity equations, which are solved by a linear complementarity solver in the discrete-time domain. The costate variable information is also evaluated by the proposed method. The parametric variational algorithm proposed in this paper is suitable for both time-invariant and time-varying systems. Two numerical examples are used to test the validity of the proposed method. The proposed algorithm is used to astrodynamics to solve a practical optimal control problem for rendezvousing spacecrafts with a finite low thrust. The numerical simulations show that the parametric variational algorithm is ef- fective for LQ optimal control problems with control inequality constraints.

A New Type of Solution Method for the Generalized Linear Complementarity Problem over a Polyhedral Cone

This paper addresses the generalized linear complementarity problem （GLCP） over a polyhedral cone. To solve the problem, we first equivalently convert the problem into an affine variational inequalities problem over a closed polyhedral cone, and then propose a new type of method to solve the GLCP based on the error bound estimation. The global and R-linear convergence rate is established. The numerical experiments show the efficiency of the method.

A New Conjugate Gradient Projection Method for Solving Stochastic Generalized Linear Complementarity Problems

In this paper, a class of the stochastic generalized linear complementarity problems with finitely many elements is proposed for the first time. Based on the Fischer-Burmeister function, a new conjugate gradient projection method is given for solving the stochastic generalized linear complementarity problems. The global convergence of the conjugate gradient projection method is proved and the related numerical results are also reported.

A linear complementarity model for multibody systems with frictional unilateral and bilateral constraints

The Lagrange-I equations and measure differential equations for multibody systems with unilateral and bilateral constraints are constructed. For bilateral constraints, frictional forces and their impulses contain the products of the filled-in relay function induced by Coulomb friction and the absolute values of normal constraint reactions. With the time-stepping impulse-velocity scheme, the measure differential equations are discretized. The equations of horizontal linear complementarity problems （HLCPs）, which are used to compute the impulses, are constructed by decomposing the absolute function and the filled-in relay function. These HLCP equations degenerate into equations of LCPs for frictional unilateral constraints, or HLCPs for frictional bilateral constraints. Finally, a numerical simulation for multibody systems with both unilateral and bilateral constraints is presented.

A Potential-Reduction Algorithm for Linear Complementarity Problems

Feasible-interior-point algorithms start from a strictly feasible interior point, but infeassible-interior-point algorithms just need to start from an arbitrary positive point, we give a potential reduction algorithm from an infeasible-starting-point for a class of non-monotone linear complementarity problem. Its polynomial complexity is analyzed. After finite iterations the algorithm produces an approximate solution of the problem or shows that there is no feasible optimal solution in a large region. Key words linear complementarity problems - infeasible-starting-point - P-matrix - potential function CLC number O 221 Foundation item: Supported by the National Natural Science Foundation of China (70371032) and the Doctoral Educational Foundation of China of the Ministry of Education (20020486035)Biography: Wang Yan-jin (1976-), male, Ph. D candidate, research direction: optimal theory and method.

GLOBAL LINEAR AND QUADRATIC ONE-STEP SMOOTHING NEWTON METHOD FOR VERTICAL LINEAR COMPLEMENTARITY PROBLEMS

A one_step smoothing Newton method is proposed for solving the vertical linear complementarity problem based on the so_called aggregation function. The proposed algorithm has the following good features: (ⅰ) It solves only one linear system of equations and does only one line search at each iteration; (ⅱ) It is well_defined for the vertical linear complementarity problem with vertical block P 0 matrix and any accumulation point of iteration sequence is its solution.Moreover, the iteration sequence is bounded for the vertical linear complementarity problem with vertical block P 0+R 0 matrix; (ⅲ) It has both global linear and local quadratic convergence without strict complementarity. Many existing smoothing Newton methods do not have the property (ⅲ).

New Mehrotra's second order predictor-corrector algorithm for P_*(κ) linear complementarity problems

It has been shown in various papers that most interior-point algorithms for linear optimization and their analysis can be generalized to P_＊（κ） linear complementarity problems.This paper presents an extension of the recent variant of Mehrotra＇s second order algorithm for linear optimijation.It is shown that the iteration-complexity bound of the algorithm is O（4κ ＋ 3）√14κ ＋ 5 nlog（x0）Ts0/ε,which is similar to that of the corresponding algorithm for linear optimization.

Interior-point algorithm based on general kernel function for monotone linear complementarity problem

A polynomial interior-point algorithm is presented for monotone linear complementarity problem （MLCP） based on a class of kernel functions with the general barrier term, which are called general kernel functions. Under the mild conditions for the barrier term, the complexity bound of algorithm in terms of such kernel function and its derivatives is obtained. The approach is actually an extension of the existing work which only used the specific kernel functions for the MLCP.

CONVERGENCE OF A MODIFIED SLP ALGORITHM FOR THE EXTENDED LINEAR COMPLEMENTARITY PROBLEM

A modified sequential linear programming algorithm is presented, whose subproblem is always solvable, for the extended linear complementarity problem (XLCP), the global convergence of the algorithm under assumption of X-row sufficiency or X-colunm monotonicity is proved. As a result, a sufficient condition for existence and boundedness of solution to the XLCP are obtained.

The Smoothing Newton Method for Solving the Extended Linear Complementarity Problem

The extended linear complementarity problem(denoted by ELCP) can be reformulated as the solution of a nonsmooth system of equations. By the symmetrically perturbed CHKS smoothing function, the ELCP is approximated by a family of parameterized smooth equations. A one-step smoothing Newton method is designed for solving the ELCP. The proposed algorithm is proved to be globally convergent under suitable assumptions.

The Nonlinear Lopsided HSS-Like Modulus-Based Matrix Splitting Iteration Method for Linear Complementarity Problems with Positive-Definite Matrices

In this paper,by means of constructing the linear complementarity problems into the corresponding absolute value equation,we raise an iteration method,called as the nonlinear lopsided HSS-like modulus-based matrix splitting iteration method,for solving the linear complementarity problems whose coefficient matrix in R^(n×n)is large sparse and positive definite.From the convergence analysis,it is appreciable to see that the proposed method will converge to its accurate solution under appropriate conditions.Numerical examples demonstrate that the presented method precede to other methods in practical implementation.

Homotopy Continuation Method for Linear Complementarity Problems

It is well known that a linear complementarity problem (LCP) can be formulated as a system of nonsmooth equations F(x) = 0, where F is a map from Rninto itself. Using the aggregate function, we construct a smooth Newton homotopy H(x,t) = 0. Under certain assumptions, we prove the existence of a smooth path defined by the Newton homotopy which leads to a solution of the original problem, and study limiting properties of the homotopy path.

A CLASS OF GENERALIZED MULTISPLITTING RELAXATION METHODS FOR LINEAR COMPLEMENTARITY PROBLEMS

Abstract In this paper,a class of generalized parallel matrix multisplitting relaxation methods for solving linear complementarity problems on the high speed multiprocessor systems is set up.This class of methods not only includes all the existing relaxation methods for the linear complementarity problems,but also yields a lot of novel ones in the sense of multisplitting.We establish the convergence theories of this class of generalized parallel multisplitting relaxation methods under the condition that the system matrix is an H matrix with positive diagonal elements.

LCP method for a planar passive dynamic walker based on an event-driven scheme

The main purpose of this paper is to present a linear complementarity problem （LCP） method for a planar passive dynamic walker with round feet based on an event-driven scheme. The passive dynamic walker is treated as a planar multi-rigid-body system. The dynamic equations of the passive dynamic walker are obtained by using Lagrange＇s equations of the second kind. The normal forces and frictional forces acting on the feet of the passive walker are described based on a modified Hertz contact model and Coulomb＇s law of dry friction. The state transition problem of stick-slip between feet and floor is formulated as an LCR which is solved with an event-driven scheme. Finally, to validate the methodology, four gaits of the walker are simulated： the stance leg neither slips nor bounces; the stance leg slips without bouncing; the stance leg bounces without slipping; the walker stands after walking several steps.

Comparison and analysis of two Coulomb friction models on the dynamic behavior of slider-crank mechanism with a revolute clearance joint

The objective of this study is to investigate the effects of the Coulomb dry friction model versus the modified Coulomb friction model on the dynamic behavior of the slider-crank mechanism with a revolute clearance joint. The normal and tangential forces acting on the contact points between the journal and the bearing are described by using a Hertzian-based contact force model and the Coulomb friction models, respectively.The dynamic equations of the mechanism are derived based on the Lagrange equations of the first kind and the Baumgarte stabilization method. The frictional force is solved via the linear complementarity problem(LCP) algorithm and the trial-and-error algorithm.Finally, three numerical examples are given to show the influence of the two Coulomb friction models on the dynamic behavior of the mechanism. Numerical results show that due to the stick friction, the slider-crank mechanism may exhibit stick-slip motion and can balance at some special positions, while the mechanism with ideal joints cannot.

NON-INTERIOR SMOOTHING ALGORITHM FOR FRICTIONAL CONTACT PROBLEMS

A new algorithm for solving the three-dimensional elastic contact problem with friction is presented. The algorithm is a non-interior smoothing algorithm based on an NCP-function. The parametric variational principle and parametric quadratic programming method were applied to the analysis of three-dimensional frictional contact problem. The solution of the contact problem was finally reduced to a linear complementarity problem, which was reformulated as a system of nonsmooth equations via an NCP-function. A smoothing approximation to the nonsmooth equations was given by the aggregate function. A Newton method was used to solve the resulting smoothing nonlinear equations. The algorithm presented is easy to understand and implement. The reliability and efficiency of this algorithm are demonstrated both by the numerical experiments of LCP in mathematical way and the examples of contact problems in mechanics.

Robust H∞ output tracking control of uncertain networked control systems

This paper investigates a problem of robust output tracking control of networked control systems(NCSs) with network-induced delays, packet dropouts, parameter uncertainties and external disturbances. Firstly, an augmented model of time-delay system is proposed for networked tracking control systems. Then, considering the piecewise differentiable characteristic of time-delay, the criterion to output tracking performance analysis and controller design are derived by using an approach of free weighting matrix, reciprocally convex and cone complementarity linearization(CCL). Finally, simulation results of numerical examples show the effectiveness of the proposed approach, and illustrate the advantages of the proposed criteria which outperform previous criteria in the literature.

Some Remarks On E＇-matrices

The class of E＇-matrices introduced in [4] is a subclass of Q0-matrices whose proper principal submatrices are Q--matrices. In this paper, we investigate the relation of the class E＇ to other known matrix classes and obtain a series of necessary and sufficient conditions for E＇-matrices.

H-infinity output feedback control for descriptor systems with delayed states

This paper considers the H-infinity dynamic output feedback control for descriptor systems with delay in states. The controller is a descriptor system without delay. Several equivalent sufficient conditions for the existence of one descriptor dynamic controller without impulsive models are given. Furthermore the explicit expression of the desired controller is obtained. The detailed design of the controller is presented using the cone complementarity linearization iterative algorithm and the LMI method. A ntumerical example is shown to illustrate the designed method.