Tổng quan luận án

Luận án tiến sĩ toán học của nghiên cứu sinh Mai Nam Phong với đề tài "Nghiên cứu tính chất nghiệm của một số dạng phương trình và hệ phương trình sai phân phi tuyến" thuộc chuyên ngành Toán giải tích (mã số 62.01), được thực hiện tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội dưới sự hướng dẫn khoa học của PGS. Vũ Văn Khương và PGS. Đặng Đình Châu.

Phương trình sai phân giữ vai trò then chốt trong lý thuyết hệ động lực rời rạc. Các mô hình sai phân xuất hiện tự nhiên khi rời rạc hóa hoặc lấy nghiệm bằng số của các phương trình vi phân mô tả các quy luật sinh học, vật lý, kỹ thuật và kinh tế. Một số mô hình thực tế tiêu biểu được văn bản viện dẫn bao gồm:

  • Mô hình sinh trưởng của một loại cây hàng năm [33]: $$x_{n+1} = \frac{\lambda x_n}{(1 + a x_n)^p + b \lambda x_n}, \quad n = 0, 1, 2, \dots \quad (1)$$ với $a, b, p \in (0, \infty)$, $\lambda \in [1, +\infty)$ và $x_0 > 0$.
  • Mô hình sản xuất tế bào máu Mackey-Glass [33]: $$x_{n+1} = \alpha x_n + \frac{\beta}{x_{n-k}^p}, \quad n = 0, 1, 2, \dots \quad (2)$$ với $\alpha \in [0, 1)$, $p, \beta \in (0, \infty)$, $k \in \mathbb{Z}^+$ và các giá trị ban đầu $x_{-k}, \dots, x_0 > 0$.
  • Mô hình tương tác vật chủ – ký sinh của R. May [36]: $$x_{n+1} = \frac{\alpha x_n}{1 + \beta y_n}, \quad y_{n+1} = \frac{\beta x_n y_n}{1 + \beta y_n}, \quad n = 0, 1, 2, \dots \quad (3)$$ với $\alpha, \beta \in (0, \infty)$ và $x_0, y_0 > 0$.
  • Mô hình dân số loài bọ cánh cứng [34]: $$x_{n+1} = a x_n + b x_{n-2} e^{-c_1 x_n - c_2 x_{n-2}}, \quad n = 0, 1, 2, \dots \quad (9)$$ với $a \in (0, 1)$, $b \in (0, \infty)$, $c_1, c_2 \in [0, \infty)$, $c_1 + c_2 > 0$ và $x_{-2}, x_{-1}, x_0 > 0$.
  • Mô hình dân số loài muỗi [26]: $$x_{n+1} = (a x_n + b x_{n-1} e^{-x_{n-1}})e^{-x_n}, \quad n = 0, 1, 2, \dots \quad (10)$$ với $a \in (0, 1)$, $b \in [0, \infty)$ và $x_{-1}, x_0 > 0$.
  • Mô hình loài ruồi xanh do Nicholson đề xuất [33]: $$x_{n+1} = (1 - \alpha) x_n + \beta x_{n-k} e^{-\gamma x_{n-k}}, \quad n = 0, 1, 2, \dots \quad (11)$$ với $\alpha \in (0, 1)$, $\beta \in (\alpha, \infty)$, $\gamma \in (0, \infty)$, $k \in \mathbb{Z}^+$ và $x_{-k}, \dots, x_0 > 0$.

Trong nghiên cứu định tính phương trình sai phân, việc tìm kiếm công thức nghiệm tường minh thường không thực hiện được hoặc biểu thức nghiệm quá phức tạp. Do đó, yêu cầu nghiên cứu trực tiếp tính chất định tính của nghiệm (như tính ổn định, tính bị chặn, tính dao động, khoảng bất biến và dáng điệu tiệm cận) đóng vai trò căn bản.

Mục tiêu cụ thể của luận án được đánh số và xác định gồm bốn nội dung:

  1. Xây dựng dạng tiệm cận của nghiệm dương không dao động của hai dạng phương trình sai phân hữu tỷ.
  2. Tính ổn định tiệm cận toàn cục của điểm cân bằng một dạng phương trình sai phân hữu tỷ.
  3. Tính bị chặn, khoảng bất biến của nghiệm và tính ổn định của nghiệm cân bằng dương của hai dạng hệ hai phương trình sai phân phi tuyến.
  4. Tốc độ hội tụ nghiệm của một dạng hệ hai phương trình sai phân phi tuyến.

