Tổng quan luận án

Luận án tiến sĩ toán học với đề tài "Sự ổn định của một vài lớp phương trình sai phân dạng phân thức" do nghiên cứu sinh Trần Hồng Thái 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 (năm 2019), thuộc chuyên ngành Toán giải tích (mã số: 9460101.02), dưới sự hướng dẫn khoa học của PGS. Vũ Văn Khương và TS. Lê Đình Định.

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

Phương trình sai phân là công cụ toán học mô tả các quá trình biến đổi trạng thái trong các hệ động lực rời rạc thuộc nhiều lĩnh vực khoa học tự nhiên và kinh tế - xã hội như dự báo dân số, mô hình sản sinh tế bào máu, mô hình phát triển thực vật, quan hệ ký sinh - vật chủ hay tương tác sinh thái côn trùng. Đối với các phương trình sai phân tuyến tính, lý thuyết giải và tìm nghiệm tường minh đã được hoàn thiện tương đối đầy đủ. Tuy nhiên, đối với phương trình sai phân phi tuyến và hệ phương trình sai phân phi tuyến, do cấu trúc biểu thức rất đa dạng nên chưa có một phương pháp tổng quát thống nhất để giải quyết tất cả các dạng.

Đặc biệt, lớp phương trình và hệ phương trình sai phân phi tuyến cấp lớn hơn một dạng phân thức thu hút sự quan tâm lớn vì thế hệ thứ $(n+1)$ phụ thuộc phi tuyến vào nhiều thế hệ trước đó. Nhiều mô hình cụ thể đã được các nhà toán học quốc tế khảo sát ở các cấp đơn giản, nhưng vẫn còn nhiều lớp phương trình sai phân dạng phân thức cấp cao với các số hạng bậc hai, bậc ba hoặc chứa các hàm tổng quát chưa được phân tích đầy đủ về tính ổn định của điểm cân bằng, tính bị chặn và dạng tiệm cận của nghiệm. Luận án được thực hiện nhằm giải quyết khoảng trống lý thuyết đối với các lớp phương trình và hệ phương trình sai phân dạng phân thức cụ thể này.

Mục tiêu và nhiệm vụ nghiên cứu

Mục đích chính của luận án là góp phần làm phong phú thêm lý thuyết định tính của phương trình sai phân thông qua việc nghiên cứu tính chất định tính của nghiệm và điểm cân bằng cho một số lớp phương trình và hệ phương trình sai phân dạng phân thức. Luận án tập trung giải quyết ba nhiệm vụ cụ thể:

  1. Nghiên cứu sự ổn định của nghiệm của hai lớp phương trình sai phân cấp hai dạng phân thức và một lớp phương trình sai phân cấp bốn dạng phân thức.
  2. Xây dựng dạng tiệm cận của nghiệm chứa hàm logarit tự nhiên của một dạng phương trình sai phân cấp ba dạng phân thức.
  3. Phân tích sự ổn định của điểm cân bằng của hai lớp hệ phương trình sai phân dạng phân thức cấp một và cấp hai.

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

  • Đối tượng nghiên cứu: Điểm cân bằng, tính ổn định địa phương, tính ổn định tiệm cận toàn cục, tính chất yên ngựa, tính bị chặn và biểu diễn tiệm cận nghiệm của các lớp phương trình sai phân phân thức cấp cao (cấp 2, cấp 3, cấp 4) và các hệ phương trình sai phân phân thức cấp 1, cấp 2.
  • Phạm vi nghiên cứu: Luận án giới hạn trong phạm vi các phương trình và hệ phương trình sai phân xác định trên không gian thực với các tham số và điều kiện ban đầu dương hoặc không âm.

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

Luận án điểm lại hệ thống các công trình nghiên cứu trong nước và quốc tế về phương trình sai phân phi tuyến và phương trình sai phân dạng phân thức:

