07.06.2013 Views

Riepilogo di tutto Bode con un PDF chiarissimo

Riepilogo di tutto Bode con un PDF chiarissimo

Riepilogo di tutto Bode con un PDF chiarissimo

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

0.0. 3.2 1<br />

Diagrammi <strong>di</strong> <strong>Bode</strong><br />

• Possibili rappresentazioni grafiche della f<strong>un</strong>zione <strong>di</strong> risposta armonica F (ω) =<br />

G(jω) sono: i Diagrammi <strong>di</strong> <strong>Bode</strong>, i Diagrammi <strong>di</strong> Nyquist e i Diagrammi<br />

<strong>di</strong> Nichols.<br />

• I Diagrammi <strong>di</strong> <strong>Bode</strong> sono due:<br />

1) il <strong>di</strong>agramma delle ampiezze rappresenta α = ln |G(jω)| in f<strong>un</strong>zione<br />

del ln ω;<br />

2) il <strong>di</strong>agramma delle fasi rappresenta β = arg G(jω) in f<strong>un</strong>zione del<br />

ln ω;<br />

• Esempio:<br />

ln G(jω) = ln<br />

Diagramma delle ampiezze<br />

G(jω) =<br />

<br />

100<br />

<br />

1 + j ω 10<br />

Diagramma delle fasi<br />

<br />

<br />

j arg G(jω)<br />

|G(jω)|e<br />

= ln |G(jω)| + j arg G(jω)<br />

1 + j ω 80<br />

<br />

<br />

1 + j ω<br />

1000<br />

<br />

= ln |G(jω)| + ln[e j arg G(jω) ]<br />

R. Zanasi, R. Morselli - Controlli Automatici - 2005/06 3. ANALISI ARMONICA


3.2. DIAGRAMMI DI BODE 3.2 2<br />

• Per i calcoli si utilizzano i logaritmi naturali. Peraltro, <strong>un</strong> cambiamento <strong>di</strong><br />

base dei logaritmi equivale ad <strong>un</strong> semplice cambiamento <strong>di</strong> scala.<br />

• Si possono utilizzare due tipi <strong>di</strong>versi <strong>di</strong> carta millimetrata:<br />

a) carta <strong>con</strong> doppia scala logaritmica per le ampiezze e carta semilogaritmica<br />

per le fasi;<br />

b) carta semilogaritmica sia per le ampiezze sia per le fasi. In questo caso<br />

la scala delle ampiezze è graduata in decibel: B =20 log 10 A.<br />

Caso a) Caso b)<br />

• I vantaggi che si hanno impiegando <strong>un</strong>a scala logaritmica sono:<br />

i) è possibile avere <strong>un</strong>a rappresentazione dettagliata <strong>di</strong> grandezze che<br />

variano in campi notevolmente estesi;<br />

ii) i <strong>di</strong>agrammi <strong>di</strong> <strong>Bode</strong> <strong>di</strong> sistemi in cascata si ottengono come somma<br />

dei <strong>di</strong>agrammi <strong>di</strong> <strong>Bode</strong> dei singoli sottosistemi;<br />

iii) i <strong>di</strong>agrammi <strong>di</strong> <strong>Bode</strong> <strong>di</strong> <strong>un</strong>a f<strong>un</strong>zione data in forma fattorizzata si<br />

ottengono come somma dei <strong>di</strong>agrammi elementari dei singoli fattori.<br />

R. Zanasi, R. Morselli - Controlli Automatici - 2005/06 3. ANALISI ARMONICA


3.2. DIAGRAMMI DI BODE 3.2 3<br />

Conversione in db e viceversa<br />

• Il decibel è <strong>un</strong>’<strong>un</strong>ità logaritmica <strong>con</strong>venzionale che normalmente si impiega<br />

per esprimere il guadagno <strong>di</strong> amplificatori.<br />

B(db) = 20 log 10 A<br />

• Per la <strong>con</strong>versione si può utilizzare il seguente <strong>di</strong>agramma:<br />

Posto A nella forma:<br />

A = r · 10 n<br />

<strong>con</strong> 1 ≤ r < 10,<br />

il valore <strong>di</strong> A in decibel è:<br />

B = 20 n + s db<br />

dove s si ricava dal <strong>di</strong>agramma a<br />

fianco: 0 ≤ s < 20.<br />

• Alc<strong>un</strong>e <strong>con</strong>versioni <strong>di</strong> uso frequente:<br />

A B = 20 log 10 A<br />

1 0<br />

√ 2 3<br />

2 6<br />

5 14<br />

10 20<br />

20 26<br />

50 34<br />

100 40<br />

1000 60<br />

10000 80<br />

0.5 −6<br />

0.1 −20<br />

0.01 −40<br />

– Esempio 1:<br />

– Esempio 2:<br />

A = 24<br />

A = 2.4 · 10 1<br />

B 20 + 8 = 28<br />

A = 0.56<br />

A = 5.6 · 10 −1<br />

B −20 + 15 = −5<br />

– Ogni 6 db il valore <strong>di</strong> A raddoppia;<br />

– Ogni 20 db il valore <strong>di</strong> A è<br />

moltiplicato per 10;<br />

R. Zanasi, R. Morselli - Controlli Automatici - 2005/06 3. ANALISI ARMONICA


