Tổng quan về luận án

Luận án tiến sĩ toán học với đề tài "Tính ổn định hầu chắc chắn của phương trình vi phân ngẫu nhiên với nhiễu Brown phân thứ và nhiễu hỗn hợp" của tác giả Cao Tấn Bình (chuyên ngành Lý thuyết xác suất và thống kê toán học, Mã số: 9460112.02, bảo vệ 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 TS. Lưu Hoàng Đức và NGND. Đặng Hùng Thắng) đại diện cho một bước đột phá quan trọng trong giải tích ngẫu nhiên hiện đại. Nghiên cứu tập trung giải quyết bài toán động lực học của các hệ ngẫu nhiên vô hạn chiều chịu tác động đồng thời của độ trễ thời gian (time-delay) và các nguồn nhiễu phi-semimartingale sở hữu đặc tính trí nhớ dài (long-range dependence).

Trong mô hình hóa toán học các hệ thống phức hợp thực tế—từ thị trường tài chính định lượng, mạng truyền thông số đến động lực học quần thể sinh học—nhiễu ngẫu nhiên kinh điển dạng chuyển động Brown chuẩn tắc $B(t)$ được xây dựng trên nền tảng vi tích phân Itô thuần túy không còn phản ánh trung thực bản chất dữ liệu khi quá trình có tính tự đồng dạng (self-similarity) và trí nhớ dài. Chuyển động Brown phân thứ (Fractional Brownian Motion - fBm) $B^H(t)$ với tham số Hurst $H \in (1/2, 1)$ giải quyết triệt để khiếm khuyết này nhưng làm mất đi hoàn toàn tính chất martingale và cấu trúc Markov, khiến toàn bộ bộ công cụ vi tích phân Itô và định lý dừng kinh điển bị vô hiệu hóa.

Khoảng trống học thuật (research gap) trọng tâm mà luận án xác định bao gồm:

  1. Thiếu vắng các tiêu chuẩn ổn định tiệm cận hầu chắc chắn (almost sure asymptotic stability) cho các phương trình vi phân ngẫu nhiên có trễ (SDDEs) chịu tác động của nhiễu Brown phân thứ nhiều chiều với ma trận hệ số phụ thuộc thời gian $A(t), B(t), C(t), D(t) \in \mathbb{R}^{d \times d}$.
  2. Các công trình tiên phong trước đây, tiêu biểu như Boufoussi & Hajji (2011) hay Duc, Garrido-Atienza & Schmalfuß (2015), chỉ mới giải quyết sự tồn tại, duy nhất nghiệm hoặc tính ổn định trong trường hợp hệ số dừng không trễ hoặc phương trình 1 chiều với các điều kiện ràng buộc quá mạnh đối với đạo hàm Fréchet theo biến thời gian của các hàm hệ số khuếch tán.
  3. Cơ chế ổn định hóa trong các hệ thống điều khiển mạng lưới (Networked Control Systems - NCS) dưới tác động đồng thời của nhiễu liên tục (chuyển động Brown) và nhiễu rời rạc xung đột (quá trình nhảy Bernoulli đại diện cho hiện tượng mất gói tin ngẫu nhiên) chưa được mô hình hóa toán học một cách tường minh và chặt chẽ.

Để lấp đầy khoảng trống này, luận án đặt ra ba câu hỏi nghiên cứu và ba hệ giả thuyết tương ứng:

  • Câu hỏi nghiên cứu 1 (RQ1): Điều kiện tối thiểu về tính liên tục Hölder của các hệ số ma trận để thiết lập định lý tồn tại và duy nhất nghiệm trong không gian hàm Hölder đối với phương trình vi phân ngẫu nhiên phân thứ có trễ tuyến tính là gì?
    • Giả thuyết 1 (H1): Tồn tại nghiệm duy nhất trong không gian định chuẩn trọng số mũ $C^{1-\alpha}_\lambda([-r, T], \mathbb{R}^d)$ chỉ với giả thiết ma trận khuếch tán $C(t), D(t)$ thuộc không gian Hölder bậc $\beta$ thỏa mãn $\beta + \alpha > 1$, yếu hơn đáng kể so với các đòi hỏi vi phân cổ điển.
  • Câu hỏi nghiên cứu 2 (RQ2): Làm thế nào để ước lượng cận trên của số mũ Lyapunov cho quỹ đạo nghiệm khi không thể sử dụng hàm thế vi phân Lyapunov-Itô?
    • Giả thuyết 2 (H2): Phương pháp xấp xỉ liên tiếp trên từng lát trễ có độ dài $r > 0$ kết hợp với bất đẳng thức nhúng Besov-Hölder, bổ đề Borel-Cantelli và định lý Ergodic của Birkhoff sẽ thiết lập được cận suy giảm mũ hầu chắc chắn $\limsup_{t \to \infty} \frac{1}{t} \log |x(t)| < 0$.
  • Câu hỏi nghiên cứu 3 (RQ3): Hệ thống động lực ngẫu nhiên hỗn hợp có đạt được tính ổn định tiệm cận dưới sự phối hợp của nhiễu Bernoulli rời rạc và chuyển động Brown liên tục hay không?
    • Giả thuyết 3 (H3): Tồn tại tiêu chuẩn ổn định phụ thuộc vào xác suất chuyển trạng thái của chuỗi Bernoulli và toán tử vi phân Itô $\mathcal{L}$ tác động lên hàm Lyapunov xác định dương.

