Operations Research/Linear Programming
21/09/2020This article is copied from wikibooks for educational purposes, see the original article
Linear Programming (LP) is a mathematical modelling technique useful for allocation of limited resources such as material, machines etc to several competing activities such as projects, services etc. A typical linear programming problem consists of a linear objective function which is to be maximized or minimized subject to a finite number of linear constraints. (By a linear function we mean a function of the form where are all variables.)
The founders of LP are George B. Dantzig, who published the simplex method in 1947, John von Neumann, who developed the theory of the duality in the same year, and Leonid Kantorovich, a Russian mathematician who used similar techniques in economics before Dantzig and won the Nobel prize in 1975 in economics. The linear programming problem was first shown to be solvable in polynomial time by Leonid Khachiyan in 1979, but a larger major theoretical and practical breakthrough in the field came in 1984 when Narendra Karmarkar introduced a new interior point method for solving linear programming problems.
The Two Variable LP model
Two variable LP models don't usually exist in practice but their study can provide us with an example as to how a general model can be handled.
We will explain the model by an example.
Consider a chemical company which produces two salts: and . Suppose an acid and a base are the only two chemicals required by the company to produce these two salts. Further suppose that by past experience the owner of the company has obtained the following data:
To make unit of the salt , units of the acid and unit of the base is required. To make unit of the salt , units of the acid and units of the base are required. At most units of the acid and units of the base are available daily. unit of the salt gets him a profit of (in whatever currency) and unit of the salt Y gets him a profit of . In addition due to a market survey he knows that the daily demand for the salt cannot exceed that for by more than unit, and that the maximum daily demand for the salt is units.
The company owner wants the best product mix. That is to say, he wants to know the amount of and he should make daily so that he gets the highest profit.
To formulate the LP model for this problem we first need to identify the decision variables. These are the variables which will represent the entities about which we have to actually make a decision about. In our problem it is clear that the entities are the salts and and so the variables should represent the amount they should be each made. So let,
= Units produced daily of salt
= Units produced daily of salt
The next step is to decide what should be the objective function for our model. It is clear from the name itself that the objective function represents our objective i.e. maximizing our total profit. Since the total profit is usually the sum of the profits obtained from and it should be equal to . This is assuming that if unit of (or ) gives a profit of (or ) then (or ) units would give a profit of (or ). Hence if represents the total daily profit the objective is:
Maximize =
Now we consider the constraints that restrict the usage of the acids and the salts. We know that the total amount of acid which we use in the day can't exceed . Since and are reasonably the amounts of the acid used in making the salt and we may conclude that our constraint is:
Similarly, for the base: Now the market dictates that the excess of the daily production of Y over X which is can't exceed which means that:
Also the maximum daily demand for the salt is units so that:
An implicit assumption is that the variables x{\mathclap{1}} and x{\mathclap{2}} can't assume negative values (Why?). The nonnegativity restrictions, account for this requirement.
The complete model can be now stated as:
Maximize =
subject to,
Any values of and that satisfy all the five constraints constitute a feasible solution. Otherwise the solution is infeasible. For example, the solution and is feasible because on substitution in the constraints none of the inequalities are violated. On the other hand the solution, the solution and is infeasible because it does not satisfy constraint (1), namely which is larger then the right hand side ().
Our goal is to find the optimum solution, i.e. one which maximizes the objective function besides being feasible. The procedure to do that will be discussed in the succeeding sections.