Tổng quan luận án

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" 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 (chuyên ngành Toán ứng dụng, mã số: 9 46 01 12, bảo vệ năm 2024).

Tính cấp thiết và khoảng trống nghiên cứu

Lý thuyết điểm bất động của ánh xạ không giãn và các mở rộng của nó giữ vị trí trung tâm trong giải tích phi tuyến, tối ưu hóa và giải phương trình vi phân. Lý thuyết này liên hệ mật thiết với các bài toán thực tế như bài toán chấp nhận lồi, bài toán chấp nhận tách đa tập (MSSFP) và bài toán trùng tách đa tập (MSSEP) – vốn phát sinh trực tiếp từ xử lý tín hiệu số, khôi phục hình ảnh cắt lớp và điều biến liều xạ trị ung thư.

Tuy nhiên, các phương pháp giải các bài toán trên trong không gian vô hạn chiều còn tồn tại các hạn chế kỹ thuật rõ rệt:

  • Nhiều phương pháp lặp cổ điển (như Krasnosel’skii-Mann, Ishikawa) chỉ đảm bảo tính hội tụ yếu trong không gian vô hạn chiều.
  • Trong bài toán MSSFP, việc xác định bước lặp $\mu_k$ của các thuật toán chiếu gradient truyền thống thường phụ thuộc vào việc tính chuẩn ma trận/toán tử chuyển $|A|$. Mặc dù đã có các hướng tiếp cận tự thích nghi hoặc hiệu chỉnh nhằm loại bỏ tham số này, một số công trình trước đó (như nghiên cứu của Tian và Zhang) vẫn chưa hoàn thiện phần chứng minh toán học về sự phân kỳ của chuỗi tham số bước lặp $\sum_{k=1}^\infty \mu_k \epsilon_k = +\infty$.
  • Đối với bài toán MSSFP và MSSEP với các họ vô hạn đếm được các tập con lồi đóng, các phương pháp của Xu hay Chen và cộng sự đòi hỏi việc tính toán trên các chuỗi/tổng vô hạn ở từng bước lặp, gây trở ngại lớn khi đưa vào tính toán số thực tế.
  • Trong bài toán bất đẳng thức biến phân (VIP) trên không gian Banach, các thuật toán giải trên tập điểm bất động chung của họ vô hạn ánh xạ không giãn còn hạn chế về khả năng đảm bảo tính hội tụ mạnh.

Mục tiêu nghiên cứu

Luận án xác định 03 mục tiêu nghiên cứu cụ thể:

  1. Mục tiêu 1: Đề xuất phương pháp hiệu chỉnh lặp mới xấp xỉ nghiệm bài toán chấp nhận tách đa tập (MSSFP) trong không gian Hilbert thực, trong đó tham số bước lặp $\mu_k$ được lựa chọn hoàn toàn độc lập với chuẩn của toán tử chuyển $|A|$.
  2. Mục tiêu 2: Đề xuất phương pháp hiệu chỉnh lặp mới xấp xỉ nghiệm bài toán trùng tách đa tập (MSSEP) trong không gian Hilbert thực, giải quyết triệt để vấn đề tính toán bằng cách chỉ sử dụng các tổng hữu hạn tại mỗi bước lặp thay cho các tổng vô hạn.
  3. Mục tiêu 3: Đề xuất phương pháp lặp lai ghép giữa phương pháp đường dốc nhất (steepest descent) và phương pháp lặp Ishikawa nhằm tìm xấp xỉ nghiệm 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 một họ vô hạn ánh xạ không giãn trong không gian Banach, thiết lập định lý hội tụ mạnh cho dãy lặp.

Đối tượng và phạm vi nghiên cứu

  • Đối tượng nghiên cứu: Các bài toán tối ưu và giải tích phi tuyến gồm: Bài toán chấp nhận tách đa tập (MSSFP), bài toán trùng tách đa tập (MSSEP/ASEP), bài toán bất đẳng thức biến phân (VIP) và bài toán tìm điểm bất động của ánh xạ không giãn/họ ánh xạ không giãn.
  • Phạm vi nghiên cứu:
    • Về mặt không gian: Không gian Hilbert thực và không gian Banach thực (bao gồm các không gian phản xạ, lồi đều, trơn đều, $q$-trơn đều, có chuẩn khả vi Gâteaux hoặc Fréchet) ở cả trường hợp hữu hạn chiều và vô hạn chiều.
    • Về mặt toán tử: Các ánh xạ tuyến tính bị chặn, ánh xạ không giãn, ánh xạ co, ánh xạ trung bình, ánh xạ không giãn chặt, toán tử đơn điệu, đơn điệu cực đại, $\lambda$-đơn điệu mạnh và $\mu$-ngược đơn điệu mạnh.