Các công trình và tác giả tiền đề

  • Mô hình ứng dụng thực tế: Luận án viện dẫn mô hình sản sinh tế bào máu do Mackey và Glass đề xuất có dạng $x_{n+1} = \frac{\beta x_n}{1 + x_{n-k}^p}$; mô hình cây hàng năm $x_{n+1} = \frac{\lambda x_n}{(1 + ax_n)^p + b\lambda x_n}$; mô hình vật chủ - ký sinh của R. May với hệ $x_{n+1} = \frac{\alpha x_n}{1 + \beta y_n}, y_{n+1} = \frac{\beta x_n y_n}{1 + \beta y_n}$; mô hình tương tác giữa sâu đục cành và cây nho ở đồng bằng Texas $x_{n+1} = \frac{\alpha x_n}{e^{y_n} + \beta y_n}, y_{n+1} = \gamma(x_n + 1)y_n$.
  • Chuyên khảo quốc tế nền tảng: Chuyên khảo năm 2002 của M.R.S. Kulenović và G. Ladas tổng hợp các kết quả về tính giới nội, tính ổn định toàn cục và tính tuần hoàn của lớp phương trình sai phân phân thức cấp hai có dạng: $$x_{n+1} = \frac{\alpha + \beta x_n + \gamma x_{n-1}}{A + Bx_n + C x_{n-1}}$$ Chuyên khảo năm 2008 của E. Camouzis và G. Ladas nghiên cứu phương trình sai phân phân thức cấp ba: $$x_{n+1} = \frac{\alpha + \beta x_n + \gamma x_{n-1} + \delta x_{n-2}}{A + Bx_n + C x_{n-1} + D x_{n-2}}$$
  • Các nghiên cứu mở rộng về phương trình phân thức:
    • Amleh và cộng sự (1999) nghiên cứu phương trình $x_{n+1} = \alpha + \frac{x_{n-1}}{x_n}$.
    • Kosmala và cộng sự nghiên cứu phương trình $y_{n+1} = \frac{p + y_{n-1}}{qy_n + y_{n-1}}$.
    • L. Berg (2002, 2003) nghiên cứu sự tồn tại và dạng tiệm cận của nghiệm hội tụ về 0 của phương trình $x_{n+1} = \frac{x_{n-1}}{1 + x_n}$ và $x_n = \frac{x_{n-3}}{1 + x_{n-1}x_{n-2}}$.
    • G. Ladas, Nesemann, X.H. Zhu, L.X. Li và R.P. Agarwal 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 cấp ba, cấp bốn có tham số $a, b$.
  • Các tác giả trong nước: Phạm Kỳ Anh, Đặng Vũ Giang, Đinh Công Hướng, Vũ Văn Khương với các công trình về phương trình sai phân phi tuyến và lý thuyết định tính hệ động lực rời rạc.

Khoảng trống luận án lựa chọn giải quyết

Mặc dù các lớp phương trình sai phân tuyến tính bậc một của tử số và mẫu số đã được khảo sát tương đối nhiều, các lớp phương trình chứa số hạng bậc hai ở tử số và mẫu số như $x_{n+1} = ax_{n-1} - \frac{bx_n^2 + cx_{n-1}^2}{1 - dx_n + ex_{n-1}}$, phương trình có hai thành phần phân thức $x_{n+1} = \frac{\alpha x_{n-1}}{\beta x_n + x_{n-1}} + \frac{x_{n-1}}{x_n}$, phương trình cấp ba phi tuyến bậc ba có dạng tiệm cận chứa số hạng logarit tự nhiên, và phương trình cấp bốn chứa các hàm khả vi tổng quát $f, g$ vẫn chưa có lời giải hoàn chỉnh. Luận án định vị nghiên cứu vào việc phân tích định tính các lớp phương trình cụ thể này cùng các hệ phương trình sai phân cấp một và cấp hai tương ứng.

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

Khung lý thuyết và khái niệm cơ bản

