Tổng quan luận án

Luận án tiến sĩ toán học của tác giả Hoàng Thế Tuấn (hoàn thành năm 2016 tại Viện Toán học thuộc Viện Hàn lâm Khoa học và Công nghệ Việt Nam, dưới sự hướng dẫn khoa học của TSKH. Đoàn Thái Sơn và GS. Nguyễn Đình Công) tập trung nghiên cứu lý thuyết định tính của phương trình vi phân phân thứ. Trong khoa học tự nhiên và kỹ thuật, các quá trình tiến hóa có tính chất phụ thuộc vào lịch sử (bộ nhớ dài) thường không thể mô tả trọn vẹn bằng phép tính vi tích phân cổ điển địa phương. Giải tích phân thứ cung cấp công cụ mô phỏng hữu hiệu cho các hiện tượng này. Tuy nhiên, phương trình vi phân phân thứ không sinh ra các toán tử có tính chất nửa nhóm, dẫn đến việc không thể thiết lập hệ động lực theo nghĩa cổ điển và không thể áp dụng trực tiếp các phương pháp giải tích đã được hoàn thiện trong lý thuyết phương trình vi phân thường. Khoảng trống nghiên cứu căn bản mà tác giả xác định là: các bài toán định tính cốt lõi của phương trình vi phân phân thứ—bao gồm lý thuyết số mũ Lyapunov, tính ổn định và không ổn định thông qua phương pháp tuyến tính hóa, cũng như sự tồn tại của đa tạp ổn định quanh điểm cân bằng—gần như chưa được giải quyết một cách hoàn chỉnh và chính xác về mặt toán học.

Mục tiêu và các nhiệm vụ nghiên cứu của luận án được xác định theo bốn nội dung chính:

  1. Nghiên cứu đặc trưng tiệm cận nghiệm của phương trình vi phân phân thứ tuyến tính có hệ số biến thiên; làm rõ giới hạn của số mũ Lyapunov cổ điển và thiết lập một khái niệm số mũ Lyapunov mới phù hợp với cấu trúc phân thứ.
  2. Mô tả cấu trúc phổ của số mũ Lyapunov mới và thiết lập tiêu chuẩn ổn định, ổn định tiệm cận cho nghiệm tầm thường của phương trình tuyến tính thông qua phổ này; tính toán giải tích tường minh cho hệ hai chiều.
  3. Chứng minh tính ổn định tiệm cận và tính không ổn định của điểm cân bằng đối với lớp phương trình vi phân phân thứ phi tuyến thông qua phổ của ma trận tuyến tính hóa, giải quyết giả thuyết đặt ra trong tài liệu chuyên ngành và khắc phục các sai lầm toán học trong các công bố trước đó.
  4. Xây dựng toán tử kiểu Lyapunov–Perron thích hợp nhằm thiết lập định lý về sự tồn tại và cấu trúc của đa tạp ổn định địa phương quanh điểm cân bằng hyperbolic trong không gian hữu hạn chiều tổng quát.

Đối tượng và phạm vi nghiên cứu:

  • Đối tượng nghiên cứu: Phương trình vi phân phân thứ với toán tử đạo hàm phân thứ Caputo cấp $\alpha \in (0, 1)$, bao gồm cả hệ tuyến tính thuần nhất hệ số biến thiên bị chặn, hệ tuyến tính hệ số hằng và hệ phi tuyến nhiễu địa phương.
  • Phạm vi không gian: Các không gian Euclide hữu hạn chiều $\mathbb{R}^d$, không gian phức $\mathbb{C}^d$, mặt cầu đơn vị $S^{d-1}$ và không gian Banach các hàm liên tục, bị chặn $C_\infty(\mathbb{R}_{\ge 0}; X)$ (với $X = \mathbb{R}^d$ hoặc $\mathbb{C}^d$).
  • Phạm vi thời gian: Biến thời gian $t$ xác định trên nửa trục thực không âm $\mathbb{R}_{\ge 0} = [0, \infty)$ hoặc đoạn hữu hạn $[0, T]$.

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

Luận án điểm lại quá trình phát triển của giải tích phân thứ từ nguồn gốc lịch sử đến các nghiên cứu hiện đại. Về mặt phép tính vi phân, hai khái niệm phổ biến nhất là đạo hàm Riemann–Liouville (xây dựng vào nửa đầu thế kỷ 19 bởi Abel, Riemann và Liouville) và đạo hàm Caputo (được M. Caputo đề xuất năm 1969 nhằm giải quyết bài toán cơ học chất nhớt). Tác giả phân tích rõ: đạo hàm Riemann–Liouville gặp hạn chế khi ứng dụng thực tế do các điều kiện ban đầu trong bài toán Cauchy không có ý nghĩa vật lý trực quan, trong khi đạo hàm Caputo cho phép sử dụng điều kiện ban đầu chứa đạo hàm nguyên cổ điển, tương thích với các phép đo thực nghiệm.

