Transformacja Laplace’a

Definicja podstawowa

Transformacją Laplace’a (przekształceniem Laplace'a) L {\displaystyle {\mathcal {L}}} nazywamy przekształcenie całkowe, które funkcji f ( t ) {\displaystyle f(t)} zmiennej rzeczywistej t {\displaystyle t} przyporządkowuje funkcję zespoloną F ( s ) {\displaystyle F(s)} zmiennej zespolonej s {\displaystyle s} za pomocą całki (zwanej całką Laplace'a)[1]

F ( s ) = 0 e s t f ( t ) d t {\displaystyle F(s)=\int \limits _{0}^{\infty }e^{-st}f(t)\,dt}

Zbiór X {\displaystyle X} wszystkich funkcji, dla których powyższa całka jest zbieżna (tj. ma skończoną wartość w punktach s {\displaystyle s} przynajmniej pewnej części płaszczyzny zespolonej), tworzy pewną przestrzeń funkcyjną.

Transformatą Laplace'a funkcji rzeczywistej f {\displaystyle f} nazywamy funkcję zespoloną F L ( f ) {\displaystyle F\equiv {\mathcal {L}}(f)} , przyporządkowaną jej za pomocą powyższej całki. Możemy więc zapisać

L : X f L ( f ) {\displaystyle {\mathcal {L}}:X\ni f\to {\mathcal {L}}(f)}

Przyjmuje się różne symbole na oznaczanie transformaty Laplace'a, wskazujące na punkt s {\displaystyle s} , do którego odnosi się wartość transformaty lub wskazujące na funkcję pierwotną f {\displaystyle f} czy też jej wartość f ( t ) {\displaystyle f(t)} :

{ L f } ( s ) {\displaystyle \{{\mathcal {L}}f\}(s)} , F ( s ) , L { f ( t ) } {\displaystyle F(s),{\mathcal {L}}\{f(t)\}} , L [ f ( t ) ] {\displaystyle {\mathcal {L}}[f(t)]}

Transformata Laplace'a jest szczególnie przydatna w rozwiązywaniu liniowych równań różniczkowych zwyczajnych, takich jak te powstające w analizie obwodów elektronicznych.

Matematykiem, który zdefiniował transformację Laplace’a i od którego nazwiska wzięła ona nazwę był Pierre Simon de Laplace.

Warunki dostateczne zbieżności całki Laplace’a

Funkcje, dla których całka Laplace'a jest na pewno zbieżna - przynajmniej w pewnej części płaszczyzny zespolonej, nazywa się oryginałami laplace'owskimi (krótko: oryginałami). Poza oryginałami istnieją inne funkcje, dla których całka Laplace'a jest zbieżna, ale nie da się ustalić wspólnych dla nich kryteriów i każdy przypadek wymaga osobnego zbadania.

Df. Oryginałem jest funkcja o dziedzinie rzeczywistej i wartościach w zbiorze liczb rzeczywistych, mająca następujące właściwości:[2]

  1. f ( t ) = 0  dla  t < 0 {\displaystyle f(t)=0{\text{ dla }}t<0}
  2. w każdym otwartym przedziale skończonym spełniony jest pierwszy i drugi warunek Dirichleta
  3. Istnieje funkcja wykładnicza, która majoryzuje, czyli ogranicza wykładniczo funkcję f ( t ) , {\displaystyle f(t),} tj. istnieją takie liczby M > 0 , ρ 0 {\displaystyle M>0,\rho \geq 0} że zachodzi nierówność
    | f ( t ) | M e d t {\displaystyle |f(t)|\leq Me^{dt}\,} dla wszystkich t R {\displaystyle t\in R}

Zbieżność w całej płaszczyźnie całek Laplace'a z oryginałów

Półpłaszczyzna zbieżności bezwzględnej. Pas zbieżności warunkowej

Jeżeli funkcja f ( t ) {\displaystyle f(t)} jest oryginałem, to:[3]

