Tổng quan về giáo trình

Giáo trình bài giảng Fundamentals of Control Systems (Cơ sở Lý thuyết Điều khiển Tự động), cụ thể tại Chapter 3: System Dynamics (Chương 3: Động học Hệ thống), do PGS. Huỳnh Thái Hoàng (Assoc. Prof. Huynh Thai Hoang) thuộc Bộ môn Điều khiển Tự động, Khoa Điện - Điện tử, Trường Đại học Bách Khoa – Đại học Quốc gia TP.HCM (Ho Chi Minh City University of Technology) biên soạn (phiên bản cập nhật ngày 06/12/2013). Tài liệu nằm trong khối kiến thức cơ sở ngành kỹ thuật then chốt, phục vụ đào tạo các ngành Kỹ thuật Điều khiển và Tự động hóa, Điện – Điện tử và các chuyên ngành kỹ thuật liên quan.

Mục tiêu học tập của phần tài liệu này tập trung vào việc trang bị cho người học:

  • Khái niệm nền tảng về động học của các hệ thống phức tạp theo thời gian thông qua mô hình toán học.
  • Phương pháp xác định và phân tích đáp ứng miền thời gian (Time response) với các tín hiệu thử tiêu chuẩn: hàm xung Dirac $\delta(t)$ và hàm nấc nhảy đơn vị $1(t)$.
  • Cơ sở lý thuyết và phương pháp đánh giá đáp ứng miền tần số (Frequency response) thông qua tỷ số giữa tín hiệu ra xác lập và tín hiệu vào dạng hình sin.
  • Kỹ năng xác định đặc tính động học của 7 khâu cơ bản và hệ thống điều khiển tổng quát thông qua các công cụ đồ thị chuẩn tắc: biểu đồ Bode (Bode diagram) và biểu đồ Nyquist (Nyquist plot).

Cấu trúc tài liệu được thiết kế theo phương pháp diễn dịch hình thức: xuất phát từ bản chất toán học của mô hình hàm truyền (Transfer Function - TF), chuyển đổi giữa miền phức Laplace ($s$), miền thời gian ($t$) và miền tần số liên tục ($j\omega$), từ đó thiết lập các quy tắc phân tích định lượng và biểu diễn đồ thị tiệm cận. Điểm đặc thù của tài liệu là tính cô đọng, quy chuẩn hóa các bước tính toán giải tích và thuật toán dựng biểu đồ.

Nội dung kiến thức cốt lõi

                             HỆ THỐNG ĐỘNG HỌC G(s)
  ĐÁP ỨNG THỜI GIAN                                          ĐÁP ỨNG TẦN SỐ
  (Time Response)                                         (Frequency Response)
Hàm trọng lượng  Hàm quá độ                             Biểu đồ Bode      Đồ thị Nyquist
g(t) = L⁻¹{G(s)} h(t) = L⁻¹{G(s)/s}                     L(ω) & φ(ω)       Tọa độ cực G(jω)
                                                             DỰ TRỮ ỔN ĐỊNH
                                                            GM & PM (Độ đo)

Các chương và chủ đề chính

