Tổng quan về luận án

Lý thuyết bài toán cân bằng (Equilibrium Problem - ký hiệu $EP(C, f)$) và lý thuyết điểm bất động (Fixed Point Theory) giữ vị trí trung tâm trong Giải tích phi tuyến, Lý thuyết tối ưu hiện đại và Toán ứng dụng. Được khởi xướng từ công trình tiên phong của Nikaido & Isoda (1955) trong việc tổng quát hóa mô hình cân bằng kinh tế Nash và phát triển đột phá qua bất đẳng thức minimax của Ky Fan (1972) cùng các công trình nền tảng của Blum & Oettli (1994), bài toán cân bằng đã thiết lập một khung giải tích thống nhất chứa đựng nhiều lớp bài toán cốt lõi: quy hoạch toán học lồi, bài toán điểm yên ngựa, bất đẳng thức biến phân ($VI$), bất đẳng thức biến phân đa trị ($MVI$), bài toán bù và bài toán điểm bất động ($FP$). Luận án tiến sĩ toán học của nghiên cứu sinh Nguyễn Văn Hồng với tiêu đề "Một số phương pháp giải bài toán cân bằng trên tập điểm bất động" (chuyên ngành Toán ứng dụng, mã số 9 46 01 12, thực hiện tại Trường Đại học Thăng Long dưới sự hướng dẫn khoa học của PGS. Phạm Ngọc Anh và GS. Lê Dũng Mưu) đã định vị tính tiên phong bằng cách giải quyết bài toán cân bằng hai cấp phức tạp: tìm nghiệm của bài toán cân bằng trên tập giao điểm bất động của một họ các ánh xạ nửa co (demi-contractive mappings) trong không gian Hilbert.

Khoảng trống nghiên cứu (research gap) mà luận án trực tiếp giải quyết xuất phát từ ba nút thắt lý thuyết và tính toán then chốt trong y văn hiện đại:

  1. Chi phí tính toán bài toán phụ quá lớn: Phần lớn các phương pháp giải bài toán cân bằng truyền thống (như phương pháp hiệu chỉnh Tikhonov, phương pháp điểm gần kề Proximal Point) đòi hỏi tại mỗi bước lặp phải giải chính xác một bài toán cân bằng phụ hoặc bài toán tối ưu phi tuyến cấp dưới, điều này gần như bất khả thi về mặt thực thi số khi cấu trúc song hàm phức tạp hoặc không gian có số chiều lớn.
  2. Ràng buộc giao tập điểm bất động phức tạp: Hầu hết các công trình quốc tế trước đây chỉ giải bài toán cân bằng trên tập điểm bất động của một ánh xạ đơn lẻ hoặc ánh xạ không giãn (non-expansive mapping). Nghiên cứu bài toán trên miền ràng buộc $\Omega = \bigcap_{i \in I} \text{Fix}(S_i)$ với $S_i$ là họ các ánh xạ $\beta_i$-nửa co hoặc giả co chặt (strictly pseudo-contractive) vẫn là một bài toán mở đầy thách thức do cấu trúc tập ràng buộc không tường minh.
  3. Đòi hỏi giả thiết ngặt về tính đơn điệu: Các thuật toán hiện hữu đòi hỏi song hàm $f$ phải thỏa mãn tính đơn điệu mạnh hoặc liên tục Lipschitz trên toàn bộ không gian, làm hẹp đáng kể phạm vi ứng dụng thực tế.

Luận án thiết lập hệ thống câu hỏi nghiên cứu và giả thuyết khoa học chặt chẽ:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để xây dựng lược đồ lặp chiếu xấp xỉ song song giải bài toán cân bằng trên tập điểm bất động mà không cần giải chính xác bài toán cân bằng phụ ở từng bước lặp?
  • Câu hỏi nghiên cứu 2 (RQ2): Bằng cách nào có thể kết hợp kỹ thuật song song với phép lặp điểm cố định kiểu Mann và phương pháp dưới đạo hàm để bảo đảm sự hội tụ mạnh trong không gian Hilbert vô hạn chiều khi tập ràng buộc là giao của một họ ánh xạ nửa co?
  • Câu hỏi nghiên cứu 3 (RQ3): Liệu việc tích hợp kỹ thuật ngoại suy quán tính (inertial extrapolation) vào nguyên lý bài toán phụ (Auxiliary Problem Principle - APP) có thể tăng tốc độ hội tụ và giảm thời gian CPU thực thi thuật toán một cách vượt trội hay không?
  • Giả thuyết khoa học 1 (H1): Dãy lặp sinh bởi lược đồ chiếu dưới đạo hàm xấp xỉ song song sẽ hội tụ mạnh về nghiệm duy nhất $x^* \in \Omega$ của bài toán $FEP(\Omega, f)$ dưới điều kiện song hàm $f$ đơn điệu mạnh $\beta$ và tập dưới vi phân xấp xỉ $\partial_{2,\epsilon} f(x,x)$ thỏa mãn điều kiện liên tục Lipschitz kiểu Hausdorff.
  • Giả thuyết khoa học 2 (H2): Kỹ thuật trích chọn song song theo quy tắc cực đại hóa độ lệch $i_0 \in \arg\max { |y_i^k - x^k| : i \in I }$ cho phép xử lý đồng thời họ hữu hạn hoặc vô hạn đếm được các ánh xạ $\beta_i$-nửa co mà vẫn bảo toàn tính ổn định tiệm cận của dãy lặp.
  • Giả thuyết khoa học 3 (H3): Thuật toán dưới đạo hàm quán tính kết hợp nguyên lý bài toán phụ giúp cải thiện từ $30%$ đến trên $60%$ số bước lặp (iterations) và thời gian thực thi (CPU-times) so với các thuật toán chiếu tuần tự cổ điển trên cùng một hệ thống thử nghiệm số.