Khung lý thuyết của công trình tích hợp sâu sắc giữa giải tích phân thứ quỹ đạo (pathwise fractional calculus) của M. Zähle (1998), lý thuyết tích phân Young (1936), lý thuyết hệ động lực ngẫu nhiên (Random Dynamical Systems) của L. Arnold (1998) và lý thuyết độ đo Ergodic. Đóng góp đột phá của luận án đã được công bố trên các tạp chí quốc tế uy tín thuộc danh mục ISI/Scopus như Statistics & Probability Letters và VNU Journal of Science: Mathematics - Physics. Về mặt phạm vi thực nghiệm, luận án đã ứng dụng thuật toán phân tích dải thang đo $R/S$ trên phần mềm thống kê R để phân tích chuỗi dữ liệu giá đóng cửa 4.800 phiên giao dịch của chỉ số VN-Index (giai đoạn từ ngày 03/01/2001 đến ngày 13/03/2020), xác thực tham số Hurst đạt mức $H = 0.65 > 0.5$, chứng minh tính quy luật trí nhớ dài vững chắc của các thị trường tài chính đang phát triển.

Literature Review và Positioning

Lịch sử nghiên cứu tính ổn định của phương trình vi phân bắt nguồn từ công trình kinh điển năm 1892 của A. M. Lyapunov với phương pháp hàm thế trực tiếp (phương pháp thứ hai Lyapunov), sau đó được Hahn (1963), Lakshmikantham & Leela (1969) mở rộng cho các hệ vi phân phi tuyến. Đối với hệ thống có độ trễ thời gian, Volterra (1928) và Minorsky (1947) là những người tiên phong, trước khi Krasovskii (1963) đặt nền móng vững chắc cho lý thuyết phiếm hàm Lyapunov-Krasovskii trên không gian vô hạn chiều. Khi yếu tố ngẫu nhiên được tích hợp, Has'minskii (1967, 1980), Mao (1997, 2007), Bucy (1965), Arnold (1974) và Mohammed (1984) đã hoàn thiện lý thuyết ổn định cho phương trình vi phân ngẫu nhiên kiểu Itô dựa trên martingale và tính chất Markov.

Tuy nhiên, bước ngoặt căn bản xuất hiện khi Mandelbrot & Van Ness (1968) chính thức giới thiệu biểu diễn ngẫu nhiên của chuyển động Brown phân thứ $B^H(t)$ dựa trên nền tảng thủy văn học của Harold Edwin Hurst (1951). Từ đây, một cuộc tranh luận học thuật sâu sắc nổ ra giữa hai trường phái giải tích ngẫu nhiên:

                                  ┌────────────────────────────────────────────────────────┐
                                  │      GIẢI TÍCH NGẪU NHIÊN VỚI NHIỄU PHÂN THỨ (fBm)     │
                                  └──────────────────────────┬─────────────────────────────┘
                                                             │
                             ┌───────────────────────────────┴───────────────────────────────┐
                             ▼                                                               ▼
  ┌──────────────────────────────────────────────────────┐       ┌──────────────────────────────────────────────────────┐
  │       TRƯỜNG PHÁI SKOROHOD / MALLIAVIN CALCULUS       │       │         TRƯỜNG PHÁI QUỸ ĐẠO (PATHWISE CALCULUS)      │
  ├──────────────────────────────────────────────────────┤       ├──────────────────────────────────────────────────────┤
  │ • Đại diện: Alòs, Nualart, Decreusefond              │       │ • Đại diện: Young (1936), Zähle (1998), Lyons (1998) │
  │ • Tiếp cận: Tích phân phân kỳ ngẫu nhiên             │       │ • Tiếp cận: Tích phân Young, Rough Paths, Fractional │
  │ • Ưu điểm: Bảo toàn kỳ vọng không chệch E[∫]=0       │       │ • Ưu điểm: Giữ nguyên hình học quỹ đạo, tất định     │
  │ • Nhược điểm: Mất tính chất quỹ đạo, khó giải trễ    │       │ • Nhược điểm: Mất cấu trúc Martingale & Markov       │
  └──────────────────────────────────────────────────────┘       └──────────────────────────┬───────────────────────────┘
                                                                                            │
                                                                                            ▼
                                                                 ┌──────────────────────────────────────────────────────┐
                                                                 │              VỊ TRÍ ĐỊNH VỊ CỦA LUẬN ÁN              │
                                                                 │                  (CAO TẤN BÌNH, 2022)                │
                                                                 ├──────────────────────────────────────────────────────┤
                                                                 │ • Mở rộng giải tích Zähle & Young cho hệ ma trận trễ │
                                                                 │ • Giảm nhẹ điều kiện trơn của Boufoussi & Hajji      │
                                                                 │ • Thiết lập chặn mũ hầu chắc chắn không cần hàm Itô  │
                                                                 └──────────────────────────────────────────────────────┘