Półpłaszczyzny zbieżności-rozbieżności transformaty Laplace'a na płaszczyźnie zespolonej
  1. istnieje taka liczba x z {\displaystyle x_{z}} , że całka Laplace'a jest zbieżna dla punktów s {\displaystyle s} takich, że R e s > x z {\displaystyle Re\,s>x_{z}} ; tym samym prosta x = x z {\displaystyle x=x_{z}} dzieli płaszczyznę zespoloną S na półpłaszczyznę zbieżności ( R e s > x z {\displaystyle Re\,s>x_{z}} ) i półpłaszczyznę rozbieżności ( R e s < x z {\displaystyle Re\,s<x_{z}} ). Na prostej x = x z {\displaystyle x=x_{z}} całka Laplace'a może być zbieżna lub rozbieżna.
  2. istnieje taka liczba x a {\displaystyle x_{a}} , że całka Laplace'a jest bezwzględnie zbieżna dla punktów s {\displaystyle s} takich, że R e s > x a {\displaystyle Re\,s>x_{a}} ; tym samym prosta x = x a {\displaystyle x=x_{a}} dzieli płaszczyznę zespoloną S {\displaystyle S} na półpłaszczyznę bezwzględnej zbieżności ( R e s > x a {\displaystyle Re\,s>x_{a}} ) i półpłaszczyznę, w której całka Laplace'a nie jest bezwzględnie zbieżna ( R e s < x a {\displaystyle Re\,s<x_{a}} ). Na prostej x = x a {\displaystyle x=x_{a}} całka Laplace'a może być bezwzględnie zbieżna lub nie być bezwzględnie zbieżna.
  3. zawsze jest x z < x a {\displaystyle x_{z}<x_{a}} lub x z = x a {\displaystyle x_{z}=x_{a}} ; gdy x z < x a {\displaystyle x_{z}<x_{a}} , to w każdym punkcie pasa nieskończonego x z < R e s < x a {\displaystyle x_{z}<Re\,s<x_{a}} całka Laplace'a jest zbieżna warunkowo; pas ten nazywa się pasem warunkowej zbieżności.

Zbieżność transformat Laplace'a z oryginałów i innych funkcji w całej płaszczyźnie[4]

Obliczając transformatę Laplace'a nie należy przywiązywać dużego znaczenia do tego, jaka jest półpłaszczyzna zbieżności. Ważny jest jedynie fakt, że półpłaszczyzna rozbieżności w ogóle istnieje. Bowiem w takiej sytuacji transformatę można rozpatrywać jako zbieżną na całej płaszczyźnie zespolonej.

Np. całka Laplace'a funkcji skokowej Heaviside'a jest zbieżna tylko na półpłaszczyźnie R e s > 0 {\displaystyle Re\,s>0} , natomiast transformatę 1 / s {\displaystyle 1/s} tej funkcji można rozpatrywać jako funkcję holomorficzną określoną na całej płaszczyźnie zespolonej, przy czym w punkcie s = 0 {\displaystyle s=0} funkcja ta ma biegun pojedynczy.

Wnioski:

  1. Na pewno dla funkcji, które są oryginałami, istnieje półpłaszczyzna zbieżności, a więc i ich transformaty Laplace'a można rozpatrywać jako zbieżne na całej płaszczyźnie zespolonej - dlatego badanie, czy funkcja jest oryginałem wystarczy, by stosować do niej cały formalizm transformaty Laplace'a. Stąd też wynika znaczenie oryginałów.
  2. Powyższy wniosek dotyczy wszystkich innych funkcji, które mają półpłaszczyznę zbieżności.

Przykłady obliczeń

Funkcja Heaviside’a; 1 ( 0 ) = 0 . 5 {\displaystyle \mathbf {1} (0)=0{.}5}

Przykład 1. Transformata funkcji skokowej

Oblicz transformatę Laplace'a funkcji skokowej Heaviside’a 1 ( t ) {\displaystyle \mathbf {1} (t)} (tzw. funkcji jednostkowej), która przyjmuje wartość 0 {\displaystyle 0} dla ujemnych chwil czasu, wartość 0.5 {\displaystyle 0.5} dla chwili t = 0.5 {\displaystyle t=0.5} , zaś wartość 1 {\displaystyle 1} dla dodatnich chwil czasu, tj.

