Tổng quan luận án

Phương trình vi phân ngẫu nhiên được đề xuất bởi Kiyosi Itô và được ứng dụng trong toán tài chính, vật lý, sinh học, lý thuyết lọc và tối ưu hóa điều khiển. Trong các mô hình thực tế, chẳng hạn như lý thuyết định giá cổ phiếu của Black – Scholes và Merton hay mô hình độ biến động ngẫu nhiên, các hệ số thường không thỏa mãn điều kiện Lipschitz toàn cục hoặc xuất hiện kỳ dị tại biên. Khi tính toán các đại lượng dạng kỳ vọng $\mathbb{E}[f(X)]$, việc tìm nghiệm giải tích chính xác chỉ thực hiện được trong số ít trường hợp, đòi hỏi phải áp dụng các phương pháp số và kỹ thuật Monte Carlo.

Phương pháp Monte Carlo đa cấp do Giles đề xuất giúp giảm khối lượng tính toán, nhưng đòi hỏi nghiệm xấp xỉ phải đạt tốc độ hội tụ mạnh trong không gian $L^p$. Tuy nhiên, lược đồ Euler-Maruyama cổ điển (PEM) đối với phương trình có hệ số tăng trên tuyến tính thường bị phân kỳ mô men, hoặc không bảo toàn được các tính chất định tính quan trọng như tính ổn định tiệm cận, tính không âm và tính không va chạm của nghiệm đúng. Các nghiên cứu giải số cho phương trình có hệ số không chính qui (hệ số tăng trên tuyến tính, chỉ thỏa mãn điều kiện liên tục Hölder hoặc nổ tại biên) vẫn còn hạn chế.

Mục tiêu nghiên cứu của luận án bao gồm:

  1. Thiết lập các định lý về sự tồn tại và duy nhất nghiệm cho các lớp phương trình vi phân ngẫu nhiên có hệ số không chính qui.
  2. Xây dựng và đánh giá tốc độ hội tụ mạnh của các lược đồ xấp xỉ Euler-Maruyama khống chế cho phương trình có hệ số dịch chuyển tăng trên tuyến tính và hệ số khuếch tán liên tục Hölder hoặc Hölder địa phương.
  3. Xây dựng lược đồ Euler-Maruyama cải tiến đảm bảo sự hội tụ, tính không âm và tính ổn định mô men trong không gian $L^p$ cho phương trình có hệ số dịch chuyển Lipschitz một phía với hệ số âm.
  4. Xây dựng lược đồ Milstein nửa ẩn cho hệ phương trình vi phân ngẫu nhiên dạng Dyson mô tả hệ điểm không va chạm và đánh giá tốc độ hội tụ.

Đối tượng nghiên cứu là phương trình vi phân ngẫu nhiên có dạng: $$dX_t = b(X_t)dt + \sigma(X_t)dW_t, \quad X_0 = x_0$$ với hệ số dịch chuyển $b(x)$ và hệ số khuếch tán $\sigma(x)$ thuộc các trường hợp:

  • $b(x)$ liên tục Lipschitz địa phương; $\sigma(x)$ liên tục $(1/2 + \alpha)$-Hölder địa phương.
  • $b(x)$ tăng trên tuyến tính, Lipschitz một phía; $\sigma(x)$ liên tục $(1/2 + \alpha)$-Hölder hoặc $(1/2 + \alpha)$-Hölder địa phương.
  • $b(x)$ Lipschitz và Lipschitz một phía với hệ số âm; $\sigma(x)$ liên tục $(1/2 + \alpha)$-Hölder địa phương.
  • $b(x)$ không liên tục hoặc có điểm nổ tại biên (như hệ Dyson); $\sigma(x)$ khả vi liên tục đến cấp hai.

Phạm vi nghiên cứu thuộc sự giao thoa giữa Giải tích ngẫu nhiên, Phương trình vi phân và Giải tích số, tập trung vào việc phân tích tính chất định lượng và định tính của nghiệm đúng và nghiệm xấp xỉ rời rạc trên trục thời gian thực.

Tổng quan tài liệu và vị trí của luận án

