Luận án tiến sĩ: Một số bài toán không chỉnh với đạo hàm bậc không nguyên - Trần Thị Khiếu (ĐH KHTN TP.HCM)

Lý thuyết vi tích phân phân số cung cấp công cụ toán học mạnh mẽ để mô hình hóa các quá trình vật lý có tính nhớ. Đạo hàm Caputo là một trong những định nghĩa toán tử phân số phổ biến nhất hiện nay. Điểm nổi bật của đạo hàm Caputo nằm ở điều kiện ban đầu. Các điều kiện biên và giá trị ban đầu trong phương trình chứa đạo hàm Caputo giữ nguyên dạng vi phân cấp nguyên cổ điển. Đặc tính này mang ý nghĩa vật lý rõ ràng và tạo thuận lợi lớn cho việc thu thập dữ liệu thực nghiệm. Toán tử Caputo tác động lên hàm số thông qua việc lấy đạo hàm nguyên trước rồi tích phân phân số cấp tương ứng. Cấu trúc toán tử này giúp mô tả chính xác quá trình suy giảm năng lượng và lan truyền sóng trong môi trường không đồng nhất. Việc phân tích tính chất giải tích của đạo hàm Caputo tạo tiền đề phân tích sự hội tụ của nghiệm.

Bên cạnh toán tử Caputo, đạo hàm Riemann-Liouville giữ vị trí trung tâm trong giải tích toán học bậc phân số. Đạo hàm Riemann-Liouville thực hiện phép tích phân phân số trước khi áp dụng toán tử vi phân cấp nguyên. Trong lý thuyết phương trình vi phân bậc phân số, hàm Mittag-Leffler đóng vai trò tương tự như hàm mũ trong giải tích cổ điển. Hàm Mittag-Leffler một tham số và hai tham số xuất hiện tự nhiên trong nghiệm giải tích của các phương trình vi phân phân số tuyến tính. Hàm số này thể hiện đặc tính phân rã chậm theo quy luật lũy thừa. Tính chất tiệm cận và chặn trên của hàm Mittag-Leffler là công cụ then chốt để đánh giá năng lượng và ước lượng sai số của nghiệm xấp xỉ. Việc kết hợp đạo hàm Riemann-Liouville với hàm Mittag-Leffler mang lại nền tảng giải tích vững chắc cho các bài toán truyền nhiệt và khuếch tán.

Theo quan niệm của Jacques Hadamard, một bài toán toán lý được coi là đặt chỉnh nếu thỏa mãn ba điều kiện: nghiệm tồn tại, nghiệm duy nhất và nghiệm phụ thuộc liên tục vào dữ liệu đầu vào. Khi một trong ba điều kiện này bị vi phạm, bài toán trở thành bài toán không chỉnh. Đa số các mô hình ngược thời gian hoặc xác định nguồn nhiệt đều rơi vào trạng thái không chỉnh theo nghĩa Hadamard. Nhiễu nhỏ trong dữ liệu đo đạc thực tế có thể dẫn đến sai số vô cùng lớn ở nghiệm tìm kiếm. Để khắc phục hiện tượng mất ổn định này, các nhà toán học phát triển lý thuyết chỉnh hóa. Mục tiêu chính là xây dựng một họ toán tử xấp xỉ liên tục phụ thuộc vào tham số chỉnh hóa. Giải pháp chỉnh hóa giúp khôi phục tính ổn định nghiệm và kiểm soát sai số trong phạm vi cho phép.

Bản chất không ổn định của bài toán ngược thời gian bắt nguồn từ đặc tính phân rã phổ của toán tử truyền tải nhiệt. Trong phương trình vi phân phân số, các tần số cao trong dữ liệu ban đầu bị triệt tiêu rất nhanh theo hàm phân rã bậc phân số. Khi tiến hành giải ngược từ thời điểm cuối về thời điểm ban đầu, toán tử ngược tương ứng sẽ khuếch đại vô hạn các thành phần tần số cao. Một lượng nhiễu đo đạc rất nhỏ cũng khiến nghiệm toán học bùng nổ và hoàn toàn mất ý nghĩa vật lý. Do đó, bài toán không thỏa mãn điều kiện ổn định của tính đặt chỉnh Hadamard. Các ví dụ giải tích cụ thể chứng minh rằng chuỗi Fourier biểu diễn nghiệm hình thức phân kỳ nếu dữ liệu đầu vào chứa sai số dù là vô cùng bé.

Để ổn định hóa nghiệm của bài toán tuyến tính, phương pháp hàm lọc bậc không nguyên được thiết lập và ứng dụng linh hoạt. Toán tử dạng hỗn hợp trong mô hình khuếch tán bao gồm sự kết hợp giữa toán tử Laplace cổ điển và toán tử vi phân cấp phân số. Kỹ thuật hàm lọc hoạt động bằng cách triệt tiêu hoặc làm suy giảm có kiểm soát các thành phần tần số cao gây mất ổn định trong khai triển phổ. Hàm lọc được thiết kế dựa trên các hàm trơn đặc biệt và phụ thuộc chặt chẽ vào tham số chỉnh hóa. Khi tham số chỉnh hóa tiến về không cùng với mức độ nhiễu suy giảm, nghiệm chỉnh hóa sẽ tiệm cận về nghiệm thực của bài toán ban đầu. Phương pháp này đảm bảo tính khả thi trong tính toán số và bảo toàn cấu trúc nghiệm vật lý.

