A multifractal model is developed to connect the Lagrangian multifractal dimensions with their Eulerian counterparts. We propose that the characteristic time scale of a Lagrangian quantity should be the Lagrangian tim...A multifractal model is developed to connect the Lagrangian multifractal dimensions with their Eulerian counterparts. We propose that the characteristic time scale of a Lagrangian quantity should be the Lagrangian time scale, and it should not be the Eulerian time scale which was widely used in previous studies on Lagrangian statistics. Using the present model, we can obtain the scaling exponents of Lagrangian velocity structure functions from the existing data or models of scaling exponents of Eulerian velocity structure functions. This model is validated by comparing its prediction with the results of experiments, direct numerical simulations, and the previous theoretical models. The comparison shows that the proposed model can better predict the scaling exponents of Lagrangian velocity structure functions, especially for orders larger than 6.展开更多
In this paper, a modified version of the Classical Lagrange Multiplier method is developed for convex quadratic optimization problems. The method, which is evolved from the first order derivative test for optimality o...In this paper, a modified version of the Classical Lagrange Multiplier method is developed for convex quadratic optimization problems. The method, which is evolved from the first order derivative test for optimality of the Lagrangian function with respect to the primary variables of the problem, decomposes the solution process into two independent ones, in which the primary variables are solved for independently, and then the secondary variables, which are the Lagrange multipliers, are solved for, afterward. This is an innovation that leads to solving independently two simpler systems of equations involving the primary variables only, on one hand, and the secondary ones on the other. Solutions obtained for small sized problems (as preliminary test of the method) demonstrate that the new method is generally effective in producing the required solutions.展开更多
In this paper, based on the theorem of the high-order velocity energy, integration and variation principle, the high-order Hamilton's principle of general holonomic systems is given. Then, three-order Lagrangian equa...In this paper, based on the theorem of the high-order velocity energy, integration and variation principle, the high-order Hamilton's principle of general holonomic systems is given. Then, three-order Lagrangian equations and four-order Lagrangian equations are obtained from the high-order Hamilton's principle. Finally, the Hamilton's principle of high-order Lagrangian function is given.展开更多
An exact augmented Lagrangian function for the nonlinear nonconvex programming problems with inequality constraints was discussed. Under suitable hypotheses, the relationship was established between the local unconstr...An exact augmented Lagrangian function for the nonlinear nonconvex programming problems with inequality constraints was discussed. Under suitable hypotheses, the relationship was established between the local unconstrained minimizers of the augmented Lagrangian function on the space of problem variables and the local minimizers of the original constrained problem. Furthermore, under some assumptions, the relationship was also established between the global solutions of the augmented Lagrangian function on some compact subset of the space of problem variables and the global solutions of the constrained problem. Therefore, f^om the theoretical point of view, a solution of the inequality constrained problem and the corresponding values of the Lagrange multipliers can be found by the well-known method of multipliers which resort to the unconstrained minimization of the augmented Lagrangian function presented.展开更多
A fuzzy adaptive control method is proposed for a flexible robot manipulator. Due to the structure characteristics of the flexible manipulator, the vibration modes must be controlled to realize the high-precision tip ...A fuzzy adaptive control method is proposed for a flexible robot manipulator. Due to the structure characteristics of the flexible manipulator, the vibration modes must be controlled to realize the high-precision tip position. The Lagrangian principle is utilized to model the dynamic function of the single-degree flexible manipulator incorporating the assumed modes method. Simulation results of the fuzzy adaptive control method in the location control and the trajectory tracking with different tip disturbances are presented and compared with the results of the classic PD control. It shows that the controller can obtain the stable and robust performance.展开更多
A new method for estimating the n (50 or 100) -year return-period waveheight, namely, the extreme waveheight expected to occur in n years, is presented on the basis of the maximum entropy principle. The main p...A new method for estimating the n (50 or 100) -year return-period waveheight, namely, the extreme waveheight expected to occur in n years, is presented on the basis of the maximum entropy principle. The main points of the method are as follows: (1) based on the Hamiltonian principle, a maximum entropy probability density function for the extreme waveheight H, f(H)=αHγ e -βH4 is derived from a Lagrangian function subject to some necessary and rational constraints; (2) the parameters α, β, and γ in the function are expressed in terms of the mean , variance V= (H-)2 and bias B= (H-)3 ; and (3) with , V and B estimated from observed data, the n -year return-period wave height H n is computed in accordance with the formula 11-F(H n)=n , where F(H n) is defined as F(H n)=∫ H n 0f(H) d H. Examples of estimating the 50 and 100-year return period waveheights by the present method and by some currently used method from observed data acquired from two hydrographic stations are given. A comparison of the estimated results shows that the present method is superior to the others.展开更多
Repeater optimization is the key for SOC (System on Chip) interconnect delay design. This paper proposes a novel optimal model for minimizing power and area overhead of repeaters while meeting the target performance...Repeater optimization is the key for SOC (System on Chip) interconnect delay design. This paper proposes a novel optimal model for minimizing power and area overhead of repeaters while meeting the target performance of on-chip interconnect lines. It also presents Lagrangian function to find the number of repeaters and their sizes required for minimizing area and power overhead with target delay constraint. Based on the 65 nanometre CMOS technology, the computed results of the intermediate and global lines show that the proposed model can significantly reduce area and power of interconnected lines, and the better performance will be achieved with the longer line. The results compared with the reference paper demonstrate the validity of this model. It can be integrated into repeater design methodology and CAD (computer aided design) tool for interconnect planning in nanometre SOC.展开更多
With the help of the time-dependent Lagrangian for a damped harmonic oscillator, the quantization of mesoscopic RLC circuit in the context of a number-phase quantization scheme is realized and the corresponding Hamilt...With the help of the time-dependent Lagrangian for a damped harmonic oscillator, the quantization of mesoscopic RLC circuit in the context of a number-phase quantization scheme is realized and the corresponding Hamiltonian operator is obtained. Then the evolution of the charge number and phase difference across the capacity are obtained. It is shown that the number-phase analysis is useful to tackle the quantization of some mesoscopic circuits and dynamical equations of the corresponding operators.展开更多
In this paper,general interpolating isogeometric boundary node method(IIBNM)and isogeometric boundary element method(IBEM)based on parameter space are proposed for 2D elasticity problems.In both methods,the integral c...In this paper,general interpolating isogeometric boundary node method(IIBNM)and isogeometric boundary element method(IBEM)based on parameter space are proposed for 2D elasticity problems.In both methods,the integral cells and elements are defined in parameter space,which can reproduce the geometry exactly at all the stages.In IIBNM,the improved interpolating moving leastsquare method(IIMLS)is applied for field approximation and the shape functions have the delta function property.The Lagrangian basis functions are used for field approximation in IBEM.Thus,the boundary conditions can be imposed directly in both methods.The shape functions are defined in 1D parameter space and no curve length needs to be computed.Besides,most methods for the treatment of the singular integrals in the boundary element method can be applied in IIBNM and IBEM directly.Numerical examples have demonstrated the accuracy of the proposed methods.展开更多
Based on the covariant Lagrangian function and Euler-Lagrange equation,a set of classical fluid equations for strong EM wave-spin plasma interaction is derived.Analysis shows that the relativistic effects may affect t...Based on the covariant Lagrangian function and Euler-Lagrange equation,a set of classical fluid equations for strong EM wave-spin plasma interaction is derived.Analysis shows that the relativistic effects may affect the interaction processes by three factors:the relativistic factor,the time component of four-spin,and the velocity-field coupling.This set of equations can be used to discuss the collective spin effects of relativistic electrons in classical regime,such as astrophysics,high-energy laser-plasma systems and so on.As an example,the spin induced ponderomotive force in the interaction of strong EM wave and magnetized plasma is investigated.Results show that the time component of four-spin,which approaches to zero in nonrelativistic situations,can increase the spin-ponderomotive force obviously in relativistic situation.展开更多
Self-adjoint theorem is introduced to match the corresponding functional of the given differential equations,and then Noether's theorem is used to determine the extended conservation laws of the original equations. F...Self-adjoint theorem is introduced to match the corresponding functional of the given differential equations,and then Noether's theorem is used to determine the extended conservation laws of the original equations. Finally, as the application of the method, the conservation laws of Drinfel'd-Sokolov-Wilson equation and Benjamin-Bona-Mahony equation are constructed.展开更多
A kind of direct methods is presented for the solution of optimal control problems with state constraints. These methods are sequential quadratic programming methods. At every iteration a quadratic programming which i...A kind of direct methods is presented for the solution of optimal control problems with state constraints. These methods are sequential quadratic programming methods. At every iteration a quadratic programming which is obtained by quadratic approximation to Lagrangian function and linear approximations to constraints is solved to get a search direction for a merit function. The merit function is formulated by augmenting the Lagrangian function with a penalty term. A line search is carried out along the search direction to determine a step length such that the merit function is decreased. The methods presented in this paper include continuous sequential quadratic programming methods and discreate sequential quadratic programming methods.展开更多
The paper concerns the problem how to purchase the reinsurance in order to make the insurer and the reinsurance company's total risk to be least under the expected value principle. When the insurer and reinsurance co...The paper concerns the problem how to purchase the reinsurance in order to make the insurer and the reinsurance company's total risk to be least under the expected value principle. When the insurer and reinsurance company take arbitrary risk measures, sufficient con- ditions for optimality of reinsurance contract are given within the restricted class of admissible contracts. Further, the explicit forms of optimal reinsurance contract under several special risk measures are given, and the method to decide parameters as well.展开更多
Non-intrusive load monitoring is a technique for monitoring the operating conditions of electrical appliances by collecting the aggregated electrical information at the household power inlet.Despite several studies on...Non-intrusive load monitoring is a technique for monitoring the operating conditions of electrical appliances by collecting the aggregated electrical information at the household power inlet.Despite several studies on the mining of unique load characteristics,few studies have extensively considered the high computational burden and sample training.Based on lowfrequency sampling data,a non-intrusive load monitoring algorithm utilizing the graph total variation(GTV)is proposed in this study.The algorithm can effectively depict the load state without the need for prior training.First,the combined Kmeans clustering algorithm and graph signals are used to build concise and accurate graph structures as load models.The GTV representing the internal structure of the graph signal is introduced as the optimization model and solved using the augmented Lagrangian iterative algorithm.The introduction of the difference operator reduces the computing cost and addresses the inaccurate reconstruction of the graph signal.With low-frequency sampling data,the algorithm only requires a little prior data and no training,thereby reducing the computing cost.Experiments conducted using the reference energy disaggregation dataset and almanac of minutely power dataset demonstrated the stable superiority of the algorithm and its low computational burden.展开更多
An augmented Lagrangian trust region method with a bi=object strategy is proposed for solving nonlinear equality constrained optimization, which falls in between penalty-type methods and penalty-free ones. At each ite...An augmented Lagrangian trust region method with a bi=object strategy is proposed for solving nonlinear equality constrained optimization, which falls in between penalty-type methods and penalty-free ones. At each iteration, a trial step is computed by minimizing a quadratic approximation model to the augmented Lagrangian function within a trust region. The model is a standard trust region subproblem for unconstrained optimization and hence can efficiently be solved by many existing methods. To choose the penalty parameter, an auxiliary trust region subproblem is introduced related to the constraint violation. It turns out that the penalty parameter need not be monotonically increasing and will not tend to infinity. A bi-object strategy, which is related to the objective function and the measure of constraint violation, is utilized to decide whether the trial step will be accepted or not. Global convergence of the method is established under mild assumptions. Numerical experiments are made, which illustrate the efficiency of the algorithm on various difficult situations.展开更多
This paper proposes a nonmonotone line search filter method with reduced Hessian updating for solving nonlinear equality constrained optimization.In order to deal with large scale problems,a reduced Hessian matrix is ...This paper proposes a nonmonotone line search filter method with reduced Hessian updating for solving nonlinear equality constrained optimization.In order to deal with large scale problems,a reduced Hessian matrix is approximated by BFGS updates.The new method assures global convergence without using a merit function.By Lagrangian function in the filter and nonmonotone scheme,the authors prove that the method can overcome Maratos effect without using second order correction step so that the locally superlinear convergence is achieved.The primary numerical experiments are reported to show effectiveness of the proposed algorithm.展开更多
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 sol...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.展开更多
Firstly,this paper proposes a generalized log-determinant optimization model with the purpose of estimating the high-dimensional sparse inverse covariance matrices.Under the normality assumption,the zero components in...Firstly,this paper proposes a generalized log-determinant optimization model with the purpose of estimating the high-dimensional sparse inverse covariance matrices.Under the normality assumption,the zero components in the inverse covariance matrices represent the conditional independence between pairs of variables given all the other variables.The generalized model considered in this study,because of the setting of the eigenvalue bounded constraints,covers a large number of existing estimators as special cases.Secondly,rather than directly tracking the challenging optimization problem,this paper uses a couple of alternating direction methods of multipliers(ADMM)to solve its dual model where 5 separable structures are contained.The first implemented algorithm is based on a single Gauss–Seidel iteration,but it does not necessarily converge theoretically.In contrast,the second algorithm employs the symmetric Gauss–Seidel(sGS)based ADMM which is equivalent to the 2-block iterative scheme from the latest sGS decomposition theorem.Finally,we do numerical simulations using the synthetic data and the real data set which show that both algorithms are very effective in estimating high-dimensional sparse inverse covariance matrix.展开更多
In this paper we consider the Cauchy problem for a semi-linear system of wave equations with Hamilton structure. We prove the existence of global smooth solution of the system for subcritical case by using conservatio...In this paper we consider the Cauchy problem for a semi-linear system of wave equations with Hamilton structure. We prove the existence of global smooth solution of the system for subcritical case by using conservation of energy and Strichartz's estimate. On the basis of Morawetz-Pohozev identity, we obtain the same result for the critical case.展开更多
基金supported by the National Natural Science Foundation of China(11072247,11021262,and 11232011)National Natural Science Associate Foundation of China(NSAF)(U1230126)973 program of China(2013CB834100)
文摘A multifractal model is developed to connect the Lagrangian multifractal dimensions with their Eulerian counterparts. We propose that the characteristic time scale of a Lagrangian quantity should be the Lagrangian time scale, and it should not be the Eulerian time scale which was widely used in previous studies on Lagrangian statistics. Using the present model, we can obtain the scaling exponents of Lagrangian velocity structure functions from the existing data or models of scaling exponents of Eulerian velocity structure functions. This model is validated by comparing its prediction with the results of experiments, direct numerical simulations, and the previous theoretical models. The comparison shows that the proposed model can better predict the scaling exponents of Lagrangian velocity structure functions, especially for orders larger than 6.
文摘In this paper, a modified version of the Classical Lagrange Multiplier method is developed for convex quadratic optimization problems. The method, which is evolved from the first order derivative test for optimality of the Lagrangian function with respect to the primary variables of the problem, decomposes the solution process into two independent ones, in which the primary variables are solved for independently, and then the secondary variables, which are the Lagrange multipliers, are solved for, afterward. This is an innovation that leads to solving independently two simpler systems of equations involving the primary variables only, on one hand, and the secondary ones on the other. Solutions obtained for small sized problems (as preliminary test of the method) demonstrate that the new method is generally effective in producing the required solutions.
基金the Natural Science Foundation of Jiangxi Provincethe Foundation of Education Department of Jiangxi Province under Grant No.[2007]136
文摘In this paper, based on the theorem of the high-order velocity energy, integration and variation principle, the high-order Hamilton's principle of general holonomic systems is given. Then, three-order Lagrangian equations and four-order Lagrangian equations are obtained from the high-order Hamilton's principle. Finally, the Hamilton's principle of high-order Lagrangian function is given.
文摘An exact augmented Lagrangian function for the nonlinear nonconvex programming problems with inequality constraints was discussed. Under suitable hypotheses, the relationship was established between the local unconstrained minimizers of the augmented Lagrangian function on the space of problem variables and the local minimizers of the original constrained problem. Furthermore, under some assumptions, the relationship was also established between the global solutions of the augmented Lagrangian function on some compact subset of the space of problem variables and the global solutions of the constrained problem. Therefore, f^om the theoretical point of view, a solution of the inequality constrained problem and the corresponding values of the Lagrange multipliers can be found by the well-known method of multipliers which resort to the unconstrained minimization of the augmented Lagrangian function presented.
文摘A fuzzy adaptive control method is proposed for a flexible robot manipulator. Due to the structure characteristics of the flexible manipulator, the vibration modes must be controlled to realize the high-precision tip position. The Lagrangian principle is utilized to model the dynamic function of the single-degree flexible manipulator incorporating the assumed modes method. Simulation results of the fuzzy adaptive control method in the location control and the trajectory tracking with different tip disturbances are presented and compared with the results of the classic PD control. It shows that the controller can obtain the stable and robust performance.
基金ThisworkisfinanciallysupportedbythePh.D.FoundationoftheMinistryoftheEducationofChina (No .2 0 0 0 4 2 30 8)
文摘A new method for estimating the n (50 or 100) -year return-period waveheight, namely, the extreme waveheight expected to occur in n years, is presented on the basis of the maximum entropy principle. The main points of the method are as follows: (1) based on the Hamiltonian principle, a maximum entropy probability density function for the extreme waveheight H, f(H)=αHγ e -βH4 is derived from a Lagrangian function subject to some necessary and rational constraints; (2) the parameters α, β, and γ in the function are expressed in terms of the mean , variance V= (H-)2 and bias B= (H-)3 ; and (3) with , V and B estimated from observed data, the n -year return-period wave height H n is computed in accordance with the formula 11-F(H n)=n , where F(H n) is defined as F(H n)=∫ H n 0f(H) d H. Examples of estimating the 50 and 100-year return period waveheights by the present method and by some currently used method from observed data acquired from two hydrographic stations are given. A comparison of the estimated results shows that the present method is superior to the others.
基金supported by the National Natural Science Foundation of the China (Grant Nos 60676009 and 60776034)the Doctor Foundation of Ministry of Education of China (Grant No 20050701015)the National Outstanding Young Scientist Foundation of China (Grant No 60725415)
文摘Repeater optimization is the key for SOC (System on Chip) interconnect delay design. This paper proposes a novel optimal model for minimizing power and area overhead of repeaters while meeting the target performance of on-chip interconnect lines. It also presents Lagrangian function to find the number of repeaters and their sizes required for minimizing area and power overhead with target delay constraint. Based on the 65 nanometre CMOS technology, the computed results of the intermediate and global lines show that the proposed model can significantly reduce area and power of interconnected lines, and the better performance will be achieved with the longer line. The results compared with the reference paper demonstrate the validity of this model. It can be integrated into repeater design methodology and CAD (computer aided design) tool for interconnect planning in nanometre SOC.
文摘With the help of the time-dependent Lagrangian for a damped harmonic oscillator, the quantization of mesoscopic RLC circuit in the context of a number-phase quantization scheme is realized and the corresponding Hamiltonian operator is obtained. Then the evolution of the charge number and phase difference across the capacity are obtained. It is shown that the number-phase analysis is useful to tackle the quantization of some mesoscopic circuits and dynamical equations of the corresponding operators.
基金The research for this paper was supported by(1)the National Natural Science Foundation of China(Grants Nos.51708429,51708428)the Open Projects Foundation(Grant No.2017-04-GF)of State Key Laboratory for Health and Safety of Bridge Structures+1 种基金Wuhan Institute of Technology Science Found(Grant No.K201734)the science and technology projects of Wuhan Urban and Rural Construction Bureau(Grants Nos.201831,201919).
文摘In this paper,general interpolating isogeometric boundary node method(IIBNM)and isogeometric boundary element method(IBEM)based on parameter space are proposed for 2D elasticity problems.In both methods,the integral cells and elements are defined in parameter space,which can reproduce the geometry exactly at all the stages.In IIBNM,the improved interpolating moving leastsquare method(IIMLS)is applied for field approximation and the shape functions have the delta function property.The Lagrangian basis functions are used for field approximation in IBEM.Thus,the boundary conditions can be imposed directly in both methods.The shape functions are defined in 1D parameter space and no curve length needs to be computed.Besides,most methods for the treatment of the singular integrals in the boundary element method can be applied in IIBNM and IBEM directly.Numerical examples have demonstrated the accuracy of the proposed methods.
基金supported by National Natural Science Foundation of China(No.12065011)Science and Technology Research Project of Jiangxi Provincial Department of Education(No.GJJ170642)。
文摘Based on the covariant Lagrangian function and Euler-Lagrange equation,a set of classical fluid equations for strong EM wave-spin plasma interaction is derived.Analysis shows that the relativistic effects may affect the interaction processes by three factors:the relativistic factor,the time component of four-spin,and the velocity-field coupling.This set of equations can be used to discuss the collective spin effects of relativistic electrons in classical regime,such as astrophysics,high-energy laser-plasma systems and so on.As an example,the spin induced ponderomotive force in the interaction of strong EM wave and magnetized plasma is investigated.Results show that the time component of four-spin,which approaches to zero in nonrelativistic situations,can increase the spin-ponderomotive force obviously in relativistic situation.
基金Supported by "Math + X" Fund of Dalian University of Technology, Science Foundation of Dalian University of Technology under Grant No. SFDUT0808the National Key Basic Research Development of China under Grant No. 2004CB318000
文摘Self-adjoint theorem is introduced to match the corresponding functional of the given differential equations,and then Noether's theorem is used to determine the extended conservation laws of the original equations. Finally, as the application of the method, the conservation laws of Drinfel'd-Sokolov-Wilson equation and Benjamin-Bona-Mahony equation are constructed.
文摘A kind of direct methods is presented for the solution of optimal control problems with state constraints. These methods are sequential quadratic programming methods. At every iteration a quadratic programming which is obtained by quadratic approximation to Lagrangian function and linear approximations to constraints is solved to get a search direction for a merit function. The merit function is formulated by augmenting the Lagrangian function with a penalty term. A line search is carried out along the search direction to determine a step length such that the merit function is decreased. The methods presented in this paper include continuous sequential quadratic programming methods and discreate sequential quadratic programming methods.
文摘The paper concerns the problem how to purchase the reinsurance in order to make the insurer and the reinsurance company's total risk to be least under the expected value principle. When the insurer and reinsurance company take arbitrary risk measures, sufficient con- ditions for optimality of reinsurance contract are given within the restricted class of admissible contracts. Further, the explicit forms of optimal reinsurance contract under several special risk measures are given, and the method to decide parameters as well.
基金supported by National Natural Science Foundation of China(No.52107117)。
文摘Non-intrusive load monitoring is a technique for monitoring the operating conditions of electrical appliances by collecting the aggregated electrical information at the household power inlet.Despite several studies on the mining of unique load characteristics,few studies have extensively considered the high computational burden and sample training.Based on lowfrequency sampling data,a non-intrusive load monitoring algorithm utilizing the graph total variation(GTV)is proposed in this study.The algorithm can effectively depict the load state without the need for prior training.First,the combined Kmeans clustering algorithm and graph signals are used to build concise and accurate graph structures as load models.The GTV representing the internal structure of the graph signal is introduced as the optimization model and solved using the augmented Lagrangian iterative algorithm.The introduction of the difference operator reduces the computing cost and addresses the inaccurate reconstruction of the graph signal.With low-frequency sampling data,the algorithm only requires a little prior data and no training,thereby reducing the computing cost.Experiments conducted using the reference energy disaggregation dataset and almanac of minutely power dataset demonstrated the stable superiority of the algorithm and its low computational burden.
文摘An augmented Lagrangian trust region method with a bi=object strategy is proposed for solving nonlinear equality constrained optimization, which falls in between penalty-type methods and penalty-free ones. At each iteration, a trial step is computed by minimizing a quadratic approximation model to the augmented Lagrangian function within a trust region. The model is a standard trust region subproblem for unconstrained optimization and hence can efficiently be solved by many existing methods. To choose the penalty parameter, an auxiliary trust region subproblem is introduced related to the constraint violation. It turns out that the penalty parameter need not be monotonically increasing and will not tend to infinity. A bi-object strategy, which is related to the objective function and the measure of constraint violation, is utilized to decide whether the trial step will be accepted or not. Global convergence of the method is established under mild assumptions. Numerical experiments are made, which illustrate the efficiency of the algorithm on various difficult situations.
基金supported by the National Science Foundation of China under Grant No.10871130the Ph.D Foundation under Grant No.20093127110005+1 种基金the Shanghai Leading Academic Discipline Project under Grant No.S30405the Innovation Program of Shanghai Municipal Education Commission under Grant No.12YZ174
文摘This paper proposes a nonmonotone line search filter method with reduced Hessian updating for solving nonlinear equality constrained optimization.In order to deal with large scale problems,a reduced Hessian matrix is approximated by BFGS updates.The new method assures global convergence without using a merit function.By Lagrangian function in the filter and nonmonotone scheme,the authors prove that the method can overcome Maratos effect without using second order correction step so that the locally superlinear convergence is achieved.The primary numerical experiments are reported to show effectiveness of the proposed algorithm.
基金This research was supported by the National Natural Science Foundation of China(Nos.71471112 and 71871140).
文摘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.
基金the National Natural Science Foundation of China(No.11971149).
文摘Firstly,this paper proposes a generalized log-determinant optimization model with the purpose of estimating the high-dimensional sparse inverse covariance matrices.Under the normality assumption,the zero components in the inverse covariance matrices represent the conditional independence between pairs of variables given all the other variables.The generalized model considered in this study,because of the setting of the eigenvalue bounded constraints,covers a large number of existing estimators as special cases.Secondly,rather than directly tracking the challenging optimization problem,this paper uses a couple of alternating direction methods of multipliers(ADMM)to solve its dual model where 5 separable structures are contained.The first implemented algorithm is based on a single Gauss–Seidel iteration,but it does not necessarily converge theoretically.In contrast,the second algorithm employs the symmetric Gauss–Seidel(sGS)based ADMM which is equivalent to the 2-block iterative scheme from the latest sGS decomposition theorem.Finally,we do numerical simulations using the synthetic data and the real data set which show that both algorithms are very effective in estimating high-dimensional sparse inverse covariance matrix.
文摘In this paper we consider the Cauchy problem for a semi-linear system of wave equations with Hamilton structure. We prove the existence of global smooth solution of the system for subcritical case by using conservation of energy and Strichartz's estimate. On the basis of Morawetz-Pohozev identity, we obtain the same result for the critical case.