Về tài liệu tổng quan và chuyên khảo, luận án ghi nhận cuốn sách đầu tiên hệ thống hóa ứng dụng của giải tích phân thứ là của K. B. Oldham và J. Spanier (1974) [31], theo sau là các chuyên khảo cơ bản của S. G. Samko, A. A. Kilbas, O. I. Marichev, cũng như các công trình của K. Diethelm [18]. Trong lĩnh vực động lực học và số mũ Lyapunov, C. Chen [9] từng xây dựng một dạng số mũ Lyapunov để khảo sát trạng thái hỗn loạn của hệ phân thứ Chen.

Đối với bài toán ổn định của phương trình phi tuyến bằng phương pháp tuyến tính hóa, luận án chỉ ra các mốc nghiên cứu từ công trình khởi đầu của E. Ahmed, A. M. A. El-Saka [2]. Tiếp sau đó, nhiều nhóm tác giả đã nỗ lực chứng minh định lý tuyến tính hóa, tiêu biểu như X. Wen và cộng sự [40], D. Qian và cộng sự [35], L. Chen và cộng sự [9], R. Zhang và cộng sự [42]. Tuy nhiên, qua khảo cứu tài liệu, tác giả luận án khẳng định các chứng minh đã công bố của những nhóm tác giả này đều chứa đựng các khiếm khuyết toán học nghiêm trọng (chẳng hạn như việc áp dụng sai các bất đẳng thức vi tích phân phân thứ hoặc giả định ngầm tính chất nửa nhóm vốn không tồn tại). Về bài toán không ổn định, J.-C. Trigeassou, N. Maamri, J. Sabatier, A. Oustaloup (gọi tắt là nhóm A–M–M trong giả thuyết ổn định) đã đề xuất giả thuyết về mối liên hệ giữa phổ tuyến tính hóa và tính ổn định, không ổn định của điểm cân bằng hệ một chiều, nhưng thiếu chứng minh giải tích tổng quát; ngoài ra chỉ có một số ít nghiên cứu kế tiếp như của C. Tuan và cộng sự [13].

Vị trí của luận án được xác lập tại điểm giao thoa giữa giải tích phân thứ và lý thuyết định tính phương trình vi phân: giải quyết triệt để các khoảng trống về mặt chứng minh toán học giải tích chặt chẽ cho phương pháp tuyến tính hóa (cả chiều ổn định lẫn không ổn định), định nghĩa một đại lượng phổ hữu hiệu thay thế số mũ Lyapunov cổ điển, và thiết lập cơ sở giải tích cho lý thuyết đa tạp bất biến (đa tạp ổn định) vốn chưa từng được giải quyết trên không gian hữu hạn chiều tổng quát.

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

Luận án vận dụng các cơ sở lý thuyết và công cụ giải tích toán học thuần túy:

  • Toán tử giải tích phân thứ: Tích phân phân thứ Riemann–Liouville cấp $\alpha > 0$ ký hiệu là $I_{a+}^\alpha x(t) := \frac{1}{\Gamma(\alpha)} \int_a^t (t-\tau)^{\alpha-1} x(\tau) d\tau$, với $\Gamma(\alpha)$ là hàm Gamma; đạo hàm phân thứ Riemann–Liouville cấp $\alpha$: $D_{a+}^\alpha x(t) := D^m I_{a+}^{m-\alpha} x(t)$; và đạo hàm phân thứ Caputo cấp $\alpha \in (0, 1)$: $^C D_{a+}^\alpha x(t) := I_{a+}^{1-\alpha} D x(t)$. Luận án nhấn mạnh đạo hàm phân thứ không có tính chất nửa nhóm, tức $D_{a+}^{\alpha_1} D_{a+}^{\alpha_2} \neq D_{a+}^{\alpha_1+\alpha_2}$.
  • Hàm Mittag-Leffler hai tham số: Định nghĩa bởi $E_{\alpha,\beta}(z) := \sum_{k=0}^\infty \frac{z^k}{\Gamma(\alpha k + \beta)}$ ($z, \beta \in \mathbb{C}, \alpha > 0$), và hàm một tham số $E_\alpha(z) := E_{\alpha,1}(z)$. Hàm này giữ vai trò hạt nhân trong nghiệm của phương trình vi phân phân thứ, tương tự vai trò của hàm mũ trong phương trình vi phân cổ điển.
  • Phân chia miền phức theo góc argument: Luận án xác định hai miền hình quạt quan trọng trên mặt phẳng phức:
    • Miền ổn định: $\Lambda_\alpha^s := { \lambda \in \mathbb{C} \setminus {0} : |\arg(\lambda)| > \frac{\alpha\pi}{2} }$.
    • Miền không ổn định: $\Lambda_\alpha^u := { \lambda \in \mathbb{C} \setminus {0} : |\arg(\lambda)| < \frac{\alpha\pi}{2} }$.
  • Bất đẳng thức Gronwall suy rộng: Sử dụng công cụ ước lượng nghiệm dạng tích phân phân thứ chặn trên bởi hàm Mittag-Leffler: nếu $|\delta(t)| \le K + \frac{L}{\Gamma(\alpha)} \int_0^t (t-\tau)^{\alpha-1} |\delta(\tau)| d\tau$ thì $|\delta(t)| \le K E_\alpha(Lt^\alpha)$.
  • Phép biến đổi Laplace và công thức biến thiên hằng số: Áp dụng biến đổi Laplace $\mathcal{L}{f(t)}(s) = \int_0^\infty e^{-st} f(t) dt$ để chuyển phương trình vi phân phân thứ sang không gian ảnh, từ đó thiết lập biểu diễn nghiệm thông qua công thức biến thiên hằng số với tích chập chứa hàm Mittag-Leffler.
  • Phương pháp nghiên cứu thực tế: Luận án thuộc hướng nghiên cứu lý thuyết định tính thuần túy trong toán học giải tích. Tác giả kết hợp phương pháp ước lượng tiệm cận giải tích phức (tích phân trên chu tuyến Hankel $\gamma(\varepsilon, \theta)$), phương pháp toán tử và điểm bất động (toán tử co Lyapunov–Perron trên không gian Banach hàm liên tục bị chặn $C_\infty(\mathbb{R}_{\ge 0}; X)$ và không gian định chuẩn có trọng $|\cdot|_w$), phương pháp đại số tuyến tính (phân rã Jordan, kỹ thuật làm nhỏ đường chéo phụ của ma trận khối Jordan), và phương pháp phản chứng giải tích.
  • Nguồn tư liệu: Toàn bộ các mệnh đề, định lý, công thức giải tích được xây dựng và chứng minh suy diễn toán học từ các bổ đề cơ bản của giải tích phân thứ, không sử dụng dữ liệu thực nghiệm hay thống kê số liệu mô phỏng.

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

