Tiêu chuẩn ổn định Routh–Hurwitz

Trong lý thuyết hệ thống điều khiển, tiêu chuẩn ổn định Routh-Hurwitz là một kiểm tra toán học là một điều kiện cần và đủ cho sự ổn định của một hệ thống điều khiển tuyến tính thời gian bất biến (LTI). Kiểm tra Routh là một thuật toán đệ quy hiệu quả mà nhà toán học người Anh Edward John Routh đề xuất vào năm 1876 để xác định xem tất cả các nghiệm của đa thức đặc trưng của một hệ thống tuyến tính có các phần thực âm hay không. Nhà toán học người Đức Adolf Hurwitz cũng độc lập đề xuất vào năm 1895 để sắp xếp các hệ số của đa thức thành một ma trận vuông, gọi là ma trận Hurwitz, và cho thấy rằng đa thức này là ổn định nếu và chỉ nếu dãy tất cả các định thức của các ma trận con chính của nó đều dương. Hai cách thức tiếp cận trên là tương đương nhau, với kiểm tra Routh cung cấp một cách hiệu quả hơn để tính toán các định thức Hurwitz so với tính toán chúng trực tiếp. Một đa thức thỏa mãn tiêu chuẩn Routh-Hurwitz được gọi là một đa thức Hurwitz.

Điều quan trọng của tiêu chuẩn này là các nghiệm p của phương trình đặc trưng của một hệ thống tuyến tính với các phần thực âm đại diện cho các nghiệm e của hệ thống đó là ổn định (bị chặn). Vì thế tiêu chuẩn này cung cấp một cách để xác định xem các phương trình chuyển động của một hệ thống tuyến tính có lời giải nào ổn định không mà thôi, mà không cần phải giải hệ thống này một cách trực tiếp. Đối với các hệ thống rời rạc, các kiểm tra sự ổn định tương ứng có thể được thực hiện bởi tiêu chuẩn Schur-Cohn, kiểm tra Jury và các kiểm tra Bistritz. Với sự ra đời của máy vi tính, tiêu chuẩn này đã trở nên ít sử dụng rộng rãi, như một cách để thay thế việc giải đa thức bằng số học, lấy xấp xỉ nghiệm trực tiếp.

Kiểm tra Routh có thể được bắt nguồn bằng cách sử dụng giải thuật Euclid và định lý Sturm trong việc đánh giá chỉ số Cauchy. Hurwitz đưa ra các điều kiện của mình theo một cách khác.

== Sử dụng giải thuật Euclid ==
Tiêu chuẩn này liên quan đến định lý Routh-Hurwitz. Thật vậy, từ phát biểu của định lý này, chúng ta có formula_1 trong đó:
- "p" là số nghiệm của đa thức "ƒ"("z") với Phần Thực âm;
- "q" là số nghiệm của đa thức "ƒ"("z") với Phần Thực dương (nên nhớ rằng "ƒ" được cho là không có nghiệm nào nằm trên trục ảo);
- "w"("x") là số biến của chuỗi Sturm suy rộng đạt được từ formula_2 và formula_3 (bởi Giải thuật Euclidean nối tiếp) trong đó formula_4 đại diện cho một số thực "y".
Bằng định lý cơ bản của đại số, mỗi đa thức bậc "n" phải có "n" nghiệm trong mặt phẳng phức (ví dụ, cho một "ƒ" không có nghiệm nào nằm trên trục ảo, "p" + "q" = "n"). Do đó, ta có điều kiện là "ƒ" là một đa thức ổn định (Hurwitz) nếu và chỉ nếu" p" − "q" = "n" (chứng minh được đưa ra dưới đây). Sử dụng định lý Routh–Hurwitz, chúng ta có thể thay thế điều kiện trên "p" và "q" bởi một điều kiện trên chuỗi Sturm suy rộng, sẽ cho lần lượt cho một điều kiện trên cá hệ số của "ƒ".

== Sử dụng ma trận ==
Cho "f"("z") là một đa thức phức. Chu trình này theo như sau:
1. Tính toán các đa thức formula_2 và formula_3 nhưformula_4 trong đó" y" là một số thực.
2. Tính toán ma trận liên quan đến formula_2 và formula_3.
3. Sắp xếp lại mỗi hàng theo cách tương tự trong đó một hàng lẽ và hàng theo sau có cùng zero ở đầu.
4. Tính toán mỗi định thức con chính của ma trận đó.
5. Nếu ít nhất một trong số các định thức con đó là âm (hoặc zero), thì đa thức "f "là không ổn định.

