2025-11-24T10:10:18.436123

Computing the Hard Scaled Relative Graph of LTI Systems

Krebbekx, Baron-Prada, Tóth et al.
Scaled Relative Graphs (SRGs) provide a novel graphical frequency-domain method for the analysis of nonlinear systems, where Linear Time-Invariant (LTI) systems are the fundamental building block. To analyze feedback loops with unstable LTI components, the hard SRG is required, since it aptly captures the input/output behavior on the extended $L_2$ space. In this paper, we develop a systematic computational method to exactly compute the hard SRG of LTI systems, which may be unstable and contain integrators. We also study its connection to the Nyquist criterion, including the multivariable case, and demonstrate our method on several examples.
academic

Calcolo del Grafico Relativo Scalato Hard dei Sistemi LTI

Informazioni Fondamentali

  • ID Articolo: 2511.17297
  • Titolo: Computing the Hard Scaled Relative Graph of LTI Systems
  • Autori: Julius P. J. Krebbekx, Eder Baron-Prada, Roland Tóth, Amritam Das
  • Istituzioni: Eindhoven University of Technology, Austrian Institute of Technology, ETH Zurich, HUN-REN Institute for Computer Science and Control
  • Classificazione: eess.SY (Sistemi e Controllo), cs.SY, math.OC
  • Data di Sottomissione: 21 novembre 2025
  • Link Articolo: https://arxiv.org/abs/2511.17297

Riassunto

I Grafici Relativi Scalati (SRG) forniscono un metodo grafico innovativo nel dominio della frequenza per l'analisi di sistemi non lineari, dove i sistemi lineari tempo-invarianti (LTI) costituiscono i blocchi costruttivi fondamentali. Per analizzare anelli di retroazione contenenti componenti LTI instabili, è necessario utilizzare l'SRG hard, poiché è in grado di catturare adeguatamente il comportamento ingresso/uscita sullo spazio L₂ esteso. Questo articolo sviluppa un metodo di calcolo sistematizzato per calcolare con precisione l'SRG hard di sistemi LTI, che possono essere instabili e contenere integratori. L'articolo inoltre indaga i collegamenti con il criterio di Nyquist, incluso il caso multivariabile, e dimostra il metodo attraverso molteplici esempi.

Contesto di Ricerca e Motivazione

Problema di Ricerca

Questo articolo si propone di risolvere il problema centrale nell'analisi di sistemi non lineari: come calcolare sistematicamente il Grafico Relativo Scalato (SRG) di sistemi di retroazione contenenti componenti LTI instabili.

Importanza del Problema

  1. Limitazioni dei metodi classici: Il diagramma di Nyquist è la pietra miliare dell'analisi dei sistemi LTI, ma è difficile generalizzarlo sistematicamente all'analisi di sistemi non lineari
  2. Esigenze pratiche: Le applicazioni industriali richiedono frequentemente l'analisi di sistemi di retroazione contenenti componenti instabili
  3. Completezza teorica: Il framework SRG esistente manca di un metodo sistematizzato per gestire sistemi instabili

Limitazioni dei Metodi Esistenti

  1. SRG Soft: Inizialmente proposto in 3, applicabile solo a sistemi SISO stabili, con problemi per sistemi instabili
  2. SRG Esteso: Risolve il problema dei sistemi SISO instabili in 11 aggiungendo informazioni dal criterio di Nyquist, ma limitato al caso monovariabile
  3. SRG Hard: Definito in 12, teoricamente in grado di gestire sistemi instabili, ma manca di un metodo di calcolo sistematizzato

Motivazione della Ricerca

Poiché il principale vantaggio dell'SRG hard risiede nel gestire sistemi instabili negli anelli di retroazione, lo sviluppo di metodi per calcolare l'SRG hard di sistemi LTI instabili avanzerebbe significativamente l'analisi SRG di sistemi non lineari.