Luận án sử dụng khung lý thuyết định tính phương trình sai phân trong không gian véc tơ thực $\mathbb{R}^k$:

  • Phương trình sai phân cấp $(k+1)$: Dạng tổng quát $x_{n+1} = F(x_n, x_{n-1}, \dots, x_{n-k})$ với ánh xạ liên tục $F: J^{k+1} \to J$, $J \subset \mathbb{R}$. Nghiệm cân bằng $x_n = \bar{x}, \forall n \ge -k$ thỏa mãn $\bar{x} = F(\bar{x}, \bar{x}, \dots, \bar{x})$.
  • Các tiêu chuẩn ổn định của điểm cân bằng $\bar{x}$:
    • Ổn định địa phương: $\forall \varepsilon > 0, \exists \delta > 0$ sao cho $\sum_{i=0}^k |x_{-i} - \bar{x}| < \delta \Rightarrow |x_n - \bar{x}| < \varepsilon, \forall n \ge 1$.
    • Hút địa phương: $\exists \gamma > 0$ sao cho $\sum_{i=0}^k |x_{-i} - \bar{x}| < \gamma \Rightarrow \lim_{n\to\infty} x_n = \bar{x}$.
    • Ổn định tiệm cận địa phương: Vừa ổn định địa phương vừa hút địa phương.
    • Hút toàn cục: Mọi nghiệm ${x_n}{n=-k}^\infty$ đều có $\lim{n\to\infty} x_n = \bar{x}$.
    • Ổn định tiệm cận toàn cục: Vừa ổn định địa phương vừa hút toàn cục.
    • Không ổn định: Không thỏa mãn điều kiện ổn định địa phương.

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

Luận án áp dụng kết hợp các phương pháp giải tích định tính thuần túy trong toán học:

  1. Phương pháp tuyến tính hóa: Xây dựng phương trình sai phân tuyến tính hóa và phương trình đặc trưng $\lambda^{k+1} - q_0\lambda^k - \dots - q_k = 0$ tại lân cận điểm cân bằng $\bar{x}$ dựa trên các đạo hàm riêng $q_i = \frac{\partial F}{\partial u_i}(\bar{x}, \dots, \bar{x})$. Đánh giá vị trí các nghiệm đặc trưng $\lambda$ so với hình tròn đơn vị $|\lambda| < 1$.
  2. Phương pháp ma trận Jacobi: Áp dụng cho các hệ phương trình sai phân phi tuyến $X_{n+1} = J_F X_n$, khảo sát phổ giá trị riêng của ma trận Jacobi $J_F(\bar{x}, \bar{y})$ để kết luận tính ổn định tiệm cận hoặc không ổn định của điểm cân bằng hệ.
  3. Nguyên lý so sánh: Sử dụng hàm tăng đơn điệu $F$ và kỹ thuật kẹp dãy $x_n \le y_n \le z_n$ để chứng minh tính bị chặn dưới, bị chặn trên của nghiệm.
  4. Định lý hội tụ: Sử dụng định lý hội tụ của El-Metwally, Grove, Ladas, Voulov (2000) dựa trên tính đơn điệu của hàm số $g(z_1, \dots, z_{k+1})$ và tính duy nhất của nghiệm hệ cực trị $m = M$ để khẳng định tính hút toàn cục.
  5. Phương pháp giải tích tiệm cận (L. Berg, Stević): Xấp xỉ phương trình sai phân rời rạc bằng phương trình vi phân liên tục tương ứng qua khai triển vi phân hình thức, thiết lập thang tiệm cận có chứa hàm logarit $\ln n$, sau đó sử dụng định lý Stević về hàm đơn điệu để chứng minh tính hợp lệ và chặn sai số của dạng tiệm cận.
  6. Phương pháp mô phỏng số: Sử dụng phần mềm MATLAB để tính toán lặp số và vẽ đồ thị quỹ đạo nghiệm nhằm kiểm chứng trực quan các kết quả giải tích.

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.

  • Phương trình sai phân cấp cao: Trình bày định lý ổn định tuyến tính hóa và bổ đề về điều kiện đủ để tất cả các nghiệm của phương trình đặc trưng bậc tùy ý nằm trong hình tròn đơn vị: nếu $\sum_{i=0}^k |q_i| < 1$ thì mọi nghiệm đều có môđun nhỏ hơn 1.
  • Định lý so sánh: Cho $F: I^{k+1} \to I$ tăng theo tất cả các biến, nếu ba dãy thỏa mãn $x_{n+1} \le F(x_n, \dots, x_{n-k})$, $y_{n+1} = F(y_n, \dots, y_{n-k})$, $z_{n+1} \ge F(z_n, \dots, z_{n-k})$ và $x_n \le y_n \le z_n$ với $-k \le n \le 0$ thì $x_n \le y_n \le z_n, \forall n > 0$.
  • Hệ phương trình sai phân phi tuyến: Định nghĩa điểm cân bằng và tính ổn định cho hệ tổng quát $Y_{n+1} = G(n, Y_n)$. Thiết lập ma trận Jacobi cho hệ cấp một hai chiều: $$J_F(x, y) = \begin{pmatrix} \frac{\partial f}{\partial x}(x, y) & \frac{\partial f}{\partial y}(x, y) \ \frac{\partial g}{\partial x}(x, y) & \frac{\partial g}{\partial y}(x, y) \end{pmatrix}$$ và hệ cấp hai hai chiều $x_{n+1} = f(x_n, x_{n-1}, y_n, y_{n-1}), y_{n+1} = g(x_n, x_{n-1}, y_n, y_{n-1})$ với ma trận Jacobi cấp $4 \times 4$: $$J_F(x, y) = \begin{pmatrix} \frac{\partial f}{\partial x_1} & \frac{\partial f}{\partial x_2} & \frac{\partial f}{\partial y_1} & \frac{\partial f}{\partial y_2} \ 1 & 0 & 0 & 0 \ \frac{\partial g}{\partial x_1} & \frac{\partial g}{\partial x_2} & \frac{\partial g}{\partial y_1} & \frac{\partial g}{\partial y_2} \ 0 & 0 & 1 & 0 \end{pmatrix}$$

