Table of Contents
Introduzione al Criterio Routh-Hurwitz
Il criterio Routh-Hurwitz è un punto di riferimento della teoria del controllo classico, fornendo un metodo sistematico per valutare la stabilità dei sistemi lineari di time-invariant (LTI) senza calcolare esplicitamente le radici dell'equazione caratteristica. Sviluppato indipendentemente da Edward John Routh nel 1874 e Adolf Hurwitz nel 1895, questa tecnica algebrica rimane ampiamente utilizzata nell'ingegneria di controllo, perché offre un modo rapido e affidabile per determinare se un sistema risponderà in una stabilità
Analizzando i coefficienti del polinomio caratteristico di un sistema, il criterio Routh-Hurwitz identifica il numero di poli (radici dell'equazione caratteristica) che si trovano nella metà destra del piano complesso. I polacchi con parti reali positive causano un comportamento instabile, mentre i poli con parti reali negative contribuiscono a risposte stabili e convergenti dettagliate. Il metodo non richiede la risoluzione di polinomi di alto grado, rendendolo pratico per l'analisi manuale e i dettaglianti.
Questo articolo spiegherà prima la base matematica del criterio Routh-Hurwitz, poi ti guiderà attraverso la costruzione dell'array Routh con un esempio passo per passo. Dopo di che, esploreremo come il criterio predice risposta del sistema e comportamento transitorio, comprese le caratteristiche chiave come overshoot e impostare il tempo. Un esempio numerico pratico illustra il processo, e discuteremo i limiti del metodo e gli strumenti complementari.
Qual è il criterio Routh-Hurwitz?
Il criterio Routh-Hurwitz è una condizione necessaria e sufficiente per la stabilità nei sistemi LTI. Funziona sull'equazione caratteristica del sistema, che deriva dal denominatore della funzione di trasferimento a ciclo chiuso dopo averlo impostato a zero. Per un sistema descritto dalla funzione di trasferimento , l'equazione caratteristica è , dove è un polinomio in s di grado
Il criterio afferma che per il sistema essere stabile (cioè tutti i poli hanno parti reali negative), tutti i coefficienti ai] devono essere nonzero e hanno lo stesso segno (tipalmente positivo), e il array Routh costruito da questi coefficienti non deve avere modifiche di segno nella prima colonna.
È importante notare che il criterio Routh-Hurwitz è puramente algebrico, non dà le posizioni dei poli reali. Invece, fornisce un conteggio dei poli destro-half (RHP) poli. Ciò lo rende inestimabile quando il polinomio caratteristico è di alto ordine (ad esempio, quarto, quinto o superiore) dove il ritrovamento delle radici è computazionalmente intensivo o analiticamente infesabile.
Il criterio è direttamente collegato alla matrice Hurwitz, una matrice quadrata costituita dai coefficienti del polinomio caratteristico. Tutti i principali minori della matrice Hurwitz devono essere positivi per la stabilità, questa è la condizione di Hurwitz. Il metodo tabulare di Routh semplifica questo processo di controllo determinante nella più semplice array Routh.
Costruendo il Routh Array: Step-by-Step
Costruire la linea Routh è una procedura sistematica, che la illustreremo usando un esempio concreto polinomiale:
Questo polinomio ha radici note a s = -1, -2, -3, -4 (tutti negativi), quindi il sistema è stabile. Lo useremo per verificare che il array Routh non dia nessun cambiamento di segno.
Passo 1: Scrivere il polinomio in forma standard
[LT][[[[]]][[[LT]]]][[[[[[[[]]]]]]]][[[[LT]]]]]]][[[[[[[LT]]]]]]][[[[[[[[[[[[f]]]]]]]]]]]]]]][[[LT]]]]]]]]]][[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[FLT]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]
Fase 2: inizializzare le prime due righe
La prima riga della serie Routh contiene i coefficienti di termini anche indicizzati (a partire da sn]): riga s4: [1, 35, 24]. La seconda riga contiene i coefficienti di termini indici dispari: file s3
]4: 1, 35, 24
3: 10, 50, 0
Passo 3: Compute righe successive utilizzando la regola Rhombus
Ogni elemento nella riga successiva viene calcolato utilizzando le due righe sopra di esso. Per un elemento in riga s[2], prima colonna, la formula è:
Ma è più facile applicare la formula standard: per riga k, colonna j (a partire dalla colonna 1), l'elemento è:
Facciamolo manualmente per chiarezza.
[LT] [[7]] [[7]] [[6]] [[6]]] [[6]] [[6]]] [[6]]][[6]]][[6]] [[6]] [[6]]]][[[6]]]][[[6]]]]][[[7]]]]][[[6]]]]]]]]
Così la riga s2: 30, 24, 0.
[FLT] [[FLT]]]
]] Usare le righe s[3 e s2]
R(s]
Così la riga s1: 42, 0.
[6][7] [[7]] [[7]]] [[6]]] [[6]]] [[[7]]]][[[7]]]][[[[7]]]][[[[6]]]]]] [[[7]]]]] [[[[7]]]]]]]] [[[FLT]]]]]]]] [[[[[[[[[[[[[[[[[[[7]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]][[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[
Quindi l'array completo di Routh è:
| s4 | 1 | 35 | 24 |
|---|---|---|---|
| s3 | 10 | 50 | 0 |
| s2 | 30 | 24 | 0 |
| s1 | 42 | 0 | 0 |
| s0 | 24 | 0 | 0 |
Passo 4: Esaminare la prima colonna
La prima colonna contiene: 1, 10, 30, 42, 24. Tutti sono positivi, e non ci sono cambiamenti di segno. Pertanto il sistema è stabile, che corrisponde alle nostre radici conosciute.
Casi speciali
Due casi speciali[] possono sorgere durante la costruzione di array Routh:
- Zero nella prima colonna:[] Se il primo elemento di una riga è zero, ma la riga non è tutti gli zeri, sostituire lo zero con un piccolo epsilon (ε), quindi procedere. Il segno di ε determina i cambiamenti dei segni. Ciò indica stabilità marginale o instabilità a seconda del resto dell'array.
- Flotta di zeri: Questo accade quando ci sono radici simmetriche (ad esempio, coppie di radici immaginarie o coppie simmetriche sull'asse reale). In questo caso, formare un polinomio ausiliario dalla fila sopra la riga zero e differenziarlo per ottenere coefficienti per la fila mancante. Ciò indica la presenza di un paio di radici puramente immaginarie o di una coppia di radici reali.
Questi casi speciali estendono l'applicabilità del criterio ai sistemi con pali sull'asse immaginario o modelli simmetrici.
Predivisione della risposta del sistema e del comportamento transitorio
Il criterio Routh-Hurwitz non è limitato ad un verdetto binario stabile/non stabile. Rivelando il numero di poli con parti reali positive, informa direttamente il comportamento transitorio del sistema. Ogni polo nella metà destra del piano contribuisce a una modalità di crescita esponenziale, causando la risposta transitoria alla diverge.
Tuttavia, il criterio non dà le posizioni esatte del polo, quindi non può esattamente quantificare il tempo di risoluzione, di impostazione o di aumento.
- Se la prima colonna non ha modifiche di segno:[ Tutti i poli hanno parti reali negative. La risposta transitoria sarà stabile, alla fine impostata a uno stato costante. Ma potrebbe esserci un decadimento lento o oscillatorio (quasi i poli immaginari) che il criterio non rileva, solo il segno della parte reale è garantito.
- Se ci sono cambiamenti di segno:[] Il numero di segni cambia è uguale al numero di poli RHP. Il sistema è instabile, e la risposta transitoria crescerà in modo incontaminato. Il numero di pali RHP indica la natura dell'instabilità: un singolo polo positivo porta alla crescita esponenziale monotonica, mentre un paio di pali RHP complessi produce oscillazioni.
- Stabilità marginale:[] Se l'equazione caratteristica ha radici sull'asse immaginario (poco pali immaginari o multipli sull'asse), l'array di Routh produrrà una fila di zeri. Ciò indica che il sistema è sul limite della stabilità, e la risposta transitoria può mostrare oscillazioni sostenute (se semplici pali immaginari) o oscillazioni rampa-up (se ripetute poli immaginari).
Per un sistema di secondo ordine, il criterio riduce a richiedere che tutti i coefficienti siano positivi (che è necessario e sufficiente per la stabilità). Per i sistemi di ordine superiore, l'array fornisce un indicatore di margine di stabilità — più vicino gli elementi di primo livello sono a zero, più il sistema è di instabilità. Questo è particolarmente utile quando si regola un loop.
Rapporto con caratteristiche transitorie
La risposta transitoria di un sistema stabile è determinata principalmente dai poli dominanti, quelli più vicini all'asse immaginario. Il criterio Routh-Hurwitz non individua direttamente questi poli, ma può aiutare a sintonizzare. Ad esempio, se il criterio mostra che l'aumento di un guadagno del controller alla fine porterà a un cambiamento di segno nella prima colonna, segna il margine di guadagno.
Le metriche transitorie chiave che sono influenzate dalla stabilità includono:
- Overshoot[] – Il valore massimo della risposta rispetto al valore finale. La risoluzione è generalmente causata da complessi poli coniugati con un basso rapporto di smorzamento. Mentre il criterio Routh-Hurwitz non calcola l'umidità, un sistema con una fila di zero (stabilità marginale) corrisponde spesso ad un rapporto di smorzamento vicino a zero, che implica un alto livello di risoluzione.
- Tempo di regolazione[ – Il tempo necessario per la risposta a rimanere entro una certa percentuale (ad esempio, 2%) del valore finale. Il tempo di impostazione è inversamente proporzionale alla parte reale del polo dominante. L'array Routh non può dare questa parte reale, ma se l'array mostra che tutti i poli sono nella metà piano sinistro, il tempo di chiusura dell'asse finito è molto lungo.
- Tempo di recupero[] – Il tempo di passare dal 10% al 90% del valore finale. Questo è legato alla frequenza naturale dei poli complessi dominanti. Ancora una volta, il criterio Routh-Hurwitz non fornisce alcun valore numerico diretto ma indica se il sistema è abbastanza stabile per avere un tempo di salita ben definito.
- Oscillazioni[ – La frequenza delle oscillazioni transitorie è determinata dalle parti immaginarie dei poli complessi. Il criterio non rivela parti immaginarie, ma la presenza di una fila di zero indica pali puramente immaginari, che producono oscillazioni sostenute nella risposta transitoria.
Così, il criterio Routh-Hurwitz agisce come un gatekeeper: assicura la stabilità come prerequisito per un'analisi significativa dei transienti utilizzando altri strumenti come locus root, diagrammi Bode o simulazioni di risposta passo.
Esempio pratico: Sistema di Terze Ordini
Applichiamo il criterio Routh-Hurwitz ad un esempio più complesso con parametri sconosciuti. Considerate la funzione di trasferimento a ciclo chiuso [] con feedback di unità. L'equazione caratteristica è:
Vogliamo determinare la gamma di K per la quale il sistema è stabile. Costruisci l'array Routh per il polinomio .
[LT] [[FLT]] [[FLT]]] [[FLT]]] [[[FLT]]]]] [[[FLT]]]]] [[[[FLT]]]] [[[FLT]]]]] [[[FLT]]]]] [[[[[FLT]]]]]]]] [[FLT]]]]]] [[[
| s3 | 1 | 10 |
|---|---|---|
| s2 | 7 | K |
| s1 | (K-70)/7 | 0 |
| s0 | K | 0 |
Per la stabilità, tutti gli elementi di prima categoria devono essere positivi: 1 > 0, 7 > 0, (K-70)/7 > 0, K > 0. Questo dà K > 70 e K > 0, quindi K > 70. Il limite inferiore è K > 70. Se K = 70, la riga s1 diventa zero, che è un caso specifico indica
Implicazioni di risposta transitorie per l'esempio
Per K appena sopra i 70, il sistema è stabile ma i poli dominanti sono molto vicini all'asse immaginario. La risposta transitoria mostrerà molto lento decadimento e grande overshoot. Il tempo di impostazione sarà lungo. Come K aumenta ulteriormente, i poli si muovono verso sinistra (generalmente), migliorando i tempi di smorzamento e di settling. Tuttavia, se K è troppo alto, l'array Routh per questo particolare sistema non mostra ulteriori cambiamenti di stabilità - la stabilità persiste.
Il criterio fornisce quindi informazioni critiche: la gamma di guadagno per la stabilità e la vicinanza all'instabilità. Valutando gli elementi di prima colonna, gli ingegneri possono calcolare i margini di guadagno e i controllori di progettazione che evitano la stabilità marginale.
Limitazioni del criterio Routh-Hurwitz
Nonostante il suo potere, il criterio Routh-Hurwitz ha limitazioni importanti che ogni ingegnere di controllo deve comprendere:
- Applicabile solo ai sistemi LTI:[ Non può essere utilizzato direttamente per sistemi di tempo, non lineari o dati campionati (anche se esistono estensioni per sistemi discreti-time tramite trasformazioni bilineari).
- Non quantifica le prestazioni transitorie:[] Il criterio non fornisce informazioni su overshoot, picco tempo, tempo di salita o tempo di impostazione esatto.
- Richiede un polinomio caratteristico noto:[ I coefficienti devono essere conosciuti esattamente. In pratica, le incertezze nei parametri possono richiedere un'analisi robusta della stabilità (ad esempio, utilizzando il teorema di Kharitonov).
- Diventa ingombrante per sistemi di ordine molto elevato:[ I calcoli manuali per i polinomi di decimo ordine sono di errore-prone, anche se il software (MATLAB, Python) automatizza la costruzione dell'array.
- Non si può gestire i ritardi di tempo direttamente:[ I sistemi con ritardi di tempo hanno equazioni caratteristiche trascendentali – il array Routh non può essere applicato direttamente a meno che non si utilizzino le approssimazioni Padé.
Per superare queste limitazioni, gli ingegneri usano comunemente il criterio Routh-Hurwitz in combinazione con altri metodi. Ad esempio, le trame di root mostrano come i poli a ciclo chiuso si muovono con guadagno, fornendo informazioni esatte. Il criterio Nyquist gestisce ritardi di tempo e non linearità. Tuttavia, il array Routh spesso funge da primo controllo rapido, soprattutto durante la progettazione preliminare o la risoluzione dei problemi.
Strumenti complementari
Il criterio Routh-Hurwitz è “necessario e sufficiente” per la stabilità dei polinomi LTI, ma per l’analisi delle prestazioni, gli ingegneri si affidano a:
- Root Locus:[] Trama la traiettoria dei poli come parametro (ad esempio, guadagno) varia.
- Trama di boda e diagrammi di Nyquist:[ Fornire margine di guadagno e margine di fase, che sono indicatori di stabilità relativa e caratteristiche di risposta transitoria.
- Step Response Simulation:[] Fornisce le metriche transitorie effettive (overshoot, tempo di configurazione) per un determinato sistema.
- Metodi di stato-spazio:[] Gli autovalori della matrice di stato (che sono i poli) possono essere calcolati numericamente, e il metodo di Lyapunov estende l'analisi di stabilità ai sistemi non lineari.
Il criterio Routh-Hurwitz occupa quindi un punto dolce storico e pratico: è abbastanza matematico da essere rigoroso, ma abbastanza semplice da fare su una lavagna bianca durante una recensione di design.
Conclusioni
Il criterio Routh-Hurwitz offre un modo chiaro e algebrico per valutare la stabilità del sistema analizzando il polinomio caratteristico senza risolvere le radici. Attraverso la costruzione del array Routh, gli ingegneri possono determinare il numero di poli instabili, comprendere la sensibilità del sistema ai cambiamenti dei parametri, e identificare i confini di simulazione della stabilità che influiscono direttamente sul comportamento transitorio.
Il software Routh-Hurwitz fornisce inoltre un controllo di alta prestazione, che consente di valutare rapidamente se un sistema di controllo si comporta in modo affidabile in condizioni transitorie.