Tổng quan về luận án

Nghiên cứu của tiến sĩ Trần Việt Anh (2018) với đề tài "Phương pháp giải bất đẳng thức biến phân trên tập nghiệm của bài toán chấp nhận tách suy rộng" thuộc chuyên ngành Toán giải tích (Mã số: 62460102), dưới sự hướng dẫn khoa học của GS. Lê Dũng Mưu và PGS. Nguyễn Hữu Điển tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội, là một công trình khoa học tiên phong thuộc lĩnh vực giải tích phi tuyến và tối ưu hóa hiện đại.

Bối cảnh khoa học của đề tài bắt nguồn từ sự phát triển sâu rộng của lý thuyết bất đẳng thức biến phân (Variational Inequality Problem - VIP), được khởi xướng bởi Philip Hartman và Guido Stampacchia (1966), sau đó được hệ thống hóa toàn diện trong không gian vô hạn chiều bởi David Kinderlehrer và Guido Stampacchia (1980) cùng Claudio Baiocchi và Antonio Capelo (1984). Song song với đó, bài toán chấp nhận tách (Split Feasibility Problem - SFP) do Yair Censor và Tommy Elfving (1994) đề xuất, cùng thuật toán CQ nổi tiếng của Charles Byrne (2002) và các mở rộng trong không gian Hilbert vô hạn chiều của Hong-Kun Xu (2010), đã trở thành công cụ toán học nền tảng mô hình hóa các bài toán thực tế như liệu pháp xạ trị điều chỉnh cường độ (Intensity-Modulated Radiation Therapy - IMRT), xử lý tín hiệu số và bài toán khôi phục ảnh y khoa.

flowchart TD
    A["Bài toán chấp nhận tách suy rộng (GSFP)"] --> B["Bài toán điểm bất động tách (SFPP)"]
    A --> C["Bài toán BĐT biến phân tách (SVIP / BSVIP)"]
    A --> D["Bài toán cân bằng tách (SEP)"]
    A --> E["Chấp nhận tách đa tập hợp (MSSFP)"]
    
    B --> F["VIP trên tập nghiệm SFPP<br/>(Krasnoselskii-Mann + Chiếu)"]
    C --> G["Bất đẳng thức biến phân hai cấp<br/>(Dưới đạo hàm tăng cường)"]
    D --> H["Nghiệm chuẩn nhỏ nhất của SEP<br/>(Đạo hàm tăng cường liên không gian)"]
    E --> I["VIP trên tập nghiệm MSSFP<br/>(Thuật toán song song + CQ)"]
    
    F --> J["Hội tụ mạnh trong không gian Hilbert vô hạn chiều"]
    G --> J
    H --> J
    I --> J

Khoảng trống nghiên cứu then chốt (Research Gap) được luận án định vị là: thuật toán CQ và các phương pháp giải SFP truyền thống đòi hỏi phải tính toán tường minh phép chiếu trực giao lên các tập lồi ràng buộc $C$ và $Q$. Tuy nhiên, trong thực tế toán học, khi $C$ và $Q$ xuất hiện dưới dạng tập ẩn phức tạp—chẳng hạn như tập điểm bất động của ánh xạ không giãn $\text{Fix}(T)$, tập nghiệm của bài toán bất đẳng thức biến phân $\text{Sol}(C, F)$, hay tập nghiệm của bài toán cân bằng $\text{Sol}(C, f)$—thì phép chiếu $P_C$ và $P_Q$ không thể xác định được một cách trực tiếp. Các dạng bài toán này tạo thành lớp bài toán chấp nhận tách suy rộng (Generalized Split Feasibility Problem - GSFP). Hơn nữa, hầu hết các phương pháp trước đó chỉ đạt được sự hội tụ yếu hoặc đòi hỏi các điều kiện ngặt nghèo như tính đơn điệu mạnh hay đơn điệu mạnh ngược của toán tử, gây bế tắc khi xử lý các toán tử chỉ thỏa mãn tính giả đơn điệu (pseudomonotone) và liên tục Lipschitz.

Luận án thiết lập 4 câu hỏi nghiên cứu và 4 giả thuyết khoa học tương ứng:

  • Câu hỏi nghiên cứu 1 ($RQ_1$): Làm thế nào để thiết lập thuật toán lặp hội tụ mạnh giải VIP đơn điệu mạnh trên tập nghiệm của bài toán điểm bất động tách (SFPP) mà không cần tính toán phép chiếu trực tiếp lên tập nghiệm ẩn?
  • Câu hỏi nghiên cứu 2 ($RQ_2$): Có thể xây dựng cơ chế chiếu dưới đạo hàm tăng cường giải bài toán bất đẳng thức biến phân hai cấp (BVIP) và bất đẳng thức biến phân tách hai cấp (BSVIP) khi các toán tử cấp dưới chỉ giả đơn điệu và liên tục Lipschitz hay không?
  • Câu hỏi nghiên cứu 3 ($RQ_3$): Làm thế nào để giải bài toán VIP trên tập nghiệm của bài toán chấp nhận tách đa tập hợp (MSSFP) với cấu trúc tính toán song song tối ưu?
  • Câu hỏi nghiên cứu 4 ($RQ_4$): Thuật toán nào cho phép tìm nghiệm có chuẩn nhỏ nhất (minimum-norm solution) của bài toán cân bằng tách (SEP) với các song hàm chỉ thỏa mãn tính giả đơn điệu và điều kiện kiểu Lipschitz?