Chương 2: Sự ổn định của điểm cân bằng của một số dạng phương trình sai phân cấp cao dạng phân thức

Chương 2 bao gồm ba nội dung nghiên cứu chính tương ứng với các lớp phương trình sai phân cấp hai, cấp ba và cấp bốn.

1. Sự ổn định tiệm cận của nghiệm của hai lớp phương trình sai phân cấp hai dạng phân thức

  • Lớp phương trình thứ nhất (2.1): $$x_{n+1} = ax_{n-1} - \frac{bx_n^2 + cx_{n-1}^2}{1 - dx_n + ex_{n-1}}, \quad n = 0, 1, 2, \dots$$ với $a, b, c, d, e > 0$ và giá trị ban đầu $x_{-1}, x_0 > 0$. Phương trình có hai điểm cân bằng là $\bar{x}_1 = 0$ và $\bar{x}_2 = \frac{a-1}{b+c-(e-d)(a-1)}$. Điểm cân bằng $\bar{x}_2$ dương khi $a > 1$ và $e < \frac{b+c}{a-1} + d$.

    • Định lý 2.1: Khi $a \ge 2, c \le \min(ae, 3b+d)$ và $e < \frac{b+c}{a-1}$, điểm cân bằng dương $\bar{x}_2$ là không ổn định. Chứng minh dựa trên việc xác định phương trình đặc trưng $\lambda^2 + p\lambda + q = 0$ với: $$p = \frac{(a-1)[2b + d(a-1)]}{b+c} > 0, \quad q = -\frac{ab + e^2 + (c-ae)(2-a)}{b+c} < 0$$ và chỉ ra phương trình có nghiệm $\lambda_2 > 1$.
    • Định lý 2.2: Khi $a < 1$, điểm cân bằng không $\bar{x}_1 = 0$ là ổn định tiệm cận địa phương do phương trình đặc trưng $\lambda^2 - a = 0$ có $|p| = 0 < 1 - a < 2$.

  • Lớp phương trình thứ hai (2.2): $$x_{n+1} = \frac{\alpha x_{n-1}}{\beta x_n + x_{n-1}} + \frac{x_{n-1}}{x_n}, \quad n = 0, 1, 2, \dots$$ với $\alpha, \beta > 0$ và $x_{-1}, x_0 > 0$. Phương trình có điểm cân bằng duy nhất $\bar{x} = \frac{\alpha}{\beta+1} + 1$. Phương trình đặc trưng có dạng: $$\lambda^2 + \frac{\alpha\beta + (\beta+1)^2}{(\beta+1)(\alpha+\beta+1)}\lambda - \frac{\alpha\beta + (\beta+1)^2}{(\beta+1)(\alpha+\beta+1)} = 0$$

    • Định lý 2.3:
      1. Điểm cân bằng $\bar{x}$ là ổn định tiệm cận địa phương nếu $0 < \beta < 1$ và $\alpha > \frac{(\beta+1)^2}{1-\beta}$.
      2. Điểm cân bằng $\bar{x}$ là điểm yên ngựa (và do đó không ổn định) nếu $\alpha < \frac{(\beta+1)^2}{1-\beta}$ hoặc khi $\beta > \frac{-(\alpha+2) + \sqrt{\alpha^2+8\alpha}}{2}$.
