Luận án: Phương pháp giải bài toán cân bằng và bất đẳng thức biến phân hai cấp - Trần Thị Hoàng Anh

Tổng quan về luận án

Luận án tiến sĩ toán học với đề tài "Phương pháp giải một vài lớp bài toán cân bằng và bất đẳng thức biến phân hai cấp" của tác giả Trần Thị Hoàng Anh (chuyên ngành Toán ứng dụng, mã số: 9 46 01 12, hoàn thành năm 2019 tại Học viện Kỹ thuật Quân sự dưới sự hướng dẫn khoa học của PGS. Phạm Ngọc Anh và GS. Phạm Thế Long) là một công trình học thuật tiên phong trong lĩnh vực giải tích lồi, lý thuyết toán tử đơn điệu và tối ưu hóa phi tuyến hiện đại.

Bối cảnh khoa học của nghiên cứu bắt nguồn từ nhu cầu cấp thiết trong việc giải quyết các bài toán tối ưu phân cấp (bilevel optimization), nơi miền ràng buộc của bài toán cấp trên (upper-level problem) không được biểu diễn tường minh mà là tập nghiệm ngầm của một bài toán cấp dưới (lower-level problem), cụ thể là bài toán bất đẳng thức biến phân (Variational Inequality - $VI$) hoặc bài toán cân bằng (Equilibrium Problem - $EP$). Các mô hình này đóng vai trò nền tảng trong cân bằng mạng giao thông Wardrop, định tuyến mạng viễn thông, bài toán cân bằng kinh tế Nash-Cournot và xử lý tín hiệu/hình ảnh số.

Khoảng trống nghiên cứu (research gap) then chốt được xác định là: Hầu hết các phương pháp giải bài toán hai cấp truyền thống—chẳng hạn như phương pháp hiệu chỉnh Tikhonov hoặc các lược đồ điểm gần kề lồng nhau kiểu Moudafi—đều đòi hỏi phải giải chính xác các bài toán phụ đơn điệu ở mỗi bước lặp, hoặc yêu cầu điều kiện đơn điệu mạnh rất ngặt nghèo của các song hàm và toán tử giá. Điều này dẫn đến sự bùng nổ khối lượng tính toán trên máy tính, và dãy nghiệm xấp xỉ dễ bị phân kỳ khi các bài toán phụ chỉ giải được gần đúng.

Để giải quyết triệt để vấn đề này, luận án tập trung trả lời 4 câu hỏi nghiên cứu và kiểm chứng 4 giả thuyết khoa học:

• Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để mở rộng kỹ thuật ánh xạ nghiệm từ bài toán bất đẳng thức biến phân sang bài toán cân bằng tổng quát trong không gian $\mathbb{R}^n$ mà vẫn bảo toàn tính chất tựa co hoặc tựa không giãn?
• Câu hỏi nghiên cứu 2 (RQ2): Có thể thiết lập một thuật toán một vòng lặp (single-loop projection) giải bài toán bất đẳng thức biến phân hai cấp $BVI(F, G, C)$ với độ phức tạp tính toán thấp hơn các thuật toán hai vòng lặp của Kopelevich hay Pham Ngoc Anh - Nguyen Van Muu không?
• Câu hỏi nghiên cứu 3 (RQ3): Thiết kế thuật toán chiếu dưới đạo hàm kết hợp xấp xỉ gắn kết (viscosity approximation) như thế nào để đạt sự hội tụ mạnh cho bài toán $VIEP(F, f, C)$ khi song hàm $f$ chỉ thỏa mãn tính giả đơn điệu và điều kiện kiểu Lipschitz?
• Câu hỏi nghiên cứu 4 (RQ4): Bằng cách nào xây dựng được thuật toán chiếu giải bài toán cân bằng $EP(f, C)$ và bài toán $EVIP(g, F, C)$ khi giả thiết đơn điệu được nới lỏng về tính chất para-đơn điệu?

Hệ thống giả thuyết nghiên cứu tương ứng:

• Giả thuyết 1 (H1): Ánh xạ nghiệm dạng Fukushima mở rộng cho bài toán cân bằng $EP(f, C)$ là tựa co với hệ số co tường minh $\delta < 1$ dưới điều kiện $\gamma$-giả đơn điệu mạnh và Lipschitz-type, và là tựa không giãn dưới điều kiện giả đơn điệu.
• Giả thuyết 2 (H2): Thuật toán chiếu kết hợp nguyên lý điểm bất động Banach cho phép dãy lặp hội tụ mạnh về nghiệm duy nhất của $BVI(F, G, C)$ mà chỉ cần 01 phép chiếu metric và 01 lần tính giá trị toán tử mỗi bước lặp.
• Giả thuyết 3 (H3): Kỹ thuật dưới vi phân chéo $\epsilon$-xấp xỉ $\partial_2^\epsilon f(x,x)$ triệt tiêu được sai số tích lũy của bài toán phụ, đảm bảo dãy lặp hội tụ mạnh tới nghiệm duy nhất của $VIEP(F, f, C)$.
• Giả thuyết 4 (H4): Phương pháp chiếu siêu phẳng phân tách kết hợp kỹ thuật lai ghép đảm bảo sự hội tụ của thuật toán cho bài toán cân bằng dưới điều kiện para-đơn điệu mà không cần tính toán ma trận Jacobi.

Khung lý thuyết (theoretical framework) của luận án tích hợp Giải tích biến phân (Variational Analysis), Lý thuyết điểm bất động (Fixed Point Theory), và Lý thuyết toán tử đơn điệu (Monotone Operator Theory). Phạm vi nghiên cứu bao quát không gian Hilbert thực vô hạn chiều $H$ và không gian Euclide hữu hạn chiều $\mathbb{R}^n$. Toàn bộ thuật toán đề xuất được lập trình, kiểm chứng thực nghiệm trên phần mềm Matlab R2013a chạy trên phần cứng Laptop CPU Intel Core i3-3110M @ 2.40GHz, 4GB RAM, đo lường chi tiết qua số bước lặp ($Iter. (k)$) và thời gian thực thi CPU ($CPU-times/s$).

Literature Review và Positioning

Lý thuyết bất đẳng thức biến phân được khởi xướng năm 1966 bởi nhà toán học Ý G. Stampacchia khi nghiên cứu bài toán điều khiển tối ưu và bài toán biên cho phương trình đạo hàm riêng elliptic. Đến năm 1994, E. Blum và W. Oettli đã thống nhất bài toán bất đẳng thức biến phân, bài toán quy hoạch tối ưu lồi, bài toán bù phi tuyến ($CP$) và bài toán cân bằng Nash trong lý thuyết trò chơi không hợp tác thành một mô hình tổng quát mang tên Bài toán cân bằng ($EP$).

Trong tiến trình phát triển thuật toán cho bài toán đơn cấp và hai cấp, y văn thế giới hình thành hai trường phái tiếp cận chủ đạo:

1. Trường phái ánh xạ điểm gần kề và hiệu chỉnh Tikhonov: Tiêu biểu là các công trình của R. T. Rockafellar (1976) với thuật toán Proximal Point Algorithm (PPA), A. Moudafi (2010) với phương pháp hiệu chỉnh Tikhonov cho bài toán cân bằng hai cấp $BEP(f, g, C)$. Trường phái này đưa bài toán hai cấp về dãy các bài toán đơn điệu đơn cấp $EP(h_\epsilon, C)$ với $h_\epsilon(x, y) = f(x, y) + \epsilon g(x, y)$. Tuy nhiên, điểm hạn chế chí mạng là mỗi bước lặp đòi hỏi tìm nghiệm giải tích chính xác của bài toán phụ—một nhiệm vụ bất khả thi trong tính toán số thực tế.
2. Trường phái đạo hàm tăng cường (Extragradient) và phép chiếu: Khởi xướng bởi G. Kopelevich (1976) cho bài toán điểm yên ngựa và $VI$, sau đó được I. Yamada (2001) phát triển cho bài toán $VIFIX$, L. D. Muu (2003, 2008), P. N. Anh và N. V. Muu (2014) mở rộng cho bài toán cân bằng kết hợp tìm điểm bất động của ánh xạ không giãn. Nhược điểm của các thuật toán hai cấp dạng này là cấu trúc vòng lặp lồng nhau (nested loops), đòi hỏi tính toán phép chiếu metric nhiều lần và bước lặp trong phải thỏa mãn điều kiện dừng ngặt.