Đối tượng nghiên cứu của luận án là các lớp phương trình sai phân hữu tỷ cấp cao và hệ hai phương trình sai phân phi tuyến ôtônôm. Phạm vi nghiên cứu tập trung vào không gian giải tích thực $\mathbb{R}$ và $\mathbb{R}^k$, khảo sát nghiệm thực với các điều kiện ban đầu dương, phân tích điểm cân bằng dương, tính ổn định cục bộ, tính ổn định tiệm cận toàn cục, cấu trúc nửa chu kỳ dao động và tốc độ hội tụ của dãy nghiệm khi $n \to \infty$.

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

Nghiên cứu định tính phương trình sai phân phi tuyến đã phát triển mạnh từ thập niên 90 của thế kỷ XX và những năm đầu thế kỷ XXI. Luận án điểm lại các công trình tiêu biểu trong và ngoài nước:

Trong nước, các hướng nghiên cứu về phương trình sai phân phi tuyến được ghi nhận qua các công trình của Vũ Ngọc Phát, Nguyễn Sinh Bảy [4, 5, 51], Đặng Vũ Giang, Đinh Công Hướng [24, 25, 28], Vũ Văn Khương [30, 31, 32].

Ở ngoài nước, các hệ thống kết quả nền tảng bao gồm:

  • M. R. S. Kulenović và G. Ladas (2001) trong chuyên khảo [35] đã tổng hợp các kết quả về tính bị chặn, tính ổn định toàn cục và tính tuần hoàn của phương trình sai phân hữu tỷ cấp hai: $$x_{n+1} = \frac{\alpha + \beta x_n + \gamma x_{n-1}}{A + B x_n + C x_{n-1}}, \quad n = 0, 1, 2, \dots \quad (4)$$ cùng các trường hợp riêng khi triệt tiêu các hệ số: phương trình $x_{n+1} = \frac{\alpha + \beta x_n}{A + B x_n}$ (khi $\gamma = C = 0$), $x_{n+1} = \frac{\beta x_n}{A + C x_{n-1}}$ (khi $\alpha = \gamma = B = 0$), và $x_{n+1} = \frac{\alpha + \beta x_n}{C x_{n-1}}$ (khi $\gamma = A = B = 0$).
  • Ladas (2001) [35] cũng trình bày các kết quả về tính bị chặn, tính ổn định và tính tuần hoàn của lớp phương trình sai phân hữu tỷ cấp ba: $$x_{n+1} = \frac{\alpha + \beta x_n + \gamma x_{n-1} + \delta x_{n-2}}{A + B x_n + C x_{n-1} + D x_{n-2}}, \quad n = 0, 1, 2, \dots \quad (8)$$
  • R. DeVault, G. Ladas, S. W. Schultz [17] nghiên cứu tính giới nội của phương trình: $$x_{n+1} = \frac{A}{x_n^p} + \frac{B}{x_{n-1}^q}, \quad n = 0, 1, \dots \quad (2.1)$$
  • E. Camouzis, G. Ladas đề xuất nghiên cứu tính ổn định và giới nội của phương trình: $$x_{n+1} = \frac{A}{x_n x_{n-1}} + \frac{1}{x_{n-3} x_{n-4}}, \quad n = 0, 1, 2, \dots \quad (2.2)$$
  • S. Stević [59] nghiên cứu tính giới nội và bền vững của phương trình: $$x_{n+1} = \sum_{i=0}^k \frac{\alpha_i}{x_{n-i}^{p_i}}, \quad n = 0, 1, 2, \dots \quad (2.3)$$
  • L. Berg [6], S. Stević [59], L. Berenhaut [10, 11] xây dựng phương pháp đánh giá tiệm cận nghiệm của phương trình sai phân.
  • X. Zhu [43, 44, 45, 46, 47] và Nesemann [48] khảo sát tính ổn định tiệm cận toàn cục của các phương trình sai phân hữu tỷ dạng đối xứng và nghịch đảo: $$x_{n+1} = \frac{x_n + x_{n-1} x_{n-2}}{x_n x_{n-1} + x_{n-2}} \quad (2.24)$$ $$x_{n+1} = \frac{x_{n-1} + x_n x_{n-2}}{x_n x_{n-1} + x_{n-2}} \quad (2.25)$$ $$x_{n+1} = \frac{x_n x_{n-1} + x_{n-2} + a}{x_n + x_{n-1} x_{n-2} + a} \quad (2.26)$$ $$x_{n+1} = \frac{x_{n-1} + x_n x_{n-2} + a}{x_n x_{n-1} + x_{n-2} + a} \quad (2.27)$$ $$x_{n+1} = \frac{x_n x_{n-2} + 1}{x_n + x_{n-2}} \quad (2.28)$$ $$x_{n+1} = \frac{x_{n-1} x_{n-2} + 1}{x_{n-1} + x_{n-2}} \quad (2.29)$$

