Chương 3: Động Học Hệ Thống - Bài Giảng Lý Thuyết Điều Khiển Tự Động

• Đề tài: Fundctrlsys Chapter 3

Tóm tắt

I. Khám phá nguyên tắc cốt lõi của động học hệ thống 2024

Chương 3 về động học hệ thống (System Dynamics) là một phần nền tảng trong lĩnh vực lý thuyết điều khiển tự động. Nội dung này tập trung vào việc nghiên cứu và thấu hiểu hành vi của các hệ thống phức tạp theo thời gian. Một nguyên tắc cơ bản được nêu rõ: các hệ thống được mô tả bởi những mô hình toán học tương tự sẽ biểu hiện các đáp ứng động học tương tự. Để phân tích các đáp ứng này, các chuyên gia thường sử dụng các tín hiệu đầu vào cơ bản như tín hiệu xung Dirac (impulse), tín hiệu nấc (step), hoặc tín hiệu hình sin. Việc phân tích này được chia thành hai lĩnh vực chính: đáp ứng thời gian (time response), xem xét hành vi của hệ thống trong miền thời gian, và đáp ứng tần số (frequency response), đánh giá cách hệ thống phản ứng với các tần số khác nhau. Hiểu rõ các khái niệm này là chìa khóa để thiết kế và tối ưu hóa các hệ thống điều khiển trong thực tế, từ kỹ thuật cơ khí, điện tử đến các quy trình công nghiệp phức tạp. Tài liệu của PGS. Huỳnh Thái Hoàng đã hệ thống hóa kiến thức này một cách bài bản.

1.1. Định nghĩa động học hệ thống System Dynamics là gì

Động học hệ thống là một phương pháp luận và mô phỏng toán học để định hình, thấu hiểu và thảo luận các vấn đề phức tạp. Trọng tâm của nó là nghiên cứu để hiểu rõ hành vi của các hệ thống theo thời gian. Điều này bao gồm việc phân tích mối quan hệ nhân quả, các vòng lặp phản hồi và độ trễ thời gian ảnh hưởng đến toàn bộ hệ thống. Trong lý thuyết điều khiển, động học hệ thống tập trung vào việc mô tả sự thay đổi trạng thái của hệ thống khi có tín hiệu tác động từ bên ngoài. Các mô hình toán học, đặc biệt là hàm truyền đạt (transfer function) trong miền Laplace, là công cụ chính để biểu diễn các đặc tính động học này. Một hệ thống được xem là ổn định nếu đầu ra của nó hội tụ về một giá trị hữu hạn khi thời gian tiến đến vô cùng. Ngược lại, hệ thống mất ổn định khi đầu ra tăng vô hạn. Phân tích động học giúp dự đoán các hành vi này trước khi triển khai hệ thống trong thực tế, cho phép các kỹ sư hiệu chỉnh thiết kế để đạt được hiệu suất mong muốn và đảm bảo an toàn vận hành.

1.2. Vai trò của tín hiệu vào cơ bản trong phân tích hệ thống

Để nghiên cứu đáp ứng động học của một hệ thống, việc lựa chọn tín hiệu đầu vào (input signal) phù hợp là cực kỳ quan trọng. Các tín hiệu cơ bản thường được sử dụng vì tính đơn giản và khả năng cung cấp thông tin toàn diện về đặc tính của hệ thống. Ba tín hiệu phổ biến nhất bao gồm:

1. Tín hiệu xung Dirac (Dirac impulse signal): Tín hiệu này có biên độ vô cùng lớn trong một khoảng thời gian vô cùng nhỏ, nhưng diện tích dưới đồ thị bằng 1. Đáp ứng của hệ thống với tín hiệu này, gọi là đáp ứng xung, chính là biến đổi Laplace ngược của hàm truyền đạt G(s).
2. Tín hiệu nấc (Step signal): Đây là tín hiệu thay đổi đột ngột từ 0 lên 1 tại thời điểm t=0 và giữ nguyên giá trị đó. Đáp ứng nấc cho thấy hệ thống chuyển từ trạng thái nghỉ sang trạng thái xác lập mới như thế nào, qua đó thể hiện các đặc tính như thời gian quá độ, độ vọt lố.
3. Tín hiệu hình sin (Sinusoidal signal): Tín hiệu này được dùng để phân tích đáp ứng tần số, cho biết hệ thống khuếch đại hay làm suy giảm biên độ và gây ra độ lệch pha như thế nào ở các tần số khác nhau. Các tín hiệu này đóng vai trò như những "bài kiểm tra tiêu chuẩn" giúp bóc tách và hiểu rõ từng khía cạnh trong hành vi của một hệ thống tuyến tính.

II. Hướng dẫn phân tích đáp ứng thời gian và đáp ứng tần số

Phân tích đáp ứng thời gian và đáp ứng tần số là hai phương pháp trụ cột để đánh giá động học hệ thống. Đáp ứng thời gian mô tả hành vi của tín hiệu đầu ra y(t) khi có một tín hiệu đầu vào u(t) cụ thể, thường là tín hiệu xung hoặc nấc. Thông qua đáp ứng xung g(t), ta có thể xác định được hàm trọng lượng của hệ thống. Trong khi đó, đáp ứng nấc h(t), hay còn gọi là hàm quá độ, cho thấy quá trình hệ thống đạt đến trạng thái ổn định mới. Ngược lại, đáp ứng tần số cung cấp một góc nhìn khác, tập trung vào cách hệ thống phản ứng ở trạng thái ổn định khi đầu vào là một tín hiệu hình sin. Nó được định nghĩa là tỷ số giữa tín hiệu ra và tín hiệu vào ở trạng thái ổn định, và được chứng minh bằng chính hàm truyền đạt G(s) khi thay s bằng jω. Phân tích này bao gồm đáp ứng biên độ (Magnitude response) và đáp ứng pha (Phase response), giúp đánh giá độ khuếch đại và độ trễ pha của hệ thống theo từng tần số.

2.1. Phân tích đáp ứng xung Impulse Response và ý nghĩa

Đáp ứng xung, ký hiệu là g(t), là hành vi của một hệ thống khi chịu tác động của tín hiệu xung Dirac. Về mặt toán học, khi tín hiệu vào U(s) = 1 (ảnh Laplace của xung Dirac), tín hiệu ra Y(s) sẽ bằng chính hàm truyền đạt G(s). Do đó, đáp ứng xung chính là biến đổi Laplace ngược của hàm truyền đạt: g(t) = L⁻¹{G(s)}. Điều này mang một ý nghĩa vật lý quan trọng: g(t) đại diện cho "bộ nhớ" hoặc "trọng số" của hệ thống tại mỗi thời điểm. Nó cho thấy cách hệ thống "nhớ" và phản ứng lại với một tác động tức thời. Một tính chất hữu ích khác là khả năng tính toán đáp ứng của hệ thống với bất kỳ tín hiệu đầu vào tùy ý nào thông qua phép tích chập (convolution) giữa hàm trọng lượng g(t) và tín hiệu đầu vào u(t). Trong thực tế, một đáp ứng xung tắt dần về 0 theo thời gian là một dấu hiệu cho thấy hệ thống ổn định.

2.2. Tìm hiểu đáp ứng nấc Step Response của hệ thống

Đáp ứng nấc, ký hiệu là h(t), mô tả hành vi của hệ thống khi đầu vào là một tín hiệu nấc đơn vị (U(s) = 1/s). Tín hiệu ra trong miền Laplace được tính bằng công thức Y(s) = U(s) * G(s) = G(s)/s. Theo tính chất của biến đổi Laplace, đáp ứng nấc là tích phân của đáp ứng xung: h(t) = ∫g(τ)dτ. Vì lý do này, h(t) còn được gọi là hàm quá độ (transient function). Đồ thị của đáp ứng nấc cung cấp những thông tin trực quan và quý giá về chất lượng của hệ thống điều khiển, bao gồm:

