Tổng quan về luận án

Lý thuyết giải tích biến phân (Variational Analysis) hiện đại đóng vai trò là nền tảng toán học cốt lõi cho lý thuyết tối ưu, điều khiển tối ưu và giải tích phi tuyến. Được khởi xướng từ đầu những năm 1960 với các nghiên cứu của Jean Jacques Moreau về dưới vi phân cho hàm lồi, tiếp nối bởi bước đột phá của Boris S. Mordukhovich vào giữa thập niên 1970 khi mở rộng dưới vi phân cho hàm không lồi, và được hệ thống hóa toàn diện trong chuyên khảo kinh điển của R. Tyrrell Rockafellar và Roger J-B. Wets (1998), giải tích biến phân đã phát triển qua hai giai đoạn lớn: giải tích biến phân bậc nhất (first-order variational analysis) đã đạt độ hoàn thiện cao, trong khi giải tích biến phân bậc hai (second-order variational analysis) đang là mũi nhọn nghiên cứu đương đại nhằm giải quyết các hệ biến phân phi tuyến phức tạp.

Trong bức tranh khoa học đó, luận án tiến sĩ toán học chuyên ngành Toán giải tích / Lý thuyết tối ưu của tác giả Lê Văn Hiển mang tựa đề "Một số kết quả về tính dưới chính quy mêtric trong giải tích biến phân và ứng dụng" (thực hiện tại Trường Đại học Vinh) đã giải quyết những khoảng trống lý thuyết trọng yếu về điều kiện chuẩn hóa ràng buộc yếu nhất cho các cấu trúc vi phân suy rộng bậc hai.

                  ┌────────────────────────────────────────────────────────┐
                  │              GIẢI TÍCH BIẾN PHÂN BẬC HAI                │
                  └───────────────────────────┬────────────────────────────┘
                                              │
                     ┌────────────────────────┴────────────────────────┐
                     ▼                                                 ▼
     ┌───────────────────────────────┐                 ┌───────────────────────────────┐
     │ ĐẠO HÀM ĐỒ THỊ NÓN PHÁP TUYẾN │                 │     TÍNH ỔN ĐỊNH XIÊN         │
     │      D N_Γ (x̄ | x̄*)           │                 │       (TILT STABILITY)        │
     └───────────────┬───────────────┘                 └───────────────┬───────────────┘
                     │                                                 │
                     │  Loại bỏ BEPP / CRCQ                            │  Đặc trưng qua đạo hàm
                     │  Chỉ cần MSCQ thuần túy                         │  đồ thị dưới gradient
                     ▼                                                 ▼
     ┌───────────────────────────────┐                 ┌───────────────────────────────┐
     │ ĐỊNH LÝ 2.10 (HỢP NHẤT)       │                 │ CẢI THIỆN ĐIỀU KIỆN ĐỦ SSOSC  │
     │ D N_Γ = ∇²(λᵀq)v + N_K(v)     │                 │ DƯỚI MSCQ (THAY VÌ MFCQ+CRCQ) │
     └───────────────┬───────────────┘                 └───────────────┬───────────────┘
                     │                                                 │
                     └────────────────────────┬────────────────────────┘
                                              ▼
                             ┌───────────────────────────────────┐
                             │ CÔNG TRÌNH CÔNG BỐ QUỐC TẾ:       │
                             │ 02 BÀI BÁO SIAM J. OPTIMIZATION   │
                             └───────────────────────────────────┘

Khoảng trống nghiên cứu (Research Gap) trung tâm xuất phát từ việc tính toán đạo hàm đồ thị (graphical derivative - được Jean-Pierre Aubin đề xuất năm 1981) của ánh xạ nón pháp tuyến $N_\Gamma$ và kiểm soát tính ổn định xiên (tilt stability - khái niệm do René Poliquin và R. T. Rockafellar giới thiệu năm 1998) trước đây luôn đòi hỏi các giả thiết chính quy rất mạnh. Cụ thể, các công trình quốc tế kinh điển của Henrion, Outrata và Surowiez (2013) hay Helmut Gfrerer và Jiri Outrata (2014), Gfrerer và Mordukhovich (2015) đều phải dựa trên giả thiết ánh xạ ràng buộc thỏa mãn tính chính quy mêtric (metric regularity), chuẩn hóa ràng buộc Mangasarian-Fromovitz (MFCQ), chuẩn hóa hạng hằng (CRCQ), hoặc phải bổ sung tính chất điểm cực biên bị chặn (Bounded Extremal Point Property - BEPP). Câu hỏi khoa học nền tảng được đặt ra: Liệu có thể loại bỏ hoàn toàn các giả thiết bổ sung khắt khe này để tính đạo hàm đồ thị của ánh xạ nón pháp tuyến và thiết lập tiêu chuẩn ổn định xiên chỉ dưới duy nhất điều kiện chuẩn hóa dưới chính quy mêtric (Metric Subregularity Constraint Qualification - MSCQ)?

Luận án tập trung giải quyết 2 câu hỏi nghiên cứu và 2 giả thuyết khoa học cụ thể:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để thiết lập công thức hiển biểu diễn đồ thị của đạo hàm đồ thị ánh xạ nón pháp tuyến $DN_\Gamma(\bar{x}|\bar{x}^*)$ khi tập ràng buộc $\Gamma = {x \in \mathbb{R}^n \mid q(x) \in \Theta}$ với $\Theta$ lồi đa diện chỉ thỏa mãn điều kiện MSCQ?
  • Giả thuyết khoa học 1 (H1): Với $\Theta$ là tập lồi đa diện, cấu trúc nón tiếp tuyến $T_{\text{gph} N_\Gamma}(\bar{x}, \bar{x}^*)$ hoàn toàn xác định được thông qua tập nhân tử tối ưu Lagrange $\Lambda(v)$ của bài toán quy hoạch tuyến tính $LP(v)$ và nón pháp tuyến của nón tới hạn $N_K(v)$ mà không cần đến tính chất BEPP hay tính chính quy mêtric đều quanh điểm khảo sát.
  • Câu hỏi nghiên cứu 2 (RQ2): Có thể sử dụng đạo hàm đồ thị của ánh xạ dưới vi phân $D(\partial f)(\bar{x}|\bar{v})$ thay thế cho dưới vi phân qua giới hạn bậc hai $\partial^2 f(\bar{x}|\bar{v})$ để đưa ra đặc trưng bậc hai đầy đủ cho cực tiểu địa phương ổn định xiên trong bài toán tối ưu không ràng buộc và quy hoạch phi tuyến dưới điều kiện MSCQ hay không?
  • Giả thuyết khoa học 2 (H2): Đạo hàm đồ thị dưới gradient cho phép thiết lập các điều kiện đủ bậc hai mạnh (SSOSC) cho tính ổn định xiên trong quy hoạch phi tuyến chỉ với MSCQ, vượt qua ranh giới của các giả thiết MFCQ và CRCQ đồng thời.