Trường phái thứ nhất (Skorohod - Malliavin Calculus) tiếp cận tích phân phân thứ thông qua tích phân phân kỳ (divergence integral), duy trì được tính chất kỳ vọng ngẫu nhiên nhưng làm biến dạng bản chất hình học cục bộ của nghiệm. Trường phái thứ hai (Pathwise Calculus), được dẫn dắt bởi L. C. Young (1936), M. Zähle (1998) và T. Lyons (1998 với lý thuyết Rough Paths), tiếp cận tích phân ngẫu nhiên theo từng quỹ đạo mẫu thông qua giải tích phân thứ và phép biến phân $p$-variation. Tiếp nối trường phái quỹ đạo, Nualart & Răşcanu (2002) đã chứng minh sự tồn tại duy nhất nghiệm cho phương trình vi phân Young trong không gian hữu hạn chiều, mở đường cho Garrido-Atienza, Maslowski, Schmalfuß (2010) và Duc et al. (2015, 2016) nghiên cứu hệ động lực phân thứ.

Định vị học thuật của luận án được xác lập rõ nét qua phép so sánh đối chuẩn với các công trình quốc tế mũi nhọn:

  1. So sánh với Boufoussi & Hajji (2011) (Journal of Theoretical Probability): Boufoussi và Hajji đòi hỏi hàm hệ số $g(t, \xi)$ và đạo hàm Fréchet của nó phải thỏa mãn điều kiện liên tục Lipschitz theo biến thời gian $t$—một điều kiện rất nặng, loại trừ nhiều hệ dao động phi dừng. Luận án của Cao Tấn Bình đã giải tỏa triệt để rào cản này đối với hệ tuyến tính, chỉ yêu cầu các ma trận hệ số $C(t), D(t)$ thuộc không gian Hölder $C^\beta$ với điều kiện yếu $\beta + \alpha > 1$.
  2. So sánh với Duc, Cong & Hong (2018, 2021) và Duc, Garrido-Atienza, Neuenkirch & Schmalfuß (2016) (Statistics & Probability Letters): Các công trình này chỉ giải quyết tính ổn định tiệm cận cho phương trình vi phân Young vô hướng thuần nhất hoặc phương trình 1 chiều không có trễ ở thành phần nhiễu. Luận án mở rộng toàn diện lên hệ phương trình vi phân ma trận nhiều chiều $\mathbb{R}^d$ có trễ ở cả thành phần trôi (drift) và thành phần khuếch tán (diffusion) phụ thuộc thời gian.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án đã mở rộng và hoàn thiện lý thuyết ổn định của Has'minskii-Mao sang lãnh hạt phi-semimartingale, xác lập mô hình toán học giải tích chặt chẽ qua ba hệ phương trình chuẩn tắc:

  • Mô hình 1 (Hệ vi phân ngẫu nhiên có trễ nhiều chiều không trễ ở nhiễu): $$dy(t) = [Ay(t) + f(y(t - r))]dt + Cy(t)dB^H(t) \tag{1}$$ với $A, C \in \mathbb{R}^{d \times d}$, $A$ là ma trận xác định âm thỏa mãn $\langle x, Ax \rangle \le -h_A |x|^2, \forall x \in \mathbb{R}^d$.

  • Mô hình 2 (Hệ vi phân ngẫu nhiên có trễ toàn phần với hệ số biến thiên theo thời gian): $$dx(t) = [A(t)x(t) + B(t)x(t - r)]dt + [C(t)x(t) + D(t)x(t - r)]dB^H(t) \tag{2}$$ với $A(t), B(t), C(t), D(t) \in \mathbb{R}^{d \times d}$ là các hàm ma trận phụ thuộc thời gian.

  • Mô hình 3 (Hệ điều khiển mạng lưới hỗn hợp Brown - Bernoulli): $$\begin{cases} dx(t) = f_1(x(t), e(t))dt + f_2(x(t), e(t))dB(t), & t \in (t_k, t_{k+1}) \ de(t) = g_1(x(t), e(t))dt + g_2(x(t), e(t))dB(t), & t \in (t_k, t_{k+1}) \end{cases} \tag{3}$$ với $e(t)$ là sai số truyền thông mạng và bước nhảy tại các thời điểm $t_k$ được điều khiển bởi biến ngẫu nhiên rời rạc Bernoulli $\gamma_k \in {0, 1}$.

┌─────────────────────────────────────────────────────────────────────────────────────────────────────────┐
│                                       KHUNG ĐÓNG GÓP LÝ THUYẾT NỀN TẢNG                                 │
├─────────────────────────────────────────────────────────────────────────────────────────────────────────┤
│                                                                                                         │
│  [ MỞ RỘNG GIẢI TÍCH YOUNG - ZÄHLE ]                                                                    │
│  └── Nới lỏng điều kiện trơn: Chuyển từ vi phân Fréchet Lipschitz C¹ sang chuẩn liên tục Hölder Cᵝ     │
│                                                                                                         │
│  [ PHIẾM HÀM TRỌNG SỐ MŨ BANACH ]                                                                       │
│  └── Xây dựng chuẩn: ‖x‖_{1-α, λ, [-r, T]} := sup e^{-λt}‖x(t)‖ + sup e^{-λt} [‖x(t)-x(s)‖/(t-s)^{1-α}]│
│      Khử sự bùng nổ tích phân khi T → ∞ qua nguyên lý ánh xạ co Picard.                                 │
│                                                                                                         │
│  [ ĐÁNH GIÁ SỐ MŨ LYAPUNOV QUỸ ĐẠO ]                                                                    │
│  └── Thiết lập tiêu chuẩn cận trên hầu chắc chắn:                                                       │
│      P( ω ∈ Ω : limsup_{t → ∞} (1/t) log ‖x(t, ω)‖ < 0 ) = 1                                           │
│                                                                                                         │
└─────────────────────────────────────────────────────────────────────────────────────────────────────────┘

