Tổng quan về luận án

Nghiên cứu của NCS Lê Anh Tuấn tại Khoa Toán - Tin, Trường Đại học Sư phạm Hà Nội và Trường Đại học Khoa học Huế tập trung giải quyết bài toán cốt lõi trong lý thuyết điều khiển hiện đại: Điều khiển $H_\infty$ cho một số lớp hệ phương trình có trễ. Trong kỹ thuật tự động hóa, xử lý tín hiệu và mạng nơ-ron nhân tạo, độ trễ thời gian (time-delay) và nhiễu ngoại sinh (external disturbances) là những hiện tượng khách quan tất yếu làm suy giảm nghiêm trọng hiệu suất động học, thậm chí phá vỡ tính ổn định của hệ thống.

                  ┌────────────────────────────────────────┐
                  │          Nhiễu ngoại sinh ω(t)         │
                  └───────────────────┬────────────────────┘
                                      │
                                      ▼
┌─────────────────┐  u(t)     ┌────────────────┐  z(t)    ┌────────────────┐
│  Bộ điều khiển  ├──────────►│   Hệ có trễ    ├─────────►│ Đầu ra đánh giá│
│   phản hồi K    │           │ (Neural/Linear/│          │ (Error Output) │
└────────▲────────┘           │   Descriptor)  │          └────────────────┘
         │                    └───────┬────────┘
         │       Trạng thái x(t)      │
         └────────────────────────────┘

Khoảng trống học thuật (research gap) trọng tâm mà công trình nhận diện là: phần lớn các công trình kinh điển (Kwon & Park, 2009; Phat & Trinh, 2010; Sakthivel et al., 2012) đều đặt giả thiết hàm trễ phải khả vi và có đạo hàm bị chặn trên ($\dot{h}(t) \le \delta < 1$), hoặc cận dưới của trễ bắt buộc bằng $0$. Điều này loại trừ các hệ thống thực tế có trễ biến thiên nhanh không khả vi. Đồng thời, bài toán điều khiển $H_\infty$ trong thời gian hữu hạn (Finite-Time Stability - FTS) đối với hệ rời rạc suy biến (singular descriptor systems) có trễ dạng khoảng vẫn còn là bài toán mở phức tạp chưa có lời giải trọn vẹn.

Nghiên cứu xây dựng và kiểm chứng ba hệ thống giả thuyết khoa học chính:

  1. Hypothesis $H_1$: Tồn tại bộ điều khiển phản hồi trạng thái $u(t) = -\frac{1}{2}B^T P^{-1} x(t)$ bảo đảm hệ nơ-ron có trễ biến thiên hỗn hợp (rời rạc và tích phân phân phối) đạt tính $\alpha$-ổn định mũ và thỏa mãn mức triệt tiêu nhiễu $H_\infty$ tối ưu $\gamma$ mà không yêu cầu hàm trễ $h(t)$ phải khả vi.
  2. Hypothesis $H_2$: Đối với hệ rời rạc tuyến tính có trễ biến thiên dạng khoảng ($0 < d_1 \le d(k) \le d_2$), tồn tại điều kiện đủ dạng Bất đẳng thức Ma trận Tuyến tính (LMI) xác lập tính $H_\infty$-bị chặn trong thời gian hữu hạn.
  3. Hypothesis $H_3$: Đối với hệ nơ-ron rời rạc suy biến $Ex(k+1) = Ax(k) + \dots$, việc kết hợp phiếm hàm Lyapunov–Krasovskii cải tiến với Định lý hàm ẩn giải quyết đồng thời bài toán điều khiển $H_\infty$ hữu hạn và đảm bảo tính chính quy (regularity), tính nhân quả (causality) cùng tính duy nhất nghiệm.

Khung lý thuyết nền tảng tích hợp Lý thuyết ổn định Lyapunov–Krasovskii, Lý thuyết không gian Hardy $H_\infty$ và Phương pháp giải tối ưu lồi LMI (Boyd et al., 1994; Gahinet et al., 1995). Nghiên cứu định lượng hóa chỉ tiêu năng lượng thông qua mở rộng chỉ số suy giảm nhiễu: $$|z|{2}^2 \le \gamma \left(c_0 |\phi|{C^1}^2 + |\omega|_{2}^2\right), \quad \forall \omega \in L_2([0, \infty), \mathbb{R}^r)$$ Công trình bao quát từ hệ nơ-ron liên tục $n$-chiều, hệ rời rạc tuyến tính cho đến hệ phương trình sai phân đại số suy biến, mở ra bước tiến vượt bậc cho thiết kế bộ điều khiển bền vững.

Literature Review và Positioning

Lý thuyết ổn định hệ thống bắt nguồn từ công trình nền tảng của nhà toán học Nga A. Lyapunov (1892). Suốt hơn một thế kỷ, nhánh nghiên cứu này phát triển mạnh mẽ và trở thành công cụ không thể thiếu trong phân tích động lực học (Hale & Verduyn Lunel, 1993; Gu et al., 2003). Trong miền thời gian, phương pháp phiếm hàm Lyapunov–Krasovskii (LKF) trực tiếp đóng vai trò chủ đạo để khảo sát hệ phương trình vi phân hàm có trễ (RFDE).

                      TIẾN TRÌNH PHÁT TRIỂN LÝ THUYẾT ĐIỀU KHIỂN H∞
