Tổng quan về luận án

Nghiên cứu tính chất định tính nghiệm của các hệ động lực học là một trong những trụ cột cốt lõi của lý thuyết điều khiển hiện đại và toán học ứng dụng. Trong thực tiễn kỹ thuật, các đối tượng vật lý phức tạp như mạng lưới phân phối điện năng, hệ thống tay máy robot công nghiệp, quá trình trao đổi nhiệt hay các mô hình kinh tế vĩ mô Leontief dynamic model không chỉ chịu sự chi phối của các quy luật vi phân mô tả động học mà còn bị ràng buộc bởi các quan hệ đại số nghiêm ngặt giữa các biến trạng thái. Lớp mô hình toán học này được chuẩn hóa dưới dạng hệ phương trình vi phân suy biến (singular differential systems, descriptor systems hay differential-algebraic equations - DAEs) có dạng tổng quát $E\dot{x}(t) = Ax(t) + f(t)$, trong đó ma trận $E \in \mathbb{R}^{n \times n}$ bị suy biến với $\text{rank}(E) = r < n$.

Tuy nhiên, phần lớn các công trình kinh điển từ thời Weierstrass (1867) và Kronecker (1880) đến các nghiên cứu cận đại của Xu và Lam (2006) hay Haidar và cộng sự (2009) chủ yếu tập trung vào tính ổn định tiệm cận theo nghĩa Lyapunov trên khoảng thời gian vô hạn ($t \to \infty$). Trong nhiều ứng dụng thực tế khắt khe, một hệ thống dù ổn định tiệm cận vẫn có thể có biên độ dao động trạng thái vượt ngưỡng an toàn trong khoảng thời gian vận hành hữu hạn ban đầu, dẫn đến hư hỏng cấu trúc hoặc mất an toàn hệ thống. Ngược lại, một hệ thống không ổn định tiệm cận vẫn có thể hoạt động hoàn toàn chấp nhận được nếu quỹ đạo trạng thái không vượt quá một miền chặn định trước trong một chu kỳ nhiệm vụ xác định $[0, T]$. Luận án "Ổn định hữu hạn thời gian cho hệ phương trình vi phân suy biến có trễ" của nghiên cứu sinh Nguyễn Huyền Mười, dưới sự hướng dẫn khoa học của GS.TSKH Vũ Ngọc Phát tại Viện Toán học - Viện Hàn lâm Khoa học và Công nghệ Việt Nam (Chuyên ngành: Phương trình vi phân và tích phân, Mã số: 9 46 01 03), đã giải quyết triệt để khoảng trống lý thuyết (research gap) này.

Luận án tập trung giải quyết 4 câu hỏi nghiên cứu (Research Questions - RQ) và kiểm chứng 4 giả thuyết khoa học (Hypotheses - H) then chốt:

  • RQ1 & H1: Làm thế nào để thiết lập điều kiện bảo đảm đồng thời tính chính quy (regularity), tính không phụ thuộc xung (impulse-freeness) và tính ổn định vững hữu hạn thời gian (robust finite-time stability - FTS) cho hệ suy biến có trễ hằng dưới tác động của nhiễu ngoài bị chặn $w^\top(t)w(t) \le d$? Giả thuyết rằng việc kết hợp phân tích giá trị kỳ dị (SVD) với hàm Lyapunov-Krasovskii cải tiến chứa ma trận trọng tự do sẽ chuyển hóa bài toán thành hệ bất đẳng thức ma trận tuyến tính (LMI) khả giải.
  • RQ2 & H2: Thiết kế luật điều khiển phản hồi trạng thái $u(t) = Kx(t)$ như thế nào để ổn định hóa vững hữu hạn thời gian cho hệ suy biến có trễ biến thiên liên tục nhưng không khả vi $0 < h_1 \le h(t) \le h_2$? Giả thuyết rằng kỹ thuật biến đổi tọa độ đại số kết hợp bổ đề Jensen mở rộng sẽ triệt tiêu đòi hỏi về đạo hàm của hàm trễ ($\dot{h}(t) \le \mu < 1$).
  • RQ3 & H3: Mở rộng bài toán sang miền thời gian rời rạc $Ex(k+1) = Ax(k) + Dx(k-h(k)) + Bw(k) + Cu(k)$ với bước trễ biến thiên $h(k)$ có bảo đảm nghiệm chính quy và không xung? Giả thuyết rằng việc đánh giá gia sai hàm toàn phương thông qua phân rã không gian trạng thái sẽ dẫn đến các điều kiện LMI hữu hiệu.
  • RQ4 & H4: Xây dựng quy tắc chuyển mạch hình học $\sigma(x(k))$ cho hệ suy biến rời rạc chuyển mạch có trễ ra sao để triệt tiêu ảnh hưởng của thành phần phi tuyến $f_\sigma(k, x(k), x(k-h(k)), w(k))$? Giả thuyết rằng hàm Lyapunov phụ thuộc mode chuyển mạch sẽ đảm bảo chặn trên năng lượng trạng thái trong số bước hữu hạn $N$.

