Tổng quan luận án

Luận án tiến sĩ kỹ thuật chuyên ngành Kỹ thuật cơ khí (mã số: 9520103) của tác giả Võ Lâm Chương, do PGS. Trương Nguyễn Luân Vũ và GS. MoonYong Lee hướng dẫn khoa học, được thực hiện tại Trường Đại học Sư phạm Kỹ thuật Thành phố Hồ Chí Minh (tháng 06/2024). Đề tài tập trung nghiên cứu: "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".

Tính cấp thiết và khoảng trống nghiên cứu

Trong các ngành công nghiệp chế biến, hóa chất, lọc - hóa dầu và năng lượng, phần lớn các quá trình sản xuất là hệ thống đa biến phức tạp (MIMO), có sự tương tác chéo mạnh giữa các vòng điều khiển và tồn tại các khâu trễ thời gian vật lý đáng kể. Việc hiệu chỉnh độc lập từng vòng điều khiển đơn lẻ gặp nhiều hạn chế do tác động tương hỗ giữa các biến vào - ra. Dù các kỹ thuật điều khiển hiện đại như điều khiển dự báo mô hình (MPC) hay điều khiển mờ/nơ-ron đã phát triển, bộ điều khiển tỉ lệ - tích phân - vi phân (PID) truyền thống vẫn chiếm khoảng 90% ứng dụng công nghiệp nhờ cấu trúc đơn giản và độ tin cậy vận hành.

Tuy nhiên, các bộ điều khiển PID bậc nguyên cổ điển bị giới hạn về bậc tự do khi phải giải quyết đồng thời bài toán bám giá trị đặt (servomechanism), triệt tiêu nhiễu quá trình (regulator) và đảm bảo tính ổn định bền vững khi có sai số mô hình. Tính toán phân số (fractional calculus) với tích phân và đạo hàm bậc thực mở rộng bộ điều khiển PID thành bộ điều khiển PID bậc phân số (FOPID/FOPI), bổ sung thêm các bậc phân số $\lambda$ và $\mu$, tạo thêm độ linh hoạt khi thiết kế. Phần lớn các nghiên cứu về FOPID trước đây chỉ dừng lại ở hệ đơn biến (SISO) hoặc hệ hai ngõ vào - hai ngõ ra (TITO 2×2) với cấu trúc đơn giản. Khoảng trống nghiên cứu lớn là việc thiếu các cấu trúc điều khiển phân số toàn diện kết hợp phân ly hệ thống, giải quyết triệt để khâu trễ thời gian, tối ưu hóa giảm bậc mô hình phân ly cho hệ đa biến bậc cao ($3\times3$, $4\times4$), và kiểm chứng tính ổn định bền vững cũng như tính khả thi trên mô hình thực nghiệm.

Mục tiêu và nhiệm vụ nghiên cứu

Luận án xác định 4 nhiệm vụ nghiên cứu cụ thể được đánh số như sau:

  1. Tìm hiểu kỹ thuật phân ly cho hệ đa biến; đề xuất giải pháp cải tiến phương pháp tính toán để thuận tiện hơn trong việc tính toán thiết kế hệ phân ly (sử dụng thuật toán tối ưu hóa bầy đàn - PSO để rút gọn và xấp xỉ các hàm truyền thành phần).
  2. Đề xuất cấu trúc điều khiển mới cho hệ đa biến nhằm cải thiện đáp ứng của hệ khi giá trị đặt thay đổi và khi bị ảnh hưởng bởi nhiễu quá trình, kết hợp loại bỏ ảnh hưởng của đặc tính trễ bằng bộ dự báo Smith (SDSP/F-SDSP); đồng thời đánh giá ổn định bền vững của cấu trúc điều khiển đề xuất.
  3. Nghiên cứu bộ điều khiển PID phân số dựa trên nền tảng toán học tính toán phân số; đề xuất các quy luật hiệu chỉnh thông số mới cho hệ đa biến: phương pháp giải tích dựa trên cấu trúc mô hình nội (IMC) cho hệ $2\times2$ và phương pháp tối ưu hóa bầy đàn đa mục tiêu (MOPSO) cho hệ bậc cao ($3\times3$, $4\times4$).
  4. Xây dựng mô hình thực nghiệm hệ bốn bồn nước liên kết (quadruple tank); đề xuất phương pháp nhận dạng hệ đa biến mở rộng từ phương pháp bình phương tối thiểu (LS) kết hợp phân ly ma trận (MFD) để phục vụ thiết kế bộ điều khiển; kiểm chứng kết quả điều khiển thời gian thực trên phần mềm Matlab (Real-Time Windows Target).