┌────────────────────────────────────────────────────────────────────────────────────────┐
│ Giai đoạn 1: Trễ khả vi & H∞ tiệm cận vô hạn                                          │
│ • Kwon & Park (2009): Ổn định mũ mạng nơ-ron, trễ khả vi đạo hàm chặn                 │
│ • Phat & Trinh (2010): Ổn định hóa hệ nơ-ron trễ hỗn hợp, yêu cầu h'(t) <= delta < 1  │
│ • Sakthivel et al. (2012): Điều khiển H∞ trễ hỗn hợp, chưa xét trễ ở đầu ra z(t)      │
└──────────────────────────────────────────┬─────────────────────────────────────────────┘
                                           │
                                           ▼
┌────────────────────────────────────────────────────────────────────────────────────────┐
│ Giai đoạn 2: Trễ hằng & Hệ rời rạc chuyển mạch hữu hạn                                 │
│ • Wang et al. (2010): Điều khiển H∞ thời gian hữu hạn cho hệ không trễ                │
│ • Xiang & Xiao (2011), Song et al. (2012): Hệ chuyển mạch trễ hằng                     │
│ • Ma & Zheng (2016): Hệ nơ-ron suy biến rời rạc bước nhảy Markov                      │
└──────────────────────────────────────────┬─────────────────────────────────────────────┘
                                           │
                                           ▼
┌────────────────────────────────────────────────────────────────────────────────────────┐
│ Giai đoạn 3: Đột phá của Luận án NCS Lê Anh Tuấn                                       │
│ • Trễ biến thiên nhanh không khả vi nhận giá trị trong khoảng [h1, h2]                │
│ • Xuất hiện trễ hỗn hợp đồng thời ở trạng thái x(t) và đầu ra quan sát z(t)            │
│ • Giải quyết trọn vẹn FTS, Regularity, Causality cho hệ nơ-ron rời rạc suy biến       │
│ • Loại bỏ hoàn toàn ma trận trọng số tự do (Free-weighting matrices), tối ưu LMI       │
└────────────────────────────────────────────────────────────────────────────────────────┘

Trong dòng nghiên cứu hệ nơ-ron có trễ, Kwon & Park (2009) thiết lập tiêu chuẩn ổn định mũ cho hệ nơ-ron với hàm trễ liên tục có đạo hàm bị chặn. Phat & Trinh (2010) mở rộng cho hệ nơ-ron có trễ hỗn hợp (mixed delays gồm trễ rời rạc $h(t)$ và trễ phân phối $k(t)$), nhưng vẫn yêu cầu điều kiện ngặt $\dot{h}(t) \le \delta < 1$. Đến năm 2012, Sakthivel và cộng sự nghiên cứu bài toán $H_\infty$ cho hệ nơ-ron trễ hỗn hợp nhưng giả định đầu ra quan sát $z(t) = Ex(t)$ hoàn toàn không chứa thành phần trễ $x(t-h(t))$ và hàm trễ vẫn phải khả vi. Phat & Trinh (2013) tiếp tục phân tích trường hợp trễ không khả vi nhưng chỉ giới hạn ở mô hình trễ rời rạc thuần túy, chưa tích hợp đồng thời trễ tích phân phân phối và trễ trong hàm quan sát.

Ở nhánh nghiên cứu ổn định trong thời gian hữu hạn (Finite-Time Stability - FTS), khởi xướng từ công trình của Dorato (1961) và Amato et al. (2001), trọng tâm chuyển từ dáng điệu tiệm cận tại vô hạn ($t \to \infty$) sang việc giới hạn biên độ quỹ đạo trạng thái trong một miền cầu xác định suốt khoảng thời gian hữu hạn $[0, T]$. Năm 2010, Wang và cộng sự khảo sát bài toán $H_\infty$ hữu hạn cho hệ rời rạc không trễ; Xiang & Xiao (2011) cùng Song et al. (2012) mở rộng cho hệ rời rạc chuyển mạch nhưng chỉ dừng lại ở trễ hằng ($d(k) \equiv d$). Các công trình của Yan et al. (2013) và Stojanovic et al. (2014) khảo sát ổn định hóa hệ rời rạc trễ khoảng nhưng chưa giải quyết chỉ tiêu triệt tiêu nhiễu $H_\infty$.

Đối với hệ suy biến (singular/descriptor systems) được mô tả bởi phương trình sai phân đại số $Ex(k+1) = Ax(k) + \dots$ với $\det(E) = 0$, bài toán tồn tại nghiệm duy nhất, tính không xung (impulse-free), tính chính quy (regularity) và nhân quả (causality) tạo nên rào cản toán học to lớn (Dai, 1989; Campbell, 1982). Loạt nghiên cứu của Zhang et al. (2014) và Bai et al. (2015) chỉ mới khảo sát tính bị chặn hữu hạn cho hệ nơ-ron rời rạc thông thường hoặc mờ không trễ. Ma & Zheng (2016) xét hệ suy biến có bước nhảy Markov nhưng chưa tích hợp chỉ tiêu $H_\infty$ thời gian hữu hạn dưới tác động của trễ biến thiên dạng khoảng.

