התמרת לפלס

התמרת לפלס היא כלי מתמטי שהשימוש בו מקל מאוד על ניתוח ההתנהגות של מערכות ליניאריות ללא תלות בזמן, כגון מעגלים חשמליים ומערכות מכניות ואופטיות. ההתמרה קרויה על-שמו של פייר-סימון לפלס.

==הגדרה==
את התהליך בו מבצעים התמרת לפלס לפונקציה <math>\ f</math> מקובל לסמן <math>\ \mathcal{L}(f)</math>. אם <math>\ f(t)</math> היא פונקציה ממשית, נהוג לסמן את ההתמרה שלה ב-<math>\ F(s)</math>, והיא מוגדרת לפי האינטגרל המסוים <math>F(s)
 = \mathcal{L} \left\{f(t)\right\} =\int_{0}^\infty e^{-st} f(t)\,dt</math>.

התמרת לפלס היא התמרה אינטגרלית: פונקציה המקבלת פונקציה (ממשית או מרוכבת) ומחזירה פונקציה מאותו סוג. ההתמרה היא העתקה ליניארית ממרחב הפונקציות הממשיות לעצמו (למעשה, ההעתקה מוגדרת רק בתנאים מסוימים, כפי שיתואר בהמשך). ייחודה הוא בכך שהיא מקיימת את הזהות <math>\ \mathcal{L}(f')=s\cdot \mathcal{L}(f)-f(0)</math>, וכך היא הופכת נגזרת של פונקציה לכפל של הפונקציה במשתנה, תכונה המקלה על ניתוח של מערכות דינמיות ליניאריות הקבועות בזמן.

ההתמרה שימושית במיוחד בפתרון של משוואות דיפרנציאליות. בהנדסת חשמל מקובל לומר שהתמרת לפלס מעבירה ממישור הזמן למישור התדר, אף על פי שלמישור לפלס אין כלל משמעות של תדר. ההתמרה שמעבירה ממישור הזמן למישור התדר היא התמרת פורייה (למעשה, אם s הוא מדומה טהור, התמרת לפלס זהה להתמרת פורייה). מישור לפלס הוא מישור מדומה ללא משמעות פיזיקלית פשוטה והוא משמש לפתרון משוואות דיפרנציאליות.

==היסטוריה==
התמרות אינטגרליות הוצגו לראשונה על ידי לאונרד אוילר שפתר בעזרת ההתמרות משוואות דיפרנציאליות רגילות ליניאריות מסדר שני. לפלס עצמו מזכיר בכתביו את עבודתו של אוילר על התמרות אינטגרליות. היה זה שפיצר שהצמיד את שמו של לפלס לביטוי <math> F(s) = \int_a^b e^{st} f(t) dt </math> שאוילר השתמש בו.

בסוף המאה ה-19 הרחיבו המתמטיקאים פואנקרה, פינקרלה , אבל ופיקאר את התמרת לפלס גם למישור המרוכב ולשני משתנים. במאה ה-20 שכללו החוקרים הארי בייטמן , ברנשטיין, דוטש , הביסייד וברומיץ' את התמרת לפלס.

== הרחבה למספרים מרוכבים ==
כאמור, התמרת לפלס פועלת גם במישור המרוכב, ובמקום הדרישה <math>s>0</math> מגיעה הדרישה <math>Re(s)>0</math>. נראה זאת בדוגמה הבאה:
: <math>\mathcal{L}(e^{i\omega t} ) = \int_0^\infty {e^{ - st} e^{i\omega t} dt = } \lim_{M \to \infty } |_0^M = {1 \over {s - i\omega }}</math>
השוויון האחרון מתקיים מכיוון ש:
: <math>\lim_{M \to \infty } |e^{i\omega M} e^{-sM}| = \lim_{M \to \infty } e^{-Re(s)M} = 0 \quad, \quad Re(s) > 0 </math>
באופן דומה ניתן להראות כי: <math>\mathcal{L}(e^{ - i\omega t} ) = \frac {1}{s + i\omega}</math>.
לכן, על פי הליניאריות של ההתמרה נקבל כי:
: <math>\mathcal{L}(\sin \omega t) = \frac{\mathcal{L}(e^{i\omega t} ) - \mathcal{L}(e^{- i\omega t} )}{2i} = \mathcal{L}(\frac {e^{i\omega t} - e^{- i\omega t} }{2i}) = \frac {\omega}{s^2 + \omega ^2 }</math>, ובדומה: <math>\mathcal{L}(\cos \omega t) = \frac{s}{s^2 + \omega ^2 }</math>.

ניתן לראות שאם s הוא מדומה (כלומר, <math>Re(s)=0</math>) התמרת לפלס הופכת להתמרת פורייה.

==סדר מעריכי של פונקציה==
נשאלת השאלה לאילו פונקציות קיימת התמרת לפלס. על פי ההגדרה, אם האינטגרל מתכנס אז ההתמרה קיימת, אבל קיים גם מדד עבור הפונקציה עצמה – הסדר המעריכי שלה:

לפונקציה <math>\ f(t)</math> יש סדר מעריכי בגובה <math>\ \alpha</math> אם קיימים קבועים <math>\ M>0, \alpha, t_0 \ge 0</math> כך שלכל <math>\ t \ge t_0</math> מתקיים <math>\ |f(t)| \le M e^{\alpha t}</math>.

לדוגמה, הסדר המעריכי של <math>\ e^{at}</math> הוא <math>\ a</math> (<math>\ a</math> יכול להיות שלילי), הסדר המעריכי של <math>\ \sin {at}</math> הוא 0 (כי נוכל לבחור <math>\ M</math> מספיק גדול), ולפולינומים יש סדר מעריכי כלשהו גדול מ-0 (שוב, כי נוכל לבחור <math>\ M</math> מספיק גדול, ובאינסוף אקספוננט "מנצח" פולינום). ל־<math>\ e^{t^2}</math> לעומת זאת, לא קיים סדר מעריכי כי לא ניתן לחסום אותה עם אקספוננט.

ואז אם <math>\ f(t)</math> רציפה למקוטעין על חצי-הישר הימני, ומסדר מעריכי <math>\ \alpha</math> אז קיימת לה התמרת לפלס לכל <math>\ Re(s) > \alpha</math> (במקרה הממשי – <math>\ s>\alpha</math>). כלומר, האקספוננט <math>\ e^{-st}</math> צריך "למשוך למטה" חזק יותר. בעוד שלכל פונקציה המקיימת את התנאים האמורים קיימת התמרת לפלס, קיימות פונקציות שאינן מקיימות תנאים אלו, ובכל זאת קיימת להן התמרה. לדוגמה, לפונקציה <math>\ f(t) = t e^{t^2} \sin{e^{t^2}}</math> לא קיים סדר מעריכי אך קיימת לה התמרת לפלס לכל <math>\ Re(s) > 0</math>.

== תכונות ==
בהינתן שתי פונקציות, <math>\ f(t)</math> ו־<math>\ g(t)</math>, וההתמרות לפלס שלהן הן <math>\ F(s)</math> ו-<math>\ G(s)</math> בהתאמה,
: <math> F(s) = \mathcal{L} \{ f(t) \} </math>
: <math> G(s) = \mathcal{L} \{ g(t) \} </math>
מתקיימות התכונות הבאות:
* כללי:
** ניתן להראות שהתמרת לפלס של כל פונקציה רציפה למקוטעין ובעלת סדר מעריכי, מתכנסת בהחלט ובמידה שווה.
** אם <math>\ f(t)</math> רציפה למקוטעין על חצי הישר הימני ובעלת סדר מעריכי <math>\ \alpha</math> אז <math>\ \lim_{Re(s) \to \infty}F(s)=\lim_{Re(s) \to \infty}\mathcal{L}(f(t))=0</math>. רואים בטבלת ההתמרות כי מעלת המכנה תמיד גבוהה ממעלת המונה.
* ליניאריות:
: <math>\mathcal{L}\left\{a f(t) + b g(t) \right\}
 = a F(s) + b G(s) </math>
* גזירות:
: <math>\mathcal{L}\{f'\}
 = s \mathcal{L}\{f\} - f(0)</math>
: <math>\mathcal{L}\{f\}
 = s^2 \mathcal{L}\{f\} - s f(0) - f'(0)</math>
: <math>\mathcal{L}\left\{ f^{(n)} \right\}
 = s^n \mathcal{L}\{f\} - s^{n - 1} f(0) - \cdots - f^{(n - 1)}(0)</math>
* קונבולוציה:
: <math>\mathcal{L}\{f(t) * g(t)\}
 = F(s) \cdot G(s) </math>
* אינטגרציה וקונבולוציה:
: <math>\mathcal{L}\left\{ \int_0^t f(\tau)\, d\tau \right\}
 = \mathcal{L}\left\{ u(t) * f(t)\right\} = {1 \over s} F(s) </math>
* איפנון:
: <math>\mathcal{L}\left\{ e^{at} f(t) \right\}
 = F(s - a)</math>
: <math>\mathcal{L}^{-1} \left\{ F(s - a) \right\}
 = e^{at} f(t)</math>
* הזזה:
: <math>\mathcal{L}\left\{ f(t - a) u(t - a) \right\}
 = e^{-as} F(s)</math>
: <math>\mathcal{L}^{-1} \left\{ e^{-as} F(s) \right\}
 = f(t - a) u(t - a)</math>
: כאן <math>\ u(t)</math> היא פונקציית מדרגה.

== התמרת לפלס של טורי חזקות ==
נגדיר טור חזקות: <math>\ f(t)=\sum_{k=0}^n a_kt^k</math>. על ידי שימוש בתכונת הליניאריות ובטבלת ההתמרות, נקבל: <math>\ \mathcal{L}(f(t))=\sum_{k=0}^n a_k\mathcal{L}(t^k)=\sum_{k=0}^n \frac{a_kk!}{s^{k+1}}</math>. לטורי חזקות אינסופיים, לעומת זאת, לא תמיד קיימת התמרת לפלס. התמרת לפלס עבור טור החזקות <math>\ f(t)=\sum_{k=0}^\infty a_kt^k</math> קיימת אם הטור מתכנס לכל <math>\ t \ge 0</math>, וקיים <math>\ N</math> טבעי כך שלכל <math>\ n \ge N</math> ולכל <math>\ \alpha >0, M>0</math> מתקיים <math>\ |a_k|\le \frac{M\alpha^k}{k!}</math>. במקרה זה נקבל: <math>\ \mathcal{L}(f(t))=\sum_{k=0}^\infty a_k\mathcal{L}(t^k)=\sum_{k=0}^\infty \frac{a_kk!}{s^{k+1}}</math>.

==משפט הערך ההתחלתי==

אם קיימת התמרת לפלס של הפונקציות <math>\ f(t), f'(t)</math> אז מתקיים:
: <math>\ \lim\limits_{t\to 0^+} f(t)= \lim\limits_{s\to\infty} sF(s)</math>
משפט זה נוח במיוחד כאשר מעוניינים בערך ההתחלתי של הפונקציה במישור הזמן, משום שהוא אינו מחייב את חישוב של התמרת לפלס ההפוכה של <math>F</math>.

==משפט הערך הסופי==

אם קיימת התמרת לפלס של הפונקציות <math>\ f(t), f'(t)</math> וכל הקטבים של <math>F</math> נמצאים בחצי המישור השמאלי, ולכל היותר ישנו קוטב בודד בראשית, אז מתקיים:
: <math>\ \lim\limits_{t\to\infty} f(t)= \lim\limits_{s\to 0} sF(s)</math>
משפט זה שימושי כאשר מעוניינים בהתנהגות המערכת לאחר ריסון תנאי ההתחלה.

== טבלת התמרות לפלס שימושיות ==
{| width="720" border="2" cellspacing="0" cellpadding="6"
|- style="background-color: #AAEECC; "
! מס'
! align="center"| פונקציה
! align="center"| מרחב הזמן <br /><math>x(t) = \mathcal{L}^{-1} \left\{ X(s) \right\}</math>
! align="center"| מרחב התדר <br /><math>X(s) = \mathcal{L}\left\{ x(t) \right\}</math>
! align="center"| תחום התכנסות<br />למערכת סיבתית
|-
| 1
| פונקציית הלם
| align="center"|<math> \delta(t-\alpha) \ </math>, <math> \delta(t) \ </math>
| align="center"|<math> e^{-\alpha s} </math>, <math> 1 \ </math>
| align="center"|לכל <math> s \ </math>
|-
| 2
| פונקציית מדרגה
| align="center"|<math> 1 \cdot u(t-\alpha) \ </math>
| align="center"|<math> { e^{-\alpha s} \over s } </math>
| align="center"| <math> \mathrm{Re}(s) > 0 \, </math>
|-
| 3
| דעיכה מעריכית
| align="center"|<math> e^{-\alpha t} \cdot u(t) \ </math>
| align="center"|<math> { 1 \over s+\alpha } </math>
| align="center"| <math> \mathrm{Re}(s) > - \alpha \ </math>
|-
| 3
| דעיכה מעריכית מוכפלת ב־t
| align="center"|<math> te^{-\alpha t} \cdot u(t) \ </math>
| align="center"|<math> { 1 \over (s+\alpha)^2 } </math>
| align="center"| <math> \mathrm{Re}(s) > - \alpha \ </math>
|-
| 4
| קירוב מעריכי
| align="center"|<math>( 1-e^{-\alpha t}) \cdot u(t) \ </math>
| align="center"|<math>\frac{\alpha}{s(s+\alpha)} </math>
| align="center"| <math> \mathrm{Re}(s) > 0\ </math>
|-
| 5
| סינוס
| align="center"| <math> \sin(\omega t) \cdot u(t) \ </math>
| align="center"| <math> { \omega \over s^2 + \omega^2 } </math>
| align="center"| <math> \mathrm{Re}(s) > 0 \ </math>
|-
| 6
| קוסינוס
| align="center"| <math> \cos(\omega t) \cdot u(t) \ </math>
| align="center"| <math> { s \over s^2 + \omega^2 } </math>
| align="center"| <math> \mathrm{Re}(s) > 0 \ </math>
|-
| 7
| סינוס היפרבולי
| align="center"| <math> \sinh(\alpha t) \cdot u(t) \ </math>
| align="center"| <math> { \alpha \over s^2 - \alpha^2 } </math>
| align="center"| <math> \mathrm{Re}(s) > | \alpha | \ </math>
|-
| 8
| קוסינוס היפרבולי
| align="center"| <math> \cosh(\alpha t) \cdot u(t) \ </math>
| align="center"| <math> { s \over s^2 - \alpha^2 } </math>
| align="center"| <math> \mathrm{Re}(s) > | \alpha | \ </math>
|-
| 9
| גל סינוס בדעיכה מעריכית
| align="center"| <math>e^{-\alpha t} \sin(\omega t) \cdot u(t) \ </math>
| align="center"| <math> { \omega \over (s+\alpha )^2 + \omega^2 } </math>
| align="center"| <math> \mathrm{Re}(s) > -\alpha \ </math>
|-
| 10
| גל קוסינוס בדעיכה מעריכית
| align="center"| <math>e^{-\alpha t} \cos(\omega t) \cdot u(t) \ </math>
| align="center"| <math> { s+\alpha \over (s+\alpha )^2 + \omega^2 } </math>
| align="center"| <math> \mathrm{Re}(s) > -\alpha \ </math>
|-
| 11
| חזקת n
| align="center"| <math>t^n \cdot u(t) </math>
| align="center"| <math> { n! \over s^{n+1} } </math>
| align="center"| <math> \mathrm{Re}(s) > 0 \, </math>
|-
| 12
| השורש ה-n-י
| align="center"| <math> \sqrt[n]{t} \cdot u(t) </math>
| align="center"| <math> s^{-(n+1)/n} \cdot \Gamma\left(1+\frac{1}{n}\right)</math>
| align="center"| <math> \mathrm{Re}(s) > 0 \, </math>
|-
| 13
| לוגריתם טבעי
| align="center"| <math> \ln \left ( { t \over t_0 } \right ) \cdot u(t) </math>
| align="center"| <math> - { t_0 \over s} \ [ \ \ln(t_0 s)+\gamma \ ] </math>
| align="center"| <math> \mathrm{Re}(s) > 0 \, </math>
|-
| 14
| פונקציית בסל <br /> מהסוג הראשון, <br /> מסדר n
| align="center"| <math> J_n( \omega t) \cdot u(t)</math>
| align="center"|<math>\frac{ \omega^n \left(s+\sqrt{s^2+ \omega^2}\right)^{-n}}{\sqrt{s^2 + \omega^2}}</math>
| align="center"| <math> \mathrm{Re}(s) > 0 \, </math> <br /> <math> (n > -1) \, </math>
|-
| 15
| פונקציית בסל מתואמת <br /> מהסוג הראשון, <br /> מסדר n
| align="center"| <math>I_n(\omega t) \cdot u(t)</math>
| align="center"| <math> \frac{ \omega^n \left(s+\sqrt{s^2-\omega^2}\right)^{-n}}{\sqrt{s^2-\omega^2}} </math>
| align="center"| <math> \mathrm{Re}(s) > | \omega | \, </math>
|-
| 16
| פונקציית בסל מהסוג השני, <br /> מסדר 0
| align="center"| <math> Y_0(\alpha t) \cdot u(t)</math>
| align="center"|
| align="center"|
|-
| 17
| פונקציית בסל מתואמת מהסוג השני, <br /> מסדר 0
| align="center"| <math> K_0(\alpha t) \cdot u(t)</math>
| align="center"|
| align="center"|
|-
| 18
| פונקציית השגיאה
| align="center"| <math> \mathrm{erf}(t) \cdot u(t) </math>
| align="center"| <math> {e^{s^2/4} \operatorname{erfc} \left(s/2\right) \over s}</math>
| align="center"| <math> \mathrm{Re}(s) > 0 \, </math>
|-
| colspan="5"|
{|
|-
| colspan="2" | הערות:
|-
| valign="top"|
* <math> u(t) \, </math> מייצג את פונקציית המדרגה.
* <math> \delta(t) \, </math> מייצג את פונקציית דלתא של דיראק.
* <math> \Gamma (z) \, </math> מייצג פונקציית גמא.
* <math> \gamma \, </math> הוא קבוע אוילר.
| valign="top"|
* <math>t \, </math>, מספר ממשי, בדרך-כלל מייצג זמן, <br />אף על פי שיכול לייצג כל מימד בלתי-תלוי.
* <math>s \, </math> הוא מספר מרוכב המסמל את התדירות הזוויתית.
* <math> \beta \,</math>, ו־<math>\omega \,</math> הם מספרים ממשיים.
* <math> n \, </math> הוא מספר שלם.
|-
| colspan="2"|
* מערכת סיבתית היא מערכת שבה התגובה להלם <math>h(t)</math> היא אפס לכל זמן <math>t<0</math> (כלומר <math>h(t)=h(t) \cdot u(t)</math>).
|}
|}

== שימושים ==
* בפתרון משוואות דיפרנציאליות בעזרת התמרת לפלס, מנצלים את הזהות <math>\ \mathcal{L}(f')=s\cdot \mathcal{L}(f)-f(0)</math>, שמקשרת בין התמרת לפלס של הנגזרת לבין התמרת לפלס של הפונקציה. בהינתן משוואה דיפרנציאלית רגילה הכוללת נגזרות מסדרים שונים, ניתן לבצע התמרת לפלס על המשוואה, לבודד את הביטוי להתמרת לפלס של הפונקציה הנעלמת – ואז לבצע התמרת לפלס הפוכה ולמצוא את פתרון המשוואה.
* בתורת הבקרה, כאשר מאפיינים מערכת, מופיעים לעיתים קרובות ביטויים המערבים נגזרות. מאחר שהתמרת לפלס הופכת גזירה למכפלה ואינטגרציה לחלוקה, הטיפול בביטויים כאלה הוא נוח, וכאשר מאפייני המערכת אינם תלויים בזמן, התמרת לפלס של משוואות התנועה נותנת פולינומים, אשר מאפשרים אנליזה קלה של המצב (ע"ע: פונקציית תמסורת).

== ראו גם ==
* התמרת פורייה

== לקריאה נוספת ==
* Schiff, Joel L., The Laplace transform: theory and applications, Springer-Verlag, 1999, New York.

==קישורים חיצוניים==

* [http://wims.unice.fr/wims/wims.cgi?session=2KADF698D7.2&+lang=en&+module=tool%2Fanalysis%2Ffourierlaplace.cn מחשבון להתמרת לפלס]
* 
* 
* 

קטגוריה:אנליזה פונקציונלית
קטגוריה:התמרת פורייה