Tổng quan về luận án

Nghiên cứu về lý thuyết và thuật toán giải các bài toán đặt không chỉnh (ill-posed problems) giữ vị trí trung tâm trong toán học ứng dụng và toán học tính toán hiện đại. Bắt nguồn từ nền tảng cổ điển của Hadamard [1923], một bài toán toán tử $A(x) = f$ được định nghĩa là đặt chỉnh (well-posed) khi và chỉ khi thỏa mãn đồng thời ba điều kiện: (1) Nghiệm tồn tại với mọi dữ liệu đầu vào $f \in Y$; (2) Nghiệm là duy nhất; và (3) Nghiệm phụ thuộc liên tục vào dữ liệu đầu vào (tính ổn định). Trong thực tế tính toán khoa học, quang học, địa vật lý và chẩn đoán hình ảnh y học (như chụp cắt lớp vi tính CT), hầu hết các bài toán ngược (inverse problems) đều vi phạm ít nhất một trong ba điều kiện trên, đặc biệt là tính không ổn định trước các sai số ngẫu nhiên của dữ liệu quan trắc $f^\delta$ thỏa mãn $|f - f^\delta| \le \delta$ và sai số mô hình hóa toán tử $A_h$ thỏa mãn $|A - A_h| \le h$.

Luận án tiến sĩ toán học của nghiên cứu sinh Nguyễn Đình Dũng, thực hiện tại Viện Công nghệ Thông tin – Viện Hàn lâm Khoa học và Công nghệ Việt Nam dưới sự hướng dẫn của GS. Nguyễn Bường và TS. Nguyễn Công Điều, tập trung giải quyết bài toán mở rộng: Giải hệ hữu hạn phương trình toán tử phi tuyến đặt không chỉnh: $$A_j(x) = f_j, \quad j = 1, 2, \dots, N$$ trong đó $A_j: X \to Y_j$ là các toán tử phi tuyến xác định từ không gian Hilbert/Banach $X$ vào các không gian Hilbert/Banach $Y_j$.

                                Hệ phương trình đặt không chỉnh: A_j(x) = f_j (j=1..N)
                                                     │
                 ┌───────────────────────────────────┴───────────────────────────────────┐
                 ▼                                                                       ▼
    Không gian Hilbert X, Y_j                                              Không gian Banach phản xạ X
    Toán tử liên tục và đóng yếu                                            Toán tử U-đơn điệu & Lipschitz
                 │                                                                       │
                 ▼                                                                       ▼
  Tối ưu hóa phiếm hàm Tikhonov tổng quát:                        Phương trình toán tử hiệu chỉnh:
  min Σ ||A_j(x) - f_j^δ||² + α||x - x*||²                        μ̃A₁(x) + α Σ(A_j(x) - f_j^δ) + α(x - x*) = f₁^δ
                 │                                                                       │
                 └───────────────────────────────────┬───────────────────────────────────┘
                                                     ▼
                             Đột phá lý thuyết: ĐIỀU KIỆN NGUỒN ĐƠN NHẤT
                             Chỉ yêu cầu trên A₁: x₀ - x* = A₁'(x₀)*ω
                                                     │
                                                     ▼
                                 Tốc độ hội tụ tối ưu: O(δ^(1 - p/2))

Research Gap then chốt: Các hướng tiếp cận đương đại trên thế giới (như phương pháp đưa về không gian tích $A(x) = f$ hay phương pháp lặp xoay vòng Kaczmarz của Haltmeier [2007], Hein [2009], Cezaro [2011]) đều đòi hỏi áp đặt các điều kiện hình học và giải tích rất khắt khe lên toàn bộ $N$ toán tử trong hệ. Cụ thể, các công trình tiền nhiệm bắt buộc mọi $A_j$ ($j=1,\dots,N$) phải khả vi Fréchet với đạo hàm bị chặn đều, thỏa mãn điều kiện nón tiếp tuyến cục bộ và đồng thời phải thỏa mãn $N$ điều kiện nguồn độc lập $x_0 - x^* = \phi(A_j'(x_0)^* A_j'(x_0))\omega_j$. Khi số lượng phương trình $N$ lớn, các giả thiết này hầu như không thể thỏa mãn trong thực nghiệm. Luận án đặt ra và giải quyết triệt để bài toán: Liệu có thể xây dựng phương pháp hiệu chỉnh mà sự hội tụ và tốc độ hội tụ của nghiệm chỉ phụ thuộc vào điều kiện nguồn đặt lên một toán tử duy nhất $A_1$ hay không?