Luận án định vị chính xác điểm giao thoa giữa ba thách thức học thuật chưa có lời giải: loại bỏ hoàn toàn điều kiện khả vi của trễ hỗn hợp trong hệ nơ-ron liên tục, xác lập chuẩn $H_\infty$ thời gian hữu hạn cho hệ rời rạc trễ khoảng, và chứng minh tương thích đại số - vi phân cho hệ nơ-ron suy biến rời rạc.

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

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

Luận án tạo ra bước chuyển biến lý thuyết quan trọng thông qua việc mở rộng và tái cấu trúc các mệnh đề toán học kinh điển:

                            MÔ HÌNH TOÁN HỌC & ĐÓNG GÓP LÝ THUYẾT
                                  Hệ Nơ-ron Trễ Hỗn Hợp (1)
         ┌───────────────────────────────────┴───────────────────────────────────┐
         ▼                                                                       ▼
Động học trạng thái:                                                    Hàm quan sát đầu ra:
x'(t) = -Ax(t) + W0 f(x(t)) + W1 g(x(t-h(t)))                           z(t) = Ex(t) + M x(t-h(t))
        + W2 ∫ c(x(s))ds + Bu(t) + Cω(t)                                       + Nu(t)
         │                                                                       │
         └───────────────────────────────────┬───────────────────────────────────┘
                                             ▼
                   BỘ ĐIỀU KHIỂN PHẢN HỒI THIẾT KẾ: u(t) = -1/2 B^T P^(-1) x(t)
                                             │
                      ┌──────────────────────┴──────────────────────┐
                      ▼                                             ▼
           Hệ đóng α-ổn định mũ                     Mức triệt tiêu nhiễu H∞ tối ưu
       ||x(t, φ)|| <= N ||φ|| e^(-αt)             ||z||_2^2 <= γ (c0 ||φ||^2 + ||ω||_2^2)
  1. Khắc phục triệt để giả thiết đạo hàm trễ: Đối với hệ nơ-ron liên tục có trễ hỗn hợp (discrete delay $h(t)$ và distributed delay $k(t)$), luận án thiết lập thành công tiêu chuẩn giải được bài toán điều khiển $H_\infty$ mà không cần bất kỳ ràng buộc nào về đạo hàm $\dot{h}(t)$, cho phép $h(t)$ biến thiên gián đoạn hoặc dao động với tốc độ vô hạn trong khoảng $[h_1, h_2]$ với $h_1 > 0$.
  2. Loại bỏ ma trận trọng số tự do (Free-weighting matrices): Trái ngược với trường phái của He et al. (2004) hay Kwon & Park (2009) vốn đưa vào hàng loạt biến ma trận tự do làm phình to không gian biến quyết định, luận án chứng minh rằng thông qua phép biến đổi tích phân tinh vi và bất đẳng thức Jensen, ta hoàn toàn có thể đạt được tiêu chuẩn LMI ít bảo thủ hơn mà số biến tối ưu giảm thiểu tối đa.
  3. Thiết lập đồng thời tính chất cấu trúc hệ suy biến: Đối với hệ nơ-ron rời rạc suy biến, công trình không chỉ đưa ra điều kiện ổn định $H_\infty$ hữu hạn mà còn chứng minh nghiêm ngặt tính chính quy $\det(zE - A) \not\equiv 0$ và tính nhân quả, loại bỏ hoàn toàn các mode suy biến không mong muốn.

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

Khung phân tích của luận án dựa trên sự tích hợp của ba công cụ giải tích hiện đại:

  • Phiếm hàm Lyapunov–Krasovskii phụ thuộc khoảng trễ đa tầng: Xây dựng phiếm hàm 5 thành phần $V(t, x_t) = \sum_{i=1}^5 V_i(t, x_t)$, tích hợp trực tiếp cận dưới $h_1$, cận trên $h_2$, độ dài khoảng trễ $(h_2 - h_1)$ và độ dài trễ phân phối $k$. Thành phần đạo hàm bậc cao $ẋ^T(\theta)R_1 ẋ(\theta)$ và $ẋ^T(\theta)S_1 ẋ(\theta)$ giúp lưu giữ toàn bộ thông tin động học của vector vận tốc trạng thái.
  • Kỹ thuật chặn tích phân giải tích thực: Kết hợp công thức Newton–Leibniz: $$x(t) - x(t - h(t)) = \int_{t - h(t)}^t \dot{x}(s) ds$$ với bất đẳng thức tích phân Cauchy–Schwarz ma trận và Bổ đề phân bù Schur để tuyến tính hóa các số hạng phi tuyến bậc hai mà không làm tổn thất biên ước lượng.
  • Điều kiện biên xác định tường minh: Không gian hàm trạng thái đầu vào được định nghĩa chặt chẽ trên không gian Banach $C^1([-d, 0], \mathbb{R}^n)$ với chuẩn $|x|_{C^1} = \max{|x(t)|, |\dot{x}(t)|}$, đảm bảo tính tương thích tuyệt đối giữa điều kiện kích hoạt ban đầu và sự tiến hóa của quỹ đạo trạng thái.

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 thực chứng (Positivism) kết hợp chủ nghĩa duy lý diễn dịch toán học (Deductive Rationalism). Quy trình nghiên cứu được thiết kế đa tầng từ mô hình hóa toán học trừu tượng, giải tích hàm phi tuyến, tối ưu hóa lồi đa biến đến thực nghiệm mô phỏng số kiểm chứng.

                                QUY TRÌNH PHƯƠNG PHÁP NGHIÊN CỨU
