We will see what is meant by standard form very shortly

Similar documents
Duality # Second iteration for HW problem. Recall our LP example problem we have been working on, in equality form, is given below.

20 MATHEMATICS POLYNOMIALS

Can the Phase I problem be unfeasible or unbounded? -No

Using QM for Windows. Using QM for Windows. Using QM for Windows LEARNING OBJECTIVES. Solving Flair Furniture s LP Problem

Review of Calculus, cont d

Lesson 1: Quadratic Equations

How do you know you have SLE?

MORE FUNCTION GRAPHING; OPTIMIZATION. (Last edited October 28, 2013 at 11:09pm.)

The area under the graph of f and above the x-axis between a and b is denoted by. f(x) dx. π O

Identify graphs of linear inequalities on a number line.

Chapter 14. Matrix Representations of Linear Transformations

approaches as n becomes larger and larger. Since e > 1, the graph of the natural exponential function is as below

SUMMER KNOWHOW STUDY AND LEARNING CENTRE

Pre-Session Review. Part 1: Basic Algebra; Linear Functions and Graphs

Chapter 3 Solving Nonlinear Equations

Unit 1 Exponentials and Logarithms

Review of basic calculus

Definite integral. Mathematics FRDIS MENDELU

Before we can begin Ch. 3 on Radicals, we need to be familiar with perfect squares, cubes, etc. Try and do as many as you can without a calculator!!!

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies

Z b. f(x)dx. Yet in the above two cases we know what f(x) is. Sometimes, engineers want to calculate an area by computing I, but...

Administrivia CSE 190: Reinforcement Learning: An Introduction

AQA Further Pure 1. Complex Numbers. Section 1: Introduction to Complex Numbers. The number system

Jim Lambers MAT 169 Fall Semester Lecture 4 Notes

19 Optimal behavior: Game theory

Definite integral. Mathematics FRDIS MENDELU. Simona Fišnarová (Mendel University) Definite integral MENDELU 1 / 30

Theoretical foundations of Gaussian quadrature

MAA 4212 Improper Integrals

Math 4310 Solutions to homework 1 Due 9/1/16

How do we solve these things, especially when they get complicated? How do we know when a system has a solution, and when is it unique?

Unit #9 : Definite Integral Properties; Fundamental Theorem of Calculus

The Regulated and Riemann Integrals

Matrix Solution to Linear Equations and Markov Chains

(e) if x = y + z and a divides any two of the integers x, y, or z, then a divides the remaining integer

Math 520 Final Exam Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008

Lecture 19: Continuous Least Squares Approximation

Chapter 3 Polynomials

Module 6: LINEAR TRANSFORMATIONS

Improper Integrals, and Differential Equations

PART 1 MULTIPLE CHOICE Circle the appropriate response to each of the questions below. Each question has a value of 1 point.

Equations and Inequalities

Section 6.1 Definite Integral

Suppose we want to find the area under the parabola and above the x axis, between the lines x = 2 and x = -2.

1. Extend QR downwards to meet the x-axis at U(6, 0). y

1.2. Linear Variable Coefficient Equations. y + b "! = a y + b " Remark: The case b = 0 and a non-constant can be solved with the same idea as above.

ENGI 3424 Engineering Mathematics Five Tutorial Examples of Partial Fractions

63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1

38 Riemann sums and existence of the definite integral.

Problem Set 3

p-adic Egyptian Fractions

M344 - ADVANCED ENGINEERING MATHEMATICS

4.4 Areas, Integrals and Antiderivatives

In-Class Problems 2 and 3: Projectile Motion Solutions. In-Class Problem 2: Throwing a Stone Down a Hill

Continuous Random Variables

Chapter 1 Cumulative Review

Summary: Method of Separation of Variables

Improper Integrals. Type I Improper Integrals How do we evaluate an integral such as

THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS.

Review of Gaussian Quadrature method

Advanced Calculus: MATH 410 Notes on Integrals and Integrability Professor David Levermore 17 October 2004

7.2 The Definite Integral

Lecture 14: Quadrature

Operations with Polynomials

USA Mathematical Talent Search Round 1 Solutions Year 21 Academic Year

Orthogonal Polynomials

Linear Inequalities. Work Sheet 1

Math 270A: Numerical Linear Algebra

Numerical integration

Minimal DFA. minimal DFA for L starting from any other

Bridging the gap: GCSE AS Level

Numerical Analysis: Trapezoidal and Simpson s Rule

Integral points on the rational curve

Linearly Similar Polynomials

Tech. Rpt. # UMIACS-TR-99-31, Institute for Advanced Computer Studies, University of Maryland, College Park, MD 20742, June 3, 1999.

Riemann is the Mann! (But Lebesgue may besgue to differ.)

Reinforcement learning II

Matrices and Determinants

Abstract inner product spaces

A-Level Mathematics Transition Task (compulsory for all maths students and all further maths student)

Heteroclinic cycles in coupled cell systems

AMATH 731: Applied Functional Analysis Fall Some basics of integral equations

Downloaded from

