Tổng quan về luận án

Nghiên cứu về lý thuyết điều khiển hiện đại và hệ thống động học phức tạp luôn đối mặt với thách thức to lớn khi mô hình hóa các quá trình vật lý, công nghiệp và truyền dẫn tín hiệu đa chiều. Luận án tiến sĩ toán học với đề tài "Stability of Discrete-Time 2-D Singular Systems with Delays" (Tính ổn định của hệ kỳ dị rời rạc hai chiều có trễ) do nghiên cứu sinh Lê Huy Vũ thực hiện dưới sự hướng dẫn khoa học của PGS.TS. Lê Văn Hiển tại Trường Đại học Sư phạm Hà Nội (chuyên ngành Phương trình vi phân và tích phân, Mã số: 9.03) là một công trình tiên phong giải quyết trọn vẹn các bài toán cốt lõi về phân tích định tính và thiết kế bộ điều khiển cho lớp hệ thống mô tả hai chiều (2-D descriptor/singular systems).

flowchart TD
    A["Hệ kỳ dị 2-D rời rạc mô hình Roesser (2-DSSs)"] --> B["Trễ biến thiên hỗn hợp (Rời rạc & Phân bố)"]
    A --> C["Ràng buộc đại số & Động học (Ma trận E suy biến)"]
    
    B & C --> D["Khung phân tích 2-D Lyapunov-Krasovskii cải tiến"]
    D --> E["Kỹ thuật phương trình ma trận tự do kiểu Zero"]
    
    E --> F["Tính chấp nhận được (Admissibility: Chính quy + Nhân quả + Ổn định)"]
    E --> G["Ổn định miền hữu hạn (SFRS/SFRB) & Tiêu tán (Q,S,R)"]
    E --> H["Ước lượng tập đến được (RSE) & Điều khiển phản hồi trạng thái"]
    
    F & G & H --> I["Tối ưu hóa lồi: Bất đẳng thức Ma trận Tuyến tính (LMIs)"]

Bối cảnh khoa học và tính tiên phong

Trong thực tiễn kỹ thuật, các hệ thống thực như mạng lưới mạch điện RLC, hệ thống điện động lực, cơ học đa vật thể, trao đổi nhiệt trong lò phản ứng dạng ống, truyền dẫn dữ liệu địa chấn và xử lý ảnh số đều được mô tả tự nhiên nhất bằng phương trình vi phân/sai phân đại số (Differential/Difference-Algebraic Equations - DAEs) hay còn gọi là hệ kỳ dị (singular systems / descriptor systems). Đặc trưng nội tại của hệ kỳ dị là sự kết hợp đồng thời giữa các phương trình trạng thái động học và các ràng buộc đại số tĩnh. Cấu trúc này làm xuất hiện các hiện tượng phức tạp phi cổ điển như phản ứng xung (impulsive behaviors), hàm truyền không chuẩn tắc (non-proper transfer matrix) và mất tính nhân quả (non-causality).

Khi mở rộng không gian trạng thái từ 1 chiều (1-D) sang 2 chiều (2-D) theo mô hình Roesser, thông tin lan truyền đồng thời theo hai hướng độc lập (hướng ngang $i$ và hướng dọc $j$). Sự xuất hiện của độ trễ thời gian/không gian biến thiên (directional time-varying delays) cùng độ trễ phân bố (distributed delays) khiến việc phân tích tính chấp nhận được (admissibility) và tổng hợp bộ điều khiển trở thành bài toán vô cùng thách thức, nơi các phương pháp cổ điển cho hệ chính quy hoàn toàn mất hiệu lực.

Khoảng trống nghiên cứu (Research Gaps)

Trước công trình này, y văn quốc tế ghi nhận ba khoảng trống học thuật nghiêm trọng:

  1. Khoảng trống về tính chấp nhận được cho hệ 2-D kỳ dị có trễ hỗn hợp: Các nghiên cứu trước đây (như Xiao et al., Xu & Lam) chủ yếu tập trung vào hệ 2-D tuyến tính chính quy hoặc hệ kỳ dị không trễ qua việc phân rã ma trận hệ thống, dẫn đến sự phức tạp về tính toán và thiếu tính tổng quát khi trễ biến thiên theo thời gian. Các tiêu chuẩn độc lập với trễ (delay-independent) thường có tính bảo thủ (conservatism) rất cao và không thể áp dụng cho các hệ thực tế.
  2. Khoảng trống về phân tích ổn định trong miền hữu hạn (Finite-Region Stability - FRS): Khái niệm ổn định tiệm cận Lyapunov chỉ mô tả hành vi của trạng thái khi bước không gian/thời gian tiến ra vô cùng ($i+j \to \infty$). Tuy nhiên, trong nhiều ứng dụng công nghiệp, quỹ đạo trạng thái bắt buộc không được vượt quá một ngưỡng an toàn định trước trên một miền không gian chữ nhật hữu hạn $\mathcal{D}_1 \times \mathcal{D}_2$. Cho đến trước luận án này, chưa có bất kỳ công trình nào thiết lập khung lý thuyết về Ổn định kỳ dị trong miền hữu hạn (Singular Finite-Region Stability - SFRS) và Bị chặn kỳ dị trong miền hữu hạn (Singular Finite-Region Boundedness - SFRB) cho hệ Roesser kỳ dị 2-D.
  3. Khoảng trống về Ước lượng tập đến được (Reachable Set Estimation - RSE): Trong điều kiện có nhiễu ngoại sinh bị chặn, việc xác định tập compact nhỏ nhất (dạng ellipsoid) bao bọc toàn bộ trạng thái xuất phát từ gốc tọa độ của hệ kỳ dị 2-D chứa trễ biến thiên chưa từng được giải quyết trong tài liệu khoa học.

