Salta al contenuto
Note per Studenti Esercizio 22 · sintesi di tipo 2 con anticipo di fase (tema d'esame gennaio 2023)

Esercizio 22sintesi di tipo 2 con anticipo di fase (tema d'esame gennaio 2023)

Esame
In questa pagina 3

Testo (tema d'esame gennaio 2023, primo appello, esercizio 3 punto i; stesso tipo di esercizio nei temi di febbraio 2020 e luglio 2021 con G(s)=s+10sG(s)=\frac{s+10}s). Dato il processo di funzione di trasferimento G(s)=s+100s,G(s)=\frac{s+100}s, è richiesto il progetto di un controllore stabilizzante C(s)C(s) che attribuisca al sistema retroazionato W=CG1+CGW=\dfrac{CG}{1+CG} tipo 22 ed errore a regime alla rampa parabolica erp(3)e_{rp}^{(3)} pari a circa 0,10{,}1, e alla funzione di trasferimento in catena aperta C(s)G(s)C(s)G(s) pulsazione di attraversamento ωA≃ωA∗=10\omega_A\simeq\omega_A^*=10 rad/s e margine di fase mφ≥mφ∗=80∘m_\varphi\ge m_\varphi^*=80^\circ.


Teoria usata: Problema di regolazione e sintesi per tentativiProblema di regolazione: data $G$ propria, senza poli in $\mathrm{Re}>0$, con $G(0)\ne0$ (cioè $G=\frac{K_B(G)}{s^r}G_1$, $G_1(0)=1$), progettare $C$ propria tale che $W=\frac{CG}{1+CG}$ sia BIBO stabile, di tipo $k$ con $|e_{rp}^{(k+1)}|\le e^$, con $\omega_A\simeq\omega_A^$ e $m_\varphi\ge m_\varphi^$. Sintesi per tentativi: Passo 1 $C'=\frac{K_B(C')}{s^\alpha}$ ($\alpha=\max(k-r,0)$, $|K_B(C')|$ minimo che soddisfa l'errore, segno di $K_B(G)$); Passo 2 $C''$ con $K_B(C'')=1$ che porta $\omega_A$ e $m_\varphi$ ai valori voluti in quattro casi ($\omega_A\gtrless\omega_A^$, $m_\varphi(\omega_A^)\gtrless m_\varphi^$): guadagno, rete anticipatrice, attenuatrice, a sella. Il criterio di Bode garantisce la stabilità.Problema di regolazione e sintesi per tentativi →, Reti anticipatriciLa rete anticipatrice $C_{ant}(s)=\frac{1+s/z}{1+s/p}=\frac{1+sT}{1+\alpha Ts}$ con $0<z<p$ ($\alpha=\frac zp<1$) ha $K_B=1$, alza il modulo di $20\log\frac pz$ dB tra $z$ e $p$ e aggiunge fase (massimo $\arcsin\frac{1-\alpha}{1+\alpha}$ in $\sqrt{zp}$). Si usa nei casi ②, ④ (parte anticipatrice) della sintesi: con $\Phi\le90^\circ$ si mette uno zero $\frac M{20}$ decadi prima di $\omega_A^$; con $90^\circ<\Phi\le180^\circ$ due zeri a $\frac M{40}$ decadi; il polo (o i poli) almeno una o due decadi dopo $\omega_A^$ per rendere la rete propria.Reti anticipatrici →, Controllori PIDIl controllore PID calcola $u=K_pe+K_d\dot e+K_i\int e$, cioè $C_{PID}(s)=K_p+K_ds+\frac{K_i}s=\frac{K_i}s\Big(1+\frac{K_p}{K_i}s+\frac{K_d}{K_i}s^2\Big)$ (due zeri e un polo nell'origine: non proprio). Nella sintesi per tentativi $\frac{K_i}s$ alza il tipo (Passo 1) e i due zeri sono una doppia azione anticipatrice (Passo 2): $PI$ = un integratore + uno zero (anticipo fino a $90^\circ$), $PID$ = integratore + due zeri (anticipo fino a $180^\circ$), $PD$ = zero senza integratore. Gli zeri si mettono $\frac M{20}$ (PI) o $\frac M{40}$ decadi ciascuno (PID) prima di $\omega_A^*$.Controllori PID →, Tipo e errore a regime di un sistema retroazionatoSe $W=\frac{\tilde G}{1+\tilde G}$ (con $\tilde G$ propria) è propria, BIBO stabile e $W(0)\ne0$, il tipo di $W$ coincide con il numero $h$ di poli nell'origine di $\tilde G$ (integratori in catena aperta). L'errore a regime corrispondente è $e_{rp}^{(h+1)}=\frac1{K_B(\tilde G)}$ per $h\ge1$ e $e_{rp}^{(1)}=\frac1{1+K_B(\tilde G)}$ per $h=0$. Per ridurre l'errore si alza il guadagno di Bode o si aggiungono integratori, compatibilmente con la stabilità.Tipo e errore a regime di un sistema retroazionato →, Pulsazione di attraversamento, margine di fase e criterio di BodeLa pulsazione di attraversamento $\omega_A$ è la $\omega>0$ in cui $|G(j\omega)|=1$ (0 dB); il margine di fase è $m_\varphi=180^\circ+\arg G(j\omega_A)$: quanta fase si può ancora perdere prima di arrivare a $-180^\circ$. Criterio di Bode: se $G$ è strettamente propria, ha guadagno di Bode $K_B>0$, nessun polo con $\mathrm{Re}>0$ e $\omega_A$ è ben definita (esiste ed è unica), allora $W=\frac G{1+G}$ è BIBO stabile $\iff m_\varphi>0$. Più $m_\varphi$ è grande, più il sistema è lontano dall'instabilità e meno oscilla.Pulsazione di attraversamento, margine di fase e criterio di Bode →.

Ipotesi sul processo. G(s)=s+100s=100 (1+s100)sG(s)=\dfrac{s+100}s=\dfrac{100\,(1+\frac s{100})}{s}: KB(G)=100>0K_B(G)=100>0, r=1r=1 (un polo nell'origine), nessuno zero in 00 e nessun polo instabile: si può applicare la sintesi.

Passo 1: tipo 2 ed errore 0,10{,}1

Il tipo richiesto è k=2k=2 e GG ha già r=1r=1 integratore: servono α=k−r=1\alpha=k-r=1 integratori nel controllore. L'errore alla rampa parabolica con tipo 22 è ∣erp(3)∣=1∣KB(C′)KB(G)∣≤0,1⇒∣KB(C′)∣≥1100⋅0,1=0,1|e_{rp}^{(3)}|=\frac1{|K_B(C')K_B(G)|}\le0{,}1\Rightarrow|K_B(C')|\ge\frac1{100\cdot0{,}1}=0{,}1. Con segno di KB(G)K_B(G) (positivo): C′(s)=0,1s.C'(s)=\frac{0{,}1}s. Allora C′G=0,1 (s+100)s2=10(1+s100)s2C'G=\dfrac{0{,}1\,(s+100)}{s^2}=\dfrac{10\left(1+\frac s{100}\right)}{s^2}: due poli nell'origine, KB=10K_B=10.

Passo 2: ωA\omega_A e margine di fase

In ω=ωA∗=10\omega=\omega_A^*=10: modulo 20log⁡[101+0,01100]=−19,9620\log\left[\frac{10\sqrt{1+0{,}01}}{100}\right]=-19{,}96 dB, fase −180∘+arctan⁡0,1=−174,3∘-180^\circ+\arctan0{,}1=-174{,}3^\circ, cioè mφ(ωA∗)=5,7∘m_\varphi(\omega_A^*)=5{,}7^\circ. Per ω<10\omega<10 il modulo è maggiore di 00 dB e in ω=10\omega=10 è già sotto (la pulsazione di attraversamento di C′GC'G è circa 3,23{,}2 rad/s <ωA∗<\omega_A^*): caso ② con M=19,96≈20M=19{,}96\approx20 dB da guadagnare e Φ=80∘−5,7∘=74,3∘<90∘\Phi=80^\circ-5{,}7^\circ=74{,}3^\circ<90^\circ: un solo zero, collocato M20=1\frac M{20}=1 decade prima di ωA∗\omega_A^*, cioè in ω=1\omega=1 rad/s: C′′=1+sC''=1+s. Il prodotto C′C′′=0,1(1+s)sC'C''=\frac{0{,}1(1+s)}s è già proprio (grado 1/11/1); la catena aperta CGCG ha però grado relativo 00 (bipropria, CG→0,1CG\to0{,}1 a frequenza infinita): per applicare il criterio di Bode (che richiede G~\tilde G strettamente propria) si aggiunge un polo lontano, a 10310^3 rad/s (22 decadi dopo ωA∗\omega_A^*, costo di fase 0,6∘0{,}6^\circ).

Controllore. C(s)=0,1 (1+s)s(1+s1000)(un PI con filtro:  Kp=0,1, Ki=0,1).C(s)=\frac{0{,}1\,(1+s)}{s\left(1+\frac s{1000}\right)}\qquad\Big(\text{un PI con filtro: }\ K_p=0{,}1,\ K_i=0{,}1\Big). Verifica. CG=0,1 (1+s)(s+100)s2(1+s1000)CG=\dfrac{0{,}1\,(1+s)(s+100)}{s^2\left(1+\frac s{1000}\right)}. In ω=10\omega=10: modulo 0,1⋅10,05⋅100,5100=1,010{,}1\cdot\frac{10{,}05\cdot100{,}5}{100}=1{,}01 (ωA=10,10\omega_A=10{,}10 rad/s) e fase −180∘+arctan⁡10+arctan⁡0,1−arctan⁡0,01=−180∘+84,29∘+5,71∘−0,57∘=−90,57∘-180^\circ+\arctan10+\arctan0{,}1-\arctan0{,}01=-180^\circ+84{,}29^\circ+5{,}71^\circ-0{,}57^\circ=-90{,}57^\circ: mφ=89,4∘≥80∘m_\varphi=89{,}4^\circ\ge80^\circ ✓. Tipo 22 (due integratori in CGCG) ed erp(3)=1KB(CG)=10,1⋅100=0,1e_{rp}^{(3)}=\frac1{K_B(CG)}=\frac1{0{,}1\cdot100}=0{,}1 ✓. KB(CG)=10>0K_B(CG)=10>0, poli in 00 (doppio) e −1000-1000 (nessuno in Re>0\mathrm{Re}>0), strettamente propria, ωA\omega_A unica: BIBO stabile per il criterio di Bode; poli di WW: −1090,7-1090{,}7, −8,12-8{,}12, −1,13-1{,}13 (verificato).

Errori comuni

  • Imporre tipo 22 con un guadagno qualunque: il tipo viene dal polo nell'origine in C′C' (qui un solo integratore in più, perché GG ne ha già uno).
  • Calcolare l'errore alla rampa parabolica con la formula del tipo 00 (11+KB\frac1{1+K_B}): per il tipo 22 è 1KB(CG)\frac1{K_B(CG)} con KB(CG)=KB(C′)KB(G)K_B(CG)=K_B(C')K_B(G).
  • Dimenticare lo zero: con il solo C′=0,1sC'=\frac{0{,}1}s la fase è −174∘-174^\circ in ωA∗\omega_A^*, margine di 5,7∘5{,}7^\circ, ben sotto 80∘80^\circ.
  • Dimenticare che il criterio di Bode richiede catena aperta strettamente propria: qui serve il polo lontano (altrimenti CG(∞)=0,1≠0CG(\infty)=0{,}1\ne0).
  • Confondere il PI con un PD: qui lo zero è in −1-1 e c'è il polo nell'origine (il tipo aumenta).

Versione ripasso

Testo. G=s+100sG=\frac{s+100}s: CC con tipo 2, e(3)≃0,1e^{(3)}\simeq0{,}1, ωA∗=10\omega_A^*=10, m∗≥80∘m^*\ge80^\circ (gennaio 2023).

Teoria collegata