Sufficiency Conditions for Equality Constrained Nonlinear Programming: Proof & Sensitivity, Slides of Computer Science

An in-depth analysis of sufficiency conditions for equality constrained nonlinear programming problems. It covers the second order sufficiency conditions, convexification using augmented lagrangians, and proof of the sufficiency conditions. Additionally, the document discusses sensitivity analysis to understand how the minimum changes when constraints are modified.

Typology: Slides

2012/2013

Uploaded on 03/27/2013

ekana
ekana 🇮🇳

4

(44)

370 documents

1 / 8

Toggle sidebar

This page cannot be seen from the preview

Don't miss anything!

bg1
NONLINEAR PROGRAMMING
LECTURE 12: SUFFICIENCY CONDITIONS
LECTURE OUTLINE
EqualityConstrainedProblems/SufciencyCon-
ditions
Convexication Using Augmented Lagrangians
Proof of the Sufciency Conditions
Sensitivity
Equality constrained problem
minimize f (x)
subject to hi(x)=0, i =1,...,m.
where f : n →, hi : n →, are continuously
differentiable. To obtain sufciency conditions, as-
sume that f and hi are twice continuously differen-
tiable.
Docsity.com
pf3
pf4
pf5
pf8

Partial preview of the text

Download Sufficiency Conditions for Equality Constrained Nonlinear Programming: Proof & Sensitivity and more Slides Computer Science in PDF only on Docsity!

NONLINEAR PROGRAMMING

LECTURE 12: SUFFICIENCY CONDITIONS

LECTURE OUTLINE

• Equality Constrained Problems/Sufficiency Con-

ditions

• Convexification Using Augmented Lagrangians

• Proof of the Sufficiency Conditions

• Sensitivity

Equality constrained problem

minimize f (x)

subject to h

i

(x) = 0, i = 1,... , m.

where f : 

n 

→ , h

i

n 

→ , are continuously

differentiable. To obtain sufficiency conditions, as-

sume that f and h

i

are twice continuously differen-

tiable.

SUFFICIENCY CONDITIONS

Second Order Sufficiency Conditions: Let x

∗ ∈ 

n

and λ

∗ ∈ 

m

satisfy

∗ ∗ ∗ ∗

∇ x

L(x , λ ) = 0, ∇ λ

L(x , λ ) = 0,

y

2

xx

L(x

, λ

)y > 0 , ∀ y = 0 with ∇h(x

)

y = 0.

Then x

is a strict local minimum.

Example: Minimize −(x

1

x 2

  • x 2

x 3

  • x 1

x 3

) subject to

∗ ∗

x 1

  • x 2

  • x 3

= 3. We have that x

∗ = x 2

= x 3

= 1 and

1

λ

= 2 satisfy the 1st order conditions. Also

∗ ∗

2

xx

L(x , λ ) = − 1 0 − 1.

We have for all y = 0 with ∇h(x

∗ )

y = 0 or y

1

  • y 2

y 3

∗ ∗

y

2

xx

L(x , λ )y = −y 1

(y 2

  • y 3

) − y 2

(y 1

  • y 3

) − y 3

(y 1

  • y 2

2 2 2

= y 1

  • y 2

  • y 3

Hence, x

is a strict local minimum.

PROOF OF SUFFICIENCY CONDITIONS

Consider the augmented Lagrangian function

c

Lc(x, λ) = f (x) + λ

h(x) + ‖h(x)‖

2

,

where c is a scalar. We have

∇xLc(x, λ) = ∇xL(x, λ),

xx

Lc(x, λ) = ∇

2 ˜ ∇

2

xx

L(x, λ) + c∇h(x)∇h(x)

where

λ = λ + ch(x). If (x

∗ , λ

) satisfy the suff. con-

ditions, we have using the lemma,

, λ

) = 0, ∇

2

xx

Lc(x

, λ

∇xLc(x ) > 0 ,

for suff. large c. Hence for some γ > 0 , � > 0 ,

) +

γ

‖x − x

) ≥ L c

(x

, λ

∗ ∗

2

L , if ‖x − x ‖ < �.

c

(x, λ

Since L

c

(x, λ

) = f (x) when h(x) = 0,

) +

γ

‖x − x

∗ ∗

2

f (x) ≥ f (x , if h(x) = 0, ‖x − x ‖ < �.

SENSITIVITY - GRAPHICAL DERIVATION

∇f(x

)

x

  • ∆x

x

∆x

a a'x = b + ∆b

a'x = b

Sensitivity theorem for the problem min a

� x=b

f (x). If b is

changed to b+∆b, the minimum x will change to x

∗ +∆x.

Since b + ∆b = a

′ (x

  • ∆x) = a

′ x + a

′ ∆x = b + a

′ ∆x, we

have a

′ ∆x = ∆b. Using the condition ∇f (x

∗ ) = −λ

∗ a,

∗ ∗ ∗ ∆cost = f (x + ∆x) − f (x ) = ∇f (x )

′ ∆x + o(‖∆x‖)

= −λ a

∆x + o(‖∆x‖)

Thus ∆cost = −λ

∗ ∆b + o(‖∆x‖), so up to first order

∆cost ∗

λ = −.

∆b

For multiple constraints a

i

x = b i

, i = 1,... , n, we have

m

∆cost = − λ i

∆bi + o(‖∆x‖).

i=

EXAMPLE

p(u)

-1 0

u slope ∇p(0) = - λ

= -

Illustration of the primal function p(u) = f x(u)

for the two-dimensional problem

2 2

minimize f (x) =

1 x 1

− x 2

− x 2 2

subject to h(x) = x 2

Here,

p(u) = min f (x) = −

1

2

u

2

− u

h(x)=u

and λ

∗ = −∇p(0) = 1, consistently with the sensitivity

theorem.

• Need for regularity of x

: Change constraint to

h(x) = x

2

2

= 0. Then p(u) = −u/ 2 − u for u ≥ 0 and

is undefined for u < 0.

PROOF OUTLINE OF SENSITIVITY THEOREM

Apply implicit function theorem to the system

∇f (x) + ∇h(x)λ = 0, h(x) = u.

For u = 0 the system has the solution (x

∗ , λ

), and

the corresponding (n + m) × (n + m) Jacobian

m

i=

λ i

∗ ∇

2 h i

(x

∗ ) ∇h(x

∗ ∇ )

2 f (x

∗ ) +

J =

∇h(x

∗ )

′ 0

is shown nonsingular using the sufficiency con-

ditions. Hence, for all u in some open sphere S

centered at u = 0, there exist x(u) and λ(u) such

that x(0) = x

, the functions x(·) and λ(·)

are continuously differentiable, and

∇f x(u) + ∇h x(u) λ(u) = 0, h x(u) = u.

For u close to u = 0, using the sufficiency condi-

tions, x(u) and λ(u) are a local minimum-Lagrange

multiplier pair for the problem min

h(x)=u

f (x).

To derive ∇p(u), differentiate h x(u) = u, to

obtain I = ∇x(u)∇h x(u) , and combine with the re-

lations ∇x(u)∇f x(u) + ∇x(u)∇h x(u) λ(u) = 0 and

∇p(u) = ∇ f x(u) = ∇x(u)∇f x(u). u