Tổng quan tài liệu và vị trí của luận án

Luận án đã tổng thuật hệ thống các công trình nghiên cứu trong nước và quốc tế liên quan trực tiếp đến từng nhóm bài toán:

1. Bài toán điểm bất động và các sơ đồ lặp cơ sở

  • Krasnosel’skii (1955) và Mann (1953): Thiết lập phương pháp lặp Mann $x_{k+1} = \alpha_k x_k + (1-\alpha_k)Tx_k$. Genel và Lindenstrass (1975) đã chỉ ra rằng trong không gian vô hạn chiều, phương pháp Mann chỉ cho kết quả hội tụ yếu. Reich (1979) mở rộng sự hội tụ yếu sang không gian Banach lồi đều.
  • Ishikawa (1974): Đề xuất sơ đồ lặp hai bước Ishikawa nhằm mở rộng phương pháp Mann cho lớp ánh xạ giả co Lipschitz, tuy nhiên thuật toán vẫn chỉ đạt được sự hội tụ yếu.
  • Browder (1967) và Halpern (1967): Halpern đề xuất phương pháp lặp $x_{k+1} = \alpha_k u + (1-\alpha_k)Tx_k$ chứng minh sự hội tụ mạnh về hình chiếu của $u$ lên tập điểm bất động. Lions (1977), Wittmann (1992), Xu (2002) lần lượt cải tiến và làm suy giảm các điều kiện trên dãy tham số ${\alpha_k}$.
  • Moudafi (2000) và Xu (2004): Giới thiệu và phát triển phương pháp xấp xỉ mềm (viscosity approximation method) bằng cách thay điểm neo $u$ bằng một ánh xạ co $f(x)$, thiết lập tính hội tụ mạnh trong không gian Hilbert.
  • Các phương pháp lai ghép: Kim và Xu (2005), Yao và cộng sự (2007), Qin và cộng sự (2009), Li (2010), Shehu (2012) đề xuất các sơ đồ kết hợp giữa Halpern, Ishikawa, Mann với phương pháp xấp xỉ mềm hoặc phương pháp đường dốc nhất (Yamada, Ceng và cộng sự 2008).

2. Bài toán chấp nhận tách (SFP) và chấp nhận tách đa tập (MSSFP)

  • Censor và Elfving (1994): Lần đầu tiên phát biểu bài toán SFP và đề xuất thuật toán giải dựa trên phương pháp chiếu gradient trong không gian Euclid hữu hạn chiều cho ứng dụng xử lý ảnh và xạ trị. Hạn chế cốt lõi là bước lặp phụ thuộc vào hệ số Lipschitz của gradient, đòi hỏi tính toán chuẩn ma trận $|A|$.
  • Zhao và Yang (2005, 2013): Đề xuất kỹ thuật tìm kiếm theo tia kiểu Armijo và phương pháp tự thích nghi để tránh việc ước lượng chuẩn toán tử.
  • Xu (2006, 2010): Đưa ra các sơ đồ lặp kế tiếp, đồng thời và tuần hoàn; đề xuất phương pháp hiệu chỉnh lặp kiểu Bruck-Bakushinsky nhằm thu được sự hội tụ mạnh. Hạn chế là bước lặp $\mu_k \in (0, \epsilon_k/(|A|^2+\epsilon_k))$ vẫn phụ thuộc $|A|$, và khi mở rộng cho họ vô hạn tập hợp thì phải tính toán trên chuỗi vô hạn.
  • Tian và Zhang (2013): Đề xuất bước lặp tự thích nghi loại bỏ $|A|$, nhưng phần chứng minh toán học chưa chứng minh được chuỗi tích $\sum \mu_k \epsilon_k$ phân kỳ khi hàm mục tiêu tiến về 0.
  • Nguyễn Bường và cộng sự (2018): Mở rộng giải thuật cho họ vô hạn tập hợp bằng cách rút gọn về tổng hữu hạn tại mỗi bước, song vẫn chưa loại bỏ được chuẩn $|A|$ khỏi công thức bước lặp.

3. Bài toán trùng tách (SEP) và trùng tách đa tập (MSSEP)

  • Byrne và Moudafi (2013): Lần đầu tiên nghiên cứu bài toán SEP trong không gian Hilbert hữu hạn chiều thông qua việc cực tiểu hóa phiếm hàm khoảng cách $f(x,y) = \frac{1}{2}|Ax - By|^2$ trên tập tích $C \times Q$.
  • Chen và cộng sự (2013, 2014, 2016), Shi và cộng sự (2014), Tian và cộng sự (2014): Đề xuất phương pháp hiệu chỉnh Tikhonov cho bài toán trùng tách xấp xỉ (ASEP) và MSSEP cho họ vô hạn tập hợp. Tồn tại chính là việc mỗi bước lặp đòi hỏi tính toán chuỗi vô hạn của các phép chiếu metric.

