Esercizio 6diagramma di Bode e di Nyquist di una funzione con zeri immaginari e poli instabili (tema d'esame luglio 2025)
In questa pagina 5
Testo (tema d'esame luglio 2025, terzo appello, esercizio 1; lo stesso esercizio è nel tema di gennaio 2023). Data la funzione di trasferimento i) si determini il diagramma di Bode (modulo e fase) della risposta in frequenza del sistema; ii) si determini il diagramma di Nyquist di per , calcolando eventuali asintoti ed intersezioni con gli assi; iii) si studi attraverso il criterio di Nyquist la stabilità BIBO del sistema retroazionato al variare di , . Si verifichi infine che lo studio con Routh conduce agli stessi risultati.
Teoria usata: Forma di Bode e diagrammi di Bode dei termini elementariLa forma di Bode scrive $W(s)=\frac{K_B}{s^h}\frac{\prod(1+T_i's)\prod(1+2\xi_i'\frac s{\omega_{ni}'}+\frac{s^2}{\omega_{ni}'^2})}{\prod(1+T_is)\prod(1+2\xi_i\frac s{\omega_{ni}}+\frac{s^2}{\omega_{ni}^2})}$ con $K_B=\lim_{s\to0}s^hW(s)$. Il diagramma di Bode (modulo in dB e fase, su asse $\omega$ logaritmico) si ottiene sommando i contributi dei termini elementari: costante, $s^{\mp1}$ ($\mp20$ dB/dec, $\mp90^\circ$), binomio ($\pm20$ dB/dec dopo $1/|T|$, fase $\pm90^\circ$, opposta se instabile), trinomio ($\pm40$ dB/dec dopo $\omega_n$, fase $\pm180^\circ$, picco $20\log2\xi$).Forma di Bode e diagrammi di Bode dei termini elementari →, Tracciare i diagrammi di Bode di una funzione di trasferimentoPer tracciare Bode di $G(s)$: (1) forma di Bode ($K_B$, $h$, termini); (2) elenco ordinato delle pulsazioni di spezzamento; (3) modulo asintotico: parte con pendenza $-20h$ dB/dec passando per $20\log|K_B|$ in $\omega=1$ e cambia pendenza di $\pm20$ ($\pm40$) a ogni spezzamento; (4) fase: $\arg K_B-90^\circ h$ iniziale più i contributi di zeri e poli (opposti se instabili), con salti di $\pm180^\circ$ per radici immaginarie; (5) correzioni (picchi dei trinomi, 3 dB agli spezzamenti).Tracciare i diagrammi di Bode di una funzione di trasferimento →, Diagrammi di Nyquist - tracciamento, asintoti e intersezioniIl diagramma di Nyquist di $G(s)$ è la curva ${G(j\omega):\omega\in\mathbb R}$ nel piano complesso, percorsa per $\omega$ crescente; per $\omega<0$ è il simmetrico rispetto all'asse reale di quella per $\omega>0$. Si traccia combinando modulo e fase di Bode: per $\omega\to0^+$ parte da $K_B$ (se $h=0$) o da infinito con angolo $\arg K_B-90^\circ h$ (asintoto verticale $\mathrm{Re}=K_B(\sum T_i'-\sum T_i)$ per $h=1$), per $\omega\to\infty$ arriva nell'origine con angolo $-90^\circ(n-m)$. Le intersezioni con gli assi si trovano da $\mathrm{Im},G(j\omega)=0$ e $\mathrm{Re},G(j\omega)=0$.Diagrammi di Nyquist - tracciamento, asintoti e intersezioni →, Criterio di Nyquist - poli sull'asse immaginario e guadagno variabileSe $G$ ha poli sull'asse immaginario il diagramma va all'infinito e si chiude al finito con archi di raggio infinito percorsi in senso orario: un polo di molteplicità $\mu$ in $j\omega_0$ ($\omega_0\ge0$) aggiunge un arco di ampiezza $\mu\pi$ tra $G(j(\omega_0-\epsilon))$ e $G(j(\omega_0+\epsilon))$ (per $\omega_0=0$ gli estremi sono i rami $\omega\to0^\mp$). Poi vale ancora $N=n_{G+}-n_{W+}$, con $n_{G+}$ che non conta i poli immaginari. Per $W=\frac{KG}{1+KG}$, $K\in\mathbb R\setminus{0}$, il punto critico è $-\frac1K$: si conta $N$ in ciascun intervallo dell'asse reale tra le intersezioni del diagramma; quando $-\frac1K$ è sul diagramma $W$ ha poli immaginari (caso critico).Criterio di Nyquist - poli sull'asse immaginario e guadagno variabile →, Criterio di Routh-HurwitzIl criterio di Routh stabilisce se un polinomio $P(s)=p_ns^n+\dots+p_0$ è di Hurwitz senza calcolarne le radici. Si costruisce la tabella (righe $n,n-1,\dots,0$) con $\text{nuovo}=\frac{\text{pivot}\cdot a-\text{prec.}\cdot b}{\text{pivot}}$; se la tabella va a compimento, il numero di variazioni di segno nella prima colonna è il numero di radici con $\mathrm{Re}>0$, le permanenze il numero con $\mathrm{Re}<0$. Elemento nullo in prima colonna: $P$ non è Hurwitz; riga nulla: radici simmetriche (spesso immaginarie) date dal polinomio ausiliario. Con parametro $K$ nei coefficienti dà gli intervalli di stabilità.Criterio di Routh-Hurwitz →.
(i) Diagrammi di Bode
Forma di Bode. è già nella forma di Bode: ( dB), (un polo nell'origine). Gli altri termini:
- numeratore : coppia di zeri immaginari , trinomio con e ;
- denominatore e : due poli reali instabili in e in .
Modulo asintotico. Parte con pendenza dB/dec e passa per dB in (perché ). In la coppia di zeri dà dB/dec e il polo in dà : la pendenza passa da a . Il modulo resta quindi piatto a dB fino a ; in il secondo polo riporta la pendenza a dB/dec. Il diagramma reale ha in un picco di antirisonanza infinito (gli zeri sono sull'asse immaginario, ).
Fase. Parte da (polo nell'origine). Gli zeri immaginari danno un salto di in . I poli instabili danno ciascuno una variazione positiva di (modulo uguale a quello dei poli stabili, fase opposta): il primo tra e , il secondo tra e . Quindi la fase sale da (valore esatto in e in ) fino a circa , subisce il salto e riparte da circa , per poi salire fino a (valori esatti in , in , in ).
| (rad/s) | |||||||
|---|---|---|---|---|---|---|---|
| (dB) | |||||||
| fase |
Grafico interattivo: Modulo asintotico (dB) di G in funzione di log10(ω): −20 dB/dec fino a ω = 1, piatto a 0 dB fino a ω = 10, poi −20 dB/dec (in ω = 1 il modulo reale ha un'antirisonanza infinita)
Grafico interattivo: Fase (gradi) di G in funzione di log10(ω): da −90° con un salto di +180° in ω = 1 (zeri immaginari) e poi fino a 270° (due poli instabili)
(ii) Diagramma di Nyquist
Per : . Razionalizzando si ottiene (La stessa cosa in forma equivalente: , .)
Comportamento agli estremi.
- : e (, : parte da ). Asintoto verticale (conferma: con ).
- : : arriva nell'origine con fase , cioè tangente all'asse immaginario negativo, dal terzo quadrante.
Intersezioni con gli assi. La parte reale si annulla solo per e la parte immaginaria si annulla per e per :
- : : il diagramma passa per l'origine ().
- : e . Intersezione con il semiasse reale negativo in .
Tracciato. Per : , (quarto quadrante): il diagramma scende da e arriva all'origine in con pendenza circa (rapporto ). Per : , (secondo quadrante): ruota in senso antiorario e taglia l'asse reale negativo in . Per : , (terzo quadrante): torna all'origine con tangente verticale. Il tratto per è il simmetrico rispetto all'asse reale.
(iii) Stabilità di con Nyquist
Dati. Poli di : (semplice), e : (il polo nell'origine non conta). Il polo nell'origine () richiede la chiusura al finito con una semicirconferenza di raggio infinito percorsa in senso orario (Criterio di Nyquist - poli sull'asse immaginario e guadagno variabileSe $G$ ha poli sull'asse immaginario il diagramma va all'infinito e si chiude al finito con archi di raggio infinito percorsi in senso orario: un polo di molteplicità $\mu$ in $j\omega_0$ ($\omega_0\ge0$) aggiunge un arco di ampiezza $\mu\pi$ tra $G(j(\omega_0-\epsilon))$ e $G(j(\omega_0+\epsilon))$ (per $\omega_0=0$ gli estremi sono i rami $\omega\to0^\mp$). Poi vale ancora $N=n_{G+}-n_{W+}$, con $n_{G+}$ che non conta i poli immaginari. Per $W=\frac{KG}{1+KG}$, $K\in\mathbb R\setminus{0}$, il punto critico è $-\frac1K$: si conta $N$ in ciascun intervallo dell'asse reale tra le intersezioni del diagramma; quando $-\frac1K$ è sul diagramma $W$ ha poli immaginari (caso critico).Criterio di Nyquist - poli sull'asse immaginario e guadagno variabile →): il ramo finisce a (con ) e il ramo parte da : l'arco passa per sull'asse reale. Il punto critico è .
Attraversamenti dell'asse reale del diagramma chiuso (verso: dal basso verso l'alto, dall'alto verso il basso):
- (per ): due volte, verso il basso ( ciascuna);
- (): due volte, verso l'alto ( ciascuna); in più il diagramma chiuso passa per l'origine anche nel punto di giunzione tra (arriva dal terzo quadrante, ) e (riparte nel secondo quadrante, ): un terzo attraversamento verso l'alto ();
- all'infinito positivo (arco di chiusura): una volta verso il basso ().
Contando con la semiretta a destra del punto critico (, quindi ):
| Posizione di | |||
|---|---|---|---|
Per il diagramma passa per il punto critico: ha due poli immaginari (). Conclusione: è BIBO stabile se e solo se . Per ha due poli instabili, per ne ha tre.
Verifica con Routh
Il polinomio caratteristico è con e ; moltiplicando per : . Tabella:
- : prima colonna : nessuna variazione, stabile ✓.
- : e (quindi il terzo elemento è positivo): segni : due variazioni, 2 poli instabili ✓.
- : ma : segni : ancora due variazioni (2 poli instabili) ✓. In il secondo elemento si annulla ma, come si vede da Nyquist, in quel punto non succede nulla di speciale (per appena maggiore o minore di le variazioni sono sempre ).
- : , il terzo elemento è positivo (numeratore e denominatore negativi) e : segni , tre variazioni, 3 poli instabili ✓.
- : la riga si annulla; il polinomio ausiliario dà ✓.
Controllo numerico dei poli: : e (2 instabili); : e (stabile); : e (stabile); : e (3 instabili). I due metodi danno gli stessi intervalli di stabilità e lo stesso numero di poli instabili, come richiesto dal testo.
Errori comuni
- Dimenticare la chiusura all'infinito per il polo nell'origine, o disegnarla antioraria.
- Contare nei il polo nell'origine: i poli sull'asse immaginario non si contano. Qui .
- Attribuire fase negativa ai poli instabili: sono ciascuno.
- Dimenticare il salto di degli zeri immaginari nella fase, o l'antirisonanza infinita nel modulo.
- Confondere gli intervalli di (, , ) con quelli di : l'ordine si inverte.
Versione ripasso
Testo. : Bode; Nyquist (asintoti, intersezioni); stabilità di , ; verifica con Routh (luglio 2025, gennaio 2023).
- Bode (Tracciare i diagrammi di Bode di una funzione di trasferimentoPer tracciare Bode di $G(s)$: (1) forma di Bode ($K_B$, $h$, termini); (2) elenco ordinato delle pulsazioni di spezzamento; (3) modulo asintotico: parte con pendenza $-20h$ dB/dec passando per $20\log|K_B|$ in $\omega=1$ e cambia pendenza di $\pm20$ ($\pm40$) a ogni spezzamento; (4) fase: $\arg K_B-90^\circ h$ iniziale più i contributi di zeri e poli (opposti se instabili), con salti di $\pm180^\circ$ per radici immaginarie; (5) correzioni (picchi dei trinomi, 3 dB agli spezzamenti).Tracciare i diagrammi di Bode di una funzione di trasferimento →): , ; zeri (, ); poli instabili in , . Modulo: dB/dec, in piatto a dB fino a , poi ; antirisonanza infinita in . Fase: , salto in , due poli instabili ciascuno: fino a .
- Nyquist (Diagrammi di Nyquist - tracciamento, asintoti e intersezioniIl diagramma di Nyquist di $G(s)$ è la curva ${G(j\omega):\omega\in\mathbb R}$ nel piano complesso, percorsa per $\omega$ crescente; per $\omega<0$ è il simmetrico rispetto all'asse reale di quella per $\omega>0$. Si traccia combinando modulo e fase di Bode: per $\omega\to0^+$ parte da $K_B$ (se $h=0$) o da infinito con angolo $\arg K_B-90^\circ h$ (asintoto verticale $\mathrm{Re}=K_B(\sum T_i'-\sum T_i)$ per $h=1$), per $\omega\to\infty$ arriva nell'origine con angolo $-90^\circ(n-m)$. Le intersezioni con gli assi si trovano da $\mathrm{Im},G(j\omega)=0$ e $\mathrm{Re},G(j\omega)=0$.Diagrammi di Nyquist - tracciamento, asintoti e intersezioni →): , ; asintoto ; passa per in ; taglia l'asse reale in (); arriva nell'origine con tangente .
- Criterio (Criterio di Nyquist - poli sull'asse immaginario e guadagno variabileSe $G$ ha poli sull'asse immaginario il diagramma va all'infinito e si chiude al finito con archi di raggio infinito percorsi in senso orario: un polo di molteplicità $\mu$ in $j\omega_0$ ($\omega_0\ge0$) aggiunge un arco di ampiezza $\mu\pi$ tra $G(j(\omega_0-\epsilon))$ e $G(j(\omega_0+\epsilon))$ (per $\omega_0=0$ gli estremi sono i rami $\omega\to0^\mp$). Poi vale ancora $N=n_{G+}-n_{W+}$, con $n_{G+}$ che non conta i poli immaginari. Per $W=\frac{KG}{1+KG}$, $K\in\mathbb R\setminus{0}$, il punto critico è $-\frac1K$: si conta $N$ in ciascun intervallo dell'asse reale tra le intersezioni del diagramma; quando $-\frac1K$ è sul diagramma $W$ ha poli immaginari (caso critico).Criterio di Nyquist - poli sull'asse immaginario e guadagno variabile →): ; chiusura: semicirconferenza oraria per il polo in ; con semiretta a destra: : , ; : , ; : , . Stabile .
- Routh (Criterio di Routh-HurwitzIl criterio di Routh stabilisce se un polinomio $P(s)=p_ns^n+\dots+p_0$ è di Hurwitz senza calcolarne le radici. Si costruisce la tabella (righe $n,n-1,\dots,0$) con $\text{nuovo}=\frac{\text{pivot}\cdot a-\text{prec.}\cdot b}{\text{pivot}}$; se la tabella va a compimento, il numero di variazioni di segno nella prima colonna è il numero di radici con $\mathrm{Re}>0$, le permanenze il numero con $\mathrm{Re}<0$. Elemento nullo in prima colonna: $P$ non è Hurwitz; riga nulla: radici simmetriche (spesso immaginarie) date dal polinomio ausiliario. Con parametro $K$ nei coefficienti dà gli intervalli di stabilità.Criterio di Routh-Hurwitz →): , prima colonna : stesse conclusioni; : .
- Errori: chiusura dimenticata/antioraria; polo in in ; fase dei poli instabili; salto degli zeri immaginari; inversione degli intervalli di .