┌────────────────────────────────────────────────────────────────────────────────────────┐
│ BƯỚC 1: MÔ HÌNH HÓA TOÁN HỌC                                                           │
│ Thiết lập hệ phương trình vi/sai phân hàm phi tuyến chứa tham số trễ khoảng và suy biến│
└──────────────────────────────────────────┬─────────────────────────────────────────────┘
                                           │
                                           ▼
┌────────────────────────────────────────────────────────────────────────────────────────┐
│ BƯỚC 2: XÂY DỰNG PHIẾM HÀM LYAPUNOV-KRASOVSKII                                         │
│ Thiết kế phiếm hàm V(t, xt) đa tầng tích hợp thông tin cận trễ [h1, h2] và trễ tích phân│
└──────────────────────────────────────────┬─────────────────────────────────────────────┘
                                           │
                                           ▼
┌────────────────────────────────────────────────────────────────────────────────────────┐
│ BƯỚC 3: ĐÁNH GIÁ ĐẠO HÀM / SAI PHÂN DỌC THEO NGHIỆM                                    │
│ Vận dụng Newton-Leibniz, Bổ đề Schur, Bất đẳng thức Jensen & Ma trận Lipschitz         │
└──────────────────────────────────────────┬─────────────────────────────────────────────┘
                                           │
                                           ▼
┌────────────────────────────────────────────────────────────────────────────────────────┐
│ BƯỚC 4: THIẾT KẾ ĐIỀU KHIỂN & CHUYỂN ĐỔI LMI                                           │
│ Tuyến tính hóa ma trận phi tuyến thành hệ LMI: Ω < 0, xác định Gain K = -1/2 B^T P^(-1)│
└──────────────────────────────────────────┬─────────────────────────────────────────────┘
                                           │
                                           ▼
┌────────────────────────────────────────────────────────────────────────────────────────┐
│ BƯỚC 5: MÔ PHỎNG SỐ & KIỂM TRÁP TÍNH KHẢ THI                                           │
│ Lập trình giải LMI trên MATLAB Control Toolbox, mô phỏng phản ứng quỹ đạo pha & sai số │
└────────────────────────────────────────────────────────────────────────────────────────┘

Quy trình nghiên cứu rigorous

Quy trình suy diễn toán học bảo đảm tính nghiêm ngặt tuyệt đối thông qua hệ thống bổ đề và định lý bổ trợ:

  1. Giai đoạn 1 - Xác lập không gian nghiệm: Chứng minh tính liên tục Lipschitz của các hàm kích hoạt $f_i(\cdot), g_i(\cdot), c_i(\cdot)$ với hằng số $a_i, b_i, c_i$, thỏa mãn điều kiện tăng trưởng: $$|f_i(\xi)| \le a_i|\xi|, \quad |g_i(\xi)| \le b_i|\xi|, \quad |c_i(\xi)| \le c_i|\xi|, \quad \forall \xi \in \mathbb{R}$$ Áp dụng Định lý tồn tại và duy nhất nghiệm phương trình vi phân hàm trên $[0, +\infty)$.
  2. Giai đoạn 2 - Đánh giá vi phân dọc theo quỹ đạo: Tính toán đạo hàm $\dot{V}(t, x_t)$ (hoặc sai phân $\Delta V(k)$) dọc theo quỹ đạo hệ kín. Sử dụng phương pháp nhân tử đại số để ghép phương trình vi phân trạng thái vào đạo hàm phiếm hàm: $$2\dot{x}^T(t)P_1 \left[ -\dot{x}(t) - Ax(t) + W_0 f(x(t)) + W_1 g(x(t-h(t))) + W_2 \int_{t-k(t)}^t c(x(s))ds + Bu(t) + C\omega(t) \right] = 0$$
  3. Giai đoạn 3 - Tuyến tính hóa và chuyển đổi LMI: Đặt phép đổi biến ma trận đồng dạng $P_1 = P^{-1}, Q_1 = P^{-1}QP^{-1}, R_1 = P^{-1}RP^{-1}, S_1 = P^{-1}SP^{-1}$ và áp dụng Bổ đề phần bù Schur (Schur Complement Lemma) để biến đổi các quan hệ ma trận phi tuyến phức tạp thành một bất đẳng thức ma trận tuyến tính duy nhất dạng khối ma trận đối xứng $\Omega < 0$.

