משוואות דיפרנציאליות רגילות/פתרון בעזרת פונקצית גרין

פתרון מד"ר לינארית עם מקדמים קבועים ותנאי התחלה

בפרק על פתרון מד"ר ע"י התמרת לפלס ראינו כי הפתרון הכללי של המד"ר

   $\ L_{n}[f(t)]=c_{n}f^{(n)}(t)+c_{n-1}f^{(n-1)}(t)+\cdots +c_{1}f'(t)+c_{0}f(t)=r(t)$

נתון על-ידי:

   $\ f(t)={\mathcal {L}}^{-1}\left[{\frac {R(s)+\mathrm {I} _{n}(s)}{c_{n}s^{n}+c_{n-1}s^{n-1}+\cdots +c_{1}s+c_{0}}}\right]$

כאשר L_n הוא אופרטור דיפרנציאלי לינארי ו-I_n הוא פולינום תנאי התחלה.

פונקצית גרין המתאימה למד"ר לינארית כנ"ל מוגדרת על ידי:

   $\ G(t)={\mathcal {L}}^{-1}\left[{\frac {1}{c_{n}s^{n}+c_{n-1}s^{n-1}+\cdots +c_{1}s+c_{0}}}\right]$

כלומר פונקצית גרין היא:

1. התמרת לפלס הפוכה של אחד-חלקי הפולינום האופייני של המד"ר;
2. פונקצית התמסורת של המערכת;
3. פתרון המד"ר עבור האילוץ המסוים $\ r(t)=\delta (t)$, כאשר δ היא פונקצית הלם.

ניתן כעת לנסח את פתרון המד"ר במונחים של פונקצית גרין (המתאימה לאופרטור L):

   $\ F(s)={\mathcal {L}}[G(t)]\cdot [R(s)+I_{n}(s)]$

תנאי התחלה הומוגניים

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

   $\ f(t)=\int \limits _{0}^{t}r(\tau )G(t-\tau )\mathrm {d} \tau =r(t)*G(t)$

כלומר פתרון המד"ר הוא הקונבולוציה של פונק' גרין G עם איבר האילוץ r.

ראינו לעיל כי בהינתן המד"ר מסדר n, $\ Lf(t)=r(t)$, עם ת"ה הומוגניים, ניתן להסתפק בפתרון הבעיה:

   $\ {\begin{cases}LG(t)=\delta (t)\\G^{(k)}(0)=0,\ \forall 0\leq k<n\end{cases}}$

ולבצע קונבולוציה בין G לבין r.

חשוב להבחין כי השימוש בפונקצית הלם מכניס אי-רציפות לפתרון. נווכח בכך באמצעות הפעלת $\ \lim _{\epsilon \to 0^{+}}\int _{0}^{\epsilon }\mathrm {d} t$ על המשואה עבור G:

   $\ \lim _{\epsilon \to 0^{+}}\int _{0}^{\epsilon }LG\mathrm {d} t=\overbrace {\lim _{\epsilon \to 0^{+}}\int _{0}^{\epsilon }\delta (t)\mathrm {d} t} ^{\equiv 1}$

1. אגף ימין שווה ל-1 ע"פ ההגדרה של פונק' הלם.
2. מכיון ש-G חסומה, $\ \lim _{\epsilon \to 0^{+}}\int _{0}^{\epsilon }G(t)\mathrm {d} t=0$ וגם $\ \lim _{\epsilon \to 0^{+}}\int _{0}^{\epsilon }G'(t)\mathrm {d} t=0$.
3. לכן, בכתיב של סכום, מתקבל הקשר:

   $\ \sum _{k=2}^{n}\lim _{\epsilon \to 0^{+}}\int _{0}^{\epsilon }G^{(k)}(t)\mathrm {d} t\quad \Rightarrow \quad \sum _{k=1}^{n-1}\left.G^{(k)}(t)\right|_{0}^{0^{+}}=1$

   בדרך כלל אי-הרציפות תבוא לידי ביטוי בנגזרת ה-(n-1) בלבד.

סיכום

1. ניתן לחשב את פונקצית גרין של אופרטור דיפרנציאלי לינארי L על-ידי התמרת לפלס הפוכה של 1-חלקי הפולינום האופייני של האופרטור.
2. לאופרטורים שונים יהיו פונקציות גרין שונות.
3. אם נוכל לחשב את פונקצית גרין של האופרטור, נוכל לקבל את פתרון המד"ר המתאימה, באופן מיידי, באמצעות נוסחת הקונבולוציה.
4. בדרך כלל אי-הרציפות של פונקצית גרין תבוא לידי ביטוי בנגזרת ה-(n-1) בלבד.

תנאי התחלה אי-הומוגניים

פרק זה לוקה בחסר. אתם מוזמנים לתרום לוויקיספר ולהשלים אותו. ראו פירוט בדף השיחה.

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

פתרון מד"ר לינארית עם תנאי שפה הומוגניים

פרק זה לוקה בחסר. אתם מוזמנים לתרום לוויקיספר ולהשלים אותו. ראו פירוט בדף השיחה.

פתרון מד"ר לינארית עם תנאי שפה אי-הומוגניים

פרק זה לוקה בחסר. אתם מוזמנים לתרום לוויקיספר ולהשלים אותו. ראו פירוט בדף השיחה.

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

• אינדקס של פונקציות גרין למשואות שונות:
  • אונ' נברסקה