Đánh giá sai số và tốc độ hội tụ là tiêu chuẩn học thuật quan trọng để xác nhận hiệu quả của phương pháp chỉnh hóa. Dưới các điều kiện tiên nghiệm khác nhau về độ trơn của nghiệm chính xác, luận án thiết lập các chặn sai số tối ưu giữa nghiệm chỉnh hóa và nghiệm thực. Với điều kiện nguồn yếu trong không gian Hilbert, tốc độ hội tụ đạt được dạng Logarithm. Đây là tốc độ hội tụ chậm nhưng phản ánh đúng bản chất không chỉnh nghiêm trọng của bài toán ngược thời gian bậc phân số. Khi nghiệm chính xác thỏa mãn điều kiện trơn cao hơn theo thang không gian Sobolev, phương pháp chỉnh hóa mang lại tốc độ hội tụ tối ưu dạng Holder. Tốc độ dạng lũy thừa này khẳng định tính ổn định vững chắc và độ tin cậy của thuật toán xấp xỉ được đề xuất.

Phương trình khuếch tán bậc phân số đóng vai trò trung tâm trong mô tả hiện tượng siêu khuếch tán và dưới khuếch tán. Toán tử Laplace bậc không nguyên biểu diễn các tương tác không gian phi cục bộ tầm xa của các hạt chuyển động. Sự phối hợp giữa đạo hàm thời gian cấp phân số và toán tử vi phân không gian phân số tạo nên mô hình toán học tổng quát. Cấu trúc phổ của toán tử này có các giá trị riêng phân bố dày đặc, làm gia tăng tính nhạy cảm của hệ đối với nhiễu dữ liệu. Bài toán ngược thời gian liên kết với phương trình khuếch tán phân số này đòi hỏi kỹ thuật phân tích giải tích nâng cao. Việc hiểu rõ tính chất triệt tiêu của các toán tử nửa nhóm sinh bởi toán tử Laplace phân số là chìa khóa để giải quyết bài toán.

Điều kiện Lipschitz toàn cục giả định rằng hàm nguồn phi tuyến có độ biến thiên bị chặn bởi một hằng số xác định trên toàn bộ không gian hàm. Đây là giả thiết toán học quan trọng giúp kiểm soát sự phi tuyến của mô hình. Bài toán ngược phi tuyến được chuyển hóa thành phương trình tích phân phi tuyến tương đương với toán tử tích phân chứa hàm Mittag-Leffler. Kỹ thuật chỉnh hóa sử dụng toán tử lọc thích hợp để thay thế các hạt nhân kỳ dị và bất ổn định. Tham số chỉnh hóa được lựa chọn tự động theo nguyên lý bù sai số tiên nghiệm hoặc nguyên lý Morley. Phương pháp này biến đổi bài toán không chỉnh ban đầu thành một phương trình toán tử chỉnh đặt, đảm bảo sự tồn tại duy nhất của nghiệm xấp xỉ trong không gian Banach.

Quy trình xây dựng nghiệm chỉnh hóa dựa trên việc áp dụng nguyên lý ánh xạ co cho toán tử tích phân chỉnh hóa. Khi tham số chỉnh hóa được cố định, ánh xạ tích phân trở thành một phép co nghiêm ngặt trong không gian các hàm liên tục theo thời gian. Sự tồn tại và duy nhất của nghiệm chỉnh hóa được chứng minh chặt chẽ thông qua dãy xấp xỉ liên tiếp Picard. Luận án thiết lập các ước lượng sai số chi tiết giữa nghiệm xấp xỉ nhận được từ dữ liệu nhiễu và nghiệm chính xác của bài toán gốc. Tốc độ hội tụ dạng Holder được xác lập rõ ràng khi mức độ nhiễu tiến về không. Kết quả phân tích khẳng định phương pháp hàm lọc duy trì sự ổn định bền vững ngay cả dưới tác động của các thành phần phi tuyến phức tạp.

Tính phi tuyến Lipschitz địa phương phá vỡ tính co toàn cục của các toán tử tích phân thông thường. Khi giải bài toán ngược thời gian, năng lượng của hệ có xu hướng tăng nhanh khi đi ngược chiều thời gian, dễ khiến nghiệm vượt ra ngoài miền kiểm soát Lipschitz. Do đó, các kỹ thuật ánh xạ co cổ điển không thể áp dụng trực tiếp trên toàn bộ miền thời gian. Ngoài ra, việc dữ liệu quan sát bị nhiễu có thể đẩy nghiệm xấp xỉ ra xa miền ổn định địa phương, gây ra sự phân kỳ nghiêm trọng trong quá trình lặp. Sự kết hợp giữa tính không chỉnh cố hữu của bài toán ngược và tính không ổn định của phi tuyến địa phương đòi hỏi phải có kỹ thuật giới hạn biên độ nghiệm một cách chủ động và kiểm soát được miền nghiệm.

