Ly thuyêt Qui hoach tuyên tinh
CH NG II: C S QUI HO CH L IƯƠ Ơ
§ 1i toán qui ho ch l i c b n ơ
1.1 Trình bày bài toán qui ho ch l i c b n ơ
Cho Γ⊆ Rn, l i, f hàm l i xác đ nh trên Γ, gi, i =1,2,…,m, m hàm lõm
cũng xác đ nh trên Γ, bi , i = 1,2,…,m, m s th c. D th y r ng, n u t p ế
h p
X = {xRn/ gi(x) bi, i = 1,2,…,m; xΓ} (1.1)
không r ng thì nó là t p h p l i, đóng trong R n (Bài t p 42).
Đ cho g n, hi u g(x) = (g 1(x), g2(x), …, gm(x)) vect m hàmơ
s , b = (b1, b2, …, bm) Rm. Khi đó t p t p X đ c vi t d i d ng ượ ế ướ
X = { xRn / g(x) b, xΓ} (1.2)
trong đó hi u g(x) b gi(x) bi, i = 1,2, …, m. Bài toán c c tr sau
đây đ c g i là bài toán qui hoach lôi c b n.ượ ơ
Tìm x0
X, đ cho
0
x X
f (x ) f (x)
min
=
(1.3)
thuy t các ph ng pháp gi i bài toán (1.3) đ c g i ế ươ ượ
thuy t qui ho ch l i. Trong ch ng này ch trình bày ph n thuy t làmế ươ ế
c s cho vi c xây d ng các thu t toán. Các ph ng pháp đ gi i bài toánơ ươ
(1.3) s đ c trình bày trong chuyên đ Qui ho ch phi tuy n h c kỳ 7. ượ ế
1.2 Đi u ki n chính qui
Đn 1.1: T p h p X xác đ nh theo (1.1) ho c (1.2) đ c g i th a ượ mãn
đi u ki n cính qui , n uế
i, 1 i m, xi: gi(xi) > bi(1.4)
Ngoài ra, t p X cũng đ c coi tho mãn ượ đi u ki n chính qui
Slater, n u ế
x0: gi(x0) > bi i, 1 i m (1.5)
D th y r ng hai đi u ki n chính qui (1.4) (1.5) t ng đ ng ươ ươ
nhau; t c n u X th a mãn đi u ki n chính qui (1.4) thì cũng th a đi u ế
ki n chính qui (1.5) và ng c l i (Bài t p 43). ượ
Lê Văn Phi, Khoa Toan Thông kê, Tr ng Đai hoc Kinh tê Tp. Hô Chi Minhươ 43
Ly thuyêt Qui hoach tuyên tinh
1.3 Hàm Lagrange
ng v i m i x ∈Γ, xét vect m chi u h(x) = b – g(x), ơ
t c là
hi(x) = b – gi(x), i = 1,2,…,m
Trên Γ× {yRm/ y 0} ta đ nh nghĩa hàm s
L(x,y) = f(x) + y, h(x) = f(x) +
m
i i i
i 1
y [b g (x)]
=
(1.6)
Hàm L(x,y) đ c đ nh nghĩa theo (1.6) đ c g i là ượ ượ hàm Lagrange ng v i
bài toán qui hoach lôi (1.3).
§ 2c đi u ki n t i u c a bài toán ư qui ho
2.1 Đi m yên ng a
Cho S Rn, T Rm và ϕ(s,t) là hàm s xác đ nh trên S ×T.
Đn 2.1: Đi m (s*, t*)
S
×
T đ c g i là đi m yên ng a (Hình ) c a hàmượ
ϕ
(s, t) trên mi n S
×
T, n uế
sS, tT : ϕ(s*, t) ϕ(s*, t*) ϕ(s, t*)
(2.1)
Hinh 2.1 Hinh 2.2
hình 2.1, đi m (s*, t*) = (0; 1) là đi m yên ng a c a m ϕ(s,t) = s + t. Vì,
n u gi nguyên t* = 1 và thay đ i s (s ch có th tăng t 0 lên 1) thì giá tr m s sế
tăng t 1 đ n 2. Trong khi đó, n u gi nguyên s* = 0 và thay đ i t (t ch có th gi m ế ế
t 1 đ n 0) thì giá tr hàm s gi m t 1 đ n 0. Trong khi đi m (1;0) không ph i ế ế
đi m yên ng a. ví d này cũng cho th y r ng, đi m yên ng a không ph i là c c
Lê Văn Phi, Khoa Toan Thông kê, Tr ng Đai hoc Kinh tê Tp. Hô Chi Minhươ 44
ϕ
ϕ(s, t*)
ϕ(s*, t*)
ϕ(s*, t)
t*
s* t
(s*, t*)
S
S = T = [0;1], ϕ(s, t) = s + t
ϕ
2
(1;1)
(0;1)
(1;0)
1 t
1
s
Ly thuyêt Qui hoach tuyên tinh
tr c a hàm Lagrange. (Giá tr c c ti u c a ϕ(s,t) = s + t là 0, t i (0;0), giá tr c c đ i
là 2 t i đi m (1;1).
D dàng ch ng minh đ c quan h sau đây: ượ
s S s S
t T t T
(s, t)
sup sup
inf inf
φ φ
(2.2)
Đlí 2.1: Đ cho hàm
ϕ
(s,t) đi m yên ng a trên mi n S
×
T thì đi u ki n
c n đ t n t i các minmax
s St T
(s, t)
maxinf
φ
s S t T
(s, t)
sup
min
φ
;
đ ng th i chúng b ng nhau .
Ch ng minh:
a) Đi u ki n c n: Gi s (s*,t*) đi m yên ng a c a hàm ϕ(s,t) trên
mi n S×T. Khi đó theo (2.1)
t T
(s*, t) (s*, t*)
sup
φ φ
và do đó
s S t T
(s, t) (s*, t*)
sup
inf
φ φ
T ng t ,ươ
s S
t T
(s*, t*) (s, t)
supinf
φ φ
Suy ra b t đ ng th c kép
s S s S s S
t T t T t T
(s, t) (s, t*) (s*, t*) (s, t*) (s, t)
sup sup sup
inf inf inf
φ φ φ φ φ
T (2.2) suy ra các d u b t đ ng th c bên trong ph i tr thành đ ng th c;
t c là
s S s S
t T
(s, t) (s, t*) (s*, t*)
supinf inf
φ = φ = φ
Hay
s St T
(s, t) (s*, t*)
maxinf
φ = φ
s S t T
(s, t*) (s*, t*)
sup
min
φ = φ
b) Đi u ki n đ : Gi s các giá tr minmax trong đ nh t n t i b ng
nhau, t c là s*S, t*T, sao cho
s S s St T
(s, t) (s, t*)
maxinf inf
φ = φ
s S t T t T
(s, t*) (s*, t)
sup sup
min
φ = φ
Suy ra
s S t T
(s, t*) (s*, t*) (s*, t)
sup
inf
φ φ φ
. Theo gi thi t các minmax b ng ế
nhau, nên các b t đ ng th c ph i tr thành đ ng th c. T c là
s S t T
(s, t*) (s*, t*) (s*, t)
sup
inf
φ = φ = φ
T đây suy ra sS, tT: ϕ(s*,t) ϕ(s*,t*) ϕ(s,t*). T c c p (s*,t*)
là đi m yên ng a c hàm ϕ(s,t) trên mi n S×T. ª1
1Trong ch ng minh trên các bi n s, t bi n thiên đ c l p, nên n u (s*,t*) (s’,t’) nh ng ế ế ế
đi m yên ng a c a ϕ(s,t) thì các c p (s*,t’) và (s’,t*) cũng là đi m yên ng a.
Lê Văn Phi, Khoa Toan Thông kê, Tr ng Đai hoc Kinh tê Tp. Hô Chi Minhươ 45
Ly thuyêt Qui hoach tuyên tinh
2.2 Đi m yên ng a c a hàm Lagrange L(x,y) đi u ki n t i
u ư
Theo ph n 1.3, hàm Lagrange c a bài toán qui hoach lôi c b n (1.3) ơ
xác đ nh trên mi n Γ× {yRm/ y 0} có d ng
L(x,y) = f(x) + y, h(x) = f(x) +
m
i i i
i 1
y [b g (x)]
=
(1.6)
Theo Đn 2.1, thì c p (x*, y*) là đi m n ng a c a L( x,y) trên mi n
Γ× {yRm/ y 0} khi
x∈Γ, y 0: L(x*,y) L(x*,y*) L(x,y*) (2.3)
Do L(x,y) l i theo x tuy n tính theo ếy, đ ng th i các t p Γ {yRm/ y
0} nh ng t p đóng, nên Đlí 2.2 tr thành: Đ cho hàm Lagrange
L(x,y) có đi m yên ng a trên mi n Γ× {yRm/ y 0} thì đi u ki n c n
đ là các minmax
xy 0
L(x, y)
maxmin
∈Γ
x y 0
L(x, y)
max
min
∈Γ
(2.4)
t n t i và b ng nhau.
Quan h gi a đi m yên ng a c a hàm L( x,y) l i gi i c a bài
toán qui hoach l i c b n (1.3) đ c cho b i đ nh lý sau đây. ơ ượ
Đlí 2.2: N u (ếx*, y*) đi m yên ng a c a hàm L( x,y) trên mi n
Γ×
{y
Rm, y
0} thì x* là l i gi i c a bài toán (1.3).
Ch ng minh: Gi s ( x*, y*) đi m yên ng a c a hàm L( x,y) trên mi n
Γ× {yRm, y 0}. T c là x∈Γ, y 0 :
f(x*) + y, h(x*) f(x*) + y*, h(x*) f(x) + y*, h(x) (2.5)
T đây suy ra:
a) y 0 : y, h(x*) y*, h(x*) (2.6)
D th y r ng đi u ày ch xãy ra khi h( x*) 0. T c b g(x*) 0; hay
g(x*) b. K t h p v i gi thi t ế ế x*∈Γ, suy ra x*X.
b) Do (2.6) đúng v i m i y 0 nên cũng đúng khi y = 0. T (2.6) suy ra
y*, h(x*) 0. M t khác, do h(x*) 0, nên ta cũng b t đ ng th c
ng c l i. ượ Suy ra y*, h(x*) = 0.
c) T đây và t (2.5) suy ra x∈Γ : f(x*) f(x) + y*, h(x). Đ c bi t khi
xX, thì h(x) 0, nên y*, h(x) 0. Do đó f(x*) f(x), xX. T c
x* là l i gi i c a bài toán (1.3).ª
Lê Văn Phi, Khoa Toan Thông kê, Tr ng Đai hoc Kinh tê Tp. Hô Chi Minhươ 46
Ly thuyêt Qui hoach tuyên tinh
Đ nh 2.2 cho bi t đi u ki n đ đ m t vect x* l i gi i c a bài toán ế ơ
QH l i c b n (1.3). T đ nh này suy ra, đ gi i bài toán (1.3) ch c n ơ
tìm m t đi m yên ng a c a hàm Lagrange t ng ng L(x,y). ươ
2.3 Đ nh lý Kuhn- Tucker
Đlý 2.3: Gi s r ng t p X th a mãn đi u ki n chính qui. Khi đó đ cho
đi m x*
∈Γ
l i gi i c a bài toán QH l i c b n (1.3) thì đi u ơ
ki n c n đ t n t i đi m y*
0 đ cho c p ( x*, y*)
đi m yên ng a c a hàm Lagrange L( x,y) trên mi n
Γ×
{y
Rm, y
0}.
Ch ng minh: Đi u ki n đ suy t đ nh 2.2 không c n đi u ki n
chính qui. Ta ch ng minh đi u ki n c n . Gi s x* là l i gi i c a bài toán
QH l i c b n (1.3). ơ
1) Xét tâp P = {(z0, z) R (m+1)/ z0 f(x*), z 0}. V i xΓ, đ t
S(x) = {(z0, z) R (m+1)/ z0 f(x), z b – g(x)}
Và đ t S =
x
S(x)
∈Γ
U
. Rõ ràng P, S đ u là t p l i trong R (m+1).
2) Đ t P0 = {(z0, z) R (m+1)/ z0 < f(x*), z < 0}. Khi y P0 = intP do đó P0
cũng t p l i. Gi s P 0 S . T c th tìm th y (z’ 0, z’) P0
S v i z’0 < f(x*), z < 0 z’0 f(x’), z b g(x’) ng v i ít nh t m t
x∈Γ. Nh v y, ư x X f(x’) < f(x*). Đi u này trái v i gi thi t r ng ế x*
là l i gi i c a (1.3).
3) Do P, S l i intP S = , nên theo đ nh tách 3 ch ng I, s t i t i ươ
m t siêu ph ng H α tách P v i S. T c là t n t i (u 0,u)R (m+1), (u0,u) 0,
đ cho
(z0, z)S (w0, w) P: u0z0 + u, z u0w0 + u, w ,
Hay u0 (z0 - w0) + u, z - w 0. (2.7)
Do (z0, z)S và (w0, w) P0 thì z0 f(x) f(x*) > w0, nên z0–w0 > 0. T
đây d dàng suy ra r ng, đ (2.7) đúng (z0, z)S và (w0, w) P0 thì u0
0. T ng t , cũng suy ra r ng uươ i 0, i = 1,2,…,m (bài t p…).
4) Ch n xΓ tùy ý, z0 = f(x), z = b – g(x) w0 = f(x*), w = 0. T (2.7) suy
ra:
xΓ: u0 f(x) + u, b – g(x) u0 f(x*).
(2.8)
N u uế0 = 0, ta có
xΓ: u, b g(x) 0. (2.9)
Do u 0, b – g(x) 0, xX, nên trên X: u, b – g(x) = 0. Hay
m
i i i
i 1
u (b g (x)) 0
= =
, xX
Hay i = 1,2,….,m: ui(bi – gi (x)) = 0, xX. (2.10)
Lê Văn Phi, Khoa Toan Thông kê, Tr ng Đai hoc Kinh tê Tp. Hô Chi Minhươ 47