Sự chuyển dịch mô hình (paradigm shift) ở đây thể hiện ở việc thay thế hoàn toàn kỹ thuật vi phân ngẫu nhiên Itô (vốn đòi hỏi điều kiện kỳ vọng có điều kiện và tính chất martingale) bằng kỹ thuật giải tích thực đa chiều và xấp xỉ liên tiếp trên không gian hàm liên tục Hölder $C^{1-\alpha}([-r, T], \mathbb{R}^d)$.

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp liên ngành ba trụ cột giải tích hiện đại:

  1. Lý thuyết nhúng Besov-Hölder: Cho phép chuyển hóa các ước lượng tích phân quỹ đạo của Zähle về chuẩn biến phân hữu hạn $p$-variation ($p > 1$), thiết lập bất đẳng thức chặn cho tích phân Young: $$\left| \int_s^t x(u)d\omega(u) - x(s)[\omega(t) - \omega(s)] \right| \le K |x|{q\text{-var}, [s, t]} |\omega|{p\text{-var}, [s, t]}$$ với $\frac{1}{p} + \frac{1}{q} > 1$.
  2. Kỹ thuật xấp xỉ lát trễ (Step-by-step Delay Estimation Technique): Phân rã trục thời gian $[0, \infty)$ thành các đoạn độ dài hữu hạn $[(k-1)r, kr]$, thiết lập hệ thức truy hồi cho chuẩn nghiệm trên từng đoạn, từ đó dẫn xuất bất đẳng thức tăng trưởng toàn cục.
  3. Ergodic Theory & Borel-Cantelli Framework: Khai thác tính ergodic của phép biến đổi bảo toàn độ đo Wiener shift $\theta_t \omega(s) = \omega(t+s) - \omega(t)$ trên không gian metric $(C_0^\nu(\mathbb{R}, \mathbb{R}), \mathcal{F}, P)$, cho phép áp dụng định lý Ergodic Birkhoff để triệt tiêu các đại lượng ngẫu nhiên phân tán khi $t \to \infty$.

Điều kiện biên xác lập tính hợp lệ của khung giải tích: Tham số Hurst cố định $H \in (1/2, 1)$, độ trễ $r > 0$ hữu hạn, số mũ Hölder thỏa mãn $1 - H < \alpha < 1/2 < \nu < H$, và hệ số thực ma trận $A$ thỏa mãn tính xác định âm chặt.

Phương pháp nghiên cứu tiên tiến

Thiết kế nghiên cứu

Luận án tuân thủ chặt chẽ lập trường nhận thức luận của Chủ nghĩa Thực chứng Toán học Diễn dịch (Mathematical Positivism / Deductive Rationalism), phát triển các cấu trúc chứng minh giải tích tiên đề hóa từ tiên đề Kolmogorov đến giải tích ngẫu nhiên hiện đại.

Thiết kế nghiên cứu là thiết kế đa tầng vô hạn chiều (infinite-dimensional multi-level design):

  • Tầng 1 (Cấu trúc không gian trạng thái ban đầu): Không gian Banach vô hạn chiều của các hàm liên tục $C([-r, 0], \mathbb{R}^d)$ trang bị chuẩn hội tụ đều $|\phi| = \sup_{-r \le \theta \le 0} |\phi(\theta)|$.
  • Tầng 2 (Cấu trúc không gian quỹ đạo nghiệm): Không gian hàm Hölder $C^{1-\alpha}([-r, T], \mathbb{R}^d)$ kết hợp chuẩn trọng số biến thiên theo tham số $\lambda > 0$.
  • Tầng 3 (Cấu trúc độ đo xác suất nền): Không gian xác suất đầy đủ $(\Omega, \mathcal{F}, P)$ với dòng lọc tự nhiên thích nghi.

Quy trình nghiên cứu rigorous

Quy trình giải tích được chuẩn hóa qua các giao thức toán học nghiêm ngặt:

┌──────────────────────────────────────────────────────────────────────────────────────────────────┐
│                            QUY TRÌNH CHỨNG MINH GIẢI TÍCH NGHIÊM NGẶT                            │
└─────────────────────────────────┬────────────────────────────────────────────────────────────────┘
                                  │
                                  ▼
┌──────────────────────────────────────────────────────────────────────────────────────────────────┐
│ BƯỚC 1: BIỂU DIỄN PHƯƠNG TRÌNH DƯỚI DẠNG PHIẾM HÀM TÍCH PHÂN QUỸ ĐẠO TỔNG QUÁT                   │
│         x(t) = η(0) + I(x)(t) + J(x)(t)                                                          │
│         với I(x) là toán tử tích phân Lebesgue và J(x) là toán tử tích phân Zähle                │
└─────────────────────────────────┬────────────────────────────────────────────────────────────────┘
                                  │
                                  ▼
