S. Karademir, N
Total Page:16
File Type:pdf, Size:1020Kb
ESSAYS ON INTEGER PROGRAMMING IN MILITARY AND POWER MANAGEMENT APPLICATIONS by Serdar Karademir BSc, Middle East Technical University, 2006 MSc, Middle East Technical University, 2008 Submitted to the Graduate Faculty of the Swanson School of Engineering in partial fulfillment of the requirements for the degree of Doctor of Philosophy University of Pittsburgh 2013 UNIVERSITY OF PITTSBURGH SWANSON SCHOOL OF ENGINEERING This dissertation was presented by Serdar Karademir It was defended on July 1st, 2013 and approved by Oleg A. Prokopyev, Ph.D., Associate Professor, Department of Industrial Engineering Nan Kong, Ph.D., Assistant Professor , Biomedical Engineering, Purdue University Jayant Rajgopal, Ph.D., Associate Professor, Department of Industrial Engineering Andrew J. Schaefer, Ph.D., Professor, Department of Industrial Engineering Dissertation Director: Oleg A. Prokopyev, Ph.D., Associate Professor, Department of Industrial Engineering ii Copyright c by Serdar Karademir 2013 iii ESSAYS ON INTEGER PROGRAMMING IN MILITARY AND POWER MANAGEMENT APPLICATIONS Serdar Karademir, PhD University of Pittsburgh, 2013 This dissertation presents three essays on important problems motivated by military and power management applications. The array antenna design problem deals with optimal arrangements of substructures called subarrays. The considered class of the stochastic as- signment problem addresses uncertainty of assignment weights over time. The well-studied deterministic counterpart of the problem has many applications including some classes of the weapon-target assignment. The speed scaling problem is of minimizing energy consumption of parallel processors in a data warehouse environment. We study each problem to discover its underlying structure and formulate tailored mathematical models. Exact, approximate, and heuristic solution approaches employing advanced optimization techniques are proposed. They are validated through simulations and their superiority is demonstrated through exten- sive computational experiments. Novelty of the developed methods and their methodological contribution to the field of Operations Research is discussed through out the dissertation. iv TABLE OF CONTENTS PREFACE ......................................... xi 1.0 INTRODUCTION .................................1 1.1 PHASED ARRAY ANTENNA DESIGN...................1 1.2 STOCHASTIC ASSIGNMENT.........................3 1.3 SPEED SCALING OF PARALLEL PROCESSORS.............4 1.4 OVERVIEW OF THE DISSERTATION....................6 2.0 IRREGULAR POLYOMINO TILING VIA INTEGER PROGRAM- MING WITH APPLICATION IN PHASED ARRAY ANTENNA DE- SIGN .........................................7 2.1 INTRODUCTION................................7 2.2 MODEL FORMULATION........................... 11 2.2.1 Nonlinear Exact Set Covering Model.................. 11 2.2.2 Linearized Model............................. 16 2.3 COLUMN GENERATION APPROACH................... 18 2.3.1 A New Branching Strategy........................ 18 2.3.2 Delayed Column Generation (DCG)................... 23 2.3.3 A New Lower-Bounding Scheme..................... 25 2.4 HEURISTIC APPROACHES.......................... 28 2.4.1 Construction Heuristics.......................... 28 2.4.1.1 Zoom-in Algorithm (ZiA):.................... 28 2.4.1.2 Meta-rectangle Tiling Algorithm (MrTA):........... 35 2.4.2 Improvement Heuristics......................... 37 v 2.5 APPROXIMATION BOUNDS......................... 48 2.6 COMPUTATIONAL RESULTS........................ 52 2.6.1 Results for Exact Approaches...................... 52 2.6.2 Results for Heuristic Approaches.................... 56 2.7 PHASED ARRAY ANTENNA SIMULATIONS............... 58 2.8 CONCLUSION................................. 60 2.9 ACKNOWLEDGMENT............................. 74 3.0 ON GREEDY APPROXIMATION ALGORITHMS FOR A CLASS OF TWO-STAGE STOCHASTIC ASSIGNMENT PROBLEMS .... 75 3.1 INTRODUCTION................................ 75 3.2 GREEDY APPROXIMATION ALGORITHMS................ 79 3.2.1 Basic Greedy Approach.......................... 79 3.2.2 Greedy Approach of Escoffier et al. [50]................. 80 3.3 NECESSARY OPTIMALITY CONDITION................. 81 3.4 ENHANCED GREEDY APPROACH..................... 83 3.4.1 Improving the first-stage assignment (EGA-I)............. 84 3.4.2 Improving the second-stage assignment (EGA-II)........... 89 3.4.3 Improving EGA with Local Search (EGA LS)............. 92 3.4.4 Analytical Observations......................... 96 3.5 COMPUTATIONAL EXPERIMENTS..................... 99 3.5.1 Setup.................................... 99 3.5.2 Results and Discussion.......................... 102 3.6 CONCLUSION................................. 109 3.7 ACKNOWLEDGMENTS............................ 109 4.0 ON SPEED SCALING VIA INTEGER PROGRAMMING ....... 110 4.1 INTRODUCTION................................ 110 4.2 OPTIMALITY CONDITIONS......................... 113 4.3 POLYNOMIALLY SOLVABLE CASES AND COMPLEXITY....... 115 4.4 A GREEDY APPROXIMATION ALGORITHM............... 117 4.5 FPTAS FOR SPECIAL CASE......................... 118 vi 4.5.1 A Dynamic Programming Approach................... 118 4.5.2 FPTAS................................... 119 4.6 AN OUTER APPROXIMATION ALGORITHM............... 121 4.7 COMPUTATIONAL EXPERIMENTS..................... 122 4.7.1 Implementation and Setup........................ 122 4.7.2 Results and Discussion.......................... 124 4.8 ACKNOWLEDGEMENT............................ 126 5.0 CONCLUSION ................................... 130 BIBLIOGRAPHY .................................... 132 vii LIST OF TABLES 1 Computational results for tetromino family: exact approaches......... 54 2 Computational results for pentomino family: exact approaches......... 55 3 Computational results for octomino family: exact approaches.......... 56 4 Computational results for octomino family: ZiA................. 57 5 Computational results for pentomino and octomino families: MrTA...... 58 6 Results for uncorrelated instances......................... 104 7 Results for correlated instances.......................... 105 8 Results for pairwise-correlated instances..................... 106 9 Results for split-like instances........................... 107 10 Results for interleaved-like instances........................ 108 11 Results for CPLEX vs. OA implementation.................... 125 12 Results for OA implementation.......................... 127 13 Results comparing OA and DP.......................... 128 viii LIST OF FIGURES 1 Monomino (F1), domino (F2), tromino (F3), tetromino (F4), and pentomino (F5) families.....................................8 2 Array and two subarray examples: rectangular and L-shaped..........9 3 Radiation patterns for different time delay control scenarios........... 10 1 4 Two polyominoes from pentomino family, F5. Polyomino fpq is shown with its rectangle hull.................................... 12 5 Minimizing versus maximizing entropy on a 20 × 20 board using pentomino family........................................ 15 6 The bijections we use and example fractional domino tilings.......... 21 7 F-, T-, and W-pentominoes at 5 × 5 zoom level................. 30 8 Obtaining a 250 × 250 tiling of pentominoes using after two successive appli- cations of zoom-in heuristic. Solution time is less than one second....... 31 9 Zooming-in the point-up octomino. The initial 12 × 16 tiling zoomed-in twice at 12 × 12 zoom level to obtain a 1728 × 2304 tiling with 0:69% optimality gap. 33 10 Creating meta-rectangles.............................. 36 11 An example of a meta-rectangle generated by Algorithm3 using initial set of rectangles f(2; 8); (4; 16); (16; 17)g and its final octomino tiling......... 38 12 Dimensions of meta-rectangles generated by MrTA and the initialization step of the tiling algorithm............................... 39 13 Initialization in Procedure Retile:A! do not cover, B! retile, and C! penalize. 40 14 Randomizing an octomino meta-rectangle tiling................. 44 15 The subset of octomino family, F8, used in experiments............. 53 ix 16 Summary of array antenna simulations results: time delay control at element level (E), 2 × 4 rectangular subarray level (R), partially optimized octomino subarray level (P), and completely optimized octomino subarray level (P+). 61 17 Partially and completely optimized polyomino tilings for 64 × 64 array size.. 62 18 Radiation patterns for 64 × 64 array size with time delay control at 2 × 4 rect- angular subarray level (R), partially optimized octomino subarray level (P), and completely optimized octomino subarray level (P+)............ 64 19 Partially and completely optimized polyomino tilings for 128 × 128 array size. 66 20 Radiation patterns for 128 × 128 array size with time delay control at element level (E), 2 × 4 rectangular subarray level (R), partially optimized octomino subarray level (P), and completely optimized octomino subarray level (P+). 68 21 Partially and completely optimized polyomino tilings for 256 × 256 array size. 70 22 Radiation patterns for 256×256 array size with time delay control at 2×4 rect- angular subarray level (R), partially optimized octomino subarray level (P), and completely optimized octomino subarray level (P+)............ 72 23 A counterexample to show that the necessary optimality condition given by Proposition 11 is not sufficient for optimality. Only arcs with nonzero weight are shown...................................... 82 24 The neighborhood N1................................ 94 25 The neighborhood N2...............................