Nghiên cứu về sự tồn tại và duy nhất nghiệm của phương trình vi phân ngẫu nhiên bắt đầu từ phương pháp xấp xỉ Picard và nguyên lý ánh xạ co trên không gian $L^2$ của Itô khi các hệ số thỏa mãn điều kiện Lipschitz toàn cục. Đối với các hệ số không Lipschitz, Stroock và Varadhan đã giải quyết bài toán martingale với hệ số dịch chuyển đo được, bị chặn và hệ số khuếch tán liên tục, elliptic đều. Krylov mở rộng cho trường hợp hệ số khuếch tán đo được, và Veretennikov chứng minh sự tồn tại, duy nhất nghiệm cho hệ nhiều chiều với hệ số dịch chuyển bị chặn và hệ số khuếch tán Lipschitz, elliptic đều. Engelbert và Schmidt đưa ra điều kiện cần và đủ cho phương trình một chiều qua phương pháp đổi thời gian. Yamada và Watanabe chứng minh mối liên hệ giữa tính duy nhất theo quỹ đạo và tính duy nhất theo phân phối, đồng thời thiết lập định lý tồn tại nghiệm mạnh duy nhất khi hệ số dịch chuyển Lipschitz và hệ số khuếch tán liên tục $1/2$-Hölder. Girsanov chỉ ra phản ví dụ khi hệ số khuếch tán liên tục $\alpha$-Hölder với $\alpha < 1/2$, và Tanaka đưa ra ví dụ về phương trình có hệ số khuếch tán là hàm dấu không có nghiệm mạnh.

Về phương diện giải số, Maruyama đề xuất lược đồ Euler-Maruyama đạt tốc độ hội tụ mạnh $1/2$ trong $L^2$ với hệ số Lipschitz toàn cục. Milstein thiết lập lược đồ bậc cao đạt tốc độ hội tụ $1$ trong $L^2$ cho phương trình có hệ số khả vi liên tục cấp hai. Khi các hệ số không chính qui, việc giải số gặp nhiều thách thức:

  • Hệ số không trơn: Hairer, Hutzenthaler và Jentzen chỉ ra ví dụ với $b$ bị chặn, trơn nhưng đạo hàm biến động cao khiến Euler-Maruyama không hội tụ với bất kỳ tốc độ đa thức nào. Hefter và Jentzen chứng minh tốc độ hội tụ của mọi phương pháp rời rạc đều rất chậm đối với quá trình Cox–Ingersoll–Ross khi tham số $\sigma$ lớn. Gyöngy và Rásonyi sử dụng kỹ thuật Yamada-Watanabe chứng minh lược đồ Euler đạt tốc độ hội tụ $\alpha$ trong $L^1$ khi hệ số khuếch tán liên tục $(1/2+\alpha)$-Hölder và hệ số dịch chuyển Lipschitz.
  • Hệ số tăng trên tuyến tính: Hutzenthaler, Jentzen và Kloeden chứng minh mô men của nghiệm xấp xỉ Euler-Maruyama cổ điển bị phân kỳ. Để khắc phục, Hu cùng Higham, Mao và Stuart đề xuất lược đồ Euler-Maruyama ẩn. Sau đó, Hutzenthaler, Jentzen, Kloeden và Sabanis phát triển lược đồ Euler-Maruyama khống chế (tamed Euler-Maruyama - TEM) dạng hiển; Wang và Gan đề xuất Milstein khống chế; Mao thiết lập Euler-Maruyama chặt.
  • Hệ số nổ tại biên: Điển hình là quá trình Bessel và chuyển động Dyson-Brown. Bru, Cépa và Lépingle, Graczyk và Malecki, Rogers và Shi, Naganuma và Taguchi nghiên cứu tính chất nghiệm hệ Dyson. Li và Menon đề xuất lược đồ Euler-Maruyama khống chế nhưng không giữ được nghiệm trong miền xác định $\Delta_d$. Jentzen, Kloeden và Neuenkirch xây dựng lược đồ chiếu nhưng không áp dụng được khi hệ số nổ trên biên. Gần đây, Ngô Hoàng Long và Taguchi đưa ra lược đồ Euler-Maruyama nửa ẩn (SIEM) bảo toàn miền giá trị $\Delta_d$ với tốc độ hội tụ $1/2$.

