Tổng quan luận án

Luận án tiến sĩ thuộc chuyên ngành Toán giải tích (chuyên sâu về lý thuyết định tính phương trình vi phân và ứng dụng), được thực hiện tại Bộ môn Giải tích, Khoa Toán - Tin, Trường Đại học Sư phạm Hà Nội dưới sự hướng dẫn khoa học của PGS.TS. Lê Văn Hiện và TS. Trịnh Tuấn Anh.

Tính cấp thiết và khoảng trống nghiên cứu

Lý thuyết định tính phương trình vi phân và lý thuyết ổn định nghiệm là hướng nghiên cứu trọng tâm trong lý thuyết điều khiển hệ thống, xuất hiện trong nhiều bài toán thực tiễn từ công nghệ thông tin, trí tuệ nhân tạo đến sinh thái học quần thể. Trong thực tế, các hệ động lực thường chứa yếu tố trễ truyền tải tín hiệu hoặc trễ sinh học. Đối với các hệ dừng có trễ hằng số, phương pháp hàm Lyapunov và bất đẳng thức ma trận tuyến tính (LMI) là công cụ phổ biến. Tuy nhiên, nhiều hệ thống thực tế có cấu trúc phi tuyến không dừng phức tạp:

  • Trong các mạng nơron nhân tạo nhiều tầng, thời gian trễ truyền tín hiệu thường tỷ lệ với thời gian hiện tại ($q_{ij}t$ với $0 < q_{ij} < 1$). Đại lượng này tương ứng với hàm trễ biến thiên không bị chặn $\tau_{ij}(t) = (1 - q_{ij})t \to \infty$ khi $t \to \infty$. Các phương pháp truyền thống (như đổi biến thời gian $t = e^\xi$) chỉ hiệu quả với hệ dừng và không áp dụng được cho hệ có hệ số biến thiên theo thời gian.
  • Khái niệm ổn định theo Lyapunov (LS) chỉ phản ánh dáng điệu tiệm cận của nghiệm tại vô hạn, trong khi các bài toán kỹ thuật thực tế đòi hỏi kiểm soát quỹ đạo trạng thái trong một khoảng thời gian hữu hạn xác định trước (khái niệm ổn định trong thời gian hữu hạn - FTS). LS và FTS là hai khái niệm độc lập. Trước công trình của tác giả, chưa có nghiên cứu nào công bố về tính ổn định hữu hạn của mạng nơron Hopfield với hệ số biến thiên và trễ tỉ lệ.
  • Vấn đề phân tích tính tiêu hao (dissipativity) của mạng nơron không dừng chứa trễ tỉ lệ khi hệ số tự phản hồi suy biến ($\inf_{t \ge 0} a_i(t) = 0$) vẫn là một thách thức toán học chưa được giải quyết bằng bất đẳng thức Halanay thông thường.
  • Trong sinh thái học quần thể, mô hình Nicholson truyền thống thường chỉ xét tốc độ suy thoái dân số tuyến tính (chỉ phù hợp với mật độ thấp). Đối với các quần thể mật độ cao hoặc mô hình bảo tồn biển, tốc độ suy giảm cá thể phụ thuộc phi tuyến vào mật độ (chẳng hạn dạng $D(t, N) = a(t) - b(t)e^{-N}$). Dáng điệu tiệm cận và tính hút của nghiệm tuần hoàn đối với lớp phương trình này chưa được khảo sát đầy đủ.

Mục tiêu nghiên cứu

Luận án xác định 4 mục tiêu nghiên cứu cụ thể:

  1. Thiết lập các điều kiện đủ thông qua tính chất phổ của $M$-ma trận đảm bảo tính ổn định trong thời gian hữu hạn và tính đồng bộ nghiệm với tốc độ lũy thừa của mô hình mạng nơron Hopfield với hệ số biến thiên và trễ tỉ lệ.
  2. Phân tích tính tiêu hao toàn cục và xác định cấu trúc tập hấp thụ dạng mũ suy rộng cho lớp hệ phương trình vi phân mô tả mạng nơron Hopfield không dừng chứa trễ tỉ lệ trong cả hai trường hợp hệ số tự phản hồi chính quy và hệ số tự phản hồi suy biến.
  3. Khảo sát tính tồn tại toàn cục, tính bền vững đều và tính tiêu hao đều của nghiệm dương đối với mô hình Nicholson phi tuyến chứa đa trễ biến thiên và hàm suy giảm phi tuyến.
  4. Thiết lập điều kiện đảm bảo sự tồn tại duy nhất của nghiệm tuần hoàn dương có tính chất hút toàn cục đối với mô hình Nicholson có trễ phi tuyến, đồng thời mở rộng cho tính hút toàn cục của điểm cân bằng dương trong trường hợp hệ số hằng số.

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

  • Đối tượng nghiên cứu: Các lớp phương trình vi phân phi tuyến có trễ, bao gồm:
    • Hệ phương trình vi phân mô tả mạng nơron Hopfield không dừng với đa trễ tỉ lệ không đồng nhất.
    • Phương trình vi phân phi tuyến mô tả mô hình động lực học quần thể Nicholson với trễ biến thiên và tốc độ suy giảm mật độ phi tuyến.
  • Phạm vi không gian và thời gian: Không gian Euclide hữu hạn chiều $\mathbb{R}^n$ (chuẩn cực đại $|\cdot|_\infty$) và không gian hàm liên tục $C([a, b], \mathbb{R}^n)$; khảo sát trên khoảng thời gian hữu hạn $[0, T]$ (đối với bài toán ổn định hữu hạn) và khoảng thời gian vô hạn $[0, \infty)$ (đối với bài toán dáng điệu tiệm cận, tính tiêu hao và tính hút toàn cục).

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 trong văn bản