Đối tượng và phạm vi nghiên cứu

  • Đối tượng nghiên cứu: Các quá trình công nghiệp đa biến vuông ($n \times n$) có trễ thời gian, kỹ thuật phân ly đơn giản hóa, bộ dự báo Smith, và các bộ điều khiển tỉ lệ - tích phân - vi phân bậc phân số (FOPI/FOPID).
  • Phạm vi nghiên cứu:
    • Về mặt lý thuyết và mô phỏng: Khảo sát các hệ đa biến vuông bậc thấp ($2\times2$) gồm mô hình Tháp chưng cất Vinante & Luyben (VL) và Bộ tách dầu nặng (Heavy Oil Fractionator); hệ đa biến bậc cao ($3\times3$ và $4\times4$) gồm Tháp chưng cất Ogunnaike & Ray (OR) và Hệ thống điều hòa nhiệt độ HVAC 4 phòng.
    • Về mặt thực nghiệm: Giới hạn trên hệ thực nghiệm bốn bồn nước liên kết (quadruple tank) tương đương hệ $2\times2$ tại phòng thí nghiệm do điều kiện kinh phí và trang thiết bị.

Tổng quan tài liệu và vị trí của luận án

Lịch sử nghiên cứu tính toán phân số và điều khiển phân số

Toán học về tính toán phân số khởi nguồn từ trao đổi giữa Leibniz và L'Hôpital vào năm 1695 về đạo hàm bậc $1/2$, tiếp tục phát triển qua các công trình của Abel (1823, bài toán tautochrone), Liouville, Riemann, Boole, Heaviside, Weyl, Hardy, Riesz. Trong lý thuyết điều khiển, Bode và Tustin (1945) là những người đầu tiên áp dụng đặc tính tần số bậc tổng quát. Từ thập niên 1990, Podlubny đã đề xuất bộ điều khiển $PI^\lambda D^\mu$ (FOPID), thúc đẩy các nghiên cứu về mô hình hóa và điều khiển phân số:

  • Nhóm nghiên cứu của Chen và cộng sự: Đóng góp các công trình nền tảng về FOPID, công cụ F-MIGO hiệu chỉnh bộ điều khiển FOPI cho hệ FOPDT (First Order Plus Delay Time), mở rộng ứng dụng trên phương tiện bay không người lái (UAV), hệ servo ổ cứng, khảo sát ổn định và nhận dạng hệ phi tuyến phân số.
  • Nhóm nghiên cứu CRONE (Oustaloup, Trigeassou, Malti): Phát triển CRONE Toolbox trên Matlab/Simulink, nghiên cứu tính ổn định bền vững, phương pháp nhận dạng sử dụng mô hình sai phân phân số, mô hình Hammerstein phân số và chuỗi Volterra phân số.
  • Các tác giả khác: Matignon và cộng sự nghiên cứu điều kiện ổn định của hệ tuyến tính phân số; các nghiên cứu mở rộng ổn định cho hệ phân số có trễ và hệ trạng thái rời rạc.

Tình hình nghiên cứu điều khiển hệ đa biến

Đối với hệ đa biến (MIMO), văn bản tổng hợp hai hướng tiếp cận chính:

  1. Điều khiển không tập trung (Decentralized/Multi-loop Control): Sử dụng các công cụ đo lường mức độ tương tác như dãy độ lợi tương đối RGA (Bristol, 1966), chỉ số Niederlinski (NI), RGA động (DRGA), RNGA. 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) cho hệ $n\times n$. Hướng này chỉ đạt hiệu quả khi mức độ tương tác chéo giữa các vòng lặp ở mức thấp.
  2. Điều khiển tập trung và điều khiển phân ly (Decoupling Control):
    • Điều khiển tập trung thuần nhất: Cần $n^2$ bộ điều khiển cho hệ $n\times n$, điển hình là MPC nhưng MPC thường chỉ dùng ở tầng điều khiển cấp cao do gánh nặng tính toán và trễ thời gian.
    • Điều khiển phân ly: Sử dụng ma trận phân ly để chuyển hệ MIMO thành các kênh SISO độc lập, sử dụng $n$ bộ điều khiển đơn biến. Ba kỹ thuật phân ly chính gồm: phân ly lý tưởng (gặp vấn đề về tính khả thi thực tế và nhạy cảm sai số), phân ly nghịch (Garrido và cộng sự, 2014 - cấu trúc ma trận phức tạp, khó mở rộng cho hệ bậc cao hơn $2\times2$) và phân ly đơn giản hóa (Vu và Lee, 2013 - ma trận phân ly có các phần tử đường chéo bằng 1, ổn định bền vững và khả thi cho hệ $n\times n$).
    • Bù thời gian trễ: Bộ dự báo Smith (SP/DTC) và các biến thể như bộ dự báo Smith có lọc (FSP) cho hệ SISO và MIMO. Các nghiên cứu kết hợp FSP với phân ly nghịch của Garrido và cộng sự chỉ phù hợp cho hệ $2\times2$.