Lớp phương trình Điểm cân bằng $\bar{x}$ Điều kiện tham số Tính chất ổn định
Phương trình (2.1) $\bar{x} = 0$ $a < 1$ Ổn định tiệm cận địa phương
Phương trình (2.1) $\bar{x} = \frac{a-1}{b+c-(e-d)(a-1)}$ $a \ge 2, c \le \min(ae, 3b+d), e < \frac{b+c}{a-1}$ Không ổn định
Phương trình (2.2) $\bar{x} = \frac{\alpha}{\beta+1} + 1$ $0 < \beta < 1, \alpha > \frac{(\beta+1)^2}{1-\beta}$ Ổn định tiệm cận địa phương
Phương trình (2.2) $\bar{x} = \frac{\alpha}{\beta+1} + 1$ $\alpha < \frac{(\beta+1)^2}{1-\beta}$ hoặc $\beta > \frac{-(\alpha+2)+\sqrt{\alpha^2+8\alpha}}{2}$ Điểm yên ngựa (không ổn định)

2. Dạng tiệm cận của nghiệm của phương trình sai phân cấp ba dạng phân thức

Khảo sát phương trình sai phân cấp ba (2.25): $$x_n = \frac{x_{n-3} - (x_n + x_{n-1})^3}{1 + x_n x_{n-1} + x_n x_{n-2} + x_{n-1} x_{n-2}}, \quad n = 0, 1, 2, \dots$$ với $x_{-3}, x_{-2}, x_{-1} \in (0, \infty)$. Điểm cân bằng của phương trình là $\bar{x} = 0$.

Bằng cách thiết lập phương trình vi phân xấp xỉ liên tục bậc một $11x^3 = -3x'$ cho nghiệm tiệm cận cấp một $x \sim \sqrt{\frac{3}{22t}}$, sau đó khai triển xấp xỉ Taylor cấp cao hơn dẫn đến phương trình vi phân $x^2(22x + 63x') = -6x'$, giải ra quan hệ $3\ln x - \frac{63}{x^2} = -22t + C$. Tác giả đề xuất dạng tiệm cận năm số hạng của nghiệm: $$\phi_n = \frac{1}{\sqrt{n}}\left(a + \frac{b\ln n + c}{n} + \frac{d\ln^2 n + e\ln n}{n^2}\right)$$

Bằng phương pháp đồng nhất hệ số đạo hàm tiệm cận, luận án xác định chính xác các hệ số: $$a = \frac{\sqrt{66}}{22}, \quad b = \frac{63\sqrt{66}}{8 \times 22^2}, \quad c = \frac{63\sqrt{66}}{8 \times 22^2}, \quad d = \frac{3 \times 63^2\sqrt{66}}{128 \times 22^3}, \quad e = \frac{63^2\sqrt{66}}{8 \times 22^3}$$

  • Định lý 2.4: Phương trình (2.25) tồn tại nghiệm có dạng tiệm cận $\phi_n$ thỏa mãn với mọi $\varepsilon > 0$, tồn tại $n(\varepsilon) \ge 0$ để: $$\phi_n - \varepsilon \frac{\ln n}{n^2\sqrt{n}} \le x_n \le \phi_n + \varepsilon \frac{\ln n}{n^2\sqrt{n}}, \quad \forall n \ge n(\varepsilon)$$ và nghiệm này hội tụ về 0 khi $n \to \infty$.

3. Sự ổn định toàn cục của nghiệm của phương trình sai phân cấp bốn dạng phân thức