• Thời gian xác lập (Settling time): Thời gian cần thiết để đáp ứng đạt và duy trì trong một phạm vi nhỏ quanh giá trị xác lập.
• Độ vọt lố (Overshoot): Mức độ đáp ứng vượt qua giá trị xác lập cuối cùng.
• Thời gian lên (Rise time): Thời gian để đáp ứng đi từ 10% đến 90% giá trị cuối cùng. Các thông số này là tiêu chí quan trọng để đánh giá hiệu suất và độ ổn định của một hệ thống điều khiển trong các ứng dụng thực tế.

2.3. Khái niệm và định nghĩa chi tiết về đáp ứng tần số

Đáp ứng tần số (Frequency Response) là một khái niệm dùng để quan sát đáp ứng của một hệ thống tuyến tính ở trạng thái ổn định khi tín hiệu đầu vào là một tín hiệu hình sin. Một đặc điểm quan trọng của hệ thống tuyến tính là: nếu đầu vào là hình sin, đầu ra ở trạng thái ổn định cũng sẽ là hình sin có cùng tần số, nhưng có thể khác về biên độ và pha. Đáp ứng tần số được định nghĩa là tỷ số phức giữa tín hiệu ra và tín hiệu vào: Frequency Response = Y(jω) / U(jω). Người ta đã chứng minh được rằng giá trị này chính bằng hàm truyền đạt G(s) khi thay thế s bằng jω, tức là G(jω). Vì G(jω) là một hàm phức, nó có thể được biểu diễn dưới dạng G(jω) = M(ω)e^(jφ(ω)), trong đó:

• M(ω) là đáp ứng biên độ, cho biết độ lợi (gain) của hệ thống theo tần số.
• φ(ω) là đáp ứng pha, cho biết độ lệch pha giữa tín hiệu ra và vào theo tần số. Phân tích này rất quan trọng trong việc thiết kế bộ lọc, phân tích độ ổn định và hiệu chỉnh bộ điều khiển.

III. Phương pháp biểu diễn đáp ứng tần số Biểu đồ Bode Nyquist

Để trực quan hóa đáp ứng tần số, các kỹ sư sử dụng các công cụ đồ họa mạnh mẽ, trong đó phổ biến nhất là biểu đồ Bode (Bode diagram) và biểu đồ Nyquist (Nyquist plot). Biểu đồ Bode là một công cụ cực kỳ hiệu quả, bao gồm hai đồ thị riêng biệt trên trục tần số logarit. Đồ thị thứ nhất, Bode magnitude plot, biểu diễn độ lợi L(ω) (tính bằng decibel) theo tần số. Đồ thị thứ hai, Bode phase plot, biểu diễn góc pha φ(ω) theo tần số. Ưu điểm của biểu đồ Bode là khả năng cộng gộp đặc tính của các khâu nối tiếp nhau, giúp đơn giản hóa việc phân tích các hệ thống phức tạp. Trong khi đó, biểu đồ Nyquist biểu diễn đáp ứng tần số G(jω) trên một mặt phẳng phức (hệ tọa độ cực). Trục hoành là phần thực và trục tung là phần ảo khi tần số ω thay đổi từ 0 đến vô cùng. Cả hai loại biểu đồ này đều cung cấp thông tin quan trọng về độ ổn định của hệ thống thông qua các chỉ số như biên độ dự trữ và biên pha dự trữ.

3.1. Kỹ thuật vẽ và đọc hiểu biểu đồ Bode Bode Diagram

