Production Capacity Modeling of Alternative, Nonidentical, Flexible Machines
|
|
- Blanche Tucker
- 5 years ago
- Views:
Transcription
1 The International Journal of Flexible Manufacturing Systems, 14, , 2002 c 2002 Kluwer Academic Publishers Manufactured in The Netherlands Production Capacity Modeling of Alternative, Nonidentical, Flexible Machines GEORGE LIBEROPOULOS Department of Mechanical and Industrial Engineering, University of Thessaly, Pedion Areos, GR Volos, Greece glib@mieuthgr Abstract Analyzing the production capacity of a flexible manufacturing system consisting of a number of alternative, nonidentical, flexible machines, where each machine is capable of producing several different part types simultaneously (by flexibly allocating its production capacity among these part types), is not a trivial task The production capacity set of such a system is naturally expressed in terms of the machine-specific production rates of all part types In this paper we also express it in terms of the total production rates of all part types over all machines More specifically, we express the capacity set as the convex hull of a set of points corresponding to all possible assignments of machines to part types, where in each assignment each machine allocates all its capacity to only one part type First, we show that within each subset of assignments having a given number of machines assigned to each part type, there is a unique assignment that corresponds to an extreme point of the capacity set Then, we propose a procedure for generating all the extreme points and facets of the capacity set Numerical experience shows that when the number of part types is less than four, the size of the capacity set (measured in terms of the number of variables times the number of constraints) is smaller, if the capacity set is expressed in terms of the total production rates of all part types over all machines than if it is expressed in terms of the machine-specific production rates of all part types When the number of part types is four or more, however, the opposite is true Key Words: alternative, nonidentical machines, production capacity, production flow control 1 Introduction and motivation In the early 1980s, Kimemia and Gershwin (1983) formulated a paradigm for the optimal control of production flow in a flexible manufacturing system Their paradigm received considerable attention in the years that followed and spawned a large number of follow-on papers Some of the material in these papers is included in the books by Gershwin (1994), Sethi and Zhang (1994), and Maimon, Khmelnitsky, and Kogan (1998) Kimemia and Gershwin s paradigm concerns a flexible manufacturing system consisting of M machine groups Machine group m has a number of identical machines, where each machine is perfectly flexible and can allocate all or part of its production capacity to N different part types simultaneously The machines are unreliable, and the main difficulty the control system faces is to meet production requirements while the machines fail and are repaired at random times The system is modeled as a continuous-time, mixed-state (having both discrete and continuous elements) stochastic dynamic system The discrete constituent of the state, α, is the vector of machine states, where each machine may be either up or down The continuous portion of the state, x, is the vector of surpluses, ie, the cumulative differences between production and demand for all part types The objective is to minimize
2 346 LIBEROPOULOS these differences The production rate vector, u (the control), is constrained to be within a capacity set, (α), which is determined by the set of operational machines The production surplus vector satisfies the differential equation ẋ(t) = u(t) d, where d is a vector of constant demand rates A feedback control law that determines the current production rate as a function of the current machine state and current production surplus, u(x, α), is sought Kimemia and Gershwin (1983) formulated this problem as a stochastic optimal control problem and set out to solve it using dynamic programming They assumed that the machine state α(t) is a continuous-time Markov process and showed that the optimal control of the system is described by a set of stochastic partial differential (Bellman) equations involving the optimal cost-to-go or value function J (x, α) More specifically, they showed that the optimal feedback law u (x, α) is the solution to the linear program (LP), min x J (x, α) u subject to u (α) In other words, the optimal production rate vector u (x, α) is the control u in the production capacity set (α) that minimizes the inner product x J (x, α) u and thereby causes the steepest descent in J (x, α) Unfortunately, an exact solution of the Bellman equations for the J (x, α) is practically infeasible for systems producing more than one part type Caramanis and Liberopoulos (1992) proposed quadratic approximations of the value functions J (x, α) and used simulation and infinitesimal perturbation analysis to optimize the parameters of these functions (see also Liberopoulos and Caramanis, 1994, 1995) The simulation involves solving a linear program (LP) repeatedly in order to determine the optimal production rate vector u(x(t), α(t)), which changes dynamically in time, since x(t) and α(t) change dynamically in time This LP has the form, min c(x(t), α(t)) u, subject to u (α(t)), (1) where c(x(t), α(t)) is the gradient of the quadratic approximation of J (x(t), α(t)) In order to solve LP (1), (α) must be precisely defined To precisely define (α), consider the machine state α = (α 1,α 2,,α M ), where α m is the number of operational machines at machine group m Let ynm k be the rate at which machine group m performs operation k on type n parts (measured in parts per unit time) Let τnm k be the time required to complete the operation It then follows that ynm k τ nm k α m, for all m (2) n k The product ynm k τ nm k is the proportion of time in which all of the machines of machine group m perform operation k on part type n The left-hand side of inequality (2) is therefore the total amount of work performed on all part types on machine group m per unit time The inequality follows because the total amount of work per unit time performed on machine group m cannot exceed the time that the machines of machine group m are available By conservation of the flow, the rate at which any given operation at any given machine group is performed on part type n is the same as the rate at which any other operation at any other
3 PRODUCTION CAPACITY MODELING 347 machine group is performed on the same part type, and is equal to the total rate at which operation k is performed on part type n, u n, that is, u n = m ynm k, for all k and n (3) The production capacity set (α) isdefined to be the set of all production rate vectors u = (u 1, u 2,,u N ) such that there exist feasible production rates ynm k satisfying (2) and (3) Note that the constraint set of LP (1) is expressed in terms of the machine-specific production rates ynm k, whereas the objective function of LP (1) is expressed in terms of the total production rates u n (α) is the projection of a polyhedral set into a lower dimensional subspace It is therefore convex and polyhedral This set is a representation of the production capacity of the system A question that arises naturally is, Would it be computationally more efficient for solving LP (1) to express (α) in terms of the variables u n instead of the variables ynm k? This question is particularly important in view of the fact that LP (1) needs to be solved repeatedly during each simulation run (since x(t) and α(t) change dynamically in time) and this provided the motivation for writing this article More specifically, our goal in this article is to determine the size of (α), when (α) is expressed in terms of the variables u n, and compare it with the size of (α), when (α) is expressed in terms of the variables ynm k Once the size of (α) is known, it can be used as a surrogate for the computational efficiency of solving LP (1), in the sense that the smaller the size of (α), the bigger the computational efficiency of solving LP (1) The size of (α) is measured in terms of the number of variables multiplied by the number of constraints The problem of modeling the capacity of alternative, nonidentical resources arises in almost all manufacturing environments where process technology is evolving, as well as in many other application areas Surprisingly, there is not much published research on this subject Leachman and Carmon (1992) proposed a procedure to generate the capacity set of alternative machine types assuming that processing times among alternative machine types are identical or proportional across the operations that they can perform Hung and Wang (1997) extended this procedure and applied it to the bin allocation planning problem (a material planning problem) in the semiconductor industry They called their procedure acceptable material set generation procedure Both procedures were devised within typical, finite-horizon, production or material planning problem formulations, but assumed that processing times among alternative machines are identical (or proportional) In this article we assume that processing times among alternative machines are different The remaining part of this article is organized as follows In Section 2 we formulate the problem of determining the size of (α) for a particular machine state α, where a given number of machines is up We therefore omit the dependence on α and call this set, for notational simplicity In Section 3 we propose an algorithm for determining the number of extreme points and the number of nonredundant facets of, as expressed in terms of the total production rates u n In Section 4 we report on the numerical experience for systems with up to four part types and four machines, and in Section 5 we draw conclusions
4 348 LIBEROPOULOS 2 Problem formulation In the rest of this article, to simplify matters, we concentrate on a flexible manufacturing system consisting of M machine groups, where each machine group has only one machine This is without loss of generality because one can always view multiple identical machines as a single machine with a given capacity The machines are different from each other, but every machine can allocate all or part of its production capacity to N different part types We assume that each part type needs only one operation that can be performed on any of the machines In short, we consider a flexible manufacturing system consisting of M alternative, nonidentical, flexible machines, where each machine can allocate all or part of its production capacity to N different part types simultaneously Such a system was considered in Maimon and Gershwin (1988) Note that when parts visit several machines in series, the sequencing and mixing of part types should be taken into account when defining the production capacity set, as was pointed out by Lasserre (1992) For ease of exposition, we concentrate on a particular machine state, namely the state where all machines are up, that is, we assume that α m = 1, for m = 1, 2,,M This is without loss of generality because any other state in which some of the machines are up and the rest are down is equivalent to the state where all machines are up in a system with fewer machines For the state where all machines are up, the production capacity set, henceforth referred to as for notational simplicity, is defined to be the set of all production rate vectors u = (u 1, u 2,,u N ) such that there exist feasible machine-specific production rates y nm satisfying, N y nm τ nm 1, for m = 1, 2,,M, (4) n=1 u n = M y nm, for n = 1, 2,,N (5) m=1 Our goal is to determine the size of, as expressed in terms of the total production rates u n, and to compare it to the size of, as expressed in terms of machine-specific production rates y nm The size of, as expressed in terms of the production rates y nm,is (NM)M = NM 2, where NM is the number of the y nm variables and M is the number of constraints (4) The size of, as expressed in terms of the total production rates u n,isnf, where N is the number of the u n variables and F is the number of nonredundant facets of, as expressed in terms of u n We need to find F To illustrate this, consider a simple example of two machines and two part types, that is, M = 2 and N = 2, with parameters, τ 11 = 1/2,τ 21 = 1, for the first machine, and τ 12 = 1,τ 22 = 2/3, for the second machine Then, can be expressed in terms of machinespecific production rates y nm by the constraints 1/2y 11 + y 21 1 and y /3y 22 1 In this case, the size of is = 8 can also be expressed in terms of the total production rates u n by the constraints 1/5u 1 + 2/5u 2 1 and 1/3u 1 + 2/9u 2 1 In this case, the number of nonredundant facets of is 2, and therefore the size of is 2 2 = 4 For a graphical representation of these constraints see Figure 1
5 PRODUCTION CAPACITY MODELING 349 Figure 1, as expressed in terms of the total production rates u n, for a case where M = 2, N = 2 3 Extreme points and facets of Ω The production capacity set, as expressed in terms of the total production rates u n,isa convex polyhedron One way to describe a convex polyhedron is to characterize it as the convex hull of a finite set of points To find a set of points, expressed in terms of the rates u n, whose convex hull is, we consider all possible assignments of M machines to N part types, such that in each assignment, each machine allocates all its capacity to only one part type The number of such assignments is N M Let S ={1, 2,,N M } be a set of indexes to all such assignments, and let s be an element of S, that is, s S For every s S, let znm s be a variable that is equal to one, if machine m is assigned to part type n in assignment s, and zero, otherwise Because for any s S, each machine is assigned to only one part type, it follows that N n=1 zs nm = 1, for m = 1, 2,,M and all s S Also, for every s S, let pn s be the maximum total production rate of part type n in assignment s, that is, pn s = M m=1 zs nm /τ nm, for n = 1, 2,,N, and let p s = (p1 s, ps 2,,ps N ) Then, the production capacity set, as expressed in terms of the rates u n, is the convex hull of points 0 and p s, for all s S Note, however, that not all points p s are extreme points (or vertices) of To obtain a nonredundant polyhedral description of, we need to identify only those points p s that are extreme points of Before setting out to identify the extreme points of, as expressed in terms of the total production rates u n, let us first try to count them To count the number of extreme points of, as expressed in terms of the total production rates u n, consider all partitions of S into subsets S r, such that each subset S r contains all assignments that have a given number of machines assigned to each part type To define each subset S r S, let in r be the number of machines that are assigned to part type n in every s S r, and let i r = (i1 r, ir 2,,ir N ) Clearly, the ir n satisfy N in r = M and 0 ir n M, for n = 1, 2,,N n=1 With the above definition in mind, S r can be defined in terms of the in r as follows: { S r = s : M znm s = ir n }, n = 1, 2,,N m=1
6 350 LIBEROPOULOS The number of elements of each subset S r is equal to the number of ways of partitioning a population of M machines into N subpopulations of sizes in r, n = 1, 2,,N, and is given by M!/ N n=1 ir n! To illustrate this, let us revisit the example at the end of Section 2 The capacity constraints for this example are shown in Figure 1, where m denotes the production capacity set of machine m and is expressed in terms of the machine-specific production rates y nm In the example, S is partitioned into three subsets, S 1 ={1}, S 2 ={3}, and S 3 ={2, 4} Subset S 1 contains the only assignment where both (all) machines are assigned to part type 1, that is, i 1 = (2, 0) Subset S 2 contains the only assignment where both (all) machines are assigned to part type 2, that is, i 2 = (0, 2) Subset S 3 contains the only two assignments where one machine is assigned to part type 1 and the other machine is assigned to part type 2, that is, i 3 = (1, 1) We now have the following result Proposition 1 In every subset S r S there is a unique assignment, s r, such that ps r is an extreme point of Proof: Consider a particular subset of assignments, say subset S r1, and consider a particular assignment in this subset, say s 1, that is, s 1 S r1 Suppose that a different assignment, say s 2, can be obtained from assignment s 1 by changing the individual assignment of only one machine, say machine m 1, from part type n 1 to some other part type, say part type n 2 Clearly, s 2 belongs to a different subset, say S r2, that is, s 2 S r2 S r1 Subsets S r1 and S r2 are related by i r 2 n 1 = i r 1 n 1 1, i r 2 n 2 = i r 1 n 2 + 1, and i r 2 n = ir 1 n, for all n n 1, n 2 Now suppose that a new assignment, say s 1, such that s 1 s 1, can be obtained from assignment s 2 by changing the individual assignment of only one machine, say machine m 2, where m 2 m 1, from part type n 2 back to part type n 1 Clearly, s 1 S r 1 Point p s 1 may be either an interior point (see Figure 2(a)), or an exterior point (see Figure 2(b)) with respect to the supporting hyperplane defined by the edge connecting points p s 1 and p s 2 Inthefirst case, p s 1 cannot be an exterior point of In the second case, p s 1 cannot be an extreme point of Therefore, only one of the points s 1 and s 1 may be an extreme point of This argument Figure 2 Projection of production rate vectors p s 1, p s 2, and p s 1 on the u n1 u n2 plane
7 PRODUCTION CAPACITY MODELING 351 can be repeated to show that in subset S r1 S there is a unique assignment, s1 S r 1, such that p s 1 is an extreme point of Proposition 1 states that in every subset S r, there is a unique assignment, s r, which corresponds to an extreme point of All other assignments in S r correspond to interior points of To illustrate Proposition 1 consider again the example of Figure 1 In this example, subset S 3 contains the only two assignments (s = 2 and s = 4) where one machine is assigned to part type 1 and the other machine is assigned to part type 2, that is, i 3 = (1, 1) Of these two assignments only one corresponds to an extreme point of, namely assignment s = 2 The following result is a consequence of Proposition 1 Corollary 1 The number of extreme points of (excluding the origin 0) is ( M+N 1 N 1 ) Proof: By Proposition 1, the number of extreme points of is equal to the number of subsets S r that partition S This number is equal to the number of ways in which M identical items (corresponding to the M machines) can be placed in N in-line compartments (corresponding to the N part types) that are separated with N 1 partitions, where each compartment may contain all, some, or none of the items Viewed differently, this is equal to the number of combinations of N 1 partitions and M items in a row Corollary gives the number of extreme points of, as expressed in terms of the rates u n For example, if N = 2, the number of extreme points of is equal to M + 1 (see Figures 1 and 11) If N = 3, the number of extreme points of is equal to (M + 2)(M + 1)/2 (see Figures 4 7) What remains to be done is to generate these points and the facets of that they define An efficient way to screen out many of the points p s that are not extreme points of is given by Proposition 2 that follows To state Proposition 2 we first need the following definitions Consider a particular assignment s Of all machines that are assigned to part type n 1 in assignment s, let m s n 1 n 2 be that machine which, when switched from part type n 1 to part type n 2, yields the largest relative gain in production capacity, that is, m s n 1 n 2 [ ] = arg max τn1 m/τ n2 m (6) m:zn s 1 m =1 Let σ (s, n 1, n 2 ) be the assignment that is obtained from assignment s by changing the individual assignment of machine m s n 1 n 2 from part type n 1 to part type n 2, where m s n 1 n 2 is given by (6) Let P s be the set of all part types that have at least one machine assigned to them in assignment s, that is, P s ={n 1 : zn s 1 m = 1, for some m} We now have the following result Proposition 2 If m σ (s,n 1,n 2 ) n 2 n 1 m s n 1 n 2 or equivalently σ (σ (s, n 1, n 2 ), n 2, n 1 ) s, for any n 1 P s and any n 2 {1, 2,,N}, then p s is not an extreme point of Proof: The proof is similar to that of Proposition 1 Assignment σ (s, n 1, n 2 ) is obtained from assignment s if the individual assignment of machine m s n 1 n 2 is changed from part
8 352 LIBEROPOULOS Figure 3 Projection of production rate vectors p s, p σ (s,n 1,n 2 ), and p σ (σ (s,n 1,n 2 ),n 2,n 1 ) on the u n1 u n2 plane type n 1 to part type n 2 so that the largest relative gain in production capacity is obtained Assignment σ (σ (s, n 1, n 2 ), n 2, n 1 ) is obtained from assignment σ (s, n 1, n 2 ) if the individual assignment of machine m σ (s,n 1,n 2 ) n 2 n 1 is changed from part type n 2 back to part type n 1 so that the largest relative gain in production capacity is obtained If m σ (s,n 1,n 2 ) n 2 n 1 m s n 1 n 2, then σ (σ (s, n 1, n 2 ), n 2, n 1 ) s and clearly p σ (σ (s,n 1,n 2 ),n 2,n 1 ) is an exterior point with respect to the supporting hyperplane defined by the edge connecting points p s and p σ (s,n 1,n 2 ) (see Figure 3) This means that point p s is not an extreme point of Proposition 2 provides a sufficient but not necessary condition for screening out points p s that are not extreme points of In other words, if a point p s satisfies this condition, it can be eliminated from the list of candidate extreme points of If it does not satisfy this condition, however, it is not necessarily an extreme point of Numerical experience, however, has shown that in a typical case there are only a few points p s that do not satisfy this condition and are not extreme points of Figure 4 shows an example of the production capacity set that results after eliminating all points that satisfy the condition stated in Proposition 2, for a case where M = 4 and N = 3 Note that two of the points that are not eliminated, namely A and B, are interior points of Figure 4, as expressed in terms of the total production rates u n, for a case where M = 4, N = 3
9 PRODUCTION CAPACITY MODELING 353 Adefinite but computationally more demanding way of testing if a point p s is an extreme point of is to solve the LP: max w s subject to M M y nm = w s znm s /τ nm( = w s pn) s, for n = 1, 2,,N, m=1 m=1 N y nm τ nm 1, n=1 for m = 1, 2,,M, (7) where the znm s, n = 1, 2,,N, m = 1, 2,,M, are given 0 and 1 coefficients for each p s Let w s, ynm, n = 1, 2,,N, m = 1, 2,,M, be the optimal solution of LP (7) The optimal solution w s is the maximum scalar by which vector p s may be multiplied without violating the capacity constraints (4) and consequently the production capacity set Clearly, if w s = 1, then p s is an extreme point of Ifw s > 1, then p s is an interior point of A graphical representation of LP (7) is shown in Figure 1 for a case where M = 2 and N = 2, and consequently S ={1, 2, 3, 4} Note that point p 4 is not an extreme point of because there exists a scalar w 4, such that w 4 > 1 and w 4 p 4 Points p 1, p 2, and p 3,on the other hand, are all extreme points of Before proceeding any further with our analysis, the following parenthetical comment is in order Certainly, the solution w s = 1 and y nm = znm s /τ nm, for all n and m, is a feasible solution of LP (7) Call this solution χ s for short One way of testing if point p s is an extreme point of, which would be more efficient than solving LP (7) from scratch, would be to examine the sign of the reduced cost coefficients of LP (7) corresponding to solution χ s More specifically, if none of these reduced cost coefficients is negative, then solution χ s is the optimal solution, and therefore p s is an extreme point of Otherwise, χ s is not the optimal solution, and therefore p s is not an extreme point of Unfortunately, however, these reduced cost coefficients cannot be computed because not all of the M + N basic variables corresponding to solution χ s are known Instead, only M +1 of the basic variables are known, namely w s and those y nm for which the corresponding znm s variables are equal to one The remaining N 1 basic variables, whose value is zero, cannot be determined Therefore, in order to test whether a point p s is definitely an extreme point of or not, LP (7) must be solved from scratch Based on our previous analysis, we propose the following procedure for generating the production capacity set, as expressed in terms of the total production rates u n 31 Procedure for generating 1 Construct a directed graph G = (S, E), where the set of nodes S is the set of indexes to all assignments and the set of arcs E is the set of ordered pairs (s,σ(s, n 1, n 2 )), for all n 1 P s, n 2 {1, 2,,N}, and s S 2 Drop the nodes (and the arcs that are attached to them) for which the sufficient condition stated in Proposition 2 holds The remaining graph is no longer a directed graph because
10 354 LIBEROPOULOS both (s,σ(s, n 1, n 2 )) and (σ (s, n 1, n 2 ), s) are arcs of the graph, for all n 1 P s, n 2 {1, 2,,N}, and s S 3 Drop the nodes (and the arcs that are attached to them) for which the solution to LP (7), w s, is such that w s > 1 4 The dropped nodes represent all the interior points of Call the resulting graph G G represents the extreme points of and the edges between them 5 Identify the facets of from graph G For systems that produce N part types, the facets of are (N 1)-dimensional and are defined as the convex hull of at least N points or, alternatively, of at least N(N 2)-dimensional faces 4 Numerical results We used the procedure described in the previous section to construct the production capacity set for cases of up to N = 4 part types and M = 4 machines Figures 4 9 show the resulting graph G for these cases Tables 1 3 show the number of facets of, as expressed in terms of the total production rates u n, that are in the positive orthant, for the cases where N = 2, 3, and 4, shown in Figures 5 10 They also show the size of, as expressed in terms of the rates y nm and u n, respectively, for the same cases, where the following notation is used: C number of extreme points of (excluding point 0), as expressed in terms of the total production rates u n ; F( j) number of facets of, as expressed in terms of the rates u n that are defined by j (N 2)-dimensional hyperplanes and are in the positive orthant; F total number of facets of, as expressed in terms of the rates u n, that are in the positive orthant; Figure 5 Graph G for a case where N = 3, M = 2
11 PRODUCTION CAPACITY MODELING 355 Figure 6 Graph G for a case where N = 3, M = 3 Figure 7 Graph G for a case where N = 3, M = 4 S 1 size of (variables constraints), as expressed in terms of the machine-specific production rates y nm ; S 2 size of (variables constraints), as expressed in terms of the total production rates u n, where S 1 = (NM)M = NM 2 and S 2 = NF Based on the results reported in Tables 1 3, we make the following observations When N = 2, the structure of, as expressed in terms of the total production rates u n,is simple Namely, has M one-dimensional facets, each of which is defined as the convex hull of two zero-dimensional points In fact, can be constructed trivially as follows: Start from the point where all machines are assigned to one part type, say part type 2, that is, (u 1, u 2 ) = (0, M m=1 1/τ 2m) Then, lay the machine-specific production capacity sets in
12 356 LIBEROPOULOS Figure 8 Graph G for a case where N = 4, M = 2 Figure 9 Graph G for a case where N = 4, M = 3 order of decreasing-slope (see Figure 11) Clearly, in this case, S 1 = 2M 2, S 2 = 2M, and therefore S 1 S 2 When N > 2, however, the number of facets of, as expressed in terms of the total production rates u n, is greater than M, for M > 1 More specifically, has M facets that are defined as the convex hull of exactly N (N 2)-dimensional faces, but it also has other facets that are defined as the convex hull of more than N (N 2)-dimensional faces For example, when N = 3 and M = 2, consists of two triangles and one parallelogram (in R 3, a triangle is a two-dimensional facet defined as the convex hull of three one-dimensional faces (edges), and a parallelogram is a two-dimensional facet defined as the convex hull of four one-dimensional faces (edges)) (see Figure 5)
13 PRODUCTION CAPACITY MODELING 357 Table 1 Number of facets of, as expressed in terms of the rates u n, and size of, as expressed in terms of the rates y nm and u n, respectively, for N = 2 M C F(2) F S 1 S M M + 1 M M 2M 2 2M Table 2 Number of facets of, as expressed in terms of the rates u n, and size of, as expressed in terms of the rates y nm and u n, respectively, for N = 3 M C F(3) F(4) F S 1 S M (M+1)(M+2) 2 M (M 1)M 2 M(M+1) 2 3M 2 3M(M+1) 2 Table 3 Number of facets of, as expressed in terms of the rates u n, and size of, as expressed in terms of the rates y nm and u n, respectively, for N = 4 M C F(4) F(5) F(6) F S 1 S M (M+1)(M+2)(M+3) 2 M (M 1)M (M 2)(M 1)M 6 M(M+1)(M+2) 6 4M 2 4M(M+1)(M+2) 6 The results in Tables 1 3 reveal a pattern, which, if extrapolated to the general case of N part types and M machines, leads to the following conjecture Conjecture 1 The total number of facets of in the positive orthant, F, as expressed in terms of total production rates u n, is ( M+N 2 N 1 )
14 358 LIBEROPOULOS Figure 10 Graph G for a case where N = 4, M = 4 Figure 11, as expressed in terms of the total production rates u n, for a case where N = 2 By Conjecture 1, the size of, as expressed in terms of total production rates u n,is ( ) M + N 2 N(M + N 2)! S 2 = N = N 1 (N 1)!(M 1)! Given that that the size of expressed in terms of the rates y nm is S 1 = NM 2,we obtain the following result Corollary 2 S 1 S 2, for N < 4, and S 1 S 2, for N 4 Corollary 2 states that when the number of part types is less than four, the size of is smaller, if is expressed in terms of the total production rates of all part types over all machines, u n, than if it is expressed in terms of the machine-specific production rates of all part types, y nm When the number of part types is four or more, however, the opposite is true To elaborate this, suppose that N = 3, M = 10 Then, by Conjecture 1, F = 55 and S 2 = 3 55 = 165, whereas S 1 = (3 10) 10 = 300, that is, S 1 is approximately
15 PRODUCTION CAPACITY MODELING 359 twice as big as S 2 In this case it is preferable to express in terms of the total production rates u n Now suppose that N = 10, M = 10 Then, by Conjecture 1, F = 48,620 and S 2 = 10 48,620 = 486,200, whereas S 1 = (10 10) 10 = 1000, that is, S 2 is approximately 500 times larger than S 1! In this case it is preferable to express in terms the machine-specific production rates y nm 5 Conclusions The production capacity set,, ofm alternative, nonidentical, flexible machines, each capable of producing N different part types, was expressed in terms of the total production rates over all machines u n, as the convex hull of the set of points corresponding to all possible assignments of machines to part types First, it was shown that within each subset of assignments having the same number of machines assigned to each part type, there is a unique assignment that corresponds to an extreme point of Then, a procedure for generating all the extreme points and facets of was proposed The topology and geometry of, as expressed in terms of the total production rates u n, is not trivial, as the two-part type case might have misled us to believe Numerical experience showed that if the number of part types is less than four, the size of is smaller, if is expressed in terms of the total production rates u n than in terms of the production rates y nm If the number of part types is four or more, however, the size of is larger if is expressed in terms of the rates u n than in terms of the rates y nm References Caramanis, M and Liberopoulos, G, Perturbation Analysis for the Design of Flexible Manufacturing System Flow Controllers, Operations Research, Vol 40, No 6, pp (November December 1992) Gershwin, SB, Manufacturing Systems Engineering, Prentice-Hall, Englewood Cliffs, NJ (1994) Hung, Y-F and Wang, Q-Z, A New Formulation Technique for Alternative Material Planning An Approach for Semiconductor Bin Allocation Planning, Computers and Industrial Engineering, Vol 32, No 2, pp (April 1997) Kimemia, J and Gershwin, S B, An Algorithm for the Computer Control of a Flexible Manufacturing System, IIE Transactions, Vol 15, No 4, pp (December 1983) Lasserre, J B, New Capacity Sets for the Hedging Point Strategy, International Journal of Production Research, Vol 30, No 12, pp (December 1992) Leachman, R C and Carmon, T F, On Capacity Modeling for Production Planning with Alternative Machine Types, IIE Transactions, Vol 24, No 4, pp (September 1992) Liberopoulos, G and Caramanis, M, Infinitesimal Perturbation Analysis for Second Derivative Estimation and Design of Manufacturing Flow Controllers, Journal of Optimization Theory and Applications, Vol 81, No 2, pp (May 1994) Liberopoulos, G and Caramanis, M, Dynamics and Design of a Class of Parameterized Manufacturing Flow Controllers, IEEE Transactions on Automatic Control, Vol 40, No 6, pp (June 1995) Maimon, O Z and Gershwin, S B, Dynamic Scheduling and Routing for Flexible Manufacturing Systems that Have Unreliable Machines, Operations Research, Vol 36, No 2, pp (March April 1988) Maimon, O Z, Khmelnitsky, E, and Kogan, K, Optimal Flow Control in Manufacturing Systems: Production, Planning and Scheduling, Kluwer Academic Publishers, Dordrecht, The Netherlands (1998) Sethi, S P and Zhang, Q, Hierarchical Decision Making in Stochastic Manufacturing Systems, Birkhauser Boston, Cambridge, MA (1994)
Optimal Control of a Failure Prone Manufacturing System With the Possibility of Temporary Increase in Production Capacity Λ Lei Huang Jian-Qiang Hu Pi
Optimal Control of a Failure Prone Manufacturing System With the Possibility of Temporary Increase in Production Capacity Lei Huang Jian-Qiang Hu Pirooz Vakili Manufacturing Engineering Department Boston
More informationAbstrct. In this paper, we consider the problem of optimal flow control for a production system with one machine which is subject to failures and prod
Monotonicity of Optimal Flow Control for Failure Prone Production Systems 1 Jian-Qiang Hu 2 and Dong Xiang 3 Communicated By W.B. Gong 1 This work is partially supported by the national Science Foundation
More informationGLOBAL ANALYSIS OF PIECEWISE LINEAR SYSTEMS USING IMPACT MAPS AND QUADRATIC SURFACE LYAPUNOV FUNCTIONS
GLOBAL ANALYSIS OF PIECEWISE LINEAR SYSTEMS USING IMPACT MAPS AND QUADRATIC SURFACE LYAPUNOV FUNCTIONS Jorge M. Gonçalves, Alexandre Megretski y, Munther A. Dahleh y California Institute of Technology
More informationChapter 1. Preliminaries
Introduction This dissertation is a reading of chapter 4 in part I of the book : Integer and Combinatorial Optimization by George L. Nemhauser & Laurence A. Wolsey. The chapter elaborates links between
More informationApplications. Stephen J. Stoyan, Maged M. Dessouky*, and Xiaoqing Wang
Introduction to Large-Scale Linear Programming and Applications Stephen J. Stoyan, Maged M. Dessouky*, and Xiaoqing Wang Daniel J. Epstein Department of Industrial and Systems Engineering, University of
More informationA Parametric Simplex Algorithm for Linear Vector Optimization Problems
A Parametric Simplex Algorithm for Linear Vector Optimization Problems Birgit Rudloff Firdevs Ulus Robert Vanderbei July 9, 2015 Abstract In this paper, a parametric simplex algorithm for solving linear
More informationarxiv: v1 [cs.cc] 5 Dec 2018
Consistency for 0 1 Programming Danial Davarnia 1 and J. N. Hooker 2 1 Iowa state University davarnia@iastate.edu 2 Carnegie Mellon University jh38@andrew.cmu.edu arxiv:1812.02215v1 [cs.cc] 5 Dec 2018
More informationSupporting hyperplanes
Supporting hyperplanes General approach when using Lagrangian methods Lecture 1 homework Shadow prices A linear programming problem The simplex tableau Simple example with cycling The pivot rule being
More informationControl of manufacturing systems using state feedback and linear programming
Proceedings of the 44th IEEE Conference on Decision and Control, and the European Control Conference 2005 Seville, Spain, December 12-15, 2005 WeA15.2 Control of manufacturing systems using state feedback
More informationA Hierarchy of Suboptimal Policies for the Multi-period, Multi-echelon, Robust Inventory Problem
A Hierarchy of Suboptimal Policies for the Multi-period, Multi-echelon, Robust Inventory Problem Dimitris J. Bertsimas Dan A. Iancu Pablo A. Parrilo Sloan School of Management and Operations Research Center,
More informationCombinatorial optimization problems
Combinatorial optimization problems Heuristic Algorithms Giovanni Righini University of Milan Department of Computer Science (Crema) Optimization In general an optimization problem can be formulated as:
More informationBranch-and-cut Approaches for Chance-constrained Formulations of Reliable Network Design Problems
Branch-and-cut Approaches for Chance-constrained Formulations of Reliable Network Design Problems Yongjia Song James R. Luedtke August 9, 2012 Abstract We study solution approaches for the design of reliably
More informationMIT Manufacturing Systems Analysis Lecture 14-16
MIT 2.852 Manufacturing Systems Analysis Lecture 14-16 Line Optimization Stanley B. Gershwin Spring, 2007 Copyright c 2007 Stanley B. Gershwin. Line Design Given a process, find the best set of machines
More informationIntroduction to Operations Research Prof. G. Srinivasan Department of Management Studies Indian Institute of Technology, Madras
Introduction to Operations Research Prof. G. Srinivasan Department of Management Studies Indian Institute of Technology, Madras Module - 03 Simplex Algorithm Lecture 15 Infeasibility In this class, we
More informationBCOL RESEARCH REPORT 07.04
BCOL RESEARCH REPORT 07.04 Industrial Engineering & Operations Research University of California, Berkeley, CA 94720-1777 LIFTING FOR CONIC MIXED-INTEGER PROGRAMMING ALPER ATAMTÜRK AND VISHNU NARAYANAN
More informationLinear Programming Methods
Chapter 11 Linear Programming Methods 1 In this chapter we consider the linear programming approach to dynamic programming. First, Bellman s equation can be reformulated as a linear program whose solution
More information3.10 Lagrangian relaxation
3.10 Lagrangian relaxation Consider a generic ILP problem min {c t x : Ax b, Dx d, x Z n } with integer coefficients. Suppose Dx d are the complicating constraints. Often the linear relaxation and the
More informationA Shadow Simplex Method for Infinite Linear Programs
A Shadow Simplex Method for Infinite Linear Programs Archis Ghate The University of Washington Seattle, WA 98195 Dushyant Sharma The University of Michigan Ann Arbor, MI 48109 May 25, 2009 Robert L. Smith
More informationInteger Programming ISE 418. Lecture 12. Dr. Ted Ralphs
Integer Programming ISE 418 Lecture 12 Dr. Ted Ralphs ISE 418 Lecture 12 1 Reading for This Lecture Nemhauser and Wolsey Sections II.2.1 Wolsey Chapter 9 ISE 418 Lecture 12 2 Generating Stronger Valid
More informationType 1. Type 1 Type 2 Type 2 M 4 M 1 B 41 B 71 B 31 B 51 B 32 B 11 B 12 B 22 B 61 M 3 M 2 B 21
Design and Operation of Manufacturing Systems The Control-Point Policy by Stanley B. Gershwin Institute Massachusetts Technology of Cambridge, 02139 USA Massachusetts I I E Annual Conference Orlando, Florida,
More informationMIT Manufacturing Systems Analysis Lectures 19 21
MIT 2.852 Manufacturing Systems Analysis Lectures 19 21 Scheduling: Real-Time Control of Manufacturing Systems Stanley B. Gershwin Spring, 2007 Copyright c 2007 Stanley B. Gershwin. Definitions Events
More informationMath 341: Convex Geometry. Xi Chen
Math 341: Convex Geometry Xi Chen 479 Central Academic Building, University of Alberta, Edmonton, Alberta T6G 2G1, CANADA E-mail address: xichen@math.ualberta.ca CHAPTER 1 Basics 1. Euclidean Geometry
More informationLP Definition and Introduction to Graphical Solution Active Learning Module 2
LP Definition and Introduction to Graphical Solution Active Learning Module 2 J. René Villalobos and Gary L. Hogg Arizona State University Paul M. Griffin Georgia Institute of Technology Background Material
More informationOn the Approximate Linear Programming Approach for Network Revenue Management Problems
On the Approximate Linear Programming Approach for Network Revenue Management Problems Chaoxu Tong School of Operations Research and Information Engineering, Cornell University, Ithaca, New York 14853,
More information3.3 Easy ILP problems and totally unimodular matrices
3.3 Easy ILP problems and totally unimodular matrices Consider a generic ILP problem expressed in standard form where A Z m n with n m, and b Z m. min{c t x : Ax = b, x Z n +} (1) P(b) = {x R n : Ax =
More informationAssignment 1: From the Definition of Convexity to Helley Theorem
Assignment 1: From the Definition of Convexity to Helley Theorem Exercise 1 Mark in the following list the sets which are convex: 1. {x R 2 : x 1 + i 2 x 2 1, i = 1,..., 10} 2. {x R 2 : x 2 1 + 2ix 1x
More informationA Necessary Condition for Learning from Positive Examples
Machine Learning, 5, 101-113 (1990) 1990 Kluwer Academic Publishers. Manufactured in The Netherlands. A Necessary Condition for Learning from Positive Examples HAIM SHVAYTSER* (HAIM%SARNOFF@PRINCETON.EDU)
More informationSection Notes 9. Midterm 2 Review. Applied Math / Engineering Sciences 121. Week of December 3, 2018
Section Notes 9 Midterm 2 Review Applied Math / Engineering Sciences 121 Week of December 3, 2018 The following list of topics is an overview of the material that was covered in the lectures and sections
More informationOn the complexity of maximizing the minimum Shannon capacity in wireless networks by joint channel assignment and power allocation
On the complexity of maximizing the minimum Shannon capacity in wireless networks by joint channel assignment and power allocation Mikael Fallgren Royal Institute of Technology December, 2009 Abstract
More informationIntroduction to the Simplex Algorithm Active Learning Module 3
Introduction to the Simplex Algorithm Active Learning Module 3 J. René Villalobos and Gary L. Hogg Arizona State University Paul M. Griffin Georgia Institute of Technology Background Material Almost any
More informationNetwork Flows. 6. Lagrangian Relaxation. Programming. Fall 2010 Instructor: Dr. Masoud Yaghini
In the name of God Network Flows 6. Lagrangian Relaxation 6.3 Lagrangian Relaxation and Integer Programming Fall 2010 Instructor: Dr. Masoud Yaghini Integer Programming Outline Branch-and-Bound Technique
More informationLectures 6, 7 and part of 8
Lectures 6, 7 and part of 8 Uriel Feige April 26, May 3, May 10, 2015 1 Linear programming duality 1.1 The diet problem revisited Recall the diet problem from Lecture 1. There are n foods, m nutrients,
More informationOptimality Conditions for Constrained Optimization
72 CHAPTER 7 Optimality Conditions for Constrained Optimization 1. First Order Conditions In this section we consider first order optimality conditions for the constrained problem P : minimize f 0 (x)
More informationSingle-part-type, multiple stage systems
MIT 2.853/2.854 Introduction to Manufacturing Systems Single-part-type, multiple stage systems Stanley B. Gershwin Laboratory for Manufacturing and Productivity Massachusetts Institute of Technology Single-stage,
More informationScientific Computing: Optimization
Scientific Computing: Optimization Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 Course MATH-GA.2043 or CSCI-GA.2112, Spring 2012 March 8th, 2011 A. Donev (Courant Institute) Lecture
More informationDS-GA 1002 Lecture notes 0 Fall Linear Algebra. These notes provide a review of basic concepts in linear algebra.
DS-GA 1002 Lecture notes 0 Fall 2016 Linear Algebra These notes provide a review of basic concepts in linear algebra. 1 Vector spaces You are no doubt familiar with vectors in R 2 or R 3, i.e. [ ] 1.1
More informationTMA 4180 Optimeringsteori KARUSH-KUHN-TUCKER THEOREM
TMA 4180 Optimeringsteori KARUSH-KUHN-TUCKER THEOREM H. E. Krogstad, IMF, Spring 2012 Karush-Kuhn-Tucker (KKT) Theorem is the most central theorem in constrained optimization, and since the proof is scattered
More information...(iii), x 2 Example 7: Geetha Perfume Company produces both perfumes and body spray from two flower extracts F 1. The following data is provided:
The LP formulation is Linear Programming: Graphical Method Maximize, Z = 2x + 7x 2 Subject to constraints, 2x + x 2 200...(i) x 75...(ii) x 2 00...(iii) where x, x 2 ³ 0 Example 7: Geetha Perfume Company
More informationSemidefinite and Second Order Cone Programming Seminar Fall 2012 Project: Robust Optimization and its Application of Robust Portfolio Optimization
Semidefinite and Second Order Cone Programming Seminar Fall 2012 Project: Robust Optimization and its Application of Robust Portfolio Optimization Instructor: Farid Alizadeh Author: Ai Kagawa 12/12/2012
More informationIBM Research Report. Stochasic Unit Committment Problem. Julio Goez Lehigh University. James Luedtke University of Wisconsin
RC24713 (W0812-119) December 23, 2008 Mathematics IBM Research Report Stochasic Unit Committment Problem Julio Goez Lehigh University James Luedtke University of Wisconsin Deepak Rajan IBM Research Division
More informationMAT-INF4110/MAT-INF9110 Mathematical optimization
MAT-INF4110/MAT-INF9110 Mathematical optimization Geir Dahl August 20, 2013 Convexity Part IV Chapter 4 Representation of convex sets different representations of convex sets, boundary polyhedra and polytopes:
More informationAn Integer Cutting-Plane Procedure for the Dantzig-Wolfe Decomposition: Theory
An Integer Cutting-Plane Procedure for the Dantzig-Wolfe Decomposition: Theory by Troels Martin Range Discussion Papers on Business and Economics No. 10/2006 FURTHER INFORMATION Department of Business
More informationIE 5531: Engineering Optimization I
IE 5531: Engineering Optimization I Lecture 15: Nonlinear optimization Prof. John Gunnar Carlsson November 1, 2010 Prof. John Gunnar Carlsson IE 5531: Engineering Optimization I November 1, 2010 1 / 24
More informationThe continuous knapsack set
The continuous knapsack set Sanjeeb Dash IBM Research sanjeebd@us.ibm.com Oktay Günlük IBM Research gunluk@us.ibm.com January 31, 2014 Laurence Wolsey Core laurence.wolsey@uclouvain.be Abstract We study
More informationTopic 6: Projected Dynamical Systems
Topic 6: Projected Dynamical Systems John F. Smith Memorial Professor and Director Virtual Center for Supernetworks Isenberg School of Management University of Massachusetts Amherst, Massachusetts 01003
More informationSolving the MWT. Recall the ILP for the MWT. We can obtain a solution to the MWT problem by solving the following ILP:
Solving the MWT Recall the ILP for the MWT. We can obtain a solution to the MWT problem by solving the following ILP: max subject to e i E ω i x i e i C E x i {0, 1} x i C E 1 for all critical mixed cycles
More informationInformation in Aloha Networks
Achieving Proportional Fairness using Local Information in Aloha Networks Koushik Kar, Saswati Sarkar, Leandros Tassiulas Abstract We address the problem of attaining proportionally fair rates using Aloha
More informationFinite Dimensional Optimization Part III: Convex Optimization 1
John Nachbar Washington University March 21, 2017 Finite Dimensional Optimization Part III: Convex Optimization 1 1 Saddle points and KKT. These notes cover another important approach to optimization,
More informationThe Simplex Method: An Example
The Simplex Method: An Example Our first step is to introduce one more new variable, which we denote by z. The variable z is define to be equal to 4x 1 +3x 2. Doing this will allow us to have a unified
More informationProjection in Logic, CP, and Optimization
Projection in Logic, CP, and Optimization John Hooker Carnegie Mellon University Workshop on Logic and Search Melbourne, 2017 Projection as a Unifying Concept Projection is a fundamental concept in logic,
More informationNo EFFICIENT LINE SEARCHING FOR CONVEX FUNCTIONS. By E. den Boef, D. den Hertog. May 2004 ISSN
No. 4 5 EFFICIENT LINE SEARCHING FOR CONVEX FUNCTIONS y E. den oef, D. den Hertog May 4 ISSN 94-785 Efficient Line Searching for Convex Functions Edgar den oef Dick den Hertog 3 Philips Research Laboratories,
More informationLinear Programming and the Simplex method
Linear Programming and the Simplex method Harald Enzinger, Michael Rath Signal Processing and Speech Communication Laboratory Jan 9, 2012 Harald Enzinger, Michael Rath Jan 9, 2012 page 1/37 Outline Introduction
More informationFacets for Node-Capacitated Multicut Polytopes from Path-Block Cycles with Two Common Nodes
Facets for Node-Capacitated Multicut Polytopes from Path-Block Cycles with Two Common Nodes Michael M. Sørensen July 2016 Abstract Path-block-cycle inequalities are valid, and sometimes facet-defining,
More informationA Review of Linear Programming
A Review of Linear Programming Instructor: Farid Alizadeh IEOR 4600y Spring 2001 February 14, 2001 1 Overview In this note we review the basic properties of linear programming including the primal simplex
More informationVector Spaces. Addition : R n R n R n Scalar multiplication : R R n R n.
Vector Spaces Definition: The usual addition and scalar multiplication of n-tuples x = (x 1,..., x n ) R n (also called vectors) are the addition and scalar multiplication operations defined component-wise:
More informationSome Background Math Notes on Limsups, Sets, and Convexity
EE599 STOCHASTIC NETWORK OPTIMIZATION, MICHAEL J. NEELY, FALL 2008 1 Some Background Math Notes on Limsups, Sets, and Convexity I. LIMITS Let f(t) be a real valued function of time. Suppose f(t) converges
More informationResearch Article Convex Polyhedron Method to Stability of Continuous Systems with Two Additive Time-Varying Delay Components
Applied Mathematics Volume 202, Article ID 689820, 3 pages doi:0.55/202/689820 Research Article Convex Polyhedron Method to Stability of Continuous Systems with Two Additive Time-Varying Delay Components
More informationChapter 3: Discrete Optimization Integer Programming
Chapter 3: Discrete Optimization Integer Programming Edoardo Amaldi DEIB Politecnico di Milano edoardo.amaldi@polimi.it Website: http://home.deib.polimi.it/amaldi/opt-16-17.shtml Academic year 2016-17
More informationLecture 9 Monotone VIs/CPs Properties of cones and some existence results. October 6, 2008
Lecture 9 Monotone VIs/CPs Properties of cones and some existence results October 6, 2008 Outline Properties of cones Existence results for monotone CPs/VIs Polyhedrality of solution sets Game theory:
More informationA Lower Bound on the Split Rank of Intersection Cuts
A Lower Bound on the Split Rank of Intersection Cuts Santanu S. Dey H. Milton Stewart School of Industrial and Systems Engineering, Georgia Institute of Technology. 200 Outline Introduction: split rank,
More informationSOME RESOURCE ALLOCATION PROBLEMS
SOME RESOURCE ALLOCATION PROBLEMS A Dissertation Presented to the Faculty of the Graduate School of Cornell University in Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy
More informationOptimization Problems with Probabilistic Constraints
Optimization Problems with Probabilistic Constraints R. Henrion Weierstrass Institute Berlin 10 th International Conference on Stochastic Programming University of Arizona, Tucson Recommended Reading A.
More informationSubject: Optimal Control Assignment-1 (Related to Lecture notes 1-10)
Subject: Optimal Control Assignment- (Related to Lecture notes -). Design a oil mug, shown in fig., to hold as much oil possible. The height and radius of the mug should not be more than 6cm. The mug must
More informationJørgen Tind, Department of Statistics and Operations Research, University of Copenhagen, Universitetsparken 5, 2100 Copenhagen O, Denmark.
DUALITY THEORY Jørgen Tind, Department of Statistics and Operations Research, University of Copenhagen, Universitetsparken 5, 2100 Copenhagen O, Denmark. Keywords: Duality, Saddle point, Complementary
More informationPartial job order for solving the two-machine flow-shop minimum-length problem with uncertain processing times
Preprints of the 13th IFAC Symposium on Information Control Problems in Manufacturing, Moscow, Russia, June 3-5, 2009 Fr-A2.3 Partial job order for solving the two-machine flow-shop minimum-length problem
More informationLinear Programming Inverse Projection Theory Chapter 3
1 Linear Programming Inverse Projection Theory Chapter 3 University of Chicago Booth School of Business Kipp Martin September 26, 2017 2 Where We Are Headed We want to solve problems with special structure!
More informationLinear programs, convex polyhedra, extreme points
MVE165/MMG631 Extreme points of convex polyhedra; reformulations; basic feasible solutions; the simplex method Ann-Brith Strömberg 2015 03 27 Linear programs, convex polyhedra, extreme points A linear
More informationDistributed Real-Time Control Systems. Lecture Distributed Control Linear Programming
Distributed Real-Time Control Systems Lecture 13-14 Distributed Control Linear Programming 1 Linear Programs Optimize a linear function subject to a set of linear (affine) constraints. Many problems can
More information3.7 Strong valid inequalities for structured ILP problems
3.7 Strong valid inequalities for structured ILP problems By studying the problem structure, we can derive strong valid inequalities yielding better approximations of conv(x ) and hence tighter bounds.
More informationLinear Algebra Review: Linear Independence. IE418 Integer Programming. Linear Algebra Review: Subspaces. Linear Algebra Review: Affine Independence
Linear Algebra Review: Linear Independence IE418: Integer Programming Department of Industrial and Systems Engineering Lehigh University 21st March 2005 A finite collection of vectors x 1,..., x k R n
More informationTHE EXISTENCE AND USEFULNESS OF EQUALITY CUTS IN THE MULTI-DEMAND MULTIDIMENSIONAL KNAPSACK PROBLEM LEVI DELISSA. B.S., Kansas State University, 2014
THE EXISTENCE AND USEFULNESS OF EQUALITY CUTS IN THE MULTI-DEMAND MULTIDIMENSIONAL KNAPSACK PROBLEM by LEVI DELISSA B.S., Kansas State University, 2014 A THESIS submitted in partial fulfillment of the
More informationMathematical Foundations -1- Constrained Optimization. Constrained Optimization. An intuitive approach 2. First Order Conditions (FOC) 7
Mathematical Foundations -- Constrained Optimization Constrained Optimization An intuitive approach First Order Conditions (FOC) 7 Constraint qualifications 9 Formal statement of the FOC for a maximum
More informationNecessary and Sufficient Conditions for Reachability on a Simplex
Necessary and Sufficient Conditions for Reachability on a Simplex Bartek Roszak a, Mireille E. Broucke a a Edward S. Rogers Sr. Department of Electrical and Computer Engineering, University of Toronto,
More informationA Generalized Uncertainty Principle and Sparse Representation in Pairs of Bases
2558 IEEE TRANSACTIONS ON INFORMATION THEORY, VOL 48, NO 9, SEPTEMBER 2002 A Generalized Uncertainty Principle Sparse Representation in Pairs of Bases Michael Elad Alfred M Bruckstein Abstract An elementary
More informationRandomized Coordinate Descent Methods on Optimization Problems with Linearly Coupled Constraints
Randomized Coordinate Descent Methods on Optimization Problems with Linearly Coupled Constraints By I. Necoara, Y. Nesterov, and F. Glineur Lijun Xu Optimization Group Meeting November 27, 2012 Outline
More informationDO NOT OPEN THIS QUESTION BOOKLET UNTIL YOU ARE TOLD TO DO SO
QUESTION BOOKLET EECS 227A Fall 2009 Midterm Tuesday, Ocotober 20, 11:10-12:30pm DO NOT OPEN THIS QUESTION BOOKLET UNTIL YOU ARE TOLD TO DO SO You have 80 minutes to complete the midterm. The midterm consists
More information56:171 Operations Research Midterm Exam - October 26, 1989 Instructor: D.L. Bricker
56:171 Operations Research Midterm Exam - October 26, 1989 Instructor: D.L. Bricker Answer all of Part One and two (of the four) problems of Part Two Problem: 1 2 3 4 5 6 7 8 TOTAL Possible: 16 12 20 10
More informationLifting for conic mixed-integer programming
Math. Program., Ser. A DOI 1.17/s117-9-282-9 FULL LENGTH PAPER Lifting for conic mixed-integer programming Alper Atamtürk Vishnu Narayanan Received: 13 March 28 / Accepted: 28 January 29 The Author(s)
More informationDuality Theory of Constrained Optimization
Duality Theory of Constrained Optimization Robert M. Freund April, 2014 c 2014 Massachusetts Institute of Technology. All rights reserved. 1 2 1 The Practical Importance of Duality Duality is pervasive
More informationJordan Journal of Mathematics and Statistics (JJMS) 8(3), 2015, pp THE NORM OF CERTAIN MATRIX OPERATORS ON NEW DIFFERENCE SEQUENCE SPACES
Jordan Journal of Mathematics and Statistics (JJMS) 8(3), 2015, pp 223-237 THE NORM OF CERTAIN MATRIX OPERATORS ON NEW DIFFERENCE SEQUENCE SPACES H. ROOPAEI (1) AND D. FOROUTANNIA (2) Abstract. The purpose
More informationIntroduction to Mathematical Programming IE406. Lecture 21. Dr. Ted Ralphs
Introduction to Mathematical Programming IE406 Lecture 21 Dr. Ted Ralphs IE406 Lecture 21 1 Reading for This Lecture Bertsimas Sections 10.2, 10.3, 11.1, 11.2 IE406 Lecture 21 2 Branch and Bound Branch
More information3.7 Cutting plane methods
3.7 Cutting plane methods Generic ILP problem min{ c t x : x X = {x Z n + : Ax b} } with m n matrix A and n 1 vector b of rationals. According to Meyer s theorem: There exists an ideal formulation: conv(x
More informationA Simpler and Tighter Redundant Klee-Minty Construction
A Simpler and Tighter Redundant Klee-Minty Construction Eissa Nematollahi Tamás Terlaky October 19, 2006 Abstract By introducing redundant Klee-Minty examples, we have previously shown that the central
More informationMATHEMATICAL PROGRAMMING I
MATHEMATICAL PROGRAMMING I Books There is no single course text, but there are many useful books, some more mathematical, others written at a more applied level. A selection is as follows: Bazaraa, Jarvis
More informationOn the Polyhedral Structure of a Multi Item Production Planning Model with Setup Times
CORE DISCUSSION PAPER 2000/52 On the Polyhedral Structure of a Multi Item Production Planning Model with Setup Times Andrew J. Miller 1, George L. Nemhauser 2, and Martin W.P. Savelsbergh 2 November 2000
More informationSome Properties of the Augmented Lagrangian in Cone Constrained Optimization
MATHEMATICS OF OPERATIONS RESEARCH Vol. 29, No. 3, August 2004, pp. 479 491 issn 0364-765X eissn 1526-5471 04 2903 0479 informs doi 10.1287/moor.1040.0103 2004 INFORMS Some Properties of the Augmented
More informationA polynomial time algorithm to solve the single-item capacitated lot sizing problem with minimum order quantities and concave costs
A polynomial time algorithm to solve the single-item capacitated lot sizing problem with minimum order quantities and concave costs Bertrand Hellion, Bernard Penz, Fabien Mangione aboratoire G-SCOP n o
More informationRegular finite Markov chains with interval probabilities
5th International Symposium on Imprecise Probability: Theories and Applications, Prague, Czech Republic, 2007 Regular finite Markov chains with interval probabilities Damjan Škulj Faculty of Social Sciences
More informationOptimization of Quadratic Forms: NP Hard Problems : Neural Networks
1 Optimization of Quadratic Forms: NP Hard Problems : Neural Networks Garimella Rama Murthy, Associate Professor, International Institute of Information Technology, Gachibowli, HYDERABAD, AP, INDIA ABSTRACT
More informationSingle-part-type, multiple stage systems. Lecturer: Stanley B. Gershwin
Single-part-type, multiple stage systems Lecturer: Stanley B. Gershwin Flow Line... also known as a Production or Transfer Line. M 1 B 1 M 2 B 2 M 3 B 3 M 4 B 4 M 5 B 5 M 6 Machine Buffer Machines are
More informationChapter 2: Linear Programming Basics. (Bertsimas & Tsitsiklis, Chapter 1)
Chapter 2: Linear Programming Basics (Bertsimas & Tsitsiklis, Chapter 1) 33 Example of a Linear Program Remarks. minimize 2x 1 x 2 + 4x 3 subject to x 1 + x 2 + x 4 2 3x 2 x 3 = 5 x 3 + x 4 3 x 1 0 x 3
More informationMultidisciplinary System Design Optimization (MSDO)
Multidisciplinary System Design Optimization (MSDO) Numerical Optimization II Lecture 8 Karen Willcox 1 Massachusetts Institute of Technology - Prof. de Weck and Prof. Willcox Today s Topics Sequential
More informationInformation Structures Preserved Under Nonlinear Time-Varying Feedback
Information Structures Preserved Under Nonlinear Time-Varying Feedback Michael Rotkowitz Electrical Engineering Royal Institute of Technology (KTH) SE-100 44 Stockholm, Sweden Email: michael.rotkowitz@ee.kth.se
More information5 Handling Constraints
5 Handling Constraints Engineering design optimization problems are very rarely unconstrained. Moreover, the constraints that appear in these problems are typically nonlinear. This motivates our interest
More information15-850: Advanced Algorithms CMU, Fall 2018 HW #4 (out October 17, 2018) Due: October 28, 2018
15-850: Advanced Algorithms CMU, Fall 2018 HW #4 (out October 17, 2018) Due: October 28, 2018 Usual rules. :) Exercises 1. Lots of Flows. Suppose you wanted to find an approximate solution to the following
More informationSummary of the simplex method
MVE165/MMG631,Linear and integer optimization with applications The simplex method: degeneracy; unbounded solutions; starting solutions; infeasibility; alternative optimal solutions Ann-Brith Strömberg
More informationLP Relaxations of Mixed Integer Programs
LP Relaxations of Mixed Integer Programs John E. Mitchell Department of Mathematical Sciences RPI, Troy, NY 12180 USA February 2015 Mitchell LP Relaxations 1 / 29 LP Relaxations LP relaxations We want
More informationRevisiting the Hamiltonian p-median problem: a new formulation on directed graphs and a branch-and-cut algorithm
Revisiting the Hamiltonian p-median problem: a new formulation on directed graphs and a branch-and-cut algorithm Tolga Bektaş 1, Luís Gouveia 2, Daniel Santos 2 1 Centre for Operational Research, Management
More informationComputer Science 385 Analysis of Algorithms Siena College Spring Topic Notes: Limitations of Algorithms
Computer Science 385 Analysis of Algorithms Siena College Spring 2011 Topic Notes: Limitations of Algorithms We conclude with a discussion of the limitations of the power of algorithms. That is, what kinds
More informationColumn Generation in Integer Programming with Applications in Multicriteria Optimization
Column Generation in Integer Programming with Applications in Multicriteria Optimization Matthias Ehrgott Department of Engineering Science The University of Auckland, New Zealand email: m.ehrgott@auckland.ac.nz
More information