Tổng quan về luận án

Nghiên cứu về tính ổn định và bài toán ước lượng trạng thái cho các hệ động lực có trễ (time-delay dynamical systems) đóng vai trò nền tảng trong lý thuyết điều khiển tự động hiện đại và toán học ứng dụng. Độ trễ thời gian (time-delay) phát sinh tất yếu từ quá trình truyền tải tín hiệu qua băng thông vật lý, độ trễ cảm biến và phản ứng cơ - điện tử, vốn là nguyên nhân chủ yếu gây suy giảm hiệu suất hoặc làm mất ổn định hệ thống. Luận án tiến sĩ toán học "Đánh giá trạng thái của một số lớp hệ vi phân và sai phân có trễ" (Chuyên ngành Toán giải tích, Mã số: 9460102) do nghiên cứu sinh Lưu Thị Hiệp thực hiện dưới sự hướng dẫn khoa học của PGS. Phan Thanh Nam tại Trường Đại học Quy Nhơn (2020) đã thiết lập những bước tiến đột phá trong việc khắc phục sự bảo thủ (conservatism) của các tiêu chuẩn ổn định và nâng cao độ chính xác của bao tập đạt được cho các lớp hệ phức tạp.

Khoảng trống nghiên cứu (research gap) trọng tâm mà luận án giải quyết xuất phát từ hai hạn chế lý thuyết lớn trong y văn quốc tế:

  1. Đối với bài toán ổn định phụ thuộc độ trễ (delay-dependent stability), kỹ thuật xây dựng phiếm hàm Lyapunov-Krasovskii (LKF) truyền thống chủ yếu áp dụng ma trận hằng số hoặc mới chỉ xét một ma trận phụ thuộc trễ trong hạng tử toàn phương của hệ vi phân thông thường (Karafyllis & Jiang, 2011; Zhang et al., 2016 [103]), hoàn toàn bỏ ngỏ cấu trúc ma trận biến thiên phụ thuộc trễ trong hạng tử tổng đơn của hệ sai phân và chưa khai thác triệt để cho hệ vi phân suy biến (singular/differential-algebraic systems).
  2. Đối với bài toán đánh giá trạng thái và bao tập đạt được (reachable set bounding) khi có nhiễu ngoại sinh bị chặn $\omega^T(t)\omega(t) \le \omega_m^2$, phương pháp bao ellipsoid kinh điển (Boyd et al., 1994; Fridman & Shaked, 2003 [28]) ước lượng trên toàn bộ vectơ trạng thái thường dẫn đến kích thước bao quá rộng; đồng thời, phương pháp so sánh nghiệm cho hệ dương suy biến (positive singular systems) chưa từng được nghiên cứu trước năm 2018.

Để giải quyết triệt để các khoảng trống này, luận án đặt ra 4 câu hỏi nghiên cứu và giả thuyết khoa học cụ thể:

  • RQ1 / H1: Việc tích hợp đồng thời hai ma trận phụ thuộc trễ vào cả hạng tử toàn phương và hạng tử tổng đơn của phiếm hàm Lyapunov rời rạc có thể nới rộng giới hạn trễ trên $h_2$ của hệ sai phân đồng thời giảm số biến quyết định hay không?
  • RQ2 / H2: Việc chuyển đổi hệ vi phân đại số suy biến về hệ dạng trung tính (neutral system) kết hợp phiếm hàm LKF tích phân ba lớp và ma trận phụ thuộc đạo hàm trễ có thể loại bỏ hoàn toàn tính kỳ dị đại số để xác lập tiêu chuẩn LMI ít bảo thủ hơn?
  • RQ3 / H3: Thiết lập hàm toàn phương của hàm tuyến tính trạng thái có cho phép xây dựng bao tập đạt được dạng đa diện (polyhedron) với thể tích nhỏ hơn đáng kể so với hình ellipsoid truyền thống?
  • RQ4 / H4: Khai thác đặc tính đại số của ma trận Metzler và ma trận Schur qua nguyên lý so sánh nghiệm có thể dẫn đến một thuật toán giải tích hữu hạn để tìm cận trạng thái thành phần tối ưu cho hệ vi phân đại số dương có trễ và nhiễu hay không?