Hướng nghiên cứu trước Tác giả tiêu biểu / Năm Hạn chế ghi nhận trong văn bản
Thiết kế FOPID cho hệ SISO Podlubny (1999), Chen và cộng sự Chưa giải quyết bài toán tương tác chéo của hệ đa biến công nghiệp.
Phân ly đơn giản hóa cho hệ $n\times n$ Vu và Lee (2013) Hàm truyền ma trận phân ly $D(s)$ và hệ sau phân ly $Q(s)$ có bậc rất cao, tính toán giải tích phức tạp.
Phân ly nghịch kết hợp IMC/FOPID Garrido và cộng sự (2014), các công bố 2018-2020 Khó xác định cấu trúc ma trận phân ly khi hệ có bậc cao hơn $2\times2$; phương pháp thiết kế phụ thuộc cấu trúc feedforward.
FOPID dùng thuật toán tiến hóa Nhóm nghiên cứu dùng GA, Bat Algorithm Cấu trúc bộ điều khiển phức tạp ($n^2$ bộ điều khiển), chỉ thử nghiệm mô phỏng trên hệ $2\times2$, thời gian tính toán tăng mạnh khi bậc hệ thống tăng.

Khoảng trống luận án giải quyết

Luận án giải quyết đồng thời các vấn đề tồn tại: kết hợp kỹ thuật phân ly đơn giản hóa với bộ dự báo Smith (F-SDSP) cho hệ đa biến có trễ; ứng dụng giải thuật PSO để rút gọn và xấp xỉ các hàm truyền thành phần có bậc nguyên/phân số; xây dựng quy tắc chỉnh định FOPID bằng IMC cho hệ $2\times2$ và MOPSO (dựa trên hàm mục tiêu kép và tiêu chuẩn $M_s$) cho hệ $3\times3$, $4\times4$; phân tích ổn định bền vững theo cấu trúc $M-\Delta$ với giá trị suy biến có cấu trúc (SSV); và tiến hành nhận dạng tham số thực nghiệm bằng phương pháp LS mở rộng kết hợp MFD trên hệ bốn bồn nước liên kết.

Cơ sở lý thuyết và phương pháp nghiên cứu

Cơ sở lý thuyết tính toán phân số

Luận án sử dụng định nghĩa tích phân và đạo hàm bậc phân số của Riemann-Liouville (R-L):

  • Tích phân phân số bậc $\alpha > 0$: $$_0D_t^{-\alpha} f(t) = \frac{1}{\Gamma(\alpha)} \int_0^t (t-\tau)^{\alpha-1} f(\tau) d\tau$$
  • Đạo hàm phân số bậc $\alpha > 0$ ($n-1 < \alpha < n$): $$_0D_t^\alpha f(t) = \frac{d^n}{dt^n} \left[ _0D_t^{-(n-\alpha)} f(t) \right] = \frac{1}{\Gamma(n-\alpha)} \frac{d^n}{dt^n} \int_0^t (t-\tau)^{n-\alpha-1} f(\tau) d\tau$$
  • Biến đổi Laplace với điều kiện đầu triệt tiêu: $$\mathcal{L}\left{ _0D_t^\alpha f(t) \right} = s^\alpha F(s), \quad \mathcal{L}\left{ _0D_t^{-\alpha} f(t) \right} = s^{-\alpha} F(s)$$
  • Phương trình sai phân vi phân tuyến tính bậc phân số (FODE): $$\sum_{i=0}^n a_i , 0D_t^{\alpha_i} y(t) = \sum{j=0}^m b_j , _0D_t^{\beta_j} u(t)$$ Vai trò mô tả động lực học của bậc phân số được làm rõ qua phương trình Langevin phân số (mô tả chuyển động Brownian có tính đến ma sát nhớt do gia tốc) và phương trình Van der Pol phân số (FO-VdP, biểu diễn sự linh hoạt của đặc tính pha tự dao động khi thay đổi bậc đạo hàm $\alpha$).