Chương 1: Giới thiệu tóm tắt về phương trình vi phân phân thứ

Chương 1 hệ thống hóa các kiến thức nền tảng làm tiền đề cho toàn bộ các chứng minh ở các chương sau:

  • Giải tích phân thứ: Trình bày chi tiết định nghĩa tích phân và đạo hàm phân thứ Riemann–Liouville và Caputo; chỉ ra tính khả tích của tích phân phân thứ trên $L^1[a,b]$; xác lập mối quan hệ giữa đạo hàm Caputo và Riemann–Liouville thông qua khai triển chuỗi chứa giá trị đạo hàm tại điểm đầu mút: $$^C D_{a+}^\alpha x(t) = D_{a+}^\alpha x(t) - \sum_{j=0}^{m-1} \frac{(t-a)^{j-\alpha}}{\Gamma(j-\alpha+1)} x^{(j)}(a)$$
  • Bản chất phi địa phương của nghiệm: Phân tích phương trình tích phân tương đương của bài toán Cauchy $^C D_{0+}^\alpha x(t) = f(t, x(t)), x(0) = x_0$: $$x(t) = x_0 + \frac{1}{\Gamma(\alpha)} \int_0^t (t-\tau)^{\alpha-1} f(\tau, x(\tau)) d\tau$$ Công thức cho thấy giá trị nghiệm tại thời điểm tương lai $t$ phụ thuộc vào toàn bộ lịch sử trạng thái của hệ trên đoạn $[0, t]$, thể hiện bản chất không có tính chất địa phương theo thời gian. Tác giả chứng minh định lý tồn tại duy nhất nghiệm địa phương trên $[0, T^*]$ (với hàm $f$ liên tục Lipschitz địa phương) và nghiệm toàn cục trên $\mathbb{R}_{\ge 0}$ (khi $f$ thỏa mãn điều kiện Lipschitz toàn cục theo biến không gian).
  • Hàm Mittag-Leffler: Giới thiệu biểu diễn tích phân của hàm $E_{\alpha,\beta}(z)$ trên các cung chu tuyến $\gamma(\varepsilon, \theta)$ và thiết lập khai triển tiệm cận bậc $p$:
    • Trong miền $|\arg(z)| \le \frac{\alpha\pi}{2}$: hàm tăng trưởng theo cấp số mũ kết hợp chuỗi lũy thừa nghịch đảo: $$E_{\alpha,\beta}(z) = \frac{1}{\alpha} z^{\frac{1-\beta}{\alpha}} \exp(z^{1/\alpha}) - \sum_{k=1}^p \frac{z^{-k}}{\Gamma(\beta - \alpha k)} + O(|z|^{-1-p})$$
    • Trong miền ổn định $\Lambda_\alpha^s$: triệt tiêu phần hàm mũ, hàm suy giảm theo quy tắc hàm lũy thừa âm: $$E_{\alpha,\beta}(z) = - \sum_{k=1}^p \frac{z^{-k}}{\Gamma(\beta - \alpha k)} + O(|z|^{-1-p})$$
  • Công thức biến thiên hằng số: Sử dụng biến đổi Laplace đối với toán tử Caputo $\mathcal{L}{^C D_{0+}^\alpha f(t)}(s) = s^\alpha \mathcal{L}{f(t)}(s) - s^{\alpha-1} f(0)$, tác giả chứng minh nghiệm của phương trình $^C D_{0+}^\alpha x(t) = Ax(t) + f(x(t))$ biểu diễn tường minh dưới dạng: $$\phi(t, x_0) = E_\alpha(t^\alpha A)x_0 + \int_0^t (t-\tau)^{\alpha-1} E_{\alpha,\alpha}((t-\tau)^\alpha A) f(\phi(\tau, x_0)) d\tau$$

