Luận án Tiến sĩ Khuất Thị Bình: Phương pháp lặp cho bài toán chấp nhận tách

Tổng quan về luận án

Bối cảnh khoa học của công trình nghiên cứu bắt nguồn từ lý thuyết điểm bất động của ánh xạ không giãn, lý thuyết toán tử đơn điệu và giải tích lồi hiện đại – những trụ cột toán học nền tảng cho việc giải quyết các bài toán tối ưu hóa, phương trình vi phân phi tuyến và bài toán ngược (inverse problems). Luận án tiến sĩ toán học với đề tài "Một số phương pháp lặp cho bài toán chấp nhận tách và các bài toán liên quan" (Chuyên ngành: Toán ứng dụng, Mã số: 9 46 01 12) do nghiên cứu sinh Khuất Thị Bình thực hiện dưới sự hướng dẫn khoa học của GS.TS. Nguyễn Bường tại Học viện Khoa học và Công nghệ (Viện Hàn lâm Khoa học và Công nghệ Việt Nam), đã thiết lập các bước tiến mang tính tiên phong trong việc phát triển các thuật toán hiệu chỉnh lặp tối ưu hóa.

Khoảng trống nghiên cứu (Research Gap) then chốt mà luận án trực tiếp giải quyết xuất phát từ các giới hạn kỹ thuật nội tại của các phương pháp lặp kinh điển:

1. Trong bài toán chấp nhận tách đa tập (MSSFP) và bài toán chấp nhận tách (SFP) khởi xướng bởi Censor và Elfving (1994), Byrne (2002), các thuật toán chiếu gradient truyền thống và hiệu chỉnh lặp kiểu Bruck–Bakushinsky (Xu, 2010) đòi hỏi kích thước bước lặp $\mu_k$ phải phụ thuộc chặt chẽ vào chuẩn phổ của toán tử chuyển tuyến tính $|A|$, vốn cực kỳ tốn kém hoặc bất khả thi để ước lượng chính xác trong không gian vô hạn chiều.
2. Trong bài toán trùng tách đa tập (MSSEP) do Chen và cộng sự (2013) nghiên cứu, các thuật toán hiện hành bắt buộc phải tính toán trên các chuỗi tổng vô hạn $\sum_{i=1}^\infty$ và $\sum_{j=1}^\infty$ ở từng bước lặp khi xét họ vô hạn đếm được các tập con lồi đóng, dẫn đến tắc nghẽn tính toán (computational intractability) trên thực tế.
3. Trong bài toán bất đẳng thức biến phân (VIP) trên không gian Banach trơn đều, các phương pháp lặp kinh điển như Krasnosel'skii–Mann hay Ishikawa chỉ bảo đảm hội tụ yếu, trong khi phương pháp xấp xỉ mềm Halpern lại gặp khó khăn khi mở rộng sang cấu trúc toán tử tổng quát và các họ ánh xạ vô hạn.

Các câu hỏi và giả thuyết nghiên cứu được thiết lập chuẩn xác:

• RQ1: Làm thế nào để thiết lập thuật toán hiệu chỉnh lặp cho bài toán MSSFP trong không gian Hilbert sao cho bước lặp $\mu_k$ hoàn toàn độc lập với chuẩn toán tử $|A|$ mà vẫn đảm bảo tính hội tụ mạnh?
  • H1: Tồn tại sơ đồ hiệu chỉnh kiểu Lavrentiev $F^k u^k + \epsilon_k(u^k - x^+) = 0$ kết hợp dãy bước lặp tự thích nghi suy giảm cho phép dãy lặp hội tụ mạnh về nghiệm có khoảng cách cực tiểu tới điểm mốc $x^+$.
• RQ2: Có thể khử bỏ phép lấy tổng vô hạn ở mỗi bước lặp trong bài toán MSSEP với họ tập chỉ số vô hạn đếm được hay không?
  • H2: Sơ đồ xấp xỉ hữu hạn lũy tiến (progressive finite-sum approximation) $U^k = \frac{1}{\bar{\alpha}k}\sum{i=1}^k \alpha_i P_{C_i}$ và $V^k = \frac{1}{\bar{\beta}k}\sum{j=1}^k \beta_j P_{Q_j}$ kết hợp hiệu chỉnh Bakushinsky–Bruck sẽ bảo toàn sự hội tụ mạnh về nghiệm có chuẩn cực tiểu.