Cấu trúc bộ điều khiển PID bậc phân số (FOPID)

Cấu trúc hàm truyền của bộ điều khiển FOPID ($PI^\lambda D^\mu$) có dạng: $$G_c(s) = K_c \left( 1 + \frac{1}{\tau_I s^\lambda} + \tau_D s^\mu \right) = K_c + \frac{K_I}{s^\lambda} + K_D s^\mu$$ Trong đó: $K_c$ là hệ số tỉ lệ, $K_I$ ($\tau_I$) là hệ số (thời gian) tích phân, $K_D$ ($\tau_D$) là hệ số (thời gian) vi phân, $\lambda > 0$ là bậc tích phân, $\mu > 0$ là bậc vi phân.

Kỹ thuật phân ly đơn giản hóa và Bộ dự báo Smith

  • Phân ly đơn giản hóa: Hệ đa biến có ma trận hàm truyền $G(s)$ kích thước $n\times n$. 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)$ ($i \neq j$) được thiết kế để triệt tiêu tương tác chéo, biến đổi đối tượng thành ma trận đường chéo $Q(s) = G(s)D(s) = \text{diag}(q_{11}(s), q_{22}(s), \dots, q_{nn}(s))$.
  • Bộ dự báo Smith đa biến: Loại bỏ khâu trễ $\theta$ ra khỏi vòng hồi tiếp kín, cho phép thiết kế bộ điều khiển cho mô hình không trễ $Q_o(s)$ thay vì $Q(s)$.

Phương pháp nhận dạng hệ thống

  • Mở rộng phương pháp bình phương tối thiểu (LS) cho hệ MIMO thông qua mô tả phân ly ma trận (MFD - Matrix Fraction Description), chuyển đổi bài toán nhận dạng hệ MIMO thành các bài toán MISO độc lập.
  • Sử dụng tín hiệu đầu vào kích thích và đo lường ngõ ra trên mô hình phần cứng bốn bồn nước để xây dựng mô hình hàm truyền $2\times2$.

Tiêu chí đánh giá chất lượng và ổn định bền vững

  • Chỉ số đánh giá sai số:
    • Tích phân trị tuyệt đối sai số: $\text{IAE} = \int_0^\infty |e(t)| dt$
    • Tích phân sai số theo trọng số thời gian: $\text{ITAE} = \int_0^\infty t |e(t)| dt$
    • Tổng biến thiên tín hiệu điều khiển: $\text{TV} = \sum_{k=1}^\infty |u(k+1) - u(k)|$
  • Hàm độ nhạy cực đại ($M_s$): $$M_s = \max_\omega |S(j\omega)| = \max_\omega \left| \frac{1}{1 + G(j\omega)G_c(j\omega)} \right|$$
  • Khung phân tích ổn định bền vững: Sử dụng cấu trúc $M-\Delta$ với mô hình bất định sai số nhân ngõ ra (multiplicative output uncertainty). Ổn định bền vững được đánh giá thông qua giá trị suy biến có cấu trúc (SSV, ký hiệu $\mu(M)$ hay $\mu$-synthesis). Điều kiện ổn định bền vững là: $$\mu(M(j\omega)) < 1, \quad \forall \omega$$

Nội dung chính theo từng chương

Chương 1. Tổng quan

Chương 1 giới thiệu bối cảnh nghiên cứu, tính cấp thiết của việc áp dụng tính toán phân số vào lĩnh vực điều khiển hệ đa biến công nghiệp. Tác giả tổng quan chi tiết các hướng nghiên cứu về tính toán phân số, các công trình của Podlubny, Chen, Oustaloup (nhóm CRONE) và các phương pháp điều khiển hệ đa biến (điều khiển phi tập trung RGA/RNGA, điều khiển tập trung MPC, phân ly lý tưởng, phân ly nghịch, phân ly đơn giản hóa). Chương này xác định các rào cản tính toán trong việc rút gọn hàm truyền phân ly và sự thiếu vắng các quy tắc hiệu chỉnh FOPID cho hệ đa biến bậc cao, từ đó xác lập mục tiêu, phạm vi và phương pháp luận của luận án.

Chương 2. Cơ sở lý thuyết