Về tính ổn định tiệm cận, Higham, Mao và Yuan chỉ ra Euler-Maruyama cổ điển có thể mất tính ổn định mô men của nghiệm đúng. Gou, Mao và Yue, Liu và Mao, Szpruch và Zhang đã nghiên cứu tính $V$-khả tích và ổn định tiệm cận cho các biến thể Euler-Maruyama nhưng chủ yếu áp dụng cho hệ số Lipschitz địa phương.

Luận án lựa chọn giải quyết các khoảng trống: xây dựng và chứng minh tốc độ hội tụ của lược đồ khống chế cho phương trình có hệ số dịch chuyển tăng trên tuyến tính kết hợp hệ số khuếch tán liên tục Hölder; thiết lập lược đồ cải tiến bảo toàn tính không âm và tính ổn định mô men; xây dựng lược đồ bậc cao Milstein nửa ẩn (SIM) cho hệ hạt không va chạm dạng Dyson.

Cơ sở lý thuyết và phương pháp nghiên cứu

Luận án sử dụng các cơ sở lý thuyết và công cụ giải tích ngẫu nhiên sau:

  • Không gian xác suất đầy đủ $(\Omega, \mathcal{F}, \mathbb{P})$ với họ lọc ${\mathcal{F}t}{t \ge 0}$ thỏa mãn các điều kiện thông thường; chuyển động Brown $m$-chiều $W_t = (W_t^1, \dots, W_t^m)^\top$.
  • Tích phân ngẫu nhiên Itô, công thức Itô cho các hàm thuộc không gian $C^{1,2}(\mathbb{R}_+ \times \mathbb{R}; \mathbb{R})$, tính chất martingale bình phương khả tích.
  • Bất đẳng thức Burkholder-Davis-Gundy (bất đẳng thức BDG), bất đẳng thức Gronwall, bất đẳng thức Young.
  • Kỹ thuật xấp xỉ hàm trị tuyệt đối Yamada-Watanabe: Sử dụng họ hàm trơn $\phi_{\delta\varepsilon}(x) = \int_0^{|x|} \int_0^y \psi_{\delta\varepsilon}(z)dzdy$ với hàm $\psi_{\delta\varepsilon}(z)$ liên tục, có giá $\text{supp},\psi_{\delta\varepsilon} \subset [\varepsilon/\delta, \varepsilon]$, thỏa mãn $\int_{\varepsilon/\delta}^\varepsilon \psi_{\delta\varepsilon}(z)dz = 1$ và $0 \le \psi_{\delta\varepsilon}(z) \le \frac{2}{z \log \delta}$. Hàm $\phi_{\delta\varepsilon}$ thỏa mãn các tính chất:
    • $(\text{YW}1)\ \phi'{\delta\varepsilon}(x) = \frac{x}{|x|} \phi'{\delta\varepsilon}(|x|)$
    • $(\text{YW}2)\ 0 \le |\phi'_{\delta\varepsilon}(x)| \le 1$
    • $(\text{YW}3)\ \phi_{\delta\varepsilon}(x) \le |x| \le \varepsilon + \phi_{\delta\varepsilon}(x)$
    • $(\text{YW}4)\ \phi''{\delta\varepsilon}(|x|) = \psi{\delta\varepsilon}(|x|) \le \frac{2\delta}{\varepsilon \log \delta}$
    • $(\text{YW}5)\ \phi''{\delta\varepsilon}(|x|) = \psi{\delta\varepsilon}(|x|) \le \frac{2}{|x|\log\delta} \mathbb{I}_{[\varepsilon/\delta, \varepsilon]}(|x|)$
  • Phương pháp Monte Carlo đa cấp (MLMC) của Giles kết hợp với định lý cơ bản về sai số một bước và sai số toàn cục.
  • Kỹ thuật thời điểm dừng $\theta_R$ để địa phương hóa và nguyên lý đánh giá số mũ Lyapunov nhằm phân tích tính ổn định tiệm cận.

Phương pháp nghiên cứu kết hợp giữa chứng minh giải tích giải tích ngẫu nhiên chặt chẽ và thực nghiệm mô phỏng số trên máy tính nhằm kiểm chứng tốc độ hội tụ của các thuật toán đề xuất (PEM, BEM, TEM, SIEM, SIM).

Nội dung chính theo từng chương

