Tổng quan luận án
Luận án tiến sĩ toán học với đề tài "Một số vấn đề về phương trình vi phân phân thứ Caputo ngẫu nhiên" do nghiên cứu sinh Phan Thị Hương thực hiện tại Đại học Kỹ thuật Lê Quý Đôn, thuộc chuyên ngành Toán ứng dụng (mã số: 9 46 01 12), dưới sự hướng dẫn khoa học của PGS. Đoàn Thái Sơn và TS. Tạ Ngọc Ánh (bảo vệ năm 2020).
Tính cấp thiết và khoảng trống nghiên cứu
Phép tính vi phân, tích phân và phương trình vi phân thường là công cụ truyền thống để mô tả các quá trình tiến hóa. Tuy nhiên, nhiều hiện tượng tự nhiên và xã hội có đặc tính phụ thuộc vào toàn bộ lịch sử quá khứ (hiệu ứng trí nhớ). Giải tích phân thứ cho phép mô hình hóa các quá trình có trí nhớ hiệu quả hơn phép tính cổ điển. Trong các dạng đạo hàm phân thứ, đạo hàm Caputo (do M. Caputo đề xuất năm 1969) có ưu thế thực tiễn do các bài toán giá trị ban đầu gắn với đạo hàm này sử dụng các điều kiện ban đầu chứa đạo hàm bậc nguyên, mang ý nghĩa vật lý rõ ràng so với đạo hàm Riemann-Liouville.
Lý thuyết phương trình vi phân phân thứ Caputo ngẫu nhiên được hình thành từ sự kết hợp giữa giải tích phân thứ Caputo và giải tích ngẫu nhiên, giúp mô tả các hệ thống vừa có tính nhớ vừa chịu tác động của nhiễu ngẫu nhiên. Mặc dù phương trình vi phân phân thứ tất định đã được nghiên cứu sâu rộng, các công trình về phương trình vi phân phân thứ Caputo ngẫu nhiên còn tương đối ít và tồn tại nhiều khoảng trống:
- Nghiên cứu về nghiệm cổ điển (classical solutions) mới chỉ được thực hiện trên một số khoảng tham số hẹp hoặc gặp sai sót trong việc thác triển nghiệm địa phương ra toàn cục do đặc tính phụ thuộc quá khứ không cho phép dịch chuyển thời gian như phương trình vi phân thường.
- Nghiên cứu về nghiệm nhẹ (mild solutions) đòi hỏi các điều kiện hệ số tương đối chặt.
- Chưa có công thức biến thiên hằng số cho phương trình vi phân phân thứ Caputo ngẫu nhiên để liên kết nghiệm của phương trình phi tuyến/không thuần nhất với phương trình tuyến tính thuần nhất tương ứng.
- Chưa có các đánh giá về sự phân tách tiệm cận giữa các nghiệm và tính chất của số mũ Lyapunov bình phương trung bình.
- Về giải số, các công bố trước đó về phương trình Volterra ngẫu nhiên với nhân kỳ dị chưa đưa ra được tốc độ hội tụ hiển của lược đồ rời rạc hóa.
Mục tiêu nghiên cứu
Mục tiêu tổng quát của luận án được chia thành hai chủ điểm chính:
- (i) Nghiên cứu một số tính chất định tính của nghiệm phương trình vi phân phân thứ Caputo ngẫu nhiên (sự tồn tại, duy nhất, tính phụ thuộc liên tục vào điều kiện ban đầu, công thức biến thiên hằng số, sự phân tách tiệm cận và số mũ Lyapunov).
- (ii) Xây dựng phương pháp giải số và đánh giá sai số, tốc độ hội tụ của nghiệm số đối với phương trình vi phân phân thứ Caputo ngẫu nhiên.
Đối tượng và phạm vi nghiên cứu
Luận án xác định các nội dung cụ thể trong đối tượng và phạm vi nghiên cứu bao gồm:
- Nội dung 1: Sự tồn tại và duy nhất nghiệm cổ điển và nghiệm nhẹ của phương trình vi phân phân thứ Caputo ngẫu nhiên.
- Nội dung 2: Công thức biến thiên hằng số cho phương trình vi phân phân thứ Caputo ngẫu nhiên.
- Nội dung 3: Một số tính chất định tính của nghiệm (sự phụ thuộc liên tục vào điều kiện ban đầu, cận dưới cho sự phân tách tiệm cận giữa hai nghiệm phân biệt, tính không âm của số mũ Lyapunov bình phương trung bình).
- Nội dung 4: Xây dựng lược đồ số kiểu Euler-Maruyama và lược đồ Euler-Maruyama mũ cho phương trình vi phân phân thứ Caputo ngẫu nhiên, đánh giá tốc độ hội tụ và tính ổn định.
- Phạm vi toán học: Phương trình vi phân phân thứ Caputo ngẫu nhiên bậc $\alpha \in \left(\frac{1}{2}, 1\right)$ trên không gian Euclide $\mathbb{R}^d$, xét trên đoạn thời gian hữu hạn $[0, T]$ và khảo sát dáng điệu tiệm cận khi $t \to \infty$. Nhiễu ngẫu nhiên được điều khiển bởi chuyển động Brown một chiều tiêu chuẩn $(W_t)_{t \ge 0}$ trên không gian xác suất đầy đủ $(\Omega, \mathcal{F}, \mathbb{P})$ trang bị bộ lọc làm rộng $(\mathcal{F}t){t \ge 0}$.
Tổng quan tài liệu và vị trí của luận án
Các hướng nghiên cứu trước được điểm lại
Trong phần mở đầu và tổng quan, tác giả điểm lại các mốc phát triển chính của giải tích phân thứ và phương trình vi phân ngẫu nhiên:
- Khái niệm tích phân và đạo hàm phân thứ Riemann-Liouville được phát triển bởi Abel, Riemann và Liouville trong nửa đầu thế kỷ 19, nhưng gặp hạn chế do điều kiện ban đầu không có ý nghĩa vật lý.
- M. Caputo (1969) đề xuất đạo hàm phân thứ Caputo nhằm giải quyết bài toán nhớt, cho phép sử dụng điều kiện ban đầu có ý nghĩa vật lý trực tiếp.
- K. Spenier cùng nhiều chuyên khảo quốc tế (như bộ sách 8 cuốn năm 2019 về giải tích phân thứ, giải số và ứng dụng) đã hệ thống hóa lý thuyết giải tích phân thứ tất định.
- Về phương trình vi phân ngẫu nhiên bậc nguyên: Lý thuyết giải tích ngẫu nhiên Itô, định lý biểu diễn Itô, định lý tồn tại duy nhất nghiệm theo điều kiện Lipschitz và tăng trưởng không quá tuyến tính, cùng phương pháp xấp xỉ Euler-Maruyama đã được chuẩn hóa.
- Về phương trình vi phân phân thứ Caputo ngẫu nhiên:
- Công trình của Wang nghiên cứu phương trình tích phân ngẫu nhiên nhưng chưa khảo sát trường hợp tham số $\alpha_2 \in \left(\frac{1}{4}, \frac{1}{2}\right)$, tương ứng với bậc phân thứ $\alpha \in \left(\frac{1}{2}, \frac{3}{4}\right)$.
- Công trình của Y. Wang và các cộng sự (2016) nghiên cứu phương trình vi phân phân thứ Caputo ngẫu nhiên trên không gian Banach, đã chứng minh được sự tồn tại và duy nhất nghiệm địa phương trên khoảng $[0, T_a]$, nhưng việc dùng phương pháp xấp xỉ liên tiếp để kéo dài nghiệm ra $[0, \infty)$ bị lỗi về mặt toán học. Do tính nhớ của đạo hàm phân thứ, việc tịnh tiến thời gian từ $[T^, T^+\delta)$ về $[0, T^*+\delta)$ không bảo toàn phương trình.
- Về giải số: Nghiên cứu của tác giả về phương trình Volterra ngẫu nhiên với nhân kỳ dị đã thiết lập lược đồ kiểu Euler nhưng chưa đưa ra được tốc độ hội tụ hiển.
Khoảng trống luận án lựa chọn giải quyết
Luận án tập trung giải quyết trực tiếp các hạn chế trên:
- Thiết lập chứng minh chặt chẽ về sự tồn tại và duy nhất nghiệm cổ điển toàn cục trên $[0, T]$ cho mọi $\alpha \in \left(\frac{1}{2}, 1\right)$ bằng cách xây dựng chuẩn có trọng số thích hợp thay vì dùng phương pháp kéo dài nghiệm địa phương.
- Giảm nhẹ các điều kiện biên và điều kiện hệ số để chứng minh sự tồn tại và duy nhất nghiệm nhẹ.
- Xây dựng công thức biến thiên hằng số cho phương trình vi phân phân thứ Caputo ngẫu nhiên dựa trên Định lý biểu diễn Itô.
- Đánh giá tốc độ phân tách tiệm cận giữa hai nghiệm phân biệt và suy ra tính không âm của số mũ Lyapunov bình phương trung bình.
- Thiết lập lược đồ số kiểu Euler-Maruyama xử lý triệt để tính kỳ dị của nhân $(t-s)^{\alpha-1}$ tại $s=t$, đồng thời chứng minh tốc độ hội tụ hiển của lược đồ này.
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 của hai nhánh toán học: Giải tích ngẫu nhiên và Giải tích phân thứ.
| Khái niệm / Công cụ lý thuyết | Định nghĩa / Biểu diễn toán học |
|---|---|
| Không gian xác suất và bộ lọc | Không gian xác suất đầy đủ $(\Omega, \mathcal{F}, \mathbb{P})$ với bộ lọc làm rộng $(\mathcal{F}t){t \ge 0}$ liên tục phải, sinh bởi chuyển động Brown $(W_t)_{t \ge 0}$ và các tập có độ đo 0. |
| Không gian $L^2$ và chuẩn bình phương trung bình | $X_t := L^2(\Omega, \mathcal{F}t, \mathbb{P})$ gồm các biến ngẫu nhiên nhận giá trị trong $\mathbb{R}^d$ với chuẩn $|f|{ms} := \sqrt{\mathbb{E}|f|^2} = \sqrt{\sum_{i=1}^d \mathbb{E}( |
| Định lý biểu diễn Itô | Với $f \in L^2(\Omega, \mathcal{F}_T, \mathbb{P})$, tồn tại duy nhất quá trình $\Xi \in M^2([0, T], \mathbb{R}^d)$ sao cho $f = \mathbb{E}(f) + \int_0^T \Xi(t) dW_t$. |
| Tích phân phân thứ Riemann–Liouville | $I_{0+}^\alpha x(t) := \frac{1}{\Gamma(\alpha)} \int_0^t (t-\tau)^{\alpha-1} x(\tau) d\tau$, với $\alpha \in (0, 1]$. |
| Đạo hàm phân thứ Caputo | $^C D_{0+}^\alpha x(t) := I_{0+}^{1-\alpha} \frac{dx(t)}{dt} = \frac{1}{\Gamma(1-\alpha)} \int_0^t (t-\tau)^{-\alpha} x'(\tau) d\tau$, với $\alpha \in (0, 1]$. |
| Hàm Mittag-Leffler một và hai tham số | $E_\alpha(M) := \sum_{k=0}^\infty \frac{M^k}{\Gamma(\alpha k + 1)}$; $E_{\alpha, \beta}(M) := \sum_{k=0}^\infty \frac{M^k}{\Gamma(\alpha k + \beta)}$ với $M \in \mathbb{R}^{d \times d}$. |
Phương pháp nghiên cứu thực tế
Luận án áp dụng các phương pháp phân tích định tính toán học và giải tích số chuyên sâu:
- Phương pháp điểm bất động Banach với chuẩn trọng số: Để chứng minh sự tồn tại và duy nhất nghiệm cổ điển và nghiệm nhẹ mà không cần chia nhỏ khoảng thời gian để kéo dài nghiệm, tác giả trang bị cho không gian Banach $H^2([0, T])$ chuẩn có trọng số dạng hàm Mittag-Leffler: $$|X|\gamma := \operatorname{ess,sup}{t \in [0, T]} \sqrt{\frac{\mathbb{E}(|X(t)|^2)}{E_{2\alpha-1}(\gamma t^{2\alpha-1})}}$$ với $\gamma > 0$ được chọn thích hợp để toán tử tích phân tương ứng trở thành một ánh xạ co.
- Kỹ thuật ước lượng tích phân và tính đẳng cự Itô: Khai thác tính chất đẳng cự Itô $\mathbb{E}\left|\int_0^t g(s) dW_s\right|^2 = \int_0^t \mathbb{E}|g(s)|^2 ds$ kết hợp bất đẳng thức Hölder để kiểm soát các số hạng ngẫu nhiên có nhân kỳ dị.
- Bất đẳng thức Gronwall cho phương trình vi phân phân thứ: Sử dụng để đánh giá sự phụ thuộc liên tục của nghiệm vào giá trị ban đầu và chặn trên của nghiệm xấp xỉ số.
- Phương pháp chứng minh phản chứng kết hợp phân rã tích phân tiệm cận: Dùng để thiết lập cận dưới cho tốc độ phân tách giữa hai nghiệm phân biệt khi $t \to \infty$.
- Kỹ thuật rời rạc hóa tránh điểm kỳ dị: Trong giải số, biến đổi nhân kỳ dị $(t-s)^{\alpha-1}$ thông qua phép chia lưới thời gian $t_k = \frac{kT}{n}$ và các toán tử trễ $\tau_n(s)$, $\rho_n(t)$ để xây dựng lược đồ số hiện tính được từng bước.
Nội dung chính theo từng chương
Chương 1: Một số kiến thức chuẩn bị
Chương 1 hệ thống hóa các khái niệm và kết quả nền tảng của giải tích ngẫu nhiên và giải tích phân thứ làm cơ sở cho các chương tiếp theo.
-
Giải tích ngẫu nhiên:
- Trình bày định nghĩa chuyển động Brown một chiều tiêu chuẩn $(W_t)_{t \ge 0}$, các tính chất quỹ đạo mẫu: liên tục Hölder địa phương cấp $\delta \in \left(0, \frac{1}{2}\right)$, không đâu khả vi và có biến phân vô hạn trên mọi khoảng con.
- Xây dựng tích phân ngẫu nhiên Itô $\int_0^T f(s) dW_s$ cho lớp hàm đơn giản $M_0([0, T], \mathbb{R})$ và mở rộng lên không gian $M^2([0, T], \mathbb{R})$, chứng minh tính đẳng cự Itô và tính tuyến tính.
- Nhắc lại Định lý biểu diễn Itô cho các biến ngẫu nhiên khả tích bình phương $\mathcal{F}_T$-đo được.
- Trình bày bài toán phương trình vi phân ngẫu nhiên $dX(t) = b(t, X(t))dt + \sigma(t, X(t))dW_t$, định lý tồn tại duy nhất nghiệm dưới điều kiện Lipschitz (N1) và điều kiện tăng trưởng không quá tuyến tính (N2).
- Nhắc lại lược đồ số Euler-Maruyama cho phương trình vi phân ngẫu nhiên bậc nguyên với đánh giá sai số bình phương trung bình: $$\mathbb{E} \sup_{0 \le t \le T} |X_n(t) - X(t)|^2 \le \frac{C_3}{n}$$
-
Giải tích phân thứ:
- Định nghĩa tích phân phân thứ Riemann-Liouville cấp $\alpha \in (0, 1]$ và đạo hàm phân thứ Caputo cấp $\alpha \in (0, 1]$. Khẳng định tính chất đạo hàm Caputo là nghịch đảo trái của tích phân Riemann-Liouville: $^C D_{0+}^\alpha I_{0+}^\alpha x(t) = x(t)$.
- Khảo sát phương trình vi phân phân thứ Caputo $^C D_{0+}^\alpha x(t) = f(t, x(t))$ và dạng tích phân tương đương $x(t) = x_0 + \frac{1}{\Gamma(\alpha)} \int_0^t (t-\tau)^{\alpha-1} f(\tau, x(\tau)) d\tau$, nhấn mạnh bản chất phụ thuộc toàn bộ quá khứ $[0, t)$.
- Đưa ra nghiệm giải tích của phương trình tuyến tính thuần nhất $^C D_{0+}^\alpha x(t) = Ax(t)$ dưới dạng $x(t) = E_\alpha(t^\alpha A)x_0$, và công thức biến thiên hằng số cho phương trình tất định $^C D_{0+}^\alpha x(t) = Ax(t) + f(x(t))$ biểu diễn qua hàm Mittag-Leffler hai tham số $E_{\alpha, \alpha}((t-\tau)^\alpha A)$.
Chương 2: Một số tính chất của nghiệm phương trình vi phân phân thứ Caputo ngẫu nhiên
Chương 2 là phần đóng góp lý thuyết trọng tâm về định tính nghiệm cho phương trình vi phân phân thứ Caputo ngẫu nhiên bậc $\alpha \in \left(\frac{1}{2}, 1\right)$ trên $[0, T]$.
-
Sự tồn tại và duy nhất nghiệm cổ điển: Xét phương trình: $$^C D_{0+}^\alpha X(t) = b(t, X(t)) + \sigma(t, X(t)) \frac{dW_t}{dt}, \quad X(0) = \eta \in X_0$$ dưới các giả thiết:
- $(H1)$ Điều kiện Lipschitz: $|b(t, x) - b(t, y)| \le L|x-y|$, $|\sigma(t, x) - \sigma(t, y)| \le L|x-y|$.
- $(H2)$ Tính bị chặn: $\int_0^T |b(t, 0)|^2 dt < \infty$, $\operatorname{ess,sup}_{t \in [0, T]} |\sigma(t, 0)| < \infty$.
Tác giả định nghĩa toán tử tích phân $T_\eta$ trên không gian Banach $H^2([0, T])$: $$T_\eta \xi(t) := \eta + \frac{1}{\Gamma(\alpha)} \int_0^t (t-\tau)^{\alpha-1} b(\tau, \xi(\tau)) d\tau + \frac{1}{\Gamma(\alpha)} \int_0^t (t-\tau)^{\alpha-1} \sigma(\tau, \xi(\tau)) dW_\tau$$ Bằng cách trang bị chuẩn trọng số $| \cdot |\gamma$ với hàm trọng số $E{2\alpha-1}(\gamma t^{2\alpha-1})$ và chọn hằng số $\gamma > \frac{3L^2(T+1)\Gamma(2\alpha-1)}{\Gamma(\alpha)^2}$, tác giả chứng minh $T_\eta$ là ánh xạ co, từ đó suy ra sự tồn tại và duy nhất nghiệm cổ điển toàn cục.
-
Sự phụ thuộc liên tục vào giá trị ban đầu: Với hai điều kiện ban đầu $\eta, \zeta \in X_0$, khoảng cách giữa hai nghiệm cổ điển $\varphi(t, \eta)$ và $\varphi(t, \zeta)$ thỏa mãn: $$|\varphi(t, \zeta) - \varphi(t, \eta)|{ms}^2 \le 3|\eta - \zeta|{ms}^2 E_{2\alpha-1}\left(\frac{(3T+3)L^2\Gamma(2\alpha-1)}{\Gamma(\alpha)^2} t^{2\alpha-1}\right)$$ Bất đẳng thức này chứng minh tính liên tục Lipschitz của nghiệm theo giá trị ban đầu trên đoạn thời gian hữu hạn.
-
Sự tồn tại và duy nhất nghiệm nhẹ: Xét phương trình bán tuyến tính $^C D_{0+}^\alpha X(t) = AX(t) + b(t, X(t)) + \sigma(t, X(t)) \frac{dW_t}{dt}$. Nghiệm nhẹ $Y(t)$ được định nghĩa thỏa mãn phương trình tích phân: $$Y(t) = E_\alpha(t^\alpha A)\eta + \int_0^t (t-\tau)^{\alpha-1} E_{\alpha, \alpha}((t-\tau)^\alpha A) b(\tau, Y(\tau)) d\tau + \int_0^t (t-\tau)^{\alpha-1} E_{\alpha, \alpha}((t-\tau)^\alpha A) \sigma(\tau, Y(\tau)) dW_\tau$$ Sử dụng toán tử $J_\eta$ và chuẩn trọng số $| \cdot |{\gamma_1}$ với $\gamma_1 > 2L^2 M_T^2 (T+1)\Gamma(2\alpha-1)$ (trong đó $M_T := \max{t \in [0, T]} |E_{\alpha, \alpha}(t^\alpha A)|$), tác giả chứng minh sự tồn tại và duy nhất nghiệm nhẹ dưới các điều kiện $(H1), (H2)$ yếu hơn các công bố trước đó.
-
Công thức biến thiên hằng số: Tác giả chứng minh nghiệm cổ điển $\varphi(t, \eta)$ và nghiệm nhẹ $\psi(t, \eta)$ trùng nhau hầu chắc chắn trên $[0, T]$. Ý tưởng chứng minh dựa trên việc áp dụng Định lý biểu diễn Itô cho phần tử ngẫu nhiên tổng quát, thiết lập ước lượng tích vô hướng: $$\mathbb{E} \left\langle \varphi(t, \eta) - \psi(t, \eta), C + \int_0^T \Xi(\tau) dW_\tau \right\rangle = 0$$ Từ đó suy ra công thức nghiệm tường minh cho hệ tuyến tính không thuần nhất $^C D_{0+}^\alpha X(t) = AX(t) + b(t) + \sigma(t) \frac{dW_t}{dt}$.
-
Cận dưới phân tách tiệm cận và số mũ Lyapunov:
- Đối với hai nghiệm phân biệt bất kỳ $\varphi(\cdot, \eta)$ và $\varphi(\cdot, \zeta)$ ($\eta \ne \zeta$), tác giả chứng minh bằng phương pháp phản chứng rằng với mọi $\epsilon > 0$: $$\limsup_{t \to \infty} t^{\frac{2\alpha}{1-\alpha} + \epsilon} |\varphi(t, \eta) - \varphi(t, \zeta)|_{ms} = \infty$$ Nghĩa là khoảng cách giữa hai nghiệm không thể tiến về 0 nhanh hơn tốc độ đa thức cấp $\frac{2\alpha}{1-\alpha}$.
- Ứng dụng kết quả trên cho phương trình song tuyến tính bị chặn $^C D_{0+}^\alpha X(t) = A(t)X(t) + B(t)X(t)\frac{dW_t}{dt}$, tác giả chứng minh số mũ Lyapunov bình phương trung bình của mọi nghiệm không tầm thường luôn không âm: $$\lambda_{ms}(\Phi(\cdot, \eta)) := \limsup_{t \to \infty} \frac{1}{t} \log |\Phi(t, \eta)|_{ms} \ge 0, \quad \forall \eta \ne 0$$
Chương 3: Lược đồ số kiểu Euler-Maruyama cho phương trình vi phân phân thứ Caputo ngẫu nhiên
Chương 3 tập trung xây dựng và phân tích lý thuyết giải tích số cho phương trình vi phân phân thứ Caputo ngẫu nhiên bậc $\alpha \in \left(\frac{1}{2}, 1\right)$.
-
Xây dựng lược đồ rời rạc hóa: Do nhân $(t-s)^{\alpha-1}$ có điểm kỳ dị tại $s = t$, tác giả đưa ra lược đồ rời rạc hóa trên lưới thời gian $t_k = \frac{kT}{n}$ ($k=0, 1, \dots, n$). Định nghĩa các hàm bước $\tau_n(s) := t_k$ khi $s \in [t_k, t_{k+1})$ và $\rho_n(t) := t_{k+1}$ khi $t \in (t_k, t_{k+1}]$. Nghiệm xấp xỉ số $X^{(n)}(t, \eta)$ được xác định bởi: $$X^{(n)}(t, \eta) = \eta + \frac{1}{\Gamma(\alpha)} \int_0^t \frac{b(\tau_n(s), X^{(n)}(\tau_n(s), \eta))}{(t-s)^{1-\alpha}} ds + \frac{1}{\Gamma(\alpha)} \int_0^t \frac{\sigma(\tau_n(s), X^{(n)}(\tau_n(s), \eta))}{(\rho_n(t) - \tau_n(s))^{1-\alpha}} dW_s$$ Lược đồ này cho phép tính toán truy hồi giá trị nghiệm tại các nút lưới mà không gặp phải tích phân suy rộng tại điểm kỳ dị.
-
Đánh giá tốc độ hội tụ: Tác giả thiết lập các bổ đề kỹ thuật chứng minh:
- Tính bị chặn đều bình phương trung bình: $\sup_{0 \le t \le T} |X^{(n)}(t, \eta)|{ms}^2 \le C_1$ với hằng số $C_1 = (1 + 3|\eta|{ms}^2) E_{2\alpha-1}\left(\frac{(6T+6)K^2\Gamma(2\alpha-1)}{\Gamma(\alpha)^2} T^{2\alpha-1}\right)$.
- Đánh giá khoảng cách giữa hai thời điểm: $|X^{(n)}(t, \eta) - X^{(n)}(\tilde{t}, \eta)|_{ms}^2 \le \frac{C_2}{n^{2\alpha-1}} + C_3 |t-\tilde{t}|^{2\alpha-1}$.
- Định lý hội tụ chính: Tồn tại hằng số $C$ độc lập với $n$ sao cho sai số bình phương trung bình giữa nghiệm số và nghiệm chính xác thỏa mãn: $$\sup_{0 \le t \le T} |X^{(n)}(t, \eta) - X(t, \eta)|_{ms}^2 \le \frac{C}{n^{2\alpha-1}}$$
- Kết quả này xác lập tốc độ hội tụ hiển cấp $\alpha - \frac{1}{2}$ theo bước lưới thời gian cho lược đồ số.
-
Lược đồ Euler-Maruyama mũ cho phương trình tuyến tính một chiều: Khảo sát phương trình tuyến tính một chiều, tác giả thiết lập lược đồ dạng mũ tương ứng, chứng minh tốc độ hội tụ và phân tích tính ổn định tiệm cận của nghiệm số.
Kết quả và những đóng góp mới
Đóng góp về mặt lý luận
- Khắc phục hạn chế trong lý thuyết tồn tại nghiệm: Chứng minh sự tồn tại và duy nhất nghiệm cổ điển toàn cục trên $[0, T]$ cho phương trình vi phân phân thứ Caputo ngẫu nhiên bậc $\alpha \in \left(\frac{1}{2}, 1\right)$, khắc phục triệt để lỗi logic về kéo dài nghiệm địa phương trong các công trình quốc tế trước đó (như Y. Wang et al., 2016).
- Giảm nhẹ điều kiện cho nghiệm nhẹ: Thiết lập định lý tồn tại và duy nhất nghiệm nhẹ đối với phương trình bán tuyến tính dưới các điều kiện Lipschitz và tăng trưởng yếu hơn so với các tài liệu đã công bố.
- Thiết lập công thức biến thiên hằng số ngẫu nhiên: Đưa ra biểu diễn tường minh của nghiệm cổ điển qua toán tử Mittag-Leffler một tham số $E_\alpha$ và hai tham số $E_{\alpha, \alpha}$, tạo công cụ giải tích hữu hiệu cho việc nghiên cứu hệ phi tuyến và hệ không thuần nhất.
- Phát hiện quy luật phân tách tiệm cận: Chứng minh khoảng cách giữa hai nghiệm phân biệt tiến về 0 không nhanh hơn tốc độ đa thức cấp $\frac{2\alpha}{1-\alpha}$, đồng thời khẳng định số mũ Lyapunov bình phương trung bình của nghiệm không tầm thường trong hệ song tuyến tính bị chặn luôn không âm ($\lambda_{ms} \ge 0$).
Đóng góp về mặt giải số và thực tiễn tính toán
- Đề xuất lược đồ số xử lý nhân kỳ dị: Xây dựng thành công lược đồ số kiểu Euler-Maruyama cho phương trình vi phân phân thứ Caputo ngẫu nhiên, giải quyết trở ngại tính toán do tính kỳ dị của nhân $(t-s)^{\alpha-1}$.
- Xác lập tốc độ hội tụ hiển: Chứng minh sai số bình phương trung bình của lược đồ số đạt cấp $\mathcal{O}\left(n^{-(2\alpha-1)}\right)$, cung cấp cơ sở toán học tường minh cho việc lập trình giải số các mô hình ngẫu nhiên có nhớ trong thực tế.
Bảng tổng hợp các kết quả chính của luận án
| Nội dung nghiên cứu | Đối tượng toán học | Kết quả chính đạt được | Công cụ / Kỹ thuật chứng minh |
|---|---|---|---|
| Nghiệm cổ điển | $^C D_{0+}^\alpha X(t) = b(t, X) + \sigma(t, X)\frac{dW_t}{dt}$ | Tồn tại, duy nhất toàn cục trên $[0, T]$ | Không gian Banach với chuẩn trọng số $E_{2\alpha-1}(\gamma t^{2\alpha-1})$, Điểm bất động Banach |
| Độ nhạy điều kiện ban đầu | $|\varphi(t, \zeta) - \varphi(t, \eta)|_{ms}^2$ | Phụ thuộc liên tục Lipschitz vào giá trị ban đầu | Bất đẳng thức Gronwall phân thứ |
| Nghiệm nhẹ & Biến thiên hằng số | Phương trình bán tuyến tính với toán tử ma trận $A$ | Biểu diễn nghiệm qua $E_\alpha, E_{\alpha, \alpha}$; chứng minh $\varphi(t) \equiv \psi(t)$ | Định lý biểu diễn Itô, chuẩn trọng số $|\cdot|_{\gamma_1}$ |
| Dáng điệu tiệm cận | Hai nghiệm phân biệt $\varphi(t, \eta), \varphi(t, \zeta)$ | $\limsup_{t \to \infty} t^{\frac{2\alpha}{1-\alpha}+\epsilon}|\Delta\varphi|{ms} = \infty$; $\lambda{ms} \ge 0$ | Phương pháp phản chứng, phân rã tích phân tiệm cận |
| Giải số Euler-Maruyama | Lược đồ số rời rạc hóa $X^{(n)}(t, \eta)$ | Tốc độ hội tụ hiển: $\mathcal{O}\left(n^{-(2\alpha-1)}\right)$ | Đánh giá kỳ vọng bình phương, xử lý dịch chuyển nhân kỳ dị |
Các kết quả trên đã được công bố trong 03 bài báo khoa học trên các tạp chí quốc tế chuyên ngành (bao gồm Stochastic Analysis and Applications, Statistics & Probability Letters, và Journal of Computational and Applied Mathematics).
Hạn chế và hướng nghiên cứu tiếp
Hạn chế được ghi nhận trong văn bản
- Luận án giới hạn phạm vi nghiên cứu ở bậc phân thứ $\alpha \in \left(\frac{1}{2}, 1\right)$. Điều kiện này xuất phát từ yêu cầu khả tích bình phương của nhân ngẫu nhiên Itô $\int_0^t (t-\tau)^{2\alpha-2} d\tau < \infty \Leftrightarrow 2\alpha - 2 > -1 \Leftrightarrow \alpha > \frac{1}{2}$. Trường hợp $\alpha \in \left(0, \frac{1}{2}\right]$ chưa được giải quyết trong khuôn khổ phương pháp này.
- Luận án mới giải quyết bài toán đặt chỉnh ở hai khía cạnh: sự tồn tại duy nhất nghiệm và sự phụ thuộc liên tục vào điều kiện ban đầu; tính liên tục của nghiệm theo thời gian $t$ và theo bậc phân thứ $\alpha$ chưa được hoàn thiện trọn vẹn trong văn bản.
Hướng nghiên cứu tiếp theo
- Tiếp tục nghiên cứu tính liên tục của nghiệm cổ điển theo thời gian và theo tham số bậc phân thứ $\alpha$ để hoàn thiện trọn vẹn bài toán đặt chỉnh (well-posedness).
- Mở rộng lý thuyết và phương pháp số cho các phương trình vi phân phân thứ ngẫu nhiên có cấu trúc phức tạp hơn hoặc chịu tác động của các quá trình nhiễu khác ngoài chuyển động Brown tiêu chuẩn (như chuyển động Brown phân thứ, nhiễu Lévy).
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, Lý thuyết xác suất và Giải tích số:
- Đối với nghiên cứu sinh và học viên cao học ngành Toán giải tích, Toán ứng dụng, Xác suất - Thống kê: Cung cấp kỹ thuật xây dựng không gian Banach với chuẩn trọng số Mittag-Leffler, kỹ thuật sử dụng Định lý biểu diễn Itô để thiết lập công thức biến thiên hằng số, và phương pháp phân tích tiệm cận cho các hệ phương trình có nhớ.
- Đối với giảng viên và nhà nghiên cứu giải tích số: Cung cấp phương pháp luận chi tiết để thiết kế và đánh giá sai số cho các lược đồ số giải phương trình vi phân ngẫu nhiên có nhân tích phân kỳ dị.
- Đối với các nhà nghiên cứu trong Cơ học, Vật lý kỹ thuật, Điều khiển học: Cung cấp công cụ toán học tường minh để mô hình hóa và giải số các quá trình tiến hóa có tính nhớ và chịu tác động của nhiễu ngẫu nhiên.
Câu hỏi thường gặp
1. Tại sao luận án lựa chọn đạo hàm phân thứ Caputo thay vì đạo hàm Riemann-Liouville?
Đạo hàm phân thứ Riemann-Liouville dẫn đến các bài toán giá trị ban đầu với điều kiện ban đầu chứa các giới hạn tích phân phân thứ tại $t=0$, vốn không có ý nghĩa vật lý cụ thể. Ngược lại, bài toán giá trị ban đầu đối với đạo hàm phân thứ Caputo sử dụng các điều kiện ban đầu chứa đạo hàm bậc nguyên thông thường ($x(0) = x_0$), có ý nghĩa vật lý rõ ràng và thuận lợi cho việc mô hình hóa các quá trình thực tế.
2. Tại sao điều kiện bậc phân thứ trong luận án lại yêu cầu $\alpha \in \left(\frac{1}{2}, 1\right)$?
Khi tính kỳ vọng bình phương trung bình của tích phân ngẫu nhiên Itô chứa nhân phân thứ $\int_0^t (t-\tau)^{\alpha-1} \sigma(\tau, X(\tau)) dW_\tau$, tính đẳng cự Itô dẫn đến việc lấy tích phân của hàm $(t-\tau)^{2\alpha-2}$. Tích phân này hội tụ khi và chỉ khi $2\alpha - 2 > -1$, tương đương $\alpha > \frac{1}{2}$. Do đó, khoảng tham số $\alpha \in \left(\frac{1}{2}, 1\right)$ là điều kiện tự nhiên để tích phân ngẫu nhiên xác định tốt trong không gian $L^2$.
3. Hạn chế trong chứng minh nghiệm toàn cục của các nghiên cứu trước (như Y. Wang et al., 2016) là gì và luận án khắc phục ra sao?
Nghiên cứu của Y. Wang et al. (2016) đã chứng minh nghiệm địa phương trên $[0, T_a]$ rồi dùng phép tịnh tiến thời gian để kéo dài nghiệm ra toàn khoảng $[0, \infty)$ tương tự như phương trình vi phân thường. Tuy nhiên, do phương trình vi phân phân thứ có tính phụ thuộc toàn bộ lịch sử quá khứ, nghiệm trên đoạn $[T^, T^+\delta)$ phụ thuộc vào giá trị trên $[0, T^*)$, nên phép dịch chuyển thời gian không bảo toàn phương trình. Luận án khắc phục hạn chế này bằng cách chứng minh trực tiếp sự tồn tại duy nhất nghiệm toàn cục trên toàn đoạn $[0, T]$ thông qua việc thiết lập chuẩn trọng số Mittag-Leffler trên không gian Banach $H^2([0, T])$.
4. Công thức biến thiên hằng số cho phương trình vi phân phân thứ Caputo ngẫu nhiên được chứng minh bằng công cụ nào?
Công thức biến thiên hằng số được chứng minh bằng cách chỉ ra nghiệm cổ điển $\varphi(t, \eta)$ trùng với nghiệm nhẹ $\psi(t, \eta)$. Tác giả sử dụng Định lý biểu diễn Itô để biểu diễn phần tử ngẫu nhiên tổng quát dưới dạng tích phân đối với chuyển động Brown, sau đó thiết lập ước lượng tích vô hướng $\mathbb{E}\left\langle \varphi(t, \eta) - \psi(t, \eta), C + \int_0^T \Xi(\tau) dW_\tau \right\rangle = 0$ với mọi hằng số $C$ và quá trình $\Xi \in M^2([0, T], \mathbb{R}^d)$.
5. Lược đồ Euler-Maruyama trong luận án đã xử lý khó khăn gì từ nhân tích phân $(t-s)^{\alpha-1}$?
Nhân $(t-s)^{\alpha-1}$ trở thành vô hạn (kỳ dị) tại $s = t$ khi $\alpha \in \left(\frac{1}{2}, 1\right)$, khiến các lược đồ Euler thông thường không áp dụng trực tiếp được. Luận án giải quyết bằng cách đưa vào các hàm bước $\tau_n(s)$ và $\rho_n(t)$ trên phân hoạch lưới thời gian, xấp xỉ số hạng khuếch tán ngẫu nhiên qua mẫu số $(\rho_n(t) - \tau_n(s))^{1-\alpha}$, đảm bảo mẫu số luôn dương và tách rời điểm kỳ dị, từ đó thu được lược đồ tính toán tường minh từng bước với tốc độ hội tụ hiển $\mathcal{O}\left(n^{-(2\alpha-1)}\right)$.
Kết luận
Luận án tiến sĩ của Phan Thị Hương đã giải quyết một cách có hệ thống các bài toán nền tảng về lý thuyết định tính và giải tích số cho phương trình vi phân phân thứ Caputo ngẫu nhiên bậc $\alpha \in \left(\frac{1}{2}, 1\right)$. Bằng việc vận dụng chuẩn trọng số Mittag-Leffler và Định lý biểu diễn Itô, tác giả đã thiết lập các định lý về sự tồn tại, duy nhất nghiệm cổ điển, nghiệm nhẹ, xây dựng công thức biến thiên hằng số, xác định cận dưới phân tách tiệm cận và tính không âm của số mũ Lyapunov. Đồng thời, luận án đã đề xuất lược đồ số kiểu Euler-Maruyama giải quyết tính kỳ dị của nhân phân thứ và chứng minh được tốc độ hội tụ hiển của phương pháp. Các kết quả này đóng góp trực tiếp vào sự phát triển của chuyên ngành Toán ứng dụng và giải tích ngẫu nhiên.