Biểu đồ Bode bao gồm hai thành phần: đồ thị biên độ và đồ thị pha. Đồ thị biên độ biểu diễn độ lợi L(ω) = 20log|G(jω)| theo log(ω). Đồ thị pha biểu diễn góc pha φ(ω) = arg(G(jω)) theo log(ω). Một trong những lợi ích lớn nhất của biểu đồ Bode là khả năng vẽ xấp xỉ bằng các đường tiệm cận, giúp đơn giản hóa quá trình phân tích. Quá trình này bắt đầu bằng việc xác định các tần số gãy (corner frequency). Tại mỗi tần số gãy, độ dốc của đường tiệm cận trên đồ thị biên độ sẽ thay đổi. Ví dụ, một khâu trễ bậc nhất sẽ làm giảm độ dốc đi 20dB/decade. Dựa vào hình dạng của biểu đồ Bode, ta có thể suy ra các đặc tính quan trọng của hệ thống. Tần số cắt biên (gain crossover frequency) ωc là tần số mà tại đó độ lợi bằng 0dB. Tần số cắt pha (phase crossover frequency) ω-π là tần số mà tại đó pha bằng -180 độ. Đây là những thông số then chốt để xác định độ ổn định.

3.2. Xác định biên độ dự trữ và biên pha dự trữ qua đồ thị

Biên độ dự trữ (Gain Margin - GM) và biên pha dự trữ (Phase Margin - PM) là hai chỉ số quan trọng nhất để đánh giá độ ổn định tương đối của một hệ thống vòng kín từ đáp ứng tần số của hệ hở.

• Biên độ dự trữ (GM) được định nghĩa là GM = 1 / |G(jω-π)|, trong đó ω-π là tần số cắt pha. Về mặt vật lý, GM cho biết hệ thống có thể tăng thêm bao nhiêu độ lợi trước khi trở nên mất ổn định. Trên biểu đồ Bode, GM [dB] = -L(ω-π). Một giá trị GM dương (dB) cho thấy hệ thống ổn định.
• Biên pha dự trữ (PM) được định nghĩa là ΦM = 180° + φ(ωc), trong đó ωc là tần số cắt biên. PM cho biết hệ thống có thể chịu thêm bao nhiêu độ trễ pha trước khi mất ổn định. Một giá trị PM dương là cần thiết cho sự ổn định. Các giá trị GM và PM càng lớn, hệ thống càng ổn định nhưng đáp ứng có thể chậm hơn. Các kỹ sư thường đặt mục tiêu thiết kế với GM > 6dB và PM trong khoảng 30°-60°.

IV. Bí quyết phân tích động học các khâu cơ bản trong hệ thống

Một hệ thống điều khiển phức tạp thường được cấu thành từ việc kết nối nhiều thành phần cơ bản. Việc hiểu rõ động học của từng khâu cơ bản này là nền tảng để phân tích toàn bộ hệ thống. Các khâu động học cơ bản bao gồm khâu khuếch đại tỷ lệ, khâu tích phân và khâu vi phân, cùng với các khâu bậc cao hơn như khâu trễ bậc nhất và khâu dao động bậc hai. Mỗi khâu có một hàm truyền đạt G(s) đặc trưng, từ đó suy ra được đáp ứng thời gian và đáp ứng tần số riêng. Ví dụ, khâu khuếch đại tỷ lệ G(s) = K chỉ làm thay đổi biên độ tín hiệu mà không gây trễ pha. Khâu tích phân G(s) = 1/s gây ra độ trễ pha -90 độ và có độ lợi giảm 20dB/decade. Ngược lại, khâu vi phân G(s) = s tạo ra độ vượt trước pha +90 độ. Bằng cách phân tích từng thành phần riêng lẻ, chúng ta có thể tổng hợp chúng lại để hiểu hành vi tổng thể của hệ thống điều khiển.

4.1. Động học khâu khuếch đại khâu tích phân và vi phân

Ba khâu động học cơ bản nhất là nền tảng của nhiều bộ điều khiển, đặc biệt là bộ điều khiển PID.

