In de analyse van systemen met tijdsvertraging komt vaak de vraag naar voren hoe de discretisatie van de systeemoperatoren het beste kan worden uitgevoerd zonder de nauwkeurigheid te verliezen, maar de rekenlast aanzienlijk te verlichten. Traditionele benaderingen van spectral discretisatie vereisen dat alle systeemvariabelen over het vertragingstijdinterval worden gediscretiseerd, wat leidt tot matrices van enorme omvang. Deze benadering is computationally intensief en vereist aanzienlijke tijd en geheugen, vooral bij grote systemen. Het gebruik van de PSD-methode biedt echter een efficiëntere aanpak die de matrixdimensies aanzienlijk verlaagt zonder in te boeten op nauwkeurigheid.

De kern van de PSD-methode ligt in het idee dat alleen de vertraagde systeemvariabelen gediscretiseerd hoeven te worden, terwijl de tijdsafhankelijke variabelen die geen vertraging vertonen, niet betrokken hoeven te worden bij de discretisatie. Dit heeft als voordeel dat de dimensionale vergroting van de matrices die de systeemoperatoren benaderen aanzienlijk wordt verminderd. Door deze aanpak wordt de hoeveelheid verwerkte data drastisch verminderd, wat niet alleen het geheugenverbruik optimaliseert, maar ook de rekentijd van de simulaties verkort.

Bij de toepassing van PSD wordt het systeem opgedeeld in twee delen: de vertraging-vrije variabelen, die niet betrokken zijn bij de dynamica van het systeem, en de vertraagde variabelen, die een directe invloed hebben op het systeemgedrag. Deze verdeling maakt het mogelijk om de discrete matrices, die de systeemoperatoren representeren, met lagere orde te benaderen dan in traditionele methoden. De voordelen van deze benadering zijn vooral duidelijk wanneer het systeem veel vertraging-vrije variabelen bevat die niet bijdragen aan de dynamiek van het systeem.

Verder is de applicatie van numerieke methoden zoals de pseudospectrale methode (PS) en impliciete Runge-Kutta-methoden voor de discretisatie van de operatoren .A en .T(h) een essentieel onderdeel van de PSD-benadering. De analytische en voortdurend reguliere aard van de karakteristieke vergelijkingen van deze operatoren maakt de pseudospectrale methode bijzonder krachtig bij het benaderen van deze operatoren met een nauwkeurigheid die bekend staat als "spectrale nauwkeurigheid." Dit betekent dat de fout tussen de werkelijke en benaderde eigenwaarden exponentieel snel afneemt naarmate het aantal discretisatiepunten toeneemt.

Het idee van PSD kan worden uitgebreid naar eigenwaarde-analysemethoden, zoals het Partial Infinitesimal Generator Discretization (PIGD), waarbij de oorspronkelijke Cauchy-problemen worden herschreven en geëvalueerd op een set van discrete punten over het tijdsvertraginginterval. Dit maakt het mogelijk om een generaliseerde eigenwaardeprobleem te formuleren, wat uiteindelijk leidt tot een standaard eigenwaardeprobleem voor de gediscretiseerde matrix die de systeemoperator .A benadert. De resultaten van deze methode bieden een nauwkeurige benadering van de werkelijke eigenwaarden van het tijdvertraging-systeem, zonder de zware belasting van een volledige discretisatie van alle systeemvariabelen.

Daarnaast wordt de methode van Partial Solution Operator Discretization (PSOD) geïntroduceerd om de systeemoperator .T(h) te discretiseren. Deze benadering vereist dat alleen de vertragingstermijnen op discrete tijdstappen worden geëvalueerd, terwijl de berekening van de niet-vertraagde variabelen wordt verminderd. Dit vermindert de noodzaak om de gehele tijdsvertragingsegmenten over meerdere discrete punten te evalueren, zoals het geval zou zijn bij een traditionele spectrale discretisatie.

Bij het toepassen van de PSD-methode wordt het systeem geherstructureerd om de matrixcomponenten te scheiden volgens de vertraging-vrije en vertraagde variabelen. Dit maakt de matrixoperaties eenvoudiger en vermindert de complexiteit van de berekeningen. De gewijzigde matrixstructuren maken het mogelijk om eigenwaarden sneller en efficiënter te berekenen zonder in te boeten op de nauwkeurigheid van de simulatie, wat essentieel is voor het werken met grote, complexe tijdsvertraging-systemen.