Khoảng trống nghiên cứu luận án lựa chọn giải quyết: Phần lớn các công trình trước đây chỉ dừng lại ở việc chứng minh tính giới nội, tính bền vững, tính tuần hoàn hoặc tính ổn định đơn thuần, chưa mô tả cấu trúc giải tích tường minh của nghiệm khi tiến dần về điểm cân bằng (dạng tiệm cận nghiệm). Đồng thời, các kết quả ổn định toàn cục của X. Zhu và Nesemann chủ yếu xét ở phương trình cấp ba; luận án đã mở rộng bài toán lên phương trình sai phân hữu tỷ cấp bốn chứa tham số $a \ge 0$ và khảo sát các hệ hai phương trình sai phân phi tuyến cấp cao.

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

Khung lý thuyết và khái niệm sử dụng

Luận án sử dụng hệ thống khái niệm chuẩn mực trong lý thuyết phương trình sai phân cấp $(k+1)$ dạng $x_{n+1} = F(x_n, x_{n-1}, \dots, x_{n-k})$ (1.1) và hệ phương trình sai phân $Y_{n+1} = G(n, Y_n)$ (1.4):

  • Điểm cân bằng: Giá trị $\bar{x}$ thỏa mãn $\bar{x} = F(\bar{x}, \dots, \bar{x})$ hoặc $\bar{Y} = G(n, \bar{Y})$.
  • Các tiêu chuẩn ổn định: Ổn định địa phương, hút địa phương, ổn định tiệm cận địa phương, hút toàn cục, và ổn định tiệm cận toàn cục (theo Elaydi [14], Kocic & Ladas [33], Kulenović & Ladas [35, 36]).
  • Phương pháp tuyến tính hóa: Hàm $F$ khả vi liên tục tại lân cận $\bar{x}$, phương trình tuyến tính hóa có dạng $y_{n+1} = \sum_{i=0}^k q_i y_{n-i}$ với $q_i = \frac{\partial F}{\partial u_i}(\bar{x}, \dots, \bar{x})$. Phương trình đặc trưng tương ứng: $$\lambda^{k+1} - q_0 \lambda^k - q_1 \lambda^{k-1} - \dots - q_k = 0 \quad (1.3)$$ Điều kiện đủ để mọi nghiệm nằm trong đĩa đơn vị mở là $\sum_{i=0}^k |q_i| < 1$.
  • Định lý Schur-Cohn cho hệ hai phương trình $x_{n+1} = f(x_n, y_n), y_{n+1} = g(x_n, y_n)$ (1.5): Điểm cân bằng $(\bar{x}, \bar{y})$ ổn định tiệm cận địa phương khi và chỉ khi mọi nghiệm của phương trình đặc trưng $\lambda^2 - \operatorname{tr}(J_F)\lambda + \det(J_F) = 0$ nằm trong đĩa đơn vị, tương đương bất đẳng thức: $$|\operatorname{tr} J_F(\bar{x}, \bar{y})| < 1 + \det J_F(\bar{x}, \bar{y}) < 2$$ trong đó $J_F$ là ma trận Jacobian tại điểm cân bằng.
  • Khái niệm về dao động và nửa chu kỳ:
    • Nửa chu kỳ dương (nửa vòng dương): Xâu số hạng liên tiếp ${x_l, \dots, x_m}$ thỏa mãn $x_i \ge \bar{x}$.
    • Nửa chu kỳ âm (nửa vòng âm): Xâu số hạng liên tiếp ${x_l, \dots, x_m}$ thỏa mãn $x_i < \bar{x}$.
    • Nghiệm không dao động: Tồn tại $N$ để $x_n \ge \bar{x}, \forall n \ge N$ hoặc $x_n < \bar{x}, \forall n \ge N$.
    • Nghiệm dao động ngặt: Với mỗi $n_0 \ge 0$, tồn tại $n_1, n_2 \ge n_0$ sao cho $(x_{n_1} - \bar{x})(x_{n_2} - \bar{x}) < 0$.
    • Nghiệm bị chặn và bền vững: Tồn tại $P, Q > 0$ sao cho $P \le x_n \le Q, \forall n \ge -k$.