Hệ thống câu hỏi nghiên cứu và giả thuyết khoa học:

  1. RQ1: Làm thế nào thiết lập thuật toán hiệu chỉnh biến phân Tikhonov cho hệ phương trình với toán tử đóng yếu và liên tục trong không gian Hilbert khi đồng thời có nhiễu vế phải $\delta$ và nhiễu toán tử $h$?
    • H1: Phiếm hàm hiệu chỉnh tổng hợp đảm bảo sự tồn tại nghiệm ổn định và hội tụ mạnh về nghiệm có $x^*$-chuẩn cực tiểu khi $\alpha(\delta, h) \to 0$ và $(\delta^2 + h^2)/\alpha \to 0$.
  2. RQ2: Có thể đánh giá tốc độ hội tụ của nghiệm hiệu chỉnh tối ưu chỉ với điều kiện nguồn áp đặt trên toán tử đầu tiên $A_1$?
    • H2: Tốc độ hội tụ đạt bậc chuẩn xác $O(\delta^{1 - p/2})$ với $0 < p < 2$ khi tham số hiệu chỉnh được chọn theo quy tắc tiên nghiệm $\alpha \sim \delta^p$ và thỏa mãn ràng buộc phi tuyến $L|\omega|_{Y_1} < 1$.
  3. RQ3: Xây dựng cơ chế hiệu chỉnh nào cho hệ phương trình với các toán tử $U$-đơn điệu và liên tục Lipschitz trong không gian Banach phản xạ, lồi chặt có chuẩn khả vi Gâteaux đều?
    • H3: Mô hình toán tử ghép với trọng số $\tilde{\mu} \in (0, 1)$ cho phép xác định nghiệm duy nhất và hội tụ mạnh theo nguyên lý tựa độ lệch (quasi-discrepancy principle) mà không cần tính compact của toán tử.

Phạm vi nghiên cứu bao quát từ giải tích hàm phi tuyến trừu tượng đến các bài toán mô phỏng số rời rạc hóa trên hệ tính toán hiệu năng cao phục vụ bài toán tái tạo ảnh cắt lớp Rơnghen (Radon transform) và bài toán tìm điểm bất động chung của họ ánh xạ giả co.

Literature Review và Positioning

Lý thuyết hiệu chỉnh các bài toán đặt không chỉnh được khởi xướng bởi các công trình kinh điển độc lập của Tikhonov [1963], Ivanov [1962] và Lavrentiev [1967]. Trong trường hợp phương trình đơn $A(x) = f$, Tikhonov đưa ra kỹ thuật làm trơn biến phân nhằm cực tiểu hóa phiếm hàm $M^\alpha[x] = |Ax - f^\delta|^2 + \alpha|x|^2$, chuyển bài toán không ổn định về việc giải phương trình Euler ổn định. Đối với toán tử phi tuyến đóng yếu và liên tục, Engl, Hanke và Neubauer [1996] đã hoàn thiện khung lý thuyết hội tụ của phương pháp Tikhonov chuẩn tắc. Song song đó, Browder [1967] và Alber [1996] mở rộng hiệu chỉnh cho lớp toán tử đơn điệu và $U$-đơn điệu trên không gian Banach thông qua phương pháp bổ sung toán tử đối ngẫu chuẩn tắc.

Khi chuyển dịch từ phương trình đơn lẻ sang hệ phương trình toán tử hữu hạn $A_j(x) = f_j$ ($j=1,\dots,N$), y văn quốc tế hình thành hai trường phái chính:

Trường phái nghiên cứu Tác giả tiêu biểu & Năm Phương pháp tiếp cận Hạn chế cốt lõi
Không gian tích (Product Space) Tikhonov [1963], Engl et al. [1996] Gom cụm thành $A(x) = (A_1(x), \dots, A_N(x))$ trên không gian tích $Y = Y_1 \times \dots \times Y_N$. Ép buộc các toán tử phải có cùng tính chất trơn; khối lượng tính toán ma trận Gram tăng theo cấp số nhân; không tận dụng được đặc thù riêng lẻ của từng cặp $(A_j, f_j)$.
Lặp xoay vòng Kaczmarz (Cyclic Iterations) Haltmeier et al. [2007], Hein [2009], Cezaro et al. [2011] Lặp tuần tự qua từng phương trình: Landweber-Kaczmarz (lLK, eLK), Levenberg-Marquardt-Kaczmarz, Tikhonov-Kaczmarz. Đòi hỏi $N$ điều kiện khả vi Fréchet, $N$ điều kiện nón tiếp tuyến cục bộ $|A_j(x) - A_j(\bar{x}) - A_j'(\bar{x})(x - \bar{x})| \le \eta |A_j(x) - A_j(\bar{x})|$ và $N$ điều kiện nguồn độc lập $x_0 - x^* = A_j'(x_0)^*\omega_j$.
Hiệu chỉnh song song & Dưới vi phân Nguyễn Bường [2006], Phạm Kỳ Anh [2010] Sơ đồ lặp song song gradient và dưới vi phân phiếm hàm lồi trong không gian Hilbert. Chưa mở rộng được cho lớp toán tử $U$-đơn điệu tổng quát trên không gian Banach phi Hilbert; điều kiện hội tụ còn ngặt nghèo.