Naast de theoretische voordelen van de PSD-methode is het ook belangrijk om te begrijpen hoe de keuze van discretisatiepunten en de resolutie van het tijdsinterval invloed kan hebben op de nauwkeurigheid en stabiliteit van de resultaten. Bij het toepassen van deze benadering is het cruciaal om een balans te vinden tussen de rekenkundige efficiëntie en de gewenste nauwkeurigheid. Het is vaak noodzakelijk om te experimenteren met verschillende instellingen voor de discretisatiepunten om de meest geschikte configuratie voor het specifieke systeem te vinden.

Hoe kan rotatie- en vermenigvuldigingspreconditionering de eigenwaarde-berekening voor tijdvertragingssystemen verbeteren?

De rotatie- en vermenigvuldigingspreconditionering biedt een geavanceerde manier om de berekening van de eigenwaarden van tijdvertragingssystemen te verbeteren, waarbij met name de eigenschommelingsmodi met lage demping snel en nauwkeurig kunnen worden bepaald. Deze techniek is van cruciaal belang voor systemen die dynamische stabiliteit vertonen en vereist efficiënte methoden voor de eigenschapsanalyse.

Het belangrijkste aspect van deze benadering is het gebruik van rotatie in de complexe vlakken van de eigenwaarden. Dit zorgt ervoor dat de dicht bij de eenheidsomloop gelegen eigenwaarden uitgerekt worden, wat de snelheid van de algoritmes die eigenwaarden berekenen, bevordert. Wanneer deze rotatie samen met een vermenigvuldiging van de eigenwaarden door een constante wordt toegepast, zoals in de rotatie- en vermenigvuldigingspreconditionering, worden de eigenwaarden in het z-vlak samengeperst. Dit kan echter de convergentiesnelheid van algoritmes die afhankelijk zijn van Krylov-subruimtes verlagen, omdat de eigenwaarden nu dichter bij elkaar liggen en dus moeilijker te onderscheiden zijn.

Voor de implementatie van deze preconditionering moeten twee hoofdopties in overweging worden genomen. De eerste optie houdt in dat de tijdstaplengte hh constant wordt gehouden, terwijl de tijdvertragingen τi\tau_i in de formule voor eigenwaarde-verandering door een vermenigvuldigingsfactor α\alpha worden aangepast. De tweede optie vergroot de tijdstaplengte hh met een factor α\alpha, terwijl de tijdvertragingen τi\tau_i hetzelfde blijven.

Het kiezen van de juiste waarden voor de parameters α\alpha en θ\theta is essentieel voor het succes van de rotatie- en vermenigvuldigingspreconditionering. De keuze van α\alpha beïnvloedt direct de snelheid van de convergentie van de eigenwaarde-algoritmes. Waarden tussen 2 en 3 voor α\alpha blijken doorgaans voldoende te zijn om de convergentie te verbeteren. Toch mag deze keuze niet in strijd zijn met de beperking αhτmax\alpha h \leq \tau_{\text{max}}, waar τmax\tau_{\text{max}} de maximale tijdvertraging in het systeem is. Wat betreft de parameter θ\theta, deze moet afhankelijk van de toepassing worden gekozen. Bij systemen waar veel elektromechanische oscillatiemodi met lage demping nodig zijn, wordt een grotere waarde voor θ\theta aanbevolen, bijvoorbeeld θ=5.74\theta = 5.74^\circ, wat correspondeert met een dempingspercentage van 10%. Voor systemen waar slechts enkele kritieke eigenwaarden snel moeten worden berekend, is een kleinere waarde voor θ\theta beter, bijvoorbeeld θ=1.72\theta = 1.72^\circ voor een dempingspercentage van 3%.

De keuze van de juiste rotatiehoek en vermenigvuldigingsfactor is dus niet alleen afhankelijk van de wens om snel resultaten te krijgen, maar moet ook rekening houden met de stabiliteit en de nauwkeurigheid van de berekeningen, die cruciaal zijn voor het betrouwbaar voorspellen van de systeemrespons en stabiliteit op lange termijn.