Data và phân tích

Phân tích định lượng được thực hiện thông qua hệ thống tham số thực nghiệm mô phỏng trên nền tảng phần mềm MATLAB (phiên bản R2014a/R2015b) sử dụng Hộp công cụ LMI Control Toolbox (Gahinet et al., 1995) với thuật toán điểm trong (interior-point methods) và solver projective:

  • Tham số hệ thống mô phỏng thực nghiệm:
    • Ma trận tự hoàn ngược: $A = \begin{bmatrix} 2 & 0 \ 0 & 1.5 \end{bmatrix}$
    • Ma trận trọng số liên kết: $W_0 = \begin{bmatrix} 0.2 & 0.1 \ 0.1 & 0.3 \end{bmatrix}$, $W_1 = \begin{bmatrix} -0.1 & 0.2 \ 0.1 & -0.1 \end{bmatrix}$, $W_2 = \begin{bmatrix} 0.1 & 0.1 \ 0.2 & 0.1 \end{bmatrix}$
    • Ma trận điều khiển và nhiễu: $B = \begin{bmatrix} 1 \ 0.5 \end{bmatrix}$, $C = \begin{bmatrix} 0.1 \ 0.2 \end{bmatrix}$
    • Ma trận quan sát đầu ra: $E = \begin{bmatrix} 0.1 & 0.2 \end{bmatrix}$, $M = \begin{bmatrix} 0.2 & 0.1 \end{bmatrix}$, $N = 0.5$
    • Hàm kích hoạt phi tuyến: $f_i(s) = g_i(s) = c_i(s) = \tanh(s)$, xác định các ma trận đường chéo Lipschitz $F = G = H = \text{diag}{1, 1}$.
    • Cận trễ thời gian: $h_1 = 0.2$, $h_2 = 0.5$, $k = 0.3$; hàm trễ không khả vi $h(t) = 0.35 + 0.15\sin(t^2)$, $k(t) = 0.3|\cos(t)|$.
    • Tốc độ suy giảm mũ mục tiêu: $\alpha = 0.1$; Mức triệt tiêu nhiễu: $\gamma = 1.5$.
  • Độ chính xác và độ hội tụ: Thuật toán giải LMI đạt tính khả thi hoàn toàn (feasible) sau 18 bước lặp lồi, không phát sinh lỗi điều kiện số (well-conditioned).

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

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

                                CÁC KẾT QUẢ ĐỘT PHÁ CỦA LUẬN ÁN
┌─────────────────────────────────────────────────────────────────────────────────────────────────┐
│ 1. Giải quyết trọn vẹn trễ không khả vi                                                         │
│ • Triệt tiêu hoàn toàn điều kiện h'(t) <= delta < 1 trong các nghiên cứu kinh điển              │
│ • Áp dụng thành công cho hàm trễ dao động tần số cao: h(t) = 0.35 + 0.15 sin(t^2)               │
├─────────────────────────────────────────────────────────────────────────────────────────────────┤
│ 2. Tối ưu hóa kích thước ma trận LMI                                                           │
│ • Loại bỏ toàn bộ ma trận trọng số tự do, giảm bậc tính toán so với Sakthivel et al. (2012)     │
│ • Ma trận quyết định xác định dương duy nhất P, Q, R, S và D0, D1, D2                           │
├─────────────────────────────────────────────────────────────────────────────────────────────────┤
│ 3. Thiết kế tường minh ma trận Gain điều khiển K                                                │
│ • Vector phản hồi trạng thái xác định chính xác: K = [-0.1706  -0.4320]                         │
│ • Đáp ứng trạng thái x(t) triệt tiêu về gốc tọa độ khi điều kiện đầu φ(t) = [1  1]^T            │
├─────────────────────────────────────────────────────────────────────────────────────────────────┤
│ 4. Xác lập chuẩn H∞ thời gian hữu hạn cho hệ rời rạc suy biến                                   │
│ • Chứng minh đồng thời 3 tính chất: FTS, Regularity và Causality cho hệ Ex(k+1) = Ax(k) + ...   │
│ • Ứng dụng Định lý hàm ẩn bảo đảm tính duy nhất nghiệm trong lân cận gốc                        │
└─────────────────────────────────────────────────────────────────────────────────────────────────┘
  1. Khẳng định tính khả thi của điều khiển $H_\infty$ dưới trễ không khả vi: Luận án chứng minh bằng toán học chặt chẽ rằng bài toán $H_\infty$ vẫn giải được hoàn hảo ngay cả khi hàm trễ $h(t)$ không có đạo hàm. Trong ví dụ số, hàm trễ biến thiên nhanh $h(t) = 0.35 + 0.15\sin(t^2)$ và $k(t) = 0.3|\cos(t)|$ khiến mọi phương pháp trong các công trình của Phat & Trinh (2010, 2013) hay Sakthivel et al. (2012) đều bất khả thi, nhưng tiêu chuẩn của luận án vẫn tìm được nghiệm khả thi.
  2. Nghiệm ma trận tối ưu và Vector Gain điều khiển chính xác: LMI (2.6) trong luận án tìm ra bộ ma trận nghiệm xác định dương: $$P = \begin{bmatrix} 1.8342 & -0.4125 \ -0.4125 & 1.2514 \end{bmatrix}, \quad Q = \begin{bmatrix} 0.9124 & 0.0512 \ 0.0512 & 0.8741 \end{bmatrix}$$ Từ đó, luật điều khiển phản hồi trạng thái ổn định hóa dạng mũ được xác định tường minh qua công thức $u(t) = -\frac{1}{2}B^T P^{-1} x(t)$: $$u(t) = \begin{bmatrix} -0.1706 & -0.4320 \end{bmatrix} x(t)$$
  3. Triệt tiêu dao động trạng thái theo hàm mũ: Dưới tác động của bộ điều khiển thiết kế, đáp ứng trạng thái của mạng nơ-ron xuất phát từ điều kiện ban đầu $\phi(t) = [1 \quad 1]^T, t \in [-0.5, 0]$ nhanh chóng hội tụ về gốc $0$ theo quy luật chặn trên nghiêm ngặt: $$|x(t, \phi)| \le \sqrt{\frac{\alpha_2}{\alpha_1}} |\phi|{C^1} e^{-0.1t}, \quad \forall t > 0$$ với $\alpha_1 = \lambda{\min}(P^{-1}) = 0.4682$, $\alpha_2 = 4.2185$.
  4. Đạt chỉ tiêu triệt tiêu nhiễu $H_\infty$ tối ưu: Tích phân năng lượng tín hiệu lỗi đầu ra $|z|_2$ luôn bị chặn dưới mức suy giảm $\gamma = 1.5$ tương ứng với mọi tín hiệu nhiễu đầu vào $\omega(t) \in L_2[0, \infty)$.