• Khâu khuếch đại tỷ lệ (Proportional Gain): Có hàm truyền đạt G(s) = K. Đáp ứng tần số của nó là một hằng số G(jω) = K, nghĩa là độ lợi L(ω) = 20log(K) không đổi và pha φ(ω) = 0° trên mọi tần số. Khâu này chỉ đơn giản là nhân tín hiệu đầu vào với một hệ số K.
• Khâu tích phân (Integral Factor): Có hàm truyền đạt G(s) = 1/s. Đáp ứng tần số là G(jω) = 1/(jω) = -j/ω. Độ lợi L(ω) = -20log(ω), là một đường thẳng có độ dốc -20dB/decade. Pha φ(ω) = -90° không đổi. Khâu này có tác dụng triệt tiêu sai số xác lập.
• Khâu vi phân (Derivative Factor): Có hàm truyền đạt G(s) = s. Đáp ứng tần số là G(jω) = jω. Độ lợi L(ω) = 20log(ω), là một đường thẳng có độ dốc +20dB/decade. Pha φ(ω) = +90° không đổi. Khâu này giúp cải thiện đáp ứng quá độ và tăng tính ổn định của hệ thống bằng cách dự đoán xu hướng thay đổi của sai số.

4.2. Đặc tính của khâu trễ bậc nhất First order lag factor

Khâu trễ bậc nhất là một trong những thành phần phổ biến nhất trong các hệ thống vật lý, mô tả các quá trình có quán tính như mạch RC, hệ thống nhiệt. Hàm truyền đạt của nó có dạng: G(s) = 1 / (Ts + 1), trong đó T là hằng số thời gian.

• Đáp ứng thời gian: Đáp ứng nấc của khâu này là h(t) = 1 - e^(-t/T), một đường cong tiệm cận đến giá trị 1. Hằng số thời gian T quyết định tốc độ của đáp ứng; T càng nhỏ, hệ thống đáp ứng càng nhanh.
• Đáp ứng tần số: G(jω) = 1 / (Tjω + 1). Biểu đồ Bode của nó có tần số gãy tại ω = 1/T. Với tần số ω << 1/T, độ lợi xấp xỉ 0dB. Với tần số ω >> 1/T, độ lợi giảm với độ dốc -20dB/decade. Góc pha thay đổi từ 0° ở tần số thấp đến -90° ở tần số cao, và bằng -45° tại tần số gãy. Khâu này hoạt động như một bộ lọc thông thấp, cho các tín hiệu tần số thấp đi qua và làm suy giảm các tín hiệu tần số cao.

4.3. Phân tích khâu dao động bậc hai Second order oscillating

Khâu dao động bậc hai mô tả các hệ thống có khả năng dao động, ví dụ như hệ con lắc lò xo có giảm chấn hoặc mạch RLC. Hàm truyền đạt chuẩn của nó là G(s) = 1 / (T²s² + 2ξTs + 1), với 0 < ξ < 1.

• Đáp ứng thời gian: Đáp ứng nấc của khâu này thể hiện dao động tắt dần. Hành vi dao động phụ thuộc mạnh vào hệ số tắt dần ξ (damping ratio). Nếu ξ gần 0, hệ thống có dao động lớn và độ vọt lố cao. Nếu ξ gần 1, đáp ứng mượt hơn và ít dao động hơn.
• Đáp ứng tần số: Biểu đồ Bode của khâu này phức tạp hơn. Tần số gãy là ω = 1/T. Với ω >> 1/T, độ dốc của đồ thị biên độ là -40dB/decade. Pha thay đổi từ 0° đến -180°. Một đặc điểm đáng chú ý là hiện tượng cộng hưởng có thể xảy ra ở gần tần số gãy nếu hệ số tắt dần ξ nhỏ, gây ra một đỉnh nhọn trên đồ thị biên độ. Hiện tượng này chỉ ra rằng hệ thống khuếch đại mạnh các tín hiệu có tần số gần với tần số dao động tự nhiên của nó.

V. Ứng dụng thực tiễn trong phân tích động học hệ thống điều khiển