Contributi Fondamentali

  1. Metodo di Calcolo Sistematizzato: Propone un algoritmo non approssimato (Teorema 2) che può calcolare con precisione l'SRG hard di sistemi LTI, inclusi sistemi instabili e sistemi con integratori
  2. Prova di Equivalenza Teorica: Dimostra che per sistemi SISO, l'SRG hard è completamente equivalente all'SRG esteso (Teorema 3)
  3. Estensione MIMO: Estende il metodo al caso multivariabile (MIMO), fornendo un framework di calcolo completo per sistemi square
  4. Collegamento con il Criterio di Nyquist: Mostra come l'SRG hard funziona come metodo alternativo al criterio MIMO di Nyquist e al criterio di Nyquist Generalizzato (GNC)
  5. Algoritmo Pratico: Basato sulla rappresentazione della funzione di trasferimento, costruisce il confine attraverso intersezioni nel piano complesso e rimozione di dischi centrati sull'asse reale

Dettagli del Metodo

Definizione del Compito

Data un sistema LTI G: L²ₑᵖ → L²ₑᵍ con funzione di trasferimento G(s) ∈ Rᵍˣᵖ(s), calcolare il suo SRG hard:

SRGe(G)=⋃T>0SRG(G∣T)\text{SRG}_e(G) = \bigcup_{T>0} \text{SRG}(G|_T)

dove G|_T è l'operatore di troncamento di G sull'intervallo di tempo 0,T.

Framework Teorico Centrale

1. Rappresentazione della Convessità Iperbolica

La chiusura di qualsiasi insieme h-convesso C può essere rappresentata come:

cl C=⋂α∈R(DRα(α)∖Drα(α))=:G({rα},{Rα})\text{cl}\,C = \bigcap_{\alpha \in \mathbb{R}} \left(D_{R_\alpha}(\alpha) \setminus D_{r_\alpha}(\alpha)\right) =: \mathcal{G}(\{r_\alpha\}, \{R_\alpha\})

dove Dr(x)D_r(x) denota il disco chiuso centrato in x con raggio r.

2. Definizione del Sistema Traslato

Per α ∈ ℝ, definire il sistema traslato:

Gα=(G0(n−q)×p)−(αI0(n−p)×p)G_\alpha = \begin{pmatrix} G \\ 0_{(n-q)\times p} \end{pmatrix} - \begin{pmatrix} \alpha I \\ 0_{(n-p)\times p} \end{pmatrix}

dove n = max{p,q}.

Componenti Fondamentali dell'Algoritmo

Proposizione 2: Calcolo del Raggio di Guadagno Massimo Rα

