Salta al contenuto
Note per Studenti Esercizio 20 · sintesi con rete a sella e PID di tipo 1 su due poli reali (temi d'esame luglio 2024, settembre 2023 e gennaio 2021)

Esercizio 20sintesi con rete a sella e PID di tipo 1 su due poli reali (temi d'esame luglio 2024, settembre 2023 e gennaio 2021)

Esame
In questa pagina 3

Testo (temi d'esame luglio 2024, terzo appello, esercizio 3; stesso testo nei temi di settembre 2023 e di gennaio 2021). Dato il sistema di funzione di trasferimento G(s)=1(1+s)(1+10s),G(s)=\frac1{(1+s)(1+10s)}, è richiesto i) il progetto di un controllore stabilizzante C1(s)C_1(s) proprio che garantisca che il sistema retroazionato sia di tipo 00, con errore di regime permanente al gradino erp(1)≃10−3e_{rp}^{(1)}\simeq10^{-3}, mentre il sistema in catena aperta abbia pulsazione di attraversamento ωA≃1\omega_A\simeq1 rad/s e margine di fase mφ≃90∘m_\varphi\simeq90^\circ; ii) il progetto di un controllore stabilizzante C2(s)C_2(s) di tipo PID (eventualmente solo P, PI o PD) che garantisca che il sistema retroazionato sia di tipo 11, con errore di regime permanente alla rampa lineare erp(2)≃0,01e_{rp}^{(2)}\simeq0{,}01, mentre il sistema in catena aperta abbia ωA≃1000\omega_A\simeq1000 rad/s e mφ≃90∘m_\varphi\simeq90^\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 →, Reti attenuatrici e reti a sellaLa rete attenuatrice $C_{att}=\frac{1+s/z}{1+s/p}$ con $0<p<z$ ha $K_B=1$ e abbassa il modulo di $20\log\frac zp$ dB oltre $z$ senza cambiare il guadagno in continua (a costo di un po' di fase persa, circa $\frac z{\omega_A^}$ rad). Caso ③: polo a $\frac M{20}$ decadi prima dello zero, zero $\ge1$ decade prima di $\omega_A^$. La rete a sella $C_{sella}=C_{att}C_{ant}$ ($p_1<z_1\le z_2<p_2$) abbassa il modulo e alza la fase: si progetta prima la parte anticipatrice (fase, con un guadagno extra $\Delta$), poi quella attenuatrice sufficientemente a sinistra per togliere $M+\Delta$ dB senza toccare la fase (caso ④).Reti attenuatrici e reti a sella →, 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 →, 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(0)=1≠0G(0)=1\ne0, nessun polo instabile (poli −1-1 e −0,1-0{,}1): KB(G)=1K_B(G)=1, r=0r=0.

Punto (i): tipo 0, erp(1)≃10−3e_{rp}^{(1)}\simeq10^{-3}, ωA∗=1\omega_A^*=1, mφ∗=90∘m_\varphi^*=90^\circ

Passo 1. Tipo 00: 11+KB(C′)≤10−3⇒KB(C′)≥1KB(G)(110−3−1)=999\frac1{1+K_B(C')}\le10^{-3}\Rightarrow K_B(C')\ge\frac1{K_B(G)}\left(\frac1{10^{-3}}-1\right)=999. Si sceglie C′=1000C'=1000 (errore 11001=9,99⋅10−4\frac1{1001}=9{,}99\cdot10^{-4}).

Passo 2. C′G=1000(1+s)(1+10s)C'G=\dfrac{1000}{(1+s)(1+10s)}. In ω=ωA∗=1\omega=\omega_A^*=1: modulo 20log⁡10002101=36,9520\log\frac{1000}{\sqrt2\sqrt{101}}=36{,}95 dB (quindi la pulsazione di attraversamento di C′GC'G è maggiore di 11: vale circa 1010 rad/s) e fase −arctan⁡1−arctan⁡10=−45∘−84,3∘=−129,3∘-\arctan1-\arctan10=-45^\circ-84{,}3^\circ=-129{,}3^\circ, cioè mφ(ωA∗)=50,7∘<90∘m_\varphi(\omega_A^*)=50{,}7^\circ<90^\circ. Caso ④: abbassare il modulo di M=36,95M=36{,}95 dB e alzare la fase di Φ=90∘−50,7∘=39,3∘\Phi=90^\circ-50{,}7^\circ=39{,}3^\circ: rete a sella.

Rete anticipatrice (per la fase): zero in 11 (cancella il polo di GG in −1-1, stabile) e polo in 100100: Cant=1+s1+s100C_{ant}=\dfrac{1+s}{1+\frac s{100}}. In ω=1\omega=1 aggiunge Δ=3,0\Delta=3{,}0 dB e +44,4∘+44{,}4^\circ (45∘−0,57∘45^\circ-0{,}57^\circ), più dei 39,3∘39{,}3^\circ necessari.

Rete attenuatrice che deve togliere M+Δ=36,95+3,0≈40M+\Delta=36{,}95+3{,}0\approx40 dB: zero in 0,10{,}1 (una decade prima di ωA∗\omega_A^*, e cancella il polo di GG in −0,1-0{,}1) e polo due decadi prima dello zero (4040 dB =2=2 decadi), cioè in 10−310^{-3}: Catt=1+10s1+1000sC_{att}=\dfrac{1+10s}{1+1000s} (perdita di fase arctan⁡10−arctan⁡1000=−5,65∘\arctan10-\arctan1000=-5{,}65^\circ).

Controllore (rete a sella) proprio: C1(s)=1000 (1+s)(1+10s)(1+s100)(1+1000s).C_1(s)=1000\,\frac{(1+s)(1+10s)}{\left(1+\frac s{100}\right)(1+1000s)}. Verifica. Le due cancellazioni zero/polo coinvolgono poli stabili di GG e lasciano C1G=1000(1+s100)(1+1000s).C_1G=\frac{1000}{\left(1+\frac s{100}\right)(1+1000s)}. In ω=1\omega=1: modulo 10001+10−41+106≈1,000\frac{1000}{\sqrt{1+10^{-4}}\sqrt{1+10^6}}\approx1{,}000 (ωA=0,9999\omega_A=0{,}9999 rad/s) e fase −arctan⁡(0,01)−arctan⁡(1000)=−0,57∘−89,94∘=−90,5∘-\arctan(0{,}01)-\arctan(1000)=-0{,}57^\circ-89{,}94^\circ=-90{,}5^\circ: mφ=89,5∘≃90∘m_\varphi=89{,}5^\circ\simeq90^\circ ✓. Tipo 00, erp(1)=11001≃10−3e_{rp}^{(1)}=\frac1{1001}\simeq10^{-3} ✓. KB>0K_B>0, poli stabili, strettamente propria, ωA\omega_A unica: BIBO stabile per il criterio di Bode (poli di WW: −98,99-98{,}99 e −1,011-1{,}011, più i due modi cancellati −1-1 e −0,1-0{,}1).

Punto (ii): tipo 1, erp(2)≃0,01e_{rp}^{(2)}\simeq0{,}01, ωA∗=1000\omega_A^*=1000, mφ∗=90∘m_\varphi^*=90^\circ

Passo 1. Tipo 11, r=0r=0: C′=KB(C′)sC'=\frac{K_B(C')}s con ∣erp(2)∣=1KB(C′)KB(G)≤0,01⇒KB(C′)=100|e_{rp}^{(2)}|=\frac1{K_B(C')K_B(G)}\le0{,}01\Rightarrow K_B(C')=100: C′=100sC'=\frac{100}s.

Passo 2. C′G=100s(1+s)(1+10s)C'G=\dfrac{100}{s(1+s)(1+10s)}. In ω=1000\omega=1000: modulo ≈100⋅11000⋅1000⋅104=10−8\approx100\cdot\frac1{1000\cdot1000\cdot10^4}=10^{-8}, cioè −160-160 dB (M=160M=160 dB da guadagnare), fase −90∘−arctan⁡1000−arctan⁡104=−269,94∘-90^\circ-\arctan1000-\arctan10^4=-269{,}94^\circ, cioè mφ(ωA∗)=−89,94∘m_\varphi(\omega_A^*)=-89{,}94^\circ: Φ=90∘+89,94∘=179,94∘>90∘\Phi=90^\circ+89{,}94^\circ=179{,}94^\circ>90^\circ: caso ② con due zeri ciascuno a M40=4\frac M{40}=4 decadi prima di ωA∗\omega_A^*, cioè in ω=0,1\omega=0{,}1 rad/s. Poiché il PID non è necessariamente proprio, non si aggiungono poli: C2(s)=100 (1+10s)2s=100+2000 s+104s2s=100⋅1s+2000+104s.C_2(s)=\frac{100\,(1+10s)^2}s=\frac{100+2000\,s+10^4s^2}s=100\cdot\frac1s+2000+10^4s. PID con Ki=100K_i=100, Kp=2000K_p=2000, Kd=104K_d=10^4 (zeri reali coincidenti in −0,1-0{,}1: Kp2=4⋅106=4KiKdK_p^2=4\cdot10^6=4K_iK_d, caso limite).

Verifica. Uno dei due zeri cancella il polo di GG in −0,1-0{,}1 (stabile): C2G=100 (1+10s)s (1+s)C_2G=\dfrac{100\,(1+10s)}{s\,(1+s)}. In ω=1000\omega=1000: modulo 100⋅1041000⋅1000=1,000\frac{100\cdot10^4}{1000\cdot1000}=1{,}000 (ωA≃1000\omega_A\simeq1000), fase −90∘+arctan⁡(104)−arctan⁡(1000)=−90∘+89,994∘−89,943∘=−89,95∘-90^\circ+\arctan(10^4)-\arctan(1000)=-90^\circ+89{,}994^\circ-89{,}943^\circ=-89{,}95^\circ: mφ=90,05∘≥90∘m_\varphi=90{,}05^\circ\ge90^\circ ✓. Tipo 11, erp(2)=1KB(C2G)=1100e_{rp}^{(2)}=\frac1{K_B(C_2G)}=\frac1{100} ✓. BIBO stabile per Bode (poli di WW: −1000,9-1000{,}9 e −0,0999-0{,}0999, più il modo cancellato −0,1-0{,}1).

Errori comuni

  • Nel punto (i) usare una sola rete anticipatrice: il modulo in ω∗\omega^* (+37+37 dB) andrebbe abbassato, non alzato; serve la sella.
  • Dimenticare di sommare alla quota da attenuare anche il guadagno Δ\Delta della parte anticipatrice (36,95+3,0≃4036{,}95+3{,}0\simeq40 dB).
  • Mettere un solo zero nel PID del punto (ii): Φ≈180∘\Phi\approx180^\circ richiede due zeri.
  • Dimenticare che il PID ha zeri e non poli: non è proprio (va bene se il testo non lo chiede, come qui: "eventualmente P, PI, PD").
  • Cancellare con gli zeri della rete un polo instabile del processo: qui i poli cancellati (−1-1 e −0,1-0{,}1) sono stabili.

Versione ripasso

Testo. G=1(1+s)(1+10s)G=\frac1{(1+s)(1+10s)}: (i) proprio, tipo 0, e(1)≃10−3e^{(1)}\simeq10^{-3}, ωA∗=1\omega_A^*=1, m∗=90∘m^*=90^\circ; (ii) PID tipo 1, e(2)≃0,01e^{(2)}\simeq0{,}01, ωA∗=1000\omega_A^*=1000, m∗=90∘m^*=90^\circ (luglio 2024, settembre 2023, gennaio 2021).

Teoria collegata