Khảo sát lớp phương trình sai phân cấp bốn tổng quát (2.38): $$x_{n+1} = \frac{g(x_n, x_{n-1}, x_{n-2}, x_{n-3})x_{n-3}^{\alpha+1} + f(x_n, x_{n-1}, x_{n-2}, x_{n-3})}{g(x_n, x_{n-1}, x_{n-2}, x_{n-3})x_{n-3}^\alpha + f(x_n, x_{n-1}, x_{n-2}, x_{n-3})}$$ với $\alpha \in [0, \infty)$ và $f, g: (0, \infty)^4 \to (0, \infty)$ là các hàm dương khả vi. Phương trình có duy nhất điểm cân bằng $\bar{x} = 1$.

  • Bổ đề 2.1: Với mọi nghiệm dương ${x_n}{n=-3}^\infty$, ta có: $$(x{n+1} - 1)(x_{n-3} - 1) \ge 0 \quad \text{và} \quad (x_{n+1} - x_{n-3})(1 - x_{n-3}) \ge 0$$
  • Định lý 2.6: Điểm cân bằng $\bar{x} = 1$ là ổn định tiệm cận toàn cục. Phương trình tuyến tính hóa tại $\bar{x} = 1$ cho $y_{n+1} = 0$, chứng minh tính ổn định địa phương. Tính hút toàn cục được chứng minh bằng cách chia dãy nghiệm thành 4 dãy con ${x_{4n+k}}_{n=0}^\infty$ ($k = 0, 1, 2, 3$). Các dãy con này đều đơn điệu, bị chặn và cùng hội tụ về giới hạn $a_k = 1$.

Chương 3: Sự ổn định tiệm cận của điểm cân bằng của một số hệ phương trình sai phân dạng phân thức

Chương 3 tập trung mở rộng nghiên cứu tính ổn định sang các hệ phương trình sai phân phân thức nhiều chiều:

  • Mục 3.1: Hệ phương trình sai phân cấp một dạng phân thức: Thiết lập bài toán cho hệ hai phương trình sai phân phân thức cấp một; chứng minh tính bị chặn của nghiệm; phân tích phương trình đặc trưng của ma trận Jacobi tại điểm cân bằng dương để tìm điều kiện ổn định tiệm cận địa phương và minh họa số trên MATLAB.
  • Mục 3.2: Hệ phương trình sai phân cấp hai dạng phân thức: Thiết lập bài toán cho hệ phương trình sai phân phân thức cấp hai; phân tích tính bị chặn của nghiệm; sử dụng ma trận Jacobi cấp $4 \times 4$ để khảo sát tính ổn định của điểm cân bằng dương và minh họa số bằng đồ thị MATLAB.

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 khoa học lý thuyết và thực tiễn:

Đóng góp về lý luận

  1. Thiết lập các điều kiện đủ tường minh về tham số cho sự ổn định tiệm cận địa phương của điểm cân bằng không, cũng như tính không ổn định của điểm cân bằng dương đối với lớp phương trình sai phân cấp hai dạng $x_{n+1} = ax_{n-1} - \frac{bx_n^2 + cx_{n-1}^2}{1 - dx_n + ex_{n-1}}$.
  2. Xác định chính xác miền tham số $(\alpha, \beta)$ để điểm cân bằng của phương trình $x_{n+1} = \frac{\alpha x_{n-1}}{\beta x_n + x_{n-1}} + \frac{x_{n-1}}{x_n}$ đạt trạng thái ổn định tiệm cận địa phương hoặc trở thành điểm yên ngựa.
  3. Xây dựng thành công công thức khai triển tiệm cận bậc cao năm số hạng chứa hàm logarit tự nhiên $\ln n$ và $\ln^2 n$ cho nghiệm hội tụ về không của phương trình sai phân cấp ba phi tuyến bậc ba, đồng thời chứng minh chặt chẽ sự tồn tại của nghiệm thực tế kẹp trong lân cận tiệm cận đó bằng định lý Stević.
  4. Mở rộng kết quả ổn định tiệm cận toàn cục của điểm cân bằng $\bar{x} = 1$ cho lớp phương trình sai phân phân thức cấp bốn chứa hai hàm số dương khả vi tùy ý $f, g$ và tham số mũ $\alpha \ge 0$, tổng quát hóa các kết quả trước đó của G. Ladas, Nesemann, Zhu, Li và Agarwal.
  5. Thiết lập điều kiện bị chặn và ổn định tiệm cận của điểm cân bằng dương cho hai lớp hệ phương trình sai phân dạng phân thức cấp một và cấp hai.

