Tổng quan luận án

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

Trong mô hình hóa toán học đối với các hiện tượng thực tế trên thang thời gian rời rạc, các hệ thống thường được biểu diễn dưới dạng phương trình sai phân. Trong nhiều lĩnh vực ứng dụng như điều khiển, kinh tế học (tiêu biểu là mô hình Leontief động đa ngành) và sinh học, các hệ thống phát sinh ở dạng ẩn mà số hạng bậc cao nhất không thể giải tường minh, dẫn đến phương trình sai phân suy biến (singular difference equations - SDEs), hay còn gọi là phương trình sai phân đại số hoặc hệ thống mô tả rời rạc (discrete-time descriptor systems). Khác với phương trình sai phân tường minh thông thường, ma trận hệ số đứng trước bước thời gian kế tiếp $E_n$ là ma trận suy biến ($\text{rank} E_n = r < d$), gây ra sự xuất hiện đồng thời của thành phần động lực học và các ràng buộc đại số, khiến việc phân tích tính giải được và tính ổn định trở nên phức tạp.

Khoảng trống nghiên cứu được tác giả xác định cụ thể như sau:

  • Các kết quả nghiên cứu trước đây về phương trình sai phân suy biến chủ yếu tập trung vào trường hợp hệ dừng (hệ số hằng) thông qua dạng chuẩn tắc Weierstrass-Kronecker hoặc giải thức ma trận nghịch đảo Drazin (S. Luenberger, S. L. Campbell).
  • Đối với hệ suy biến có hệ số biến thiên theo thời gian (time-varying), các kết quả còn rất hạn chế. Các nghiên cứu hiện có chủ yếu tiếp cận theo phương pháp hàm Lyapunov cho phương trình bán tuyến tính (Phạm Kỳ Anh và cộng sự) hoặc chỉ xét hệ cấp một với các dạng nhiễu hẹp, chưa khảo sát toàn diện ảnh hưởng của nhiễu hai phía (xuất hiện ở cả hai vế).
  • Định lý kiểu Bohl-Perron thiết lập mối quan hệ giữa tính bị chặn của nghiệm hệ không thuần nhất với tính ổn định mũ/ổn định đều của hệ thuần nhất tương ứng vốn đã phát triển hoàn chỉnh trên phương trình vi phân và sai phân thông thường (Bohl, Perron, Daleckii, Krein, Aulbach, Nguyễn Văn Minh, Pituk), nhưng chưa được mở rộng đầy đủ cho các hệ phương trình sai phân suy biến tuyến tính cấp một và cấp hai có hệ số biến thiên.
  • Khái niệm số mũ Bohl (Bohl exponent) và tính vững của số mũ Bohl đối với hệ rời rạc suy biến chưa từng được nghiên cứu trong các y văn trước đó.
  • Đối với phương trình sai phân suy biến cấp hai (và bậc cao) có hệ số biến thiên, các vấn đề về tính giải được của bài toán giá trị ban đầu, tính tương thích của điều kiện ban đầu, tiêu chuẩn ổn định mũ và bán kính ổn định vững dưới tác động của nhiễu có cấu trúc vẫn là một bài toán mở chưa được giải quyết.

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. Phát triển phương pháp tiếp cận bằng phép chiếu (projector-based approach) để phân tích tính ổn định và tính ổn định dưới tác động của nhiễu một phía và nhiễu hai phía cho phương trình sai phân suy biến tuyến tính cấp một có hệ số biến thiên.
  2. Thiết lập các định lý dạng Bohl-Perron cho phương trình sai phân suy biến tuyến tính cấp một và cấp hai; xây dựng khái niệm số mũ Bohl và khảo sát tính chất nhạy của số mũ Bohl đối với nhiễu hệ số.
  3. Nghiên cứu tính giải được của bài toán giá trị ban đầu và xác định cấu trúc tường minh của tập điều kiện ban đầu tương thích cho phương trình sai phân suy biến tuyến tính cấp hai có hệ số biến thiên thông qua chỉ số lạ (strangeness index).
  4. Thiết lập tiêu chuẩn ổn định mũ thông qua nguyên lý so sánh và ước lượng cận trên tường minh cho nhiễu có cấu trúc để hệ sai phân suy biến cấp hai bảo toàn chỉ số và tính ổn định mũ (ổn định vững).

Câu hỏi nghiên cứu