Chương 1: Sơ lược về giải số phương trình vi phân ngẫu nhiên

Chương 1 hệ thống hóa các khái niệm và kết quả nền tảng về phương trình vi phân ngẫu nhiên Itô:

  • Định nghĩa nghiệm mạnh, điều kiện tồn tại duy nhất nghiệm theo tiêu chuẩn Lipschitz và tăng tuyến tính toàn cục, tính chất dòng của nghiệm $X(t; t_0, x_0) = X(t; s, X(s; t_0, x_0))$.
  • Tính bị chặn mô men bậc $p \ge 2$: $\mathbb{E}[|X_t - X_s|^p] \le C(t-s)^{p/2}$, tính liên tục theo thời gian và tính ổn định tiệm cận qua số mũ Lyapunov $\limsup_{t\to\infty} \frac{1}{t} \log |X_t|$.
  • Cơ sở phương pháp Monte Carlo đa cấp: ước lượng $\mathbb{E}[f(X)]$ qua tổng các cấp $\sum_{j=1}^L \frac{1}{n_j} \sum_{i=1}^{n_j} (S_j - S_{j-1})_i$, phân tích hàm chi phí tối ưu phụ thuộc vào tốc độ hội tụ mạnh $r_s$ và hội tụ yếu $r_w$.
  • Định lý cơ bản về sự hội tụ theo trung bình: nếu sai số một bước của kỳ vọng đạt cỡ $h^{p_1}$ và sai số trung bình bậc hai đạt cỡ $h^{p_2}$ thì tốc độ hội tụ mạnh toàn cục là $r_s = p_2 - 1/2$.
  • Cấu trúc và tốc độ hội tụ của lược đồ Euler ($r_s = 1/2$ cho nhiễu nhân, $r_s = 1$ cho nhiễu cộng tính) và lược đồ Milstein ($r_s = 1$).
  • Tổng thuật các kết quả về phương trình với hệ số không chính qui: lớp hệ số không đủ trơn (Hölder), hệ số tăng trên tuyến tính, hệ số suy biến hoặc nổ tại biên và các nghiên cứu về tính ổn định của nghiệm rời rạc.

Chương 2: Lược đồ Euler-Maruyama khống chế cho phương trình vi phân ngẫu nhiên