1 ( t ) = { 0   dla  t < 0 1 2 dla  t = 0 1 dla  t > 0 {\displaystyle \mathbf {1} (t)={\begin{cases}0\ \quad \quad {\text{dla }}\,t<0\\{\frac {1}{2}}\quad \quad {\text{dla }}\,t=0\\1\quad \quad \,{\text{dla }}\,t>0\end{cases}}}

Funkcja ta odgrywa szczególną rolę w teorii transformaty Laplace'a.

Rozwiązanie: Obliczamy:[3]

L { 1 ( t ) } = 0 1 ( t ) e s t d t = lim T + 0 T e s t d t = lim T + ( e s t s | 0 T ) = {\displaystyle {\mathcal {L}}\{\mathbf {1} (t)\}=\int \limits _{0}^{\infty }\mathbf {1} (t)e^{-st}\,dt=\lim _{T\to +\infty }\int \limits _{0}^{T}e^{-st}\,dt=\lim _{T\to +\infty }{\Bigg (}-{\frac {e^{-st}}{s}}{\Bigg |}_{0}^{T}{\Bigg )}=}
= lim T + ( 1 s 1 s e s T ) = lim T + 1 s [ 1 e x T ( cos ω T i sin ω T ) ] {\displaystyle =\lim _{T\to +\infty }{\Big (}{\frac {1}{s}}-{\frac {1}{s}}e^{-sT}{\Big )}=\lim _{T\to +\infty }{\frac {1}{s}}{\big [}1-e^{-xT}(\cos \omega T-i\sin \omega T){\big ]}}

przy czym w ostatnim kroku obliczeń podstawiono s = x + i ω {\displaystyle s=x+i\omega } oraz skorzystano z wzoru Eulera ( e i ω T = cos ω T i sin ω T {\displaystyle e^{-i\omega T}=\cos \omega T-i\sin \omega T} ). Widać, że obliczana całka jest dla x > 0 {\displaystyle x>0} zbieżna i równa 1 s {\displaystyle {\frac {1}{s}}} , zaś rozbieżna dla x 0 {\displaystyle x\leq 0} , a zatem x z = 0 {\displaystyle x_{z}=0} . Ponadto x a = 0 {\displaystyle x_{a}=0} , ponieważ

0 + | e s t | d t = 0 + e x t x d t = { 1 x dla  x > 0 +   dla  x 0 {\displaystyle \int \limits _{0}^{+\infty }{\big |}e^{-st}{\big |}\,dt=\int \limits _{0}^{+\infty }e^{-xt}x\,dt={\begin{cases}{\frac {1}{x}}\quad \quad \,{\text{dla }}x>0\\+\infty \ \quad {\text{dla }}x\leq 0\end{cases}}}

Ostatecznie transformata Laplace'a funkcji jednostkowej wynosi:

L { 1 ( t ) } = 1 s {\displaystyle {\mathcal {L}}\{\mathbf {1} (t)\}={\frac {1}{s}}}

Uwaga:

W Tablicy transformat (dalej) podana jest wartość transformaty L { a } = a 1 s {\displaystyle {\mathcal {L}}\{a\}=a{\frac {1}{s}}} ; nie jest to jednak transformata z funkcji stałej f ( t ) = a {\displaystyle f(t)=a} ; w zapisie tym milcząco zakłada się, że funkcja pierwotna ma postać f ( t ) = a 1 ( t ) {\displaystyle f(t)=a\cdot \mathbf {1} (t)} , co zapewnia spełnienie warunku, że funkcja jest oryginałem laplace'owskimi, a to gwarantuje istnienie transformaty.

Przykład 2. Transformata funkcji eksponencjalnej

Niech b {\displaystyle b} będzie dowolną liczbą zespoloną:[5]

L { e a t 1 ( t ) } = 0 + e a t e s t d t = 0 + e ( a + s ) t d t {\displaystyle {\mathcal {L}}\{e^{-at}\mathbf {1} (t)\}=\int \limits _{0}^{+\infty }e^{-at}e^{-st}\,dt=\int \limits _{0}^{+\infty }e^{-(a+s)t}\,dt}

Rozwiązanie: Suma a + s {\displaystyle a+s} jest stałą liczbą (podobnie jak wcześniej s {\displaystyle s} ) . Porównując tę całkę z obliczoną w przykładzie 1 możemy natychmiast napisać wynik

L { e a t 1 ( t ) } = 1 s + a {\displaystyle {\mathcal {L}}\{e^{-at}\mathbf {1} (t)\}={\frac {1}{s+a}}}

Interpretacja i związek z transformatą Fouriera

Postacie funkcji R e ( e s t ) {\displaystyle Re\,(e^{-st})} dla różnych zespolonych częstotliwości s = σ + i ω {\displaystyle s{=}\sigma {+}i\omega } , co można zapisać w postaci e σ t cos ( ω t ) . {\displaystyle e^{-\sigma t}\cos(\omega t).} σ = 0 {\displaystyle \sigma {=}0} zawiera tylko cosinusy. Dodatnie σ {\displaystyle \sigma } zawiera cosinusy tłumione. Ujemne σ {\displaystyle \sigma } zawiera cosinusy rosnące eksponencjalnie.

Za pomocą całki Laplace'a funkcja rzeczywista f ( t ) {\displaystyle f(t)} , reprezentująca zmieniający się w czasie t {\displaystyle t} sygnał (np. pole elektryczne fali, przychodzącej do odbiornika), jest transformowana na płaszczyznę S ( płaszczyznę zespoloną); dokonuje się to poprzez scałkowanie iloczynu funkcji f ( t ) {\displaystyle f(t)} z wyrażeniami typu e s t {\displaystyle e^{-st}} dla czasu t {\displaystyle t} od {\displaystyle -\infty } do + {\displaystyle +\infty } :

f ( t ) e s t d t ; {\displaystyle \int _{-\infty }^{\infty }f(t)e^{-st}\,dt\,;}

przy tym s {\displaystyle s} jest liczbą zespoloną, stałą w procesie obliczania całki.

Aby zrozumieć, co otrzymuje się w wyniku takiego działania, trzeba najpierw poznać działanie transformacji Fouriera (por. też analiza Fouriera), wyrażonej za pomocą analogicznej całki:

f ^ ( ω ) = f ( t ) e i ω t d t {\displaystyle {\hat {f}}(\omega )=\int \limits _{-\infty }^{\infty }f(t)e^{-i\omega t}dt}

W całce Fouriera funkcje harmoniczne e i ω t = cos ω t i sin ω t , ω R {\displaystyle e^{-i\omega t}=\cos \omega t-i\sin \omega t,\omega \in R} mnożone są przez sygnał f ( t ) {\displaystyle f(t)} ; wynikowa całka dostarcza informacji nt. zawartości poszczególnych harmonicznych, wchodzących w skład sygnału (dokonuje rozkładu sygnału na jego widmo).

Transformacja Laplace’a wykonuje podobne działanie, ale o bardziej ogólnym charakterze: funkcje e s t {\displaystyle e^{-st}} w całce Laplace'a zależą od zmiennej zespolonej s = x + i ω {\displaystyle s=x+i\omega } , czyli de facto mają postać e s t = e x t e i ω t {\displaystyle e^{-st}=e^{-xt}\cdot e^{-i\omega t}} ; dzięki temu pozwalają dokonać nie tylko analizy zawartości harmonicznych e i ω t {\displaystyle e^{-i\omega t}} w sygnale f ( t ) {\displaystyle f(t)} , ale również efektów zaniku sygnału w czasie i przestrzeni, poprzez funkcję e x t . {\displaystyle e^{-xt}.} Na przykład krzywa sinusoidalna tłumiona może być modelowana za pomocą transformacji Laplace'a. Transformacja Fouriera stanowi więc szczególny przypadek przekształcenia Laplace’a dla s = i ω . {\displaystyle s=i\omega .}

Podobnie uogólnieniem dyskretnej transformaty Fouriera stanowi transformata Z, z którą powiązana jest transformata Laplace’a (zob. metoda Tustina).

Własności

Poniżej zestawiono właściwości transformaty Laplace'a[6]

Liniowość

L { a f ( t ) + b g ( t ) } = a L { f ( t ) } + b L { g ( t ) } {\displaystyle {\mathcal {L}}\{af(t)+bg(t)\}=a{\mathcal {L}}\{f(t)\}+b{\mathcal {L}}\{g(t)\}}

Transformata pochodnej

(a) pierwsza pochodna

L { f } = s L { f } f ( 0 + ) , {\displaystyle {\mathcal {L}}\{f'\}=s{\mathcal {L}}\{f\}-f(0^{+}),}
gdzie:
f ( 0 + ) {\displaystyle f(0^{+})} - granica prawostronna funkcji f ( t ) {\displaystyle f(t)} w chwili t = 0 {\displaystyle t=0} ; jest to np. warunek początkowy dla funkcji f ( t ) {\displaystyle f(t)} będącej szukanym rozwiązaniem równania różniczkowego funkcji

(b) druga pochodna

L { f } = s 2 L { f } s f ( 0 + ) f ( 0 + ) {\displaystyle {\mathcal {L}}\{f''\}=s^{2}{\mathcal {L}}\{f\}-sf(0^{+})-f'(0^{+})} ,
gdzie:
f ( 0 + ) {\displaystyle f(0^{+})} , f ( 0 + ) {\displaystyle f'(0^{+})} - warunki początkowe dla funkcji f ( t ) {\displaystyle f(t)} w chwili t = 0 {\displaystyle t=0} ( f ( 0 + ) {\displaystyle f'(0^{+})} oznacza granicę prawostronną pochodnej funkcji f ( t ) {\displaystyle f(t)} w chwili t = 0 {\displaystyle t=0} )

(c) n-ta pochodna

L { f ( n ) } = s n L { f } k = 1 n s n k f ( k 1 ) ( 0 + ) = s n L { f } s n 1 f ( 0 + ) s n 2 f ( 0 + ) f ( n 1 ) ( 0 + ) {\displaystyle {\mathcal {L}}\{f^{(n)}\}=s^{n}{\mathcal {L}}\{f\}-\sum _{k=1}^{n}s^{n-k}f^{(k-1)}(0^{+})=s^{n}{\mathcal {L}}\{f\}-s^{n-1}f(0^{+})-s^{n-2}f'(0^{+})-\ldots -f^{(n-1)}(0^{+})} ,
gdzie f ( 0 + ) , f ( 0 + ) , , f ( n 1 ) ( 0 + ) {\displaystyle f(0^{+}),f'(0^{+}),\ldots ,f^{(n-1)}(0^{+})} - warunki początkowe funkcji f ( t ) {\displaystyle f(t)} w chwili t = 0 {\displaystyle t=0}

Pochodna n-tego rzędu transformaty

F ( n ) ( s ) = ( 1 ) n L { t n f ( t ) } {\displaystyle F^{(n)}(s)=(-1)^{n}{\mathcal {L}}\{t^{n}f(t)\}}

Transformata całki

L { 0 t f ( τ ) d τ } = 1 s F ( s ) {\displaystyle {\mathcal {L}}\left\{\int \limits _{0}^{t}f(\tau )d\tau \right\}={\frac {1}{s}}F(s)}
gdzie F ( s ) = L { f ( t } ( s ) {\displaystyle F(s)={\mathcal {L}}\{f(t\}(s)} jest transformatą oryginału

Oznacza to, że: Całkowanie oryginału f(t) w dziedzinie czasu można zastąpić transformatą z oryginału, podzieloną przez zmienną s w dziedzinie częstotliwości

Całka transformaty

L { f ( t ) t } = s F ( σ ) d σ {\displaystyle {\mathcal {L}}\left\{{\frac {f(t)}{t}}\right\}=\int \limits _{s}^{\infty }F(\sigma )\,d\sigma }

Przesunięcie w dziedzinie transformaty

L { e a t f ( t ) } = F ( s a ) {\displaystyle {\mathcal {L}}\{e^{at}f(t)\}=F(s-a)}
L 1 { F ( s a ) } = e a t f ( t ) {\displaystyle {\mathcal {L}}^{-1}\{F(s-a)\}=e^{at}f(t)}

Transformata funkcji z przesunięciem

L { f ( t a ) 1 ( t a ) } = e a s F ( s ) {\displaystyle {\mathcal {L}}\{f(t-a)1(t-a)\}=e^{-as}F(s)}
L 1 { e a s F ( s ) } = f ( t a ) 1 ( t a ) , {\displaystyle {\mathcal {L}}^{-1}\{e^{-as}F(s)\}=f(t-a)1(t-a),}
gdzie 1 ( t ) {\displaystyle 1(t)} oznacza skok jednostkowy.

Splot jednostronny

L { 0 t f ( u ) g ( t u ) d u } = L { f g } = L { f } L { g } {\displaystyle {\mathcal {L}}\left\{\int \limits _{0}^{t}f(u)\cdot g(t-u)\,du\right\}={\mathcal {L}}\{f*g\}={\mathcal {L}}\{f\}{\mathcal {L}}\{g\}}
Jest to tzw. twierdzenie Borela o splocie.

Transformata funkcji okresowej o okresie T

L { f } = 1 1 e s T 0 T e s t f ( t ) d t {\displaystyle {\mathcal {L}}\{f\}={\frac {1}{1-e^{-sT}}}\int \limits _{0}^{T}e^{-st}f(t)\,dt}

Własności graniczne

lim t 0 f ( t ) = lim s s F ( s ) {\displaystyle \lim _{t\to 0}f(t)=\lim _{s\to \infty }sF(s)}
lim t f ( t ) = lim s 0 s F ( s ) {\displaystyle \lim _{t\to \infty }f(t)=\lim _{s\to 0}sF(s)}

Tabela transformat Laplace’a

Poniżej zestawiono transformaty Laplace'a częściej spotykanych funkcji[7]. Zakładamy, że funkcje, na których dokonuje się transformacji Laplace'a są oryginałami (wg. definicji wyżej podanej), w szczególności f ( t ) = 0 {\displaystyle f(t)=0} dla t < 0 {\displaystyle t<0} , mimo że jawnie nie zostało to zapisane w poniższych wzorach. (Tylko jednak przy takim założeniu transformata odwrotna zastosowana do transformaty funkcji pierwotnej zwróci funkcję pierwotną.)

Oznaczenia: a {\displaystyle a} - stała liczba

L { a } = a 1 s {\displaystyle {\mathcal {L}}\{a\}=a{\frac {1}{s}}}
L { a t } = a 1 s 2 {\displaystyle {\mathcal {L}}\{at\}=a{\frac {1}{s^{2}}}}
L { a t n } = a n ! s n + 1 d l a n = 0 , 1 , 2 , 3 , {\displaystyle {\mathcal {L}}\{at^{n}\}=a{\frac {n!}{s^{n+1}}}\qquad dla\quad n=0,1,2,3,\dots }
L { e a t } = 1 s a {\displaystyle {\mathcal {L}}\{e^{at}\}={\frac {1}{s-a}}}
L { sin ( a t + b ) } = a cos b + s sin b s 2 + a 2 {\displaystyle {\mathcal {L}}\{\sin(at+b)\}={\frac {a\cdot \cos b+s\cdot \sin b}{s^{2}+a^{2}}}}
L { sin ( a t ) } = a s 2 + a 2 {\displaystyle {\mathcal {L}}\{\sin(at)\}={\frac {a}{s^{2}+a^{2}}}}
L { cos ( a t ) } = s s 2 + a 2 {\displaystyle {\mathcal {L}}\{\cos(at)\}={\frac {s}{s^{2}+a^{2}}}}
L { sinh ( a t + b ) } = 1 2 ( a s ) e b + 1 2 ( a + s ) e b s 2 a 2 {\displaystyle {\mathcal {L}}\{\sinh(at+b)\}={\frac {{\frac {1}{2}}(a-s)e^{-b}+{\frac {1}{2}}(a+s)e^{b}}{s^{2}-a^{2}}}}
L { sinh ( a t ) } = a s 2 a 2 {\displaystyle {\mathcal {L}}\{\sinh(at)\}={\frac {a}{s^{2}-a^{2}}}}
L { cosh ( a t ) } = s s 2 a 2 {\displaystyle {\mathcal {L}}\{\cosh(at)\}={\frac {s}{s^{2}-a^{2}}}}
L { t n e a t } = n ! ( s a ) n + 1 {\displaystyle {\mathcal {L}}\{t^{n}e^{at}\}={\frac {n!}{(s-a)^{n+1}}}}
L { e a t sin ( b t ) } = b ( s a ) 2 + b 2 {\displaystyle {\mathcal {L}}\{e^{at}\sin(bt)\}={\frac {b}{(s-a)^{2}+b^{2}}}}
L { t 2 b sin ( b t ) } = s ( s 2 + b 2 ) 2 {\displaystyle {\mathcal {L}}\left\{{\frac {t}{2b}}\sin(bt)\right\}={\frac {s}{(s^{2}+b^{2})^{2}}}}
L { e a t cos ( b t ) } = s a ( s a ) 2 + b 2 {\displaystyle {\mathcal {L}}\{e^{at}\cos(bt)\}={\frac {s-a}{(s-a)^{2}+b^{2}}}}
L { ln ( a t ) } = γ + ln ( s ) ln ( a ) s , {\displaystyle {\mathcal {L}}\{\ln(at)\}=-{\frac {\gamma +\ln(s)-\ln(a)}{s}},}

gdzie γ 0 , 5772156649 {\displaystyle \gamma \approx 0,5772156649} stała Eulera.

Transformata odwrotna Laplace’a

Transformatą odwrotną funkcji C s F ( s ) C {\displaystyle \mathbb {C} \ni s\to F(s)\in \mathbb {C} } nazywamy taką funkcję R + = ( 0 , + ) t f ( t ) R , {\displaystyle \mathbb {R_{+}} =(\,0\,,\,+\infty \,)\ni t\to f(t)\in \mathbb {R} ,} której transformatą jest F : {\displaystyle F{:}}

L 1 { F } = f , {\displaystyle {\mathcal {L}}^{-1}\{F\}=f,} jeżeli F = L { f } . {\displaystyle F={\mathcal {L}}\{f\}.}

Transformata odwrotna dana jest wzorem[1]

f ( t ) = L 1 { F } ( t ) = 1 2 π i t i t + i e s t F ( s ) d s {\displaystyle f(t)={\mathcal {L}}^{-1}\{\mathrm {F} \}(t)={\frac {1}{2\pi {\rm {i}}}}\int _{t-{\rm {i}}\cdot \infty }^{t+{\rm {i}}\cdot \infty }{\rm {e}}^{st}\mathrm {F} (s)\,{\rm {d}}s}

Rozwiązywanie równań różniczkowych za pomocą transformaty Laplace'a

Równania różniczkowe liniowe o stałych współczynnikach można rozwiązać wykorzystując transformatę Laplace'a.[8] Rozważmy równanie różniczkowe rzędu n {\displaystyle n} postaci

x ( n ) + a 1 x ( n 1 ) + a 2 x ( n 2 ) + a n 1 x ( 1 ) + a n x = f ( t ) {\displaystyle x^{(n)}+a_{1}x^{(n-1)}+a_{2}x^{(n-2)}+\dots a_{n-1}x^{(1)}+a_{n}x=f(t)}

gdzie:

  • x ( t ) {\displaystyle x(t)} - szukana funkcja zmiennej t {\displaystyle t} , określona przez warunki początkowe:
    x ( 0 ) = c 1 , d x d t ( 0 ) = c 2 , , d x n 1 d t x n 1 ( 0 ) = c n {\displaystyle x(0)=c_{1},\,\,{\frac {dx}{dt}}(0)=c_{2},\dots ,{\frac {dx^{n-1}}{dtx^{n-1}}}(0)=c_{n}}
  • x ( 1 ) , , x ( n 2 ) , x ( n 1 ) , x ( n ) {\displaystyle x^{(1)},\dots ,x^{(n-2)},x^{(n-1)},x^{(n)}} - pochodne funkcji x ( t ) {\displaystyle x(t)} po zmiennej t {\displaystyle t} rzędu 1 , , n 2 , n 1 , n {\displaystyle 1,\dots ,n-2,n-1,n}
  • a 1 , a 2 , , a n 1 , a n {\displaystyle a_{1},a_{2},\dots ,a_{n-1},a_{n}} - zadane stałe współczynniki
  • f ( t ) {\displaystyle f(t)} - zadana funkcja zmiennej t {\displaystyle t} , która posiada transformatę Laplace'a

Metoda szukania rozwiązania:

(1) Wykonujemy transformatę Laplace'a na obu stronach równania różniczkowego; korzystając z liniowości transformaty otrzymuje się:

L { x ( n ) } + a 1 L { x ( n 1 ) } + a 2 L { x ( n 2 ) } + a n 1 L { x ( 1 ) } + a n L { x } = L { f ( t ) } {\displaystyle {\mathcal {L}}\{x^{(n)}\}+a_{1}{\mathcal {L}}\{x^{(n-1)}\}+a_{2}{\mathcal {L}}\{x^{(n-2)}\}+\dots a_{n-1}{\mathcal {L}}\{x^{(1)}\}+a_{n}{\mathcal {L}}\{x\}={\mathcal {L}}\{f(t)\}}

(2) Korzystając ze wzorów na transformaty pochodnej funkcji poszczególnych rzędów po prostych rachunkach otrzymuje się:

L { x } = L { f ( t ) } + c 1 s n 1 + ( c 2 + a 1 c 1 ) s n 2 + + ( c n + a 1 c n 1 + + a n ) s 0 s n + a 1 s n 1 + a 2 s n 2 + + a n {\displaystyle {\mathcal {L}}\{x\}={\frac {{\mathcal {L}}\{f(t)\}+c_{1}s^{n-1}+(c_{2}+a_{1}c_{1})s^{n-2}+\dots +(c_{n}+a_{1}c_{n-1}+\dots +a_{n})s^{0}}{s^{n}+a_{1}s^{n-1}+a_{2}s^{n-2}+\dots +a_{n}}}}
- powyższy wzór przedstawia transformatę Laplace'a szukanej funkcji x ( t ) {\displaystyle x(t)} .

(3) Transformatę Laplace'a rozkłada się na ułamki proste i odczytuje się z tabeli transformat funkcje pierwotne poszczególnych ułamków - w ten sposób otrzymuje się rozwiązanie.

Uwaga:

Ostatni krok jest analogiczny do sposobu znanego z teorii całek: funkcję pierwotną znajduje się najprościej na podstawie tabeli, zestawiającej funkcje podcałkowe i odpowiadające im funkcje pierwotne.

Przykład

Znaleźć funkcję x ( t ) {\displaystyle x(t)} spełniającą równanie[9]

d x d t + x = sin t {\displaystyle {\frac {dx}{dt}}+x=\sin t}

z warunkiem początkowym x ( 0 ) = 0 {\displaystyle x(0)=0}

Rozwiązanie: Znajdujemy transformatę Laplace'a obu stron równania

L { x ( 1 ) } + L { x } = L { sin ( t ) } {\displaystyle {\mathcal {L}}\{x^{(1)}\}+{\mathcal {L}}\{x\}={\mathcal {L}}\{\sin(t)\}}

i korzystając z Tabeli transformat dla pochodnej funkcji i funkcji sin(t) otrzymamy:

s L { x } + L { x } = 1 s 2 + 1 {\displaystyle s{\mathcal {L}}\{x\}+{\mathcal {L}}\{x\}={\frac {1}{s^{2}+1}}}

Z powyższego równania obliczamy L { x } {\displaystyle {\mathcal {L}}\{x\}} , wykonując proste przekształcenia algebraiczne (i w tym tkwi siła tej metody, że równanie różniczkowe jest zamieniane z pomocą transformaty Laplace'a na proste równanie algebraiczne)

L { x } = 1 ( s + 1 ) ( s 2 + 1 ) {\displaystyle {\mathcal {L}}\{x\}={\frac {1}{(s+1)(s^{2}+1)}}}

Po rozbiciu na ułamki proste mamy:

L { x } = 1 2 1 s + 1 1 2 s s 2 + 1 + 1 2 1 s 2 + 1 {\displaystyle {\mathcal {L}}\{x\}={\frac {1}{2}}{\frac {1}{s+1}}-{\frac {1}{2}}{\frac {s}{s^{2}+1}}+{\frac {1}{2}}{\frac {1}{s^{2}+1}}}

Aby znaleźć funkcję pierwotną należałoby wykonać transformatę odwrotną Laplace'a na powyższym równaniu. Jednak najprościej jest skorzystać z tabeli transformat, gdzie odnajdujemy oryginały na podstawie transformat, jakimi są poszczególne ułamki proste. Otrzymujemy:

x ( t ) = 1 2 e t 1 2 cos t + 1 2 sin t {\displaystyle x(t)={\frac {1}{2}}e^{-t}-{\frac {1}{2}}\cos t+{\frac {1}{2}}\sin t}

Powyższe rozwiązanie spełnia warunki początkowe x ( 0 ) = 0 {\displaystyle x(0)=0} , co łatwo sprawdzić, podstawiając do rozwiązania t = 0. {\displaystyle t=0.}

Zastosowanie

 Osobny artykuł: Funkcja przejścia.

Transformata Laplace’a posiada kilka własności, które czynią ją szczególnie użyteczną w analizie liniowych układów dynamicznych. W inżynierii i fizyce jako narzędzie analizy graficznej wykorzystywana jest płaszczyzna S. Na płaszczyźnie S, mnożenie przez s {\displaystyle s} daje efekt różniczkowania (zob. człon różniczkujący), dzielenie przez s {\displaystyle s} daje efekt całkowania (zob. człon całkujący). Za pomocą transformaty Laplace'a można efektywnie rozwiązywać równania różniczkowe liniowe o stałych współczynnikach.[6]

Analiza pierwiastków zespolonych równania na płaszczyźnie S i przedstawienie ich na wykresie Arganda, może ujawnić informacje na temat charakterystyk częstotliwościowych i na temat stabilności układu (przebieg rzeczywistej funkcji czasu).

Kodowanie oznaczenia

W Unikodzie symbol transformaty Laplace’a ma postać:

Znak Unicode Kod HTML Nazwa unikodowa Nazwa polska
U+2112 &#x2112; lub &#8466; SCRIPT CAPITAL L pisana wielka litera L

W LaTeX-u używa się znacznika:

Znak LaTeX
L {\displaystyle {\mathcal {L}}} \mathcal L

Zobacz też

Przypisy

Bibliografia

  • Włodzimierz Krysicki, Lech Włodarski: Analiza matematyczna w zadaniach. Część II. Warszawa: Państwowe Wydawnictwo Naukowe, 1987, s. 378–385. ISBN 83-01-02440-2.
  • B. Piłat, M. J. Wasilewski, Tablice całek, Wydawnictwo Naukowo-Techniczne, Warszawa 1985, s. 98-102 (Tabela: Oryginały i transformaty Laplace'a) oraz s. 107 (Całkowanie równań liniowych metodą przekształcenia Laplace'a). ISBN 978-83-20-40-525-5
  • T. Trajdos, Matematyka dla inżynierów, PWN, Warszawa 1974, s. 567-608, s. 617-618. ISBN 978-83-20-40-152-3
  • W. Żakowski, W. Leksiński, Matematyka cz. IV, Wydawnictwo Naukowo-Techniczne, Warszawa 1978, s. 371-409. ISBN 978-83-01-19359-1

Linki zewnętrzne

  • Materiały dydaktyczne DSP AGH. dsp.agh.edu.pl. [zarchiwizowane z tego adresu (2013-12-17)].
  • What does the Laplace Transform really tell us? A visual explanation (plus applications) - wspaniałe objaśnienie transformacji Laplace'a, ilustrowane wykresami 3D
  • p
  • d
  • e
transformacje całkowe
inne transformacje
w rachunku prawdopodobieństwa
Kontrola autorytatywna (transformacja całkowa):
  • LCCN: sh85074666
  • NDL: 00567362
  • BnF: 119531733
  • BNCF: 35184
  • NKC: ph117703
  • J9U: 987007555498105171
  • Britannica: topic/Laplace-transform
  • SNL: Laplace-transformasjon