Tương ứng là các giả thuyết khoa học:

  • Giả thuyết 1 ($H_1$): Sự kết hợp giữa phương pháp chiếu gradient hiệu chỉnh và kỹ thuật lặp Krasnoselskii-Mann sẽ bảo toàn tính hội tụ mạnh của dãy lặp về nghiệm duy nhất của VIP trên tập nghiệm SFPP trong không gian Hilbert vô hạn chiều.
  • Giả thuyết 2 ($H_2$): Phép chiếu lên nửa không gian biến thiên (half-space) trong kỹ thuật dưới đạo hàm tăng cường sẽ triệt tiêu đòi hỏi về tính đơn điệu mạnh của toán tử cấp dưới trong bài toán BSVIP.
  • Giả thuyết 3 ($H_3$): Cấu trúc lặp song song kết hợp toán tử liên hợp $A^*$ sẽ khắc phục hoàn toàn hiện tượng nghẽn tính toán khi số lượng tập lồi $M, N$ trong MSSFP tăng cao.
  • Giả thuyết 4 ($H_4$): Toán tử đơn giá $T_f$ xây dựng từ song hàm cân bằng giả đơn điệu $f$ thỏa mãn tính tựa không giãn và nguyên lý bán đóng (demiclosedness principle), cho phép xây dựng thuật toán tìm nghiệm chuẩn nhỏ nhất của SEP.

Khung lý thuyết của công trình tích hợp lý thuyết toán tử đơn điệu, giải tích lồi hiện đại và giải tích phi tuyến. Đóng góp đột phá của luận án đã được công bố trên 6 bài báo khoa học thuộc danh mục ISI/Scopus quốc tế uy tín, được kiểm chứng thực nghiệm số trên phần mềm MATLAB R2012a (máy tính Intel Core i3-3217U CPU @ 1.80GHz, 2 GB RAM) với sai số dừng $\varepsilon \le 10^{-6}$.

Literature Review và Positioning

Lịch sử phát triển của lý thuyết tối ưu hóa và giải tích phi tuyến ghi nhận các nhánh nghiên cứu chính:

                  ┌───────────────────────────────────────────────────────────┐
                  │                 LÝ THUYẾT BẤT ĐẲNG THỨC BIẾN PHÂN          │
                  │   Hartman & Stampacchia (1966) -> Kinderlehrer (1980)    │
                  └─────────────────────────────┬─────────────────────────────┘
                                                │
         ┌──────────────────────────────────────┴──────────────────────────────────────┐
         ▼                                                                             ▼
┌─────────────────────────────────┐                                   ┌─────────────────────────────────┐
│     PHƯƠNG PHÁP ĐẠO HÀM TĂNG CƯỜNG    │                                   │     BÀI TOÁN CHẤP NHẬN TÁCH (SFP)   │
│ • Levitin & Polyak (1966)       │                                   │ • Censor & Elfving (1994)       │
│ • Korpelevich (1976)            │                                   │ • Byrne (2002): Thuật toán CQ   │
│ • Censor et al. (2011): Subgrad │                                   │ • Xu (2010): Vô hạn chiều       │
│ • Kraikaew & Saejung (2014)     │                                   │ • Ceng et al. (2012): Min-norm  │
└────────────────┬────────────────┘                                   └────────────────┬────────────────┘
                 │                                                                     │
                 └──────────────────────────────────┬──────────────────────────────────┘
                                                    │
                                                    ▼
                  ┌───────────────────────────────────────────────────────────┐
                  │             ĐỊNH VỊ ĐỘT PHÁ CỦA LUẬN ÁN (2018)             │
                  │  Giải VIP/BSVIP/SEP/MSSFP trên tập nghiệm ẩn suy rộng      │
                  │  Toán tử giả đơn điệu + Lipschitz, Hội tụ MẠNH vô hạn chiều│
                  └───────────────────────────────────────────────────────────┘
  1. Dòng nghiên cứu về VIP và phương pháp lặp chiếu: Khởi đầu từ phương pháp gradient chiếu của Alan A. Levitin và Boris T. Polyak (1966), thuật toán đòi hỏi toán tử $F$ phải đơn điệu mạnh hoặc đơn điệu mạnh ngược. Để nới lỏng giả thiết, Galina M. Korpelevich (1976) đề xuất phương pháp đạo hàm tăng cường (Extragradient method) cho toán tử đơn điệu và liên tục Lipschitz nhưng đòi hỏi 2 phép chiếu lên tập $C$ trong mỗi bước lặp. Năm 2011, Yair Censor cùng cộng sự đã tạo bước ngoặt khi đề xuất phương pháp dưới đạo hàm tăng cường (Subgradient Extragradient method), thay thế phép chiếu thứ hai bằng phép chiếu lên nửa không gian $T_k = {\omega \in H : \langle x^k - \tau F(x^k) - y^k, \omega - y^k \rangle \le 0}$, giảm thiểu chi phí tính toán nhưng chỉ đạt hội tụ yếu. Đến năm 2014, Rapeepan Kraikaew và Satit Saejung kết hợp phương pháp này với kỹ thuật Halpern để thu được dãy hội tụ mạnh.
  2. Dòng nghiên cứu về bài toán chấp nhận tách (SFP): Sau đề xuất nền tảng của Yair Censor và Tommy Elfving (1994), Charles Byrne (2002) phát triển thuật toán CQ: $$x^{k+1} = P_C(x^k + \gamma A^T (P_Q - I)Ax^k)$$ với $\gamma \in (0, 2/L)$. Hong-Kun Xu (2010) mở rộng thuật toán CQ sang không gian Hilbert vô hạn chiều dùng toán tử liên hợp $A^*$. Tuy nhiên, hạn chế cốt tử của thuật toán CQ là đòi hỏi tính được tường minh $P_C$ và $P_Q$.
  3. Dòng nghiên cứu về bài toán cân bằng (EP) và bài toán điểm bất động: Bất đẳng thức cân bằng được Hukukane Nikaido và Kazuo Isoda đề xuất năm 1955 nhằm tổng quát hóa cân bằng Nash. Ky Fan (1972) phát triển định lý minimax tồn tại nghiệm. Lê Dũng Mưu và Werner Oettli (1992), Eugen Blum và Werner Oettli (1994) đã định hình thuật ngữ "Bài toán cân bằng" (Equilibrium Problem - EP). Patrick Louis Combettes và Sever Adrian Hirstoaga (2005) đã mở đường kết nối bài toán cân bằng đơn điệu với lý thuyết điểm bất động của ánh xạ không giãn thông qua toán tử giải quy $T_r^f$.
  4. Dòng nghiên cứu về bài toán tìm nghiệm chuẩn nhỏ nhất: Lu-Chuan Ceng cùng các đồng nghiệp (2012) đề xuất thuật toán tìm nghiệm có chuẩn nhỏ nhất của bài toán chấp nhận tách thông qua toán tử hiệu chỉnh $\nabla f_{\alpha_k} := A^*(I - P_Q)A + \alpha_k I$. Tại Việt Nam, các hướng nghiên cứu của Phạm Kỳ Anh, Nguyễn Bường (2014), Phan Quốc Khánh, Lê Dũng Mưu và Nguyễn Thị Thu Thủy đã đạt được nhiều kết quả quan trọng về các bài toán biến phân đa cấp.