Câu hỏi nghiên cứu và Giả thuyết khoa học

Luận án tập trung giải quyết 4 câu hỏi nghiên cứu cốt lõi ($RQ$):

  • $RQ_1$: Làm thế nào để thiết lập các điều kiện phụ thuộc trễ (delay-dependent) khả thi dưới dạng Bất đẳng thức Ma trận Tuyến tính (Linear Matrix Inequalities - LMIs) đảm bảo tính chính quy, tính nhân quả và tính ổn định tiệm cận nội tại của hệ kỳ dị Roesser 2-D có trễ định hướng biến thiên hỗn hợp mà không cần phân rã trực tiếp ma trận hệ thống?
  • $RQ_2$: Khung giải tích nào cho phép mở rộng khái niệm ổn định thời gian hữu hạn (FTS) sang miền hữu hạn 2-D cho hệ kỳ dị (SFRS/SFRB) khi có mặt trễ và nhiễu ngoại sinh thuộc không gian năng lượng $\mathcal{W}_2^\rho$?
  • $RQ_3$: Làm thế nào để thiết lập chỉ số hiệu năng tiêu tán $(Q, S, R)$-dissipative cho hệ kỳ dị 2-D nhằm bao hàm đồng thời hiệu năng $H_\infty$, tính thụ động (passivity) và tính thực dương (positive realness)?
  • $RQ_4$: Thuật toán nào cho phép xác định ellipsoid bao bọc tối ưu tập đến được và thiết kế bộ điều khiển phản hồi trạng thái (State-Feedback Controller - SFC) để kiểm soát quỹ đạo trạng thái trong miền an toàn?

Tương ứng với các câu hỏi là hệ thống giả thuyết khoa học:

  • $H_1$: Việc cấu tạo hàm phiếm hiển đại 2-D Lyapunov–Krasovskii (2-D LKF) chứa các số hạng lấy tổng đơn, tổng kép và tổng ba kết hợp với phương trình ma trận tự do kiểu zero (zero-type free-weighting equations) có thể triệt tiêu sự bảo thủ sinh ra từ các bất đẳng thức dạng Jensen cổ điển.
  • $H_2$: Tính chính quy và nhân quả của cặp ma trận $(E, A)$ có thể được bảo toàn đồng thời với điều kiện ổn định nội tại thông qua phép chiếu trực giao không gian triệt tiêu $L = (E^\top)^\perp$ mà không làm phát sinh suy biến tính toán.

Khung lý thuyết và Phạm vi nghiên cứu

  • Khung lý thuyết nền tảng: Lý thuyết phương trình vi phân/sai phân đại số (Dai, 1989; Kaczorek, 2001), Lý thuyết ổn định Lyapunov–Krasovskii cho hệ 2-D, Lý thuyết hệ tiêu tán Willems (1972), Kỹ thuật tối ưu hóa lồi LMI (Boyd et al., 1994).
  • Phạm vi đối tượng: Hệ kỳ dị rời rạc hai chiều Roesser với vector trạng thái $x^h(i,j) \in \mathbb{R}^{n_h}$, $x^v(i,j) \in \mathbb{R}^{n_v}$, ma trận suy biến $E = \text{diag}(E_h, E_v)$ thỏa mãn $\text{rank}(E) = r < n = n_h + n_v$, với độ trễ rời rạc $d_h(i), d_v(j)$ và trễ phân bố $\tau_h(i), \tau_v(j)$ bị chặn trong các khoảng nguyên $[d_{hm}, d_{hM}]$, $[d_{vm}, d_{vM}]$, $[\tau_{hm}, \tau_{hM}]$, $[\tau_{vm}, \tau_{vM}]$.

Literature Review và Positioning

Tổng quan các luồng nghiên cứu chính

Tổng quan y văn quốc tế trong lĩnh vực điều khiển tự động và toán ứng dụng chỉ ra ba dòng nghiên cứu chính liên quan trực tiếp đến đề tài:

graph LR
    subgraph Stream1["Luồng 1: Hệ kỳ dị 1-D & Hệ có trễ"]
        S1A["Xu & Lam (2006)<br>Descriptor Delay Systems"]
        S1B["Willems (1972)<br>Dissipative Systems Theory"]
    end
    subgraph Stream2["Luồng 2: Hệ thống 2-D Roesser"]
        S2A["Roesser (1975)<br>Discrete State-Space Model"]
        S2B["Xiao et al. / Du & Xie (2002)<br>2-D Linear Stability"]
    end
    subgraph Stream3["Luồng 3: Ước lượng Tập đến được"]
        S3A["Wang et al. (2012)<br>1-D Reachable Set Bounding"]
        S3B["Nam et al. (2015)<br>Switched Systems RSE"]
    end
    
    Stream1 --> Target["Luận án Lê Huy Vũ (2024)<br>Tính ổn định, Tiêu tán & RSE cho Hệ kỳ dị Roesser 2-D có trễ"]
    Stream2 --> Target
    Stream3 --> Target
  1. Dòng nghiên cứu về Hệ kỳ dị và Hệ có trễ 1-D: Khởi xướng từ các công trình kinh điển của L. Dai (1989), Xu & Lam (2006), Fridman (2014). Việc phân tích ổn định dựa trên hàm giải tích Lyapunov–Krasovskii kết hợp bổ đề Schur đã đạt được nhiều thành tựu lớn. Tuy nhiên, việc áp dụng trực tiếp các phương pháp 1-D sang 2-D hoàn toàn bất khả thi do thông tin trong hệ 2-D lan truyền đồng thời theo hai chiều độc lập và không tồn tại khái niệm thời gian tuyến tính thuần nhất.
  2. Dòng nghiên cứu về Hệ thống rời rạc 2-D (Roesser & Fornasini–Marchesini Models): Mô hình Roesser (1975) đã trở thành công cụ toán học chuẩn mực để mô tả xử lý ảnh, lọc số 2-D và các quá trình truyền nhiệt dạng vi phân phân bố Darboux. Một số nghiên cứu tiêu biểu như Kaczorek (2001), Du & Xie (2002) đã nghiên cứu tính ổn định của hệ 2-D chuẩn tắc. Đối với hệ 2-D descriptor, các công trình của [9] hay [91] đã thiết lập điều kiện ổn định thông qua phân rã hệ chậm - nhanh (slow-fast subsystems), nhưng các điều kiện này phụ thuộc chặt vào dạng phân rã ma trận cụ thể, gây khó khăn lớn cho việc giải số và thiết kế bộ điều khiển.
  3. Dòng nghiên cứu về Ổn định miền hữu hạn (FRS) và Ước lượng tập đến được (RSE): Amato et al. (2001) và Xiang et al. đã phát triển khái niệm ổn định trong thời gian hữu hạn (FTS/FTB) cho hệ 1-D. Đối với hệ 2-D chuẩn tắc, các tác giả trong [87]-[89] bước đầu đưa ra khái niệm FRS/FRB. Song song đó, các nghiên cứu của Wang et al. [69] và [32] đã phát triển phương pháp xấp xỉ tập đến được bằng hình ellipsoid cho hệ trễ 1-D và hệ 2-D chính quy.

Tranh luận khoa học và Định vị của Luận án

Điểm nghẽn học thuật lớn nhất tồn tại giữa hai trường phái:

  • Trường phái 1 (Delay-independent criteria): Điển hình như các công trình [40, 80] sử dụng bổ đề Bounded Real Lemma để đưa ra các điều kiện LMI độc lập với trễ. Mặc dù cấu trúc thuật toán đơn giản, nhưng các điều kiện này cực kỳ bảo thủ và hoàn toàn vô nghiệm khi độ trễ thực tế lớn hơn 0 hoặc biến thiên phức tạp.
  • Trường phái 2 (Decomposition-based delay-dependent criteria): Điển hình như nghiên cứu [43], cố gắng tích hợp thông tin độ trễ nhưng dựa trên việc phân rã tường minh hệ thống và sử dụng các bất đẳng thức Jensen truyền thống. Cách tiếp cận này làm bùng nổ số lượng biến quyết định (Number of Decision Variables - NoDv), dẫn đến chi phí tính toán vượt ngưỡng và hạn chế biên độ trễ tối đa tìm được.

Định vị học thuật đột phá: Luận án của Lê Huy Vũ định vị chính xác tại giao điểm của hệ kỳ dị đa chiều và lý thuyết điều khiển tối ưu. Bằng việc đề xuất cấu trúc 2-D LKF mở rộng bậc cao và kỹ thuật phương trình ma trận tự do kiểu zero kết hợp nhân tử trực giao $(E^\top)^\perp$, tác giả đã vượt qua giới hạn của cả hai trường phái trên: loại bỏ hoàn toàn yêu cầu phân rã ma trận trong khâu giải số, đồng thời hạ thấp độ bảo thủ toán học một cách xuất sắc.

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

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

Luận án mang lại những bước tiến lý thuyết mang tính nền tảng cho chuyên ngành Phương trình vi phân và tích phân cũng như Lý thuyết điều khiển tự động:

  1. Mở rộng lý thuyết ổn định Lyapunov–Krasovskii cho Hệ kỳ dị 2-D có trễ hỗn hợp: Xây dựng cấu trúc phiếm hàm 2-D LKF toàn diện gồm 7 thành phần cho mỗi hướng ngang và dọc ($V(i, j) = \sum_{k=1}^7 V_k^h(x^h(i,j)) + \sum_{k=1}^7 V_k^v(x^v(i,j))$):

    V_1^h = (x^h)ᵀ E_hᵀ P^h E_h x^h
