Het Routh .Hurwitz stabiliteitscriterium is een hoeksteen van de klassieke controle theorie, het verstrekken van een directe algebraïsche methode om te bepalen of een lineair tijd-invariant systeem stabiel is zonder op te lossen voor zijn polen. Voor integer-order systemen .Waar de karakteristieke polynomial heeft alleen gehele krachten van de Laplace variabele s] .De test is zowel efficiënt als elegant . Echter , een groeiend aantal van de reële toepassingen , van viscoelastische demping tot elektrochemische impedantie spectroscopie , worden beschreven door fractionele-order differentiaalvergelijkingen . Deze systemen introduceren derivaten van niet-integreerde orde , en hun karakteristieke vergelijkingen bevatten fractionele krachten van , waardoor de traditionele Routh array niet direct kan worden samengesteld . Engineers en onderzoekers moeten daarom de onderliggende filosofie van de Routh .Hurwitz analyse aanpassen om fractionele-orde dynamicatica te hanteren . Dit artikel legt uit hoe deze ideeën, die praktische methoden die de eenvoud van het oorspronkelijke

Het klassieke criterium van Routh

Voordat u extensies onderzoekt, is het essentieel om de klassieke methode te herdefiniëren.

Het criterium Routh . Hurwitz construeren een tabel (de Routh-array) uit de coëfficiënten. De stabiliteit voorwaarde is dat elk element in de eerste kolom van de array moet hetzelfde teken (meestal positief). Als er een teken verandert, het aantal teken verandert gelijk aan het aantal rechts-half-vlak-polen. De schoonheid van de test ligt in zijn zuiver algebraïsche aard: geen wortel-bepaling, geen complexe analyse. Voor integer-order systemen, dit werkt perfect omdat polynomialen hebben een eindig aantal wortels, en de linker-halfvlak-conditie is ondubbelzinnig.

Het criterium geldt ook voor systemen met vertraging (via de Padé benadering) en discrete-tijdsystemen (met behulp van de bilineaire transformatie), maar al deze uitbreidingen zijn afhankelijk van het polynomiale kader. Fractional-order systemen breken dat kader, wat een nieuw perspectief vereist.

Fractational-Order Systems: Een kort overzicht

Een fractioneel-orde systeem is een systeem waarvan de dynamiek wordt beschreven door een fractionele differentiaalvergelijking van de vorm

