Table of Contents
Introduzione all'analisi della stabilità nei sistemi discreti
L'analisi della stabilità costituisce la colonna portante dell'ingegneria del sistema di controllo, assicurando che un sistema risponda prevedibilmente e rimanga limitato alle normali condizioni operative. Per i sistemi discreti-tempo—comune nel controllo digitale, nell'elaborazione del segnale e nelle applicazioni incorporate—il criterio di stabilità si sposta dal piano continuo s-plano sinistro al cerchio dell'unità nel piano z-plano.
Questo articolo fornisce un confronto completo di Routh-Hurwitz e Jury’s test di stabilità, esplorando le loro basi matematiche, i passaggi procedurali, gli adattamenti pratici e i punti di forza relativi.
Stabilità del sistema discreto: Criteri del cerchio dell'unità
Prima di immergersi nei metodi specifici, è essenziale ricordare la condizione fondamentale di stabilità per i sistemi discreti-time. Un sistema discorsi lineare tempo-invariante descritto dalla funzione di trasferimento H(z) = ]]N(z)]]] / ]D(z) è solo
Mentre gli algoritmi di ricerca delle radici possono calcolare direttamente le posizioni dei poli, diventano ingombranti per i polinomi di alto ordine. Quindi, i criteri algebrici come Routh-Hurwitz e Jury’ la prova offre controlli efficienti e non iterative. La sfida chiave è che Jury’s test funziona direttamente sui coefficienti di trasformazione D(znomi)[FFFFFurwiur]
Criteri di Routh-Hurwitz e la sua adattamento ai sistemi discreti
Routh-Hurwitz classico per sistemi continui
Il criterio Routh-Hurwitz fornisce condizioni necessarie e sufficienti per la stabilità dei sistemi lineari continui.
[LT] [[FLT]] [[FLT]]]] ]] ]] ]] ]]]][FLT]]]][FLT][F82]
Il sistema è stabile se tutti gli elementi nella prima colonna dell'array sono positivi e nonzero. Un cambio di segno indica le radici nel piano di metà destra e quindi l'instabilità. Questo approccio evita il calcolo esplicito della radice ed è diretto per i polinomi fino a ordini moderati.
Trasformazione bilineare: conversione dello z-Plane al s-Plane
Per applicare Routh-Hurwitz a sistemi discreti, gli ingegneri utilizzano la trasformazione bilineare (chiamata anche Tustin’s transformation). Questa mappatura riguarda la variabile [z] a una nuova variabile complessa ]w]] (spesso indicata come ]] [[
z[] = (1 + wT[/2] / (1 – ]wT/2]
[FLT]] [[FLT]]] [[FLT]]] = 2, la forma più semplice ]]] [ = (1 + ]w]]]]]] [[[]]]]]]]]] [[FLT[FLT]]]]]]]]]]]]]]]]]]]] [[[[[[[[[[[[[[[[[FLT]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]] [[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[
Passi per l'applicazione di Routh-Hurwitz a sistemi discreti
- Ottenga il polinomio caratteristico discreto[[ []D(z)] dal sistema & n. 8217; la funzione di trasferimento.
- Sostituire la trasformazione bilineare[[] in []D(z)] per ottenere un polinomio in []w]]. Questo passaggio può comportare l'espansione e la semplificazione algebrica, soprattutto per i sistemi di ordine più alto.
- Applicare il criterio Routh-Hurwitz[ al polinomio trasformato in [w[]. Costruisci il array Routh ed esamina i segni delle voci di primo livello.
- Conclusione di stabilità:[] Se tutte le voci di primo livello sono positive, il sistema discreto è stabile (tutti i poli all'interno del cerchio dell'unità).
Mentre questo approccio funziona, introduce una maggiore complessità algebrica. La trasformazione bilineare può aumentare significativamente l'ordine polinomiale se non gestito con attenzione. Inoltre, la trasformazione stessa può diventare una fonte di errore numerico, soprattutto per sistemi di alto ordine o quando i periodi di campionamento non sono normalizzati correttamente. Nonostante questi svantaggi, il metodo Routh-Hurwitz rimane un'opzione praticabile per gli ingegneri già familiari con la sua applicazione continua e che preferiscono un unico quadro.
Limitazioni dell'adattamento Routh-Hurwitz
- L'onere algebrico:[ La trasformazione bilineare espande il polinomio, potenzialmente introducendo errori di arrotondamento e facendo calcoli manuali noiosi.
- Sensibilità numerica:[ Per i polinomi di alto ordine, le magnitudine del coefficiente possono variare ampiamente, portando a array ill-condizionato.
- Casi speciali:[] Il metodo richiede la manipolazione di polinomi ausiliari per le radici sull'asse immaginario (nell' []w[]]-plano) che corrispondono ai poli sul cerchio unitario – una condizione di stabilità marginale che è spesso indesiderato.
- Inefficienza per sistemi puramente discreti:[] Poiché il metodo non è nativo del dominio discreto, molti libri di testo di controllo raccomandano Jury’s test come l'alternativa più diretta e meno incline all'errore.
Tuttavia, il criterio Routh-Hurwitz viene talvolta insegnato nei corsi di controllo come ponte tra analisi continua e discreta, e può essere applicato utilizzando strumenti di calcolo simbolici per mitigare le difficoltà algebriche.
Test di stabilità della giuria: progettato per il dominio discreto
Origini e Principio
Giury’s test di stabilità, introdotto da Eliahu I. Giuria negli anni '60, è la controparte discreta del criterio Routh-Hurwitz. Esamina direttamente i coefficienti del polinomio caratteristico nomi]D(z)] senza richiedere trasformazioni di dominio.
Considerare il polinomio caratteristico minico dell'ordine n:
D(z) = a0zn + a1zn–1 + … + an–1z + an = 0, with a0 > 0.
Jury’s test fornisce una procedura sistematica per determinare se tutte le radici si trovano all'interno del cerchio dell'unità formando una tabella con [n[] + 1 righe. Le righe vengono generate sequenzialimente utilizzando calcoli determinanti che coinvolgono i coefficienti.
Procedura per la costruzione della tabella di giuria
- Imposta i coefficienti[] nella prima riga (row 0) in ordine naturale: a]0, ]]1[FLT]
- [[LT]][[FLT]]][[FLT]]]][[[FLT]]]]][[FLT:]]][[[FLT:]]]][[] [[FLT:[FLT]]]][[[[[FLT]]]]]][[[FLT]]]][[[[FLT]]]]]]]]]]][[[[[[[[[[[[FLT]]]]]]]]]]]]]]]]]]][[[[[[[[[[[[[[[[[[[[[[FLT]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]][[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[
- [LT][LT][LT][[[[[]]]][LT][[[[[[]]]]]][[[[[[[[[]]]]]]]][[[[[LT]]]]][[[[[[[[[]]]]]]]]]]]]]][[[[[[F]]]]]]]]]]]]]]]]]]]]]]]]][[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[F]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]
- Continua fino a quando non rimane un solo elemento non zero[[ nella prima colonna, o fino al crollo della tabella. La riga finale fornisce l'ultima condizione.
- Verificare le condizioni:[] Per la stabilità, il seguente deve contenere: [
- ]D](1) > 0
- ] [ ] ] []]D[[](–1) > 0
- Per la tavola della Giuria, il primo elemento di ogni fila numerata (crescite 1, 3, 5, ...) deve essere maggiore di zero.
Queste condizioni sono necessarie e sufficienti per tutte le radici per soddisfare |[]z| < 1. In pratica, il test viene implementato nei sistemi di algebra del computer, ma il processo manuale è gestibile per ordini fino a circa 4 o 5.
Esempio: Sistema di Seconda Ordina
[FLT]] [[FLT]]]] ]] ]]]][FLT]][FLT]][F[F[FLT]][FLT]][[[FLT]]][[[[FLT][[[[FLT]]]]]][[[[[[FLT]]]]][[[[F[F[F[[[F]]]]]]]][[[FLT]]]]]]][[[FLT]]]]]]][[[F[F[[[[[[[[[[[FLT]]]]]]]]]]]]]][FLT]]]]]]]]]]]][[[[[[[[[[[[[[
- Row 0: 1, ]a][1, ]]a[]]2]
- Raccolta 1:] ]]]], ]]]]]]1, 1
- [LT] [FLT] [[[6]]] [[f]]] [[f]]]] [[fl]]]] [[f]]]][[[f]]]][[[fl]]]][[[f]]]]]]][[f]]]]][[[f]]]]]][FLT]][f]][f]]][[[f]]]]]][[[[[[[[f]]]]]]]]]]]]]]]]]]]][[[[[[[[[[[f]]]]]]]]]]]]]]]]]][[[[[[[[[[[[[[[[[[[f]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]
Condizioni di stabilità per il secondo ordine:
- D[](1) = 1 + ]a]1 + ]a[]2]]] > 0
- ] ] ] ]]] ]]] ]]]]]]]] + [FLT:]]]]]]]]][FLT:[FLT]]]]][FLT:[FLT]]]]]][FLT][FLT]][FLT:[FLT:[FLT]]]]]]]]][FLT:[FLT:[FLT]]]][FLT:[FLT]]]]]][F[FLT]]]]]]][FLT][FLT]][
- a][2[]] | < 1
Queste sono le condizioni familiari per un sistema stabile di secondo ordine discreto. Si noti che Jury’s test produce direttamente queste semplici disuguaglianze senza alcuna trasformazione.
Vantaggi della Jury’s Test
- Applicazione diretta:[]] Funziona sul z]-dominio polinomiale; non è necessario mappatura.
- Set completo di condizioni:[ Fornisce sia criteri di stabilità necessari che sufficienti.
- Casi speciali di manegge:[]] Include controlli per le radici sul cerchio unitario (stabilità marginale) attraverso il []D[[]]](1) e [D[](–1) condizioni.
- Adatto all'implementazione digitale: La costruzione della tabella è algoritmica e facilmente programmata.
Confronto testa a testa di Routh-Hurwitz e Test di stabilità della giuria
| Aspect | Routh-Hurwitz (Adapted) | Jury’s Stability Test |
|---|---|---|
| Domain | Continuous (s-plane) originally; adapted via bilinear transform to w-plane | Discrete (z-plane) directly |
| Stability region | Left-half plane (after transformation, maps to inside unit circle) | Inside unit circle (|z|<1) |
| Procedure | Construct Routh array from w-polynomial; check sign of first column | Construct Jury table from z-polynomial; check series of inequalities |
| Complexity | Moderate after transformation; transformation adds extra algebraic steps | Low to moderate; table construction straightforward for low to medium orders |
| Numerical stability | Can suffer from coefficient scaling after bilinear transform | Generally good; determinants may cause swelling for high orders |
| Handling of marginal stability | Auxiliary polynomial required; roots on imaginary axis in w-plane | Direct checks through D(1) and D(-1) conditions |
| Educational accessibility | Students often learn continuous version first; adaptation extends existing knowledge | Requires understanding of unit circle concept; but method is self-contained |
| Software implementation | Easily coded via array construction; bilinear transform may require symbolic toolbox | Easily coded via recursive table; native to z-domain |
| Typical use case | When system is already in s-domain or hybrid continuous/discrete | When characteristic polynomial is directly from a discrete transfer function |
Punti di confronto dettagliati
1. Applicabilità e Dominio Fit
La prova è la scelta naturale per sistemi discreti puri, soprattutto nelle applicazioni di controllo digitale dove vengono campionati l'uscita del controller e l'ingresso delle piante. Routh-Hurwitz richiede una trasformazione che aggiunge complessità e può oscurare l'interpretazione fisica dei poli rispetto al cerchio dell'unità.
2. Complesso algebrico
Per ordini fino a 4 e 5, entrambi i metodi sono gestibili a mano. Tuttavia, Jury’s test evita l'espansione bilineare, che può produrre termini con potenze miste e grandi coefficienti. Ad esempio, un terzo ordine D(z)] trasformati a ]w produce un coefficiente di terzo ordine in sé stesso
3. Robustezza numerica
Routh-Hurwitz nel w[]]-dominio può soffrire di malcondizionamenti se la trasformazione bilineare crea coefficienti con magnitudine molto diverse. Jury’s test affronta anche problemi numerici per polinomi di alto ordine (ad esempio, ordine > 10) a causa di sottratti ripetuti e calcoli determinanti, ma questi sono meno gravi.
4. Capacità di determinare il numero di polacchi non stabili
Routh-Hurwitz, in virtù della Routh array’s significa cambiamenti, può indicare direttamente il numero di poli nella destra-half [w[]]-plano, che corrisponde al numero di poli al di fuori del cerchio dell'unità.
5. Gestione di Coefficienti Polinomiali Speciali
Quando un sistema discreto ha zero coefficienti (poteri di licenziamento), entrambi i metodi richiedono cautela. Routh-Hurwitz tratta zero voci nella matrice utilizzando piccole approssimazioni di epsilon. Jury’s test ha bisogno di gestire gli zeri durante la costruzione della tabella, ma l'algoritmo rimane valido fino al coefficiente di riferimento a][FLT: stabilità dis]]]
Applicazioni pratiche e quando scegliere ogni metodo
Utilizzare la prova della giuria quando
- Lavori esclusivamente nel dominio discreto (z-plane).
- Vuoi un test a una fase senza trasformazioni di dominio.
- Stai analizzando un controller digitale o un sistema di dati campionati con un periodo di campionamento noto.
- È necessario condizioni chiare per la stabilità marginale (radici sul cerchio dell'unità).
- L'ordine polinomiale è moderato (≤ 8 per i controlli manuali; più alto per il software).
Utilizzare Routh-Hurwitz (via Bilinear Transform) Quando
- Il tuo sistema di controllo è già nel dominio continuo e si sta discostando in seguito.
- Preferite un quadro di analisi unificato per sistemi continui e discreti.
- Hai software (ad esempio, MATLAB, Python con SymPy) per gestire automaticamente la trasformazione bilineare.
- È necessario contare il numero esatto di pali instabili (tramite le modifiche del array Routh).
- State insegnando il legame tra criteri di stabilità continui e discreti.
Esempi reali-mondo
Esempio: Regolatore digitale della velocità del motore. Un microcontroller implementa un controller PI per un motore DC. Il polinomio caratteristico a ciclo chiuso deriva dal dominio z. Un ingegnere utilizzerebbe Jury’s test per verificare la stabilità per vari tassi di campionamento e guadagni proporzionali.
Esempio: Progettazione del filtro adattivo. Nel processo di elaborazione del segnale, i filtri IIR richiedono pali all'interno del cerchio dell'unità per evitare l'autooscillazione.
Esempio: Sistema ibrido continuo-discreto. Un impianto continuo controllato da un compensatore digitale richiede discretizzazione tramite una presa di ordine zero. L'equazione caratteristica generale nel domini z può essere analizzata sia da Jury’s test direttamente o utilizzando la trasformazione bilineare per applicare Routh-Hurwitz. Molti ingegneri di controllo preferiscono il soggiorno Jury&82#
Attuazione e automazione del software
Entrambi i metodi sono disponibili nelle caselle di controllo del sistema.
- MATLAB:[] La funzione (o script personalizzati) può costruire l'array Routh, mentre è disponibile nella Toolbox Robotics o può essere programmato. La funzione esegue la discorsizzazione e la trasformazione bilineare se necessario.
- Python (SciPy and control library): La libreria include per Jury’s test. Per Routh-Hurwitz, viene utilizzata la funzione o la costruzione manuale dell'array.
- Strumenti simbolici:[ SymPy e Mathematica permettono la manipolazione algebrica per la trasformazione bilineare e la successiva generazione di array.
Quando si utilizza il software, il passaggio di trasformazione per Routh-Hurwitz è automatizzato, riducendo l'errore manuale. Tuttavia, per la prototipazione rapida dei sistemi discreti, Jury’s test rimane la scelta più diretta.
Limitazioni e Caveats
Entrambi diventano inflessibili per i polinomi di alto ordine (ordine > 10) a causa dell'accumulo di errori numerici. Per tali casi, la radice-finanza computazionale (ad esempio, metodi di autovalore) è spesso più affidabile. Inoltre, entrambi i test assumono la caratteristica polinomiale è noto esattamente; modellare incertezze e variazioni dei parametri richiedono una solida analisi della stabilità al di là di questi criteri algebrici.
Un altro importante avvertimento: l'adattamento Routh-Hurwitz tramite la trasformazione bilineare è valido solo per sistemi lineari invarianti, non si estende direttamente ai sistemi discorsivi o non lineari.
Per la Jury’s test, la condizione D](1) > 0 e (–1)]]]]][
Conclusioni
La scelta tra Routh-Hurwitz e Jury’ la prova di stabilità per sistemi discreti dipende dal contesto, dalla familiarità e dalle risorse computazionali. Jury’s test è il metodo più diretto e costruito appositamente per polinomi a tempo discreto, offrendo condizioni chiare ed evitando la trasformazione bilinea del rinforzo.
Per la maggior parte del moderno design di controllo digitale e l'analisi del sistema discreta, Jury’ la prova di stabilità è l'approccio raccomandato. Tuttavia, la competenza in entrambe le tecniche assicura che l'ingegnere può adattarsi a varie dichiarazioni di problemi e sfruttare il software disponibile per verificare i risultati.