Journal of Optimization Theory and Applications最新文献

In this study, we propose an algorithmic framework for solving a class of optimal control problems. This class is associated with the minimum-time interception of moving target problems, where a plant with a given state equation must approach a moving target whose trajectory is known a priori. Our framework employs an analytical description of the distance from an arbitrary point to the reachable set of the plant. The proposed algorithm is always convergent and cannot be improved without losing the guarantee of its convergence to the correct solution for arbitrary Lipschitz continuous trajectories of the moving target. In practice, it is difficult to obtain an analytical description of the distance to the reachable set for the sophisticated state equation of the plant. Nevertheless, it was shown that the distance can be obtained for some widely used models, such as the Dubins car, in an explicit form. Finally, we illustrate the generality and effectiveness of the proposed framework for simple motions and the Dubins model.

This work introduces a new Riemannian optimization method for registering open parameterized surfaces with a constrained global optimization approach. The proposed formulation leads to a rigorous theoretic foundation and guarantees the existence and the uniqueness of a global solution. We also propose a new Bayesian clustering approach where local distributions of surfaces are modeled with spherical Gaussian processes. The maximization of the posterior density is performed with Hamiltonian dynamics which provide a natural and computationally efficient spherical Hamiltonian Monte Carlo sampling. Experimental results demonstrate the efficiency of the proposed method.

Many tasks in real life scenarios can be naturally formulated as nonconvex optimization problems. Unfortunately, to date, the iterative numerical methods to find even only the local minima of these nonconvex cost functions are extremely slow and strongly affected by the initialization chosen. We devise a predictor–corrector strategy that efficiently computes locally optimal solutions to these problems. An initialization-free convex minimization allows to predict a global good preliminary candidate, which is then corrected by solving a parameter-free nonconvex minimization. A simple algorithm, such as alternating direction method of multipliers works surprisingly well in producing good solutions. This strategy is applied to the challenging problem of decomposing a 1D signal into semantically distinct components mathematically identified by smooth, piecewise-constant, oscillatory structured and unstructured (noise) parts.

In this paper, we investigate bounds of solution set of the polynomial complementarity problem. When a polynomial complementarity problem has a solution, we propose a lower bound of solution norm by entries of coefficient tensors of the polynomial. We prove that the proposing lower bound is larger than some existing lower bounds appeared in tensor complementarity problems and polynomial complementarity problems. When the solution set of a polynomial complementarity problem is nonempty, and the coefficient tensor of the leading term of the polynomial is an (R_0)-tensor, we propose a new upper bound of solution norm of the polynomial complementarity problem by a quantity defining by an optimization problem. Furthermore, we prove that when coefficient tensors of the polynomial are partially symmetric, the proposing lower bound formula with respect to tensor tuples reaches the maximum value, and the proposing upper bound formula with respect to tensor tuples reaches the minimum value. Finally, by using such partial symmetry, we obtain bounds of solution norm by coefficients of the polynomial.

Taking inspiration from what is commonly done in single-objective optimization, most local algorithms proposed for multiobjective optimization extend the classical iterative scalar methods and produce sequences of points able to converge to single efficient points. Recently, a growing number of local algorithms that build sequences of sets has been devised, following the real nature of multiobjective optimization, where the aim is that of approximating the efficient set. This calls for a new analysis of the necessary optimality conditions for multiobjective optimization. We explore conditions for sets of points that share the same features of the necessary optimality conditions for single-objective optimization. On the one hand, from a theoretical point of view, these conditions define properties that are necessarily satisfied by the (weakly) efficient set. On the other hand, from an algorithmic point of view, any set that does not satisfy the proposed conditions can be easily improved by using first-order information on some objective functions. We analyse both the unconstrained and the constrained case, giving some examples.

We consider how to solve a class of non-Lipschitz mathematical programs with equilibrium constraints (MPEC) where the objective function involves a non-Lipschitz sparsity-inducing function and other functions are smooth. Solving the non-Lipschitz MPEC is highly challenging since the standard constraint qualifications fail due to the existence of equilibrium constraints and the subdifferential of the objective function is unbounded due to the existence of the non-Lipschitz function. On the one hand, for tackling the non-Lipschitzness of the objective function, we introduce a novel class of locally Lipschitz approximation functions that consolidate and unify a diverse range of existing smoothing techniques for the non-Lipschitz function. On the other hand, we use the Kanzow and Schwartz regularization scheme to approximate the equilibrium constraints since this regularization can preserve certain perpendicular structure as in equilibrium constraints, which can induce better convergence results. Then an approximation method is proposed for solving the non-Lipschitz MPEC and its convergence is established under weak conditions. In contrast with existing results, the proposed method can converge to a better stationary point under weaker qualification conditions. Finally, a computational study on the sparse solutions of linear complementarity problems is presented. The numerical results demonstrate the effectiveness of the proposed method.

In this paper, we investigate the asymptotic properties of a particular class of state-dependent sweeping processes. While extensive research has been conducted on the existence and uniqueness of solutions for sweeping processes, there is a scarcity of studies addressing their behavior in the limit of large time. Additionally, we introduce novel algorithms designed for the resolution of quasi-variational inequalities. As a result, we introduce a new derivative-free algorithm to find zeros of nonsmooth Lipschitz continuous mappings with a linear convergence rate. This algorithm can be effectively used in nonsmooth and nonconvex optimization problems that do not possess necessarily second-order differentiability conditions of the data.

This paper studies the convergence of the proximal point method for quasiconvex functions in finite dimensional complete Riemannian manifolds. We prove initially that, in the general case, when the objective function is proper and lower semicontinuous, each accumulation point of the sequence generated by the method, if it exists, is a limiting critical point of the function. Then, under the assumptions that the sectional curvature of the manifold is bounded above by some non negative constant and the objective function is quasiconvex we analyze two cases. When the constant is zero, the global convergence of the algorithm to a limiting critical point is assured and if it is positive, we prove the local convergence for a class of quasiconvex functions, which includes Lipschitz functions.

Much is known about when a locally optimal solution depends in a single-valued Lipschitz continuous way on the problem’s parameters, including tilt perturbations. Much less is known, however, about when that solution and a uniquely determined multiplier vector associated with it exhibit that dependence as a primal–dual pair. In classical nonlinear programming, such advantageous behavior is tied to the combination of the standard strong second-order sufficient condition (SSOC) for local optimality and the linear independent gradient condition (LIGC) on the active constraint gradients. But although second-order sufficient conditons have successfully been extended far beyond nonlinear programming, insights into what should replace constraint gradient independence as the extended dual counterpart have been lacking. The exact answer is provided here for a wide range of optimization problems in finite dimensions. Behind it are advances in how coderivatives and strict graphical derivatives can be deployed. New results about strong metric regularity in solving variational inequalities and generalized equations are obtained from that as well.

Many solution algorithms for multiobjective optimization problems are based on scalarization methods that transform the multiobjective problem into a scalar-valued optimization problem. In this article, we study the theory of weighted (p)-norm scalarizations. These methods minimize the distance induced by a weighted (p)-norm between the image of a feasible solution and a given reference point. We provide a comprehensive theory of the set of eligible weights and, in particular, analyze the topological structure of the normalized weight set. This set is composed of connected subsets, called weight set components which are in a one-to-one relation with the set of optimal images of the corresponding weighted (p)-norm scalarization. Our work generalizes and complements existing results for the weighted sum and the weighted Tchebycheff scalarization and provides new insights into the structure of the set of all Pareto optimal solutions.