Khung lý thuyết của luận án tích hợp Giải tích hàm phi tuyến, Giải tích lồi hiện đại, Lý thuyết toán tử đơn điệu, Định lý Ky Fan Minimax, Nguyên lý ánh xạ co Banach và Kỹ thuật ngoại suy Nesterov-Polyak. Đóng góp đột phá của luận án đã được lượng hóa qua 04 công trình khoa học chất lượng cao (gồm 01 bài báo trên tạp chí hệ ISI/SCI, 02 bài báo trên tạp chí hệ SCIE và 01 bài gửi đăng tạp chí SCIE), cùng các thực nghiệm mô phỏng số chi tiết trên môi trường máy tính hiệu năng cao Intel Core i9-9900KS CPU @ 4.00GHz và Intel Core i3-3110M running MATLAB. Phạm vi nghiên cứu bao quát từ không gian Euclide hữu hạn chiều $\mathbb{R}^n$ ($n = 5, 10, \dots$) đến không gian Hilbert vô hạn chiều thực $H$.

Literature Review và Positioning

Khung tổng quan y văn của luận án được cấu trúc qua ba dòng chảy học thuật chủ đạo:

                      ┌────────────────────────────────────────┐
                      │    DÒNG 1: LÝ THUYẾT BÀI TOÁN CÂN BẰNG │
                      │  Nikaido-Isoda (1955), Ky Fan (1972),   │
                      │       Blum & Oettli (1994)             │
                      └──────────────────┬─────────────────────┘
                                         │
                      ┌──────────────────┴─────────────────────┐
                      │     DÒNG 2: LÝ THUYẾT ĐIỂM BẤT ĐỘNG    │
                      │ Banach (1922), Mann (1953), Halpern    │
                      │  (1967), Yamada (2001), Cegielski      │
                      └──────────────────┬─────────────────────┘
                                         │
                      ┌──────────────────┴─────────────────────┐
                      │   DÒNG 3: BÀI TOÁN TỐI ƯU HAI CẤP      │
                      │    Moudafi (2000), Santos (2011),      │
                      │   Anh & Muu (2014), Anh et al. (2021)  │
                      └──────────────────┬─────────────────────┘
                                         │
                                         ▼
                      ┌────────────────────────────────────────┐
                      │      VỊ TRÍ ĐỘT PHÁ CỦA LUẬN ÁN        │
                      │ 4 Thuật toán chiếu song song xấp xỉ &  │
                      │  dưới đạo hàm quán tính trên giao họ   │
                      │      toán tử nửa co trong không gian H │
                      └────────────────────────────────────────┘
  1. Lý thuyết bài toán cân bằng thuần túy: Khởi đầu từ bài toán bất đẳng thức minimax của Ky Fan (1972), Brezis et al. (1972) và định nghĩa tường minh về bài toán cân bằng của Muu & Oettli (1992), Blum & Oettli (1994). Các mở rộng tiếp theo của Bigi et al. (2013), Iusem et al. (2003) chủ yếu tập trung vào các điều kiện tồn tại nghiệm định tính và phương pháp hàm phạt, phương pháp điểm gần kề.
  2. Lý thuyết điểm bất động cho toán tử phi tuyến: Phát triển từ nguyên lý ánh xạ co Banach (1922) sang các mở rộng cho ánh xạ không giãn và ánh xạ nửa co bởi Mann (1953), Halpern (1967), Caristi (1976), Cegielski (2012) và phương pháp hướng giảm lai ghép của Yamada (2001).
  3. Giao thoa giữa bài toán cân bằng và điểm bất động: Các công trình của Moudafi (2000), Combettes (2001), Takahashi (2007), Anh & Muu (2014), Duc et al. (2016), Hai (2014) và gần đây là Anh et al. (2021) nghiên cứu bài toán cân bằng trên tập nghiệm của bất đẳng thức biến phân hoặc điểm bất động.

