Post on 28-Dec-2015
1
Ch
apte
r 8
Chapter 8
Nonlinear Programming with Constraints
2
Ch
apte
r 8
3
Ch
apte
r 8
4
Ch
apte
r 8
Methods for Solving NLP Problems
5
Ch
apte
r 8
*)(**)( xhxf (a) Where 414.1/1* is called the Lagrange multiplier for the constraint h = 0
)21,21(* x ; see Fig. E 8.1a
6
Ch
apte
r 8
7
Ch
apte
r 8
8
Ch
apte
r 8
9
Ch
apte
r 8
10
Ch
apte
r 8
11
Ch
apte
r 8
12
Ch
apte
r 8
Note that there are n + m equations in the n + m unknowns x and λ
13
Ch
apte
r 8
14
Ch
apte
r 8
15
Ch
apte
r 8
16
Ch
apte
r 8
17
1 2Minimize : ( )f x xx
2 21 2Subject to : ( ) 25g x x x
By the Lagrange multiplier method.
Solution: The Lagrange function is
2 21 2 1 2( , ) ( 25)L u x x x x x u
The necessary conditions for a stationary point are
2 11
2 0L
x uxx
1 22
2 0L
x uxx
2 21 2 25
Lx x
u
2 21 2(25 ) 0u x x
Ch
apte
r 8
18
Ch
apte
r 8
19
Ch
apte
r 8
20
Ch
apte
r 8
21
Ch
apte
r 8
Penalty functions for handling equality constraints
22
Ch
apte
r 8
23
Ch
apte
r 8
for handling inequality constraints
Note g must be >0 ; r 0
24
Ch
apte
r 8
The logarithmic barrier function formulation for m constraints is
25
Ch
apte
r 8
26
Ch
apte
r 8
27
Ch
apte
r 8
28
Ch
apte
r 8
29
Ch
apte
r 8
Use xc = 2 yc = 2 for linearization
(step bounds)
30
Ch
apte
r 8
31
Ch
apte
r 8
32
Ch
apte
r 8
33
Ch
apte
r 8
Quadratic Programming (QP)
34
Ch
apte
r 8
8.3 QUADRATIC PROGRAMMING
35
Use of Quadratic Programming to Design Multivariable Controllers
(Model Predictive Control)
• Targets (set points) selected by real-time optimization software based on current operating and economic conditions
• Minimize square of deviations between predicted future outputs and specific reference trajectory to new targets using QP
• Framework handles multiple input, multiple output (MIMO) control problems with constraints on manipulated and controlled variables. Dynamics obtained from transfer function model.
Ch
apte
r 8
36
Successive Quadratic Programming
• Considered by some to be the best general nonlinear programming algorithm
• Repetitively approximates nonlinear objective function with quadratic function and nonlinear constraints with linear constraints
• Uses line search rather than QP step for each iteration• Inequality constrained Quadratic Programming (IQP)
keeps all inequality constraints• Equality constrained Quadratic Programming (EQP) only
keeps equality constraints by utilizing and active set strategy
• SQP is an Infeasible Path method
Ch
apte
r 8
37
Ch
apte
r 8
38
Ch
apte
r 8
solve for ,x
39
Ch
apte
r 8
Generalized Reduced Gradient (GRG)
40
Ch
apte
r 8
41
Ch
apte
r 8
42
Ch
apte
r 8
43
Ch
apte
r 8
44
Ch
apte
r 8
45
Ch
apte
r 8
46
Ch
apte
r 8
47
Ch
apte
r 8
48
Ch
apte
r 8
49
Ch
apte
r 8
50
Ch
apte
r 8
51
• sequential simplex• conjugate gradient• Newton’s method• Quasi-Newton
Ch
apte
r 8
52
Ch
apte
r 8
53
Ch
apte
r 8
54
Ch
apte
r 8
55
Ch
apte
r 8
56
Ch
apte
r 8
57
Ch
apte
r 8
58
Ch
apte
r 8
59
Ch
apte
r 8
60
Ch
apte
r 8
61
Ch
apte
r 8
62
Ch
apte
r 1
we have to supplement blast furnace gas with fuel oil, but we want to minimize the purchase of fuel oil.
Ch
apte
r 8
63
Ch
apte
r 1
define: X1 = amount of fuel oil used in generator 1X2 = amount of fuel oil used in generator 2X3 = amount of BFG used in generator 1X4 = amount of BFG used in generator 2P1 = mw output of generator 1P2 = mw output of generator 2
range of operation of generator 1 and generator 218 ≤ P1 ≤ 3014 ≤ P2 ≤ 25
Fuel effects in the generators are additive (can operateon either BFG or fuel oil)
Ch
apte
r 8
64
Ch
apte
r 1
10 units of BFG are available (on the average):
1 unit BFG = Btu equivalent of 1 ton/hr. fuel oil.
We need 50 mw power at all times.
5021 PP
Experimental data needed?
Ch
apte
r 8
65
Ch
apte
r 1
Mathematical Statement
1 2
fuel oil to gen 1 fuel oil to gen 2
Min f x x
a. operating ranges 18 ≤ P1 ≤ 30& requirements 14 ≤ P2 ≤ 25
P1 + P2 = 50
b. availability of x3 + x4
blast furnace gas
c. operating P11 (x1) fuel oilcharacteristics P12 (x3) BFG
P21 (x2) fuel oil P22 (x4)
BFG
gen 1
gen 2
2221212111 PPPPPP
Ch
apte
r 8
66
Ch
apte
r 1
BFG addmay
wesince sufficientnot isequation this
00145.15186.4609.1 21111 XXP
22212
422
221
12111
312
BFG )(
oil fuel )(
1generator in reqt.)(
PPP
XP
XP
PPP
XBFGP
fcn of burners, heat transfer characteristics (convex functions)
Ch
apte
r 8
67
Ch
apte
r 1
Solution
NLP 4 ineq. const.piece-wise LP 6 eq. const.
1
2
3.05
30
20
optf
P
P
No fuel oil is used in generator 1.In generator 2, fuel oil provides 58%of the power (rest is BFG).
heat transfer characteristics may change, or BFG mayvary w.r.t. time (on-line solution)
Ch
apte
r 8