Trong y văn tồn tại hai luồng tranh luận học thuật lớn:

  • Luồng quan điểm thứ nhất (đại diện bởi các phương pháp gradient cổ điển) cho rằng để đạt được hội tụ mạnh bắt buộc toán tử phải thỏa mãn tính đơn điệu mạnh hoặc đơn điệu mạnh ngược.
  • Luồng quan điểm thứ hai (tiêu biểu bởi Censor et al., Kraikaew & Saejung) tìm cách nới lỏng sang toán tử giả đơn điệu nhưng vấp phải trở ngại lớn về độ phức tạp khi giải các bài toán biến phân hai cấp với ràng buộc tách (BSVIP) và bài toán cân bằng tách (SEP).

Luận án của Trần Việt Anh đã định vị chính xác điểm nghẽn này: giải quyết bài toán VIP hai cấp mà bài toán cấp trên là đơn điệu mạnh, bài toán cấp dưới là bài toán biến phân tách hoặc bài toán cân bằng tách với toán tử chỉ cần giả đơn điệu và liên tục Lipschitz, thiết lập thuật toán đầu tiên trên thế giới giải quyết trọn vẹn lớp bài toán này với chứng minh hội tụ mạnh trong không gian Hilbert vô hạn chiều.

Tiêu chí so sánh Thuật toán CQ của Byrne (2002) / Xu (2010) Thuật toán Min-norm của Lu-Chuan Ceng et al. (2012) Thuật toán BSVIP & Min-norm SEP của Trần Việt Anh (2018)
Không gian khảo sát Hilbert hữu hạn / vô hạn chiều Hilbert thực Không gian Hilbert thực vô hạn chiều $H_1, H_2$
Bản chất tập ràng buộc Tập lồi đóng $C, Q$ tường minh Tập lồi đóng $C, Q$ tường minh Tập nghiệm ẩn: $\text{Fix}(T), \text{Sol}(C, F_1), \text{Sol}(C, f)$
Tính chất toán tử Toán tử tuyến tính bị chặn $A$ Ánh xạ không giãn, đơn điệu Toán tử giả đơn điệu, liên tục Lipschitz
Bản chất hội tụ Hội tụ yếu (Weak convergence) Hội tụ mạnh (Strong convergence) Hội tụ mạnh (Strong convergence)
Yêu cầu phép chiếu 02 phép chiếu trực giao $P_C, P_Q$ Chiếu trực giao $P_C, P_Q$ Chiếu nửa không gian biến thiên (Subgradient)
Khả năng giải BSVIP / SEP Không khả thi Không áp dụng được cho BSVIP/SEP Giải quyết triệt để BSVIP và tìm Min-norm SEP

Đó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 lý thuyết toán tử đơn điệu và giải tích phi tuyến qua các đóng góp lý thuyết nền tảng:

                            ┌──────────────────────────────────────────────────────────┐
                            │               KHUNG PHÂN TÍCH TÍCH HỢP                   │
                            └────────────────────────────┬─────────────────────────────┘
                                                         │
         ┌───────────────────────────────────────────────┼──────────────────────────────────────────────┐
         ▼                                               ▼                                              ▼
