תורת הבקרה/התמרת לפלס

מתוך ויקיספר, אוסף ספרי הלימוד והמדריכים החופשי.

כללית, התמרת לפלס היא התמרה אינטגרלית לינארית אשר הופכת נגזרת למכפלה ואינטגרציה לחלוקה, וכל התמרה היא מצורה של מנת פולינומים מוכללים^[1].

כאשר מבצעים התמרת לפלס על מישור הזמן (במשתנה ‎‎t) עוברים למה שנקרא "מישור התדר" (במשתנה s) אשר נקרא גם "מישור לפלס". נהוג לכתוב פונקציות במישור הזמן באות קטנה ופונקציות במישור התדר באות גדולה.

[עריכה] שימוש בתורת הבקרה

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

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

[עריכה] הגדרה

$\ F(s) = \mathcal{L} \left\{f(t)\right\} =\int_0^{\infty} e^{-st} f(t)\,dt$

[עריכה] דוגמאות

(להשלים) ק'טרי'בע הנ'קרהנ רקקרערקקיימנ

[עריכה] מטלאב

syms a t
laplace(exp(a*t))

[עריכה] מייפל

y:=exp(a*t)
laplace(y(t),t,s)

[עריכה] התמרת לפלס הפוכה

$\ f(t)= \mathcal{L}^{-1} \left\{F(s)\right\} = {1\over 2\pi j} \int_{c-j\infty}^{c+j\infty} e^{st} F(s)\,ds$

כאשר $\ t>0,\ c>c_0$ ו-c₀ הוא הסדר המעריכי של f.

[עריכה] דוגמאות

(להשלים)

[עריכה] מטלאב

syms a s
ilaplace(1/(s+a))

[עריכה] מייפל

Y:=1/(s+a)
invlaplace(Y(s),s,t)

[עריכה] שיטות למציאת התמרת לפלס הפוכה

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

[עריכה] קטבים בודדים

כאשר F היא מנת פולינומים מצומצמים עם קטבים בודדים (ממשיים או מרוכבים), נפרק את F לשברים חלקיים:

$\ F(s)={B(s)\over A(s)}= \frac{r_1}{s-p_1}+\frac{r_2}{s-p_2}+...+\frac{r_n}{s-p_n}$

כאשר $\ r_i, p_j$ הם הקטבים והשאריות, בהתאמה. מתקבל:

$\ f(t)= \mathcal{L}^{-1} \left\{F(s)\right\}= u(t) \cdot \sum_{k=1}^n r_ke^{p_kt}$

כאשר u היא פונקצית מדרגה. מתכונות הפונקציות המרוכבות ניתן להראות:

$\ r_k= \frac{B(p_k)}{A'(s)|_{s=p_k}}$

זכרו כי נוסחה זו נכונה לקטבים בודדים בלבד!

[עריכה] קטבים מרוכבים צמודים

במקרה זה, יופיע עבור כל זוג קטבים כאלה ביטוי מהצורה:

$\ F_k(s)= \frac{r_k}{s-p_k}+\frac{r_{k+1}}{s-p_{k+1}}$

כאשר F_k הוא פרקציה של הפונקציה F. נניח כי:

$\ r_{k,k+1}=a\pm jb,\ p_{k,k+1}= \alpha\pm j\beta$

אז מתקבל:

$\ f_k(t)= r_ke^{p_kt}+ r_{k+1}e^{p_{k+1}t}= 2|r_k|e^{\alpha t} \cos(\beta t + \arg(r_k))$

כאשר השיוויון האחרון נובע מתכונות המספרים המרוכבים.

[עריכה] קוטב מסדר n

במקרה זה, מפרקים לשברים חלקיים:

$\ F_k(s)= \frac{B_k(s)}{(s-p_k)^n}= \frac{r_{k_1}}{s-p_k}+ \frac{r_{k_2}}{(s-p_k)^2}+...+ \frac{r_{k_n}}{(s-p_k)^n}$

כאשר F_k הוא פרקציה של הפונקציה F, ו-B_k היא פרקציה של הפונקציה B.

מכאן מבצעים את התמרת לפלס הפוכה באמצעות טבלאות.

[עריכה] מטלאב

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

B=[2 3]; A=[1 3 1];
[r,p,k]=residue(B,A)

$\ r=\binom{1}{1},\ p=\binom{-2.6180}{-0.382},\ k=[\ ]$

[עריכה] מייפל

פירוק לשברים חלקיים:

convert(1/(s^3(s+1)(s+2)),parfrac,s);

[עריכה] משפט הערך ההתחלתי

אם קיימת התמרת לפלס של הפונקציות $\ f(t), f'(t)$ אז מתקיים:

$\ \lim\limits_{t\to 0^+} f(t)= \lim\limits_{s\to\infty} sF(s)$

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

[עריכה] משפט הערך הסופי

אם קיימת התמרת לפלס של הפונקציות $\ f(t), f'(t)$ וכל הקטבים של F נמצאים בחצי המישור השמאלי, ולכל היותר ישנו קוטב בודד בראשית, אז מתקיים:

$\ \lim\limits_{t\to\infty} f(t)= \lim\limits_{s\to 0} sF(s)$

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

[עריכה] קישורים חיצוניים

ערך בוויקיפדיה: התמרת לפלס

[עריכה] הערות

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

[0] מנת פולינומים מוכללים, כלומר תיתכן חזקה שאינה מספר טבעי.

[1]

תורת הבקרה/התמרת לפלס

מתוך ויקיספר, אוסף ספרי הלימוד והמדריכים החופשי.

תוכן עניינים

[עריכה] שימוש בתורת הבקרה

[עריכה] הגדרה

[עריכה] דוגמאות

[עריכה] מטלאב

[עריכה] מייפל

[עריכה] התמרת לפלס הפוכה

[עריכה] דוגמאות

[עריכה] מטלאב

[עריכה] מייפל

[עריכה] שיטות למציאת התמרת לפלס הפוכה

[עריכה] קטבים בודדים

[עריכה] קטבים מרוכבים צמודים

[עריכה] קוטב מסדר n

[עריכה] מטלאב

[עריכה] מייפל

[עריכה] משפט הערך ההתחלתי

[עריכה] משפט הערך הסופי

[עריכה] קישורים חיצוניים

[עריכה] הערות

צפיות

כלים אישיים

ניווט

קהילה

חיפוש

תיבת כלים

הדפסה/ייצוא