Luận án tập trung giải quyết các câu hỏi khoa học:

  1. Trong điều kiện nào thì phương trình sai phân suy biến tuyến tính cấp một có hệ số biến thiên bảo toàn được tính ổn định đều và ổn định mũ khi ma trận hệ số chịu tác động của nhiễu một phía hoặc hai phía?
  2. Có thể đặc trưng hóa tính ổn định mũ và ổn định đều của hệ sai phân suy biến thuần nhất thông qua tính bị chặn của nghiệm trong các không gian dãy $l_p$ của hệ không thuần nhất tương ứng (nguyên lý Bohl-Perron) hay không?
  3. Khái niệm số mũ Bohl được định nghĩa như thế nào cho hệ sai phân suy biến, và số mũ này có tính chất bền vững (robustness) trước các nhiễu cho phép hay không?
  4. Làm thế nào để giải trừ liên kết (decouple) một phương trình sai phân suy biến tuyến tính cấp hai có hệ số biến thiên về dạng chính quy không có chỉ số lạ, và làm thế nào để xác định bán kính ổn định vững dưới tác động của nhiễu có cấu trúc?

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

  • Đối tượng nghiên cứu:
    • Phương trình sai phân suy biến tuyến tính cấp một có hệ số biến thiên: $E_n y(n+1) = A_n y(n) + q_n$, $n \in \mathbb{N}(n_0)$, với $E_n, A_n \in \mathbb{K}^{d \times d}$, $\text{rank} E_n = r < d$.
    • Phương trình sai phân suy biến tuyến tính cấp hai có hệ số biến thiên: $A_n x(n+2) + B_n x(n+1) + C_n x(n) = f_n$, $n \in \mathbb{N}(n_0)$, với $A_n, B_n, C_n \in \mathbb{C}^{d \times d}$, $\text{rank} A_n \equiv d_1 < d$.
  • Phạm vi nghiên cứu:
    • Về cấu trúc toán học: Giới hạn ở các hệ có chỉ số kéo theo bằng 1 (tractability index-1) đối với hệ cấp một và hệ không có chỉ số lạ (strangeness-free, $\mu = 0$) đối với hệ cấp hai.
    • Về không gian và thời gian: Thời gian rời rạc $n \in \mathbb{N}(n_0) = {n \in \mathbb{Z} : n \ge n_0}$, không gian trạng thái hữu hạn chiều $\mathbb{K}^d$ ($\mathbb{R}^d$ hoặc $\mathbb{C}^d$).

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 mà tác giả điểm lại

Tác giả đã hệ thống hóa các hướng nghiên cứu nền tảng liên quan trực tiếp đến đề tài:

  1. Nghiên cứu hệ sai phân suy biến dừng và giải thức giải tích:

    • S. Luenberger (1977, 1979) [52, 54] là những công trình đầu tiên nghiên cứu hệ sai phân suy biến thông qua dạng phân tích Weierstrass-Kronecker cho chùm ma trận chính quy ${E, A}$, ứng dụng trực tiếp vào mô hình cân đối liên ngành Leontief động $x(k) = Ax(k) + B[x(k+1) - x(k)] + d(k)$.
    • S. L. Campbell (1980, 1982) [17] phát triển phương pháp giải phương trình sai phân suy biến bằng nghịch đảo Drazin cho hệ dừng và đưa ra điều kiện giải được cho hệ không dừng dưới giả thiết hạng ma trận không đổi.

  2. Phương pháp phép chiếu và chỉ số kéo theo:

    • R. März và các cộng sự (từ những năm 1980) [48] xây dựng lý thuyết giải tích dựa trên toán tử chiếu (projector-based analysis) cho phương trình vi phân - đại số (DAEs).
    • Phạm Kỳ Anh, Vũ Hoàng Linh và các cộng sự [3, 4, 5, 6] mở rộng khái niệm chỉ số kéo theo cấp 1 (tractability index-1) cho phương trình sai phân đại số phi tuyến [5], thiết lập lý thuyết Floquet và định lý giản lược Lyapunov cho hệ tuần hoàn [6], và áp dụng phương pháp hàm Lyapunov cho hệ sai phân tựa tuyến tính [4].
    • Ngô Thị Thanh Nga (2015) [59] bước đầu nghiên cứu hành vi tiệm cận của nghiệm phương trình sai phân suy biến chịu tác động của nhiễu.

  3. Lý thuyết Bohl-Perron và tính bị chặn của nghiệm:

    • P. Bohl (1904) và O. Perron (1930) khởi xướng định lý Bohl-Perron cho phương trình vi phân thường: tính chất "đầu vào bị chặn - trạng thái bị chặn" (bounded input - bounded state) tương đương với tính ổn định mũ của hệ thuần nhất.
    • Daleckii và Krein (1974) [22] tổng quát hóa định lý Bohl-Perron trên không gian Banach.
    • B. Aulbach và Nguyễn Văn Minh (1996) [8], cùng với P. Pituk (2000) [70], mở rộng định lý Bohl-Perron cho phương trình sai phân chính quy trong không gian Banach với các lớp không gian dãy $l_p$ ($1 \le p \le \infty$). E. Braverman (2004, 2007) [9, 12] tiếp tục mở rộng cho phương trình sai phân có trễ.
    • Nguyễn Hữu Dư và Vũ Hoàng Linh (2010, 2011) [10, 19] thiết lập định lý Bohl-Perron cho phương trình vi phân - đại số (DAEs).

  4. Lý thuyết ổn định vững và bán kính ổn định:

    • D. Hinrichsen và A. J. Pritchard (1986, 1989) [40, 41] xây dựng lý thuyết bán kính ổn định (stability radius) cho hệ phương trình vi phân liên tục chịu nhiễu có cấu trúc $\dot{x} = (A + B\Delta C)x$. L. Qiu và các cộng sự (1995) thiết lập công thức tính bán kính ổn định thực.
    • R. Byers, N. K. Nichols (1993) [16] và L. Qiu, E. J. Davison (1998) [71] tính bán kính ổn định phức cho hệ DAEs dừng chỉ số 1: $E\dot{x} = Ax$.
    • Nguyễn Hữu Dư và Vũ Hoàng Linh (2006, 2013) [24, 25] nghiên cứu bán kính ổn định cho hệ DAEs suy biến kỳ dị; B. Jacob (2001) [43] áp dụng lý thuyết toán tử tìm công thức giải tích cho hệ biến thiên theo thời gian; Nguyễn Hữu Dư, Vũ Hoàng Linh và Ngô Thị Thanh Nga (2013) [61] mở rộng cho hệ sai phân suy biến biến thiên.
    • V. Mehrmann và Đỗ Đức Thuận (2014) [56] nghiên cứu bán kính ổn định của hệ sai phân suy biến bậc cao hệ số hằng dựa trên kỹ thuật giải trừ liên kết ma trận.

  5. Phương pháp tiếp cận chỉ số lạ (Strangeness-index):

    • P. Kunkel và V. Mehrmann (từ những năm 1990) phát triển lý thuyết chỉ số lạ cho DAEs.
    • T. Brüll (2006) [13] chuyển giao toàn diện lý thuyết chỉ số lạ sang hệ thống mô tả rời rạc (discrete-time descriptor systems), phân loại các bất biến cục bộ và quy trình hồi quy rút gọn hệ thống.

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