┌─────────────────────────────────┐             ┌─────────────────────────────────┐            ┌─────────────────────────────────┐
│     GIẢI TÍCH LỒI NÂNG CAO      │             │    TOÁN TỬ KHÔNG GIÃN & LẶP     │            │    CHIẾU DƯỚI ĐẠO HÀM TĂNG CƯỜNG│
│ • Định lý Moreau-Rockafellar    │             │ • Kỹ thuật Krasnoselskii-Mann   │            │ • Nửa không gian biến thiên     │
│ • Dưới vi phân & Nón pháp tuyến │             │ • Nguyên lý bán đóng (Demiclose)│            │   $T_k = \{\omega: \dots \le 0\}$│
│ • Song hàm cân bằng $f(x, y)$   │             │ • Phép chiếu metric $P_C, P_Q$  │            │ • Cắt siêu phẳng siêu tốc       │
└────────────────┬────────────────┘             └────────────────┬────────────────┘            └────────────────┬────────────────┘
                 │                                               │                                              │
                 └───────────────────────────────────────────────┼──────────────────────────────────────────────┘
                                                                 │
                                                                 ▼
                            ┌──────────────────────────────────────────────────────────┐
                            │           THUẬT TOÁN HỘI TỤ MẠNH ĐA KHÔNG GIAN           │
                            │              $x_{k+1} = \alpha_k x_k + (1-\alpha_k)T(z_k)$ │
                            └──────────────────────────────────────────────────────────┘
  1. Xây dựng toán tử tựa không giãn $T_f$ từ song hàm giả đơn điệu: Đối với song hàm cân bằng $f: H \times H \to \mathbb{R} \cup {+\infty}$ thỏa mãn các điều kiện Blum-Oettli mở rộng và điều kiện kiểu Lipschitz với các hằng số $c_1, c_2 > 0$, luận án đã xây dựng ánh xạ đơn trị $T_f: C \to C$: $$T_f(x) := \text{argmin}\left{\lambda f(s(x), y) + \frac{1}{2}|y - x|^2 : y \in C\right}$$ trong đó: $$s(x) := \text{argmin}\left{\lambda f(x, y) + \frac{1}{2}|y - x|^2 : y \in C\right}$$ với $0 < \lambda < \min\left{\frac{1}{2c_1}, \frac{1}{2c_2}\right}$. Luận án chứng minh rằng $T_f$ là ánh xạ tựa không giãn (quasi-nonexpansive), thỏa mãn $\text{Fix}(T_f) = \text{Sol}(C, f)$ và đặc biệt chứng minh $T_f$ thỏa mãn nguyên lý bán đóng (demiclosedness principle)—kết quả chưa từng có trong các công trình của Patrick Louis Combettes và Sever Adrian Hirstoaga (2005).

  2. Thiết lập mô hình bất đẳng thức biến phân hai cấp tách (BSVIP): Định nghĩa bài toán tìm $x^* \in \Omega$ sao cho $\langle F(x^), x - x^ \rangle \ge 0, \forall x \in \Omega$, với tập ràng buộc $\Omega = {x^* \in \text{Sol}(C, F_1) : Ax^* \in \text{Sol}(Q, F_2)}$. Đây là bước chuyển đổi hệ hình toán học từ bài toán tối ưu trên tập lồi tường minh sang tối ưu trên đa tạp nghiệm của hệ phương trình toán tử phi tuyến phân tách liên không gian.

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

Khung phân tích tích hợp chặt chẽ 3 trụ cột toán học:

  • Trụ cột 1 - Lý thuyết dưới vi phân Moreau-Rockafellar: Khai thác tính chất $\partial(f_1 + f_2)(x) = \partial f_1(x) + \partial f_2(x)$ và đặc trưng nón pháp tuyến ngoài $N_C(x) = \partial \delta_C(x)$ để giải tích hóa các điều kiện tối ưu bậc một của các bài toán quy nạp hàm chỉ.
  • Trụ cột 2 - Lý thuyết toán tử không giãn và kỹ thuật lặp Krasnoselskii-Mann: Sử dụng bổ đề lồi dãy Banach và bổ đề Opial: $$\liminf_{i \to \infty} |y^{k_i} - y| < \liminf_{i \to \infty} |y^{k_i} - T(y)| \quad (\text{khi } y \ne T(y))$$ để thiết lập tính đóng yếu của tập nghiệm và triệt tiêu sai số phân kỳ.
  • Trụ cột 3 - Kỹ thuật xấp xỉ liên tục liên không gian: Kết nối không gian Hilbert $H_1$ và $H_2$ thông qua toán tử tuyến tính bị chặn $A$ và toán tử liên hợp $A^$, đảm bảo $|A^| = |A|$.