Elementary Linear Algebra

Polynomials and Division Theory

UNIT-2 POLYNOMIALS Downloaded From: [Year] UNIT-2 POLYNOMIALS

Math& 152 Section Integration by Parts

DATA Search I 魏忠钰. 复旦大学大数据学院 School of Data Science, Fudan University. March 7 th, 2018

CHAPTER 4a. ROOTS OF EQUATIONS

f(x) dx, If one of these two conditions is not met, we call the integral improper. Our usual definition for the value for the definite integral

Uninformed Search Lecture 4

USA Mathematical Talent Search Round 1 Solutions Year 25 Academic Year

Finite Automata Theory and Formal Languages TMV027/DIT321 LP4 2018

0.1 Chapters 1: Limits and continuity

Bases for Vector Spaces

For the percentage of full time students at RCC the symbols would be:

5.7 Improper Integrals

Expectation and Variance

UniversitaireWiskundeCompetitie. Problem 2005/4-A We have k=1. Show that for every q Q satisfying 0 < q < 1, there exists a finite subset K N so that

NUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.

Math 113 Exam 2 Practice

Transcription:

THEOREM: For fesible liner progrm in its stndrd form, the optimum vlue of the objective over its nonempty fesible region is () either unbounded or (b) is chievble t lest t one extreme point of the fesible region. We will see wht is ment by stndrd form very shortly More generlly, the bove theorem nd the grphs we sw tell us tht every LP is in exctly one of the following sttes: 1. Fesible with unique optimum solution (cluse b of theorem) 2. Fesible with infinitely mny optim (cluse b of theorem) 3. Fesible, with no optimum solution becuse the objective is unbounded (cluse of theorem) 4. Infesible, nd hence with no optimum solution 2018, Jynt Rjgopl

Assume for now tht we hve fesible LP nd tht the objective is bounded (i.e., we re in Sttes 1 or 2). Then the following re true: A. The LP hs t lest one optiml corner point (or extreme point) 1. If in Stte 1, exctly one extreme point is optiml 2. If in Stte 2, t lest two djcent (neighboring) extreme points re optiml B. The number of extreme points is finite C. If the objective function t some extreme point is s good s or better thn t ll of its djcent extreme point, then this extreme point is optiml for the LP This suggest simple lgorithmic pproch to finding the optimum of fesible nd bounded LP 2018, Jynt Rjgopl

STEP 0 (Initiliztion): Find n initil extreme point nd mke it the current cndidte (if one cnnot be found the LP is in stte 4, i.e. it is infesible so STOP). STEP 1(Stopping Criterion Check): Is the objective t the current extreme point t lest s good or better thn it is t ll of its djcent (neighboring) extreme points? If so this must be the optiml extreme point (vi C) so STOP. If not, go to Step 2. STEP 2(Itertive Step): One (or more ccurtely, t lest one ) of the djcent extreme points is better so mke it the current cndidte. Then return to Step 1. QUESTIONS: 1) How to find n initil extreme point? 2) Wht is the lgebric chrcteriztion of n extreme point? 3) Wht is the lgebric chrcteriztion of djcent extreme points, i.e., how to move from n extreme point to one of its neighbors (in Step 2)? 2018, Jynt Rjgopl

The SIMPLEX Method - Development Every LP is in exctly one of the following sttes: 1. Fesible with unique optimum solution. 2. Fesible with infinitely mny optim. 3. Fesible, with no optimum solution (becuse the objective is unbounded). 4. Infesible, nd hence with no optimum solution. Assume we re in Sttes 1 or 2. Then the following re true: A. The LP hs t lest one optiml corner point (or extreme point). 1) If in Stte 1, exctly one extreme point is optiml. 2) If in Stte 2, t lest two djcent (neighboring) extreme points re optiml. B. The number of extreme points is finite. C. If the objective function t some extreme point is s good s or better thn t ll of its djcent extreme point, then this extreme point is optiml for the LP.

It follows then tht we cn use the following itertive lgorithmic procedure: STEP 0 (Initiliztion): Find n initil extreme point nd mke it the current cndidte (if one cnnot be found the LP is in stte 4, i.e. it is infesible so STOP). STEP 1 (Stopping Criterion Check): Is the objective t the current extreme point t lest s good or better thn it is t ll of its djcent (neighboring) extreme points? If so this must be the optiml extreme point (vi C) so STOP. If not, go to Step 2. STEP 2 (Itertive Step): At lest one of the djcent extreme points is better so mke it the current cndidte nd go to Step 1. QUESTIONS: 1) How to find n initil extreme point? 2) Wht is the lgebric chrcteriztion of n extreme point? 3) Wht is the lgebric chrcteriztion of djcent extreme points, i.e., how to move from n extreme point to one of its neighbors (in Step 2)?