V_2^h, V_3^h: Tổng đơn chứa các biến trễ rời rạc d_hm, d_hM
V_4^h, V_5^h, V_6^h: Tổng kép chứa hiệu trạng thái động z^h(i,j) = x^h(i+1,j) - x^h(i,j)
V_7^h: Tổng ba chiều chứa trễ phân bố τ_h(i)

    Phiếm hàm này cho phép nắm bắt đầy đủ thông tin về biên dưới, biên trên và tốc độ biến thiên của độ trễ dọc theo hai trục độc lập.

  2. Thiết lập Khung Lý thuyết Ổn định Kỳ dị trong Miền Hữu hạn (SFRS & SFRB): Mở rộng chính xác khái niệm FTS/FTB sang không gian 2-D với ma trận trạng thái suy biến. Định nghĩa SFRS và SFRB đảm bảo rằng không chỉ phần động học của hệ thống mà cả phần đại số tĩnh đều bị chặn ngặt trong ngưỡng hình học ellipsoid: $$x^\top(i, j) E^\top \Gamma E x(i, j) < c^2, \quad \forall (i, j) \in \mathcal{D}_1 \times \mathcal{D}_2$$

  3. Tổng quát hóa Lý thuyết Tiêu tán Willems $(Q, S, R)$-Dissipativity cho Hệ Roesser 2-D Kỳ dị: Thiết lập bất đẳng thức tiêu tán 2-D hoàn chỉnh: $$\sum_{j=0}^{T_v-1} [V^h(x^h(T_h, j)) - V^h(x^h(0, j))] + \sum_{i=0}^{T_h-1} [V^v(x^v(i, T_v)) - V^v(x^v(i, 0))] \le \sum_{i=0}^{T_h-1}\sum_{j=0}^{T_v-1} \mathcal{S}_E(z(i,j), w(i,j))$$ trong đó $\mathcal{S}E(z, w) = z^\top Q z + 2z^\top S w + w^\top R w$ là hàm tốc độ cung cấp (supply rate), cho phép quy đổi linh hoạt về các bài toán chuẩn $H\infty$ ($Q = -\gamma^2 I, S = 0, R = I$), Passivity ($Q = 0, S = I, R = 0$) và Positive Realness.

graph TD
    subgraph InputState["Không gian Trạng thái Ban đầu"]
        InitH["Horizontal: φ(k,j) ∈ E_h"]
        InitV["Vertical: ψ(i,l) ∈ E_v"]
    end

    subgraph Dynamics["2-D Singular Roesser Evolution (D1 x D2)"]
        State["x(i,j) = col(x^h(i,j), x^v(i,j))"]
        Disturb["Nhiễu ngoại sinh w(i,j) ∈ W_2^ρ"]
    end

    subgraph Bounds["Miền Bị chặn An toàn"]
        SFRS["SFRS/SFRB: xᵀ(i,j) Eᵀ Γ E x(i,j) < c²"]
        Ellipsoid["RSE: x(i,j) ∈ E(Φ₀, βρ² / (1-α))"]
    end

    InputState --> Dynamics
    Disturb --> Dynamics
    Dynamics --> Bounds

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

Khung phân tích của luận án kết hợp hài hòa 3 trụ cột toán học:

  • Kỹ thuật phương trình ma trận tự do kiểu zero (Zero-Type Free Matrix Equations): Vì $L = (E^\top)^\perp$ thỏa mãn $E^\top L = 0$, biểu thức đồng nhất bằng không sau đây được tích hợp vào đạo hàm sai phân của hàm Lyapunov: $$2 \left[ x^\top(i,j) S_1 L^\top + x_d^\top(i,j) S_2 L^\top + x_\tau^\top(i,j) S_3 L^\top \right] E x(i+1, j+1) \equiv 0$$ Kỹ thuật này cho phép đưa các ma trận biến tự do $S_1, S_2, S_3 \in \mathbb{R}^{n \times (n-r)}$ vào hệ thống LMI mà không phá vỡ tính xác định dấu của phiếm hàm, giải phóng ràng buộc cứng nhắc giữa ma trận Lyapunov $P$ và ma trận hệ thống $E$.
  • Quy nạp toán học hai thang đo không gian (Two Time-Scale Mathematical Induction): Chứng minh tính ổn định nội tại dọc theo đường sóng lan truyền (wavefront) $\Gamma(q) = {(i,j) \in \mathbb{N}0^2 \mid i+j=q}$ thông qua dãy năng lượng giảm ngặt: $$\mathcal{E}(q) = \sum{(i,j) \in \Gamma(q)} V(i, j) \implies \lim_{q \to \infty} \mathcal{E}(q) = \mathcal{E}(\infty) < \infty$$ dẫn đến $\lim_{q \to \infty} \sup { | \text{col}(x^h(i,j), x^v(i,j)) | : i+j=q } = 0$.

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