3.2. DIAGRAMMI DI BODE 3.2 4<br />

• Si <strong>con</strong>sideri <strong>un</strong>a generica f<strong>un</strong>zione <strong>di</strong> trasferimento:<br />

G(s) = K1<br />

s m + bm−1 s m−1 + . . . + b1 s + b0<br />

s h (s n−h + an−1 s n−h−1 + . . . + ah+1 s + ah)<br />

• Il fattore s h corrisponde ad <strong>un</strong> eventuale polo nell’origine avente or<strong>di</strong>ne<br />

<strong>di</strong> molteplicità h: se la f<strong>un</strong>zione <strong>di</strong> trasferimento non presenta poli<br />

nell’origine, è h=0.<br />

• Forma fattorizzata a poli e zeri:<br />

G(s) = K1<br />

• Forma fattorizzata a costanti <strong>di</strong> tempo:<br />

G(s) = K<br />

(s − z1) (s − z2) . . . (s − zm)<br />

s h (s − ph+1) (s − ph+2) . . . (s − pn)<br />

1+τ ′ 1s 1+τ ′ 2s . . .<br />

s h 1+τ1s 1+τ2s . . .<br />

<br />

1+2δ ′ 1 s<br />

ω ′ +<br />

n1<br />

s2<br />

′ 2<br />

ω n1<br />

<br />

1+2δ1 s s2<br />

ωn1<br />

+<br />

ω 2 n1<br />

<br />

. . .<br />

<br />

. . .<br />

• La f<strong>un</strong>zione <strong>di</strong> risposta armonica si ottiene ponendo s=jω:<br />

<br />

′ ′ 1+jωτ 1 1+jωτ 2 . . . 1+2δ<br />

G(jω) = K<br />

′ jω<br />

1ω<br />

′ −<br />

n1<br />

ω2<br />

<br />

ω<br />

(jω) h <br />

1+jωτ1 1+jωτ2 . . .<br />

′ 2<br />

n1<br />

<br />

1+2δ1 jω ω2<br />

ωn1<br />

−<br />

ω 2 n1<br />

. . .<br />

<br />

. . .<br />

• Il <strong>di</strong>agramma <strong>di</strong> <strong>Bode</strong> della f<strong>un</strong>zione G(jω) si ottiene come somma dei<br />

<strong>di</strong>agrammi <strong>di</strong> <strong>Bode</strong> delle seguenti f<strong>un</strong>zioni elementari:<br />

K<br />

(j ω) −h<br />

1 + j ω τ ±1<br />

<br />

1 + j 2 δ ω ωn<br />

ω2 −<br />

ω 2 ±1<br />

n<br />

R. Zanasi, R. Morselli - Controlli Automatici - 2005/06 3. ANALISI ARMONICA


3.2. DIAGRAMMI DI BODE 3.2 5<br />

• Guadagno costante: G(jω) = K.<br />

- I <strong>di</strong>agrammi dei moduli e<br />

della fasi sono in<strong>di</strong>pendenti da<br />

ω, cioè sono costanti;<br />

- Se K > 0, il <strong>di</strong>agrammi delle<br />

fasi è β = 0;<br />

- Se K < 0, il <strong>di</strong>agramma<br />

delle fasi è β = −π.<br />

• Poli nell’origine: G(jω) = (j ω) −h .<br />

- I <strong>di</strong>agrammi riportati a fianco<br />

sono relativi al caso h = 2;<br />

- In scala (ln, ln), il <strong>di</strong>agramma<br />

delle ampiezze è <strong>un</strong>a<br />

retta passante per il p<strong>un</strong>to<br />

(ω, α) = (1, 0 db) e <strong>di</strong> inclinazione<br />

−h;<br />

- Il <strong>di</strong>agramma delle fasi è costante,<br />

identicamente uguale<br />

a −h π<br />

2 ;<br />

• Diagrammi <strong>di</strong> <strong>Bode</strong> delle ampiezze e delle fasi:<br />

ln G(jω) = α + j β = ln 1 π<br />

− j h<br />

ωh 2<br />

= −h ln ω − j h π<br />

2<br />

R. Zanasi, R. Morselli - Controlli Automatici - 2005/06 3. ANALISI ARMONICA


3.2. DIAGRAMMI DI BODE 3.2 6<br />

• Poli reali: G(jω) = 1 + j ω τ −1 =<br />

ln G(jω) = α + j β = ln<br />

• Diagrammi asintotici <strong>di</strong> <strong>Bode</strong>:<br />

Diagramma delle ampiezze:<br />

a) per ω ≪ 1/τ, si ottiene<br />

α0, cioè il <strong>di</strong>agramma<br />

tende a coincidere <strong>con</strong> l’asse<br />

delle ascisse.<br />

b) per ω ≫1/τ, si ottiene<br />

α ln 1<br />

ωτ<br />

= ln 1<br />

τ<br />

− ln ω<br />

cioè il <strong>di</strong>agramma tende alla<br />

retta passante per il p<strong>un</strong>to<br />

ln ω = ln (1/τ) e <strong>di</strong><br />

inclinazione −1;<br />

1<br />

1 + jωτ<br />

1<br />

√ 1 + ω 2 τ 2<br />

+ j (−arctan ωτ)<br />