Khung lý thuyết của luận án được xây dựng vững chắc trên nền tảng của Lý thuyết Vi phân suy rộng Mordukhovich, Lý thuyết Tính chính quy mêtric và dưới chính quy mêtric của Ioffe-Dontchev-Rockafellar, cùng Khung hàm chính quy gần kề (prox-regularity) của Poliquin-Rockafellar. Đóng góp đột phá của luận án mang tính lượng hóa rõ nét: hợp nhất các kết quả độc lập trước đó của Dontchev-Rockafellar (1996), Henrion et al. (2013), Gfrerer-Mordukhovich (2015) thành một định lý tổng quát duy nhất; chứng minh trực tiếp điều kiện đủ bậc hai mạnh SSOSC bảo toàn tính ổn định xiên dưới MSCQ; và thành quả nghiên cứu đã được công bố tại 02 bài báo khoa học trên tạp chí chuyên ngành xếp hạng hàng đầu thế giới trong lĩnh vực tối ưu: SIAM Journal on Optimization (SIOPT).

Phạm vi nghiên cứu được xác lập chuẩn xác trong không gian Euclid hữu hạn chiều $\mathbb{R}^n$, khảo sát các ánh xạ đa trị $F: \mathbb{R}^n \rightrightarrows \mathbb{R}^m$, tập ràng buộc phi tuyến $\Gamma$ với $\Theta$ lồi đa diện, các hàm mục tiêu $f: \mathbb{R}^n \to \overline{\mathbb{R}}$ chính quy gần kề và liên tục dưới vi phân, cùng lớp bài toán quy hoạch toàn phương có một ràng buộc toàn phương (QCQP).


Literature Review và Positioning

Lịch sử phát triển của giải tích biến phân bậc hai và lý thuyết ổn định tối ưu ghi nhận ba dòng chảy học thuật chủ lưu tương tác chặt chẽ với nhau:

DÒNG CHẢY 1: Lý thuyết Chính quy Mêtric & Chuẩn hóa ràng buộc
Lyusternik (1934) ──► Borwein (1986) ──► Ioffe (1979, 2008) ──► Dontchev & Rockafellar (2004)
[Định lý Ánh xạ mở]   [Metric Regularity] [MSCQ / Dưới chính quy]  [Hệ thống hóa Khái niệm]

DÒNG CHẢY 2: Đạo hàm đồ thị của Ánh xạ Nón pháp tuyến
Aubin (1981) ──► Dontchev-Rockafellar (1996) ──► Henrion et al. (2013) ──► Gfrerer-Mordukhovich (2015)
[Contingent Der.]  [Tập lồi đa diện Γ]           [Giả thiết Metric Reg.]     [MSCQ + giả thiết BEPP]
                                                                                   │
                                                         LUẬN ÁN LÊ VĂN HIỂN ◄─────┘
                                                         [Bỏ hoàn toàn BEPP, chỉ cần MSCQ]

DÒNG CHẢY 3: Lý thuyết Ổn định xiên (Tilt Stability)
Poliquin-Rockafellar (1998) ──► Mordukhovich-Rockafellar (2012) ──► Mordukhovich-Nghia (2015)
[Đặc trưng qua ∂²f]            [Đặc trưng SSOSC dưới LICQ]         [Giới thiệu USOSC dưới MFCQ+CRCQ]
                                                                                   │
                                                         LUẬN ÁN LÊ VĂN HIỂN ◄─────┘
                                                         [Đặc trưng qua D(∂f), SSOSC dưới MSCQ]

Dòng chảy thứ nhất tập trung vào lý thuyết chính quy mêtric và chuẩn hóa ràng buộc. Khởi đầu từ Định lý Lyusternik-Graves và định lý ánh xạ mở Banach-Schauder, Jonathan Borwein (1986) cùng Borwein & Zhuang (1988) đã định hình khái niệm chính quy mêtric gắn liền với tính chất mở tuyến tính và tính chất Aubin của ánh xạ ngược. Alexander D. Ioffe (1979) đề xuất khái niệm điểm chính quy, đặt nền móng cho khái niệm dưới chính quy mêtric (metric subregularity), tương đương với tính chất tĩnh lặng (calmness) của ánh xạ ngược, trước khi Asen L. Dontchev và R. T. Rockafellar (2004) chuẩn hóa thuật ngữ này. Ioffe và Outrata (2008) chỉ ra rằng điều kiện chuẩn hóa trong hệ vi phân suy rộng đối ngẫu thực chất tương đương với tính chính quy mêtric của ánh xạ thích hợp.

