Think you've seen it all? Think again. Outside those doors, we might see anything. We could find new worlds, terrifying monsters, impossible things. And if you come with me... nothing will ever be the same again!
- The Doctor, Doctor Who
Látum f vera fall sem skilgreint er á menginu \mathbb{R}_+=\{t\in \mathbb{R}; t\geq 0\} með gildi í {\mathbb{C}} og gerum ráð fyrir að f sé heildanlegt á sérhverju lokuðu og takmörkuðu bili [0,b]. Laplace–mynd f, sem við táknum með {\cal L} f eða {\cal L}\{f\}, er skilgreind með formúlunni
{\cal L} f(s)=\int_0^ \infty e^{-st}f(t)\, dt.\label{7.1.1}
Skilgreiningarmengi fallsins {\cal L} f samanstendur af öllum tvinntölum s þannig að heildið í hægri hliðinni sé samleitið. Laplace-ummyndun er vörpunin {\cal L} sem úthlutar falli f Laplace-mynd sinni {\cal L} f.
Attention!
Stundum er skilgreind ,,tvíhliða Laplace-mynd‘‘ falls þar sem heildið er reiknað frá -\infty til \infty og þá gæti maður kallað okkar Laplace-mynd ,,einhliða Laplace-mynd‘‘.
Við segjum að fallið f:\mathbb{R}_+\to {\mathbb{C}} sé af veldisvísisgerð ef til eru jákvæðir fastar M og c þannig að
|f(t)|\leq Me^{c t}, \qquad t\in \mathbb{R}_+.\label{7.1.2}
Ef f er þar að auki heildanlegt á sérhverju takmörkuðu bili [0,b], þá er {\cal L} f skilgreint fyrir öll s\in {\mathbb{C}} með \operatorname{Re\, } s>c. Við fáum að auki vaxtartakmarkanir á {\cal L} f,
|{\cal L} f(s) |\leq
\int_0^\infty e^{-\operatorname{Re\, } st} Me^{c t} \, dt =
\dfrac M{\operatorname{Re\, }\, s-c}, \qquad \operatorname{Re\, }\, s>c.
Laplace-ummyndunin er línuleg vörpun, en það þýðir að
{\cal L}\{\alpha f+\beta g\}(s)=\alpha{\cal L}\{f\}(s)+\beta{\cal L}\{g\}(s)
ef f og g eru föll af veldisvísisgerð, \alpha og \beta eru tvinntölur og s\in {\mathbb{C}} liggur í skilgreiningarmengi fallanna {\cal L}\{f\} og {\cal L}\{g\}.
Gerum ráð fyrir að föllin f,g\in C(\mathbb{R}_+) séu bæði af veldisvísisgerð og að til sé fasti c þannig að
{\cal L} f(s)={\cal L} g(s), \qquad s\in {\mathbb{C}}, \quad \operatorname{Re\, } s\geq c.
Þá er f(t)=g(t) fyrir öll t\in \mathbb{R}_+.
- Ef \alpha\in \mathbb{R} og \alpha>-1, þá er
{\cal L}\{t^\alpha\}(s)=\dfrac {\Gamma(\alpha+1)}{s^{\alpha+1}},
og sér í lagi gildir ef \alpha er heiltala að
{\cal L}\{t^\alpha\}(s)=\dfrac {\alpha!}{s^{\alpha+1}}.
- Fyrir sérhvert \alpha\in {\mathbb{C}} gildir
{\cal L}\{e^{\alpha t}\}(s)=\dfrac{1}{s-\alpha}.
{\cal L}\{e^{\alpha t}f\}(s) = {\cal L}\{f\}(s-\alpha).
\begin{aligned}
{\cal L}\{\cos\beta t\}(s) &=
\frac 12 {\cal L}\{e^{i\beta t}\}(s) +\frac 12{\cal L}\{e^{-i\beta t}\}(s)\\
&=\frac 12\left[\dfrac 1{s-i\beta}+\dfrac 1{s+i\beta}\right]
=\dfrac s{s^2+\beta^2},\\
{\cal L}\{\sin\beta t\}(s) &=
\frac 1{2i}{\cal L}\{e^{i\beta t}\}(s) -\frac 1{2i}{\cal L}\{e^{-i\beta t}\}(s)\\
&=\frac 1{2i}\left[\dfrac 1{s-i\beta}-\dfrac 1{s+i\beta}\right]
=\dfrac {\beta}{s^2+\beta^2},\\
{\cal L}\{\cosh \beta t\}(s) &=
\frac 12 {\cal L}\{e^{\beta t}\}(s) +\frac 12{\cal L}\{e^{-\beta t}\}(s)\\
&=\frac 12\left[\dfrac 1{s-\beta}+\dfrac 1{s+\beta}\right]
=\dfrac s{s^2-\beta^2},\\
{\cal L}\{\sinh \beta t\}(s) &=
\frac 1{2}{\cal L}\{e^{\beta t}\}(s) -\frac 1{2}{\cal L}\{e^{-i\beta t}\}(s)\\
&=\frac 1{2}\left[\dfrac 1{s-\beta}-\dfrac 1{s+\beta}\right]
=\dfrac \beta{s^2-\beta^2},\\
{\cal L}\{e^{\alpha t}t^{\beta}\}(s)
&=\dfrac{\Gamma(\beta+1)}{(s-\alpha)^{\beta+1}},\\
{\cal L}\{e^{\alpha t}\cos \beta t\}(s)
&=\dfrac{s-\alpha}{(s-\alpha)^2+\beta^2},\\
{\cal L}\{e^{\alpha t}\sin \beta t\}(s)
&=\dfrac{\beta}{(s-\alpha)^2+\beta^2},\\
{\cal L}\{e^{\alpha t}\cosh \beta t\}(s)
&=\dfrac{s-\alpha}{(s-\alpha)^2-\beta^2},\\
{\cal L}\{e^{\alpha t}\sinh \beta t\}(s)
&=\dfrac{\beta}{(s-\alpha)^2-\beta^2}.\end{aligned}
Ef u\in C^ m(\mathbb{R}_+) og u, u', u'', \dots, u^{(m-1)}, eru af veldisvísisgerð, þá er {\cal L}\{u^{(m)}\}(s) skilgreint fyrir öll s\in {\mathbb{C}} með \operatorname{Re\, } s nógu stórt og
{\cal L}\{u^{(m)}\}(s)=s^
m{\cal L}\{u\}(s)-s^{m-1}u(0)-\cdots
-su^{(m-2)}(0)-u^{(m-1)}(0).
Sér í lagi gildir að ef U(s)={\cal L}\{u(t)\}(s), þá er
{\cal L}\{u'\}(s) = sU(s)-u(0),\qquad\mbox{og}\qquad
{\cal L}\{u''\}(s) =s^2U(s)-su(0)-u'(0).
Leysa á upphafsgildisverkefni af taginu:
a_mu^{(m)}+\cdots +a_1u'+a_0u=f(t), \qquad u(0)=b_0,\ u'(0)=b_1,\ldots,u^{(m-1)}(0)=b_{m-1}.
Skref 1: Reiknið Laplace-mynd hvorrar hliðar fyrir sig. Setning 10.1.7. gefur aðferð til að reikna Laplace-mynd vinstri hliðarinnar.
Skref 2: Notið jöfnuna sem kemur út úr Skrefi 1 til að fá formúlu fyrir U(s)={\cal L}\{u(t)\}(s).
Skref 3: Notið töflu eða andhverfa Laplace-ummyndun til að finna samfellt fall u(t) sem hefur U(s) sem fundið var í Skrefi 2 sem Laplace-mynd. Gætuð t.d. þurft að nota stofnbrotaliðun.
Skref 4: Fallið u(t) sem fannst í Skrefi 3 er lausn upphafsgildisverkefnisins.
Fyrir vigurgilt fall u(t)=\big(u_1(t), \ldots, u_m(t)\big) má skilgreina Laplace-mynd þannig að tekin er Laplace-mynd í hverju hniti fyrir sig,
{\cal L}\{u(t)\}(s)=\big({\cal L}\{u_1(t)\}(s), \ldots,{\cal L}\{u_m(t)\}(s)\big).
Laplace-mynd fylkjagilds falls er reiknuð á sama hátt.
Um sérhvert m\times m fylki A gildir
{\cal L}\{e^{tA}\}(s) = (sI-A)^{-1}.
Gerum ráð fyrir að P(D)=a_mD^ m+\cdots+a_1D+a_0 sé línulegur afleiðuvirki með fastastuðla. Látum g\in C^{\infty}(\mathbb{R}) vera fallið sem uppfyllir
P(D)g=0,\qquad \mbox{me\eth }
g(0)=g'(0)=\cdots=g^{(m-2)}(0)=0,\mbox{ og }
g^{(m-1)}(0)=1/a_m.
Þá er G(t,\tau)=g(t-\tau) Green-fall virkjans P(D).
Ef a er einhver punktur þá hefur upphafsgildisverkefnið
P(D)u=f(t),
með
u(a)=u'(a)=\cdots=u^{(m-1)}(a)=0,
ótvírætt ákvarðaða lausn u_p\in C^m(I) sem gefin er með formúlunni
u_p(t) = \int_a^ t G(t,\tau)f(\tau) \, d\tau, \qquad t\in I,
og G(t,\tau) er Green-fall virkjans P(D).
Með sama táknmál og hér að ofan gildir að
{\cal L}\{g\}(s)=\frac{1}{P(s)}.
Einnig gildir
{\cal L}\{u_p\}(s)={\cal L}\left\{\int_0^tg(t-\tau)f(\tau)\, d\tau\right\}(s)=
{\cal L}\{g\}(s){\cal L}\{f\}(s).
Ef f og g eru föll af veldisvísisgerð og heildanleg á sérhverju bili [0,b], þá er
{\cal L}\left\{\int_0^tf(t-\tau)g(\tau)\, d\tau\right\}(s)=
{\cal L}\{f\}(s){\cal L}\{g\}(s).
Ef f er af veldisvísisgerð og heildanlegt á sérhverju bili [0,b], þá er
{\cal L}\left\{\int_0^t f(\tau) \, d\tau\right\}(s) = \dfrac 1s
{\cal L}\{f\}(s).
Látum f:\mathbb{R}_+\to {\mathbb{C}} vera fall af veldisvísigerð sem er heildanlegt á sérhverju bili [0,b]. Þá er {\cal L}f fágað á menginu \{s\in {\mathbb{C}}\mid \operatorname{Re\, } s>c\} (þar sem c er fastinn úr 10.1.2.) og
\frac{d^k}{ds^k}{\cal L}\{f\}(s)=(-1)^k{\cal L}\{t^kf(t)\}(s),\qquad \operatorname{Re\, } s>c.
Heaviside-fallið H(x) er skilgreint þannig að H(t)=1 ef t\geq 0 og H(t)=0 ef t<0.
Fallið H_a(t)=H(t-a) er þannig að H_a(t)=1 ef t\geq a og H(t)=0 ef t<a.
Fyrir a\geq 0 er
{\cal L}\{H_a\}(s)=\frac{e^{-as}}{s}.
Látum f:\mathbb{R}_+\to {\mathbb{C}} vera fall af veldisvísisgerð. Þá gildir um sérhvert a\geq 0 að
{\cal L}\{H(t-a)f(t-a)\}(s) = e^{-as}{\cal L}\{f\}(s).
þar sem fallið t\mapsto H(t-a)f(t-a) tekur gildið 0 fyrir öll t<a.
Látum a vera gefna rauntölu. Skilgreinum fall f_\epsilon(t) þannig að f_\epsilon(t)=1/\epsilon ef a<t<a+\epsilon en f(t)=0 fyrir öll önnur gildi á t. Athugið að ef a>0 þá er
\int_0^\infty f_\epsilon(t)\,dt=1.
Skilgreinum nú \delta_a sem "markgildið" \lim_{\epsilon\to 0}f_\epsilon. Sérstaklega skilgreinum við Dirac delta fallið sem \delta=\delta_0. Athugið að \delta_a(t)=\delta(t-a).
Attention!
Athugið að \delta_a er ekki venjulegt fall heldur útvíkkuð gerð af falli sem er kölluð dreififall.
Hægt er að reikna heildi af \delta_a og við fáum að
\int_c^d \delta_a(t)\,dt=\left\{\begin{array}{ll}
1 & \mbox{ef }c\leq a\leq d,\\
0 & \mbox{annars.}
\end{array}\right.
Auðvelt er að sjá að ef c\leq a\leq d þá er \int_c^d f(t)\delta_a(t)\,dt=f(a). Sérstaklega gildir að
{\cal L}\{\delta_a\}(s)=e^{-as}.
Ef f:\mathbb{R}_+\to {\mathbb{C}} er samfellt deildanlegt á köflum og uppfyllir |f(t)|\leq Me^{ct}, t\in \mathbb{R}_+, þar sem M og c eru jákvæðir fastar, þá gildir um sérhvert \xi>c og sérhvert t>0 að
\begin{aligned}
\tfrac 12(f(t+)+f(t-)) &= \lim_{R\to +\infty} \dfrac 1{2\pi}
\int_{-R}^{+R}e^{(\xi+i\eta)t}{\cal L} f(\xi+i\eta) \, d\eta\\
&= \lim_{R\to +\infty} \dfrac 1{2\pi i}
\int_{\xi-iR}^{\xi+iR}e^{\zeta t}{\cal L} f(\zeta) \, d\zeta,\nonumber\end{aligned}
þar sem \int_{\xi-iR}^{\xi+iR} táknar að heildað sé eftir línustrikinu með upphafspunktinn \xi-iR og lokapunktinn \xi+iR. Ef {\cal L} f(\xi+i\eta) er í L^ 1(\mathbb{R}) sem fall af \eta, þá er f samfellt í t og
\begin{aligned}
f(t)&= \dfrac 1{2\pi}
\int_{-\infty}^{+\infty}e^{(\xi+i\eta)t}{\cal L} f(\xi+i\eta) \, d\eta\\
&= \dfrac 1{2\pi i}
\int_{\xi-i\infty}^{\xi+i\infty}e^{\zeta t}{\cal L} f(\zeta) \,
d\zeta,\nonumber\end{aligned}
þar sem \int_{\xi-i\infty}^{\xi+i\infty} táknar að heildað sé eftir línunni \{\xi+i\eta; \eta\in \mathbb{R}\} í stefnu vaxandi \eta.
Látum f:\mathbb{R}_+\to {\mathbb{C}} vera samfellt deildanlegt á köflum og af veldisvísisgerð, með |f(t)|\leq Me^{ct}, t>0, og gerum ráð fyrir að hægt sé að framlengja {\cal L} f yfir í fágað fall á {\mathbb{C}}\setminus A, þar sem A er endanlegt mengi. Ef \xi>c, M_r táknar hálfhringinn sem stikaður er með \gamma_r(\theta)=\xi+ire^{i\theta}, \theta\in [0,\pi] og
\max_{\zeta\in M_r}|{\cal L} f(\zeta)|\to 0, \qquad r\to +\infty,
þá er
\frac 12(f(t+)+f(t-))=
\sum_{\alpha\in A}\operatorname{Res}(e^{\zeta t}{\cal L} f(\zeta),\alpha).
Ef f er samfellt í t, þá gildir
f(t)= \sum_{\alpha\in A}\operatorname{Res}(e^{\zeta t}{\cal L} f(\zeta),\alpha).
Notum áfram táknmálið í 10.2.3. Táknum með A núllstöðvamengi margliðunnar P(\zeta). Þá er
g(t)= \sum\limits_{\alpha\in A}
\operatorname{Res}\bigg( \dfrac {e^{t\zeta}}{P(\zeta)},\alpha\bigg).
If it’s time to go, remember what you’re leaving. Remember the best.
- The Doctor, Doctor Who