Vị trí và khoảng trống luận án giải quyết

Luận án định vị trọng tâm vào việc khắc phục đồng thời các hạn chế nêu trên:

  1. Loại bỏ hoàn toàn sự phụ thuộc vào chuẩn $|A|$ trong công thức tham số bước lặp của phương pháp hiệu chỉnh lặp cho bài toán MSSFP.
  2. Thiết kế giải thuật lặp giải bài toán MSSEP với họ vô hạn tập hợp sao cho mỗi bước lặp chỉ xử lý tổng hữu hạn có chỉ số tăng dần theo số bước lặp, loại bỏ hoàn toàn việc tính chuỗi vô hạn.
  3. Xây dựng sơ đồ lai ghép Ishikawa - đường dốc nhất cho bài toán VIP trên không gian Banach, thiết lập hội tụ mạnh cho họ vô hạn ánh xạ không giãn.

Cơ sở lý thuyết và phương pháp nghiên cứu

1. Khung lý thuyết và khái niệm sử dụng

Luận án vận dụng các khái niệm và định lý cơ bản trong giải tích hàm và giải tích phi tuyến:

  • Tính chất toán tử trong không gian Hilbert $H$: Ánh xạ $L$-Lipschitz, ánh xạ không giãn ($L=1$), ánh xạ co ($L<1$), ánh xạ $\lambda$-đơn điệu mạnh, ánh xạ $\mu$-ngược đơn điệu mạnh, ánh xạ $\alpha$-trung bình ($T = (1-\alpha)I + \alpha T'$ với $T'$ không giãn), ánh xạ không giãn chặt ($T = \frac{1}{2}(I+V)$).
  • Phép chiếu metric $P_C$: Tính chất của phép chiếu trên tập con lồi đóng: $\langle x - P_C x, y - P_C x \rangle \le 0, \forall y \in C$ và $|u - P_C u|^2 + |P_C u - z|^2 \le |u - z|^2, \forall z \in C$.
  • Không gian Banach $E$ và ánh xạ đối ngẫu:
    • Môđun trơn $\rho_E(\tau) = \sup{\frac{1}{2}(|x+y| + |x-y|) - 1 : |x| \le 1, |y| \le \tau}$. Không gian $E$ là $q$-trơn đều nếu $\rho_E(\tau) \le c\tau^q$ ($c>0, 1<q\le 2$), tiêu biểu là các không gian $L^p, l^p, W_m^p$ với $1<p<\infty$.
    • Ánh xạ đối ngẫu chuẩn tắc $j: E \to 2^{E^}$ xác định bởi $j(x) = {x^ \in E^* : \langle x, x^* \rangle = |x|^2, |x^| = |x|}$. Trong không gian Hilbert, $j \equiv I$. Ánh xạ đối ngẫu tổng quát $j_q(x) = {x^ \in E^* : \langle x, x^* \rangle = |x|^q, |x^*| = |x|^{q-1}}$.
  • Lý thuyết hiệu chỉnh: Phương pháp hiệu chỉnh Lavrentiev ($Ax + \epsilon x = f$), phương pháp hiệu chỉnh Tikhonov, và phương pháp hiệu chỉnh lặp kiểu Bruck - Bakushinsky.

2. Phương pháp nghiên cứu

  • Phương pháp phân tích toán học thuần túy: Thiết lập các đánh giá bất đẳng thức biến phân, sử dụng các bổ đề giải tích dãy thực (bổ đề Xu, bổ đề Maingé, bổ đề dãy số thực không âm) để chứng minh tính bị chặn và tính hội tụ mạnh của dãy nghiệm xấp xỉ trong cấu trúc không gian vô hạn chiều.
  • Phương pháp xấp xỉ hữu hạn hóa: Kỹ thuật cắt cụt chỉ số tổng vô hạn thành tổng hữu hạn có trọng số điều chỉnh $\bar{\gamma}k = \sum{i=1}^k \gamma_i$ và $\bar{\lambda}k = \sum{j=1}^k \lambda_j$ tại mỗi bước lặp $k$.
  • Phương pháp thực nghiệm tính toán số: Lập trình giải thuật trên máy tính, chạy thử nghiệm trên các bài toán mẫu (test problems) với các ma trận và tập ràng buộc cụ thể để đo lường sai số $|x_{k+1} - x_k|$, vẽ đồ thị hội tụ và đối sánh hiệu năng.

3. Nguồn tư liệu và dữ liệu mô hình hóa

