Trasformata di Laplace Numeri complessi Un numero complesso z, costituito da una parte reale ed una immaginaria, è scritto come: z  x  jy, x, y   Re[ z ] j: Im[ z ] j 2  1 Z y  O y z  x  jy x z *  x  jy Z* Forma cartesiana Numeri complessi z  zz *  x  jy x  jy  Modulo  x 2  jxy  jxy  y 2  x 2  y 2 y tg  x y   arctg  2k , k  0 intero x Fase Numeri complessi y jy z  x  jy z x  z cos   x x z *  x  jy y  z sin  z  z cos   j sin   e j  cos   j sin  Ricordando le formule di Eulero: e j  e  j  2 cos  e j  e  j  2 j sin  e  j  cos   j sin  e j  e  j cos   2 e j  e  j sin   2j Forma trigonometrica Numeri complessi y z z  x  jy  z  z e j x Forma esponenziale z *  x  jy z *  x  jy  z cos   j z sin   z e  j z  z  Forma polare Numeri complessi 1 1 x  jy z  x  jy x  jy x  jy  x  jy x y  2   j x  y2 x2  y2 x2  y2 y  2 y x  y2  z   arctg  arctg x x x2  y2 z  x x2 2  y2   x 2 y2 2  y2  2 x2  y2  2  2 x y 1 x2  y2 La trasformata di Laplace è un operatore che associa ad una funzione del tempo f(t) definita per t≥0 una funzione F(s) a valori complessi definita per valori della variabile complessa s. f (t ), t  0 F (s) s    j L’utilizzo delle trasformate di Laplace consente di semplificare notevolmente i calcoli nella risoluzione di equazioni differenziali: operazioni di derivazione ed integrazione nel dominio del tempo corrispondono ad operazioni di tipo algebrico nel dominio delle trasformate. Problema differenziabile L Problema algebrico Soluzione del problema differenziabile 1 L Soluzione del problema algebrico  F ( s)   f (t )e  st dt 0 Qualsiasi funzione f(t), per cui esiste un valore della variabile s tale che l’integrale è finito, si dice trasformabile secondo Laplace. L’insieme di tutti i valori complessi s per cui esiste, finito, l’integrale e quindi la funzione F(s), viene detto dominio di convergenza, ed è rappresentato da un semipiano del piano s, posto a destra di una retta parallela all’asse immaginario, di equazione Re[s]=σ0. Tale retta viene denominata asse di convergenza ed il valore σ0 ascissa di convergenza. Gradino  T L[1]   e  st dt  lim  e  st dt  0 T  0 T  e  e  sT 1 1  lim    lim    T  T  s s s  s 0  st Re[ s ]  0 e sT  e   j T  e T e  jT  e T cos(T )  j sin( T ) Im[s] k L[ k ]  s Re[s] Proprietà di linearità L f1 (t )  F1 (s) L f 2 (t )  F2 (s) Lc1 f1 (t )  c2 f 2 (t )  c1F1 (s)  c2 F2 (s) Rampa unitaria  T L[t ]   te st dt  lim  te st dt  T  0  te st  lim  T   s  te  lim  T   s T 0  st T 0 0 e  st   dt   s  0 T  fg ' dt  fg   f ' gdt T 1  st   2e   s  0  Te  sT 1  sT 1  1  lim   2 e  2 2 T  s s s  s  Re[ s ]  0 Esponenziale  T L[e at ]   e at e  st dt  lim  e ( a  s )t dt  0 T  0 T  1 ( a  s )t   lim  e   T  a  s  0 1 ( a  s )T 1 1  lim e   T  a  s as sa Re[ s]  Re[ a] Cosinusoide  e jt  e  jt  Lcost   L   2   1  L e jt  L e  jt  2 1 1 1  s     2  s  j s  j  s 2   2     Re[ s ]  0 Sinusoide  e jt  e  jt  Lsin t   L   2j   1  L e jt  L e  jt  2j      1  1 1     2 j  s  j s  j  s 2   2 Re[ s ]  0 Traslazione   L e kt f (t )  F (s  k )   0 0 kt  st ( s  k ) t e f ( t ) e dt  f ( t ) e dt F ( s  k )   L f (t  k )  e  ks F ( s)   f (t  k )e  st dt  f (t  k )  0 , tk u t k 0    k  f (u )e  s (u  k ) du  e  sk  f (u )e  su du e  sk F ( s ) 0 f (t ) Impulso A t0 A f (t )  , 0  t  t0 t0  0, t  0, t  t0 0 t0 L’area sottesa vale A A A f (t )  1(t )  1(t  t0 ) t0 t0 A 1 A 1  st0 F (s)   e  t0 s t0 s  A1  1  e  st0 t0 s  Funzione impulsiva A g (t )  lim , 0  t  t0 t0 0 t 0  0, t  0, t  t0   A1 G ( s )  lim 1  e  st0  t 0 0 t s 0 A  lim t0 0   d 1  e  st0 dt0 Ase  st0  lim A t0 0 d s t0 s  dt0 Impulso di Dirac 1  (t )  lim , 0  t  t0 t0 0 t 0  0, t  0, t  t0 (s)  1 Impulso di Dirac Risposta all’impulso u (t ) y(t )   t   h(t ) y (t )  t Ogni segnale x(t) può essere espresso come convoluzione con l’impulso di Dirac x(t )  x(t )   (t )    (t ) x(t   )d 0 Dim:   1 1 x(t   )d  lim  x(t   )d 0 lim  0   0  0 Per il teorema del valor medio:   0,  :   x(t  )d  xt    0 1 lim xt     lim xt     x(t )  0   0 0   Impulso di Dirac u (t ) y (t )  Il segnale in uscita può essere calcolato attraverso la convoluzione del segnale di ingresso con la risposta impulsiva. L’uscita del sistema all’ingresso x(t) sarà del tipo: t  y (t )  x(t )      (t ) x(t   )d  0  t y (t )   x(t   ) (t )d  0 t   h(t ) x(t   )d  h(t )  (t ) 0 Impulso di Dirac Problemi La risposta impulsiva di un sistema può essere ricavata applicando in ingresso un segnale che approssimi l’impulso di Dirac e misurando l’uscita corrispondente. L’impulso di Dirac è un’astrazione matematica che può solo essere approssimata. In molti casi non è possibile né conveniente applicare al sistema una sollecitazione impulsiva per non danneggiare il sistema a causa dell’elevata ampiezza dell’impulso. Esercizio   n! L t  n 1 s Sapendo che calcolare n  L t 2e3t  Esercizio Calcolare  L e 2t sin( 4t )  Esercizio Calcolare   2  2  cos t  3 , t  3   L   2   0, t  3  Esercizio Calcolare Lt sin t  Teorema della derivata   L f (t )  sF (s)  f (0) '    fg ' dt  fg   f ' gdt  L f ' (t )   f ' (t )e  st dt  0 T   T  st  st  s  f (t )e dt    lim  f (t )e 0 T  0   T    st  sT  lim  f (T )e  f (0)  s  f (t )e dt   T  0     s  f (t )e  st dt  f (0)  sF ( s )  f (0) 0 Si è sfruttato il fatto che f(t) è di ordine esponenziale per t che tende all’infinito Teorema della derivata   L f '' (t )  s 2 F ( s)  sf (0)  f ' (0)     L f '' (t )  sL f ' (t )  f ' (0)  ssF ( s )  f (0)  f ' (0)  s 2 F ( s )  sf (0)  f ' (0)   L f ''' (t )  s 3 F ( s)  s 2 f (0)  sf ' (0)  f '' (0) Teorema dell’integrale t g ' (t )  f (u), g (t )   f (u )du g (0)  0 0   L g ' (t )  F ( s) sLg (t )  F ( s) F (s) Lg (t )  s t  F (s) L   f (u )du   s 0  Teorema del valore finale lim f (t )  lim sF ( s ) t  Nell’ipotesi che tale limite esista s 0   Dal teorema della derivata si ha: f ' (t )e  st dt  sF ( s)  f (0) 0  da cui lim s 0  f ' (t )e  st dt  lim sF ( s)  f (0) s 0 0 Eseguendo il limite sotto il segno di integrale il che è lecito per l’analiticità della funzione:   0 f ' (t )dt  lim sF ( s)  f (0) s 0 e quindi lim f (t )  f (0)  lim sF ( s )  f (0) t  s 0 Teorema del valore iniziale lim f (t )  lim sF ( s ) t 0 s   lim s   f ' (t )e  st dt  lim sF ( s )  f (0) s  0  lim sF ( s )  f (0) 0 s  lim f (t )  lim sF ( s ) t 0 s  Integrale di convoluzione   L   f1   f 2 t   d   F1 ( s ) F2 ( s ) 0  Utilità x  3x  2 x  0, x(0)  a, x (0)  b Lx(t )  X ( s ) Lx (t )  sX ( s )  x(0) Lx(t )  s 2 X ( s )  sx (0)  x ' (0) s 2 X (s)  as  b  3sX (s)  a   2 X (s)  0   X ( s) s 2  3s  2  as  b  3a as  b  3a 2a  b a  b X ( s)    s  1s  2 s  1 s  2 Utilità as  b  3a 2a  b a  b X ( s)    s  1s  2 s  1 s  2  2a  b  1  a  b  x(t )  L X ( s )  L  L    s 1   s  2  1 1  1  1  1  x(t )  (2a  b) L1   ( a  b ) L  s  2   s  1 x(t )  (2a  b)e t  (a  b)e 2t Problema differenziabile L Problema algebrico Soluzione del problema differenziabile 1 L Soluzione del problema algebrico Tecniche di antitrasformazione Frazione razionale propria F ( s)  Il denominatore di F(s) ha: n radici distinte radici con molteplicità maggiore di 1 radici complesse coniugate N ( s) D( s ) Tecniche di antitrasformazione Il denominatore di F(s) ha: n radici distinte N ( s) F (s)  s  p1 s  p2  s  pn  F (s)  Rn R1 R2   s  p1  s  p2  s  pn  POLI RESIDUI Tecniche di antitrasformazione Il denominatore di F(s) ha: n radici distinte Rn N (s) R1 R2    s  p1 s  p2 s  pn  s  p1  s  p2  s  pn  Calcoliamo R1 s  p1 N ( s) s  p1 R1  s  p1 R2    s  p1 Rn  s  p1 s  p2 s  pn  s  p1  s  p2  s  pn  s  p1 R2    s  p1 Rn N ( s)  R1  s  p2 s  pn  s  p2  s  pn  N ( s) lim  R1 s  p1 s  p  s  p  2 n Tecniche di antitrasformazione Il denominatore di F(s) ha: n radici distinte N ( s) F (s)  s  p1 s  p2  s  pn  Rn R1 R2 F (s)    s  p1  s  p2  s  pn  Rk  lim F (s)s  pk , s  pk k  1,2,, n Tecniche di antitrasformazione Il denominatore di F(s) ha: n radici distinte  1  1  1  1  1  L F ( s )  R1 L     R2 L      Rn L        s  p s  p s  p 1  2  n     1 1 n   Rk e pk t k 1 Esercizio s2 1 F (s)  3 s s Tecniche di antitrasformazione Il denominatore di F(s) ha: radici con molteplicità maggiore di 1 N ( s) F ( s)   s  pk   Rk Rk1 Rk 2 F (s)      1 s  pk  s  pk  s  pk  Tecniche di antitrasformazione Il denominatore di F(s) ha: radici con molteplicità maggiore di 1 Rk Rk1 Rk 2 N (s) F ( s)        1 s  pk  s  pk  s  pk  s  pk  Calcoliamo Rk1    s  p Rk s  pk  N (s)  s  pk  Rk1  s  pk  Rk 2    k s  pk  s  pk  s  pk  s  pk  1    Rk1  lim F ( s )s  pk   s  pk Tecniche di antitrasformazione Il denominatore di F(s) ha: radici con molteplicità maggiore di 1 Calcoliamo Rk2 N ( s)  Rk1  s  pk Rk 2  s  pk  Rk 3  s  pk  Rk 4    s  pk  2 3  1 Rk dN ( s ) 2  2  Rk 2  2s  pk Rk 3  3s  pk  Rk 4   (   1)s  pk  Rk ds  d F ( s )s  pk  Rk 2  lim s  pk ds   Tecniche di antitrasformazione Il denominatore di F(s) ha: radici con molteplicità maggiore di 1 Calcoliamo Rk3 dN ( s ) 2  2  Rk 2  2s  pk Rk 3  3s  pk  Rk 4   (   1)s  pk  Rk ds d 2 N ( s)  3      2 R  6 s  p R   (   1 )(   2 ) s  p Rk k3 k k4 k 2 ds  d 2 F ( s )s  pk  1 Rk 3  lim 2 s  pk ds 2   Tecniche di antitrasformazione Il denominatore di F(s) ha: radici con molteplicità maggiore di 1 Calcoliamo Rkj   d j 1 F ( s)s  pk  1 Rkj  lim , j 1 s  p ( j  1)! k ds Ricordando che   L t n e at   j  1,2,,  n! s  a n1  1  1 k 1 pt L   t e k   s  p   (k  1)! 1 Esercizio s 3  5s F ( s)  s  1s  23 Tecniche di antitrasformazione Il denominatore di F(s) ha: radici complesse coniugate F ( s)  F ( s)  N ( s) s  p  s  p*   p  a  jb p *  a  jb N ( s) R1 R2   s  p  s  p* s  p s  p*   N ( s) N ( p) N (a  jb) R1  lim    u  jv * s p s  p* p p j 2b N ( s) N ( p* ) N (a  jb) * R2  lim*  *   R1  u  jv s p s  p p p  j 2b Tecniche di antitrasformazione Il denominatore di F(s) ha: radici complesse coniugate u  jv u  jv F ( s)   s    j  s    j  L1 F ( s )  u 2  v 2 e jarctg v u e   j t  u 2  v 2 e  jarctg v u v v    j  t  arctg   j  t  arctg   u u  u 2  v 2 et e  e      j  t  arctg uv   j  t  arctg uv      e   1 2 2 t  e L F ( s )  2 u  v e   2     v  L F ( s )  2 u  v e cos t  arctg  u  1 2 2 t e   j t  Esercizio Funzione di trasferimento Funzione di trasferimento Funzione di trasferimento