Chương 2 trình bày có hệ thống các nền tảng toán học và kỹ thuật điều khiển:

  • Định nghĩa Riemann-Liouville, biến đổi Laplace cho toán tử phân số, phương trình vi phân sai phân phân số (FODE).
  • Minh chứng tính đa dạng và linh hoạt của bậc phân số thông qua phương trình chuyển động Langevin phân số và phương trình tự dao động Van der Pol phân số (FO-VdP).
  • Khảo sát đặc tính đáp ứng bậc thang của khâu tích phân phân số $s^{-\lambda}$ và vi phân phân số $s^\mu$, xấp xỉ Oustaloup trong miền tần số.
  • Trình bày kỹ thuật phân ly đơn giản hóa cho hệ đa biến, cấu trúc bộ dự báo Smith đơn biến và đa biến.
  • Cơ sở lý thuyết nhận dạng hệ thống tuyến tính rời rạc bằng phương pháp bình phương tối thiểu (LS) và mở rộng cho hệ đa biến qua mô hình MFD.

Chương 3. Đánh giá hệ thống và ổn định bền vững

Chương 3 xây dựng khung đánh giá định lượng cho các thuật toán điều khiển:

  • Định nghĩa các chỉ số chất lượng miền thời gian: IAE, ITAE và tổng biến thiên tín hiệu điều khiển TV.
  • Phân tích ý nghĩa hình học của hàm độ nhạy $S(s)$ và hàm độ nhạy bổ sung $T(s)$, xác lập giới hạn của hàm độ nhạy cực đại $M_s$ để đảm bảo độ dự trữ ổn định.
  • Thiết lập mô hình bất định sai số nhân ở ngõ ra và chuyển đổi hệ kín về cấu trúc chuẩn $M-\Delta$.
  • Trình bày định lý độ lợi nhỏ (Small Gain Theorem), điều kiện ổn định nội, và công cụ giá trị suy biến có cấu trúc (SSV - $\mu$) để kiểm chứng tính ổn định bền vững của hệ đa biến có bậc nguyên và bậc phân số.

Chương 4. Các phương pháp thiết kế đề xuất

Chương này là trọng tâm đóng góp lý thuyết và thuật toán của luận án:

  • Cấu trúc điều khiển F-SDSP: Kết hợp bộ phân ly đơn giản hóa $D(s)$ và bộ dự báo Smith đa biến, triệt tiêu tương tác chéo và loại bỏ khâu trễ ra khỏi vòng hồi tiếp kín.
  • Rút gọn mô hình bằng PSO: Sử dụng giải thuật tối ưu hóa bầy đàn để xấp xỉ các phần tử hàm truyền phức tạp của $D(s)$ và $Q(s)$ về dạng hàm bậc nhất có trễ (FOPDT) hoặc bậc hai có trễ (SOPDT) có bậc nguyên hoặc bậc phân số, giảm thiểu sai số tính toán giải tích.
  • Quy tắc chỉnh định cho hệ $2\times2$ dựa trên IMC: Thiết lập cấu trúc điều khiển mô hình nội, suy dẫn thông số bộ điều khiển FOPI/FOPID thông qua việc lựa chọn thời hằng bộ lọc mong muốn $\tau_c$ để cân bằng giữa bám tín hiệu đặt và kháng nhiễu.
  • Quy tắc chỉnh định cho hệ bậc cao ($3\times3$, $4\times4$) bằng MOPSO:
    • Thiết lập bài toán tối ưu hóa đa mục tiêu với 2 hàm mục tiêu: cực tiểu hóa sai số bám điểm đặt ($J_r = \text{ITAE}_r$) và cực tiểu hóa sai số khi có nhiễu tác động ($J_d = \text{ITAE}_d$).
    • Thuật toán MOPSO tìm kiếm tập nghiệm tối ưu Pareto.
    • Lựa chọn điểm tối ưu cuối cùng trên đường biên Pareto thông qua ràng buộc hàm độ nhạy cực đại $M_s$ quy định cho từng vòng lặp.

Chương 5. Kết quả mô phỏng và thực nghiệm

Chương 5 trình bày toàn bộ kết quả kiểm chứng trên các mô hình chuẩn và hệ thống phần cứng thực tế:

1. Mô phỏng hệ TITO ($2\times2$)

  • Tháp chưng cất Vinante & Luyben (VL): So sánh F-SDSP với các phương pháp của Luyben, Vu & Lee, Garrido. Kết quả chỉ ra phương pháp đề xuất cải thiện đáng kể các chỉ số IAE, ITAE và TV; biểu đồ SSV ($\mu < 1$) khẳng định tính ổn định bền vững.
  • Bộ tách dầu nặng (Heavy Oil Fractionator): Mô hình có đặc tính trễ lớn và tương tác mạnh. Bộ điều khiển đề xuất duy trì đáp ứng quá độ ổn định, triệt tiêu can nhiễu giữa 2 kênh và tối ưu hóa năng lượng điều khiển thể hiện qua chỉ số TV.