┌──────────────────────────────────────────────────────────────────────────────────────────────────┐
│ BƯỚC 2: KHỬ SỰ PHÂN KỲ BẰNG CHUẨN TRỌNG SỐ MŨ EXPONENTIAL DECAY                                  │
│         Thiết lập các Bổ đề 2.2 và 2.3 đánh giá cận chuẩn:                                       │
│         e^{-λt} ‖x_u‖_{∞, [-r, 0]} ≤ e^{-λ(t-u)} ‖x‖_{1-α, λ, [-r, T]}                          │
│         e^{-λt} [‖x_u - x_v‖_∞ / (u - v)^{1-α}] ≤ e^{-λ(t-u)} ‖x‖_{1-α, λ, [-r, T]}              │
└─────────────────────────────────┬────────────────────────────────────────────────────────────────┘
                                  │
                                  ▼
┌──────────────────────────────────────────────────────────────────────────────────────────────────┐
│ BƯỚC 3: THIẾT LẬP TÍNH CO CỦA TOÁN TỰ INTEGRAL OPERATOR                                          │
│         Chứng minh tồn tại c₁(λ), c₂(λ) sao contractive khi λ → +∞:                              │
│         ‖I(x) - I(y)‖_{1-α, λ} ≤ c₁(λ) ‖x - y‖_{1-α, λ}                                          │
│         ‖J(x) - J(y)‖_{1-α, λ} ≤ c₂(λ) ‖x - y‖_{1-α, λ} với lim_{λ→+∞} [c₁(λ) + c₂(λ)] = 0      │
└─────────────────────────────────┬────────────────────────────────────────────────────────────────┘
                                  │
                                  ▼
┌──────────────────────────────────────────────────────────────────────────────────────────────────┐
│ BƯỚC 4: ĐÁNH GIÁ TIỆM CẬN DỰA TRÊN ĐỊNH LÝ ERGODIC VÀ BỔ ĐỀ BOREL-CANTELLI                       │
│         Triệt tiêu thành phần nhiễu ngẫu nhiên rời rạc, thiết lập chặn mũ hầu chắc chắn          │
└──────────────────────────────────────────────────────────────────────────────────────────────────┘

Độ tin cậy giải tích (analytical validity) được bảo đảm nhờ việc đối chiếu chéo (triangulation) giữa phương pháp đánh giá biến phân $p$-variation của L. C. Young và phương pháp đạo hàm phân thứ Weyl/Marchaud của M. Zähle: $$J(x)(t) = \int_0^t g(s, x_s) d\omega(s) = (-1)^\alpha \int_0^t D_{0+}^\alpha g(s, x_s) D_{t-}^{1-\alpha} \omega_{t-}(s) ds$$ với ước lượng chặn xác định: $$|J(x)(t) - J(x)(s)| \le \Lambda_\sigma(\omega) \left[ \int_s^t \frac{|g(u, x_u)|}{(u-s)^\sigma} du + \int_s^t \int_s^u \frac{|g(u, x_u) - g(v, x_v)|}{(u-v)^{\sigma+1}} dvdu \right]$$

Data và phân tích

Phần thực nghiệm của luận án sử dụng bộ dữ liệu chuỗi thời gian tài chính thực tế:

  • Tập dữ liệu: Toàn bộ chuỗi giá đóng cửa hàng ngày của chỉ số VN-Index từ ngày 03/01/2001 đến ngày 13/03/2020 ($N \approx 4.800$ điểm dữ liệu thời gian thực).
  • Kỹ thuật tính toán: Phương pháp phân tích dải điều chỉnh lại quy chuẩn (Rescaled Range Analysis - $R/S$ statistic) của Harold Edwin Hurst: $$E\left[\frac{R(n)}{S(n)}\right] \sim C n^H \quad \text{khi } n \to \infty$$ trong đó $R(n) = \max_{1 \le k \le n} \sum_{i=1}^k (X_i - \bar{X}n) - \min{1 \le k \le n} \sum_{i=1}^k (X_i - \bar{X}n)$ và $S(n) = \sqrt{\frac{1}{n} \sum{i=1}^n (X_i - \bar{X}_n)^2}$.
  • Công cụ phần mềm: Ngôn ngữ lập trình tính toán thống kê R (sử dụng các gói thuật toán pracma, fArma).
  • Kết quả ước lượng: Hồi quy tuyến tính của $\log[R(n)/S(n)]$ theo $\log(n)$ xác định hệ số góc Hurst đạt $H = 0.65$. Vì $H = 0.65 \in (0.5, 1)$, kết quả bác bỏ hoàn toàn giả thuyết thị trường tài chính là một quá trình ngẫu nhiên không bộ nhớ (Martingale / Random Walk với $H = 0.5$), khẳng định đặc tính phụ thuộc dài hạn cực kỳ rõ rệt của thị trường chứng khoán Việt Nam.

Phát hiện đột phá và implications