• L’approssimazione asintotica del <strong>di</strong>agramma delle ampiezze è pertanto<br />

costituita dalle due semirette<br />

<br />

0 per ln ω ≤ ln 1<br />

τ ,<br />

α =<br />

ln 1 τ − ln ω per ln ω ≥ ln 1 τ ,<br />

• Il massimo errore della rappresentazione asintotica si ha per ω =1/τ e vale<br />

ln √ 2 (circa 3db).<br />

• La pendenza −1, su carta semilogaritmica <strong>di</strong>venta −20 db/decade;<br />

R. Zanasi, R. Morselli - Controlli Automatici - 2005/06 3. ANALISI ARMONICA


3.2. DIAGRAMMI DI BODE 3.2 7<br />

• Diagramma asintotico delle fasi:<br />

La spezzata si ottiene<br />

collegando i due<br />

asintoti β = 0 e β =<br />

−π/2 <strong>con</strong> la tangente<br />

al <strong>di</strong>agramma effettivo<br />

nel p<strong>un</strong>to corrispondente<br />

alla pulsazione<br />

<strong>di</strong> rottura<br />

ω0 = 1/τ, p<strong>un</strong>to in<br />

cui è β =−π/4.<br />

• Essendo β = −arctan ωτ , si può scrivere<br />

<br />

d β <br />

<br />

d ln ω<br />

ω=ω0<br />

= d β<br />

<br />

d ω <br />

<br />

d ω d ln ω<br />

ω=ω0<br />

= − ω0 τ<br />

1 + ω2 = −1<br />

0 τ 2 2<br />

• Le pulsazioni ωa e ωb si determinano utilizzando la relazione<br />

da cui si ottiene<br />

ln ω0<br />

ωa<br />

= ln ωb<br />

π/4<br />

ln ω0 − ln ωa<br />

ω0<br />

= π<br />

2<br />

=<br />

→<br />

π/4<br />

ln ωb − ln ω0<br />

ω0<br />

ωa<br />

= ωb<br />

ω0<br />

= 1<br />

2<br />

= e π 2 = 4, 81<br />

• I <strong>di</strong>agrammi <strong>di</strong> <strong>Bode</strong> della f<strong>un</strong>zione G(jω)=1+jωτ si ottengono ribaltando<br />

attorno all’asse delle ascisse quelli della f<strong>un</strong>zione G(jω)=(1+jωτ) −1 .<br />

• Quando la costante <strong>di</strong> tempo τ è negativa, il <strong>di</strong>agramma delle ampiezze<br />

risulta immutato, mentre il <strong>di</strong>agramma delle fasi risulta ribaltato rispetto<br />

all’asse delle ascisse.<br />

R. Zanasi, R. Morselli - Controlli Automatici - 2005/06 3. ANALISI ARMONICA


3.2. DIAGRAMMI DI BODE 3.2 8<br />

• Esempio. Si <strong>con</strong>sideri la seguente f<strong>un</strong>zione <strong>di</strong> risposta armonica:<br />

• Diagrammi <strong>di</strong> <strong>Bode</strong>:<br />

G(jω) =<br />

5, 6 (1 + j ω 0, 5)<br />

(1 + j ω 4) (1 + j ω 0, 25) (1 + j ω 0, 125)<br />

R. Zanasi, R. Morselli - Controlli Automatici - 2005/06 3. ANALISI ARMONICA


3.2. DIAGRAMMI DI BODE 3.2 9<br />

• Poli complessi <strong>con</strong>iugati (0≤δ


3.2. DIAGRAMMI DI BODE 3.2 10<br />

• Il <strong>di</strong>agramma effettivo si può <strong>di</strong>scostare sensibilmente da quello asintotico:<br />

in particolare, per δ =0 e in corrispondenza della pulsazione <strong>di</strong> rottura ωn,<br />

lo scostamento è infinito.<br />

• Il <strong>di</strong>agramma delle ampiezze ha le seguenti proprietà:<br />

1) per 0≤δ ≤1/ √ 2, presenta <strong>un</strong> massimo;<br />

2) per 0≤δ ≤1/2, interseca l’asse a destra del p<strong>un</strong>to ω =ωn;<br />

3) per 1/2≤δ ≤1/ √ 2, interseca l’asse a sinistra del p<strong>un</strong>to ω =ωn;<br />

4) per 1/ √ 2≤δ ≤1, non interseca l’asse delle ascisse ed è pertanto tutta<br />

al <strong>di</strong> sotto della sua approssimazione asintotica.<br />

• Pulsazione <strong>di</strong> risonanza ωR. Posto u=ω/ωn, il massimo dell’ampiezza<br />

corrisponde ad <strong>un</strong> minimo della f<strong>un</strong>zione<br />

(1 − u 2 ) 2 + 4 δ 2 u 2<br />

Derivando e uguagliando a zero la derivata, si ottiene<br />

−4 (1 − u 2 ) u + 8 δ 2 u = 0<br />

Trascurando la soluzione nulla si ottiene<br />

uR = 1−2 δ2 →<br />

√<br />

ωR = ωn 1 − 2 δ2 • Picco <strong>di</strong> risonanza MR: si calcola come modulo della f<strong>un</strong>zione <strong>di</strong> risposta<br />