Điều kiện biên xác định tường minh: bước điều khiển $\delta \in \left(0, \frac{1}{|A|^2 + 1}\right)$, hệ số đơn điệu mạnh $\mu \in \left(0, \frac{2\beta}{L^2}\right)$, dãy số ${\alpha_k} \subset (0, 1)$ thỏa mãn $\lim_{k \to \infty} \alpha_k = \alpha \in (0, 1)$, và dãy suy giảm ${\lambda_k} \subset (0, 1)$ thỏa mãn $\lim_{k \to \infty} \lambda_k = 0, \sum_{k=0}^{\infty} \lambda_k(1-\alpha_k) = \infty$.

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ủ triết lý duy thực phản biện (Critical Realism) và chủ nghĩa thực chứng diễn dịch (Deductive Positivism) trong toán học thuần túy. Thiết kế nghiên cứu đa tầng kết hợp chứng minh giải tích giải tích trừu tượng trong không gian Hilbert vô hạn chiều với kiểm chứng thực nghiệm mô phỏng số.

                               ┌───────────────────────────────────────────────────────────┐
                               │                 SƠ ĐỒ THUẬT TOÁN BỐN BƯỚC                 │
                               └─────────────────────────────┬─────────────────────────────┘
                                                             │
                 ┌───────────────────────────────────────────┴───────────────────────────────────────────┐
                 │                                                                                       │
                 ▼                                                                                       ▼
┌─────────────────────────────────┐                                                     ┌─────────────────────────────────┐
│     BƯỚC 1: CHIẾU TRÊN $H_2$    │                                                     │    BƯỚC 2: KÉO NGƯỢC VỀ $H_1$   │
│   $u_k = P_Q(Ax_k)$             │ ──────────────────────────────────────────────────> │ $y_k = P_C(x_k + \delta A^*(Su_k - Ax_k))$ │
└─────────────────────────────────┘                                                     └────────────────┬────────────────┘
                                                                                                         │
                 ┌───────────────────────────────────────────────────────────────────────────────────────┘
                 │
                 ▼                                                                                       ▼
┌─────────────────────────────────┐                                                     ┌─────────────────────────────────┐
│     BƯỚC 3: CHIẾU GRADIENT VIP  │                                                     │   BƯỚC 4: TỔ HỢP LẶP PHI TUYẾN  │
│ $z_k = P_C(y_k - \lambda_k \mu F(y_k))$ │ ──────────────────────────────────────────> │ $x_{k+1} = \alpha_k x_k + (1-\alpha_k)T(z_k)$ │
└─────────────────────────────────┘                                                     └─────────────────────────────────┘

Quy trình nghiên cứu rigorous

Quy trình lặp tổng quát giải VIP trên tập nghiệm của bài toán điểm bất động tách (SFPP) được thực thi nghiêm ngặt qua 4 bước:

  1. Bước 1 (Xử lý không gian ảnh $H_2$): Tính hình chiếu metric trên tập lồi $Q$: $$u_k = P_Q(Ax_k)$$
  2. Bước 2 (Kéo ngược thông tin liên hợp về không gian nguồn $H_1$): $$y_k = P_C\big(x_k + \delta A^*(Su_k - Ax_k)\big)$$
  3. Bước 3 (Thực hiện bước suy giảm giải VIP cấp trên): $$z_k = P_C\big(y_k - \lambda_k \mu F(y_k)\big)$$
  4. Bước 4 (Tổ hợp lặp Krasnoselskii-Mann với toán tử không giãn $T$): $$x_{k+1} = \alpha_k x_k + (1 - \alpha_k)T(z_k)$$

Quy trình chứng minh hội tụ 7 bước mẫu mực:

  • Bước 1: Thiết lập bất đẳng thức suy giảm cơ sở: $$|y_k - x^|^2 \le |x_k - x^|^2 - \delta(1 - \delta|A|^2)|Su_k - Ax_k|^2 - \delta|u_k - Ax_k|^2$$
  • Bước 2: Chứng minh tính bị chặn (boundedness) của các dãy ${x_k}, {y_k}, {z_k}$ và ${F(y_k)}$.
  • Bước 3: Chứng minh tính tiệm cận của khoảng cách giữa hai bước lặp: $$|y_{k+1} - y_k|^2 \le |x_{k+1} - x_k|^2 - \delta(1 - \delta|A|^2)|Su_{k+1} - Su_k + Ax_k - Ax_{k+1}|^2$$
  • Bước 4: Áp dụng bổ đề dãy suy rộng chứng minh: $$\lim_{k \to \infty} |x_{k+1} - x_k| = 0, \quad \lim_{k \to \infty} |x_k - T(y_k)| = 0$$
  • Bước 5: Thiết lập tính triệt tiêu vi sai toán tử: $$\lim_{k \to \infty} |y_k - T(y_k)| = 0, \quad \lim_{k \to \infty} |u_k - Su_k| = 0, \quad \lim_{k \to \infty} |Ax_k - u_k| = 0$$
  • Bước 6: Chứng minh giới hạn trên biến phân qua Bổ đề Opial: $$\limsup_{k \to \infty} \langle F(x^), x^ - y_k + \lambda_k \mu F(y_k) \rangle \le 0$$
  • Bước 7: Áp dụng bổ đề bất đẳng thức giải tích thực chứng minh dãy ${x_k}$ hội tụ mạnh $x_k \to x^* = P_\Omega(x_0)$.

Data và phân tích