Luận án tập trung xử lý dứt điểm các giới hạn lý thuyết:

  • Xây dựng cặp toán tử chiếu chuẩn tắc duy nhất nhằm khử liên kết phương trình sai phân suy biến cấp một có hệ số biến thiên, khắc phục việc mất tính duy nhất khi chọn phép chiếu thông thường.
  • Thiết lập tiêu chuẩn ổn định khi hệ chịu đồng thời cả nhiễu ở ma trận trạng thái $A_n$ lẫn ma trận dẫn đầu $E_n$ (nhiễu hai phía).
  • Chuyển tiếp định lý Bohl-Perron sang phương trình sai phân suy biến, xây dựng lý thuyết số mũ Bohl đầu tiên cho hệ rời rạc suy biến.
  • Khảo sát toàn diện phương trình sai phân suy biến cấp hai có hệ số biến thiên bằng cách kết hợp dạng chuẩn tắc không có chỉ số lạ của T. Brüll với nguyên lý so sánh ma trận để định lượng bán kính ổn định vững.

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 vận dụng các cơ sở lý thuyết toán học chuyên sâu:

  1. Lý thuyết chỉ số kéo theo (Tractability index-1):

    • Đặt $N_n := \ker E_n$, $Q_n$ là toán tử chiếu lên $N_n$, $P_n := I - Q_n$.
    • Giả thiết tồn tại dãy toán tử khả nghịch $T_n \in GL(\mathbb{K}^d)$ sao cho $T_n|{N_n}$ là đẳng cấu từ $N_n$ lên $N{n-1}$.
    • Các ma trận và không gian bổ trợ: $$G_n := E_n - A_n T_n Q_n, \quad S_n := {z \in \mathbb{K}^d : A_n z \in \text{im} E_n}$$
    • Hệ có chỉ số kéo theo bằng 1 khi và chỉ khi $\text{rank} E_n = r$ không đổi và $N_{n-1} \cap S_n = {0}$, tương đương với tính khả nghịch của ma trận $G_n$.
    • Cặp phép chiếu chuẩn tắc duy nhất (canonical projector pair): $$\widetilde{Q}_{n-1} := -T_n Q_n G_n^{-1} A_n, \quad \widetilde{P}_n := I - \widetilde{Q}_n$$

  2. Lý thuyết chỉ số lạ (Strangeness-index):

    • Áp dụng phân tích chùm ma trận cho cặp $(E_k, A_k)$ với các bất biến cục bộ của T. Brüll: hạng tiến $r_f = \text{rank}(E)$, hạng lùi $r_b = \text{rank}(A)$, $h_f = \text{rank}(Z^H A)$ (với $Z$ là cơ sở của $\ker E^H$), chỉ số lạ $s = h_f - \min(h_f, n - r_f)$.
    • Quá trình rút gọn đệ quy xác định chỉ số lạ $\mu = \min{i \in \mathbb{N} : s_i = 0}$. Trường hợp $\mu = 0$ là hệ không có chỉ số lạ (strangeness-free).

  3. Lý thuyết số mũ Bohl (Bohl exponent):

    • Mở rộng khái niệm số mũ Bohl từ hệ phương trình vi phân và sai phân thông thường của P. Bohl [15] và Daleckii - Krein [22] sang hệ sai phân suy biến để đo lường cận trên chính xác của tốc độ suy giảm nghiệm.

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