Chương 2 tập trung nghiên cứu lược đồ Euler-Maruyama khống chế (TEM) cho phương trình vi phân ngẫu nhiên một chiều: $$X_t = x_0 + \int_0^t b(s, X_s)ds + \int_0^t \sigma(s, X_s)dW_s$$ Tác giả thiết lập hệ điều kiện giả thiết:

  • (A1) $x b(t,x) \vee |\sigma(t,x)|^2 \le L(1+|x|^2)$; hoặc (A1') $x b(s,x) + \frac{p_0-1}{2}|\sigma(s,x)|^2 \le L(1+|x|^2)$.
  • (A2) $b(t,x)$ thỏa mãn điều kiện Lipschitz một phía: $(x-y)(b(t,x)-b(t,y)) \le L|x-y|^2$.
  • (A3) $b(t,x)$ thỏa mãn điều kiện Lipschitz địa phương tăng đa thức: $|b(t,x)-b(t,y)| \le L(1+|x|^l+|y|^l)|x-y|$ và $|b(t,x)| \le L(1+|x|^{l+1})$; hoặc (A3') liên tục Lipschitz địa phương và bị chặn địa phương bởi $L_R$.
  • (A4) $\sigma(t,x)$ liên tục $(1/2+\alpha)$-Hölder: $|\sigma(t,x)-\sigma(t,y)| \le L|x-y|^{1/2+\alpha}$ với $\alpha \in (0, 1/2]$; hoặc (A4') Hölder địa phương tăng đa thức; hoặc (A4'') Hölder địa phương bị chặn trên tập compact bởi $L_R$.
  • (A5) Điều kiện liên tục Hölder theo biến thời gian với số mũ $\beta \in (0, 1]$.

Dưới các điều kiện (A1'), (A3') và (A4''), luận án chứng minh định lý về sự tồn tại và duy nhất nghiệm mạnh của phương trình thông qua việc xây dựng dãy cắt $(b_N, \sigma_N)$ kết hợp hàm xấp xỉ Yamada-Watanabe $\phi_{\delta\varepsilon}$ và bổ đề Borel-Cantelli.

Luận án xây dựng lược đồ Euler-Maruyama khống chế liên tục: $$X_t^{(n)} = x_0 + \int_0^t b_n(s, X_{\eta_n(s)}^{(n)})ds + \int_0^t \sigma(s, X_{\eta_n(s)}^{(n)})dW_s$$ trong đó $b_n(t,x) = \frac{b(t,x)}{1 + n^{-\lambda}|b(t,x)|}$ với tham số $\lambda \in [\alpha, 1/2]$, $\eta_n(t) = t_k = \frac{kT}{n}$ khi $t \in [t_k, t_{k+1})$.

Kết quả đánh giá tốc độ hội tụ mạnh của lược đồ khống chế: $$\sup_{\tau \in \mathcal{T}} \mathbb{E}\left[|X_\tau - X_\tau^{(n)}|\right] \le \begin{cases} \frac{C}{n^\alpha} & \text{nếu } \alpha \in (0, 1/2] \ \frac{C}{\log n} & \text{nếu } \alpha = 0 \end{cases}$$ với $\mathcal{T}$ là tập các thời điểm dừng $\tau \le T$.

Bảng tổng hợp các điều kiện và tốc độ hội tụ tương ứng trong Chương 2:

Nhóm điều kiện hệ số dịch chuyển $b$ Nhóm điều kiện hệ số khuếch tán $\sigma$ Tốc độ hội tụ mạnh trong $L^1$
Tăng trên tuyến tính, Lipschitz một phía (A1, A2, A3) Liên tục $(1/2+\alpha)$-Hölder toàn cục (A4) Cỡ $n^{-\alpha}$ với $\alpha \in (0, 1/2]$
Tăng trên tuyến tính, Lipschitz một phía (A1, A2, A3) Liên tục $1/2$-Hölder toàn cục ($\alpha = 0$) Cỡ $(\log n)^{-1}$
Tăng đa thức, Lipschitz địa phương (A1, A2, A3) Liên tục $(1/2+\alpha)$-Hölder địa phương (A4') Xác định theo bậc Hölder $\alpha$

Chương 3: Sự hội tụ, tính không âm và ổn định của lược đồ Euler-Maruyama cải tiến

Chương 3 khảo sát lớp phương trình vi phân ngẫu nhiên có hệ số dịch chuyển Lipschitz một phía với hằng số âm và hệ số khuếch tán liên tục Hölder địa phương:

  • Mở rộng kỹ thuật xấp xỉ Yamada-Watanabe để xử lý tính chất không âm của nghiệm.
  • Xây dựng lược đồ Euler-Maruyama cải tiến kết hợp giữa kỹ thuật khống chế và phép chiếu/chặn bảo toàn miền giá trị.
  • Chứng minh nghiệm xấp xỉ bảo toàn tính không âm của nghiệm đúng hầu chắc chắn: nếu $X_0 \ge 0$ thì $X_k \ge 0$ với mọi bước rời rạc $k$.
  • Đánh giá tốc độ hội tụ mạnh của lược đồ cải tiến, đạt mức tương đương với trường hợp hệ số chính qui.
  • Chứng minh tính ổn định mũ của mô men trong không gian $L^p$ ($\mathbb{E}[|X_k|^p] \le C e^{-\gamma t_k}$) khi bước nhảy thời gian $h$ cố định.
  • Thực nghiệm số so sánh giữa lược đồ Euler-Maruyama gốc (PEM), lược đồ Euler-Maruyama nghịch (BEM) và lược đồ cải tiến không âm khống chế (TEM), minh họa qua đồ thị sai số theo thang đo $\log_2 - \log_2$ và mô phỏng quỹ đạo nghiệm.

Chương 4: Lược đồ Milstein nửa ẩn cho hệ điểm không va chạm

Chương 4 nghiên cứu hệ phương trình vi phân ngẫu nhiên nhiều chiều dạng Dyson mô tả hệ $d$ hạt tương tác không va chạm trên miền $\Delta_d = {x \in \mathbb{R}^d: x_1 < x_2 < \dots < x_d}$: $$X_i(t) = X_i(0) + \int_0^t \left[\sum_{j \neq i} \frac{\gamma_{ij}}{X_i(s) - X_j(s)} + b_i(X_i(s))\right]ds + \int_0^t \sigma_i(X_i(s))dW_i(s), \quad 1 \le i \le d$$

  • Phân tích tính chất kỳ dị nổ của hệ số dịch chuyển khi các hạt tiến lại gần nhau ($X_i(s) - X_j(s) \to 0$).
  • Thiết lập lược đồ Euler-Maruyama nửa ẩn (SIEM) và phát triển lên lược đồ xấp xỉ Milstein nửa ẩn (SIM) nhằm tăng bậc chính xác.
  • Thiết lập biểu diễn giải tích của sai số rời rạc và các đánh giá chặn mô men cho thành phần nửa ẩn.
  • Chứng minh nghiệm xấp xỉ của lược đồ Milstein nửa ẩn không bao giờ va chạm (bảo toàn miền $\Delta_d$) và đánh giá tốc độ hội tụ mạnh của lược đồ.
  • Trình bày ví dụ mô phỏng số so sánh sai số trung bình bình phương $mse_E(k)$ của lược đồ Euler nửa ẩn và $mse_M(k)$ của lược đồ Milstein nửa ẩn theo thang đo $\log_2$.

Kết quả và những đóng góp mới

Các kết quả và đóng góp khoa học mới của luận án gồm:

  1. Chứng minh định lý tồn tại và duy nhất nghiệm mạnh cho phương trình vi phân ngẫu nhiên có hệ số dịch chuyển liên tục Lipschitz địa phương tăng đa thức và hệ số khuếch tán liên tục $(1/2+\alpha)$-Hölder địa phương.
  2. Xây dựng lược đồ Euler-Maruyama khống chế (TEM) dạng hiển cho phương trình có hệ số dịch chuyển tăng trên tuyến tính và hệ số khuếch tán liên tục Hölder; xác định chính xác tốc độ hội tụ mạnh của lược đồ đạt cỡ $n^{-\alpha}$ khi $\alpha \in (0, 1/2]$ và cỡ $(\log n)^{-1}$ khi $\alpha = 0$.
  3. Đề xuất lược đồ Euler-Maruyama cải tiến bảo toàn tính không âm của nghiệm, đạt tốc độ hội tụ phù hợp và chứng minh tính ổn định mũ theo chuẩn trong không gian $L^p$ với bước nhảy thời gian cố định.
  4. Xây dựng lược đồ Milstein nửa ẩn (SIM) cho hệ điểm không va chạm dạng Dyson, chứng minh tính bảo toàn miền hình học $\Delta_d$ của nghiệm xấp xỉ và thiết lập đánh giá tốc độ hội tụ mạnh bậc cao cho hệ phương trình có hệ số nổ tại biên.

Hạn chế và hướng nghiên cứu tiếp

Hạn chế và phạm vi mở rộng của luận án:

  • Các đánh giá tốc độ hội tụ của lược đồ Euler-Maruyama khống chế trong Chương 2 và Chương 3 chủ yếu được thiết lập chi tiết cho không gian một chiều $\mathbb{R}$. Việc mở rộng các kỹ thuật xấp xỉ Yamada-Watanabe đa chiều cho phương trình có hệ số khuếch tán Hölder tổng quát trên $\mathbb{R}^d$ vẫn là bài toán phức tạp.
  • Đối với hệ điểm không va chạm dạng Dyson trong Chương 4, lược đồ Milstein nửa ẩn đòi hỏi hệ số khuếch tán $\sigma$ phải thỏa mãn điều kiện trơn (khả vi liên tục đến cấp hai). Hướng nghiên cứu tiếp theo là khảo sát các lược đồ số bậc cao khi hệ số khuếch tán $\sigma$ không trơn hoặc có cấu trúc ma trận tương tác phức tạp hơn.

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 Lý thuyết xác suất và Thống kê toán học, Toán giải tích và Giải tích số.
  • Các nhà nghiên cứu trong lĩnh vực mô hình hóa toán tài chính (mô hình biến động ngẫu nhiên, mô hình Cox-Ingersoll-Ross), vật lý thống kê (hệ ma trận ngẫu nhiên, hệ hạt tương tác Dyson) cần công cụ mô phỏng số chính xác và kiểm soát sai số thuật toán.
  • Bạn đọc quan tâm đến kỹ thuật ước lượng giải tích ngẫu nhiên, phương pháp xấp xỉ hàm Yamada-Watanabe và phương pháp Monte Carlo đa cấp.

Câu hỏi thường gặp

1. Vì sao lược đồ Euler-Maruyama cổ điển (PEM) bị phân kỳ đối với phương trình có hệ số dịch chuyển tăng trên tuyến tính?

Khi hệ số dịch chuyển $b(x)$ tăng nhanh hơn tuyến tính (ví dụ $b(x) = ax - bx^3$ với $b>0$), ở các bước rời rạc có bước nhảy Brown lớn, giá trị nghiệm xấp xỉ của lược đồ PEM bị phóng đại mất kiểm soát, dẫn đến mô men $\mathbb{E}[|X_k|^p]$ phân kỳ ra vô cùng ở thời gian hữu hạn, mặc dù nghiệm đúng vẫn có mô men hữu hạn.

2. Bản chất kỹ thuật của lược đồ Euler-Maruyama khống chế (TEM) là gì?

Lược đồ TEM thay thế hệ số dịch chuyển $b(t, x)$ bằng hệ số đã được khống chế $b_n(t, x) = \frac{b(t, x)}{1 + n^{-\lambda}|b(t, x)|}$ với tham số $\lambda \in [\alpha, 1/2]$. Phép biến đổi này giữ cho hệ số dịch chuyển bị chặn ở mỗi bước chia lưới $n$, ngăn chặn sự bùng nổ mô men nhưng vẫn bảo đảm tiệm cận về hàm gốc khi $n \to \infty$.

3. Kỹ thuật xấp xỉ Yamada-Watanabe đóng vai trò gì trong các chứng minh của luận án?

Hàm trị tuyệt đối $|x|$ không khả vi cấp hai tại điểm $0$ nên không thể áp dụng trực tiếp công thức Itô. Kỹ thuật Yamada-Watanabe xây dựng họ hàm trơn $\phi_{\delta\varepsilon}(x) \in C^2$ xấp xỉ đều hàm $|x|$, cho phép áp dụng công thức Itô và triệt tiêu thành phần kỳ dị của hệ số khuếch tán liên tục $(1/2+\alpha)$-Hölder thông qua các đánh giá chặn đạo hàm cấp hai $\phi''_{\delta\varepsilon}$.

4. Ưu điểm nổi bật của lược đồ Milstein nửa ẩn (SIM) đối với hệ hạt Dyson là gì?

Hệ hạt Dyson có hệ số dịch chuyển chứa thành phần $\frac{\gamma_{ij}}{X_i - X_j}$ nổ vô cùng tại biên miền $\Delta_d$. Lược đồ Milstein nửa ẩn vừa giải quyết tính kỳ dị này để đảm bảo các hạt xấp xỉ không bao giờ va chạm vào nhau ($X_1(t) < X_2(t) < \dots < X_d(t)$ hầu chắc chắn), vừa tận dụng thành phần đạo hàm của hệ số khuếch tán để nâng cao bậc hội tụ mạnh so với lược đồ Euler nửa ẩn thông thường.

Kết luận

Luận án của nghiên cứu sinh Lương Đức Trọng đã giải quyết có hệ thống bài toán giải số cho các lớp phương trình vi phân ngẫu nhiên có hệ số không chính qui. Tác giả đã chứng minh sự tồn tại duy nhất nghiệm, thiết lập tốc độ hội tụ mạnh cho lược đồ Euler-Maruyama khống chế dưới điều kiện hệ số dịch chuyển tăng trên tuyến tính và hệ số khuếch tán liên tục Hölder, đồng thời phát triển các lược đồ cải tiến bảo toàn tính không âm và lược đồ Milstein nửa ẩn bảo toàn tính không va chạm cho hệ Dyson. Các kết quả đạt được đóng góp hoàn thiện cơ sở lý thuyết giải số ngẫu nhiên và cung cấp các thuật toán mô phỏng tin cậy cho các mô hình ứng dụng trong vật lý và toán tài chính.