armonica in corrispondenza della pulsazione ωR:<br />

MR =<br />

1<br />

(1 − 1 − 2 δ 2 ) 2 + 4 δ 2 (1 − 2 δ 2 )<br />

MR =<br />

1<br />

2 δ √ 1 − δ 2<br />

R. Zanasi, R. Morselli - Controlli Automatici - 2005/06 3. ANALISI ARMONICA


3.2. DIAGRAMMI DI BODE 3.2 11<br />

• Diagramma <strong>di</strong> <strong>Bode</strong> della fasi:<br />

• Approssimazione asintotica del <strong>di</strong>agramma delle fasi si ottiene <strong>con</strong>gi<strong>un</strong>gendo<br />

gli asintoti β = 0 e β = −π <strong>con</strong> <strong>un</strong> segmento inclinato <strong>di</strong> pendenza<br />

opport<strong>un</strong>a.<br />

• Essendo β = −arctan<br />

<br />

d β <br />

<br />

d ln ω<br />

=<br />

ω=ωn<br />

d β<br />

d u<br />

2 δ u<br />

1−u 2 , si deduce:<br />

<br />

d u <br />

<br />

1<br />

d ln ω = − <br />

u=1 1 + 2 δ u<br />

1 − u 2<br />

2 2 δ (1 + u 2 ) u<br />

(1 − u 2 ) 2<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

u=1<br />

= − 1<br />

δ<br />

• Le pulsazioni ωa e ωb sono legate alla pulsazione ωn dalla relazione:<br />

π/2<br />

ln ωn − ln ωa<br />

=<br />

dalla quale si ottiene<br />

ωn<br />

ωa<br />

= ωb<br />

ωn<br />

π/2<br />

ln ωb − ln ωn<br />

= e π 2 δ = 4, 81 δ<br />

= 1<br />

δ<br />

R. Zanasi, R. Morselli - Controlli Automatici - 2005/06 3. ANALISI ARMONICA


3.2. DIAGRAMMI DI BODE 3.2 12<br />

• Diagrammi delle ampiezze e delle fasi in scala semilogaritmica:<br />

R. Zanasi, R. Morselli - Controlli Automatici - 2005/06 3. ANALISI ARMONICA


3.2. DIAGRAMMI DI BODE 3.2 13<br />

Graficazione “qualitativa” dei <strong>di</strong>agrammi <strong>di</strong> <strong>Bode</strong><br />

Si procede alla graficazione del <strong>di</strong>agramma asintotico <strong>di</strong> <strong>Bode</strong> della f<strong>un</strong>zione<br />

<strong>di</strong> risposta armonica assegnata. Due possibili meto<strong>di</strong>.<br />

Primo metodo: somma dei singoli <strong>con</strong>tributi<br />

a) La f<strong>un</strong>zione G(s) viene messa nella forma “a costanti <strong>di</strong> tempo”:<br />

G(s) =<br />

10(s − 1)<br />

s(s + 1)(s 2 + 8s + 25)<br />

→ G(s) = −10<br />

25<br />

(1 − s)<br />

s(1 + s)(1 + 8s<br />

25<br />

b) Si tracciano i <strong>di</strong>agrammi asintotici <strong>di</strong> <strong>Bode</strong> delle singole componenti:<br />

K = − 10<br />

25 , G1(s) = (1 − s), G2(s) = 1<br />

s , G3(s) =<br />

1<br />

(1 + s) , G4(s) =<br />

1<br />

(1 + 8s s2 + 25 25 )<br />

+ s2<br />

25 )<br />

c) Si sommano i singoli <strong>con</strong>tributi per ottenere il <strong>di</strong>agramma asintotico della<br />

f<strong>un</strong>zione G(s).<br />

- Il <strong>con</strong>tributo del termine K è costante: |K| = −7.96 db e arg K = −π.<br />

- Lo zero instabile (1 − s), e il polo stabile (1 + s) −1 agis<strong>con</strong>o alla pulsazione<br />

ω = 1 e fornis<strong>con</strong>o due <strong>con</strong>tribuiti uguali e <strong>con</strong>trari nel <strong>di</strong>agramma delle<br />

ampiezze. Il loro <strong>con</strong>tributo nel <strong>di</strong>agramma delle fasi si somma: l’ampiezza<br />

complessiva per ω → ∞ è −π.<br />

- La coppia <strong>di</strong> poli complessi <strong>con</strong>iugati (1+ 8s<br />

25<br />

+ s2<br />

25 )−1 determina sul <strong>di</strong>agramma<br />

asintotico delle ampiezze <strong>un</strong>a attenuazione <strong>di</strong> −40 db/dec a partire dalla pulsazione<br />

ωn = 5. Il <strong>con</strong>tributo al <strong>di</strong>agramma delle fasi è negativo <strong>di</strong> ampiezza<br />

complessiva −π al variare <strong>di</strong> ω. Le pulsazioni alle quali si ha <strong>un</strong> cambiamento<br />

<strong>di</strong> pendenza del <strong>di</strong>agramma asintotico <strong>di</strong> <strong>Bode</strong> delle fasi sono le seguenti<br />

ωa = 1<br />

4.81 , ωb = 4.81, ¯ωa = ωn<br />

4.81 δ , ¯ωb = ωn4.81 δ<br />