Bảng so sánh định vị học thuật của luận án với các công trình quốc tế tiêu biểu:

Luận án đã định vị chính xác khoảng trống khoa học: Thiết lập các thuật toán chiếu cải tiến và chiếu dưới đạo hàm có cấu trúc một vòng lặp trực tiếp, triệt tiêu hoàn toàn sự phụ thuộc vào các bài toán phụ giải chính xác, đồng thời hạ thấp các giả thiết đơn điệu từ "đơn điệu mạnh" xuống "giả đơn điệu" và "para-đơn điệu".

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

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

Luận án mở rộng trực tiếp các định lý nền tảng của giải tích lồi và lý thuyết điểm bất động:

1. Mở rộng nguyên lý ánh xạ nghiệm Fukushima (1992): Tác giả đã tổng quát hóa cấu trúc ánh xạ nghiệm từ không gian bài toán $VI$ sang bài toán cân bằng $EP(f, C)$ trong $\mathbb{R}^n$. Bằng việc áp dụng định lý Moreau-Rockafellar về dưới vi phân của tổng hai hàm lồi, luận án chứng minh rằng nghiệm của $EP(f, C)$ trùng khớp với tập điểm bất động $Fix(S)$ của ánh xạ: $$S(x) = \text{argmin}\left{ f(x, y) + \frac{1}{2}\langle y - x, G(y - x)\rangle : y \in C \right}$$ với $G$ là ma trận đối xứng xác định dương cấp $n$.
2. Chứng minh tính tựa co (Quasi-contraction) tường minh: Khi song hàm $f$ là $\gamma$-giả đơn điệu mạnh và thỏa mãn điều kiện kiểu Lipschitz với hằng số $c_1, c_2 > 0$, với tham số điều chỉnh $\lambda \in (0, \frac{1}{2c_1})$ và $\gamma > c_2$, ánh xạ nghiệm $S(x)$ đạt tính chất tựa co với hệ số co chính xác: $$\delta = \frac{1}{\sqrt{1 + 2\lambda(\gamma - c_2)}} \in (0, 1)$$ kéo theo bất đẳng thức co ngặt: $|S(x) - x^| \le \delta |x - x^|, \forall x \in C, x^* \in \text{Sol}(f, C)$.
3. Chứng minh tính tựa không giãn (Quasi-nonexpansiveness): Dưới giả thiết $f$ chỉ là giả đơn điệu thông thường, bằng cách xây dựng ánh xạ nghiệm hai tầng $z_x = \text{argmin}{\lambda f(x, y) + \frac{1}{2}|y - x|^2 : y \in C}$, luận án chứng minh thành công tính chất tựa không giãn: $|S(x) - x^| \le |x - x^|$ với $\lambda \in (\frac{1}{2c_1}, \frac{1}{2c_2})$.

graph TD
    A["Bài toán cân bằng EP(f, C)"] --> B["Ánh xạ nghiệm S(x) dạng Fukushima cải tiến"]
    B --> C{"Giả thiết đơn điệu của song hàm f"}
    C -->|gamma-giả đơn điệu mạnh + Lipschitz-type| D["Tính tựa co: ||S(x) - x*|| <= delta ||x - x*|| với delta < 1"]
    C -->|Giả đơn điệu + Lipschitz-type| E["Tính tựa không giãn: ||S(x) - x*|| <= ||x - x*||"]
    D --> F["Thiết lập thuật toán lặp điểm bất động đơn giản"]
    E --> F
    F --> G["Hội tụ mạnh về nghiệm duy nhất x* trong Hilbert H"]

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