2. Mô phỏng hệ đa biến bậc cao ($3\times3$ và $4\times4$)

  • Tháp chưng cất Ogunnaike & Ray (OR - $3\times3$): Xấp xỉ $q_{11}, q_{22}, q_{33}$ bằng PSO, áp dụng MOPSO tìm tập nghiệm Pareto theo $J_r$ và $J_d$. Hệ thống đáp ứng độc lập tại 3 vòng điều khiển khi thay đổi giá trị đặt và chịu nhiễu tải. Đồ thị SSV nằm hoàn toàn dưới ngưỡng 1.
  • Hệ thống điều hòa không khí HVAC 4 phòng ($4\times4$): Điều khiển nhiệt độ độc lập cho 4 vùng không gian tương tác nhiệt phức tạp. Cả 4 vòng điều khiển đạt được chất lượng bám điểm đặt nhanh, không vọt lố và đảm bảo ổn định bền vững.

3. Thực nghiệm trên hệ bốn bồn nước liên kết (Quadruple Tank)

  • Thiết lập phần cứng, thu thập dữ liệu vào - ra (điện áp bơm - mức nước $h_1, h_2$).
  • Nhận dạng hệ thống bằng thuật toán LS mở rộng với MFD, thu được ma trận hàm truyền $2\times2$.
  • Dùng PSO xấp xỉ hàm truyền phân ly $q_{11}(s), q_{22}(s)$ về dạng bậc phân số.
  • Thiết kế bộ điều khiển FOPID và nhúng vào mô hình Simulink chạy thời gian thực (Real-Time Windows Target).
  • Kết quả thực nghiệm khẳng định mức chất lỏng bám chính xác giá trị đặt, triệt tiêu tương tác giữa hai bồn và hệ thống vận hành ổn định.
Đối tượng nghiên cứu Kích thước Phương pháp đề xuất áp dụng Phương pháp so sánh đối chứng Tiêu chí kiểm chứng
Tháp Vinante & Luyben (VL) $2\times2$ F-SDSP + IMC FOPI/FOPID Luyben, Vu - Lee, Garrido IAE, ITAE, TV, biểu đồ SSV ($\mu$)
Bộ tách dầu nặng (Heavy Oil) $2\times2$ F-SDSP + IMC FOPI/FOPID Các phương pháp kinh điển đã công bố IAE, ITAE, TV, biểu đồ SSV ($\mu$)
Tháp Ogunnaike & Ray (OR) $3\times3$ Phân ly + PSO rút gọn + MOPSO FOPID Các bộ điều khiển đa biến công bố trước ITAE bám ($J_r$), ITAE nhiễu ($J_d$), $M_s$, SSV
Hệ thống HVAC 4 phòng $4\times4$ Phân ly + PSO rút gọn + MOPSO FOPID Cấu trúc phi tập trung / phân ly truyền thống Đáp ứng 4 kênh nhiệt độ, TV, SSV
Hệ 4 bồn nước liên kết $2\times2$ Nhận dạng MFD-LS + PSO + FOPID thực nghiệm Bộ điều khiển PID bậc nguyên trên phần cứng thật Mức nước $h_1, h_2$, điện áp điều khiển thực tế

Chương 6. Kết luận

Tổng kết toàn bộ các kết quả lý thuyết, mô phỏng và thực nghiệm đã đạt được; nêu rõ các giới hạn còn tồn tại về quy mô hệ thống và định hướng các nghiên cứu mở rộng tiếp theo.

Kết quả và những đóng góp mới

