In this paper, on the basis of the logarithmic barrier function and KKT conditions, we propose a combined homotopy infeasible interior-point method (CHIIP) for convex nonlinear programming problems. For any convex n...In this paper, on the basis of the logarithmic barrier function and KKT conditions, we propose a combined homotopy infeasible interior-point method (CHIIP) for convex nonlinear programming problems. For any convex nonlinear programming, without strict convexity for the logarithmic barrier function, we get different solutions of the convex programming in different cases by CHIIP method.展开更多
Under the environment of electric power market, economic dispatch (ED) problem should consider network constraints, unit ramp rates, besides the basic constraints. For this problem, it is important to establish the ef...Under the environment of electric power market, economic dispatch (ED) problem should consider network constraints, unit ramp rates, besides the basic constraints. For this problem, it is important to establish the effective model and algorithm. This paper examines the decoupled conditions that affect the solution optimality to this problem. It proposes an effective model and solution method. Based on the look-ahead technique, it finds the number of time intervals to guarantee the solution optimality. Next, an efficient technique for finding the optimal solution via the interior point methods is described. Test cases, which include dispatching six units over 5 time intervals on the IEEE 30 test system with line flows and ramp constraints are presented. Results indicate that the computational effort as measured by iteration counts or execution time varies only modestly with the problem size.展开更多
This paper proposes an infeasible interior-point algorithm with full-Newton step for linear complementarity problem,which is an extension of Roos about linear optimization. The main iteration of the algorithm consists...This paper proposes an infeasible interior-point algorithm with full-Newton step for linear complementarity problem,which is an extension of Roos about linear optimization. The main iteration of the algorithm consists of a feasibility step and several centrality steps. At last,we prove that the algorithm has O(nlog n/ε) polynomial complexity,which coincides with the best known one for the infeasible interior-point algorithm at present.展开更多
In this paper,we are mainly devoted to solving fixed point problems in more general nonconvex sets via an interior point homotopy method.Under suitable conditions,a constructive proof is given to prove the existence o...In this paper,we are mainly devoted to solving fixed point problems in more general nonconvex sets via an interior point homotopy method.Under suitable conditions,a constructive proof is given to prove the existence of fixed points,which can lead to an implementable globally convergent algorithm.展开更多
A new higher-order accurate space-time discontinuous Galerkin(DG)method using the interior penalty flux and discontinuous basis functions,both in space and in time,is pre-sented and fully analyzed for the second-order...A new higher-order accurate space-time discontinuous Galerkin(DG)method using the interior penalty flux and discontinuous basis functions,both in space and in time,is pre-sented and fully analyzed for the second-order scalar wave equation.Special attention is given to the definition of the numerical fluxes since they are crucial for the stability and accuracy of the space-time DG method.The theoretical analysis shows that the DG discre-tization is stable and converges in a DG-norm on general unstructured and locally refined meshes,including local refinement in time.The space-time interior penalty DG discre-tization does not have a CFL-type restriction for stability.Optimal order of accuracy is obtained in the DG-norm if the mesh size h and the time stepΔt satisfy h≅CΔt,with C a positive constant.The optimal order of accuracy of the space-time DG discretization in the DG-norm is confirmed by calculations on several model problems.These calculations also show that for pth-order tensor product basis functions the convergence rate in the L∞and L2-norms is order p+1 for polynomial orders p=1 and p=3 and order p for polynomial order p=2.展开更多
The finite-dimensional variational inequality problem (VIP) has been studied extensively in the literature because of its successful applications in many fields such as economics, transportation, regional science and ...The finite-dimensional variational inequality problem (VIP) has been studied extensively in the literature because of its successful applications in many fields such as economics, transportation, regional science and operations research. Barker and Pang[1] have given an excellent survey of theories, methods and applications of VIPs.展开更多
This work focuses on the application of the reconstruction method of differentiated backprojection (DBP)-projection onto convex sets (POCS) in the interior problem.First,we present the definition of the interior p...This work focuses on the application of the reconstruction method of differentiated backprojection (DBP)-projection onto convex sets (POCS) in the interior problem.First,we present the definition of the interior problem and real truncated Hilbert transform,and then outline the implementation steps of DBP-POCS.After that,we introduce the middle-part known condition for region of interest (ROI) accurate reconstruction and the unique condition of the interior problem,and verify the uniqueness and stability of the interior problem accurate reconstruction through numerical experiments,and then compare the results for the interior problem in reconstruction images using filtered backprojection (FBP).In addition,the authors also design the application models of ROI reconstruction and make an initial attempt to the application of DBP-POCS method in the interior problem.展开更多
On the basis of the formulations of the logarithmic barrier function and the idea of following the path of minimizers for the logarithmic barrier family of problems the so called "centralpath" for linear pro...On the basis of the formulations of the logarithmic barrier function and the idea of following the path of minimizers for the logarithmic barrier family of problems the so called "centralpath" for linear programming, we propose a new framework of primal-dual infeasible interiorpoint method for linear programming problems. Without the strict convexity of the logarithmic barrier function, we get the following results: (a) if the homotopy parameterμcan not reach to zero,then the feasible set of these programming problems is empty; (b) if the strictly feasible set is nonempty and the solution set is bounded, then for any initial point x, we can obtain a solution of the problems by this method; (c) if the strictly feasible set is nonempty and the solution set is unbounded, then for any initial point x, we can obtain a (?)-solution; and(d) if the strictly feasible set is nonempty and the solution set is empty, then we can get the curve x(μ), which towards to the generalized solutions.展开更多
The incorrect surface current may be obtained in the vicinity of the resonant frequencies when the method of moments is used to solve either the electric or magnetic field surface integral equation. An effective met...The incorrect surface current may be obtained in the vicinity of the resonant frequencies when the method of moments is used to solve either the electric or magnetic field surface integral equation. An effective method is presented to determine the correct surface current, i.e., the correct surface current, i.e., the correct surface current is composed of the non-resonant mode current and the normalized resonant mode current multiplied by an unknown complex factor. The unknown complex factor can be obtained by employing the condition that the total field inside the conducting closed body must be zero at specified interior points. A numerical example is given for an infinitely long and perfectly conducting circular cylinder at the interior resonance, and the calculated surface currents are in good agreement with the analytical ones. The validity and accuracy of the presented method is thus verified.展开更多
Optimal adjustment algorithm for p coordinates is a generalization of the optimal pair adjustment algorithm for linear programming, which in turn is based on von Neumann’s algorithm. Its main advantages are simplicit...Optimal adjustment algorithm for p coordinates is a generalization of the optimal pair adjustment algorithm for linear programming, which in turn is based on von Neumann’s algorithm. Its main advantages are simplicity and quick progress in the early iterations. In this work, to accelerate the convergence of the interior point method, few iterations of this generalized algorithm are applied to the Mehrotra’s heuristic, which determines the starting point for the interior point method in the PCx software. Computational experiments in a set of linear programming problems have shown that this approach reduces the total number of iterations and the running time for many of them, including large-scale ones.展开更多
Transmission line manipulations in a power system are necessary for the execution of preventative or corrective main- tenance in a network, thus ensuring the stability of the system. In this study, primal-dual interio...Transmission line manipulations in a power system are necessary for the execution of preventative or corrective main- tenance in a network, thus ensuring the stability of the system. In this study, primal-dual interior-point methods are used to minimize costs and losses in the generation and transmission of the predispatch active power flow in a hydroelectric system with previously scheduled line manipulations for preventative maintenance, over a period of twenty-four hours. The matrix structure of this problem and the modification that it imposes on the system is also broached in this study. From the computational standpoint, the effort required to solve a problem with or without line manipulations is similar, and the reasons for this are also discussed in this study. Computational results sustain our findings.展开更多
Low-order wavefront error account for a large proportion of wave aberrations.A compensation method for low order aberration of projection lithography objective based on Interior Point Method is presented.Compensation ...Low-order wavefront error account for a large proportion of wave aberrations.A compensation method for low order aberration of projection lithography objective based on Interior Point Method is presented.Compensation model between wavefront error and degree of movable lens freedom is established.Converting over-determined system to underdetermined system,the compensation is solved by Interior Point Method(IPM).The presented method is compared with direct solve the over-determined system.Then,other algorithm GA,EA and PS is compared with IPM.Simulation and experimental results show that the presented compensation method can obtained compensation with less residuals compared with direct solve the over-determined system.Also,the presented compensation method can reduce computation time and obtain results with less residuals compare with AGA,EA and PS.Moreover,after compensation,RMS of wavefront error of the experimental lithography projection objective decrease from 56.05 nm to 17.88 nm.展开更多
In this paper, a high order compact difference scheme and a multigrid method are proposed for solving two-dimensional (2D) elliptic problems with variable coefficients and interior/boundary layers on nonuniform grids....In this paper, a high order compact difference scheme and a multigrid method are proposed for solving two-dimensional (2D) elliptic problems with variable coefficients and interior/boundary layers on nonuniform grids. Firstly, the original equation is transformed from the physical domain (with a nonuniform mesh) to the computational domain (with a uniform mesh) by using a coordinate transformation. Then, a fourth order compact difference scheme is proposed to solve the transformed elliptic equation on uniform girds. After that, a multigrid method is employed to solve the linear algebraic system arising from the difference equation. At last, the numerical experiments on some elliptic problems with interior/boundary layers are conducted to show high accuracy and high efficiency of the present method.展开更多
The generation expansion planning is one of complex mixed-integer optimization problems, which involves a large number of continuous or discrete decision variables and constraints. In this paper, an interior point wit...The generation expansion planning is one of complex mixed-integer optimization problems, which involves a large number of continuous or discrete decision variables and constraints. In this paper, an interior point with cutting plane (IP/CP) method is proposed to solve the mixed-integer optimization problem of the electrical power generation expansion planning. The IP/CP method could improve the overall efficiency of the solution and reduce the computational time. Proposed method is combined with the Bender's decomposition technique in order to decompose the generation expansion problem into a master investment problem and a slave operational problem. The numerical example is presented to compare with the effectiveness of the proposed algorithm.展开更多
The element-free Galerkin method is proposed to solve free vibration of rectangular plates with finite interior elastic point supports and elastically restrained edges.Based on the extended Hamilton's principle for t...The element-free Galerkin method is proposed to solve free vibration of rectangular plates with finite interior elastic point supports and elastically restrained edges.Based on the extended Hamilton's principle for the elastic dynamics system,the dimensionless equations of motion of rectangular plates with finite interior elastic point supports and the edge elastically restrained are established using the element-free Galerkin method.Through numerical calculation,curves of the natural frequency of thin plates with three edges simply supported and one edge elastically restrained,and three edges clamped and the other edge elastically restrained versus the spring constant,locations of elastic point support and the elastic stiffness of edge elastically restrained are obtained.Effects of elastic point supports and edge elastically restrained on the free vibration characteristics of the thin plates are analyzed.展开更多
It is well known that the incorrect results will be given using either the electric or magnetic field integral equation to calculate the radar cross section (RCS) of a closed body at the interior resonance. In this pa...It is well known that the incorrect results will be given using either the electric or magnetic field integral equation to calculate the radar cross section (RCS) of a closed body at the interior resonance. In this paper, an effective iterative technique is used to correct the calculated surface current density from the electric field integral equation. The radar cross section is computed for an infinite conducting circular cylinder at the interior resonance, and the obtained results are in good agreement with the analytical results. The backscattering cross section of an infinite triangular cylinder in the vicinity of a resonant frequency is also calculated. It is shown that the presence method is efficient and accurate.展开更多
This paper proposes an interior-point technique for detecting the nondominated points of multi-objective optimization problems using the direction-based cone method.Cone method decomposes the multi-objective optimizat...This paper proposes an interior-point technique for detecting the nondominated points of multi-objective optimization problems using the direction-based cone method.Cone method decomposes the multi-objective optimization problems into a set of single-objective optimization problems.For this set of problems,parametric perturbed KKT conditions are derived.Subsequently,an interior point technique is developed to solve the parametric perturbed KKT conditions.A differentiable merit function is also proposed whose stationary point satisfies the KKT conditions.Under some mild assumptions,the proposed algorithm is shown to be globally convergent.Numerical results of unconstrained and constrained multi-objective optimization test problems are presented.Also,three performance metrics(modified generational distance,hypervolume,inverted generational distance)are used on some test problems to investigate the efficiency of the proposed algorithm.We also compare the results of the proposed algorithm with the results of some other existing popular methods.展开更多
文摘In this paper, on the basis of the logarithmic barrier function and KKT conditions, we propose a combined homotopy infeasible interior-point method (CHIIP) for convex nonlinear programming problems. For any convex nonlinear programming, without strict convexity for the logarithmic barrier function, we get different solutions of the convex programming in different cases by CHIIP method.
文摘Under the environment of electric power market, economic dispatch (ED) problem should consider network constraints, unit ramp rates, besides the basic constraints. For this problem, it is important to establish the effective model and algorithm. This paper examines the decoupled conditions that affect the solution optimality to this problem. It proposes an effective model and solution method. Based on the look-ahead technique, it finds the number of time intervals to guarantee the solution optimality. Next, an efficient technique for finding the optimal solution via the interior point methods is described. Test cases, which include dispatching six units over 5 time intervals on the IEEE 30 test system with line flows and ramp constraints are presented. Results indicate that the computational effort as measured by iteration counts or execution time varies only modestly with the problem size.
基金Supported by the National Natural Science Fund Finances Projects(71071119)
文摘This paper proposes an infeasible interior-point algorithm with full-Newton step for linear complementarity problem,which is an extension of Roos about linear optimization. The main iteration of the algorithm consists of a feasibility step and several centrality steps. At last,we prove that the algorithm has O(nlog n/ε) polynomial complexity,which coincides with the best known one for the infeasible interior-point algorithm at present.
基金Supported by the NNSF of China(11026079)Supported by the Youth Backbone Teacher Foundation of Henan Province(173)
文摘In this paper,we are mainly devoted to solving fixed point problems in more general nonconvex sets via an interior point homotopy method.Under suitable conditions,a constructive proof is given to prove the existence of fixed points,which can lead to an implementable globally convergent algorithm.
文摘A new higher-order accurate space-time discontinuous Galerkin(DG)method using the interior penalty flux and discontinuous basis functions,both in space and in time,is pre-sented and fully analyzed for the second-order scalar wave equation.Special attention is given to the definition of the numerical fluxes since they are crucial for the stability and accuracy of the space-time DG method.The theoretical analysis shows that the DG discre-tization is stable and converges in a DG-norm on general unstructured and locally refined meshes,including local refinement in time.The space-time interior penalty DG discre-tization does not have a CFL-type restriction for stability.Optimal order of accuracy is obtained in the DG-norm if the mesh size h and the time stepΔt satisfy h≅CΔt,with C a positive constant.The optimal order of accuracy of the space-time DG discretization in the DG-norm is confirmed by calculations on several model problems.These calculations also show that for pth-order tensor product basis functions the convergence rate in the L∞and L2-norms is order p+1 for polynomial orders p=1 and p=3 and order p for polynomial order p=2.
基金The NNSF (10071031) of China and National 973 Project.
文摘The finite-dimensional variational inequality problem (VIP) has been studied extensively in the literature because of its successful applications in many fields such as economics, transportation, regional science and operations research. Barker and Pang[1] have given an excellent survey of theories, methods and applications of VIPs.
基金supported by the National Natural Science Foundation of China (Grant No.60872116)
文摘This work focuses on the application of the reconstruction method of differentiated backprojection (DBP)-projection onto convex sets (POCS) in the interior problem.First,we present the definition of the interior problem and real truncated Hilbert transform,and then outline the implementation steps of DBP-POCS.After that,we introduce the middle-part known condition for region of interest (ROI) accurate reconstruction and the unique condition of the interior problem,and verify the uniqueness and stability of the interior problem accurate reconstruction through numerical experiments,and then compare the results for the interior problem in reconstruction images using filtered backprojection (FBP).In addition,the authors also design the application models of ROI reconstruction and make an initial attempt to the application of DBP-POCS method in the interior problem.
文摘On the basis of the formulations of the logarithmic barrier function and the idea of following the path of minimizers for the logarithmic barrier family of problems the so called "centralpath" for linear programming, we propose a new framework of primal-dual infeasible interiorpoint method for linear programming problems. Without the strict convexity of the logarithmic barrier function, we get the following results: (a) if the homotopy parameterμcan not reach to zero,then the feasible set of these programming problems is empty; (b) if the strictly feasible set is nonempty and the solution set is bounded, then for any initial point x, we can obtain a solution of the problems by this method; (c) if the strictly feasible set is nonempty and the solution set is unbounded, then for any initial point x, we can obtain a (?)-solution; and(d) if the strictly feasible set is nonempty and the solution set is empty, then we can get the curve x(μ), which towards to the generalized solutions.
文摘The incorrect surface current may be obtained in the vicinity of the resonant frequencies when the method of moments is used to solve either the electric or magnetic field surface integral equation. An effective method is presented to determine the correct surface current, i.e., the correct surface current, i.e., the correct surface current is composed of the non-resonant mode current and the normalized resonant mode current multiplied by an unknown complex factor. The unknown complex factor can be obtained by employing the condition that the total field inside the conducting closed body must be zero at specified interior points. A numerical example is given for an infinitely long and perfectly conducting circular cylinder at the interior resonance, and the calculated surface currents are in good agreement with the analytical ones. The validity and accuracy of the presented method is thus verified.
文摘Optimal adjustment algorithm for p coordinates is a generalization of the optimal pair adjustment algorithm for linear programming, which in turn is based on von Neumann’s algorithm. Its main advantages are simplicity and quick progress in the early iterations. In this work, to accelerate the convergence of the interior point method, few iterations of this generalized algorithm are applied to the Mehrotra’s heuristic, which determines the starting point for the interior point method in the PCx software. Computational experiments in a set of linear programming problems have shown that this approach reduces the total number of iterations and the running time for many of them, including large-scale ones.
文摘Transmission line manipulations in a power system are necessary for the execution of preventative or corrective main- tenance in a network, thus ensuring the stability of the system. In this study, primal-dual interior-point methods are used to minimize costs and losses in the generation and transmission of the predispatch active power flow in a hydroelectric system with previously scheduled line manipulations for preventative maintenance, over a period of twenty-four hours. The matrix structure of this problem and the modification that it imposes on the system is also broached in this study. From the computational standpoint, the effort required to solve a problem with or without line manipulations is similar, and the reasons for this are also discussed in this study. Computational results sustain our findings.
文摘Low-order wavefront error account for a large proportion of wave aberrations.A compensation method for low order aberration of projection lithography objective based on Interior Point Method is presented.Compensation model between wavefront error and degree of movable lens freedom is established.Converting over-determined system to underdetermined system,the compensation is solved by Interior Point Method(IPM).The presented method is compared with direct solve the over-determined system.Then,other algorithm GA,EA and PS is compared with IPM.Simulation and experimental results show that the presented compensation method can obtained compensation with less residuals compared with direct solve the over-determined system.Also,the presented compensation method can reduce computation time and obtain results with less residuals compare with AGA,EA and PS.Moreover,after compensation,RMS of wavefront error of the experimental lithography projection objective decrease from 56.05 nm to 17.88 nm.
文摘In this paper, a high order compact difference scheme and a multigrid method are proposed for solving two-dimensional (2D) elliptic problems with variable coefficients and interior/boundary layers on nonuniform grids. Firstly, the original equation is transformed from the physical domain (with a nonuniform mesh) to the computational domain (with a uniform mesh) by using a coordinate transformation. Then, a fourth order compact difference scheme is proposed to solve the transformed elliptic equation on uniform girds. After that, a multigrid method is employed to solve the linear algebraic system arising from the difference equation. At last, the numerical experiments on some elliptic problems with interior/boundary layers are conducted to show high accuracy and high efficiency of the present method.
文摘The generation expansion planning is one of complex mixed-integer optimization problems, which involves a large number of continuous or discrete decision variables and constraints. In this paper, an interior point with cutting plane (IP/CP) method is proposed to solve the mixed-integer optimization problem of the electrical power generation expansion planning. The IP/CP method could improve the overall efficiency of the solution and reduce the computational time. Proposed method is combined with the Bender's decomposition technique in order to decompose the generation expansion problem into a master investment problem and a slave operational problem. The numerical example is presented to compare with the effectiveness of the proposed algorithm.
基金Project supported by the National Natural Science Foundation of China (Grant No.10872163)the Natural Science Foundation of Education Department of Shaanxi Province (Grant No.08JK394)
文摘The element-free Galerkin method is proposed to solve free vibration of rectangular plates with finite interior elastic point supports and elastically restrained edges.Based on the extended Hamilton's principle for the elastic dynamics system,the dimensionless equations of motion of rectangular plates with finite interior elastic point supports and the edge elastically restrained are established using the element-free Galerkin method.Through numerical calculation,curves of the natural frequency of thin plates with three edges simply supported and one edge elastically restrained,and three edges clamped and the other edge elastically restrained versus the spring constant,locations of elastic point support and the elastic stiffness of edge elastically restrained are obtained.Effects of elastic point supports and edge elastically restrained on the free vibration characteristics of the thin plates are analyzed.
基金This project was supported by the Foundation of MOE of China (No. 00179).
文摘It is well known that the incorrect results will be given using either the electric or magnetic field integral equation to calculate the radar cross section (RCS) of a closed body at the interior resonance. In this paper, an effective iterative technique is used to correct the calculated surface current density from the electric field integral equation. The radar cross section is computed for an infinite conducting circular cylinder at the interior resonance, and the obtained results are in good agreement with the analytical results. The backscattering cross section of an infinite triangular cylinder in the vicinity of a resonant frequency is also calculated. It is shown that the presence method is efficient and accurate.
基金financial support from Council of Scientific and Industrial Research,India through a research fellowship(File No.09/1217(0025)/2017-EMR-I)to carry out this research workDebdas Ghosh acknowledges the research grant(MTR/2021/000696)from SERB,India to carry out this research work.
文摘This paper proposes an interior-point technique for detecting the nondominated points of multi-objective optimization problems using the direction-based cone method.Cone method decomposes the multi-objective optimization problems into a set of single-objective optimization problems.For this set of problems,parametric perturbed KKT conditions are derived.Subsequently,an interior point technique is developed to solve the parametric perturbed KKT conditions.A differentiable merit function is also proposed whose stationary point satisfies the KKT conditions.Under some mild assumptions,the proposed algorithm is shown to be globally convergent.Numerical results of unconstrained and constrained multi-objective optimization test problems are presented.Also,three performance metrics(modified generational distance,hypervolume,inverted generational distance)are used on some test problems to investigate the efficiency of the proposed algorithm.We also compare the results of the proposed algorithm with the results of some other existing popular methods.