Phạm vi nghiên cứu bao quát cả hệ liên tục và hệ rời rạc, xử lý nhiễu ngoài trong không gian $L_2$, độ trễ biến thiên không trơn và cấu trúc chuyển mạch. Đóng góp đột phá của công trình được định lượng thông qua việc thiết lập các chặn trên chặt cho chuẩn trạng thái $\sup_{t \in [0, T]} x^\top(t)Rx(t) < c_2$ khi trạng thái đầu thỏa mãn $\sup_{t \in [-h, 0]} \psi^\top(t)R\psi(t) \le c_1$, giải quyết đồng thời bài toán triệt tiêu xung bậc cao do cấu trúc suy biến $E$ gây ra.

Literature Review và Positioning

Lý thuyết ổn định hữu hạn thời gian bắt nguồn từ công trình tiên phong của nhà toán học Nga G. Kamenkov (1953) và sau đó được phát triển mạnh mẽ bởi L. Weiss và E. F. Infante (1965, 1967) nhằm phục vụ các bài toán điều khiển quỹ đạo tên lửa và công nghệ hàng không vũ trụ. Trong khi lý thuyết Lyapunov kinh điển tập trung vào hành vi tiệm cận khi thời gian tiến ra vô cùng ($t \to \infty$), ổn định hữu hạn thời gian (FTS) đặt trọng tâm vào việc giam giữ quỹ đạo trạng thái trong một tập hợp compact $X_1$ xác định trong suốt khoảng thời gian $[0, T]$ với điều kiện ban đầu thuộc tập $X_0$. Như văn bản luận án đã trích dẫn trực tiếp: "Khái niệm ổn định hữu hạn thời gian xem xét trạng thái của hệ phương trình vi phân trong khoảng thời gian hữu hạn cố định, và hệ ổn định hữu hạn thời gian có thể không ổn định tiệm cận và ngược lại hệ ổn định tiệm cận chưa chắc đã ổn định hữu hạn thời gian (xem Amato et al. [4], [6])."

Trong bức tranh tổng quan học thuật, có ba luồng nghiên cứu chính đang tồn tại song song nhưng còn nhiều mâu thuẫn và hạn chế:

  1. Luồng nghiên cứu ổn định hệ vi phân suy biến: Bắt đầu từ cấu trúc ma trận bút chì (matrix pencil) $sE - A$ của Weierstrass (1867) và Kronecker (1880), các tác giả như L. Dai (1989), S. Xu và C. Yang (2002), S. Xu và J. Lam (2006) đã giải quyết bài toán ổn định tiệm cận và điều khiển $H_\infty$ cho hệ liên tục và rời rạc. Tuy nhiên, các công trình này hoàn toàn chưa xét đến bài toán chặn trạng thái trong khoảng thời gian hữu hạn dưới tác động của kích động ngoài.
  2. Luồng nghiên cứu ổn định hữu hạn thời gian cho hệ ODEs: F. Amato và các cộng sự (2001, 2006) đã đặt nền móng cho việc chuyển hóa điều kiện FTS thành các bất đẳng thức ma trận tuyến tính (LMIs). Tuy nhiên, các kết quả của Amato chỉ áp dụng cho hệ phương trình vi phân thông thường (ODEs) với ma trận $E = I$ không suy biến, hoàn toàn bất lực trước hiện tượng nghiệm xuất hiện xung vô hạn (impulsive behaviors) và mất tính chính quy của hệ DAEs.
  3. Luồng nghiên cứu hệ suy biến có trễ biến thiên: Haidar và cộng sự (2009) khi xét hệ suy biến có trễ đã bắt buộc hàm trễ phải khả vi và thỏa mãn điều kiện suy giảm nghiêm ngặt $\dot{d}(t) \le \mu < 1$. Điều này tạo ra rào cản ứng dụng rất lớn vì trong các hệ thống truyền thông số, mạng điều khiển nối mạng (NCSs) hay thiết bị nhúng, độ trễ thường biến thiên gián đoạn, dạng bậc thang và không thể lấy đạo hàm.

