Gebruik van Routh-hurwitz om systeemrespons en tijdelijk gedrag te voorspellen

Inleiding tot het criterium van Routh-Hurwitz

Het criterium Routh-Hurwitz is een hoeksteen van de klassieke controletheorie, die een systematische methode voor het beoordelen van de stabiliteit van lineaire tijd-invariant (LTI) systemen zonder expliciet de wortels van de karakteristieke vergelijking berekenen. Ontwikkeld onafhankelijk van Edward John Routh in 1874 en Adolf Hurwitz in 1895, deze algebraïsche techniek blijft wijd gebruikt in de controle engineering, omdat het biedt een snelle, betrouwbare manier om te bepalen of een systeem zal reageren op een stabiele, voorspelbare manier of vertonen oscillaties en instabiliteit. Stabiliteit is fundamenteel voor het ontwerp van een systeem een instabiel systeem kan leiden tot gevaarlijke ondoordringbaarheid, aanhoudende oscillaties, of zelfs catastrofale mislukking, vooral in toepassingen zoals lucht- en ruimtevaart, auto-besturing, chemische processen en robotica.

Door de coëfficiënten van een systeem te analyseren, wordt het aantal polen (wortels van de karakteristieke vergelijking) dat in de rechterhelft van het complexe vlak ligt, geïdentificeerd. Polen met positieve reële delen veroorzaken onstabiel gedrag, terwijl polen met negatieve reële delen bijdragen tot stabiele, convergente reacties. De methode vereist geen oplossing van polynomen van hoge mate, waardoor het praktisch is voor handmatige analyse en vroeg-fase ontwerp. Bovendien zijn inzichten in voorbijgaande gedragspatronen zoals zettingstijd, .. en oscillion frequentie helpen ingenieurs met het afstemmen van controllers om te voldoen aan prestatiespecificaties voordat ze zich verbinden aan gedetailleerde simulaties of hardware-prototypes.

Dit artikel zal eerst de wiskundige basis van het Routh-Hurwitz criterium uitleggen, en u vervolgens begeleiden door de Routh-array te bouwen met een stap-voor-stap voorbeeld. Daarna zullen we onderzoeken hoe het criterium systeemrespons en transiënt gedrag voorspelt, inclusief belangrijke kenmerken zoals overloop en afwikkelingstijd. Een praktisch numeriek voorbeeld zal het proces illustreren, en we zullen de beperkingen en aanvullende hulpmiddelen van de methode bespreken. Aan het eind zult u een solide inzicht hebben in hoe u het Routh-Hurwitz-criterium als praktisch, gezaghebbend instrument in de analyse van het controlesysteem kunt toepassen.

Wat is het criterium van Routh-Hurwitz?

Het criterium Routh-Hurwitz is een noodzakelijke en voldoende voorwaarde voor stabiliteit in LTI-systemen. Het werkt op de karakteristieke vergelijking van het systeem, die wordt afgeleid van de noemer van de functie van de overdracht van gesloten lus nadat het op nul is ingesteld. Voor een systeem beschreven door de overdrachtsfunctie is de karakteristieke vergelijking , waarbij een polynoom is in s graad n:

Het criterium stelt dat voor het systeem stabiel moet zijn (dat wil zeggen dat alle polen negatieve reële delen hebben), dat a_i niet nul moeten zijn en hetzelfde teken (typisch positief) moeten hebben, en dat de Routh-array die is samengesteld uit deze coëfficiënten geen tekenveranderingen in de eerste kolom mag hebben. Als er een teken verandert, is het aantal tekenveranderingen gelijk aan het aantal polen met positieve reële delen, wat wijst op instabiliteit.

Het is belangrijk om op te merken dat het criterium Routh-Hurwitz zuiver algebraïsch is, maar dat het niet de werkelijke poollocaties geeft. In plaats daarvan wordt een aantal palen van rechts-half-vlak (RHP) gegeven. Dit maakt het onschatbaar wanneer de karakteristieke polynomial van hoge orde is (bijvoorbeeld vierde, vijfde of hoger) waar wortelvinding computationeel of analytisch niet haalbaar is.