Các tranh luận khoa học nổi bật trong y văn xoay quanh hai trường phái:

  • Tranh luận 1: Giải chính xác (Exact) đối đầu với Giải xấp xỉ (Inexact): Trường phái cổ điển của Tikhonov và Proximal Point đòi hỏi giải chính xác nghiệm của bài toán cân bằng phụ tại mỗi bước lặp $\min { f(x_k, y) + \frac{1}{2\lambda_k} |y - x_k|^2 }$. Điều này dẫn đến sự bế tắc khi triển khai số. Trường phái xấp xỉ (Santos & Scheimberg, 2011) chủ trương sử dụng $\epsilon$-dưới vi phân nhưng mới chỉ giới hạn cho bài toán cân bằng đơn lẻ với giả thiết paramonotone.
  • Tranh luận 2: Hội tụ yếu (Weak convergence) đối đầu với Hội tụ mạnh (Strong convergence): Trong không gian Hilbert vô hạn chiều, các thuật toán kiểu Mann hoặc gradient chiếu tiêu chuẩn thường chỉ đảm bảo hội tụ yếu. Để đạt được hội tụ mạnh, bắt buộc phải sử dụng các kỹ thuật biến đổi phức tạp như phép chiếu Halpern, phép chiếu thu hẹp CQ của Takahashi, hoặc đưa thêm số hạng đơn điệu mạnh.

So sánh với hai nghiên cứu quốc tế tiêu biểu:

  • So sánh với Santos & Scheimberg (2011 - Tạp chí Optimization): Santos và Scheimberg đề xuất phương pháp chiếu dưới đạo hàm xấp xỉ ($IPSM$) giải bài toán cân bằng $EP(C, f)$. Tuy nhiên, phương pháp $IPSM$ chỉ áp dụng cho miền ràng buộc $C$ lồi đóng thông thường, giải thuật xử lý tuần tự và không thể giải quyết bài toán khi miền ràng buộc bị khống chế bởi tập nghiệm điểm bất động $\Omega = \bigcap_{i \in I} \text{Fix}(S_i)$. Luận án của Nguyễn Văn Hồng đã tổng quát hóa vượt bậc phương pháp $IPSM$ sang lược đồ chiếu song song xấp xỉ đồng thời xử lý họ toán tử nửa co.
  • So sánh với Yamada (2001): Thuật toán hướng giảm lai ghép của Yamada ($HSDM$) chỉ giải bài toán bất đẳng thức biến phân đơn trị $VI(C, F)$ trên tập điểm bất động của một ánh xạ không giãn đơn lẻ $T$. Luận án mở rộng khung phân tích từ toán tử đơn trị $F$ sang song hàm phi tuyến $f(x, y)$, nâng cấp từ ánh xạ không giãn lên họ ánh xạ $\beta_i$-nửa co tổng quát, đồng thời kết hợp kỹ thuật ngoại suy quán tính giúp triệt tiêu hiện tượng trễ dao động trong quá trình lặp.

Đó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 và làm sâu sắc thêm nền tảng lý thuyết cân bằng Blum-Oettli và lý thuyết điểm bất động thông qua các đóng góp cụ thể:

  • Mở rộng phạm vi lớp toán tử ràng buộc: Luận án mở rộng lớp toán tử ràng buộc từ ánh xạ không giãn ($L=0$) sang lớp ánh xạ $\beta_i$-nửa co ($d$-demi-contractive) thỏa mãn: $$|S(x) - x^|^2 \le |x - x^|^2 + d |x - S(x)|^2, \quad \forall x \in C, , x^* \in \text{Fix}(S)$$ và lớp ánh xạ giả co chặt (strictly pseudo-contractive).
  • Hệ thống 04 thuật toán mới với chứng minh hội tụ mạnh:
    1. Thuật toán 2.1 (Chiếu song song xấp xỉ): Kết hợp $\epsilon$-dưới đạo hàm xấp xỉ và hướng giảm lai ghép Yamada.
    2. Thuật toán 2.2 (Dưới đạo hàm song song): Kết hợp toán tử lặp điểm cố định Mann với kỹ thuật song song trên toàn không gian $H$.
    3. Thuật toán 2.3 (Chiếu đạo hàm tăng cường song song): Ứng dụng kỹ thuật đạo hàm tăng cường giải bài toán hai cấp.
    4. Thuật toán 3.1 & 3.2 (Dưới đạo hàm quán tính & Nguyên lý bài toán phụ quán tính song song): Tích hợp bước ngoại suy quán tính $w^k = x^k + \theta_k(x^k - x^{k-1})$ kết hợp nguyên lý bài toán phụ (APP).
  • Thiết lập các Định lý hội tụ mạnh đột phá: Chứng minh toán học chặt chẽ rằng toàn bộ các dãy lặp ${x^k}$ đều hội tụ mạnh trong không gian Hilbert thực $H$ về nghiệm duy nhất $x^*$ (khẳng định tại các Định lý 2.1, Định lý 2.2, Định lý 2.3, Định lý 3.1, Định lý 3.2).

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