Tài liệu tập trung triển khai toàn diện Chương 3 (System Dynamics) với 4 mảng chuyên đề mang tính tuần tự logic:

  1. Khái niệm về động học hệ thống (The concept of system dynamics): Định nghĩa động học hệ thống là việc nghiên cứu hành vi của hệ thống phức tạp biến thiên theo thời gian. Các hệ thống có mô hình toán học tương tự nhau sẽ biểu diễn các đáp ứng động học tương đương. Việc khảo sát được thực hiện qua các tín hiệu thử nghiệm đầu vào cơ bản: xung Dirac, nấc thang (step) và hàm sin.
  2. Đáp ứng miền thời gian (Time response):
    • Đáp ứng xung (Impulse response / Weighting function): Ký hiệu $g(t)$, là đáp ứng của hệ thống khi ngõ vào là xung Dirac ($U(s) = 1$). Theo đó, $Y(s) = G(s)$ và $g(t) = \mathcal{L}^{-1}{G(s)}$. Đáp ứng của hệ thống đối với một tín hiệu ngõ vào bất kỳ có thể được tính thông qua tích chập giữa hàm trọng lượng và ngõ vào.
    • Đáp ứng nấc (Step response / Transient function): Ký hiệu $h(t)$, là đáp ứng khi ngõ vào là hàm nấc đơn vị ($U(s) = 1/s$). Ta có $Y(s) = G(s)/s$ và $h(t) = \mathcal{L}^{-1}{G(s)/s} = \int_0^t g(\tau)d\tau$. Hàm quá độ là tích phân của hàm trọng lượng.
    • Ví dụ mẫu: Khảo sát hệ thống có hàm truyền $G(s) = \frac{s+1}{s(s+5)}$, giải tích tìm ra:
      • Hàm trọng lượng: $g(t) = \frac{1}{5} + \frac{4}{5}e^{-5t}$.
      • Hàm quá độ: $h(t) = \frac{1}{5}t - \frac{4}{25} + \frac{4}{25}e^{-5t}$.
  3. Đáp ứng miền tần số (Frequency response):
    • Khảo sát đáp ứng xác lập của hệ thống tuyến tính dưới tác động của tín hiệu sin $u(t) = U_m \sin(\omega t)$. Ngõ ra xác lập có dạng $y(t) = Y_m \sin(\omega t + \phi)$, cùng tần số $\omega$ nhưng khác biệt về biên độ và pha.
    • Định nghĩa: Đáp ứng tần số là tỷ số giữa ngõ ra xác lập và ngõ vào hình sin, chứng minh bằng $G(j\omega) = G(s)|_{s=j\omega} = M(\omega)e^{j\phi(\omega)}$.
    • Đáp ứng biên độ (Magnitude response): $M(\omega) = |G(j\omega)| = \sqrt{P^2(\omega) + Q^2(\omega)}$, thể hiện độ khuếch đại (gain) của hệ thống theo tần số.
    • Đáp ứng pha (Phase response): $\phi(\omega) = \arg G(j\omega) = \text{arctg}\left(\frac{Q(\omega)}{P(\omega)}\right)$, thể hiện độ lệch pha giữa tín hiệu ra và tín hiệu vào theo tần số.
    • Công cụ đồ thị:
      • Biểu đồ Bode gồm 2 đồ thị: Đồ thị biên độ $L(\omega) = 20\lg M(\omega)$ (đơn vị dB) theo trục hoành logarit và đồ thị pha $\phi(\omega)$.
      • Biểu đồ Nyquist: Đồ thị biểu diễn $G(j\omega)$ trong hệ tọa độ cực khi $\omega$ quét từ $0 \to +\infty$.
    • Tần số cắt và độ dự trữ ổn định:
      • Tần số cắt biên ($\omega_c$ - Gain crossover frequency): Tần số tại đó $M(\omega_c) = 1$ (tương đương $L(\omega_c) = 0\text{ dB}$).
      • Dự trữ biên ($GM$ - Gain Margin): $GM = \frac{1}{M(\omega_{-\pi})} \Leftrightarrow GM = -L(\omega_{-\pi})\text{ [dB]}$, thể hiện độ lợi bổ sung cần thiết tại tần số cắt pha để đưa hệ thống đến biên giới ổn định.
      • Dự trữ pha ($PM$ - Phase Margin): $PM = 180^\circ + \phi(\omega_c)$, là độ trễ pha bổ sung tại tần số cắt biên cần thiết để đưa hệ thống đến biên giới ổn định.
  4. Đặc tính động học của các khâu cơ bản (Dynamics of typical components):