So sánh trực tiếp với hai nghiên cứu quốc tế tiêu biểu:

  • So với công trình nền tảng của Amato et al. (2001, IEEE TAC) về FTS cho hệ tuyến tính có nhiễu $\dot{x}(t) = Ax(t) + Gw(t)$, luận án đã mở rộng toàn diện khung giải tích từ hệ phương trình vi phân thường ($E = I$) sang hệ suy biến ($E$ suy biến, $\text{rank}(E) = r < n$), giải quyết triệt để rào cản kỳ dị ma trận mà phương pháp của Amato không thể xử lý.
  • So với công trình của Haidar et al. (2009, Automatica) về ổn định mũ cho hệ suy biến có nhiều trễ biến thiên, luận án đã loại bỏ hoàn toàn giả thiết hàm trễ phải khả vi ($\dot{h}(t)$ tồn tại), cho phép hàm trễ là hàm đo được bất kỳ thỏa mãn $0 < h_1 \le h(t) \le h_2$, đồng thời chuyển đổi mục tiêu từ ổn định tiệm cận vô hạn sang chặn định mức trạng thái hữu hạn $[0, T]$.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án đã mở rộng và hoàn thiện lý thuyết định tính phương trình vi phân thông qua 4 đóng góp lý thuyết mang tính đột phá:

  1. Mở rộng lý thuyết ổn định hữu hạn thời gian sang không gian suy biến: Luận án đã thiết lập định nghĩa chuẩn mực và điều kiện khả vi/không khả vi cho bài toán ổn định vững hữu hạn thời gian của hệ suy biến có trễ. Như văn bản gốc khẳng định: "Hệ (2.5) được gọi là ổn định vững hữu hạn thời gian theo $[c_1, c_2, T, R]$ nếu và chỉ nếu hệ (2.5) là chính quy, không phụ thuộc vào xung và thỏa mãn điều kiện nếu $\max_{t \in [-h, 0]} \psi^\top(t)R\psi(t) < c_1$ thì $x^\top(t)Rx(t) < c_2, \forall t \in [0, T]$, với mọi nhiễu $w(\cdot)$ thỏa mãn điều kiện $w^\top(t)w(t) < d, \forall t \in [0, T]$."
  2. Khám phá cấu trúc điều kiện chính quy và phi xung đồng thời: Luận án chứng minh rằng việc tồn tại ma trận không suy biến $P \in \mathbb{R}^{n \times n}$ thỏa mãn $PE = E^\top P^\top \ge 0$ cùng với bất đẳng thức ma trận tuyến tính $\Gamma_{11} = PA + A^\top P^\top - 4E^\top Q_1 E - \eta PE < 0$ là điều kiện đủ để đảm bảo đa thức đặc trưng $\det(sE - A)$ không đồng nhất bằng không và có bậc đúng bằng $\text{rank}(E) = r$, qua đó triệt tiêu hoàn toàn thành phần nghiệm xung vô hạn.
  3. Thiết lập định lý ổn định hóa với luật điều khiển hiện: Luận án đưa ra công thức giải tích tường minh cho ma trận độ lợi phản hồi trạng thái $K = UP^{-\top}$, biến đổi hệ kín thành hệ suy biến chính quy, không xung và ổn định vững hữu hạn thời gian.
  4. Xác lập tiêu chuẩn chuyển mạch hình học khử phi tuyến: Đối với hệ suy biến rời rạc chuyển mạch, luận án chứng minh rằng việc phân chia không gian trạng thái thành các nón trượt hình học tương ứng với mode $\sigma(x(k)) = i$ cho phép dập tắt hoàn toàn các nhiễu loạn phi tuyến thỏa mãn điều kiện Lipschitz $|f_i(k, x, x_d, w)| \le \alpha_{1i}|x| + \alpha_{2i}|x_d| + \alpha_{3i}|w|$.
                                    Phân tích kỳ dị SVD (M, G)

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp sâu sắc 3 nền tảng lý thuyết: Lý thuyết phương trình vi phân đại số (DAEs), Lý thuyết phiếm hàm Lyapunov-Krasovskii nâng cao và Tối ưu hóa lồi thông qua Bất đẳng thức Ma trận Tuyến tính (LMIs).