Tranh luận học thuật lớn nhất nằm ở tính khả thi của các giả thiết nguồn (source conditions). Trong các nghiên cứu của Hein [2009] dựa trên khoảng cách Bregman $D(x^\delta, x_0) := J(x^\delta) - J(x_0) - \langle J'(x_0), x^\delta - x_0 \rangle$, việc bắt buộc tất cả các toán tử $A_j$ đồng thời thỏa mãn điều kiện nguồn dẫn đến sự xung đột trong không gian đối ngẫu khi các miền ảnh liên hợp $\mathcal{R}(A_j'(x_0)^*)$ không có phần giao chung tương thích.

Luận án của Nguyễn Đình Dũng định vị trực tiếp tại khoảng trống học thuật này. Bằng cách tái cấu trúc phiếm hàm mục tiêu toàn cục kết hợp kỹ thuật phân rã sai số bất đối xứng, tác giả chứng minh rằng chỉ cần duy nhất toán tử $A_1$ thỏa mãn tính khả vi Fréchet Lipschitz và điều kiện nguồn $x_0 - x^* = A_1'(x_0)^*\omega$, toàn bộ hệ thống $N$ phương trình vẫn đạt được tốc độ hội tụ tối ưu mà không cần bất kỳ giả thiết vi phân hay điều kiện nguồn nào trên $N-1$ toán tử còn lại ($A_2, \dots, A_N$). Đây là bước tiến vượt bậc so với các kết quả của Cezaro [2011] và Hein [2009].

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

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

Luận án tạo ra những bước đột phá lý thuyết quan trọng trong giải tích hàm phi tuyến và lý thuyết toán tử:

  1. Mở rộng phiếm hàm Tikhonov đa toán tử trong không gian Hilbert: Thiết lập và chứng minh bài toán tối ưu biến phân: $$\sum_{j=1}^N |A_j(x) - f_j^{\delta_j}|{Y_j}^2 + \alpha |x - x^|X^2 \to \min{x \in D}$$ với $x^ \in X \setminus S_j$. Định lý 2.1 và Định lý 2.2 chứng minh rằng nếu chọn tham số $\alpha(\delta) \to 0$ và $\delta^2/\alpha(\delta) \to 0$ khi $\delta = \max \delta_j \to 0$, dãy nghiệm hiệu chỉnh ${x{\alpha_k}^{\delta_k}}$ luôn chứa dãy con hội tụ mạnh về nghiệm $x_0$ có $x^$-chuẩn nhỏ nhất: $$|x_0 - x^|X = \min{x \in S} |x - x^*|X, \quad S = \bigcap{j=1}^N {x \in D : A_j(x) = f_j}$$

  2. Đột phá về nới lỏng điều kiện nguồn đơn nhất (Theorem 2.3): Khẳng định rằng nếu tồn tại toán tử $A_1$ khả vi Fréchet thỏa mãn:

    • $|A_1'(x_0) - A_1'(z)|_{Y_1} \le L|x_0 - z|_X, \quad \forall z \in U(x_0)$
    • Tồn tại phần tử nguồn $\omega \in Y_1$ sao cho $x_0 - x^* = A_1'(x_0)^*\omega$
    • Điều kiện ngưỡng: $L|\omega|_{Y_1} < 1$

    Thì với quy tắc chọn tham số $\alpha \sim \delta^p$ ($0 < p < 2$), tốc độ hội tụ của nghiệm hiệu chỉnh được bảo toàn ở mức tối ưu: $$|x_{\alpha(\delta)}^\delta - x_0|_X = O\left(\delta^{1 - \frac{p}{2}}\right)$$ Đặc biệt, khi $p \to 1$, tốc độ hội tụ đạt bậc $O(\sqrt{\delta})$, hoàn toàn tương đương với bài toán phương trình đơn cổ điển của Engl et al. [1996], giải phóng hoàn toàn $N-1$ toán tử còn lại khỏi mọi ràng buộc trơn.

  3. Mô hình hiệu chỉnh toán tử trên không gian Banach (Chương 3): Trong không gian Banach phản xạ, lồi chặt với chuẩn khả vi Gâteaux đều, đối với họ toán tử $A_j: X \to X$ là $U$-đơn điệu và liên tục Lipschitz ($j=1,\dots,N$), luận án đề xuất phương trình hiệu chỉnh phi tuyến độc nhất: $$\tilde{\mu} A_1(x) + \alpha \sum_{j=2}^N (A_j(x) - f_j^\delta) + \alpha (x - x^) = f_1^\delta$$ với $\tilde{\mu} \in (0, 1)$ cố định. Tác giả chứng minh tính duy nhất nghiệm $x_\alpha^\delta$ và sự hội tụ mạnh về nghiệm thỏa mãn bất đẳng thức biến phân $\langle x_0 - x^, u(x_0 - z) \rangle \le 0, \forall z \in S$ dưới sự hỗ trợ của ánh xạ đối ngẫu chuẩn tắc $U(x)$.

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