Khung phân tích của luận án kết hợp chặt chẽ 3 trụ cột toán học:

• Lý thuyết toán tử đơn điệu ngược ($\eta$-ism): Ánh xạ $G$ thỏa mãn $\langle G(x) - G(y), x - y \rangle \ge \eta |G(x) - G(y)|^2$, cho phép toán tử $I - \lambda G$ là không giãn khi $\lambda \in (0, 2\eta)$.
• Nguyên lý xấp xỉ gắn kết (Viscosity Approximation Method): Toán tử co $T_\lambda = I - \lambda\mu F$ với $F$ là $\beta$-đơn điệu mạnh và $L$-liên tục Lipschitz. Hằng số co ngặt được lượng hóa chính xác: $$\tau = 1 - \sqrt{1 - \mu(2\beta - \mu L^2)} \in (0, 1] \quad \text{với} \quad \mu \in \left(0, \frac{2\beta}{L^2}\right)$$
• Kỹ thuật dưới vi phân chéo $\epsilon$-xấp xỉ ($\partial_2^\epsilon f(x,x)$): Giải quyết triệt để sự không trơn của song hàm mà không cần giả thiết khả vi Frechet hay tính toán ma trận Hessian phức tạp.

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

Thiết kế nghiên cứu

Luận án tuân thủ thế giới quan thực chứng suy diễn logic-tiên đề (deductive-axiomatic positivism) của toán học giải tích hiện đại. Quy trình nghiên cứu đi từ việc mô hình hóa toán học các lớp bài toán hai cấp, phân tích tính chất tập nghiệm, thiết kế thuật toán lặp số, chứng minh hội tụ giải tích giải nghiệm lý thuyết, và kiểm thử tính toán trên các bài toán chuẩn (benchmarks).

Quy trình nghiên cứu rigorous

Quy trình thuật toán và kỹ thuật chứng minh hội tụ được chuẩn hóa qua 4 bước lặp toán học chặt chẽ:

1. Phép chiếu Metric và Siêu phẳng phân tách: Sử dụng đặc trưng hình học của nón pháp tuyến ngoài $N_C(x)$ và tính chất đơn điệu mạnh ngược của toán tử chiếu metric $Pr_C$: $$|x - Pr_C(x)|^2 \le |x - y|^2 - |y - Pr_C(x)|^2, \quad \forall x \in H, y \in C$$
2. Thiết lập điều kiện thông số lặp (CP): Các dãy tham số bước nhảy ${\rho_n}, {\beta_n}, {\epsilon_n}, {\alpha_n}$ được kiểm soát bởi hệ điều kiện Maingé nghiêm ngặt: $$\sum_{n=0}^\infty \beta_n = \infty, \quad \sum_{n=0}^\infty \beta_n^2 < \infty, \quad \lim_{n\to\infty} \alpha_n = 0, \quad \sum_{n=0}^\infty \frac{\alpha_n \beta_n}{\rho_n} = \infty$$ Quy tắc chọn tham số thực nghiệm: $\rho_n = 1, \beta_n = n^{-a}, \alpha_n = n^{-b}, \epsilon_n = n^{-1}$ với $a \in (0.5; 1)$ và $b \in (0; 1-a)$.
3. Kỹ thuật chứng minh hội tụ: Sử dụng tính chất Kadec-Klee (trong không gian Hilbert, nếu $x_k \rightharpoonup x^$ và $|x_k| \to |x^|$ thì $x_k \to x^*$), bổ đề hội tụ tựa-Fejér và bất đẳng thức giải tích thực để khẳng định sự hội tụ mạnh không phụ thuộc vào điểm xuất phát $x_0 \in C$.

Data và phân tích

Để kiểm chứng tính đúng đắn và hiệu năng thuật toán, luận án xây dựng hệ thống bài toán kiểm thử (Test Problems) phi tuyến nhiều chiều trong $\mathbb{R}^n$ ($n = 5, 10, 20, 50, 100$).

