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 hiệu chỉnh giải hệ phương trình toán tử đặt không chỉnh" do nghiên cứu sinh Nguyễn Đình Dũng thực hiện tại Viện Công nghệ Thông tin thuộc Viện Hàn lâm Khoa học và Công nghệ Việt Nam, dưới sự hướng dẫn khoa học của GS. Nguyễn Bường và TS. Nguyễn Công Điều, thuộc chuyên ngành Toán học tính toán (mã số 62.01).
Tính cấp thiết và khoảng trống nghiên cứu
Trong mô hình hóa các bài toán thực tế và ứng dụng tính toán khoa học, tồn tại một lớp lớn các bài toán mà nghiệm không ổn định; nghĩa là những biến thiên nhỏ của dữ liệu đầu vào (do sai số đo đạc, nhiễu hoặc quá trình làm tròn số trên máy tính điện tử) dẫn tới sai lệch lớn ở dữ liệu đầu ra, thậm chí làm bài toán rơi vào trạng thái vô nghiệm theo nghĩa thông thường. Lớp bài toán này được gọi là bài toán đặt không chỉnh (ill-posed problems).
Khi bài toán mở rộng từ một phương trình đơn lẻ sang hệ gồm hữu hạn $N$ phương trình toán tử đặt không chỉnh dạng $A_j(x) = f_j$ ($j = 1, 2, \dots, N$), các hướng tiếp cận truyền thống như phương pháp lặp xoay vòng Landweber-Kaczmarz, Levenberg-Marquardt-Kaczmarz hoặc đưa về phương trình toán tử trên không gian tích gặp phải những rào cản kỹ thuật lớn khi $N$ tăng cao. Cụ thể, các nghiên cứu trước đây (như công trình của M. Haltmeier, T. Hein, A. De Cezaro) đều đòi hỏi phải áp đặt đồng thời ba điều kiện khắt khe lên toàn bộ $N$ toán tử $A_j$:
- Điều kiện khả vi Fréchet với các đạo hàm Fréchet giới nội đều trong lân cận nghiệm;
- Điều kiện nón tiếp tuyến cục bộ;
- Điều kiện nguồn (source condition) trên nghiệm của hệ đối với từng toán tử thành phần.
Khoảng trống nghiên cứu được luận án nhận diện là: việc đòi hỏi điều kiện nguồn và điều kiện giải tích trên toàn bộ các toán tử thành phần làm thu hẹp phạm vi áp dụng thực tế và gây phức tạp lớn cho việc đánh giá tốc độ hội tụ. Do đó, nhu cầu khoa học đặt ra là xây dựng các lược đồ hiệu chỉnh mới cho phép nới lỏng điều kiện, trong đó sự hội tụ và tốc độ hội tụ của nghiệm hiệu chỉnh chỉ cần dựa trên điều kiện của một toán tử duy nhất trong hệ.
Mục tiêu và câu hỏi nghiên cứu
Luận án tập trung giải quyết hai mục tiêu và hai câu hỏi nghiên cứu cụ thể sau:
- Mục tiêu 1: Xây dựng phương pháp hiệu chỉnh dạng Tikhonov giải hệ phương trình toán tử phi tuyến và tuyến tính với các toán tử liên tục và đóng yếu trên không gian Hilbert (xét cả hai trường hợp: chỉ có nhiễu vế phải và có đồng thời nhiễu vế phải lẫn nhiễu toán tử), sao cho sự hội tụ và đánh giá tốc độ hội tụ của nghiệm hiệu chỉnh chỉ cần thỏa mãn điều kiện nguồn đặt lên một toán tử thành phần $A_1$.
- Mục tiêu 2: Xây dựng và khảo sát phương pháp hiệu chỉnh dạng Browder-Tikhonov giải hệ phương trình toán tử phi tuyến đối với lớp toán tử $U$-đơn điệu và liên tục Lipschitz trên không gian Banach phản xạ, lồi chặt có chuẩn khả vi Gâteaux đều; đồng thời đề xuất nguyên lý tựa độ lệch để lựa chọn tham số hiệu chỉnh và đánh giá tốc độ hội tụ của nghiệm hiệu chỉnh dựa trên một toán tử $A_1$.
Hai câu hỏi nghiên cứu được văn bản luận án xác lập trực tiếp gồm:
- Câu hỏi 1: Liệu có thể xây dựng được phương pháp hiệu chỉnh cho hệ phương trình toán tử đặt không chỉnh mà sự hội tụ cũng như đánh giá tốc độ hội tụ của nghiệm hiệu chỉnh chỉ cần dựa vào điều kiện đặt lên một toán tử hay không?
- Câu hỏi 2: Khi các toán tử $A_j$ là toán tử $U$-đơn điệu trên không gian Banach, liệu có thể xây dựng được các phương pháp hiệu chỉnh hiệu quả tương tự như trường hợp phương trình đơn lẻ hay không?
Đối tượng và phạm vi nghiên cứu
- Đối tượng nghiên cứu: Hệ hữu hạn phương trình toán tử phi tuyến và tuyến tính đặt không chỉnh $A_j(x) = f_j$ ($j = 1, \dots, N$) với các lớp toán tử liên tục, đóng yếu trên không gian Hilbert, toán tử $U$-đơn điệu và liên tục Lipschitz trên không gian Banach.
- Phạm vi nghiên cứu: Nghiên cứu lý thuyết giải tích hàm phi tuyến, lý thuyết bài toán đặt không chỉnh, phương pháp biến phân, đánh giá tiên nghiệm/hậu nghiệm và mô phỏng thực nghiệm tính toán số trên máy tính.
Tổng quan tài liệu và vị trí của luận án
Các hướng nghiên cứu tiền đề
Luận án đã điểm lại lịch sử phát triển của lý thuyết bài toán đặt không chỉnh và các công trình nền tảng:
- Khái niệm bài toán đặt chỉnh: Do J. Hadamard (1923, 1932) đề xuất khi nghiên cứu phương trình vi phân elliptic và parabolic với ba tiêu chí: tồn tại nghiệm, nghiệm duy nhất, và nghiệm phụ thuộc liên tục vào dữ liệu ban đầu.
- Phương pháp hiệu chỉnh Tikhonov: A. N. Tikhonov đặt nền móng giải bài toán không chỉnh bằng cách đưa về bài toán cực tiểu phiếm hàm làm trơn trên tập compact lồi.
- Khái niệm tựa nghiệm: V. K. Ivanov đề xuất khái niệm tựa nghiệm (quasi-solution) làm cực tiểu khoảng cách $\inf_{x \in M} \rho_Y(A(x), f)$ trên tập compact $M$ khi dữ liệu vế phải bị nhiễu rơi ra ngoài miền ảnh $A(M)$.
- Phương pháp phương trình xấp xỉ: M. M. Lavrentiev đề xuất thay thế phương trình không chỉnh bằng phương trình xấp xỉ giải được và nghiệm phụ thuộc liên tục vào vế phải.
- Toán tử liên tục, đóng yếu và toán tử đơn điệu: H. W. Engl nghiên cứu phương trình với toán tử liên tục và đóng yếu trên không gian Hilbert; Ya. I. Alber mở rộng cho toán tử đơn điệu và hemi-liên tục từ không gian Banach vào không gian liên hợp. F. E. Browder phát triển phương pháp hiệu chỉnh cho toán tử đơn điệu.
- Bài toán điểm bất động chung: J. von Neumann (1949) đề xuất phương pháp lặp phép chiếu luân phiên tìm điểm chung trên giao của các tập lồi đóng; các nghiên cứu mở rộng của Phạm Kỳ Anh (hiệu chỉnh lặp song song), Nguyễn Bường (toán tử đơn điệu, $h$-liên tục, dưới vi phân phiếm hàm lồi).
- Hệ phương trình toán tử đặt không chỉnh: M. Haltmeier (2007) đưa ra phương pháp lặp cải tiến Landweber-Kaczmarz (gồm lLK và eLK); T. Hein (2009) sử dụng khoảng cách Bregman để đánh giá tốc độ hội tụ của hệ; A. De Cezaro (2011) nghiên cứu phương pháp lặp Tikhonov-Kaczmarz cho hệ phương trình.
| Hướng nghiên cứu / Tác giả | Đối tượng toán tử | Không gian | Nhược điểm / Rào cản kỹ thuật |
|---|---|---|---|
| A. N. Tikhonov, H. W. Engl | Phương trình đơn lẻ $A(x) = f$ | Hilbert | Chưa xét trực tiếp hệ phương trình độc lập |
| M. Haltmeier (2007) | Hệ phương trình $A_j(x) = f_j$ | Hilbert | Phương pháp lặp Kaczmarz phụ thuộc thứ tự lặp |
| T. Hein (2009) | Hệ phương trình $A_j(x) = f_j$ | Hilbert | Đòi hỏi điều kiện nguồn trên tất cả $N$ toán tử |
| A. De Cezaro (2011) | Hệ phương trình $A_j(x) = f_j$ | Hilbert | Đòi hỏi điều kiện nón tiếp tuyến và nguồn trên cả $N$ toán tử |
| Luận án (Nguyễn Đình Dũng) | Hệ phương trình $A_j(x) = f_j$ | Hilbert & Banach | Chỉ cần điều kiện nguồn trên 1 toán tử $A_1$; mở rộng cho toán tử $U$-đơn điệu |
Khoảng trống luận án lựa chọn giải quyết
Luận án tập trung khắc phục trực tiếp hạn chế của phương pháp đưa về không gian tích và phương pháp lặp xoay vòng Kaczmarz: loại bỏ yêu cầu phải thỏa mãn điều kiện nguồn và điều kiện trơn trên toàn bộ $N$ toán tử, thay vào đó chỉ cần một toán tử thành phần (chẳng hạn $A_1$) thỏa mãn điều kiện khả vi Fréchet và điều kiện nguồn dạng $x_0 - x^* = A_1'(x_0)^*\omega$.
Cơ sở lý thuyết và phương pháp nghiên cứu
Cơ sở lý thuyết và các khái niệm sử dụng
Luận án xây dựng trên nền tảng giải tích hàm hiện đại với các định nghĩa và không gian hàm sau:
- Không gian Hilbert và không gian Banach: Không gian Banach phản xạ, không gian lồi chặt, không gian có chuẩn khả vi Gâteaux và khả vi Gâteaux đều.
- Toán tử đối ngẫu chuẩn tắc: Toán tử đa trị $U: X \to 2^{X^}$ xác định bởi $U(x) = {u(x) \in X^ : \langle x, u(x)\rangle = |x||u(x)|, |u(x)| = |x|}$.
- Các lớp toán tử phi tuyến:
- Toán tử đóng yếu: Nếu $x_n \rightharpoonup x$ và $A(x_n) \rightharpoonup y$ thì $A(x) = y$.
- Toán tử $U$-đơn điệu: Thỏa mãn $\langle A(x) - A(y), u(x - y)\rangle \ge 0$ với mọi $x, y \in X$.
- Toán tử $U$-đơn điệu mạnh: Tồn tại $\alpha > 0$ sao cho $\langle A(x) - A(y), u(x - y)\rangle \ge \alpha |x - y|^2$.
- Toán tử ngược $U$-đơn điệu mạnh: Tồn tại $\lambda > 0$ sao cho $\langle A(x) - A(y), u(x - y)\rangle \ge \lambda |A(x) - A(y)|^2$.
- Toán tử $m$-$U$-đơn điệu: $A$ là $U$-đơn điệu và $R(A + \lambda I) = X$ với mọi $\lambda > 0$.
- Toán tử liên tục Lipschitz: Tồn tại hằng số $L > 0$ sao cho $|A(x) - A(y)| \le L|x - y|$.
- Ánh xạ giả co: Thỏa mãn $\langle Tx - Ty, u(x - y)\rangle \le |x - y|^2$; khi đó toán tử $A := I - T$ là toán tử $U$-đơn điệu.
- Giới hạn Banach: Phiếm hàm tuyến tính liên tục $\mu$ trên không gian $\ell^\infty$ thỏa mãn $|\mu| = \mu_k(1) = 1$ và $\mu_k(a_{k+1}) = \mu_k(a_k)$.
Phương pháp nghiên cứu
- Phương pháp giải tích biến phân: Chuyển bài toán tìm nghiệm hệ phương trình về bài toán tối ưu cực tiểu hóa phiếm hàm Tikhonov tổng hợp: $$\sum_{j=1}^N |A_j(x) - f_j^\delta|_{Y_j}^2 + \alpha |x - x^*|X^2 \to \min{x \in D}$$
- Phương pháp hiệu chỉnh toán tử: Thiết lập phương trình hiệu chỉnh phi tuyến trên không gian Banach: $$A_1(x) + \alpha \sum_{j=2}^N \tilde{\mu} (A_j(x) - f_j^\delta) + \alpha(x - x^*) = f_1^\delta, \quad \tilde{\mu} \in (0, 1)$$
- Kỹ thuật đánh giá tiên nghiệm và hậu nghiệm: Sử dụng các bất đẳng thức hàm, nguyên lý tựa độ lệch (quasi-residual principle) để xác định quy tắc chọn tham số hiệu chỉnh $\alpha = \alpha(\delta)$.
- Phương pháp tính toán số: Rời rạc hóa bài toán, xây dựng thuật toán lặp và lập trình thử nghiệm số.
Nguồn tư liệu và dữ liệu
Nghiên cứu sử dụng hệ thống tài liệu giải tích hàm, lý thuyết phương trình đạo hàm riêng, mô hình toán học của bài toán chụp cắt lớp vi tính bức xạ X (nguyên lý suy giảm cường độ theo định luật Beer-Lambert và biến đổi tích phân Radon), cùng các bộ dữ liệu mô phỏng tính toán số thực nghiệm.
Nội dung chính theo từng chương
Chương 1: Hệ phương trình toán tử đặt không chỉnh
Chương 1 đóng vai trò hệ thống hóa cơ sở lý thuyết, giới thiệu bài toán nguồn gốc thực tế và khảo sát các phương pháp hiệu chỉnh hiện có.
- Khái niệm cơ bản (Mục 1.1): Trình bày các định nghĩa về tích vô hướng, không gian Hilbert, không gian Banach, tính lồi chặt, tính khả vi Gâteaux, dưới vi phân của phiếm hàm lồi, dãy cực tiểu hóa, sự hội tụ yếu và hội tụ mạnh.
- Phương pháp hiệu chỉnh Tikhonov cho phương trình đơn lẻ (Mục 1.2):
- Phân tích ví dụ kinh điển của Hadamard về bài toán Cauchy cho phương trình Laplace chứng minh tính không ổn định của nghiệm.
- Trình bày phương pháp hiệu chỉnh Tikhonov cho toán tử tuyến tính liên tục qua phiếm hàm $M^\alpha[x] = |A_h x - f_\delta|^2 + \alpha |x|^2$ và phương trình Euler $(A_h^* A_h + \alpha I)x = A_h^* f_\delta$.
- Trình bày hiệu chỉnh Tikhonov cho toán tử phi tuyến liên tục và đóng yếu, chứng minh sự tồn tại dãy con hội tụ về nghiệm có $x^$-chuẩn nhỏ nhất (Định lý 1.3) và tốc độ hội tụ $|x_\alpha^{\delta,\eta} - x_0| = O(\sqrt{\delta})$ dưới điều kiện nguồn $x_0 - x^ = A'(x_0)^*\omega$ (Định lý 1.4).
- Trình bày phương pháp hiệu chỉnh Browder-Tikhonov cho toán tử $U$-đơn điệu trên không gian Banach qua phương trình $A(x) + \alpha(x - x^*) = f_\delta$ (Định lý 1.5).
- Hệ phương trình toán tử đặt không chỉnh và mô hình dẫn xuất (Mục 1.3):
- Mô hình điểm bất động chung: Chứng minh bài toán tìm điểm bất động chung cho một họ ánh xạ giả co $T_j$ ($j = 1, \dots, N$) quy về hệ phương trình toán tử $U$-đơn điệu $A_j(x) = \theta$ với $A_j := I - T_j$.
- Mô hình chụp cắt lớp X-quang: Phân tích cơ chế suy giảm cường độ tia X qua định luật Beer: $$dI(x) = -\mu(x)I(x)dx$$ Xây dựng biểu thức hình chiếu tích phân trên hệ tọa độ quay $O\zeta\xi$: $$p(\xi, \theta) = -\ln\left(\frac{I(\xi, \theta)}{I_0}\right) = \int_{-\infty}^{+\infty} \mu(x, y) d\zeta$$ Rời rạc hóa miền tích phân dẫn đến hệ phương trình đại số đặt không chỉnh $\sum_{k=1}^{N_s} A_{jk}\mu_j = p_{i,k}$.
- Đánh giá các phương pháp trước: Phân tích phương pháp lặp xoay vòng dạng Newton, phương pháp lặp Tikhonov-Kaczmarz của Cezaro (2011) với công thức dừng lặp (1.27), chỉ ra độ phức tạp tính toán và việc phụ thuộc điều kiện trên toàn bộ $N$ toán tử.
Chương 2: Hiệu chỉnh cho hệ phương trình với toán tử liên tục và đóng yếu
Chương 2 trình bày các đóng góp lý thuyết cốt lõi của luận án đối với hệ phương trình trên không gian Hilbert.
- Trường hợp nhiễu vế phải (Mục 2.1):
- Xét hệ $A_j(x) = f_j$ ($j = 1, \dots, N$) với $A_j: D \subset X \to Y_j$ liên tục, đóng yếu và dữ liệu nhiễu $|f_j^{\delta_j} - f_j|_{Y_j} \le \delta_j$.
- Thiết lập bài toán cực tiểu hóa phiếm hàm: $$\sum_{j=1}^N |A_j(x) - f_j^{\delta_j}|_{Y_j}^2 + \alpha |x - x^*|X^2 \to \min{x \in D}$$
- Định lý 2.1 (Sự ổn định): Chứng minh với $\alpha > 0$, mọi dãy nghiệm cực tiểu ${x_k}$ tương ứng với $f_j^{\delta_{j_k}} \to f_j^{\delta_j}$ đều chứa dãy con hội tụ mạnh về nghiệm cực tiểu của bài toán.
- Định lý 2.2 (Sự hội tụ): Chứng minh khi $\alpha(\delta) \to 0$ và $\delta^2/\alpha(\delta) \to 0$ khi $\delta \to 0$, dãy nghiệm hiệu chỉnh ${x_{\alpha(\delta)}^\delta}$ hội tụ mạnh về nghiệm $x_0$ có $x^*$-chuẩn nhỏ nhất của hệ phương trình.
- Định lý 2.3 (Tốc độ hội tụ): Luận án chứng minh rằng chỉ cần giả thiết toán tử $A_1$ khả vi Fréchet, đạo hàm Fréchet liên tục Lipschitz với hằng số $L$, thỏa mãn điều kiện nguồn $x_0 - x^* = A_1'(x_0)^*\omega$ và $L|\omega|{Y_1} < 1$, khi chọn tham số $\alpha \sim \delta^p$ ($0 < p < 2$), tốc độ hội tụ đạt được đánh giá: $$|x{\alpha(\delta)}^\delta - x_0|_X = O\left(\delta^{1 - \frac{p}{2}}\right)$$
- Trường hợp nhiễu vế phải và nhiễu toán tử (Mục 2.2): Mở rộng bài toán khi toán tử $A_j$ bị xấp xỉ bởi $A_j^h$ thỏa mãn $|A_j^h(x) - A_j(x)|_{Y_j} \le h g(|x|_X)$ với $g(t)$ là hàm không âm giới nội. Thiết lập sự hội tụ và đánh giá sai số dưới tác động đồng thời của mức nhiễu dữ liệu $\delta$ và mức nhiễu toán tử $h$.
- Hệ phương trình với toán tử tuyến tính liên tục (Mục 2.3): Khảo sát trường hợp $A_j$ là các toán tử tuyến tính liên tục, dẫn xuất phương trình Euler tương ứng và quy tắc chọn tham số hiệu chỉnh tuyến tính.
- Kết quả tính toán số (Mục 2.4):
- Thiết lập quy tắc dừng lặp cho thuật toán.
- Thử nghiệm trên các hệ phương trình toán tử tuyến tính và phi tuyến đặt không chỉnh cụ thể, đối sánh giữa nghiệm giải tích và nghiệm hiệu chỉnh số để kiểm chứng tính đúng đắn của lý thuyết.
Chương 3: Hiệu chỉnh tìm nghiệm cho hệ phương trình phi tuyến với toán tử U − đơn điệu và liên tục Lipschitz trên không gian Banach
Chương 3 mở rộng phạm vi nghiên cứu sang không gian Banach phản xạ và lồi chặt có chuẩn khả vi Gâteaux đều.
- Xây dựng phương pháp hiệu chỉnh (Mục 3.1):
- Đề xuất mô hình hiệu chỉnh giải phương trình toán tử phi tuyến tổng hợp: $$A_1(x) + \alpha \sum_{j=2}^N \tilde{\mu} (A_j(x) - f_j^\delta) + \alpha(x - x^*) = f_1^\delta$$ trong đó $\tilde{\mu} \in (0, 1)$ là hằng số cố định, $\alpha > 0$ là tham số hiệu chỉnh.
- Chứng minh sự tồn tại và duy nhất nghiệm hiệu chỉnh $x_\alpha^\delta$ trên không gian Banach.
- Nguyên lý tựa độ lệch chọn tham số hiệu chỉnh (Mục 3.2): Xây dựng nguyên lý tựa độ lệch (quasi-residual principle) dạng hậu nghiệm để xác định tham số $\alpha = \alpha(\delta)$ tự động từ dữ liệu nhiễu mà không cần biết trước nghiệm chính xác.
- Đánh giá tốc độ hội tụ (Mục 3.3): Chứng minh tính hội tụ mạnh của nghiệm hiệu chỉnh $x_\alpha^\delta \to x_0$ khi mức nhiễu $\delta \to 0$, đồng thời thiết lập đánh giá tốc độ hội tụ chỉ dựa trên điều kiện đặt lên một toán tử $A_1$.
- Kết quả tính toán số (Mục 3.4): Trình bày các ví dụ số trên không gian Banach, mô phỏng quá trình lặp và đánh giá sai số thực tế giữa nghiệm hiệu chỉnh và nghiệm giải tích.
Kết quả và những đóng góp mới
Đóng góp mới về mặt lý luận khoa học
- Nới lỏng điều kiện nguồn cho hệ toán tử trên không gian Hilbert: Luận án đã giải quyết câu hỏi nghiên cứu thứ nhất bằng việc xây dự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, đóng yếu; chứng minh nghiệm hiệu chỉnh hội tụ về nghiệm có $x^*$-chuẩn nhỏ nhất. Điểm mới căn bản là tốc độ hội tụ $|x_\alpha^\delta - x_0|_X = O(\delta^{1 - p/2})$ được thiết lập chỉ dựa trên điều kiện nguồn và tính trơn của một toán tử $A_1$, hoàn toàn giải phóng $N-1$ toán tử còn lại khỏi các ràng buộc ngặt nghèo về điều kiện nguồn và điều kiện nón tiếp tuyến.
- Mở rộng lý thuyết hiệu chỉnh trên không gian Banach: Luận án đã giải quyết câu hỏi nghiên cứu thứ hai khi đề xuất phương pháp hiệu chỉnh dạng Browder-Tikhonov cho hệ phương trình với các toán tử $U$-đơn điệu và liên tục Lipschitz trên không gian Banach phản xạ, lồi chặt có chuẩn khả vi Gâteaux đều; xây dựng nguyên lý tựa độ lệch chọn tham số hiệu chỉnh và chứng minh sự hội tụ mạnh cùng tốc độ hội tụ của phương pháp.
- Tổng quát hóa cho trường hợp đa nhiễu: Thiết lập hoàn chỉnh cơ sở lý thuyết hội tụ cho hệ phương trình khi có đồng thời cả nhiễu vế phải $\delta$ lẫn nhiễu xấp xỉ toán tử $h$.
Kết quả và giải pháp thực tiễn
- Thuật toán tính toán số: Cung cấp quy tắc dừng lặp và thuật toán giải số cho cả hệ phương trình toán tử tuyến tính và phi tuyến đặt không chỉnh.
- Khả năng ứng dụng mô hình: Kết nối lý thuyết hiệu chỉnh hệ phương trình với bài toán thực tế về tái tạo hình ảnh cấu trúc vật thể trong chụp cắt lớp vi tính X-quang và bài toán tìm điểm bất động chung của họ ánh xạ giả co.
Hạn chế và hướng nghiên cứu tiếp
Hạn chế được ghi nhận từ phạm vi nghiên cứu
- Các kết quả đánh giá tốc độ hội tụ vẫn đòi hỏi toán tử đại diện $A_1$ phải thỏa mãn tính khả vi Fréchet và điều kiện nguồn cục bộ; chưa mở rộng được cho trường hợp toàn bộ các toán tử trong hệ đều không khả vi Fréchet.
- Không gian khảo sát trong Chương 3 giới hạn ở lớp không gian Banach phản xạ, lồi chặt có chuẩn khả vi Gâteaux đều, chưa áp dụng trực tiếp cho các không gian Banach tổng quát bất kỳ.
Hướng nghiên cứu tiếp
- Mở rộng phương pháp hiệu chỉnh cho hệ phương trình chứa các toán tử đa trị hoặc toán tử không có tính chất đơn điệu trên không gian Banach.
- Nghiên cứu các điều kiện nguồn tổng quát (chẳng hạn điều kiện nguồn dạng hàm biến phân) thay thế cho điều kiện nguồn đạo hàm Fréchet.
- Áp dụng các thuật toán hiệu chỉnh của luận án vào việc giải quyết trực tiếp các bài toán ngược quy mô lớn trong địa chấn học và xử lý ảnh y tế nhiệt học.
Giá trị tham khảo
Luận án là tài liệu chuyên khảo có giá trị học thuật cao trong chuyên ngành Toán học tính toán và Giải tích hàm ứng dụng:
- Đối với nghiên cứu sinh và học viên cao học ngành Toán: Cung cấp các kỹ thuật giải tích hàm phi tuyến, phương pháp đánh giá bất đẳng thức biến phân và kỹ thuật chứng minh sự hội tụ của các bài toán đặt không chỉnh. Đặc biệt, Định lý 2.3 và kỹ thuật xử lý phương trình hiệu chỉnh trên không gian Banach ở Chương 3 là mẫu mực phân tích lý thuyết có thể kế thừa.
- Đối với giảng viên và nhà nghiên cứu: Luận án cung cấp bức tranh tổng quan và phân loại hệ thống các phương pháp hiệu chỉnh hiện đại (Tikhonov, Landweber-Kaczmarz, Levenberg-Marquardt, Browder-Tikhonov); là tài liệu tham khảo tốt cho các bài giảng chuyên đề về giải tích số và bài toán ngược.
- Đối với kỹ sư tính toán và xử lý ảnh: Mục 1.3 của Chương 1 cung cấp mô hình toán học chi tiết về bài toán tái tạo hình ảnh chụp cắt lớp vi tính dựa trên định luật suy giảm tia X và hệ phương trình tích phân rời rạc, hỗ trợ việc thiết kế thuật toán khôi phục dữ liệu ảnh khi có nhiễu.
Câu hỏi thường gặp
1. Khái niệm bài toán đặt không chỉnh theo Hadamard được định nghĩa như thế nào trong luận án?
Theo J. Hadamard, một bài toán tìm nghiệm $A(x) = f$ ($A: X \to Y$) được gọi là đặt chỉnh nếu thỏa mãn đồng thời ba điều kiện: (1) Phương trình có nghiệm $x_0$ với mọi $f \in Y$; (2) Nghiệm $x_0$ là duy nhất; (3) Nghiệm $x_0$ phụ thuộc liên tục vào dữ liệu $f$. Nếu vi phạm ít nhất một trong ba điều kiện này, bài toán được gọi là bài toán đặt không chỉnh.
2. Sự khác biệt căn bản giữa phương pháp của luận án và phương pháp lặp Kaczmarz của A. De Cezaro (2011) là gì?
Phương pháp lặp Tikhonov-Kaczmarz của Cezaro đòi hỏi toàn bộ $N$ toán tử $A_j$ trong hệ phải thỏa mãn điều kiện khả vi Fréchet, điều kiện nón tiếp tuyến cục bộ và điều kiện nguồn. Ngược lại, phương pháp của luận án đưa về bài toán tối ưu phiếm hàm tổng hợp, chứng minh sự hội tụ và tốc độ hội tụ của nghiệm hiệu chỉnh mà chỉ cần đặt điều kiện nguồn và tính trơn lên duy nhất một toán tử $A_1$.
3. Tốc độ hội tụ của nghiệm hiệu chỉnh trong Chương 2 đạt được mức độ nào?
Trong Định lý 2.3, khi toán tử $A_1$ khả vi Fréchet có đạo hàm thỏa mãn điều kiện Lipschitz với hằng số $L$, tồn tại phần tử nguồn $\omega \in Y_1$ sao cho $x_0 - x^* = A_1'(x_0)^*\omega$ với $L|\omega|{Y_1} < 1$, khi chọn tham số hiệu chỉnh $\alpha \sim \delta^p$ ($0 < p < 2$), tốc độ hội tụ của nghiệm hiệu chỉnh $x{\alpha(\delta)}^\delta$ về nghiệm chính xác $x_0$ đạt mức: $$|x_{\alpha(\delta)}^\delta - x_0|_X = O\left(\delta^{1 - \frac{p}{2}}\right)$$
4. Bài toán thực tế nào được trình bày chi tiết trong Chương 1 để minh họa cho nguồn gốc hệ phương trình toán tử đặt không chỉnh?
Luận án trình bày chi tiết mô hình vật lý của bài toán chụp cắt lớp vi tính X-quang. Dựa trên định luật Beer về độ suy giảm cường độ tia X $dI(x) = -\mu(x)I(x)dx$, tích phân hàm hệ số hấp thụ tuyến tính $\mu(x, y)$ dọc theo đường truyền tia được xác định qua hàm hình chiếu $p(\xi, \theta) = -\ln(I/I_0)$. Rời rạc hóa tích phân này trên nhiều góc quét dẫn trực tiếp đến hệ phương trình đại số đặt không chỉnh $\sum_{k=1}^{N_s} A_{jk}\mu_j = p_{i,k}$.
5. Dạng phương trình hiệu chỉnh được luận án đề xuất cho hệ toán tử U-đơn điệu trên không gian Banach ở Chương 3 là gì?
Luận án đề xuất phương trình hiệu chỉnh dạng: $$A_1(x) + \alpha \sum_{j=2}^N \tilde{\mu} (A_j(x) - f_j^\delta) + \alpha(x - x^) = f_1^\delta$$ trong đó $X$ là không gian Banach phản xạ, lồi chặt có chuẩn khả vi Gâteaux đều, $A_j$ là các toán tử $U$-đơn điệu liên tục Lipschitz, $x^$ là phần tử ban đầu, $\tilde{\mu} \in (0, 1)$ là hằng số cố định, và $\alpha > 0$ là tham số hiệu chỉnh được xác định qua nguyên lý tựa độ lệch.
Kết luận
Luận án tiến sĩ của tác giả Nguyễn Đình Dũng đã giải quyết có hệ thống các vấn đề lý thuyết và tính toán then chốt trong việc tìm nghiệm xấp xỉ ổn định cho hệ phương trình toán tử đặt không chỉnh trên cả không gian Hilbert và không gian Banach. Bằng việc xây dựng các lược đồ hiệu chỉnh Tikhonov và Browder-Tikhonov cải tiến, công trình đã nới lỏng thành công các điều kiện nguồn ngặt nghèo trước đây, chứng minh sự hội tụ và thiết lập đánh giá tốc độ hội tụ chỉ dựa trên một toán tử thành phần. Các kết quả nghiên cứu đã được công bố trên các tạp chí chuyên ngành toán học trong nước, quốc tế và được kiểm chứng qua các thử nghiệm tính toán số.