Khâu động học Hàm truyền $G(s)$ Đáp ứng thời gian ($g(t), h(t)$) Đáp ứng tần số ($M(\omega), \phi(\omega)$) Tiệm cận biểu đồ Bode
Khâu khuếch đại (Proportional) $K$ $g(t) = K\delta(t)$
$h(t) = K \cdot 1(t)$		 $M(\omega) = K \Rightarrow L(\omega) = 20\lg K$
$\phi(\omega) = 0^\circ$		 Đường nằm ngang $20\lg K$
Khâu tích phân (Integral) $\frac{1}{s}$ $g(t) = 1(t)$
$h(t) = t \cdot 1(t)$		 $M(\omega) = \frac{1}{\omega} \Rightarrow L(\omega) = -20\lg\omega$
$\phi(\omega) = -90^\circ$		 Đường thẳng dốc $-20\text{ dB/dec}$ qua $(1, 0\text{ dB})$
Khâu vi phân (Derivative) $s$ $g(t) = \delta'(t)$
$h(t) = \delta(t)$		 $M(\omega) = \omega \Rightarrow L(\omega) = 20\lg\omega$
$\phi(\omega) = +90^\circ$		 Đường thẳng dốc $+20\text{ dB/dec}$ qua $(1, 0\text{ dB})$
Quán tính bậc nhất (1st-order lag) $\frac{1}{Ts+1}$ $g(t) = \frac{1}{T}e^{-\frac{t}{T}}1(t)$
$h(t) = (1-e^{-\frac{t}{T}})1(t)$		 $M(\omega) = \frac{1}{\sqrt{1+T^2\omega^2}}$
$\phi(\omega) = -\text{arctg}(T\omega)$		 $\omega < \frac{1}{T}$: trục hoành ($0\text{ dB}$)
$\omega > \frac{1}{T}$: dốc $-20\text{ dB/dec}$
Vi phân bậc nhất (1st-order lead) $Ts+1$ $g(t) = T\delta'(t) + \delta(t)$
$h(t) = T\delta(t) + 1(t)$		 $M(\omega) = \sqrt{1+T^2\omega^2}$
$\phi(\omega) = \text{arctg}(T\omega)$		 $\omega < \frac{1}{T}$: trục hoành ($0\text{ dB}$)
$\omega > \frac{1}{T}$: dốc $+20\text{ dB/dec}$
Dao động bậc hai (2nd-order oscillating) $\frac{1}{T^2s^2 + 2\xi Ts + 1}$
($0 < \xi < 1$)		 $g(t) = \frac{\omega_n}{\sqrt{1-\xi^2}}e^{-\xi\omega_n t}\sin(\omega_n\sqrt{1-\xi^2}t)$
$h(t) = 1 - \frac{e^{-\xi\omega_n t}}{\sqrt{1-\xi^2}}\sin(\dots)$		 $M(\omega) = \frac{1}{\sqrt{(1-T^2\omega^2)^2 + 4\xi^2 T^2\omega^2}}$
$\phi(\omega) = -\text{arctg}\left(\frac{2\xi T\omega}{1-T^2\omega^2}\right)$		 $\omega < \frac{1}{T}$: trục hoành ($0\text{ dB}$)
$\omega > \frac{1}{T}$: dốc $-40\text{ dB/dec}$
Khâu trễ (Time delay) $e^{-Ts}$ $g(t) = \delta(t-T)$
$h(t) = 1(t-T)$		 $M(\omega) = 1 \Rightarrow L(\omega) = 0\text{ dB}$
$\phi(\omega) = -\omega T$		 Biên độ không đổi ($0\text{ dB}$), pha giảm tuyến tính
  1. Động học của hệ thống điều khiển (Dynamics of control systems):
    • Xét hàm truyền phân thức hữu tỉ: $G(s) = \frac{b_0 s^m + b_1 s^{m-1} + \dots + b_m}{a_0 s^n + a_1 s^{n-1} + \dots + a_n}$.
    • Biến đổi Laplace của hàm quá độ: $H(s) = \frac{G(s)}{s} = \frac{b_0 s^m + \dots + b_m}{s(a_0 s^n + \dots + a_n)}$.
    • Các định lý giới hạn về đáp ứng thời gian:
      • Nếu $G(s)$ không chứa khâu tích phân/vi phân lý tưởng: $g(\infty) = \lim_{s\to 0} s G(s) = 0$; $h(\infty) = \lim_{s\to 0} s H(s) = \frac{b_m}{a_n} \neq 0$.
      • Nếu $G(s)$ chứa khâu tích phân lý tưởng ($a_n = 0$): $g(\infty) = \frac{b_m}{a_{n-1}} \neq 0$; $h(\infty) \to \infty$.
      • Nếu $G(s)$ chứa khâu vi phân lý tưởng ($b_m = 0$): $g(\infty) = 0$; $h(\infty) = 0$.
      • Tính chất strictly proper ($m < n$): $g(0) = \lim_{s\to\infty} s G(s) = 0$.
      • Tính chất proper ($m \le n$): $h(0) = \lim_{s\to\infty} s H(s) = 0$ (khi $m < n$).
    • Hệ thống gồm các khâu mắc nối tiếp $G(s) = \prod_{i=1}^l G_i(s)$:
      • $L(\omega) = \sum_{i=1}^l L_i(\omega)$ và $\phi(\omega) = \sum_{i=1}^l \phi_i(\omega)$.
      • Biểu đồ Bode tổng bằng tổng các biểu đồ Bode thành phần.
    • Quy trình 5 bước dựng biểu đồ Bode tiệm cận cho hàm truyền dạng chuẩn $G(s) = K s^\nu \prod G_i(s)$:
      • Bước 1: Xác định tất cả các tần số gãy $\omega_i = 1/T_i$, sắp xếp tăng dần $\omega_1 < \omega_2 < \omega_3 \dots$
      • Bước 2: Xác định điểm tựa $A$ có tọa độ $(\omega_0, L(\omega_0) = 20\lg K - 20\nu\lg\omega_0)$ với $\omega_0 < \omega_1$ (nếu $\omega_1 > 1$ có thể chọn $\omega_0 = 1$).
      • Bước 3: Qua điểm $A$, vẽ đường tiệm cận đầu tiên có độ dốc $-20\nu\text{ dB/dec}$ (nếu có khâu tích phân) hoặc $+20\nu\text{ dB/dec}$ (nếu có khâu vi phân), kéo dài đến tần số gãy đầu tiên.
      • Bước 4: Tại mỗi tần số gãy $\omega_i$, cộng thêm vào độ dốc hiện tại:
        • $-20\text{ dB/dec} \times \mu_i$ đối với khâu quán tính bậc nhất bội $\mu_i$.
        • $+20\text{ dB/dec} \times \mu_i$ đối với khâu vi phân bậc nhất bội $\mu_i$.
        • $-40\text{ dB/dec} \times \mu_i$ đối với khâu dao động bậc hai bội $\mu_i$.
        • $+40\text{ dB/dec} \times \mu_i$ đối với khâu vi phân bậc hai bội $\mu_i$.
      • Bước 5: Lặp lại Bước 4 cho đến khi qua hết tần số gãy cuối cùng.

