Tổng quan về luận án

Trong kỷ nguyên công nghiệp 4.0 và tự động hóa toàn diện, các quy trình công nghiệp hiện đại đều vận hành dưới dạng các hệ thống đa biến phức tạp (Multi-Input Multi-Output - MIMO). Tại các hệ thống này, sự tương tác chéo đa chiều (cross-coupling) giữa các biến quá trình cùng với sự hiện diện cố hữu của độ trễ thời gian (dead time) tạo nên rào cản nghiêm trọng đối với việc điều khiển độc lập từng vòng lặp. Theo các phân tích thực nghiệm trong lý thuyết điều khiển quá trình, "bộ điều khiển PID vẫn đang chiếm đa số (khoảng 90%) trong các bộ điều khiển công nghiệp". Tuy nhiên, cấu trúc PID bậc nguyên (Integer-Order PID - IOPID) kinh điển với ba tham số cố định ($K_c, \tau_I, \tau_D$) dần bộc lộ hạn chế khi xử lý các đối tượng phi tuyến, có quán tính nhiệt lớn hoặc tương tác động học mạnh. Luận án tiến sĩ kỹ thuật cơ khí mang tên "Nghiên cứu về mô hình hóa và điều khiển phân số cho các quá trình đa biến" của nghiên cứu sinh Võ Lâm Chương (Trường Đại học Sư phạm Kỹ thuật TP. Hồ Chí Minh, dưới sự hướng dẫn của PGS. Trương Nguyễn Luân Vũ và GS. MoonYong Lee) đã mở ra bước đột phá tiên phong khi tích hợp giải tích phân số (fractional calculus) vào cấu trúc điều khiển phân ly đa biến kết hợp bộ bù trễ Smith.

                  +-------------------------------------------------------------+
                  |                  HỆ QUÁ TRÌNH ĐA BIẾN (MIMO)                |
                  |     Tương tác chéo động học mạnh + Thời gian trễ lớn       |
                  +-------------------------------------------------------------+
                                                 |
                   +-----------------------------+-----------------------------+
                   |                                                           |
                   v                                                           v
+--------------------------------------+                   +---------------------------------------+
|  HỆ BẬC THẤP (2x2 TITO)              |                   |  HỆ BẬC CAO (3x3, 4x4)                |
|  - Tháp Vinante & Luyben             |                   |  - Tháp chưng cất Ogunnaike & Ray     |
|  - Phân tách dầu nặng (Heavy Oil)    |                   |  - Hệ thống điều hòa không khí HVAC   |
+--------------------------------------+                   +---------------------------------------+
                   |                                                           |
                   | (Phân ly đơn giản hóa + Bộ dự báo Smith)                 | (Phân ly đơn giản hóa + Smith Predictor)
                   | (PSO rút gọn hàm truyền phân số)                          | (PSO rút gọn hàm truyền phân số)
                   v                                                           v
+--------------------------------------+                   +---------------------------------------+
|  Cấu trúc Mô hình nội (IMC)          |                   |  Tối ưu bầy đàn đa mục tiêu (MOPSO)   |
|  - Tinh chỉnh thời hằng lọc mong muốn|                   |  - Tối ưu đồng thời: bám quỹ đạo Jr   |
|  - Cân bằng Servomechanism/Regulator |                   |    & triệt tiêu nhiễu Jd              |
|  - Đảm bảo độ nhạy cực đại Ms        |                   |  - Lựa chọn nghiệm trên Pareto Front  |
+--------------------------------------+                   +---------------------------------------+
                   |                                                           |
                   +-----------------------------+-----------------------------+
                                                 |
                                                 v
                  +-------------------------------------------------------------+
                  |               KIỂM CHỨNG & ĐÁNH GIÁ ĐỘ BỀN VỮNG             |
                  |  - Phân tích cấu trúc M-Δ & Giá trị suy biến SSV (μ < 1)    |
                  |  - Thực nghiệm thời gian thực: Hệ 4 bồn nước (Quad-Tank)   |
                  +-------------------------------------------------------------+

Khoảng trống nghiên cứu (research gap) trọng tâm mà luận án xác định xuất phát từ thực tế: phần lớn các nghiên cứu về bộ điều khiển tỉ lệ - tích phân - vi phân bậc phân số (Fractional-Order PID - FOPID, $PI^\lambda D^\mu$) từ trước đến nay chỉ tập trung trên hệ thống đơn biến một ngõ vào - một ngõ ra (Single-Input Single-Output - SISO). Việc mở rộng FOPID sang hệ đa biến MIMO gặp phải ba thách thức học thuật cốt lõi: (1) Sự bùng nổ khối lượng tính toán khi giải tích toán học ma trận phân ly và ma trận tương đương sau phân ly, (2) Khó khăn trong việc thiết lập quy luật chỉnh định tường minh cho 5 thông số ($K_c, \tau_I, \tau_D, \lambda, \mu$) trên từng vòng lặp khi bậc của hệ thống tăng cao, và (3) Đảm bảo tính ổn định bền vững (robust stability) trước sai số mô hình hóa không thể tránh khỏi trong môi trường thực tế.

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

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để loại bỏ hoàn toàn ảnh hưởng của ma trận thời gian trễ trong hệ đa biến mà không làm tăng độ phức tạp cấu trúc của bộ điều khiển?
  • Câu hỏi nghiên cứu 2 (RQ2): Giải pháp nào tối ưu hóa việc xấp xỉ rút gọn các hàm truyền tương đương bậc cao phức tạp sau phân ly về dạng hàm truyền phân số chuẩn tắc?
  • Câu hỏi nghiên cứu 3 (RQ3): Thiết lập quy tắc chỉnh định thông số FOPID như thế nào để cân bằng tối ưu giữa chất lượng đáp ứng bám giá trị đặt (servomechanism) và khả năng kháng nhiễu tải (regulator) đồng thời thỏa mãn biên độ ổn định bền vững?
  • Giả thuyết khoa học 1 (H1): Việc kết hợp kỹ thuật phân ly đơn giản hóa (Simplified Decoupling) với bộ dự báo Smith (Smith Predictor) cho hệ đa biến sẽ triệt tiêu hoàn toàn thành phần trễ pha khỏi mẫu số phương trình đặc trưng vòng kín.
  • Giả thuyết khoa học 2 (H2): Ứng dụng giải thuật tối ưu hóa bầy đàn (Particle Swarm Optimization - PSO) sẽ giảm thiểu sai số xấp xỉ hàm truyền tương đương vượt trội so với các phương pháp khớp hệ số cổ điển.
  • Giả thuyết khoa học 3 (H3): Bộ điều khiển FOPID được chỉnh định bằng cấu trúc mô hình nội (IMC) hoặc tối ưu bầy đàn đa mục tiêu (MOPSO) sẽ cải thiện đáng kể chỉ số tích phân sai số tuyệt đối ($IAE$), sai số tích lũy nhân trọng số thời gian ($ITAE$) và tổng biến thiên tín hiệu điều khiển ($TV$) so với IOPID truyền thống.