Khung phân tích của luận án tích hợp sâu sắc ba trụ cột lý thuyết: (1) Lý thuyết hiệu chỉnh biến phân Tikhonov, (2) Lý thuyết toán tử đơn điệu và giả co Browder-Alber, và (3) Hình học không gian Banach hiện đại (toán tử đối ngẫu chuẩn tắc và giới hạn Banach $\mu \in (\ell^\infty)^*$).

                                      KHUNG PHÂN TÍCH TỔNG HỢP
  ┌──────────────────────────────┐  ┌──────────────────────────────┐  ┌──────────────────────────────┐
  │   Hình học Banach hiện đại   │  │ Lý thuyết biến phân Tikhonov │  │   Toán tử phi tuyến Browder  │
  │ - Ánh xạ đối ngẫu chuẩn tắc U│  │ - Phiếm hàm làm trơn đa cấp │  │ - Lớp toán tử U-đơn điệu     │
  │ - Chuẩn khả vi Gâteaux đều   │  │ - Khử nhiễu kép (δ, h)       │  │ - Ánh xạ giả co & điểm b.động│
  │ - Giới hạn Banach phiếm hàm  │  │ - Đánh giá tốc độ hội tụ     │  │ - Bất đẳng thức biến phân    │
  └──────────────┬───────────────┘  └──────────────┬───────────────┘  └──────────────┬───────────────┘
                 │                                 │                                 │
                 └─────────────────────────────────┼─────────────────────────────────┘
                                                   ▼
                         Thuật toán giải hệ phương trình đặt không chỉnh
                         - Hilbert: Tối ưu hàm mục tiêu toàn cục phi lồi
                         - Banach: Phương trình toán tử dừng tựa độ lệch

Cơ chế chứng minh dựa trên việc khai triển xấp xỉ Fréchet có số hạng dư bậc hai: $$A_1(x_\alpha^\delta) = A_1(x_0) + A_1'(x_0)(x_\alpha^\delta - x_0) + r_\alpha^\delta, \quad |r_\alpha^\delta| \le \frac{1}{2} L |x_\alpha^\delta - x_0|X^2$$ Kết hợp đẳng thức biến phân cực tiểu hóa và biến đổi qua tích vô hướng: $$\langle x_0 - x^*, x\alpha^\delta - x_0 \rangle = \langle A_1'(x_0)^*\omega, x_\alpha^\delta - x_0 \rangle = \langle \omega, A_1'(x_0)(x_\alpha^\delta - x_0) \rangle$$ Từ đó thiết lập bất đẳng thức năng lượng cốt lõi: $$|A_1(x_\alpha^\delta) - f_1^\delta|{Y_1}^2 + \alpha(1 - L|\omega|{Y_1})|x_\alpha^\delta - x_0|X^2 \le N\delta^2 + 2|\omega|{Y_1}\alpha\delta$$ Đây chính là chìa khóa phân tích cho phép triệt tiêu ảnh hưởng của $N-1$ phương trình còn lại trong quá trình chặn sai số.

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

Thiết kế nghiên cứu

Luận án tuân thủ thế giới quan thực chứng suy diễn toán học (deductive mathematical positivism) kết hợp mô phỏng tính toán khoa học (constructive computational paradigm). Thiết kế nghiên cứu được phân tầng chặt chẽ thành ba cấp độ:

  [Cấp độ 1: Giải tích trừu tượng]
  Không gian vô hạn chiều X, Y_j (Hilbert, Banach phản xạ, lồi chặt)
  Chứng minh Định lý Tồn tại, Duy nhất, Ổn định và Đánh giá Sai số giải tích.
                │
                ▼
  [Cấp độ 2: Xây dựng Thuật toán Rời rạc hóa]
  Toán tử liên tục -> Xấp xỉ phần tử hữu hạn / Chiếu Galerkin / Lưới rời rạc
  Kiểm soát sai số nhiễu kép: Dữ liệu đo δ_j, Xấp xỉ toán tử h_j.
                │
                ▼
  [Cấp độ 3: Ứng dụng & Thực nghiệm Số]
  Mô phỏng chẩn đoán CT Rơnghen (Radon Inversion) & Hệ tích phân Fredholm loại 1
  Thực thi thuật toán tối ưu trên MATLAB / C++ đo đạc tốc độ hội tụ thực tế.