Phương pháp nghiên cứu thực tế

Luận án sử dụng phương pháp nghiên cứu định tính giải tích thuần túy, bao gồm:

  1. Phương pháp xấp xỉ tiệm cận của Berg – Stević: Dự đoán dạng giải tích tiệm cận $\phi_n = \bar{x} + a t_0^n + b t_0^{2n}$, thiết lập đa thức đặc trưng để tìm nghiệm riêng $t_0 \in (0, 1)$, sau đó xây dựng hai dãy kẹp ${y_n}$ và ${z_n}$ thỏa mãn $y_n \le x_n \le z_n$. Sử dụng Định lý 2.1 về ánh xạ đơn điệu nhiều biến trên nón dương để khẳng định sự tồn tại nghiệm tiệm cận.
  2. Phương pháp phân tích cấu trúc nửa chu kỳ và dãy con: Thiết lập các đồng nhất thức sai phân cho hiệu $(x_{n+1}-1)$, $(x_{n+1}-x_{n-1})$, $(x_{n+1}-x_{n-3})$, phân loại các điều kiện ban đầu thành 7 trường hợp độc lập để xác định quy luật xuất hiện tuần hoàn của các nửa chu kỳ. Kết hợp tính đơn điệu và bị chặn của các dãy con theo bước nhảy chu kỳ để lấy giới hạn, chứng minh điểm cân bằng là điểm hút toàn cục.
  3. Phương pháp Jacobian và không gian pha rời rạc: Xác định tập bất biến $I \times J$, tính toán các đạo hàm riêng của toán tử chuyển tiếp, khảo sát tính chất phổ của ma trận Jacobian để suy ra tính ổn định và tốc độ hội tụ.

Nguồn tư liệu là các công trình khoa học, chuyên khảo và bài báo chuyên ngành Toán giải tích trong danh mục tài liệu tham khảo của luận án.

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

Chương 1. Kiến thức chuẩn bị

Chương 1 hệ thống hóa các định nghĩa, bổ đề và định lý nền tảng làm công cụ cho các chương sau:

  • Trình bày định nghĩa phương trình sai phân cấp $(k+1)$, nghiệm cân bằng và điểm cân bằng.
  • Nêu chi tiết 6 trạng thái ổn định của điểm cân bằng: ổn định địa phương, hút địa phương, ổn định tiệm cận địa phương, hút toàn cục, ổn định tiệm cận toàn cục, không ổn định.
  • Trình bày Định lý ổn định tuyến tính hóa và điều kiện tổng hệ số đạo hàm riêng bé hơn 1.
  • Định nghĩa khái niệm nửa chu kỳ dương, nửa chu kỳ âm, độ dài nửa chu kỳ, nghiệm dao động, nghiệm dao động ngặt, nghiệm không dao động, nghiệm bị chặn và bền vững.
  • Thiết lập cơ sở lý thuyết cho hệ phương trình sai phân tổng quát $Y_{n+1} = G(n, Y_n)$ và hệ hai phương trình sai phân ôtônôm (1.5), định nghĩa tập bất biến $I \times J$, ma trận Jacobian $J_F$ và định lý Schur-Cohn.

Chương 2. Ba dạng phương trình sai phân hữu tỷ

Chương 2 tập trung khảo sát ba dạng phương trình sai phân hữu tỷ cấp cao:

1. Dáng điệu tiệm cận nghiệm của hai dạng phương trình sai phân hữu tỷ (Mục 2.1)

  • Phương trình dạng thứ nhất (2.4): $$x_{n-k} = \sum_{i=1}^{k-1} \frac{A_i}{x_{n-i}} - x_n - 1$$ với $A_i \in [0, \infty)$, $\sum_{i=1}^{k-1} A_i > 1$, giá trị ban đầu $x_{-k}, \dots, x_{-1} > 0$.

    • Điểm cân bằng dương: $\bar{x} = \sqrt{\sum_{i=1}^{k-1} A_i - 1}$.
    • Phương trình đặc trưng: $p(t) = \bar{x}^2 t^k + A_1 t^{k-1} + \dots + A_{k-1} - 1 = 0$ (2.9). Vì $p(0) = A_{k-1} - 1 < 0$, $p(1) = 2\bar{x}^2 > 0$ và $p'(t) > 0, \forall t \in (0, 1]$, tồn tại duy nhất nghiệm $t_0 \in (0, 1)$.
    • Định lý 2.2: Phương trình (2.4) có nghiệm dương không dao động hội tụ đơn điệu về $\bar{x}$ với dạng biểu diễn tiệm cận: $$x_n = \bar{x} + a t_0^n + b t_0^{2n} + o(t_0^{2n})$$ trong đó hệ số $b$ được xác định tường minh: $$b = \frac{a^2 t_0^{2k}}{\bar{x} p(t_0^2)} \left[ \sum_{i=1}^{k-1} A_i t_0^{-2i} - \left(\sum_{i=1}^{k-1} A_i t_0^{-i} + \bar{x}^2\right)^2 \right]$$

  • Phương trình dạng thứ hai (2.5): $$x_{n+1} = \frac{A_1}{x_n x_{n-1} x_{n-2}} + \frac{A_2}{x_{n-1} x_{n-2} x_{n-3}} + \frac{A_3}{x_{n-2} x_{n-3} x_{n-4}} - \frac{1}{x_{n-5}^3}$$ với $A_1, A_2, A_3 \in [0, \infty)$ và $A = A_1 + A_2 + A_3 - 1 > 0$.

    • Điểm cân bằng dương: $\bar{x} = \sqrt[4]{A}$.
    • Phương trình đặc trưng tương ứng với phương trình tuyến tính hóa (2.17): $$p(t) = \bar{x}^4 t^6 + A_1 t^5 + (A_1 + A_2)t^4 + (A_1 + A_2 + A_3)t^3 + (A_2 + A_3)t^2 + A_3 t - 3 = 0 \quad (2.18)$$ Do $p(0) = -3 < 0$, $p(1) = 4A > 0$ và $p'(t) > 0, \forall t \in (0, 1)$, tồn tại nghiệm duy nhất $t_0 \in (0, 1)$.
    • Định lý 2.3: Tồn tại nghiệm dương không dao động hội tụ đơn điệu về điểm cân bằng dương $\bar{x} = \sqrt[4]{A}$ với dạng tiệm cận $x_n = \bar{x} + a t_0^n + b t_0^{2n} + o(t_0^{2n})$.
    • Định lý 2.4: Mở rộng kết quả cho phương trình hữu tỷ tổng quát cấp $(k+m-2)$: $$x_{n+1} = \sum_{i=1}^{k-1} \frac{A_i}{\prod_{j=0}^{m-1} x_{n-i-j+1}} - \frac{1}{x_{n-k-m+2}^m}$$ hội tụ về điểm cân bằng dương $\bar{x} = \sqrt[m+1]{\sum_{i=1}^{k-1} A_i - 1}$.

2. Tính ổn định của một dạng phương trình sai phân hữu tỷ (Mục 2.2)

  • Khảo sát phương trình sai phân hữu tỷ cấp bốn (2.30): $$x_{n+1} = \frac{x_{n-1} x_{n-3} + x_{n-1}^2 + a}{x_{n-1}^2 x_{n-3} + x_{n-1} + a}, \quad n = 0, 1, 2, \dots$$ với $a \in [0, \infty)$ và các giá trị ban đầu $x_{-3}, x_{-2}, x_{-1}, x_0 \in (0, \infty)$.
    • Điểm cân bằng dương duy nhất: $\bar{x} = 1$.
    • Bổ đề 2.1: Nghiệm dương ${x_n}$ là tầm thường về cuối (nghĩa là $x_n = 1, \forall n \ge n_0$) khi và chỉ khi đồng thời: $$(x_{-3} - 1)(x_{-1} - 1) = 0 \quad (2.33) \quad \text{và} \quad (x_{-2} - 1)(x_0 - 1) = 0 \quad (2.34)$$
    • Bổ đề 2.2: Khi các giá trị ban đầu $x_i \ne 1$ ($i = -3, -2, -1, 0$), các hệ thức sau nghiệm đúng với mọi $n \ge 0$:
      • (a) $(x_{n+1} - 1)(x_{n-1} - 1)(x_{n-3} - 1) < 0$
      • (b) $(x_{n+1} - x_{n-1})(x_{n-1} - 1) < 0$
      • (c) $(x_{n+1} - x_{n-3})(x_{n-3} - 1) < 0$
      • (d) $(x_{n+1} - \frac{1}{x_{n-1}})(x_{n-1} - 1) > 0$
      • (e) $(x_{n+1} - \frac{1}{x_{n-3}})(x_{n-3} - 1) > 0$
    • Định lý 2.5: Phân tích 7 trường hợp giá trị ban đầu của nghiệm dao động ngặt, chứng minh quy luật xuất hiện của độ dài các nửa chu kỳ dương ($+$) và âm ($-$) có tính tuần hoàn thuộc một trong các dạng:
      • Dạng 1: $\dots, 1^-, 2^+, 1^-, 2^+, \dots$ (từ Trường hợp 1: $x_{p-3}>1, x_{p-2}<1, x_{p-1}>1, x_p>1$)
      • Dạng 2: $\dots, 1^+, 1^-, 1^+, 3^-, 1^+, 1^-, 1^+, 3^-, \dots$ (từ Trường hợp 2 & 3)
      • Dạng 3: $\dots, 4^+, 2^-, 4^+, 2^-, \dots$ (từ Trường hợp 4: $x_{p-3}>1, x_{p-2}<1, x_{p-1}<1, x_p<1$)
      • Dạng 4: $\dots, 5^+, 1^-, 5^+, 1^-, \dots$ (từ các Trường hợp 5, 6, 7 có chứa số hạng bằng 1)
    • Định lý 2.6: Điểm cân bằng $\bar{x} = 1$ của phương trình (2.30) là ổn định tiệm cận toàn cục trên $(0, \infty)$.