• Môi trường thực thi: Ngôn ngữ lập trình Matlab R2013a, phần cứng Laptop Intel Core i3-3110M CPU @ 2.40GHz, 4GB RAM.
• Tiêu chuẩn dừng: $|x_{k+1} - x_k| \le Tol = 10^{-4}$ hoặc $10^{-6}$.
• Chỉ số đo lường hiệu năng: Số bước lặp thuật toán ($Iter. (k)$) và tổng thời gian CPU tiêu tốn ($CPU-times/s$).

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

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

Luận án đã công bố 5 phát hiện cốt lõi mang tính đột phá về mặt thuật toán:

graph LR
    subgraph "4 Thuật toán Đột phá của Luận án"
        T1["Thuật toán 1: Chiếu 1 bước cho BVI(F, G, C)<br/>- 1 phép chiếu Pr_C<br/>- G: eta-ism, F: beta-đơn điệu mạnh"]
        T2["Thuật toán 2: Chiếu dưới đạo hàm cho VIEP(F, f, C)<br/>- Dưới vi phân chéo d_2^eps f(x,x)<br/>- f: Giả đơn điệu, Lipschitz-type"]
        T3["Thuật toán 3: Chiếu dưới đạo hàm cho EVIP(g, F, C)<br/>- F: Para-đơn điệu, đóng yếu<br/>- g: Đơn điệu mạnh"]
        T4["Thuật toán 4: Chiếu cải tiến cho EP(f, C)<br/>- f: Para-đơn điệu<br/>- Chiếu siêu phẳng phân tách"]
    end

1. Thuật toán chiếu một bước giải $BVI(F, G, C)$ (Chương 2): Khắc phục hoàn toàn cấu trúc hai vòng lặp cồng kềnh của thuật toán Pham Ngoc Anh - Nguyen Van Muu (2014). Tại mỗi bước lặp $k$, thuật toán chỉ thực hiện duy nhất: $$y^k = Pr_C(x^k - \lambda G(x^k)), \quad x^{k+1} = \alpha_k (x^k - \lambda\mu F(x^k)) + (1 - \alpha_k) y^k$$ Dãy ${x^k}$ được chứng minh hội tụ mạnh tới nghiệm duy nhất $x^* \in \Omega$ mà không cần bất kỳ vòng lặp phụ nào.
2. Thuật toán chiếu - dưới đạo hàm giải $VIEP(F, f, C)$ (Chương 3): Xử lý thành công lớp bài toán bất đẳng thức biến phân trên tập nghiệm của bài toán cân bằng khi song hàm $f$ chỉ là giả đơn điệu (thay vì đơn điệu mạnh) và thỏa mãn điều kiện kiểu Lipschitz. Thuật toán kết hợp véc-tơ dưới đạo hàm $g^n \in \partial_2^{\epsilon_n} f(x^n, x^n)$ và hướng hiệu chỉnh $u^n = g^n + \alpha_n F(x^n)$, đạt sự hội tụ mạnh giải tích.
3. Thuật toán chiếu - dưới đạo hàm giải $EVIP(g, F, C)$ (Chương 4): Giải quyết trọn vẹn bài toán cân bằng trên tập nghiệm của bài toán bất đẳng thức biến phân dưới điều kiện ánh xạ $F$ chỉ thỏa mãn tính chất para-đơn điệu và đóng yếu trên $C$.
4. Thuật toán kiểu chiếu giải bài toán cân bằng $EP(f, C)$ (Chương 5): Đề xuất cơ chế chiếu siêu phẳng lai ghép cho phép giải bài toán cân bằng đơn cấp với song hàm para-đơn điệu mà không cần đến tính liên tục Lipschitz của toán tử gradient.
5. Dữ liệu thực nghiệm định lượng vượt trội: Thử nghiệm số trên Matlab R2013a chứng minh thuật toán của luận án giảm từ 35% đến 60% thời gian thực thi CPU và giảm từ 40% đến 50% số bước lặp so với các thuật toán đạo hàm tăng cường cổ điển của Kopelevich và thuật toán hiệu chỉnh Tikhonov của Moudafi trên cùng một bài toán kiểm thử kích thước $n = 100$.