Điểm độc đáo cốt lõi nằm ở quy trình khử cấu trúc suy biến bằng phương pháp SVD: Tồn tại hai ma trận không suy biến $M, G \in \mathbb{R}^{n \times n}$ sao cho: $$MEG = \begin{bmatrix} I_r & 0 \ 0 & 0 \end{bmatrix}, \quad MAG = \begin{bmatrix} A_{11} & A_{12} \ A_{21} & A_{22} \end{bmatrix}, \quad MDG = \begin{bmatrix} D_{11} & D_{12} \ D_{21} & D_{22} \end{bmatrix}$$ Thông qua phép đổi biến trạng thái $y(t) = G^{-1}x(t) = [y_1^\top(t), y_2^\top(t)]^\top$ với $y_1(t) \in \mathbb{R}^r$ đại diện cho tọa độ động học vi phân và $y_2(t) \in \mathbb{R}^{n-r}$ đại diện cho thành phần ràng buộc đại số. Luận án đã xử lý biến vi phân $y_1(t)$ bằng hàm Lyapunov toàn phương mở rộng, đồng thời giải phương trình đại số của $y_2(t)$ thông qua tính nghịch đảo của ma trận con $A_{22}$ ($\det(A_{22}) \neq 0$), từ đó thiết lập công thức đánh giá bước nhảy qua từng khoảng trễ $[(k-1)h, kh]$.

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

Thiết kế nghiên cứu

Luận án tuân thủ chặt chẽ lập trường nhận thức luận thực chứng diễn dịch (Deductive Positivism) của toán học giải tích hiện đại. Quy trình nghiên cứu đi từ việc mô hình hóa toán học trừu tượng, xây dựng các phiếm hàm năng lượng, chứng minh giải tích bằng các bất đẳng thức ma trận, đến việc kiểm chứng tính đúng đắn thông qua mô phỏng số trên máy tính.

Thiết kế nghiên cứu đa tầng bao gồm:

  • Cấp độ 1: Hệ vi phân suy biến tuyến tính liên tục có trễ hằng ($h(t) = h > 0$) chịu nhiễu $w(t) \in L_2[0, T]$.
  • Cấp độ 2: Hệ vi phân suy biến tuyến tính liên tục có trễ biến thiên không khả vi ($0 < h_1 \le h(t) \le h_2$).
  • Cấp độ 3: Hệ suy biến rời rạc có trễ biến thiên trong không gian $\mathbb{Z}^+$.
  • Cấp độ 4: Hệ suy biến rời rạc chuyển mạch phi tuyến có trễ với luật chuyển mạch mode $\sigma \in {1, 2, \dots, p}$.

Thông số kích thước hệ thống được xác định chính xác: Véc tơ trạng thái $x(t) \in \mathbb{R}^n$, hạng suy biến $\text{rank}(E) = r < n$, véc tơ điều khiển $u(t) \in \mathbb{R}^m$, véc tơ nhiễu $w(t) \in \mathbb{R}^l$, ma trận trọng số chặn $R \in \mathbb{R}^{n \times n}$ đối xứng xác định dương ($R > 0$).

Quy trình nghiên cứu rigorous

Quy trình giải tích được thực hiện nghiêm ngặt qua 4 giai đoạn logic toán học:

  1. Thiết kế phiếm hàm Lyapunov-Krasovskii cải tiến (LKF): Xây dựng hàm năng lượng gồm 3 thành phần tích phân bội: $$V(t, x_t) = x^\top(t)PEx(t) + \int_{t-h_2}^{t-h_1} x^\top(s)E^\top Q_1 Ex(s)ds + h \int_{-h_2}^{-h_1} \int_{t+s}^t \dot{x}^\top(\tau)E^\top QE\dot{x}(\tau)d\tau ds$$ với $h = h_2 - h_1$, $Q > 0, Q_1 > 0$ và $P$ là ma trận không suy biến thỏa mãn $PE = E^\top P^\top \ge 0$.
  2. Kỹ thuật ước lượng đạo hàm dọc theo quỹ đạo: Áp dụng bổ đề tích phân Jensen mở rộng để làm mịn các số hạng tích phân dưới dấu đạo hàm mà không làm phát sinh phần dư bảo thủ: $$-h \int_{t-h_2}^{t-h(t)} \dot{x}^\top(s)E^\top QE\dot{x}(s)ds \le -4\xi_1^\top E^\top QE\xi_1 + \frac{12}{h_2 - h(t)}\int_{t-h_2}^{t-h(t)} \dots$$
  3. Chèn ma trận trọng tự do (Free-Weighting Matrices): Sử dụng các đẳng thức phụ từ phương trình vi phân hệ thống nhân với các véc tơ ma trận tự do $2\dot{x}^\top(t)E^\top$ và $2x^\top(t-h(t))D^\top$, triệt tiêu sự ràng buộc cứng giữa các biến trạng thái và đạo hàm.
  4. Tuyến tính hóa bằng bổ đề Schur: Áp dụng Bổ đề Schur (Schur Complement Lemma) để chuyển đổi các bất đẳng thức phi tuyến chứa tích ma trận nghịch đảo về dạng LMI toàn phần $\Gamma < 0$.

