ĐẠI HỌC THÁI NGUN
TRƯỜNG ĐẠI HỌC KHOA HỌC
TRỊNH DUY BÌNH
ĐIỀU KIỆN TỐI ƯU CHO NGHIỆM
HỮU HIỆU CỦA BÀI TỐN TỐI ƯU
ĐA MỤC TIÊU VỚI CÁC HÀM ỔN
ĐỊNH
LUẬN VĂN THẠC SĨ TỐN HỌC
Thái Ngun - 2013
Số hóa bởi Trung tâm Học liệu />ĐẠI HỌC THÁI NGUN
TRƯỜNG ĐẠI HỌC KHOA HỌC
TRINH DUY BÌNH
ĐIỀU KIỆN TỐI ƯU CHO NGHIỆM
HỮU HIỆU CỦA BÀI TỐN TỐI ƯU
ĐA MỤC TIÊU VỚI CÁC HÀM ỔN
ĐỊNH
Chun ngành : TỐN ỨNG DỤNG
Mã số : 60 46 01 12
LUẬN VĂN THẠC SĨ TỐN HỌC
Người hướng dẫn khoa học:
PGS.TS ĐỖ VĂN LƯU
Thái Ngun - 2013
Số hóa bởi Trung tâm Học liệu />Mục lục
Mở đầu 1
1 Hàm ổn định và đạo hàm tiếp liên 3
1.1 Hàm ổn định và đạo hàm tiếp liên . . . . . . . . . . . 3
1.2 Các quy tắc tính đạo hàm tiếp liên . . . . . . . . . . 11
2 Điều kiện tối ưu 18
2.1 Jacobian suy rộng Clarke . . . . . . . . . . . . . . . . 18
2.2 Các hàm vững . . . . . . . . . . . . . . . . . . . . . . 21
2.3 Các điều kiện tối ưu . . . . . . . . . . . . . . . . . . 30
2.4 Các quy tắc nhân tử Lagrange . . . . . . . . . . . . . 37
Kết luận 41
Tài liệu tham khảo 42
ii
Số hóa bởi Trung tâm Học liệu />Mở đầu
Lý thuyết các điều kiện tối ưu vectơ là một phần quan trọng
của lý thuyết tối ưu hóa. Người ta thiết lập các điều kiện tối ưu cho
các bài tốn tối ưu khơng trơn với các hàm Lipschitz địa phương
dưới ngơn ngữ các dưới vi phân khác nhau, chẳng hạn dưới vi phân
hàm lồi, các dưới vi phân Clarke, Michel Penot, Mordukhovich . . .
Lớp các hàm ổn định tại mỗi điểm của một tập rộng hơn lớp các
hàm Lipschitz địa phương trên tập đó. Bài tốn tối ưu vectơ với các
hàm ổn định được Jimenez - Novo [6] nghiên cứu và thiết lập các
điều kiện tối ưu dưới ngơn ngữ đạo hàm tiếp liên . Jimenez - Novo
[6] cũng chỉ ra đạo hàm tiếp liên có nhiều tính chất phong phú trong
lớp các hàm ổn định. Đây là vấn đề thời sự được nhiều tác giả quan
tâm nghiên cứu. Chính vì thế em chọn đề tài: " Điều kiện tối ưu cho
nghiệm hữu hiệu của bài tốn tối ưu đa mục tiêu với các hàm ổn
định "
Luận văn trình bày lý thuyết đạo hàm tiếp liên cho hàm ổn định,
hàm vững và các điều kiện tối ưu của Jimenez - Novo [6] cho bài
tốn tối ưu đa mục tiêu tổng qt trong khơng gian định chuẩn và
bài tốn tối ưu đa mục tiêu với ràng buộc nón và ràng buộc đẳng
thức trong khơng gian hữu hạn chiều.
Luận văn bao gồm phần mở đầu, hai chương, kết luận và danh
mục các tài liệu tham khảo.
Chương 1. Hàm ổn định và đạo hàm tiếp liên
Trình bày các kết quả nghiên cứu của Jimenez - Novo ([6], 2008)
về hàm ổn định và đạo hàm tiếp liên của các hàm ổn định, các quy
tắc tính đạo hàm tiếp liên trong lớp các hàm ổn định bao gồm: quy
1
Số hóa bởi Trung tâm Học liệu />tắc hàm hợp, đạo hàm tiếp liên của một tổng, một tích, một thương
hai hàm và max của một số hữu hạn hàm ổn định.
Chương 2. Điều kiện tối ưu
Trình bày các kết quả của Jimenez - Novo [6] về hàm vững, mối
quan hệ với hàm ổn định, hàm khả vi Hadamard, hàm Lipschitz địa
phương, đạo hàm tiếp liên của hàm vững, và các điều kiện tối ưu
cho bài tốn tối ưu vectơ tổng qt với các hàm ổn định và hàm
vững dưới ngơn ngữ đạo hàm tiếp liên. Chương 2 cũng trình bày các
điều kiện cần và đủ tối ưu dưới dạng các quy tắc nhân tử Lagrange
cho bài tốn tối ưu vectơ với ràng buộc nón và ràng buộc đẳng thức
trong khơng gian hữu hạn chiều, dưới ngơn ngữ đạo hàm tiếp liên.
Luận văn này được hồn thành tại trường Đại học Khoa học, Đại
học Thái Ngun dưới sự hướng dẫn tận tình của PGS TS. Đỗ Văn
Lưu. Em xin bày tỏ lòng biết ơn chân thành và sâu sắc về sự tận
tâm và nhiệt tình của Thầy trong suốt q trình em thực hiện luận
văn.
Em xin chân thành cảm ơn Ban Giám hiệu, phòng Đào tạo Khoa
học và Quan hệ quốc tế, Khoa Tốn - Tin trường Đại học Khoa học,
Đại học Thái Ngun đã quan tâm và giúp đỡ em trong suốt thời
gian học tập tại trường.
Em cũng xin gửi lời cảm ơn đến gia đình, bạn bè và các đồng
nghiệp đã động viên, giúp đỡ em trong q trình học tập của mình.
Do thời gian và kiến thức còn hạn chế nên luận văn khơng tránh
khỏi những thiếu sót. Em rất mong nhận được sự góp ý của các
thầy cơ để luận văn được hồn thiện hơn. Em xin chân thành cảm ơn!
Thái Ngun, ngày 15 tháng 7 năm 2013
Tác giả
Trịnh Duy Bình
2
Số hóa bởi Trung tâm Học liệu />Chương 1
Hàm ổn định và đạo hàm tiếp liên
Chương 1 trình bày các kết quả nghiên cứu của Jimenez - Novo
([6], 2008) về hàm ổn định và đạo hàm tiếp liên của các hàm ổn
định, các quy tắc tính đạo hàm tiếp liên trong lớp các hàm ổn định
bao gồm: quy tắc hàm hợp, đạo hàm tiếp liên của một tổng, một
tích, một thương hai hàm và max của một số hữu hạn hàm ổn định.
1.1 Hàm ổn định và đạo hàm tiếp liên
Giả sử X, Y là các khơng gian định chuẩn, M là một tập con của
X. Ký hiệu B (x
0
, δ) là hình cầu mở tâm x
0
, bán kính δ.
Kí hiệu intM, clM, coM, coneM tương ứng là phần trong của M,
bao đóng của M, bao lồi của M, nón sinh bởi M, và
Cone
+
M = {αx : α > 0, x ∈ M} .
Nón D ⊂ Y được gọi là nón nhọn nếu D ∩ (−D) = {0} .Trong
chương này ta giả thiết D là nón lồi đóng nhọn có phần trong khác
rỗng. Nón D sinh ra thứ tự trong khơng gian Y.
Định nghĩa 1.1.1.
Giả sử M ⊂ X và x
0
∈ X.
(a) Nón tiếp tuyến của M tại x
0
là
T (M, x
0
) =
v ∈ X : ∃t
n
→ 0
+
, ∃v
n
→ v
sao cho x
0
+ t
n
v
n
∈ M, ∀n ∈ N
;
3
Số hóa bởi Trung tâm Học liệu />(b) Nón các phương đạt được của M tại x
0
là
A(M, x
0
) =
v ∈ X : ∀t
n
→ 0
+
, ∃v
n
→ v
sao cho x
0
+ t
n
v
n
∈ M, ∀n ∈ N
;
(c) Nón tiếp tuyến phần trong của M tại điểm x
0
là
IT (M, x
0
) =
v ∈ X : ∃δ > 0 sao cho x
0
+ tu ∈ M
∀t ∈ (0, δ] , ∀u ∈ B (v, δ)
;
(d) Nón tiếp tuyến phần trong dãy của M tại x
0
là
IT
s
(M, x
0
) =
v ∈ X : ∃δ > 0, ∃t
n
→ 0
+
sao cho
x
0
+ t
n
u ∈ M ∀n ∈ N, ∀u ∈ B (v, δ)
.
Chú ý rằng ta ln có A (M, x
0
) ⊂ T (M, x
0
) .Tập hợp M được
gọi là khả dẫn xuất (derivable) tại điểm x
0
(xem [10]) nếu
A(M, x
0
) = T (M, x
0
).
Ta nói rằng hàm f : X → Y là khả dẫn xuất đồ thị (graph -
derivable) tại x
0
theo phương v ∈ X nếu với ∀y ∈ Y ,
(v, y) ∈ T (graphf, (x
0
, f(x
0
))) ⇒ (v, y) ∈ A (graphf, (x
0
, f(x
0
))),
trong đó graphf = {(x, y) ∈ X × Y : y = f (x)}.
Hàm f là khả dẫn xuất đồ thị tại x
0
nếu f là khả dẫn xuất đồ
thị tại x
0
theo mọi phương ∀v ∈ X khi và chỉ khi
T (graphf, (x
0
, f(x
0
))) = A (graphf, (x
0
, f(x
0
)))
tức là tập hợp graphf là khả dẫn xuất tại điểm (x
0
, f(x
0
)).
Sau đây là một ví dụ về hàm khả dẫn xuất đồ thị.
Ví dụ 1.1.1
Giả sử f : R → R được cho bởi f(x) = x sin
1
x
nếu x = 0 và
f(0) = 0. Khi đó f là khả dẫn xuất đồ thị tại 0.
Đạo hàm Hadamard của f : X → Y tại x
0
theo phương v ∈ X là
df(x
0
, v) = lim
(t,u)→(0
+
,v)
f(x
0
+ tu) − f(x
0
)
t
.
4
Số hóa bởi Trung tâm Học liệu />f là hàm khả vi Hadamard tại x
0
theo phương v nếu df(x
0
, v) tồn
tại, và f là hàm khả vi Hadamard tại x
0
nếu df (x
0
, v) tồn tại với
mọi v ∈ X. Nếu f khả vi Frechet tại x
0
, đạo hàm Frechet của nó
được kí hiệu bởi ∇f(x
0
).
Nếu hàm f : X → R thì đạo hàm Hadamard dưới và trên của f tại
x
0
được định nghĩa bởi
df (x
0
, v) = lim
(t,u)→(0
+
,v)
inf (f (x
0
+ tu) − f (x
0
))
t
tương ứng df (x
0
, v) = lim
(t,u)→(0
+
,v)
sup (f (x
0
+ tu) − f (x
0
))
t
.
Định nghĩa 1.1.2. [1]
Ta nói rằng f : X → Y là một hàm ổn định tại x
0
∈ X nếu tồn
tại một miền lân cận U của x
0
và k > 0 sao cho
f(x) − f(x
0
) ≤ kx − x
0
, ∀x ∈ U.
Nếu f(x) − f(x
) ≤ kx − x
, ∀x, x
∈ U, ta sẽ nói rằng f là
hàm Lipschitz địa phương tại x
0
.
Nếu với mỗi x
0
∈ X, tồn tại một lân cận U của x
0
sao cho f là
hàm Lipschitz địa phương tại x
0
thì ta nói rằng f là hàm Lipschitz
địa phương trên X.
Rõ ràng rằng hàm f Lipschitz địa phương trên X kéo theo f là
hàm ổn định tại mỗi x
0
∈ X, nhưng điều ngược lại khơng đúng.
Ví dụ hàm f trong Ví dụ (1.1.1) là ổn định với mỗi x
0
∈ R nhưng f
khơng là hàm Lipschitz địa phương tại 0.
Chúng ta sẽ sử dụng đạo hàm theo phương sau đây cho ánh xạ
đa trị (xem [8]).
Định nghĩa 1.1.3.
Giả sử f : X → Y và x
0
, v ∈ X. Đạo hàm tiếp liên của f tại x
0
theo phương v là tập hợp sau đây trong Y:
∂
∗
f (x
0
) v =
y ∈ Y : ∃ (t
n
, v
n
) → (0
+
, v) sao cho
lim
n→∞
(f (x
0
+ t
n
v
n
) − f (x
0
)) /t
n
= y
.
5
Số hóa bởi Trung tâm Học liệu />Chú ý rằng ∂
∗
f(x
0
)v được gọi là đạo hàm tiếp liên của ánh xạ đa
trị x → {f (x)} . Do đó ∂
∗
f(x
0
)v là một tập hợp đóng và ánh xạ đa
trị v → ∂
∗
f (x
0
) v thuần nhất dương.
Nếu f khả vi Hadamard tại x
0
thì
∂
∗
f(x
0
)v = {df(x, v)} ∀v ∈ X.
Cho f : X → Y và g : X → Z, trong đó Z là một khơng gian
định chuẩn khác.
Hàm (f, g) : X → Y × Z được xác định bởi (f, g)(x) = (f(x), g(x)).
Trong khơng gian Y × Z, ta lấy (y, z) = y + z.
Định lý sau (cách chứng minh rất đơn giản) cho ta một số tính
chất của hàm (f, g) liên quan đến tính ổn định và đạo hàm tiếp liên.
Mệnh đề 1.1.1. [6]
(i) f và g ổn định tại x
0
⇔ (f, g) ổn định tại x
0
.
(ii) Với mỗi v ∈ X, ∂
∗
(f, g) (x
0
) v ⊂ ∂
∗
f (x
0
) v × ∂
∗
g (x
0
) v. Nếu f
hoặc g khả vi Hadamard tại x
0
theo phương v thì ta có dấu bằng.
Bao hàm thức trong (ii) có thể đúng như ví dụ 1.1.2 dưới đây.
Nhận xét 1.1.1.
(1) Rõ ràng là nếu f : X → R
p
là ổn định taị x
0
thì
∂
∗
f (x
0
) v = φ, ∀v ∈ X.
Thật vậy với mỗi dãy t
n
→ 0
+
và v
n
→ v, (với n đủ lớn) ta có
f(x
0
+ t
n
v
n
) − f(x
0
)
t
n
≤ kv
n
≤ k(v + 1).
bởi vì y
n
= (f (x
0
+ t
n
v
n
) − f (x
0
)) /t
n
là một dãy bị chặn, tồn
tại một dãy con (ta vẫn kí hiệu như cũ) sao cho
y
n
→ y ∈ ∂
∗
f(x
0
)v.
(2) Nếu f : X → Y là ổn định taị x
0
với hằng số k trên một lân
cận của x
0
thì ∂
∗
f(x
0
)v ≤ kv. Hơn nữa, nếu Y = R
p
, thì
∂
∗
f(x
0
)v là một tập hợp compact với mọi v ∈ X.
6
Số hóa bởi Trung tâm Học liệu />(3) Nếu Y = R và df(x
0
, v), df(x
0
, v) hữu hạn thì
{df(x
0
, v), df(x
0
, v)} ⊂ ∂
∗
f(x
0
)v ⊂ [df(x
0
, v), df(x
0
, v)].
Chiều ngược lại của bao hàm thức cũng đúng khi f liên tục
trong một lân cận của x
0
.
(4) Nếu Y = R
p
thì
∂
∗
f (x
0
) v ⊂
p
i=1
∂
∗
f
i
(x
0
) v ⊂
p
i=1
df
i
(x
0
, v) , df
i
(x
0
, v)
.
∂
∗
f(x
0
)v là một tập hợp lồi khi Y = R và f liên tục trong
một lân cận của x
0
. Điều này khơng đúng cho khơng gian định
chuẩn tùy ý Y , thậm chí nếu f là hàm Lipschitz địa phương tại
x
0
.
Ví dụ 1.1.2
(a) Giả sử A = {2
−n
: n ∈ N} ∪{0} và f : R → R được xác định bởi
f(x) = dist(x, A). Khi đó f là một hàm Lipschitz tồn cục
trên R hằng số 1. Nhưng f khơng khả vi Hadamard tại x
0
= 0
theo phương v = 1.
Thật vậy, ta có
0 ≤ f(x) ≤
x
3
, ∀x ∈
0,
1
2
,
f(2
−n
) = 0 và f(s
n
) =
s
n
3
, ∀n ∈ N, trong đó s
n
là điểm giữa
của đoạn [2
−n−1
, 2
−n
], tức là s
n
= 3 · 2
−n−2
.
Do đó df(x
0
, v) = 0 và df(x
0
, v) =
1
3
. Vì thế, ∂
∗
f(x
0
)v = [0,
1
3
].
(b) Giả sử
B = {0}
∞
n=1
7
8
· 2
−n
,
5
4
· 2
−n
và chúng ta xác định g : R → R bởi g(x) = 2dist(x, B). Khi đó,
g là Lipschitz tồn cục và thỏa mãn
0 ≤ g (x) ≤ f (x) ≤ x/3, ∀v ∈ [0, 1/2]
7
Số hóa bởi Trung tâm Học liệu />g = 0 trên B, và g (s
n
) = s
n
/3 ∀n ∈ N.
Các đạo hàm Hadamard dưới và trên theo phương v = 1 là
dg (x
0
, v) = 0 và dg (x
0
, v) =
1
3
.
Do đó, ∂
∗
g(x
0
)v =
0,
1
3
.
(c) Giả sử h = (f, g). Khi đó, h là Lipschitz tồn cục và ta có
∂
∗
h(x
0
)v =
co{(0, 0), (
1
5
, 0)}
∪
co{(
1
5
, 0), (
1
3
,
1
3
)}
.
Tập này khơng lồi trong R
2
.
Những kết quả sau đây cho ta những tính chất của hàm ổn định
và đạo hàm tiếp liên.
Mệnh đề 1.1.2.
f : X → Y là ổn định tại x
0
khi và chỉ khi df(x
0
, 0) = 0.
Chứng minh.
(=⇒) Ta suy ra trực tiếp từ những định nghĩa bởi vì với δ > 0 nào
đó,
f (x
0
+ tu) − f (x
0
)
t
≤ k u , ∀t ∈ (0, δ) , ∀u ∈ B (0.1) .
(⇐= ) Giả sử > 0 cố định. Theo giả thiết ∃δ
1
, δ
2
> 0 sao cho
f(x
0
+ tu) − f(x
0
)
t
< ∀t ∈ (0, δ
1
), u ∈ B(v, δ
2
). (1.1)
Chọn a cố định, η ∈ (0, δ
2
) và δ = η · δ
1
. Khi đó,
B(0, δ) \ {0} ⊂ (0, δ
1
) · S
η
⊂ (0, δ
1
) · B(0, δ
2
),
trong đó S
η
= {x ∈ X : x = η}. Thật vậy, cho w ∈ B(0, δ) \ {0},
ta xác định t = w/η và u = t
−1
w. Khi đó, t < δ
1
và u = η.
Vì vậy, w = t · u ∈ (0, δ
1
) · δ
η
. Bao hàm thức thứ hai là hiển nhiên.
Nếu 0 < x − x
0
< δ, nó sẽ là w = x − x
0
= t · u
với t ∈ (0, δ
1
), u = η
8
Số hóa bởi Trung tâm Học liệu />Do đó,
f (x) − f (x
0
)
x − x
0
=
f (x
0
+ w) − f (x
0
)
w
=
f (x
0
+ tu) − f (x
0
)
t
·
1
u
< ε
1
n
= k
(bất đẳng thức cuối đúng do phương trình (1.1)).
Mệnh đề 1.1.3.
Giả sử f : X → Y, x
0
, v ∈ X và y ∈ Y
(i) Nếu df(x
0
, v) = y thì ∂
∗
f(x
0
)v = {y}.
(ii) Nếu Y là hữu hạn chiều, f ổn định tại x
0
và ∂
∗
f(x
0
)v = {y} thì
df(x
0
, v) = y.
Chứng minh.
Phần (i) là hiển nhiên. Ta chứng minh phần (ii)
Ta chứng minh rằng nếu
y
n
:= (f(x
0
+ t
n
v
n
) − f(x
0
))/t
n
→ y, ∀(t
n
, v
n
) → (0
+
, v)
thì df(x
0
, v) = y. Bởi vì f ổn định tại x
0
dãy (y
n
) bị chặn,
tức là y
n
∈ clB(0, r), ∀n ∈ N với r > 0, nào đó. Với mỗi k ∈ N ta
xác định
n(k) = max{n ∈ N : y
n
∈ B(y,
1
k
)}.
Max là đạt được, bởi vì nếu khơng thì có vơ hạn y
n
∈ B(y,
1
k
). Do
đó, y
∈ B(y,
1
k
) (vì clB(0, r) \ B(y, 1/k) là một tập compact).
Nhưng y
∈ ∂
∗
f (x
0
) v = {y}. Ta đi đến mâu thuẫn. Hiển nhiên,
y
n
→ y.
Nếu f chỉ liên tục tại x
0
thì Mệnh đề 1.1.3(ii) là sai. Ví dụ sau
minh họa điều đó.
9
Số hóa bởi Trung tâm Học liệu />Ví dụ 1.1.3
Cho f : R → R được cho bởi
f(x) =
1
n
, nếu x =
1
n
, n ∈ N \ {0}
| x |, nếu x =
1
n
, n ∈ N \ {0}
Ta có ∂
∗
f(0)1 = {1}, nhưng f khơng khả vi Hadamard tại 0 theo
phương 1.
Kết quả tiếp theo là hệ quả trực tiếp của Mệnh đề 1.1.2 và 1.1.3(i).
Hệ quả 1.1.1.
Khi f : X → Y là ổn định tại x
0
thì ∂
∗
f(x
0
)0 = {0}.
Điều ngược lại cũng đúng đối với hàm liên tục khi Y hữu hạn
chiều.
Mệnh đề 1.1.4.
Giả sử f : X → R
p
liên tục tại x
0
. Nếu ∂
∗
f(x
0
)0 = {0} thì f ổn
định tại x
0
.
Chứng minh.
Giả sử f khơng ổn định tại x
0
. Khi đó, tồn tại dãy x
n
→ x
0
sao
cho
f(x
n
− f(x
0
) > nx
n
− x
0
, ∀n ∈ N. (1.2)
Đặt t
n
= f(x
n
) − f(x
0
) > 0 và v
n
= (x
n
− x
0
)/t
n
. Vì f liên tục tại
x
0
, ta có t
n
→ 0
+
, và từ (1.2), ta suy ra v
n
→ 0.
Lấy dãy con nếu cần thiết, ta có thể giả sử rằng
y
n
:= (f(x
0
+t
n
v
n
)−f(x
0
))/t
n
hội tụ tới y nào đó y ∈ R
p
với y = 1
bởi vì y
n
= 1. Do đó, y ∈ ∂
∗
f(x
0
)0. Điều này mâu thuẫn với giả
thiết.
Nếu ∂
∗
f(x
0
)0 = {0} thì f khơng ổn định tại x
0
do hệ quả 1.1.1
nhưng f có thể liên tục tại x
0
.
Chẳng hạn f : R → R, f(x) =
| x |, x
0
= 0. Lấy t
n
= 1/n và
v
n
= 1/n, ta có lim
x→∞
f (t
n
v
n
) /t
n
= 1 ∈ ∂
∗
f (x
0
) 0.
10
Số hóa bởi Trung tâm Học liệu />Hệ quả 1.1.2.
f : X → R
p
liên tục tại x
0
và ∂
∗
f (x
0
) 0 = {0} và tương đương với
df (x
0
, 0) = 0.
Hệ quả 1.1.2 được suy ra từ các Mệnh đề 1.1.4 và 1.1.2.
1.2 Các quy tắc tính đạo hàm tiếp liên
Tiếp theo chúng ta sẽ trình bày một số quy tắc tính tốn đạo hàm
tiếp liên. Giả sử rằng Z là một khơng gian định chuẩn.
Định lý 1.2.1. (Quy tắc hàm hợp)
Giả sử f : X → Y và g : Y → Z khả vi Hadamard tại f(x
0
).
Khi đó,
dg(f(x
0
), ·)(∂
∗
f(x
0
)v) ⊂ ∂
∗
(g ◦ f)(x
0
)v, ∀v ∈ X (1.3)
Hơn nữa, nếu Y hữu hạn chiều và f ổn định tại x
0
thì g ◦ f ổn
định tại x
0
và có đẳng thức trong (1.3).
Chứng minh.
Giả sử y ∈ ∂
∗
f(x
0
)v. Khi đó, ∃(t
n
, v
n
) → (0
+
, v) sao cho
y
n
:= (f (x
0
+ t
n
v
n
) − f (x
0
)) /t
n
→ y.
Bởi vì g khả vi Hadamard tại f(x
0
), cho nên dg(f(x
0
), y) tồn tại và
ta có
g(f(x
0
) + tu) = g(f(x
0
)) + tdg(f(x
0
), y) + tα(t, u),
trong đó
lim
(t,u)→(0
+
,y)
α (t, u) = 0.
Đặc biệt, với t = t
n
và u = y
n
thì f(x
0
+ t
n
v
n
) = f(x
0
) + t
n
y
n
, và ta
có
g (f (x
0
+ t
n
v
n
)) = g (f (x
0
) + t
n
y
n
)
= g (f (x
0
)) + t
n
dg (f (x
0
) , y) + t
n
α (t
n
, y
n
) .
11
Số hóa bởi Trung tâm Học liệu />Do đó,
lim
n→∞
(g ◦ f) (x
0
+ t
n
v
n
) − (g ◦ f) (x
0
)
t
n
= dg (f (x
0
) , y) .
Điều này kéo theo
dg(f(x
0
), y) ∈ ∂
∗
(g ◦ f)(x
0
)v.
Tiếp theo ta chứng minh tính ổn định của g ◦ f, dễ thấy vì g cũng
ổn định tại f(x
0
) do Mệnh đề 1.1.2.
Giả sử z ∈ ∂
∗
(g ◦ f)(x
0
)v. Khi đó tồn tại (t
n
, v
n
) → (0
+
, v)
sao cho
g(f(x
0
+ t
n
v
n
)) − g(f(x
0
))
t
n
→ z. (1.4)
Do f ổn định tại x
0
và Y hữu hạn chiều, ta có thể giả sử rằng:
(lấy một dãy con nếu cần)
y
n
:=
f(x
0
+ t
n
v
n
) − f(x
0
)
t
n
→ y. (1.5)
Do đó y ∈ ∂
∗
f(x
0
)v. Do g khả vi Hadamard tại f (x
0
),
và f(x
0
+ t
n
v
n
) = f(x
0
) + t
n
y
n
và sử dụng (1.4) và (1.5) ta suy ra
dg (f (x
0
) , y) = lim
n→∞
g (f (x
0
) + t
n
y
n
) − g (f (x
0
))
t
n
= lim
x→∞
g (f (x
0
+ t
n
v
n
)) − g (f (x
0
))
t
n
= z.
Vì vậy, z ∈ dg(f(x
0
), ·)(∂
∗
f(x
0
)v).
Kết quả sau đây suy ra trực tiếp từ Định lý 1.2.1, khi lấy g = ψ.
Chú ý rằng dψ(y, u) = ψ
0
(u), ∀y, u ∈ Y .
Hệ quả 1.2.1.
Giả sử f : X → Y , và ψ : Y → Z là một hàm affine liên tục được
cho bởi ψ(y) = ψ
0
(y) + b với ψ
0
: Y → Z là hàm tuyến tính liên tục
và b ∈ Z. Khi đó,
ψ
0
(∂
∗
f(x
0
)v) ⊂ ∂
∗
(ψ ◦ f)(x
0
)v, ∀v ∈ X. (1.6)
Hơn nữa, nếu Y hữu hạn chiều và f ổn định tại x
0
thì ψ ◦ f là ổn
định tại x
0
và đẳng thức xảy ra trong (1.6).
12
Số hóa bởi Trung tâm Học liệu />Chú ý rằng khi Y hữu hạn chiều và ψ : Y → Z là một hàm affine
thì ψ liên tục trên Y .
Mệnh đề 1.2.1. (Quy tắc tổng)
Giả sử f, g là hai hàm từ X vào Y . Khi đó,
{y + z : (y, z) ∈ ∂
∗
(f, g) (x
0
) v} ⊂ ∂
∗
(f + g) (x
0
) v.
Nếu Y hữu hạn chiều và f (hoặc g) ổn định tại x
0
thì
∂
∗
(f + g)(x
0
)v ⊂ ∂
∗
f(x
0
)v + ∂
∗
g(x
0
)v. (1.7)
Nếu f (hoặc g) khả vi Hadamard tại x
0
theo phương v thì đẳng thức
xảy ra trong (1.7).
Chứng minh.
Phần thứ nhất suy ra từ Hệ quả 1.2.1 với (f, g) : X → Y ×Y thay
cho f và ψ : Y × Y → Y được cho bởi cơng thức ψ(y, z) = y + z.
Phần thứ hai suy ra từ định nghĩa và chú ý đến Nhận xét 1.2.1(1)
và phần thứ ba suy ra từ Định nghĩa.
Các hàm f và g trong Ví dụ 1.2.1(c) cho thấy bao hàm thức (1.7)
có thể chặt.
Để minh họa quy tắc tổng này, ta xét h : X × Y → R
p
được cho bởi h(x, y) = f(x) + g(y),
trong đó f : X → R
p
, g : Y → R
p
và (x
0
, y
0
), (v, u) ∈ X × Y .
Nếu f ổn định tại x
0
thì
∂
∗
h(x
0
, y
0
)(v, u) ⊂ ∂
∗
f(x
0
)v + ∂
∗
g(y
0
)u,
và đẳng thức xảy ra nếu f khả vi Hadamard tại x
0
theo phương v.
Để chứng minh điều này, ta xét f và g từ X × Y vào R
p
f(x, y) = f(x) và g(x, y) = g(y) và khi đó áp dụng Mệnh đề 1.2.1
cho h = f + g, ta nhận được
∂
∗
f (x
0
, y
0
) (v, u) = ∂
∗
f (x
0
) v,
và
∂
∗
g (x
0
, y
0
) (v, u) = ∂
∗
g (y
0
) u
13
Số hóa bởi Trung tâm Học liệu />Mệnh đề 1.2.2. (Quy tắc tích)
Giả sử f, g là hai hàm từ X vào R. Khi đó,
{g(x
0
)y + f(x
0
)z : (y, z) ∈ ∂
∗
(f, g)(x
0
)v} ⊂ ∂
∗
(fg)(x
0
)v. (1.8)
Nếu f ổn định tại x
0
, và hoặc f(x
0
) = 0 hoặc g liên tục tại x
0
, thì
∂
∗
(fg)(x
0
)v ⊂ g(x
0
)∂
∗
f(x
0
)v + f(x
0
)∂
∗
g(x
0
)v. (1.9)
Hơn nữa, nếu f khả vi Hadamard tại x
0
theo phương v thì đẳng thức
xảy ra trong (1.8) và (1.9).
Nếu f, g ổn định tại x
0
thì fg ổn định x
0
, và đẳng thức này xảy
ra trong (1.8). Hơn nữa,
∂
∗
(fg)(x
0
)v = ∂
∗
(g(x
0
)f + f(x
0
)g)(x
0
)v. (1.10)
Chứng minh.
Giả sử h : R
2
→ R là một hàm h(y, z) = yz. Hàm này khả vi phân
Frechet, với ∇h(y
0
, z
0
)(u, w) = z
0
u + y
0
w. Chú ý rằng fg = h ◦ (f, g)
và áp dụng Định lý 1.2.1 để nhận được phần 1, tức là bao hàm thức
(1.8).
Giả sử rằng f ổn định tại x
0
.
Giả sử x
n
= x
0
+ t
n
v
n
, với (t
n
, v
n
) → (o
+
, v). Khi đó,
w
n
: =
f (x
n
) g (x
n
) − f (x
0
) g (x
0
)
t
n
= g (x
0
)
f (x
n
) − f (x
0
)
t
n
+ f (x
n
)
g (x
n
) − g (x
0
)
t
n
.
(1.11)
Nếu (w
n
) hội tu tới w ∈ ∂
∗
(fg)(x
0
)v, thì do Nhận xét 1.1.1(1)
(ta lấy thêm 1 dãy con nếu cần thiết), ta nhận được
y
n
:= (f(x
n
) − f(x
0
))/t
n
→ y ∈ ∂
∗
f(x
0
)v.
Nếu f(x
0
) = 0, do f(x
n
) → f(x
0
), ta có f(x
n
) = 0, ∀n đủ lớn và
sủ dụng (1.11), ta có
g(x
n
) − g(x
0
)
t
n
=
w
n
− g(x
0
)y
n
f(x
n
)
→
w − g(x
0
)y
f(x
0
)
∈ ∂
∗
g(x
0
)v.
14
Số hóa bởi Trung tâm Học liệu />Nếu f(x
0
) = 0 và g liên tục tại x
0
, vì y
n
= f(x
n
)/t
n
→ y, ta suy
ra
w
n
=
f(x
n
)
t
n
g(x
n
) → w = g(x
0
)y ∈ g(x
0
)∂
∗
f(x
0
)v + 0∂
∗
g(x
0
)v.
Đối với phần thứ ba, khi f khả vi Hadamard tại x
0
theo phương
v, từ Mệnh đề 1.1.1(ii) ta có
∂
∗
(f, g)(x
0
)v = ∂
∗
f(x
0
)v × ∂
∗
g(x
0
)v.
Do đó,
{g(x
0
)y + f(x
0
)z : (y, z) ∈ ∂
∗
(f, g)(x
0
)v}
= g(x
0
)∂
∗
f(x
0
)v + f(x
0
)∂
∗
g(x
0
)v.
Vì thế mà ta có đẳng thức trong (1.8) và (1.9).
Đối với phần cuối, đẳng thức trong (1.8) và tính ổn định của fg
được chứng minh như trong phần đầu tiên bằng cách áp dụng phần
thứ hai của Định lý 1.2.1 vì (f, g) ổn định bởi mệnh đề 1.1.1(i).
Để chứng minh đẳng thức (1.10) ta hãy xét phiếm hàm tuyến tính
ψ : R
2
→ R được định nghĩa bởi : ψ(y, z) = g(x
0
)y + f(x
0
)z.
Sử dụng Hệ quả 1.2.1, ta nhận được
ψ(∂
∗
(f, g)(x
0
)v) = ∂
∗
(ψ ◦ (f, g))(x
0
)v
= ∂
∗
(g(x
0
)f + f(x
0
)g)(x
0
)v.
Tuy nhiên
ψ (∂
∗
(f, g) (x
0
) v) = {g (x
0
) y + f (x
0
) z : (y, z) ∈ ∂
∗
(f, g) (x
0
) v}
= ∂
∗
(fg) (x
0
) v
như chúng ta vừa chỉ ra và chứng minh được hồn thành.
Ví dụ sau đây cho thấy bao hàm thức trong Mệnh đề 1.2.2 có thể
là chặt.
Ví dụ 1.2.1
Giả sử f : R → R, g : R → R, x
0
= 0 và v = 1. Để cho ngắn gọn,
15
Số hóa bởi Trung tâm Học liệu />ta xét
A = {g(x
0
)y + f(x
0
)z : (y, z) ∈ ∂
∗
(f, g)(x
0
)v},
B = ∂
∗
(fg) (x
0
) v,
C = g (x
0
) ∂
∗
f (x
0
) v + f (x
0
) ∂
∗
g (x
0
) v,
D = ∂
∗
(g (x
0
) f + f (x
0
) g) (x
0
) v.
(a) f(x) = x, g(x) = sin(1/x), nếu x = 0 và g(0) = 1. Khi đó f khả
vi Hadamard tại x
0
nhưng khơng liên tục tại x
0
,
A = C = D = {1} và B = [−1, 1].
Các bao hàm thức (1.9) và (1.10) khơng thỏa mãn và (1.8) chặt.
(b) f(x) = x sin(1/x), nếu x = 0, f(0) = 0; g(x) = sin(1/x),
nếu x = 0, g(0) = 1. Khi đó, f ổn định tại x
0
, nhưng khơng
khả vi Hadamard tại x
0
theo phương v, A = {1}, B = [0, 1] và
C = D = [−1, 1]. (1.9) thỏa mãn với bao hàm thức chặt, nhưng
(1.10) lại khơng thỏa mãn.
(c) f(x) = 1 +x sin(1/x), g(x) = 1 −x sin(1/x), và f(0) = g(0) = 0.
Khi đó, f và g ổn định tại x
0
,
∂
∗
(f, g)(x
0
)v = {(y, −y) : −1 ≤ y ≤ 1},
∂
∗
f(x
0
)v = ∂
∗
g(x
0
)v = [−1, 1], A = B = D = {0} và C = [- 2,2];
(1.8) và (1.10) thỏa mãn đẳng thức, và (1.9) thỏa mãn với bao
hàm thức chặt.
Tương tự vào Mệnh đề 1.2.2, ta có mệnh đề sau đây.
Mệnh đề 1.2.3. (Quy tắc thương)
Giả sử f, g là hai hàm từ X vào R với g(x
0
) = 0.
Khi đó,
g(x
0
)y − f(x
0
)z
g(x
0
)
2
: (y, z) ∈ ∂
∗
(f, g)(x
0
)v
⊂ ∂
∗
(f/g)(x
0
)v. (1.12)
Nếu g ổn định tại x
0
(hoặc f ổn định tại x
0
, và hoặc f(x
0
) = 0 hoặc
g liên tục tại x
0
) thì
∂
∗
(f/g)(x
0
)v ⊂
g(x
0
)∂
∗
f(x
0
)v − f(x
0
)∂
∗
g(x
0
)v
g(x
0
)
2
. (1.13)
16
Số hóa bởi Trung tâm Học liệu />Hơn nữa, nếu g khả vi Hadamard tại x
0
theo phương v thì ta có đẳng
thức trong (1.12) và (1.13).
Nếu f và g ổn định tại x
0
thì f/g ổn định tại x
0
, và ta có đẳng
thức trong (1.12).
Hơn nữa,
∂
∗
(f/g)(x
0
)v =
1
g(x
0
)
2
∂
∗
(g(x
0
)f − f(x
0
)g)(x
0
)v.
Mệnh đề 1.2.4. (Quy tắc hàm max)
Giả sử f = (f
1
, f
2
, . . . , f
p
) : X → R
p
và g : X → R được xác định
bởi
g(x) = max {f
i
(x) : i = 1, 2, . . . , p} .
Khi đó,
α ∈ R : ∃y = (y
1
, y
2
, . . . , y
p
) ∈ ∂
∗
f (x
0
) v.
sao cho α = max {y
i
: i = I (f (x
0
))}
⊂ ∂
∗
g (x
0
) v,
trong đó I (f (x
0
)) = {i ∈ {1, 2, . . . , p} : g (x
0
) = f
i
(x
0
)}. Nếu f ổn
định tại x
0
, thì g ổn định tại x
0
và đẳng thức này xảy ra.
Chứng minh.
Ta xác định hàm h : R
p
→ R trong đó
h (y
1
, y
2
, . . . , y
p
) = max {y
i
: i = 1, 2, . . . , p}.
Hiển nhiên rằng g = h ◦ f. Hàm h khả vi Hadamard tại y
0
∈ R
p
với
dh (y
0
, u) = max {u
i
: i ∈ I (y
0
)}
ta áp dụng Định lý 1.2.1 và nhận được điều phải chứng minh.
17
Số hóa bởi Trung tâm Học liệu />Chương 2
Điều kiện tối ưu
Chương 2 trình bày các kết quả của Jimenez - Novo [6] về hàm
vững, mối quan hệ với hàm ổn định, hàm khả vi Hadamard, hàm
Lipschitz địa phương, đạo hàm tiếp liên của hàm vững, và các điều
kiện tối ưu cho bài tốn tối ưu vectơ tổng qt với các hàm ổn định
và hàm vững dưới ngơn ngữ đạo hàm tiếp liên. Trong chương này,
chúng tơi cũng trình bày các điều kiện cần và đủ tối ưu dưới dạng
các quy tắc nhân tử Lagrange cho bài tốn tối ưu vectơ với ràng
buộc nón và ràng buộc đẳng thức trong khơng gian hữu hạn chiều,
dưới ngơn ngữ đạo hàm tiếp liên.
2.1 Jacobian suy rộng Clarke
Cho hàm vectơ f : R
n
→ R
m
. Khi đó hàm f có dạng:
f (x) = (f
1
(x) , . . . , f
m
(x))
Giả sử mỗi hàm f
i
là Lipschitz địa phương tại x. Khi đó, hàm f là
Lipschitz địa phương tại x.
Định lí Rademacher khẳng định rằng hàm f khả vi hầu khắp nơi
trong một lân cận của x (mà trên đó f Lipschitz), tức là mỗi f
i
khả
vi hầu khắp nơi trong lân cận đó của x.
Ký hiệu Ω
f
là tập các điểm mà hàm f khơng khả vi. Nếu tại x
hàm f có đạo hàm, thì ta ký hiệu Jf (x) là m×n- ma trận Jacobian
thơng thường của f.
18
Số hóa bởi Trung tâm Học liệu />Định nghĩa 2.1.1.
Jacobian suy rộng của hàm vectơ f tại x, ký hiệu ∂f (x), được
định nghĩa như sau:
∂f (x) := co {lim Jf (x
i
) : x
i
→ x, x
i
/∈ Ω
f
}
Ta xét khơng gian các m × n- ma trận R
m×n
với chuẩn:
M
m×n
=
m
i=1
r
i
2
1/2
trong đó r
i
là hàng thứ i của M. Ký hiệu B
m×n
là hình cầu đơn vị
mở trong R
m×n
.
Định lý 2.1.1. [1]
Giả sử mỗi hàm f
i
Lipschitz địa phương tại x (i = 1, . . . , m)
khi đó,
(i) ∂f (x) = ∅, lồi, compắc trong R
m×n
;
(ii) Ánh xạ đa trị ∂f là đóng tại x, tức là nếu
x
i
→ x, Z
i
∈ ∂f (x
i
) , Z
i
→ Z, thì Z ∈ ∂f (x) ;
(iii) ∂f nửa liên tục trên tại x: ∀ε > 0, ∃δ > 0, ∀x ∈ x + δB
∂f (x) ⊂ ∂f (x) + εB
m×n
;
(iv) Nếu f
i
Lipschitz địa phương tại x với hằng số K
i
, thì
f Lipschitz địa phương tại x với hằng số K = (K
1
, . . . , K
m
)
và ∂f (x) ⊂ KB
m×n
;
(v) ∂f (x) ⊂ ∂f
1
(x) × ∂f
2
(x) × . . . × ∂f
m
(x), trong đó
∂f
1
(x) × . . . × ∂f
m
(x)
là tập tất cả m × n- ma trận có hàng thứ i thuộc ∂f
i
(x) , ∀i.
Nếu m = 1 thì ∂f (x) = ∂f
1
(x) (tức là Jacobian suy rộng trùng với
gradient suy rộng).
19
Số hóa bởi Trung tâm Học liệu />Định lý 2.1.2. (Định lí về giá trị trung bình cho hàm vectơ) [1]
Giả sử f là hàm vectơ Lipschitz trên tập lồi mở U ⊂ R
n
; x, y ∈ U
Khi đó,
f (y) − f (x) ∈ co ∂f ([x, y]) (y − x) (2.1)
Vế phải của (2.1) có nghĩa là bao lồi của tất cả các điểm dạng
Z (y − x),
trong đó Z ∈ ∂f (u), với u nào đó thuộc [x, y] và
[co ∂f ([x, y])] (y − x) = co [∂f ([x, y]) (y − x)]
Định lý 2.1.3. (Jacobian suy rộng của hàm hợp)[1]
Giả sử h : R
n
→ R
m
Lipschitz địa phương tại x, g : R
m
→ R
Lipschitz địa phương tại h (x) , f = g ◦ h
Khi đó, f cũng Lipschitz địa phương tại x, và
∂f (x) ⊂ co {∂g (h (x)) ∂h (x)} . (2.2)
Hơn nữa, nếu g khả vi chặt Hadamard tại h (x), thì (2.2) có dấu
bằng, và phép tốn co có thể bỏ được.
Bây giờ ta sẽ so sánh đạo hàm liên tục với Jacobian suy rộng. Giả
sử f : R
n
→ R
p
là hàm Lipschitz địa phương tại x
0
.
Quan hệ giữa Jacobian suy rộng Clarke và ∂
∗
f (x
0
) v được đưa ra
trong Mệnh đề 2.1.1. Trước hết ta cần một kết quả sau.
Bổ đề 2.1.1.
Giả sử f : R
n
→ R
p
là một Lipschitz địa phương tại x
0
và v ∈ R
n
.
Nếu lim
n→∞
(x
n
− x
0
) /t
n
= v với t
n
→ 0
+
, thì (ta chọn một dãy con
nếu cần thiết)
lim
n→∞
(f (x
n
) − f (x
0
)) t
n
= Av với A nào đó ∈ ∂f (x
0
) .
Chứng minh.
Giả sử v
n
được xác định bởi đẳng thức x
n
= x
0
+ t
n
v
n
thì v
n
→ v.
Sử dụng định lý giá trị trung bình (Định lý 2.1.2) ta nhận được
f (x
n
) − f (x
0
) = A
n
t
n
v
n
20
Số hóa bởi Trung tâm Học liệu />với A
n
∈ co {∂f (x) : x ∈ [x
0
, x
n
]}. Do ánh xạ đa trị x → ∂f (x) là
nửa liên tục trên, lồi và giá trị compact, ta có dãy A
n
hội tụ tới
A ∈ ∂f (x
0
).
Mệnh đề 2.1.1.
Nếu f : R
n
→ R
p
là một Lipschitz địa phương tại x
0
thì
∂
∗
f (x
0
) v ⊂ ∂f (x
0
) v ∀v ∈ R
n
Chứng minh của mệnh đề này suy trực tiếp từ định nghĩa của
∂
∗
f (x
0
) v và Bổ đề 2.1.1. Dễ thấy rằng bao hàm thức trên nói chung
là chặt. Chẳng hạn, cho f : R → R được cho bởi f (x) = |x|, ta có
∂
∗
f (0) 1 = {1} , ∂f (0) = [−1, 1]
và
∂f (0) 1 = [−1, 1] .
Lợi thế của ∂f (x
0
) v là tính lồi, trong khi đó ∂
∗
f (x
0
) v nói chung
khơng lồi (xem Ví dụ 1.1.2).
2.2 Các hàm vững
Định nghĩa 2.2.1. [9]
Hàm f : X → Y được gọi là vững tại x
0
∈ X theo phương v ∈ X,
hay nói gọn là f vững tại (x
0
, v) nếu
lim
(t,u)→(0
+
,v)
f (x
0
+ tu) − f (x
0
+ tv)
t
= 0.
Đẳng thức trên tương đương với
lim
(t,v
,v
)→(0
+
,v,v)
f (x
0
+ tv
) − f (x
0
+ tv
)
t
= 0.
f được gọi là vững tại x
0
nếu f vững tại x
0
theo mọi phương.
Đặc biệt, f vững tại (x
0
, 0) khi và chỉ khi f ổn định tại x
0
(Mệnh
đề 1.1.2).
Nhắc lại rằng f vững tại x
0
khi f là Lipschitz địa phương tại x
0
.
21
Số hóa bởi Trung tâm Học liệu />Tương tự, f vững tại x
0
khi f khả vi Hadamard tại x
0
. Ba khái
niệm này hồn tồn khác nhau, lớp các hàm vững rộng hơn lớp của
Lipschitz hoặc khả vi Hadamard (xem Ví dụ 2.2.1).
Tính vững tại (x
0
, v) với v = 0 khơng phụ thuộc vào giá trị của f
tại x
0
, tức là nếu f vững tai (x
0
, v) với v = 0, g (x) = f (x) , ∀x = x
0
và g (x
0
) = y với y tùy ý y ∈ Y thì g vững tại (x
0
, v).
Nếu f : X → Y và g : X → Z. Rõ ràng là f và g vững tại (x
0
, v)
nếu và chỉ nếu (f, g) vững tại (x
0
, v).
Đối với các hàm giá trị trong R, ta biết rằng tổng, cận trên đúng
và cận dưới đúng của một số hạng hữu hạn của các hàm vững là
hàm vững. Ta hãy xem các tính chất của một tích, một thương và
hợp của các hàm vững.
Mệnh đề 2.2.1.
Giả sử f, g là hai hàm từ X vào R. Nếu f và g là các hàm vững
tại (x
0
, v) và bị chặn trên
Γ (x
0
, v, δ) = {x
0
+ tu : t ∈ (0, δ) , u ∈ B (v, δ)}
với δ > 0 nào đó thì fg vững tại (x
0
, v).
Chứng minh.
Ta có đồng nhất thức
f (x
0
+ tu) g (x
0
+ tu) − f (x
0
+ tv) g (x
0
+ tv)
t
=
f (x
0
+ tu) − f (x
0
+ tv)
t
g (x
0
+ tu)
+f (x
0
+ tv)
g (x
0
+ tu) − g (x
0
+ tv)
t
và sau đó chuyển qua giới hạn khi (t, u) → (0
+
, 0), ta được điều phải
chứng minh.
Chứng minh tương tự ta nhận được mệnh đề sau đây.
22
Số hóa bởi Trung tâm Học liệu />