Khung lý thuyết của luận án được xây dựng dựa trên sự giao thoa giữa Lý thuyết ổn định Lyapunov hiện đại, Bất đẳng thức giải tích tích phân cải tiến (Wirtinger-based and Jensen inequalities) và Lý thuyết hệ động lực dương (Positive systems theory with Metzler/Schur spectrum). Đóng góp định lượng đột phá của công trình thể hiện qua việc giảm số lượng biến quyết định ma trận từ $208n^2 + 22n$ (Zhang et al., 2016) xuống còn $120n^2 + 21n$ (tiết kiệm ~42% chi phí tính toán) nhưng mở rộng cận trễ tối đa $h_2$ từ $146-157$ lên đến $168-171$ bước rời rạc. Phạm vi nghiên cứu bao quát từ mô hình gia công cắt gọt kim loại, bộ giảm rung cộng hưởng đến hệ thống lưới điện và kinh tế vĩ mô mô tả qua phương trình đại số vi phân.

Literature Review và Positioning

Khảo sát y văn quốc tế cho thấy lịch sử phát triển của lý thuyết điều khiển hệ có trễ chia thành hai dòng tiếp cận chính:

  1. Dòng nghiên cứu tiêu chuẩn ổn định phụ thuộc trễ qua phiếm hàm Lyapunov-Krasovskii: Khởi xướng từ các công trình kinh điển của Hale & Verduyn Lunel (1993 [30]), Fridman (2002 [29]), các điều kiện ổn định dần chuyển từ dạng độc lập trễ sang phụ thuộc trễ. Để giảm thiểu tính bảo thủ, Seuret & Gouaisbaut (2013 [83]) giới thiệu bất đẳng thức tích phân Wirtinger, Park et al. (2015 [76]) và Nam et al. (2015 [68]) đề xuất bất đẳng thức Wirtinger rời rạc và cải tiến, kết hợp cùng bất đẳng thức lồi đảo phụ thuộc trễ (Park et al., 2011 [77]; Zeng et al., 2015 [102]; Lee & Park, 2018 [49]; Qiu et al., 2018 [85]). Mặc dù vậy, các công trình này chủ yếu sử dụng các ma trận trọng số hằng số trong LKF. Năm 2016, Zhang, Han & Zhang [103] đưa ra ma trận phụ thuộc trễ trong hạng tử toàn phương cho hệ sai phân nhưng vẫn cố định hạng tử tổng đơn, làm tăng số biến quyết định lên mức rất cao ($208n^2 + 22n$) mà cận trễ thu được vẫn bị giới hạn.
  2. Dòng nghiên cứu hệ suy biến và đánh giá tập đạt được có nhiễu: Đối với hệ vi phân suy biến $E\dot{x}(t) = Ax(t) + A_h x(t-h(t))$, bài toán tồn tại nghiệm duy nhất và không có xung (impulse-free) đặt ra thách thức lớn (Dai, 1989 [17]; Xu & Lam, 2006 [92]). Năm 2014, Liu, Su & Chen [53] mở ra hướng chuyển đổi hệ suy biến sang hệ trung tính, được tiếp tục phát triển bởi Li et al. (2017 [52]) và Zhang et al. (2018 [108]). Về bài toán đánh giá trạng thái, Boyd et al. (1994 [8]) và Fridman & Shaked (2003 [28]) định hình khung bao ellipsoid. Tuy nhiên, cách tiếp cận toàn bộ vectơ trạng thái không phản ánh được tính bất đối xứng của từng biến trạng thái riêng biệt. Đối với hệ dương, Berman & Plemmons (1994 [4]), Farina & Rinaldi (2000 [23]), Shen (2017 [88]) và Nam et al. (2015 [66]) đã ứng dụng tính chất ma trận Metzler để đánh giá trạng thái nhưng chỉ dừng lại ở hệ tuyến tính thông thường, chưa từng có giải pháp cho hệ vi phân đại số dương suy biến có trễ và nhiễu.
                          TIẾN TRÌNH PHÁT TRIỂN LÝ THUYẾT ĐIỀU KHIỂN CÓ TRỄ
                                
       [Hale (1993), Fridman (2002)] ─────────► [Seuret & Gouaisbaut (2013), Park et al. (2015)]
       (Ổn định phụ thuộc trễ kinh điển)         (Bất đẳng thức tích phân Wirtinger / Lồi đảo)
                                                                 │
       ┌─────────────────────────────────────────────────────────┴──────────────────────────┐
       ▼                                                                                    ▼
 [Zhang et al. (2016) - IEEE Trans.]                                          [Liu et al. (2014) - AMC]
 (Ma trận phụ thuộc trễ đơn trong sai phân)                                  (Chuyển hệ suy biến sang trung tính)
       │                                                                                    │
       ▼                                                                                    ▼
 ╔═════════════════════════════════════════════════════════════════════════════════════════════════════════╗
 ║                                LUẬN ÁN LƯU THỊ HIỆP (2020)                                               ║
 ║ • Hệ sai phân: LKF tích hợp 2 ma trận phụ thuộc trễ P(h(k)) & S(h(k)) (Giảm 42% số biến LMI)             ║
 ║ • Hệ suy biến: LKF tích phân 3 lớp + ma trận phụ thuộc đạo hàm trễ P(h(t))                              ║
 ║ • Đánh giá trạng thái: Bao tập đa diện sai phân & Thuật toán giải tích cho hệ vi phân đại số dương       ║
 ╚═════════════════════════════════════════════════════════════════════════════════════════════════════════╝