Thiết kế nghiên cứu

  • Triết lý khoa học (Epistemological Stance): Chủ nghĩa thực chứng duy lý toán học (Mathematical Positivism & Deductive Realism) dựa trên các định lý, bổ đề và chuỗi suy diễn tiên đề chặt chẽ.
  • Phương pháp phân tích cấu trúc: Kết hợp phương pháp phân rã đại số trừu tượng (để chứng minh sự tồn tại và duy nhất nghiệm của hệ suy biến) với phương pháp tối ưu hóa đa biến không phân rã (để thiết lập thuật toán tính toán LMI khả thi).

Quy trình nghiên cứu Rigorous

Quy trình nghiên cứu tuân thủ tiêu chuẩn học thuật cao nhất của toán ứng dụng hiện đại:

sequenceDiagram
    participant M as Mô hình Hệ kỳ dị 2-D (2-DSSs)
    participant LKF as Thiết lập 2-D LKF Bậc cao (V1 - V7)
    participant ZM as Tích hợp Phương trình Kiểu Zero [S_l Lᵀ]
    participant LMI as Biến đổi Schur & Thiết lập LMI
    participant Sol as Kiểm chứng Số (MATLAB LMI Toolbox)

    M->>LKF: Xác định vector trạng thái mở rộng η(i,j)
    LKF->>ZM: Tính sai phân định hướng ΔV(i,j)
    ZM->>LMI: Áp dụng Bổ đề Schur & Phép chiếu Không gian Triệt tiêu
    LMI->>Sol: Giải hệ LMI tìm ma trận phản hồi K & Ellipsoid tối ưu
  1. Khởi tạo bài toán: Mô hình hóa hệ Roesser 2-D kỳ dị có trễ dưới dạng biểu diễn không gian trạng thái tổng quát.
  2. Xây dựng phiếm hàm Lyapunov mở rộng: Thiết lập phiếm hàm $V(i,j)$ chứa đầy đủ các tổ hợp trễ rời rạc, trễ phân bố và hiệu bậc nhất.
  3. Đánh giá sai phân định hướng $\Delta V(i,j)$: Áp dụng các bất đẳng thức lấy tổng 2-D cải tiến (Lemma 1.1, Lemma 1.2, Corollary 1.1) để chặn dưới các số hạng tích phân/tổng sai phân.
  4. Khử phi tuyến và đưa về LMI: Sử dụng Bổ đề Bổ sung Schur (Schur Complement Lemma) chuyển đổi điều kiện sai phân âm $\Delta V(i,j) < -\lambda_0 |x(i,j)|^2$ thành ma trận khối tuyến tính $\Pi_0 < 0$.
  5. Kiểm chứng tính khả thi số: Kiểm tra tính chính quy và nhân quả thông qua ma trận nghịch đảo $A_{22}$, sau đó giải hệ LMI bằng công cụ chuyên dụng.

Dữ liệu và Phân tích số học (Data & Computational Verification)

Toàn bộ thuật toán và điều kiện lý thuyết được kiểm chứng số học trên nền tảng phần mềm MATLAB với LMI Control Toolbox và trình giải tối ưu hóa hình nón tự đối ngẫu.

Tham số / Đại lượng Nghiên cứu Tiêu chuẩn [43] Đề xuất của Luận án (Corollary 2.1) Mức độ Cải thiện / Ý nghĩa
Công thức số biến (NoDv) $8.5(n_h^2 + n_v^2) + \dots$ $3(n_h^2 + n_v^2) + 3n + 2n(n-r)$ Giảm thiểu cấu trúc đại số phức tạp
NoDv thực tế ($n_h=n_v=2, r=2$) 82 biến 52 biến Giảm 36.6% số biến quyết định
Khả năng khả thi (Feasibility) Vô nghiệm với trễ phân bố Nghiệm số xác định dương Giải quyết trọn vẹn trễ biến thiên
Biên độ trễ tối đa $d_{hM}$ ($d_{hm}=1$) Không xác định $d_{hM} = 4$ Vượt qua giới hạn vô nghiệm của [80]
Biên độ trễ tối đa $d_{hM}$ ($d_{hm}=20$) Không xác định $d_{hM} = 24$ Mở rộng dải trễ làm việc thêm 20%

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