Luận án tự tổng kết các đóng góp khoa học và thực tiễn chính như sau:

  1. Đề xuất cấu trúc điều khiển phân ly F-SDSP: Kết hợp kỹ thuật phân ly đơn giản hóa và bộ dự báo Smith cho hệ đa biến có trễ, giải quyết đồng thời bài toán tương tác chéo giữa các vòng điều khiển và loại bỏ ảnh hưởng của thời gian trễ ra khỏi đặc tính động học vòng kín.
  2. Ứng dụng giải thuật PSO trong rút gọn và xấp xỉ mô hình: Đề xuất quy trình sử dụng thuật toán tối ưu hóa bầy đàn để xấp xỉ các phần tử của ma trận phân ly $D(s)$ và ma trận sau phân ly $Q(s)$ về các dạng chuẩn tắc (bậc nguyên hoặc phân số), giải quyết triệt để sự phức tạp trong tính toán giải tích khi bậc hệ thống tăng cao.
  3. Phát triển các quy luật hiệu chỉnh thông số FOPID đa biến:
    • Đối với hệ bậc thấp ($2\times2$): Xây dựng quy luật giải tích dựa trên cấu trúc mô hình nội (IMC), giảm số lượng tham số cần tinh chỉnh về thời hằng bộ lọc mong muốn $\tau_c$.
    • Đối với hệ bậc cao ($3\times3$ và $4\times4$): Ứng dụng giải thuật MOPSO tối ưu hóa đồng thời chất lượng bám điểm đặt ($J_r$) và kháng nhiễu ($J_d$), sử dụng biên Pareto kết hợp ràng buộc độ nhạy cực đại $M_s$ để xác định bộ thông số điều khiển tối ưu.
  4. Mở rộng phân tích ổn định bền vững cho hệ phân số đa biến: Áp dụng cấu trúc $M-\Delta$ với mô hình sai số nhân ngõ ra và công cụ giá trị suy biến có cấu trúc (SSV / $\mu$) để kiểm chứng định lượng tính ổn định bền vững của hệ thống điều khiển phân số đa biến khi có sai lệch mô hình.
  5. Đóng góp thực nghiệm và phương pháp nhận dạng: Mở rộng phương pháp bình phương tối thiểu (LS) kết hợp mô tả phân ly ma trận (MFD) để nhận dạng hệ MIMO; chế tạo và kiểm chứng thành công bộ điều khiển FOPID thời gian thực trên mô hình phần cứng bốn bồn nước liên kết qua Matlab Real-Time Windows Target.

Hạn chế và hướng nghiên cứu tiếp

Theo công bố của tác giả trong văn bản, các hạn chế và hướng phát triển được xác định bao gồm:

  • Hạn chế về quy mô hệ thống: Nghiên cứu lý thuyết và thuật toán giới hạn trên các hệ đa biến vuông ($n\times n$); phần mô phỏng dừng lại ở hệ $4\times4$ và phần thực nghiệm chỉ thực hiện trên hệ $2\times2$ do hạn chế về kinh phí và thiết bị chế tạo.
  • Hạn chế về cấu trúc mô hình nhận dạng: Quá trình nhận dạng thực nghiệm sử dụng mô hình tuyến tính quanh điểm làm việc, chưa xét đến tính phi tuyến mạnh khi dải hoạt động thay đổi rộng.
  • Hướng nghiên cứu tiếp theo:
    • Mở rộng thuật toán phân ly và điều khiển phân số cho các hệ thống đa biến không vuông ($m \times n$).
    • Nghiên cứu các giải thuật nhận dạng trực tuyến và điều khiển thích nghi bậc phân số cho các đối tượng có tham số biến đổi theo thời gian.
    • Ứng dụng cấu trúc F-SDSP và MOPSO FOPID trên các quy trình công nghiệp quy mô lớn hơn trong thực tế.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên sâu trong lĩnh vực Kỹ thuật điều khiển và Tự động hóa, Kỹ thuật Cơ điện tử và Kỹ thuật Cơ khí:

  • Đối với nghiên cứu sinh và học viên cao học: Cung cấp tài liệu hệ thống về toán học tính toán phân số (toán tử Riemann-Liouville, biến đổi Laplace, hàm Mittag-Leffler, phương trình vi phân FODE), cấu trúc bộ điều khiển FOPID và phương pháp phân tích ổn định bền vững $M-\Delta$ / SSV.
  • Đối với giảng viên và nhà nghiên cứu: Cung cấp phương pháp luận kết hợp giữa giải thuật tiến hóa (PSO, MOPSO) và các cấu trúc điều khiển kinh điển (IMC, Smith Predictor, Decoupling) để giải quyết bài toán điều khiển quá trình đa biến có trễ.
  • Đối với kỹ sư điều khiển quá trình: Cung cấp quy trình nhận dạng thực tế hệ đa biến bằng phương pháp MFD-LS và giải pháp triển khai thuật toán điều khiển thời gian thực trên phần mềm Matlab Real-Time Windows Target kết nối phần cứng.

Câu hỏi thường gặp

1. Cấu trúc F-SDSP trong luận án giải quyết vấn đề gì của hệ đa biến có trễ?