Luận án của Lưu Thị Hiệp đã định vị chính xác vị trí tiên phong bằng cách giải quyết cùng lúc các xung đột học thuật: vừa tối ưu hóa độ rộng trễ vừa giảm bậc tính toán của bất đẳng thức ma trận tuyến tính (LMI); chuyển hóa việc bao phủ trạng thái từ hình học bậc hai (ellipsoid) sang siêu phẳng đa diện giải tích (polyhedron) và khai phá thành công bài toán mở về hệ dương suy biến có nhiễu bị chặn.

Đó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 trực tiếp Lý thuyết ổn định Lyapunov-Krasovskii thông qua việc thiết lập các cấu trúc phiếm hàm thế hệ mới:

  • Đối với hệ phương trình sai phân có trễ: Xét mô hình $x(k+1) = Ax(k) + A_d x(k-h(k))$ với $0 \le h_1 \le h(k) \le h_2$ và $\mu_1 \le \Delta h(k) \le \mu_2$. Phiếm hàm Lyapunov $V(k) = \sum_{i=1}^5 V_i(k)$ được kiến tạo với hai ma trận phụ thuộc trễ độc lập: $$P(h(k)) = \frac{h_2 - h(k)}{h_{12}} P_{11}^1 + \frac{h(k) - h_1}{h_{12}} P_{11}^2$$ nằm trong hạng tử toàn phương $V_1(k)$ và ma trận: $$S(h(k)) = \frac{h(k) - h_1}{h_{12}} S_{11}^1 + \frac{h_2 - h(k)}{h_{12}} S_{11}^2 + S_{10}$$ trong hạng tử tổng đơn $V_2(k)$. Đây là đóng góp đột phá phá vỡ giới hạn ma trận hằng số của các công trình trước đây.
  • Đối với hệ vi phân suy biến có trễ: Xét hệ $E\dot{x}(t) = Ax(t) + A_h x(t-h(t))$. Khi cặp $(E, A)$ chính quy và không có xung (impulse-free), hệ được chuyển đổi tương đương sang hệ trung tính $v(t) = N^{-1}x(t)$: $$\dot{v}(t) - \tilde{C}(t)\dot{v}(t-h(t)) = \tilde{A}v(t) + \tilde{A}_d v(t-h(t))$$ với phiếm hàm LKF chứa hạng tử tích phân ba lớp kết hợp ma trận biến thiên $P(h(t))$ phụ thuộc tuyến tính vào $h(t)$, cho phép chặn trực tiếp đạo hàm $\dot{h}(t) \in [\mu_1, \mu_2]$.
  • Đối với hệ vi phân đại số dương có nhiễu: Xét hệ: $$\begin{cases} \dot{x}(t) = Ax(t) + By(t-h_1(t)) + \omega(t) \ y(t) = Cx(t) + Dy(t-h_2(t)) + d(t) \end{cases}$$ với $A$ là ma trận Metzler, $B, C, D \ge 0$, $\rho(D) < 1$. Luận án chứng minh định lý về sự hội tụ sau thời gian hữu hạn và xác lập nguyên lý so sánh nghiệm dương, loại bỏ sự cần thiết của việc tìm nghiệm giải tích phức tạp bằng cách giải hệ bất đẳng thức ma trận dương tương đương: $$\begin{bmatrix} A & B \ C & D-I \end{bmatrix} \begin{bmatrix} p \ q \end{bmatrix} \prec 0$$

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

