Sửa đổi Phân loại các bài toán quy hoạch toán học

{{mới}}
Có nhiều cách để phân loại các bài toán [[quy hoạch toán học]]:
* Lồi đối lập với Không lồi
* Trơn đối lập với Không trơn
* Tuyến tính đối lập với Phi tuyến.

==Lồi và không lồi==

===Tập lồi===

Chúng ta nói rằng <math>D \subset \mathbb{R}^n</math> là một tập lồi nếu

:<math>(1 - t)x + ty \in D \text{ với mọi } x \in D, y \in D \text{ và } t \in (0, 1)</math>.

===Bao lồi===

Tập lồi nhỏ nhất chứa <math>\Omega \subset \mathbb{R}^n</math> được gọi là ''bao lồi'' của <math>\Omega</math> và được ký hiệu bởi <math>\text{co}\Omega</math>.

===Hàm lồi===

Một hàm <math>f : R^n \to \mathbb{\bar{R}}</math> được gọi là ''lồi'' nếu ''tập trên đồ thị'' của nó,

{{NumBlk|:|<math>\text{epi}f := \{(x, \alpha) : x \in \mathbb{R}, \alpha \in \mathbb{R}, \alpha \ge f(x)</math>|{{EquationRef|1}}}}

là một tập lồi trong không gian tích <math>\mathbb{R}^n \times \mathbb{R}</math>.

===Hàm chính thường===

Hàm <math>f</math> được gọi là chính thường nếu <math>f(x) < +\infty</math> với ít nhất một phần tử <math>x \in \mathbb{R}^n</math> và <math>f(x) > -\infty</math> với mọi <math>x \in \mathbb{R}^n</math>. Hàm <math>f : \mathbb{R}^n \to \mathbb{\bar{R}}</math> được gọi là ''lõm'' nếu hàm <math>-f</math> được xác định bởi công thức <math>(-f)(x) = -f(x)</math> là lồi.

Theo các quy ước thường dùng (xem Rockafellar (1970), tr. 24),

:<math>\alpha + (+\infty) = (+\infty) + \alpha = +\infty \text{ với } -\infty < \alpha \le +\infty</math>

:<math>\alpha + (-\infty) = (-\infty) + \alpha = -\infty \text{ với } -\infty \le \alpha < +\infty</math>

:<math>\alpha(+\infty) = (+\infty)\alpha = +\infty, \alpha(-\infty) = (-\infty)\alpha = -\infty, \text{ với } 0 < \alpha \le +\infty</math>

:<math>\alpha(+\infty) = (+\infty)\alpha = -\infty, \alpha(-\infty) = (-\infty)\alpha = +\infty, \text{ với } -\infty \le \alpha < 0</math>

:<math>0(+\infty) = (+\infty)0 = 0 = 0(-\infty) = (-\infty)0,\ -(-\infty) = +\infty,\ \text{inf}\empty = +\infty,\ \text{sup}\empty = - \infty</math>.

Các tổ hợp <math>(+\infty) + (-\infty)</math> và <math>(-\infty) + (+\infty)</math> là vô nghĩa và sẽ được tránh sử dụng.

===Bất đẳng thức Jensen===

Từ định nghĩa hàm lồi và công thức ({{EquationNote|1}}) ta suy ra rằng hàm số <math>f : \mathbb{R}^n \to \mathbb{R} \cup\{ +\infty \}</math> là lồi khi và chỉ khi

{{NumBlk|:|<math>f((1 - t)x + ty) \le (1 - t)f(x) + tf(y), \forall x, y \in \mathbb{R}^n, \forall t \in (0, 1).</math>|{{EquationRef|2}}}}

Tổng quát hơn, một hàm <math>f : \mathbb{R}^n \to \mathbb{R} \cup \{+\infty\}</math> là lồi khi và chỉ khi

:<math>f(\lambda _1 x_1 + ... + \lambda _k x_k) \le \lambda _1 f(x_1) + ... + \lambda _k f(x_k) \qquad \text{(bất đẳng thức Jensen)}</math>

với mọi <math>x_1, ... , x_k \in \mathbb{R}^n</math> và <math>\lambda _1 \ge 0, ... , \lambda _kk \ge 0, \lambda _1 + ... + \lambda _k = 1</math>. (Xem Rockafellar (1970), Định lý 4.3.)

===Hàm lồi chặt===

Giả sử rằng <math>D \subset \mathbb{R}^n</math> là một tập lồi. Nếu bất đẳng thức trong ({{EquationNote|2}}) đúng với mọi <math>x, y \in D</math> và với mọi <math>t \in (0, 1)</math>, thì ta nói <math>f</math> là lồi trên <math>D</math>. Nếu bất đẳng thức trong ({{EquationNote|2}}) nghiệm đúng như một bất đẳng thức chặt với mọi <math>x, y \in D</math> mà <math>x \ne y</math> và với mọi <math>t \in (0, 1)</math>, thì ta nói <math>f</math> là lồi chặt trên <math>D</math>.

===Bài toán quy hoạch lồi===

Ta nói rằng <math>(P)</math> là bài toán quy hoạch lồi nếu <math>D</math> là tập lồi và <math>f</math> là hàm lồi. Nếu <math>(P)</math> bài toán quy hoạch lồi, thì

{{NumBlk|:|<math>\text{Sol}(P) = \text{loc}(P).</math>|{{EquationRef|3}}}}

===Bài toán quy hoạch không lồi===

Nếu <math>D</math> là tập không lồi hoặc <math>f</math> là hàm không lồi, thì ta nói <math>(P)</math> là ''bài toán quy hoạch không lồi.''

Xét bài toán

{{NumBlk|:|<math>\text{min}\{f(x) = (x_1 - c_1)^2 + (x_2 - c_2)^2 : x \in D\},</math>|{{EquationRef|4}}}}

ở đó <math>D = {x = (x_1, x_2) \in \mathbb{R}^2 : x_1 \ge 0} \cup {x = (x_1, x_2) : x_2 \ge 0}</math> và <math>c = (c_1, c_2) = (-2, -1)</math>. Nhận xét rằng <math>f</math> là lồi, trong khi <math>D</math> là không lồi. Rõ ràng rằng ({{EquationNote|4}}) là tương đương với bài toán sau đây:

{{NumBlk|:|<math>\text{min}\{\lVert x - c \rVert : x \in D\}.</math>|{{EquationRef|5}}}}

Có thể thấy rằng tập nghiệm của ({{EquationNote|4}}) và ({{EquationNote|5}}) chỉ gồm một điểm <math>(-2, 0)</math>, còn tập nghiệm địa phương gồm hai điểm: <math>(-2, 0)</math> và <math>(0, -1)</math>. Điều này chứng tỏ rằng đẳng thức ({{EquationNote|3}}) nói chung không đúng với các bài toán quy hoạch không lồi.

Đặt <math>f_1(x) = -x + 2, f_2(x) = x + \frac{3}{2}</math> với mọi <math>x \in \mathbb{R}</math>. Định nghĩa <math>f(x) = \text{min}{f_1(x), f_2(x)}</math> và chọn <math>D = [0, 2] \subset R</math>. Với hàm <math>f</math> và tập <math>D</math> đó, chúng ta có

:<math>\text{Sol}(P) =\{2\},\quad \text{loc}(P) = \{0, 2\};</math>

tức là đẳng thức ({{EquationNote|3}}) không đúng cho bài toán trong ví dụ này. Ở đây, <math>f</math> là hàm không lồi và <math>D</math> là tập lồi.

Các hàm lồi có nhiều tính chất thú vị. Chẳng hạn, hàm lồi là liên tục tại mỗi điểm trong của miền hữu hiệu của nó và nó là khả vi theo hướng tại mỗi điểm thuộc miền đó.

===Miền hữu hiệu===

Đối với mỗi hàm số <math>f : \mathbb{R}^n \to \bar{\mathbb{R}}</math>, tập hợp

:<math>\text{dom}f := \{x \in \mathbb{R}^n: -\infty < f(x) < +\infty\}</math>

được gọi là ''miền hữu hiệu'' của f.

===Đạo hàm theo hướng===

Đối với một điểm <math>\bar{x} \in \text{dom}f</math> và một véctơ <math>v \in \mathbb{R}^n</math>, nếu giới hạn

:<math>f'(\bar{x}; v) := \lim_{t \downarrow 0} \frac{f(\bar{x} + tv) - f(\bar{x})}{t}</math>

(có thể nhận các giá trị <math>+\infty</math> và <math>-\infty</math>) tồn tại, thì <math>f</math> được gọi là khả vi theo hướng tại <math>\bar{x}</math> theo hướng <math>v</math> và giá trị <math>f'(\bar{x}; v)</math> được gọi là đạo hàm theo hướng của <math>f</math> tại <math>\bar{x}</math> theo hướng <math>v</math>. Nếu <math>f'(\bar{x}; v)</math> tồn tại với mọi <math>v \in \mathbb{R}^n</math>, thì <math>f</math> được gọi là ''khả vi theo hướng'' tại <math>\bar{x}</math>.

Trong hai định lý tiếp theo đây, <math>f : \mathbb{R}n \to \mathbb{R} \cup \{+\infty\}</math> là một hàm lồi và chính thường.

'''Định lý 0.0.1.''' (Xem Rockafellar (1970), Định lý 10.1) Nếu <math>\bar{x} \in \mathbb{R}^n</math> và <math>\delta > 0</math> là điểm và số thực sao cho hình cầu mở <math>B(\bar{x}, \delta)</math> chứa trong <math>\text{dom}f</math>, thì phần hạn chế của <math>f</math> trên <math>B(\bar{x}, \delta)</math> là một hàm số thực liên tục.

'''Định lý 0.0.2.''' (Xem Rockafellar (1970), Định lý 23.1) Nếu <math>\bar{x} \in \text{dom}f</math>, thì với mỗi <math>v \in \mathbb{R}^n</math> giới hạn

:<math>f'(\bar{x};\ v) := \lim_{t \downarrow 0} \frac{f(\bar{x} + tv) - f(\bar{x})}{t}</math>

tồn tại, và ta có

:<math>f'(\bar{x};\ v) = \inf_{t > 0} \frac{f(\bar{x} + tv) - f(\bar{x})}{t}.</math>

===Nón pháp tuyến===

''Nón pháp tuyến'' <math>N_D(\bar{x})</math> của tập lồi <math>D \subset \mathbb{R}^n</math> tại một điểm <math>\bar{x} \in \mathbb{R}^n</math> được cho bởi công thức

<math>
N_D(\bar{x}) =
\begin{cases}
\{ x^* \in R^n : \langle x^*, x - \bar{x} \rangle \le 0 \text{ với mọi } x \in D\} & \text{nếu } \bar{x} \in D\\
\empty & \text{nếu } \bar{x} \notin D.
\end{cases}
</math>

===Dưới vi phân===

Dưới vi phân <math>\partial f(\bar{x})</math> của một hàm lồi <math>f : R^n \to \mathbb{R}</math> tại một điểm <math>\bar{x} \in \mathbb{R}n</math> được định nghĩa bằng công thức

:<math>\partial f(\bar{x}) = \{ x^* \in \mathbb{R}^n : f(\bar{x}) + \langle x^*, x - \bar{x} \rangle \le f(x) \text{ với mọi } x \in \mathbb{R}^n</math>

===Phần trong tương đối===

Tập con <math>M \subset \mathbb{R}^n</math> được gọi là một ''tập affine'' (a-phin) nếu <math>tx + (1 - t)y \in M</math> với mọi <math>x \in M, y \in M</math> và <math>t \in R</math>. Đối với một <math>D \subset \mathbb{R}^n</math>, ''bao affine'' aff<math>D of D</math> là tập affine nhỏ nhất chứa <math>D</math>. Phần trong tương đối của <math>D</math> được xác định bởi công thức

:<math>\text{ri}D := \{ x \in D : \exists \partial > 0 \text{ sao cho } B(x, \partial) \cap \text{aff}D \subset D \}</math>

Định lý sau đây mô tả mối quan hệ giữa đạo hàm theo hướng và dưới vi phân của các hàm lồi.

'''Định lý 0.0.3.''' (Xem Rockafellar (1970), Định lý 23.4) Cho <math>f</math> là một hàm lồi trên <math>R^n</math>. Nếu <math>x \notin \text{dom}f</math>, thì <math>\partial f(x)</math> là rỗng. Nếu <math>x \in \text{ri}(\text{dom}f)</math>, thì <math>\partial f(x)</math> là khác rỗng và

:<math>f'(x;v) = \sup \{ \langle x^*, v \rangle : x^* \in \partial f(x) \}, \quad \forall v \in R^n. </math>

Ngoài ra, <math>\partial f(x)</math> là tập khác rỗng và giới nội khi và chỉ khi

:<math>x \in \text{int}(\text{dom} f);</math>

trong trường hợp đó, <math>f'(x; v)</math> là hữu hạn với mỗi <math>v \in \mathbb{R}^n</math>.

Kết quả sau đây được gọi là Định lý Moreau-Rockafellar.

'''Định lý 0.0.4.''' (Xem Rockafellar (1970), Định lý 23.8) Cho <math>f = f_1 + ... + f_k</math>, ở đó <math>f_1 + ... + f_k</math> là các hàm lồi, chính thường trên <math>\mathbb{R}^n</math>. Nếu

:<math>\bigcap_{i = 1}^{k}\text{ri}(\text{dom} f_i) \ne \empty,</math>

thì

:<math>\partial f(x) = \partial f_1 (x) + ... + \partial f_k (x), \quad \forall x \in R^n.</math>

Điều kiện cần và đủ tối ưu cho các bài toán quy hoạch lồi được phát biểu như sau.

'''Định lý 0.0.5.''' (Xem Rockafellar (1970), Định lý 27.4) Giả sử rằng <math>f</math> là một hàm lồi, chính thường trên <math>\mathbb{R}^n</math> và <math>D \subset \mathbb{R}^n</math> là một tập lồi. Nếu bao hàm thức

{{NumBlk|:|<math>0 \in \partial f(\bar{x}) + N_D(\bar{x})</math>|{{EquationRef|6}}}}

nghiệm đúng với <math>\bar{x} \in \mathbb{R}^n</math>, thì <math>\bar{x}</math> là nghiệm của <math>(P)</math>. Ngược lại, nếu

{{NumBlk|:|<math>\text{ri}(\text{dom} f) \cap \text{ri}D \ne \empty ,</math>|{{EquationRef|7}}}}

thì ({{EquationNote|6}}) là điều kiện cần và đủ cho <math>\bar{x} \in \mathbb{R}^n</math> là nghiệm của <math>(P)</math>. Nói riêng ra, nếu <math>D = \mathbb{R}^n</math>, thì <math>\bar{x}</math> là một nghiệm của <math>(P)</math> khi và chỉ khi <math>0 \in \partial f(\bar{x})</math>.

Bao hàm thức ({{EquationNote|6}}) có nghĩa là tồn tại <math>x^* \in \partial f (\bar{x})</math> và <math>u^* \in N_D(\bar{x})</math> sao cho 0 = x<sup>∗</sup> + u<sup>∗</sup>. Nhận xét rằng ({{EquationNote|7}}) là một điều kiện chính quy cho bài toán quy hạch lồi có dạng <math>(P)</math>.

Định lý 0.0.5 là công cụ hiệu quả để giải nhiều bài toán quy hoạch lồi. Chúng ta có thể minh họa điều đó bằng cách xét ví dụ sau đây.

===Điểm Fermat===

Cho <math>A, B, C</math> là ba điểm trong không gian hai chiều <math>\mathbb{R}^2</math> với các tọa độ tương ứng là

:<math>a = (a_1, a_2), b = (b_1, b_2), c = (c_1, c_2).</math>

Giả sử rằng không tồn tại đường thẳng nào chứa tất cả ba điểm đó. Bài toán đặt ra là tìm một điểm <math>M</math> trong <math>\mathbb{R}^2</math> với các tọa độ <math>\bar{x} = (\bar{x}_1, \bar{x}_2)</math> sao cho tổng khoảng cách từ <math>M</math> tới <math>A, B</math> và <math>C</math> là tối thiểu. Điều đó có nghĩa rằng <math>\bar{x}</math> là một nghiệm của bài toán quy hoạch lồi không có ràng buộc sau:

{{NumBlk|:|<math>\text{min} \{ f(x) := \lVert x - a \rVert + \lVert x - b \rVert + \lVert x - c \rVert : x \in \mathbb{R}^2 \}.</math>|{{EquationRef|8}}}}

Sử dụng định lý Weierstrass và tính lồi chặt của <math>f(x)</math> trên <math>\mathbb{R}^2</math>, ta có thể chứng tỏ rằng ({{EquationNote|8}}) có duy nhất nghiệm; xem Lee, Tam, Yen (2005), tr. 11–13. Để ý rằng <math>f = f_1 + f_2 + f_3</math>, ở đó <math>f_1(x) = \lVert x - a \rVert, f_2(x) = \lVert x - b \rVert, f_3(x) = \lVert x - c \rVert</math>. Theo Định lý 0.0.5, <math>\bar{x}</math> là nghiệm của ({{EquationNote|8}}) khi và chỉ khi <math>0 \in \partial f(\bar{x})</math>. Vì <math>\text{dom} f_i = \mathbb{R}^2 (i = 1, 2 ,3)</math>, sử dụng Định lý 0.0.4 ta có thể viết bao hàm thức cuối dưới dạng

:<math>0 \in \partial f_1 (\bar{x}) + f_2 (\bar{x}) + f_3 (\bar{x})</math>.

Tiếp theo, bằng cách tính toán các dưới vi phân (xem Lee, Tam, Yen (2005), tr. 10), ta có thể xác định được điểm <math>\bar{x}</math> là nghiệm duy nhất của bao hàm thức này. Trong ngôn ngữ của Hình học Euclide, ta có các kết luận sau:

# Nếu một trong ba góc của tam giác <math>ABC</math>, ví dụ như góc <math>\hat{A}</math>, là lớn hơn hoặc bằng 120&deg;, thì <math>M \equiv A</math> là nghiệm duy nhất của bài toán đang được xét.
# Nếu tất cả các góc của tam giác <math>ABC</math> đều nhỏ hơn 120&deg;, thì nghiệm duy nhất của bài toán là điểm <math>M</math> nhìn các cạnh <math>BC</math>, <math>AC</math> và <math>AB</math> của tam giác <math>ABC</math> dưới cùng một góc 120&deg;(Điểm <math>M</math> đặc biệt này được gọi là điểm Fermat hay điểm Torricelli (xem Weisstein (1999)). Có thể chứng tỏ rằng điểm Fermat thuộc phần trong của tam giác <math>ABC</math>.)

Trong bài toán <math>(P)</math>, nếu <math>D</math> là tập nhiệm của một hệ các phương trình và bất phương trình, thì điều kiện tối ưu bậc nhất có thể viết được ở dạng có sử dụng các nhân tử Lagrange.

Chúng ta hãy xét bài toán <math>(P)</math> dưới các giả thiết <math>f : \mathbb{R}^n \to \mathbb{R}</math> là một hàm lồi,

{{NumBlk|:|<math>D = \{x \in \mathbb{R}^n : g_1(x) \le, ..., g_m(x) \le 0, h_1(x) = 0, ..., h_s(x) = 0\}</math>|{{EquationRef|9}}}}

ở đó <math>g_i: \mathbb{R}^n \to \mathbb{R}, i = 1, . . . , m</math>, là các hàm lồi, <math>h_j: \mathbb{R}^n \to \mathbb{R}, j = 1, . . . , s</math> là các hàm affine, nghĩa là tồn tại <math>a_j \in \mathbb{R}^n</math> và <math>\alpha _j \in \mathbb{R}</math> sao cho <math>h_j(x) = \langle a_j, x \rangle + \alpha _j</math> với mọi <math>x \in \mathbb{R}^n</math>. Chúng ta chấp nhận rằng các ràng buộc bất đẳng thức (tương ứng, các ràng buộc bất đẳng thức) có thể không có mặt trong ({{EquationNote|9}}). Để cho gọn, chúng ta sử dụng cách viết hình thức <math>m = 0</math> (tương ứng, <math>s = 0</math>) để chỉ rằng các ràng buộc bất đẳng thức (tương ứng, các ràng buộc bất đẳng thức) không có mặt trong ({{EquationNote|9}}).

'''Định lý 0.0.6.''' (Định lý Kuhn-Tucker cho các bài toán quy hoạch lồi; xem Rockafellar (1970), tr. 283) Giả sử rằng <math>(P)</math> là bài toán quy hoạch lồi với <math>D</math> được cho bởi ({{EquationNote|9}}). Giả sử rằng các giả thiết đặt lên <math>f, g_i (i = 1, ... , m)</math> và <math>h_j (j = 1, ... s)</math> như đã nói ở trên được thỏa mãn. Giả sử rằng tồn tại véctơ <math>z \in \mathbb{R}^n</math> sao cho

{{NumBlk|:|<math>g_i(z) < 0 \text{ với } i = 1, ..., m \quad \text{và} \quad h_j (z) = 0 \text{ với } j = 1, ..., s</math>|{{EquationRef|10}}}}

Khi đó, <math>\bar{x}</math> là một nghiệm của <math>(P)</math> khi và chỉ khi tồn tại <math>m + s</math> số thực <math>\lambda _1, ... , \lambda _m, \mu _1, ... , \mu _s,</math> được gọi là các nhân tử Langrange tương ứng với <math>\bar{x}</math>, sao cho các điều skiện Kuhn-Tucker sau đây được thỏa mãn:

# <math>\lambda _i \ge 0, g_i(\bar{x}) \le 0 \text{ và } \lambda _i f _i(\bar{x}) = 0 \text{ với } i = 1, ... , m,</math>
# <math>h _j (\bar{x}) = 0 \text{ với } j = 1, ... , s,</math>
# <math>0 \in \partial f(\bar{x}) + \sum_{i = 1}^{m} \lambda _i \partial g _i(\bar{x}) + \sum_{j = 1}^{s} \mu _j a_j.</math>

Nhận xét rằng ({{EquationNote|10}}) là một điều kiện ''chuẩn hóa ràng buộc'' (còn được gọi là ''điều kiện chính quy ràng buộc'', hay đơn giản là ''điều kiện chính quy'') cho các bài toán hoạch lồi. Nếu <math>s = 0</math> thì nó trở thành

:<math> \exists z \in \mathbb{R}^n \text{ s.t. } g_i(z) < 0 \text{ với } i = 1, ..., m \quad \text{(Điều kiện Slater)}</math>

Nếu <math>m = 0</math> thì ({{EquationNote|10}}) tương đương với đòi hỏi rằng <math>D</math> là khác rỗng. Trên thực tế, trong trường hợp đó, có thể bỏ qua điều kiện ({{EquationNote|10}}) trong phát biểu của Định lý 0.0.6.

==Trơn và không trơn==
===Định nghĩa===
Để cho gọn, nếu <math>f : \mathbb{R}^n \to R</math> là hàm khả vi Fréchet liên tục, thì ta nói rằng <math>f</math> là một <math>\text{hàm } C^1</math> (tức là hàm thuộc lớp <math>C^1</math>).

Ta gọi <math>(P)</math> là bài toán quy hoạch toán học trơn nếu <math>f : \mathbb{R}^n \to \mathbb{R}</math> là hàm <math>C^1</math> và <math>D</math> biểu diễn được dưới dạng ({{EquationNote|9}}), ở đó <math>g_i: \mathbb{R}^n \to \mathbb{R} (i = 1, . . . , m)</math> và <math>h_j : \mathbb{R}^n \to \mathbb{R} (j = 1, . . . , s)</math> là các hàm <math>C^1</math>. Ngược lại, <math>(P)</math> được gọi là bài toán quy hoạch toán học không trơn.

===Hàm Lipschitz địa phương===

Một hàm <math>f : \mathbb{R}^n \to \mathbb{R}</math> được gọi là Lipschitz địa phương tại <math>\bar{x} \in \mathbb{R}^n</math> nếu tồn tại hằng số <math>l \ge 0</math> và lân cận <math>U</math> của <math>\bar{x}</math> sao cho

{{NumBlk|:|<math>\lVert f(x') - f(x) \rVert \le l \lVert x' - x \rVert \quad \text{ với mọi } x, x' \in U</math>|{{EquationRef|11}}}}

Nếu <math>f</math> là Lipschitz địa phương tại mỗi điểm thuộc <math>\mathbb{R}^n</math>, thì <math>f</math> được gọi là hàm số Lipschitz địa phương trên <math>\mathbb{R}^n</math>.

Nếu bất đẳng thức trong ({{EquationNote|11}}) đúng với mọi <math>x, x' /in C</math>, ở đó <math>C</math> là một tập con của <math>\mathbb{R}^n</math>, thì ta nói <math>f</math> là Lipschitz trên <math>C</math> với hệ số Lipschitz <math>l</math>.

===Đạo hàm theo hướng suy rộng theo nghĩa Clarke===

Nếu <math>f</math> là Lipschitz địa phương tại <math>\bar{x}</math>, thì đạo hàm theo hướng suy rộng theo nghĩa Clarke của <math>f</math> tại <math>\bar{x}</math> theo hướng <math>v \in \mathbb{R}^n</math> được định nghĩa bằng công thức

<math>f^0(\bar{x}; v) := \lim_{x \to \bar{x}} \sup_{t \downarrow 0} \frac{f(x + tv) - f(x)}{t}</math>

<math>= sup \Bigl\{ \xi \in \mathbb{R} : \forall \text{ các dãy } x_k \to \bar{x} \text{ và } t_k \to 0+ \text{ sao cho } \xi = \lim_{k \to + \infty} \frac{f(x_k + t_kv) - f(x_k)}{t_k} \Bigr\}.</math>

===Dưới vi phân Clarke===

''Dưới vi phân Clarke'' của <math>f</math> tại <math>\bar{x}</math> được cho bởi công thức

<math>\partial f (\bar{x}) := \{ x^* \in \mathbb{R}^n : f^0 (\bar{x}; v) \ge \langle x^*, v \rangle \text{ với mọi } v \in \mathbb{R}^n \}</math>

Định lý 0.0.7. ( Xem Clarke (1983), các Mệnh đề 2.1.2, 2.2.4, 2.2.6 và 2.2.7)Cho f : Rn → R là một hàm số thực. Khi đó, các khẳng định sau đây nghiệm đúng:

(a) Nếu f là Lipschitz địa phương tại x¯ ∈ Rn, thì

f0(¯x; v) = max{hx∗, vi : x∗ ∈ ∂f(¯x)}

với mọi v ∈ Rn.

(b) Nếu f là hàm C1, thì f là hàm số Lipschitz địa phương và ∂f(¯x) = {∇f(¯x)}, f0(¯x; v) = h∇f(¯x), vi với mọi x¯ ∈ Rn và v ∈ Rn.

(c) Nếu f là lồi, thì f là hàm số Lipschitz địa phương và, với mỗi x¯ ∈ Rn, dưới vi phân Clarke ∂f(¯x) trùng với dưới vi phân của f tại x¯theo nghĩa Giải tích lồi, tức là dưới vi phân được định nghĩa bởi . Ngoài ra, f0 (¯x; v) = f0 (¯x; v) với mỗi v ∈ Rn.

Liên quan đến khẳng định (c) ở trên, chúng ta lưu ý rằng, đạo hàm theo hướng f0(¯x; v) tồn tại (xem Định lý 0.0.2).

===Nón tiếp tuyến Clarke===

Cho C ⊂ Rn là tập con khác rỗng. Nón tiếp tuyến Clarke TC(x) của C tại x ∈ C là tập hợp tất cả các véctơ v ∈ Rn thỏa mãn d0 C(x; v) = 0, ở đó d0 C(x; v) ký hiệu đạo hàm theo hướng suy rộng theo nghĩa Clarke của hàm số Lipschitz

dC(z) := inf{ky − zk : y ∈ C}

tại x theo hướng v.

===Nón pháp tuyến Clarke===

Nón pháp tuyến Clarke NC(x) của C tại x được định nghĩa là nón đối ngẫu của TC(x), tức là

NC(x) = {x∗ ∈ Rn: hx∗, vi ≤ 0 với mọi v ∈ TC(x)}.

Định lý 0.0.8. (Xem Clarke (1983), các Mệnh đề 2.4.3, 2.4.4 và 2.4.5) Với mỗi tập con khác rỗng C ⊂ Rn và với mỗi điểm x ∈ C, các khẳng định sau nghiệm đúng:

(a) NC(x) = n∪t≥0 t∂dC(x)o.

(b) Nếu C là lồi, thì NC(x) trùng với nón pháp tuyến của C tại x được định nghĩa bởi công thức , và TC(x) trùng với bao đóng tôpô của hình nón

cone(C − x) := {tz : t ≥ 0, z ∈ C − x}.

(c) Bao hàm thức v ∈ TC(x) nghiệm đúng khi và chỉ khi với mỗi dãy điểm xk trong C hội tụ đến x và dãy số tk trong (0, +∞) hội tụ đến 0, tồn tại dãy véctơ vk trong Rn hội tụ đến v sao cho xk + tkvk ∈ C với mọi k.

Xét bài toán (P) dưới các giả thiết f : Rn → R là hàm Lipschitz địa phương,

D = {x ∈ C : g1(x) ≤ 0, . . . , gm(x) ≤ 0, h1(x) = 0, . . . , hs(x) = 0}, (12)

với C ⊂ Rn là tập khác rỗng, gi: Rn → R (i = 1, . . . , m) và hj: Rn → R(j = 1, . . . , s) là các hàm Lipschitz địa phương.

Định lý 0.0.9. (Xem Clarke (1983), Định lý 6.1.1 và Nhận xét 6.1.2) Nếu x¯ là nghiệm địa phương của (P), thì tồn tại m + s + 1 số thực λ0 ≥ 0, λ1 ≥ 0, . . . , λm ≥ 0, µ1, . . . , µs, không đồng thời bằng 0, sao cho

0 ∈ λ0∂f(¯x) +Xm i=1 λi∂gi(¯x) +Xs j=1 µj∂hj (¯x) + NC(¯x) (13)

và

λigi(¯x) = 0 với mọi i = 1, 2, . . . , m. (14)

===Điều kiện cần tối ưu bậc nhất ở dạng Fritz-John===

Định lý trên phát biểu điều kiện cần tối ưu bậc nhất ở dạng Fritz-John cho một lớp bài toán không trơn. Dưới những điều kiện chính quy ràng buộc thích hợp, nhân tử λ0 tương ứng với hàm mục tiêu f là dương. Trong trường hợp đó, bằng cách chia cả hai vế của bao hàm thức trong (13) và các đẳng thức trong (14) cho λ0, và đặt λei = λi/λ0 cho mỗi i = 1, . . . , m, µej = µj/λ0 cho mỗi j = 1, . . . , s, chúng ta thu được

0 ∈ ∂f(¯x) +Xmi=1λei∂gi(¯x) +X s j=1 µej∂hj (¯x) + NC(¯x) (15)

và

λeigi(¯x) = 0 với mọi i = 1, 2, . . . , m. (16)

Tương tự như trong trường hợp các bài toán quy hoạch lồi (Định lý 1.6), nếu x¯ ∈ D và (15), (16) nhiệm đúng, thì các số λe1 ≥ 0, . . . , λem ≥ 0, µe1 ∈ R, . . . , µes ∈ R được gọi là các nhân tử Lagrange tương ứng với x¯.

Hai quy tắc nhân tử Lagrange sau đây suy ra từ Định lý 0.0.9 (xem Clarke(1983), tr. 234–236).

Hệ quả 0.0.1. Nếu x¯ là một nghiệm địa phương của (P) và nếu điều kiện chuẩn hóa ràng buộc


h0 ∈Xmi=1λi∂gi(¯x) +Xsj=1µj∂hj (¯x) + NC(¯x),

λ1 ≥ 0, . . . , λm ≥ 0, µ1 ∈ R, . . . , µs ∈ R;

λigi(¯x) = 0 với i = 1, . . . , mi

=⇒[λ1 = · · · = λm = 0, µ1 = · · · = µs = 0i]

được thỏa mãn, thì tồn tại các nhân tử Lagrange λ1 ≥ 0, . . . , λm ≥ 0, µ1 ∈ R, . . . , µs ∈ R sao cho λigi(¯x) = 0 với i = 1, 2, . . . , m, và

0 ∈ ∂f(¯x) +X m i=1 λi∂gi(¯x) +X s j=1 µj∂hj (¯x) + NC(¯x).

Hệ quả 0.0.2. Giả sử rằng x¯ là một nghiệm địa phương của bài toán quy hoạch trơn (P), ở đó D được cho bởi công thức (9). Nếu điều kiện chuẩn hóa ràng buộc Mangasarian-Fromovitz (viết tắt: (MFCQ))

Các véctơ {∇hj (¯x) : j = 1, . . . , s} là độc lập tuyến tính,

và tồn tại v ∈ Rn sao cho h∇hj (¯x), vi = 0

với mỗi j = 1, . . . , s, và h∇gi(¯x), vi < 0

với mỗi i = 1, . . . , m thỏa mãn gi(¯x) = 0

được thỏa mãn, thì tồn tại các nhân tử Lagrange λ1 ≥ 0, . . . , λm ≥ 0, µ1 ∈ R, . . . , µs ∈ R sao cho λigi(¯x) = 0 với i = 1, 2, . . . , m, và

0 = ∇f(¯x) +X m i=1 λi∇gi(¯x) +X s j=1 µj∇hj (¯x).

Sử dụng Định lý 0.0.9 ta dễ dàng chứng minh được quy tắc nhân tử Lagrange cho các bài toán quy hoạch lồi đã được phát biểu trong Định lý 0.0.6 (xem Lee, Tam, Yen (2005), tr. 18–19).

==Tuyến tính và phi tuyến==

===Tập lồi đa diện===

Tập hợp D ⊂ Rn được gọi là tập lồi đa diện nếu như ta có thể biểu diễn D dưới dạng giao của một số hữu hạn các nửa không gian đóng của R n; nghĩa là tồn tại các véctơ khác không a1, . . . , am ∈ Rn và các số thực β1, . . . , βm sao cho

D = {x ∈ Rn : hai , xi ≥ βi với i = 1, . . . , m}. (17)

Nói một cách khác, D là tập nghiệm của một hệ gồm hữu hạn các bất đẳng thức tuyến tính. (Chúng ta quy ước rằng giao của một họ rỗng của các nửa không gian đóng của Rn là Rn. Vì thế, D = Rn cũng là một tập lồi đa diện.)

===Điểm cực biên===

Một điểm x ∈ D được gọi là điểm cực biên của D nếu như không thể nào biểu diễn x dưới dạng x = (1 − t)y + tz, ở đó y ∈ D, z ∈ D, y 6= z, và t ∈ (0, 1). Tập hợp tất cả các điểm cực biên của D được ký hiệu bởi extrD. Ký hiệu bởi A ma trận cấp m × n với các phần tử aij (i = 1, . . . , m, j = 1, . . . , n), ở đó aij là thành phần thứ j của véctơ ai. Đặt b = (β1, . . . , βm) ∈ Rm. Khi đó, ta có thể viết lại (17) như sau:

D = {x ∈ Rn: Ax ≥ b}.

Từ đây về sau, đối với hai véctơ tùy ý y = (y1, . . . , ym) ∈ Rm và z = (z1, . . . , zm) ∈ R m, ta viết y ≥ z nếu yi ≥ zi với i = 1, . . . , m. Ta sẽ viết y > z nếu yi > zi với mọi i = 1, . . . , m. Vì

{x ∈ Rn: Ax = b} = {x ∈ Rn: Ax ≥ b, (−A)x ≥ −b},

ta suy ra rằng {x ∈ Rn: Ax = b} là một tập lồi đa diện.

===Bài toán quy hoạch tuyến tính===

Bài toán (P) được gọi là bài toán quy hoạch tuyến tính nếu D là tập lồi đa diện và f(x) là phiếm hàm tuyến tính. Ngược lại, (P) được gọi là bài toán quy hoạch phi tuyến.

Có ba dạng điển hình của bài toán quy hoạch tuyến:

min{f(x) = hc, xi : x ∈ Rn, Ax ≥ b},

min{f(x) = hc, xi : x ∈ Rn, Ax = b, x ≥ 0},

min{f(x) = hc, xi : x ∈ Rn, Ax ≥ b, Cx = d};

chúng được gọi tương ứng là dạng chuẩn, dạng chính tắc, và dạng tổng quát. Ở
đây, A ∈ Rm×n, C ∈ Rs×n là các ma trận cho trước, c ∈ Rn, b ∈ Rm và d ∈ Rs là các véctơ cho trước.

Xét bài toán quy hoạch tuyến tính dạng chuẩn:

min nx1 +12x2 : x = (x1, x2), x1 + x2 ≥ 1, x1 ≥ 0, x2 ≥ 0o.

Dễ kiểm tra rằng Sol(P) = {(0, 1)}. Nhận xét rằng tập ràng buộc

D = {x ∈ R2: x1 + x2 ≥ 1, x1 ≥ 0, x2 ≥ 0}

có hai điểm cực biên, cụ thể là extr = {(1, 0), (0, 1)}. Một trong hai điểm đó là nghiệm của bài toán được xét.

===Bài toán đối ngẫu===

Bài toán đối ngẫu của các bài toán quy hoạch tuyến tính dạng chuẩn, dạng chính tắc, và dạng tổng quát, tương ứng là:

max{hb, yi : y ∈ Rm, ATy = c, y ≥ 0},

max{hb, yi : y ∈ Rm, ATy ≤ c},

max{hb, yi + hd, zi : (y, z) ∈ Rm × Rs, ATy + CTz = c, y ≥ 0}.

Vì bài toán quy hoạch tuyến tính cũng là một bài toán quy hoạch lồi, nên nó có tất cả các tính chất của bài toán quy hạch lồi. Ngoài ra, bài toán quy hoạch tuyến tính còn có những tính chất đặc biệt khác.

Định lý 0.0.10. (Xem Dantzig (1963)) Cho (P) là bài toán quy hoạch tuyến tính ở một trong các dạng điển hình. Các tính chất sau đây nghiệm đúng:

(i) Nếu tập ràng buộc là khác rỗng và nếu v(P) > −∞, thì Sol(P) là tập lồi đa diện khác rỗng.

(ii) Nếu cả hai tập extrD và Sol(P) đều khác rỗng, thì giao extrD ∩ Sol(P) cũng là khác rỗng.

(iii) Nếu rankA = n và tập D := {x ∈ Rn: Ax = b, x ≥ 0} là khác rỗng, thì D phải có ít nhất một điểm cực biên.

(iv) Giá trị tối ưu v(P) của (P) và giá trị tối ưu v(P0) of của bài toán đối ngẫu(P0) của (P) là bằng nhau, nếu như tập ràng buộc của ít nhất là một trong hai bài toán là khác rỗng.

==Tham khảo==
*Lee, Tam, Yen (2005), Chương 1.