6 results
Search Results
Now showing 1 - 6 of 6
- ItemOn M-stationary points for a stochastic equilibrium problem under equilibrium constraints in electricity spot market modeling(Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2007) Henrion, René; Römisch, WernerModeling several competitive leaders and followers acting in an electricity market leads to coupled systems of mathematical programs with equilibrium constraints, called equilibrium problems with equilibrium constraints (EPECs). We consider a simplified model for competition in electricity markets under uncertainty of demand in an electricity network as a (stochastic) multi-leader-follower game. First order necessary conditions are developed for the corresponding stochastic EPEC based on a result of Outrata [17]. For applying the general result an explicit representation of the co-derivative of the normal cone mapping to a polyhedron is derived (Proposition 3.2). Later the co-derivative formula is used for verifying constraint qualifications and for identifying M-stationary solutions of the stochastic EPEC if the demand is represented by a finite number of scenarios.
- ItemAre Quasi-Monte Carlo algorithms efficient for two-stage stochastic programs?(Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2016) Heitsch, Holger; Leövey, Hernan; Römisch, WernerQuasi-Monte Carlo algorithms are studied for designing discrete approximations of two-stage linear stochastic programs with random right-hand side and continuous probability distribution. The latter should allow for a transformation to a distribution with independent marginals. The twostage integrands are piecewise linear, but neither smooth nor lie in the function spaces considered for QMC error analysis. We show that under some weak geometric condition on the two-stage model all terms of their ANOVA decomposition, except the one of highest order, are continuously differentiable and that first and second order ANOVA terms have mixed first order partial derivatives and belong to L2. Hence, randomly shifted lattice rules (SLR) may achieve the optimal rate of convergence O(n-1+delta) with 2 (0; 1 2 ] and a constant not depending on the dimension if the effective superposition dimension is at most two. We discuss effective dimensions and dimension reduction for two-stage integrands. The geometric condition is shown to be satisfied almost everywhere if the underlying probability distribution is normal and principal component analysis (PCA) is used for transforming the covariance matrix. Numerical experiments for a large scale two-stage stochastic production planning model with normal demand show that indeed convergence rates close to the optimal are achieved when using SLR and randomly scrambled Sobol’ point sets accompanied with PCA for dimension reduction.
- ItemConditioning of linear-quadratic two-stage stochastic optimization problems(Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2013) Emich, Konstantin; Henrion, René; Römisch, WernerIn this paper a condition number for linear-quadratic two-stage stochastic optimization problems is introduced as the Lipschitz modulus of the multifunction assigning to a (discrete) probability distribution the solution set of the problem. Being the outer norm of the Mordukhovich coderivative of this multifunction, the condition number can be estimated from above explicitly in terms of the problem data by applying appropriate calculus rules. Here, a chain rule for the extended partial second-order subdifferential recently proved by Mordukhovich and Rockafellar plays a crucial role. The obtained results are illustrated for the example of two-stage stochastic optimization problems with simple recourse.
- ItemDiscrepancy distances and scenario reduction in two-stage stochastic integer programming(Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2007) Henrion, René; Küchler, Christian; Römisch, WernerPolyhedral discrepancies are relevant for the quantitative stability of mixed-integer two-stage and chance constrained stochastic programs. We study the problem of optimal scenario reduction for a discrete probability distribution with respect to certain polyhedral discrepancies and develop algorithms for determining the optimally reduced distribution approximately. Encouraging numerical experience for optimal scenario reduction is provided.
- ItemProblem-based optimal scenario generation and reduction in stochastic programming(Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2018) Henrion, Réne; Römisch, WernerScenarios are indispensable ingredients for the numerical solution of stochastic programs. Earlier approaches to optimal scenario generation and reduction are based on stability arguments involving distances of probability measures. In this paper we review those ideas and suggest to make use of stability estimates based only on problem specific data. For linear two-stage stochastic programs we show that the problem-based approach to optimal scenario generation can be reformulated as best approximation problem for the expected recourse function which in turn can be rewritten as a generalized semi-infinite program. We show that the latter is convex if either right-hand sides or costs are random and can be transformed into a semi-infinite program in a number of cases. We also consider problem-based optimal scenario reduction for two-stage models and optimal scenario generation for chance constrained programs. Finally, we discuss problem-based scenario generation for the classical newsvendor problem.
- ItemScenario reduction in stochastic programming with respect to discrepancy distances(Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2006) Henrion, René; Küchler, Christian; Römisch, WernerDiscrete approximations to chance constraints and mixed-integertwo-stage stochastic programs require moderately sized scenario sets. The relevant distances of (multivariate) probability distributions for deriving quantitative stability results for such stochastic programs are B-discrepancies, where the class B of Borel sets depends on their structural properties. Hence, the optimal scenario reduction problem for such models is stated with respect to B-discrepancies. In this paper, upper and lower bounds, and some explicit solutions for optimal scenario reduction problems are derived. In addition, we develop heuristic algorithms for determining nearly optimally reduced probability measures, discuss the case of the cell discrepancy (or Kolmogorov metric) in some detail and provide some numerical experience.