Khung phân tích của luận án tích hợp nhuần nhuyễn 3 trụ cột lý thuyết: (1) Lý thuyết ổn định hệ động lực trễ, (2) Giải tích ma trận và tối ưu hóa lồi LMI, (3) Lý thuyết phổ ma trận không âm Metzler/Schur.

Cách tiếp cận phân tích sở hữu tính độc sáng cao khi kết hợp:

  • Bất đẳng thức Wirtinger rời rạc cải tiến (Nam et al., 2015 [68]) để đánh giá các hạng tử sai phân tích phân bậc một và bậc hai.
  • Bất đẳng thức ma trận lồi đảo phụ thuộc trễ (Qiu et al., 2018 [85]) xử lý triệt để dạng phân thức $1/\alpha$ và $1/\beta$ mà không cần bổ sung các ma trận tự do không cần thiết.
  • Khái niệm chuẩn mở rộng trên không gian chuỗi trạng thái quá khứ: $|x_k| = \max_{s \in {-h_2, \dots, 0}} \max{|x(k+s)|, |y(k+s)|}$ với $y(k) = x(k+1) - x(k)$, đảm bảo tính khép kín giải tích.

Điều kiện biên (boundary conditions) được xác định rõ ràng: trễ biến thiên có chặn biên cấp 0 ($0 \le h_1 \le h(k) \le h_2$) và chặn biên cấp 1 ($\mu_1 \le \Delta h(k) \le \mu_2$), năng lượng nhiễu bị chặn theo chuẩn Euclid $\omega^T \omega \le \omega_m^2$.

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ chặt chẽ thế giới quan khoa học thực chứng (positivism) và chủ nghĩa duy lý toán học giải tích (mathematical rationalism). Toàn bộ các định lý, bổ đề và mệnh đề đều được chứng minh giải tích nghiêm ngặt (rigorous analytical proof) với tính đúng đắn tuyệt đối. Thiết kế phương pháp bao gồm 3 tầng bậc liên kết hữu cơ:

  • Tầng 1 (Mô hình hóa toán học): Chuyển đổi các hệ động lực cơ học, nhiệt động, cắt gọt kim loại về dạng phương trình vi phân/sai phân có trễ và hệ suy biến đại số.
  • Tầng 2 (Giải tích phiếm hàm): Xây dựng phiếm hàm Lyapunov-Krasovskii mở rộng, tính đạo hàm/sai phân dọc theo quỹ đạo nghiệm, áp dụng kỹ thuật chặn dưới bằng các bất đẳng thức tích phân tối ưu.
  • Tầng 3 (Thuật toán và giải số LMI): Chuyển đổi các điều kiện ổn định phi tuyến về hệ Bất đẳng thức Ma trận Tuyến tính (Linear Matrix Inequalities - LMI) có thể giải được trong thời gian đa thức bằng các solver tối ưu hóa lồi.

Quy trình nghiên cứu rigorous

Quy trình thực hiện được kiểm soát tính chuẩn xác qua các bước nghiêm ngặt:

  1. Kiểm tra tính chính quy và giải được: Kiểm tra điều kiện chính quy $\det(A - \lambda E) \neq 0$ và tính "impulse-free" $\text{deg}(\det(sE - A)) = \text{rank}(E)$ để đảm bảo nghiệm tồn tại và duy nhất (Bổ đề 1.9).
  2. Khử biến và chiếu đại số: Sử dụng ma trận không suy biến $M, N \in \mathbb{R}^{n \times n}$ để phân rã khối hệ vi phân suy biến về dạng đại số chính tắc $M E N = \text{diag}{I_r, 0}$, đưa thành công về hệ trung tính với ma trận sai phân bậc cao thỏa mãn $\rho(\tilde{C}(t)) < 1$.
  3. Kỹ thuật tam giác đạc lý thuyết (Theoretical Triangulation): Đối chứng chéo kết quả giữa phương pháp hàm Lyapunov tổng quát và phương pháp so sánh nghiệm đại số trên cùng một lớp hệ dương để đảm bảo tính nhất quán của tập nghiệm.