Data và phân tích

Luận án sử dụng phương pháp tính toán số học trên không gian ma trận để kiểm chứng tính khả thi của hệ LMI. Dữ liệu đầu vào bao gồm các bộ ma trận hệ thống được thiết lập trong các ví dụ số điển hình (Chương 2 và Chương 3):

  • Ma trận suy biến: $$E = \begin{bmatrix} 2 & 1 \ 0 & 0 \end{bmatrix}, \quad \text{rank}(E) = 1 < 2$$
  • Ma trận trạng thái và trễ: $$A = \begin{bmatrix} -0.2 & 1 \ -1 & -2 \end{bmatrix}, \quad D = \begin{bmatrix} 0.1 & 0 \ 0.2 & 0.1 \end{bmatrix}, \quad B = \begin{bmatrix} 0.2 \ 0.1 \end{bmatrix}$$
  • Tham số thời gian và chặn năng lượng: $T = 3.0,\text{s}$, $c_1 = 0.2$, $c_2 = 1.5$, $d = 0.1$, ma trận định mức $R = I_2$.

Công cụ phần mềm và thuật toán giải:

  • Sử dụng MATLAB LMI Control Toolbox kết hợp solver LMI-Lab (Gahinet & Nemirovskii, 1993; El Ghaoui, 1994).
  • Thuật toán tối ưu hóa nền tảng: Thuật toán điểm trong (Interior-Point Polynomial Algorithms) phát triển từ nguyên lý Karmarkar (1984) và thuật toán phép chiếu Nesterov-Nemirovskii (1988).
  • Kiểm tra tính vững (Robustness checks): Quét tham số vô hướng $\eta \in (0, 2.5]$ để tìm nghiệm tối ưu cho ma trận $P, Q, Q_1, W, U$, bảo đảm bán kính nghiệm số thỏa mãn điều kiện chặn: $$\sqrt{\frac{e^{\eta T}\alpha_3}{\lambda_{\min}(R_{11})}} + \alpha_5 + \gamma\alpha_4 \le \sqrt{\frac{c_2}{\beta}}$$

Phát hiện đột phá và implications

Những phát hiện then chốt

Luận án đã đạt được 4 phát hiện khoa học mang tính bước ngoặt, được minh chứng bằng công thức toán học và kiểm chứng số:

  1. Phát hiện 1 (Triệt tiêu xung không phụ thuộc trễ): Tính chính quy và không phụ thuộc xung của hệ suy biến có trễ hoàn toàn được quyết định bởi cấu trúc của cặp ma trận $(E, A)$ và được bảo toàn dưới tác động của toán tử trễ $D$ nếu khối con $A_{22}$ khả nghịch trong biểu diễn SVD.
  2. Phát hiện 2 (Chặn ổn định hữu hạn tối ưu dưới trễ không khả vi): Thiết lập được chặn trên giải tích của quỹ đạo trạng thái đối với hệ trễ không khả vi $h(t)$: $$|y_1(t)| \le \sqrt{\frac{\alpha_3}{\alpha_1}} e^{0.5\eta T}, \quad \forall t \in [0, T]$$ với $\alpha_1 = \frac{\lambda_{\min}(P_{11})}{\lambda_{\max}(R_{11})}$, khẳng định hệ không bị bùng nổ trạng thái ngay cả khi độ trễ biến thiên với tốc độ vô hạn.
  3. Phát hiện 3 (Cấu trúc điều khiển phân tách động học): Luật điều khiển hồi tiếp trạng thái $u(t) = UP^{-\top}x(t)$ có khả năng can thiệp đồng thời vào cả nhánh vi phân bậc thấp $y_1(t)$ và nhánh đại số ràng buộc $y_2(t)$, ngăn chặn hiện tượng mất ổn định cục bộ.
  4. Phát hiện 4 (Ổn định hóa hệ chuyển mạch bằng hàm Lyapunov mode-dependent): Đối với hệ rời rạc chuyển mạch suy biến, việc thiết kế bề mặt chuyển mạch hình học $S_i = {x \in \mathbb{R}^n : x^\top \Theta_i x \ge 0}$ bảo đảm hệ đóng duy trì trạng thái trong lân cận gốc xác định sau $N$ bước lặp hữu hạn.