• RQ3: Làm thế nào để lai ghép phương pháp đường dốc nhất với phương pháp Ishikawa nhằm giải VIP trên tập điểm bất động chung của họ vô hạn ánh xạ không giãn trong không gian Banach?
  • H3: Cấu trúc lai ghép giữa phương pháp gradient dốc nhất và sơ đồ hai bước Ishikawa sẽ nâng cấp bậc hội tụ từ hội tụ yếu lên hội tụ mạnh trong không gian Banach $q$-trơn đều.

Khung lý thuyết của luận án tích hợp: Lý thuyết toán tử đơn điệu cực đại, Lý thuyết hiệu chỉnh Tikhonov–Lavrentiev cho bài toán đặt không chỉnh, Lý thuyết phép chiếu metric trên tập lồi đóng, và Lý thuyết hình học không gian Banach (ánh xạ đối ngẫu chuẩn tắc $j$). Đóng góp đột phá của công trình được định lượng hóa thông qua việc giảm độ phức tạp tính toán từ vô hạn $O(\infty)$ xuống $O(k)$ tại mỗi bước lặp, xóa bỏ $100%$ sự phụ thuộc vào việc tính toán chuẩn ma trận $|A|$, và chứng minh giải tích thành công sự hội tụ mạnh trong không gian vô hạn chiều. Phạm vi nghiên cứu bao quát các không gian hàm vô hạn chiều Hilbert và Banach phản xạ, với ứng dụng thực tiễn được mô hình hóa trong bài toán tái tạo hình ảnh chụp cắt lớp vi tính (CT) và tối ưu hóa liều xạ trị điều biến cường độ (IMRT).

Literature Review và Positioning

Lịch sử phát triển của các phương pháp lặp giải bài toán tách và điểm bất động được định hình bởi ba dòng nghiên cứu chính:

[Dòng 1: Bài toán Điểm bất động & VIP]
   ├── Krasnosel'skii (1955) & Mann (1953) [Hội tụ yếu]
   ├── Ishikawa (1974) [Mở rộng 2 bước, hội tụ yếu]
   └── Halpern (1967) & Moudafi (2000) [Xấp xỉ mềm, hội tụ mạnh]

[Dòng 2: Bài toán Tách SFP & MSSFP]
   ├── Censor & Elfving (1994) [SFP trong R^n, phụ thuộc ||A||]
   ├── Byrne (2002) [CQ Algorithm, bước lặp 0 < gamma < 2/||A||^2]
   ├── Xu H.K. (2006, 2010) [Hiệu chỉnh Bruck-Bakushinsky, phụ thuộc ||A||]
   └── Tian & Zhang (2017) [Tự thích nghi, chưa chứng minh trọn vẹn chuỗi phân kỳ]

[Dòng 3: Bài toán Trùng tách SEP & MSSEP]
   ├── Byrne & Moudafi (2013) [SSEA đồng thời, hội tụ yếu]
   └── Chen et al. (2013, 2014) [Hiệu chỉnh Tikhonov, vướng tổng vô hạn ở mỗi bước k]

Cuộc tranh luận học thuật sâu sắc tập trung vào hai mâu thuẫn lớn:

1. Hội tụ yếu đối kháng Hội tụ mạnh: Trong không gian vô hạn chiều, các thuật toán Mann và Ishikawa nguyên bản chỉ bảo đảm hội tụ yếu. Genel và Lindenstrass (1975) đã đưa ra phản ví dụ kinh điển chứng minh phương pháp Mann không thể hội tụ mạnh nếu không có các điều kiện thu hẹp không gian. Để đạt được hội tụ mạnh, các nhà toán học quốc tế như Wittmann (1992), Xu (2002, 2004) phải áp dụng các điều kiện kiểm soát bước lặp ngặt nghèo trên chuỗi tham số hiệu chỉnh $\sum |\alpha_{k+1} - \alpha_k| < \infty$.
2. Ước lượng chuẩn toán tử đối kháng Thuật toán tự thích nghi: Để đảm bảo tính không giãn của toán tử chuyển tiếp, các thuật toán của Byrne (2002) và Xu (2010) buộc bước lặp $\mu_k \in (0, 2/|A|^2)$. Tian và Zhang (2017) cố gắng khắc phục bằng bước lặp thích nghi $\mu_k = \gamma_k f(x_k)/|A^*(I-P_Q)Ax_k|^2$, nhưng bằng chứng hội tụ bị khiếm khuyết nghiêm trọng do chưa chứng minh được $\sum_{k=1}^\infty \mu_k \epsilon_k = +\infty$ khi $f(x_k) \to 0$.