Implications đa chiều

• Ý nghĩa lý thuyết: Làm phong phú kho tàng lý thuyết tối ưu hóa không trơn và giải tích biến phân, thiết lập cầu nối toán học vững chắc giữa bài toán cân bằng và lý thuyết điểm bất động của ánh xạ co/tựa không giãn.
• Ý nghĩa phương pháp luận: Cung cấp khung thuật toán tổng quát có thể áp dụng trực tiếp cho các bài toán phân cấp trong không gian vô hạn chiều (như phương trình đạo hàm riêng và điều khiển tối ưu phân bố).
• Ứng dụng thực tiễn:
  • Giao thông vận tải: Tối ưu hóa lưu lượng dòng xe trên mạng lưới giao thông đô thị dựa trên nguyên lý cân bằng Wardrop hai cấp (cấp trên quản lý đèn tín hiệu/thu phí, cấp dưới mô hình hóa hành vi người tham gia giao thông).
  • Viễn thông & Mạng máy tính: Giải quyết bài toán định tuyến tối ưu cân bằng tải (traffic engineering) và phân bổ tài nguyên băng thông công bằng trong mạng không dây thế hệ mới.
  • Kinh tế học & Quản trị: Xác định điểm cân bằng Nash-Cournot trong thị trường năng lượng cạnh tranh không hoàn hảo có ràng buộc chính sách từ chính phủ.

Limitations và Future Research

Mặc dù đạt được những bước tiến đột phá, luận án cũng chỉ rõ các giới hạn nội tại một cách khách quan:

1. Sự phụ thuộc vào hằng số Lipschitz: Các thuật toán ở Chương 2 và Chương 3 vẫn yêu cầu biết trước các hằng số Lipschitz $L$ hoặc hằng số Lipschitz-type $c_1, c_2$ để chọn độ dài bước nhảy $\lambda, \mu$. Trong thực tế, việc ước lượng chính xác các hằng số này đối với các bài toán phi tuyến phức tạp là rất khó khăn.
2. Ràng buộc về tính lồi: Miền ràng buộc $C$ bắt buộc phải là tập lồi đóng và các song hàm phải có tính chất lồi theo biến thứ hai. Nghiên cứu chưa bao quát các tập ràng buộc không lồi hoặc các bài toán cân bằng tựa lồi tổng quát.
3. Tính toán dưới vi phân chéo: Việc xác định chính xác phần tử dưới vi phân $g^n \in \partial_2^{\epsilon_n} f(x^n, x^n)$ đòi hỏi cấu trúc của song hàm $f$ phải đủ tường minh.

Chương trình nghiên cứu 10 năm tiếp theo (Future Research Agenda):

• Hướng 1: Xây dựng các thuật toán tự thích nghi bước nhảy kiểu Armijo hoặc thuật toán không đòi hỏi hằng số Lipschitz (gradient-free / linesearch-free algorithms) cho bài toán cân bằng hai cấp.
• Hướng 2: Mở rộng phương pháp giải cho bài toán bất đẳng thức biến phân đa trị $MVI(T, C)$ và bài toán điểm bất động đa trị $MFP(F, C)$.
• Hướng 3: Nghiên cứu bài toán cân bằng hai cấp với cấu trúc không lồi (non-convex bilevel equilibrium) ứng dụng trong học sâu (bilevel deep learning / meta-learning).
• Hướng 4: Phát triển các biến thể thuật toán song song và phân tán (distributed/asynchronous algorithms) để xử lý dữ liệu lớn trên nền tảng điện toán đám mây.

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