Để vượt qua giới hạn của tính chất địa phương, kỹ thuật xấp xỉ bằng dãy hàm cắt cụt được phát triển. Hàm nguồn Lipschitz địa phương được xấp xỉ bởi một dãy các hàm nhân tạo thỏa mãn điều kiện Lipschitz toàn cục trên toàn không gian. Mỗi hàm trong dãy xấp xỉ trùng khớp với hàm ban đầu bên trong một hình cầu bán kính xác định và bị chặn trơn tru ở bên ngoài. Thông qua phép cắt cụt này, bài toán cục bộ được chuyển đổi thành một họ bài toán toàn cục có thể giải được bằng phương pháp hàm lọc. Bán kính cắt cụt được lựa chọn cẩn thận phụ thuộc vào mức độ trơn của nghiệm ban đầu và mức độ nhiễu. Kỹ thuật này giúp kiểm soát chặt chẽ biên độ của nghiệm trong suốt quá trình chỉnh hóa.

Nghiệm chỉnh hóa của bài toán phi tuyến Lipschitz địa phương được xác định như nghiệm của phương trình tích phân đã qua xử lý cắt cụt và làm trơn phổ. Bằng cách thiết lập mối quan hệ đồng thời giữa tham số chỉnh hóa, tham số cắt cụt và mức độ nhiễu dữ liệu, tính ổn định nghiệm được khôi phục hoàn toàn. Luận án chứng minh rằng nghiệm chỉnh hóa hội tụ mạnh về nghiệm chính xác trong chuẩn không gian thích hợp khi mức nhiễu tiến dần về không. Đánh giá sai số tiệm cận đạt được ở dạng Logarithm hoặc Holder tùy thuộc vào không gian tiên nghiệm của nghiệm thực. Kết quả này đánh dấu bước tiến quan trọng trong lý thuyết giải tích số cho các phương trình đạo hàm bậc không nguyên với cấu trúc phi tuyến phức tạp.

Thuật toán tính toán số cho bài toán tuyến tính sử dụng phương pháp sai phân hữu hạn kết hợp với khai triển hàm trực giao. Không gian bài toán được rời rạc hóa thành lưới điểm đều, trong khi đạo hàm thời gian bậc phân số được xấp xỉ qua công thức sai phân vi phân Caputo. Toán tử Laplace dạng hỗn hợp được rời rạc hóa bằng ma trận vi phân đối xứng. Để kiểm tra tính ổn định, dữ liệu cuối thời gian được cộng thêm nhiễu trắng nhân tạo với độ lớn dao động từ 1% đến 10%. Thuật toán áp dụng phương pháp hàm lọc để lọc bỏ nhiễu và khôi phục phân bố trạng thái ban đầu. Các bảng số liệu và đồ thị sai số khẳng định nghiệm số xấp xỉ bám sát nghiệm giải tích chính xác, loại bỏ hoàn toàn hiện tượng dao động bùng nổ của bài toán không chỉnh.

Đối với bài toán phi tuyến có hàm nguồn, thuật toán lặp Picard kết hợp kỹ thuật cắt cụt được cài đặt để giải phương trình tích phân phi tuyến rời rạc. Tại mỗi bước lặp, phương pháp hàm lọc bậc không nguyên được áp dụng để giữ cho chuỗi nghiệm luôn nằm trong miền ổn định. Quá trình tính toán số được tiến hành trên hai lớp mô hình: mô hình có toán tử Laplace bậc không nguyên và mô hình toán tử hỗn hợp. Các ví dụ số với hàm nguồn dạng đa thức và hàm mũ minh họa khả năng khôi phục biên dạng nhiệt ban đầu rất mượt mà. Đồ thị so sánh giữa nghiệm giải tích, nghiệm thực nghiệm khi có nhiễu và nghiệm chỉnh hóa cho thấy sai số suy giảm đều đặn theo đúng tốc độ hội tụ lý thuyết đã dự báo.

Luận án thực hiện so sánh đối chuẩn giữa phương pháp chỉnh hóa Tikhonov cải biên và phương pháp hàm lọc bậc không nguyên mới đề xuất. Phương pháp Tikhonov bổ sung số hạng phạt vào hàm mục tiêu bình phương tối thiểu để kiểm soát chuẩn của nghiệm. Trong khi đó, phương pháp hàm lọc can thiệp trực tiếp vào phổ giá trị riêng của toán tử vi phân phân số. Kết quả thực nghiệm số cho thấy phương pháp hàm lọc đạt độ chính xác cao hơn tại các vùng biên và có thời gian thực thi thuật toán ngắn hơn. Sự kết hợp linh hoạt giữa hai phương pháp này mang lại giải pháp toàn diện cho việc giải quyết các bài toán ngược thời gian trong phương trình khuếch tán bậc phân số, mở ra hướng ứng dụng rộng rãi trong kỹ thuật và vật lý hiện đại.