Khung phân tích của luận án là sự hợp nhất của 3 trụ cột toán học:

  1. Giải tích lồi và dưới vi phân xấp xỉ: Khai thác khái niệm $\epsilon$-dưới vi phân chéo $\partial_{2,\epsilon} f(x,x) = {w \in H : f(x,y) + \epsilon \ge \langle w, y-x \rangle, \forall y \in C}$ liên tục Lipschitz theo khoảng cách Hausdorff $\rho(A, B) = \max {\sup_{a \in A} \inf_{b \in B} |a-b|, \sup_{b \in B} \inf_{a \in A} |a-b|}$.
  2. Nguyên lý điều khiển bước lặp và tham số đa chiều: Thiết lập các bất đẳng thức ràng buộc tham số tối ưu: $$\tau \in (0, \beta), \quad 0 < \tau_k \le \gamma_k < \min \left{ \frac{2\beta}{L^2}, \frac{2(\beta-\tau)}{L^2-\tau^2}, \frac{1}{\tau} \right}, \quad 0 < a \le \alpha_{k,i} < \min \left{ \frac{1-\beta_i}{2} : i \in I \right}$$
  3. Quy tắc lựa chọn song song cực đại (Parallel argmax selection rule): $$y^k := y_{i_0}^k \quad \text{với} \quad i_0 \in \arg\max { |y_i^k - x^k| : i \in I }$$ Quy tắc này giải phóng thuật toán khỏi sự phụ thuộc vào cấu trúc thứ tự tuần tự, cho phép phân bổ tính toán song song độc lập trên các nhân vi xử lý.
                      KHUNG PHÂN TÍCH TỔNG QUÁT CỦA LUẬN ÁN
  ┌─────────────────────────────────────────────────────────────────────────┐
  │                         ĐIỂM XUẤT PHÁT: x^k                             │
  └────────────────────────────────────┬────────────────────────────────────┘
                                       │
                    ┌──────────────────┴──────────────────┐
                    ▼                                     ▼
      ┌───────────────────────────┐         ┌───────────────────────────┐
      │  NHÁNH 1: CHIẾU SONG SONG │         │  NHÁNH 2: NGOẠI SUY QUÁN  │
      │       (Parallel Branch)   │         │     TÍNH (Inertial Branch)│
      │   y_i^k = (1-α)x^k + αS_i │         │ w^k = x^k + θ(x^k - x^{k-1})│
      │ i_0 = argmax ||y_i^k-x^k||│         │                           │
      └─────────────┬─────────────┘         └─────────────┬─────────────┘
                    │                                     │
                    └──────────────────┬──────────────────┘
                                       │
                                       ▼
      ┌─────────────────────────────────────────────────────────────────┐
      │ BƯỚC CẬP NHẬT DƯỚI ĐẠO HÀM XẤP XỈ / NGUYÊN LÝ BÀI TOÁN PHỤ (APP)│
      │     u^k ∈ ∂_{2,τ_k} f(y^k, y^k),  x^{k+1} = P_C(y^k - γ_k u^k)   │
      └────────────────────────────────┬────────────────────────────────┘
                                       │
                                       ▼
      ┌─────────────────────────────────────────────────────────────────┐
      │     HỘI TỤ MẠNH: ||x^k - x*|| → 0 trong không gian Hilbert H    │
      └─────────────────────────────────────────────────────────────────┘

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

Thiết kế nghiên cứu

  • Hệ hình nghiên cứu (Research Paradigm): Luận án vận hành hoàn toàn trên hệ hình thực chứng diễn dịch (Deductive Positivism) của Toán học thuần túy và Toán ứng dụng, sử dụng phương pháp chứng minh tiên đề giải tích chặt chẽ kết hợp kiểm chứng mô phỏng số (numerical experimentation).
  • Không gian trạng thái: Không gian Hilbert thực $H$ tổng quát trang bị tích vô hướng $\langle \cdot, \cdot \rangle$ và chuẩn cảm sinh $|u| = \sqrt{\langle u, u \rangle}$, bao quát cả không gian Euclide hữu hạn chiều $\mathbb{R}^n$ và không gian hàm vô hạn chiều $L^2, \ell^2$.
  • Cơ chế hội tụ: Sử dụng bổ đề kỹ thuật Maingé (2008) cho các dãy số không đơn điệu nhằm triệt tiêu đòi hỏi về tính đơn điệu của dãy khoảng cách $|x^k - x^*|$, kết hợp tính chất tựa không giãn của toán tử $S_\alpha = (1-\alpha)I + \alpha S$.

Quy trình nghiên cứu rigorous

Quy trình giải tích và tính toán được chuẩn hóa qua 4 giai đoạn logic:

┌─────────────────────────────────────────────────────────────────────────────┐
│ GIAI ĐOẠN 1: MÔ HÌNH HÓA TOÁN HỌC                                           │
│ Biểu diễn bài toán thực tế dưới dạng FEP(Ω, f) với Ω = ∩ Fix(S_i)           │
└──────────────────────────────────────┬──────────────────────────────────────┘
                                       │
                                       ▼
┌─────────────────────────────────────────────────────────────────────────────┐
│ GIAI ĐOẠN 2: THIẾT KẾ THUẬT TOÁN & BỘ ĐIỀU KHIỂN BƯỚC LẶP                    │
│ Xây dựng công thức bước lặp, ràng buộc tham số τ, γ_k, α_{k,i}, θ_k         │
└──────────────────────────────────────┬──────────────────────────────────────┘
                                       │
                                       ▼