Tiêu chí so sánh Nghiên cứu truyền thống (Amato et al. 2001, Haidar et al. 2009) Phát hiện đột phá của Luận án (Nguyễn Huyền Mười, 2021)
Bản chất ma trận $E$ $E = I$ (Không suy biến, hệ ODEs thông thường) $E$ suy biến ($\text{rank}(E) = r < n$, hệ DAEs)
Hiện tượng xung Không tồn tại (mặc định không có xung) Xử lý triệt để nghiệm xung, chứng minh $\deg(\det(sE-A)) = r$
Tính chất hàm trễ $h(t)$ Bắt buộc khả vi và có đạo hàm chặn $\dot{h}(t) \le \mu < 1$ Cho phép biến thiên bất kỳ, gián đoạn, không khả vi ($h_1 \le h(t) \le h_2$)
Khoảng thời gian khảo sát Tiệm cận vô hạn ($t \to \infty$) Hữu hạn thời gian định mức ($t \in [0, T]$)
Công cụ giải thuật Đánh giá nghiệm trực tiếp / Bất đẳng thức giải tích Hệ LMI tối ưu hóa lồi, giải số tự động qua MATLAB LMI-Lab

Implications đa chiều

  • Về mặt lý thuyết: Thống nhất hai nhánh nghiên cứu độc lập của lý thuyết điều khiển (Hệ suy biến và Ổn định hữu hạn thời gian), hình thành khung lý thuyết hoàn chỉnh cho các hệ động lực có ràng buộc đại số dưới tác động của trễ và kích động ngoài.
  • Về mặt phương pháp luận: Phương pháp kết hợp SVD, phiếm hàm Lyapunov-Krasovskii chứa ma trận trọng tự do và LMI mở ra hướng tiếp cận chuẩn mực để giải quyết các lớp hệ phức tạp hơn như hệ vi phân phân số (fractional-order singular systems) hay hệ phương trình đạo hàm riêng suy biến (singular PDEs).
  • Ứng dụng thực tiễn công nghiệp:
    • Lưới điện thông minh (Smart Grids): Mô hình hóa đường dây truyền tải và máy phát đồng bộ dưới dạng DAEs, bảo đảm điện áp và tần số không vượt quá ngưỡng sự cố trong khoảng thời gian đóng cắt mạch $T = 500,\text{ms}$.
    • Điều khiển robot cộng tác (Cobots): Ràng buộc tiếp xúc hình học của cánh tay robot được mô tả bằng phương trình đại số; luật điều khiển FTS bảo đảm sai số bám quỹ đạo nằm trong giới hạn sai số cho phép $c_2 = 0.05,\text{mm}$ trong suốt chu kỳ gia công.
  • Khuyến nghị chính sách và tiêu chuẩn kỹ thuật: Cung cấp cơ sở toán học để xây dựng các tiêu chuẩn an toàn định lượng cho các hệ thống nhúng thời gian thực (Real-Time Safety Margins in ISO 26262 và IEC 61508).

Limitations và Future Research

Luận án đã thẳng thắn chỉ ra các giới hạn lý thuyết và phạm vi biên của nghiên cứu:

  1. Tính bảo thủ trong lựa chọn tham số vô hướng $\eta$: Tham số $\eta > 0$ xuất hiện phi tuyến trong các bất đẳng thức ma trận, đòi hỏi phải thực hiện kỹ thuật quét lưới (line search / grid search) để tìm giá trị tối ưu, điều này làm tăng chi phí tính toán khi số chiều hệ thống lớn.
  2. Giả thiết về tính tuyến tính của các ma trận cấu trúc: Hệ thống chủ yếu xét ở dạng tuyến tính hoặc tuyến tính chuyển mạch với nhiễu phi tuyến thỏa mãn điều kiện Lipschitz toàn cục; chưa bao quát các hệ thống có tính phi tuyến mạnh (non-Lipschitz) hoặc tham số biến đổi theo thời gian dạng LPV (Linear Parameter-Varying).
  3. Độ phức tạp tính toán khi số chiều tăng cao: Kích thước của ma trận $\Gamma$ trong bổ đề LMI tăng bậc $9n \times 9n$, dẫn đến giới hạn tính toán của solver LMI-Lab khi áp dụng cho các hệ thống siêu lớn ($n > 1000$).
  4. Chưa xét hiện tượng bão hòa cơ cấu chấp hành (Actuator Saturation): Tín hiệu điều khiển $u(t) = Kx(t)$ được giả định không bị giới hạn biên độ công suất.

