非光滑方程组牛顿法的全局收敛性分析

通过定义一种新的*-微分,本文给出了局部Lipschitz非光滑方程组的牛顿法,并对其全局收敛性进行了研究.该牛顿法结合了非光滑方程组的局部收敛性和全局收敛性.最后,我们把这种牛顿法应用到非光滑函数的光滑复合方程组问题上,得到了较好的收敛性.     

凸约束非光滑方程组一个新的谱梯度投影算法

基于寻找分离超平面的三种经典线搜索技术,本文提出了一种自适应线搜索技术.结合谱梯度投影法,提出了凸约束非光滑单调方程组的一个谱梯度投影算法.该算法不需要计算和存储任何矩阵,因而适合求解大规模非光滑的非线性单调方程组.在较弱的条件下,证明了方法的全局收敛性,并分析了算法的收敛率.数值试验结果表明算法是有效的和鲁棒的.     

一个求解广义圆锥互补问题的无导数光滑算法

本文研究了一个求解广义圆锥互补问题的无导数光滑算法.利用光滑函数将广义圆锥互补问题等价转化成一个光滑方程组,然后再利用牛顿法求解此方程组.该算法采用了一种新的非单调无导数线搜索技术,并且在适当条件下具有全局和局部二次收敛性质.数值实验结果表明算法是非常有效的.     

非光滑非线性互补问题的牛顿法

研究了非光滑的非线性互补问题. 首先将非光滑的非线性互补问题转化为一个非光滑方程组,然后用牛顿法求解这个非光滑方程组. 在该牛顿法中,每次迭代只需一个原始函数B-微分中的一个元素. 最后证明了该牛顿法的超线性收敛性.     

一类非线性二阶锥规划的非光滑牛顿法

本文研究了非线性二阶锥规划问题.利用投影映射将非线性二阶锥规划问题的KKT最优性条件转化成非光滑方程组,获得了一个修正的中心路径非光滑牛顿法.在适当的条件下保证方程组的B-次微分在任意点都可逆,并且证明算法具有全局收敛性.     

二阶锥约束随机变分不等式问题的数值方法研究

本文研究二阶锥约束随机变分不等式(SOCCSVI)问题,运用样本均值近似(SAA)方法结合光滑Fischer-Burmeister互补函数来求解该问题.首先,将SOCCSVI问题的Karush-Kuhn-Tucker系统转化为与之等价的方程组,并证明了该方程组的雅可比矩阵的非奇异性.其次,构造了光滑牛顿算法求解该方程组.最后,文章给出了两个数值实验证明了算法的有效性.     

求解非线性互补问题的一类光滑Broyden-like方法

非线性互补问题可以等价地转换为光滑方程组来求解. 基于一种新的非单调线搜索准则, 提出了求解非线性互补问题等价光滑方程组的一类新的非单调光滑 Broyden-like 算法.在适当的假设条件下, 证明了该算法的全局收敛性与局部超线性收敛性. 数值实验表明所提出的算法是有效的.     

求非线性规划的非可行滤子无二次子规划方法

本文研究了不等式约束的非线性规划问题.利用带滤子的无二次子规划(QP-free)非可行域方法,构造一个等价于原约束问题的一阶KKT条件的非光滑方程组,给出解这个方程组的迭代算法,并获得算法的全局收敛性.     

依赖凝聚函数求解非线性互补问题的一种微分方程方法

利用凝聚函数一致逼近非光滑极大值函数的性质,将非线性互补问题转化为参数化光滑方程组.然后,对此方程组给出了一种微分方程解法,并且证明了非线性互补问题的解是微分方程系统的渐进稳定平衡点.在适当的假设条件下,证明了所给出的算法具有二次收敛速度.数值结果表明了此算法的有效性.     

求解非线性方程组的一种新的全局收敛的Levenberg-Marquardt算法

本文提出了求解非线性方程组的一种新的全局收敛的Levenberg-Marquardt算法,即μk=ακ(θ||F_k|| (1-θ)||J_k~TF_k||),θ∈[0,1],其中ακ利用信赖域技巧来修正.在不必假设雅可比矩阵非奇异的局部误差界条件下,证明了该算法是全局收敛和局部二次收敛的.数值试验表明该算法能有效地求解奇异非线性方程组问题.     

A distributionally robust optimization model for batch nonlinear switched time-delay system considering uncertain output measurements

In this paper, we consider a nonlinear switched time-delay (NSTD) system with unknown switching times and unknown system parameters, where the output measurement is uncertain. This system is the underling dynamical system for the batch process of glycerol bioconversion to 1,3-propanediol induced by Klebsiella pneumoniae. The uncertain output measurement is regarded as a stochastic vector (whose components are stochastic variables) and the only information about its distribution is the first-order moment. The objective of this paper is to identify the unknown quantities of the NSTD system. For this, a distributionally robust optimization problem (a bi-level optimization problem) governed by the NSTD system is proposed, where the relative error under the environment of uncertain output measurements is involved in the cost functional. The bi-level optimization problem is transformed into a single-level optimization problem with non-smooth term through the application of duality theory in probability space. By applying the smoothing technique, the non-smooth term is approximated by a smooth term and the convergence of the approximation is established. Then, the gradients of the cost functional with respect to switching times and system parameters are derived. A hybrid optimization algorithm is developed to solve the transformed problem. Finally, we verify the obtained switching times and system parameters, as well as the effectiveness of the proposed algorithm, by solving this distributionally robust optimization problem.     