Chương 2: Số mũ Lyapunov cho phương trình vi phân phân thứ

Chương 2 giải quyết bài toán định lượng tốc độ tiệm cận của phương trình phân thứ tuyến tính thuần nhất hệ số biến thiên $^C D_{0+}^\alpha x(t) = A(t)x(t)$ với $A(t)$ liên tục và bị chặn ($M := \sup_{t\ge 0} |A(t)| < \infty$):

  • Giới hạn của số mũ Lyapunov cổ điển: Số mũ Lyapunov cổ điển của hàm $f$ được xác định bởi $\chi(f) := \limsup_{t\to\infty} \frac{1}{t} \log |f(t)|$. Định lý 2.2 chứng minh rằng: với mọi nghiệm không tầm thường ($x_0 \neq 0$), số mũ Lyapunov cổ điển luôn thỏa mãn: $$\chi(\phi(\cdot, x_0)) \ge 0$$ Bằng phương pháp phản chứng, nếu $\chi < 0$, nghiệm suy giảm nhanh hơn hàm mũ, dẫn tới tích phân trong biểu diễn nghiệm hội tụ về 0 và mâu thuẫn với $|x_0| \neq 0$. Vì $\chi$ không thể nhận giá trị âm ngay cả khi nghiệm suy giảm tiệm cận về 0 theo hàm lũy thừa, số mũ Lyapunov cổ điển mất khả năng phân loại tính ổn định của hệ vi phân phân thứ.
  • Thiết lập số mũ Lyapunov phân thứ: Nhận thấy hàm Mittag-Leffler thực $E_\alpha: \mathbb{R} \to \mathbb{R}{\ge 0}$ đơn điệu tăng ngặt, có hàm ngược $\log\alpha^M: \mathbb{R}{>0} \to \mathbb{R}$, tác giả định nghĩa số mũ Lyapunov phân thứ cấp $\alpha$: $$\chi\alpha(f) := \limsup_{t\to\infty} \frac{1}{t^\alpha} \log_\alpha^M(|f(t)|)$$ Đồng thời thiết lập liên hệ giải tích giữa $\chi(f)$ và $\chi_\alpha(f)$:
    • Nếu $\chi(f) > 0$ thì $\chi_\alpha(f) = (\chi(f))^\alpha > 0$.
    • Nếu $\lim\sup_{t\to\infty} t^\alpha |f(t)| < \infty$ thì $\chi_\alpha(f) = -\frac{1}{\lim\sup_{t\to\infty} t^\alpha |f(t)| \cdot \Gamma(1-\alpha)} < 0$.
  • Cấu trúc phổ Lyapunov phân thứ: Phổ $\Sigma_\alpha := { \chi_\alpha(\phi(\cdot, x_0)) : x_0 \in \mathbb{R}^d \setminus {0} }$ bị chặn trên bởi $M$ ($\Sigma_\alpha \subset (-\infty, M]$). Phần phổ không âm $\Sigma_\alpha \cap \mathbb{R}{\ge 0}$ chứa tối đa $d$ giá trị phân biệt $\lambda_j < \lambda{j-1} < \dots < \lambda_1$. Tập $S := { x_0 \in \mathbb{R}^d : \chi_\alpha(\phi(\cdot, x_0)) < 0 }$ và các tập mức $E_i := { x_0 \in \mathbb{R}^d : \chi_\alpha(\phi(\cdot, x_0)) \le \lambda_i }$ tạo thành một dãy không gian con tuyến tính lồng nhau: $S \subsetneq E_j \subsetneq E_{j-1} \subsetneq \dots \subsetneq E_1 = \mathbb{R}^d$.
  • Liên hệ với tính ổn định: Nghiệm tầm thường ổn định khi và chỉ khi $\Sigma_\alpha \subset (-\infty, 0]$; và ổn định tiệm cận khi $\Sigma_\alpha \subset (-\infty, 0)$.
  • Cấu trúc phổ trên mặt cầu đơn vị $S^{d-1}$: Tập $\Lambda_\alpha := { \chi_\alpha(\phi(\cdot, x_0)) : x_0 \in S^{d-1} }$ được chứng minh là tập hợp gồm một đoạn liên tục âm ghép với các điểm rời rạc không âm: $\Lambda_\alpha = [a, b] \cup {\lambda_1, \dots, \lambda_j}$ với $b < 0$ (nếu $S \neq \emptyset$).
  • Tính toán tường minh cho hệ hai chiều ($d=2$): Luận án đưa ra công thức giải tích chính xác số mũ $\chi_\alpha$ cho hệ tuyến tính hai chiều hệ số hằng theo ba dạng chuẩn của ma trận hệ số $A$:
Dạng ma trận hệ số $A$ Điều kiện tham số và tọa độ điểm đầu $(x_1, x_2)$ Số mũ Lyapunov phân thứ $\chi_\alpha(\phi(\cdot, x))$
$\begin{pmatrix} \lambda & 0 \ 0 & \gamma \end{pmatrix}$ ($\lambda \ge \gamma$) $x_1 x_2 \neq 0, \lambda < 0$ $\max \left{ \frac{\lambda}{|x_1|}, \frac{\gamma}{|x_2|} \right}$
$\begin{pmatrix} \lambda & 0 \ 0 & \gamma \end{pmatrix}$ ($\lambda \ge \gamma$) $\lambda \ge 0$ $\lambda$
$\begin{pmatrix} \lambda & 1 \ 0 & \lambda \end{pmatrix}$ $\lambda < 0, x_1 x_2 \neq 0$ $\max \left{ \frac{-\lambda^2}{|x_2 - \lambda x_1|}, \frac{\lambda}{|x_2|} \right}$
$\begin{pmatrix} \lambda & 1 \ 0 & \lambda \end{pmatrix}$ $\lambda \ge 0$ $\lambda$
$\begin{pmatrix} r\cos\psi & -r\sin\psi \ r\sin\psi & r\cos\psi \end{pmatrix}$ $|\psi| < \frac{\alpha\pi}{2}$ $r \left( \cos\frac{\psi}{\alpha} \right)^\alpha$
$\begin{pmatrix} r\cos\psi & -r\sin\psi \ r\sin\psi & r\cos\psi \end{pmatrix}$ $|\psi| = \frac{\alpha\pi}{2}$ $0$
$\begin{pmatrix} r\cos\psi & -r\sin\psi \ r\sin\psi & r\cos\psi \end{pmatrix}$ $|\psi| > \frac{\alpha\pi}{2}$ $-\frac{r}{\max{|x_1\cos\psi + x_2\sin\psi|, |x_2\cos\psi - x_1\sin\psi|}}$

Chương 3: Tính ổn định của phương trình vi phân phân thứ

Chương 3 khảo sát hệ phi tuyến tự trị $^C D_{0+}^\alpha x(t) = Ax(t) + f(x(t))$ với $f(0)=0$ và thỏa mãn tính chất Lipschitz địa phương suy biến tại gốc: $\lim_{r\to 0} \ell_f(r) = 0$:

  • Chứng minh tính ổn định tiệm cận: Thiết lập định lý: nếu toàn bộ phổ ma trận tuyến tính hóa nằm trong miền ổn định ($\sigma(A) \subset \Lambda_\alpha^s$), tức mọi giá trị riêng $\lambda \in \sigma(A)$ đều thỏa mãn $|\arg(\lambda)| > \frac{\alpha\pi}{2}$, thì điểm cân bằng $x^* = 0$ là ổn định tiệm cận. Để chứng minh, tác giả thiết lập toán tử Lyapunov–Perron: $$\mathcal{T}(\xi)(t) = E_\alpha(t^\alpha A)x_0 + \int_0^t (t-\tau)^{\alpha-1} E_{\alpha,\alpha}((t-\tau)^\alpha A) f(\xi(\tau)) d\tau$$ Chứng minh rằng trên hình cầu đóng bán kính $r$ của không gian Banach $C_\infty(\mathbb{R}{\ge 0}; \mathbb{R}^d)$, toán tử $\mathcal{T}$ là ánh xạ co khi $|x_0|$ và $r$ đủ nhỏ nhờ các đánh giá tích phân tiệm cận của hàm Mittag-Leffler trong miền $\Lambda\alpha^s$.
  • Phân tích sai sót của các công trình trước: Luận án dành mục 3.2.3 thảo luận chi tiết các sai lầm toán học của các tác giả trước:
    • Nhóm Qian và cộng sự [35]: Áp dụng sai bất đẳng thức tích phân phân thứ khi hoán đổi thứ tự lấy tích phân và chuẩn mà không kiểm soát tính kỳ dị của nhân $(t-\tau)^{\alpha-1}$.
    • Nhóm Wen và cộng sự [40]: Nhầm lẫn giữa tính chất của nghiệm hệ phương trình vi phân nguyên với hệ phân thứ khi giả định nghiệm suy giảm mũ thay vì suy giảm lũy thừa ($t^{-\alpha}$).
    • Nhóm Chen và cộng sự [9], Zhang và cộng sự [42]: Giả định tính chất nửa nhóm cho toán tử tiến hóa phân thứ, dẫn đến việc thiết lập các phương trình tích phân sai bản chất.
  • Chứng minh tính không ổn định: Thiết lập định lý khẳng định: nếu phổ ma trận $A$ chứa ít nhất một giá trị riêng nằm trong miền không ổn định ($\sigma(A) \cap \Lambda_\alpha^u \neq \emptyset$, tức tồn tại $\lambda$ có $|\arg(\lambda)| < \frac{\alpha\pi}{2}$), thì điểm cân bằng $x^* = 0$ là không ổn định. Chứng minh được thực hiện bằng cách phân rã không gian theo các không gian con bất biến, xây dựng toán tử tích phân Lyapunov–Perron trên không gian các hàm bị chặn dưới chuẩn có trọng số $|\xi|w := \sup{t\ge 0} e^{-\gamma t} |\xi(t)|$ (với $\gamma > 0$ phù hợp), chứng minh tính co và chỉ ra sự tồn tại của quỹ đạo nghiệm thoát khỏi lân cận của 0. Kết quả này chứng minh trọn vẹn Giả thuyết A–M–M cho hệ $d$-chiều tổng quát.

