Control Del Aire Progrmacion Lineal

Download as pdf or txt
Download as pdf or txt
You are on page 1of 8

Operations Research Fall 2012

Review 2: Controlling Air Pollution

Background: The three main types of pollutants in this airshed are


particulate matter (sulfur oxides, and hydrocarbons).
Objectives:
The new standards require that the company reduce its annual emission of
these pollutants by the amounts 60, 150 and 125 respectively.
The steelworks has two primary sources of pollution, namely, the blast
furnaces for making pig iron and the open-hearth furnaces for changing iron
into steel.
Table 3.13 shows how much emission (in millions of pounds per year) can
be eliminated from each type of furnace by fully using any abatement
method to its technological limit.

Zhou, Jian: [email protected] School of Management, Shanghai University 48 / 56


Operations Research Fall 2012

Review 2: Controlling Air Pollution


The NORI and LEETS CO.

Zhou, Jian: [email protected] School of Management, Shanghai University 49 / 56


Operations Research Fall 2012

Review 2: Controlling Air Pollution

How to solve the problem?

Step 1: List four basic elements

X Decision variables:
Fraction of the maximum feasible use of an abatement method (so xj
cannot exceed 1), xj (j = 1, 2, · · · , 6)
X Objective function:
Z –Minimum total annual cost
X Constraints:
Reduction requirements, fraction constraints and nonnegative
constraints
X Parameters:
······
Zhou, Jian: [email protected] School of Management, Shanghai University 50 / 56
Operations Research Fall 2012

Review 2: Controlling Air Pollution

Step 2: Give the mathematical model

minZ = 8x1 + 10x2 + 7x3 + 6x4 + 11x5 + 9x6







subject to:







12x1 + 9x2 + 25x3 + 20x4 + 17x5 + 13x6 > 60




 35x1 + 42x2 + 18x3 + 31x4 + 56x5 + 49x6 > 150





 37x1 + 53x2 + 28x3 + 24x4 + 29x5 + 20x6 > 125


xj 6 1, forj = 1, 2, · · · , 6 xj > 0, forj = 1, 2, · · · , 6

Zhou, Jian: [email protected] School of Management, Shanghai University 51 / 56


Operations Research Fall 2012
Step 4:Results by Matlab:
Step 3: Solve the model by software

c = [8, 10, 7, 6, 11, 9];


a = [−12, −9, −25, −20, −17, −13;
−35, −42, −18, −31, −56, −49;
−37, −53, −28, −24, −29, −20];
b = [−60, −150, −125];
lb = [0, 0, 0, 0, 0, 0];
ub = [1, 1, 1, 1, 1, 1];
Aeq = [ ];
Beq = [ ];
[x, fval ] = linprog (c, a, b, Aeq, Beq, lb, ub)

Zhou, Jian: [email protected] School of Management, Shanghai University 52 / 56


Operations Research Fall 2012

Solve the Problem by Excel

H3 = SUMPRODUCT (B3 : G 3 ∗ B11 : G 11)


H4 = SUMPRODUCT (B4 : G 4 ∗ B11 : G 11)
······

Zhou, Jian: [email protected] School of Management, Shanghai University 53 / 56


Operations Research Fall 2012

Solve the Problem by Lindo

Lindo Code:

min 8x1 + 10x2 + 7x3 + 6x4 + 11x5 + 9x6


subject to

12x1 + 9x2 + 25x3 + 20x4 + 17x5 + 13x6 <= 60


35x1 + 42x2 + 18x3 + 31x4 + 56x5 + 49x6 <= 150
37x1 + 53x2 + 28x3 + 24x4 + 29x5 + 20x6 <= 125

Zhou, Jian: [email protected] School of Management, Shanghai University 54 / 56


Operations Research Fall 2012

Results by Lindo

Zhou, Jian: [email protected] School of Management, Shanghai University 55 / 56

You might also like