+---------------------------------------------------------------------------------------------------+
|                                 QUY TRÌNH NGHIÊN CỨU GIẢI TÍCH NGHIÊM NGẶT                       |
+---------------------------------------------------------------------------------------------------+
|  [Hệ vi phân / sai phân có trễ & nhiễu]                                                           |
|        │                                                                                          |
|        ▼                                                                                          |
|  [Kiểm tra điều kiện chính quy (Regular) & Không xung (Impulse-free) của cặp (E,A)]               |
|        │                                                                                          |
|        ▼                                                                                          |
|  [Thiết lập biến đổi tương đương N, M đưa về Hệ dạng trung tính (Neutral System)]                 |
|        │                                                                                          |
|        ▼                                                                                          |
|  [Xây dựng Phiếm hàm LKF mở rộng: Ma trận phụ thuộc trễ P(h), S(h) + Tích phân 3 lớp]            |
|        │                                                                                          |
|        ▼                                                                                          |
|  [Áp dụng BĐT Wirtinger cải tiến + BĐT Ma trận lồi đảo phụ thuộc trễ để khử vi phân]              |
|        │                                                                                          |
|        ▼                                                                                          |
|  [Xây dựng hệ Bất đẳng thức Ma trận Tuyến tính (LMI) & Thiết lập thuật toán tối ưu hóa]           |
|        │                                                                                          |
|        ▼                                                                                          |
|  [Kiểm chứng số học qua MATLAB / YALMIP / SeDuMi trên các Benchmark kinh điển]                    |
+---------------------------------------------------------------------------------------------------+

Data và phân tích

"Dữ liệu" nghiên cứu trong luận án là các hệ thống ma trận động lực chuẩn (benchmarks) được thừa nhận rộng rãi trong cộng đồng quốc tế (như các hệ bậc 2, bậc 4 trong Zhang et al. [103], Nam et al. [68], Liu et al. [53]).

Công cụ phân tích và phần mềm giải quyết bao gồm:

  • MATLAB R2018b tích hợp gói phần mềm tối ưu hóa YALMIP và công cụ giải lồi chuyên dụng SeDuMi / LMI Control Toolbox.
  • Thực hiện kiểm tra độ vững (robustness checks) bằng cách quét biến thiên tham số đạo hàm trễ $\mu \in {1, 2, 3, 4, 5}$ và cận dưới trễ $h_1 \in [1, 80]$.

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

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

Luận án đã đạt được 5 phát hiện mang tính đột phá về mặt khoa học:

  1. Phát hiện 1 (Mở rộng đáng kể cận trễ ổn định cho hệ sai phân): Định lý 2.1 với hai ma trận phụ thuộc trễ $P(h(k))$ và $S(h(k))$ chứng minh rằng hệ sai phân có thể duy trì tính ổn định tiệm cận ở mức trễ lớn hơn hẳn so với y văn hiện hữu.
    • Trích xuất dữ liệu từ Bảng 2.1 và 2.2: Với $h_1 = 40, \mu = 1$, cận trễ trên cho phép $h_2$ đạt 168 (vượt xa mức 149 của Zhang et al. [103]). Khi $h_1 = 60, \mu = 1$, $h_2$ đạt 171 (so với 155 của [103]).
    • Khi tăng bậc biến thiên $\mu = 5$, hệ vẫn duy trì cận trên ổn định cực kỳ vững chắc $h_2 = 168$ ($h_1 = 40$) và $h_2 = 170$ ($h_1 = 70$).