Chương 3. Hai dạng hệ phương trình sai phân phi tuyến

Chương 3 trình bày các kết quả nghiên cứu định tính đối với hai dạng hệ hai phương trình sai phân phi tuyến:

  • Mục 3.1: Nghiên cứu tính bị chặn, khoảng bất biến và tính ổn định tiệm cận của nghiệm đối với dạng hệ phương trình sai phân phi tuyến thứ nhất; chứng minh dáng điệu tiệm cận của nghiệm thông qua ma trận Jacobian và đưa ra các ví dụ minh họa bằng số.
  • Mục 3.2: Nghiên cứu dạng hệ hai phương trình sai phân phi tuyến thứ hai; chứng minh tính bị chặn, tính ổn định của điểm cân bằng dương, xác định dáng điệu toàn cục của nghiệm, đánh giá tốc độ hội tụ của nghiệm về điểm cân bằng dương và trình bày các ví dụ số kiểm chứng.

Bảng tổng hợp các dạng phương trình và hệ phương trình trong luận án

Dạng phương trình / Hệ phương trình Cấp / Bậc Điểm cân bằng $\bar{x}$ hoặc $(\bar{x}, \bar{y})$ Tính chất định tính chính được chứng minh
$x_{n-k} = \sum_{i=1}^{k-1} \frac{A_i}{x_{n-i}} - x_n - 1$ Cấp $k$ $\bar{x} = \sqrt{\sum_{i=1}^{k-1} A_i - 1}$ Tồn tại nghiệm dương không dao động hội tụ về $\bar{x}$ với tiệm cận $x_n = \bar{x} + a t_0^n + b t_0^{2n} + o(t_0^{2n})$
$x_{n+1} = \sum_{i=1}^3 \frac{A_i}{\prod x} - \frac{1}{x_{n-5}^3}$ Cấp 6 $\bar{x} = \sqrt[4]{A_1 + A_2 + A_3 - 1}$ Nghiệm không dao động hội tụ đơn điệu về $\bar{x}$ với tiệm cận bậc hai theo nghiệm đặc trưng $t_0 \in (0, 1)$
$x_{n+1} = \frac{x_{n-1} x_{n-3} + x_{n-1}^2 + a}{x_{n-1}^2 x_{n-3} + x_{n-1} + a}$ Cấp 4 $\bar{x} = 1$ Ổn định tiệm cận toàn cục; độ dài các nửa chu kỳ xuất hiện tuần hoàn theo 4 dạng quy luật
Hệ hai phương trình sai phân phi tuyến (Mục 3.1) Hệ cấp cao $(\bar{x}, \bar{y})$ dương Xác định tập bất biến, tính bị chặn, ổn định tiệm cận địa phương
Hệ hai phương trình sai phân phi tuyến (Mục 3.2) Hệ cấp cao $(\bar{x}, \bar{y})$ dương Tính bị chặn, ổn định toàn cục và ước lượng tốc độ hội tụ nghiệm

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