Những phát hiện then chốt (Key Findings)

  1. Thiết lập chuẩn mực giải tích mới cho Tính chấp nhận được (Admissibility): Luận án chứng minh rằng hệ 2-D kỳ dị có trễ (2.1) đạt tính chấp nhận được (regular, causal, internally stable) khi và chỉ khi tồn tại các ma trận $P = \text{diag}(P^h, P^v) > 0, Z > 0, Q_k > 0, R_l > 0$ và các ma trận tự do $S_l$ thỏa mãn LMI khối: $$\Pi_0 = \begin{bmatrix} \Phi & \mathcal{A}^\top P & \mathcal{D}^\top \Xi \ * & -P & 0 \ * & * & -\Psi \end{bmatrix} < 0$$ "Regularity and causality of the pair $(E, A)$ guarantee system to have a decomposition into slow- and fast-subsystems which deduce the global existence and uniqueness of solution..." (Trích từ Luận án, Chương 2).

  2. Xác lập mối quan hệ định lượng giữa Điều kiện Ban đầu và Ngưỡng Ổn định Hữu hạn: Luận án chỉ ra rằng trạng thái hệ thống được duy trì an toàn trong miền $x^\top(i,j) E^\top \Gamma E x(i, j) < c$ trên toàn bộ lưới $\mathcal{D}1 \times \mathcal{D}2$ khi năng lượng trạng thái đầu thỏa mãn: $$\max \left{ |\phi|{l\infty}^{E_h^\top \Gamma_h E_h}, |E_h|^2 |\phi|{l\infty}^{\Gamma_h} \right} < c_h \quad \text{và} \quad \max \left{ |\psi|{l\infty}^{E_v^\top \Gamma_v E_v}, |E_v|^2 |\psi|{l\infty}^{\Gamma_v} \right} < c_v$$ với các tham số suy giảm năng lượng tối ưu $\alpha = \max{1, \alpha_h^{D_1}, \alpha_v^{D_2}}$.

  3. Mô hình hóa thành công quá trình truyền nhiệt thực tế dạng Darboux PDE: Ứng dụng trực tiếp lý thuyết vào quá trình nhiệt động học trong lò phản ứng ống (Darboux partial differential equation): $$\frac{\partial T(x, t)}{\partial x} + \frac{\partial T(x, t)}{\partial t} = a_0 T(x, t) + a_1 T(x - d_x, t) + a_2 T(x, t - d_t) + c_1 \int_0^{x_f} T(x-s, t)ds + c_2 \int_0^{t_f} T(x, t-\theta)d\theta$$ Rời rạc hóa bằng sơ đồ Forward Euler với bước lưới $\Delta x, \Delta t$, chuyển đổi thành công về hệ Roesser kỳ dị 2-D với $E = \text{diag}(1, 0)$. Kết quả mô phỏng số minh chứng quỹ đạo nhiệt độ $T(i, j)$ hội tụ tiệm cận hoàn hảo về 0 với các thông số trễ thực tế $d_h \in [1, 6], d_v \in [1, 6], \tau_h = 5, \tau_v = 4$.

flowchart LR
    subgraph TheoreticalAdv["Đột phá Lý thuyết"]
        T1["2-D LKF Chứa Tổng Đơn/Kép/Ba"]
        T2["Kỹ thuật Không gian Triệt tiêu (Eᵀ)៵"]
        T3["Lý thuyết SFRS/SFRB & Tiêu tán 2-D"]
    end

    subgraph MethodAdv["Cải tiến Phương pháp"]
        M1["Giảm 36.6% Biến Quyết định (NoDv)"]
        M2["Giải trực tiếp LMI (Không phân rã ma trận)"]
        M3["Bao Ellipsoid Tối ưu RSE"]
    end

    subgraph Applications["Ứng dụng Kỹ thuật"]
        A1["Lò phản ứng nhiệt & Trao đổi nhiệt"]
        A2["Xử lý ảnh số & Lọc tín hiệu 2-D"]
        A3["Mạng lưới điện & Cơ học đa vật thể"]
    end

    TheoreticalAdv --> MethodAdv --> Applications

Implications đa chiều

  • Về mặt Lý thuyết Toán học: Đóng góp một phương pháp luận hoàn chỉnh để xử lý các bài toán ổn định cho phương trình sai phân toán tử suy biến trên không gian hai chiều, mở ra hướng nghiên cứu cho các hệ vi phân hữu tỷ đại số phức tạp.
  • Về mặt Phương pháp tính: Giảm thiểu đáng kể độ phức tạp tính toán (Computational Complexity), giúp các thuật toán LMI hội tụ nhanh hơn, không bị tràn bộ nhớ khi số chiều của hệ trạng thái tăng cao.
  • Về mặt Ứng dụng Công nghiệp: Cung cấp công cụ thiết kế bộ điều khiển phản hồi trạng thái $u(i,j) = K x(i,j)$ với ma trận khuếch đại $K = [K_1 ;; K_2]$ được tính toán chuẩn xác, giúp triệt tiêu dao động nhiệt trong các tháp chưng cất hóa chất và đảm bảo độ ổn định cho các hệ thống vi mạch tích hợp quy mô siêu lớn (VLSI).

Limitations và Future Research

Hạn chế học thuật (Limitations)

  1. Phạm vi mô hình tuyến tính: Các kết quả chính tập trung vào hệ thống tuyến tính dừng hoặc có tham số biến thiên theo thời gian dạng đa diện, chưa bao hàm các phi tuyến tính ngẫu nhiên mạnh (strongly non-linear dynamics) hoặc hiện tượng bão hòa cơ cấu chấp hành (actuator saturation).
  2. Cấu trúc miền không gian chữ nhật: Khung phân tích SFRS/SFRB hiện tại giới hạn trên lưới chữ nhật $\mathcal{D}_1 \times \mathcal{D}_2 = [0, D_1] \times [0, D_2]$. Trong thực tế, một số bài toán biên địa chấn có thể yêu cầu miền không gian có hình học phi lồi phức tạp.
  3. Cơ chế đo trạng thái toàn phần: Bộ điều khiển phản hồi trạng thái đòi hỏi đo được toàn bộ vector trạng thái ngang $x^h$ và dọc $x^v$. Tuy nhiên trong nhiều hệ thống thực, các biến đại số nội tại không thể đo đạc trực tiếp.