Luận án sử dụng hai mô hình bài toán ứng dụng thực tế để làm cơ sở dẫn xuất dữ liệu:

  1. Bài toán chụp cắt lớp vi tính X-quang (CT Scanner): Dựa trên định luật hấp thụ tổng quát Ber $dI(x) = -\mu(x)I(x)dx$, tích phân đường của hệ số suy giảm tuyến tính $p(v,\theta) = -\ln(I/I_0) = \int_{-\infty}^{+\infty} \mu(x,y)du$. Quá trình rời rạc hóa dẫn về hệ phương trình đại số tuyến tính $Ax = b$ với nhiễu $|w| \le \epsilon$, chuyển hóa thành bài toán chấp nhận tách kết hợp hiệu chỉnh ràng buộc chuẩn $l_1$: $\min_{x \in E^n} \frac{1}{2}|Ax - b|_2^2$ với $|x|_1 \le t$.
  2. Bài toán tối ưu hóa điều biến liều xạ trị (IMRT): Phân chia thể tích cơ thể thành $m$ phần tử thể tích (voxel) và $n$ chùm tia nhỏ (beamlet). Ma trận liều hấp thụ $H = (d_{ij}){m \times n}$. Các ràng buộc tập lồi gồm: vùng cường độ chùm tia không âm $X+ = {x \in E^n : x_j \ge 0}$, giới hạn liều cực đại tại các tổ chức rủi ro OAR ($H_{\max,t}$), giới hạn liều cực tiểu tại thể tích đích PTV ($H_{\min,t}$) và hàm liều tương đương $E_t(h) = (\frac{1}{N_t}\sum_{i \in S_t} h_i^{\sigma_t})^{1/\sigma_t}$, tạo thành bài toán MSSFP dạng tìm $p^* \in X_+ \cap \bigcap_l X_l$ sao cho $Hp^* \in \bigcap_v H_v$.

Nội dung chính theo từng chương

CẤU TRÚC NỘI DUNG LUẬN ÁN:

Chương 1: Một số kiến thức bổ trợ

Chương 1 hệ thống hóa các định nghĩa toán tử, các bổ đề giải tích cơ sở và mô hình hóa toán học các ứng dụng thực tế.

Nội dung chi tiết:

  • Toán tử và phép chiếu: Trình bày cấu trúc không gian Hilbert $H$, không gian Banach $E$, các lớp toán tử không giãn, trung bình, ngược đơn điệu mạnh, ánh xạ đối ngẫu $j, j_q$.
  • Các thuật toán điểm bất động cổ điển: Phân tích cơ chế và điều kiện hội tụ của phương pháp lặp Mann (1.4), Ishikawa (1.5), Halpern (1.7), xấp xỉ mềm Moudafi (1.10), cùng các sơ đồ cải biên kết hợp giữa Mann - Halpern (Kim & Xu 1.12), Krasnosel'skii-Mann biến tính (Yao và cộng sự 1.13), và Ishikawa - Halpern (Qin và cộng sự 1.14).
  • Mô hình bài toán MSSFP và MSSEP:
    • Bài toán MSSFP: Tìm $x \in C := \bigcap_{i \in J_1} C_i$ sao cho $Ax \in Q := \bigcap_{j \in J_2} Q_j$.
    • Bài toán MSSEP: Tìm $z = [x, y]$ với $x \in C := \bigcap_{i \in J_1} C_i, y \in Q := \bigcap_{j \in J_2} Q_j$ sao cho $Ax = By$.
    • Phân tích giải thuật tự thích nghi của Tian - Zhang (1.18), giải thuật hiệu chỉnh của Nguyễn Bường và cộng sự (1.20) và phương pháp hiệu chỉnh Tikhonov cho ASEP của Chen và cộng sự (1.31).
  • Phân tích toán học hai bài toán ứng dụng:
    • Thiết lập phương trình chiếu Radon rời rạc trong chụp cắt lớp X-quang $A x = b$ và bài toán tối ưu với ràng buộc chuẩn $l_1$.
    • Thiết lập mô hình bài toán MSSFP trong tối ưu hóa xạ trị IMRT dựa trên ma trận phân bố liều $H = (d_{ij})$ và các ràng buộc phi tuyến lồi $H_{\max,t}, H_{\min,t}, \Omega_t, \Gamma_t$.

Chương 2: Phương pháp hiệu chỉnh lặp xấp xỉ nghiệm bài toán chấp nhận tách và trùng tách đa tập

Chương 2 tập trung giải quyết Mục tiêu 1 và Mục tiêu 2 của luận án thông qua hai phương pháp hiệu chỉnh lặp mới.

1. Bài toán chấp nhận tách đa tập (MSSFP) – Phương pháp hiệu chỉnh kiểu Lavrentiev