Hoewel de rotatie- en vermenigvuldigingspreconditionering aanzienlijke voordelen biedt voor de eigenwaarde-analyse van tijdvertragingssystemen, heeft het ook zijn beperkingen. In tegenstelling tot de PIGD-methoden met de verschuivings-inversie transformatie, die eigenwaarden dicht bij vooraf gedefinieerde verschuivingen snel kan berekenen, kunnen de PSOD-methoden met rotatie- en vermenigvuldigingspreconditionering langzamer zijn en vereisen ze meer iteraties, wat de rekenintensiteit van het proces vergroot. De eigenwaarden worden namelijk geconcentreerd nabij de eenheidsomloop, wat een uitdaging vormt voor de snelheid van convergentie.

Daarom is het essentieel om een balans te vinden tussen de rekenkracht die nodig is voor het toepassen van deze preconditionering en de vereiste nauwkeurigheid van de eigenwaardeanalyse. De gekozen benadering moet altijd worden afgestemd op de specifieke eisen van de toepassing en de systemen die geanalyseerd worden.

In de praktijk zal het toepassen van de rotatie- en vermenigvuldigingspreconditionering in combinatie met zorgvuldig gekozen parameters α\alpha en θ\theta niet alleen de nauwkeurigheid van de berekeningen verbeteren, maar ook de betrouwbaarheid van de stabiliteitsanalyse verhogen. Hierdoor kunnen grote en complexe tijdvertragingssystemen met vertrouwen worden geëvalueerd, wat essentieel is voor de ontwikkeling van geavanceerde controle- en simulatiemethoden voor dynamische systemen.

Hoe de PS-Methode bijdraagt aan de Discretisatie van Tijdvertragingen in Grootschalige Systemen

In het kader van grootschalige tijdsvertraging systemen wordt de schatting van toestandsvariabelen vaak uitgevoerd door integratie van de tijdsdiscretisatie. Dit proces kan als volgt worden geformuleerd: de schattingen van de toestanden van het systeem, zoals x^1,1,k\hat{x}_1,1,k en y^1,1,k\hat{y}_1,1,k, worden bepaald over het tijdsinterval [0,h][0, h]. Het proces wordt vaak geformuleerd door gebruik te maken van integralen en de parameters ϕx,1,0\phi_x,1,0 en ϕy,1,0\phi_y,1,0 als initiële toestanden, waaruit verdere schattingen kunnen worden berekend.

Het proces zelf kan wiskundig als volgt worden uitgedrukt:

x^1,1,k=x^(h+θM,1,k),k=0,1,,M1\hat{x}_{1,1,k} = \hat{x}(h + \theta M,1,k), \quad k = 0, 1, \ldots, M-1
y^1,1,k=y^(h+θM,1,k),k=0,1,,M1\hat{y}_{1,1,k} = \hat{y}(h + \theta M,1,k), \quad k = 0, 1, \ldots, M-1

Hierbij wordt de schatting van de toestanden x^1,1,k\hat{x}_1,1,k en y^1,1,k\hat{y}_1,1,k benaderd door een integraal over een bepaald segment, waarbij de toestand ϕx,1,0\phi_x,1,0 en ϕy,1,0\phi_y,1,0 dienen als uitgangspunt. Dit resulteert in de volgende uitdrukkingen voor de toestanden:

j=1N0h+θM,1,kzjds\sum_{j=1}^{N} \int_0^{h+\theta M,1,k} z_j ds

en

j=1N0h+θM,1,kwjds.\sum_{j=1}^{N} \int_0^{h+\theta M,1,k} w_j ds.

De discretisatie van het tijdsegment kan verder geformaliseerd worden door te werken met de zogenaamde PS-partiële discretisatie, die de relatie tussen x^1,1,k\hat{x}_1,1,k en y^1,1,k\hat{y}_1,1,k tijdens het interval [h,0][-h, 0] beschrijft. Dit zorgt ervoor dat de schattingen x^1,1,k\hat{x}_1,1,k en y^1,1,k\hat{y}_1,1,k onafhankelijk blijven van de vertragingsovergang, waarbij de berekeningen van de toestanden niet worden beïnvloed door eerdere vertragingen of waarden die buiten het tijdsinterval [τmax,0)[-\tau_{\text{max}}, 0) vallen.

Het proces kan verder worden geoptimaliseerd door gebruik te maken van geavanceerde discretisatie technieken, zoals de Lagrange-interpolatiemethode, die wordt toegepast om de toestanden van het systeem over subintervallen te schatten. Dit wordt met name belangrijk wanneer we werken met kleinere subintervallen die de nauwkeurigheid van de schattingen kunnen verbeteren.

Wanneer we deze techniek toepassen, kan de PS-partiële discretisatie van het verschuivingssegment worden uitgedrukt als volgt:

x^1,i,k=x^(h+θM,i,k),i=2,3,,Q,k=0,1,,M.\hat{x}_{1,i,k} = \hat{x}(h + \theta M,i,k), \quad i = 2, 3, \ldots, Q, \quad k = 0, 1, \ldots, M.

Voor i=Qi = Q, waar het subinterval korter wordt dan de voorgaande, is er een noodzaak om de schattingen x^1,Q,k\hat{x}_{1,Q,k} en y^1,Q,k\hat{y}_{1,Q,k} te berekenen door middel van de Lagrange-interpolatie, waarbij de waarden van de toestanden worden berekend door een gewogen som van de waarden op voorgaande tijdstappen.

De PS-methode biedt vervolgens een uitdrukking voor de partiële discretisatie van het verschuivingssegment T(h)T(h) van het systeem:

T^M=[0IMn2U],\hat{T}_M = \left[\begin{array}{ccc} 0 & I_{Mn2} & U \\
\end{array}\right],

waarbij UU de gewogen matrix is die de interpolatiewaarden van de toestanden vertegenwoordigt.

Een essentieel onderdeel van het proces van discretisatie is de berekening van de eigenwaarden van grootschalige tijdvertraging systemen. Dit wordt gedaan door de zogenaamde "onvolledige coördinaterotatie" toe te passen, waarmee eigenwaarden met dempingsverhoudingen onder een bepaalde drempel ζ\zeta kunnen worden gevangen. Deze rotatie wordt uitgevoerd door de toestanden te transformeren met een rotatiehoek van θ=arcsinζ\theta = \arcsin \zeta in tegenwijzerzin. Deze techniek helpt om de efficiëntie van het algoritme te verbeteren door de dispersie van de eigenwaarden te verhogen, wat resulteert in een versnelde convergentie van het iteratieve oplossingsproces.

Een belangrijk aspect van de efficiëntie van de berekeningen is de toepassing van preconditioneringstechnieken. Dit kan worden gedaan door de zogenaamde rotatie- en vermenigvuldigingspreconditionering toe te passen, wat de spreiding van de eigenwaarden verder verbetert. Dit leidt tot een snellere convergentie bij het oplossen van grootschalige systemen, waarbij de berekeningen van de toestanden in de tijdsegmenten sneller kunnen worden uitgevoerd.

De rotatie- en vermenigvuldigingspreconditionering wordt als volgt geïmplementeerd: de vertraging τi\tau_i wordt geherstructureerd naar 1/α1/\alpha van de oorspronkelijke waarde, waarbij de vertragingintervallen opnieuw worden verdeeld in kleinere subintervallen. Dit zorgt ervoor dat de iteraties sneller convergeren door de structuur van de systeemtoestanden te optimaliseren.

Wat belangrijk is voor de lezer is het besef dat het gebruik van de PS-partiële discretisatie in tijdsvertraging systemen een krachtige methode is om de schattingen van de toestanden over verschillende tijdsintervallen te verbeteren. Het biedt een fundamentele benadering voor het oplossen van complexe systemen met vertragingen, waarbij traditionele methoden niet effectief zijn. Bovendien speelt de toepassing van preconditionering een cruciale rol in de snelheid en efficiëntie van de oplossing van dergelijke grootschalige systemen, vooral wanneer de dimensies van het systeem aanzienlijk groot zijn. De juiste implementatie van deze technieken kan de rekentijd aanzienlijk verkorten, wat van groot belang is voor toepassingen in de engineering en wetenschappen.

Wat is het effect van preconditionering in de discretisatie van tijdvertraging systemen?

In de context van tijdvertraging systemen en hun benadering via discretisatie, is het cruciaal om te begrijpen hoe verschillende implementaties van preconditionering de matrixstructuur beïnvloeden en uiteindelijk de oplossing van het systeem kunnen verbeteren. De volgende uitleg betreft de specifieke effecten van de PSOD-PS methode, waarbij de matrixrepresentaties van het systeem worden geoptimaliseerd door gebruik te maken van twee verschillende implementaties van rotatie- en vermenigvuldigingspreconditionering.