Dòng chảy thứ hai liên quan trực tiếp đến việc tính toán đạo hàm đồ thị của ánh xạ nón pháp tuyến $N_\Gamma$. Năm 1996, Dontchev và Rockafellar thiết lập công thức mô tả đồ thị $DN_\Gamma$ cho tập $\Gamma$ lồi đa diện. Đến năm 2013, René Henrion, Jiri Outrata và Thomas Surowiez tính được $DN_\Gamma$ khi ánh xạ $M_q(x) := q(x) - \Theta$ là chính quy mêtric quanh điểm khảo sát; công thức được đơn giản hóa nếu thỏa mãn cả MFCQ và chuẩn hóa ràng buộc hạng hằng (CRCQ). Năm 2014, Gfrerer và Outrata đã có bước tiến lớn khi đưa ra lược đồ chứng minh trực tiếp đạo hàm đồ thị khi thay chính quy mêtric bằng tính dưới chính quy mêtric kết hợp với tính chính quy mêtric đều lân cận. Năm 2015, Helmut Gfrerer và Boris S. Mordukhovich nới lỏng tiếp điều kiện chính quy mêtric đều thành tính chất điểm cực biên bị chặn (BEPP).

Dòng chảy thứ ba là lý thuyết ổn định xiên. Poliquin và Rockafellar (1998) thiết lập mối liên hệ giữa ổn định xiên với tính xác định dương của dưới vi phân bậc hai qua giới hạn $\partial^2 f$. Mordukhovich và Rockafellar (2012) chứng minh rằng dưới chuẩn hóa độc lập tuyến tính (LICQ), điều kiện đủ bậc hai mạnh (SSOSC) là điều kiện cần và đủ cho tính ổn định xiên. Mordukhovich và Outrata (2012) nới lỏng điều kiện sang MFCQ kết hợp CRCQ thì SSOSC là điều kiện đủ. Đến năm 2015, Mordukhovich và Nguyễn Mậu Nam Nghĩa chỉ ra SSOSC không phải điều kiện cần dưới MFCQ + CRCQ và giới thiệu điều kiện đủ bậc hai đều (USOSC).

Về mặt học thuật, tồn tại cuộc tranh luận sâu sắc: Liệu có thể đạt được công thức giải tích chuẩn tắc cho $DN_\Gamma$ và kiểm soát tính ổn định xiên mà không cần viện dẫn tới tính chất hình học phức tạp như BEPP hay các điều kiện hạng hằng ngặt nghèo như CRCQ?

Vị trí của luận án trong y văn quốc tế được xác lập thông qua việc so sánh đối sánh trực tiếp với các công trình chuẩn mực:

  1. So với Henrion, Outrata & Surowiez (2013, SIAM J. Optim.): Nghiên cứu của Henrion et al. đòi hỏi $M_q$ phải chính quy mêtric trên toàn lân cận hoặc thỏa mãn đồng thời MFCQ và CRCQ. Luận án của Lê Văn Hiển đã hạ bậc giả thiết xuống mức chỉ cần $M_q$ dưới chính quy mêtric (MSCQ) tại duy nhất điểm khảo sát $\bar{x}$.
  2. So với Gfrerer & Mordukhovich (2015, SIAM J. Optim.): Công trình của Gfrerer và Mordukhovich yêu cầu đồng thời MSCQ và giả thiết hình học phụ trợ BEPP. Luận án chứng minh rằng khi $\Theta$ là tập lồi đa diện, giả thiết BEPP hoàn toàn dư thừa và có thể triệt tiêu.
  3. So với Mordukhovich & Outrata (2012, Math. Program.)Mordukhovich & Nghia (2015, Math. Program.): Luận án đã mở rộng phạm vi áp dụng của điều kiện đủ bậc hai mạnh SSOSC cho tính ổn định xiên từ không gian thỏa mãn MFCQ + CRCQ sang không gian chỉ cần thỏa mãn MSCQ – một bước tiến nhảy vọt về tính bao quát của lý thuyết cực tiểu địa phương ổn định.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án tạo ra bước phát triển mới cho hệ thống lý thuyết giải tích biến phân bậc hai của Boris S. Mordukhovich và R. T. Rockafellar qua ba đóng góp cốt lõi:

Thứ nhất, hoàn thiện biểu diễn giải tích nón tiếp tuyến đồ thị nón pháp tuyến $T_{\text{gph} N_\Gamma}(\bar{x}, \bar{x}^)$. Luận án chứng minh công thức đóng: $$T_{\text{gph} N_\Gamma}(\bar{x}, \bar{x}^) = \left{ (v, v^) \in \mathbb{R}^n \times \mathbb{R}^n ;\middle|; \exists \lambda \in \Lambda(v) : v^ \in \nabla^2(\lambda^T q)(\bar{x})v + N_K(v) \right}$$ qua đó xác lập công thức tính đạo hàm đồ thị: $$DN_\Gamma(\bar{x}|\bar{x}^)(v) = \nabla^2(\lambda^T q)(\bar{x})v \Big|{\lambda \in \Lambda(v)} + N_K(v)$$ trong đó $\Lambda(v)$ là tập nghiệm tối ưu của bài toán quy hoạch tuyến tính $LP(v): \min { -v^T \nabla^2(\lambda^T q)(\bar{x})v \mid \lambda \in \Lambda }$, $K = T\Gamma(\bar{x}) \cap {\bar{x}^}^\perp$ là nón tới hạn, và $N_K(v)$ được biểu diễn hiển: $$N_K(v) = \left{ \nabla q(\bar{x})^T \mu ;\middle|; \mu^T \nabla q(\bar{x})v = 0, ; \mu \in T_{N_\Theta(\bar{y})}(\lambda) \right} = \sum_{i \in I_q(\bar{x}) \cap v^\perp} t_i \nabla q_i(\bar{x}) - t_0 \bar{x}^* \quad (t_0, t_i \ge 0)$$

Thứ hai, thay đổi căn bản cách tiếp cận nghiên cứu tính ổn định xiên (Tilt Stability Paradigm Shift). Thay vì dựa vào cấu trúc dưới vi phân qua giới hạn bậc hai $\partial^2 f(\bar{x}|\bar{v})$ vốn cực kỳ khó tính toán tường minh khi có ràng buộc phi tuyến, luận án tiên phong sử dụng đạo hàm đồ thị dưới gradient $D(\partial f)(\bar{x}|\bar{v})$. Luận án chứng minh rằng với hàm $f$ chính quy gần kề và liên tục dưới vi phân, tính ổn định xiên tương đương với tính xác định dương của toán tử đạo hàm đồ thị dưới gradient: $$\langle v^, v \rangle > 0, \quad \forall v \neq 0, ; v^ \in D(\partial f)(\bar{x}|0)(v)$$