Luận án điểm lại các nhánh nghiên cứu lý thuyết và thực nghiệm liên quan trực tiếp đến đề tài:

Nhóm vấn đề Các tác giả và công trình được trích dẫn Nội dung trọng tâm đã công bố
Lý thuyết ổn định và PTVP có trễ Hale & Verduyn Lunel [12], Kolmanovskii & Myshkis [21], Smith [31], Wu [38], Lakshmikantham et al. [1] Nền tảng phương trình vi phân hàm, phương pháp hàm Lyapunov-Krasovskii và dáng điệu tiệm cận của nghiệm.
Ổn định thời gian hữu hạn (FTS) Kamenkov (1953) [20], Amato et al. (2001) [2], Hien & Trinh (2016) [16] Nguồn gốc FTS trong cơ học điều khiển; phân tích tính độc lập giữa FTS và ổn định Lyapunov; thiết lập điều kiện LMI cho hệ tuyến tính có trễ biến thiên.
Mạng nơron và trễ tỉ lệ Hopfield [28], Roska & Chua [34], Liao & Wang [40], Zhou [41], Chen et al. [42], Xu et al. [44], Zhang et al. [45], Berezansky et al. [8] Động lực học mạng nơron; trễ tỉ lệ trong truyền tín hiệu nhiều tầng; phương pháp đổi biến $t = e^\xi$ cho hệ dừng; điều kiện ổn định đa thức cho hệ hằng số.
Tính tiêu hao và BĐT Halanay Hale (1988) [32], Halanay (1966) [19], Wen et al. [36], Liu [37], Zhou et al. (2014) [43] Khái niệm tập hấp thụ và tính tiêu hao toàn cục; ứng dụng bất đẳng thức Halanay và cải biên; phương pháp đổi biến chuyển trễ tỉ lệ về trễ hằng cho hệ dừng.
Mô hình Nicholson có trễ Nicholson (1954) [29], Berezansky et al. (2010, 2011) [3, 4], Faria (2017) [6], Liu [22], Saker & Agarwal [26], Berezansky et al. [10], Wang [35] Mô hình sinh trưởng ve châu Úc; dáng điệu tiệm cận của nghiệm tuần hoàn; mô hình có số hạng thu hoạch; khảo sát mô hình với tốc độ suy giảm phi tuyến Type-I ($a - be^{-N}$) và Type-II ($\frac{aN}{b+N}$).

Khoảng trống nghiên cứu được luận án giải quyết

  1. Các kết quả ổn định và đồng bộ hóa trước đây (như trong [42, 44, 45]) chỉ áp dụng cho mạng nơron có trọng số kết nối hằng số và trễ bị chặn, hoặc dựa vào phép đổi biến trạng thái chỉ dùng được cho hệ dừng. Luận án giải quyết trường hợp hệ số biến thiên theo thời gian chứa trễ tỉ lệ không đồng nhất và không bị chặn.
  2. Phương pháp hàm Lyapunov toàn phương trong [42] đưa ra các điều kiện ngặt. Luận án sử dụng kỹ thuật so sánh kết hợp $M$-ma trận để nới lỏng các điều kiện này.
  3. Bất đẳng thức Halanay kinh điển yêu cầu hệ số tiêu hao phải âm đều ($\alpha(t) \le -\sigma < 0$). Luận án xây dựng kỹ thuật đánh giá cho phép hệ số tự phản hồi suy biến về 0 trên vô hạn ($\inf_{t \ge 0} a_i(t) = 0$).
  4. Các nghiên cứu mô hình Nicholson phần lớn giả định tỉ lệ tử vong tuyến tính. Luận án thiết lập kết quả hoàn chỉnh về sự tồn tại duy nhất và tính hút toàn cục của nghiệm tuần hoàn cho mô hình có tốc độ suy giảm phụ thuộc phi tuyến vào mật độ ($D(t, N) = a(t) - b(t)e^{-N}$).

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