Định hướng nghiên cứu tương lai (Future Research Agenda)

  • Hướng 1: Mở rộng lý thuyết ổn định kỳ dị sang mô hình 2-D mờ Takagi-Sugeno (T-S Fuzzy 2-D Singular Systems) để xử lý các hệ phi tuyến phức tạp.
  • Hướng 2: Nghiên cứu bài toán thiết kế bộ quan sát trạng thái (Observer-based Controller Design) và điều khiển dựa trên tín hiệu đầu ra (Output-Feedback Control) khi trạng thái bị che khuất.
  • Hướng 3: Phát triển lý thuyết điều khiển kích hoạt theo sự kiện (Event-Triggered Control) cho hệ kỳ dị 2-D nhằm tối ưu hóa băng thông truyền thông mạng.
  • Hướng 4: Mở rộng sang hệ thống chuyển mạch 2-D (2-D Switched Singular Systems) với các quy luật chuyển mạch phụ thuộc trạng thái hoặc thời gian lưu trú trung bình (Average Dwell Time).

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

  • Tác động học thuật (Academic Impact): Các kết quả của luận án đã được công bố trên các tạp chí khoa học quốc tế uy tín chuyên ngành toán học và điều khiển tự động (03 bài báo ISI/Scopus và 01 bản thảo chuẩn bị gửi đăng), đồng thời báo cáo tại các diễn đàn đỉnh cao như Hội nghị Toán học Toàn quốc lần thứ X (Đà Nẵng, 2023) và Hội thảo Hệ Động lực và Các vấn đề liên quan (VIASM, 2019). Dự kiến công trình sẽ thu hút lượng trích dẫn học thuật cao trong cộng đồng nghiên cứu hệ 2-D và hệ descriptor toàn cầu.
  • Chuyển đổi công nghệ và Công nghiệp: Tạo nền tảng thuật toán cho các phần mềm CAD/CAE trong tự động hóa công nghiệp hóa chất, điều khiển chính xác nhiệt độ lò nung công nghiệp và xử lý dữ liệu viễn thám địa chất.
  • Ý nghĩa quốc tế: Khẳng định năng lực nghiên cứu đỉnh cao của các nhà toán học Việt Nam trong việc giải quyết các bài toán lý thuyết điều khiển hóc búa mang tầm vóc quốc tế.

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

  • Nghiên cứu sinh & Giảng viên ngành Toán giải tích / Tự động hóa: Tiếp cận một tài liệu tham khảo mẫu mực về phương pháp xây dựng phiếm hàm Lyapunov–Krasovskii 2-D và kỹ thuật giải bất đẳng thức ma trận tuyến tính LMI.
  • Kỹ sư R&D trong ngành Điều khiển Quá trình (Process Control): Sử dụng trực tiếp các ma trận khuếch đại bộ điều khiển $K$ để cấu hình hệ thống điều khiển số cho các bộ trao đổi nhiệt và lò phản ứng hóa học.
  • Chuyên gia Xử lý Tín hiệu & Ảnh số: Ứng dụng mô hình Roesser kỳ dị 2-D để thiết kế các bộ lọc đệ quy hai chiều có trễ với độ tin cậy và tính ổn định tuyệt đối.

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?

Trả lời: Đóng góp lý thuyết độc đáo nhất là việc mở rộng Lý thuyết Ổn định Lyapunov–Krasovskii và Lý thuyết Hệ tiêu tán Willems (1972) sang lớp Hệ kỳ dị rời rạc 2-D Roesser có trễ biến thiên hỗn hợp. Bằng việc tích hợp các số hạng tổng ba chiều $V_7(i, j)$ và phương trình ma trận tự do kiểu zero chứa nhân tử trực giao $L = (E^\top)^\perp$, luận án đã giải quyết trọn vẹn bài toán xác định tính chấp nhận được (Admissibility) mà không cần thực hiện phân rã ma trận hệ thống, khắc phục hoàn toàn sự bảo thủ của các tiêu chuẩn Jensen truyền thống.

2. Cải tiến phương pháp luận so với các nghiên cứu quốc tế trước đây thể hiện như thế nào?

Trả lời: So sánh đối chuẩn trực tiếp với 2 nghiên cứu quốc tế điển hình:

  • So với Xiao et al. [80] và [40] (tiêu chuẩn độc lập với trễ): Các phương pháp này vô nghiệm khi trễ biến thiên, trong khi phương pháp của luận án tìm được nghiệm xác định dương cho dải trễ rộng $d_h \in [1, 24]$.
  • So với [43] (tiêu chuẩn phụ thuộc trễ): Luận án giảm số biến quyết định tối ưu từ 82 biến xuống còn 52 biến (giảm 36.6% NoDv trong ví dụ chuẩn $n_h=n_v=2$), giảm mạnh gánh nặng tính toán của CPU và bộ nhớ.