Luận án thuần túy sử dụng các phương pháp phân tích toán học lý thuyết (định tính):

  • Kỹ thuật khử liên kết giải tích (Decoupling technique): Sử dụng ma trận $G_n^{-1}$ và cặp toán tử chiếu chuẩn tắc ${\widetilde{P}n, \widetilde{Q}n}$ để phân rã không gian nghiệm thành hai thành phần trực giao: thành phần động học $\widetilde{P}{n-1} y(n)$ (thỏa mãn một phương trình sai phân chính quy) và thành phần đại số $\widetilde{Q}{n-1} y(n)$ (xác định đơn trị qua số hạng không thuần nhất).
  • Phương pháp toán tử Cauchy: Xây dựng toán tử tiến hóa Cauchy $\Phi(n, m)$ tương thích với toán tử chiếu chuẩn tắc để biểu diễn công thức biến thiên hằng số cho nghiệm của hệ suy biến.
  • Nguyên lý so sánh (Comparison principle): Thiết lập bất đẳng thức ma trận và hệ so sánh dương để đánh giá chuẩn vector nghiệm của hệ cấp hai thông qua nghiệm của phương trình sai phân vô hướng hoặc phương trình sai phân cấp một tương đương.
  • Phương pháp giải tích hàm và lý thuyết toán tử: Vận dụng nguyên lý bị chặn đều (Uniform Boundedness Principle), bất đẳng thức Young dạng rời rạc cho phép cuộn dãy số, và bất đẳng thức Gronwall rời rạc để chặn nghiệm.
  • Phương pháp phân tích ma trận: Ứng dụng phân tích giá trị kỳ dị (Singular Value Decomposition - SVD), nghịch đảo Drazin, và các phép biến đổi tương đương toàn cục bảo toàn cấu trúc ma trận.

Nguồn tư liệu và dữ liệu nghiên cứu

Luận án sử dụng ngữ liệu và mô hình cấu trúc toán học trừu tượng:

  • Dữ liệu tham chiếu mô hình kinh tế vĩ mô: Mô hình cân đối động Leontief [52, 54] $Bx(k+1) = (I - A + B)x(k) - d(k)$ với ma trận vốn $B$ suy biến.
  • Hệ thống các định lý, bổ đề và bài báo chuyên ngành xuất bản quốc tế về DAEs và SDEs.

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

Chapter 1: Preliminaries

Chương này trình bày các định nghĩa, khái niệm nền tảng và các kết quả bổ trợ được sử dụng xuyên suốt luận án.

Nội dung chi tiết:

  • Tiếp cận hệ sai phân suy biến tuyến tính bằng chỉ số kéo theo (Mục 1.1): Khảo sát hệ cấp một $E_n y(n+1) = A_n y(n) + q_n$. Thiết lập Lemma 1.1 chứng minh sự tương đương giữa tính khả nghịch của $G_n = E_n - A_n T_n Q_n$ và điều kiện tách không gian $N_{n-1} \oplus S_n = \mathbb{K}^d$. Định nghĩa hệ có chỉ số kéo theo cấp 1. Xây dựng cặp toán tử chiếu chuẩn tắc duy nhất $\widetilde{P}_n, \widetilde{Q}_n$ (Lemma 1.3 và 1.4) và khử liên kết hệ thành: $$\begin{cases} \widetilde{P}_n y(n+1) = \widetilde{P}_n \widetilde{G}n^{-1} A_n \widetilde{P}{n-1} y(n) + \widetilde{P}_n \widetilde{G}n^{-1} q_n \ \widetilde{Q}{n-1} y(n) = T_n \widetilde{Q}n \widetilde{G}n^{-1} q_n \end{cases}$$ Chứng minh bài toán Cauchy với điều kiện ban đầu tương thích $P{n_0-1}(y(n_0) - y_0) = 0$ có nghiệm duy nhất được biểu diễn qua toán tử Cauchy $\Phi(n, m) = \prod{k=n-1}^m \widetilde{P}_k \widetilde{G}_k^{-1} A_k$.

  • Tiếp cận hệ sai phân suy biến tuyến tính bằng chỉ số lạ (Mục 1.2): Trình bày lý thuyết của T. Brüll về quan hệ tương đương toàn cục giữa các cặp dãy ma trận ${(E_k, A_k)} \sim {(\widetilde{E}_k, \widetilde{A}_k)}$ qua các dãy ma trận khả nghịch $P_k, Q_k$. Liệt kê hệ thống bất biến địa phương $(r_f, r_b, h_f, h_b, c, a, s, d, u, v)$, quy trình hồi quy ma trận và định nghĩa chỉ số lạ $\mu$ theo công thức (1.22). Chứng minh tính tương đương giữa hai cách tiếp cận: hệ có chỉ số kéo theo cấp 1 tương đương với hệ không có chỉ số lạ ($\mu = 0$). Đồng thời, tóm tắt các kết quả của Mehrmann và Thuận (2014) [56] về bán kính ổn định của hệ sai phân suy biến dừng cấp hai.

  • Các kết quả bổ trợ khác (Mục 1.3): Thiết lập bổ đề Gronwall rời rạc (Lemma 1.1), bổ đề nghịch đảo ma trận Sherman-Morrison-Woodbury dạng đại số (Lemma 1.2: $(I+FN)^{-1} = I - F(I+NF)^{-1}N$), nguyên lý bị chặn đều Banach-Steinhaus, bất đẳng thức Young rời rạc cho tích chập dãy số (1.35), và phân tích giá trị kỳ dị (SVD).

Chapter 2: Singular systems of first-order difference equations

Chương này tập trung giải quyết toàn diện bài toán ổn định, ổn định dưới tác động của nhiễu, các định lý dạng Bohl-Perron và số mũ Bohl cho hệ sai phân suy biến cấp một.