Để giải bài toán MSSFP với $J_1 = \mathbb{N}^+, J_2 = \mathbb{N}^+$, tác giả đề xuất phương trình hiệu chỉnh kiểu Lavrentiev: $$F^k u^k + \epsilon_k (u^k - x^+) = 0 \quad (2.3)$$ Trong đó $x^+ \in H_1$ là điểm dự đoán ban đầu, $\epsilon_k > 0$ là tham số hiệu chỉnh, toán tử $F^k$ được định nghĩa: $$F^k = I - U^k - A^*(I - V^k)A$$ với các toán tử trung bình lấy trên tổng hữu hạn: $$U^k = \frac{1}{\bar{\gamma}k}\sum{i=1}^k \gamma_i P_{C_i}, \quad V^k = \frac{1}{\bar{\lambda}k}\sum{j=1}^k \lambda_j P_{Q_j}, \quad \bar{\gamma}k = \sum{i=1}^k \gamma_i, \quad \bar{\lambda}k = \sum{j=1}^k \lambda_j$$

Từ phương trình hiệu chỉnh, sơ đồ lặp hiện được xây dựng: $$z^{k+1} = \left(I - \mu_k (F^k + \epsilon_k (I - x^+))\right) z^k, \quad k \ge 1 \quad (2.11)$$

Các điều kiện áp đặt lên hệ tham số:

  • (a) $\mu_k, \epsilon_k \in (0, 1)$, $\lim_{k\to\infty} \mu_k/\epsilon_k = 0$, $\lim_{k\to\infty} \epsilon_k = 0$, $\epsilon_{k+1} < \epsilon_k$, và $\sum_{k=1}^\infty \mu_k \epsilon_k = \infty$.
  • (b) $\gamma_i > 0, \lambda_j > 0$ thỏa mãn $\sum_{i=1}^\infty \gamma_i = 1$ và $\sum_{j=1}^\infty \lambda_j = 1$.
  • Dãy tham số mẫu được chứng minh thỏa mãn toàn bộ giả thiết: $\mu_k = \frac{1}{(k+1)^a}$, $\epsilon_k = \frac{1}{(k+1)^b}$ với $0 < b < a$ và $a + 2b < 1$; $\gamma_i = \lambda_i = \frac{1}{i(i+1)}$.

Kết quả chứng minh (Định lý 2.1 & Định lý 2.2):

  • Dãy nghiệm hiệu chỉnh ${u^k}$ duy nhất và hội tụ mạnh về $p^* = P_\Gamma x^+$, tức là nghiệm của bài toán MSSFP thỏa mãn: $$|p^* - x^+| \le |p - x^+| \quad \forall p \in \Gamma$$
  • Dãy lặp ${z^k}$ sinh bởi thuật toán (2.11) hội tụ mạnh về $p^*$ khi $k \to \infty$. Bước lặp $\mu_k$ hoàn toàn độc lập với chuẩn $|A|$.

2. Bài toán trùng tách đa tập (MSSEP) – Phương pháp hiệu chỉnh kiểu Bakushinsky - Bruck

Tác giả mở rộng nguyên lý hiệu chỉnh kết hợp kỹ thuật tính toán tổng hữu hạn cho bài toán MSSEP: tìm $z = [x,y]$ với $x \in \bigcap_{i=1}^\infty C_i, y \in \bigcap_{j=1}^\infty Q_j$ sao cho $Ax = By$.

Cơ chế thuật toán:

  • Chuyển bài toán về dạng cực tiểu hàm phi tuyến trên không gian tích $H = H_1 \times H_2$.
  • Áp dụng phép chiếu metric trên các tổng hữu hạn $U^k$ và $V^k$ tại bước thứ $k$.
  • Chứng minh tính bị chặn và sự hội tụ mạnh của dãy nghiệm lặp về nghiệm có chuẩn cực tiểu của bài toán MSSEP mà không cần tính toán chuỗi vô hạn.
  • Cung cấp các ví dụ số kiểm chứng tính đúng đắn và tốc độ hội tụ của thuật toán.

Chương 3: Phương pháp lai ghép đường dốc nhất với phương pháp Ishikawa xấp xỉ nghiệm bài toán bất đẳng thức biến phân

Chương 3 giải quyết Mục tiêu 3 của luận án trong không gian Banach.