Positioning của luận án nằm ở giao điểm giải quyết triệt để hai nút thắt trên. So sánh với các công trình quốc tế tiêu biểu:

• So với Xu (2010, Inverse Problems): Thuật toán của Xu yêu cầu điều kiện $0 < \mu_k \le \epsilon_k/(|A|^2 + \epsilon_k)$, nghĩa là phải tính chính xác $|A|$. Luận án đã giải phóng hoàn toàn ràng buộc này bằng kỹ thuật hiệu chỉnh Lavrentiev mới, trong đó tham số $\mu_k$ chỉ cần thỏa mãn các điều kiện đại số phân kỳ độc lập: $\lim_{k\to\infty} \mu_k/\epsilon_k = 0$ và $\sum_{k=1}^\infty \mu_k \epsilon_k = \infty$.
• So với Chen et al. (2013, Fixed Point Theory and Applications): Phương pháp của Chen giải MSSEP yêu cầu tính toán toàn bộ tổng vô hạn $\sum_{j=1}^\infty \beta_j P_{Q_j}$ ngay từ bước lặp đầu tiên ($k=1$). Luận án giải quyết triệt để bằng cấu trúc xấp xỉ tổng hữu hạn mở rộng lũy tiến $k$, loại bỏ độ phức tạp vô hạn mà vẫn duy trì định lý hội tụ mạnh nghiêm ngặt.

Đó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à thách thức các lý thuyết nền tảng của Lavrentiev (1967), Bakushinsky (1967), Bruck (1975), Ishikawa (1974) và Goldstein (1964) thông qua 4 định lý hội tụ đột phá:

"Phương pháp lặp Krasnosel’skii–Mann và Ishikawa cho kết quả hội tụ yếu, trong khi phương pháp lặp Halpern và phương pháp xấp xỉ mềm cho kết quả hội tụ mạnh trong không gian vô hạn chiều." (Trích dẫn từ phần Mở đầu luận án)

Xuất phát từ nhận định cốt lõi đó, tác giả đã phát triển các mô hình lý thuyết mới:

Mô hình 1: Hiệu chỉnh Lavrentiev cho MSSFP vô hạn đếm được Xét phương trình toán tử hiệu chỉnh: $$F^k u^k + \epsilon_k (u^k - x^+) = 0$$ với $F^k := I - U^k - A^*(I - V^k)A$, trong đó $U^k = \frac{1}{\bar{\alpha}k}\sum{i=1}^k \alpha_i P_{C_i}$ và $V^k = \frac{1}{\bar{\beta}k}\sum{j=1}^k \beta_j P_{Q_j}$. Dãy lặp hiện được xây dựng: $$z^{k+1} = (I - \mu_k(F^k + \epsilon_k(I - x^+)))z^k, \quad k \ge 1$$ Luận án chứng minh:

"Dãy lặp ${z^k}$ xác định bởi $z^{k+1} = (I - \mu_k(F^k + \epsilon_k(I - x^+)))z^k, k \ge 1$ hội tụ mạnh tới $p^$ thỏa mãn $|p^ - x^+| \le |p - x^+|, \forall p \in \Gamma$." (Trích dẫn Chương 2)

Mô hình 2: Hiệu chỉnh Bakushinsky–Bruck tổng hữu hạn cho MSSEP Giải quyết bài toán tìm $z = [x, y]$ sao cho $x \in C = \cap_{i=1}^\infty C_i, y \in Q = \cap_{j=1}^\infty Q_j$ thỏa mãn $Ax = By$. Sơ đồ lặp: $$z^{k+1} = P_S^k \left( (1 - \epsilon_k \mu_k) z^k - \mu_k G^* G z^k \right)$$ trong đó phép chiếu $P_S^k$ chỉ thực hiện trên tích trực tiếp các tổng hữu hạn $U^k \times V^k$ với $G = [A \quad -B]$.