Chương 4: Đa tạp ổn định của phương trình vi phân phân thứ

Chương 4 mở rộng lý thuyết hệ động lực cho phương trình phân thứ phi tuyến bằng việc chứng minh sự tồn tại của đa tạp ổn định:

  • Thiết lập bài toán: Điểm cân bằng $x^* = 0$ được gọi là hyperbolic nếu ma trận tuyến tính hóa $A$ không có giá trị riêng nào nằm trên đường biên $|\arg(\lambda)| = \frac{\alpha\pi}{2}$, tức $\sigma(A) \subset \Lambda_\alpha^s \cup \Lambda_\alpha^u$. Không gian pha $\mathbb{R}^d$ phân rã thành tổng trực tiếp của không gian con ổn định và không gian con không ổn định: $\mathbb{R}^d = X^s \oplus X^u$.
  • Kỹ thuật giải tích: Tác giả áp dụng phép biến đổi ma trận để đưa $A$ về dạng chuẩn Jordan, sau đó sử dụng kỹ thuật làm nhỏ đường chéo phụ của các khối Jordan: với mỗi $\varepsilon > 0$, tồn tại phép đổi biến tuyến tính biến ma trận khối Jordan $J_k(\lambda)$ thành ma trận có các phần tử trên đường chéo phụ bằng $\varepsilon$. Kỹ thuật này cho phép kiểm soát hoàn toàn các số hạng phi chéo khi đánh giá chuẩn của hàm ma trận Mittag-Leffler $E_\alpha(t^\alpha A)$.
  • Xây dựng toán tử Lyapunov–Perron: Chiếu phương trình lên hai không gian con $X^s$ và $X^u$ thông qua các toán tử chiếu $P^s$ và $P^u$. Toán tử Lyapunov–Perron được xây dựng dưới dạng tích phân hai phần: phần ổn định lấy tích phân tiến từ $0$ đến $t$, phần không ổn định lấy tích phân lùi từ $t$ đến $\infty$: $$\mathcal{T}{x_0^s}(\xi)(t) = E\alpha(t^\alpha A)x_0^s + \int_0^t (t-\tau)^{\alpha-1} P^s E_{\alpha,\alpha}((t-\tau)^\alpha A) f(\xi(\tau)) d\tau - \int_t^\infty (\tau-t)^{\alpha-1} P^u E_{\alpha,\alpha}((\tau-t)^\alpha A) f(\xi(\tau)) d\tau$$
  • Định lý đa tạp ổn định: Chứng minh rằng toán tử $\mathcal{T}{x_0^s}$ là ánh xạ co trên lân cận của không gian Banach. Điểm bất động duy nhất xác định đồ thị của ánh xạ Lipschitz $\psi: X^s \cap B(0, \delta) \to X^u$. Đa tạp ổn định địa phương $W^s(U) = { x_0 \in U : \phi(t, x_0) \in U, \forall t \ge 0 \text{ và } \lim{t\to\infty} \phi(t, x_0) = 0 }$ tồn tại, là một đa tạp Lipschitz tiếp xúc với không gian con ổn định $X^s$ tại gốc tọa độ. Chương kết thúc bằng ví dụ minh họa tính toán đa tạp ổn định cho hệ phi tuyến phẳng hai chiều.

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

Luận án tự xác định các đóng góp mới về mặt lý luận giải tích và thực tiễn phương pháp toán học như sau:

  1. Về lý thuyết số mũ Lyapunov cho hệ phân thứ: Chứng minh được tính chất không âm của số mũ Lyapunov cổ điển đối với mọi nghiệm không tầm thường của phương trình vi phân phân thứ tuyến tính có hệ số liên tục và bị chặn. Phát hiện này giải thích nguyên nhân lý thuyết số mũ cổ điển không thể dùng làm công cụ đo tính ổn định tiệm cận cho hệ phân thứ.
  2. Đề xuất số mũ Lyapunov phân thứ: Định nghĩa thành công khái niệm số mũ Lyapunov phân thứ $\chi_\alpha(f)$ dựa trên ánh xạ ngược của hàm Mittag-Leffler thực $\log_\alpha^M$. Khái niệm này cho phép phản ánh chính xác tốc độ suy giảm theo hàm lũy thừa và đặc trưng hóa hoàn toàn tính ổn định, ổn định tiệm cận của nghiệm qua phổ $\Sigma_\alpha$. Đồng thời, cung cấp thuật toán giải tích tính toán tường minh số mũ này cho hệ hai chiều tổng quát.
  3. Chứng minh hoàn chỉnh định lý ổn định và không ổn định tuyến tính hóa: Thiết lập chứng minh giải tích chuẩn xác đầu tiên cho tính ổn định tiệm cận (khi $\sigma(A) \subset \Lambda_\alpha^s$) và tính không ổn định (khi phổ cắt miền $\Lambda_\alpha^u$) đối với phương trình vi phân phân thứ phi tuyến cấp $\alpha \in (0, 1)$ trong $\mathbb{R}^d$. Khẳng định tính đúng đắn của Giả thuyết A–M–M trên không gian đa chiều, đồng thời vạch rõ và khắc phục các sai lầm toán học trong các công trình quốc tế đã xuất bản của các nhóm tác giả trước đó.
  4. Thiết lập định lý đa tạp ổn định: Xây dựng thành công toán tử Lyapunov–Perron thích hợp kết hợp với kỹ thuật làm nhỏ đường chéo phụ ma trận Jordan, chứng minh sự tồn tại và tính chất trơn Lipschitz của đa tạp ổn định địa phương gần điểm cân bằng hyperbolic cho hệ phân thứ phi tuyến trong không gian hữu hạn chiều tổng quát.