\bar{\sigma}(G_\alpha) & \text{se } G \in RH_\infty^{q\times p} \\ \infty & \text{altrimenti} \end{cases}$$ **Idea Chiave**: - Per sistemi stabili, utilizzare la norma H∞ (valore singolare massimo) - Per sistemi instabili o non regolari, Rα = ∞ #### Proposizione 3: Calcolo del Raggio di Guadagno Minimo rα $$r_\alpha = \begin{cases} \underline{\sigma}(G_\alpha) & \text{se } G_\alpha \text{ è a fase minima} \\ 0 & \text{altrimenti} \end{cases}$$ **Idea Chiave**: - Attraverso il calcolo del sistema inverso $G_\alpha^{-1}$ e applicazione della Proposizione 2 - Solo quando Gα è a fase minima (tutti gli zeri di trasmissione nel semipiano sinistro aperto), rα è non nullo - Utilizzo della forma di Smith-McMillan: i poli di Gα sono gli zeri di $G_\alpha^{-1}$, e viceversa ### Algoritmo Principale (Teorema 2) **Per sistemi LTI square** G: L²ₑᵖ → L²ₑᵖ con funzione di trasferimento G(s) ∈ Rᵖˣᵖ(s): $$\text{SRG}_e(G) = \mathcal{G}(\{r_\alpha\}, \{R_\alpha\})$$ **Passi dell'Algoritmo**: 1. Per ogni α ∈ ℝ, costruire Gα = G - αI 2. Utilizzare la Proposizione 2 per calcolare Rα 3. Utilizzare la Proposizione 3 per calcolare rα 4. Calcolare l'intersezione di tutte le regioni anulari ### Punti di Innovazione Tecnica #### 1. Gestione dello Spazio L₂ Esteso - I metodi tradizionali operano nello spazio L₂, incapaci di gestire segnali divergenti - L'SRG hard è definito nello spazio esteso L²ₑ, che include segnali periodici e divergenti - La gestione dell'instabilità avviene considerando segnali troncati per tutti i T > 0 #### 2. Rilevamento della Fase Minima - Innovazione chiave: ridurre il calcolo di rα al rilevamento se Gα è a fase minima - Fase minima ⟺ tutti gli zeri di trasmissione nel semipiano sinistro aperto - Per sistemi non a fase minima, rα = 0, riflettendo i potenziali problemi nella retroazione #### 3. Unificazione con il Criterio di Nyquist Per sistemi SISO, è provato che: $$\text{SRG}_e(G) = \text{SRG}'(G)$$ dove SRG'(G) è l'SRG esteso, che include informazioni dall'SRG soft e dal criterio di Nyquist. #### 4. Gestione di Casi Speciali - **Integratori**: Portano automaticamente a Rα = ∞ - **Sistemi non regolari**: Rilevati attraverso la forma di Smith-McMillan - **Ritardi temporali**: Portano a rα = 0, Rα rimane inalterato ## Configurazione Sperimentale ### Sistemi di Esempio #### Sistemi SISO 1. **G₁(s) = 1/(s² + s + 1)**: Sistema del secondo ordine stabile 2. **G₂(s) = 1/(s(s+1))**: Sistema con integratore 3. **G₃(s) = 3/((s-2)(s/10+1))**: Sistema instabile (polo nel semipiano destro) #### Sistemi MIMO 1. **G₄(s)**: Sistema 2×2 stabile $$G_4(s) = \begin{pmatrix} \frac{1}{s+1} & \frac{1}{s+2} \\ \frac{1}{s+4} & \frac{1}{s+3} \end{pmatrix}$$ 2. **G₅(s)**: Sistema 2×2 instabile $$G_5(s) = \begin{pmatrix} \frac{s+7}{s-1} & \frac{s-5}{(s+2)^2} \\ \frac{1}{(s+4)^3} & \frac{s}{(s+3)^2} \end{pmatrix}$$ 3. **G₆(s)**: Sistema 3×3 con integratore ### Metodo di Visualizzazione - Tracciare il confine dell'SRG hard nel piano complesso (area grigia) - Sovrapporre il confine dell'SRG soft (area grigio scuro) - Contrassegnare punti critici come il punto -1 (per l'analisi di stabilità) ## Risultati Sperimentali ### Risultati Principali #### Risultati dei Sistemi SISO 1. **G₁**: L'SRG hard coincide con l'SRG soft (poiché il sistema è stabile), formando una regione limitata 2. **G₂**: A causa dell'integratore, l'SRG hard si estende all'infinito, includendo il semiasse reale [0,∞) 3. **G₃**: L'SRG hard del sistema instabile include regioni aggiuntive, riflettendo le informazioni di avvolgimento di Nyquist **Verifica dell'Equivalenza dell'SRG Esteso**: - Per G₂, $\mathcal{N}_{G_2} \cap \mathbb{R} = [0,\infty)$ - Questo corrisponde completamente all'SRG hard nella figura 2b - Verifica la correttezza del Teorema 3 #### Risultati dei Sistemi MIMO **G₄ (Sistema 2×2 stabile)**: - Il confine dell'SRG hard e dell'SRG soft sono molto vicini - Entrambi sono separati dal punto -1, indicando stabilità della retroazione **G₅ (Sistema 2×2 instabile)**: - L'SRG hard mostra **due buchi** (figura 2e) - Questi buchi corrispondono a regioni di configurazione di anello chiuso stabile - Il buco più piccolo può essere utilizzato per il progetto del guadagno **G₆ (Sistema 3×3 con integratore)**: - L'SRG hard è parzialmente illimitato (a causa dell'integratore) - Il confine dell'SRG soft tende a ±j∞ nella direzione Im ### Applicazione dell'Analisi di Stabilità #### Esempio di Progetto del Guadagno (Basato su G₅) L'articolo presenta tre strategie di selezione del guadagno: 1. **Scalatura Pura**: k₁ = 0, k₂ = 1/5 - Intervallo [-6, 0) ⊈ SRGₑ(G₅) - Dopo la scalatura il punto -1 è separato dall'SRG 2. **Inversione e Scalatura**: k₁ = 0, k₂ = -2 - Invertire il piccolo buco destro sul punto -1 - Sfruttare le proprietà geometriche dell'SRG hard 3. **Traslazione**: k₁ = -1.5, k₂ = 1 - Traslare il piccolo buco a sinistra in modo che includa il punto -1 - Garantire la stabilità dell'anello chiuso Ogni caso soddisfa la condizione di stabilità: $$\text{dist}(-1, \text{SRG}_e(k_1I + k_2G_5)) \geq r > 0$$ ### Scoperte Chiave 1. **Precisione**: Per sistemi SISO e sistemi MIMO a matrice normale, il confine è preciso (non approssimato) 2. **Stabilità Numerica**: Il confine numerico risultante da un insieme finito di α è sempre un'approssimazione esterna, non sottostima mai l'SRG 3. **Intuizione Geometrica**: L'SRG hard fornisce un'interpretazione geometrica intuitiva, visualizzando il margine di stabilità 4. **Modularità**: Una volta calcolato l'SRG hard di G, può essere riutilizzato per diversi progetti di guadagno senza ricalcolo ## Lavori Correlati ### Sviluppo del Framework SRG 1. **SRG Originale** [2,3]: Proposto da Ryu et al. e Chaffey et al., basato su SRG soft nello spazio L₂ 2. **SRG Esteso** [4,11]: Esteso ai sistemi SISO instabili aggiungendo informazioni dal criterio di Nyquist 3. **SRG Hard** [12]: Proposto da Chen, definito su L₂[0,T], teoricamente in grado di gestire sistemi instabili ### Campi di Applicazione dell'SRG 1. **Analisi Non Lineare**: Criterio del cerchio generalizzato [4], teorema del piccolo guadagno [3] 2. **Controllo con Reset**: Analisi [8] e progetto [9] 3. **Sistemi Decentralizzati**: Condizioni di stabilità [10] 4. **Caratteristiche nel Dominio della Frequenza**: Diagrammi di Bode non lineari e definizioni di larghezza di banda [5] 5. **Analisi di Fase**: Concetti di anticipo/ritardo di fase [6] ### Relazione con i Metodi Classici 1. **Criterio di Nyquist**: L'SRG hard nel caso SISO è equivalente all'SRG esteso, che include informazioni di Nyquist 2. **Criterio MIMO di Nyquist** [13]: L'SRG hard fornisce un metodo alternativo senza necessità di calcolare il numero di avvolgimenti 3. **Criterio di Nyquist Generalizzato (GNC)** [14]: L'SRG hard evita il problema dell'assemblaggio manuale delle tracce caratteristiche ### Vantaggi di Questo Articolo 1. **Sistematizzazione**: Fornisce un algoritmo esplicito, non un metodo euristico 2. **Universalità**: Gestisce sistemi stabili, instabili e con integratori 3. **Precisione**: Metodo non approssimato, fornisce confini precisi 4. **Scalabilità**: Si estende naturalmente al caso MIMO ## Conclusioni e Discussione ### Conclusioni Principali 1. **Completezza dell'Algoritmo**: Sviluppa il primo metodo sistematizzato per calcolare con precisione l'SRG hard di sistemi LTI (inclusi sistemi instabili) 2. **Unificazione Teorica**: Dimostra l'equivalenza dell'SRG hard e dell'SRG esteso nel caso SISO, unificando i due framework 3. **Valore Pratico**: L'SRG hard fornisce un'alternativa pratica al criterio MIMO di Nyquist e al GNC, con i seguenti vantaggi: - Non è necessario calcolare il numero di avvolgimenti - L'SRG può essere riutilizzato quando il guadagno cambia - Fornisce direttamente i margini di guadagno L₂ e di stabilità 4. **Intuizione Geometrica**: La costruzione geometrica basata sull'intersezione di dischi fornisce una visualizzazione chiara e una comprensione intuitiva ### Limitazioni #### 1. Dipendenza dal Modello - Il metodo attuale si basa sulla rappresentazione della funzione di trasferimento, richiedendo un modello parametrico del sistema - Il vantaggio del criterio di Nyquist è la capacità di utilizzare direttamente i dati di risposta in frequenza misurati - **Esigenza Futura**: Calcolare rα e Rα direttamente dai dati ingresso/uscita #### 2. Sfida del Rilevamento della Fase Minima - La Proposizione 2 richiede solo il rilevamento dell'esistenza di poli instabili (non il conteggio) - La Proposizione 3 richiede il rilevamento degli zeri a fase non minima - **Difficoltà Nota**: Il rilevamento degli zeri a fase non minima dai dati è problematico [22] #### 3. Sistemi Non Square - La teoria completa attuale si applica solo a sistemi square (p = q) - Caso Wide (p > q): Gα ha sempre un nucleo algebrico non banale, portando a rα = 0 - Caso Tall (p < q): È possibile utilizzare la pseudoinversa sinistra, ma richiede molti dettagli tecnici aggiuntivi #### 4. Ipotesi di Buona Posizione - L'analisi di stabilità dell'SRG hard richiede l'assunzione della buona posizione del sistema di retroazione - L'SRG soft e l'SRG esteso garantiscono la buona posizione oltre ai margini di guadagno - Questo è uno svantaggio dell'SRG hard rispetto all'SRG esteso ### Direzioni Future #### 1. Metodi Guidati dai Dati Calcolare l'SRG hard direttamente da dati di risposta in frequenza o dati ingresso/uscita, in particolare: - Rilevare poli instabili dai dati - Identificare zeri a fase non minima dai dati #### 2. Estensione a Sistemi Non Square - Perfezionare la teoria per sistemi tall (p < q) - Sviluppare metodi alternativi per sistemi wide #### 3. Applicazioni a Sistemi Non Lineari - Applicare i risultati all'analisi di sistemi MIMO non lineari [7] - Combinare con il framework di analisi della retroazione non lineare dell'SRG hard in [12] #### 4. Ottimizzazione Computazionale - Sviluppare algoritmi numerici efficienti - Affrontare la complessità computazionale per sistemi MIMO ad alta dimensione - Selezione adattiva dell'insieme α per bilanciare precisione e costo computazionale #### 5. Analisi di Robustezza - Estendere a sistemi incerti - Integrare con la teoria del controllo robusto ## Valutazione Approfondita ### Punti di Forza #### 1. Rigore Teorico - **Prove Matematiche Complete**: Tutti i risultati principali (Proposizioni 1-3, Teoremi 2-3) hanno prove rigorose - **Fondamento Geometrico**: Basato su una solida base teorica della convessità iperbolica - **Prova di Equivalenza**: Il Teorema 3 stabilisce un collegamento esplicito con i metodi classici #### 2. Innovazione del Metodo - **Framework Unificato**: Primo metodo unificato per gestire sistemi LTI stabili e instabili - **Trucco Elegante del Sistema Inverso**: Calcolare rα attraverso l'inversione di Gα è una soluzione elegante - **Costruzione Geometrica**: Il metodo di intersezione dei dischi è sia intuitivo che facile da implementare #### 3. Valore Pratico - **Algoritmo Implementabile**: Il metodo si basa su strumenti standard dei sistemi lineari (valori singolari, rilevamento della fase minima) - **Progetto Modulare**: L'SRG calcolato può essere riutilizzato per diversi progetti di guadagno - **Visualizzazione Intuitiva**: La rappresentazione geometrica nel piano complesso facilita la comprensione e l'applicazione #### 4. Adeguatezza degli Esempi - Copre sistemi SISO e MIMO, stabili e instabili, con integratori - Dimostra praticamente l'applicazione al progetto del guadagno - Visualizzazione chiara ed efficace ### Insufficienze #### 1. Limitazioni Teoriche - **Restrizione Square**: I risultati principali si applicano solo a sistemi square, limitando l'applicabilità - **Buona Posizione**: Non può garantire la buona posizione come l'SRG esteso - **Ipotesi di Matrice Normale**: Per sistemi MIMO generali, il confine potrebbe non essere preciso #### 2. Sfide Pratiche - **Dipendenza dal Modello**: Richiede un modello di funzione di trasferimento preciso - **Rilevamento della Fase Minima**: Il rilevamento dai dati è problematico - **Complessità Computazionale**: Per sistemi MIMO ad alta dimensione, è necessario calcolare su un continuo di α #### 3. Insufficienza Sperimentale - **Mancanza di Analisi del Tempo di Calcolo**: Non sono riportati i tempi di calcolo dell'algoritmo - **Discussione sulla Precisione Numerica**: L'errore numerico dovuto a un insieme finito di α non è quantificato - **Confronto Pratico con GNC**: Solo confronto teorico, mancano casi pratici di confronto #### 4. Problemi di Presentazione - **Simboli Densi**: L'abbondanza di simboli matematici potrebbe influire sulla leggibilità - **Mancanza di Pseudocodice**: Sebbene il metodo sia chiaro, manca una descrizione dell'algoritmo direttamente implementabile - **Spiegazione Grafica**: Alcuni grafici (come gli artefatti numerici nella figura 2b) richiedono maggiore spiegazione ### Valutazione dell'Impatto #### Contributo al Campo 1. **Completezza Teorica**: Colma il vuoto nel framework SRG nel gestire sistemi instabili 2. **Progresso Metodologico**: Fornisce un nuovo strumento per l'analisi di sistemi non lineari 3. **Modernizzazione della Teoria Classica**: Fornisce un'interpretazione geometrica moderna del criterio di Nyquist #### Valore Pratico 1. **Progetto di Controllo**: Applicabile direttamente al progetto di sistemi di retroazione con componenti instabili 2. **Analisi di Stabilità**: Fornisce una visualizzazione del margine di stabilità più intuitiva rispetto ai metodi tradizionali 3. **Valore Didattico**: Il metodo geometrico aiuta a comprendere il comportamento dei sistemi di retroazione #### Riproducibilità - **Riproducibilità Teorica**: Tutte le prove sono completamente fornite nell'appendice - **Fattibilità di Implementazione**: Basato su strumenti standard, facile da implementare - **Chiarezza degli Esempi**: Fornisce funzioni di trasferimento concrete, facilitando la verifica - **Mancanza di Codice**: Non è fornita un'implementazione open source (questo è un peccato) ### Scenari di Applicabilità #### Applicazioni Più Adatte 1. **Progetto di Sistemi di Retroazione Complessi**: Specialmente sistemi contenenti componenti instabili 2. **Analisi di Robustezza**: Quando è necessario visualizzare il margine di stabilità 3. **Sistemi Multivariabili**: Regolazione del guadagno e analisi di stabilità di sistemi MIMO 4. **Sistemi Non Lineari**: Come blocco costruttivo fondamentale per l'analisi SRG non lineare #### Scenari Meno Adatti 1. **Puramente Guidato dai Dati**: Quando non è possibile ottenere un modello preciso 2. **Sistemi Non Square**: Quando p ≠ q il metodo è incompleto 3. **Applicazioni in Tempo Reale**: Se la complessità computazionale è un vincolo critico 4. **Sistemi Estremamente Ad Alta Dimensione**: Quando la dimensione è molto elevata il calcolo potrebbe diventare difficile ### Valutazione Complessiva Questo è un **articolo teorico di controllo di alta qualità** che fornisce contributi importanti al framework del Grafico Relativo Scalato. I principali vantaggi sono: - Fornisce il primo metodo sistematizzato per gestire l'SRG hard di sistemi LTI instabili - La teoria è rigorosa e le prove sono complete - Stabilisce un collegamento chiaro con il criterio classico di Nyquist - Pone le basi per l'analisi di sistemi non lineari Principali spazi di miglioramento: - Estensione a sistemi non square - Sviluppo di versioni guidate dai dati - Fornitura di implementazione open source e più casi pratici **Indice di Raccomandazione**: ⭐⭐⭐⭐☆ (4/5) Questo lavoro ha un valore di riferimento importante per i ricercatori che lavorano nel controllo non lineare, controllo robusto e analisi di sistemi, in particolare per coloro che si interessano del metodo SRG o che hanno bisogno di gestire sistemi instabili. ## Riferimenti (Selezionati) 1. **[2] Ryu et al., 2022**: "Scaled relative graphs: Nonexpansive operators via 2D Euclidean geometry" - Definizione matematica originale dell'SRG 2. **[3] Chaffey et al., 2023**: "Graphical Nonlinear System Analysis" - Introduzione dell'SRG nella teoria del controllo 3. **[11] Krebbekx et al., 2025**: "Scaled Relative Graph Analysis of General Interconnections" - Definizione dell'SRG esteso 4. **[12] Chen, 2025**: "Soft and Hard Scaled Relative Graphs" - Prima definizione dell'SRG hard 5. **[13] MacFarlane & Postlethwaite, 1977**: Letteratura classica sul criterio MIMO di Nyquist 6. **[14] Desoer & Wang, 1980**: Criterio di Nyquist Generalizzato (GNC)