MESM543
This presentation is the property of its rightful owner.
Sponsored Links
1 / 87

MESM543 PowerPoint PPT Presentation


  • 83 Views
  • Uploaded on
  • Presentation posted in: General

MESM543. Operations Research OPTIMIZATION. Operations Research (OR) is the field of how to form mathematical models of complex management decision problems and how to analyze the models to gain insight about possible solutions.

Download Presentation

MESM543

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.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.


- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -

Presentation Transcript


Mesm543

MESM543

Operations Research OPTIMIZATION

Operations Research (OR) is the field of how to form mathematical models of complex management decision problems and how to analyze the models to gain insight about possible solutions.


Mesm543

Applications grouped by type of organizational client

  • Business

  • Government and Non-Profit

  • Health Care

  • Military


Mesm543

Applications grouped by function

  • Planning, Strategic Decision-Making

  • Production

  • Distribution, Logistics, Transportation

  • Supply Chain Management

  • Marketing Engineering

  • Financial Engineering


Mesm543

Build Your Knowledge to increase your success in practice

  • Linear Programming

  • Non-linear Programming

  • Dynamic Programming

  • Markov Decision Processes

  • Multiple Criteria Decision Making

  • Queuing Models

  • General Simulation

Decisions


Mesm543

Production System

7


Mesm543

Production Systems

8


Mesm543

OR Process

Assessment

Real world problem

Real world solution

Abstraction

Interpretation

Analysis

Model

Model solution

9


Mesm543

Operations Research Techniques

  • Linear Programming

  • Dynamic Programming

  • Integer Programming

  • Nonlinear programming

  • Goal Programming

  • Network Programming

10

26.10.2014


Mesm543

The principal phases for implementing OR in practice

  • Definition of the problem

    • Decision alternative definition

    • Objective determination

    • Operational limitations specification

  • Construction of the model

    Entails translating the problem definition into mathematical relationships

  • Solution of the model

    It entails the use of well-defined optimization algoritms

  • Validation of the model

    Does the model provide a reasonable prediction of the system’s behavior?

  • Implementation of the solution

    Translation of the model’s results into operating instructions


Mesm543

Linear Programming (LP)

An important topic of Deterministic Operations Research

  • Agenda

  • Modeling problems

  • Examples of models and some classical problems

  • Graphical interpretation of LP

  • Solving LP by Simplex using MS Excel

  • Some theoretical ideas behind LP and Simplex


Mesm543

Example 1: Product Mix Problem

Fertilizer manufacturing company, 2 types of fertilizer

Type A: high phosphorus

Type B: low phosphorus


Mesm543

Product Mix Problem: Modeling

Step 1. The decision variables

Daily production ofType A: x tons

Type B: y tons

Step 2. The objective function (maximize profit)

z =15x + 10y


Mesm543

Product Mix Problem: Modeling..

Step 3. The constraints

Limited supply of raw materials per day:

Urea:2x + y ≤ 1500

Potash:x + y ≤ 1200

Rock Phosphate:x ≤ 500


Mesm543

Product Mix Problem: Complete model

Maximizez( x, y) = 15 x + 10y

subject to2x + y ≤ 1500

x + y ≤ 1200

x ≤ 500

x ≥ 0,

y ≥ 0

Interesting Aspects:

Linearity, Inequalities

Feasible solutions: (0, 0), (1, 1), …

Infeasible solutions: (600, 500), …


Mesm543

Background: Petroleum Refinery

Example 2. Blending Problem

Three types of petrol (minimum Octane rating: 85, 90, 95)

Four types of oils (Octane rating: 68, 86, 91, 99)

Blending oils  petrol, with proportional Octane rating

Objective: best product mix [how much of each petrol, oil to sell]


Mesm543

Example 2. Blending Problem, the data


Mesm543

68x11 + 86x21 + 91x31 + 99x41

Its Octane Rating:

≥ 95,

x11 + x21 + x31 + x41

Blending Problem: Modeling

Step 1. The decision variables

xij = barrels/day of oil i( i = 1, 2, 3, or 4)

to make petrol j (j = 1, 2, or 3)

Total premium petrol per day = x11 + x21 + x31 + x41

68x11 + 86x21 + 91x31 + 99x41 - 95(x11 + x21 + x31 + x41) ≥ 0.


Mesm543

premium

45.15(x11 + x21 + x31 + x41) +

42.95(x12 + x22 + x32 + x42) +

40.99(x13 + x23 + x33 + x43) +

36.85 (4000 – (x11 + x12 + x13)) +

36.85 (5050 – (x21 + x22 + x23)) +

38.95 (7100 –(x31 + x32 + x33)) +

38.95 (4300 – (x41 + x42 + x43))

super

regular

Oil 1

Oil 2

Oil 3

Oil 4

Blending Problem: Modeling..

Step 2.The objective function

Maximize profit  Maximize revenue


Mesm543

Blending Problem: Modeling...

Step 3.The constraints

(a) The OcR constraints:

68x11 + 86x21 + 91x31 + 99x41 - 95(x11 + x21 + x31 + x41) ≥ 0

