I. Tổng quan về phân tích ổn định hệ thống trong lý thuyết điều khiển
Phân tích ổn định hệ thống là một nhiệm vụ nền tảng trong kỹ thuật điều khiển tự động. Một hệ thống chỉ hữu ích khi nó hoạt động một cách ổn định. Khái niệm ổn định cơ bản nhất là ổn định BIBO (Bounded-Input, Bounded-Output). Theo định nghĩa này, một hệ thống được coi là ổn định nếu mọi tín hiệu vào bị chặn (bounded input) đều tạo ra một tín hiệu ra bị chặn (bounded output). Nếu một tín hiệu vào bị chặn có thể gây ra một tín hiệu ra không bị chặn, hệ thống được coi là không ổn định. Trạng thái nằm giữa ổn định và không ổn định được gọi là ổn định tới hạn (marginally stable). Việc phân tích này phụ thuộc mật thiết vào hàm truyền đạt (transfer function) của hệ thống, được biểu diễn dưới dạng G(s) = Y(s)/U(s). Mẫu số của hàm truyền đạt, khi cho bằng 0, tạo thành phương trình đặc tính. Nghiệm của phương trình này, được gọi là các cực của hệ thống, quyết định trực tiếp đến tính ổn định. Về mặt toán học, tính ổn định của một hệ thống tuyến tính bất biến theo thời gian (LTI) hoàn toàn được quyết định bởi vị trí các cực và zero trên mặt phẳng phức s. Một hệ thống là ổn định nếu và chỉ nếu tất cả các cực của nó nằm ở nửa bên trái của mặt phẳng phức. Nếu có bất kỳ cực nào nằm ở nửa bên phải, hệ thống sẽ không ổn định. Nếu có cực nằm trên trục ảo và không có cực lặp, hệ thống sẽ ổn định tới hạn. Do đó, mục tiêu chính của việc phân tích ổn định hệ thống là xác định vị trí của các cực này mà không cần phải giải tường minh phương trình đặc tính, vốn có thể rất phức tạp đối với các hệ bậc cao.
1.1. Khái niệm ổn định BIBO và vai trò của hàm truyền đạt
Khái niệm ổn định BIBO (Bounded-Input, Bounded-Output) là tiêu chuẩn cốt lõi để đánh giá tính ổn định. Một hệ thống ổn định BIBO đảm bảo rằng với bất kỳ tín hiệu vào đơn vị (unit step input) hoặc tín hiệu bị chặn nào khác, đầu ra sẽ không tăng đến vô hạn. Hàm truyền đạt G(s) đóng vai trò trung tâm trong việc phân tích này. Đa thức mẫu số của G(s), hay còn gọi là đa thức đặc tính A(s), chứa đựng toàn bộ thông tin về động học nội tại và tính ổn định của hệ thống. Các nghiệm của phương trình đặc tính A(s) = 0 chính là các cực của hệ thống. Vị trí của các cực này trên mặt phẳng phức sẽ quyết định hành vi của hệ thống theo thời gian, bao gồm cả đáp ứng quá độ và đáp ứng trạng thái xác lập.
1.2. Mối quan hệ giữa vị trí các cực và zero và tính ổn định
Vị trí của các cực trên mặt phẳng phức là yếu tố quyết định tính ổn định. Theo tài liệu của Huỳnh Thái Hoàng, một hệ thống được xác định là ổn định nếu tất cả các cực của nó có phần thực âm (nằm ở nửa trái mặt phẳng s). Nếu có ít nhất một cực có phần thực dương (nằm ở nửa phải mặt phẳng s), hệ thống sẽ không ổn định. Các cực nằm trên trục ảo (phần thực bằng 0) dẫn đến trạng thái ổn định tới hạn, hệ thống sẽ dao động với biên độ không đổi. Mặc dù các zero không trực tiếp quyết định tính ổn định, chúng ảnh hưởng đến hình dạng của đáp ứng thời gian và có thể làm triệt tiêu hoặc khuếch đại các mode dao động tương ứng với các cực.
II. Các thách thức khi xác định ổn định từ phương trình đặc tính
Thách thức lớn nhất trong phân tích ổn định hệ thống là việc xác định vị trí các cực mà không cần giải phương trình đặc tính bậc cao. Đối với các hệ bậc nhất hoặc hệ bậc hai, việc tìm nghiệm là tương đối đơn giản. Tuy nhiên, với các hệ thống thực tế có bậc từ 3 trở lên, việc giải tường minh phương trình đại số để tìm chính xác giá trị các cực là không khả thi và tốn kém về mặt tính toán. Một phương pháp tiếp cận ban đầu là kiểm tra điều kiện cần. Điều kiện cần để một hệ thống ổn định là tất cả các hệ số của phương trình đặc tính phải tồn tại và cùng dấu (thường là dương). Nếu bất kỳ hệ số nào bằng 0 hoặc trái dấu, hệ thống chắc chắn có ít nhất một cực ở nửa phải mặt phẳng phức hoặc trên trục ảo, và do đó không ổn định. Tuy nhiên, đây chỉ là điều kiện cần chứ không phải điều kiện đủ. Một hệ thống có thể có tất cả các hệ số dương nhưng vẫn không ổn định. Ví dụ, phương trình s³ + 4s² + 5s + 2 = 0 có tất cả hệ số dương, nhưng việc kiểm tra sâu hơn là cần thiết để kết luận. Do đó, các kỹ sư cần những công cụ mạnh mẽ hơn, các tiêu chuẩn đại số cho phép kết luận về tính ổn định mà không cần tìm nghiệm trực tiếp. Các tiêu chuẩn này cung cấp một câu trả lời "có/không" về tính ổn định và thậm chí còn cho biết số lượng cực không ổn định.
2.1. Điều kiện cần Hạn chế của việc kiểm tra hệ số đa thức
Việc kiểm tra các hệ số của đa thức đặc tính là bước đầu tiên và nhanh nhất. Một đa thức đặc tính A(s) = a₀sⁿ + a₁sⁿ⁻¹ + ... + aₙ = 0 phải có tất cả các hệ số aᵢ > 0 để hệ thống có thể ổn định. Nếu điều kiện này không được thỏa mãn, có thể kết luận ngay hệ thống không ổn định. Ví dụ, phương trình s³ + 3s² - 2s + 1 = 0 chắc chắn không ổn định vì có hệ số âm. Tuy nhiên, nếu tất cả các hệ số đều dương, không thể kết luận ngay là hệ thống ổn định. Đây là một cạm bẫy phổ biến và đòi hỏi phải áp dụng các tiêu chuẩn đầy đủ hơn như Routh-Hurwitz.
2.2. Sự cần thiết của các tiêu chuẩn đại số để phân tích hệ bậc cao
Đối với các hệ thống có bậc n ≥ 3, việc giải phương trình đặc tính là bất khả thi. Các tiêu chuẩn đại số như Routh-Hurwitz và Hurwitz cung cấp một phương pháp hệ thống để xác định sự tồn tại của các cực có phần thực dương. Các phương pháp này chỉ dựa trên các hệ số của phương trình đặc tính và không yêu cầu tính toán nghiệm. Chúng không chỉ trả lời câu hỏi hệ thống có ổn định hay không, mà còn cho biết chính xác có bao nhiêu cực nằm ở nửa phải mặt phẳng phức, cung cấp thông tin quý giá cho việc thiết kế và hiệu chỉnh bộ điều khiển.
III. Hướng dẫn áp dụng tiêu chuẩn ổn định Routh Hurwitz chi tiết nhất
Một trong những công cụ đại số mạnh mẽ nhất để phân tích ổn định hệ thống là tiêu chuẩn ổn định Routh-Hurwitz. Tiêu chuẩn này cung cấp một điều kiện cần và đủ để xác định tính ổn định của một hệ thống tuyến tính mà không cần phải tìm các cực của phương trình đặc tính. Phương pháp này dựa trên việc xây dựng một bảng gọi là bảng Routh từ các hệ số của đa thức đặc tính a₀sⁿ + a₁sⁿ⁻¹ + ... + aₙ = 0. Quy trình bắt đầu bằng việc điền hai hàng đầu tiên của bảng với các hệ số chẵn và lẻ của đa thức. Các hàng tiếp theo được tính toán một cách đệ quy dựa trên hai hàng ngay phía trên nó. Sau khi hoàn thành bảng Routh, tính ổn định của hệ thống được xác định bằng cách kiểm tra các phần tử trong cột đầu tiên. Điều kiện cần và đủ để hệ thống ổn định là tất cả các phần tử trong cột đầu tiên của bảng Routh phải có cùng dấu (dương). Nếu có bất kỳ sự đổi dấu nào trong cột đầu tiên, hệ thống sẽ không ổn định. Hơn nữa, số lần đổi dấu trong cột đầu tiên chính bằng số lượng cực nằm ở nửa bên phải mặt phẳng phức. Phương pháp này cực kỳ hữu ích trong việc phân tích ảnh hưởng của một tham số (ví dụ như độ lợi K) đến tính ổn định của hệ thống.
3.1. Quy tắc xây dựng bảng Routh từ phương trình đặc tính
Để xây dựng bảng Routh, hàng thứ nhất bao gồm các hệ số có chỉ số chẵn (a₀, a₂, a₄, ...). Hàng thứ hai bao gồm các hệ số có chỉ số lẻ (a₁, a₃, a₅, ...). Các phần tử của các hàng tiếp theo (từ hàng thứ ba trở đi) được tính theo một công thức xác định. Ví dụ, phần tử đầu tiên của hàng thứ ba được tính bằng: b₁ = (a₁a₂ - a₀a₃)/a₁. Quá trình này được lặp lại cho đến khi tất cả các hàng được điền đủ, với tổng số n+1 hàng cho một đa thức bậc n. Sự chính xác trong việc xây dựng bảng là rất quan trọng để có kết luận đúng đắn.
3.2. Cách diễn giải kết quả và xử lý các trường hợp đặc biệt
Diễn giải bảng Routh rất trực tiếp: không có sự đổi dấu nào ở cột đầu tiên có nghĩa là hệ thống ổn định. Số lần đổi dấu cho biết số cực không ổn định. Tuy nhiên, có hai trường hợp đặc biệt cần xử lý. Trường hợp thứ nhất là khi một phần tử ở cột đầu tiên bằng 0 nhưng các phần tử khác trên cùng hàng không bằng 0. Trong tình huống này, số 0 được thay thế bằng một số dương rất nhỏ (epsilon ε) và tiếp tục tính toán. Dấu của các phần tử bên dưới sẽ phụ thuộc vào dấu của ε. Trường hợp thứ hai là khi cả một hàng chứa toàn số 0. Điều này chỉ ra sự tồn tại của các cặp cực đối xứng qua gốc tọa độ. Một đa thức phụ được hình thành từ hàng ngay trên hàng số 0, và hàng số 0 được thay thế bằng các hệ số của đạo hàm đa thức phụ đó.
IV. Phương pháp phân tích đáp ứng thời gian cho hệ bậc nhất và bậc hai
Việc phân tích đáp ứng thời gian là một cách tiếp cận quan trọng khác để hiểu về hành vi của hệ thống, đặc biệt là đối với các hệ bậc nhất và hệ bậc hai, vốn là nền tảng cho việc mô hình hóa nhiều hệ thống phức tạp hơn. Đáp ứng thời gian của một hệ thống bao gồm hai thành phần chính: đáp ứng quá độ (transient response) và đáp ứng trạng thái xác lập (steady-state response). Đáp ứng quá độ là hành vi của hệ thống khi nó chuyển từ trạng thái ban đầu đến trạng thái cuối cùng. Đáp ứng trạng thái xác lập là hành vi của đầu ra khi thời gian tiến đến vô cùng. Đối với một hệ bậc nhất với hàm truyền G(s) = K/(τs+1), đáp ứng với tín hiệu vào bậc thang đơn vị được đặc trưng bởi hằng số thời gian τ. Thời gian xác lập (settling time) thường được coi là khoảng 4τ. Đối với một hệ bậc hai chuẩn G(s) = ωₙ² / (s² + 2ζωₙs + ωₙ²), hành vi của nó được quyết định bởi hai tham số quan trọng: tần số dao động tự nhiên (natural frequency) ωₙ và tỷ số cản (damping ratio) ζ. Các chỉ tiêu chất lượng này mô tả hệ thống phản ứng nhanh như thế nào và mức độ dao động của nó.
4.1. Đặc điểm đáp ứng quá độ của hệ thống điều khiển
Đáp ứng quá độ của một hệ thống được đánh giá thông qua các chỉ tiêu chất lượng quan trọng. Thời gian lên (rise time) là thời gian để đáp ứng đi từ 10% đến 90% giá trị cuối cùng. Độ vọt lố (percent overshoot) là giá trị đỉnh lớn nhất của đường cong đáp ứng được đo từ giá trị xác lập, thể hiện mức độ dao động của hệ thống. Thời gian xác lập (settling time) là thời gian cần thiết để đáp ứng đạt đến và duy trì trong một phạm vi nhỏ (thường là 2% hoặc 5%) quanh giá trị xác lập. Các chỉ tiêu này liên quan trực tiếp đến vị trí các cực của hệ bậc hai, cụ thể là tỷ số cản và tần số dao động tự nhiên.
4.2. Phân tích sai số xác lập và đáp ứng trạng thái ổn định
Đáp ứng trạng thái xác lập mô tả hành vi dài hạn của hệ thống. Một chỉ tiêu quan trọng trong giai đoạn này là sai số xác lập (steady-state error), là sự khác biệt giữa tín hiệu vào (mong muốn) và tín hiệu ra (thực tế) khi thời gian tiến đến vô cùng. Sai số này phụ thuộc vào loại của hệ thống (số lượng cực tại gốc tọa độ) và dạng của tín hiệu vào (bậc thang, dốc, parabol). Việc phân tích sai số xác lập rất quan trọng trong các ứng dụng điều khiển vị trí hoặc theo dõi quỹ đạo, nơi độ chính xác cao là yêu cầu bắt buộc.
V. Ứng dụng Quỹ đạo nghiệm số để đánh giá ổn định và thiết kế bộ điều khiển
Phương pháp Quỹ đạo nghiệm số (Root Locus) là một công cụ đồ họa mạnh mẽ cho phép phân tích ổn định hệ thống và thiết kế bộ điều khiển. Quỹ đạo nghiệm số là tập hợp tất cả các vị trí có thể có của các cực vòng kín khi một tham số hệ thống, thường là độ lợi K, thay đổi từ 0 đến vô cùng. Bằng cách vẽ đồ thị này, kỹ sư có thể hình dung được các cực di chuyển như thế nào trên mặt phẳng phức và từ đó xác định phạm vi giá trị của K để hệ thống duy trì ổn định. Một hệ thống sẽ trở nên không ổn định khi bất kỳ nhánh nào của quỹ đạo nghiệm số đi vào nửa bên phải của mặt phẳng phức. Điểm mà quỹ đạo cắt trục ảo tương ứng với giá trị độ lợi tới hạn K_cr, tại đó hệ thống bắt đầu dao động. Phương pháp này không chỉ cho biết hệ thống ổn định hay không, mà còn cung cấp cái nhìn sâu sắc về hành vi động học của hệ thống. Ví dụ, nó cho thấy hệ thống sẽ trở nên dao động hơn hay ít hơn khi K thay đổi, thông qua việc quan sát sự di chuyển của các cực đối với trục thực. Việc xây dựng quỹ đạo nghiệm số tuân theo một tập hợp các quy tắc dựa trên điều kiện pha và điều kiện biên độ của phương trình đặc tính, giúp vẽ phác thảo quỹ đạo một cách nhanh chóng mà không cần tính toán phức tạp.
5.1. Nguyên tắc và các quy tắc vẽ quỹ đạo nghiệm số
Việc vẽ quỹ đạo nghiệm số dựa trên phương trình đặc tính 1 + K*G(s)H(s) = 0. Một điểm s trên mặt phẳng phức thuộc quỹ đạo nghiệm số nếu nó thỏa mãn điều kiện pha: góc của G(s)H(s) phải bằng (2l+1)180°. Các quy tắc chính bao gồm: quỹ đạo bắt đầu từ các cực và kết thúc tại các zero của vòng hở; số nhánh bằng số cực; các đoạn tiệm cận cho các nhánh đi ra vô cùng; các điểm tách/nhập trên trục thực; và các góc xuất phát/góc tới tại các cực/zero phức. Những quy tắc này cho phép phác họa nhanh chóng hành vi của hệ thống.
5.2. Diễn giải quỹ đạo nghiệm số để xác định khoảng ổn định
Sau khi vẽ, việc diễn giải quỹ đạo nghiệm số là chìa khóa. Phần quỹ đạo nằm ở nửa trái mặt phẳng phức tương ứng với các giá trị của K làm cho hệ thống ổn định. Điểm giao của quỹ đạo với trục ảo (nếu có) là cực kỳ quan trọng. Tại điểm này, hệ thống ở trạng thái ổn định tới hạn. Giá trị của K tại giao điểm này được gọi là độ lợi tới hạn. Bất kỳ giá trị K nào lớn hơn độ lợi tới hạn sẽ làm cho hệ thống không ổn định. Bằng cách này, quỹ đạo nghiệm số cung cấp một công cụ thiết kế trực quan để chọn độ lợi K nhằm đạt được cả tính ổn định và các chỉ tiêu chất lượng mong muốn như độ vọt lố và thời gian xác lập.