Mô hình 3: Lai ghép Đường dốc nhất – Ishikawa cho VIP trong Banach Trên không gian Banach trơn đều $E$, giải bài toán bất đẳng thức biến phân $\langle Fp^, j(p^ - p) \rangle \le 0, \forall p \in \text{Fix}(T)$ với toán tử $F$ là $\gamma$-$j$-đơn điệu mạnh và $\mu$-giả co chặt ($\gamma + \mu > 1$): $$\begin{cases} y^k = (1 - \beta_k) x^k + \beta_k T_k x^k \ x^{k+1} = (1 - \alpha_k) y^k + \alpha_k (I - \lambda_k F) y^k \end{cases}$$ Sơ đồ này tạo ra bước chuyển đổi hệ hình (paradigm shift), khẳng định rằng sự kết hợp lai ghép có thể chuyển hóa hoàn toàn đặc tính hội tụ yếu của cấu trúc Ishikawa thành hội tụ mạnh trong không gian Banach.

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

Khung phân tích của luận án tích hợp liên hoàn ba phân nhánh toán học thuần túy:

1. Toán tử không giãn chặt và Phép chiếu Metric: Tận dụng tính chất $\frac{1}{2}$-ngược đơn điệu mạnh của $I - P_{C_i}$ và tính chất $\frac{1}{2|A|^2}$-ngược đơn điệu mạnh của $A^*(I - P_{Q_j})A$ để thiết lập tính đơn điệu của toán tử tổng hợp $F^k$.
2. Lý thuyết Hiệu chỉnh (Regularization Theory): Áp dụng tham số nhiễu $\epsilon_k \to 0$ để biến bài toán đặt không chỉnh (ill-posed) thành dãy các bài toán đơn điệu mạnh đặt chỉnh (well-posed), có nghiệm duy nhất $u^k$.
3. Giải tích Tiệm cận Dãy số (Asymptotic Sequence Analysis): Vận dụng hệ thống bổ đề điều khiển bước lặp phân kỳ $\sum \mu_k \epsilon_k = \infty$ và kỹ thuật phân tích dãy con không đơn điệu (Maingé Lemma) để khẳng định tính hội tụ mạnh.

Điều kiện biên (Boundary Conditions):

• Không gian $H_1, H_2, H_3$ là các không gian Hilbert thực vô hạn chiều; $E$ là không gian Banach phản xạ, $q$-trơn đều có chuẩn khả vi Gâteaux liên tục.
• Tập nghiệm $\Gamma \neq \emptyset$.
• Dãy trọng số thỏa mãn chuẩn hóa: $\alpha_i > 0, \sum_{i=1}^\infty \alpha_i = 1$; $\beta_j > 0, \sum_{j=1}^\infty \beta_j = 1$.

graph TD
    A["Họ vô hạn tập lồi đóng {C_i}, {Q_j}"] --> B["Xấp xỉ Tổng hữu hạn U^k, V^k"]
    B --> C["Toán tử Đơn điệu F^k = I - U^k - A*(I - V^k)A"]
    C --> D["Hiệu chỉnh Lavrentiev / Bakushinsky-Bruck: F^k + epsilon_k(I - x^+)"]
    D --> E["Sơ đồ lặp tự thích nghi không phụ thuộc ||A||"]
    E --> F["Dãy lặp {z^k} Hội tụ mạnh về Nghiệm tối ưu p* in Gamma"]

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ủ nghiêm ngặt quan điểm nhận thức luận duy lý - diễn dịch hình thức (Rationalist-Deductive Epistemology) của toán học giải tích hàm hiện đại, kết hợp phương pháp luận kiểm chứng số học (Numerical Constructive Verification).

Thiết kế nghiên cứu đa tầng bao gồm:

• Tầng 1 (Trừu tượng hóa toán học): Thiết lập không gian hàm vô hạn chiều (Hilbert $H$ và Banach $q$-trơn đều $E$ như không gian $L^p, \ell^p, W^{m,p}$ với $1 < p < \infty$).
• Tầng 2 (Rời rạc hóa cấu trúc): Chuyển đổi các bài toán tối ưu hóa liên tục chứa toán tử truyền qua phép rời rạc hóa hệ toán tử sang bài toán chấp nhận tách nhiều tập.
• Tầng 3 (Thuật toán hóa): Xây dựng thuật toán lặp phân kỳ có kiểm soát sai số, giải phóng điều kiện tiên nghiệm về chuẩn.
• Tầng 4 (Mô phỏng thực nghiệm số): Thiết kế các kịch bản thử nghiệm số trên máy tính nhằm đánh giá tốc độ hội tụ và sai số chuẩn.

Quy trình nghiên cứu rigorous

Quy trình giải tích logic trải qua 5 bước nghiêm ngặt:

1. Chứng minh tính đặt chỉnh của bài toán phụ hiệu chỉnh: Chứng minh toán tử $F^k + \epsilon_k (I - x^+)$ là $\epsilon_k$-đơn điệu mạnh, từ đó suy ra sự tồn tại và duy nhất của nghiệm $u^k$ với mỗi $k \ge 1$.
2. Đánh giá tính bị chặn (Boundedness): Chứng minh dãy nghiệm xấp xỉ ${u^k}$ và dãy lặp ${z^k}$ luôn nằm trong một hình cầu đóng bị chặn $S(0, M_1)$ trong không gian vô hạn chiều với: $$\sup_{i,j,k \ge 1} \left{ |u^k|, |P_{C_i} u^{k-1}|, |A^* P_{Q_j} A u^{k-1}| \right} \le M_1$$
3. Chứng minh triệt tiêu tiệm cận: Chứng minh $\lim_{k\to\infty} |(I - P_{C_i})u^k| = 0$ và $\lim_{k\to\infty} |(I - P_{Q_j})Au^k| = 0$ với mọi chỉ số $i, j \ge 1$.
4. Đánh giá bước nhảy liên tiếp: Thiết lập bất đẳng thức sai phân: $$|u^k - u^{k-1}| \le \frac{2M_1}{\epsilon_k} \left( \frac{\alpha_k}{\bar{\alpha}_k} + \frac{\beta_k}{\bar{\beta}k} \right) + \frac{|\epsilon{k-1} - \epsilon_k|}{\epsilon_k} (M_1 + |x^+|) =: d_k$$
5. Chứng minh Hội tụ mạnh (Strong Convergence): Vận dụng tính chất không gian phản xạ, nguyên lý compact yếu và toán tử đơn điệu cực đại để chứng minh toàn bộ dãy lặp hội tụ theo chuẩn mạnh về $p^* = P_\Gamma(x^+)$.

Data và phân tích

Mặc dù là công trình toán lý thuyết thuần túy, luận án thực hiện các thực nghiệm số học chính xác trên môi trường tính toán khoa học (MATLAB R2020b / Python SciPy Optimization framework):

• Tham số thực nghiệm chuẩn hóa:
  • Dãy bước lặp: $\mu_k = \frac{1}{(k+1)^a}$, tham số hiệu chỉnh: $\epsilon_k = \frac{1}{(k+1)^b}$ với các số mũ được kiểm soát $0 < b < a$ và $a + 2b < 1$ (ví dụ chọn $a = 0.5, b = 0.2$).
  • Dãy trọng số tập hợp: $\alpha_i = \beta_i = \frac{1}{i(i+1)}$, đảm bảo $\sum_{i=1}^\infty \alpha_i = 1$ và $\bar{\alpha}k = \sum{i=1}^k \alpha_i = 1 - \frac{1}{k+1} \to 1$.
• Độ phức tạp tính toán: Ở bước lặp thứ $k$, thuật toán chỉ yêu cầu tính $k$ phép chiếu metric trên $C_i$ và $k$ phép chiếu trên $Q_j$, thay vì tính tổng vô hạn.
• Tiêu chuẩn dừng thuật toán: Sai số tương đối chuẩn Euclid $|z^{k+1} - z^k|_2 < 10^{-6}$ hoặc $|F^k z^k|_2 < 10^{-5}$.

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

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

Luận án xác lập 4 phát hiện toán học mang tính đột phá:

Kết quả phản trực giác (Counter-intuitive result) nổi bật là: Việc cắt cụt chuỗi vô hạn thành các tổng hữu hạn tăng dần theo bước lặp $k$ không hề làm suy giảm vận tốc hội tụ hay gây lệch nghiệm tiệm cận, mà trái lại còn tạo ra tính ổn định số học vượt trội, ngăn ngừa hiện tượng tích tụ sai số làm tròn (round-off error accumulation).

Implications đa chiều

• Về mặt Lý thuyết: Đóng góp trực tiếp vào giải tích phi tuyến và lý thuyết xấp xỉ nghiệm, mở rộng phổ áp dụng của phương pháp hiệu chỉnh Lavrentiev từ phương trình toán tử tuyến tính sang hệ bất đẳng thức bao hàm thức và bài toán chấp nhận tách.
• Về mặt Phương pháp luận: Cung cấp một framework mẫu mực cho việc chuyển hóa các thuật toán hội tụ yếu thành hội tụ mạnh, và từ thuật toán "lý thuyết vô hạn" thành thuật toán "tính toán hữu hạn khả thi".
• Về mặt Ứng dụng Y tế & Xử lý tín hiệu:
  1. Khôi phục cấu trúc ảnh chụp cắt lớp vi tính (CT): Mô hình hóa bài toán tích phân Radon rời rạc $Ax = b$ với ràng buộc lồi $|x|_1 \le t$ (nén tín hiệu/Lasso). Thuật toán cho phép tái tạo ma trận hệ số hấp thụ tuyến tính $\mu(x,y)$ từ tập hình chiếu $p(u, \theta)$ với độ chính xác cao và triệt tiêu nhiễu đo đạc $w$ ($|w| \le \varepsilon$).
  2. Tối ưu hóa liều xạ trị IMRT: Giải bài toán tìm vector cường độ chùm tia nhỏ $p^* \in X_+ \cap \bigcap_l X_l$ sao cho liều hấp thụ $h^* = Hp^* \in \bigcap_v H_v$. Đảm bảo liều lượng xạ trị tại khối u (PTV) đạt ngưỡng tiêu diệt tế bào ác tính ($E_t(h) \ge E_t^{\min}$) trong khi liều lượng tại các cơ quan lành có nguy cơ (OAR) được kiểm soát nghiêm ngặt dưới ngưỡng an toàn ($E_t(h) \le E_t^{\max}$).

Limitations và Future Research

Luận án thẳng thắn chỉ ra các giới hạn nghiên cứu:

1. Ràng buộc Tính Lồi Chặt chẽ: Toàn bộ khung lý thuyết xây dựng trên giả định các tập $C_i, Q_j$ là tập con lồi đóng. Đối với các bài toán tối ưu hóa không lồi (non-convex split feasibility problems) nảy sinh trong học sâu (Deep Learning), thuật toán chưa thể áp dụng trực tiếp.
2. Chi phí tính toán phép chiếu $P_C, P_Q$: Thuật toán giả định rằng phép chiếu metric lên từng tập thành phần có thể tính được dạng hiện (explicit form). Khi các tập $C_i$ có cấu trúc hình học phi tuyến phức tạp, việc tính $P_{C_i}$ đòi hỏi các thuật toán con (sub-algorithms) giải bài toán tối ưu phụ.
3. Tốc độ hội tụ tiệm cận (Convergence Rate): Do sử dụng bước lặp suy giảm để đảm bảo hội tụ mạnh trong không gian vô hạn chiều, tốc độ hội tụ thực tế đạt bậc $O(1/\sqrt{k})$ hoặc $O(1/k^b)$, chậm hơn so với các phương pháp gia tốc bậc hai.

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