dove δ = 0.8 è il coefficiente <strong>di</strong> smorzamento della coppia <strong>di</strong> poli complessi<br />

<strong>con</strong>iugati.<br />

- La <strong>di</strong>fficoltà nell’utilizzare questo metodo sta nel fatto che la somma dei<br />

singoli <strong>con</strong>tributi non è sempre agevole.<br />

R. Zanasi, R. Morselli - Controlli Automatici - 2005/06 3. ANALISI ARMONICA


3.2. DIAGRAMMI DI BODE 3.2 14<br />

Diagrammi asintotici <strong>di</strong> <strong>Bode</strong> delle ampiezze e delle fasi della f<strong>un</strong>zione G(s)<br />

α = α |G(jω)|(db) = |G(jω)|<br />

40<br />

|G2|<br />

20<br />

0<br />

|K|<br />

-20<br />

-40<br />

-60<br />

-80<br />

ϕ0 = − 3π<br />

2<br />

−2π<br />

− 5π<br />

2<br />

−3π<br />

− 7π<br />

2<br />

arg |G(jω)|<br />

π<br />

π<br />

2<br />

0<br />

− π<br />

2<br />

−π<br />

−1<br />

0.1 1 5 10 100<br />

β<br />

B<br />

0.1 ωa 1 ¯ωa 5 10 100<br />

ωa<br />

G1, G3<br />

G4<br />

¯ωa<br />

ωb<br />

¯ω<br />

ωb<br />

A<br />

−2<br />

−3<br />

G4<br />

¯ωb<br />

¯ωb<br />

+1<br />

C<br />

|G1|<br />

ω [log 10]<br />

|G3|<br />

|G4|<br />

ω [log 10]<br />

G2<br />

K<br />

K, G2<br />

G1, G3<br />

ϕ∞ = − 7π<br />

2<br />

R. Zanasi, R. Morselli - Controlli Automatici - 2005/06 3. ANALISI ARMONICA


3.2. DIAGRAMMI DI BODE 3.2 15<br />

Se<strong>con</strong>do metodo: graficazione “rapida”<br />

Diagramma delle ampiezze<br />

a) Si in<strong>di</strong>vidua nella f<strong>un</strong>zione G(s) tutte le pulsazioni in corrispondenza delle<br />

quali si ha <strong>un</strong> cambiamento <strong>di</strong> pendenza. Tali pulsazioni coincidono <strong>con</strong><br />

gli zeri reali, i poli reali e <strong>con</strong> le pulsazioni naturali ωn delle coppie <strong>di</strong> poli<br />

e zeri complessi <strong>con</strong>iugati della f<strong>un</strong>zione G(s). Nel caso in esame si ha<br />

ω = 1 e ω = 5. Tali pulsazioni vengono or<strong>di</strong>nate in or<strong>di</strong>ne crescente <strong>di</strong><br />

modulo.<br />

b) Tenendo <strong>con</strong>to del fatto che gli zeri (reali o complessi <strong>con</strong>iugati) determinano<br />

<strong>un</strong> incremento <strong>di</strong> pendenza (rispettivamente <strong>di</strong> +1 e <strong>di</strong> +2) e che,<br />

viceversa, i poli (reali o complessi <strong>con</strong>iugati) determinano <strong>un</strong> decremento<br />

della pendenza del <strong>di</strong>agramma asintotico (rispettivamente <strong>di</strong> -1 e <strong>di</strong> -2), è<br />

chiaro che la “forma” del <strong>di</strong>agramma asintotico è già nota a priori prima<br />

<strong>di</strong> iniziare la graficazione.<br />

Nel caso in esame, per esempio, si ha:<br />

-1<br />

β<br />

(o,x)<br />

<br />

x<br />

x<br />

1 5<br />

In corrispondenza della pulsazione ω = 1 non si ha cambiamento si pendenza perchè<br />

a questa pulsazione agis<strong>con</strong>o <strong>con</strong>temporaneamente sia <strong>un</strong> polo che <strong>un</strong>o zero.<br />

c) Si determina la posizione “verticale” del <strong>di</strong>agramma asintotico <strong>di</strong> <strong>Bode</strong>.<br />

- Se la f<strong>un</strong>zione G(s) è <strong>di</strong> tipo 0, il posizionamento verticale è automaticamente<br />

determinato dal calcolo del guadagno statico G(0).<br />

- Se il sistema è <strong>di</strong> tipo 1, o in generale <strong>di</strong> tipo h, il posizionamento verticale<br />

può avvenire, per esempio, calcolando la posizione del p<strong>un</strong>to A in<br />

corrispondenza della pulsazione ¯ω alla quale si ha il primo cambiamento<br />

<strong>di</strong> pendenza. Siano (¯ω, β) le coor<strong>di</strong>nate del p<strong>un</strong>to A. Se il sistema è <strong>di</strong><br />

R. Zanasi, R. Morselli - Controlli Automatici - 2005/06 3. ANALISI ARMONICA<br />

-1<br />

¯ω<br />

A<br />

-3<br />

ω


3.2. DIAGRAMMI DI BODE 3.2 16<br />

d)<br />

tipo h, il valore della coor<strong>di</strong>nata β si determina in base alla formula:<br />

<br />

<br />