Cấu trúc F-SDSP (Simplified Decoupling Smith Predictor) kết hợp kỹ thuật phân ly đơn giản hóa với bộ dự báo Smith. Ma trận phân ly đơn giản hóa biến đổi ma trận truyền dẫn đầy đủ của hệ MIMO thành ma trận đường chéo $Q(s)$ (loại bỏ tương tác chéo giữa các vòng), trong khi bộ dự báo Smith đưa các khâu trễ thời gian ra ngoài vòng hồi tiếp kín, cho phép thiết kế các bộ điều khiển đơn biến độc lập trên các kênh không trễ.

2. Thuật toán PSO được sử dụng vào khâu nào trong quy trình thiết kế của tác giả?

Thuật toán PSO (Particle Swarm Optimization) đơn mục tiêu được tác giả sử dụng vào khâu rút gọn và xấp xỉ mô hình: thay vì giải tích toán học phức tạp để tìm hàm truyền của ma trận phân ly $D(s)$ và ma trận hệ sau phân ly $Q(s)$, PSO tối ưu hóa sai lệch để đưa các hàm truyền bậc cao về dạng chuẩn FOPDT hoặc SOPDT có bậc nguyên hoặc bậc phân số.

3. Phương pháp chỉnh định bộ điều khiển FOPID cho hệ $3\times3$ và $4\times4$ khác hệ $2\times2$ như thế nào?

Đối với hệ $2\times2$, tác giả sử dụng phương pháp giải tích dựa trên cấu trúc mô hình nội (IMC) với tham số điều chỉnh là thời hằng $\tau_c$. Đối với hệ bậc cao ($3\times3$ và $4\times4$), do tính toán giải tích không còn khả thi, tác giả sử dụng thuật toán tối ưu hóa bầy đàn đa mục tiêu (MOPSO) tối ưu hóa đồng thời hai hàm mục tiêu $J_r$ (bám điểm đặt) và $J_d$ (kháng nhiễu), sau đó chọn điểm nghiệm trên đường biên Pareto thỏa mãn ràng buộc hàm độ nhạy cực đại $M_s$.

4. Tính ổn định bền vững của các bộ điều khiển đề xuất được kiểm chứng bằng công cụ nào?

Tính ổn định bền vững được kiểm chứng thông qua cấu trúc bất định sai số nhân ngõ ra đưa về dạng chuẩn $M-\Delta$. Đại lượng đánh giá là giá trị suy biến có cấu trúc (SSV, ký hiệu $\mu(M)$). Hệ thống được chứng minh đạt ổn định bền vững khi đường đặc tính biên độ của SSV thỏa mãn điều kiện $\mu(M(j\omega)) < 1$ trên toàn bộ dải tần số khảo sát.

5. Quá trình thực nghiệm trên hệ bốn bồn nước được tiến hành như thế nào?

Tác giả thu thập dữ liệu đáp ứng mức nước theo điện áp bơm, sử dụng kỹ thuật phân ly ma trận (MFD) kết hợp phương pháp bình phương tối thiểu (LS) để nhận dạng mô hình hàm truyền $2\times2$. Sau đó, PSO được áp dụng để xấp xỉ hàm truyền phân ly về dạng bậc phân số, tính toán các thông số bộ điều khiển FOPID và nạp vào sơ đồ Simulink điều khiển trực tiếp hệ thống phần cứng ở chế độ thời gian thực thông qua Matlab Real-Time Windows Target.

Kết luận

Luận án tiến sĩ của tác giả Võ Lâm Chương đã giải quyết trọn vẹn bài toán mô hình hóa và điều khiển phân số cho các quá trình đa biến vuông có trễ thời gian. Bằng việc kết hợp cấu trúc phân ly đơn giản hóa, bộ dự báo Smith, giải thuật tiến hóa (PSO, MOPSO) và kỹ thuật điều khiển mô hình nội (IMC), luận án đã thiết lập các quy tắc chỉnh định tham số FOPID hiệu quả cho cả hệ bậc thấp ($2\times2$) lẫn hệ bậc cao ($3\times3$, $4\times4$). Tính đúng đắn, chất lượng đáp ứng động học và độ ổn định bền vững của các phương pháp đề xuất đã được chứng minh chặt chẽ qua các mô hình mô phỏng chuẩn (Vinante & Luyben, Heavy Oil Fractionator, Ogunnaike & Ray, HVAC) và hệ thống thực nghiệm bốn bồn nước liên kết thời gian thực.