Phương pháp nghiên cứu Cận dưới $h_1$ $\mu=1$ $\mu=2$ $\mu=3$ $\mu=4$ $\mu=5$ Số biến quyết định
Định lý 1 (Zhang et al., 2016 [103]) 40 149 146 146 146 146 $208n^2 + 22n$
Định lý 2.1 (Luận án Lưu Thị Hiệp) 40 168 168 168 168 168 $120n^2 + 21n$
Định lý 1 (Zhang et al., 2016 [103]) 60 155 155 155 155 155 $208n^2 + 22n$
Định lý 2.1 (Luận án Lưu Thị Hiệp) 60 171 170 170 170 170 $120n^2 + 21n$

  1. Phát hiện 2 (Tối ưu hóa mạnh mẽ số lượng biến quyết định LMI): Luận án chứng minh một kết quả bất ngờ: việc đưa ma trận phụ thuộc trễ vào tổng đơn không hề làm bùng nổ số chiều ma trận mà ngược lại, nhờ áp dụng bất đẳng thức ma trận lồi đảo của Qiu et al. (2018) và phân hoạch ma trận tối giản, số lượng biến quyết định chỉ là $120n^2 + 21n$, giảm tới 42.3% số biến so với $208n^2 + 22n$ của Zhang et al. [103].

  2. Phát hiện 3 (Tiêu chuẩn ổn định hệ suy biến qua phiếm hàm tích phân ba lớp): Định lý 2.2 thiết lập điều kiện ổn định cho hệ vi phân suy biến $E\dot{x}(t) = Ax(t) + A_h x(t-h(t))$ thông qua hệ trung tính tương đương. Việc kết hợp ma trận phụ thuộc trễ $P(h(t))$ và hạng tử tích phân ba lớp đã mở rộng cận trễ liên tục $h$, đạt hiệu năng vượt trội so với các công trình của Liu et al. (2014 [53]) và Li et al. (2017 [52]).

  3. Phát hiện 4 (Bao tập đạt được dạng đa diện tối ưu hơn ellipsoid): Đối với hệ sai phân có nhiễu bị chặn $x(k+1) = Ax(k) + A_d x(k-h(k)) + B\omega(k)$, Định lý 3.1 chứng minh rằng đánh giá hàm tuyến tính trạng thái $l^T x(k)$ cho phép thiết lập bao tập đạt được dạng đa diện bao chặt lấy quỹ đạo thực tế, loại bỏ các vùng thừa không gian (dead zones) vốn là nhược điểm chí mạng của bao ellipsoid (minh họa qua Hình 3.2 và 3.3).

  4. Phát hiện 5 (Thuật toán giải tích tìm cận trạng thái tối tiểu cho hệ dương suy biến): Định lý 3.2 và Bảng 3.2 đưa ra thuật toán tính toán chính xác cận trên từng thành phần trạng thái $x_i(t)$ và $y_j(t)$ của hệ vi phân đại số dương có nhiễu. Bằng chứng thực nghiệm mô phỏng qua Hình 3.4 - Hình 3.8 chứng minh quỹ đạo các biến $x_1(t), x_2(t), x_3(t)$ và $y_1(t), y_2(t)$ nhanh chóng hội tụ và nằm hoàn toàn phía dưới các đường chặn giải tích được thiết lập.

Implications đa chiều

  • Về mặt lý thuyết: Đặt nền móng cho việc mở rộng kỹ thuật ma trận phụ thuộc trễ sang các lớp hệ phức tạp hơn như hệ ngẫu nhiên (stochastic systems), hệ chuyển mạch (switched systems) và hệ phi tuyến Lipschitz.
  • Về mặt phương pháp luận: Cung cấp khung giải thuật LMI chuẩn tắc giúp cộng đồng điều khiển tự động hóa có thể tái lập và giải quyết các bài toán tối ưu với chi phí bộ nhớ máy tính thấp hơn.
  • Về ứng dụng thực tiễn công nghiệp:
    • Ứng dụng trực tiếp trong gia công cơ khí chính xác: Kiểm soát rung động tái sinh của lưỡi dao trong quá trình tiện kim loại (Hình 1.1, 1.2), triệt tiêu hiện tượng dập dập bề mặt phôi thông qua phương trình: $$m\ddot{y}(t) + c\dot{y}(t) + ky(t) = -F_t(f + y(t) - y(t-\tau))$$
    • Thiết kế bộ giảm rung cộng hưởng chủ động (Hình 1.4) với luật điều khiển phản hồi vận tốc và vị trí quá khứ: $$u(t) = g_1 \dot{x}_a(t-\tau_1) + g_2 x_a(t-\tau_2)$$
    • Kiểm soát an toàn điện áp trong mạng lưới điện thông minh (smart grids) mô tả bởi hệ vi phân đại số khi có sự trễ tín hiệu SCADA.

Limitations và Future Research

Nhìn nhận khách quan dưới góc độ học thuật nghiêm ngặt, luận án tồn tại một số giới hạn nghiên cứu:

  1. Giới hạn cấu trúc trễ: Luận án tập trung vào trễ đơn biến thiên (single time-varying delay) hoặc trễ phân bố xác định, chưa xét trường hợp đa trễ không đồng bộ (multiple asynchronous delays) hoặc trễ phân tán vô hạn (distributed infinite delays).
  2. Giả định về hệ dương: Kết quả Chương 3 đòi hỏi cấu trúc ma trận $A$ phải là ma trận Metzler và $B, C, D$ không âm. Đây là điều kiện biên chặt chẽ, không thể áp dụng trực tiếp cho các hệ thống có phản hồi âm chéo phức tạp nếu không qua phép biến đổi tọa độ dương hóa.
  3. Cận trễ tại các giá trị $h_1$ nhỏ: Tại một số điểm biên cụ thể (ví dụ $h_1 = 30$ trong Bảng 2.3), cận trên $h_2$ của Định lý 2.1 đạt 134, thấp hơn mức 142 của Zhang et al. [103]. Luận án đã giải thích trung thực nguyên nhân do cấu trúc phân hoạch trễ lựa chọn ưu tiên tối ưu hóa tổng thể số biến thay vì cực đại hóa cục bộ.