┌─────────────────────────────────────────────────────────────────────────────┐
│ GIAI ĐOẠN 3: CHỨNG MINH HỘI TỤ ĐỊNH LÝ TOÁN HỌC                              │
│ Chứng minh tính bị chặn → Điểm tụ yếu x̄ ∈ Ω → Khẳng định hội tụ mạnh ||x^k-x*||→ 0│
└──────────────────────────────────────┬──────────────────────────────────────┘
                                       │
                                       ▼
┌─────────────────────────────────────────────────────────────────────────────┐
│ GIAI ĐOẠN 4: THỰC NGHIỆM SỐ & BENCHMARK ĐỐI CHỨNG                            │
│ Chạy code MATLAB giải bài toán kích thước lớn, phân tích Iterations & CPU-times │
└─────────────────────────────────────────────────────────────────────────────┘

Tính hợp lệ nội tại (Internal Validity) được kiểm soát thông qua tính đúng đắn tuyệt đối của các phép biến đổi đại số tích vô hướng: $$|u \pm v|^2 = |u|^2 \pm 2\langle u, v \rangle + |v|^2, \quad |u+v|^2 \le |u|^2 + 2\langle v, u+v \rangle$$ và các tính chất hình học của phép chiếu metric $P_C(x) = \arg\min_{y \in C} |y - x|^2$: $$\langle x - P_C(x), P_C(x) - y \rangle \ge 0, \quad \forall y \in C$$

Data và phân tích số liệu thực nghiệm

Dữ liệu mô phỏng được thiết kế theo chuẩn benchmark quốc tế cho bài toán cân bằng phi tuyến trong $\mathbb{R}^5$:

  • Song hàm cân bằng $f(x, y) = \langle F(x) + Qy + q, y - x \rangle$
  • Ma trận $Q = AA^T + B + D$ được xây dựng từ:
    • Ma trận vuông $A \in \mathbb{R}^{5 \times 5}$: $$A = \begin{pmatrix} -2 & 1 & 0 & 1 & -1 \ 1 & 2 & 1 & 0 & 2 \ 0 & 1 & 3 & 1 & 2 \ 0 & 1 & 3 & 1 & 0 \ 2 & 0 & 1 & -1 & 3 \end{pmatrix}$$
    • Ma trận phản xứng $B \in \mathbb{R}^{5 \times 5}$ ($B = -B^T$): $$B = \begin{pmatrix} 0 & 1 & 2 & 1 & -1 \ -1 & 3 & 2 & 0 & 2 \ -2 & -2 & 1 & 1 & -3 \ -1 & 0 & -1 & 1 & 0 \ 1 & -2 & 3 & 0 & 2 \end{pmatrix}$$
    • Ma trận đường chéo dương $D = \text{diag}(5, 3, 12, 15, 22)$
    • Vectơ hằng số $q = (2, 3, -4, 1, 5)^T$
  • Các chỉ số ma trận định lượng: Giá trị riêng nhỏ nhất $\lambda_{\min}(Q) = 10.2313$, chuẩn ma trận $|Q| = 58.0016$.
  • Ánh xạ phi tuyến $F(x) = (\eta x_1 + \eta x_2 + \sin x_1, -\eta x_1 + \eta x_2 + \sin x_2, (\eta-1)x_3, (\eta-1)x_4, (\eta-1)x_5)^T$ với $\eta > 1 + |Q|$. Hằng số liên tục Lipschitz của $F$ được chứng minh tường minh là $L = \sqrt{2(2\eta^2 + 2\eta + 1)}$ và $F$ là $(\eta-1)$-đơn điệu mạnh. Do đó song hàm $f$ đơn điệu mạnh với hằng số $\beta = \eta + 9.2313$.
  • Hai ánh xạ nửa co phi tuyến phức tạp: $$S_1(x) = \left(x_1, \frac{1}{3}\sin x_2, \frac{1}{3}x_3, \frac{1}{2}x_4, \sin^3 x_5\right)^T, \quad S_2(x) = \left(\frac{1}{4}x_1, x_2, \sin x_3, \sin^2 x_4, x_5\right)^T$$
  • Môi trường phần mềm và phần cứng: Chương trình thực thi trên MATLAB R2016a / R2013a; phần cứng máy tính CPU Intel Core i9-9900KS @ 4.00GHz và Laptop Intel Core i3-3110M CPU @ 2.40GHz 4GB RAM. Phép chiếu lên tập đa giác lồi $C$ được giải tối ưu bằng hàm solver fmincon. Tiêu chuẩn dừng thuật toán: sai số tuyệt đối $\max { |y^k - x^k|, |x^{k+1} - x^k| } \le \epsilon = 10^{-4}$ hoặc $10^{-6}$.

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

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

  1. Xác lập sự hội tụ mạnh không điều kiện compact: Luận án chứng minh thành công sự hội tụ mạnh của toàn bộ 4 thuật toán trong không gian Hilbert thực vô hạn chiều mà không cần giả thiết tập ràng buộc $C$ hay tập điểm bất động $\Omega$ phải compact. Đây là bước tiến vượt bậc so với định lý Ky Fan cổ điển vốn phụ thuộc chặt vào tính compact.
  2. Hiệu quả vượt trội của quy tắc trích chọn song song cực đại: Cơ chế $i_0 \in \arg\max_{i \in I} |y_i^k - x^k|$ giúp giảm thiểu số lượng phép chiếu đơn lẻ, cho phép thuật toán đạt độ chính xác $\epsilon = 10^{-6}$ chỉ sau vài chục bước lặp (so với hàng trăm bước lặp của các thuật toán chiếu luân phiên tuần tự).
  3. Hiện tượng gia tốc quán tính (Inertial Acceleration Phenomenon): Việc bổ sung số hạng quán tính $\theta_k(x^k - x^{k-1})$ tạo ra động lực vượt qua các thung lũng hẹp của hàm mục tiêu, giúp tốc độ giảm sai số $|x^k - x^*|$ đạt mức độ dốc cao hơn rõ rệt mà không gây mất ổn định nghiệm.
  4. Triệt tiêu nút thắt tính toán bài toán phụ: Bằng cách sử dụng phần tử $\epsilon$-dưới đạo hàm $u^k \in \partial_{2,\tau_k} f(y^k, y^k)$, chi phí cho mỗi bước lặp chuyển thành một phép chiếu trực tiếp lên tập lồi thay vì phải giải phương trình phi tuyến phức tạp, tiết kiệm hơn $50%$ thời gian CPU trên mỗi vòng lặp.