Đóng góp về thực tiễn và ứng dụng

  • Cung cấp các công thức giải tích trực tiếp giúp xác định hành vi dài hạn của các hệ động lực rời rạc dạng phân thức mà không cần giải phương trình tường minh.
  • Xây dựng các thuật toán tính toán lặp và kịch bản mô phỏng số trên phần mềm MATLAB để trực quan hóa quỹ đạo nghiệm của phương trình và hệ phương trình sai phân.
  • Các kết quả công bố trong 06 công trình khoa học (05 bài báo viết chung và 01 bài báo độc lập của tác giả) được báo cáo tại các seminar chuyên ngành của Trường Đại học Khoa học Tự nhiên – ĐHQG Hà Nội và Trường Đại học Sư phạm Kỹ thuật Hưng Yên.

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

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

  • Hạn chế: Luận án tập trung khảo sát các lớp phương trình và hệ phương trình cụ thể với giả thiết các tham số và điều kiện ban đầu là số thực dương hoặc không âm; chưa khảo sát trường hợp nghiệm phức hoặc tham số đổi dấu; chưa mở rộng phân tích các hiện tượng động lực học phi tuyến phức tạp khác như tính phân nhánh (bifurcation), tính tuần hoàn tổng quát hay tính hỗn loạn (chaos).
  • Hướng nghiên cứu tiếp:
    • Tiếp tục mở rộng nghiên cứu tính ổn định, tính bị chặn và tính tuần hoàn cho các lớp phương trình sai phân dạng phân thức có cấp cao hơn ($k > 4$).
    • Mở rộng khảo sát các hệ phương trình sai phân phân thức $m$ chiều ($m \ge 3$) với cấu trúc trễ phức tạp.
    • Nghiên cứu các lớp phương trình sai phân ngẫu nhiên hoặc phương trình sai phân chứa tham số biến thiên theo thời gian $n$.

Giá trị tham khảo

Luận án có giá trị tham khảo chuyên môn trong các lĩnh vực:

  • Nghiên cứu sinh và học viên cao học chuyên ngành Toán giải tích, Phương trình vi phân và Hệ động lực: Tham khảo phương pháp nghiên cứu định tính, kỹ thuật tuyến tính hóa, kỹ thuật phân tích ma trận Jacobi cho hệ sai phân nhiều chiều và phương pháp xây dựng khai triển tiệm cận nghiệm rời rạc của L. Berg và Stević.
  • Giảng viên và nhà nghiên cứu toán học: Sử dụng làm tài liệu tham khảo cho các bài giảng chuyên đề về phương trình sai phân phi tuyến, lý thuyết ổn định Lyapunov rời rạc và các mô hình toán sinh học.
  • Các nhà mô hình hóa ứng dụng: Tham khảo cấu trúc các mô hình phân thức và điều kiện ổn định cân bằng để áp dụng vào việc mô phỏng rời rạc hóa các bài toán sinh thái học, dịch tễ học và kinh tế học.
  • Nội dung nổi bật đáng tham khảo: Kỹ thuật tìm hệ số tiệm cận chứa logarit ở Mục 2.2 và phương pháp phân rã dãy con bốn thành phần để chứng minh ổn định toàn cục ở Mục 2.3.

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

1. Luận án nghiên cứu những lớp phương trình sai phân cụ thể nào trong Chương 2?