Luận án đã đạt được các đóng góp mới về mặt lý thuyết giải tích và phương pháp nghiên cứu hệ động lực rời rạc:

  1. Xây dựng dạng tiệm cận giải tích tường minh: Thiết lập thành công công thức xấp xỉ tiệm cận bậc hai cho nghiệm dương không dao động của hai lớp phương trình sai phân hữu tỷ cấp cao (2.4) và (2.5), xác định chính xác số mũ hội tụ thông qua nghiệm dương duy nhất $t_0 \in (0, 1)$ của đa thức đặc trưng liên kết.
  2. Chứng minh tính ổn định tiệm cận toàn cục cho phương trình hữu tỷ cấp bốn chứa tham số: Giải quyết trọn vẹn bài toán ổn định toàn cục của điểm cân bằng $\bar{x} = 1$ cho phương trình (2.30) với tham số $a \ge 0$. Đóng góp này làm sáng tỏ cơ chế dao động thông qua việc phân loại 7 trường hợp giá trị ban đầu và chứng minh tính tuần hoàn của độ dài các nửa chu kỳ.
  3. Phát triển kết quả định tính cho hệ phương trình sai phân phi tuyến: Mở rộng các tiêu chuẩn về tính bị chặn, khoảng bất biến, tính ổn định và tốc độ hội tụ nghiệm cho hai dạng hệ hai phương trình sai phân phi tuyến cấp cao.

Các kết quả khoa học của luận án đã được công bố trong 05 bài báo khoa học (trong đó có 01 bài báo đứng tên độc lập của tác giả và 04 bài báo viết chung cùng cán bộ hướng dẫn).

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

Căn cứ theo nội dung và phạm vi nghiên cứu được trình bày trong văn bản:

  • Hạn chế:
    • Các kết quả tiệm cận và ổn định chủ yếu áp dụng cho phương trình và hệ phương trình ôtônôm (hệ số hằng), chưa xét đến trường hợp các hệ số biến thiên theo bước thời gian $n$ (non-autonomous).
    • Nghiên cứu dạng tiệm cận chỉ tập trung vào lớp nghiệm dương không dao động; đối với các nghiệm dao động ngặt, luận án chỉ phân tích cấu trúc nửa chu kỳ chứ chưa thiết lập biểu thức tiệm cận giải tích tương đương.
    • Hệ phương trình sai phân mới khảo sát ở không gian hai chiều ($k=2$).
  • Hướng nghiên cứu tiếp theo:
    • Mở rộng nghiên cứu dạng tiệm cận và tính ổn định cho các phương trình sai phân hữu tỷ và hệ phương trình sai phân có hệ số thay đổi theo thời gian.
    • Nghiên cứu tính chất định tính của hệ phương trình sai phân phi tuyến với số chiều tổng quát ($k \ge 3$).
    • Khảo sát các hiện tượng phân nhánh (bifurcation) và tính hỗn loạn (chaos) trong các mô hình rời rạc tương ứng.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên khảo hữu ích cho:

  • Nghiên cứu sinh, học viên cao học và giảng viên chuyên ngành Toán giải tích, Phương trình vi phân và Hệ động lực: Cung cấp kỹ thuật chi tiết về phương pháp xấp xỉ tiệm cận Berg – Stević và phương pháp phân tích chu kỳ nghiệm.
  • Các nhà nghiên cứu mô hình hóa toán học trong Sinh học, Y học và Kinh tế: Tham khảo các kỹ thuật khảo sát tính ổn định và dáng điệu tiệm cận của các phương trình sai phân hữu tỷ xuất phát từ các mô hình thực tế như mô hình dân số muỗi, bọ cánh cứng, ruồi xanh Nicholson và mô hình tế bào máu Mackey-Glass.
  • Phần có giá trị tham khảo học thuật cao: Mục 2.1 (kỹ thuật xây dựng dãy kẹp để tìm tiệm cận nghiệm) và Mục 2.2 (kỹ thuật phân loại 7 trường hợp nửa chu kỳ để chứng minh ổn định tiệm cận toàn cục).

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