• Hướng 1: Tích hợp kỹ thuật gia tốc quán tính (Inertial extrapolation techniques kiểu Nesterov / FISTA) vào sơ đồ hiệu chỉnh Lavrentiev và Ishikawa để nâng tốc độ hội tụ lên $O(1/k^2)$.
• Hướng 2: Mở rộng thuật toán cho bài toán chấp nhận tách với cấu trúc tập không lồi (Non-convex split equality problems) thông qua phép chiếu gần kề (Proximal operators) và kỹ thuật xấp xỉ đa tạp.
• Hướng 3: Nghiên cứu các phiên bản thuật toán ngẫu nhiên (Stochastic / Block-coordinate iterative algorithms), trong đó ở mỗi bước chỉ chọn ngẫu nhiên một tập con các chỉ số $i \in J_1, j \in J_2$, phục vụ xử lý dữ liệu lớn (Big Data optimization).
• Hướng 4: Ứng dụng thuật toán mở mạng (Deep algorithm unrolling) để thiết kế các kiến trúc mạng nơ-ron sâu giải bài toán ngược y sinh với tốc độ thời gian thực.

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

• Ảnh hưởng Học thuật: Đề tài mở ra hướng tiếp cận mới trong chuyên ngành Giải tích số và Tối ưu hóa, là tài liệu tham khảo nền tảng với tiềm năng trích dẫn cao trong các tạp chí ISI/Scopus chuyên ngành uy tín (Journal of Optimization Theory and Applications, Optimization, Numerical Algorithms, Fixed Point Theory and Algorithms).
• Chuyển đổi Công nghệ R&D Y tế: Cung cấp thuật toán lõi cho các phần mềm lập kế hoạch xạ trị (Treatment Planning Systems - TPS) trong máy gia tốc y tế tuyến tính (LINAC), giúp tối ưu hóa chùm tia xạ trị IMRT chính xác hơn, bảo vệ mô lành của bệnh nhân ung thư.
• Tác động Công nghiệp Số: Tối ưu hóa thuật toán tái tạo ảnh CT công nghiệp và CT y tế, giảm thiểu thời gian quét của đầu thu/nguồn phát X-quang, từ đó giảm liều lượng bức xạ hấp thụ của bệnh nhân mà vẫn bảo đảm độ phân giải hình ảnh chẩn đoán.
• Ý nghĩa Quốc tế: Đóng góp giải pháp học thuật của Việt Nam vào mạng lưới nghiên cứu quốc tế về giải tích phi tuyến và bài toán ngược.

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

• Nghiên cứu sinh & Nhà toán học trẻ: Tiếp cận một hệ thống kỹ thuật chứng minh giải tích chuẩn mực về hội tụ mạnh trong không gian vô hạn chiều, làm chủ phương pháp hiệu chỉnh toán tử.
• Các nhà nghiên cứu Tối ưu hóa & Giải tích hàm: Khai thác khung thuật toán tự thích nghi không phụ thuộc chuẩn $|A|$ và kỹ thuật tổng hữu hạn lũy tiến để phát triển cho các lớp bài toán biến phân phức tạp hơn.
• Kỹ sư R&D Xử lý tín hiệu & Hình ảnh Y tế: Ứng dụng trực tiếp thuật toán lặp để giải bài toán khôi phục tín hiệu thưa (Sparse signal recovery), khử nhiễu ảnh và tái tạo ảnh chụp cắt lớp.
• Chuyên gia Vật lý Y khoa & Bác sĩ Xạ trị: Hưởng lợi từ mô hình phân bổ liều IMRT chính xác, giúp cân bằng tối ưu giữa việc tiêu diệt khối u và bảo tồn tế bào khỏe mạnh.

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 phát triển thành công Phương pháp hiệu chỉnh Lavrentiev tự thích nghi cho bài toán MSSFP trong không gian Hilbert vô hạn chiều. Luận án đã mở rộng lý thuyết hiệu chỉnh Lavrentiev kinh điển (vốn chỉ áp dụng cho phương trình toán tử tuyến tính đơn điệu không âm $Ax + \epsilon x = f$) sang hệ toán tử phi tuyến phức hợp $F^k = I - U^k - A^*(I - V^k)A$ trên các họ vô hạn tập lồi đóng, đồng thời giải phóng hoàn toàn điều kiện bước lặp khỏi chuẩn phổ $|A|$.

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 ở điểm nào?

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

• So với Xu H.K. (2010): Loại bỏ triệt để điều kiện $0 < \mu_k \le \epsilon_k/(|A|^2 + \epsilon_k)$, giúp thuật toán vận hành mà không cần biết trước thông tin về toán tử $A$.
• So với Chen et al. (2013): Luận án đã thay thế việc lấy tổng vô hạn ở từng bước lặp bằng các toán tử trung bình xấp xỉ hữu hạn mở rộng lũy tiến $U^k = \frac{1}{\bar{\alpha}k}\sum{i=1}^k \alpha_i P_{C_i}$, giảm độ phức tạp từ $\infty$ xuống hữu hạn mà vẫn chứng minh được tính hội tụ mạnh.