• Tác động học thuật: Các kết quả cốt lõi của luận án đã được công bố trên 05 bài báo khoa học tại các tạp chí quốc tế uy tín thuộc danh mục ISI/Scopus (như Journal of Global Optimization, Numerical Algorithms, Optimization Letters, Applicable Analysis). Công trình được cộng đồng nghiên cứu giải tích biến phân trong nước và quốc tế trích dẫn, mở ra hướng tiếp cận mới trong việc xử lý các bài toán hai cấp bằng phép chiếu đơn.
• Tác động ngành và công nghệ: Cung cấp thuật toán tối ưu hóa nền tảng cho các phần mềm mô phỏng giao thông thông minh (ITS), hệ thống quản lý lưới điện thông minh (Smart Grid) trong việc cân bằng cung cầu phụ tải điện thời gian thực.
• Chính sách và kinh tế - xã hội: Giúp các nhà hoạch định chính sách kinh tế lượng hóa chính xác các kịch bản cân bằng thuế - thị trường và chính sách giá carbon thông qua mô hình cân bằng tổng thể có ràng buộc phân cấp.

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

• Nghiên cứu sinh và học viên cao học ngành Toán ứng dụng / Khoa học máy tính: Tiếp cận các kỹ thuật chứng minh hội tụ hiện đại, phương pháp giải tích lồi tiên tiến và nguồn tài liệu tham khảo chuẩn mực về bài toán hai cấp.
• Các nhà toán học và chuyên gia tối ưu hóa: Kế thừa các định lý về tính tựa co của ánh xạ nghiệm và kỹ thuật dưới vi phân chéo để phát triển các lớp thuật toán mới cho bài toán cân bằng ngẫu nhiên.
• Kỹ sư R&D trong ngành Giao thông, Viễn thông và AI: Khai thác mã nguồn thuật toán để tích hợp vào các module giải bài toán tối ưu hóa phân bổ tài nguyên mạng và cân bằng hệ thống quy mô lớn.
• Chuyên gia phân tích chính sách kinh tế: Ứng dụng mô hình cân bằng hai cấp để mô phỏng tương tác chiến lược giữa cơ quan quản lý nhà nước và các tác nhân thị trường.

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 của ai?

Đóng góp độc đáo nhất là việc mở rộng nguyên lý ánh xạ nghiệm của K. Fukushima (1992) từ bài toán bất đẳng thức biến phân sang bài toán cân bằng tổng quát $EP(f, C)$ trong không gian $\mathbb{R}^n$, đồng thời chứng minh được tính tựa co ngặt (quasi-contraction) với hệ số co tường minh: $$\delta = \frac{1}{\sqrt{1 + 2\lambda(\gamma - c_2)}} < 1$$ dưới điều kiện song hàm $\gamma$-giả đơn điệu mạnh và thỏa mãn tính chất kiểu Lipschitz ($c_1, c_2$). Kết quả này mở rộng trực tiếp lý thuyết toán tử không giãn của Fukushima và định lý Moreau-Rockafellar về dưới vi phân của hàm lồi.

2. Điểm cải tiến phương pháp luận cốt lõi so với 02 nghiên cứu quốc tế tiền nhiệm là gì?

So với thuật toán của G. Kopelevich (1976) (đòi hỏi 02 phép chiếu metric $Pr_C$ mỗi bước) và thuật toán của P. N. Anh & N. V. Muu (2014) (cấu trúc 02 vòng lặp lồng nhau với vòng trong phải giải xấp xỉ bài toán $VI$), thuật toán chiếu của luận án (Chương 2) chỉ sử dụng 01 phép chiếu metric duy nhất và tính toán 01 lần giá trị toán tử mỗi bước lặp, triệt tiêu hoàn toàn vòng lặp trong, giúp giảm hơn 50% thời gian thực thi CPU trên máy tính.

3. Phát hiện bất ngờ nhất (counter-intuitive) có dữ liệu thực nghiệm chứng minh là gì?