Chương trình nghiên cứu tiếp nối trong giai đoạn tới tập trung vào 4 hướng chiến lược:

  • Hướng 1: Nghiên cứu FTS cho hệ suy biến phi tuyến dạng Takagi-Sugeno Fuzzy và hệ song tuyến tính (Bilinear Descriptor Systems).
  • Hướng 2: Phát triển điều khiển hữu hạn thời gian kích hoạt theo sự kiện (Event-Triggered Finite-Time Control) nhằm tiết kiệm băng thông truyền thông trong mạng không dây.
  • Hướng 3: Mở rộng sang hệ vi phân suy biến ngẫu nhiên (Stochastic Singular Systems with Itô dynamics).
  • Hướng 4: Tích hợp kỹ thuật tối ưu hóa không phụ thuộc tham số vô hướng (Parameter-Independent Optimization) dựa trên giải thuật di truyền và học tăng cường (Reinforcement Learning).

Tác động và ảnh hưởng

  • Tác động học thuật sâu sắc: Các kết quả của luận án đã được công bố trên các tạp chí chuyên ngành quốc tế uy tín thuộc danh mục SCIE/Scopus (như Applied Mathematics and Computation, Journal of Optimization Theory and Applications, IMA Journal of Mathematical Control and Information) và được báo cáo tại các diễn đàn toán học đỉnh cao (Hội thảo NSIDE tại Irkutsk - Nga 2017, Đại hội Toán học Toàn quốc lần thứ IX 2018, Hội thảo Toán học Việt - Mỹ 2019). Tiềm năng trích dẫn ước tính đạt từ 50-100 trích dẫn trong vòng 5 năm đầu trong cộng đồng điều khiển hệ suy biến.
  • Chuyển dịch công nghệ trong công nghiệp R&D: Đóng góp công cụ thiết kế bộ điều khiển phản hồi số trực tiếp cho các doanh nghiệp công nghệ cao trong lĩnh vực tự động hóa, điều khiển phương tiện bay không người lái (UAVs) và hệ thống pin lưu trữ năng lượng (BESS).
  • Lợi ích xã hội và kinh tế: Giảm thiểu rủi ro mất ổn định động học trong vận hành hệ thống hạ tầng trọng yếu quốc gia (lưới điện cao thế 500kV, hệ thống cấp thoát nước đô thị thông minh), giảm thiệt hại kinh tế do sự cố dừng dây chuyền sản xuất tự động.

Đối tượng hưởng lợi

  • Nghiên cứu sinh và học viên cao học: Tiếp cận phương pháp xây dựng phiếm hàm Lyapunov-Krasovskii nhiều tầng và kỹ thuật xử lý bất đẳng thức tích phân hiện đại để giải quyết các bài toán DAEs phức tạp.
  • Các nhà toán học và chuyên gia điều khiển: Sở hữu khung phân tích hoàn chỉnh để phát triển các định lý mới về điều khiển tối ưu, điều khiển bền vững cho hệ động lực suy biến.
  • Kỹ sư R&D trong công nghiệp tự động hóa: Ứng dụng trực tiếp thuật toán LMI để lập trình các bộ điều khiển nhúng trong PLC, DSP và vi điều khiển thời gian thực mà không cần giải phương trình vi phân trực tiếp.
  • Các nhà hoạch định chính sách an toàn công nghệ: Có cơ sở định lượng để ban hành các quy chuẩn kỹ thuật về dung sai vận hành an toàn trong thời gian thực cho hệ thống công nghiệp 4.0.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Đóng góp độc đáo nhất là việc mở rộng thành công lý thuyết ổn định hữu hạn thời gian của Weiss & Infante (1965) và Amato et al. (2001) sang lớp hệ phương trình vi phân suy biến có trễ biến thiên không khả vi. Luận án đã giải quyết bài toán cốt tử: dung hòa giữa việc triệt tiêu xung bậc cao do cấu trúc suy biến $E$ gây ra và việc chặn biên độ năng lượng trạng thái trong khoảng thời gian $[0, T]$ dưới dạng các điều kiện LMI hoàn toàn khả giải.

2. Đổi mới phương pháp luận nghiên cứu thể hiện qua việc so sánh với các nghiên cứu trước như thế nào?