Nội dung chi tiết:

  • Các khái niệm ổn định (Mục 2.1): Định nghĩa hình thức cho nghiệm tầm thường $x = 0$ của hệ thuần nhất: ổn định, ổn định đều, ổn định tiệm cận và ổn định mũ. Định lý 2.1 và 2.2 chứng minh rằng khi dãy toán tử chiếu chuẩn tắc ${\widetilde{P}_m}$ bị chặn:

    • Hệ ổn định đều khi và chỉ khi $|\Phi(n, m)| \le C, \forall n \ge m \ge n_0$.
    • Hệ ổn định mũ khi và chỉ khi $|\Phi(n, m)| \le K \omega^{n-m}$ với $K > 0$ và $0 < \omega < 1$. Tác giả xây dựng Ví dụ 2.1 cụ thể trong $\mathbb{R}^3$ có chứa tham số $\alpha > 0$ để minh họa: hệ ổn định mũ khi $\alpha < 1$, ổn định đều khi $\alpha = 1$, và không ổn định khi $\alpha > 1$.

  • Tính ổn định của phương trình chịu nhiễu (Mục 2.2):

    • Nhiễu một phía: Xét hệ $E_n y(n+1) = (A_n + B_n)y(n)$. Luận án đề xuất 4 giả thiết đối với ma trận nhiễu $B_n$:
      • Giả thiết 1: $\sup_{n \in \mathbb{N}(n_0)} |T_n \widetilde{Q}_n \widetilde{G}_n^{-1} B_n| < 1$.
      • Giả thiết 2: $\sum_{n=n_0}^\infty |\widetilde{P}_n \widetilde{G}_n^{-1} B_n| < \infty$.
      • Giả thiết 3: Tồn tại $\delta > 0, N > n_0$ sao cho $|\widetilde{P}_n \widetilde{G}_n^{-1} B_n| \le \delta, \forall n \ge N$.
      • Giả thiết 4: $|\widetilde{P}_n \widetilde{G}_n^{-1} B_n| \to 0$ khi $n \to \infty$. Định lý 2.3 chứng minh hệ chịu nhiễu bảo toàn tính ổn định đều dưới Giả thiết 1 và 2. Định lý 2.4 và Hệ quả 2.1 chứng minh hệ bảo toàn tính ổn định mũ dưới Giả thiết 1 và Giả thiết 3 (hoặc 4).
    • Nhiễu hai phía: Khảo sát hệ $(E_n + F_n)y(n+1) = (A_n + B_n)y(n)$. Luận án chỉ ra một ví dụ phản chứng (2.8): một nhiễu $F_n$ cực nhỏ làm thay đổi hạng có thể khiến nghiệm bùng nổ ra vô hạn ($y_2(n) = (n!)^2 y_2(0) \to \infty$). Từ đó, tác giả đưa ra Giả thiết 5 yêu cầu $F_n$ phải thuộc lớp nhiễu có cấu trúc cho phép thỏa mãn $\sup_{n} |F_n G_n^{-1}| < 1$, đảm bảo ma trận suy biến giữ nguyên chỉ số 1.

  • Định lý ổn định kiểu Bohl-Perron (Mục 2.3): Nghiên cứu mối quan hệ giữa đầu vào $q_n$ và nghiệm $y(n)$. Chứng minh 3 định lý cơ bản: nếu với mọi dãy vế phải $q \in l_p$ mà nghiệm $y$ tương ứng thuộc $l_q$ ($1 \le p \le q \le \infty$), thì hệ thuần nhất tương ứng phải đạt tính ổn định mũ hoặc ổn định đều.

  • Số mũ Bohl và tính ổn định mũ (Mục 2.4): Định nghĩa số mũ Bohl trên hệ sai phân suy biến thông qua toán tử Cauchy $\Phi(n, m)$. Chứng minh mối liên hệ chặt chẽ giữa dấu âm của số mũ Bohl cực đại với tính ổn định mũ của hệ. Thiết lập định lý về tính vững (robustness) của số mũ Bohl khi ma trận hệ số chịu các nhiễu cho phép nhỏ.

  • Trường hợp hàm toán tử chiếu chuẩn tắc không bị chặn (Mục 2.5): Khảo sát trường hợp suy biến phức tạp khi $|\widetilde{P}_n| \to \infty$. Phân tích sự phá vỡ các điều kiện cần trong Định lý 2.1, 2.2 và đề xuất các không gian nghiệm có trọng số để bảo toàn đặc trưng tiệm cận.

Chapter 3: Singular systems of second-order difference equations

Chương này trình bày các phát triển hoàn toàn mới dành riêng cho hệ phương trình sai phân suy biến tuyến tính cấp hai có hệ số biến thiên: $$A_n x(n+2) + B_n x(n+1) + C_n x(n) = f_n, \quad n \ge n_0$$ với $A_n$ suy biến ($\text{rank} A_n \equiv d_1 < d$).