Kiến thức nền tảng được xây dựng

Nội dung xây dựng cơ sở toán học chặt chẽ dựa trên:

  • Lý thuyết hàm phức và phép biến đổi Laplace thuận/nghịch.
  • Phân tích đa thức và nghiệm của phương trình đặc trưng.
  • Mối quan hệ đối ngẫu giữa phân tích quá độ miền thời gian và phổ tần số miền xác lập.

Kỹ năng phát triển

  • Kỹ năng phân tích giải tích đáp ứng động học của hệ thống thông qua phân tích phân thức tối giản.
  • Kỹ năng phác thảo đồ thị tiệm cận Bode và đồ thị Nyquist theo quy trình chuẩn.
  • Kỹ năng xác định và biện luận định lượng các chỉ số dự trữ ổn định ($GM, PM$) từ biểu đồ tần số.

Phương pháp giảng dạy và học tập

Tài liệu sử dụng phương pháp sư phạm cấu trúc hóa dạng mô-đun toán học kỹ thuật, kết hợp giữa chứng minh công thức và hình học trực quan:

Cách tiếp cận sư phạm

Bài giảng áp dụng tiến trình sư phạm 4 bước cho từng đối tượng khảo sát:

  1. Xác lập mô hình toán học (Hàm truyền $G(s)$).
  2. Khai triển giải tích miền thời gian ($g(t), h(t)$ qua biến đổi Laplace ngược).
  3. Khai triển giải tích miền tần số ($M(\omega), \phi(\omega)$ qua phép thế $s = j\omega$).
  4. Đồ thị hóa thông qua đường cong lý thuyết và các đường tiệm cận xấp xỉ.