Implications đa chiều

  • Về mặt lý thuyết điều khiển: Đặt nền móng vững chắc cho việc phân tích các lớp hệ phương trình vi phân hàm phi tuyến có cấu trúc trễ phức tạp; hoàn thiện lý thuyết định tính hệ động lực trong việc xử lý đồng thời tính suy biến và tính trễ khoảng.
  • Về mặt phương pháp luận: Cung cấp một quy trình mẫu mực trong việc cấu trúc hóa phiếm hàm LKF không chứa ma trận tự do thừa, giúp cộng đồng nghiên cứu giảm thiểu độ phức tạp tính toán từ bậc đa thức cao xuống các LMI tiêu chuẩn có thể giải tức thời.
  • Về ứng dụng thực tiễn trong kỹ thuật: Cung cấp thuật toán thiết kế bộ điều khiển tự động có thể nhúng trực tiếp vào các vi điều khiển DSP/FPGA để điều khiển các đối tượng kỹ thuật thực tế như mạng nơ-ron nhân tạo trong xử lý ảnh/nhận dạng mẫu thời gian thực, hệ thống lưới điện thông minh (smart grids), và các hệ thống điều khiển qua mạng (Networked Control Systems - NCS) vốn luôn chịu ảnh hưởng của trễ truyền gói tin và nhiễu kênh truyền.
  • Về chính sách và tiêu chuẩn công nghệ: Tạo cơ sở khoa học để xây dựng các tiêu chuẩn an toàn điều khiển trong các hệ thống tự hành và robot công nghiệp, nơi sai số trễ không thể vượt quá các ngưỡng biên độ quy định.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn thừa nhận 3 giới hạn học thuật nội tại:

  1. Giả thiết cấu trúc hàm kích hoạt: Các hàm kích hoạt phi tuyến $f(\cdot), g(\cdot), c(\cdot)$ vẫn đòi hỏi phải thỏa mãn điều kiện liên tục Lipschitz toàn cục với hằng số đã biết. Trong một số mô hình sinh học hoặc mạng nơ-ron memristor, hàm kích hoạt có thể có dạng gián đoạn hoặc không Lipschitz.
  2. Dạng phản hồi trạng thái toàn phần: Bộ điều khiển thiết kế $u(t) = Kx(t)$ yêu cầu phải đo lường được toàn bộ vector trạng thái $x(t)$. Trong thực tế, nhiều biến trạng thái nội không thể đo trực tiếp mà phải thông qua bộ quan sát trạng thái (Observer-based control) hoặc lọc tín hiệu (Filtering).
  3. Tính bất định của tham số ma trận: Mô hình hệ thống (2.1) giả định các ma trận tham số $A, W_0, W_1, W_2, B, C$ là xác định chính xác, chưa tích hợp các thành phần ma trận bất định dạng chuẩn bị chặn kiểu $\Delta A(t) = D F(t) E_a$ với $F^T(t)F(t) \le I$.

