A VU-decomposition method for a second-order cone programming problem

A vu-decomposition method for solving a second-order cone problem is presented in this paper. It is first transformed into a nonlinear programming problem. Then, the structure of the Clarke subdifferential corresponding to the penalty function and some results of itsvu-decomposition are given. Under a certain condition, a twice continuously differentiable trajectory is computed to produce a second-order expansion of the objective function. A conceptual algorithm for solving this problem with a superlinear convergence rate is given.

Anti-interference beam pattern design based on second-order cone programming optimization

When signal-to-interference ratio is low, the energy of strong interference leaked from the side lobe of beam pattern will infect the detection of weak target. Therefore, the beam pattern needs to be optimized. The existing Dolph-Chebyshev weighting method can get the lowest side lobe level under given main lobe width, but for the other non-uniform circular array and nonlinear array, the low side lobe pattern needs to be designed specially. The second order cone programming optimization （SOCP） algorithm proposed in the paper transforms the optimization of the beam pattern into a standard convex optimization problem. Thus there is a paradigm to follow for any array formation, which not only achieves the purpose of Dolph-Chebyshev weighting, but also solves the problem of the increased side lobe when the signal is at end fire direction The simulation proves that the SOCP algorithm can detect the weak target better than the conventional beam forming.

Improved design of reconfigurable frequency response masking filters based on second-order cone programming

In order to improve the design results for the reconfigurable frequency response masking FRM filters an improved design method based on second-order cone programming SOCP is proposed.Unlike traditional methods that separately design the proposed method takes all the desired designing modes into consideration when designing all the subfilters. First an initial solution is obtained by separately designing the subfilters and then the initial solution is updated by iteratively solving a SOCP problem. The proposed method is evaluated on a design example and simulation results demonstrate that jointly designing all the subfilters can obtain significantly lower minimax approximation errors compared to the conventional design method.

Two new predictor-corrector algorithms for second-order cone programming

Based on the ideas of infeasible interior-point methods and predictor-corrector algorithms, two interior-point predictor-corrector algorithms for the second-order cone programming （SOCP） are presented. The two algorithms use the Newton direction and the Euler direction as the predictor directions, respectively. The corrector directions belong to the category of the Alizadeh-Haeberly-Overton （AHO） directions. These algorithms are suitable to the cases of feasible and infeasible interior iterative points. A simpler neighborhood of the central path for the SOCP is proposed, which is the pivotal difference from other interior-point predictor-corrector algorithms. Under some assumptions, the algorithms possess the global, linear, and quadratic convergence. The complexity bound O（rln（εo/ε）） is obtained, where r denotes the number of the second-order cones in the SOCP problem. The numerical results show that the proposed algorithms are effective.

Robust Blind Separation for MIMO Systems against Channel Mismatch Using Second-Order Cone Programming

To improve the deteriorated capacity gain and source recovery performance due to channel mismatch problem,this paper reports a research about blind separation method against channel mismatch in multiple-input multiple-output(MIMO) systems.The channel mismatch problem can be described as a channel with bounded fluctuant errors due to channel distortion or channel estimation errors.The problem of blind signal separation/extraction with channel mismatch is formulated as a cost function of blind source separation(BSS) subject to the second-order cone constraint,which can be called as second-order cone programing optimization problem.Then the resulting cost function is solved by approximate negentropy maximization using quasi-Newton iterative methods for blind separation/extraction source signals.Theoretical analysis demonstrates that the proposed algorithm has low computational complexity and improved performance advantages.Simulation results verify that the capacity gain and bit error rate(BER) performance of the proposed blind separation method is superior to those of the existing methods in MIMO systems with channel mismatch problem.

AN INFEASIBLE-INTERIOR-POINT PREDICTOR-CORRECTOR ALGORITHM FOR THE SECOND-ORDER CONE PROGRAM

A globally convergent infeasible-interior-point predictor-corrector algorithm is presented for the second-order cone programming （SOCP） by using the Alizadeh- Haeberly-Overton （AHO） search direction. This algorithm does not require the feasibility of the initial points and iteration points. Under suitable assumptions, it is shown that the algorithm can find an -approximate solution of an SOCP in at most O（√n ln（ε0/ε）） iterations. The iteration-complexity bound of our algorithm is almost the same as the best known bound of feasible interior point algorithms for the SOCP.

SMOOTHING NEWTON ALGORITHM FOR THE CIRCULAR CONE PROGRAMMING WITH A NONMONOTONE LINE SEARCH

In this paper, we present a nonmonotone smoothing Newton algorithm for solving the circular cone programming(CCP) problem in which a linear function is minimized or maximized over the intersection of an affine space with the circular cone. Based on the relationship between the circular cone and the second-order cone(SOC), we reformulate the CCP problem as the second-order cone problem(SOCP). By extending the nonmonotone line search for unconstrained optimization to the CCP, a nonmonotone smoothing Newton method is proposed for solving the CCP. Under suitable assumptions, the proposed algorithm is shown to be globally and locally quadratically convergent. Some preliminary numerical results indicate the effectiveness of the proposed algorithm for solving the CCP.