Implications đa chiều

                       CÁC HỆ GIÁ TRỊ VÀ ỨNG DỤNG ĐA CHIỀU
  ┌───────────────────────────────────┬───────────────────────────────────┐
  │      Ý NGHĨA LÝ THUYẾT & PHƯƠNG PHÁP│       Ý NGHĨA THỰC TIỄN & CHÍNH SÁCH│
  ├───────────────────────────────────┼───────────────────────────────────┤
  │• Thống nhất lý thuyết tối ưu hai  │• Mô hình cân bằng thị trường điện │
  │  cấp, biến phân và điểm bất động. │  thông minh đa tác tử (Smart Grid).│
  │• Cung cấp khung thuật toán song   │• Phân bổ lưu lượng giao thông đô  │
  │  song hóa cho tính toán phân tán. │  thị cực đại hóa phúc lợi xã hội. │
  │• Xóa bỏ đòi hỏi giải chính xác    │• Giải bài toán Nash đa hãng sản   │
  │  bài toán phụ cấp dưới.           │  xuất với hàm Nikaido-Isoda.      │
  └───────────────────────────────────┴───────────────────────────────────┘
  • Đóng góp lý thuyết: Mở rộng biên giới của Lý thuyết toán tử đơn điệu và Giải tích lồi, cung cấp công cụ mạnh mẽ để nghiên cứu các bài toán đặt không chỉnh (ill-posed problems) trong kinh tế học toán lượng và tối ưu hóa đa mục tiêu.
  • Đổi mới phương pháp luận: Thiết lập quy trình chuẩn tắc kết hợp giữa giải tích xấp xỉ dưới vi phân và ngoại suy quán tính có thể chuyển giao trực tiếp cho các lớp bài toán khác như bài toán nứt gãy cơ học, bài toán khôi phục tín hiệu/hình ảnh số và học máy phân tán (Distributed Machine Learning).
  • Ứng dụng thực tiễn: Mô hình hóa hoàn hảo bài toán cân bằng thị trường cạnh tranh Nash của $N$ hãng sản xuất thông qua hàm Nikaido-Isoda $\Phi(x, y) = \sum_{i=1}^N [f_i(x_1, \dots, x_i, \dots, x_N) - f_i(x_1, \dots, y_i, \dots, x_N)]$, hỗ trợ các cơ quan điều tiết thị trường thiết lập chính sách giá trần, thuế carbon và kiểm soát độc quyền dựa trên bằng chứng toán học lượng hóa.

Limitations và Future Research

Nhìn nhận một cách khách quan và khoa học, luận án tồn tại các giới hạn biên sau:

  1. Giả thiết về tính đơn điệu mạnh của song hàm: Sự hội tụ mạnh của các thuật toán trong luận án vẫn dựa trên giả thiết song hàm $f$ là $\beta$-đơn điệu mạnh hoặc giả đơn điệu mạnh. Đối với các bài toán cân bằng tổng quát chỉ thỏa mãn tính đơn điệu yếu, tựa đơn điệu (quasimonotone) hoặc không đơn điệu, thuật toán chưa áp dụng trực tiếp được.
  2. Yêu cầu ước lượng trước các hằng số cấu trúc: Việc xác định kích thước bước lặp $\gamma_k$ đòi hỏi phải biết trước hoặc ước lượng cận trên của hằng số Lipschitz $L$ và hằng số đơn điệu $\beta$. Trong thực tế tính toán quy mô siêu lớn, việc tính toán chuẩn ma trận hoặc hằng số Lipschitz có thể tốn kém chi phí.
  3. Giới hạn trong không gian Hilbert: Mặc dù không gian Hilbert vô hạn chiều là rất tổng quát, nhưng nhiều bài toán điều khiển tối ưu và phương trình đạo hàm riêng tự nhiên thuộc về các không gian Banach phản xạ không có tích vô hướng (như không gian Sobolev $W^{1,p}$ với $p \ne 2$).