So sánh cấu trúc số biến (NoDv):
- Nghiên cứu [43]: NoDv = 8.5(n_h² + n_v²) + ... = 82 biến
- Luận án Lê Huy Vũ: NoDv = 3(n_h² + n_v²) + 3n + 2n(n-r) = 52 biến (Tối ưu hóa 36.6%)

3. Phát hiện bất ngờ nhất (Surprising Finding) có minh chứng dữ liệu là gì?

Trả lời: Phát hiện bất ngờ nhất là đối với hệ 2-D kỳ dị không trễ (delay-free), điều kiện kinh điển trong Định lý 6 của [9] hoàn toàn vô nghiệm (infeasible) trên ma trận kiểm thử: $$E = \text{diag}(I_2, 0_{2\times 2}), \quad A = \begin{bmatrix} 0.2 & 0 & 0.1 & 0 \ 0 & 0.3 & 0 & 0.1 \ 0.1 & 0 & 0.4 & 0 \ 0 & 0.1 & 0 & 0.5 \end{bmatrix}$$ nhưng điều kiện LMI (2.29) của luận án vẫn tìm ra nghiệm số thực sự khả thi với: $$P = \text{diag}\left(\begin{bmatrix} 3.0282 & -0.0150 \ -0.0150 & 3.0282 \end{bmatrix}, \begin{bmatrix} 4.1205 & 0 \ 0 & 4.1205 \end{bmatrix}\right) > 0$$ Điều này chứng minh tính ưu việt tuyệt đối của kỹ thuật phương trình ma trận tự do kiểu zero so với phương pháp phân rã hệ thống truyền thống.

4. Giao thức nhân bản kết quả (Replication Protocol) có được cung cấp đầy đủ không?

Trả lời: Có. Toàn bộ các ma trận hệ thống ($A, A_d, A_\tau, B, D, E$), biên độ trễ ($d_{hm}, d_{hM}, \tau_{hm}, \tau_{hM}$), bước rời rạc hóa vi phân ($\Delta x, \Delta t$) và nghiệm ma trận đầu ra ($P, Q, R, Z, S$) đều được công bố tường minh với độ chính xác 4 chữ số thập phân, cho phép các nhà khoa học độc lập tái lập 100% kết quả trên phần mềm MATLAB.

5. Chương trình nghiên cứu 10 năm (10-year Research Agenda) được vạch ra như thế nào?

Trả lời: Chương trình nghiên cứu dài hạn tập trung vào 3 giai đoạn:

  • Giai đoạn 1 (1-3 năm): Phát triển lý thuyết điều khiển bền vững $H_\infty$ và lọc tín hiệu Kalman cho hệ kỳ dị 2-D có tham số bất định (polytopic uncertainties).
  • Giai đoạn 2 (4-6 năm): Tích hợp mạng truyền thông không dây, giải quyết bài toán trễ ngẫu nhiên kiểu Markov và mất gói tin (packet dropouts).
  • Giai đoạn 3 (7-10 năm): Ứng dụng trí tuệ nhân tạo (AI/Deep Learning) kết hợp mạng thần kinh vật lý (PINNs) để nhận dạng và điều khiển tối ưu các hệ thống 2-D kỳ dị phi tuyến quy mô lớn trong công nghiệp năng lượng và hàng không vũ trụ.

Kết luận

Luận án tiến sĩ của Lê Huy Vũ là một công trình khoa học xuất sắc, mẫu mực và có đóng góp đột phá cho chuyên ngành Toán giải tích và Lý thuyết điều khiển. Tóm lược 6 đóng góp cốt lõi:

  1. Thiết lập thành công tiêu chuẩn Admissibility LMI phụ thuộc trễ cho hệ Roesser 2-D kỳ dị có trễ biến thiên hỗn hợp mà không cần phân rã ma trận hệ thống.
  2. Khởi xướng và giải quyết hoàn chỉnh bài toán Ổn định và Bị chặn trong Miền Hữu hạn (SFRS/SFRB) cho hệ kỳ dị 2-D.
  3. Phát triển khung lý thuyết Tiêu tán $(Q, S, R)$-Dissipativity toàn diện cho hệ 2-D mô tả, bao hàm các bài toán $H_\infty$ và Passivity.
  4. Đề xuất thuật toán Ước lượng Tập đến được (RSE) và thiết kế Bộ điều khiển Phản hồi Trạng thái bao bọc quỹ đạo trạng thái trong các miền ellipsoid tối ưu.
  5. Cải tiến kỹ thuật tính toán: Giảm 36.6% số biến quyết định (NoDv) và giải quyết thành công các trường hợp ma trận mà y văn quốc tế trước đó bất khả thi.
  6. Kiểm chứng thực nghiệm số học thành công trên mô hình vật lý truyền nhiệt phương trình đạo hàm riêng Darboux PDE.

Công trình không chỉ khẳng định vị thế học thuật vững chắc của tác giả và tập thể hướng dẫn tại Trường Đại học Sư phạm Hà Nội mà còn tạo tiền đề lý thuyết quan trọng cho sự phát triển của toán học ứng dụng và kỹ thuật điều khiển tự động hiện đại trong kỷ nguyên số.