Nội dung chi tiết:

  • Bài toán giá trị ban đầu và tính tương thích (Mục 3.1): Dưới giả thiết hệ không có chỉ số lạ ($\mu = 0$), tác giả sử dụng các phép biến đổi ma trận unita từng bước để đưa bộ ma trận $(A_n, B_n, C_n)$ về dạng tam giác khối chuẩn tắc: $$d_1 + d_2 + d_3 = d, \quad d_1 = \text{rank}(A_n)$$ Từ đó, xây dựng tường minh không gian các điều kiện ban đầu tương thích ${x(n_0), x(n_0+1)}$ để bài toán Cauchy có nghiệm duy nhất, đồng thời thiết lập mối liên hệ đẳng cấu giữa tập nghiệm của hệ suy biến cấp hai ban đầu với tập nghiệm của một hệ phương trình sai phân chính quy cấp hai thu gọn.

  • Khái niệm và tiêu chuẩn ổn định mũ (Mục 3.2 & 3.3): Đưa ra định nghĩa ổn định mũ cho hệ sai phân cấp hai thuần nhất. Bằng cách áp dụng nguyên lý so sánh (chuyển đổi hệ cấp hai sang hệ phương trình sai phân vector cấp một tương đương kết hợp với chặn ma trận chuẩn), tác giả chứng minh các tiêu chuẩn tường minh cho tính ổn định mũ thông qua bán kính phổ của ma trận suy rộng kết hợp.

  • Định lý Bohl-Perron cho hệ cấp hai (Mục 3.4): Mở rộng định lý Bohl-Perron: chứng minh rằng tính chất $l_p$-đầu vào dẫn đến $l_\infty$-trạng thái của hệ không thuần nhất bảo toàn tính ổn định mũ của hệ thuần nhất cấp hai.

  • Tính ổn định vững dưới tác động của nhiễu có cấu trúc (Mục 3.5): Khảo sát hệ cấp hai bị nhiễu: $$(A_n + D_{1,n} \Delta_{1,n} E_n)x(n+2) + (B_n + D_{2,n} \Delta_{2,n} E_n)x(n+1) + (C_n + D_{3,n} \Delta_{3,n} E_n)x(n) = 0$$ Xác định cấu trúc ma trận nhiễu cho phép (allowable perturbations) để bảo toàn chỉ số lạ $\mu = 0$ và kích thước các khối $d_1, d_2, d_3$. Bằng cách sử dụng nguyên lý so sánh, luận án thiết lập công thức tính bán kính ổn định vững và đưa ra cận trên tường minh cho chuẩn của khối nhiễu $\Delta = [\Delta_1, \Delta_2, \Delta_3]^T$ đảm bảo hệ bị nhiễu vẫn duy trì tính ổn định mũ.

Bảng tổng hợp đối chiếu phương pháp và kết quả giữa hai lớp phương trình

Tiêu chí phân tích Phương trình sai phân suy biến cấp một (Chương 2) Phương trình sai phân suy biến cấp hai (Chương 3)
Dạng phương trình $E_n y(n+1) = A_n y(n) + q_n$ $A_n x(n+2) + B_n x(n+1) + C_n x(n) = f_n$
Công cụ giải tích chính Phép chiếu chuẩn tắc, chỉ số kéo theo 1 (Tractability index-1) Chỉ số lạ bằng 0 (Strangeness-free form, $\mu = 0$)
Điều kiện ma trận dẫn đầu $\text{rank} E_n = r < d$, ma trận $G_n$ khả nghịch $\text{rank} A_n \equiv d_1 < d$, phân rã khối $d_1 + d_2 + d_3 = d$
Phương pháp khử liên kết Tách trực giao qua cặp phép chiếu $\widetilde{P}n, \widetilde{Q}{n-1}$ Phép biến đổi tương đương toàn cục về dạng khối tam giác
Biểu diễn nghiệm Toán tử tiến hóa Cauchy $\Phi(n, m)$ và biến thiên hằng số Chuyển đổi về hệ sai phân chính quy thu gọn tương đương
Công cụ đánh giá ổn định Đánh giá chuẩn toán tử Cauchy, bất đẳng thức Gronwall rời rạc Nguyên lý so sánh ma trận (Comparison principle)
Công cụ đặc trưng hóa nghiệm Số mũ Bohl $\kappa_B(A)$, định lý dạng Bohl-Perron trên $l_p$ Định lý dạng Bohl-Perron cấp hai, bán kính phổ ma trận khối
Phạm vi khảo sát nhiễu Nhiễu một phía $B_n$ và nhiễu hai phía $(F_n, B_n)$ Nhiễu có cấu trúc giới hạn $D_i \Delta_i E$ ($i=1, 2, 3$)

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

Căn cứ theo phần tự công bố của luận án, các đóng góp mới về mặt lý luận và khoa học gồm:

  1. Về phương trình sai phân suy biến cấp một:

    • Xây dựng thành công phương pháp phép chiếu chuẩn tắc duy nhất, loại bỏ sự phụ thuộc vào việc lựa chọn phép chiếu tùy ý trong các nghiên cứu trước đây.
    • Thiết lập các điều kiện đủ tường minh (Giả thiết 1 đến 5) để nghiệm của hệ sai phân suy biến bảo toàn tính ổn định đều và ổn định mũ trước tác động của các lớp nhiễu một phía và nhiễu hai phía có cấu trúc.
    • Mở rộng thành công 3 định lý kiểu Bohl-Perron từ phương trình sai phân chính quy sang hệ sai phân suy biến biến thiên theo thời gian.
    • Lần đầu tiên đưa ra định nghĩa số mũ Bohl cho hệ sai phân suy biến rời rạc, chứng minh mối liên hệ tương đương giữa số mũ Bohl âm với tính ổn định mũ, và chứng minh tính bền vững của số mũ Bohl trước nhiễu ma trận.

  2. Về phương trình sai phân suy biến cấp hai:

    • Xây dựng quy trình đại số tường minh để xác định tập điều kiện ban đầu tương thích ${x_0, x_1}$ cho bài toán Cauchy dưới giả thiết không có chỉ số lạ.
    • Thiết lập mối quan hệ tương đương một - một giữa tập nghiệm của hệ suy biến cấp hai và tập nghiệm của hệ sai phân chính quy thu gọn.
    • Đề xuất kỹ thuật nguyên lý so sánh mới, cho phép thiết lập các tiêu chuẩn ổn định mũ không phụ thuộc vào phổ ma trận thời gian thực.
    • Thiết lập công thức ước lượng cận trên tường minh cho ma trận nhiễu có cấu trúc, cung cấp tiêu chuẩn định lượng để hệ duy trì tính giải được và tính ổn định mũ (bán kính ổn định vững).

  3. Kiến nghị và giải pháp đề xuất từ luận án:

    • Đề xuất áp dụng kỹ thuật khử liên kết chuẩn tắc và nguyên lý so sánh vào việc tính toán số trị ổn định trong các bài toán điều khiển tự động rời rạc có ràng buộc đại số.
    • Khuyến nghị sử dụng cận nhiễu đã tính toán để thiết kế bộ điều khiển bền vững cho các mô hình kinh tế động học Leontief khi ma trận vốn $B$ có các phần tử chịu sai số ngẫu nhiên hoặc biến động tham số.

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