Toàn bộ các thuật toán đề xuất đều được hiện thực hóa và kiểm thử số độc lập:

  • Cấu hình phần cứng và môi trường: Máy tính xách tay CPU Intel(R) Core(TM) i3-3217U @ 1.80GHz, bộ nhớ trong 2 GB RAM, hệ điều hành Windows/Linux.
  • Nền tảng phần mềm: MATLAB R2012a, sử dụng các hàm tối ưu hóa lồi chuyên dụng và các phép toán ma trận đại số tuyến tính.
  • Tiêu chuẩn dừng (Stopping criterion): Thuật toán tự động dừng khi sai số tương đối giữa hai bước lặp liên tiếp đạt ngưỡng dung sai: $$E_k = |x_{k+1} - x_k| < \varepsilon \quad (\text{với } \varepsilon = 10^{-4}, 10^{-6})$$
  • Kiểm định tính vững (Robustness checks): Các thử nghiệm số được tiến hành với nhiều điểm xuất phát ngẫu nhiên $x_0 \in C$, ma trận toán tử ngẫu nhiên $A$, các hằng số bước $\delta, \lambda_k, \alpha_k$ khác nhau nhằm kiểm tra tốc độ hội tụ CPU time (giây) và số bước lặp (Number of iterations - Iter). Kết quả nghiệm tính toán được so sánh trực tiếp với nghiệm giải tích chính xác (Exact solution) để xác nhận sai số hội tụ tuyệt đối.

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

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

  1. Khám phá nguyên lý bán đóng cho lớp song hàm cân bằng giả đơn điệu: Chứng minh hoàn chỉnh rằng toán tử lặp $T_f$ thỏa mãn nguyên lý bán đóng: nếu dãy ${x_n} \subset C$ hội tụ yếu $x_n \rightharpoonup x$ và $|T_f(x_n) - x_n| \to 0$ thì $x \in \text{Fix}(T_f) = \text{Sol}(C, f)$. Đây là tiền đề lý thuyết quyết định mở đường cho việc giải bài toán cân bằng phi tuyến không compact.
  2. Hội tụ mạnh không điều kiện trên tập nghiệm SFPP: Khác với thuật toán CQ chỉ đạt hội tụ yếu trong không gian Hilbert vô hạn chiều, thuật toán đề xuất đạt chuẩn hội tụ mạnh với tốc độ suy giảm sai số tiệm cận bậc nhất $O(1/k)$.
  3. Giải thuật đầu tiên cho bài toán BSVIP: Đề xuất và chứng minh thành công thuật toán dưới đạo hàm tăng cường giải bài toán bất đẳng thức biến phân hai cấp với ràng buộc là tập nghiệm của bài toán biến phân tách, trong đó các ánh xạ giá $F_1, F_2$ chỉ cần thỏa mãn tính giả đơn điệu và liên tục Lipschitz.
  4. Giải pháp đột phá tìm nghiệm chuẩn nhỏ nhất của SEP: Thiết lập phương pháp đạo hàm tăng cường liên kết hai không gian thông qua toán tử tuyến tính bị chặn $A$ và toán tử liên hợp $A^*$, giải quyết triệt để bài toán tìm nghiệm có năng lượng tối thiểu (minimum-norm) của bài toán cân bằng tách.
                                  ┌───────────────────────────────────────────────────────────┐
                                  │                HỆ GIÁ TRỊ TÁC ĐỘNG ĐA CHIỀU               │
                                  └─────────────────────────────┬─────────────────────────────┘
                                                                │
                 ┌──────────────────────────────────────────────┼──────────────────────────────────────────────┐
                 ▼                                              ▼                                              ▼
┌─────────────────────────────────┐            ┌─────────────────────────────────┐            ┌─────────────────────────────────┐
│        Ý NGHĨA LÝ THUYẾT        │            │       Ý NGHĨA PHƯƠNG PHÁP       │            │       ỨNG DỤNG THỰC TIỄN        │
│ • Hoàn thiện giải tích đơn điệu │            │ • Thay thế chiếu tập ẩn bằng    │            │ • Liệu pháp xạ trị IMRT         │
│ • Cầu nối VIP - Điểm bất động   │            │   chiếu nửa không gian          │            │ • Khôi phục ảnh y khoa, MRI     │
│ • Mở rộng không gian vô hạn     │            │ • Lặp song song cho MSSFP       │            │ • Tối ưu hóa mạng viễn thông    │
└─────────────────────────────────┘            └─────────────────────────────────┘            └─────────────────────────────────┘

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp khung giải tích toán học thống nhất bao bọc các lớp bài toán tối ưu hóa, bài toán bù (Complementarity problems), bài toán điểm bất động Brouwer/Kakutani và cân bằng mạng lưới giao thông Nash-Wardrop trong cùng một mô hình tổng quát.
  • Về mặt phương pháp luận: Xóa bỏ hoàn toàn gánh nặng tính toán của phép chiếu trực giao lên tập lồi phức tạp bằng kỹ thuật chiếu trên nửa không gian biến thiên $T_k$, mở rộng khả năng áp dụng thuật toán lên các không gian hàm vô hạn chiều như $L^2(\Omega)$ và không gian Sobolev $H^1(\Omega)$.
  • Về mặt ứng dụng thực tiễn:
    • Liệu pháp xạ trị IMRT: Cho phép xác định liều lượng chiếu xạ tối ưu lên khối u (tập $Q$) đồng thời bảo vệ nghiêm ngặt các mô lành xung quanh (tập $C$), trong đó $A$ mô hình hóa ma trận hấp thụ bức xạ.
    • Xử lý tín hiệu và ảnh: Giải quyết bài toán khôi phục ảnh mờ/nhiễu (Image deblurring/denoising) và nén dữ liệu cảm biến (Compressed sensing) với số lượng ràng buộc đa tập hợp khổng lồ.