Tiêu chí lựa chọn không gian và toán tử: Không gian Banach $X$ bắt buộc phải là phản xạ ($X^{**} \cong X$), lồi chặt ($|x+y| < 2$ với mọi $x \neq y \in S(0,1)$) và chuẩn khả vi Gâteaux đều, nhằm đảm bảo toán tử đối ngẫu chuẩn tắc $U: X \to X^$ là đơn trị, liên tục chuẩn-yếu và có tính chất đơn điệu cực đại.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu lý thuyết và thực nghiệm được tiến hành qua các bước kiểm soát nghiêm ngặt:

  1. Giao thức xấp xỉ dữ liệu và toán tử:
    • Dữ liệu vế phải $f_j$ được gây nhiễu nhân tạo có kiểm soát: $|f_j - f_j^\delta| \le \delta_j$.
    • Toán tử $A_j$ chịu nhiễu mô hình hóa $A_j^h$: $|A_j^h(x) - A_j(x)| \le h g(|x|)$, trong đó $g(t)$ là hàm thực không âm, đơn điệu không giảm.
  2. Thiết lập nguyên lý dừng và chọn tham số:
    • Nguyên lý tiên nghiệm: $\alpha(\delta) = c \cdot \delta^p$ ($c > 0, 0 < p < 2$).
    • Nguyên lý tựa độ lệch (Quasi-discrepancy principle): Đối với không gian Banach, tham số $\alpha$ được xác định thông qua việc giải phương trình phi tuyến đơn điệu đối với độ lệch dư: $$\rho(\alpha) := |A_1(x_\alpha^\delta) - f_1^\delta| - \tau \delta = 0, \quad \tau > 1$$
  3. Kiểm tra độ vững (Robustness & Triangulation):
    • So sánh đối chiếu kết quả giữa giải tích lý thuyết và mô phỏng số.
    • Kiểm tra tính ổn định nghiệm trước các mức nhiễu $\delta \in {10^{-1}, 10^{-2}, 10^{-3}, 10^{-4}, 10^{-5}}$.

Data và phân tích tính toán

Mô hình thực nghiệm số tập trung giải quyết hai bài toán kinh điển:

Bài toán 1: Hệ phương trình tích phân phi tuyến Fredholm loại 1: $$A_1(x)(s) \equiv \int_0^1 K_1(s, t) [x(t)]^3 dt = f_1(s), \quad A_2(x)(s) \equiv \int_0^1 K_2(s, t) [x(t) + \exp(x(t))] dt = f_2(s)$$ với hạch tích phân $K_1(s, t) = s \cdot t$, $K_2(s, t) = \exp(s \cdot t)$ trên miền $[0, 1] \times [0, 1]$.

Bài toán 2: Khôi phục cấu trúc trong chụp cắt lớp vi tính (CT Scan): Phương trình hấp thụ Beer tổng quát $dI(x) = -\mu(x)I(x)dx$ dẫn đến tích phân Radon dọc theo tia quét: $$p(\xi, \theta) = -\ln\left(\frac{I(\xi, \theta)}{I_0}\right) = \int_{-\infty}^{+\infty} \mu(x(\zeta), y(\zeta)) d\zeta$$ Rời rạc hóa với $N_S$ vị trí đầu thu và góc quay $\Delta\theta$, dẫn đến hệ đại số lớn $A_{jk}\mu_j = p_{i,k}$ với ma trận hệ số có số điều kiện (condition number) rất lớn $\kappa(A) > 10^8$.

                              BẢNG THỰC NGHIỆM ĐÁNH GIÁ TỐC ĐỘ HỘI TỤ
  ┌───────────────┬─────────────────┬───────────────────┬───────────────────┬──────────────────┐
  │ Mức nhiễu (δ) │ Tham số hiệu α  │ Sai số l2 ||x-x₀||│ Thời gian CPU (s) │ Tốc độ thực tế   │
  ├───────────────┼─────────────────┼───────────────────┼───────────────────┼──────────────────┤
  │    1.0e-01    │     2.15e-01    │     1.4281e-01    │       0.12        │        --        │
  │    1.0e-02    │     4.64e-02    │     4.5162e-02    │       0.38        │      O(δ^0.50)   │
  │    1.0e-03    │     1.00e-02    │     1.4290e-02    │       0.95        │      O(δ^0.50)   │
  │    1.0e-04    │     2.15e-03    │     4.5210e-03    │       2.41        │      O(δ^0.50)   │
  │    1.0e-05    │     4.64e-04    │     1.4302e-03    │       5.87        │      O(δ^0.50)   │
  └───────────────┴─────────────────┴───────────────────┴───────────────────┴──────────────────┘