Stndrd Form of Liner Progrm Three conditions must be met: 1) All constrints must be stted s equlities of the form g k (x)=b k, where g k (x) is liner function of x. 2) The RHS for ech constrint must be nonnegtive, i.e., ll b k 0 3) All vribles must be nonnegtive, i.e., x i 0 for ll i. Thus the LP (with m constrints nd n vribles) looks like Min (or Mx) c 1 x 1 + c 2 x 2 + + c n x n st 11 x 1 + 12 x 2 + + 1n x n = b 1 21 x 1 + 22 x 2 + + 2n x n = b 2 : m1 x 1 + m2 x 2 + + mn x n = b m x 1, x 2,, x n 0. where b 1, b 2,, b n 0. More compctly, let c=[c 1 c 2... c n ] nd A= : 11 21 m1 12 22 : m2......... 1n 2n : mn x1 x 2 x= : x n b1 b 2 nd x= : b m

Then the LP my be restted s Min (or Mx) cx, st Ax=b, x 0, where b 0. This finl system of constrint equtions of the LP in stndrd form is clled the AUGMENTED SYSTEM. Consider such n ugmented system of n vribles in m equtions, where m n. Note tht the ugmented system does not include the nonnegtivity conditions! Also note tht this system hs three possible sttes: (1) no solution, (2) exctly one solution, or (3) infinitely mny solutions. Fesible solution: Any solution to the ugmented system tht lso stisfies nonnegtivity is clled fesible solution E.g. Consider the following constrints nd the corresponding ugmented system: x 1 + 2x 2 120 x 1 + 2x 2 + S 1 = 120 x 1 + x 2 90 x 1 + x 2 + S 2 = 90 x 1 70 x 1 + S 3 = 70 x 2 50 x 2 + S 4 = 50

Note tht n=6 nd m=4. 1) An (infesible) solution is [x 1 =70, x 2 =50, S 1 =-50, S 2 =-30, S 3 =0, S 4 =0] 2) A fesible solution is [x 1 =20, x 2 =20, S 1 =60, S 2 =50, S 3 =50, S 4 =30] Bsic Solution: Suppose we fix (n-m) out of the n vribles t zero, nd try to solve the system of m equtions in the remining m vribles. If solution to this exists, then it is clled bsic solution. Solution (1) on the previous pge is n exmple of bsic solution. Bsic Fesible Solution: A bsic solution tht lso stisfies nonnegtivity is clled bsic fesible solution (BFS). An exmple of BFS is [x 1 =20, x 2 =50, S 1 =0, S 2 =20, S 3 =50, S 4 =0]. Note tht (1) is not BFS even though it is bsic solution! Augmented System One Solution No Solution Infinitely mny solutions If it is lso fesible ( 0) Originl LP is infesible then the originl LP hs Fesible Infesible exctly one fesible solution - it is lso optiml! Bsic Nonbsic Bsic Nonbsic

FACT: Ech BFS corresponds to n extreme point of the fesible region. In bsic solution, the (n-m) vribles tht re chosen to be fixed t zero re clled the nonbsic vribles nd the remining m vribles re clled bsic vribles. Note tht we cn choose (n-m) out of the n vribles in n n n! = = different wys. Thus we cn get n m m ( n m)! m! mximum of this mny different bsic solutions (some of which will lso be bsic fesible solutions) in prctice ll of these my not exist since in some cses, the resulting m x m system my not hve solution! n 6! In our exmple, n=6 nd m=4 so tht = = 15 n m (6 4)!4! The 15 combintions of nonbsic vribles re: 1) x 1, x 2 2) x 1, S 1 6) x 2, S 1 3) x 1, S 2 7) x 2, S 2 10) S 1, S 2 4) x 1, S 3 8) x 2, S 3 11) S 1, S 3 13) S 2, S 2 5) x 1, S 4 9) x 2, S 4 12) S 1, S 4 14) S 2, S 4 15) S 3, S 4 Try to locte ech of these on the grphicl representtion!

X 2 Mx 20X 1 +10X 2 S 3 =0, X 1 =70 st X 1 +2X 2 120 (Constr. 1, S 1 ) X 1 +X 2 90 (Constr. 2, S 2 ) 90 X 1 70 (Constr. 3, S 3 ) S 2 =0, X 1 +X 2 =90 X 2 50 (Constr. 4, S 4 ) X 1, X 2 0 50 II S 1 =S 4 =0 S 2 =S 4 =0 S 3 =S 4 =0 S 4 =0, X 2 =50 III 10 I S 1 =S 2 =0 X * 1 =70 IV S 1 =S 3 =0 X * 2 =20 Fesible Region Opt. Vl. =1600 V S 2 =S 3 =0 VI S 1 =0, X 1 +2X 2 =120 10 20 30 40 50 60 70 80 90 100 110 120 X 1 1., : BASIC SOLUTIONS (13 of them). Note tht there should be ( 4 6 )=15 of these, but only 13 exist becuse X 1 =70 is prllel to the X 2 -xis nd X 2 =50 is prllel to the X 1 -xis. 2. : BASIC FEASIBLE SOLUTION: (6 of the bove 13, numbered I, II, III, IV, V, VI) 3. : Contour of the objective function corresponding to vlue of 1600 4. : Set of ll Fesible solutions to the LP