A New Method Combining Interior and Exterior Approaches for Linear Programming

In this paper we present a new method combining interior and exterior approaches to solve linear programming problems. With the assumption that a feasible interior solution to the input system is known, this algorithm uses it and appropriate constraints of the system to construct a sequence of the so called station cones whose vertices tend very fast to the solution to be found. The computational experiments show that the number of iterations of the new algorithm is significantly smaller than that of the second phase of the simplex method. Additionally, when the number of variables and constraints of the problem increase, the number of iterations of the new algorithm increase in a slower manner than that of the simplex method.

Polynomial Complexity Bounds of Mehrotra-type Predictor-corrector Algorithms for Linear Programming over Symmetric Cones

We establish polynomial complexity corrector algorithms for linear programming over bounds of the Mehrotra-type predictor- symmetric cones. We first slightly modify the maximum step size in the predictor step of the safeguard based Mehrotra-type algorithm for linear programming, that was proposed by Salahi et al. Then, using the machinery of Euclidean Jordan algebras, we extend the modified algorithm to symmetric cones. Based on the Nesterov-Todd direction, we obtain O（r log ε1） iteration complexity bound of this algorithm, where r is the rank of the Jordan algebras and ε is the required precision. We also present a new variant of Mehrotra-type algorithm using a new adaptive updating scheme of centering parameter and show that this algorithm enjoys the same order of complexity bound as the safeguard algorithm. We illustrate the numerical behaviour of the methods on some small examples.

Station Cone Algorithm for Linear Programming

Recently we have proposed anew method combininginterior and exterior approaches to solve linear programming problems. This method uses an interior point, and from there connected to the vertex of the so called station cone which is also a solution of the dual problem. This allows us to determine the entering vector and the new station cone. Here in this paper, we present a new modified algorithm for the case, when at each iteration we determine a new interior point. The new building interior point moves toward the optimal vertex. Thanks to the shortened from both inside and outside, the new version allows to find quicker the optimal solution. The computational experiments show that the number of iterations of the new modified algorithm is significantly smaller than that of the second phase of the dual simplex method.

A stable implicit nodal integration-based particle finite element method(N-PFEM)for modelling saturated soil dynamics

In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.

Day-ahead Optimization Schedule for Gas-electric Integrated Energy System Based on Second-order Cone Programming

This paper proposes an optimal day-ahead opti-mization schedule for gas-electric integrated energy system(IES)considering the bi-directional energy flow.The hourly topology of electric power system(EPS),natural gas system(NGS),energy hubs(EH)integrated power to gas(P2G)unit,are modeled to minimize the day-ahead operation cost of IES.Then,a second-order cone programming(SOCP)method is utilized to solve the optimization problem,which is actually a mixed integer nonconvex and nonlinear programming issue.Besides,cutting planes are added to ensure the exactness of the global optimal solution.Finally,simulation results demonstrate that the proposed optimization schedule can provide a safe,effective and economical day-ahead scheduling scheme for gas-electric IES.

Optimal design and verification of temporal and spatial filters using second-order cone programming approach

Temporal filters and spatial filters are widely used in many areas of signal processing. A number of optimal design criteria to these problems are available in the literature. Various computational techniques are also presented to optimize these criteria chosen. There are many drawbacks in these methods. In this paper, we introduce a unified framework for optimal design of temporal and spatial filters. Most of the optimal design problems of FIR filters and beamformers are included in the framework. It is shown that all the design problems can be reformulated as convex optimization form as the second-order cone programming （SOCP） and solved efficiently via the well-established interior point methods. The main advantage of our SOCP approach as compared with earlier approaches is that it can include most of the existing methods as its special cases, which leads to more flexible designs. Furthermore, the SOCP approach can optimize multiple required performance measures, which is the drawback of earlier approaches. The SOCP approach is also developed to optimally design temporal and spatial two-dimensional filter and spatial matrix filter. Numerical results demonstrate the effectiveness of the proposed approach.

Nonsingularity in second-order cone programming via the smoothing metric projector

Based on the differential properties of the smoothing metric projector onto the second-order cone,we prove that,for a locally optimal solution to a nonlinear second-order cone programming problem,the nonsingularity of the Clarke's generalized Jacobian of the smoothing Karush-Kuhn-Tucker system,constructed by the smoothing metric projector,is equivalent to the strong second-order sufficient condition and constraint nondegeneracy,which is in turn equivalent to the strong regularity of the Karush-Kuhn-Tucker point.Moreover,this nonsingularity property guarantees the quadratic convergence of the corresponding smoothing Newton method for solving a Karush-Kuhn-Tucker point.Interestingly,the analysis does not need the strict complementarity condition.

Exact Computable Representation of Some Second-Order Cone Constrained Quadratic Programming Problems