Bài tập và ví dụ minh họa

Tài liệu tích hợp các ví dụ giải mẫu cụ thể từng bước tính toán. Điển hình là bài toán giải tích hàm phân thức: $$G(s) = \frac{s+1}{s(s+5)}$$ Các bước phân tích mẫu phân tách: $$\mathcal{L}^{-1}\left{\frac{s+1}{s(s+5)}\right} = \mathcal{L}^{-1}\left{\frac{1}{5s} + \frac{4}{5(s+5)}\right}$$ cho thấy rõ cách thức chuyển giao từ lý thuyết sang giải pháp kỹ thuật giải tích.

Hướng dẫn tự học

Người học được định hướng nắm vững quy tắc giải tích 5 bước để dựng nhanh biểu đồ Bode biên độ mà không phụ thuộc hoàn toàn vào phần mềm tính toán, đồng thời luyện tập phân tích các điểm kỳ dị, tần số cắt và độ lệch pha trên đồ thị Nyquist.

Điểm nổi bật và cập nhật

  • Hệ thống hóa toàn diện 7 khâu cơ bản: Phân loại và lập bảng đối sánh chi tiết toàn bộ các khâu động học kinh điển trong lý thuyết điều khiển (Khuếch đại, Tích phân, Vi phân, Quán tính bậc nhất, Vi phân bậc nhất, Dao động bậc hai, Trễ).
  • Quy chuẩn hóa 5 bước vẽ tiệm cận Bode: Cung cấp thuật toán rõ ràng, hỗ trợ việc phác thảo nhanh đường đặc tính biên độ tần số logarit của hệ nối tiếp bất kỳ thông qua các tần số gãy $\omega_i = 1/T_i$ và độ dốc chuẩn ($\pm 20\text{ dB/dec}, \pm 40\text{ dB/dec}$).
  • Khảo sát giải tích các định lý giá trị biên: Trình bày chi tiết các hệ quả toán học của định lý giá trị cuối và giá trị đầu đối với hàm quá độ $h(t)$ và hàm trọng lượng $g(t)$ ứng với từng bậc của đa thức tử ($m$) và đa thức mẫu ($n$), phân định rõ ràng giữa hệ proper ($m \le n$) và strictly proper ($m < n$).
  • Mối liên hệ hình học - giải tích: Kết nối trực tiếp giữa công thức giải tích $G(j\omega)$ với độ đo vật lý của dự trữ biên ($GM$) và dự trữ pha ($PM$).