Nội dung chi tiết:

  • Phát biểu bài toán: Tìm $p^* \in C := \bigcap_{i=1}^\infty Fix(T_i)$ sao cho: $$\langle F p^, j(p^ - p) \rangle \le 0 \quad \forall p \in C$$ trong đó $E$ là không gian Banach thực trơn đều, $F: E \to E$ là toán tử $\kappa$-$j$-đơn điệu mạnh và $\mu$-giả co chặt ($\kappa + \mu > 1$), ${T_i}$ là họ vô hạn các ánh xạ không giãn trên $E$.
  • Cấu trúc thuật toán lai ghép: Kết hợp kỹ thuật xấp xỉ hai bước của Ishikawa với hướng giảm dốc nhất của toán tử $F$:
    • Bước phụ: tính toán hướng lặp trung gian dựa trên tổ hợp lồi có trọng số của ánh xạ không giãn.
    • Bước chính: điều chỉnh xấp xỉ nghiệm dọc theo hướng tác động của toán tử $F$ kết hợp với dãy hệ số co dần về 0.
  • Chứng minh toán học: Sử dụng các tính chất hình học của không gian Banach $q$-trơn đều và ánh xạ đối ngẫu $j$, tác giả chứng minh dãy lặp hội tụ mạnh về nghiệm duy nhất $p^*$ của bài toán VIP.
  • Ví dụ số minh họa: Luận án xây dựng các ví dụ số cụ thể kiểm tra sự hội tụ của dãy lặp trong không gian hàm, minh họa tốc độ suy giảm sai số qua các bước lặp.

Bảng tổng hợp so sánh các phương pháp lặp trong luận án

Tiêu chí Phương pháp MSSFP (Mục 2.1) Phương pháp MSSEP (Mục 2.2) Phương pháp VIP lai ghép (Chương 3)
Không gian nghiên cứu Hilbert thực $H_1, H_2$ (vô hạn chiều) Hilbert thực $H_1, H_2, H_3$ (vô hạn chiều) Banach thực $E$ ($q$-trơn đều, phản xạ)
Lớp bài toán Chấp nhận tách đa tập ($\bigcap C_i, \bigcap Q_j$) Trùng tách đa tập ($Ax = By$) Bất đẳng thức biến phân trên $\bigcap Fix(T_i)$
Nguyên lý cốt lõi Hiệu chỉnh Lavrentiev + Chiếu hữu hạn Hiệu chỉnh Bakushinsky - Bruck Lai ghép Ishikawa + Đường dốc nhất
Loại hội tụ Hội tụ mạnh (Strong convergence) Hội tụ mạnh (Strong convergence) Hội tụ mạnh (Strong convergence)
Sự phụ thuộc chuẩn $|A|$ Không phụ thuộc (đã loại bỏ) Phụ thuộc tham số hiệu chỉnh Không áp dụng toán tử chuyển
Xử lý họ tập vô hạn Tổng hữu hạn $\bar{\gamma}_k, \bar{\lambda}_k$ Tổng hữu hạn tại mỗi bước lặp Họ vô hạn ánh xạ không giãn
Nghiệm giới hạn Nghiệm $p^*$ gần điểm dự đoán $x^+$ nhất Nghiệm có chuẩn cực tiểu Nghiệm duy nhất của bài toán VIP

Kết quả và những đóng góp mới

Luận án đạt được các kết quả khoa học và đóng góp mới sau:

  1. Đóng góp lý thuyết cho bài toán chấp nhận tách đa tập (MSSFP):

    • Xây dựng thành công phương pháp hiệu chỉnh lặp kiểu Lavrentiev cho bài toán MSSFP trong không gian Hilbert thực vô hạn chiều với họ vô hạn đếm được các tập con lồi đóng.
    • Loại bỏ hoàn toàn sự phụ thuộc vào chuẩn của toán tử chuyển $|A|$ trong công thức xác định bước lặp $\mu_k$.
    • Chứng minh trọn vẹn định lý hội tụ mạnh của dãy lặp về nghiệm gần điểm dự đoán $x^+$ nhất, khắc phục lỗ hổng chứng minh trong công trình trước đây của Tian và Zhang.

  2. Đóng góp thuật toán cho bài toán trùng tách đa tập (MSSEP):

    • Đề xuất phương pháp hiệu chỉnh lặp kiểu Bakushinsky - Bruck giải bài toán MSSEP trong không gian Hilbert với các họ tập vô hạn đếm được.
    • Giải quyết triệt để khó khăn trong thực hành tính toán bằng cách thay thế chuỗi vô hạn của các phép chiếu metric bằng các tổng hữu hạn có chỉ số tăng dần theo bước lặp $k$.
    • Thiết lập các điều kiện tham số tường minh đảm bảo tính hội tụ mạnh của dãy lặp về nghiệm có chuẩn cực tiểu.

  3. Đóng góp cho bài toán bất đẳng thức biến phân (VIP) trong không gian Banach:

    • Xây dựng thành công sơ đồ lặp lai ghép giữa phương pháp Ishikawa và phương pháp đường dốc nhất để giải bài toán VIP trên tập điểm bất động chung của một họ vô hạn ánh xạ không giãn.
    • Nâng bậc hội tụ từ hội tụ yếu (của phương pháp Ishikawa truyền thống) lên hội tụ mạnh trong không gian Banach trơn đều.

  4. Đóng góp thực tiễn và tính toán:

    • Xây dựng các ví dụ số minh họa cụ thể cho từng thuật toán đề xuất ở Chương 2 và Chương 3, kiểm chứng bằng thực nghiệm số về tốc độ hội tụ và tính ổn định của giải thuật.
    • Hệ thống hóa mô hình toán học ứng dụng của bài toán MSSFP trong kỹ thuật 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ị IMRT.