Solving the quadratically constrained quadratic programming(QCQP)problem is in general NP-hard.Only a few subclasses of the QCQP problem are known to be polynomial-time solvable.Recently,the QCQP problem with a nonconvex quadratic objective function over one ball and two parallel linear constraints is proven to have an exact computable representation,which reformulates the original problem as a linear semidefinite program with additional linear and second-order cone constraints.In this paper,we provide exact computable representations for some more subclasses of the QCQP problem,in particular,the subclass with one secondorder cone constraint and two special linear constraints.

Unified convergence analysis of a second-order method of multipliers for nonlinear conic programming

In this paper,we accomplish the unified convergence analysis of a second-order method of multipliers(i.e.,a second-order augmented Lagrangian method)for solving the conventional nonlinear conic optimization problems.Specifically,the algorithm that we investigate incorporates a specially designed nonsmooth(generalized)Newton step to furnish a second-order update rule for the multipliers.We first show in a unified fashion that under a few abstract assumptions,the proposed method is locally convergent and possesses a(nonasymptotic)superlinear convergence rate,even though the penalty parameter is fixed and/or the strict complementarity fails.Subsequently,we demonstrate that for the three typical scenarios,i.e.,the classic nonlinear programming,the nonlinear second-order cone programming and the nonlinear semidefinite programming,these abstract assumptions are nothing but exactly the implications of the iconic sufficient conditions that are assumed for establishing the Q-linear convergence rates of the method of multipliers without assuming the strict complementarity.

The use of the node-based smoothed finite element method to estimate static and seismic bearing capacities of shallow strip footings

The node-based smoothed finite element method(NS-FEM)is shortly presented for calculations of the static and seismic bearing capacities of shallow strip footings.A series of computations has been performed to assess variations in seismic bearing capacity factors with both horizontal and vertical seismic accelerations.Numerical results obtained agree very well with those using the slip-line method,revealing that the magnitude of the seismic bearing capacity is highly dependent upon the combinations of various directions of both components of the seismic acceleration.An upward vertical seismic acceleration reduces the seismic bearing capacity compared to the downward vertical seismic acceleration in calculations.In addition,particular emphasis is placed on a separate estimation of the effects of soil and superstructure inertia on each seismic bearing capacity component.While the effect of inertia forces arising in the soil on the seismic bearing capacity is non-trivial,and the superstructure inertia is the major contributor to reductions in the seismic bearing capacity.Both tables and charts are given for practical application to the seismic design of the foundations.

On Iteration Complexity of a First-Order Primal-Dual Method for Nonlinear Convex Cone Programming

Nonlinear convex cone programming(NCCP)models have found many practical applications.In this paper,we introduce a flexible first-order primal-dual algorithm,called the variant auxiliary problem principle(VAPP),for solving NCCP problems when the objective function and constraints are convex but may be nonsmooth.At each iteration,VAPP generates a nonlinear approximation of the primal augmented Lagrangian model.The approximation incorporates both linearization and a distance-like proximal term,and then the iterations of VAPP are shown to possess a decomposition property for NCCP.Motivated by recent applications in big data analytics,there has been a growing interest in the convergence rate analysis of algorithms with parallel computing capabilities for large scale optimization problems.We establish O(1/t)convergence rate towards primal optimality,feasibility and dual optimality.By adaptively setting parameters at different iterations,we show an O(1/t2)rate for the strongly convex case.Finally,we discuss some issues in the implementation of VAPP.

A Modified and Simplified Full Nesterov–Todd Step O(N)Infeasible Interior-Point Method for Second-Order Cone Optimization

We present a modified and simplified version of an infeasible interior-point method for second-order cone optimization published in 2013(Zangiabadi et al.in J Optim Theory Appl,2013).In the earlier version,each iteration consisted of one socalled feasibility step and a few centering steps.Here,each iteration consists of only a feasibility step.Thus,the new algorithm improves the number of iterations and the improvement is due to a lemma which gives an upper bound for the proximity after the feasibility step.The complexity result coincides with the best-known iteration bound for infeasible interior-point methods.

极限下限分析的区域光滑径向点插值法

极限分析的高效数值计算方法在结构设计和安全评定中具有非常重要的作用.为了更加有效地求解极限分析问题,将区域光滑径向点插值法与二阶锥规划相结合,提出了理想弹塑性结构极限下限分析的一种新方法.将问题域离散为简单的三角形背景单元,每个单元进一步划分成若干个光滑域.为了将复杂的域积分转化为简单的边界积分,并且避免计算形函数的导数,采用广义梯度光滑技术对每个光滑域进行应变光滑处理.由于径向点插值法构造的形函数满足Kronecker delta性质,本质边界条件可以直接施加.依据下限定理,在满足以等效积分弱形式表达的自平衡应力场平衡条件的基础上,用二阶锥规划成功构建了极限分析下限法的计算模型,从而可方便地通过基于原始-对偶内点法的数学规划求解器MOSEK直接求解该问题.数值算例结果表明,本文所提方法有效地克服了维数障碍问题,具有较高的计算精度,并且计算结果对网格畸变十分不敏感.