Thứ ba, thiết lập hệ điều kiện đủ bậc hai tối ưu cho quy hoạch phi tuyến dưới chuẩn hóa yếu MSCQ. Luận án chứng minh mệnh đề lý thuyết nền tảng: Nếu $\bar{x}$ là điểm dừng thỏa mãn MSCQ và thỏa mãn điều kiện đủ bậc hai mạnh SSOSC: $$v^T \nabla_{xx}^2 L(\bar{x}, \lambda) v > 0, \quad \forall v \in K(\bar{x}, \bar{x}^) \setminus {0}, ; \forall \lambda \in \Lambda(\bar{x}, \bar{x}^)$$ thì $\bar{x}$ chắc chắn là điểm cực tiểu địa phương ổn định xiên.

                    ┌──────────────────────────────────────────────────┐
                    │ TẬP HỢP CÁC CHUẨN HÓA RÀNG BUỘC (CQ TAXONOMY)    │
                    └────────────────────────┬─────────────────────────┘
                                             │
                       ┌─────────────────────┴─────────────────────┐
                       ▼                                           ▼
             ┌───────────────────┐                       ┌───────────────────┐
             │       LICQ        │                       │       CRCQ        │
             └─────────┬─────────┘                       └─────────┬─────────┘
                       │                                           │
                       ▼                                           ▼
             ┌───────────────────┐                       ┌───────────────────┐
             │       MFCQ        │                       │       CPLD        │
             └─────────┬─────────┘                       └─────────┬─────────┘
                       │                                           │
                       └─────────────────────┬─────────────────────┘
                                             │
                                             ▼
                             ┌───────────────────────────────┐
                             │       MSCQ (DƯỚI CHÍNH QUY)   │  ◄── Điểm tựa của luận án
                             └───────────────┬───────────────┘
                                             │
                                             ▼
                             ┌───────────────────────────────┐
                             │ TẬP ĐIỂM CHẤP NHẬN ĐƯỢC Γ     │
                             └───────────────────────────────┘

Khung phân tích độc đáo

Khung phân tích của luận án là sự tích hợp hữu cơ giữa ba trụ cột lý thuyết lớn:

  1. Lý thuyết Điểm tựa Mêtric (Metric Regularity Theory) kế thừa từ Alexander Ioffe, Asen Dontchev và R. T. Rockafellar;
  2. Lý thuyết Giải tích Đa trị và Đối đạo hàm của Boris S. Mordukhovich;
  3. Hình học Biến phân và Hàm chính quy gần kề của René Poliquin và R. T. Rockafellar.