β = <br />

s<br />

<br />

hG(s) ¯ω h<br />

<br />

<br />

<br />

<br />

s=0<br />

cioè si sostituisce ¯ω h al posto degli h poli nell’origine, mentre in tutti gli<br />

altri termini <strong>di</strong> G(s) si mette s = 0.<br />

Nel caso in esame si ha h = 1 ed ¯ω = 5 per cui<br />

<br />

<br />

β = <br />

s G(s) <br />

<br />

5 = <br />

10(s − 1)<br />

<br />

s=0 5(s + 1)(s2 <br />

<br />

<br />

+ 8s + 25) <br />

s=0<br />

= 2<br />

25<br />

= −21.94 db<br />

È ora possibile tracciare il <strong>di</strong>agramma asintotico complessivo tracciando, a<br />

partire da A, i vari tratti della “spezzata”, ogn<strong>un</strong>o <strong>con</strong> la propria pendenza.<br />

Nel caso in esame, per esempio, il tratto <strong>di</strong> spezzata che precede il p<strong>un</strong>to A si<br />

determina in<strong>di</strong>viduando il p<strong>un</strong>to B. Questo p<strong>un</strong>to si calcola a partire da A <strong>di</strong>minuendo<br />

la pulsazione <strong>di</strong> <strong>un</strong>a decade ed aumentando <strong>di</strong> 20 db l’ampiezza: B = (0.5, β + 20).<br />

Allo stesso modo si procede per determinare il tratto che segue il p<strong>un</strong>to A. In questo<br />

caso, essendo -3 la pendenza <strong>di</strong> questo tratto, il p<strong>un</strong>to C si determina aumentando<br />

la pulsazione <strong>di</strong> <strong>un</strong>a decade e <strong>di</strong>minuendo l’ampiezza <strong>di</strong> 60 db: C = (50, β − 60).<br />

Diagramma delle fasi<br />

Anche la graficazione del <strong>di</strong>agramma asintotico delle fasi può essere fatta<br />

più “rapidamente” se si procede nel modo seguente (si faccia riferimento al<br />

<strong>di</strong>agramma delle fasi riportato nella precedente figura):<br />

a) Si in<strong>di</strong>vidua la fase <strong>di</strong> partenza ϕ0 del <strong>di</strong>agramma asintotico delle fasi<br />

calcolando la fase iniziale della f<strong>un</strong>zione approssimante G0(s) per ω → 0.<br />

. Tale fase è<br />

Nel caso in esame, per esempio, la fase iniziale è ϕ0 = −3π 2<br />

comprensiva del segno negativo della costante K e della fase costante − π<br />

2<br />

introdotta dal polo nell’origine.<br />

b) Si prendono in <strong>con</strong>siderazione i poli e gli zeri in or<strong>di</strong>ne crescente della loro<br />

pulsazione critica (cioè il valore assoluto del corrispondente polo o zero,<br />

oppure la ωn nel caso <strong>di</strong> coppie <strong>di</strong> poli o zeri complessi <strong>con</strong>iugati). Il<br />

<strong>di</strong>agramma asintotico <strong>di</strong> ciasc<strong>un</strong> elemento viene <strong>di</strong>segnato in <strong>un</strong>a <strong>di</strong>versa<br />

R. Zanasi, R. Morselli - Controlli Automatici - 2005/06 3. ANALISI ARMONICA


3.2. DIAGRAMMI DI BODE 3.2 17<br />

“fascia” in f<strong>un</strong>zione dell’azione introdotta dai precedenti elementi, e in<br />

f<strong>un</strong>zione del fatto che l’elemento sia stabile o instabile.<br />

Nel caso in esame, per esempio, i primi due elementi da prendere in <strong>con</strong>siderazione<br />

all’aumentare <strong>di</strong> ω sono il polo stabile (s+1) −1 e lo zero instabile<br />

(s − 1). Questi due elementi agis<strong>con</strong>o <strong>con</strong>temporaneamente e ciasc<strong>un</strong>o<br />

<strong>di</strong> essi introduce <strong>un</strong>o “sfasamento”, per ω ∈ [−∞, ∞], pari a − π<br />

2 . Il<br />

<strong>con</strong>tributo complessivo <strong>di</strong> questi due elementi è <strong>un</strong> <strong>di</strong>agramma asintotico<br />

<strong>di</strong> ampiezza −π da <strong>di</strong>segnare verso il basso nella fascia − 3π<br />

2 , −5π<br />

2 . Anche<br />

la coppia <strong>di</strong> poli complessi <strong>con</strong>iugati (1 + 8s s2<br />

25 + 25 )−1 , essendo stabili,<br />

introduce <strong>un</strong>o sfasamento <strong>di</strong> ampiezza complessiva −π. Il loro <strong>con</strong>tributo<br />

al <strong>di</strong>agramma asintotico va quin<strong>di</strong> <strong>di</strong>segnato nella fascia − 5π,<br />

−7π<br />

2 2 .<br />

c) Si procede alla graficazione del <strong>di</strong>agramma asintotico complessivo “interpolando”<br />

i <strong>di</strong>agrammi asintotici delle fasi dei singoli elementi, ogn<strong>un</strong>o dei<br />