Hạn chế và hướng nghiên cứu tiếp

Hạn chế của nghiên cứu

  • Các kết quả lý thuyết trong luận án chủ yếu thiết lập trên lớp ánh xạ tuyến tính bị chặn $A$ và lớp ánh xạ phi tuyến không giãn hoặc co. Chưa mở rộng cho các lớp toán tử phi tuyến tổng quát hơn như toán tử tựa không giãn (quasi-nonexpansive), toán tử giả co (pseudocontractive) không liên tục Lipschitz hoặc các toán tử đa trị tổng quát.
  • Các thử nghiệm số trong luận án dừng lại ở mức độ các ví dụ toán học mô phỏng trên không gian mẫu, chưa triển khai thử nghiệm trên các bộ dữ liệu y tế lâm sàng thực tế quy mô lớn của chụp cắt lớp hay xạ trị IMRT.

Hướng nghiên cứu tiếp theo

  • Mở rộng phương pháp hiệu chỉnh lặp loại bỏ chuẩn toán tử chuyển cho các bài toán cân bằng tách (Split Equilibrium Problems) và bài toán bao hàm thức biến phân đơn điệu suy rộng.
  • Nghiên cứu các sơ đồ lặp quán tính (inertial extrapolation) kết hợp với kỹ thuật hiệu chỉnh Lavrentiev nhằm tăng tốc độ hội tụ của thuật toán.
  • Ứng dụng các thuật toán đề xuất vào bài toán xử lý ảnh số lớn, phân tích dữ liệu lớn và giải quyết bài toán điều biến liều xạ trị trên dữ liệu hình ảnh chụp cắt lớp bệnh nhân thực tế.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên khảo trong lĩnh vực Toán ứng dụng và Tối ưu hóa:

  • Dành cho nghiên cứu sinh, học viên cao học và giảng viên chuyên ngành Toán giải tích, Toán ứng dụng, Tối ưu hóa:
    • Tham khảo kỹ thuật xây dựng phương trình hiệu chỉnh Lavrentiev và Bakushinsky - Bruck để thu được hội tụ mạnh trong không gian vô hạn chiều.
    • Tham khảo phương pháp kỹ thuật xử lý xấp xỉ họ vô hạn các tập lồi đóng bằng chuỗi tổng hữu hạn tăng dần.
    • Tham khảo kỹ thuật phân tích bất đẳng thức trên không gian Banach $q$-trơn đều và các tính chất của ánh xạ đối ngẫu $j, j_q$.
  • Dành cho các nhà nghiên cứu trong lĩnh vực Xử lý tín hiệu, Thị giác máy tính và Vật lý y khoa (Xạ trị ung thư):
    • Tham khảo quy trình mô hình hóa bài toán ngược trong tái tạo ảnh chụp cắt lớp vi tính $Ax=b$ có nhiễu kết hợp ràng buộc chuẩn $l_1$.
    • Tham khảo mô hình toán học tối ưu hóa phân bố liều xạ trị IMRT với các ràng buộc lồi trên thể tích khối u (PTV) và cơ quan lành có nguy cơ (OAR).

Câu hỏi thường gặp

1. Điểm khác biệt cốt lõi giữa phương pháp hiệu chỉnh lặp Lavrentiev cho bài toán MSSFP trong luận án so với phương pháp của Xu (2010) là gì?

Trong phương pháp của Xu (2010), bước lặp $\mu_k$ bị giới hạn nghiêm ngặt bởi điều kiện $0 < \mu_k < \epsilon_k / (|A|^2 + \epsilon_k)$, nghĩa là bắt buộc phải tính toán hoặc ước lượng cận trên của chuẩn toán tử chuyển $|A|$. Luận án đã đề xuất phương pháp hiệu chỉnh kiểu Lavrentiev với điều kiện tham số $\mu_k, \epsilon_k \in (0,1)$, $\lim_{k\to\infty} \mu_k/\epsilon_k = 0$ và $\sum \mu_k \epsilon_k = \infty$ (ví dụ $\mu_k = 1/(k+1)^a, \epsilon_k = 1/(k+1)^b$ với $0<b<a, a+2b<1$), hoàn toàn không chứa đại lượng $|A|$ trong biểu thức xác định bước lặp.