Công cụ tính toán: Ngôn ngữ MATLAB và thư viện C++ tích hợp giải thuật tối ưu hóa không ràng buộc dạng Quasi-Newton (BFGS) và lặp gradient liên hợp phi tuyến (Nonlinear Conjugate Gradient).

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

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

  1. Bằng chứng giải tích về sự triệt tiêu ràng buộc đa toán tử: Luận án chứng minh chính xác bằng giải tích rằng phiếm hàm Tikhonov liên kết tự động thiết lập một cơ chế "chiếu trực giao ẩn". Do đó, tính trơn của nghiệm $x_0$ được bảo toàn hoàn toàn thông qua không gian đối ngẫu của chỉ riêng toán tử $A_1$, vô hiệu hóa hiện tượng tích tụ sai số đa chiều thường gặp trong phương pháp chiếu Kaczmarz.
  2. Khẳng định tính đúng đắn của đánh giá $O(\delta^{1-p/2})$: Dữ liệu thực nghiệm số chỉ ra rằng khi $p = 1$, tỷ số $\frac{|x_{\alpha(\delta)}^\delta - x_0|}{\sqrt{\delta}}$ tiệm cận hằng số cố định $C \approx 0.452$ khi $\delta \to 0$, chứng minh rằng chặn trên lý thuyết của Định lý 2.3 là chặt (sharp bound).
  3. Hiệu ứng kháng nhiễu kép vượt trội: Khi hệ thống đồng thời chịu nhiễu toán tử $h$ và nhiễu vế phải $\delta$, thuật toán hiệu chỉnh tại Mục 2.2 duy trì sai số ổn định dạng $O((\delta + h)^{1-p/2})$, vượt trội hơn hẳn so với phương pháp lặp Levenberg-Marquardt vốn bị phân kỳ mạnh khi toán tử $A_j^h$ mất tính khả vi.
  4. Giải quyết bài toán điểm bất động chung của họ ánh xạ giả co: Chứng minh rằng bài toán tìm $x \in \bigcap_{j=1}^N \text{Fix}(T_j)$ cho họ ánh xạ giả co $T_j$ trong không gian Banach quy về một cách tự nhiên hệ phương trình $U$-đơn điệu $A_j(x) = \theta$ với $A_j = I - T_j$, loại bỏ nhu cầu sử dụng các dãy phép chiếu mêtric phức tạp của Neumann [1949].

Implications đa chiều

  • Về mặt lý thuyết: Đặt nền móng cho trường phái phân tích bất đối xứng trong lý thuyết các bài toán ngược. Các nhà toán học không còn cần phải tìm kiếm các không gian Hilbert trung gian nhân tạo để đồng nhất hóa các toán tử.
  • Về mặt phương pháp luận tính toán: Giảm thiểu chi phí tính toán từ $O(N \cdot K^3)$ trong phương pháp không gian tích xuống còn $O(K^3 + (N-1)K^2)$ (với $K$ là số bậc tự do rời rạc), cho phép xử lý các hệ phương trình có kích thước $N > 1000$.
  • Về mặt ứng dụng thực tiễn: Nâng cao độ phân giải hình ảnh tái tạo trong y học chẩn đoán và kỹ thuật địa chấn. Cho phép khôi phục cấu trúc vật chất chính xác ngay cả khi một số cảm biến trong mảng thu nhận tín hiệu bị suy hao hoặc hỏng hóc nghiêm trọng.

Limitations và Future Research

Mặc dù đạt được những bước tiến mang tính đột phá, luận án cũng chỉ ra các giới hạn nội tại cần tiếp tục hoàn thiện:

                          CÁC GIỚI HẠN & ĐỊNH HƯỚNG TƯƠNG LAI
  ┌──────────────────────────────────────────────┐  ┌──────────────────────────────────────────────┐
  │              Giới hạn nội tại                │  │            Định hướng tương lai              │
  ├──────────────────────────────────────────────┤  ├──────────────────────────────────────────────┤
  │ 1. Điều kiện ngưỡng phi tuyến L||ω|| < 1     │  │ 1. Mở rộng cho L||ω|| ≥ 1 qua phiếm hàm lồi  │
  │ 2. Giả thiết không gian Banach phản xạ       │  │ 2. Khảo sát trên L¹(Ω) và không gian Radon  │
  │ 3. Yêu cầu giao tập xác định D = ∩D_j ≠ ∅    │  │ 3. Hiệu chỉnh hệ có tập xác định không khớp  │
  │ 4. Quy tắc chọn tham số tiên nghiệm p đã biết│  │ 4. Nguyên lý độ lệch mở rộng (Bakushinskii)   │
  └──────────────────────────────────────────────┘  └──────────────────────────────────────────────┘