Luận án nghiên cứu bốn lớp phương trình sai phân cụ thể trong Chương 2:

  • Lớp cấp hai chứa bậc hai: $x_{n+1} = ax_{n-1} - \frac{bx_n^2 + cx_{n-1}^2}{1 - dx_n + ex_{n-1}}$ ($a, b, c, d, e > 0$).
  • Lớp cấp hai hai phân thức: $x_{n+1} = \frac{\alpha x_{n-1}}{\beta x_n + x_{n-1}} + \frac{x_{n-1}}{x_n}$ ($\alpha, \beta > 0$).
  • Lớp cấp ba phi tuyến: $x_n = \frac{x_{n-3} - (x_n + x_{n-1})^3}{1 + x_n x_{n-1} + x_n x_{n-2} + x_{n-1} x_{n-2}}$.
  • Lớp cấp bốn tổng quát: $x_{n+1} = \frac{g(x_n, x_{n-1}, x_{n-2}, x_{n-3})x_{n-3}^{\alpha+1} + f(x_n, x_{n-1}, x_{n-2}, x_{n-3})}{g(x_n, x_{n-1}, x_{n-2}, x_{n-3})x_{n-3}^\alpha + f(x_n, x_{n-1}, x_{n-2}, x_{n-3})}$.

2. Điểm cân bằng của phương trình cấp bốn (2.38) có giá trị bằng bao nhiêu và có tính chất ổn định gì?

Phương trình cấp bốn (2.38) có duy nhất một điểm cân bằng là $\bar{x} = 1$. Điểm cân bằng này có tính chất ổn định tiệm cận toàn cục với mọi tham số $\alpha \in [0, \infty)$ và các hàm dương khả vi $f, g$.

3. Dạng tiệm cận của nghiệm phương trình cấp ba (2.25) chứa các số hạng nào?

Dạng tiệm cận năm số hạng của nghiệm phương trình (2.25) được xác định bởi công thức: $$\phi_n = \frac{1}{\sqrt{n}}\left(a + \frac{b\ln n + c}{n} + \frac{d\ln^2 n + e\ln n}{n^2}\right)$$ với các hệ số thực xác định chính xác: $a = \frac{\sqrt{66}}{22}$, $b = c = \frac{63\sqrt{66}}{8 \times 22^2}$, $d = \frac{3 \times 63^2\sqrt{66}}{128 \times 22^3}$, $e = \frac{63^2\sqrt{66}}{8 \times 22^3}$.

4. Khi nào điểm cân bằng dương của phương trình cấp hai (2.2) trở thành điểm yên ngựa?

Điểm cân bằng $\bar{x} = \frac{\alpha}{\beta+1} + 1$ của phương trình (2.2) là điểm yên ngựa (và do đó không ổn định) khi thỏa mãn một trong hai điều kiện: $\alpha < \frac{(\beta+1)^2}{1-\beta}$ (với $0 < \beta < 1$) hoặc $\beta > \frac{-(\alpha+2) + \sqrt{\alpha^2+8\alpha}}{2}$.

5. Phương pháp nào được sử dụng để chứng minh tính hút toàn cục của phương trình cấp bốn trong luận án?

Luận án sử dụng phương pháp phân tích dãy nghiệm ${x_n}{n=-3}^\infty$ thành bốn dãy con rời nhau ${x{4n}}$, ${x_{4n+1}}$, ${x_{4n+2}}$, ${x_{4n+3}}$. Dựa vào Bổ đề 2.1, chứng minh cả bốn dãy con này đều là các dãy đơn điệu và bị chặn, từ đó suy ra giới hạn hữu hạn của chúng và chứng minh cả bốn giới hạn đều bằng 1.

Kết luận

Luận án tiến sĩ của tác giả Trần Hồng Thái là công trình nghiên cứu chuyên sâu thuộc chuyên ngành Toán giải tích về lý thuyết định tính của phương trình và hệ phương trình sai phân phi tuyến dạng phân thức. Tác giả đã giải quyết trọn vẹn các bài toán về sự ổn định địa phương, tính chất yên ngựa, tính ổn định tiệm cận toàn cục và dạng biểu diễn tiệm cận logarit của nghiệm cho nhiều lớp phương trình sai phân cấp hai, cấp ba, cấp bốn và hệ phương trình sai phân tương ứng. Các kết quả nghiên cứu được chứng minh chặt chẽ bằng các công cụ giải tích toán học thuần túy, đóng góp các tiêu chuẩn định tính mới và mở rộng nhiều công trình tiền đề trong lý thuyết phương trình sai phân rời rạc.