Bij de eerste implementatie van de preconditionering worden de parameters zoals de vertragingstijd τi (i = 1, 2, ..., m), de frequentie h, en de systeemmatrixcomponenten van de matrices U, Ū en Ũ in het systeem behouden. Het doel is de discretisatie van de tijdvertragingen in subintervallen, die met de parameter Q′ en de stapgrootte αh worden herverdeeld. Deze transformatie leidt tot een herschikking van de matrixmaten en de toevoeging van nieuwe termen die verband houden met de vertragingen. De matrices worden aangepast, waarbij bepaalde submatrices constant blijven, terwijl andere herberekend worden om de impact van de vermenigvuldigingspreconditionering te integreren.

In de tweede implementatie wordt dezelfde set van matrices geherstructureerd, maar met de nadruk op het verhogen van de invloed van de vertragingstijden door ze met een factor α te vermenigvuldigen. Het belangrijkste verschil is dat, hoewel de parameters in de eerste implementatie constant blijven, ze in de tweede implementatie worden geschaald, wat leidt tot een evenredige verandering in de elementen van de matrices. Dit wordt weerspiegeld in de evolutie van de termen zoals tQ,k,j en θM,i,k, die veranderen in reactie op de schaling van de tijdvertragingen.

Beide benaderingen, ondanks hun verschillen in de manier waarop vertragingen worden behandeld, leiden tot een identieke eindresultaat wat betreft de structuur van de systeemmatrix. De effectiviteit van deze methoden wordt vergeleken door de laatste rijen van de matrices T̂M,N en TM,N te analyseren, waarbij blijkt dat ze overeenkomen, wat wijst op een consistente aanpak van de preconditionering.

Bij het werken met grootschalige tijdvertraging systemen is het essentieel te begrijpen hoe de keuze van de implementatie invloed heeft op de nauwkeurigheid en stabiliteit van de berekeningen. Terwijl de eerste implementatie meer gericht is op het behouden van de oorspronkelijke structuur van de matrix, is de tweede implementatie meer geschikt voor situaties waarin er meer controle nodig is over de vertragingstijden zelf. Deze afstemming van de preconditionering kan de convergentie van de oplossingsmethoden versnellen, vooral bij complexe systemen met lange tijdvertragingen.

De bovenstaande implementaties zijn gebaseerd op de veronderstelling dat de matrixstructuur voldoende goed is geconditioneerd voor de gekozen discretisatie. Dit is van belang voor de lezer, omdat een slecht geconditioneerde matrix kan leiden tot numerieke instabiliteiten, ongeacht de toegepaste preconditionering. Het is daarom noodzakelijk om bij de keuze van de implementatie niet alleen de theoretische voordelen van de preconditionering in overweging te nemen, maar ook de praktische effecten op de numerieke stabiliteit van het systeem.

Naast het begrijpen van de matrixtransformaties die plaatsvinden bij het gebruik van preconditionering, is het ook van belang om in gedachten te houden hoe deze methoden zich verhouden tot de algehele opzet van het systeem. Wanneer de tijdvertragingen in een systeem worden discretiseerd, kunnen de gekozen instellingen voor de vertragingstijden en de schaalfactoren (zoals α) direct van invloed zijn op de snelheid van convergentie en de nauwkeurigheid van de uiteindelijke oplossing. Dit betekent dat de keuze van de preconditionering niet los kan worden gezien van de specifieke eisen van het probleem dat men probeert op te lossen, zoals de grootte van het systeem en de gewenste precisie van de oplossing.

Hoe vertragingen in feedback- en controlekanalen de stabiliteit van krachtensystemen beïnvloeden

In moderne krachtensystemen, zoals de UHV North China-Central China interconnectie, spelen tijdvertragingen in feedback- en controlekanalen een cruciale rol in de dynamische stabiliteit van het netwerk. Het begrijpen van de effecten van deze vertragingen op de kleine signaalstabiliteit is essentieel voor het ontwerp van effectieve controlesystemen. Het is belangrijk te realiseren dat vertragingen niet alleen ontstaan door fysieke afstanden of technologie, maar ook door de complexiteit van het netwerk en de noodzaak voor gedistribueerde controlesystemen.

In de context van een UHV-systeem, bijvoorbeeld, zijn de vertragingen in de feedbacksignalen van twee controllers in de Daihai- en Donghai-krachtcentrales van het Mengxi- en Shandong-netwerk respectievelijk 120 ms en 100 ms. Deze vertragingen beïnvloeden de efficiëntie van het energiebeheer en kunnen leiden tot aanzienlijke oscillaties of zelfs destabilisatie van het systeem, vooral wanneer de netfrequentie van inter-area oscillatiemodi niet goed wordt beheerd. De noodzaak om deze vertragingen correct te modelleren en te begrijpen is dus van groot belang voor het handhaven van de stabiliteit.