68x12 + 86x22 + 91x32 + 99x42 - 90(x12 + x22 + x32 + x42) ≥ 0

68x13 + 86x23 + 91x33 + 99x43 - 85(x13 + x23 + x33 + x43) ≥ 0


Mesm543

Blending Problem: Modeling....

Step 3.The constraints..

(b) Can’t use more oil than we have:

x11 + x12 + x13 ≤ 4000

x21 + x22 + x23 ≤ 5050

x31 + x32 + x33 ≤ 7100

x41 + x42 + x43 ≤ 4300


Mesm543

Blending Problem: Modeling…..

Step 3.The constraints...

(c) The demand constraints:

x11 + x21 + x31 + x41 ≤ 10,000

x13 + x23 + x33 + x43 ≥ 15,000

(d) Allowed values of variables

xij ≥ 0 for i = 1, 2, 3, 4, and j = 1, 2, 3.


Mesm543

Blending Problem: complete model

Maximize:

45.15(x11 + x21 + x31 + x41) + 42.95(x12 + x22 + x32 + x42) + 40.99(x13 + x23 + x33 + x43) +

36.85(4000 – (x11 + x12 + x13)) + 36.85 (5050 – (x21 + x22 + x23)) +

38.95 (7100 –(x31 + x32 + x33)) + 38.95 (4300 – (x41 + x42 + x43))

Subject to:

68x11 + 86x21 + 91x31 + 99x41 - 95(x11 + x21 + x31 + x41) ≥ 0

68x12 + 86x22 + 91x32 + 99x42 - 90(x12 + x22 + x32 + x42) ≥ 0

68x13 + 86x23 + 91x33 + 99x43 - 85(x13 + x23 + x33 + x43) ≥ 0

x11 + x12 + x13 ≤ 4000

x21 + x22 + x23 ≤ 5050

x31 + x32 + x33 ≤ 7100

x41 + x42 + x43 ≤ 4300

x11 + x21 + x31 + x41 ≤ 10,000

x13 + x23 + x33 + x43 ≥ 15,000

xij ≥ 0 for I = 1, 2, 3, 4, and j = 1, 2, 3.

Octane

rating

Supply

Demand


Mesm543

Example 3: Transportation problem

Background:

Company has several factories (sinks), and

several suppliers (sources)

Objective: Minimize the cost of transportation


Mesm543

Example 3. Transportation problem, the data


Mesm543

Transportation problem: the model

Step 1. The decision variables

xij= amount of ore shipped from mine i to plant j per day.

Step 2:The objective function

Minimize the transportation costs:

Minimize: 11x11 + 8x12 + 2x13 + 7x21 + 5x22 + 4x23


Mesm543

Transportation problem: the model..

Step 3.The constraints

(a) Shipment from each mine less than daily production

x11 + x12 + x13 ≤ 800[capacity of mine 1]

x21 + x22 + x23 ≤ 300[capacity of mine 2]

(b) Demand of each plant must be met

x11 +x21≥ 400[demand at plant 1]

x12 + x22≥ 500[demand at plant 2]

x13+ x23≥ 200[demand at plant 3]

(c) Decision variables can’t be negative

xij ≥ 0, for all i= 1, 2, j = 1, 2, 3.


Mesm543

Transportation problem: historical note

Kantorovich in USSR in the 1930’s,

Koopmans in 1940’s

Dantzig in 1950’s  Simplex method

Kantorovich and Koopmans, Nobel prize (Economics) in 1975


Mesm543

Point in a 1D space:x = c

c

0

y

3

2x+3y = 9

2x+3y = 3

2

2x+3y = 6

2x+3y = 0

1

x

3

4.5

(0,0)

1.5

The Geometry of Linear Programs

Line in 2D:ax + by = c


Mesm543

z

1

Plane

x + y + z = 1

Plane:

x + y + z = 2

z

Plane:

x + y + z = 0

1

x

1

1

1

y

y

1

x

The Geometry of Linear Programs

Plane in 3D:ax + by + cz = d


Mesm543

y

3

2x+3y > 6

2

2x+3y < 6

1

x

4.5

3

1.5

(0,0)

2x+3y = 6

The Geometry of Linear Programs

Hyper-plane in n-Dimensions: a1x1 + a2x2 + … + anxn = c

??

2-D Half spaces:


Mesm543

y

y ≥0

1500

1000

500

x ≥0

x

500

Feasible set

1000

(0,0)

1500

x + y ≤ 1200

x ≤ 500

2x+y ≤ 1500

The Geometry of LP: Product Mix revisited

maxz( x, y) = 15 x + 10y

ST2x + y ≤ 1500

x + y ≤ 1200

x ≤ 500

x ≥ 0,

y ≥ 0


Mesm543

y

y ≥0

15x + 10y = 5000

1500

15x + 10y = 13,500

1000

15x + 10y = 0

(300, 900)

500

x ≥0

x

500

Feasible set

1000

(0,0)

1500

x + y ≤ 1200

x ≤ 500

2x+y ≤ 1500

The Geometry of LP: Product Mix revisited

maxz( x, y) = 15 x + 10y

ST2x + y ≤ 1500

x + y ≤ 1200