Het criterium is direct gekoppeld aan de Hurwitz matrix, een vierkante matrix gevormd uit de coëfficiënten van de karakteristieke polynomial. Alle leidende belangrijkste minderjarigen van de Hurwitz matrix moet positief zijn voor stabiliteit dit is de Hurwitz voorwaarde. Routh... › › methode vereenvoudigt dit determinant-checking proces in de meer eenvoudige Routh array.

Het opbouwen van de Routh Array: Stap-voor-stap

De bouw van de Routh-array is een systematische procedure. We zullen het illustreren met behulp van een concreet voorbeeld polynomial:

Deze veelterm heeft wortels bij s = -1, -2, -3, -4 (alle negatieve), dus het systeem is stabiel. We zullen het gebruiken om te controleren of de Routh array geen tekenwijzigingen geeft.

Stap 1: Schrijf de Polynomial in standaardvorm

Zorg ervoor dat de polynoom is gerangschikt in afnemende vermogens van s met alle coëfficiënten geschreven. Voor ons voorbeeld, de coëfficiënten zijn:[
a₄ = 1 (s^4[), a[_3[ = 10 (s^3[), a_2[] = 35 (s[^2[]), a[₁ = 50 (s[[¹), a_0[ = 24 (constant term).

Stap 2: Initialiseer de eerste twee rijen

De eerste rij van de Routh-array bevat de coëfficiënten van evengeïndexeerde termen (beginnend van s[ⁿ): rij s^4[: [1, 35, 24]. De tweede rij bevat de coëfficiënten van oneven-geindexeerde termen: rij s^3[: [10, 50, 0]. (Vul ontbrekende termen in met nullen.)

Dus:
Rij s⁴: 1, 35, 24
Rij s³: 10, 50, 0

Stap 3: Bereken volgende rijen met behulp van de Rhombus regel

Elk element in de volgende rij wordt berekend met behulp van de twee rijen erboven. Voor een element in rij s², eerste kolom, is de formule:

Maar het is makkelijker om de standaardformule toe te passen: voor rij k, kolom j (beginnend vanaf kolom 1), is het element:

Laten we het handmatig doen voor duidelijkheid.

Row s[²:[
Eerste element: gebruik van rijen s[^4[] en s[^3[.
R(s[^2[]]1) = -[ (1×50) - (10×35) ] / 10 = -[50 - 350]/10 = -[-300]/10 = 30.[
Tweede element: volgende kolom (j=2). R(s[⁴],2) = 35, R(s^{,3]=24; R(s³]]]]]] = 50][LT:[20][LT][}

Zo luidt rij s²: 30, 24, 0.

Roeien ¹:
Gebruik de rijen s^3[] en s[^2[].[
R(s^1[],1) = -[ (10×24) - (30×50) ] / 30 = -[240 - 1500]/30 = -[[SLT:1260]/30 = 42.
R(s^1[,2) = -[ (10×0) - (30×0) ] ] ] / 30 = 0.

Aldus rij s¹: 42, 0.

Roemen⁰:
Gebruik rijen s² en s[^1[]
R(s[^0[],1) = -[ (30×0) - (42×0) ] ] / 42 = 0. Maar de laatste rij is altijd de constante term van de polynomial, die 24 zou moeten zijn. Eigenlijk is er een bekende regel: het laatste element van de laatste rij is gelijk aan a₀ (de constante term). We kunnen berekenen als: R(s⁰],1) = -[ (30×24] - (42×) ] ] ] - (42 = [[-100]] ] = 24] = 24] ] ] [[ [FLT

Dus de complete Routh array is:

s⁴	1	35	24
s³	10	50	0
s²	30	24	0
s¹	42	0	0
s⁰	24	0	0

Stap 4: Onderzoek de eerste kolom

De eerste kolom bevat: 1, 10, 30, 42, 24. Allen zijn positief, en er zijn geen tekenen veranderingen. Daarom is het systeem stabiel, die overeenkomt met onze bekende wortels.

Bijzondere zaken

Twee bijzondere gevallen kunnen zich voordoen tijdens de bouw van een Routh-array:

Zero in de eerste kolom: Als het eerste element van een rij nul is maar de rij niet alle nullen is, vervang dan de nul door een klein epsilon (ε), ga dan verder. Het teken ε bepaalt de tekenwijzigingen. Dit geeft de marginale stabiliteit of instabiliteit aan, afhankelijk van de rest van de reeks.
Entire rij van nullen: Dit gebeurt wanneer er symmetrisch geplaatste wortels (bijvoorbeeld paren van imaginaire wortels of symmetrische paren rond de echte as) zijn. In dit geval, vorm een hulppolynomial van de rij boven de nul rij en differentieer het om coëfficiënten te verkrijgen voor de ontbrekende rij. Dit geeft de aanwezigheid aan van een paar puur denkbeeldige wortels of een paar echte wortels met tegengestelde tekens (d.w.z. instabiliteit).

Deze speciale gevallen breiden de toepasbaarheid van het criterium uit tot systemen met polen op de denkbeeldige as of symmetrische patronen.

Voorspelling van systeemrespons en tijdelijk gedrag

Het criterium Routh-Hurwitz is niet beperkt tot een binaire stabiele/onstabiele uitspraak. Door het aantal polen met positieve reële delen te onthullen, informeert het direct het voorbijgaande gedrag van het systeem. Elke pool in het rechter halfvlak draagt bij aan een exponentieel groeiende modus, waardoor de voorbijgaande respons uiteenvalt. Omgekeerd leiden alle polen in het linker halfvlak tot verval exponentieel, uiteindelijk tot een stabiele toestand.

Het criterium geeft echter niet de exacte poollocaties, zodat het niet precies kan kwantificeren overschrijding, afwikkelingstijd, of stijgingstijd. Niettemin, het geeft essentiële aanwijzingen:

Als de eerste kolom geen tekenwijzigingen heeft: Alle polen hebben negatieve echte delen. De voorbijgaande respons zal stabiel zijn, uiteindelijk zich vestigend op een steady state. Maar er kan een traag of oscillerend verval (bijna denkbeeldige polen) zijn dat het criterium niet alleen het teken van het echte deel opspoort, is gegarandeerd.
Als er tekenen veranderen: Het aantal tekenveranderingen is gelijk aan het aantal RHP-polen. Het systeem is onstabiel en de voorbijgaande respons zal ongebonden groeien. Het aantal RHP-polen wijst op de aard van instabiliteit: een enkele positieve pool leidt tot monotone exponentiële groei, terwijl een paar complexe RHP-polen groeiende oscillaties veroorzaken.
Marginale stabiliteit: Als de karakteristieke vergelijking wortels heeft op de denkbeeldige as (zuiver imaginaire of meerdere polen op de as), zal de Routh-array een rij nullen produceren. Dit geeft aan dat het systeem op de grens van stabiliteit staat, en de voorbijgaande respons kan aanhoudende oscillaties (indien eenvoudige denkbeeldige polen) of oplopende oscillaties (indien herhaalde denkbeeldige polen) vertonen.

Door het resultaat van Routh-Hurwitz te combineren met kennis van de systeemvolgorde, kunnen ingenieurs geïnformeerde voorspellingen doen. Voor een tweede-ordesysteem vermindert het criterium tot het vereisen van alle coëfficiënten om positief te zijn (wat zowel noodzakelijk als voldoende is voor stabiliteit).Voor hogere-orde systemen, de array biedt een stabiliteitsmarge indicator .Hoe dichter de eerste-kolom elementen zijn op nul, hoe dichter het systeem is instabiele. Dit is vooral handig bij het aanpassen van de winstparameters in een controlelus.

Relatie met voorbijgaande eigenschappen

De voorbijgaande respons van een stabiel systeem wordt voornamelijk bepaald door de dominante polen ..die het dichtst bij de denkbeeldige as. Het criterium Routh-Hurwitz niet direct lokaliseren deze polen, maar het kan helpen bij het afstemmen . Bijvoorbeeld , als het criterium toont dat het verhogen van een controller winst uiteindelijk zal leiden tot een teken verandering in de eerste kolom , markeert het de winst marge . Ingenieurs vaak dit gebruiken om proportionele winsten in P , PI , of PID controllers zonder dat het hele systeem te simuleren elke keer .

Belangrijke voorbijgaande metrieken die door stabiliteit worden beïnvloed zijn:

Overshoot . . De maximale piekwaarde van de respons ten opzichte van de uiteindelijke waarde. Overshoot wordt meestal veroorzaakt door complexe geconjugeerde polen met een lage dempingsverhouding. Terwijl het Routh-Hurwitz criterium geen demping berekent, komt een systeem met een rij nullen (marginale stabiliteit) vaak overeen met een dempingsverhouding bij nul, wat een hoge overbelasting impliceert.
Insteltijd
Begintijd
Oscillaties .. De frequentie van voorbijgaande oscillaties wordt bepaald door de denkbeeldige delen van de complexe polen. Het criterium onthult geen denkbeeldige delen, maar de aanwezigheid van een rij nullen duidt puur denkbeeldige polen aan, die duurzame oscillaties veroorzaken in de voorbijgaande respons.

Het criterium Routh-Hurwitz fungeert dus als poortwachter: het zorgt voor stabiliteit als een voorwaarde voor zinvolle transiënte analyse met behulp van andere instrumenten zoals wortellocus, Bode-ploegen of stapresponssimulaties.

Praktisch voorbeeld: Derde-ordesysteem

Laten we het Routh-Hurwitz criterium toepassen op een complexer voorbeeld met onbekende parameters. Beschouw de closed-loop transfer functie met eenheid feedback. De karakteristieke vergelijking is:

We willen het bereik van K bepalen waarvoor het systeem stabiel is. Teken de Routh-array voor polynoom .

Rijen:
s^3[: 1, 10
s²: 7, K[
s[^1[: [^{s[^]0[]: [} Eigenlijk eenvoudiger: het element van de laatste rij is de constante term K (mits de rij hierboven niet alle nullen is). Meer systematisch: rij s[^0[ is afgeleid uit rijen s²] en s¹]. De formule die wordt gebruikt: ^{(het is de constante term rechtstreeks gelijk aan de array).}

s³	1	10
s²	7	K
s¹	(K-70)/7	0
s⁰	K	0

Voor stabiliteit moeten alle elementen van de eerste kolom positief zijn: 1 > 0, 7 > 0, (K-70)/7 > 0, K > 0. Dit geeft K > 70 en K > 0, dus K > 70. De ondergrens is K > 70. Als K = 70, de rij s^1[] nul wordt, wat een speciaal geval is dat marginale stabiliteit aangeeft (zuiver denkbeeldige wortels). Dus het systeem is stabiel voor K > 70. Dit vertelt de ingenieur dat de winsten boven 70 een stabiel voorbijgaand resultaat opleveren, maar als K 70 nadert, worden de tijdelijke oscillaties onaangetast. De tijdelijke respons voor K = 200, bijvoorbeeld, zal twee complexe geconjugeerde polen in het linker halve vlak en een echte pool hebben, wat leidt tot een overbruggende maar uiteindelijk bestendig reageren. Het Routh-Hurwitz criterium hier expliciet voorspelt de stabiliteitsgrens.

Voorbijgaande responsimplicaties voor het voorbeeld

Voor K net boven de 70 is het systeem stabiel, maar de dominante polen zijn zeer dicht bij de denkbeeldige as. De voorbijgaande reactie zal zeer langzaam verval en grote overschrijding vertonen. De afwikkelingstijd zal lang zijn. Als K verder toeneemt, bewegen de polen naar links (in het algemeen), verbeteren de demping en de afwikkelingstijd. Echter, als K te hoog is, zou de Routh-array voor dit specifieke systeem geen extra tekenenveranderingen vertonen. In de praktijk kan zeer hoge K andere problemen veroorzaken zoals actuatorsaturatie of geluidsversterking, maar het Routh-Hurwitz-criterium controleert alleen de stabiliteit van de lus.

Het criterium geeft dus kritische informatie: de winst voor stabiliteit, en de nabijheid van instabiliteit. Door de eerste kolom elementen te evalueren, kunnen ingenieurs winstmarges en ontwerpcontrollers berekenen die marginale stabiliteit vermijden.

Beperkingen van het criterium van Routh-Hurwitz

Ondanks zijn kracht heeft het criterium van Routh-Hurwitz belangrijke beperkingen die elke control engineer moet begrijpen:

Alleen van toepassing op LTI-systemen: Het kan niet direct worden gebruikt voor tijdvariabel, niet-lineair of bemonsterde gegevenssystemen (hoewel er uitbreidingen bestaan voor discrete tijdsystemen via bilineaire transformaties).
Kwantificeert geen transiënte prestaties: Het criterium geeft geen informatie over overschrijding, piektijd, stijgingstijd of exacte afwikkelingstijd. Daarvoor zijn root locus, frequentierespons of tijddomeinsimulaties nodig.
Vereist een bekende karakteristieke polynoom: De coëfficiënten moeten precies bekend zijn. In de praktijk kunnen onzekerheden in parameters een robuuste stabiliteitsanalyse vereisen (bijvoorbeeld door Kharitonovs stelling te gebruiken).
Wordt omslachtig voor zeer hoge-orde systemen: Handmatige berekeningen voor tiende-orde polynomen zijn foutgevoelig, hoewel software (MATLAB, Python) automatiseert de array constructie.
Kan vertragingen niet direct verwerken: Systemen met vertraging hebben transcendentale kenmerkende vergelijkingen.De Routh-array kan niet direct worden toegepast tenzij Padé-amperamenten worden gebruikt.

Om deze beperkingen te overwinnen, gebruiken ingenieurs vaak het Routh-Hurwitz criterium in combinatie met andere methoden. Bijvoorbeeld, wortellocaties tonen hoe gesloten-luspalen bewegen met winst, het verstrekken van exacte voorbijgaande informatie. Het Nyquist criterium behandelt vertragingen en niet-lineairheden. Niettemin dient de Routh-array vaak als een snelle eerste controle . vooral tijdens het voorlopige ontwerp of probleemoplossing.

Aanvullende hulpmiddelen

Het criterium Routh-Hurwitz is ..onnodig en voldoende ..voor de stabiliteit van LTI-polynomen, maar voor prestatieanalyse, ingenieurs vertrouwen op:

Root Locus: Plots de baan van polen als een parameter (bijv., gain) varieert. Het toont direct de dempingsverhouding en natuurlijke frequentie.
Bode Plots en Nyquist Diagrams: Zorg voor winstmarge en fasemarge, die indicatoren zijn voor relatieve stabiliteit en tijdelijke responskenmerken.
Stapresponsimulatie: Geeft de werkelijke tijdelijke metriek (overshoot, afwikkelingstijd) voor een bepaald systeem.
State-Space Methoden: De eigenwaarden van de staatmatrix (die de polen zijn) kunnen numeriek worden berekend, en Lyapunovs methode breidt stabiliteitsanalyse uit tot niet-lineaire systemen.

Het criterium Routh-Hurwitz neemt dus een historische en praktische zoet plekje in beslag: het is wiskundig genoeg om rigoureus te zijn, maar toch eenvoudig genoeg om tijdens een ontwerponderzoek op een whiteboard te worden gedaan. Het blijft een essentiële vaardigheid voor elke control engineer.

Conclusie

Het Routh-Hurwitz criterium biedt een duidelijke, algebraïsche manier om de stabiliteit van het besturingssysteem te beoordelen door de karakteristieke polynomiale analyse zonder op te lossen voor wortels. Door de bouw van de Routh-array kunnen ingenieurs het aantal onstabiele polen bepalen, de gevoeligheid van het systeem voor parameterveranderingen begrijpen en stabiliteitsgrenzen identificeren die direct van invloed zijn op tijdelijk gedrag. Hoewel het geen exacte tijdelijke metrieken zoals ondoordringbare of fixatietijd biedt, dient het als een basisscreeningsinstrument dat diepere analyse informeert met behulp van root locus, frequentierespons of simulatie.

De beheersing van het Routh-Hurwitz-criterium biedt u de mogelijkheid om snel te beoordelen of een besturingssysteem zich betrouwbaar zal gedragen onder voorbijgaande omstandigheden. Of u nu een eenvoudige proportionele controller ontwerpt of een complex multi-loopsysteem afstemt, de Routh-array is een betrouwbaar startpunt. Voor verdere lezing, overwegen om het [Routh