Limitations và Future Research

Mặc dù đạt được những đóng góp đột phá, luận án cũng ghi nhận các giới hạn học thuật khách quan:

  • Giới hạn 1 (Sự phụ thuộc vào hằng số Lipschitz): Các tham số bước $\mu, \lambda, \delta$ đòi hỏi phải biết trước hoặc chặn trên được các hằng số Lipschitz $L, L_1, L_2$, các hệ số kiểu Lipschitz $c_1, c_2$ và chuẩn toán tử $|A|$. Trong thực tế quy mô lớn, việc ước lượng chính xác $|A|$ đòi hỏi chi phí tính toán trị riêng ma trận rất tốn kém.
  • Giới hạn 2 (Phạm vi không gian khảo sát): Toàn bộ lý thuyết được phát triển trên không gian Hilbert thực. Việc mở rộng sang không gian Banach tổng quát (nơi tích vô hướng được thay bằng cặp đối ngẫu và toán tử chiếu metric được thay bằng toán tử chiếu mòn/chiếu suy rộng) vẫn chưa được giải quyết.
  • Giới hạn 3 (Kiểm thử thực nghiệm): Các ví dụ số chủ yếu được thực hiện trên không gian hữu hạn chiều $\mathbb{R}^n$ với số chiều vừa phải ($n \le 1000$).

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

  1. Phát triển các thuật toán với kỹ thuật bước tự thích nghi (Self-adaptive stepsize rules / Armijo-type line search) nhằm triệt tiêu hoàn toàn sự phụ thuộc vào hằng số Lipschitz và chuẩn ma trận $|A|$.
  2. Mở rộng lý thuyết bất đẳng thức biến phân hai cấp tách lên không gian Banach phản xạ và không gian Hadamard (không gian metric có độ cong âm).
  3. Nghiên cứu bài toán biến phân ngẫu nhiên tách (Stochastic Split Variational Inequalities) ứng dụng trong học máy và tối ưu hóa dữ liệu lớn (Big Data Optimization).
  4. Tích hợp thuật toán vào các hệ thống xạ trị y tế thời gian thực và xử lý ảnh siêu phân giải.

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

  • Tác động học thuật: Các kết quả của luận án đã được công bố trên 6 công trình khoa học quốc tế thuộc hệ thống ISI/Scopus (bao gồm các tạp chí uy tín hàng đầu như Numerical Algorithms, Applied Mathematics and Computation, Journal of Global Optimization, Optimization). Công trình nhận được sự quan tâm và trích dẫn sâu rộng từ cộng đồng giải tích phi tuyến quốc tế.
  • Tác động ngành và công nghệ: Cung cấp thuật toán nền tảng cho việc phát triển các phần mềm tính toán liều lượng xạ trị y tế IMRT chính xác cao, giúp rút ngắn thời gian lập kế hoạch điều trị ung thư và giảm thiểu tối đa tổn thương mô lành.
  • Tác động xã hội và giáo dục: Đóng góp tài liệu học thuật chuyên khảo cao cấp cho công tác đào tạo sau đại học ngành Toán giải tích và Toán ứng dụng tại Việt Nam, nâng cao vị thế toán học ứng dụng Việt Nam trên trường quốc tế.

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

  • Nghiên cứu sinh và học viên cao học Toán học: Tiếp cận phương pháp luận chuẩn mực về giải tích phi tuyến, kỹ thuật chứng minh hội tụ mạnh và công cụ toán tử lặp đa cấp.
  • Các nhà khoa học thuộc lĩnh vực Tối ưu hóa: Khai thác các cấu trúc thuật toán mới để giải quyết các bài toán tối ưu đa mục tiêu, bài toán bù suy rộng và cân bằng Nash hai cấp.
  • Kỹ sư R&D trong ngành Xử lý tín hiệu và Khôi phục ảnh: Ứng dụng thuật toán lặp song song để xử lý các luồng dữ liệu phân tán quy mô lớn với độ trễ thấp.
  • Chuyên gia Vật lý y khoa và Xạ trị ung thư: Sử dụng mô hình toán học chấp nhận tách để tối ưu hóa quy trình điều biến chùm tia xạ trị lâm sà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 nào?

Đóng góp độc đáo nhất là việc xây dựng thành công toán tử đơn trị $T_f: C \to C$ từ song hàm cân bằng giả đơn điệu $f$ thỏa mãn điều kiện kiểu Lipschitz, chứng minh toán tử này là tựa không giãn và thỏa mãn nguyên lý bán đóng (demiclosedness principle) với $\text{Fix}(T_f) = \text{Sol}(C, f)$. Kết quả này mở rộng trực tiếp lý thuyết của Patrick Louis Combettes và Sever Adrian Hirstoaga (2005) (vốn chỉ áp dụng cho song hàm đơn điệu) và khắc phục thiếu sót cốt tử của các công trình quốc tế trước đó (chẳng hạn như tài liệu [5] trong luận án) khi không chứng minh được nguyên lý bán đóng cho toán tử tựa không giãn.

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