Các kiến nghị và giải pháp khoa học đề xuất: Luận án đề xuất việc nghiên cứu định tính các hệ vi phân phân thứ bắt buộc phải sử dụng các đánh giá hàm Mittag-Leffler trên chu tuyến phức thay cho các bất đẳng thức vi tích phân thông thường, và khuyến nghị áp dụng toán tử kiểu Lyapunov–Perron làm khung phân tích chuẩn tắc cho các bài toán đa tạp bất biến phân thứ.

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

Dựa trên phạm vi và các giả thiết toán học được xác định trong văn bản, luận án có các giới hạn kỹ thuật:

  • Giới hạn cấp phân thứ: Toàn bộ các kết quả chính được thiết lập riêng cho đạo hàm phân thứ Caputo cấp $\alpha \in (0, 1)$. Luận án chưa xét trường hợp cấp phân thứ lớn hơn 1 ($\alpha \ge 1$) hoặc cấp biến thiên theo thời gian.
  • Giới hạn cấu trúc hệ: Các định lý ổn định phi tuyến và đa tạp ổn định ở Chương 3 và Chương 4 chỉ áp dụng cho hệ tự trị (không phụ thuộc tường minh vào thời gian $t$) và hàm phi tuyến $f$ phải thỏa mãn điều kiện triệt tiêu đạo hàm địa phương tại gốc ($\lim_{r\to 0} \ell_f(r) = 0$). Các hệ có nhiễu tuần hoàn, nhiễu ngẫu nhiên hoặc trễ thời gian chưa nằm trong khuôn khổ giải quyết của văn bản.
  • Phạm vi không gian pha: Các kết quả đa tạp ổn định mới dừng lại ở không gian Euclide hữu hạn chiều $\mathbb{R}^d$, chưa mở rộng sang không gian vô hạn chiều (như phương trình đạo hàm riêng phân thứ).
  • Hướng nghiên cứu tiếp được gợi mở từ văn bản: Nghiên cứu sự tồn tại của đa tạp không ổn định ($W^u$), đa tạp trung tâm ($W^c$) quanh điểm cân bằng phi hyperbolic, và mở rộng bài toán số mũ Lyapunov phân thứ cho các hệ phương trình vi phân ngẫu nhiên phân thứ.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên khảo chuyên sâu, có giá trị học thuật đối với:

  • Nghiên cứu sinh và nhà nghiên cứu chuyên ngành phương trình vi phân, hệ động lực: Tham khảo trực tiếp phương pháp luận giải tích hiện đại, đặc biệt là kỹ thuật thiết lập toán tử Lyapunov–Perron trên không gian Banach hàm bị chặn và kỹ thuật biến đổi chuẩn Jordan để xử lý các toán tử không có tính chất nửa nhóm.
  • Giảng viên giảng dạy các chuyên đề toán giải tích cao học: Phần Phụ lục (các đánh giá giải tích phức của hàm Mittag-Leffler hai tham số trên các miền $\Lambda_\alpha^s$ và $\Lambda_\alpha^u$) và Chương 1 (chứng minh công thức biến thiên hằng số qua biến đổi Laplace) cung cấp hệ thống bổ đề hoàn chỉnh, chuẩn xác để đưa vào bài giảng.
  • Các nhà nghiên cứu mô hình hóa trong cơ học vật liệu, vật lý kỹ thuật, lưu biến học: Cung cấp tiêu chuẩn giải tích tường minh (tiêu chuẩn góc phổ ma trận $\frac{\alpha\pi}{2}$ và bảng tra số mũ Lyapunov hai chiều ở Mục 2.3) để kiểm tra tính ổn định tiệm cận của các trạng thái cân bằng trong mô hình vật liệu có trí nhớ hoặc hệ thống điện tử phân thứ.

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

1. Vì sao số mũ Lyapunov cổ điển không thể sử dụng để khảo sát tính ổn định tiệm cận của phương trình vi phân phân thứ?
Trả lời: Theo Định lý 2.2 của luận án, đối với phương trình vi phân phân thứ tuyến tính hệ số liên tục và bị chặn, số mũ Lyapunov cổ điển của mọi nghiệm không tầm thường luôn không âm ($\chi(\phi(\cdot, x_0)) \ge 0$). Do nghiệm của phương trình phân thứ trong miền ổn định suy giảm tiệm cận về 0 theo hàm lũy thừa ($t^{-\alpha}$) chứ không suy giảm nhanh theo hàm mũ âm ($e^{-ct}$), số mũ cổ điển (dựa trên hàm $\log$) bị triệt tiêu về 0, không thể phân biệt được trạng thái nghiệm suy giảm về 0 với trạng thái nghiệm dao động hoặc tăng trưởng chậm hơn hàm mũ.