Khung lý thuyết của công trình được xây dựng vững chắc trên nền tảng: Giải tích phân số Riemann-Liouville, lý thuyết điều khiển mô hình nội (Internal Model Control - IMC), lý thuyết phân tích độ ổn định bền vững thông qua cấu trúc $M-\Delta$ và Giá trị suy biến có cấu trúc (Structured Singular Value - SSV, ký hiệu $\mu$), cùng giải thuật tiến hóa bầy đàn đa mục tiêu (MOPSO). Phạm vi nghiên cứu bao quát từ mô phỏng số trên các mô hình chuẩn quốc tế (Benchmark models) gồm hệ 2×2 (Tháp Vinante & Luyben, Bộ tách dầu nặng Heavy Oil Fractionator), hệ 3×3 (Tháp chưng cất Ogunnaike & Ray), hệ 4×4 (Hệ thống nhiệt độ HVAC 4 phòng) cho đến thực nghiệm phần cứng thời gian thực trên hệ thống bốn bồn nước liên kết (Quadruple Tank System).

Literature Review và Positioning

Lịch sử tính toán phân số khởi nguồn từ bức thư trao đổi kinh điển năm 1695 giữa Leibniz và L'Hôpital khi bàn về ý nghĩa của đạo hàm bậc $1/2$, trong đó Leibniz đã dự báo: "Nó sẽ dẫn đến tranh cãi mà một ngày nào đó sẽ tạo ra các kết quả có ích". Trải qua hơn ba thế kỷ phát triển nền tảng toán học bởi Euler, Lagrange, Abel, Liouville và Riemann, tính toán phân số chính thức thâm nhập vào lý thuyết điều khiển tự động hiện đại nhờ công trình đột phá của Podlubny (1999) khi giới thiệu cấu trúc bộ điều khiển $PI^\lambda D^\mu$. Nhóm nghiên cứu của Chen và các cộng sự (2007, 2009) đã phát triển các công cụ xấp xỉ số học, thiết lập quy luật chỉnh định cho hệ bậc một có trễ (FOPTD) và mở rộng sang các hệ thống servo cơ điện chính xác. Song song đó, nhóm trường phái CRONE (Oustaloup, Trigeassou, Malti) tại Pháp đã khẳng định tính ưu việt của phân số trong việc thiết kế hệ thống điều khiển bền vững miền tần số và nhận dạng hệ thống phi tuyến.

                           TIẾN TRÌNH PHÁT TRIỂN LÝ THUYẾT ĐIỀU KHIỂN PHÂN SỐ
                                      VÀ VỊ TRÍ CỦA LUẬN ÁN
                                                |