x ≤ 500

x ≥ 0,

y ≥ 0

Try point: x = 0, y = 0:


Mesm543

Summary

1. LP formulations are very common in modern industry

2. Beautiful connection between Algebra and Geometry

3. Geometry not useful for > 3 variables

4. Practical problems: 1000’s of variables (see next slide)

5. Need Algebraic method !


Mesm543

Dynamic Programming – Recursive formula (forward, Backward)

A state k is the information that is needed at any stage to make an optimal decision.

A problem can be decomposed into stages, where each stage j represents an operation.

Define fj(xj) as fj(xj) = maximum return for stage i given state yj

Forward Recursive:


Mesm543

Forward vs. Backward Recursion

Backward recursion formula:


Mesm543

Budgeting Problem –FirstExample1.m

A corporation has $5 million to allocate to its three plants for possible expansion Each plant has submitted a number of proposals on how it intends to spend the money.

Each proposal gives the cost of the expansion (c) and the total revenue expected (r). The following table gives the proposals generated:

Problem formulation

Define state x[6]=(0 1 2 3 4 5); available resources

Define Stage j[3]=(1 2 3); The operation or processing of the problem takes place

Define Policy or Proposal of problem m[4]=(1 2 3 4)


Mesm543

%Dynamic Programming 1st example, CMU

% F. Nagi 5/2/12

clear

clc

Alloc = 6; %Allocation 7 million +1 for index '0'

r=[0 5 6 0;

0 8 9 12;

0 4 0 0]; %revenue of Proposal

c=[0 1 2 0;

0 2 3 4;

0 1 0 0]; %cost of proposal

x =[0 ;

1 ;

2 ;

3 ;

4 ;

5 ];

f2(1:Alloc,1:Alloc)=0; %initialze Answer function f2

x=x';

for j=1 %state

for m=1:4 %proposal

for k=1:Alloc %state

if c(j,m)<=x(k)

f1(m,k)=r(j,m);

end

end

end

end

f1=max(f1)

for j =2:3 %stage

for m=1:4 %proposal

for k=1:Alloc %states

if c(j,m)<=x(k)

f1(j,k)=max(r(j,m)+f1(j-1,x(k) c(j,m)+1), f2(j,k))

f2(j,k)=f1(j,k);

end

end

end

end

f2=max(f1)


Mesm543

WorkStation Production Problem

Ans, f1 =

0 5 6 6 6 6

0 5 8 13 14 17

0 5 9 13 17 18

Ans, f2 =

0 5 9 13 17 18


Mesm543

Problem formulation

Resource Alloc = 5; %Allocation 4 cycle time +1 for index '0'

r=[0 31 62 0 0;

0 47 0 0 0;

0 14 28 42 56]; %revenue of Products

c=[0 2 4 0 0;

0 3 0 0 0;

0 1 2 3 4]; %cycle time of products

x =[0 ;

1 ;

2 ;

3 ;

4 ]; % cycle time


Mesm543

The cycle of a Genetic Algorithms is presented below

Each cycle in Genetic Algorithms produces a new generation of possible solutions for a given problem


Mesm543

Genetic Algorithm Steps

  • The elements of the population are encoded into bit-strings, called chromosomes.

  • The performance of the strings, often called fitness, is then evaluated with the help of some functions f(x) , representing the constraints of the problem.

  • The crossover operation that recombines the bits (genes) of each two selected strings (chromosomes) . The crossover helps to span over the solution space


Mesm543

Genetic Algorithm Steps

In Mutation the bits at one or more randomly selected positions of the chromosomes are altered. The mutation process helps to overcome trapping at local maxima

Mutation of a chromosome at

the 5th bit position.

Example: The Genetic Algorithms cycle is illustrated in this example for maximizing a function f(x) = x2in the interval 0 = x = 31.

In this example the fitness function is f (x) itself. Starts with 4 initial strings. The fitness value of the strings and the percentage fitness of the total are estimated in Table A.


Mesm543

Example

Table A

  • Since fitness of the second string is large, we select 2 copies of the second string and one each for the first and fourth string in the mating pool. The selection of the partners in the mating pool is also done randomly

  • Here in table B, we selected partner of string 1 to be the 2-nd string and partner of 4-th string to be the 2nd string. The crossover points for the first-second and second-fourth strings have been selected after 0-th and 2-nd bit positions respectively in table B.

Table B


Mesm543

Example

The second generation of the population without mutation in the first generation is presented in table C.

Table C


Mesm543

gaex.m % Example in >>gatool

0-4, five bits

25 =32 (0-31)

Fitness f(x)

x -variable


Mesm543

GA- Problems

  • Use Gatool and maximize the quadratic equation

  • f(x) = x2 +4x-1 within the range

  • −5 ≤x≤0.

2. Use Gatool and maximize the function

f(x1, x2, x3)= −5 sin(x1) sin(x2) sin(x3) + (- sin(5x1) sin(5x2)sin(x3))

where 0 <= xi <= pi, for 1 <= i <= 3.

3. Create a “gatool” to minimize the function

f(x) = cos(2x) within the range 0≤x≤3.14


  • Login