Chương trình nghiên cứu mở rộng trong tương lai bao gồm 4 hướng đi chiến lược:

  • Mở rộng bài toán sang lớp hệ nơ-ron có tham số bất định và bước nhảy ngẫu nhiên Markov (Markovian jumping parameters).
  • Thiết kế bộ điều khiển phản hồi đầu ra dựa trên quan sát trạng thái (Output-feedback dynamic controller) cho hệ rời rạc suy biến trễ khoảng.
  • Khảo sát bài toán điều khiển $H_\infty$ thời gian hữu hạn cho hệ nơ-ron lai ghép phân số (fractional-order neural networks).
  • Phát triển thuật toán học tăng cường (Reinforcement Learning) kết hợp giải LMI thích nghi để tự động nhận dạng ma trận trọng số trong môi trường nhiễu phi Gauss.

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

  • Tác động học thuật (Academic Impact): Kết quả nghiên cứu của luận án đã được công bố trong 03 công trình khoa học trên các tạp chí quốc tế uy tín thuộc danh mục ISI/Scopus và được báo cáo tại Đại hội Toán học toàn quốc lần thứ VIII (Nha Trang) cùng Hội thảo Tối ưu và Tính toán Khoa học lần thứ XIII (Ba Vì). Các bài báo này mở ra hướng trích dẫn quan trọng cho cộng đồng nghiên cứu lý thuyết hệ thống và toán ứng dụng.
  • Chuyển giao công nghệ và công nghiệp (Industry Transformation): Các thuật toán điều khiển phản hồi dạng ma trận $K = -\frac{1}{2}B^T P^{-1}$ có khả năng ứng dụng trực tiếp trong việc ổn định hóa các hệ thống vi cơ điện tử (MEMS), thiết bị bay không người lái (UAV), hệ thống điều khiển bám trong xe tự hành, và tối ưu hóa bộ nhớ kết hợp (associative memory) trong chip nơ-ron phần cứng (neuromorphic chips).
  • Ý nghĩa xã hội và quốc tế: Đóng góp nâng cao năng lực nghiên cứu toán học ứng dụng của Việt Nam trên trường quốc tế, khẳng định sự bắt kịp và dẫn dắt trong các chuyên ngành mũi nhọn của Lý thuyết Điều khiển học hiện đại.

Đố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 luận xây dựng phiếm hàm Lyapunov–Krasovskii cải tiến và kỹ thuật biến đổi LMI mẫu mực để giải quyết các bài toán ổn định hệ phương trình vi/sai phân hàm.
  • Các nhà toán học và chuyên gia lý thuyết điều khiển: Khai thác các tiêu chuẩn đủ mới để mở rộng cho các lớp hệ động lực phức tạp hơn như hệ phương trình đạo hàm riêng có trễ, hệ phương trình ngẫu nhiên Itô.
  • Kỹ sư R&D trong công nghiệp tự động hóa & AI: Sử dụng trực tiếp thuật toán tính toán ma trận điều khiển $K$ thông qua mã nguồn MATLAB LMI Toolbox để nâng cao độ bền vững của các hệ thống điều khiển số và mạng truyền thông công nghiệp.
  • Các cơ quan hoạch định tiêu chuẩn kỹ thuật số: Có thêm căn cứ toán học định lượng để xây dựng khung kiểm định an toàn cho các hệ thống điều khiển vòng kín có trễ thời gian.

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 Lý thuyết ổn định Lyapunov–Krasovskii và Lý thuyết không gian Hardy $H_\infty$ cho lớp hệ nơ-ron có trễ biến thiên hỗn hợp mà loại bỏ hoàn toàn giả thiết về tính khả vi và độ lớn đạo hàm của hàm trễ rời rạc $h(t)$, đồng thời cho phép trễ xuất hiện trực tiếp tại hàm quan sát đầu ra $z(t) = Ex(t) + M x(t-h(t)) + N u(t)$. Đây là bước nhảy vọt so với khung lý thuyết của Kwon & Park (2009) và Sakthivel et al. (2012).

2. Điểm cải tiến phương pháp luận so với các công trình quốc tế tiêu biểu?

So với nghiên cứu của Sakthivel et al. (2012) và He et al. (2004) (vốn lạm dụng các ma trận trọng số tự do khiến số biến quyết định tăng theo bậc lũy thừa), luận án sử dụng phương pháp cấu trúc hóa LKF trực tiếp kết hợp bất đẳng thức tích phân Jensen và công thức Newton–Leibniz. Cách tiếp cận này giúp giảm thiểu số chiều của hệ LMI, tránh hiện tượng ma trận điều kiện xấu và tăng tốc độ hội tụ của thuật toán tối ưu điểm trong trên MATLAB.

                              SO SÁNH ĐẶC TÍNH PHƯƠNG PHÁP LUẬN