+---------------------------------------------------------------------------------------------------+
| 1695: Khởi xướng giải tích phân số (Bức thư lịch sử giữa Leibniz và L'Hôpital)                   |
| 1823 - Thế kỷ 20: Hoàn thiện toán học giải tích phân số (Abel, Liouville, Riemann, Caputo)         |
+---------------------------------------------------------------------------------------------------+
                                                |
                                                v
+---------------------------------------------------------------------------------------------------+
| 1999: Podlubny đề xuất bộ điều khiển FOPID (PI^λ D^μ) mở rộng bậc phân số                        |
| 2000 - 2015: Trường phái CRONE (Oustaloup et al.) & Nhóm Chen et al. phát triển SISO FOPID        |
+---------------------------------------------------------------------------------------------------+
                                                |
                                                v
+---------------------------------------------------------------------------------------------------+
| TIẾP CẬN ĐA BIẾN TRUYỀN THỐNG:                                                                    |
| - Phân ly nghịch (Garrido et al., 2011, 2013): Bế tắc cấu trúc khi n ≥ 3                         |
| - Phân ly đơn giản hóa (Vu & Lee, 2010): Gánh nặng giải tích rút gọn ma trận bậc cao              |
| - Điều khiển phân tán thuần túy: Thất bại khi tương tác chéo mạnh (RGA/DRGA lệch)                |
+---------------------------------------------------------------------------------------------------+
                                                |
                                                v
+===================================================================================================+
| VỊ TRÍ ĐỘT PHÁ CỦA LUẬN ÁN (Võ Lâm Chương, 2024):                                                 |
| 1. Thiết lập cấu trúc F-SDSP: Phân ly đơn giản hóa + Bộ dự báo Smith cho hệ MIMO n x n           |
| 2. Đột phá rút gọn mô hình: Tích hợp PSO tối ưu hóa hàm truyền phân số tương đương                |
| 3. Chỉnh định kép: IMC phân số (cho 2x2) & MOPSO đa mục tiêu (Jr, Jd, Ms) (cho 3x3, 4x4)          |
| 4. Khung ổn định bền vững M-Δ / SSV (μ < 1) & Thực nghiệm phần cứng Quad-Tank thời gian thực     |
+===================================================================================================+

Trong bức tranh tổng quan về điều khiển hệ đa biến MIMO, tồn tại hai nhánh tiếp cận chính:

  1. Điều khiển không tập trung (Decentralized/Multi-loop Control): Sử dụng ma trận độ lợi tương quan RGA (Bristol, 1966), chỉ số Niederlinski (NI), hoặc ma trận tương quan động DRGA để ghép cặp biến. Vu và Lee (2010) đề xuất phương pháp thiết kế độc lập dựa trên hàm truyền vòng hở tương đương (EOTF). Nhược điểm chí mạng của nhánh này là chỉ hiệu quả khi tương tác giữa các vòng lặp ở mức yếu; khi tương tác chéo tăng cao, hệ thống dễ rơi vào mất ổn định.
  2. Điều khiển tập trung và kỹ thuật phân ly (Decoupling Techniques): Để giải quyết tương tác mạnh, các kỹ thuật phân ly lý tưởng, phân ly nghịch (Inverted Decoupling) và phân ly đơn giản hóa (Simplified Decoupling) được áp dụng. Garrido và các cộng sự (2011, 2013) đã có những đóng góp xuất sắc về phân ly nghịch kết hợp IMC cho hệ 2×2. Tuy nhiên, khi mở rộng lên hệ $n \times n$ ($n \ge 3$), việc lựa chọn cấu trúc ma trận phân ly trở nên vô cùng phức tạp và không có quy luật tổng quát. Mặt khác, Vu và Lee (2010) phát triển phân ly đơn giản hóa cho hệ $n \times n$, nhưng việc rút gọn các hàm truyền phức tạp sau phân ly bằng các phép toán đại số giải tích bậc nguyên truyền thống gặp phải gánh nặng tính toán khổng lồ.

Về vấn đề bù thời gian trễ, bộ dự báo Smith (Smith Predictor - SP) kinh điển và bộ dự báo Smith có lọc (Filtered Smith Predictor - FSP) của Normey-Rico và Camacho (2007) đã thành công rực rỡ trên hệ SISO và bước đầu áp dụng trên MIMO kết hợp phân ly nghịch. Tuy nhiên, việc kết hợp giữa phân ly đơn giản hóa và bộ dự báo Smith cho hệ đa biến bậc cao sử dụng bộ điều khiển phân số hoàn toàn chưa từng được công bố trong y văn quốc tế.

Luận án của Võ Lâm Chương đã định vị chính xác vào khoảng trống học thuật này. Tác giả không chỉ thiết lập cấu trúc điều khiển phân ly đơn giản hóa kết hợp bộ dự báo Smith và điều khiển phân số (F-SDSP) cho hệ $n \times n$, mà còn giải quyết triệt để bài toán rút gọn mô hình bằng giải thuật tiến hóa PSO và xây dựng khung chỉnh định thông số toàn diện qua IMC và MOPSO.

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

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

Luận án tạo nên bước chuyển dịch mô hình (paradigm shift) từ không gian điều khiển bậc nguyên truyền thống sang không gian bậc phân số tổng quát cho hệ thống MIMO đa biến có trễ.

  • Mở rộng lý thuyết điều khiển $PI^\lambda D^\mu$ của Podlubny: Thiết lập mối quan hệ toán học giữa bậc tích phân phân số $\lambda$ và bậc đạo hàm phân số $\mu$ ($0 < \lambda, \mu < 2$) đối với khả năng bù pha và tạo đặc tính biên tần số mong muốn trong hệ đa biến. Khác với bậc nguyên chỉ có độ dốc Bode biên độ cố định $-20\text{ dB/decade}$ hoặc $+20\text{ dB/decade}$, toán tử vi tích phân phân số $s^\mu$ và $s^{-\lambda}$ cho phép tạo độ dốc tùy ý $-20\lambda\text{ dB/decade}$ và $+20\mu\text{ dB/decade}$, giúp hệ thống đạt độ dự trữ pha phẳng (isofrequency phase margin), triệt tiêu hiện tượng vọt lố và tăng cường tính bền vững.
  • Mô hình toán học động học phi tuyến tổng quát: Chứng minh vai trò của bậc phân số trong việc mô tả chính xác các phương trình động lực học phức tạp như phương trình Langevin phân số mô tả chuyển động Brownian có nhớ, phương trình Van der Pol phân số (FO-VdP) bảo toàn chu trình giới hạn tự dao động với tính linh hoạt cao, và phương trình truyền nhiệt Bernoulli phân số.
  • Mệnh đề toán học về cấu trúc triệt tiêu độ trễ (Proposition 1): Luận án chứng minh rằng thông qua việc tích hợp ma trận dự báo Smith đa biến vào cấu trúc phân ly đơn giản hóa, ma trận hàm truyền vòng kín tương đương loại bỏ hoàn toàn các số hạng trễ pha $e^{-\theta_{ij}s}$ ra khỏi phương trình đặc trưng của hệ kín, chuyển bài toán điều khiển hệ có trễ phức tạp thành bài toán điều khiển hệ không trễ tương đương.

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

Khung phân tích của luận án tích hợp 4 trụ cột lý thuyết vững chắc:

  1. Lý thuyết Giải tích phân số Riemann-Liouville: Sử dụng định nghĩa tích phân và đạo hàm phân số: $$_0D_t^{-\alpha} f(t) = \frac{1}{\Gamma(\alpha)} \int_0^t (t-\tau)^{\alpha-1} f(\tau) d\tau$$ $$_0D_t^{\alpha} f(t) = \frac{d}{dt} \left[ _0D_t^{-(1-\alpha)} f(t) \right] = \frac{1}{\Gamma(1-\alpha)} \frac{d}{dt} \int_0^t \frac{f(\tau)}{(t-\tau)^\alpha} d\tau$$ với biến đổi Laplace tương ứng $\mathcal{L}{_0D_t^{\pm \alpha} f(t)} = s^{\pm \alpha} F(s)$ dưới điều kiện biên triệt tiêu.

  2. Kỹ thuật phân ly đơn giản hóa và xấp xỉ mô hình qua PSO: "NCS đề xuất sử dụng giải thuật tối ưu hóa bầy đàn (PSO) trong việc rút gọn và đơn giản hóa các hàm truyền thành phần của ma trận phân ly cũng như ma trận tương đương của hệ sau khi phân ly." Ma trận phân ly $D(s)$ có các phần tử đường chéo $d_{ii}(s) = 1$, các phần tử ngoài đường chéo $d_{ij}(s)$ được xác định nhằm biến đổi ma trận hàm truyền thực $G(s)$ thành ma trận đường chéo $Q(s) = G(s)D(s)$. Giải thuật PSO tối ưu hóa hàm mục tiêu sai số bình phương giữa đáp ứng tần số/thời gian của hàm gốc và hàm phân số rút gọn dạng: $$q_{ii}^(s) = \frac{K_i}{\tau_i s^{\alpha_i} + 1} e^{-\theta_i s} \quad \text{hoặc} \quad q_{ii}^(s) = \frac{K_i}{\tau_{1i} s^{\alpha_i} + \tau_{2i} s^{\beta_i} + 1} e^{-\theta_i s}$$

  3. Cơ chế chỉnh định thông số đa tầng:

    • Hệ bậc thấp (2×2 TITO): Áp dụng cấu trúc IMC phân số, quy đổi hàm truyền bộ điều khiển $Q_{IMC}(s) = q_{ii-}^{-1}(s) f_i(s)$ về dạng $FOPID$ chuẩn: $$G_c(s) = K_c \left( 1 + \frac{1}{\tau_I s^\lambda} + \tau_D s^\mu \right) \left( \frac{1}{\tau_f s + 1} \right)$$
    • Hệ bậc cao (3×3, 4×4): Ứng dụng thuật toán MOPSO với hai hàm mục tiêu xung đột: $$J_r = \int_0^\infty |e_r(t)| dt \quad (IAE_{\text{setpoint}}), \qquad J_d = \int_0^\infty |e_d(t)| dt \quad (IAE_{\text{disturbance}})$$ Nghiệm tối ưu được lựa chọn trên biên Pareto (Pareto Front) kết hợp điều kiện ràng buộc hàm độ nhạy cực đại $M_s \le 1.4 - 2.0$ nhằm bảo đảm chất lượng động học và độ ổn định bền vững.

  4. Khung đánh giá ổn định bền vững $M-\Delta$ và SSV: Thiết lập ma trận truyền $M(s)$ từ sơ đồ khối hệ thống kín với sai số nhân ngõ ra $\Delta(s)$ thỏa mãn $|\Delta(s)|\infty \le 1$. Tiêu chuẩn ổn định bền vững được kiểm chứng qua đại lượng suy biến có cấu trúc: $$\mu{\mathbf{\Delta}}(M(j\omega)) < 1, \quad \forall \omega \in [0, \infty)$$

                                 KHUNG PHÂN TÍCH TỔNG THỂ CỦA LUẬN ÁN
                                                  |
+---------------------------------------------------------------------------------------------------+
| TRỤ CỘT 1: TOÁN HỌC PHÂN SỐ RIEMANN-LIOUVILLE                                                      |
| L{_0D_t^α f(t)} = s^α F(s) | L{_0D_t^-λ f(t)} = s^-λ F(s) | Xấp xỉ Oustaloup liên tục miền tần số |
+---------------------------------------------------------------------------------------------------+
                                                  |
                                                  v
+---------------------------------------------------------------------------------------------------+
| TRỤ CỘT 2: CẤU TRÚC PHÂN LY ĐƠN GIẢN HÓA & DỰ BÁO SMITH (F-SDSP)                                 |
| - Ma trận phân ly D(s) với d_ii = 1 | Ma trận sau phân ly Q(s) = G(s)D(s)                         |
| - Smith Predictor: Loại bỏ triệt để e^-θs khỏi mẫu số phương trình đặc trưng vòng kín            |
| - Rút gọn hàm truyền q_ii(s), d_ij(s) qua giải thuật bầy đàn PSO                                  |
+---------------------------------------------------------------------------------------------------+
                                                  |
                                                  v
+---------------------------------------------------------------------------------------------------+
| TRỤ CỘT 3: CƠ CHẾ CHỈNH ĐỊNH THÔNG SỐ TỐI ƯU                                                      |
| - Hệ 2x2: Cấu trúc Mô hình nội phân số (IMC) | Tinh chỉnh thời hằng mong muốn                     |
| - Hệ 3x3, 4x4: Tối ưu bầy đàn đa mục tiêu (MOPSO) | Tối ưu hóa Jr & Jd trên Pareto Front         |
+---------------------------------------------------------------------------------------------------+
                                                  |
                                                  v
+---------------------------------------------------------------------------------------------------+
| TRỤ CỘT 4: PHÂN TÍCH ỔN ĐỊNH BỀN VỮNG M-Δ & SSV                                                   |
| - Bất định sai số nhân ngõ ra Δ(s) | Đánh giá Structured Singular Value μ(M) < 1                  |
+---------------------------------------------------------------------------------------------------+

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ẽ thế giới quan thực chứng (positivism paradigm), kết hợp hài hòa giữa mô hình hóa toán học giải tích, tính toán mô phỏng tiến hóa ngẫu nhiên và kiểm chứng thực nghiệm phần cứng vật lý. Thiết kế nghiên cứu được phân cấp đa tầng:

  • Cấp độ mô hình hóa giải tích: Xây dựng cấu trúc ma trận phân ly và phương trình tương đương.
  • Cấp độ tối ưu hóa tham số: Sử dụng thuật toán heuristic PSO và MOPSO.
  • Cấp độ kiểm chứng mô phỏng: Đánh giá định lượng trên 4 mô hình đối chứng tiêu chuẩn quốc tế.
  • Cấp độ thực nghiệm thời gian thực: Triển khai giải thuật điều khiển trực tiếp trên đối tượng vật lý thông qua MATLAB Real-Time Windows Target.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu được thực hiện qua các giai đoạn nghiêm ngặt:

  1. Thu thập dữ liệu và Nhận dạng đa biến: Mở rộng phương pháp bình phương cực tiểu (Least Squares - LS) từ hệ đơn biến SISO sang hệ đa biến MIMO thông qua kỹ thuật Phân rã ma trận (Matrix Fraction Description - MFD). Trên hệ thực nghiệm bốn bồn nước, tín hiệu kích thích nhị phân giả ngẫu nhiên (PRBS) được đưa vào hai máy bơm ngõ vào, thu thập dữ liệu đáp ứng mức nước ($h_1, h_2$) với chu kỳ lấy mẫu chính xác.
  2. Thiết kế ma trận phân ly và rút gọn mô hình qua PSO: Tính toán các phần tử $d_{ij}(s)$ của ma trận $D(s)$ và các hàm truyền đường chéo $q_{ii}(s)$. Áp dụng giải thuật PSO với kích thước quần thể bầy đàn 50 cá thể, 100 thế hệ lặp để tìm bộ tham số tối ưu cho hàm truyền phân số xấp xỉ, đảm bảo tối thiểu hóa chỉ số sai số toàn phương trung bình (MSE).
  3. Tổng hợp bộ điều khiển FOPI/FOPID: Thiết lập 5 tham số cho từng vòng lặp bằng giải tích IMC (đối với hệ 2×2) và thuật toán MOPSO với lưu trữ Pareto (đối với hệ 3×3 và 4×4).
  4. Kiểm tra tính hợp thức và độ bền vững: Xây dựng sơ đồ phân tích cấu trúc $M-\Delta$ với ma trận trọng số bất định $W_m(s)$, quét trên toàn dải tần số $\omega \in [10^{-3}, 10^3]\text{ rad/s}$ để khẳng định đường cong SSV luôn nằm dưới ngưỡng giới hạn đơn vị ($\mu < 1$).

Data và phân tích

Các mô hình chuẩn đối chứng (Benchmark Datasets) được sử dụng để thẩm định bao gồm:

  • Tháp Vinante & Luyben (VL Column - 2×2 TITO): Đối tượng chưng cất phân tách methanol-nước với độ trễ lớn và tương tác mạnh. $$G_{VL}(s) = \begin{bmatrix} \frac{-2.2 e^{-s}}{7s+1} & \frac{1.3 e^{-0.3s}}{7s+1} \ \frac{-2.8 e^{-1.8s}}{9.5s+1} & \frac{4.3 e^{-0.35s}}{9.2s+1} \end{bmatrix}$$
  • Bộ tách dầu nặng (Heavy Oil Fractionator - 2×2): Hệ thống có quán tính thời gian lớn và thời gian trễ phi đối xứng. $$G_{HOF}(s) = \begin{bmatrix} \frac{4.05 e^{-27s}}{50s+1} & \frac{1.77 e^{-28s}}{60s+1} \ \frac{5.39 e^{-18s}}{50s+1} & \frac{5.72 e^{-14s}}{60s+1} \end{bmatrix}$$
  • Tháp chưng cất Ogunnaike & Ray (OR Column - 3×3): Hệ thống chưng cất ba ngõ vào - ba ngõ ra với tương tác động học đa chiều phức tạp.
  • Hệ thống điều hòa không khí HVAC 4 phòng (4×4 MIMO): Mô hình nhiệt độ không gian với 4 biến vào là lưu lượng nước lạnh và 4 biến ra là nhiệt độ phòng.
  • Hệ thực nghiệm bốn bồn nước liên kết (Quadruple Tank): Hệ thống phần cứng bao gồm 4 bồn chứa mica trong suốt, 2 máy bơm DC điều khiển bằng điện áp qua card thu thập dữ liệu (DAQ), 4 cảm biến áp suất đo mức nước chuyển đổi chuẩn tín hiệu $0-5\text{V}$, kết nối máy tính điều khiển qua phần mềm MATLAB/Simulink chạy Real-Time Windows Target.

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

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

Kết quả mô phỏng và thực nghiệm trong luận án đã mang lại những phát hiện đột phá có ý nghĩa khoa học và thực tiễn cao:

  1. Vượt trội toàn diện về chỉ số tích phân sai số: Trên tất cả các mô hình thử nghiệm, cấu trúc F-SDSP đề xuất cho kết quả cải thiện vượt bậc so với các công trình công bố quốc tế trước đây (như phương pháp phân ly nghịch IMC của Garrido et al., phương pháp phân ly đơn giản hóa của Vu & Lee, và các bộ điều khiển PID bậc nguyên kinh điển). Cụ thể, trên mô hình tháp Vinante & Luyben (VL), chỉ số $IAE$ khi thay đổi giá trị đặt và khi có nhiễu tác động giảm từ $25%$ đến $48%$, trong khi thời gian xác lập (settling time) giảm hơn $35%$.
  2. Giảm thiểu độ biến thiên tín hiệu điều khiển (Total Variation - TV): Một phát hiện phản trực giác nhưng mang giá trị kỹ thuật sâu sắc là: mặc dù bộ điều khiển phân số có bậc động học phức tạp hơn, nhưng tín hiệu điều khiển $u(t)$ tạo ra lại mượt mà hơn đáng kể. Tổng biến thiên $TV = \sum |u_{k+1} - u_k|$ của F-SDSP giảm từ $30%$ đến $55%$ so với bộ điều khiển IOPID bậc nguyên. Điều này chứng minh khâu vi phân phân số $s^\mu$ ($0 < \mu < 1$) không khuếch đại nhiễu tần số cao mạnh như khâu vi phân bậc 1 ($s^1$), giúp bảo vệ cơ cấu chấp hành (van điều khiển, máy bơm) khỏi hiện tượng mài mòn cơ khí và bão hòa công suất.
  3. Độ chính xác vượt trội của giải thuật PSO trong rút gọn hàm truyền: So sánh với phương pháp so khớp hệ số (Coefficient Matching) và phương pháp sai số dự báo (PEM), PSO giúp tìm được hàm truyền phân số tương đương $q_{ii}^*(s)$ có sai lệch pha và biên độ gần như trùng khớp hoàn toàn với hàm gốc $q_{ii}(s)$ trên toàn dải tần số làm việc, đặc biệt hiệu quả với các hệ bậc cao 3×3 và 4×4.
  4. Khẳng định độ ổn định bền vững tuyệt đối qua SSV: Đường biểu diễn giá trị suy biến có cấu trúc $\mu_{\mathbf{\Delta}}(M(j\omega))$ cho tất cả các trường hợp (Tháp VL, Heavy Oil, Tháp OR, HVAC và Hệ 4 bồn nước) đều đạt giá trị cực đại $\mu_{\max} \le 0.75 - 0.88$, nhỏ hơn rất nhiều so với ngưỡng giới hạn $1.0$. Điều này chứng minh hệ thống duy trì độ ổn định tuyệt đối ngay cả khi các tham số của đối tượng thực bị biến thiên $\pm 20%$ so với mô hình danh định.
  5. Tính khả thi tuyệt đối trong môi trường thời gian thực: Kết quả thực nghiệm trên hệ bốn bồn nước (Quadruple Tank) chứng minh thuật toán FOPID sau khi rời rạc hóa bằng phương pháp Oustaloup vận hành hoàn hảo với chu kỳ lấy mẫu thời gian thực $T_s = 0.1\text{s}$ trên nền MATLAB Real-Time Windows Target mà không xảy ra hiện tượng trễ nhịp tính toán (CPU overrun). Mức chất lỏng tại các bồn $h_1, h_2$ bám chính xác quỹ đạo đặt với độ vọt lố $< 2%$ và triệt tiêu nhiễu tải nhanh chóng trong vòng dưới 25 giây.
+---------------------------------------------------------------------------------------------------+
|                     BẢNG TỔNG HỢP SO SÁNH HIỆU NĂNG ĐIỀU KHIỂN (BENCHMARKS)                       |
+------------------------------+---------------------------+----------------------------------------+
| Tiêu chí đánh giá            | Phương pháp kinh điển     | Phương pháp đề xuất (F-SDSP Luận án)   |
|                              | (IOPID / Inverted Dec.)   |                                        |
+------------------------------+---------------------------+----------------------------------------+
| Tích phân sai số tuyệt đối   | Cơ sở chuẩn định mức      | Giảm từ 25% đến 48%                    |
| (IAE / ITAE)                 | (Baseline = 100%)         | (Đáp ứng quá độ nhanh, ít dao động)    |
+------------------------------+---------------------------+----------------------------------------+
| Tổng biến thiên tín hiệu     | Cao (Dễ rung lắc van,     | Giảm từ 30% đến 55%                    |
| điều khiển (Total Var. - TV) | bão hòa chấp hành)        | (Tín hiệu êm thuận, chống mòn cơ khí)  |
+------------------------------+---------------------------+----------------------------------------+
| Độ nhạy cực đại (Ms)         | 1.8 - 2.4 (Biên kém)      | 1.3 - 1.6 (Tối ưu theo biên Pareto)    |
+------------------------------+---------------------------+----------------------------------------+
| Đỉnh suy biến cấu trúc       | 0.95 - 1.15               | 0.75 - 0.88                            |
| (Peak SSV / μ_max)           | (Tiệm cận/vượt ngưỡng bất | (Nhỏ hơn 1.0 -> Ổn định bền vững cao)  |
|                              | định khi có sai số mô hình|                                        |
+------------------------------+---------------------------+----------------------------------------+
| Khả năng mở rộng cho MIMO    | Rất khó khăn khi n ≥ 3    | Hoàn toàn tổng quát cho hệ n x n       |
| bậc cao                      | (Bùng nổ giải tích)       | (Nhờ tối ưu hóa tiến hóa PSO/MOPSO)    |
+------------------------------+---------------------------+----------------------------------------+

Implications đa chiều

  • Về mặt học thuật: Cung cấp khung phương pháp luận hoàn chỉnh kết hợp giải tích phân số với các kỹ thuật phân ly đa biến hiện đại, mở đường cho các nghiên cứu lý thuyết mở rộng FOPID trên không gian trạng thái phân số đa chiều.
  • Về mặt kỹ thuật công nghệ: Cung cấp quy trình tự động hóa thiết kế bộ điều khiển: từ nhận dạng dữ liệu vào-ra MFD-LS, rút gọn mô hình qua PSO, đến tìm thông số FOPID tối ưu qua MOPSO.
  • Về mặt ứng dụng thực tiễn: Ứng dụng trực tiếp vào các ngành công nghiệp mũi nhọn như lọc - hóa dầu (tháp chưng cất, tháp phân tách), hóa chất, xử lý nước thải công nghiệp và hệ thống quản lý năng lượng nhiệt tòa nhà thông minh (HVAC), giúp tiết kiệm điện năng và nâng cao tuổi thọ thiết bị.

Limitations và Future Research

Nhằm duy trì tính khách quan khoa học, luận án thẳng thắn thừa nhận các giới hạn kỹ thuật:

  1. Giới hạn về cấu trúc ma trận vuông: Luận án mới chỉ tập trung giải quyết triệt để cho hệ đa biến vuông ($n \times n$ ngõ vào - ngõ ra). Đối với hệ không vuông ($m \times n$ với $m \ne n$), việc giải ma trận phân ly đòi hỏi các phép nghịch đảo giả ma trận phức tạp chưa được đề cập sâu.
  2. Giả định thời gian trễ bất biến: Các thông số thời gian trễ $\theta_{ij}$ được giả định là hằng số đã biết hoặc nhận dạng tĩnh. Trong các hệ thống dòng chảy biến thiên mạnh, độ trễ thời gian có thể thay đổi theo thời gian (time-varying delay).
  3. Chi phí tính toán ngoại tuyến của MOPSO: Thuật toán tối ưu hóa đa mục tiêu MOPSO đòi hỏi thời gian tính toán mô phỏng ban đầu tương đối lớn để bầy hạt hội tụ về đường biên Pareto.
  4. Mô hình tuyến tính quanh điểm làm việc: Nhận dạng hệ thống dựa trên mô hình tuyến tính hóa cục bộ MFD-LS; do đó, khi hệ thống vận hành trôi dạt xa khỏi điểm làm việc danh định, chất lượng điều khiển có thể bị suy giảm.

Hướng nghiên cứu phát triển tiếp theo bao gồm:

  • Mở rộng cấu trúc F-SDSP cho hệ thống đa biến không vuông ($m \times n$) và hệ có trễ biến thiên ngẫu nhiên theo thời gian.
  • Kết hợp giải thuật FOPID với cấu trúc thích nghi (Adaptive FOPID) hoặc mạng nơ-ron sâu để tự động chỉnh định trực tuyến tham số $(\lambda, \mu)$ khi đối tượng có tính phi tuyến mạnh.
  • Tích hợp điều khiển trượt bậc phân số (Fractional-Order Sliding Mode Control) nhằm nâng cao khả năng kháng nhiễu bất định tuyệt đối.
  • Hiện thực hóa thuật toán FOPID trực tiếp trên các chip vi xử lý nhúng chuyên dụng công nghiệp như DSP (Digital Signal Processor), FPGA hoặc PLC công nghiệp thay vì máy tính cá nhân.
  • Triển khai thử nghiệm quy mô pilot tại các nhà máy lọc dầu và hóa chất thực tế.

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

  • Tác động học thuật: Các kết quả của luận án đã được công bố trên các tạp chí quốc tế chuyên ngành uy tín thuộc danh mục ISI/Scopus (như hợp tác với nhóm nghiên cứu của GS. Jae Hak Jung, Đại học Yeungnam, Hàn Quốc). Công trình là tài liệu tham khảo giá trị cao cho các nghiên cứu tiếp theo về điều khiển phân số đa biến, dự kiến thu hút trích dẫn khoa học mạnh mẽ trong cộng đồng tự động hóa.
  • Chuyển đổi công nghiệp: Cung cấp giải pháp nâng cấp trực tiếp cho các hệ thống điều khiển quá trình hiện hữu mà không cần thay đổi phần cứng đắt tiền, chỉ cần thay đổi thuật toán điều khiển trong phần mềm DCS/SCADA từ IOPID sang FOPID.
  • Hiệu quả kinh tế - xã hội: Giảm dao động quá độ trong các tháp chưng cất hóa dầu giúp tăng độ tinh khiết sản phẩm thêm $1.5 - 3%$, giảm tiêu hao năng lượng nhiệt tái đun sôi; tối ưu hóa điều khiển nhiệt độ HVAC giúp tiết kiệm từ $8 - 15%$ điện năng tiêu thụ trong các tòa nhà thương mại quy mô lớn.

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

  • Nghiên cứu sinh và Giảng viên ngành Tự động hóa/Cơ điện tử: Tiếp cận được phương pháp luận kết hợp giải tích phân số với các công cụ tối ưu hóa tiến hóa hiện đại và khung phân tích ổn định bền vững $M-\Delta$.
  • Kỹ sư R&D và Tự động hóa quá trình: Sở hữu công cụ thiết kế bộ điều khiển có khả năng giải quyết triệt để bài toán tương tác chéo và thời gian trễ trong các dây chuyền sản xuất thực tế.
  • Nhà quản trị năng lượng và vận hành nhà máy: Nắm bắt giải pháp nâng cao hiệu suất vận hành, giảm chi phí bảo trì cơ cấu chấp hành nhờ tính năng triệt tiêu rung lắc tín hiệu điều khiển của FOPID.

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ì? Đó là việc thiết lập thành công cấu trúc điều khiển phân ly đơn giản hóa kết hợp bộ dự báo Smith và bộ điều khiển PID bậc phân số (F-SDSP) cho hệ đa biến tổng quát $n \times n$. Công trình đã mở rộng lý thuyết điều khiển phân số của Podlubny từ không gian đơn biến SISO sang không gian đa biến MIMO có trễ, loại bỏ hoàn toàn số hạng trễ ra khỏi phương trình đặc trưng vòng kín và chứng minh tính ổn định bền vững thông qua cấu trúc $M-\Delta$ với $\mu_{\max} < 1$.

  2. Đột phá phương pháp luận so với các nghiên cứu tiền nhiệm quốc tế là gì? So với kỹ thuật phân ly nghịch của Garrido et al. (2011, 2013) vốn bị bế tắc khi bậc hệ thống $n \ge 3$, và kỹ thuật phân ly đơn giản hóa của Vu & Lee (2010) vốn bị quá tải tính toán đại số giải tích, luận án đã tiên phong ứng dụng giải thuật tối ưu bầy đàn PSO để xấp xỉ rút gọn hàm truyền phân số tương đương, đồng thời tích hợp MOPSO giải quyết bài toán tối ưu đa mục tiêu (cân bằng giữa $IAE_{\text{setpoint}}$ và $IAE_{\text{disturbance}}$ dưới ràng buộc $M_s$) cho hệ bậc cao 3×3 và 4×4.

  3. Phát hiện thực nghiệm gây bất ngờ nhất là gì? Phát hiện về sự suy giảm mạnh mẽ của tổng biến thiên tín hiệu điều khiển ($TV$ giảm $30 - 55%$) khi sử dụng bộ điều khiển FOPID so với IOPID. Về mặt trực quan, bộ điều khiển bậc phân số có cấu trúc toán học phức tạp hơn thường bị nghi ngại sẽ gây dao động ngõ ra điều khiển; tuy nhiên, dữ liệu thực nghiệm đã chứng minh bậc đạo hàm không nguyên $\mu \in (0, 1)$ tạo ra hiệu ứng lọc thông dải tự nhiên, giảm sốc áp lực cho van và bơm điều khiển.

  4. Giao thức tái lập nghiên cứu (Replication Protocol) có được cung cấp đầy đủ không? Luận án cung cấp chi tiết toàn bộ ma trận hàm truyền của các mô hình chuẩn quốc tế (Vinante-Luyben, Heavy Oil, Ogunnaike-Ray, HVAC), thông số kích thước bầy đàn/số thế hệ của PSO/MOPSO, sơ đồ khối phần cứng hệ bốn bồn nước (kích thước bồn, dải đo cảm biến $0-5\text{V}$, điện áp bơm $0-12\text{V}$) và sơ đồ cấu hình Simulink chạy trên nền tảng Real-Time Windows Target.

  5. Chương trình nghiên cứu 10 năm tiếp theo được định hình như thế nào? Chương trình nghiên cứu dài hạn tập trung vào 3 giai đoạn: (1) 2024-2027: Phát triển lý thuyết F-SDSP cho hệ không vuông $m \times n$ và hệ có trễ biến thiên; (2) 2027-2030: Nhúng toàn diện giải thuật FOPID tự chỉnh định trên nền tảng phần cứng công nghiệp chuyên dụng (DSP/FPGA/PLC); (3) 2030-2034: Ứng dụng công nghiệp hóa quy mô lớn trong các nhà máy lọc dầu, nhà máy điện sinh học và lưới nhiệt thông minh đô thị.

Kết luận

  1. Luận án đã giải quyết trọn vẹn bài toán điều khiển các quá trình công nghiệp đa biến phức tạp có tương tác chéo mạnh và độ trễ thời gian lớn bằng việc đề xuất cấu trúc điều khiển F-SDSP tiên phong.
  2. Ứng dụng thành công giải thuật tối ưu hóa bầy đàn (PSO) trong việc xấp xỉ rút gọn các hàm truyền ma trận phân ly phức tạp về dạng hàm truyền phân số chuẩn tắc, hóa giải rào cản tính toán đại số giải tích bậc cao.
  3. Xây dựng bộ quy tắc chỉnh định thông số FOPID toàn diện: sử dụng cấu trúc mô hình nội (IMC) cho hệ 2×2 và tối ưu hóa bầy đàn đa mục tiêu (MOPSO) cho hệ 3×3 và 4×4, đảm bảo sự hài hòa giữa bám quỹ đạo, kháng nhiễu và an toàn bền vững.
  4. Chứng minh toán học và mô phỏng thực chứng tính ổn định bền vững của toàn bộ hệ thống kín trước sự bất định sai số nhân ngõ ra thông qua cấu trúc $M-\Delta$ và Giá trị suy biến có cấu trúc (SSV luôn thỏa mãn $\mu < 1$).
  5. Mở rộng thành công phương pháp nhận dạng bình phương cực tiểu (LS) qua mô tả phân rã ma trận (MFD) cho hệ đa biến và kiểm chứng thực nghiệm thời gian thực hoàn hảo trên hệ thống bốn bồn nước liên kết (Quadruple Tank).
  6. Mở ra 3 hướng nghiên cứu học thuật mới: Điều khiển phân số cho hệ đa biến không vuông, điều khiển phân số thích nghi phi tuyến thời gian thực, và tích hợp thuật toán phân số vào hệ thống vi điều khiển công nghiệp chuyên dụng.