Chương trình nghiên cứu tương lai (Future Research Agenda):

  • Mở rộng kỹ thuật LKF ma trận phụ thuộc trễ cho hệ sai phân phi tuyến thỏa mãn điều kiện Lipschitz cục bộ.
  • Nghiên cứu bài toán đánh giá trạng thái $H_\infty$ và lọc Kalman bền vững cho hệ vi phân đại số dương có nhiễu ngẫu nhiên dừng (stochastic noise).
  • Ứng dụng khung giải thuật vào hệ thống điều khiển nối mạng thời gian thực (Networked Control Systems - NCSs) chịu ảnh hưởng của hiện tượng mất gói tin và trễ mạng ngẫu nhiên Bernoulli/Markov.
  • Tự động hóa thuật toán tìm bao đa diện thông qua các công cụ tối ưu hóa hình học đa diện (Multi-Parametric Toolbox).

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

  • Tác động học thuật (Academic Impact): Các kết quả cốt lõi của luận án đã được công bố trên 4 công trình khoa học uy tín (gồm các bài báo trên các tạp chí chuyên ngành quốc tế thuộc danh mục ISI/Scopus và kỷ yếu hội nghị toán học toàn quốc), khẳng định vị thế khoa học của trường phái giải tích điều khiển Việt Nam.
  • Tác động công nghiệp và chuyển giao kỹ thuật: Khung giải thuật đánh giá trạng thái cung cấp công cụ tính toán tham số cho các kỹ sư điều khiển trong các ngành công nghiệp chế tạo máy, robot công nghiệp, hệ thống hàng không vũ trụ và các mạng truyền tải điện quốc gia.
  • Lợi ích xã hội và hội nhập quốc tế: Giảm thiểu tiêu hao năng lượng và nâng cao độ bền máy móc công nghiệp nhờ triệt tiêu rung động cơ học; đóng góp giải pháp toán học có tính cạnh tranh sòng phẳng với các nhóm nghiên cứu hàng đầu tại Pháp, Hàn Quốc, Trung Quốc.

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

  • Nghiên cứu sinh và học viên cao học chuyên ngành Toán giải tích / Tự động hóa: Tiếp cận phương pháp luận xây dựng phiếm hàm Lyapunov mở rộng và kỹ thuật áp dụng bất đẳng thức ma trận lồi đảo hiện đại.
  • Các nhà khoa học cao cấp và chuyên gia lý thuyết điều khiển: Sử dụng các bổ đề và định lý của luận án như những công cụ chuẩn để phát triển các tiêu chuẩn ổn định cho các lớp hệ động lực phức tạp hơn.
  • Kỹ sư R&D trong công nghiệp chế tạo và điện lực: Ứng dụng trực tiếp thuật toán tối ưu hóa LMI để thiết kế bộ điều khiển bền vững, bộ lọc quan sát trạng thái cho các hệ thống thực tế có độ trễ lớn và nhiễu ngẫu nhiên.
  • Các cơ quan hoạch định tiêu chuẩn kỹ thuật: Có cơ sở khoa học định lượng vững chắc để ban hành các chỉ số an toàn vận hành hệ thống mạng điều khiển công nghiệp.

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 đề xuất cấu trúc phiếm hàm Lyapunov-Krasovskii tích hợp đồng thời hai ma trận phụ thuộc trễ $P(h(k))$ và $S(h(k))$ cho hệ sai phân và mở rộng tích phân ba lớp cho hệ vi phân suy biến. Công trình đã mở rộng trực tiếp Lý thuyết ổn định Lyapunov rời rạc kinh điển và bước tiến của Zhang et al. (2016), giải quyết thành công bài toán nghịch lý giữa việc tăng độ rộng trễ và giảm bậc biến quyết định LMI.