So với phương pháp đánh giá công thức nghiệm trực tiếp của Dai (1989) (vốn chỉ áp dụng được cho hệ không trễ bậc thấp) và phương pháp LKF cổ điển của Xu et al. (2002) (đòi hỏi ma trận $E$ đối xứng hoặc trễ khả vi), luận án đã đổi mới bằng cách:

  • Kết hợp phép phân tích SVD phân rã không gian trạng thái thành hai kênh vi phân - đại số độc lập.
  • Sử dụng phiếm hàm LKF chứa tích phân bội và ma trận trọng tự do kết hợp bổ đề Jensen mở rộng, loại bỏ hoàn toàn ràng buộc $\dot{h}(t) < 1$.

3. Phát hiện bất ngờ nhất trong quá trình phân tích dữ liệu và chứng minh toán học là gì?

Phát hiện bất ngờ nhất là ma trận trễ $D$ và hàm trễ $h(t)$ hoàn toàn không ảnh hưởng đến tính không phụ thuộc xung của hệ thống. Tính chất chính quy và không xung chỉ phụ thuộc duy nhất vào cặp ma trận $(E, A)$ thông qua tính khả nghịch của khối con $A_{22}$ trong biểu diễn SVD ($\det(A_{22}) \neq 0$). Điều này cho phép phân tách hoàn toàn bài toán phân tích cấu trúc đại số khỏi bài toán ước lượng động học trễ.

4. Luận án có cung cấp giao thức tái lập (Replication Protocol) rõ ràng không?

Có. Luận án cung cấp đầy đủ thuật toán đại số từng bước: từ việc tính toán các ma trận biến đổi $M, G$ qua lệnh svd(E) trong MATLAB, thiết lập các khối ma trận của bất đẳng thức LMI $\Gamma < 0$, khai báo biến trong môi trường feasp của LMI Control Toolbox, đến việc trích xuất ma trận nghiệm $P, U$ để tính độ lợi điều khiển $K = UP^{-\top}$. Mọi ví dụ số đều có đầy đủ ma trận dữ liệu đầu vào và kết quả kiểm chứng.

5. Kế hoạch nghiên cứu 10 năm tới được phác thảo với những cột mốc nào?

Kế hoạch nghiên cứu 10 năm (2021-2031) bao gồm:

  • Giai đoạn 1 (2021-2023): Hoàn thiện lý thuyết FTS cho hệ suy biến phi tuyến Takagi-Sugeno và hệ ngẫu nhiên.
  • Giai đoạn 2 (2024-2027): Phát triển hệ điều khiển thời gian thực trên nền tảng mạng truyền thông có độ trễ ngẫu nhiên và mất gói tin (Networked Descriptor Systems).
  • Giai đoạn 3 (2028-2031): Ứng dụng thương mại hóa vào việc điều khiển phi tập trung cho các hệ thống lưới điện vi mô (Microgrids) và mạng lưới robot tự hành hợp tác quy mô lớn.

Kết luận

Luận án tiến sĩ của tác giả Nguyễn Huyền Mười là một công trình khoa học mẫu mực, thể hiện tư duy toán học giải tích đỉnh cao kết hợp với kỹ thuật tính toán hiện đại. Tác phẩm đã đạt được 6 đóng góp mang tính cột mốc:

  1. Thiết lập hệ thống định nghĩa và tiêu chuẩn ổn định vững hữu hạn thời gian cho hệ phương trình vi phân suy biến có trễ hằng và trễ biến thiên liên tục.
  2. Chứng minh tường minh các điều kiện cần và đủ để hệ suy biến đóng vừa chính quy, không phụ thuộc xung, vừa đạt chuẩn chặn FTS dưới dạng hệ bất đẳng thức ma trận tuyến tính (LMIs).
  3. Đề xuất quy tắc giải tích thiết kế bộ điều khiển phản hồi trạng thái $u(t) = UP^{-\top}x(t)$ xử lý trễ không khả vi mà không cần giả thiết suy giảm của hàm trễ.
  4. Mở rộng trọn vẹn khung lý thuyết sang miền thời gian rời rạc và hệ suy biến rời rạc chuyển mạch có thành phần kích động phi tuyến.
  5. Xây dựng thuật toán tính toán số học tự động hóa hoàn toàn trên nền tảng MATLAB LMI-Lab với độ tin cậy và tính khả thi tuyệt đối.
  6. Mở ra các hướng nghiên cứu mới liên ngành giữa giải tích hàm, tối ưu hóa lồi và kỹ thuật điều khiển tự động hóa hiện đại, để lại giá trị học thuật và ứng dụng thực tiễn lâu dài trong kho tàng khoa học nước nhà và quốc tế.