quali è stato <strong>di</strong>segnato nella fascia più opport<strong>un</strong>a. Nel caso in esame,<br />

è evidente che la determinazione dei p<strong>un</strong>ti interme<strong>di</strong> <strong>di</strong> cambiamento <strong>di</strong><br />

pendenza risulta notevolmente semplificato rispetto al caso <strong>di</strong> semplice<br />

somma dei singoli <strong>con</strong>tributi asintotici. Si noti inoltre che la fase finale<br />

ϕ∞ ottenuta è in accordo <strong>con</strong> la fase finale della f<strong>un</strong>zione approssimante<br />

G∞(s) per ω → ∞.<br />

R. Zanasi, R. Morselli - Controlli Automatici - 2005/06 3. ANALISI ARMONICA


3.4. LA FORMULA DI BODE 3.4 1<br />

La formula <strong>di</strong> <strong>Bode</strong><br />

• Una f<strong>un</strong>zione <strong>di</strong> trasferimento razionale fratta è a a fase minima se non<br />

ha né poli né zeri nel semipiano destro del piano s.<br />

• Per sistemi a fase minima, detta ωc la pulsazione in corrispondenza della<br />

quale si vuole calcolare la fase βc, vale la formula <strong>di</strong> <strong>Bode</strong>:<br />

βc = 1<br />

∞<br />

d α<br />

<br />

<br />

ln cotanh <br />

π d u u<br />

<br />

<br />

du<br />

2<br />

−∞<br />

in cui si è posto α := ln |G(jω)|, u := ln ω<br />

ωc<br />

= ln ω−ln ωc.<br />

• La fase βc in corrispondenza <strong>di</strong> <strong>un</strong>a<br />

data pulsazione ωc <strong>di</strong>pende essen-<br />

d α<br />

zialmente dalla pendenza d u del <strong>di</strong>agramma<br />

delle ampiezze nell’intorno<br />

<strong>di</strong> quella pulsazione ωc.<br />

• Esempio: βc = 0·A1 − 1·A2 − 2·A3<br />

dove:<br />

A1 = π<br />

4 −β1, A2 = β1+β2, A3 = π<br />

4 −β2<br />

R. Zanasi, R. Morselli - Controlli Automatici - 2005/06 3. ANALISI ARMONICA


3.4. LA FORMULA DI BODE 3.4 2<br />

• Significato della variabile <strong>di</strong> integrazione u: se il <strong>di</strong>agramma α è riferito ai<br />

logaritmi naturali, la variabile u non è altro che l’ascissa ln ω <strong>con</strong> l’origine<br />

traslata in ln ωc.<br />

• La <strong>con</strong><strong>di</strong>zione necessaria e sufficiente per la vali<strong>di</strong>tà della formula <strong>di</strong> <strong>Bode</strong>,<br />

cioè che la f<strong>un</strong>zione <strong>di</strong> trasferimento sia a fase minima, è sod<strong>di</strong>sfatta per<br />

la quasi totalità dei sistemi che normalmente si <strong>con</strong>siderano.<br />

• Esempio <strong>di</strong> rete elettrica a fase non minima:<br />

•<br />

• F<strong>un</strong>zione <strong>di</strong> trasferimento<br />

Vu(s) =<br />

1/Cs − R<br />

R + 1/Cs Vi(s)<br />

che, posto T :=RC, <strong>di</strong>venta<br />

G(s) = Vu(s)<br />

Vi(s)<br />

= 1 − T s<br />

1 + T s ,<br />

cioè <strong>un</strong>a f<strong>un</strong>zione non a fase<br />

non minima;<br />

• Diagrammi <strong>di</strong> <strong>Bode</strong> delle<br />

ampiezze e delle fasi: il <strong>di</strong>agramma<br />

delle ampiezze è costante<br />

α = 0 (|G(jω)| = 1),<br />

mentre il <strong>di</strong>agramma delle fasi<br />

varia gradualmente da 0 ◦ a<br />

−180 ◦ .<br />

È chiaro che applicando la formula <strong>di</strong> <strong>Bode</strong> all’esempio si sarebbe invece<br />

dedotta <strong>un</strong>a fase identicamente nulla.<br />

R. Zanasi, R. Morselli - Controlli Automatici - 2005/06 3. ANALISI ARMONICA


3.4. LA FORMULA DI BODE 3.4 3<br />

• La f<strong>un</strong>zione <strong>di</strong> trasferimento trascendente<br />

G(s) = e −t0 s<br />

che rappresenta <strong>un</strong> ritardo finito <strong>di</strong> valore t0, non è a fase minima.<br />

• Essendo<br />

G(jω) = e −jωt0 = cos ω t0 − j sen ω t0 ,<br />

la f<strong>un</strong>zione <strong>di</strong> risposta armonica ha modulo identicamente <strong>un</strong>itario e fase<br />

crescente linearmente <strong>con</strong> la frequenza.<br />

• Per ricavare i <strong>di</strong>agrammi <strong>di</strong> <strong>Bode</strong>, si scrive<br />

ln ω<br />

ln G(jω) = α + j β = 0 − j ω t0 = 0 − j t0 e<br />

relazione dalla quale si deduce che il <strong>di</strong>agramma delle fasi ha <strong>un</strong> andamento<br />

