Tổng quan luận án
Luận án Tiến sĩ Toán học (Toán ứng dụng), chuyên ngành Lý thuyết Xác suất và Thống kê toán học (mã số 9460112.02) của nghiên cứu sinh Nguyễn Văn Tân, được hoàn thành năm 2023 tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội, dưới sự hướng dẫn khoa học của PGS. Nguyễn Tiến Dũng và TS. Trần Mạnh Cường. Luận án nghiên cứu về việc ứng dụng giải tích Malliavin trong việc khảo sát tốc độ hội tụ và các tính chất mật độ nghiệm của một số lớp phương trình vi phân ngẫu nhiên (SDE).
Tính cấp thiết và khoảng trống nghiên cứu
Trong lý thuyết xác suất hiện đại, giải tích Malliavin (phép tính vi phân ngẫu nhiên vô hạn chiều trên không gian Wiener) cung cấp các kỹ thuật hữu hiệu để nghiên cứu cấu trúc phân bố xác suất và tính chính quy của các phiếm hàm ngẫu nhiên. Hai hướng nghiên cứu lý thuyết nổi bật sử dụng giải tích Malliavin gồm có:
- Chứng minh sự tồn tại, tính trơn ($C^\infty$) và thiết lập các ước lượng cận Gauss hai phía cho mật độ xác suất cùng xác suất đuôi của nghiệm phương trình vi phân ngẫu nhiên.
- Thiết lập định lý giới hạn trung tâm và đánh giá tốc độ hội tụ của nghiệm các phương trình vi phân ngẫu nhiên.
Từ góc độ văn bản của tác giả, các khoảng trống nghiên cứu cụ thể được xác định như sau:
- Đối với xấp xỉ Smoluchowski-Kramers cho phương trình vi phân ngẫu nhiên bậc hai mô tả chuyển động của hạt có khối lượng rất nhỏ $\mu \to 0$, các nghiên cứu trước đây (như công trình của Freidlin) chủ yếu dừng lại ở việc chứng minh sự hội tụ theo xác suất và hội tụ trong không gian $L^2(\Omega)$ với tốc độ cấp $\mathcal{O}(\mu)$. Tuy nhiên, các thông tin về phân bố xác suất thông qua khoảng cách Kolmogorov $\sup_{x \in \mathbb{R}} |P(X_t^\mu \le x) - P(X_t \le x)|$ vẫn chưa có ước lượng Berry-Esseen hiển.
- Đối với phương trình vi phân ngẫu nhiên có trễ và lược đồ giải số xấp xỉ Carathéodory, các đánh giá tốc độ hội tụ yếu trước đây thường đòi hỏi hàm thử phải thỏa mãn điều kiện liên tục Lipschitz. Bài toán ước lượng tốc độ hội tụ yếu đối với lớp các hàm thử chỉ cần đo được và bị chặn chưa có nhiều kết quả cụ thể.
- Đối với phương trình vi phân hàm ngẫu nhiên phân thứ (SFDE) điều khiển bởi chuyển động Brown phân thứ ($B^H$) với tham số Hurst $H > 1/2$, các nghiên cứu trước đây chỉ khảo sát trường hợp không có số hạng tích phân trễ phân bố. Việc thiết lập sự tồn tại nghiệm, ước lượng cận Gauss hai phía và tính trơn của mật độ khi đồng thời xuất hiện số hạng tích phân trễ và các dạng nhiễu cộng tính cũng như nhiễu nhân tính đòi hỏi kỹ thuật phân tích phức tạp chưa được giải quyết đầy đủ.
Mục tiêu nghiên cứu
Luận án xác định ba mục tiêu nghiên cứu cụ thể:
- Thiết lập một ước lượng Berry-Esseen hiển cho tốc độ hội tụ theo khoảng cách Kolmogorov của xấp xỉ Smoluchowski-Kramers đối với phương trình vi phân ngẫu nhiên khi khối lượng hạt $\mu \to 0$.
- Nghiên cứu sự hội tụ yếu của nghiệm phương trình vi phân ngẫu nhiên có trễ liên quan đến tham số thời gian trễ và áp dụng trực tiếp để thu được tốc độ hội tụ yếu cho lược đồ xấp xỉ Carathéodory đối với mọi hàm đo được bị chặn.
- Chứng minh tính trơn và thiết lập ước lượng Gauss cận trên, cận dưới cho mật độ xác suất của nghiệm phương trình vi phân hàm ngẫu nhiên phân thứ với $H > 1/2$ trong cả hai trường hợp nhiễu cộng tính và nhiễu nhân tính.
Đối tượng và phạm vi nghiên cứu
- Đối tượng nghiên cứu:
- Nghiệm của phương trình vi phân ngẫu nhiên bậc hai mô tả chuyển động hạt có khối lượng nhỏ và phương trình vi phân ngẫu nhiên khuếch tán giới hạn.
- Nghiệm của phương trình vi phân ngẫu nhiên có trễ thời gian và dãy nghiệm xấp xỉ Carathéodory.
- Nghiệm của phương trình vi phân hàm ngẫu nhiên điều khiển bởi chuyển động Brown phân thứ $B^H$ với tham số Hurst $H \in (1/2, 1)$.
- Phạm vi nghiên cứu: Luận án được thực hiện trên không gian xác suất đầy đủ $(\Omega, \mathcal{F}, P)$, trên đoạn thời gian hữu hạn $[0, T]$ và không gian trạng thái thực $\mathbb{R}$. Các hàm hệ số được giả định thỏa mãn tính liên tục Lipschitz, điều kiện tăng trưởng tuyến tính, các điều kiện chính quy về đạo hàm riêng cấp 1, cấp 2 hoặc khả vi vô hạn lần, cùng điều kiện suy biến đều của hệ số khuếch tán.
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 văn bản
Tác giả đã hệ thống hóa các tài liệu nền tảng và các công trình chuyên sâu liên quan trực tiếp đến đề tài:
- Lý thuyết nền tảng giải tích Malliavin: Khởi nguồn từ công trình của Malliavin (1978a, b, c), sau đó được phát triển và hoàn thiện bởi Shigekawa (1980), Bismut (1981), Stroock (1981a, b), Ikeda & Watanabe (1984), Bichteler và cộng sự (1987), Malliavin (1991), Sanz-Solé (2005), Malliavin & Thalmaier (2005), Nualart (2006), Di Nunno và cộng sự (2009), Nourdin & Peccati (2012), Ishikawa (2016).
- Xấp xỉ Smoluchowski-Kramers: Xuất phát từ các mô hình vật lý của Smoluchowski (1916) và Kramers (1940), lý thuyết động học ngẫu nhiên của Nelson (1967). Về mặt toán học ngẫu nhiên, các kết quả nền tảng được thiết lập bởi Freidlin (2004), Cerrai & Freidlin (2006), Cerrai (2014, 2017), Hottovy và cộng sự (2015), Nguyen Tien Dung & Privault (2019).
- Phương trình vi phân ngẫu nhiên có trễ và xấp xỉ Carathéodory: Xuất phát từ phương pháp giải số Carathéodory (1955) cho phương trình tất định, Bell & Mohammed (1989) lần đầu mở rộng cho SDE. Tiếp theo là các nghiên cứu của Turo (1996) cho phương trình vi phân hàm ngẫu nhiên, Mao & Liu (1998) cho phương trình tiến hóa trễ bán tuyến tính, Ferrante & Rovira (2010) cho phương trình có trễ điều khiển bởi fBm, Faizullah (2012) dưới chuyển động G-Brownian, Benabdallah & Bourza (2019) cho SDE với cận phản xạ, Mao (2008) cho các lớp SDE trễ tổng quát, Hu, Mohammed & Yan (2008) về lược đồ Euler cho phương trình trễ.
- Tích phân phân thứ và ước lượng mật độ: Định nghĩa tích phân đường ngẫu nhiên theo fBm bởi Zähle (1998), Samko và cộng sự (1993). Các công trình về ước lượng mật độ và tính trơn của nghiệm SDE phân thứ bởi Nualart & Saussereau (2009), Nourdin & Viens (2009), Baudoin & Ouyang (2014), Cass, Hairer, Litterer & Tindel (2015), Takeuchi (2007), cùng các công trình của Nguyen Tien Dung (2013, 2016, 2019), Nguyen Tien Dung và cộng sự (2019).
Vị trí và khoảng trống luận án giải quyết
Luận án kế thừa trực tiếp các kỹ thuật giải tích Malliavin kết hợp với phương pháp ước lượng tốc độ hội tụ yếu gần đây của PGS. Nguyễn Tiến Dũng để giải quyết ba vấn đề cụ thể:
- Đưa ra chặn trên tường minh cấp $\mathcal{O}(\mu^{1/2})$ cho khoảng cách Kolmogorov trong xấp xỉ Smoluchowski-Kramers, bổ sung thông tin phân bố xác suất mà các đánh giá $L^2$ trước đây chưa phản ánh.
- Đưa ra đánh giá hội tụ yếu cấp $\mathcal{O}(n^{-1/2})$ cho xấp xỉ Carathéodory trên không gian hàm đo được bị chặn, vượt qua giới hạn của các kỹ thuật giải tích số truyền thống chỉ áp dụng cho hàm Lipschitz.
- Xây dựng công thức đổi biến khử nhiễu ngẫu nhiên phân thứ để thiết lập sự tồn tại nghiệm mạnh duy nhất, chặn Gauss hai phía cấp $t^{2H}$ và tính trơn vô hạn cấp cho mật độ nghiệm của phương trình vi phân hàm ngẫu nhiên có trễ tích phân với $H > 1/2$.
Cơ sở lý thuyết và phương pháp nghiên cứu
Cơ sở lý thuyết và khung phân tích
Luận án sử dụng các cấu trúc toán học của lý thuyết giải tích ngẫu nhiên hiện đại:
- Khai triển Wiener-Itô: Không gian xác suất $(\Omega, \mathcal{F}, P)$ với dòng thông tin tự nhiên $\mathcal{F} = {\mathcal{F}t, t \in [0, T]}$ sinh bởi chuyển động Brown $B = (B_t){t \in [0, T]}$. Mọi biến ngẫu nhiên $\xi \in L^2(\Omega, \mathcal{F}T, P)$ đều có biểu diễn duy nhất $\xi = \sum{n=0}^\infty I_n(f_n)$ với $f_n \in L^2([0, T]^n)$ là các hàm đối xứng, thỏa mãn hệ thức đẳng cự $|\xi|{L^2(P)}^2 = \sum{n=0}^\infty n! |f_n|_{L^2([0, T]^n)}^2$.
- Toán tử tích phân Skorohod ($\delta$): Mở rộng của tích phân Itô cho các quá trình ngẫu nhiên $u \in \text{Dom}(\delta)$ không nhất thiết thích ứng với bộ lọc $\mathcal{F}t$, xác định qua công thức $\delta(u) = \sum{n=0}^\infty I_{n+1}(\tilde{f}_n)$ với $\tilde{f}_n$ là hàm đối xứng hóa của $f_n$.
- Toán tử đạo hàm Malliavin ($D$): Toán tử vi phân xác định trên không gian Sobolev ngẫu nhiên $\mathbb{D}^{1,2}$ bởi $D_t F = \sum_{n=1}^\infty n I_{n-1}(f_n(\cdot, t))$. Định nghĩa không gian $\mathbb{D}^{k,p}$ với chuẩn $|F|{k,p}^p = \mathbb{E}|F|^p + \sum{j=1}^k \mathbb{E}\left[ \int_{[0, T]^j} |D_{t_1,\dots,t_j}^j F|^2 dt_1 \dots dt_j \right]^{p/2}$ và không gian $\mathbb{D}^\infty = \bigcap_{k \ge 1, p > 1} \mathbb{D}^{k,p}$.
- Tính chất đối ngẫu Malliavin-Skorohod: Hệ thức cơ bản $\mathbb{E}\left[ F \int_0^T u(t)\delta B_t \right] = \mathbb{E}\left[ \int_0^T u(t) D_t F dt \right]$ với mọi $F \in \mathbb{D}^{1,2}$ và $u \in \text{Dom}(\delta)$.
- Khung phân tích mật độ Nourdin-Viens: Nếu $F \in \mathbb{D}^{1,2}$ thỏa mãn $\mathbb{E}[F] = 0$ và tồn tại các hằng số dương $c, C$ sao cho $c \le \int_0^T D_t F \cdot \mathbb{E}[D_t F | \mathcal{F}_t] dt \le C$ hầu chắc chắn, thì $F$ có mật độ xác suất $p_F(x)$ thỏa mãn ước lượng Gauss hai phía: $$\frac{\mathbb{E}|F|}{2C} \exp\left( -\frac{x^2}{2c} \right) \le p_F(x) \le \frac{\mathbb{E}|F|}{2c} \exp\left( -\frac{x^2}{2C} \right), \quad \forall x \in \mathbb{R}.$$
- Tích phân phân thứ dạng Volterra và Zähle: Chuyển động Brown phân thứ $B^H$ với tham số Hurst $H \in (1/2, 1)$ có hàm hiệp phương sai $R_H(t, s) = \frac{1}{2}(t^{2H} + s^{2H} - |t-s|^{2H})$ và biểu diễn Volterra $B_t^H = \int_0^t K_H(t, s) dB_s$, trong đó nhân $K_H(t, s) = C_H s^{1/2-H} \int_s^t (u-s)^{H-3/2} u^{H-1/2} du$ với $C_H = \left[ \frac{H(2H-1)}{\mathcal{B}(2-2H, H-1/2)} \right]^{1/2}$.
Phương pháp nghiên cứu thực tế
Luận án là công trình toán học lý thuyết định tính và giải tích hàm ngẫu nhiên. Các phương pháp suy luận và công cụ kỹ thuật bao gồm:
- Phương pháp xấp xỉ hàm trơn và đối ngẫu: Xấp xỉ hàm chỉ tiêu hoặc hàm đo được bị chặn bởi dãy hàm trơn có giá compact $g_\varepsilon$, áp dụng quy tắc xích cho đạo hàm Malliavin, kết hợp công thức đối ngẫu giữa đạo hàm Malliavin và tích phân ngẫu nhiên Skorohod để chuyển sai phân kỳ vọng về tích phân chứa đạo hàm Malliavin cấp 1 và cấp 2.
- Phương pháp đổi biến phi tuyến giải tích: Sử dụng hàm thế $F(t, x) = \int_0^x \frac{1}{\sigma(t, u)} du$ và định lý hàm ẩn để chuyển phương trình vi phân hàm ngẫu nhiên phân thứ với nhiễu nhân tính về phương trình với nhiễu cộng tính có hệ số dịch chuyển trôi thỏa mãn điều kiện Lipschitz.
- Kỹ thuật ước lượng moment và giải tích bất đẳng thức: Sử dụng hệ thống các bất đẳng thức giải tích ngẫu nhiên gồm bất đẳng thức Burkholder-Davis-Gundy (BDG), bất đẳng thức Hölder, bổ đề Gronwall, bổ đề Gronwall dạng trễ (Gronwall-type), định lý so sánh nghiệm phương trình vi phân, định lý biểu diễn Martingale Dubins-Schwarz, định lý Fernique cho quá trình Gauss và bất đẳng thức Markov.
Nguồn tư liệu và dữ liệu
Nguồn tư liệu hoàn toàn bao gồm các công trình nghiên cứu toán học lý thuyết đã công bố trên các tạp chí quốc tế chuyên ngành xác suất, các sách chuyên khảo về giải tích Malliavin và phương trình vi phân ngẫu nhiên. Dữ liệu trong luận án là các cấu trúc giải tích, công thức nghiệm, các chặn giải tích và định lý toán học được chứng minh chặt chẽ.
Nội dung chính theo từng chương
Chương 1: Giải tích Malliavin
Chương 1 tổng hợp các kiến thức cơ sở về giải tích Malliavin cần thiết cho các chương sau:
- Khai triển Wiener-Itô và Tích phân Skorohod: Định nghĩa tích phân Itô lặp $I_n(f_n) = \int_0^T \dots \int_0^{t_2} f_n(t_1, \dots, t_n) dB_{t_1} \dots dB_{t_n}$ trên tập $S_n = {(t_1, \dots, t_n) \in [0, T]^n: 0 \le t_1 < \dots < t_n \le T}$; thiết lập tính chất đối xứng hóa hàm $f_n$ và công thức đẳng cự Wiener-Itô; định nghĩa toán tử Skorohod $\delta(u) = \int_0^T u(t) \delta B_t$ và chứng minh tính mở rộng của nó đối với tích phân Itô khi $u$ là quá trình thích ứng.
- Đạo hàm Malliavin và các quy tắc vi phân: Định nghĩa đạo hàm Malliavin $D_t F$, xây dựng các không gian $\mathbb{D}^{1,2}$, $\mathbb{D}^{k,p}$, chứng minh quy tắc tích $D_t(F_1 F_2) = F_1 D_t F_2 + F_2 D_t F_1$ và quy tắc xích cho hàm Lipschitz $\varphi$: $D_t \varphi(F) = \varphi'(F) D_t F$. Trình bày công thức Clark-Ocone và định lý biểu diễn đối ngẫu.
- Tính khả vi Malliavin của nghiệm SDE: Khảo sát phương trình vi phân ngẫu nhiên chuẩn $X_t = x_0 + \int_0^t a(s, X_s)ds + \int_0^t \sigma(s, X_s)dB_s$. Dưới các điều kiện Lipschitz và tăng trưởng tuyến tính, nghiệm $X_t \in \mathbb{D}^{1,2}$ với đạo hàm Malliavin thỏa mãn phương trình vi phân ngẫu nhiên tuyến tính: $$D_\theta X_t = \sigma(\theta, X_\theta) + \int_\theta^t \tilde{a}(s) D_\theta X_s ds + \int_\theta^t \tilde{\sigma}(s) D_\theta X_s dB_s, \quad 0 \le \theta \le t \le T.$$
- Tiêu chuẩn tồn tại và tính trơn của mật độ: Trình bày tiêu chuẩn Malliavin: nếu $F \in \mathbb{D}^\infty$ và $|DF|{L^2[0, T]}^{-1} \in \bigcap{p > 1} L^p(\Omega)$ thì $F$ có mật độ xác suất thuộc lớp $C^\infty(\mathbb{R})$; thiết lập các mệnh đề đánh giá cận trên và cận dưới Gauss cho mật độ dựa trên đại lượng $\int_0^T D_t F \mathbb{E}[D_t F | \mathcal{F}_t] dt$.
Chương 2: Ước lượng tốc độ hội tụ cho xấp xỉ Smoluchowski-Kramers và nghiệm của phương trình vi phân ngẫu nhiên có trễ
Chương 2 trình bày hai kết quả lý thuyết chính của luận án:
1. Ước lượng Berry-Esseen trong xấp xỉ Smoluchowski-Kramers
Khảo sát phương trình vi phân ngẫu nhiên bậc hai với khối lượng hạt $\mu > 0$: $$\mu \ddot{X}_t^\mu = b(t, X_t^\mu) + \sigma(t, X_t^\mu)\dot{B}_t - \dot{X}_t^\mu, \quad X_0^\mu = x_0, \quad \dot{X}_0^\mu = y_0.$$ Bằng công thức tích phân từng phần, phương trình được chuyển về dạng tích phân tương đương: $$X_t^\mu = x_0 + \mu y_0(1 - e^{-t/\mu}) + \int_0^t b(s, X_s^\mu)ds + \int_0^t \sigma(s, X_s^\mu)dB_s - \mu \int_0^t e^{-\frac{t-s}{\mu}} b(s, X_s^\mu)ds - \int_0^t e^{-\frac{t-s}{\mu}} \sigma(s, X_s^\mu)dB_s.$$ Khi $\mu \to 0$, $X_t^\mu$ hội tụ đến nghiệm $X_t$ của phương trình khuếch tán $dX_t = b(t, X_t)dt + \sigma(t, X_t)dB_t$.
- Giả thiết:
- $(B_1)$: $b, \sigma: [0, T] \times \mathbb{R} \to \mathbb{R}$ thỏa mãn tính liên tục Lipschitz và tăng trưởng tuyến tính với hằng số $L > 0$.
- $(B_2)$: $b(t, \cdot), \sigma(t, \cdot) \in C^2(\mathbb{R})$ có các đạo hàm cấp một và cấp hai bị chặn.
- Điều kiện elliptic đều: $|\sigma|0 := \inf{(t, x) \in [0, T] \times \mathbb{R}} |\sigma(t, x)| > 0$.
- Các kết quả bổ trợ:
- Đánh giá moment sai phân: $\mathbb{E}|X_t^\mu - X_t|^p \le C \mu^{p/2}$ với mọi $p \ge 2, \mu \in (0, 1)$.
- Tính khả vi Malliavin và đánh giá sai phân đạo hàm Malliavin: $\mathbb{E}|DX_t^\mu - DX_t|_{L^2[0, T]}^2 \le C \mu$.
- Đánh giá moment nghịch đảo đạo hàm Malliavin: $\mathbb{E}\left[ |DX_t|_{L^2[0, T]}^{-2p} \right] \le C t^{-p}$ với mọi $p \ge 1$.
- Định lý 2.6 (Kết quả chính): Với các giả thiết trên, với mỗi $t \in (0, T]$, khoảng cách Kolmogorov thỏa mãn ước lượng Berry-Esseen hiển: $$\sup_{x \in \mathbb{R}} |P(X_t^\mu \le x) - P(X_t \le x)| \le C t^{-1/2} \sqrt{\mu}, \quad \forall \mu \in (0, 1),$$ trong đó hằng số $C > 0$ không phụ thuộc vào $t$ và $\mu$.
2. Sự hội tụ yếu của nghiệm SDE có trễ và xấp xỉ Carathéodory
Khảo sát hệ phương trình vi phân ngẫu nhiên có trễ: $$X_{\tau}(t) = \varphi(0) + \int_0^t b(s, X_{\tau}(s), X_{\tau}(s - \tau))ds + \int_0^t \sigma(s, X_{\tau}(s), X_{\tau}(s - \tau))dB_s, \quad t \in [0, T],$$ với điều kiện ban đầu $X_\tau(t) = \varphi(t)$ trên $[-\tau, 0]$, trong đó $\varphi$ là hàm tất định bị chặn.
- Giả thiết:
- $(H_1)$: $b, \sigma: [0, T] \times \mathbb{R}^2 \to \mathbb{R}$ thỏa mãn tính liên tục Lipschitz và tăng trưởng tuyến tính theo hằng số $L$.
- $(H_2)$: $b(t, \cdot, \cdot), \sigma(t, \cdot, \cdot) \in C^2(\mathbb{R}^2)$ với các đạo hàm riêng bị chặn bởi $L$.
- Điều kiện elliptic: $|\sigma(t, x, y)| \ge \sigma_0 > 0$.
- Bổ đề cơ bản (Bổ đề 2.7): Cho $F_1, F_2 \in \mathbb{D}^{1,2}$ với $|DF_1|{L^2[0, T]} > 0$ h.c. và $g \in \mathcal{G} = {g: \mathbb{R} \to \mathbb{R} \text{ đo được}, |g|\infty \le 1}$, ta có: $$|\mathbb{E} g(F_1) - \mathbb{E} g(F_2)| \le C \left( \mathbb{E}|DF_1|{L^2[0, T]}^{-2} \cdot \mathbb{E}\int_0^T \int_0^T |D\theta D_r F_1|^2 d\theta dr + (\mathbb{E}|DF_1|{L^2[0, T]}^{-4})^{1/2} \right)^{1/2} |F_1 - F_2|{L^2(\Omega)}.$$
- Định lý 2.20: Cho $0 < \tau_1 < \tau_2$. Với mọi $g \in \mathcal{G}$ và $t \in (0, T]$: $$|\mathbb{E} g(X_{\tau_1}(t)) - \mathbb{E} g(X_{\tau_2}(t))| \le C_{T, L} \left( |\tau_1 - \tau_2|^2 + t^{-1} \int_0^{t \wedge \tau_1} |\varphi(s - \tau_1) - \varphi(s - \tau_2)|^4 ds + t^{-1} \int_{t \wedge \tau_1}^{t \wedge \tau_2} |\varphi(0) - \varphi(s - \tau_2)|^4 ds \right)^{1/4}.$$
- Hệ quả 2.21: Nếu $\varphi$ liên tục Hölder bậc $\beta \in (0, 1)$, thì: $$\sup_{0 \le t \le T} |\mathbb{E} g(X_{\tau_1}(t)) - \mathbb{E} g(X_{\tau_2}(t))| \le C_{T, L} |\tau_1 - \tau_2|^{\frac{\beta}{4} \wedge \frac{1}{2}}.$$
- Hệ quả 2.22 (Ứng dụng cho xấp xỉ Carathéodory): Xét xấp xỉ Carathéodory $x^n(t) = x_0 + \int_0^t b(s, x^n(s - 1/n))ds + \int_0^t \sigma(s, x^n(s - 1/n))dB_s$. Khi đó với mọi hàm $g \in \mathcal{G}$: $$\sup_{0 \le t \le T} |\mathbb{E} g(x^n(t)) - \mathbb{E} g(x(t))| \le \frac{C}{\sqrt{n}}, \quad \forall n \ge 1.$$
Chương 3: Tính trơn và ước lượng mật độ của nghiệm phương trình vi phân hàm ngẫu nhiên phân thứ
Chương 3 nghiên cứu lớp phương trình vi phân hàm ngẫu nhiên phân thứ (SFDE) tổng quát: $$X_t = \eta(0) + \int_0^t \left( \int_{-r}^0 \rho(u - s) g(u, X_{u+s}) du + a(s, X_s) \right) ds + \int_0^t \sigma(s, X_s) dB_s^H, \quad t \in [0, T],$$ với điều kiện ban đầu $X_t = \eta(t)$ trên $[-r, 0]$, trong đó $H \in (1/2, 1)$ và tích phân theo $B^H$ được hiểu theo nghĩa Riemann-Stieltjes/Zähle.
1. Trường hợp nhiễu cộng tính ($\sigma(s, x) = \sigma_s$)
Phương trình có dạng: $$X_t = \eta(0) + \int_0^t \left( \int_{-r}^0 \rho(u - s) g(X_{u+s}) du + a(s, X_s) \right) ds + \int_0^t \sigma_s dB_s^H.$$
- Giả thiết: $(A_1)$ $\rho \in L^1([-r, 0])$; $(A_2)$ $a, g$ liên tục Lipschitz; $(A_3)$ $\sigma$ là hàm tất định bị chặn thỏa mãn $\int_0^T \left( \int_s^T |\sigma_r| K_H(r, s) dr \right)^2 ds < \infty$.
- Đạo hàm Malliavin: $X_t \in \mathbb{D}^{1,2}$ và thỏa mãn: $$D_\theta X_t = \int_\theta^t \left( \int_{-r}^0 \rho(u - s) \tilde{g}(u) D_\theta X_{u+s} du + \tilde{a}(s) D_\theta X_s \right) ds + \sigma_\theta \partial_1 K_H(t, \theta), \quad 0 \le \theta < t.$$
- Đánh giá giải tích: Tồn tại hằng số $c, C > 0$ sao cho $|D_\theta X_t| \le C K_H(t, \theta)$ và khi $\inf_{t \in [0, T]} |\sigma(t)| > 0$: $$c t^{2H} \le \int_0^t D_\theta X_t \cdot \mathbb{E}[D_\theta X_t | \mathcal{F}_\theta] d\theta \le C t^{2H}, \quad \text{h.c.}$$
- Định lý 3.11 (Ước lượng Gauss cho nhiễu cộng tính): Mật độ xác suất $p_{X_t}(x)$ của $X_t$ tồn tại và thỏa mãn: $$\frac{\mathbb{E}|X_t - \mathbb{E} X_t|}{2 C t^{2H}} \exp\left( -\frac{(x - \mathbb{E} X_t)^2}{2 c t^{2H}} \right) \le p_{X_t}(x) \le \frac{\mathbb{E}|X_t - \mathbb{E} X_t|}{2 c t^{2H}} \exp\left( -\frac{(x - \mathbb{E} X_t)^2}{2 C t^{2H}} \right), \quad \forall x \in \mathbb{R}.$$
2. Trường hợp nhiễu nhân tính ($\sigma = \sigma(t, x)$)
- Giả thiết $(A_4)$: $a, g, \sigma \in C_b^{0, 1}([0, T] \times \mathbb{R})$ và tồn tại $m > 0$ sao cho $|\sigma(t, x)| \ge m > 0$ với mọi $(t, x) \in [0, T] \times \mathbb{R}$.
- Phương pháp đổi biến phi tuyến: Đặt $F(t, x) = \int_0^x \frac{1}{\sigma(t, u)} du$. Hàm ngược $G(t, y)$ xác định trên $[0, T] \times \mathbb{R}$ thỏa mãn $F(t, G(t, y)) = y$. Xét phương trình phụ cho $Y_t$ có dạng nhiễu cộng tính: $$Y_t = y_0 + \int_0^t \left( A(s, Y_s) + \frac{1}{\sigma(s, G(s, Y_s))} \int_{s-r}^s \rho(u - s) g(u, G(u, Y_u)) du \right) ds + B_t^H,$$ trong đó $A(s, y) = \frac{\partial F}{\partial s}(s, G(s, y)) + \frac{a(s, G(s, y))}{\sigma(s, G(s, y))}$.
- Mệnh đề 3.12 (Sự tồn tại và duy nhất nghiệm mạnh): Phương trình có nghiệm mạnh duy nhất $X_t = G(t, Y_t)$, thích ứng với $\mathcal{F}_t$, và các quỹ đạo của $X_t$ liên tục Hölder bậc $(H - \varepsilon)$ với mọi $\varepsilon \in (0, H)$.
- Mệnh đề 3.13 & 3.14 (Tính khả vi Malliavin): $X_t \in \mathbb{D}^{1,2}$ với $D_\theta X_t = \sigma(t, X_t) D_\theta Y_t$ và thỏa mãn đánh giá $|D_\theta X_t| \le C K_H(t, \theta)$ h.c.
- Định lý 3.16 (Ước lượng Gauss cho nhiễu nhân tính): Mật độ $p_{X_t}(x)$ tồn tại và thỏa mãn chặn hai phía dạng Gauss: $$\frac{\mathbb{E}|X_t - \mathbb{E} X_t|}{2 C t^{2H}} \exp\left( -\frac{(x - \mathbb{E} X_t)^2}{2 c t^{2H}} \right) \le p_{X_t}(x) \le \frac{\mathbb{E}|X_t - \mathbb{E} X_t|}{2 c t^{2H}} \exp\left( -\frac{(x - \mathbb{E} X_t)^2}{2 C t^{2H}} \right), \quad \forall x \in \mathbb{R}.$$
- Định lý 3.18 (Tính trơn vô hạn cấp của mật độ): Giả sử các hệ số $a, g, \sigma$ khả vi vô hạn lần theo $x$ với các đạo hàm mọi cấp bị chặn. Khi đó:
- $X_t \in \mathbb{D}^\infty = \bigcap_{k \ge 1, p > 1} \mathbb{D}^{k, p}$.
- Đại lượng $|DX_t|{L^2[0, T]}^{-1} \in \bigcap{\alpha > 1} L^\alpha(\Omega)$ thông qua đánh giá: $$\mathbb{E}\left[ \left( \int_0^t |D_\theta X_t|^2 d\theta \right)^{-\alpha} \right] < \infty, \quad \forall \alpha > 1.$$
- Mật độ xác suất $p_{X_t}(x)$ khả vi vô hạn lần trên $\mathbb{R}$ ($p_{X_t} \in C^\infty(\mathbb{R})$).
Kết quả và những đóng góp mới
Tổng hợp các kết quả khoa học của luận án
Luận án đã đạt được các đóng góp mới về mặt lý thuyết giải tích ngẫu nhiên, được tóm tắt trong bảng sau:
| STT | Đối tượng nghiên cứu | Phương pháp / Kỹ thuật | Kết quả đạt được | Cấp độ đánh giá |
|---|---|---|---|---|
| 1 | Xấp xỉ Smoluchowski-Kramers ($\mu \to 0$) | Giải tích Malliavin, chặn moment sai phân, đạo hàm Malliavin cấp 1 và cấp 2 | Ước lượng Berry-Esseen hiển theo khoảng cách Kolmogorov | $\sup_{x \in \mathbb{R}} |P(X_t^\mu \le x) - P(X_t \le x)| \le C t^{-1/2} \sqrt{\mu}$ |
| 2 | SDE có trễ tham số $\tau_1, \tau_2$ | Toán tử Malliavin đối ngẫu trên không gian hàm đo được bị chặn $\mathcal{G}$ | Đánh giá tốc độ hội tụ yếu theo khoảng cách tham số trễ | Chặn cấp $|\tau_1 - \tau_2|^{1/2}$ (hoặc cấp $\frac{\beta}{4} \wedge \frac{1}{2}$ với dữ liệu $\beta$-Hölder) |
| 3 | Lược đồ giải số Carathéodory $x^n(t)$ | Áp dụng kết quả SDE có trễ với $\tau_1 = 0, \tau_2 = 1/n$ | Tốc độ hội tụ yếu của lược đồ Carathéodory trên hàm đo được bị chặn | $\sup_{0 \le t \le T} |\mathbb{E} g(x^n(t)) - \mathbb{E} g(x(t))| \le \frac{C}{\sqrt{n}}$ |
| 4 | SFDE phân thứ ($H > 1/2$) nhiễu cộng tính | Nhân Volterra $K_H(t, s)$, tiêu chuẩn mật độ Nourdin-Viens | Sự tồn tại và ước lượng Gauss cận trên, cận dưới hai phía cho mật độ | $p_{X_t}(x) \asymp \frac{\mathbb{E}|X_t - \mathbb{E}X_t|}{t^{2H}} \exp\left( -\frac{(x - \mathbb{E}X_t)^2}{\text{const} \cdot t^{2H}} \right)$ |
| 5 | SFDE phân thứ ($H > 1/2$) nhiễu nhân tính | Đổi biến phi tuyến qua hàm ẩn $F(t, x)$, ước lượng giải tích | Nghiệm mạnh duy nhất, liên tục Hölder $(H-\varepsilon)$, ước lượng Gauss hai phía | $p_{X_t}(x)$ tồn tại và thỏa mãn chặn Gauss tỷ lệ theo bậc $t^{2H}$ |
| 6 | Tính trơn nghiệm SFDE phân thứ | Chặn moment nghịch đảo đạo hàm Malliavin cấp cao trong $\mathbb{D}^\infty$ | Chứng minh mật độ nghiệm là hàm trơn vô hạn lần | $p_{X_t} \in C^\infty(\mathbb{R})$ |
Các công trình khoa học đã công bố của tác giả
Các kết quả chính của luận án đã được công bố trong hai công trình khoa học:
- [CT1] Nguyen Van Tan (2020), "Smoothness and Gaussian Density Estimates for Stochastic Functional Differential Equations With Fractional Noise", Statistics, Optimization & Information Computing, 8(4), pp. 883–900.
- [CT2] Nguyen Van Tan, Nguyen Tien Dung (2020), "A Berry-Esseen Bound in the Smoluchowski-Kramers Approximation", Journal of Statistical Physics, 179, pp. 671–685.
Hạn chế và hướng nghiên cứu tiếp
Hạn chế được ghi nhận trong văn bản
- Các kết quả về phương trình vi phân ngẫu nhiên có trễ trong Chương 2 mới khảo sát cho trường hợp một độ trễ cố định $\tau > 0$, chưa bao quát các hệ phương trình có nhiều độ trễ (đa trễ) hoặc độ trễ biến đổi theo thời gian $\tau(t)$.
- Kỹ thuật giải tích hàm phân thứ trong Chương 3 chỉ áp dụng cho tham số Hurst $H \in (1/2, 1)$, nơi tích phân đường xác định theo nghĩa Riemann-Stieltjes quỹ đạo. Trường hợp $H \in (0, 1/2)$ có cấu trúc số học phân kỳ phức tạp hơn chưa được khảo sát trong nghiên cứu này.
- Các điều kiện đặt lên hệ số khuếch tán $\sigma$ đều yêu cầu tính elliptic ngặt ($\inf |\sigma| > 0$) và các đạo hàm riêng bị chặn đều, chưa xét đến trường hợp suy biến tại biên hoặc hệ số tăng trưởng siêu tuyến tính.
Hướng nghiên cứu tiếp
- Mở rộng bài toán khảo sát sự hội tụ yếu trong lược đồ xấp xỉ Carathéodory cho các hệ phương trình vi phân ngẫu nhiên tổng quát hơn có cấu trúc đa trễ hoặc biến trễ.
- Nghiên cứu tính chất mật độ và tốc độ hội tụ cho các hệ phương trình vi phân ngẫu nhiên có bước nhảy (quá trình Lévy) hoặc phương trình vi phân ngẫu nhiên phân thứ có phản xạ tại biên.
Giá trị tham khảo
Luận án có giá trị tham khảo chuyên môn trong các lĩnh vực giải tích ngẫu nhiên và toán ứng dụng:
- Đối với nghiên cứu sinh và học viên cao học ngành Toán giải tích, Xác suất - Thống kê: Luận án cung cấp một tài liệu tham khảo có cấu trúc hệ thống về việc áp dụng giải tích Malliavin để xử lý các bài toán biên phân bố xác suất, kỹ thuật ước lượng moment nghịch đảo của ma trận/toán tử Malliavin, và phương pháp đổi biến phi tuyến giải tích trên không gian Wiener.
- Đối với giảng viên và nhà nghiên cứu lý thuyết xác suất: Cung cấp công thức hiển về chặn Berry-Esseen trong động học Smoluchowski-Kramers và khung lý thuyết đánh giá sai số yếu của các phương pháp giải số trên lớp hàm đo được tổng quát.
- Đối với các chuyên gia trong lĩnh vực Toán tài chính và Vật lý thống kê: Các mô hình hạt có quán tính nhỏ (Smoluchowski-Kramers) và phương trình trễ nhớ dài (fBm) cung cấp cơ sở giải tích cho các bài toán định giá tài sản ngẫu nhiên có trễ, lọc ngẫu nhiên và mô hình hóa hệ thống vật lý có tương tác trễ.
Câu hỏi thường gặp
1. Luận án đạt được tốc độ hội tụ Berry-Esseen cụ thể là bao nhiêu trong xấp xỉ Smoluchowski-Kramers?
Tốc độ hội tụ theo khoảng cách Kolmogorov được chứng minh là cấp $\mathcal{O}(\mu^{1/2})$ với chặn hiển $\sup_{x \in \mathbb{R}} |P(X_t^\mu \le x) - P(X_t \le x)| \le C t^{-1/2} \sqrt{\mu}$, trong đó $\mu$ là khối lượng hạt và $C$ là hằng số dương độc lập với $t$ và $\mu$.
2. Điểm mới căn bản trong kết quả về xấp xỉ Carathéodory của luận án là gì?
Điểm mới căn bản là ước lượng tốc độ hội tụ yếu cấp $\mathcal{O}(n^{-1/2})$ được thiết lập cho toàn bộ không gian hàm đo được bị chặn bất kỳ $\mathcal{G} = {g: \mathbb{R} \to \mathbb{R}, |g|_\infty \le 1}$, thay vì chỉ áp dụng cho lớp các hàm thử thỏa mãn tính liên tục Lipschitz như các công trình trước đây.
3. Phương pháp nào được sử dụng để xử lý nhiễu nhân tính trong phương trình vi phân hàm phân thứ?
Tác giả sử dụng phép biến đổi hàm thế $F(t, x) = \int_0^x \frac{1}{\sigma(t, u)} du$ và hàm ngược giải tích $G(t, y)$ để chuyển phương trình ban đầu có nhiễu nhân tính về phương trình vi phân hàm ngẫu nhiên phụ $Y_t$ có cấu trúc nhiễu cộng tính với hệ số trôi mới vẫn bảo toàn tính liên tục Lipschitz.
4. Điều kiện nào bảo đảm mật độ xác suất của nghiệm phương trình vi phân hàm phân thứ là khả vi vô hạn lần ($C^\infty$)?
Mật độ xác suất $p_{X_t}(x) \in C^\infty(\mathbb{R})$ khi các hệ số $a, g, \sigma$ khả vi vô hạn lần theo biến không gian với các đạo hàm mọi cấp bị chặn đều, hệ số $\sigma$ thỏa mãn điều kiện elliptic đều $\inf |\sigma| \ge m > 0$, và tham số Hurst $H \in (1/2, 1)$.
5. Vai trò của toán tử Skorohod và hệ thức đối ngẫu trong luận án là gì?
Toán tử Skorohod $\delta$ đóng vai trò là toán tử đối ngẫu liên hợp của đạo hàm Malliavin $D$. Hệ thức đối ngẫu $\mathbb{E}[F \delta(u)] = \mathbb{E}\left[\int_0^T u(t) D_t F dt\right]$ cho phép biến đổi sai phân tích phân của các hàm chỉ tiêu không trơn thành các tích phân kỳ vọng chứa đạo hàm Malliavin cấp 1 và cấp 2, từ đó thu được các đánh giá định lượng cho khoảng cách Kolmogorov và sai số hội tụ yếu.
Kết luận
Luận án của nghiên cứu sinh Nguyễn Văn Tân là một công trình nghiên cứu toán học chuyên sâu thuộc chuyên ngành Lý thuyết Xác suất và Thống kê toán học, ứng dụng giải tích Malliavin vào việc khảo sát phương trình vi phân ngẫu nhiên. Công trình đã giải quyết trọn vẹn ba bài toán lý thuyết: thiết lập ước lượng Berry-Esseen hiển cho xấp xỉ Smoluchowski-Kramers, đánh giá tốc độ hội tụ yếu của lược đồ xấp xỉ Carathéodory trên không gian hàm đo được bị chặn, và chứng minh tính trơn cùng ước lượng cận Gauss hai phía cho nghiệm phương trình vi phân hàm ngẫu nhiên phân thứ ($H > 1/2$). Các kết quả đạt được có tính chính xác cao về mặt giải tích ngẫu nhiên, đã được thẩm định và công bố trên các tạp chí khoa học chuyên ngành quốc tế.