3. Kết quả bất ngờ nhất có sự hỗ trợ của dữ liệu/chứng minh toán học là gì?

Kết quả bất ngờ nhất là: Dãy lặp ${z^k}$ sinh bởi sơ đồ hiệu chỉnh lặp Lavrentiev cấp số hữu hạn $U^k, V^k$ vẫn bảo toàn tính hội tụ mạnh trực tiếp về nghiệm $p^* = P_\Gamma(x^+)$ mà không đòi hỏi bất kỳ điều kiện compact địa phương nào trên không gian Hilbert vô hạn chiều $H_1, H_2$. Dữ liệu đánh giá bước nhảy qua bất đẳng thức sai phân xác nhận rằng số hạng dư $d_k \to 0$ khi $k \to \infty$ hoàn toàn kiểm soát được sai số tích lũy.

4. Luận án có cung cấp quy trình tái lặp (Replication Protocol) hoàn chỉnh không?

Hoàn toàn đầy đủ. Toàn bộ các công thức toán tử $U^k, V^k, F^k$, điều kiện tham số $\mu_k, \epsilon_k, \alpha_i, \beta_j$, cùng các thuật toán lặp tường minh (2.11), (1.19), (1.31) đều được mô tả chi tiết từng bước. Các nhà nghiên cứu khác có thể lập trình tái lập chính xác $100%$ các thực nghiệm số trên bất kỳ ngôn ngữ tính toán khoa học nào (MATLAB, Python, C++).

5. Khung chương trình nghiên cứu 10 năm được định hình như thế nào?

Luận án định hình lộ trình nghiên cứu mở rộng sang: (1) Thuật toán gia tốc quán tính bậc cao cho bài toán chấp nhận tách; (2) Giải bài toán tách với ràng buộc tập không lồi (Non-convex split problems); (3) Ứng dụng trong học máy tối ưu phân tán (Distributed optimization); và (4) Tích hợp vào hệ thống chẩn đoán hình ảnh y sinh thời gian thực.

Kết luận

Luận án tiến sĩ của NCS. Khuất Thị Bình đã hoàn thành xuất sắc các mục tiêu nghiên cứu đặt ra, tạo nên những dấu ấn học thuật nổi bật:

1. Xây dựng thành công phương pháp hiệu chỉnh lặp Lavrentiev mới cho bài toán chấp nhận tách đa tập (MSSFP), giải phóng hoàn toàn sự phụ thuộc của kích thước bước lặp vào chuẩn toán tử $|A|$.
2. Giải quyết triệt để bài toán trùng tách đa tập (MSSEP) trên các họ vô hạn đếm được các tập lồi đóng bằng phương pháp hiệu chỉnh lặp Bakushinsky–Bruck với tổng hữu hạn lũy tiến ở mỗi bước lặp, xóa bỏ rào cản tính toán tổng vô hạn.
3. Thiết lập sơ đồ lai ghép tiên tiến giữa phương pháp đường dốc nhất và lặp Ishikawa, đạt được sự hội tụ mạnh cho bài toán bất đẳng thức biến phân (VIP) trên tập điểm bất động chung của họ vô hạn ánh xạ không giãn trong không gian Banach trơn đều.
4. Chứng minh giải tích nghiêm ngặt toàn bộ các định lý hội tụ mạnh trong không gian vô hạn chiều, bảo đảm tính đúng đắn toán học tuyệt đối.
5. Mô hình hóa thành công ứng dụng thực tiễn trong bài toán tái tạo cấu trúc ảnh chụp cắt lớp vi tính (CT) và tối ưu hóa liều lượng chùm tia xạ trị điều biến cường độ (IMRT).
6. Mở ra 3 hướng nghiên cứu mới: Tối ưu hóa cấu trúc không lồi, thuật toán quán tính gia tốc, và tối ưu hóa ngẫu nhiên phân tán trong thời đại khoa học dữ liệu.