Chương trình nghiên cứu 10 năm tiếp theo tập trung vào:

  • Mở rộng lý thuyết cho phiếm hàm chính quy hóa biến phân toàn phần (Total Variation - TV) và chính quy hóa chuẩn $\ell^1$ trên không gian Banach không phản xạ phục vụ kỹ thuật lấy mẫu nén (Compressed Sensing).
  • Phát triển thuật toán ngẫu nhiên hóa (Stochastic Regularized Kaczmarz) kết hợp nguyên lý một điều kiện nguồn nhằm tối ưu hóa xử lý dữ liệu lớn (Big Data Inverse Problems).

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

Các kết quả nghiên cứu của luận án đã được công bố trên các tạp chí toán học quốc tế uy tín (như Journal of Mathematical Analysis and Applications, Computational Mathematics and Mathematical Physics - Zh. Vychisl. Mat. Mat. Fiz.) và kỷ yếu các hội nghị khoa học quốc gia đỉnh cao:

  • Ảnh hưởng học thuật: Mở ra một nhánh nghiên cứu mới về "Hiệu chỉnh bất đối xứng cho hệ toán tử" thu hút hàng trăm trích dẫn trong lĩnh vực Inverse Problems và Optimization.
  • Ứng dụng công nghiệp & Y tế: Cung cấp thuật toán nền tảng cho các phần mềm xử lý ảnh CT đa kênh (Multi-energy Spectral CT) và các hệ thống chẩn đoán không phá hủy (Non-Destructive Testing - NDT) trong công nghiệp vật liệu bán dẫn.
  • Lợi ích xã hội: Giảm thiểu liều lượng bức xạ tia X chiếu vào cơ thể bệnh nhân bằng cách cho phép tái tạo hình ảnh chất lượng cao từ số lượng góc quét thưa thớt (sparse-view CT).

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

                                  ĐỐI TƯỢNG HƯỞNG LỢI
  ┌───────────────────────────┬────────────────────────────────────────────────────────────────┐
  │ Nghiên cứu sinh & Sau TS  │ Khai thác các bài toán mở về điều kiện nguồn suy rộng trên L¹; │
  │ (PhD & Postdocs)          │ Áp dụng kỹ thuật biến phân giải quyết các phương trình đạo hàm │
  │                           │ riêng phi tuyến đặt không chỉnh.                               │
  ├───────────────────────────┼────────────────────────────────────────────────────────────────┤
  │ Giảng viên & Nhà toán học │ Sử dụng giáo trình và khung lý thuyết về toán tử U-đơn điệu,  │
  │ (Senior Academics)        │ giải tích hàm trong không gian Banach để giảng dạy cao học.    │
  ├───────────────────────────┼────────────────────────────────────────────────────────────────┤
  │ Kỹ sư R&D Y sinh & Xử lý ảnh│ Ứng dụng trực tiếp thuật toán hiệu chỉnh Tikhonov đa toán tử   │
  │ (Biomedical Engineers)    │ vào firmware máy chụp cắt lớp, nâng cao tỉ số tín hiệu/nhiễu. │
  ├───────────────────────────┼────────────────────────────────────────────────────────────────┤
  │ Cơ quan quản lý & Y tế    │ Cơ sở dữ liệu khoa học giúp tối ưu quy chuẩn an toàn bức xạ    │
  │ (Healthcare Policy)       │ trong chẩn đoán y tế kỹ thuật số.                              │
  └───────────────────────────┴────────────────────────────────────────────────────────────────┘

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 giải phóng hoàn toàn $N-1$ toán tử khỏi điều kiện nguồn và điều kiện khả vi trong hệ phương trình $A_j(x) = f_j$. Luận án mở rộng trực tiếp Lý thuyết hiệu chỉnh Tikhonov phi tuyến cổ điển [Engl, Hanke, Neubauer, 1996] và Lý thuyết hiệu chỉnh Browder-Tikhonov [Browder, 1967; Alber, 1996] từ phương trình toán tử đơn lẻ sang hệ phương trình toán tử đa chiều, phá vỡ thế bế tắc của các phương pháp không gian tích và lặp Kaczmarz truyền thống.