Chương trình nghiên cứu tương lai (Future Research Agenda):

  • Hướng 1: Phát triển kỹ thuật bước nhảy tự điều chỉnh (Self-adaptive / Armijo-like linesearch rule) để loại bỏ hoàn toàn sự phụ thuộc vào việc biết trước hằng số Lipschitz $L$.
  • Hướng 2: Mở rộng thuật toán cho lớp bài toán cân bằng với song hàm giả đơn điệu (pseudomonotone) hoặc phi đơn điệu bằng cách tích hợp kỹ thuật siêu mặt phẳng cắt (hyperplane cutting methods) và phép lặp Tseng.
  • Hướng 3: Nghiên cứu phiên bản ngẫu nhiên (Stochastic Parallel Inertial Subgradient Methods) để giải quyết các bài toán tối ưu dữ liệu lớn và học máy nơi gradient chỉ có thể ước lượng qua mẫu ngẫu nhiên.
  • Hướng 4: Mở rộng hệ lý thuyết từ không gian Hilbert sang không gian Banach phản xạ thực sử dụng khoảng cách Bregman (Bregman distance) và toán tử ánh xạ chuẩn tắc (normalized duality mapping).

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

  • Tác động học thuật quốc tế: Các kết quả của luận án được công bố trên 03 tạp chí uy tín thuộc danh mục Web of Science (SCI/SCIE) và báo cáo tại các diễn đàn học thuật đỉnh cao như Hội nghị Toán học Toàn quốc lần thứ X (VMC 2023 - Đà Nẵng), Hội nghị Quốc tế ISAS 2022 (ĐH Bách khoa TP.HCM), Hội thảo Tối ưu và Tính toán Khoa học lần thứ 21 (Ba Vì, 2023) và Viện Nghiên cứu Cao cấp về Toán (VIASM, 2021). Ước tính các thuật toán này sẽ thu hút chỉ số trích dẫn cao từ cộng đồng nghiên cứu tối ưu hóa phi tuyến toàn cầu.
  • Chuyển đổi công nghệ và tính toán công nghiệp: Cung cấp thuật toán tối ưu hóa lõi có khả năng tích hợp vào các phần mềm mô phỏng hệ thống năng lượng, điều độ lưới điện thông minh, và phần mềm tối ưu hóa mạng viễn thông thế hệ mới (5G/6G Resource Allocation).
  • Lợi ích xã hội và hoạch định chính sách: Cung cấp công cụ toán học định lượng hỗ trợ các nhà phân tích chính sách kinh tế vĩ mô mô phỏng hành vi thị trường, dự báo điểm cân bằng cung-cầu và đánh giá tác động của các chính sách thuế, bảo vệ môi trường.

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

  • Nghiên cứu sinh và Giảng viên ngành Toán ứng dụng: Tiếp cận hệ thống chứng minh toán học chuẩn mực, các kỹ thuật giải tích lồi hiện đại và kỹ năng thiết kế thuật toán song song.
  • Nhà nghiên cứu Vận trù học (Operations Research): Ứng dụng các thuật toán giải bài toán bù phi tuyến, bất đẳng thức biến phân hai cấp trong logistics và chuỗi cung ứng.
  • Kỹ sư R&D trong Công nghệ thông tin & Viễn thông: Sử dụng các thuật toán chiếu song song để tối ưu hóa băng thông, quản lý công suất phát và định tuyến mạng lưới.
  • Chuyên gia phân tích kinh tế & hoạch định chính sách: Áp dụng mô hình cân bằng Nash-Nikaido-Isoda để dự báo cân bằng thị trường cạnh tranh không hoàn hảo.

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 học thuyết nào?

Đóng góp lý thuyết độc đáo nhất là việc thiết lập thành công định lý hội tụ mạnh cho lớp bài toán cân bằng hai cấp trên tập giao các điểm bất động $\Omega = \bigcap_{i \in I} \text{Fix}(S_i)$ của họ toán tử $\beta_i$-nửa co trong không gian Hilbert vô hạn chiều. Luận án mở rộng trực tiếp khung lý thuyết cân bằng Blum-Oettli (1994) và nguyên lý xấp xỉ dưới vi phân của Santos & Scheimberg (2011), khắc phục triệt để hạn chế chỉ hội tụ yếu của các lược đồ Mann kinh điển bằng cách tích hợp điều kiện bước nhảy $\gamma_k < \min \left{ \frac{2\beta}{L^2}, \frac{2(\beta-\tau)}{L^2-\tau^2}, \frac{1}{\tau} \right}$ cùng quy tắc chọn cực đại khoảng cách song song.

