Materie: | Appunti |
Categoria: | Elettronica |
Download: | 363 |
Data: | 20.02.2006 |
Numero di pagine: | 12 |
Formato di file: | .doc (Microsoft Word) |
Download
Anteprima
elettronica-place-diagrammi-bode_1.zip (Dimensione: 97.22 Kb)
trucheck.it_elettronica:-la-place-e-i-diagrammi-di-bode.doc 258.5 Kb
readme.txt 59 Bytes
Testo
TRASFORMATA DI LAPLACE
E
DIAGRAMMI DI BODE
a. s. 2005 / 2006
PREMESSA
Per lo studio dei sistemi di controllo si utilizzano modelli matematici dinamici lineari. L’analisi o il progetto di un sistema di controllo comporta, in generale, la risoluzione di equazioni differenziali ordinarie a coefficienti costanti.
Lo schema generale, riportato in figura, mostra un percorso alternativo alla risoluzione di equazioni differenziali, che consente ugualmente di giungere alla soluzione del problema nel campo reale.
La trasformazione dal campo reale al campo complesso avviene tramite la TRASFORMATA DI LA PLACE ( L ) che associa ad una generica funzione del tempo f ( t ) una funzione complessa F ( s ) di variabile complessa s.
CONCETTO DI TRASFORMATA
Ogni metodo trasformazionale permette di semplificare un’operazione matematica
( per esempio: qualora risulti difficile fare il prodotto di due numeri, conviene effettuare la somma dei loro logaritmi ).
TRASFORMATA DI LA PLACE
Ogni trasformata ha il suo specifico campo di utilizzazione; la L - trasformata permette di risolvere in modo semplificato le equazioni integro – differenziali poiché le trasforma in equazioni algebriche.
La trasformata di Laplace è la trasformata di una funzione di variabile reale come ad esempio il tempo t . Tale funzione deve soddisfare alle condizioni di essere nulla nella regione negativa del tempo cioè: f ( t ) = 0 per t < 0 e definita nella regione positiva del tempo.
Si definisce trasformata di Laplace ( L ) di una funzione f ( t ), la funzione della variabile complessa S, definita per mezzo della relazione integrale:
in cui :
• f ( t ) = funzione del tempo, nulla per t < 0, definita funzione originale;
• s = variabile complessa, definita operatore di Laplace, espressa come: s = i + j ;
• F ( s ) = funzione della variabile s risultante dall'operazione di trasformazione
(L-trasformata).
Nella ricerca della funzione di trasferimento di un quadripolo contenente elementi reattivi, si giunge ad una equazione differenziale.
Applicando la L-trasformata alla precedente equazione differenziale si ottiene una equazione algebrica F ( s ), nella quale la variabile non è più t ma s.
In quanto equazione algebrica, la F ( s ) può essere trattata con le regole dell'algebra, con l'obiettivo ad esempio di ottenere il rapporto tra il segnale di uscita e quello di ingresso del quadripolo.
Determinata la funzione di trasferimento in funzione della variabile s, è possibile:
• antitrasformare, ovvero ottenere di nuovo una funzione della variabile t;
• effettuare considerazioni che, pur sviluppate analizzando la funzione in s, consentono di trarre conclusioni in merito al comportamento generale del quadripolo nel dominio del tempo.
PROPRIETÀ DELLA L - TRASFORMATA
a) Teorema della linearità
La trasformata di una combinazione lineare di funzione del tempo è uguale alla combinazione lineare delle L - trasformate delle funzioni.
In formule:
se L [ f ( t ) ] = F ( s ) e L [ g ( t ) ] = G ( s )
ed A e B sono costanti, risulta verificata l’ eguaglianza :
L [ A * f ( t ) + B * g ( t ) ] = A * F ( s ) + B * G ( s )
b) Teorema della derivata
Se L [ f ( t ) ] = F ( s ) risulta:
L [ d f ( t ) / d t ] = S * F (s) – f ( 0+ )
Nella precedente espressione , f ( 0+) rappresenta un termine che tiene conto delle condizioni iniziali, ovvero il valore posseduto da f ( t ) per t = 0.
Nel caso di derivate multiple si ha:
L [ d n f (t) / d t n ] = S n F (s) – S n - 1 * f ( 0+) – S n - 2 * f ’ ( 0+) - …..
ESEMPIO: Per la tensione ai capi di un componente induttivo si è ricavata la relazione:
VL ( t ) = L * d i ( t ) / d t
Applicando il teorema della derivata si ottiene:
VL (s) = S * L* I (s) – L * I0
Avendo indicato con I0 la corrente eventualmente circolante nell‘ induttanza all’istante iniziale.Il termine L* I0 che esprime l’energia immagazzinata inizialmente nell’ induttanza , ha le dimensioni di une tensione e può essere rappresentata con un generatore nello schema equivalente del componente L – trasformato.
c) Teorema dell’ integrale
Se L [ f ( t ) ] = F ( s ) , risulta:
L [ f ( t ) d t ] = [ F ( s ) / s ] + [ f ( 0+) / s] con f ( 0+) = ooo+ f (t) d t
Nella precedente espressione , f ( 0+) rappresenta un termine che tiene conto delle condizioni iniziali , ovvero del valore posseduto da f ( t) per t = 0.
ESEMPIO: Per la tensione ai capi di un condensatore si è ricavata la relazione :
VC ( t ) = ( 1 / C ) * i (t) d t
Applicando il teorema dell’ integrale si ottiene:
VC ( s ) = (1 / sC ) * I ( s ) + ( V0 / s )
Avendo indicato con V0 la differenza di potenziale presente ai capi del condensatore qualora sia inizialmente scarico. Il termine V0 / s può essere rappresentato con un generatore nello schema equivalente del componente L -trasformato.
d) Teorema del valore iniziale
Se L [ f ( t ) ] = F ( s ) risulta:
lim t -> 0+ f ( t ) = lim s -> > s * F (s)
Questo teorema consente di determinare il valore assunto dalla f ( t ) nell’ intorno (positivo) dell’ origine, conoscendo la L – trasformata F ( s ), ovvero senza che sia necessario antitraformare.
e) Teorema del valore finale
Se L [ f ( t ) ] = F ( s ) risulta:
lim t -> > f ( t ) = lim s -> 0 s * F (s)
Questo teorema consente di determinare il valore assunto dalla f ( t ) a regime (cioè per tt ), conoscendo la L – trasformata F ( s ), ovvero senza che sia necessario antitrasformare.
CADUTA DI TENSIONE AI CAPI DEI BIPOLI PASSIVI
Espressioni in funzione del tempo Espressioni L - trasformate
Un componente reattivo può essere schematizzato con
• Z C = 1 / s C se è un condensatore
• Z L = s L se è un induttore
a meno di eventuali condizioni iniziali schematizzate con un generatore di tensione (Vo / s) nel condensatore ed un generatore di corrente ( - L Io ) nell'induttore.
ESEMPIO 1
Determinare la funzione di trasferimento del seguente circuito:
Essendo V i (t) = V r (t) + V c (t) con V r (t) = R * i (t) e V c (t) = 1/C * / i (t) d t ,
applicando le L-trasformate con l'ipotesi che il condensatore sia inizialmente scarico, si ha:
V r (s) = R * I (s) e V c (s) = 1 / (sC) * I (s) ; allora:
ESEMPIO 2
Determinare la funzione di trasferimento del seguente circuito:
Essendo V i (t) = V r (t) + V L (t) con V r (t) = R * i (t) e V L(t) = L * d i (t) / d t ,
applicando le L - trasformate e nell'ipotesi che l'induttanza non sia inizialmente percorsa da corrente, si ha: V r (s) = R * I (s) e V L(s) = s L * I (s) ; allora:
ESEMPIO 3
Determinare la funzione di trasferimento del seguente circuito:
Poichè V i (t) = V r (t) + V L (t) + V c (t) con V r (t) = R * i (t), V c (t) = 1/C * / i (t) d t
e VL(t) = L * d i (t) / d t , applicando le L - trasformate nell'ipotesi che il condensatore sia inizialmente scarico e che l'induttanza non sia inizialmente percorsa da corrente, si ha:
V r (s) = R * I (s) , V c (s) = 1 / (sC) * I (s) e V L(s) = s L * I (s) ; per cui:
L’ANTITRASFORMAZIONE
Antitrasformare significa ricavare l’espressione della funzione del tempo corrispondente ad una funzione di s.
A ) Utilizzo delle tabelle di trasformazione
Per casi semplici si utilizza direttamente le tabelle.
ESEMPIO N. 1
posto RC = p e dividendo numeratore e denominatore per i si ha:
Dalla tabella(riga 1 e 6),se si pone a = k = 1 / / , si ha che l’antitrasformata è :
Se: R=100k0 ; C=20 F
ESEMPIO N. 2
Posto L / R = / ; dividendo numeratore e denominatore per si ha :
che si antitrasforma come nel precedente caso:
Se: R=1S ; L=100mH
B) Lo sviluppo in frazioni parziali
Nei casi precedenti si è visto l’utilizzo della corrispondenza tra la funzione
F (s) = k / (s+a) e la corrispondente funzione nel dominio del tempo f (t) = K e -at .
La F(s) considerata costituisce una espressione notevole, in quanto al numeratore presenta una costante e al denominatore un binomio nella variabile s.
Alla base del metodo di antitrasformazione vi è la possibilità di scomporre un polinomio nel quale si conoscono le radici (valori della variabile che annullano il polinomio), nel prodotto di
tanti binomi quante sono le radici stesse, per trasformarlo poi in una somma di frazioni che hanno al numeratore delle costanti ed al denominatore ciascuno dei precedenti binomi.
Se A(S) = N(S) / D(S) è l’espressione della funzione di trasferimento e (P1, P2, P3…..) sono valori di s che annullano il denominatore, essa può essere scritta come segue:
A(S) = N(S) / D(S) = N(S) / (s - P1) (s - P2) (s - P3)…..
I valori di s che annullano il polinomio al numeratore si dicono ZERI.
I valori di s che annullano il polinomio al denominatore si dicono POLI.
Si sviluppano metodi che consentano l’antitrasformazione se il grado del polinomio al numeratore è minore di quello al denominatore cioè se il numero degli zeri è minore di quello dei poli.
Effettuare lo SVILUPPO IN FRAZIONI PARZIALI significa determinare i valori da assegnare alle costanti A, B, C, …. che permettono di soddisfare l'eguaglianza:
A(S) =N(S) / D(S) =N(S) / (s - P1) (s - P2) (s - P3)…= A / (s - P1) + B / (s - P2) + C / (s - P3)
Per determinare i valori da dare alle costanti A, B, C, …. è necessario che N ( s ) risulti uguale al numeratore ottenuto ricavando il denominatore comune tra le frazioni parziali ed effettuando la somma.
Se vengono determinati A, B, C, …. è possibile antitrasformare separatamente ciascuno degli addendi che assumono l'espressione e - a t = 1 / (s + a ) ottenendo l'andamento della funzione di trasferimento nel tempo.
ESEMPIO 1
In questo esempio di antitrasformata N ( s ) = 1
Data la seguente A ( s) determiniamo la corrispondente f ( t ) cioè l'antitrasformata:
Si ricavano i poli di A ( s ) cioè le radici di s2 + 4 s + 3 : P1 = - 1 e P2 = - 3.
È possibile ora lo sviluppo in frazioni parziali:
I valori delle costanti A e B si ricavano ponendo: A ( s+3 ) + B (s+1) = 1 cioè
s ( A+B) + 3A + B = 1 per cui dobbiamo porre A+B = 0 e 3A+B = 1 affinché l'eguaglianza sia verificata.
Risolvendo il precedente sistema si ottiene : A = 0,5 e B = - 0,5.
Sostituendo i valori di A e B calcolati nello sviluppo in frazioni parziali, si ha:
Per antitrasformare si ricorre alla tabella:
• per la prima frazione parziale si ha K = 0,5 ed a = 1
• per la seconda frazione parziale si ha K = - 0,5 ed a = 3
Pertanto:
a ( t ) = L - 1 [ A ( s ) ] = 0,5 e - t - 0,5 e - 3 t
Graficamente l'andamento di a ( t ) sarà:
ESEMPIO 2
In questo esempio la A(s) ha uno zero: s = -1 e due poli: P1 = -3 e P2 = -5 .
Determinare l'antitrasformata della seguente funzione:
Il denominatore è già scomposto in due binomi, pertanto si effettua solo la scomposizione in frazioni parziali:
L’eguaglianza dei numeratori fornisce:
s+1 = A (s+5) + B (s+3) ;
s+1 = A s + 5 A + B s + 3 B ;
s+1 = s ( A+B ) + ( 5A + 3B ) ;
da cui considerando l'uguaglianza tra i due membri:
Sostituendo nello sviluppo in frazioni parziali:
Dalla tabella ricaviamo:
a ( t ) = L - 1 [ A ( s ) ] = - e - 3 t + 2 e - 5 t
CASI PARTICOLARI DELLO SVILUPPO IN FRAZIONI PARZIALI
1. POLI MULTIPLI
Se il denominatore della funzione di trasferimento si annulla per valori di “ s ” coincidenti (POLI MULTIPLI), cioè nella scomposizione in fattori del denominatore si ottiene un binomio del tipo: ( s - p ) n occorre individuare tante frazioni parziali quant’è l’ordine di molteplicità “ n ” del polo multiplo; queste frazioni presentano al numeratore una costante come di consueto, mentre i denominatori sono potenze decrescenti del binomio “ s-p ” da “n” a 1.
ESEMPIO
Antitrasformare:
Lo sviluppo in frazioni parziali sarà:
A ( s + 3 ) 2 + B ( s + 2 ) + C ( s + 3 ) ( s + 2 ) = s ( s + 1 )
s 2 (A+C) + s ( 6A+B+5C) + 9 A + 2 B + 6 C = s 2 + s
eguagliando i coefficienti di " s " di grado eguale e i termini noti, si ottiene:
Risolvendo si ottiene: A = 2; B = - 6; C = - 1; e quindi ( vedi tabella n° 9 ):
a ( t ) = L - 1 [ A ( t ) ] = 2 e - 2 t - 6 e - 3 t - e - 3 t
2. POLI COMPLESSI CONIUGATI
Se i poli di A(s) risultano complessi coniugati, il denominatore deve essere portato alla forma normale: D ( s ) = ( s + a )2 + +2 mentre il numeratore dello sviluppo in frazioni parziali assume la forma di un binomio in “s”, cioè: N(s) = As + B ; dopo di che si procede con il metodo consueto per la ricerca dei coefficienti A e B, riconducendosi poi alle relazioni n°11 e n°18 della TABELLA.
ESEMPIO
Antitrasformare:
I poli sono: P1 = - 2 – j e P2 = - 2 + j ; pertanto il denominatore avrà la forma normale:
D ( s ) = ( s + a )2 + +2 e quindi: 2a = 4 → a = 2; a2 + 2 = 5; 2 = 1 → = 1[rad/s].
Si ricordi che S costituisce il coefficiente dell’immaginario della variabile complessa “s” e rappresenta una pulsazione ( = 2 f), pertanto non vengono presi in considerazione i valori negativi.
Con i risultati precedenti:
L’eguaglianza dei numeratori fornisce: A = 0 e B = 1, da cui:
Per l’antitrasformazione si usa il n°11 della TABELLA con a = 2 e = 1 [ rad / s ] :
a ( t ) = e - 2 t sin ( t )
TABELLA
f ( t ) per t > 0
F(s)
1
K f ( t )
K F(s)
2
u ( t ) = 1
gradino unitario
1 / s
3
( t )
impulso unitario di DIRAC
1
4
t
rampa unitaria
1 / ( s2 )
5
t n
t n - 1 .
( n - 1 ) !
n intero > 0
n!
s n + 1
1 / ( s n )
6
e - a t
1 / ( s + a )
7
( 1 / a ) ( 1 - e- a t)
1 / [ s ( s + a ) ]
8
e -a t - e - b t.
b – a
1 .
(s + a) (s +b)
9
t e - a t
1 / [(s + a)2]
10
t n e - a t
1 . (t n - 1 ) e - a t
(n-1)!
con n intero e > 0
n! .
(s + a)n+1
1 / [ ( s + a ) n ]
11
e - a t sin ( s t )
.
(s + a)2 + 2
12
K sin((t + t)
con c = arctg ( / b)
K = radice quadrata ( b2++2)
s + b .
s2 + 2
13
( K t + 1) e - a t
con K = b – a
s + b .
(s + a)2
14
sin ( s t )
.
s2 + 2
15
cos ( c t )
s .
s2 + 2
16
t sin ( t )
2 2 .
(s2 + 2)2
17
t cos ( t )
s2 - 2 .
(s2 + 2)2
18
K e - a t sin ( t + t)
con c = arctg [ / ( b - a )]
K = radice quadrata[(b - a)2++2]
s + b .
(s + a)2 + 2
I DIAGRAMMI DI BODE
Un generico sistema è caratterizzato dalla sua FUNZIONE DI TRASFERIMENTO
G ( s ); per la sua struttura, essa risulta indipendente dal segnale di ingresso applicato, pertanto dipende solo dalla costituzione interna del sistema. Ponendo s = j s, la funzione di trasferimento che si ottiene è detta “risposta in frequenza” del sistema. Tale funzione ha un significato ben preciso, in base al concetto di scomposizione di un generico segnale periodico in somma di componenti sinusoidali di varie pulsazioni .
TEOREMA: un sistema lineare stazionario asintoticamente stabile, soggetto ad un ingresso sinusoidale presenta, a regime, una risposta sinusoidale, avente la stessa frequenza del segnale d’ingresso e dipendente dalla funzione di trasferimento G ( j G ).
Per rappresentare graficamente la “risposta in frequenza” di un sistema, si utilizzano i diagrammi di Bode, che riportano l’andamento del modulo e della fase di G ( j G ).
Il DIAGRAMMA DEL MODULO ha, sull’asse delle ordinate, il modulo della G (j G) espresso in decibel e, sull’asse delle ascisse, il logaritmo in base 10 della pulsazione b.
Il DIAGRAMMA DELLA FASE riporta il valore della fase di G ( j G ) in funzione sempre del logaritmo di .
Il valore del modulo della “risposta in frequenza” diventa:
| G ( j | ) | d b = 20 log 10 | G ( j | ) |
I grafici in scala logaritmica, rispetto alla scala lineare, hanno i seguenti vantaggi:
• capacità di rappresentare un campo molto esteso di valore di c;
• rappresentazione di una funzione data in forma fattorizzata come somma di diagrammi corrispondenti ad ogni singolo fattore.
Si supponga di voler rappresentare mediante i diagrammi di Bode la seguente funzione:
Raccogliendo al numeratore (-Z1) ed al denominatore (-P1):
Ponendo K = K1 ( Z1 / P1 ); z = 1 / - Z1; ;p = 1 / - P1 si ricava:
Che esprime la risposta in frequenza del sistema nella cosiddetta “forma di Bode”.
Applicando il logaritmo ad entrambi i membri si ricava:
log10 | G ( j | ) | = log10 | k | + log10 | 1+ j||z | + log10 | 1 / (1+ j//p ) |
Per ottenere il diagramma di Bode, è sufficiente tracciare i diagrammi dei termini costituenti la risposta, quindi sommare tutti i contributi.
DIAGRAMMI ELEMENTARI DI BODE
COSTANTE: G ( j C ) = K
Si esegue il modulo e lo si esprime in decibel | G ( j ) | dB = 20 lg 10 K ( esso non dipende da ); la fase ) = arctg ( 0 / K ) ; / = 0° se K > 0 ; mentre > = 180° se K < 0.
[ = w ]
ZERO NELL'ORIGINE : G ( j Z ) = j
Il modulo | G ( j ( ) | dB = 20 lg 10 rappresenta una retta avente pendenza
20 dB / decade e passante per l'origine degli assi; la fase / = acrtg ( / 0 ) = 90°
[[ = w ]
POLO NELL'ORIGINE: G ( j ) = 1 / j /
Il modulo | G ( j I ) | dB = - 20 lg 10 rappresenta una retta avente pendenza
- 20 dB / decade e passante per l'origine degli assi; la fase / = - acrtg ( a / 0 ) = - 90°
[[ = w ]
ZERO REALE: G ( j Z ) = 1 + j Z
Si esegue il modulo e lo si esprime in decibel | G ( j S ) | dB = 20 lg 10 | 1+ j | . Per comodità si preferisce una rappresentazione approssimata ( per asintoti ) della funzione ottenuta:
• per p z > 1 si ha | G ( j > ) | dB = 20 lg + 20 lg z
Il diagramma del modulo è approssimato: le due rette che lo compongono sono gli asintoti ai quali tende la funzione effettiva per valori di z sufficientemente lontani dalla pulsazione relativa allo zero. Infatti, in tale punto, si ha lo scostamento massimo del diagramma reale da quello approssimato; la funzione reale, per la pulsazione . = 1 / / z , detta pulsazione di taglio, assume il valore di 3 dB. L'errore diminuisce all'allontanarsi dalla pulsazione di taglio e, a distanza di una decade, si può ritenere trascurabile.
Per la fase P = arctg ( z / 1 ) si determina la rappresentazione grafica approssimata:
• per p z 1 ossia > >> 1 / / z si ha = 90° [ z = T z ]
POLO REALE : G ( j P ) = 1 / ( 1 + j / p )
Il modulo: | G ( j ( ) | dB = - 20 lg 10 | 1 + j p |
Anche in questo caso per comodità si fa una rappresentazione approssimata della funzione:
• per p p > 1 si ha | G ( j > ) | dB = - 20 lg - 20 lg p
Le semirette del diagramma del modulo sono gli asintoti ai quali tende la funzione per valori di a sufficientemente lontani dalla pulsazione del polo. Alla pulsazione di taglio A = 1 / / p , la funzione reale assume il valore di - 3 dB; l'errore diminuisce allontanandosi dalla pulsazione di taglio e, a distanza di una decade, si può ritenere trascurabile.
Per la fase P = - arctg ( p / 1 ) si determina la rappresentazione grafica approssimata:
• per p p 1 ossia > >> 1 / / p si ha = - 90° [ p =T p ]
1
1