2. Số mũ Lyapunov phân thứ $\chi_\alpha(f)$ được định nghĩa như thế nào và dựa trên cơ sở toán học nào?
Trả lời: Dựa trên tính chất đơn điệu tăng ngặt của hàm Mittag-Leffler thực $E_\alpha: \mathbb{R} \to \mathbb{R}{\ge 0}$ (với $\lim{x\to -\infty} E_\alpha(x) = 0$ và $\lim_{x\to\infty} E_\alpha(x) = \infty$), tác giả xây dựng hàm ngược liên tục $\log_\alpha^M: \mathbb{R}{>0} \to \mathbb{R}$. Số mũ Lyapunov phân thứ cấp $\alpha$ của hàm $f$ được định nghĩa bằng: $$\chi\alpha(f) = \limsup_{t\to\infty} \frac{1}{t^\alpha} \log_\alpha^M(|f(t)|)$$ Số mũ này nhận giá trị âm khi nghiệm suy giảm tiệm cận về 0 theo quy tắc lũy thừa và nhận giá trị dương khi nghiệm tăng trưởng theo cấp số mũ.

3. Điều kiện cần và đủ trên ma trận tuyến tính hóa để điểm cân bằng của hệ phi tuyến ổn định tiệm cận hoặc không ổn định là gì?
Trả lời: Theo kết quả Chương 3, tính ổn định của hệ phi tuyến $^C D_{0+}^\alpha x(t) = Ax(t) + f(x(t))$ (với $f(0)=0, \lim_{r\to 0}\ell_f(r)=0$) được xác định bởi vị trí phổ $\sigma(A)$ của ma trận $A$ so với góc tới hạn $\frac{\alpha\pi}{2}$:

  • Điểm cân bằng $x^* = 0$ ổn định tiệm cận nếu toàn bộ các giá trị riêng của $A$ nằm trong miền $\Lambda_\alpha^s$, tức $|\arg(\lambda)| > \frac{\alpha\pi}{2}, \forall \lambda \in \sigma(A)$.
  • Điểm cân bằng $x^* = 0$ không ổn định nếu tồn tại ít nhất một giá trị riêng của $A$ nằm trong miền $\Lambda_\alpha^u$, tức $|\arg(\lambda)| < \frac{\alpha\pi}{2}$.

4. Nghiệm của phương trình vi phân phân thứ khác biệt căn bản như thế nào so với phương trình vi phân thường về tính địa phương theo thời gian?
Trả lời: Đối với phương trình vi phân thường, để xác định nghiệm tại thời điểm $t > t_0$, người ta chỉ cần biết thông tin trạng thái địa phương tại điểm hiện tại $t_0$. Ngược lại, biểu thức tích phân tương đương của phương trình phân thứ Caputo tích phân trên toàn bộ khoảng $[0, t]$ với nhân $(t-\tau)^{\alpha-1}$. Do đó, trạng thái nghiệm tại tương lai phụ thuộc đồng thời vào thông tin hiện tại và toàn bộ lịch sử quá khứ của quá trình trên đoạn $[0, t_0)$, thể hiện bản chất có bộ nhớ và tính chất phi địa phương theo thời gian.

5. Vì sao không thể áp dụng trực tiếp các phương pháp của hệ động lực cổ điển cho phương trình vi phân phân thứ?
Trả lời: Lý do căn bản là các toán tử đạo hàm và tích phân phân thứ không thỏa mãn tính chất nửa nhóm ($D_{a+}^{\alpha_1} D_{a+}^{\alpha_2} \neq D_{a+}^{\alpha_1+\alpha_2}$). Do đó, toán tử giải của phương trình vi phân phân thứ không sinh ra một nhóm hoặc nửa nhóm chuyển pha một tham số ${\Phi_t}{t\ge 0}$ thỏa mãn $\Phi{t+s} = \Phi_t \circ \Phi_s$, làm vô hiệu hóa các công cụ cổ điển vốn dựa trên tính chất nửa nhóm của hệ động lực.

Kết luận

Luận án tiến sĩ của tác giả Hoàng Thế Tuấn là một công trình chuyên khảo giải tích toán học chặt chẽ, giải quyết trọn vẹn các bài toán định tính cơ bản của phương trình vi phân phân thứ Caputo cấp $\alpha \in (0, 1)$ trên không gian hữu hạn chiều. Tác giả đã chứng minh giới hạn của số mũ Lyapunov cổ điển, đề xuất thành công khái niệm số mũ Lyapunov phân thứ để đặc trưng hóa tính ổn định của hệ tuyến tính, đồng thời cung cấp bảng tính toán tường minh cho hệ hai chiều. Bằng phương pháp toán tử Lyapunov–Perron kết hợp kỹ thuật phân rã Jordan và tiệm cận giải tích phức của hàm Mittag-Leffler, luận án đã chứng minh chuẩn xác định lý ổn định tiệm cận, tính không ổn định thông qua tuyến tính hóa và xác lập sự tồn tại của đa tạp ổn định địa phương gần điểm cân bằng hyperbolic.