Phát hiện bất ngờ nhất là thuật toán chiếu dưới đạo hàm (Chương 3 & 4) vẫn đạt được sự hội tụ mạnh ($x^k \to x^*$) ngay cả khi song hàm cân bằng $f$ chỉ thỏa mãn tính giả đơn điệu hoặc toán tử $F$ chỉ là para-đơn điệu (những lớp điều kiện yếu hơn rất nhiều so với đơn điệu mạnh). Thực nghiệm trên Matlab R2013a cho thấy dãy lặp ${x^k}$ hội tụ đơn điệu về nghiệm chuẩn xác với sai số $|x_{k+1} - x_k| < 10^{-6}$ chỉ sau 38 bước lặp trên bài toán thử nghiệm kích thước $n = 50$.

4. Luận án có cung cấp giao thức tái lập (Replication Protocol) không?

Có. Luận án cung cấp đầy đủ:

• Công thức giải tích của các bài toán kiểm thử ($TestProb.$) với các ma trận xác định dương $G$, hàm mục tiêu và toán tử phi tuyến.
• Quy tắc thiết lập chính xác các tham số lặp: $\lambda = 0.5/c_1$, bước co $\mu \in (0, 2\beta/L^2)$, các dãy số ${\alpha_k}, {\beta_k}, {\epsilon_k}$ theo điều kiện (CP).
• Toàn bộ thông số môi trường phần mềm (Matlab R2013a) và phần cứng thử nghiệm (Intel Core i3-3110M @ 2.40GHz, 4GB RAM).

5. Lộ trình nghiên cứu 10 năm tới được phác thảo như thế nào?

Lộ trình 10 năm tập trung vào 3 trọng tâm:

1. Phát triển thuật toán tự tìm bước nhảy không phụ thuộc hằng số Lipschitz (giai đoạn 1-3 năm).
2. Mở rộng lý thuyết sang không gian Banach vô hạn chiều và bài toán cân bằng hai cấp ngẫu nhiên (stochastic bilevel equilibrium) (giai đoạn 4-6 năm).
3. Ứng dụng thuật toán phân tán trong học máy phân cấp (bilevel federated learning) và điều khiển tối ưu mạng lưới vạn vật kết nối IoT (giai đoạn 7-10 năm).

Kết luận

Luận án tiến sĩ của tác giả Trần Thị Hoàng Anh đã hoàn thành xuất sắc các mục tiêu nghiên cứu đề ra với 5 đóng góp học thuật cốt lõi:

1. Thiết lập và chứng minh hoàn chỉnh tính tựa co và tựa không giãn của ánh xạ nghiệm dạng Fukushima mở rộng cho bài toán cân bằng trong không gian Euclide $\mathbb{R}^n$.
2. Đề xuất thuật toán chiếu một bước mới cho bài toán bất đẳng thức biến phân hai cấp $BVI(F, G, C)$, chứng minh sự hội tụ mạnh trong không gian Hilbert mà không cần cấu trúc vòng lặp lồng nhau.
3. Xây dựng thuật toán chiếu - dưới đạo hàm giải bài toán bất đẳng thức biến phân trên tập nghiệm của bài toán cân bằng $VIEP(F, f, C)$ dưới điều kiện giả đơn điệu và liên tục kiểu Lipschitz.
4. Phát triển thuật toán chiếu dưới đạo hàm cho bài toán cân bằng trên tập nghiệm bài toán bất đẳng thức biến phân $EVIP(g, F, C)$ với giả thiết toán tử para-đơn điệu và đóng yếu.
5. Đề xuất thuật toán kiểu chiếu siêu phẳng cải tiến giải bài toán cân bằng đơn cấp $EP(f, C)$ dưới điều kiện song hàm para-đơn điệu.

Công trình tạo ra bước tiến quan trọng trong việc chuyển đổi mô hình từ các thuật toán xấp xỉ hai vòng lặp phức tạp sang các lược đồ chiếu đơn hiệu năng cao, mở ra các nhánh nghiên cứu liên ngành giàu tiềm năng giữa toán tối ưu, khoa học dữ liệu và kỹ thuật điều khiển tự động hiện đại.