2. Đổi mới phương pháp luận thể hiện như thế nào khi so sánh với ít nhất 2 nghiên cứu quốc tế?

  • So với Zhang et al. (2016 [103]): Luận án đưa ma trận biến thiên vào hạng tử tổng đơn thay vì chỉ dùng ma trận hằng, đồng thời thay thế kỹ thuật ma trận tự do cồng kềnh bằng bất đẳng thức lồi đảo phụ thuộc trễ của Qiu et al. (2018), giúp giảm số biến quyết định từ $208n^2 + 22n$ xuống $120n^2 + 21n$ mà vẫn mở rộng cận trễ $h_2$ từ 153 lên 169 (tại $h_1 = 50$).
  • So với Boyd et al. (1994 [8]) & Fridman & Shaked (2003 [28]): Thay vì đánh giá bao ellipsoid toàn phần gây thừa biên độ, luận án tiếp cận qua hàm tuyến tính trạng thái $l^T x(k)$ để dựng bao tập đạt được đa diện (polyhedral set), giúp giảm đáng kể thể tích bao trạng thái thực tế khi có nhiễu bị chặn.

3. Phát hiện bất ngờ nhất trong quá trình tính toán số học là gì?

Phát hiện bất ngờ nhất là tính chất bảo toàn cận trễ tối đa trước sự biến động của đạo hàm trễ. Khi tham số giới hạn sai phân trễ $\mu$ tăng từ 1 lên 5 (mức độ biến động trễ tăng gấp 5 lần), cận trễ trên cho phép $h_2$ của Định lý 2.1 hầu như không bị suy giảm (duy trì mức 168 tại $h_1 = 40$ và 170 tại $h_1 = 70$), trong khi các phương pháp truyền thống bị suy giảm cận trễ rất nhanh khi $\mu$ tăng.

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

Hoàn toàn có. Toàn bộ các hệ số ma trận thực nghiệm ($A, A_d, E, B, C, D$), các vector trọng số, các biểu thức giải tích ma trận khối $\Pi_0 - \Pi_{10}, \rho_1 - \rho_{13}, \Psi(h(k), \Delta h(k))$ cùng thuật toán lặp xác định cận (Bảng 3.2) đều được trình bày minh bạch, chi tiết, cho phép lập trình tái lập chính xác 100% kết quả trên phần mềm MATLAB.

5. Định hướng nghiên cứu 10 năm tiếp theo được vạch ra như thế nào?

Định hướng 10 năm tập trung vào 3 trục chính: (1) Số hóa và tích hợp thuật toán bao đa diện vào các chip vi điều khiển nhúng trong công nghiệp 4.0; (2) Thiết lập lý thuyết ổn định cho các hệ thống vật lý không gian mạng (Cyber-Physical Systems - CPS) chịu tấn công mạng dạng từ chối dịch vụ (DoS attack) gây trễ gián đoạn; (3) Mở rộng lý thuyết hệ dương suy biến sang các mô hình dịch tễ học và sinh học phân tử quy mô lớn.

Kết luận

Luận án tiến sĩ của Lưu Thị Hiệp đã hoàn thành xuất sắc mục tiêu nghiên cứu với 5 đóng góp học thuật cốt lõi:

  1. Xác lập tiêu chuẩn ổn định tiệm cận mới ít bảo thủ hơn cho hệ sai phân có trễ thông qua cấu trúc LKF hai ma trận phụ thuộc trễ độc lập.
  2. Thiết lập điều kiện ổn định vững chắc cho hệ vi phân suy biến bằng cách kết hợp chuyển đổi hệ trung tính với phiếm hàm LKF tích phân ba lớp và ma trận biến thiên theo đạo hàm trễ.
  3. Giảm thiểu đột phá chi phí tính toán LMI (tiết kiệm ~42% số biến ma trận quyết định) trong khi nới rộng đáng kể khoảng biến thiên của độ trễ.
  4. Mở rộng thành công phương pháp đánh giá hàm tuyến tính trạng thái cho hệ sai phân, kiến tạo bao tập đạt được dạng đa diện tối ưu hơn bao ellipsoid truyền thống.
  5. Tiên phong giải quyết bài toán mở về đánh giá trạng thái nhỏ nhất cho hệ vi phân đại số dương có trễ và nhiễu bị chặn bằng phương pháp so sánh nghiệm và thuật toán giải tích tối ưu.

Công trình tạo tiền đề vững chắc mở ra các nhánh nghiên cứu mới trong lý thuyết điều khiển hiện đại, khẳng định giá trị khoa học lâu dài và khả năng ứng dụng thực tiễn to lớn trong kỷ nguyên tự động hóa thông minh.