2. Đổi mới phương pháp luận thể hiện như thế nào khi so sánh với ít nhất 2 công trình quốc tế?

So với nghiên cứu của Santos & Scheimberg (2011), luận án đã chuyển hóa phương pháp chiếu dưới đạo hàm xấp xỉ tuần tự trên miền lồi $C$ thành phương pháp chiếu song song xấp xỉ đa toán tử trên giao tập điểm bất động. So với công trình của Yamada (2001) vốn chỉ áp dụng hướng giảm lai ghép cho bất đẳng thức biến phân đơn trị $VI(C, F)$ với một ánh xạ không giãn duy nhất, luận án đã giải quyết thành công song hàm cân bằng phi tuyến không khả vi $f(x, y)$ trên họ vô hạn các ánh xạ nửa co, đồng thời tích hợp số hạng gia tốc quán tính để tối ưu hóa thời gian thực thi.

3. Phát hiện nào gây bất ngờ nhất về mặt toán học và thực nghiệm số?

Phát hiện bất ngờ nhất là việc quy tắc lựa chọn song song $i_0 \in \arg\max_{i \in I} |y_i^k - x^k|$ không chỉ duy trì tính hội tụ toán học chặt chẽ mà còn đóng vai trò như một bộ lọc tích cực: chỉ cần cập nhật theo hướng của toán tử có độ lệch lớn nhất tại mỗi bước lặp, toàn bộ hệ thống các toán tử ràng buộc còn lại đều tự động co về tập điểm bất động chung $\Omega$ với tốc độ tiệm cận nhanh hơn việc tính toán trung bình trọng số toàn bộ các toán tử.

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

Luận án cung cấp đầy đủ và chi tiết toàn bộ tham số cấu trúc, ma trận thực nghiệm ($A, B, D \in \mathbb{R}^{5 \times 5}$ tường minh), vectơ $q$, hằng số liên tục Lipschitz $L = \sqrt{2(2\eta^2 + 2\eta + 1)}$, điều kiện tham số $\tau, \gamma_k, \alpha_{k,i}$, phần mềm MATLAB R2016a/R2013a, hàm solver fmincon, phần cứng CPU Intel Core i9-9900KS @ 4.00GHz / i3-3110M @ 2.40GHz và tiêu chuẩn dừng sai số $\epsilon \le 10^{-6}$, đảm bảo khả năng tái lập kết quả thực nghiệm chính xác $100%$.

5. Chương trình nghị sự nghiên cứu 10 năm (10-year research agenda) được vạch ra như thế nào?

Chương trình nghị sự 10 năm định hướng tập trung vào: (1) Xây dựng lý thuyết cân bằng trên đa tạp Riemannian; (2) Phát triển các giải thuật học sâu tăng cường (Reinforcement Learning) dựa trên nguyên lý cân bằng Nash đa tác tử; (3) Tích hợp thuật toán tối ưu quán tính vào các hệ thống tính toán lượng tử (Quantum Optimization Algorithms) để giải quyết các bài toán cân bằng quy mô hàng triệu biến số.

Kết luận

Luận án tiến sĩ của tác giả Nguyễn Văn Hồng là một công trình nghiên cứu khoa học xuất sắc, mẫu mực và hoàn chỉnh trong chuyên ngành Toán ứng dụng, thể hiện qua 5 kết luận cốt lõi:

  1. Hoàn thiện khung lý thuyết mới: Xây dựng thành công 04 thuật toán chiếu mở rộng và dưới đạo hàm quán tính giải quyết trọn vẹn bài toán cân bằng trên tập giao điểm bất động của họ toán tử nửa co.
  2. Chứng minh hội tụ mạnh toàn cục: Thiết lập các định lý hội tụ mạnh nghiêm ngặt trong không gian Hilbert vô hạn chiều, vượt qua rào cản hội tụ yếu của các phương pháp lặp truyền thống.
  3. Đột phá về hiệu năng tính toán: Giảm thiểu đáng kể chi phí tính toán thông qua kỹ thuật $\epsilon$-dưới đạo hàm và gia tốc quán tính, loại bỏ hoàn toàn yêu cầu giải chính xác bài toán cân bằng phụ.
  4. Kiểm chứng thực nghiệm vững chắc: Các mô phỏng số chi tiết trên phần mềm MATLAB với bộ dữ liệu ma trận đầy đủ đã chứng minh tính ưu việt vượt trội của thuật toán đề xuất về cả số bước lặp (Iterations) và thời gian CPU.
  5. Mở ra các hướng nghiên cứu liên ngành: Đặt nền móng vững chắc cho việc ứng dụng lý thuyết cân bằng hiện đại vào các bài toán kinh tế thị trường Nash, mạng lưới giao thông thông minh và tối ưu hóa hệ thống phân tán trong kỷ nguyên số.