esponenziale. Anche in questo caso l’applicazione della formula <strong>di</strong> <strong>Bode</strong><br />

avrebbe <strong>con</strong>dotto ad <strong>un</strong> risultato errato (β = 0).<br />

R. Zanasi, R. Morselli - Controlli Automatici - 2005/06 3. ANALISI ARMONICA


3.4. LA FORMULA DI BODE 3.4 4<br />

Diagrammi <strong>di</strong> <strong>Bode</strong>: riepilogo<br />

• I <strong>di</strong>agrammi <strong>di</strong> <strong>Bode</strong> si basano su alc<strong>un</strong>e proprietà dei logaritmi e sulle<br />

proprietà valide per il modulo e l’argomento <strong>di</strong> f<strong>un</strong>zioni complesse:<br />

1. |A B| = |A| |B| ⇒ log10 |A B| = log10 |A| + log10 |B|<br />

<br />

<br />

2. <br />

A<br />

<br />

|A|<br />

B<br />

=<br />

|B| ⇒ log <br />

<br />

<br />

A<br />

<br />

10 B<br />

= log10 |A| − log10 |B|<br />

3. arg(A B) = arg(A) + arg(B)<br />

<br />

A<br />

4. arg = arg(A) − arg(B)<br />

B<br />

• Queste proprietà permettono <strong>di</strong> ri<strong>con</strong>durre il problema <strong>di</strong> tracciare i <strong>di</strong>agrammi<br />

<strong>di</strong> bode <strong>di</strong> <strong>un</strong>a generica f<strong>un</strong>zione razionale fratta<br />

•<br />

G(jω) = K (jω + b)...(−ω2 + 2δbωbjω + ω 2 b )<br />

(jω) h (jω + d)...(−ω 2 + 2δdωdjω + ω 2 d )<br />

alla “somma” dei <strong>di</strong>agrammi <strong>di</strong> <strong>Bode</strong> delle seguenti f<strong>un</strong>zioni elementari:<br />

K<br />

(jω) −h<br />

(jω + a) ±1<br />

(−ω 2 + 2δaωajω + ω 2 a) ±1<br />

È possibile ricavare i <strong>di</strong>agrammi asintotici <strong>di</strong> <strong>Bode</strong> anche senza ricorrere<br />

alla somma dei singoli <strong>con</strong>tributi, ma utilizzando <strong>un</strong> metodo <strong>di</strong>retto.<br />

R. Zanasi, R. Morselli - Controlli Automatici - 2005/06 3. ANALISI ARMONICA


3.4. LA FORMULA DI BODE 3.4 5<br />

Assi nei <strong>di</strong>agrammi <strong>di</strong> <strong>Bode</strong><br />

• I Diagrammi <strong>di</strong> <strong>Bode</strong> usano l’asse orizzontale in scala logaritmica. Considerando<br />

l’asse reale R e fissata <strong>un</strong>a origine, <strong>un</strong>a pulsazione ω corrisponde<br />

ad <strong>un</strong> p<strong>un</strong>to sull’asse <strong>con</strong> coor<strong>di</strong>nata x = log 10 ω. Accanto all’asse si<br />

possono quin<strong>di</strong> in<strong>di</strong>care o i valori della coor<strong>di</strong>nata x oppure <strong>di</strong>rettamente<br />

i valori <strong>di</strong> ω; questa se<strong>con</strong>da soluzione è la più comoda.<br />

- 1<br />

0.1<br />

0 1 2<br />

1 10 100<br />

0 0.3 0.5 0.7 0.9 1 x= log <br />

1 2 3 5 8 10<br />

x= log <br />

• Per il <strong>di</strong>segno qualitativo dei <strong>di</strong>agrammi <strong>con</strong>viene memorizzare alc<strong>un</strong>i valori:<br />

log 10 2 0.3 log 10 3 0.5 log 10 5 0.7 log 10 8 0.9<br />

• L’asse verticale nei <strong>di</strong>agrammi <strong>di</strong> ampiezza è graduato in decibel (db):<br />

A| db<br />

def<br />

= 20 log10 A<br />

Con questa scala, le pendenze caratteristiche dei <strong>di</strong>agrammi <strong>di</strong> <strong>Bode</strong> sono<br />

±20 db/decade, ±40 db/decade, ecc. Per como<strong>di</strong>tà tali pendenze<br />

vengono in<strong>di</strong>cate rispettivamente <strong>con</strong> i SIMBOLI ±1, ±2, ecc.<br />

N.B.: se la scala verticale fosse semplicemente logaritmica (y = log 10 A)<br />

le pendenze caratteristiche dei <strong>di</strong>agrammi <strong>di</strong> <strong>Bode</strong> sarebbero ±1 e ±2.<br />

• L’asse verticale nei <strong>di</strong>agrammi <strong>di</strong> fase può essere graduato sia in ra<strong>di</strong>anti sia<br />

in gra<strong>di</strong>. In ogni caso il <strong>di</strong>agramma delle fasi può essere traslato verso l’alto<br />

o verso il basso <strong>di</strong> multipli interi <strong>di</strong> 2π, o <strong>di</strong> 360 o , mantenendo inalterato<br />

il suo significato.<br />

R. Zanasi, R. Morselli - Controlli Automatici - 2005/06 3. ANALISI ARMONICA

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!