So với thuật toán CQ cổ điển của Charles Byrne (2002) và Hong-Kun Xu (2010) (chỉ áp dụng cho tập lồi tường minh và hội tụ yếu), luận án đã:

  • Tích hợp kỹ thuật lặp Krasnoselskii-Mann và chiếu suy giảm để xử lý tập ràng buộc ẩn $\text{Fix}(T), \text{Fix}(S)$ thu được hội tụ mạnh.
  • Thay thế hoàn toàn phép chiếu trực giao tốn kém trên tập $C$ trong phương pháp đạo hàm tăng cường bằng phép chiếu trên nửa không gian biến thiên $T_k$, cho phép giải quyết bài toán BSVIP với toán tử giả đơn điệu mà không làm bùng nổ chi phí tính toán.

3. Phát hiện bất ngờ nhất trong quá trình nghiên cứu là gì?

Phát hiện bất ngờ nhất là tính ổn định tiệm cận của cấu trúc lặp song song khi giải bài toán VIP trên tập nghiệm chấp nhận tách đa tập hợp (MSSFP). Mặc dù số lượng tập lồi ràng buộc $M$ và $N$ trong hai không gian $H_1, H_2$ tăng lên, thuật toán song song kết hợp toán tử liên hợp $A^*$ vẫn bảo toàn bán kính hội tụ và không làm suy giảm tốc độ hội tụ thực tế của dãy nghiệm.

4. Quy trình tái lập thực nghiệm (Replication protocol) có được cung cấp đầy đủ không?

Luận án cung cấp chi tiết tuyệt đối quy trình thực nghiệm số:

  • Cấu hình: MATLAB R2012a, CPU Intel Core i3-3217U @ 1.80GHz, RAM 2GB.
  • Khởi tạo: Thiết lập ma trận toán tử $A$, các hàm giá $F, F_1, F_2$ đơn điệu mạnh/giả đơn điệu, các tham số bước chính xác: $\delta \in (0, 1/(|A|^2 + 1))$, $\mu \in (0, 2\beta/L^2)$, dãy tham số $\alpha_k = \frac{k}{2(k+1)}$, $\lambda_k = \frac{1}{k+1}$.
  • Tiêu chuẩn kiểm chứng: So sánh nghiệm lặp $x_k$ với nghiệm giải tích chính xác $x^*$, theo dõi trực tiếp sai số qua các thang đo thời gian thực CPU time và số bước lặp (Iter).

5. Chương trình nghị sự nghiên cứu 10 năm tới được định hình như thế nào?

Chương trình nghị sự 10 năm tập trung vào 3 trọng tâm:

  1. Tự thích nghi hóa hoàn toàn (Full self-adaptivity): Loại bỏ hoàn toàn các hằng số tiên nghiệm trong các bước lặp.
  2. Mở rộng cấu trúc không gian: Chuyển dịch toàn bộ hệ thống giải thuật sang không gian Banach hình học và đa tạp Riemann.
  3. Ứng dụng Trí tuệ nhân tạo: Nhúng các thuật toán tối ưu hóa chấp nhận tách vào kiến trúc mạng nơ-ron sâu (Deep Neural Networks) để giải các bài toán học tăng cường đa tác tử và tối ưu hóa phân tán quy mô siêu lớn.

Kết luận

  1. Xây dựng thành công toán tử tựa không giãn $T_f$ thỏa mãn nguyên lý bán đóng từ song hàm cân bằng giả đơn điệu, giải quyết khoảng trống lý thuyết về mối quan hệ giữa bài toán cân bằng và bài toán điểm bất động.
  2. Đề xuất và chứng minh định lý hội tụ mạnh cho thuật toán giải bài toán bất đẳng thức biến phân đơn điệu mạnh trên tập nghiệm của bài toán điểm bất động tách (SFPP) trong không gian Hilbert vô hạn chiều.
  3. Thiết lập thuật toán đầu tiên trong y văn giải bài toán bất đẳng thức biến phân hai cấp với ràng buộc là tập nghiệm của bài toán bất đẳng thức biến phân tách (BSVIP) cho lớp ánh xạ giá giả đơn điệu và liên tục Lipschitz.
  4. Phát triển thuật toán song song hiệu năng cao giải bài toán VIP trên tập nghiệm của bài toán chấp nhận tách đa tập hợp (MSSFP), cải tiến vượt bậc so với các kết quả của GS. Nguyễn Bường (2014).
  5. Đề xuất thuật toán đầu tiên giải quyết bài toán tìm nghiệm có chuẩn nhỏ nhất của bài toán cân bằng tách (SEP) với song hàm giả đơn điệu thỏa mãn điều kiện kiểu Lipschitz.
  6. Toàn bộ 6 công trình khoa học công bố trên các tạp chí quốc tế ISI/Scopus khẳng định tính tiên phong, chuẩn mực học thuật và giá trị ứng dụng thực tiễn to lớn của luận án trong tối ưu hóa hiện đại và xử lý bài toán công nghệ y sinh.