Cơ sở lý thuyết và các khái niệm sử dụng

  • Lý thuyết $M$-ma trận: Một ma trận $A = (a_{ij}) \in \mathbb{R}^{n \times n}$ là $M$-ma trận nếu $a_{ij} \le 0$ với mọi $i \ne j$ và các định thức con chính đều dương. Luận án sử dụng các tính chất tương đương:
    • $\text{Re}(\lambda_j) > 0$ với mọi giá trị riêng $\lambda_j \in \lambda(A)$.
    • Tồn tại ma trận không âm $B \ge 0$ và số thực $s > \rho(B)$ sao cho $A = sI_n - B$.
    • Tồn tại vectơ $\xi \in \mathbb{R}^n, \xi \succ 0$ sao cho $A\xi \succ 0$.
    • Nếu $A$ là $M$-ma trận không suy biến, tồn tại $\chi \in \mathbb{R}^n, \chi \succ 0$ với $|\chi|_\infty = 1$ sao cho $A\chi \succ 0$.
  • Lý thuyết ổn định Lyapunov và ổn định thời gian hữu hạn:
    • Ổn định đều và ổn định tiệm cận toàn cục: Định nghĩa thông qua hàm ban đầu $\phi \in C([ -r, 0], \mathbb{R}^n)$ và chuẩn $|\phi|C = \sup{-r \le s \le 0} |\phi(s)|$.
    • Ổn định mũ toàn cục (GES): $|x(t, \phi)| \le \beta |\phi|_C e^{-\alpha(t - t_0)}$.
    • Ổn định trong thời gian hữu hạn (FTS): Với $T > 0, r_1 < r_2$, hệ là FTS đối với $(r_1, r_2, T)$ nếu xuất phát từ $|\phi|C \le r_1$ thì quỹ đạo thỏa mãn $|x(t)|\infty < r_2$ với mọi $t \in [0, T]$.
  • Khái niệm về tính tiêu hao:
    • Tập hút toàn cục: Tập compact $\Omega \subset \mathbb{R}^n$ thỏa mãn $\lim\sup_{t \to \infty} \rho(x(t), \Omega) = 0$.
    • Tập hút mũ suy rộng: $\rho(x(t), \Omega) \le \kappa(|x_0|\infty) e^{-\sigma(t)}$ với $\lim{t \to \infty} \sigma(t) = \infty$.
    • Tính tiêu hao toàn cục: Tồn tại tập hấp thụ bị chặn $B \subset \mathbb{R}^n$ sao cho mọi quỹ đạo đi vào và nằm trong $B$ sau một thời gian hữu hạn.
  • Công cụ giải tích vi phân:
    • Đạo hàm Dini trên bên phải: $D^+ v(t) = \limsup_{h \to 0^+} \frac{v(t+h) - v(t)}{h}$.
    • Hàm dấu suy rộng: $\sigma_f(t) = 1$ nếu $f(t) > 0$ (hoặc $f(t) = 0, f'(t) > 0$); $\sigma_f(t) = 0$ nếu $f(t) = 0, f'(t) = 0$; $\sigma_f(t) = -1$ nếu $f(t) < 0$ (hoặc $f(t) = 0, f'(t) < 0$). Khi đó $D^+ |f(t)| = f'(t)\sigma_f(t)$.

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

Luận án kết hợp các phương pháp toán học giải tích:

  1. Nguyên lý so sánh và bất đẳng thức vi - tích phân: Xây dựng hệ hàm so sánh dạng mũ hoặc dạng lũy thừa, sử dụng đạo hàm Dini và lập luận phản chứng để đánh giá chặn trên cho sai sai phân các nghiệm.
  2. Kỹ thuật phổ ma trận: Chuyển các hệ điều kiện vi phân phi tuyến phức tạp về bài toán kiểm tra tính không suy biến của $M$-ma trận thông qua việc tồn tại vectơ dương riêng $\xi \succ 0$.
  3. Phương pháp phiếm hàm Lyapunov-Krasovskii: Sử dụng cho các hệ tuyến tính có trễ hỗn hợp biến thiên để đưa về bất đẳng thức ma trận tuyến tính (LMI).
  4. Mô phỏng số: Thiết lập thuật toán tính toán tốc độ đồng bộ cực đại $\sigma_{\max}$, giải tích số trên các ví dụ ma trận cấp 2 và vẽ quỹ đạo pha trạng thái để kiểm chứng lý thuyết.

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

Chương 1: Sơ bộ một số kết quả liên quan

Chương 1 hệ thống hóa các kiến thức nền tảng và kết quả bổ trợ phục vụ cho các chương sau:

  • Lý thuyết $M$-ma trận: Trình bày 5 điều kiện tương đương đặc trưng cho $M$-ma trận không suy biến (Mệnh đề 1.1) và tính chất chuẩn hóa vectơ dương $\chi \succ 0, |\chi|_\infty = 1$.
  • Hệ phương trình vi phân có trễ: Trình bày bài toán Cauchy $x'(t) = f(t, x_t)$, các định lý về sự tồn tại và duy nhất nghiệm địa phương/toàn cục dưới điều kiện Lipschitz, điều kiện tăng trưởng tuyến tính $|f(t, \phi)| \le a(t)|\phi|_C + b(t)$, và điều kiện dạng Nagumo. Trình bày các định lý ổn định nghiệm theo hàm Lyapunov $V(t, \phi)$.
  • Tính ổn định trong thời gian hữu hạn (FTS):
    • Định nghĩa hình học của FTS trên tập trạng thái ban đầu $X_0$ và tập quỹ đạo $X_t$ (hoặc dạng ellipsoid $E_R(\rho) = {x \in \mathbb{R}^n: x^\top R x < \rho}$).
    • Phân tích mối quan hệ giữa FTS và ổn định Lyapunov qua hai ví dụ phản chứng:
      1. Phương trình $x'(t) = -1.1 \frac{t+2}{t+1} x(t)$ ổn định tiệm cận theo Lyapunov nhưng không ổn định hữu hạn đối với $(r_1 = 1, r_2 = 1.05, T = 10)$.
      2. Phương trình $x'(t) = -0.1 \frac{t}{t+6} x(t)$ ổn định hữu hạn đối với $(r_1 = 1, r_2 = 1.5, T = 10)$ nhưng không ổn định tiệm cận theo Lyapunov.
    • Trình bày kết quả FTS cho hệ tuyến tính với trễ hỗn hợp biến thiên $x'(t) = Ax(t) + Dx(t - \tau(t)) + G \int_{t - \kappa(t)}^t x(s)ds$ thông qua điều kiện ma trận LMI $\Pi = \Pi_0 + \Pi_1 + \Pi_2 < 0$ (Định lý 1.1).
  • Tính tiêu hao và Bất đẳng thức Halanay: Trình bày nguyên lý thiết lập tập hấp thụ $B(0, \sqrt{-\gamma^*/(\alpha_0 + \beta_0)} + \epsilon)$ qua bất đẳng thức Halanay. Phân tích phương pháp đổi biến $y(t) = x(e^t)$ của Zhou et al. [43] và chỉ ra hạn chế không thể mở rộng cho hệ có hệ số biến thiên.
  • Bổ đề bổ trợ: Trình bày bất đẳng thức ma trận $2x^\top y \le x^\top M x + y^\top M^{-1} y$ và bổ đề về sự tồn tại của hàm tuần hoàn tiệm cận từ tích phân $\int_0^\infty |u(t + \omega) - u(t)|dt < \infty$.

Chương 2: Tính ổn định hữu hạn của lớp hệ phương trình vi phân mô tả mạng nơron Hopfield với hệ số biến thiên và trễ tỉ lệ

Chương này nghiên cứu mô hình mạng nơron Hopfield với hệ số biến thiên và đa trễ tỉ lệ không đồng nhất: $$x'i(t) = -a_i(t)x_i(t) + \sum{j=1}^n b_{ij}(t)f_j(x_j(t)) + \sum_{j=1}^n c_{ij}(t)g_j(x_j(q_{ij}t)) + I_i(t), \quad t > 0, \quad i \in [n]$$ với $q_{ij} \in (0, 1)$ và điều kiện đầu $x_i(0) = x_{0i}$.

  • Các giả thiết kỹ thuật:
    • (A2.1): Hàm kích hoạt thỏa mãn điều kiện Lipschitz suy rộng: $l_{ik}^- \le \frac{f_i(x) - f_i(y)}{x - y} \le l_{ik}^+$, đặt $L_{fi} = \max{|l_{i1}^+|, |l_{i1}^-|}, L_{gi} = \max{|l_{i2}^+|, |l_{i2}^-|}$.
    • (A2.2): $A(t) \ge A = \text{diag}{a_1, \dots, a_n} > 0$, $|b_{ij}(t)| \le b_{ij}$, $|c_{ij}(t)| \le c_{ij}$. Đặt $M = B L_f + C L_g - A$.
  • Kết quả chính về Ổn định hữu hạn (Định lý 2.1): Cho trước $0 < r_1 < r_2$ và $T > 0$. Hệ mạng nơron ổn định hữu hạn đối với $(r_1, r_2, T)$ nếu tồn tại số dương $\gamma > 0$ và vectơ $\xi \in \mathbb{R}^n, \xi \succ 0$ thỏa mãn:
    1. $(M - \gamma I)\xi \prec 0$,
    2. $e^{\gamma T} < \frac{r_2}{C(\xi) r_1}$, trong đó $C(\xi) = \frac{\xi^+}{\xi^-} = \frac{\max_{i} \xi_i}{\min_{i} \xi_i}$.
  • Kết quả về Dáng điệu tiệm cận và Sự đồng bộ nghiệm (Định lý 2.2): Nếu ma trận $-M = A - B L_f - C L_g$ là một $M$-ma trận không suy biến, thì hai nghiệm bất kỳ $x(t)$ và $x^(t)$ của hệ sẽ đồng bộ hóa với tốc độ lũy thừa: $$|x(t) - x^(t)|\infty \le \beta \frac{|x(0) - x^*(0)|\infty}{(1 + t)^{\sigma_0}}, \quad \forall t \ge 0$$ với $\beta = C(\xi)$ và hằng số $\sigma_0 > 0$ được xác định duy nhất từ phương trình $H_i(\sigma) = 0$.
  • Thuật toán xác định tốc độ đồng bộ $\sigma_{\max}$:
    1. Tìm $\xi \succ 0$ sao cho $M\xi \prec 0$.
    2. Tính $\eta = \min_{i \in [n]} { a_i \xi_i - \sum_{j=1}^n \xi_j (b_{ij} L_{fj} + c_{ij} L_{gj}) }$.
    3. Giải phương trình $\sigma \xi_i + \sum_{j=1}^n L_{gj} c_{ij} \xi_j [e^{\sigma \ln(1/q_{ij})} - 1] = \eta$ tìm nghiệm $\sigma_i^$, sau đó lấy $\sigma_{\max} = \min_{i \in [n]} \sigma_i^$.
  • Ví dụ và so sánh học thuật:
    • Ví dụ 1 (Cấp 2): Hệ có $M = \begin{pmatrix} -2 & 2 \ 2 & -2 \end{pmatrix}$. Do $\det(M) = 0$ nên $-M$ suy biến (không ổn định tiệm cận). Tuy nhiên với $\gamma = 0.005$, ma trận $\gamma I - M$ là $M$-ma trận không suy biến. Với $r_1 = 1, r_2 = 1.01$, hệ đạt FTS trên đoạn $[0, 15]$ với $T_{\max} = 18$.
    • Ví dụ 2: Phương trình $x'(t) = -(1 + e^{-t})x(t) + (1 - e^{-t})x(0.5t)$ thỏa mãn điều kiện FTS với mọi $r_1, r_2, T > 0$, nhưng quỹ đạo tiến về $x(0)e^{-2} > 0$, chứng minh hệ không ổn định tiệm cận theo Lyapunov.
    • So sánh với [42, 44, 45]: Luận án chứng minh điều kiện của Chen et al. [42] đòi hỏi bán kính phổ $\rho(K) < 1$ (bị vi phạm trong Ví dụ 3), và điều kiện LMI trong [44, 45] đòi hỏi $q > \frac{1}{(\gamma - 1)^2}$ (rất ngặt so với điều kiện của luận án chỉ cần $\gamma > 2$ với mọi $q \in (0, 1)$).

Chương 3: Tính tiêu hao của lớp hệ phương trình vi phân mô tả mạng nơron không dừng chứa trễ tỉ lệ

Chương này mở rộng nghiên cứu sang tính tiêu hao toàn cục và xác định kích thước tập hấp thụ cho hệ: $$x'i(t) = -a_i(t)x_i(t) + \sum{j=1}^n b_{ij}(t)f_j(x_j(t)) + \sum_{j=1}^n c_{ij}(t)g_j(x_j(q_{ij}t)) + I_i(t), \quad t > 0$$

  • Thiết lập sơ bộ và giả thiết:
    • (A3.1): Hàm kích hoạt thỏa mãn Lipschitz: $|f_j(a) - f_j(b)| \le L_{fj}|a - b|$, $|g_j(a) - g_j(b)| \le L_{gj}|a - b|$.
    • (A3.2): Trọng số và đầu vào bị chặn: $|b_{ij}(t)| \le b_{ij}$, $|c_{ij}(t)| \le c_{ij}$, $|I_i(t)| \le \bar{I}_i$.
    • Đặt $\gamma = \max_{i \in [n]} { \sum_{j=1}^n (b_{ij}|f_j(0)| + c_{ij}|g_j(0)|) + \bar{I}_i }$.
  • Trường hợp hệ số tự phản hồi chính quy ($a_i(t) \ge a_i > 0$):
    • Định lý 3.1: Giả sử $M = A - B L_f - C L_g$ là một $M$-ma trận không suy biến. Khi đó:
      1. Tập $\Omega = \left{ x \in \mathbb{R}^n : |x|\infty \le \frac{\gamma}{(M\chi)+} \right}$ là tập hút mũ suy rộng của hệ, với $\chi \succ 0, |\chi|_\infty = 1$ sao cho $M\chi \succ 0$.
      2. Hệ tiêu hao toàn cục, và hình cầu đóng $B\left(0, \frac{\gamma}{(M\chi)_+} + \epsilon\right)$ là tập hấp thụ với mọi $\epsilon > 0$.
    • Đánh giá quỹ đạo nghiệm thỏa mãn: $$|x(t)|\infty \le \frac{\gamma}{m^*} + \kappa{x_0} e^{-\hat{\alpha}\ln(1+t)}, \quad \forall t \ge 0$$ trong đó $m^* = \min_{i \in [n]} { a_i \chi_i - \sum_{j=1}^n (b_{ij} L_{fj} + c_{ij} L_{gj})\chi_j }$.
    • Hệ quả 3.1: Nếu $a_i > \sum_{j=1}^n (L_{fi} b_{ji} + L_{gi} c_{ji})$ với mọi $i \in [n]$, hệ tiêu hao toàn cục. Kết quả này nới lỏng trực tiếp điều kiện trong [37] vốn đòi hỏi $a_i > \sum_{j=1}^n L_{fi} b_{ji} + \frac{L_{gi}}{1 - \tau_M} \sum_{j=1}^n c_{ji}$.
  • Trường hợp hệ số tự phản hồi suy biến:
    • (A3.4): Tồn tại hàm $\phi(t) > 0$ và các hằng số $\hat{a}i > 0$ sao cho $a_i(t) \ge \hat{a}i \phi(t)$, $\sup{t \ge 0} \int_0^t \phi(s)ds < \infty$ hoặc $\lim{t \to \infty} \int_0^t \phi(s)ds = \infty$.
    • Bằng cách xây dựng hàm so sánh tích phân suy rộng, luận án chứng minh tính tiêu hao toàn cục mà không cần giả thiết $\inf_{t \ge 0} a_i(t) > 0$.

Chương 4: Tính hút toàn cục của nghiệm tuần hoàn dương của một mô hình Nicholson có trễ

Chương này nghiên cứu mô hình động lực học quần thể Nicholson có trễ với tốc độ suy giảm phi tuyến Type-I: $$N'(t) = -D(t, N(t)) + \sum_{k=1}^p \beta_k(t) N(t - \tau_k(t)) e^{-\gamma_k(t) N(t - \tau_k(t))}, \quad t \ge 0$$ với hàm tử vong phi tuyến $D(t, N) = a(t) - b(t)e^{-N}$, trong đó $a(t), b(t), \beta_k(t), \gamma_k(t), \tau_k(t)$ là các hàm $\omega$-tuần hoàn, liên tục và dương.

  • Nghiệm dương toàn cục và tính bền vững đều (Uniform Persistence):
    • Luận án chứng minh rằng với mọi điều kiện đầu liên tục dương $\phi(t) > 0, t \in [-\tau, 0]$, nghiệm $N(t, \phi)$ tồn tại duy nhất và dương trên toàn bộ $[0, \infty)$.
    • Sử dụng phương pháp so sánh vi-tích phân, luận án chỉ ra tính bị chặn trên (tiêu hao đều) và bị chặn dưới tách rời điểm 0 (bền vững đều): $$0 < m_* \le \liminf_{t \to \infty} N(t) \le \limsup_{t \to \infty} N(t) \le M^* < \infty$$
  • Sự tồn tại và tính hút toàn cục của nghiệm tuần hoàn dương:
    • Bằng nguyên lý ánh xạ co hoặc định lý điểm bất động kết hợp với tính tiêu hao đều, luận án chứng minh phương trình tồn tại duy nhất một nghiệm $\omega$-tuần hoàn dương $N^*(t)$.
    • Mọi nghiệm dương $N(t)$ khác đều bị hút về $N^(t)$ khi $t \to \infty$: $$\lim_{t \to \infty} |N(t) - N^(t)| = 0$$
  • Áp dụng cho mô hình hệ số hằng số: Khi các tham số $a, b, \beta_k, \gamma_k, \tau_k$ là hằng số, nghiệm tuần hoàn suy biến thành điểm cân bằng dương $N^* > 0$. Luận án thiết lập điều kiện đại số tường minh đảm bảo điểm cân bằng dương $N^*$ tồn tại duy nhất và có tính hút toàn cục trên miền nghiệm dương.

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

Luận án tự tổng kết các đóng góp học thuật mới về mặt lý thuyết và ứng dụng:

  1. Về lý thuyết ổn định thời gian hữu hạn và đồng bộ mạng nơron:
    • Thiết lập được các điều kiện đủ mới thông qua tính chất phổ của $M$-ma trận để đảm bảo tính ổn định trong thời gian hữu hạn (FTS) cho mạng nơron Hopfield không dừng chứa đa trễ tỉ lệ không đồng nhất.
    • Chứng minh được rằng khi $-M$ là $M$-ma trận không suy biến, các quỹ đạo trạng thái của mạng nơron sẽ đồng bộ hóa với tốc độ lũy thừa $\mathcal{O}((1+t)^{-\sigma_0})$. Đưa ra thuật toán giải tích xác định tốc độ đồng bộ tối ưu $\sigma_{\max}$.
  2. Về lý thuyết tiêu hao của hệ động lực không dừng:
    • Giải quyết bài toán phân tích tính tiêu hao toàn cục cho lớp mạng nơron chứa trễ tỉ lệ không bị chặn mà không cần sử dụng phương pháp đổi biến thời gian.
    • Xác định tường minh bán kính của tập hút mũ suy rộng $\Omega$ và tập hấp thụ $B$ trong cả hai trường hợp: hệ số tự phản hồi chính quy ($a_i(t) \ge a_i > 0$) và hệ số tự phản hồi suy biến ($\inf_{t \ge 0} a_i(t) = 0$).
  3. Về lý thuyết định tính mô hình sinh thái Nicholson phi tuyến:
    • Chứng minh tính tồn tại toàn cục, tính bền vững đều và tính tiêu hao đều của nghiệm dương đối với mô hình Nicholson có trễ với tốc độ suy giảm mật độ phi tuyến dạng $D(t, N) = a(t) - b(t)e^{-N}$.
    • Thiết lập các điều kiện đủ đảm bảo sự tồn tại, tính duy nhất và tính hút toàn cục của nghiệm tuần hoàn dương (đối với hệ tuần hoàn) và của điểm cân bằng dương (đối với hệ dừng).
  4. Công bố khoa học:
    • Các kết quả chính của luận án đã được công bố trong 03 công trình khoa học trên các tạp chí quốc tế chuyên ngành thuộc danh mục ISI.
    • Các nội dung đã được báo cáo và thảo luận tại: Xemina Bộ môn Giải tích (ĐH Sư phạm Hà Nội), Hội nghị khoa học Nghiên cứu sinh (Khoa Toán - Tin, ĐH Sư phạm Hà Nội) và Xemina Phòng Tối ưu và Điều khiển (Viện Toán học, Viện Hàn lâm Khoa học và Công nghệ Việt Nam).

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

Dựa trên phạm vi và nội dung nghiên cứu thực tế được trình bày trong văn bản:

  • Hạn chế:
    • Các kết quả về mạng nơron tập trung chủ yếu vào mô hình dạng Hopfield với trễ tỉ lệ dạng đại số $q_{ij}t$; chưa xét đến các cấu trúc mạng phức tạp hơn như mạng nơron ngẫu nhiên, mạng nơron bậc cao hoặc mạng có chứa trễ trung hòa.
    • Đối với mô hình Nicholson, luận án mới chỉ khảo sát dạng hàm suy giảm mật độ phi tuyến Type-I ($D(t, N) = a(t) - b(t)e^{-N}$); chưa mở rộng cho lớp hàm suy giảm Type-II ($D(N) = \frac{aN}{b+N}$) hoặc mô hình mạng lưới nhiều loài tương tác với trễ phân bố.
  • Hướng nghiên cứu tiếp:
    • Tiếp tục phát triển phương pháp so sánh $M$-ma trận cho các lớp phương trình vi phân hàm phi tuyến không dừng chứa đồng thời trễ tỉ lệ và trễ phân phối ngẫu nhiên.
    • Mở rộng khảo sát tính ổn định và tính hút toàn cục cho mô hình sinh thái với hàm suy thoái mật độ phi tuyến Type-II và các hệ động lực sinh học tương tác đa quần thể.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên khảo trong lĩnh vực giải tích định tính phương trình vi phân:

  • Đối tượng thụ hưởng:
    • Nghiên cứu sinh, học viên cao học và giảng viên ngành Toán giải tích, Toán ứng dụng: Tham khảo kỹ thuật áp dụng lý thuyết $M$-ma trận kết hợp đạo hàm Dini để giải quyết các phương trình vi phân không dừng chứa trễ không bị chặn.
    • Các nhà nghiên cứu trong lĩnh vực Lý thuyết điều khiển và Trí tuệ nhân tạo: Tham khảo các tiêu chuẩn ổn định thời gian hữu hạn (FTS), thuật toán ước lượng tốc độ đồng bộ lũy thừa và phương pháp ước lượng tập hấp thụ trong thiết kế mạng nơron Hopfield.
    • Các nhà nghiên cứu Toán sinh học và Động lực học quần thể: Tham khảo phương pháp đánh giá tính bền vững đều và tính hút của nghiệm tuần hoàn cho mô hình Nicholson với hàm tử vong phi tuyến.
  • Phần đáng tham khảo nhất:
    • Kỹ thuật chứng minh Định lý 2.1 và thuật toán tìm tốc độ đồng bộ $\sigma_{\max}$ trong Chương 2.
    • Kỹ thuật chặn nghiệm xây dựng tập hút mũ suy rộng cho hệ suy biến trong Chương 3.
    • Khung phân tích tính bền vững và nghiệm tuần hoàn dương của mô hình Nicholson trong Chương 4.

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

1. Tính ổn định trong thời gian hữu hạn (FTS) và tính ổn định theo Lyapunov (LS) khác nhau như thế nào?

Tính ổn định theo Lyapunov là khái niệm định tính, xác định dáng điệu tiệm cận của nghiệm khi thời gian dần ra vô hạn ($t \to \infty$). Ngược lại, tính ổn định trong thời gian hữu hạn mang tính định lượng, đòi hỏi mọi quỹ đạo xuất phát từ một miền cho trước $|x(0)| \le r_1$ không được vượt quá ngưỡng giới hạn $|x(t)| < r_2$ trong suốt khoảng thời gian định trước $[0, T]$. Hai khái niệm này độc lập: một hệ có thể đạt FTS nhưng không ổn định Lyapunov, và ngược lại, một hệ ổn định tiệm cận Lyapunov vẫn có thể vi phạm ngưỡng FTS trong khoảng thời gian ngắn.

2. Tại sao phương pháp đổi biến thời gian truyền thống không áp dụng được cho mô hình trong Chương 2 và Chương 3?

Với các hệ dừng (hệ số hằng số) có trễ tỉ lệ $x'(t) = g(x(t), x(qt))$, phép đổi biến trạng thái $t = e^\xi$ sẽ chuyển trễ tỉ lệ $qt$ thành trễ hằng số $\xi - \tau$ (với $\tau = -\ln q$). Tuy nhiên, đối với hệ có hệ số biến thiên theo thời gian $a_i(t), b_{ij}(t), c_{ij}(t)$, phép đổi biến này sẽ làm biến đổi ma trận hệ số thành dạng phức tạp phụ thuộc biến mới $e^\xi$, phá vỡ cấu trúc giải tích của hệ và không thể áp dụng các định lý Lyapunov truyền thống.

3. Điều kiện ổn định dựa trên M-ma trận của luận án có ưu điểm gì so với các công bố của Chen et al. [42] hay Zhang et al. [45]?

Điều kiện trong [42] sử dụng hàm Lyapunov toàn phương, dẫn đến yêu cầu bán kính phổ $\rho(K) < 1$, điều kiện này rất ngặt và dễ bị vi phạm ngay cả trên các hệ ma trận cấp 2 đơn giản. Các điều kiện trong [44, 45] đòi hỏi tham số trễ phải thỏa mãn $q > \frac{1}{(\gamma - 1)^2}$. Trong khi đó, cách tiếp cận bằng nguyên lý so sánh $M$-ma trận của luận án chỉ yêu cầu $-M$ là $M$-ma trận không suy biến (hoặc $\gamma I - M$ không suy biến), có hiệu lực với mọi $q \in (0, 1)$ mà không cần ràng buộc ngặt về giá trị của trễ.

4. Tại sao luận án lại khảo sát mô hình Nicholson với tốc độ suy giảm phi tuyến?

Mô hình Nicholson cổ điển giả định tỉ lệ tử vong của quần thể phụ thuộc tuyến tính vào số lượng cá thể. Giả thiết này chỉ phù hợp với các quần thể có mật độ cá thể thấp. Trong thực tế (như các khu bảo tồn biển hoặc quần thể mật độ cao), tỉ lệ tử vong thường chịu áp lực cạnh tranh gay gắt và phụ thuộc phi tuyến vào mật độ cá thể, được mô tả bởi các hàm suy giảm phi tuyến dạng $D(t, N) = a(t) - b(t)e^{-N}$ (Type-I).

5. Khái niệm tập hút mũ suy rộng trong Chương 3 được hiểu như thế nào?

Một tập compact $\Omega$ được gọi là tập hút mũ suy rộng nếu khoảng cách từ quỹ đạo nghiệm đến tập $\Omega$ thỏa mãn đánh giá suy giảm $\rho(x(t), \Omega) \le \kappa(|x_0|\infty) e^{-\sigma(t)}$, trong đó $\sigma(t)$ là một hàm không giảm thỏa mãn $\lim{t \to \infty} \sigma(t) = \infty$ (ví dụ $\sigma(t) = \hat{\alpha} \ln(1+t)$). Khái niệm này tổng quát hơn tập hút mũ thông thường (khi $\sigma(t) = \alpha t$) và cho phép thiết lập tính tiêu hao toàn cục cho các hệ động lực có tốc độ suy giảm nghiệm dạng đa thức hoặc hàm logarit.

Kết luận

Luận án tiến sĩ đã giải quyết trọn vẹn các bài toán định tính về tính ổn định trong thời gian hữu hạn, tính tiêu hao toàn cục và dáng điệu tiệm cận của nghiệm cho hai lớp hệ phương trình vi phân phi tuyến có trễ quan trọng: mạng nơron Hopfield với hệ số biến thiên chứa trễ tỉ lệ và mô hình sinh thái quần thể Nicholson với tốc độ suy giảm phi tuyến. Bằng việc kết hợp nguyên lý so sánh vi-tích phân với lý thuyết phổ của $M$-ma trận, tác giả đã khắc phục được các giới hạn của phương pháp đổi biến và bất đẳng thức Halanay kinh điển, đưa ra các tiêu chuẩn tường minh dễ kiểm tra cho các hệ động lực không dừng. Các kết quả lý thuyết được minh chứng bằng các ví dụ số học và 03 công trình công bố trên các tạp chí quốc tế uy tín trong danh mục ISI.