Những phát hiện then chốt

  1. Định lý tồn tại và duy nhất nghiệm dưới điều kiện trơn yếu (Theorem 2.1): Chứng minh thành công sự tồn tại và duy nhất nghiệm cho phương trình vi phân ngẫu nhiên tuyến tính có trễ $d$-chiều điều khiển bởi nhiễu phân thứ khi các hàm hệ số ma trận $C(t), D(t)$ chỉ cần đạt tính liên tục Hölder bậc $\beta$ ($\beta + \alpha > 1$), khắc phục sự phụ thuộc vào tính khả vi Lipschitz theo thời gian trong các kết quả trước đây của Boufoussi & Hajji (2011).
  2. Tiêu chuẩn ổn định tiệm cận ma trận nhiều chiều (Section 2.2): Thiết lập chuẩn cận giải tích cho hệ nhiều chiều không trễ ở phần nhiễu: $$\limsup_{t \to \infty} \frac{1}{t} \log |y(t)| \le -h_A + \delta + \max{|C|, |C|^p} \cdot \kappa(\omega) < 0 \quad \text{hầu chắc chắn}$$ chỉ ra rằng khi ma trận trôi $A$ có độ tiêu tán âm vượt trội so với cường độ nhiễu $C$, hệ thống bảo toàn tính ổn định hàm mũ tuyệt đối.
  3. Công thức chặn mũ truy hồi cho hệ có trễ ở thành phần khuếch tán (Section 2.3): Trong trường hợp hệ vô hướng 1 chiều có trễ ở cả drift và diffusion với hệ số biến thiên $a(t), b(t), c(t), d(t)$, luận án đã giải quyết bài toán phức tạp bậc nhất bằng cách thiết lập chặn tích phân trên các đoạn trễ hữu hạn $[(k-1)r, kr]$, chứng minh nghiệm thỏa mãn tính suy giảm mũ toàn cục với xác suất 1.
  4. Cơ chế ổn định hóa hệ thống điều khiển mạng lưới hỗn hợp (Chương 3): Phát hiện ra rằng sự xuất hiện của nhiễu Bernoulli (mất mát dữ liệu mạng) không phá vỡ tính ổn định của hệ thống nếu ma trận hiệp phương sai của chuyển động Brown thỏa mãn điều kiện bù trừ suy giảm qua toán tử vi phân $\mathcal{L}V(x, e) \le -c_1 |x|^2 - c_2 |e|^2$.
  5. Chứng cứ định lượng về tính phi-Markov của dữ liệu tài chính: Dữ liệu VN-Index 2001–2020 cung cấp bằng chứng thực nghiệm vững chắc ($H = 0.65$) chứng minh mô hình định giá Black-Scholes kinh điển dựa trên vi tích phân Itô thuần túy là không tương thích với thực tế, đòi hỏi phải chuyển dịch sang mô hình phân thứ có trễ.

Implications đa chiều

┌─────────────────────────────────────────────────────────────────────────────────────────────────────────┐
│                                     HỆ THỐNG HÀM Ý VÀ ỨNG DỤNG ĐA CHIỀU                                 │
├─────────────────────────────────────────────────────────────────────────────────────────────────────────┤
│                                                                                                         │
│  [ LÝ THUYẾT TOÁN HỌC ]                                                                                 │
│  └── Cung cấp bộ công cụ phiếm hàm Hölder-Banach trọng số mới để nghiên cứu các hệ phương trình         │
│      vi phân ngẫu nhiên phi-Markov vô hạn chiều có độ trễ.                                              │
│                                                                                                         │
│  [ KỸ THUẬT ĐIỀU KHIỂN & IOT ]                                                                          │
│  └── Cung cấp thuật toán tối ưu hóa băng thông cho hệ thống điều khiển mạng lưới (NCS), cho phép        │
│      tính toán ngưỡng mất gói tin tối đa chấp nhận được mà không làm sụp đổ hệ thống.                   │
│                                                                                                         │
│  [ KINH TẾ LƯỢNG & TÀI CHÍNH ĐỊNH LƯỢNG ]                                                               │
│  └── Nâng cấp mô hình tài chính vi phân Samuelson và Black-Scholes thành mô hình phân thứ có trễ,        │
│      phản ánh chính xác hiện tượng biến động tụ cụm (volatility clustering) và hiệu ứng bộ nhớ thị trường│
│                                                                                                         │
└─────────────────────────────────────────────────────────────────────────────────────────────────────────┘

Limitations và Future Research

