Tổng quan luận án
Luận án tiến sĩ toán học của nghiên cứu sinh Nguyễn Văn Nhân, chuyên ngành Toán Giải Tích (mã số 1.01), được thực hiện dưới sự hướng dẫn khoa học của GS. Đặng Đình Ang và GS. Bùi Doãn Khanh tại Trường Đại học Sư phạm TP. Hồ Chí Minh năm 2000. Đề tài tập trung nghiên cứu ứng dụng phương pháp mômen và các kỹ thuật chỉnh hóa toán học để giải quyết một số bài toán ngược không chỉnh (ill-posed problems) trong giải tích ứng dụng.
Tính cấp thiết và khoảng trống nghiên cứu
Trong giải tích ứng dụng và các ngành khoa học kỹ thuật, nhiều bài toán thực tế yêu cầu xác định các đặc trưng hàm số hoặc cấu trúc hình học bên trong vật thể từ dữ liệu thực nghiệm. Tuy nhiên, dữ liệu thu thập từ các thiết bị đo đạc thực tế chỉ là một tập giá trị rời rạc tại một số vị trí hữu hạn, được gọi là các mômen. Hơn nữa, dữ liệu đo đạc luôn chứa sai số ngẫu nhiên $\epsilon > 0$. Dù sai số này rất nhỏ, sự không ổn định nội tại có thể dẫn đến sai số rất lớn ở nghiệm khôi phục, hoặc khiến bài toán trở nên vô nghiệm. Theo nghĩa Hadamard, các bài toán này là bài toán không chỉnh (không tồn tại nghiệm, nghiệm không duy nhất, hoặc nghiệm không phụ thuộc liên tục vào dữ liệu ban đầu). Do đó, việc xây dựng phương pháp tiếp cận mômen kết hợp thuật toán chỉnh hóa và đánh giá sai số nghiệm có ý nghĩa khoa học và thực tiễn cơ bản.
Mục tiêu nghiên cứu
Luận án xác định 3 mục tiêu nghiên cứu cụ thể tương ứng với 3 bài toán giải tích ứng dụng:
- Xây dựng phương pháp xấp xỉ và thuật toán chỉnh hóa dựa trên đa thức nội suy Lagrange để khôi phục hàm giải tích trong không gian Hardy $H^2(U)$ trên đĩa đơn vị từ dãy mômen đo đạc bị nhiễu, đồng thời thiết lập đánh giá sai số chỉnh hóa.
- Thiết lập mô hình toán học cho bài toán xác định hình dáng dị vật hai chiều trong lòng đất từ dữ liệu gradient trọng lực trên mặt đất; chứng minh tính duy nhất nghiệm của phương trình tích phân phi tuyến tương ứng; xấp xỉ bài toán bằng phương trình tích phân tuyến tính, đưa về bài toán mômen tuyến tính và chỉnh hóa bằng phương pháp Tikhonov.
- Đề xuất phương pháp tiếp cận mômen dựa trên tính chất giải tích và điều kiện mômen của biến đổi Radon để khôi phục trường vectơ hai chiều và ba chiều; chuyển bài toán tìm nghiệm chỉnh hóa về bài toán tìm điểm bất động của toán tử co và đánh giá sai số.
Đối tượng và phạm vi nghiên cứu
- Đối tượng nghiên cứu:
- Lớp hàm giải tích thuộc không gian Hardy $H^2(U)$ trên đĩa đơn vị phức $U = {z \in \mathbb{C} : |z| < 1}$.
- Dị vật $\Omega$ có mật độ đồng nhất trong nửa mặt phẳng $\mathbb{R}^2$ dưới lòng đất, biểu diễn qua hàm biên biên độ $z = \sigma(x)$.
- Trường vectơ vận tốc dòng chảy hoặc trường phân bố ứng suất hai chiều và ba chiều thông qua các tích phân đường/mặt của biến đổi Radon.
- Các dãy số liệu mômen rời rạc thuộc không gian dãy bị chặn $\ell^\infty$.
- Phạm vi nghiên cứu: Luận án giới hạn khảo sát trong khuôn khổ giải tích hàm, lý thuyết thế vị, phương trình tích phân và giải tích toán học thuần túy; các đánh giá sai số được thực hiện dưới dạng chặn giải tích tường minh theo mức nhiễu $\epsilon$.
Tổng quan tài liệu và vị trí của luận án
Luận án dành một phần dung lượng lớn để hệ thống hóa lịch sử phát triển của bài toán mômen và các nghiên cứu liên quan đến bài toán ngược trong địa vật lý và chụp ảnh cắt lớp.
Lịch sử bài toán mômen và các hướng nghiên cứu trước
Bài toán mômen có lịch sử phát triển gần 150 năm, khởi nguồn từ cơ học và lý thuyết giải tích:
- P. L. Tchebycheff (1873): Đặt nền móng đầu tiên khi khảo sát mức độ xác định của hàm phân bố từ dãy mômen, tìm cận trên và cận dưới đúng của tích phân thông qua công cụ phân số liên tục.
- E. Heine (1861, 1878, 1881): Phát triển lý thuyết phân số liên tục và ứng dụng đa thức trực giao vào xấp xỉ toàn phương.
- A. A. Markoff (1884, 1896): Chứng minh trọn vẹn bài toán của Tchebycheff trong luận án năm 1884 và mở rộng bài toán mômen dạng $\int_a^b x^n \varphi(x) dx = \mu_n$ với $0 \le \varphi(x) \le L$.
- T. J. Stieltjes (1894–1895): Công bố công trình khảo cứu về phân số liên tục, đưa ra công cụ tích phân Stieltjes để giải quyết bài toán tìm hàm không giảm $w(x)$ trên $[0, +\infty)$ thỏa $\int_0^{+\infty} x^n dw(x) = \mu_n$; phân định bài toán mômen xác định (nghiệm duy nhất) và bài toán mômen không xác định (vô số nghiệm).
- H. Hamburger (1920, 1921): Mở rộng bài toán Stieltjes ra toàn bộ trục thực $(-\infty, +\infty)$ và thiết lập điều kiện cần và đủ cho tính xác định của nghiệm.
- R. Nevanlinna: Sử dụng lý thuyết hàm giải tích hiện đại, đề xuất khái niệm lời giải cực trị (extremal solutions).
- M. Riesz (1921, 1922, 1923): Khảo sát bài toán mômen trên cơ sở các đa thức tựa trực giao (quasi-orthogonal polynomials).
- T. Carleman (1923, 1926) và F. Hausdorff (1923): Carleman thiết lập mối liên hệ giữa bài toán mômen, hàm tựa giải tích và dạng toàn phương vô hạn biến; Hausdorff đưa ra tiêu chuẩn cho bài toán mômen trên khoảng hữu hạn và điều kiện để nghiệm có tính chất liên tục, khả vi.
- N. I. Achyeser và M. G. Krein (1934): Mở rộng lý thuyết cho bài toán mômen lượng giác $\int_{-\pi}^\pi e^{int} dw(t) = \mu_n$.
- M. H. Stone, E. K. Haviland, H. Cramér: Mở rộng bài toán mômen sang không gian Hilbert, hệ toán tử trừu tượng và không gian nhiều chiều.
Các nghiên cứu chuyên sâu về bài toán thành phần
- Bài toán khôi phục hàm giải tích: Được khảo sát trong các công trình của J. R. Cannon, P. Colli, D. Đặng Đình Ang và các cộng sự.
- Bài toán dị thường trọng lực trong địa vật lý: R. A. Smith chứng minh tính duy nhất nghiệm; M. H. P. Bott và D. W. Oldenburg tính toán cụ thể cho trường hợp hai chiều; D. Đặng Đình Ang, R. Gorenflo và D. D. Vy chứng minh định lý duy nhất nghiệm tổng quát và chỉnh hóa bài toán ba chiều bằng xấp xỉ hữu hạn chiều.
- Bài toán biến đổi Radon và trường vectơ: S. A. Johnson, Y. Takuso, và G. W. Faris khảo sát các bài toán chụp ảnh cắt lớp (tomography) để tái tạo trường vận tốc và ứng suất vật liệu.
Khoảng trống luận án lựa chọn giải quyết
- Đối với bài toán khôi phục hàm giải tích: Luận án đề xuất cấu trúc đa thức xấp xỉ trực tiếp từ hệ phương trình nội suy Lagrange trên các điểm rời rạc, thiết lập thuật toán lựa chọn số chiều $m(\epsilon)$ phụ thuộc mức nhiễu để đạt tính ổn định và xác lập bậc hội tụ sai số.
- Đối với bài toán trọng lực: Khảo sát mô hình hai chiều với dữ kiện là gradient trọng lực thay vì dị thường trọng lực; chứng minh tính duy nhất nghiệm bằng phương pháp hàm điều hòa trong nửa mặt phẳng và biến đổi bài toán xấp xỉ tuyến tính thành bài toán mômen tương đương thông qua biến đổi Laplace.
- Đối với bài toán trường vectơ: Sử dụng các điều kiện mômen của biến đổi Radon và kỹ thuật nhiễu toán tử để đưa việc chỉnh hóa về bài toán tìm điểm bất động của toán tử co.
Cơ sở lý thuyết và phương pháp nghiên cứu
Lý thuyết, khái niệm và khung phân tích
- Khái niệm bài toán không chỉnh Hadamard: Bài toán toán học được gọi là không chỉnh nếu vi phạm ít nhất một trong ba điều kiện: (i) có nghiệm, (ii) nghiệm duy nhất, (iii) nghiệm phụ thuộc liên tục vào dữ liệu ban đầu.
- Không gian Hardy $H^2(U)$: Không gian các hàm phức $F(z) = \sum_{k=0}^\infty a_k z^k$ giải tích trên đĩa đơn vị $U = {z \in \mathbb{C} : |z| < 1}$ với chuẩn xác định bởi: $$|F|{H^2}^2 = \sum{k=0}^\infty |a_k|^2 < \infty$$
- Lý thuyết thế vị và hàm điều hòa: Biểu diễn thế trọng lực $U(x,z)$ trong nửa mặt phẳng $\mathbb{R}^2$ ($-\infty < z < H, H > 0$) thỏa mãn phương trình Laplace $\Delta U = 0$ bên ngoài vật thể khối lượng.
- Biến đổi Laplace và biến đổi Radon: Ứng dụng biến đổi tích phân để phân tích tính duy nhất của nghiệm bài toán mômen và biểu diễn hình học của trường vectơ.
- Lý thuyết điểm bất động của toán tử co: Nguyên lý ánh xạ co Banach trong không gian metric đầy đủ được sử dụng để xây dựng nghiệm chỉnh hóa cho phương trình toán tử bị nhiễu.
Phương pháp nghiên cứu
Luận án sử dụng phương pháp nghiên cứu giải tích toán học thuần túy, kết hợp các kỹ thuật suy diễn định tính và định lượng giải tích:
- Phương pháp nội suy và giải tích phức: Sử dụng đa thức Lagrange, công thức đạo hàm hàm phức, các khai triển chuỗi Taylor và đẳng thức tổ hợp để đánh giá các chặn hệ số đại số.
- Phương pháp hàm điều hòa và giải tích tiệm cận: Sử dụng tính chất triệt tiêu tại vô cực của hàm điều hòa và đạo hàm riêng theo hướng để chứng minh tính duy nhất nghiệm của phương trình tích phân phi tuyến.
- Phương pháp chỉnh hóa Tikhonov và nhiễu toán tử: Biến đổi các phương trình tích phân loại một không chỉnh thành các bài toán cực tiểu hóa phiếm hàm hoặc bài toán toán tử co ổn định.
Nguồn tư liệu và dữ liệu
Dữ liệu sử dụng trong luận án là dữ liệu toán học trừu tượng:
- Dãy dữ liệu đo đạc thực nghiệm rời rạc $\mu = (\mu_n) \in \ell^\infty$ với giả thiết độ lệch chuẩn bị chặn $|\mu - \mu^0|{\ell^\infty} = \sup{n} |\mu_n - \mu_n^0| \le \epsilon$.
- Dữ liệu hàm gradient trọng lực $f_0(x) = \frac{\partial^2 U}{\partial z^2}(x,H)$ trên bề mặt $z = H$.
- Dữ liệu tích phân đường và tích phân mặt từ phép chiếu biến đổi Radon.
Nội dung chính theo từng chương
Chương 1: Tổng quan
Chương 1 trình bày tổng quan lịch sử phát triển của bài toán mômen từ nửa cuối thế kỷ 19 đến thời điểm nghiên cứu, hệ thống hóa các kết quả của Tchebycheff, Heine, Stieltjes, Markoff, Hamburger, Nevanlinna, Riesz, Carleman, Hausdorff, Achyeser, Krein, Stone, Haviland và Cramér. Chương này phân tích cơ sở cơ học của bài toán phân bố khối lượng, sự ra đời của tích phân Stieltjes, và quá trình chuyển tiếp từ bài toán mômen giải tích cổ điển sang phương pháp mômen ứng dụng trong các bài toán ngược. Đồng thời, chương 1 giới thiệu khái quát cấu trúc 3 bài toán chính sẽ được giải quyết trong các chương 2, 3 và 4.
Chương 2: Khôi phục hàm giải tích trên đĩa đơn vị từ dãy các mômen : chỉnh hóa và sai số
Chương 2 giải quyết bài toán xác định hàm $u \in H^2(U)$ thỏa mãn: $$u(z_n) = \mu_n, \quad |z_n| < 1, \quad z_n \in U, \quad n = 1, 2, \dots$$ với $(z_n)$ là dãy vô hạn các điểm trong đĩa đơn vị $U$ có khoảng cách đến biên $\partial U$ lớn hơn 0, và $(\mu_n) \in \ell^\infty$.
1. Chứng minh tính không chỉnh theo Hadamard
Tác giả chứng minh tập hợp các dãy $(\mu_n) \in \ell^\infty$ làm cho bài toán vô nghiệm là trù mật trong $\ell^\infty$. Giả sử ứng với $\mu^0 = (\mu_n^0)$, bài toán có nghiệm $u_0 \in H^2(U)$. Với $\epsilon > 0$, xét dữ liệu nhiễu $\mu = \mu^0 + \epsilon e_1$ với $e_1 = (1, 0, 0, \dots)$. Nếu tồn tại nghiệm $u \in H^2(U)$ thỏa $u(z_n) = \mu_n$, xét $v = u - u_0 \in H^2(U)$ thì $v(z_n) = 0$ với mọi $n \ge 2$. Do $(z_n)_{n \ge 2}$ có điểm tụ trong $U$, suy ra $v \equiv 0$ trên $U$, dẫn đến $v(z_1) = \epsilon = 0$ (mâu thuẫn).
2. Xây dựng đa thức xấp xỉ và Bổ đề 2.1
Từ dữ liệu đo đạc $\mu = (\mu_n)$ với $|\mu - \mu^0|{\ell^\infty} \le \epsilon$, với mỗi $m \in \mathbb{N}$, hệ số $(a{mk}){k=0}^{m-1}$ được xác định từ hệ phương trình tuyến tính: $$\sum{k=0}^{m-1} a_{mk} z_n^k = \mu_n, \quad n = 1, \dots, m$$ Đa thức xấp xỉ $P_m(z)$ được định nghĩa bởi: $$P_m(z) = \sum_{0 \le k \le m/2} a_{mk} z^k$$
Tác giả thiết lập Bổ đề 2.1 bằng quy nạp toán học để biểu diễn và chặn các hệ số của đa thức nội suy Lagrange: $$R_p(Z_s, z) = \sum_{n=1}^s z_n^p \prod_{j \in {1,\dots,s} \setminus {n}} \frac{z - z_j}{z_n - z_j} = \sum_{k=0}^{s-1} C_{s,p,k}(Z_s) z^k$$ với các chặn tổ hợp $|C_{\beta,s,p}| \le C_{p-1}^{s-1}$.
3. Định lý 2.1 và Đánh giá sai số
Giả sử $0 < \alpha < 1$ thỏa mãn $\frac{\sqrt{\alpha}}{(1-\sqrt{\alpha})^2} < 1$, dãy $(z_n) \subset U$ phân biệt thỏa $|z_n| \le \alpha$. Đặt: $$D_m = \sup_{n \in {1,\dots,m}} \prod_{j \in {1,\dots,m} \setminus {n}} |z_n - z_j|^{-1}$$ Nếu $u_0 \in H^2(U)$ và $u_0' \in H^2(U)$, tác giả chứng minh tồn tại hàm $m(\epsilon) \to \infty$ khi $\epsilon \to 0$ sao cho sai số hội tụ về 0: $$|P_{m(\epsilon)} - u_0|{H^2}^2 \le \epsilon^2 + \frac{2(1 + m(\epsilon)/2)}{\sqrt{1-\sqrt{\alpha}}(1-\alpha)^2} \left[ 2\epsilon^2 m^2(\epsilon) D_m^2 (1+\alpha)^{2m(\epsilon)-2} + 2\frac{\alpha^{m(\epsilon)}}{(1-\sqrt{\alpha})^2 (1-\alpha)^2} \sup_n |a_n|^2 \right] + \frac{4|u_0'|{H^2}^2}{m^2(\epsilon)}$$
4. Định lý 2.2: Trường hợp riêng dãy điểm $z_n = 1/n$
Khi chọn dãy điểm $z_n = 1/n$, đại lượng $D_m$ được đánh giá tường minh: $$D_m = \sup_{n \in {1,\dots,m}} \frac{n^{m-1} m!}{(n-1)!(m-n)!} \le \frac{m^{m-1} 2^{m-1}}{(m-1)!} \le \frac{m 4^{m-1}}{\sqrt{2\pi(m-1)}}$$ Khi đó, với $m(\epsilon) = [f^{-1}(1/\epsilon)]$, hàm $f(\theta) = 4(\theta+2)\theta^2(1+\alpha)^{2\theta-2}$, tác giả thiết lập đánh giá sai số đạt tốc độ logarit: $$|P_{m(\epsilon)} - u_0|{H^2} \le \left( 1 + |u_0|{H^2} + 4|u_0'|_{H^2} \right) \left( \ln \frac{1}{\epsilon} \right)^{-1}, \quad \forall \epsilon \in (0, \epsilon_0)$$
| Đại lượng / Ký hiệu | Định nghĩa / Biểu thức giải tích | Vai trò trong đánh giá sai số |
|---|---|---|
| Không gian nghiệm $H^2(U)$ | $F(z) = \sum_{k=0}^\infty a_k z^k, \sum |a_k|^2 < \infty$ | Không gian hàm giải tích cần khôi phục |
| Mức nhiễu $\epsilon$ | $\sup_n |\mu_n - \mu_n^0| \le \epsilon$ | Biên sai số đo đạc thực nghiệm |
| Đa thức xấp xỉ $P_m(z)$ | $\sum_{0 \le k \le m/2} a_{mk} z^k$ | Nghiệm chỉnh hóa dạng đa thức cắt cụt |
| Hằng số hình học $D_m$ | $\sup_n \prod_{j \ne n} |z_n - z_j|^{-1}$ | Đại lượng đặc trưng cho sự phân bố điểm đo |
| Bậc hội tụ sai số (dãy $1/n$) | $O\left( (\ln(1/\epsilon))^{-1} \right)$ | Tốc độ hội tụ của nghiệm chỉnh hóa khi $\epsilon \to 0$ |
Chương 3: Phương trình trọng lực phi tuyến, tính duy nhất nghiệm và xấp xỉ bằng các mômen tuyến tính
Chương 3 xét bài toán ngược địa vật lý: xác định biên dị vật $z = \sigma(x)$ ($0 \le x \le 1$) nằm trong mô hình nửa mặt phẳng đất $-\infty < z < H$ ($H > 0$) từ dữ liệu gradient trọng lực trên mặt đất $z = H$. Giả định độ chênh mật độ $\rho = 1$.
1. Thiết lập phương trình tích phân phi tuyến
Gradient trọng lực $f_0(x) = \frac{\partial^2 U}{\partial z^2}(x,H)$ sinh ra trên mặt đất có biểu thức: $$f_0(x) = -\int_0^1 \left[ \frac{H - \sigma(\xi)}{(x-\xi)^2 + (H-\sigma(\xi))^2} - \frac{H}{(x-\xi)^2 + H^2} \right] d\xi$$ Bài toán quy về giải phương trình tích phân phi tuyến $f(x, \sigma) = f_0(x)$ để tìm $\sigma(x)$.
2. Định lý 3.1 về tính duy nhất nghiệm
Tác giả chứng minh phương trình trên có nhiều nhất một nghiệm liên tục $\sigma = \sigma(x)$ trên $[0,1]$ thỏa mãn điều kiện biên $\sigma(0) = \sigma(1) = 0$ và $0 < \sigma(x) \le a < H$.
- Chứng minh: Giả sử tồn tại hai nghiệm $\sigma_1$ và $\sigma_2$ tương ứng với hai miền $\Omega_1, \Omega_2$. Xét hàm thế hiệu: $$F(x,z) = \frac{\partial U_1}{\partial z} - \frac{\partial U_2}{\partial z} = \int_{\mathbb{R}^2} (\chi_{\Omega_1} - \chi_{\Omega_2}) \frac{z - \zeta}{(x-\xi)^2 + (z-\zeta)^2} d\xi d\zeta$$ Hàm $F(x,z)$ là hàm điều hòa trên $\mathbb{R}^2 \setminus (\overline{\Omega}_1 \cup \overline{\Omega}_2)$. Do dữ kiện gradient trọng lực trùng nhau trên mặt đất $z = H$, ta có $\frac{\partial F}{\partial z}(x,H) = 0$. Kết hợp tính điều hòa $\Delta F = 0$ và tính triệt tiêu tiệm cận khi $|x| \to \infty$ hoặc $z \to \infty$, tác giả suy ra $F(x,z) \equiv 0$ trên toàn miền $z \ge H$. Theo nguyên lý thác triển giải tích cho hàm điều hòa, suy ra $\Omega_1 = \Omega_2$, tức $\sigma_1(x) \equiv \sigma_2(x)$.
3. Tuyến tính hóa và đưa về bài toán mômen tuyến tính
Khi độ sâu dị vật $\sigma(\xi)$ bé so với khoảng cách $H$, khai triển tuyến tính đưa phương trình về dạng tích phân tuyến tính: $$\int_0^1 \frac{(x-\xi)^2 - H^2}{[(x-\xi)^2 + H^2]^2} \sigma(\xi) d\xi = f_0(x)$$ Thông qua phép đổi biến và tích phân tích chập, bài toán được chuyển đổi về hệ phương trình mômen tuyến tính tương đương: $$\int_0^1 \frac{\xi^n}{(x + \xi + 1)^2} \sigma(\xi) d\xi = \mu_n$$ Tính duy nhất nghiệm của hệ mômen này được chứng minh chặt chẽ nhờ tính đơn định của biến đổi tích phân Laplace. Cuối cùng, tác giả áp dụng phương pháp Tikhonov để chỉnh hóa phương trình tuyến tính và thiết lập các đánh giá sai số nghiệm.
Chương 4: Khôi phục trường vectơ từ biến đổi Radon, tiếp cận mômen và chỉnh hóa
Chương 4 nghiên cứu bài toán chụp ảnh cắt lớp (tomography) nhằm khôi phục trường vectơ hai chiều và ba chiều từ các số liệu đo tích phân đường hoặc tích phân mặt (biến đổi Radon). Bài toán này xuất phát từ nhu cầu xác định trường vận tốc dòng chảy chất lỏng hoặc phân bố tensor ứng suất dư trong lòng vật liệu kim loại.
1. Mô hình hóa bài toán biến đổi Radon
Trường vectơ cần khôi phục được biểu diễn qua ánh xạ tích phân Radon. Dữ liệu thực nghiệm là các phép chiếu tích phân theo các hướng đơn vị $\mathbf{U}_n$ ($|\mathbf{U}_n| = 1$ trong $\mathbb{R}^k$). Sử dụng tính chất của hàm giải tích và điều kiện mômen của biến đổi Radon, tác giả chuyển bài toán khôi phục hàm $f: \mathbb{R}^k \to \mathbb{R}$ về hệ phương trình mômen toán tử.
2. Kỹ thuật chỉnh hóa và nguyên lý điểm bất động
Bài toán xác định hàm $f$ từ dãy mômen biến đổi Radon là bài toán ngược không chỉnh nghiêm trọng. Tác giả đề xuất kỹ thuật làm nhiễu phương trình toán tử:
- Phương trình toán tử ban đầu $A f = \mu$ được thay thế bởi phương trình hiệu chỉnh có tham số $A_\alpha f_\alpha = \mu^\epsilon$.
- Biến đổi phương trình hiệu chỉnh về dạng phương trình điểm bất động $f_\alpha = T_\alpha(f_\alpha)$, trong đó $T_\alpha$ được chứng minh là một toán tử co trên không gian hàm thích hợp khi tham số chỉnh hóa $\alpha$ được chọn tương thích với sai số $\epsilon$.
- Áp dụng nguyên lý ánh xạ co Banach, nghiệm chỉnh hóa duy nhất được xây dựng thông qua dãy lặp hội tụ.
- Luận án thiết lập các bất đẳng thức đánh giá sai số giữa nghiệm chỉnh hóa và nghiệm chính xác, sau đó mở rộng thuật toán cho trường vectơ 2D và 3D.
Kết quả và những đóng góp mới
Luận án đã đạt được các kết quả giải tích và đóng góp học thuật mới sau:
Đóng góp về mặt lý luận
- Thiết lập thuật toán chỉnh hóa giải tích trên đĩa đơn vị: Xây dựng tường minh đa thức xấp xỉ $P_m(z)$ cho hàm thuộc không gian Hardy $H^2(U)$ dựa trên hệ phương trình đại số nội suy Lagrange từ dữ liệu mômen rời rạc có nhiễu.
- Chứng minh tính giải được và đánh giá sai số tường minh: Chứng minh tính trù mật của tập dữ liệu vô nghiệm trong $\ell^\infty$ (khẳng định tính không chỉnh Hadamard); thiết lập công thức đánh giá sai số tổng quát trong Định lý 2.1 và xác lập tốc độ hội tụ cấp logarit $(\ln(1/\epsilon))^{-1}$ cho trường hợp dãy điểm $z_n = 1/n$ trong Định lý 2.2.
- Chứng minh tính duy nhất nghiệm cho bài toán gradient trọng lực 2D: Thiết lập mô hình toán học sử dụng gradient trọng lực $\frac{\partial^2 U}{\partial z^2}(x,H)$ và chứng minh định lý duy nhất nghiệm (Định lý 3.1) cho phương trình tích phân phi tuyến xác định biên dị vật $\sigma(x)$ bằng phương pháp hàm điều hòa trong nửa mặt phẳng.
- Liên kết bài toán trọng lực tuyến tính với bài toán mômen Laplace: Đưa phương trình trọng lực xấp xỉ tuyến tính về bài toán mômen tuyến tính tương đương và chứng minh tính duy nhất nghiệm qua biến đổi Laplace; áp dụng thành công phương pháp chỉnh hóa Tikhonov.
- Đề xuất phương pháp điểm bất động cho bài toán biến đổi Radon: Vận dụng điều kiện mômen của biến đổi Radon kết hợp kỹ thuật nhiễu toán tử để quy bài toán khôi phục trường vectơ không chỉnh về bài toán tìm điểm bất động của toán tử co, mở rộng kết quả cho trường vectơ 2 chiều và 3 chiều.
Đóng góp về thực tiễn và giải pháp đề xuất
- Đề xuất nguyên lý lựa chọn tham số cắt cụt số chiều $m(\epsilon)$ và tham số chỉnh hóa $\alpha(\epsilon)$ phụ thuộc trực tiếp vào biên độ nhiễu đo đạc $\epsilon$, cung cấp quy trình toán học để đảm bảo tính ổn định của nghiệm.
- Cung cấp cơ sở giải tích cho việc xử lý số liệu đo đạc gradient trọng lực trong thăm dò địa vật lý và xử lý dữ liệu chiếu cắt lớp trong đo đạc vận tốc dòng chảy, ứng suất cơ học.
Hạn chế và hướng nghiên cứu tiếp
Hạn chế
- Các kết quả về trường trọng lực trong Chương 3 chủ yếu được thiết lập cho mô hình hai chiều phẳng với giả thiết chênh lệch mật độ là hằng số ($\rho = 1$) và biên dị vật $\sigma(x)$ thỏa mãn các điều kiện biên trơn, biên độ bị chặn ($0 < \sigma(x) \le a < H$).
- Bậc hội tụ sai số của phương pháp nội suy Lagrange trong đĩa đơn vị với cách chọn điểm $z_n = 1/n$ đạt tốc độ logarit $(\ln(1/\epsilon))^{-1}$, là tốc độ hội tụ tương đối chậm khi mức nhiễu $\epsilon$ giảm.
- Các đánh giá sai số và nghiệm chỉnh hóa dừng lại ở mức độ giải tích toán học thuần túy, chưa đi kèm các bảng thử nghiệm tính toán số mô phỏng trên máy tính hay số liệu thực nghiệm hiện trường.
Hướng nghiên cứu tiếp
- Mở rộng bài toán gradient trọng lực phi tuyến cho mô hình không gian ba chiều toàn vẹn với hàm mật độ biến thiên $\rho(x,y,z)$.
- Khảo sát các quy luật phân bố điểm đo $(z_n)$ khác trên đĩa đơn vị $U$ nhằm cải thiện tốc độ hội tụ của nghiệm chỉnh hóa vượt trên cấp logarit.
- Xây dựng thuật toán số hóa và chương trình tính toán thực nghiệm cho bài toán khôi phục trường vectơ từ biến đổi Radon trong không gian ba chiều.
Giá trị tham khảo
Luận án là tài liệu tham khảo chuyên khảo 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: Tham khảo phương pháp tiếp cận mômen, kỹ thuật đánh giá hệ số đa thức Lagrange nhiều biến qua Bổ đề 2.1, và phương pháp sử dụng nguyên lý điểm bất động của toán tử co để giải bài toán ngược không chỉnh.
- Nhà nghiên cứu trong lĩnh vực Địa vật lý và Vật lý Toán: Tham khảo mô hình thế trọng lực, kỹ thuật sử dụng hàm điều hòa để chứng minh tính duy nhất nghiệm của phương trình tích phân phi tuyến và phương pháp tuyến tính hóa bài toán trọng lực.
- Chuyên gia xử lý ảnh và Kỹ thuật đo lường (Tomography): Tham khảo cách tiếp cận biến đổi Radon dựa trên điều kiện mômen để tái tạo trường vectơ hai chiều và ba chiều.
Câu hỏi thường gặp
1. Tại sao bài toán khôi phục hàm giải tích từ dãy mômen trên đĩa đơn vị lại không chỉnh theo Hadamard?
Bài toán đặt ra là tìm $u \in H^2(U)$ thỏa mãn $u(z_n) = \mu_n$. Tác giả chứng minh rằng tập hợp các dãy dữ liệu $(\mu_n) \in \ell^\infty$ làm cho bài toán vô nghiệm là trù mật trong $\ell^\infty$. Do đó, với một nghiệm $u_0$ ứng với dữ liệu chính xác $\mu^0$, chỉ cần một nhiễu vô cùng bé $\epsilon e_1$ cũng khiến bài toán mới hoàn toàn không tồn tại nghiệm trong $H^2(U)$, vi phạm điều kiện tồn tại nghiệm và tính phụ thuộc liên tục của nghiệm vào dữ liệu đo.
2. Tác giả đã giải quyết tính không chỉnh của bài toán trong Chương 2 bằng cách nào?
Tác giả thay thế việc tìm hàm giải tích vô hạn chiều bằng việc tìm đa thức xấp xỉ cắt cụt $P_m(z) = \sum_{0 \le k \le m/2} a_{mk} z^k$, trong đó hệ số $(a_{mk})$ là nghiệm của hệ phương trình nội suy Lagrange cấp $m$. Bằng cách chọn bậc đa thức $m = m(\epsilon)$ phụ thuộc nghịch đảo vào sai số đo đạc $\epsilon$ (sao cho $m(\epsilon) \to \infty$ khi $\epsilon \to 0$), nghiệm xấp xỉ $P_{m(\epsilon)}$ được chứng minh là hội tụ về nghiệm chính xác $u_0$ theo chuẩn của không gian Hardy $H^2(U)$.
3. Ưu điểm của việc sử dụng dữ kiện gradient trọng lực thay vì dị thường trọng lực trong Chương 3 là gì?
Dữ kiện gradient trọng lực $f_0(x) = \frac{\partial^2 U}{\partial z^2}(x,H)$ là đạo hàm cấp hai của thế trọng lực, phản ánh trực tiếp và nhạy bén hơn sự thay đổi hình học của biên dị vật $\sigma(x)$ so với dị thường trọng lực (đạo hàm cấp một). Về mặt toán học, biểu thức tích phân của gradient trọng lực cho phép áp dụng hiệu quả tính chất của hàm điều hòa $F(x,z)$ trong nửa mặt phẳng và điều kiện triệt tiêu tại vô cực để thiết lập chứng minh tính duy nhất nghiệm cho bài toán phi tuyến hai chiều.
4. Bổ đề 2.1 trong luận án đóng vai trò gì trong việc đánh giá sai số chỉnh hóa?
Bổ đề 2.1 thiết lập công thức tường minh và chặn trên cho các hệ số $C_{s,p,k}(Z_s)$ của đa thức nội suy Lagrange $R_p(Z_s, z)$ qua các bất đẳng thức tổ hợp $|C_{\beta,s,p}| \le C_{p-1}^{s-1}$. Bổ đề này là công cụ nền tảng để phân tách và kiểm soát sai số giữa hệ số đa thức xấp xỉ $a_{mk}$ và hệ số khai triển Taylor thực tế $a_k$ của hàm cần tìm, từ đó dẫn đến việc chứng minh Định lý 2.1 và Định lý 2.2.
5. Phương pháp toán tử co được ứng dụng như thế nào trong bài toán biến đổi Radon ở Chương 4?
Đối với bài toán khôi phục trường vectơ từ biến đổi Radon, tác giả đưa bài toán về phương trình toán tử không chỉnh dạng $A f = \mu$. Sau đó, phương trình được làm nhiễu có kiểm soát thành $A_\alpha f_\alpha = \mu^\epsilon$ và biến đổi tương đương thành phương trình dạng $f_\alpha = T_\alpha(f_\alpha)$. Tác giả chứng minh rằng với việc lựa chọn tham số $\alpha$ phù hợp, toán tử $T_\alpha$ thỏa mãn điều kiện co trong không gian metric hàm, đảm bảo sự tồn tại duy nhất của nghiệm chỉnh hóa và sự hội tụ của thuật toán lặp.
Kết luận
Luận án tiến sĩ của Nguyễn Văn Nhân đã nghiên cứu một cách có hệ thống việc áp dụng phương pháp mômen và các kỹ thuật chỉnh hóa giải tích để giải quyết ba bài toán ngược không chỉnh quan trọng: khôi phục hàm trong không gian Hardy $H^2(U)$, xác định hình dáng dị vật từ gradient trọng lực phi tuyến, và khôi phục trường vectơ từ biến đổi Radon. Công trình đã đóng góp các kết quả toán học chuẩn xác về chứng minh tính duy nhất nghiệm bằng hàm điều hòa và biến đổi Laplace, đồng thời thiết lập các đánh giá sai số chỉnh hóa tường minh thông qua đa thức nội suy Lagrange và nguyên lý ánh xạ co. Các kết quả này cung cấp cơ sở lý thuyết giải tích quan trọng cho các bài toán ngược trong vật lý toán, địa vật lý thăm dò và kỹ thuật chụp ảnh cắt lớp.