De vertragingen van feedback- en controlekanalen moeten altijd worden meegenomen in de simulaties en analyses van een systeem. De tijdvertragingen in de signalen tussen de units, zoals tussen de Daihai- en Gaoer-krachtcentrales, kunnen leiden tot oscillaties in de inter-area frequenties. Bijvoorbeeld, de rotorsnelheid van de turbines in verschillende centrales kan onderhevig zijn aan vertragingen in de communicatie van feedbacksignalen, wat de frequentie van oscillaties in het netwerk kan veranderen. Het effect van vertragingen op de kleine signaalstabiliteit kan zelfs resulteren in een verlies van stabiliteit bij hogere vermogensniveau’s.

De efficiëntie van de controle en het gedrag van de feedbacksignalen wordt vaak geanalyseerd met behulp van eigenwaarde-analysemethoden. Eigenwaarden zijn immers directe indicatoren van de stabiliteit van het systeem, waarbij negatieve realen duiden op een stabiel systeem en positieve realen wijzen op een instabiele toestand. Het nauwkeurig berekenen van deze eigenwaarden, rekening houdend met de vertragingen, kan alleen worden uitgevoerd met geavanceerde technieken zoals de PIGD-PS-methode, die speciaal ontworpen is om grote tijdvertragingen in systemen te verwerken.

De PIGD-PS-methode is bewezen zeer effectief te zijn in het verbeteren van de nauwkeurigheid van de stabiliteitsanalyse van systemen met meerdere tijdvertragingen. Door gebruik te maken van een discretisatiematrix van het infinitesimale generatorenmodel, kan deze methode de invloed van vertragingen nauwkeuriger in kaart brengen. Dit maakt het mogelijk om eigenwaarden die verband houden met de vertragingen te berekenen met behulp van technieken zoals de QR-methode en de Newton-methode. Door iteratief de waarden te verfijnen, kunnen de exacte waarden worden vastgesteld, waarmee de stabiliteit van het systeem beter kan worden voorspeld.

Het belang van een correcte analyse van tijdvertragingen kan niet worden onderschat. Vertragingen in het netwerk kunnen de dynamiek aanzienlijk veranderen, zelfs als ze klein lijken. Zo blijkt uit onderzoeken dat wanneer de tijdvertragingen in de feedbacksignalen worden verwaarloosd, de kleine signaalstabiliteit van het systeem aanzienlijk kan verslechteren. Bij grotere vertragingen, zoals in de interconnectie van het Noord- en Centraal-China-systeem, kunnen zelfs kleine vertragingen leiden tot onvoorziene oscillaties die het netwerk destabiliseren.

Daarom moet elke systeemontwerper of operator van krachtensystemen het effect van vertragingen als een integraal onderdeel van hun stabiliteitsanalyse beschouwen. Het gebruik van geavanceerde methoden zoals PIGD-PS helpt niet alleen bij het identificeren van kritieke vertragingen, maar ook bij het verbeteren van de respons van het netwerk onder variërende omstandigheden. Het begrijpen van de exacte impact van vertragingen op de stabiliteit vereist een gedetailleerde kennis van de systeemdynamiek en de tijdsafhankelijke variabelen die de oscillaties beïnvloeden.

Bij de analyse van tijdvertragingen is het belangrijk te realiseren dat een te gedetailleerde discretisatie van het systeem niet altijd noodzakelijk is voor het behoud van de stabiliteit. De keuze van de juiste graad van discretisatie kan bepalend zijn voor de efficiëntie van de berekeningen. In veel gevallen blijkt een lagere granulariteit al voldoende te zijn om de stabiliteit van het systeem adequaat te voorspellen, zoals aangetoond door simulaties met een discretisatiegraad van N = 20.

Een ander aspect dat van belang is bij de analyse van tijdvertragingen in grote systemen is het effect van pseudo-vertraagde toestandsvariabelen die kunnen ontstaan bij het gebruik van traditionele DDE-gebaseerde methoden. Deze pseudo-vertraagde variabelen kunnen de nauwkeurigheid van de berekeningen beïnvloeden en leiden tot inefficiënties. Daarom wordt aangeraden om de nieuwste methoden zoals PIGD-PS te gebruiken, die deze problemen minimaliseren en de stabiliteitsanalyse van het systeem versnellen.