Đối tượng sử dụng giáo trình

  • Sinh viên đại học: Dành cho sinh viên năm thứ 2 và năm thứ 3 thuộc các ngành Kỹ thuật Điều khiển và Tự động hóa, Điện - Điện tử, Kỹ thuật Cơ điện tử, Kỹ thuật Hệ thống Công nghiệp theo học học phần Lý thuyết Điều khiển Tự động / Cơ sở Điều khiển Tuyến tính.
  • Yêu cầu tiên quyết (Prerequisites):
    • Toán cao cấp (Giải tích vi tích phân, Đại số tuyến tính, Phương trình vi phân).
    • Phép biến đổi Laplace và ứng dụng giải phương trình vi phân tuyến tính.
    • Kiến thức cơ bản về biến đổi tín hiệu và lý thuyết mạch điện.
  • Giảng viên: Sử dụng làm khung bài giảng chuẩn cho các bài giảng lý thuyết trên lớp, định hướng cấu trúc ra đề thi tự luận và thiết kế bài tập lớn.
  • Kỹ sư và người tự học: Tài liệu tham khảo chuẩn mực để tra cứu công thức hàm truyền, đáp ứng bước nhảy, đáp ứng xung và phương pháp đánh giá ổn định hệ thống tuyến tính.

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

1. Giáo trình này phù hợp với ai?

Tài liệu được thiết kế phục vụ trực tiếp cho sinh viên đại học ngành Kỹ thuật Điều khiển và Tự động hóa, Điện – Điện tử, cùng các kỹ sư cần củng cố kiến thức nền tảng về mô hình hóa và phân tích động học hệ thống điều khiển liên tục.

2. Cần chuẩn bị kiến thức nền nào để học?

Người học cần nắm vững toán giải tích, phương trình vi phân tuyến tính, biến đổi Laplace, biến số phức và các khái niệm cơ bản về mạch điện hoặc cơ học tuyến tính.

3. Điểm khác biệt trong phương pháp trình bày của tài liệu này là gì?

Tài liệu trình bày dưới dạng bài giảng cô đọng, song song hóa phân tích giữa miền thời gian ($t$) và miền tần số ($j\omega$), chuẩn hóa quy trình 5 bước phác thảo biểu đồ Bode tiệm cận và cung cấp các hệ quả định lý giá trị biên tường minh.

4. Làm sao để tự học và áp dụng hiệu quả nội dung Chapter 3?

Người học nên tự giải lại các ví dụ mẫu phân tích biến đổi Laplace ngược của hàm phân thức, ghi nhớ bảng đặc tính 7 khâu cơ bản, sau đó thực hành từng bước dựng biểu đồ Bode tiệm cận từ việc xác định tần số gãy $\omega_i = 1/T_i$ đến tính độ dốc tổng hợp.

5. Có công cụ biểu diễn đồ thị nào được đề cập trong tài liệu?

Tài liệu giới thiệu và chuẩn hóa cách sử dụng hai công cụ phân tích tần số kinh điển là Biểu đồ Bode (gồm đồ thị biên độ logarit $L(\omega)$ và đồ thị pha $\phi(\omega)$) và Biểu đồ Nyquist (đồ thị cực của $G(j\omega)$).

Kết luận

Bài giảng Chapter 3: System Dynamics trong tập tài liệu Fundamentals of Control Systems của PGS. Huỳnh Thái Hoàng (Đại học Bách Khoa – ĐHQG-HCM) cung cấp khung kiến thức toán học và kỹ thuật chuẩn xác về động học hệ thống điều khiển tuyến tính liên tục.

Lộ trình học tập đề xuất: $$\text{Mô hình hóa } G(s) \longrightarrow \text{Khảo sát 7 khâu cơ bản} \longrightarrow \text{Phân tích đáp ứng thời gian } g(t), h(t) \longrightarrow \text{Dựng biểu đồ Bode/Nyquist} \longrightarrow \text{Đánh giá } GM, PM$$

Nội dung đóng vai trò cầu nối lý thuyết then chốt trước khi người học tiến tới các nội dung tiếp theo về đánh giá tính ổn định hệ thống (Chương 4) và thiết kế bộ điều khiển (Chương 5, 6).