Lý thuyết về động học hệ thống có ứng dụng sâu rộng trong việc phân tích và thiết kế các hệ thống điều khiển thực tế. Khi một hệ thống điều khiển được cấu thành từ nhiều khâu cơ bản mắc nối tiếp, hàm truyền đạt tổng thể sẽ là tích của các hàm truyền đạt thành phần. Điều này dẫn đến một tính chất quan trọng trên biểu đồ Bode: đồ thị Bode của toàn hệ thống bằng tổng các đồ thị Bode của từng khâu. Kỹ thuật này cho phép các kỹ sư xây dựng xấp xỉ biểu đồ Bode cho các hệ thống phức tạp một cách nhanh chóng. Bằng cách phân tích hàm truyền đạt G(s), ta có thể đưa ra những nhận xét quan trọng về hành vi của hệ thống ở trạng thái ổn định. Ví dụ, sự hiện diện của một khâu tích phân lý tưởng trong G(s) (mẫu số có s) ngụ ý rằng hàm quá độ sẽ tiến đến vô cùng, trong khi sự có mặt của khâu vi phân lý tưởng (tử số có s) cho thấy hàm quá độ sẽ tiến về 0. Những hiểu biết này rất cần thiết để đảm bảo hệ thống hoạt động đúng như thiết kế.

5.1. Cách xấp xỉ biểu đồ Bode cho hệ thống nhiều khâu

Việc vẽ xấp xỉ biểu đồ Bode cho một hệ thống phức tạp là một kỹ năng quan trọng. Giả sử hệ thống có hàm truyền đạt dạng G(s) = K * s^α * G1(s) * G2(s) * .... Quá trình vẽ bao gồm các bước sau:

1. Xác định tần số gãy: Tìm tất cả các tần số gãy ωi = 1/Ti từ các khâu bậc nhất và bậc hai, sau đó sắp xếp chúng theo thứ tự tăng dần.
2. Vẽ điểm bắt đầu: Chọn một tần số ω0 nhỏ hơn tần số gãy đầu tiên. Tại đó, vẽ điểm A có tọa độ (ω0, 20log(K) + α * 20log(ω0)).
3. Vẽ tiệm cận đầu tiên: Từ điểm A, vẽ một đường tiệm cận có độ dốc α * 20 dB/decade (dương nếu là khâu vi phân, âm nếu là khâu tích phân) kéo dài đến tần số gãy tiếp theo.
4. Thay đổi độ dốc: Tại mỗi tần số gãy ωi, cộng thêm vào độ dốc hiện tại một lượng tương ứng với khâu đó (-20dB/dec cho khâu trễ bậc nhất, +20dB/dec cho khâu vượt trước bậc nhất, -40dB/dec cho khâu dao động bậc hai, v.v.).
5. Lặp lại: Tiếp tục quá trình cho đến khi đi qua tất cả các tần số gãy. Phương pháp này cung cấp một cái nhìn tổng quan nhanh chóng và chính xác về đáp ứng tần số của hệ thống.

5.2. Nhận xét về trạng thái ổn định dựa trên hàm truyền G s

Hàm truyền đạt G(s) chứa đựng nhiều thông tin về hành vi của hệ thống ở trạng thái ổn định (t → ∞). Dựa vào định lý giá trị cuối cùng của biến đổi Laplace, ta có thể dự đoán giá trị của đáp ứng mà không cần tìm biểu thức trong miền thời gian.

• Nếu G(s) không chứa khâu tích phân hoặc khâu vi phân lý tưởng (tức là a_n ≠ 0 và b_m ≠ 0), thì đáp ứng xung g(t) sẽ tắt dần về 0, trong khi đáp ứng nấc h(t) sẽ hội tụ về một giá trị hữu hạn khác không (b_m/a_n). Đây là hành vi của một hệ thống ổn định và có độ lợi tĩnh.
• Nếu G(s) chứa một khâu tích phân lý tưởng (a_n = 0), đáp ứng nấc h(t) sẽ tiến đến vô cùng. Hệ thống này không tự ổn định ở trạng thái mới mà sẽ tiếp tục thay đổi.
• Nếu G(s) chứa một khâu vi phân lý tưởng (b_m = 0), đáp ứng nấc h(t) sẽ tiến về 0 ở trạng thái ổn định. Điều này có nghĩa là hệ thống chỉ phản ứng với sự thay đổi của tín hiệu đầu vào, không phải giá trị tuyệt đối của nó. Những nhận xét này giúp đánh giá sơ bộ tính chất của hệ thống ngay từ mô hình toán học.