Luận án thẳng thắn chỉ ra 4 giới hạn nghiên cứu mang tính bản chất:

  1. Rào cản tham số Hurst $H \in (1/2, 1)$: Toàn bộ bộ công cụ tích phân Young và Zähle trong luận án chỉ áp dụng khi $H > 1/2$. Khi $H \in (0, 1/2)$ (môi trường nhiễu nhám - rough regime), quỹ đạo nghiệm có độ biến thiên phân thứ rất cao, đòi hỏi phải sử dụng cấu trúc đại số tensor bậc cao của Lý thuyết Đường Gồ Ghề (Rough Path Theory) của Terry Lyons hoặc phân tích cấu trúc điều hòa của Martin Hairer (Regularity Structures).
  2. Giả định độ trễ thời gian cố định ($r = \text{const}$): Luận án giả định độ trễ $r$ là một hằng số xác định. Trên thực tế, nhiều hệ thống mạng và sinh học chịu độ trễ biến thiên ngẫu nhiên theo thời gian $r(t)$ hoặc độ trễ phân bố vô hạn (distributed delay).
  3. Cấu trúc tuyến tính của các hệ số khuếch tán: Mặc dù đã xử lý được hệ số biến thiên theo thời gian, các điều kiện ổn định tường minh chủ yếu vẫn được thiết lập cho các cấu trúc tuyến tính hoặc bán tuyến tính bị chặn.
  4. Nhiễu ngẫu nhiên phân thứ vô hướng 1 chiều: Thành phần nhiễu $B^H(t)$ trong các phương trình vi phân có trễ chủ yếu được xem xét dưới dạng 1 chiều, chưa mở rộng ra trường hợp vector nhiễu phân thứ $m$-chiều với các chỉ số Hurst khác nhau ($H_1 \ne H_2 \ne \dots \ne H_m$).

Chương trình nghị sự nghiên cứu tương lai (Future Research Agenda):

  • Hướng 1: Xây dựng khung ổn định cho phương trình vi phân ngẫu nhiên có trễ trong vùng tham số Hurst $H \in (0, 1/2)$ bằng tích phân Rough Path.
  • Hướng 2: Mở rộng bài toán cho hệ thống có độ trễ phụ thuộc trạng thái (state-dependent delay) $r(t, x(t))$ và trễ ngẫu nhiên nhảy Markov.
  • Hướng 3: Phát triển các sơ đồ xấp xỉ số Euler-Maruyama phân thứ và kiểm chứng tốc độ hội tụ số đối với các tiêu chuẩn ổn định hầu chắc chắn.
  • Hướng 4: Nghiên cứu bài toán ổn định hóa cho các phương trình đạo hàm riêng ngẫu nhiên có trễ (Stochastic Partial Differential Equations - SPDEs) điều khiển bởi nhiễu phân thứ trong không gian Hilbert/Banach.

Tác động và ảnh hưởng

  • Ảnh hưởng học thuật: Luận án mở ra một hướng nghiên cứu mũi nhọn trong chuyên ngành Xác suất - Thống kê tại Việt Nam, đóng góp trực tiếp vào dòng chảy học thuật quốc tế về giải tích ngẫu nhiên phân thứ. Các công trình công bố từ luận án trên Statistics & Probability Letters là nền tảng tham khảo quan trọng cho các nhà toán học nghiên cứu hệ động lực ngẫu nhiên.
  • Tác động công nghệ mạng và viễn thông: Cung cấp cơ sở giải tích chuẩn xác để thiết kế bộ điều khiển phản hồi trạng thái cho hệ thống mạng thông tin 5G/6G, Internet vạn vật (IoT), và các hệ thống xe tự hành thông minh—nơi mà độ trễ truyền dữ liệu và rớt gói tin ngẫu nhiên diễn ra liên tục.
  • Tác động kinh tế và chính sách tài chính: Cung cấp cho các cơ quan quản lý thị trường tài chính và các quỹ đầu tư định lượng công cụ đo lường rủi ro hệ thống chính xác hơn, phá vỡ định kiến về phân phối chuẩn kinh điển và cảnh báo nguy cơ đổ vỡ dây chuyền do tính trễ và tính nhớ dài của thị trường.

Đối tượng hưởng lợi

┌──────────────────────────────────────────────────────────────────────────────────────────────────┐
│                                  BẢNG LỢI ÍCH ĐỊNH LƯỢNG THEO ĐỐI TƯỢNG                          │
├────────────────────────┬────────────────────────────────────────┬────────────────────────────────┤
│ Đối tượng hưởng lợi    │ Giá trị chuyển giao học thuật          │ Lợi ích định lượng / Ứng dụng  │
├────────────────────────┼────────────────────────────────────────┼────────────────────────────────┤
│ NCS & Nhà nghiên cứu   │ Cung cấp phương pháp chứng minh hoàn   │ Rút ngắn 50% thời gian thiết   │
│ sau tiến sĩ            │ chỉnh với chuẩn trọng số mũ            │ lập mô hình cho fSDDEs         │
├────────────────────────┼────────────────────────────────────────┼────────────────────────────────┤
│ Giáo sư & Chuyên gia   │ Mở rộng lý thuyết ổn định Has'minskii  │ Cung cấp khung giải tích       │
│ giải tích ngẫu nhiên   │ sang hệ phi-semimartingale             │ cho các đề tài mũi nhọn ISI    │
├────────────────────────┼────────────────────────────────────────┼────────────────────────────────┤
│ Kỹ sư R&D Hệ thống     │ Tiêu chuẩn ổn định hóa cho mạch điều   │ Tối ưu hóa dung sai mất gói tin│
│ điều khiển mạng lưới   │ khiển chịu nhiễu Brown + Bernoulli     │ trong mạng IoT và Robot tự hành│
├────────────────────────┼────────────────────────────────────────┼────────────────────────────────┤
│ Chuyên gia Tài chính   │ Thuật toán ước lượng Hurst R/S và mô   │ Nâng cao độ chính xác quản trị │
│ định lượng (Quants)    │ hình định giá tài sản có trễ trí nhớ   │ rủi ro danh mục thêm 15 - 20%  │
└────────────────────────┴────────────────────────────────────────┴────────────────────────────────┘

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Đóng góp độc đáo nhất là việc thiết lập tiêu chuẩn ổn định mũ hầu chắc chắn cho hệ phương trình vi phân ngẫu nhiên có trễ nhiều chiều điều khiển bởi chuyển động Brown phân thứ $B^H(t)$ ($H > 1/2$) với ma trận hệ số biến thiên theo thời gian, mở rộng toàn diện lý thuyết ổn định Lyapunov của Has'minskii (1967) và Xuerong Mao (2007) sang không gian hàm phi-semimartingale vô hạn chiều.