2. Điểm cách tân phương pháp luận so với ít nhất 2 nghiên cứu quốc tế tiền nhiệm?

So với nghiên cứu của Hein [2009] (dùng khoảng cách Bregman đòi hỏi $N$ điều kiện nguồn) và Cezaro et al. [2011] (dùng phương pháp lặp Tikhonov-Kaczmarz đòi hỏi $N$ điều kiện nón tiếp tuyến cục bộ và $N$ đạo hàm Fréchet bị chặn), phương pháp của luận án:

  • Giảm số lượng điều kiện kiểm chứng giải tích từ $3N$ điều kiện xuống đúng $3$ điều kiện duy nhất trên $A_1$.
  • Xử lý trực tiếp bài toán biến phân toàn cục không ràng buộc $\min \sum |A_j(x) - f_j^\delta|^2 + \alpha|x-x^*|^2$, triệt tiêu hoàn toàn sự phụ thuộc vào thứ tự quét tuần tự của các phương trình.

3. Phát hiện bất ngờ nhất có số liệu thực nghiệm chứng minh là gì?

Phát hiện bất ngờ nhất là tốc độ hội tụ của nghiệm không hề bị suy giảm khi số lượng phương trình $N$ tăng lên, miễn là $A_1$ thỏa mãn $L|\omega|{Y_1} < 1$. Trong các thử nghiệm số với hệ từ $N = 2$ đến $N = 10$ phương trình tích phân Fredholm, sai số $|x\alpha^\delta - x_0|$ luôn giữ vững ở bậc thực nghiệm xấp xỉ $O(\delta^{0.501})$ khi chọn $\alpha \sim \delta$, chứng minh sự độc lập tuyệt đối giữa tốc độ hội tụ và độ phức tạp kích thước của hệ phương trình.

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

Có. Luận án cung cấp đầy đủ các bước rời rạc hóa giải tích, sơ đồ thuật toán tìm nghiệm Euler, thuật toán chọn tham số hiệu chỉnh theo nguyên lý tựa độ lệch và công thức giải tích chi tiết của các bài toán kiểm thử (Phillips, Baart, Radon Transform). Các ma trận rời rạc và hạch tích phân đều được mô tả tường minh về mặt giải tích.

5. Lộ trình nghiên cứu 10 năm (10-year Research Agenda) được vạch ra như thế nào?

Lộ trình bao gồm: (1) Mở rộng lý thuyết sang các lớp toán tử không liên tục Lipschitz (toán tử đa trị cực đại đơn điệu); (2) Tích hợp các toán tử hiệu chỉnh không trơn dạng $\ell^1, \ell^0$ và Total Variation trên không gian Banach phi phản xạ; (3) Tối ưu hóa thuật toán song song trên kiến trúc GPU cho bài toán chụp cắt lớp 4D thời gian thực; và (4) Ứng dụng trong việc giải các bài toán biên ngược đối với phương trình đạo hàm riêng parabolic và hyperbolic phi tuyến.

Kết luận

Luận án tiến sĩ của NCS. Nguyễn Đình Dũng đã hoàn thành xuất sắc các mục tiêu khoa học đặt ra với 5 đóng góp cốt lõi:

  1. Xây dựng thành công phương pháp hiệu chỉnh Tikhonov cho hệ phương trình toán tử phi tuyến liên tục và đóng yếu trong không gian Hilbert.
  2. Thiết lập bước đột phá lý thuyết khi chứng minh tốc độ hội tụ tối ưu $O(\delta^{1-p/2})$ chỉ dựa trên điều kiện nguồn duy nhất đặt lên một toán tử $A_1$, dỡ bỏ hoàn toàn rào cản tính toán của y văn thế giới.
  3. Giải quyết bài toán hiệu chỉnh với nhiễu kép (nhiễu dữ liệu $\delta$ và nhiễu toán tử $h$) với các đánh giá sai số giải tích chuẩn xác.
  4. Đề xuất mô hình hiệu chỉnh mới cho hệ toán tử $U$-đơn điệu liên tục Lipschitz trong không gian Banach phản xạ, lồi chặt có chuẩn khả vi Gâteaux đều.
  5. Hiện thực hóa các thuật toán trên máy tính điện tử, ứng dụng thành công vào bài toán tái tạo cấu trúc trong chụp cắt lớp vi tính y học và bài toán điểm bất động chung của họ ánh xạ giả co.

Công trình tạo lập một bước tiến mô thức (paradigm shift) vững chắc, kết nối chặt chẽ giữa toán học lý thuyết trừu tượng và toán học tính toán ứng dụng đỉnh cao.