We consider the metric projection operator from the real Hilbert space onto a strongly convex set. We prove that the restriction of this operator on the complement of some neighborhood of the strongly convex set is Li...
详细信息
We consider the metric projection operator from the real Hilbert space onto a strongly convex set. We prove that the restriction of this operator on the complement of some neighborhood of the strongly convex set is Lipschitz continuous with the Lipschitz constant strictly less than 1. This property characterizes the class of strongly convex sets and (to a certain degree) the Hilbert space. We apply the results obtained to the question concerning the rate of convergence for the gradient projection algorithm with differentiable convex function and strongly convex set. (C) 2012 Elsevier Inc. All rights reserved.
A domain decomposition method (DDM) is presented to solve the distributed optimal control problem. The optimal control problem essentially couples an elliptic partial differential equation with respect to the state va...
详细信息
A domain decomposition method (DDM) is presented to solve the distributed optimal control problem. The optimal control problem essentially couples an elliptic partial differential equation with respect to the state variable and a variational inequality with respect to the constrained control variable. The proposed algorithm, called SA-GP algorithm, consists of two iterative stages. In the inner loops, the Schwarz alternating method (SA) is applied to solve the state and co-state variables, and in the outer loops the gradient projection algorithm (GP) is adopted to obtain the control variable. Convergence of iterations depends on both the outer and the inner loops, which are coupled and affected by each other. In the classical iteration algorithms, a given tolerance would be reached after sufficiently many iteration steps, but more iterations lead to huge computational cost. For solving constrained optimal control problems, most of the computational cost is used to solve PDEs. In this paper, a proposed iterative number independent of the tolerance is used in the inner loops so as to save a lot of computational cost. The convergence rate of L-2-error of control variable is derived. Also the analysis on how to choose the proposed iteration number in the inner loops is given. Some numerical experiments are performed to verify the theoretical results. (C) 2011 Elsevier B.V. All rights reserved.
An optimal control problem governed by the first bi-harmonic equation with the integral constraint for the state and its spectral approximations based on a mixed formulation are investigated. The optimality conditions...
详细信息
An optimal control problem governed by the first bi-harmonic equation with the integral constraint for the state and its spectral approximations based on a mixed formulation are investigated. The optimality conditions of the exact and the discrete optimal control systems are derived. The a priori error estimates of high order spectral accuracy are obtained. Furthermore, a simple and efficient iterative algorithm is proposed to solve mixed discrete system. Some numerical examples are performed to verify the theoretical results.
This paper presents a new dual formulation for quadratically constrained convex programs. The special structure of the derived dual problem allows us to apply the gradient projection algorithm to produce a simple expl...
详细信息
This paper presents a new dual formulation for quadratically constrained convex programs. The special structure of the derived dual problem allows us to apply the gradient projection algorithm to produce a simple explicit method involving only elementary vector-matrix operations, proven to converge at a linear rate.
An integral state-constrained optimal control problem governed by an elliptic partial differential equation and its finite element approximation are considered. The finite element approximation is constructed on multi...
详细信息
An integral state-constrained optimal control problem governed by an elliptic partial differential equation and its finite element approximation are considered. The finite element approximation is constructed on multimeshes. An L-2-norm a priori error estimate of the finite element approximation is obtained. Further, some superconvergence results are proved. Based on these superconvergence results, almost optimal L-infinity-norm error estimates are derived. Some recovery algorithms are then proposed to produce a posteriori error estimators of gradient type. To solve the finite element system, a simple and yet efficient iterative gradient projection algorithm is proposed and its convergence rate is proved. Some numerical examples are performed to confirm theoretical analysis.
In this paper implicit function-based parameterizations for orthogonal and oblique rotation matrices are proposed. The parameterizations are used to construct Newton algorithms for minimizing differentiable rotation c...
详细信息
In this paper implicit function-based parameterizations for orthogonal and oblique rotation matrices are proposed. The parameterizations are used to construct Newton algorithms for minimizing differentiable rotation criteria applied to m factors and p variables. The speed of the new algorithms is compared to that of existing algorithms and to that of Newton algorithms based on alternative parameterizations. Several rotation criteria were examined and the algorithms were evaluated over a range of values for m. Initial guesses for Newton algorithms were improved by subconvergence iterations of the gradient projection algorithm. Simulation results suggest that no one algorithm is fastest for minimizing all criteria for all values of m. Among competing algorithms, the gradient projection algorithm alone was faster than the implicit function algorithm for minimizing a quartic criterion over oblique rotation matrices when m is large. In all other conditions, however, the implicit function algorithms were competitive with or faster than the fastest existing algorithms. The new algorithms showed the greatest advantage over other algorithms when minimizing a nonquartic component loss criterion.
The supply chain network equilibrium problem with capacity constraints (SCNE-C) is an extension of the supply chain network equilibrium problem (SCNE), which also takes into account capacity constraints which refer to...
详细信息
The supply chain network equilibrium problem with capacity constraints (SCNE-C) is an extension of the supply chain network equilibrium problem (SCNE), which also takes into account capacity constraints which refer to the maximum production capacity for a manufacturer or the maximum storage/display space for a retailer. Due to inherent link interactions in the demand functions and cost functions, the SCNE-C problem is formulated as a mathematical model using the variational inequality ( VI) approach. This VI model is characterised by the so-called Wardrop second principle ( in terms of the 'generalised' route cost). To solve the model, a path-based four-loop nested diagonalisation method, along with a supernetwork representation, is proposed and demonstrated with a few numerical examples. The obtained results fully comply with the Wardrop second principle at both retailer sector and demand markets and can provide useful route information of the product. In addition, the stricter the capacity constraints imposed, the lower the quantity demanded will be, and provided at a higher product price. The concepts developed in this paper can be extended into many other spatial price equilibrium problems.
The traffic assignment problem with nonadditive path costs was discussed, especially considering the emission effect on traditional traffic assignment model. Based on Venigalla's research, the traffic assignment m...
详细信息
The traffic assignment problem with nonadditive path costs was discussed, especially considering the emission effect on traditional traffic assignment model. Based on Venigalla's research, the traffic assignment model considering emission at different engine operating modes was completed. A path-based algorithm namely gradientprojection (GP) algorithm is used to solve the above mathematic model. The case study shows that the GP algorithm has a rapid convergence speed than FW algorithm. The results also show that GP algorithm can be widely used in solving the traffic assignment problem with nonadditive path costs.
Noisy independent component analysis (ICA) is viewed as a method of factor rotation in exploratory factor analysis (EFA). Starting from an initial EFA solution, rather than rotating the loadings towards simplicity, th...
详细信息
ISBN:
(纸本)9783540744931
Noisy independent component analysis (ICA) is viewed as a method of factor rotation in exploratory factor analysis (EFA). Starting from an initial EFA solution, rather than rotating the loadings towards simplicity, the factors are rotated orthogonally towards independence. An application to Thurstone's box problem in psychometrics is presented using a new data matrix containing measurement error. Results show that the proposed rotational approach to noisy ICA recovers the components used to generate the mixtures quite accurately and also produces simple loadings.
In this article, we study algorithms for computing the nonparametric maximum likelihood estimator (NPMLE) of the failure function with two types of censored data: doubly censored data and (type 2) interval-censored da...
详细信息
In this article, we study algorithms for computing the nonparametric maximum likelihood estimator (NPMLE) of the failure function with two types of censored data: doubly censored data and (type 2) interval-censored data. We consider two projection methods, namely the iterative convex minorant algorithm (ICM) and a generalization of the Rosen algorithm (GR) and compare these methods to the well-known EM algorithm. The comparison conducted via simulation studies shows that the hybrid algorithms that alternately use the EM and GR for doubly censored data or, alternately, use the EM and ICM for (type 2) interval-censored data appear to be much more efficient than the EM, especially in large sample situation.
暂无评论