waarbij de differentiatieorders αi en βi niet noodzakelijk gehele getallen zijn. De Laplace-transformatie van een dergelijke vergelijking (onder nul initiële voorwaarden) levert een overdrachtsfunctie op die termen bevat als s[]γ met γ-fractionele. Bijvoorbeeld een eenvoudige fractionele-orde-integrator heeft de overdrachtsfunctie 1/[]s[[α[ met 0 < α < 2. Deze systemen vertonen eigenschappen die niet in in in gehele ordesystemen zijn gezien: oneindig-dimensionale toestand-ruimte-representatie, niet-exponentieel ontspanning, en lange-afstandsgeheugen.

Fractionele calculus heeft toepassingen gevonden in de controle van flexibele structuren, modellering van biologische weefsels, batterij impedantie karakterisatie, en zelfs financiële tijd series. De noodzaak om stabiliteit voor dergelijke systemen te controleren is acuut, maar de klassieke Routh . Hurwitz tabel kan niet worden gebouwd omdat de karakteristieke vergelijking is niet een polynoom.

Waarom de klassieke Routh . Hurwitz test mislukt voor Fractational Systems

Het obstakel is fundamenteel: een fractionele-orde karakteristiek vergelijking, zoals

De krachten van s zijn geen integer veelvouden van elkaar. De Routh-array vereist integer vermogen zodat een systematische eliminatie kan worden uitgevoerd. Bovendien leiden fractionele krachten tot takafscheidingen in het complexe vlak; het concept van

Omdat het criterium van Routh .Hurwitz in wezen een algebraïsche test is voor de Hurwitz-eigendom van een polynomiale, kan het niet direct worden toegepast op een niet-polynomiale uitdrukking. In plaats daarvan zijn er wijzigingen nodig die de geest van de test behouden, die stabiliteit bepalen door een getransformeerde algebraïsche conditie te onderzoeken en daarbij rekening te houden met de fractionele aard.

Het criterium van de stabiliteit van Matignon: een directe algebraïsche vervanging

De bekendste uitbreiding is het Matignon stabiliteitscriterium (ook wel de Matignon stelling genoemd). Het biedt een noodzakelijke en voldoende voorwaarde voor een grote klasse van fractionele-ordesystemen met overeenkomstige orden. Een fractioneel systeem zou overeenkomstige orden hebben als alle fractionele krachten in de karakteristieke vergelijking gehele veelvouden van een basisorde α zijn (bv. α, 2α, 3α, ...). Veel praktische systemen, vooral die welke zijn afgeleid van fractionele PID controllers of van fractionele viscoelastische modellen, zijn evenredig.

Matignon-criterium: Voor een fractioneel-ordesysteem met karakteristieke vergelijking

waarbij α een positief reëel getal is, is het systeem stabiel indien en alleen indien

Met andere woorden, de wortels moeten zich bevinden in een gebied van het complexe vlak buiten een sector van hoek απ/2 gecentreerd op de negatieve reële as. Deze voorwaarde reduceert tot de klassieke toestand van de linker-halfvlakte wanneer α = 1 (sindsdien απ/2 = π/2), maar voor α < 1 is het stabiliteitsgebied een wig.

Het Matignon-criterium is krachtig omdat het het stabiliteitsprobleem terug transformeert in een wortellocatieprobleem voor de getransformeerde polynoom P in de variabele w = sα[[]]. Men kan de Routh

Dit leidt tot een praktische procedure:

  1. Identificeer de basisfractieorde α zodat de karakteristieke vergelijking een polynoom wordt in sα.
  2. Define w = s[α en herschrijf de karakteristieke vergelijking als P(w) = 0, een conventionele polynoom.
  3. Toepassing van de klassieke Routh
  4. Controleer de argumentvoorwaarde: Voor elke wortel wi, berekent u de overeenkomstige [[i[[[FLT:]]] = wi[[1/α[[] en controleert u of

Het Matignon-criterium wordt vaak gebruikt als het primaire instrument voor stabiliteitsanalyse van evenredige fractionele-ordesystemen. Voor niet-commensuraatsystemen (waar de fractionele orden geen veelvouden van een gemeenschappelijke basis zijn), moeten alternatieve methoden worden gebruikt, zoals de frequentie-domeinbenadering of numerieke inversie.

Alternatieve analytische benaderingen

De Mittag •Leffler-functie en de vertegenwoordiging van de staat en de ruimte

Fractional-order systemen kunnen worden weergegeven in een pseudo-state-ruimte-vorm met behulp van de fractionele afgeleide van orde α. De stabiliteitstoestand wordt een beperking van de eigenwaarden van het systeemmatrix. In het bijzonder, voor een lineaire fractionele-orde systeem beschreven door

het systeem stabiel is als en alleen als

waarbij λ(A) de eigenwaarden van de matrix zijn A. Dit is het directe multi-variabele analoog van het Matignon-criterium. Voor overeenkomstige systemen reduceert deze voorwaarde tot hetzelfde wigvormige stabiliteitsgebied. Ingenieurs kunnen de eigenwaarden berekenen en hun argumenten controleren met behulp van standaardnumerieke lineaire algebra.

Transformatie naar een Integer-Order Systeem via variabele substitutie

Een ander analytisch pad is om een nieuwe variabele z = s[α in te voeren en vervolgens de klassieke Routh

Frequentie-domeintechnieken

De stabiliteit kan ook worden beoordeeld aan de hand van de Nyquist- of Bodecriteria aangepast aan de fractionele-orde-overdrachtfuncties. De open-loop-overdrachtsfunctie G(s)[ van een fractioneel-ordesysteem is een niet-rationele functie, maar de frequentierespons kan numeriek worden geëvalueerd. Het Nyquist-stabiliteitscriterium houdt nog steeds in: het aantal encircllementen van het −1 punt bepaalt de stabiliteit van gesloten-loop. De directe toepassing van de Routh .Hurwitz-methode wordt echter vaak aanbevolen omdat het plotten voorkomt en een algebraïsche conditie biedt.

Numeriek en Computatiegereedschap

Gezien de complexiteit van fractionele-ordersystemen is de techniek sterk afhankelijk van softwaretools. De FOMCON-toolbox (Fractional-Order Modeling and Control) voor MATLAB is een populaire keuze. Het biedt functies om stabiliteitscontroles te berekenen met behulp van het Matignon-criterium, om benaderingen (bv. recursieve benadering van Oustaloup) van fractionele operators te construeren en tijd-domeinsimulaties uit te voeren.

Andere instrumenten zijn:

  • CRONE toolbox (Commande Robuste d
  • Python-bibliotheken: fractionele (voor fractionele calculus) en controle (voor algemene analyse van het controlesysteem) kunnen samen worden gebruikt om het Matignon-criterium uit te voeren.
  • Symbolische berekening in Mathematica of Maple kan algebraïsche manipulaties van fractionele krachten en wortel-vinding in het complexe vlak verwerken.

Voor een snelle controle gebruiken veel ingenieurs het volgende numerieke algoritme:

  1. Definieer de karakteristieke vergelijking als een symbolische uitdrukking in s met fractionele krachten.
  2. Converteer naar een polynoom in w = sα als het systeem evenredig is.
  3. Bereken de wortels van de polynoom in w met behulp van een numerieke oplosser.
  4. Voor elke wortel, bereken de overeenkomstige s en zijn argument.
  5. Controleer of

Dit algoritme kan in enkele regels code worden geïmplementeerd en is robuust voor de meest evenredige systemen. Voor niet-commensurate systemen moet men vertrouwen op frequentie-domein methoden of op raster gebaseerde numerieke evaluatie van de karakteristieke vergelijking in het complexe vlak.

Praktische stapsgewijze procedure voor ingenieurs

Hieronder volgt een pragmatische workflow die de klassieke Routh . Hurwitz geest combineert met de nodige fractionele wijzigingen.

  1. Behoud de karakteristieke vergelijking van de overdrachtsfunctie van het systeem’s. Voor een fractioneel-ordesysteem zal dit termen omvatten als sγ met γ non-integer.
  2. Controleer op commensurability.[ Bepaal of alle fractionele orden gehele getallen van een gemeenschappelijke basis α zijn. Bijvoorbeeld, orden 0,5, 1.5, 2.0 zijn veelvouden van 0,5; orden 0.7 en 1.2 zijn geen veelvouden van één basis (tenzij de basis 0,1 is, maar dat is zelden rationeel). Als het systeem evenredig is, ga dan door naar de polynomiale transformatie. Zo niet, ga dan naar stap 6 (frequentie-domeinmethode).
  3. Definieer de polynome in w = s[α] Schrijf de karakteristieke vergelijking als P(w) = 0[]. Bijvoorbeeld, []s[.5[] + 3 s[1.2[ + 2 is niet evenredig (2.5 en 1.2 zijn geen veelvouden van een gemeenschappelijke basis), dus deze stap zou mislukken. Een beter voorbeeld: 1.5] + 3 s0.5] + 2 + 0]] met α ]w[FLT
  4. Voeg de klassieke Routh
  5. Verifiëren van de argumentvoorwaarde. Bereken de wortels van P(w).Voor elke wortel wi[][, berekening [[FLT:]]]i[] = w]i[1/α. De complexe kracht wordt multigewaardeerd; gebruik de belangrijkste waarde. Controle [arg(si]] > απ/2. Als dit voor alle wortels geldt, is de fractie gelijk. Als de wortel de basiswaarde in strijd is met de voorwaarde, is het systeem instabiel.
  6. Voor niet-commensureersystemen, gebruik een frequentie-domeinbenadering.[ Evalueer de karakteristieke vergelijking langs de grens van de stabiliteitszone, die voor fractionele systemen een sector is.Het Nyquist-complot van de open-loop-overdrachtsfunctie kan de stabiliteit van gesloten lus aangeven. Als alternatief, gebruik je numerieke root-finding in het s]-vlak door de fractionele operator te discreteren, bijvoorbeeld door gebruik te maken van de Oustaloup-avantage. Pas vervolgens de klassieke Routh . Hurwitz-test toe op het benaderende integer-ordersysteem. Deze methode geeft een bij benadering verkregen resultaat en moet worden gevalideerd met simulaties.

Veel ingenieurs vinden dat de overeenkomstige case een groot deel van de praktische fractionele-controle ontwerpen omvat, met name fractionele PID controllers (waar de afgeleide en integrale orders zijn veelvouden van een gemeenschappelijke α, vaak 0,5 of 0,1). De procedure hierboven is daarom algemeen toepasbaar.

Illustratieve casestudy: Fractional-Order PID Controller

Overweeg een fractionele-orde PID-regelaar (PIαDβ) die een eerste-ordeinstallatie bestuurt. De totale gesloten-loop-karakteristiek vergelijking kan zijn:

Dit systeem is niet commensuraat omdat de fractionele orden 2.5, 1.2, 0.1 geen gehele veelvouden van een gemeenschappelijke basis zijn. We kunnen het Matignon criterium niet direct toepassen. In plaats daarvan benaderen we de fractionele operators met een Oustaloup filter van orde, bijvoorbeeld 5 over een frequentiebereik [10, 100] rad/s. Dit levert een high-order integer-order overdracht functie. De benadering karakteristieke polynomial zal van orde 5×5 = 25 voor elke fractionele term, maar we kunnen de orde verminderen met behulp van evenwichtige truncatie. Dan kan de klassieke Routh . Hurwitz test worden toegepast op de gereduceerde polynomial. Een stabiliteitscontrole via simulatie bevestigt het resultaat. Deze benadering methode is de feitelijke standaard in de industrie wanneer analytische tools niet beschikbaar zijn.

Als de controller daarentegen met passende orders is ontworpen (bijvoorbeeld α = 0,5, dus orden zijn 0,5, 1,0, 1.5, 2.0), wordt het systeem evenredig. Dan kunnen we het Matignon criterium direct gebruiken, zodat we een exacte noodzakelijke en voldoende conditie krijgen (modulo de branche-cut interpretatie). Dit illustreert het voordeel van het ontwerpen van fractionele controllers met passende orders indien mogelijk.

Beperkingen en doorlopend onderzoek

De uitbreiding van de analyse van Routh .Hurwitz tot fractionele-ordersystemen is geen enkele formule, maar een verzameling technieken, elk met zijn eigen domein van toepasbaarheid. Het Matignon criterium is exact voor evenredige systemen. Voor niet-commensurate systemen kan de benadering van een integer-order foutmeldingen veroorzaken, vooral als het systeem licht bevochtigd hogefrequentiemodi heeft. Onderzoekers onderzoeken directe algebraïsche criteria voor niet-commensuraatsystemen met behulp van het concept van fractionele-order Routh arrays (bijvoorbeeld, met behulp van de fractionele-order variabele transformatie []s = (z)]1/α[ en het toepassen van een aangepaste array). Een andere veelbelovende richting is het gebruik van lineaire matrix-ongelijkheiden (LMI's) aangepast aan fractionele-ordersystemen.

Bovendien is de stabiliteitsregio voor fractionele systemen met irrationele orden geen eenvoudige wig; het kan een ingewikkelder gebied worden in het Riemann-oppervlak. Praktische ingenieurs vertrouwen meestal op numerieke evaluatie van de karakteristieke vergelijkingsomvang en fase langs de stabiliteitsgrens, die computationeel intensief is maar betrouwbare resultaten oplevert voor specifieke parameterwaarden.

Conclusie

De uitbreiding van de Routh . Hurwitz-analyse tot fractionele-ordesystemen vereist dat het klassieke polynomiale kader wordt overschreden. Het Matignon-stabiliteitscriterium voor evenredige systemen brengt het probleem opnieuw om in een wortel-argument-voorwaarde die kan worden gecontroleerd met behulp van de klassieke Routh-array en numerieke wortel-bepaling. Voor niet-commensuraatsystemen blijven frequentie-domeinmethoden en numerieke benaderingen de primaire instrumenten. Door zowel de theoretische als praktische numerieke workflow te begrijpen, kunnen ingenieurs met vertrouwen stabiliteitsanalyse uitvoeren voor een breed scala aan fractionele-ordesystemen die de kloof tussen klassieke controle en nieuwe fractionele-ordetoepassingen overbruggen.

Voor verdere lezing kunnen lezers de Wikipedia-invoer op het criterium van Routh