A Hybrid Smoothing-Nonsmooth Newton-Type Algorithm Yielding an Exact Solution of the $P_0$-LCP

We propose a hybrid smoothing-nonsmooth Newton-type algorithm for solving the P0 linear complementarity problem (P0-LCP) based on the techniques used in the non-smooth Newton method and smoothing Newton method. Under some assumptions, the proposed algorithm can find an exact solution of P0-LCP in finite steps. Preliminary numerical results indicate that the proposed algorithm is promising.     

Level set based multi-scale methods for large deformation contact problems

We consider the numerical simulation of contact problems in elasticity with large deformations. The non-penetration condition is described by means of a signed distance function to the obstacle's boundary. Techniques from level set methods allow for an appropriate numerical approximation of the signed distance function preserving its non-smooth character. The emerging non-convex optimization problem subject to non-smooth inequality constraints is solved by a non-smooth multiscale SQP method in combination with a non-smooth multigrid method as interior solver. Several examples in three space dimensions including applications in biomechanics illustrate the capability of our methods.     

A periodically forced,piecewise linear system. Part I: Local singularity and grazing bifurcation

A methodology for the local singularity of non-smooth dynamical systems is systematically presented in this paper, and a periodically forced, piecewise linear system is investigated as a sample problem to demonstrate the methodology. The sliding dynamics along the separation boundary are investigated through the differential inclusion theory. For this sample problem, a perturbation method is introduced to determine the singularity of the sliding dynamics on the separation boundary. The criteria for grazing bifurcation are presented mathematically and numerically. The grazing flows are illustrated numerically. This methodology can be very easily applied to predict grazing motions in other non-smooth dynamical systems. The fragmentation of the strange attractors of chaotic motion will be presented in the second part of this work.     

A novel algorithm for area traffic capacity control with elastic travel demands

A non-linear area traffic control system with limited capacity is considered in this paper. Optimal signal settings and link capacity expansions can be determined while trip distribution and network flow are in equilibrium. This problem can be formulated as a non-linear mathematical program with equilibrium constraints. For the objective function a non-linear constrained optimization program for signal settings and link capacity expansion is determined. For the constraint set the elastic user equilibrium traffic assignment obeying Wardrop’s first principle can be formulated as a variational inequality. Since the constrained optimization problem is non-convex, only local optima can be obtained. In this paper, a novel algorithm using a non-smooth trust region approach is proposed. Numerical tests are performed using a real data city network and various example test networks in which the effectiveness and robustness of the proposed method are confirmed as compared to other well-known solution methods.     

Chaotic attractor generation and critical value analysis via switching approach

The generation of novel chaotic funnel-shaped attractors is introduced and the analysis of related critical values is given with a proposed switching method in this paper. The underlying mechanism involves a simple three-dimensional switched system and a hysteretically switching signal. Moreover, theoretic analysis is carried out to study the attractor generation and the corresponding critical values by fully utilizing the specific structure of the non-smooth system. Based on carefully derivation, the critical values and related stability regions of the created attractors are estimated explicitly, which is usually impossible for general non-smooth dynamics. In addition, it is demonstrated by simulation that various attractor patterns are generated conveniently by adjusting suitable system parameters.     

Approximation of the eigenvalues of a fourth order differential equation with non-smooth coefficients

The eigenvalues of a fourth order, generalized eigenvalue problem in one dimension, with non-smooth coefficients are approximated by a finite element method, introduced in an earlier work by the author and A. Lutoborski, in the context of a similar source problem with non-smooth coefficients. Error estimates for the approximate eigenvalues and eigenvectors are obtained, showing a better performance of this method, when applied to eigenvalue approximation, compared to a standard finite element method with arbitrary mesh.     

Non-Smooth Optimization for Robust Control of Infinite-Dimensional Systems

We use a non-smooth trust-region method for H _∞-control of infinite-dimensional systems. Our method applies in particular to distributed and boundary control of partial differential equations. It is computationally attractive as it avoids the use of system reduction or identification. For illustration the method is applied to control a reaction-convection-diffusion system, a Van de Vusse reactor, and to a cavity flow control problem.     

D’Alembert function for exact non-smooth modal analysis of the bar in unilateral contact

Non-smooth modal analysis is an extension of modal analysis to non-smooth systems, prone to unilateral contact conditions for instance. The problem of a one-dimensional bar subject to unilateral contact on its boundary has been previously investigated numerically and the corresponding spectrum of vibration could be partially explored. In the present work, the non-smooth modal analysis of the above system is reformulated as a set of functional equations through the use of both d’Alembert solution to the wave equation and the method of steps for Neutral Delay Differential Equations. The system features a strong internal resonance condition and it is established that irrational and rational periods of vibration should be carefully distinguished. For irrational periods, it was previously proven that the displacement field of the non-smooth modes of vibration is characterized with piecewise-linear functions in space and time and such a motion is unique for a prescribed energy. However, for rational periods, which are the subject of this work, new periodic solutions are found analytically. Findings consist of families of iso-periodic solutions with piecewise-smooth displacement fields in space and time and continua of piecewise-smooth periodic solutions of the same energy and frequency.