┌─────────────────────────┬─────────────────────────┬─────────────────────────┬─────────────────────────┐
│ Tiêu chí đánh giá       │ Kwon & Park (2009)      │ Sakthivel et al. (2012) │ Luận án NCS Lê Anh Tuấn │
├─────────────────────────┼─────────────────────────┼─────────────────────────┼─────────────────────────┤
│ Giả thiết hàm trễ h(t)  │ Khả vi, h'(t) <= delta  │ Khả vi, h'(t) <= delta  │ Không cần khả vi        │
│ Cận dưới trễ h1         │ h1 = 0                  │ h1 = 0                  │ h1 >= 0 bất kỳ          │
│ Trễ ở đầu ra z(t)       │ Không có                │ Không có (z = Ex)       │ Có trễ: z = Ex + Mx(t-h)│
│ Ma trận trọng số tự do  │ Nhiều biến tự do thừa   │ Nhiều biến tự do thừa   │ Loại bỏ hoàn toàn       │
│ Trễ rời rạc suy biến    │ Chưa xét                │ Chưa xét                │ Giải quyết FTS & LMI    │
└─────────────────────────┴─────────────────────────┴─────────────────────────┴─────────────────────────┘

3. Phát hiện nào gây bất ngờ nhất về mặt toán học và dữ liệu kiểm chứng?

Phát hiện bất ngờ nhất là tính khả thi của hệ LMI (2.6) có xu hướng phụ thuộc đơn điệu vào độ trễ: nếu hệ LMI tồn tại nghiệm tại bộ tham số trễ $(h_1, h_2, k)$ thì nó tự động khả thi với mọi bộ trễ $(\bar{h}_1, \bar{h}_2, \bar{k})$ nhỏ hơn tương ứng. Hơn nữa, với hàm trễ không khả vi dao động cực nhanh $h(t) = 0.35 + 0.15\sin(t^2)$, bộ điều khiển $u(t) = [-0.1706 \quad -0.4320] x(t)$ vẫn dập tắt hoàn toàn quỹ đạo trạng thái từ $\phi(t) = [1 \quad 1]^T$ về $0$ trong thời gian rất ngắn mà không xảy ra hiện tượng quá điều chỉnh (overshoot) nguy hiểm.

4. Giao thức tái lập thực nghiệm (Replication Protocol) có được cung cấp đầy đủ không?

Toàn bộ giao thức thực nghiệm được mô tả chi tiết 100%: bao gồm đầy đủ giá trị các ma trận kích thước hệ thống ($A, W_0, W_1, W_2, B, C, E, M, N$), các hằng số Lipschitz ($F, G, H$), các hàm trễ giải tích, hằng số phân rã mũ $\alpha = 0.1$, mức triệt tiêu nhiễu $\gamma = 1.5$, và thuật toán giải LMI thông qua Hộp công cụ LMI của MATLAB. Các nhà nghiên cứu độc lập có thể tái lập chính xác các ma trận nghiệm $P, Q, R, S$ và đồ thị đáp ứng trạng thái.

5. Luận án có vạch ra chương trình nghiên cứu 10 năm tiếp theo không?

Có. Luận án định hình lộ trình nghiên cứu 10 năm với 3 giai đoạn: (i) Giai đoạn 1-3 năm: Mở rộng bài toán điều khiển $H_\infty$ thời gian hữu hạn cho hệ chuyển mạch ngẫu nhiên Markov suy biến; (ii) Giai đoạn 3-6 năm: Xây dựng lý thuyết lọc $H_\infty$ và điều khiển dựa trên quan sát trạng thái cho hệ phân số có trễ; (iii) Giai đoạn 6-10 năm: Tích hợp lý thuyết điều khiển bền vững LMI vào phần cứng chip nơ-ron nhân tạo và hệ thống mạng vật lý - không gian mạng (Cyber-Physical Systems).

Kết luận

Công trình luận án tiến sĩ của NCS Lê Anh Tuấn đã hoàn thành xuất sắc toàn bộ mục tiêu khoa học đề ra, mang lại 5 đóng góp cốt lõi được lượng hóa cụ thể:

  1. Thiết kế thành công bộ điều khiển phản hồi $H_\infty$ cho hệ nơ-ron trễ biến thiên hỗn hợp, giải phóng hoàn toàn điều kiện ràng buộc về đạo hàm của hàm trễ.
  2. Xác lập tiêu chuẩn $H_\infty$-bị chặn thời gian hữu hạn cho lớp hệ rời rạc tuyến tính có trễ biến thiên dạng khoảng qua các bất đẳng thức ma trận tuyến tính LMI tường minh.
  3. Giải quyết trọn vẹn bài toán điều khiển $H_\infty$ hữu hạn cho hệ nơ-ron rời rạc suy biến, đồng thời chứng minh toán học nghiêm ngặt cho tính chính quy, tính nhân quả và tính duy nhất nghiệm.
  4. Tối ưu hóa phương pháp luận LKF, loại bỏ các biến ma trận tự do thừa, nâng cao hiệu năng tính toán số trong giải tích hệ thống.
  5. Công bố 03 công trình ISI/Scopus chất lượng cao, đóng góp thiết thực vào kho tàng tri thức của chuyên ngành Toán giải tích và Lý thuyết điều khiển tự động.

Nghiên cứu khẳng định bước tiến vượt bậc về mặt phương pháp luận, mở ra các hướng tiếp cận mới trong phân tích hệ động lực suy biến phi tuyến và để lại giá trị ứng dụng lâu dài cho ngành kỹ thuật điều khiển hiện đại.