2. Luận án giải quyết trở ngại tính toán tổng vô hạn trong bài toán trùng tách đa tập (MSSEP) như thế nào?

Đối với họ vô hạn các tập con lồi đóng ${C_i}{i=1}^\infty$ và ${Q_j}{j=1}^\infty$, các nghiên cứu trước đây (như Chen và cs 2013) yêu cầu tính toàn bộ chuỗi vô hạn các phép chiếu tại mỗi bước lặp $k$. Luận án khắc phục bằng cách sử dụng toán tử trung bình trên tổng hữu hạn: $$U^k = \frac{1}{\bar{\gamma}k}\sum{i=1}^k \gamma_i P_{C_i}, \quad V^k = \frac{1}{\bar{\lambda}k}\sum{j=1}^k \lambda_j P_{Q_j}$$ với $\bar{\gamma}k = \sum{i=1}^k \gamma_i, \bar{\lambda}k = \sum{j=1}^k \lambda_j$. Tại mỗi bước $k$, thuật toán chỉ cần tính tổng từ $1$ đến $k$, giúp khối lượng tính toán ở mỗi bước là hữu hạn và khả thi trong thực tế.

3. Tại sao tác giả kết hợp phương pháp lặp Ishikawa với phương pháp đường dốc nhất trong Chương 3?

Phương pháp lặp Ishikawa truyền thống chỉ đảm bảo sự hội tụ yếu của dãy lặp trong không gian vô hạn chiều. Bằng cách kết hợp sơ đồ hai bước của Ishikawa với hướng giảm dốc nhất của toán tử $\kappa$-$j$-đơn điệu mạnh và $\mu$-giả co chặt $F$, thuật toán tạo ra lực kéo hiệu chỉnh đủ mạnh để đảm bảo dãy lặp hội tụ mạnh về nghiệm duy nhất của 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. Bài toán chụp cắt lớp vi tính (CT) được mô hình hóa thành bài toán chấp nhận tách như thế nào trong luận án?

Dựa trên định luật hấp thụ Ber $dI(x) = -\mu(x)I(x)dx$, bài toán khôi phục hệ số suy giảm tuyến tính $\mu(x,y)$ từ các hình chiếu $p(v,\theta)$ sau khi rời rạc hóa dẫn về hệ phương trình đại số $Ax = b$. Có tính đến sai số đo lường $|w| \le \epsilon$ và điều kiện giảm nhiễu dựa trên chuẩn $l_1$, bài toán trở thành bài toán hiệu chỉnh: tìm $\min_{x \in E^n} \frac{1}{2}|Ax - b|_2^2$ thỏa mãn ràng buộc $|x|_1 \le t$. Đây chính là bài toán chấp nhận tách (SFP) với tập lồi $C = {x \in E^n : |x|_1 \le t}$ và $Q = {b}$.

5. Dãy nghiệm của phương pháp hiệu chỉnh lặp đề xuất cho MSSFP hội tụ về điểm nào trong tập nghiệm?

Theo Định lý 2.1 và Định lý 2.2 của luận án, khi $k \to \infty$, dãy nghiệm hiệu chỉnh ${u^k}$ và dãy lặp hiện ${z^k}$ hội tụ mạnh về điểm $p^* \in \Gamma$ (với $\Gamma = C \cap A^{-1}Q$ là tập nghiệm của bài toán MSSFP). Điểm $p^$ này chính là hình chiếu metric của điểm dự đoán ban đầu $x^+$ lên tập nghiệm $\Gamma$, thỏa mãn đặc trưng khoảng cách cực tiểu: $|p^ - x^+| \le |p - x^+|$ với mọi $p \in \Gamma$.

Kết luận

Luận án tiến sĩ của tác giả Khuất Thị Bình đã giải quyết ba bài toán trung tâm trong lý thuyết điểm bất động và tối ưu hóa phi tuyến: bài toán chấp nhận tách đa tập (MSSFP), bài toán trùng tách đa tập (MSSEP) và bài toán bất đẳng thức biến phân (VIP). Bằng việc vận dụng các kỹ thuật hiệu chỉnh Lavrentiev, Bakushinsky - Bruck và lai ghép đường dốc nhất - Ishikawa, công trình đã đề xuất các sơ đồ lặp mới đạt được sự hội tụ mạnh trong không gian Hilbert và không gian Banach vô hạn chiều. Các đóng góp nổi bật bao gồm việc loại bỏ sự phụ thuộc vào chuẩn toán tử chuyển $|A|$, hữu hạn hóa các phép tính tổng vô hạn tại mỗi bước lặp và thiết lập các đánh giá hội tụ toán học chặt chẽ kèm ví dụ số minh họa.