=== Ví dụ ===
- Cho formula_10 (để đơn giản chúng ta chỉ lấy các hệ số phần thực) trong đó formula_11 (để tránh một nghiệm trong zero do đó ta có thể sử dụng định lý Routh–Hurwitz). Đầu tiên, ta phải tính các đa thức thực formula_2 vàformula_3:
Chú ý rằng chúng ta phải giả sử "b" khác zero trong phép chia đầu tiên. Chuỗi Sturm tổng quát trong trường hợp này là 

=== Tiêu chuẩn Routh–Hurwitz cho các đa thức bậc hai, ba và bốn ===
Dưới đây, chúng ta giả sử rằng hệ số của bậc cao nhất (ví dụ. formula_26 trong một đa thức bậc hai) là dương. Nếu cần thiết,luôn luôn có thể đạt được điều này bằng cách nhân đa thức này với formula_27.
- Đối với một đa thức bậc hai, formula_28, tất cả các nghiệm là nằm trong mặt phẳng bên trái (và hệ thống này với phương trình đặc tính formula_29 là ổn định) nếu tất cả các hệ số thỏa mãn formula_30.
- Đối với một đa thức bậc baformula_31, tất cả các hệ số phải thỏa mãn  formula_30, và formula_33
- Đối với một đa thức bậc bốn formula_34, tất cả các hệ số phải thỏa mãn formula_30, và formula_36 và formula_37 Nói chung tiêu chuẩn ổn định Routh phát biểu rằng tất cả các thành phần Cột đầu tiên của mảng Routh là cùng dấu.
Các hệ thống đáp ứng tiêu chí nêu trên được gọi là vòng lặp kín ổn định, nếu không chúng là không ổn định vì đổi dấu trong các thành phần cột đầu tiên.

=== Ví dụ về bậc cao ===
Có thể sử dụng phương pháp bảng để xác định sự ổn định khi các nghiệm của một đa thức đặc tính bậc cao khó đạt được. Đối với đa thức "bậc thứ n"
- formula_38
bảng này có "n" + 1 hàng và có dạng sau:

trong đó các thành phần formula_39 và formula_40 có thể được tính toán như sau:
- formula_41
- formula_42
Khi hoàn thành, số dấu thay đổi trong cột đầu tiên sẽ là số các cực không-âm.

Trong cột đầu tiên, có hai lần đổi dấu (0.75 → −3, và −3 → 3), do đó có hai nghiệm không-âm trong đó hệ thống sẽ không ổn định.

Đôi khi sự hiện diện của các cực trên trục ảo tạo ra một trạng thái ổn định biên. Trong trường hợp này các hệ số của "Mảng Routh" trong toàn bộ một hàng trở thành zero và do đó đáp án xa hơn cho đa thức  này là tìm các lần đổi dấu là không thể. Thì cách tiếp cận khác sẽ được xem xét.  Hàng nằm trên hàng chứa các zero của đa thức này được gọi là "Đa thức Phụ".
- formula_43
Chúng ta có bảng sau:

Trong trường hợp đa thức Phụ này là formula_44 một lần nữa lại bằng zero. Bước tiếp theo là đạo hàm phương trình ở trên để đạt được đa thức tiếp theo. formula_45. Các hệ số của hàng chứa zero bây giờ trở thành "8" và "24". Chu trình mảng Routh được tiến hành sử dụng các giá trị này để đạt được hai điểm trên trục ảo. Hai điểm nằm trên trục ảo này là nguyên nhân chính gây ra độ ổn định biên.

== Xem thêm ==
- Kỹ thuật điều khiển
- Phái sinh của mảng Routh
- Tiêu chuẩn ổn định Nyquist
- Định lý Routh-Hurwitz
- Quỹ đạo nghiệm số
- Hàm truyền
- Tiêu chuẩn ổn định Jury
- Tiêu chuẩn ổn định Bistritz
- Định lý Kharitonov
- Tiêu chuẩn Liénard–Chipart

== Tham khảo ==
- Felix Gantmacher (J.L. Brenner translator) (1959) "Applications of the Theory of Matrices", pp 177–80, New York: Interscience.

== Liên kết ngoài ==
- Một đoạn mã MATLAB thực thi kiểm tra Routh-Hurwitz