Hạn chế được ghi nhận từ văn bản

  • Luận án chỉ giới hạn phạm vi khảo sát ở lớp hệ phương trình sai phân tuyến tính; chưa xét đến các tương tác phi tuyến bậc cao.
  • Giả thiết về chỉ số chỉ dừng lại ở mức cơ bản: chỉ số kéo theo bằng 1 đối với hệ cấp một và không có chỉ số lạ ($\mu = 0$) đối với hệ cấp hai; chưa giải quyết trường hợp hệ có chỉ số cao ($\mu \ge 1$) vốn xuất hiện các đạo hàm sai phân bậc cao của số hạng không thuần nhất.
  • Đối với phương trình cấp cao, luận án mới chỉ dừng lại ở hệ cấp hai, chưa mở rộng trực tiếp công thức giải tích cho phương trình suy biến cấp $k$ tổng quát ($k > 2$).
  • Khảo sát nhiễu hai phía ở hệ cấp một và nhiễu có cấu trúc ở hệ cấp hai phụ thuộc vào các điều kiện giả thiết hạn chế cấu trúc (allowable structured perturbations), chưa bao quát toàn bộ các dạng nhiễu tự do bất kỳ.

Hướng nghiên cứu tiếp

Văn bản luận án đề xuất các hướng phát triển tiếp theo:

  1. Mở rộng phương pháp tiếp cận chỉ số lạ và nguyên lý so sánh cho hệ phương trình sai phân suy biến bậc cao tổng quát ($k \ge 3$).
  2. Nghiên cứu tính ổn định và tính giải được cho hệ sai phân suy biến có chỉ số cao ($\mu \ge 1$).
  3. Khảo sát tính ổn định của hệ sai phân suy biến chịu tác động của trễ (delay singular difference equations) hoặc hệ chuyển mạch suy biến (switched singular systems).
  4. Nghiên cứu hệ sai phân suy biến ngẫu nhiên hoặc chịu tác động của tham số ngẫu nhiên.

Giá trị tham khảo

Luận án có giá trị khoa học và ứng dụng thực tiễn rõ ràng đối với các nhóm đối tượng:

  • 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à tích phân, Toán ứng dụng: Có thể sử dụng Chương 1 và Chương 2 như một tài liệu tham khảo chuẩn mực về phương pháp toán tử chiếu, kỹ thuật phân rã nghiệm, và lý thuyết Bohl-Perron trên thang thời gian rời rạc.
  • Giảng viên và nhà nghiên cứu lý thuyết điều khiển hệ thống: Tham khảo trực tiếp Chương 3 về các tiêu chuẩn ổn định vững và công thức tính cận trên của nhiễu ma trận cho các hệ thống mô tả rời rạc (descriptor systems).
  • Các chuyên gia mô hình hóa kinh tế vĩ mô và kinh tế lượng: Ứng dụng trực tiếp kết quả giải trừ liên kết và phân tích độ nhạy của hệ suy biến vào mô hình Leontief động đa ngành, nhằm đánh giá độ ổn định của sản lượng khi ma trận vốn đầu tư chịu biến động thực tế.

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

1. Tại sao nghiên cứu phương trình sai phân suy biến (SDEs) lại phức tạp hơn phương trình sai phân thông thường?

Trong phương trình sai phân thông thường $x(n+1) = f(x(n))$, trạng thái tiếp theo được xác định tường minh duy nhất từ trạng thái trước đó. Trong phương trình sai phân suy biến $E_n y(n+1) = A_n y(n) + q_n$, ma trận $E_n$ bị suy biến ($\ker E_n \neq {0}$). Điều này dẫn đến hai khó khăn căn bản:

  • Bài toán giá trị ban đầu không thể giải được với mọi điều kiện ban đầu tùy ý mà chỉ giải được với các điều kiện ban đầu tương thích nằm trong một không gian con xác định.
  • Nghiệm của hệ là sự đan xen giữa thành phần động lực học tiến hóa theo thời gian và thành phần đại số tức thời; nếu số hạng không thuần nhất hoặc ma trận hệ số bị nhiễu không đúng cấu trúc, hệ có thể rơi vào tình trạng vô nghiệm, có vô số nghiệm, hoặc nghiệm bùng nổ mất ổn định.

