1 / 146

The Simplex Algorithm

The Simplex Algorithm. 虞台文. 大同大學資工所 智慧型多媒體研究室. Content. Basic Feasible Solutions The Geometry of Linear Programs Moving From Bfs to Bfs Organization of a Tableau Choosing a Profitable Column Pivot Selection Anticycling Algorithm The Two-Phase Simplex Algorithm. Linear Programming.

Download Presentation

The Simplex Algorithm

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. The Simplex Algorithm 虞台文 大同大學資工所 智慧型多媒體研究室

  2. Content • Basic Feasible Solutions • The Geometry of Linear Programs • Moving From Bfs to Bfs • Organization of a Tableau • Choosing a Profitable Column • Pivot Selection • Anticycling Algorithm • The Two-Phase Simplex Algorithm

  3. Linear Programming Basic Feasible Solutions 大同大學資工所 智慧型多媒體研究室

  4. The Goal of the Simplex Algorithm Any LP can be converted into the standard form

  5. Convex Polytope defines a convex set F • Bounded • Unbounded • Empty F may be

  6. The Basic Idea of Simplex Algorithm Finding optimum by moving around the corner of the convex polytopein cost descent sense basic feasible solution at a corner

  7. The Basic Idea of Simplex Algorithm How to find an initial feasible solution? Finding optimum by moving around the corner of the convex polytopein cost descent sense How to move from corner to corner? basic feasible solution at a corner

  8. Assumption 1 Assume that A is of rank m. Independent There is a basis has an inverse.

  9. Basic Solution The basis solution corresponding to B is: The basis solution may be infeasible.

  10. x1 x2 x3 x4 x5 x6 x7 Example

  11. x1 x2 x3 x4 x5 x6 x7 Example feasible

  12. x1 x2 x3 x4 x5 x6 x7 Example infeasible

  13. Basic Feasible Solution defines an nonempty convex set F If a basic solution is in F, it is called a basic feasible solution (bfs).

  14. The Existence of BFS’s • Bounded • Unbounded • Empty F may be defines an nonempty convex set F If F is nonempty, at least one bfs.

  15. The Existence of BFS’s • Bounded • Unbounded • Empty F may be defines an nonempty convex set F How to find it? If F is nonempty, at least one bfs.

  16. Assumptions defines an nonempty convex set F • A is of rank m. • F is nonempty. • c’x is bounded below for xF.

  17. Assumptions Are all bfs’s the vertices of the convex polytope defined by F? defines an nonempty convex set F • A is of rank m. • F is nonempty. • c’x is bounded below for xF. Ensure there is a bounded solution. There exist bfs’s.

  18. The Simplex Alogrithm The Geometry of Linear Programs 大同大學資工所 智慧型多媒體研究室

  19. Linear Subspaces of Rd • S  Rdis a subspace of Rdif it is closed under vector addition and scalar multiplication. • S defined below is a subspace of Rd. a set of homogenous linear equations

  20. Dimensions

  21. an affine subspace an affine subspace Affine subspaces A translated linear subspace. E.g., a set of nonhomogenous linear equations

  22. an affine subspace Dimensions

  23. Subsets of Rd • The following subsets are notsubspace or affine subspace. A line segment The first quadrant A halfspace

  24. The dimension of any subset of Rd is the smallest dimension of any affine subspace which contains it. Dimensions • The following subsets are notsubspace or affine subspace. Dim = 2 Dim = 2 Dim = 1 A line segment The first quadrant A halfspace

  25. The Feasible Spaces of LP defines an nonempty convex set F • A is of rank m. • F is nonempty. • c’x is bounded below for xF.

  26. Fact: Halfspaces are convex. Hyperplane/Halfspace An affine subspace of Rd of dimension d1 is called a hyperplane. Ahyperplane defines two closed halfspaces:

  27. Fact: Halfspaces are convex. Convex Polytopes The intersection of a finite number of halfspaces, when it is bounded and nonempty is called a convex polytope, or simply a polytope.

  28. x3 x1 x2 Example  (0, 0, 3)  (1, 0, 3)   (0, 1, 3)  (2, 0, 2)   (2, 0, 0) (0, 0, 0)  (0, 2, 0) (2, 2, 0)

  29. More on Polytopes • Geometric views of a polytope • The convex hull of a finite set of points. • The intersection of a finite number of halfspaces, when it is bounded and nonempty, i.e., defined by • The algebraic view of a polytope Any relation? The feasible space defined by LP (in standard form).

  30. LP LP  Polytopes xn x1 x2 xnm xnm+1 Assume n  m columns m columns A can always be in this form if rank(A)=m.

  31. LP LP  Polytopes xn x1 x2 xnm xnm+1 Assume n  m columns m columns

  32. LP LP  Polytopes 0

  33. Polytopes  LP m inequalities m equalities

  34. Slack variables Polytopes  LP m inequalities m equalities

  35. Slack variables Polytopes  LP

  36. Slack variables Polytopes and LP H H I

  37. Polytopes and LP The answer

  38. Polytopes and LP The answer

  39. Polytopes & F of LP defines a polytope Are there any relations? defines a feasible set

  40. Polytopes & F of LP defines a polytope Some points in P are vertices. Are there any relations? Some points in F are bfs’s. defines a feasible set

  41. a convex polytope a feasible set of the corresponding LP is a vertex of P is a bfs. Theorem 1

  42. Theorem 1 Pf) “” See textbook.

  43. is a vertex, then it cannot be the strict combination of points of P, i.e., Theorem 1 Pf) “” Fact:

  44. is a vertex, then it cannot be the strict combination of points of P, i.e., Theorem 1 Pf) “” Fact:

  45. Let Suppose not. Then, some dj 0. Theorem 1 Pf) “” Fact: We want to show at most m nonzero xi’s. We want to show that vectors in B are linearly independent. Define B = {Aj: xj > 0, 1  j  n} >0 if  is sufficiently small.

  46. Let Suppose not. Then, some dj 0. Define Theorem 1  is sufficiently small. Pf) “” Fact: We want to show at most m nonzero xi’s. We want to show that vectors in B are linearly independent. Define B = {Aj: xj > 0, 1  j  n} >0 if  is sufficiently small.

  47. Let Theorem 1 Pf) “” Fact: We want to show at most m nonzero xi’s. Define B = {Aj: xj > 0, 1  j  n} Aj are linearly independent. |B|  m. Since rank(A) = m, we can always augment B to include m linearly independent vectors. Using B to form basic columns renders x a bfs. 

  48. Let Two different , their corresponding B and B’ must be different. Discussion If |B| < m, there may be many ways to augment to m linearly independent vectors. Two different B and B’ may corresponds to the same bfs. Pf) “” Fact: We want to show at most m nonzero xi’s. Define B = {Aj: xj > 0, 1  j  n} Aj are linearly independent. |B|  m. Since rank(A) = m, we can always augment B to include m linearly independent vectors. Using B to form basic columns renders x a bfs. 

  49. x3  (0, 0, 3)  (1, 0, 3)   (0, 1, 3)  (2, 0, 2)   (2, 0, 0) (0, 0, 0) x1  x2 (0, 2, 0) (2, 2, 0) Example

  50. Example Corresponding F of LP

More Related