1. Điểm cân bằng của phương trình sai phân hữu tỷ dạng (2.4) được xác định như thế nào và dạng tiệm cận của nghiệm có biểu thức ra sao?
Điểm cân bằng dương của phương trình (2.4) được xác định bởi công thức $\bar{x} = \sqrt{\sum_{i=1}^{k-1} A_i - 1}$ (với điều kiện $\sum_{i=1}^{k-1} A_i > 1$). Dạng tiệm cận của nghiệm dương không dao động hội tụ về $\bar{x}$ có cấu trúc $x_n = \bar{x} + a t_0^n + b t_0^{2n} + o(t_0^{2n})$, trong đó $t_0$ là nghiệm dương duy nhất thuộc khoảng $(0, 1)$ của phương trình đặc trưng $p(t) = \bar{x}^2 t^k + A_1 t^{k-1} + \dots + A_{k-1} - 1 = 0$.

2. Điều kiện cần và đủ để nghiệm của phương trình sai phân hữu tỷ cấp bốn (2.30) trở thành nghiệm tầm thường về cuối là gì?
Theo Bổ đề 2.1 của luận án, nghiệm dương ${x_n}$ của phương trình (2.30) là tầm thường về cuối (nghĩa là tồn tại $n_0$ để $x_n = 1, \forall n \ge n_0$) khi và chỉ khi các giá trị ban đầu thỏa mãn đồng thời hai điều kiện: $(x_{-3} - 1)(x_{-1} - 1) = 0$ và $(x_{-2} - 1)(x_0 - 1) = 0$.

3. Quy luật xuất hiện của các nửa chu kỳ dương và âm của nghiệm dao động ngặt trong phương trình (2.30) diễn ra như thế nào?
Dựa trên Bổ đề 2.2 và Định lý 2.5, khi phân tích 7 trường hợp quan hệ giá trị ban đầu, độ dài của các nửa chu kỳ dương ($+$) và âm ($-$) xuất hiện tuần hoàn theo một trong bốn quy luật: $\dots, 1^-, 2^+, 1^-, 2^+, \dots$; hoặc $\dots, 1^+, 1^-, 1^+, 3^-, 1^+, 1^-, 1^+, 3^-, \dots$; hoặc $\dots, 4^+, 2^-, 4^+, 2^-, \dots$; hoặc $\dots, 5^+, 1^-, 5^+, 1^-, \dots$.

4. Phương pháp nào được sử dụng để chứng minh sự tồn tại nghiệm tiệm cận trong Chương 2?
Tác giả sử dụng phương pháp xấp xỉ tiệm cận của Berg – Stević kết hợp Định lý 2.1 về hàm liên tục không giảm theo từng biến trên khoảng đóng $I = [\bar{x}, \infty)$. Kỹ thuật thực hiện bằng cách xây dựng hai dãy $y_n = \bar{x} + a t_0^n + q_1 t_0^{2n}$ và $z_n = \bar{x} + a t_0^n + q_2 t_0^{2n}$ để thiết lập bất đẳng thức kẹp $y_{n-k} \le f(n, y_{n-k+1}, \dots, z_{n+1}) \le z_{n-k}$ với $n$ đủ lớn.

5. Điều kiện Schur-Cohn để điểm cân bằng của hệ hai phương trình sai phân đạt ổn định tiệm cận địa phương là gì?
Theo Định lý ổn định tuyến tính hóa cho hệ hai phương trình trong Chương 1, điểm cân bằng $(\bar{x}, \bar{y})$ đạt ổn định tiệm cận địa phương khi và chỉ khi mọi nghiệm của phương trình đặc trưng $\lambda^2 - \operatorname{tr}(J_F)\lambda + \det(J_F) = 0$ nằm trong đĩa đơn vị, tương đương với điều kiện bất đẳng thức: $|\operatorname{tr} J_F(\bar{x}, \bar{y})| < 1 + \det J_F(\bar{x}, \bar{y}) < 2$.

Kết luận

Luận án tiến sĩ của tác giả Mai Nam Phong đã giải quyết trọn vẹn các bài toán định tính trọng tâm đối với ba dạng phương trình sai phân hữu tỷ cấp cao và hai dạng hệ hai phương trình sai phân phi tuyến. Tác giả đã thiết lập biểu thức tiệm cận bậc hai tường minh cho nghiệm dương không dao động và chứng minh tính ổn định tiệm cận toàn cục của điểm cân bằng dương thông qua việc phân tích quy luật nửa chu kỳ tuần hoàn. Công trình bổ sung các phương pháp giải tích chặt chẽ vào lý thuyết định tính hệ động lực rời rạc và cung cấp cơ sở toán học cho việc nghiên cứu các mô hình sai phân phi tuyến trong thực tiễn.