Đột phá phương pháp luận phân tích thể hiện ở việc kết hợp kỹ thuật xấp xỉ biến phân bậc hai dọc theo hướng dịch chuyển $x_t = \bar{x} + tv + \frac{1}{2}t^2 z$, khai thác tính đối ngẫu mạnh giữa quy hoạch tuyến tính tham số $LP(v)$ và $DP(v)$, cùng với việc ứng dụng Bổ đề Chặn sai số Hoffman (Hoffman's Error Bound Lemma) trên các nón lồi đa diện $\widetilde{\Theta} := \text{pos}{b_i \mid i \in I_0}$. Kỹ thuật này cho phép kiểm soát sai số $d_{\Psi_{I_0}(\bar{x}^)}(\lambda^k) \le \beta t_k (|v_k^| + c_1 c_2 |v_k|)$ mà không hề cần đến giả thiết chính quy mêtric lân cận hay tính bị chặn cực biên BEPP.

Điều kiện biên (Boundary Conditions) của khung phân tích được xác định rõ: Tập $\Theta \subset \mathbb{R}^m$ bắt buộc phải là tập lồi đa diện biểu diễn dưới dạng hệ bất đẳng thức tuyến tính hữu hạn $\Theta = {y \in \mathbb{R}^m \mid \langle b_i, y \rangle \le \alpha_i, i = 1, \dots, \ell}$, hàm ánh xạ ràng buộc $q: \mathbb{R}^n \to \mathbb{R}^m$ và hàm mục tiêu $f: \mathbb{R}^n \to \overline{\mathbb{R}}$ phải thuộc lớp khả vi liên tục hai lần ($C^2$).


Phương pháp nghiên cứu tiên tiến

Thiết kế nghiên cứu

Nghiên cứu tuân thủ lập trường triết học nhận thức luận Diễn dịch Toán học Tiên đề (Deductive Mathematical Rationalism / Formal Positivism) của Toán học thuần túy. Thiết kế nghiên cứu cấu trúc đa tầng (Multi-level Analytical Design):

  • Tầng 1 (Cấu trúc hình học bậc một): Nghiên cứu nón tiếp tuyến Bouligand-Severi $T_\Omega(\bar{x})$, nón pháp tuyến Fréchet $\widehat{N}\Omega(\bar{x})$, nón pháp tuyến Mordukhovich $N\Omega(\bar{x})$ và dưới vi phân Fréchet $\widehat{\partial}\phi(\bar{x})$.
  • Tầng 2 (Toán tử biến phân bậc hai): Đạo hàm đồ thị $DF(\bar{x}|\bar{y})$, đối đạo hàm Fréchet $\widehat{D}^* F(\bar{x}|\bar{y})$, đạo hàm đồ thị của ánh xạ dưới vi phân $D(\partial f)(\bar{x}|\bar{v})$ và đạo hàm đồ thị của ánh xạ nón pháp tuyến $DN_\Gamma(\bar{x}|\bar{x}^*)$.
  • Tầng 3 (Hệ thống tham số và tính ổn định): Phương trình suy rộng chứa tham số $0 \in \Phi(p, x) + N_\Gamma(x)$, tính chất tĩnh lặng cô lập (isolated calmness) của ánh xạ nghiệm $S(p)$, và tính ổn định xiên $\text{tilt}(f, \bar{x})$.

Quy trình nghiên cứu rigorous

Quy trình chứng minh giải tích biến phân được thực hiện qua các giao thức toán học nghiêm ngặt:

                  ┌────────────────────────────────────────────────────────┐
                  │ QUY TRÌNH CHỨNG MINH ĐỊNH LÝ ĐẠO HÀM ĐỒ THỊ 2.10       │
                  └───────────────────────────┬────────────────────────────┘
                                              │
                      ┌───────────────────────┴───────────────────────┐
                      ▼                                               ▼
      ┌───────────────────────────────┐               ┌───────────────────────────────┐
      │ BƯỚC 1: CHỨNG MINH BAO HÀM    │               │ BƯỚC 2: CHỨNG MINH BAO HÀM    │
      │ THỨ THUẬN (VẾ TRÁI ⊂ VẾ PHẢI) │               │ THỨ ĐẢO (VẾ PHẢI ⊂ VẾ TRÁI)   │
      └───────────────┬───────────────┘               └───────────────┬───────────────┘
                      │                                               │
                      │  Xây dựng bài toán tối ưu                     │  Dãy x_k* ∈ N_Γ(x_k),
                      │  nhiễu tham số (P_t)                          │  (x_k, x_k*) ──► (x̄, x̄*)
                      │                                               │
                      │  Ước lượng nghiệm xt:                         │  Tách nón lồi đa diện Θ̃
                      │  ‖x̄ + tv - xt‖ = o(t)                         │  pos{b_i | i ∈ I_0} độc lập
                      │                                               │
                      │  Khảo sát giới hạn nón                        │  Áp dụng Bổ đề Chặn sai số
                      │  pháp tuyến qua t ↓ 0                         │  Hoffman cho đa diện
                      │                                               │
                      ▼                                               ▼
      ┌───────────────────────────────┐               ┌───────────────────────────────┐
      │ KẾT HỢP HAI CHIỀU:            │               │ XỬ LÝ TRƯỜNG HỢP µ ∉ T_N(λ):  │
      │ THIẾT LẬP ĐẲNG THỨC NÓN HOÀN  │               │ HIỆU CHỈNH DÃY DỊCH CHUYỂN    │
      │ HẢO VÀ CÔNG THỨC DN_Γ HIỂN    │               │ BẰNG TỔ HỢP TUYẾN TÍNH DƯƠNG  │
      └───────────────────────────────┘               └───────────────────────────────┘
  1. Thiết lập Bổ đề Nón tới hạn (Lemma 2.1): Chứng minh biểu diễn đối ngẫu của nón tới hạn $N_K(v) = { \nabla q(\bar{x})^T \mu \mid \mu^T \nabla q(\bar{x})v = 0, \mu \in T_{N_\Theta(\bar{y})}(\lambda) }$ thông qua phân tích cực nón $N_K(v) = K^\circ \cap {v}^\perp = (\nabla q(\bar{x})^T N_\Theta(\bar{y}) + \mathbb{R}\bar{x}^*) \cap {v}^\perp$.
  2. Ước lượng tiệm cận nghiệm xấp xỉ: Xây dựng dãy nghiệm bài toán tối ưu phụ thuộc tham số $(P_t): \min { | \bar{x} + tv + \frac{1}{2}t^2 z + \alpha(\bar{x}^* + tv^) - x |^2 \mid x \in \Gamma \cap B_r(\bar{x}) }$. Sử dụng tính dưới chính quy mêtric $d_\Gamma(\bar{x} + tv + \frac{1}{2}t^2 z) \le \kappa_1 d_\Theta(q(\bar{x} + tv + \frac{1}{2}t^2 z)) \le o(t^2)$ để thu được ước lượng then chốt: $$|\bar{x} + tv - x_t|^2 \le o(t^2) + o(|x_t - \bar{x}|^2) \implies |\bar{x} + tv - x_t| = o(t)$$ từ đó khẳng định $(v, v^) \in T_{\text{gph} N_\Gamma}(\bar{x}, \bar{x}^*)$.
  3. Triệt tiêu giả thiết phụ bằng Bổ đề Hoffman: Trong chiều chứng minh ngược lại, với dãy $(x_k, x_k^) \to (\bar{x}, \bar{x}^)$ có $x_k^* \in N_\Gamma(x_k)$, áp dụng Bổ đề Hoffman cho ánh xạ $\Psi_{I_0}(x^) := { \lambda \in \mathbb{R}^m \mid \nabla q(\bar{x})^T \lambda = x^, \langle a_i, \lambda \rangle \le 0 }$ để tìm được dãy nhân tử xấp xỉ $\tilde{\lambda}^k \in \Psi_{I_0}(\bar{x}^)$ thỏa mãn $|\lambda^k - \tilde{\lambda}^k| \le \beta t_k (|v_k^| + c_1 c_2 |v_k|)$, hoàn tất chứng minh đẳng thức tập hợp tuyệt đối.

Data và phân tích

Đối tượng phân tích trong luận án là cấu trúc không gian hình học biến phân thuần túy. Các công cụ toán học và đại lượng định lượng bao gồm:

  • Môđun dưới chính quy mêtric: $\text{subreg},F(\bar{x}|\bar{y}) := \inf { \kappa \in \mathbb{R}+ \mid \exists r > 0, d{F^{-1}(\bar{y})}(x) \le \kappa d_F(x)(\bar{y}), \forall x \in B_r(\bar{x}) }$.
  • Môđun chính xác của tính ổn định xiên $\text{tilt}(f, \bar{x})$.
  • Tính xác định dương của ma trận Hessian Lagrange suy rộng $\nabla_{xx}^2 L_g(\bar{x}, \alpha, \lambda) > 0$.

Tính vững (Robustness) của các kết quả toán học được kiểm chứng qua việc xây dựng các phản ví dụ sắc bén (chẳng hạn Ví dụ 2.11 và Ví dụ 2.12), chứng minh tính đúng đắn của công thức khi điều kiện chính quy mêtric bị phá vỡ hoàn toàn ($(0,0) \notin \text{int}(q(\bar{x}) + \nabla q(\bar{x})\mathbb{R}^2 - \Theta)$) và điều kiện BEPP của Gfrerer-Mordukhovich không thỏa mãn.


Phát hiện đột phá và implications

Những phát hiện then chốt

Luận án đã đạt được 5 phát hiện khoa học mang tính đột phá:

┌────────────────────────────────────────────────────────────────────────────────────────┐
│                        5 PHÁT HIỆN ĐỘT PHÁ CỦA LUẬN ÁN                                 │
├────────────────────────────────────────────────────────────────────────────────────────┤
│ 1. CÔNG THỨC HIỂN DN_Γ DƯỚI MSCQ THUẦN TÚY                                            │
│    Hợp nhất và mở rộng hoàn toàn các kết quả của Dontchev-Rockafellar (1996),          │
│    Henrion et al. (2013), Gfrerer-Mordukhovich (2015) khi loại bỏ hoàn toàn BEPP/CRCQ. │
├────────────────────────────────────────────────────────────────────────────────────────┤
│ 2. ĐẶC TRƯNG ỔN ĐỊNH XIÊN QUA ĐẠO HÀM ĐỒ THỊ DƯỚI GRADIENT D(∂f)                      │
│    Lần đầu tiên thiết lập cầu nối trực tiếp giữa D(∂f) và tilt stability thay cho      │
│    việc sử dụng dưới vi phân bậc hai ∂²f phức tạp.                                     │
├────────────────────────────────────────────────────────────────────────────────────────┤
│ 3. NỚI LỎNG ĐIỀU KIỆN ĐỦ SSOSC DƯỚI DUY NHẤT MSCQ                                      │
│    Chứng minh SSOSC là điều kiện đủ cho ổn định xiên trong NLP dưới MSCQ, vượt trội    │
│    hẳn so với đòi hỏi MFCQ + CRCQ của Mordukhovich-Outrata (2012).                     │
├────────────────────────────────────────────────────────────────────────────────────────┤
│ 4. TIÊU CHUẨN TƯỜNG MINH CHO QUY HOẠCH TOÀN PHƯƠNG QCQP                               │
│    Xây dựng tiêu chuẩn ổn định xiên dạng đại số ma trận tường minh cho bài toán        │
│    toàn phương có một ràng buộc toàn phương thỏa mãn MSCQ.                             │
├────────────────────────────────────────────────────────────────────────────────────────┤
│ 5. ĐẶC TRƯNG TĨNH LẶNG CÔ LẬP CHO PHƯƠNG TRÌNH SUY RỘNG                                │
│    Thiết lập điều kiện cần và đủ cho tính isolated calmness của ánh xạ nghiệm          │
│    S(p) = {x | 0 ∈ Φ(p, x) + N_Γ(x)} dựa trên toán tử DN_Γ vừa tính toán.              │
└────────────────────────────────────────────────────────────────────────────────────────┘
  1. Công thức đạo hàm đồ thị tổng quát cho nón pháp tuyến: Xác lập Định lý 2.10 chứng minh $DN_\Gamma(\bar{x}|\bar{x}^*)(v) = \nabla^2(\lambda^T q)(\bar{x})v + N_K(v)$ với $\lambda \in \Lambda(v)$ cho tập lồi đa diện $\Theta$ dưới duy nhất giả thiết MSCQ.
  2. Tiêu chuẩn ổn định xiên thông qua đạo hàm đồ thị dưới gradient: Định hình một phương pháp hoàn toàn mới: Cực tiểu địa phương của hàm $f$ chính quy gần kề và liên tục dưới vi phân là ổn định xiên khi và chỉ khi đạo hàm đồ thị của ánh xạ dưới vi phân $D(\partial f)(\bar{x}|0)$ có tính xác định dương.
  3. Mở rộng điều kiện đủ bậc hai mạnh SSOSC trong quy hoạch phi tuyến: Bẻ gãy ranh giới của các giả thiết kinh điển. Trong khi Mordukhovich và Outrata (2012) bắt buộc phải có đồng thời MFCQ và CRCQ, luận án chứng minh SSOSC vẫn là điều kiện đủ bảo đảm ổn định xiên khi chỉ có MSCQ.
  4. Đặc trưng tường minh cho bài toán quy hoạch toàn phương có ràng buộc toàn phương (QCQP): Khai thác cấu trúc ma trận của bài toán $\min { \frac{1}{2}x^T A x + a^T x \mid \frac{1}{2}x^T B x + b^T x \le 0 }$ dưới MSCQ, đưa ra tiêu chuẩn ổn định xiên khả tích trực tiếp trên hệ số đại số.
  5. Đặc trưng hóa tính tĩnh lặng cô lập (Isolated Calmness) của phương trình suy rộng: Sử dụng công thức $DN_\Gamma$ để xác lập tiêu chuẩn tương đương cho tính tĩnh lặng cô lập của ánh xạ nghiệm $S(p) := { x \in \mathbb{R}^n \mid 0 \in \Phi(p, x) + N_\Gamma(x) }$ thông qua điều kiện hạch rỗng: $D\Phi(\bar{p}, \bar{x})(0, v) + DN_\Gamma(\bar{x}|-\Phi(\bar{p}, \bar{x}))(v) \subset {0} \implies v = 0$.

Implications đa chiều

  • Ý nghĩa lý thuyết: Thống nhất các nhánh nghiên cứu phân mảnh trong giải tích biến phân bậc hai; khẳng định vị thế trung tâm của tính dưới chính quy mêtric (MSCQ) như là điều kiện chuẩn hóa yếu nhất mà vẫn bảo toàn được tính toán giải tích vi phân bậc hai chính xác.
  • Ý nghĩa phương pháp luận: Khẳng định đạo hàm đồ thị dưới gradient $D(\partial f)$ là một công cụ giải tích mạnh mẽ, linh hoạt và thuận tiện hơn nhiều so với dưới vi phân bậc hai Mordukhovich $\partial^2 f$ trong việc phân tích độ nhạy và tính ổn định tối ưu.
  • Ý nghĩa thuật toán và ứng dụng: Đặt nền móng hội tụ vững chắc cho các thuật toán tối ưu số hiện đại: Thuật toán Quy hoạch Toàn phương Tuần tự (SQP), Phương pháp Nhân tử Lagrange Tăng cường (Augmented Lagrangian Methods), và các thuật toán Newton phi trơn giải phương trình suy rộng trong điều kiện suy biến ràng buộc (khi MFCQ bị vi phạm).

Limitations và Future Research

Luận án nghiêm túc chỉ ra các giới hạn lý thuyết và điều kiện biên:

  1. Giới hạn không gian: Toàn bộ các kết quả chính được thiết lập trong không gian Euclid hữu hạn chiều $\mathbb{R}^n$. Việc chuyển đổi sang không gian Banach vô hạn chiều sẽ gặp rào cản lớn do tính compact yếu của hình cầu đơn vị không còn được bảo toàn.
  2. Giới hạn tập đích $\Theta$: Công thức Định lý 2.10 đòi hỏi $\Theta$ là tập lồi đa diện để sử dụng Bổ đề Hoffman và tính hữu hạn của biểu diễn mặt nón. Kết quả chưa bao quát trực tiếp cho các nón phi đa diện quan trọng như nón nửa xác định dương (semidefinite cone $\mathcal{S}+^n$) trong Quy hoạch Bán xác định (SDP) hay nón cấp hai (second-order cone / Lorentz cone $\mathcal{K}{\text{SOC}}$).
  3. Giới hạn độ trơn: Ánh xạ ràng buộc $q$ và hàm mục tiêu $f$ bắt buộc phải thuộc lớp khả vi liên tục hai lần ($C^2$).
                ┌────────────────────────────────────────────────────────┐
                │        CHƯƠNG TRÌNH NGHIÊN CỨU 10 NĂM TỚI               │
                └───────────────────────────┬────────────────────────────┘
                                            │
           ┌────────────────────────────────┼────────────────────────────────┐
           ▼                                ▼                                ▼
┌──────────────────────┐         ┌──────────────────────┐         ┌──────────────────────┐
│  HƯỚNG 1: MỞ RỘNG    │         │  HƯỚNG 2: BẬC TRƠN   │         │  HƯỚNG 3: THUẬT TOÁN │
│  NÓN KHÔNG ĐA DIỆN   │         │  C¹,¹ VÀ HÀM GHÉP    │         │  TỐI ƯU SỐ & ALM/SQP │
├──────────────────────┤         ├──────────────────────┤         ├──────────────────────┤
│ Mở rộng Định lý 2.10 │         │ Xử lý bài toán tối ưu│         │ Đánh giá tốc độ hội tụ│
│ cho Nón SDP S₊ⁿ và   │         │ không trơn dạng hàm  │         │ siêu tuyến tính khi  │
│ Lorentz Cone K_SOC   │         │ ghép composite C¹,¹  │         │ vi phạm MFCQ         │
└──────────────────────┘         └──────────────────────┘         └──────────────────────┘

Chương trình nghiên cứu tiếp nối (Future Research Agenda):

  • Hướng 1: Mở rộng công thức tính đạo hàm đồ thị của ánh xạ nón pháp tuyến cho quy hoạch nón tổng quát (Conic Programming), đặc biệt là nón $\mathcal{S}+^n$ và $\mathcal{K}{\text{SOC}}$ dưới điều kiện MSCQ.
  • Hướng 2: Nghiên cứu tính ổn định xiên cho lớp hàm mục tiêu không trơn dạng hàm ghép $f(x) = g(h(x))$ với $g$ lồi chính thường và $h$ khả vi $C^{1,1}$.
  • Hướng 3: Thiết kế và phân tích tốc độ hội tụ siêu tuyến tính của thuật toán SQP và thuật toán Newton nửa trơn (Semismooth Newton) giải phương trình suy rộng dựa trên môđun ổn định xiên $\text{tilt}(f, \bar{x})$ đã tính toán.

Tác động và ảnh hưởng

  • Tác động học thuật quốc tế: Công trình đóng góp trực tiếp vào dòng chảy khoa học đỉnh cao thông qua 02 bài báo công bố trên tạp chí quốc tế hàng đầu danh mục ISI/SCI-Q1: SIAM Journal on Optimization (SIOPT) – tạp chí danh giá bậc nhất của Hiệp hội Toán học Ứng dụng và Công nghiệp Hoa Kỳ (SIAM). Kết quả này định vị toán học Việt Nam trên bản đồ nghiên cứu giải tích biến phân thế giới.
  • Tác động tới phần mềm tối ưu hóa: Cung cấp cơ sở toán học chuẩn xác để nâng cấp các solver tối ưu hóa phi tuyến phi lồi thương mại và mã nguồn mở (như IPOPT, KNITRO, SNOPT), giúp các bộ giải không bị dừng đột ngột (breakdown) tại các điểm tối ưu suy biến thỏa mãn MSCQ nhưng vi phạm MFCQ.
  • Tác động kỹ thuật thực tiễn: Ứng dụng trong điều khiển tối ưu thời gian thực cho hệ thống robot chuyển động phức tạp, tính toán điều độ kinh tế lưới điện thông minh (Smart Grid Optimal Power Flow), và tối ưu hóa hình học cấu trúc (Structural Topology Optimization) – nơi các ràng buộc cơ học thường xuyên vi phạm các điều kiện chuẩn hóa cổ điển.

Đối tượng hưởng lợi

  • Nghiên cứu sinh và Giảng viên ngành Toán: Tiếp cận một tài liệu chuyên khảo mẫu mực về giải tích biến phân hiện đại, nắm bắt các kỹ thuật xử lý vi phân suy rộng bậc hai phức tạp.
  • Các nhà nghiên cứu Tối ưu hóa thuật toán (R&D): Sử dụng các đặc trưng $DN_\Gamma$ và tilt stability để chứng minh sự hội tụ của các thuật toán tối ưu hóa phi trơn và tối ưu hóa hệ biến phân.
  • Kỹ sư Hệ thống và Điều khiển Tự động: Áp dụng lý thuyết ổn định xiên để thiết kế các bộ điều khiển bền vững (robust controllers) trước nhiễu tham số tuyến tính và phi tuyến.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Đóng góp độc đáo nhất là Định lý 2.10, thiết lập công thức tường minh tính đạo hàm đồ thị ánh xạ nón pháp tuyến $DN_\Gamma(\bar{x}|\bar{x}^*)(v) = \nabla^2(\lambda^T q)(\bar{x})v|_{\lambda \in \Lambda(v)} + N_K(v)$ cho tập lồi đa diện $\Theta$ chỉ dưới duy nhất điều kiện MSCQ. Định lý này mở rộng và hoàn thiện Lý thuyết Vi phân suy rộng bậc hai Mordukhovich-Rockafellar, xóa bỏ hoàn toàn giả thiết BEPP của Gfrerer-Mordukhovich (2015) và giả thiết chính quy mêtric của Henrion et al. (2013).

2. Đột phá phương pháp luận so với các nghiên cứu quốc tế tiền nhiệm thể hiện ở đâu?

So với nghiên cứu của Henrion et al. (2013) (đòi hỏi $M_q$ chính quy mêtric quanh điểm xét) và Gfrerer & Mordukhovich (2015) (đòi hỏi MSCQ + BEPP), luận án đã kết hợp lược đồ xấp xỉ biến phân trực tiếp của Gfrerer-Outrata (2014) với kỹ thuật đối ngẫu của Ioffe-Outrata (2008) và Bổ đề Chặn sai số Hoffman trên nón lồi đa diện. Đột phá này cho phép kiểm soát sai số tiệm cận bậc hai $o(t^2)$ mà không cần tính chính quy mêtric đều trên lân cận.

3. Phát hiện nào gây bất ngờ nhất về mặt toán học trong luận án?

Phát hiện bất ngờ nhất là tính đủ của điều kiện SSOSC đối với tính ổn định xiên trong quy hoạch phi tuyến dưới duy nhất chuẩn hóa MSCQ. Trước đây, y văn quốc tế (đặc biệt là Mordukhovich & Outrata, 2012) luôn mặc định rằng để SSOSC bảo đảm ổn định xiên thì cần đồng thời cả MFCQ và CRCQ. Luận án chứng minh rằng thông qua lăng kính đạo hàm đồ thị dưới gradient $D(\partial f)$, chỉ cần MSCQ – điều kiện yếu hơn rất nhiều – là đủ để thiết lập kết quả này.

4. Luận án có cung cấp giao thức tái lập/kiểm chứng toán học đầy đủ không?

Có. Toàn bộ các định lý trong Chương 2 và Chương 3 đều đi kèm chứng minh toán học giải tích biến phân tường minh từng bước (constructive analytical proofs). Mọi đánh giá sai số chuẩn Euclid, nón tiếp tuyến, và phép chiếu mêtric đều được tường minh hóa kèm theo các phản ví dụ số học 2 chiều ($n=2, m=2$) có tính toán cụ thể để độc giả kiểm chứng độc lập.

5. Kế hoạch nghiên cứu 10 năm tới được phác thảo thế nào?

Kế hoạch 10 năm tập trung vào việc: (i) Mở rộng công thức $DN_\Gamma$ sang không gian vô hạn chiều và nón phi đa diện (SDP, SOCP); (ii) Phát triển lý thuyết ổn định xiên cấp cao (higher-order tilt stability); (iii) Tích hợp các điều kiện ổn định xiên vào việc thiết kế thế hệ thuật toán tối ưu số mới chịu lỗi suy biến ràng buộc.


Kết luận

Luận án tiến sĩ của tác giả Lê Văn Hiển đã tạo nên một dấu ấn khoa học xuất sắc trong chuyên ngành Giải tích toán học và Lý thuyết Tối ưu với 5 kết luận cốt lõi:

  1. Thiết lập thành công công thức hiển của đạo hàm đồ thị ánh xạ nón pháp tuyến $DN_\Gamma(\bar{x}|\bar{x}^*)$ khi $\Theta$ là tập lồi đa diện dưới điều kiện chuẩn hóa yếu nhất là tính dưới chính quy mêtric (MSCQ).
  2. Xóa bỏ hoàn toàn giả thiết phụ BEPP và CRCQ tồn tại trong các công trình quốc tế kinh điển, hợp nhất toàn diện các kết quả tính toán $DN_\Gamma$ trong giải tích biến phân bậc hai.
  3. Đề xuất phương pháp luận mới sử dụng đạo hàm đồ thị dưới gradient $D(\partial f)$ để đặc trưng hóa trọn vẹn tính ổn định xiên (tilt stability) cho các bài toán tối ưu không ràng buộc với hàm mục tiêu chính quy gần kề.
  4. Chứng minh điều kiện đủ bậc hai mạnh SSOSC bảo đảm tính ổn định xiên trong quy hoạch phi tuyến chỉ với giả thiết MSCQ, cải tiến vượt bậc kết quả kinh điển của Mordukhovich và Outrata (2012).
  5. Xây dựng các tiêu chuẩn đại số ma trận tường minh cho tính ổn định xiên trong quy hoạch toàn phương có ràng buộc toàn phương (QCQP) và tính tĩnh lặng cô lập của nghiệm phương trình suy rộng.

Công trình không chỉ nâng cao vị thế học thuật của giải tích biến phân tại Việt Nam với 02 bài báo công bố trên tạp chí quốc tế hàng đầu SIAM Journal on Optimization, mà còn mở ra những hướng nghiên cứu chiến lược về giải tích biến phân bậc hai cho các cấu trúc nón phi đa diện và thuật toán tối ưu số hiện đại trong tương lai.