2. Ý nghĩa của cặp phép chiếu chuẩn tắc (canonical projectors) trong Chương 2 là gì?

Khi giải trừ liên kết (decoupling) hệ cấp một bằng các phép chiếu thông thường, các ma trận chiếu $P_n$ phụ thuộc vào việc lựa chọn cơ sở đại số của hạt nhân $\ker E_n$. Luận án xây dựng cặp phép chiếu chuẩn tắc: $$\widetilde{Q}_{n-1} = -T_n Q_n G_n^{-1} A_n, \quad \widetilde{P}_n = I - \widetilde{Q}n$$ Cặp phép chiếu này là duy nhất (không phụ thuộc vào việc chọn $T_n$ và $Q_n$), giúp phân tách chính xác và bất biến không gian trạng thái thành không gian con động học $S_n$ và không gian ràng buộc đại số $N{n-1}$, từ đó xác định duy nhất toán tử Cauchy $\Phi(n, m)$ cho hệ.

3. Cần những điều kiện gì để hệ suy biến cấp một bảo toàn tính ổn định mũ dưới tác động của nhiễu một phía $B_n$?

Theo Định lý 2.4 và Hệ quả 2.1, hệ suy biến bảo toàn tính ổn định mũ khi:

  • Hệ không bị nhiễu có chỉ số kéo theo bằng 1, ổn định mũ và dãy phép chiếu chuẩn tắc ${\widetilde{P}_n}$ bị chặn.
  • Ma trận nhiễu $B_n$ thỏa mãn Giả thiết 1: $\sup_{n} |T_n \widetilde{Q}_n \widetilde{G}_n^{-1} B_n| < 1$ (để đảm bảo tính khả nghịch của thành phần đại số và duy trì chỉ số 1).
  • Ma trận nhiễu $B_n$ thỏa mãn Giả thiết 3: Tồn tại $\delta > 0$ đủ nhỏ và chỉ số $N$ sao cho $|\widetilde{P}_n \widetilde{G}_n^{-1} B_n| \le \delta$ với mọi $n \ge N$ (hoặc Giả thiết 4: $|\widetilde{P}_n \widetilde{G}_n^{-1} B_n| \to 0$ khi $n \to \infty$).

4. Tại sao một nhiễu hai phía $F_n$ tùy ý dù nhỏ cũng có thể phá vỡ hoàn toàn tính ổn định của hệ?

Ví dụ (2.8) trong văn bản chỉ ra rằng nếu ma trận nhiễu $F_n$ xuất hiện ở vế trái $(E_n + F_n)y(n+1)$ làm thay đổi hạng của $E_n$ hoặc phá vỡ cấu trúc không gian triệt tiêu (ker), nó có thể biến các ràng buộc đại số tĩnh thành các phương trình động lực có nghiệm gia tăng nhanh. Cụ thể trong ví dụ, thành phần nghiệm đại số bị biến thành $y_2(n) = (n!)^2 y_2(0)$, tiến ra vô hạn khi $n \to \infty$. Do đó, để khảo sát tính ổn định, bắt buộc phải đặt ra điều kiện nhiễu có cấu trúc cho phép (allowable structured perturbation) như Giả thiết 5.

5. Dạng chuẩn tắc không có chỉ số lạ (strangeness-free) giúp ích gì cho việc giải hệ cấp hai ở Chương 3?

Dưới giả thiết chỉ số lạ $\mu = 0$, bộ ba ma trận $(A_n, B_n, C_n)$ có thể biến đổi tương đương toàn cục về dạng khối tam giác có cấu trúc xác định theo các kích thước $d_1 + d_2 + d_3 = d$. Cấu trúc này cho phép:

  • Tách rời các phương trình thuần túy đại số để kiểm tra điều kiện tương thích của dữ kiện ban đầu $x_0, x_1$.
  • Triệt tiêu các ẩn số phụ và đưa hệ về một phương trình sai phân chính quy cấp hai trên không gian con số chiều $d_1$, từ đó áp dụng trực tiếp nguyên lý so sánh để đánh giá nghiệm mà không làm biến dạng động học của hệ gốc.

Kết luận

Luận án tiến sĩ của tác giả Ngô Thị Thanh Nga đã giải quyết một cách hệ thống bài toán ổn định và ổn định vững cho hai lớp phương trình sai phân suy biến tuyến tính cấp một và cấp hai có hệ số biến thiên theo thời gian. Thông qua việc kết hợp chặt chẽ giữa phương pháp phép chiếu chuẩn tắc và lý thuyết chỉ số lạ, tác giả đã thiết lập thành công các tiêu chuẩn ổn định đều, ổn định mũ, các định lý đặc trưng hóa nghiệm dạng Bohl-Perron, cũng như lý thuyết số mũ Bohl cho hệ rời rạc suy biến. Các kết quả đạt được không chỉ lấp đầy các khoảng trống lý thuyết quan trọng trong lĩnh vực phương trình vi phân và tích phân mà còn cung cấp công cụ toán học định lượng tin cậy cho việc phân tích và điều khiển các hệ thống động lực học rời rạc có ràng buộc trong thực tế.