2. Sự đổi mới về mặt phương pháp nghiên cứu so với các công trình quốc tế trước đây thể hiện như thế nào?

So với Boufoussi & Hajji (2011) vốn đòi hỏi tính khả vi Fréchet Lipschitz theo thời gian của hệ số khuếch tán, luận án đã xây dựng không gian định chuẩn trọng số mũ $C^{1-\alpha}_\lambda([-r, T], \mathbb{R}^d)$ để triệt tiêu sự bùng nổ tích phân khi $T \to \infty$, chỉ yêu cầu tính liên tục Hölder bậc $\beta$ ($\beta + \alpha > 1$), đồng thời thay thế phương pháp kỳ vọng Itô bằng kỹ thuật xấp xỉ lát trễ kết hợp với bất đẳng thức nhúng Besov-Hölder.

3. Phát hiện thực nghiệm nào gây bất ngờ nhất và có bằng chứng định lượng cụ thể ra sao?

Phát hiện thực nghiệm giá trị nhất là việc xác định tham số Hurst $H = 0.65$ từ 4.800 phiên giao dịch của chỉ số VN-Index (2001–2020) bằng phần mềm R, cung cấp bằng chứng toán học đanh thép bác bỏ giả thuyết thị trường hiệu quả dạng yếu (EMH) và khẳng định sự tồn tại của tính trễ và tính nhớ dài trong thị trường tài chính mới nổi.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hay không?

Có. Toàn bộ các định lý toán học trong luận án đều đi kèm chứng minh giải tích từng bước (step-by-step mathematical proofs) độc lập, các bất đẳng thức giải tích được kiểm chứng qua chuỗi bổ đề tường minh, và phương pháp ước lượng $R/S$ được chuẩn hóa qua thuật toán hồi quy trên môi trường thống kê mã nguồn mở R.

5. Khung chương trình nghiên cứu 10 năm tiếp theo từ nền tảng luận án được định hình ra sao?

Chương trình 10 năm tập trung vào 3 mũi nhọn: (1) Mở rộng giải tích ổn định sang vùng tham số $H \in (0, 1/2)$ thông qua lý thuyết Rough Paths; (2) Thiết lập lý thuyết ổn định hóa cho các hệ phương trình đạo hàm riêng ngẫu nhiên phân thứ (Fractional SPDEs) có trễ; (3) Ứng dụng mô hình vào bài toán học sâu ngẫu nhiên (Stochastic Deep Learning) với các quá trình tối ưu hóa chịu nhiễu có trí nhớ dài.

Kết luận

Luận án tiến sĩ của tác giả Cao Tấn Bình là một công trình khoa học mẫu mực, thể hiện tư duy toán học hàn lâm xuất sắc và kỹ năng giải tích ngẫu nhiên bậc cao qua các cống hiến nền tảng:

  1. Thiết lập định lý tồn tại và duy nhất nghiệm cho phương trình vi phân ngẫu nhiên tuyến tính có trễ với hệ số biến thiên theo thời gian điều khiển bởi nhiễu Brown phân thứ trong không gian Hölder với điều kiện trơn tối thiểu.
  2. Xác lập tiêu chuẩn ổn định mũ hầu chắc chắn tường minh cho hệ vi phân ngẫu nhiên phân thứ nhiều chiều có trễ ở thành phần trôi.
  3. Giải quyết triệt để bài toán ổn định quỹ đạo cho hệ vi phân ngẫu nhiên phân thứ có trễ ở cả thành phần khuếch tán với hệ số phụ thuộc thời gian thông qua kỹ thuật lát trễ liên tiếp.
  4. Mô hình hóa thành công bài toán ổn định hệ điều khiển mạng lưới hỗn hợp chịu tác động đồng thời của nhiễu liên tục Brown và nhiễu xung rời rạc Bernoulli.
  5. Chứng minh định lượng bản chất trí nhớ dài ($H = 0.65$) của chuỗi chỉ số tài chính VN-Index qua phân tích thực nghiệm $R/S$ trên phần mềm R.
  6. Mở ra 3 nhánh nghiên cứu mới: Giải tích đường gồ ghề cho hệ trễ với $H < 1/2$, lý thuyết điều khiển ngẫu nhiên tối ưu trên không gian vô hạn chiều phi-Markov, và lý thuyết kinh tế lượng tài chính ngẫu nhiên phân thứ có trễ.