פונקציית ליאפונוב

פרק 1

הרעיון הגדול: להוכיח יציבות בלי לפתור את המערכת

עד כה, כשרצינו לדעת אם מערכת יציבה, הדרך הייתה "לפתור ולראות": במערכת ליניארית $\dot x=Ax$ חישבנו ערכים עצמיים ובדקנו אם כולם בחצי המישור השמאלי. אך שיטה זו נשענת על כך שאנו יודעים לפתור — או לפחות לאפיין — את הפתרון במפורש. במערכות לא־ליניאריות, $\dot x=f(x)$ , כמעט אף פעם אין פתרון מפורש. האם זה אומר שאיננו יכולים לומר דבר על יציבות?

התשובה של ליאפונוב: אפשר, ובאלגנטיות. במקום לשאול "לאן הולך המסלול?", נשאל "מה קורה לגודל סקלרי חכם אחד לאורך המסלול?". אם נצליח למצוא פונקציה $$V(x)$$ שמתנהגת כמו אנרגיה — חיובית בכל מקום פרט לשיווי המשקל, ורק יורדת בזמן — אז המסלול, יהיה אשר יהיה, מוכרח לשקוע אל נקודת המינימום של $$V$$ , כלומר אל שיווי המשקל. לא פתרנו שום דבר, ובכל זאת הוכחנו יציבות.

הרעיון בשורה אחת

אם קיימת "אנרגיה מוכללת" $$V(x)>0$$ שאינה עולה לעולם לאורך מסלולי המערכת ( $\dot V\leq 0$ ) — המערכת יציבה. אם היא ממש יורדת ( $\dot V<0$ ) — המסלולים מתכנסים לשיווי המשקל. פונקציה כזו נקראת פונקציית ליאפונוב, והשיטה כולה נקראת השיטה הישירה של ליאפונוב (direct method / second method) — "ישירה" כי היא עוקפת לחלוטין את הצורך בפתרון.

למה זה כלי כה מרכזי בבקרה מודרנית?

משלוש סיבות. ראשית, זהו הכלי הכללי היחיד ליציבות לא־ליניארית — ערכים עצמיים פשוט לא מוגדרים למערכת לא־ליניארית. שנית, גם בעולם הליניארי, גישת ליאפונוב מולידה את משוואת ליאפונוב — כלי אלגברי שבודק יציבות בלי לחשב ערכים עצמיים, ומופיע בכל פינה של הבקרה המודרנית (גרמיאנים, LQR, תצפתנים, LMI). שלישית, פונקציות ליאפונוב אינן רק כלי אנליזה אלא גם כלי תכן: שיטות שלמות של תכנון בקרים (backstepping, בקרה אדפטיבית) בנויות סביב בנייה מכוונת של פונקציית ליאפונוב.

מטרת העמוד

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

פרק 2

מפת מושגים: "ליאפונוב" בכל ההקשרים

אלכסנדר מיכאילוביץ' ליאפונוב (1857–1918), מתמטיקאי רוסי, פרסם ב־1892 את עבודת הדוקטורט שלו "הבעיה הכללית של יציבות התנועה" — ומאז שמו מוצמד לחצי תריסר מושגים שונים. סטודנטים רבים שומעים "ליאפונוב" בהרצאה ואינם בטוחים לאיזה ליאפונוב הכוונה. הנה המפה המלאה:

המושג	מה הוא	איפה פוגשים אותו
יציבות במובן ליאפונוב	ההגדרה המדויקת של יציבות: מסלול שמתחיל קרוב — נשאר קרוב (ε–δ)	הבסיס לכל דיון ביציבות; פרק 4 כאן
פונקציית ליאפונוב	פונקציה סקלרית $$V(x)$$ דמוית־אנרגיה המשמשת להוכחת יציבות	נושא העמוד הזה
השיטה הישירה (השנייה) של ליאפונוב	המשפטים שמסיקים יציבות מקיום פונקציית ליאפונוב	פרק 8 כאן
השיטה העקיפה (הראשונה) של ליאפונוב	הסקת יציבות מקומית של מערכת לא־ליניארית מתוך הליניאריזציה שלה	משפט הליניאריזציה; קורסי בקרה לא־ליניארית
משוואת ליאפונוב	המשוואה המטריצית $A^{\mathsf T}P+PA=-Q$ למערכות ליניאריות	עמוד נפרד; מוזכרת בפרק 11 כאן
מעריך (אקספוננט) ליאפונוב	קצב ההתבדרות של מסלולים סמוכים; מדד לכאוס	מערכות דינמיות וכאוס — לא בקורס בקרה סטנדרטי

איך לא ללכת לאיבוד

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

פרק 3

נקודות שיווי משקל — האובייקט שאת יציבותו בודקים

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

\dot x=f(x),\qquad x\in\mathbb{R}^n

נקודה $$x_e$$ נקראת נקודת שיווי משקל אם $$f(x_e)=0$$ — מערכת שמתחילה שם, נשארת שם לנצח (הנגזרת אפס). מערכת לא־ליניארית יכולה להחזיק כמה נקודות שיווי משקל, וכל אחת עשויה להיות יציבה או לא־יציבה בפני עצמה.

הנחה סטנדרטית: שיווי המשקל בראשית

בלי הגבלת הכלליות נניח תמיד ש־ $$x_e=0$$ . אם שיווי המשקל נמצא במקום אחר, מגדירים משתנה מוזז $$z=x-x_e$$ ומקבלים מערכת חדשה שבה שיווי המשקל הוא בראשית. זה חוסך סימונים ומאפשר לנסח את כל המשפטים ביחס ל־ $$x=0$$ .

דוגמה מוכרת: המטוטלת $\ddot\theta=-\frac{g}{\ell}\sin\theta$ מחזיקה שתי נקודות שיווי משקל מהותיות — $\theta=0$ (תלויה למטה) ו־ $\theta=\pi$ (עומדת על הראש). האינטואיציה ברורה: הראשונה יציבה, השנייה לא. תורת ליאפונוב תיתן לנו את הכלים להוכיח זאת — ולהכליל לכל מערכת.

פרק 4

הגדרות היציבות במובן ליאפונוב

לפני שנוכיח יציבות, חייבים להגדיר אותה במדויק. ליאפונוב נתן את ההגדרה שהפכה לסטנדרט — בסגנון $\varepsilon$ – $\delta$ המוכר מחדו״א. יש שלוש דרגות עיקריות, מהחלשה לחזקה:

דרגה 1 — יציבות (במובן ליאפונוב)

שיווי המשקל $$x=0$$ נקרא יציב אם לכל $\varepsilon>0$ קיים $\delta>0$ כך שכל מסלול שמתחיל במרחק קטן מ־ $\delta$ נשאר לעד במרחק קטן מ־ $\varepsilon$ :

\|x(0)\|<\delta\;\Longrightarrow\;\|x(t)\|<\varepsilon\quad\forall t\geq 0

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

דרגה 2 — יציבות אסימפטוטית

שיווי המשקל נקרא יציב אסימפטוטית אם הוא יציב, וגם מסלולים שמתחילים מספיק קרוב מתכנסים אליו:

\|x(0)\|<\delta_0\;\Longrightarrow\;\lim_{t\to\infty}x(t)=0

במילים: "מתחיל קרוב — נשאר קרוב וגם חוזר הביתה". זו הדרישה השימושית ברוב יישומי הבקרה.

דרגה 3 — יציבות אסימפטוטית גלובלית

אם ההתכנסות מתקיימת מכל תנאי התחלה ( $\delta_0=\infty$ ), היציבות נקראת גלובלית (GAS — Globally Asymptotically Stable). במילים: לא משנה מאיפה מתחילים — תמיד חוזרים הביתה.

עדינות חשובה — התכנסות לבדה אינה יציבות!

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

ובצד השני — אי־יציבות

שיווי משקל שאינו יציב נקרא לא־יציב: קיים $\varepsilon$ כך שקרוב ככל שנתחיל, יש מסלולים שבורחים מעבר לו. המטוטלת ההפוכה ( $\theta=\pi$ ) היא הדוגמה הקלאסית: כל סטייה זעירה גדלה ומפילה אותה.

דרגה	נשאר קרוב?	מתכנס?	מכל מקום?	דוגמה
יציבות	כן	לא בהכרח	—	מטוטלת ללא חיכוך
יציבות אסימפטוטית	כן	כן (מקומית)	לא בהכרח	מטוטלת עם חיכוך
יציבות גלובלית (GAS)	כן	כן	כן	$\dot x=-x$
אי־יציבות	לא	—	—	מטוטלת הפוכה

פרק 5

האינטואיציה: כדור, קערה ואנרגיה

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

האנרגיה חיובית תמיד

האנרגיה הכוללת (קינטית + פוטנציאלית) חיובית בכל מצב, ומתאפסת רק בנקודה אחת — הכדור נח בתחתית. לאנרגיה יש "מינימום גלובלי" יחיד בדיוק בשיווי המשקל.

החיכוך רק גורע ממנה

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

ולכן — הכדור שוקע לתחתית

כמות חיובית שרק יורדת חייבת להתכנס. ואם היא מתכנסת לאפס, המצב חייב להתכנס לנקודה היחידה שבה האנרגיה אפס — התחתית. הוכחנו התכנסות בלי לדעת דבר על צורת המסלול.

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

תמונת קווי הגובה

דרך מצוינת לדמיין: משטחי הרמה $$V(x)=c$$ הם "קווי גובה" סביב שיווי המשקל, כמו במפה טופוגרפית של עמק. התנאי $\dot V\le 0$ אומר שהמסלול לעולם אינו חוצה קו גובה כלפי מעלה — הוא יכול רק לרדת מקו לקו או להישאר עליו. מסלול שכלוא בתוך קו גובה סגור לא יכול לברוח החוצה — וזו בדיוק יציבות. אם הוא ממש יורד כל הזמן ( $\dot V<0$ ) — הוא שוקע אל תחתית העמק.

פרק 6

הגדרות: פונקציות מוגדרות חיובית ושלילית

כדי לנסח משפטים מדויקים, צריך שפה מדויקת ל"פונקציה שמתנהגת כמו אנרגיה". נניח $V:\mathbb{R}^n\to\mathbb{R}$ רציפה וגזירה ברציפות, עם $$V(0)=0$$ , ותהי $$D$$ סביבה של הראשית.

מוגדרת חיובית (PD)

$$V(0)=0$$ וגם $$V(x)>0$$ לכל $x\in D,\;x\neq 0$ . "גבעה" שמינימומה היחיד בראשית.

חצי־מוגדרת חיובית (PSD)

$$V(0)=0$$ וגם $V(x)\geq 0$ . מותר לה להתאפס גם מחוץ לראשית — יש לה "עמקים שטוחים".

מוגדרת / חצי־מוגדרת שלילית

$$V$$ מוגדרת (חצי־מוגדרת) שלילית אם $$-V$$ מוגדרת (חצי־מוגדרת) חיובית.

דוגמאות ממחישות ב־ $\mathbb{R}^2$

פונקציה	סיווג	למה
$$V=x_1^2+x_2^2$$	מוגדרת חיובית	מתאפסת רק בראשית
$$V=x_1^2$$	חצי־מוגדרת חיובית	מתאפסת על כל הציר $$x_1=0$$
$$V=x_1^2-x_2^2$$	לא מוגדרת (indefinite)	חיובית בכיוונים מסוימים, שלילית באחרים — "אוכף"
$$V=x_1^2+x_2^4$$	מוגדרת חיובית	לא כל PD היא ריבועית!

המקרה הריבועי — הקשר למטריצות

המשפחה החשובה ביותר של מועמדות היא הפונקציות הריבועיות:

V(x)=x^{\mathsf T}P\,x,\qquad P=P^{\mathsf T}

עבורן, מוגדרוּת הפונקציה שקולה בדיוק למוגדרוּת המטריצה: $$V$$ מוגדרת חיובית אם ורק אם $P\succ 0$ (כל הערכים העצמיים של $$P$$ חיוביים). בדיקה מעשית: קריטריון סילבסטר — כל המינורים הראשיים המובילים חיוביים. עבור $2\times2$ : $p_{11}>0$ וגם $\det P>0$ . הפונקציות הריבועיות יובילו אותנו בפרק 11 היישר אל משוואת ליאפונוב.

פרק 7

הנגזרת לאורך מסלולים — $\dot V$

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

$$V(x)$$ היא פונקציה של המצב. כשהמצב נע לפי $\dot x=f(x)$ , הערך $$V(x(t))$$ הופך לפונקציה של הזמן. נגזרתו מחושבת בכלל השרשרת:

\dot V(x)=\frac{d}{dt}V\bigl(x(t)\bigr) =\nabla V(x)\cdot\dot x =\sum_{i=1}^{n}\frac{\partial V}{\partial x_i}\,f_i(x) =\frac{\partial V}{\partial x}\,f(x)

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

הביטוי $\dot V=\nabla V\cdot f(x)$ תלוי רק בשני דברים שיש לנו ביד: הפונקציה $$V$$ שאנחנו בחרנו, ושדה הווקטורים $$f$$ שמגדיר את המערכת. הפתרון $$x(t)$$ אינו מופיע בשום מקום. אנחנו מחשבים ביטוי אלגברי במשתני המצב, ובודקים את סימנו כפונקציה של $$x$$ . זו כל הפואנטה של "השיטה הישירה" — היא עוקפת את הצורך לפתור את המד״ר.

הפירוש הגאומטרי

הגרדיאנט $\nabla V$ ניצב לקו הגובה של $$V$$ ומצביע "במעלה ההר". המכפלה $\nabla V\cdot f$ היא ההיטל של וקטור המהירות על כיוון העלייה. אם המכפלה שלילית — וקטור המהירות מצביע "במורד", כלומר פנימה דרך קו הגובה. מסלול שתמיד חוצה קווי גובה פנימה לא יכול אלא לשקוע אל המרכז. אם המכפלה אפס — המסלול משיק לקו הגובה, גולש עליו בלי לרדת.

פרק 8

השיטה הישירה: משפטי היציבות של ליאפונוב

כעת יש בידינו את כל המרכיבים, ואפשר לנסח את המשפטים המרכזיים. תהי $\dot x=f(x)$ עם שיווי משקל בראשית, ותהי $$V(x)$$ גזירה ברציפות בסביבה $$D$$ של הראשית.

משפט 1 — יציבות

אם בסביבה $$D$$ מתקיים:

V(x)>0\;\;\forall x\neq 0,\quad V(0)=0, \qquad \dot V(x)\leq 0\;\;\forall x\in D

אז שיווי המשקל $$x=0$$ יציב, ו־ $$V$$ נקראת פונקציית ליאפונוב של המערכת.

משפט 2 — יציבות אסימפטוטית

אם בנוסף הנגזרת מוגדרת שלילית (שלילית ממש מחוץ לראשית):

\dot V(x)<0\;\;\forall x\in D,\;x\neq 0

אז שיווי המשקל יציב אסימפטוטית: מסלולים שמתחילים מספיק קרוב מתכנסים לראשית.

למה המשפטים נכונים — סקיצת ההוכחה באינטואיציה

עבור משפט 1: נתון $\varepsilon$ . נבחר קו גובה $$V=c$$ שכלוא כולו בתוך כדור ברדיוס $\varepsilon$ (אפשר, כי $$V$$ רציפה וחיובית). נבחר $\delta$ קטן מספיק כך שכל מסלול המתחיל בכדור $\delta$ מתחיל עם $$V. מאחר ש־\(V$ אינה עולה לעולם, המסלול לא יוכל לחצות את קו הגובה $V=c$ החוצה — ולכן יישאר כלוא בתוך כדור $\varepsilon$. זו בדיוק הגדרת היציבות.$

עבור משפט 2: $$V(x(t))$$ היא פונקציה יורדת וחסומה מלמטה (באפס), ולכן מתכנסת לגבול כלשהו $c^\ast\geq 0$ . אילו $c^\ast>0$ , המסלול היה נשאר לנצח בטבעת שבה $\dot V$ שלילית וחסומה הרחק מאפס — אבל אז $$V$$ הייתה יורדת בקצב קבוע לפחות, וחוצה כל גבול כלפי מטה. סתירה. לכן $c^\ast=0$ , כלומר $V\to 0$ , ומכאן (מוגדרוּת חיובית) $x\to 0$ .

שימו לב לכיוון הלוגי של המשפטים

המשפטים הם תנאי מספיק, לא הכרחי: מציאת פונקציית ליאפונוב מוכיחה יציבות, אבל אי־מציאה אינה מוכיחה דבר! אם ניסיתם מועמדת $$V$$ וקיבלתם $\dot V$ שאינה בסימן הנכון — ייתכן שהמערכת יציבה ופשוט בחרתם מועמדת לא מוצלחת. נסו אחרת. (משפטי converse מבטיחים שלמערכת יציבה אסימפטוטית תמיד קיימת פונקציית ליאפונוב — אבל הם לא אומרים איך למצוא אותה. עוד על כך בפרק 14.)

פרק 9

יציבות גלובלית ואי־חסימות רדיאלית

המשפטים בפרק הקודם מקומיים — הם מבטיחים יציבות רק בסביבת הראשית. כדי לקבל יציבות גלובלית (מכל תנאי התחלה), נדרשת תוספת אחת שקל לפספס:

משפט 3 — יציבות אסימפטוטית גלובלית (ברבשין–קרסובסקי)

אם $$V$$ מוגדרת חיובית על כל $\mathbb{R}^n$ , הנגזרת $\dot V$ מוגדרת שלילית על כל $\mathbb{R}^n$ , ובנוסף $$V$$ בלתי־חסומה רדיאלית:

\|x\|\to\infty\;\Longrightarrow\;V(x)\to\infty

אז שיווי המשקל יציב אסימפטוטית גלובלית.

למה צריך את התנאי הנוסף?

התנאי (radially unbounded) מבטיח שכל קו גובה $$V=c$$ הוא קבוצה חסומה — "טבעת סגורה" שכולאת את המסלול. בלעדיו, קווי גובה עלולים "להיפתח לאינסוף": קיימות פונקציות מוגדרות חיוביות, כמו $V=\frac{x_1^2}{1+x_1^2}+x_2^2$ , שבהן קווי גובה גבוהים אינם סגורים, ומסלול יכול "לגלוש במורד $$V$$ " ובכל זאת לברוח לאינסוף בכיוון $$x_1$$ — כי בכיוון זה $$V$$ נשארת חסומה (קטנה מ־1) לעולם. אי־חסימות רדיאלית סוגרת את הפרצה: לרדת ב־ $$V$$ פירושו בהכרח להישאר באזור חסום.

בדיקה מהירה

כל פונקציה ריבועית $V=x^{\mathsf T}Px$ עם $P\succ0$ היא אוטומטית בלתי־חסומה רדיאלית (כי $V\geq\lambda_{\min}(P)\|x\|^2$ ). לכן במערכות ליניאריות התנאי הזה מתקיים "בחינם" — הוא הופך רלוונטי בעיקר במערכות לא־ליניאריות עם מועמדות לא־ריבועיות.

פרק 10

עקרון האינווריאנטיות של לה־סאל

בפועל, קורה שוב ושוב אותו תרחיש מתסכל: מוצאים $$V$$ טבעית ויפה (בדרך כלל האנרגיה), מחשבים $\dot V$ — ומקבלים ביטוי חצי־מוגדר שלילי בלבד: $\dot V\leq 0$ , עם התאפסות גם מחוץ לראשית. משפט 2 אינו חל, ומשפט 1 נותן רק "יציבות" — למרות שאנחנו בטוחים שהמערכת מתכנסת. הדוגמה הקלאסית: מטוטלת עם חיכוך, שבה $\dot V=-b\,\omega^2$ מתאפסת בכל פעם שהמהירות הזוויתית אפס — גם כשהזווית אינה אפס.

עקרון האינווריאנטיות של לה־סאל (LaSalle) מציל את המצב. האבחנה שלו: העובדה ש־ $\dot V=0$ בנקודה מסוימת אינה אומרת שהמסלול יכול להישאר שם. השאלה הנכונה היא — האם המסלול יכול לחיות לנצח בתוך הקבוצה שבה $\dot V=0$ ?

ניסוח העקרון

תהי $$V$$ מוגדרת חיובית עם $\dot V\leq 0$ בתחום חסום ואינווריאנטי $\Omega$ . נגדיר:

E=\{x\in\Omega:\;\dot V(x)=0\}

ותהי $$M$$ הקבוצה האינווריאנטית הגדולה ביותר המוכלת ב־ $$E$$ (קבוצה שמסלול המתחיל בה נשאר בה). אז כל מסלול מתכנס אל $$M$$ . בפרט, אם $M=\{0\}$ — הראשית יציבה אסימפטוטית.

איך משתמשים בזה בפועל — הבדיקה בשלושה צעדים

מצאו את $$E$$

פתרו $\dot V(x)=0$ . לדוגמה, במטוטלת: $\dot V=-b\omega^2=0\Rightarrow\omega=0$ — הקבוצה $$E$$ היא כל ציר הזוויות.

שאלו: מי יכול להישאר ב־ $$E$$ ?

הציבו את התנאי של $$E$$ בחזרה במשוואות המערכת ודרשו שהוא יישמר בזמן. במטוטלת: להישאר עם $\omega\equiv0$ דורש $\dot\omega=0$ , כלומר $\sin\theta=0$ — רק נקודות שיווי המשקל עצמן.

הסיקו

אם הקבוצה האינווריאנטית היחידה ב־ $$E$$ (בתחום הרלוונטי) היא הראשית — קיבלתם יציבות אסימפטוטית, למרות ש־ $\dot V$ רק חצי־מוגדרת.

האינטואיציה

המסלול יכול "לחלוף דרך" נקודות שבהן $\dot V=0$ — כמו מטוטלת שברגע השיא של התנופה מהירותה אפס והחיכוך רגעית אינו גורע אנרגיה. אבל היא לא יכולה להישאר שם: הכבידה מיד מאיצה אותה שוב, המהירות חוזרת להיות שונה מאפס, והאנרגיה ממשיכה לרדת. רק מי שיכול לגור לצמיתות בקבוצת ה"אפס־גריעה" ניצל מהשקיעה — ובמטוטלת, זה רק שיווי המשקל עצמו.

פרק 11

פונקציות ריבועיות ומערכות ליניאריות — הגשר למשוואת ליאפונוב

נפעיל כעת את השיטה על המקרה החשוב ביותר: מערכת ליניארית $\dot x=Ax$ עם מועמדת ריבועית $V=x^{\mathsf T}Px$ , $P=P^{\mathsf T}\succ0$ . נחשב את הנגזרת לאורך המסלולים:

\dot V=\dot x^{\mathsf T}Px+x^{\mathsf T}P\dot x =(Ax)^{\mathsf T}Px+x^{\mathsf T}P(Ax) =x^{\mathsf T}\bigl(A^{\mathsf T}P+PA\bigr)x

קיבלנו ש־ $\dot V$ היא בעצמה תבנית ריבועית, עם המטריצה $A^{\mathsf T}P+PA$ . כדי ש־ $\dot V$ תהיה מוגדרת שלילית, נדרוש שמטריצה זו תהיה שלילית מוגדרת — כלומר שתשווה למינוס מטריצה חיובית מוגדרת כלשהי $Q\succ0$ :

A^{\mathsf T}P+PA=-Q

זוהי משוואת ליאפונוב — התגלמותה האלגברית של השיטה הישירה בעולם הליניארי. במקום לחפש פונקציה, מחפשים מטריצה $$P$$ : בוחרים $Q\succ0$ (למשל $$Q=I$$ ), פותרים את המשוואה הליניארית עבור $$P$$ , ובודקים אם יצאה חיובית מוגדרת. המשפט המרכזי (שמוכח ונחקר לעומק בעמוד הנפרד על משוואת ליאפונוב) אומר:

המשפט — יציבות ליניארית בשפת ליאפונוב

המטריצה $$A$$ יציבה הורביץ (כל ערכיה העצמיים בחצי המישור השמאלי) אם ורק אם לכל $Q\succ0$ קיים פתרון יחיד $P\succ0$ למשוואת ליאפונוב. כלומר: במערכות ליניאריות, פונקציית ליאפונוב ריבועית תמיד קיימת כשהמערכת יציבה — והחיפוש אחריה הוא בעיה אלגברית פתירה לחלוטין.

סדר הפעולות הנכון — בוחרים את $$Q$$ , פותרים את $$P$$

טעות נפוצה היא לבחור $$P$$ (למשל $$P=I$$ ) ולבדוק מה יוצא $$Q$$ . הבעיה: אם יצא ש־ $$Q$$ אינה מוגדרת — לא הוכחתם כלום (אולי סתם בחרתם $$P$$ גרוע). הכיוון הנכון והמובטח הוא ההפוך: בחרו $Q\succ0$ , פתרו את $$P$$ . אם $$A$$ יציבה, מובטח ש־ $$P$$ תצא חיובית מוגדרת; אם $$P$$ לא יצאה חיובית מוגדרת — המערכת בהכרח אינה יציבה אסימפטוטית. המבחן חד־משמעי לשני הכיוונים.

מעבר לבדיקת יציבות, המבנה הריבועי־ליניארי הזה הוא אבן הפינה של הבקרה המודרנית כולה: משוואת ריקטי של LQR היא "משוואת ליאפונוב עם איבר ריבועי", הגרמיאנים של שליטוּת וצפיוּת הם פתרונות משוואת ליאפונוב, ותנאי הדטקטביליות שראינו בעמוד הקודם נוסח בדיוק כקיום פתרון לאי־שוויון ליאפונוב $(A-LC)^{\mathsf T}P+P(A-LC)\prec0$ .

פרק 12

איך מוצאים פונקציית ליאפונוב?

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

האנרגיה הפיזיקלית

אם למערכת יש מקור פיזיקלי — מכני, חשמלי, תרמי — נסו קודם את האנרגיה הכוללת. במערכות מכניות: קינטית + פוטנציאלית. במעגלים: $\tfrac12CV_C^2+\tfrac12LI_L^2$ . לרוב תקבלו $\dot V\leq0$ (חצי־מוגדרת) — ואז לה־סאל משלים את העבודה.

ריבועית פשוטה: $V=\|x\|^2$

המועמדת הזולה ביותר: $V=x_1^2+\dots+x_n^2$ . לוקח שניות לחשב את $\dot V$ , ולעיתים מפתיעות זה פשוט עובד (ראו דוגמה 1 בפרק הבא).

ריבועית כללית: $V=x^{\mathsf T}Px$

אם המערכת ליניארית (או שבודקים את הליניאריזציה) — פתרו את משוואת ליאפונוב עם $$Q=I$$ . זו דרך שיטתית לחלוטין, והפתרון $$P$$ משמש לעיתים כמועמדת טובה גם למערכת הלא־ליניארית המלאה בסביבת שיווי המשקל.

ריבועית + אינטגרל של האי־ליניאריות

כשיש איבר לא־ליניארי מהצורה $$g(x_1)$$ , נסו לכלול ב־ $$V$$ את $\int_0^{x_1}g(\sigma)\,d\sigma$ — בדיוק כמו שהאנרגיה הפוטנציאלית של המטוטלת היא אינטגרל של $\sin$ . תבניות אלה (Lur'e-type) נפוצות מאוד.

שיטות מתקדמות

שיטת הגרדיאנט המשתנה (variable gradient), פונקציות סכום־ריבועים (SOS) עם פותרים נומריים, ו־backstepping לבניית בקר ופונקציה יחד. אלה חומר לקורסי המשך — כאן חשוב רק לדעת שהן קיימות.

עצה מעשית לבחינה

בתרגילים ובבחינות, המועמדת כמעט תמיד אחת מהשלוש הראשונות: אנרגיה, $\|x\|^2$ , או ריבועית עם מקדמים שמותאמים כך שאיברים "לא נוחים" ב־ $\dot V$ יתבטלו. אם חישבתם $\dot V$ ונשאר איבר בעייתי עם סימן לא ידוע — נסו לשקלל: $$V=ax_1^2+bx_2^2$$ ובחרו $$a,b$$ כך שהאיבר המעורב ייעלם.

פרק 13

דוגמאות חישוב מלאות

שלוש דוגמאות בסדר עולה של תחכום: מערכת סקלרית שבה הליניאריזציה נכשלת וליאפונוב מצליח; מערכת מישורית עם לה־סאל; והמטוטלת עם חיכוך — הדוגמה הקנונית של השיטה כולה.

דוגמה 1 — $\dot x=-x^3$ : היכן שהליניאריזציה שותקת

מערכת סקלרית פשוטה: $\dot x=-x^3$ , שיווי משקל בראשית. נשים לב שהליניאריזציה סביב הראשית היא $\dot x\approx 0\cdot x$ — ערך עצמי אפס, על הציר המדומה, והשיטה העקיפה של ליאפונוב אינה מכריעה. ננסה את השיטה הישירה עם המועמדת הפשוטה ביותר:

V(x)=x^2\qquad\Rightarrow\qquad \dot V=\frac{\partial V}{\partial x}\,f(x)=2x\cdot(-x^3)=-2x^4

בדיקת התנאים: $$V=x^2$$ מוגדרת חיובית ובלתי־חסומה רדיאלית; $\dot V=-2x^4$ שלילית ממש לכל $x\neq0$ — מוגדרת שלילית. כל תנאי משפט 3 מתקיימים על כל הישר:

x=0\ \text{יציב אסימפטוטית גלובלית}

הלקח

זו הדגמה חיה לעוצמת השיטה הישירה: הליניאריזציה (ערך עצמי אפס) לא יכלה להכריע דבר, ואילו פונקציית ליאפונוב של שורה אחת נתנה את התוצאה החזקה ביותר האפשרית — יציבות גלובלית. אגב, ההתכנסות כאן אינה מעריכית (היא איטית, כמו $1/\sqrt{t}$ ) — ליאפונוב מוכיח התכנסות גם כשאין דעיכה אקספוננציאלית.

דוגמה 2 — מערכת מישורית: $\dot V$ חצי־מוגדרת ולה־סאל מציל

נתונה המערכת הליניארית:

\dot x_1=x_2,\qquad \dot x_2=-x_1-x_2

ננסה את המועמדת $$V=x_1^2+x_2^2$$ (מוגדרת חיובית, בלתי־חסומה רדיאלית). הנגזרת:

\dot V=2x_1\dot x_1+2x_2\dot x_2 =2x_1x_2+2x_2(-x_1-x_2) =-2x_2^2

האיברים המעורבים התבטלו ונשארנו עם $\dot V=-2x_2^2\leq0$ — אך זו רק חצי־מוגדרת שלילית: היא מתאפסת על כל הציר $$x_2=0$$ , לא רק בראשית. משפט 2 אינו חל ישירות. נפעיל את לה־סאל:

שלב 1: $E=\{\dot V=0\}=\{x_2=0\}$ — כל ציר $$x_1$$ . שלב 2: האם מסלול יכול להישאר ב־ $$E$$ ? להישאר עם $x_2\equiv0$ דורש גם $\dot x_2\equiv0$ . אבל מהמשוואה השנייה, עם $$x_2=0$$ :

\dot x_2=-x_1-x_2=-x_1=0\;\Longrightarrow\;x_1=0

שלב 3: הקבוצה האינווריאנטית הגדולה ביותר בתוך $$E$$ היא הראשית בלבד. לפי לה־סאל:

x=0\ \text{יציב אסימפטוטית גלובלית}

אימות צולב

המערכת ליניארית עם $A=\begin{bmatrix}0&1\\-1&-1\end{bmatrix}$ , שערכיה העצמיים הם $-\tfrac12\pm j\tfrac{\sqrt3}{2}$ — חלק ממשי שלילי, יציבה הורביץ. מסקנת לה־סאל מתאשרת. אגב, פתרון משוואת ליאפונוב $A^{\mathsf T}P+PA=-I$ נותן כאן $P=\begin{bmatrix}3/2&1/2\\1/2&1\end{bmatrix}\succ0$ — פונקציה ריבועית "משופרת" שעבורה $\dot V=-\|x\|^2$ מוגדרת שלילית ממש, בלי צורך בלה־סאל כלל. שתי הדרכים כשרות; לה־סאל חסך לנו את פתרון המשוואה.

דוגמה 3 — המטוטלת עם חיכוך: אנרגיה + לה־סאל

המטוטלת המרוסנת, עם $x_1=\theta$ (זווית) ו־ $x_2=\dot\theta$ (מהירות זוויתית), $a=g/\ell>0$ , $$b>0$$ מקדם החיכוך:

\dot x_1=x_2,\qquad \dot x_2=-a\sin x_1-b\,x_2

המועמדת — האנרגיה הפיזיקלית (פוטנציאלית + קינטית, בנרמול מתאים):

V(x)=a\,(1-\cos x_1)+\tfrac{1}{2}x_2^2

בתחום $|x_1|<\pi$ זו פונקציה מוגדרת חיובית (מתאפסת רק בראשית). נחשב את הנגזרת:

\dot V=a\sin x_1\cdot\dot x_1+x_2\,\dot x_2 =a x_2\sin x_1+x_2(-a\sin x_1-b x_2) =-b\,x_2^2\;\leq\;0

איברי הכבידה התבטלו במדויק (הם משמרים אנרגיה — רק ממירים פוטנציאלית לקינטית), ונשאר בדיוק ההספק שהחיכוך גורע. שוב חצי־מוגדרת — לה־סאל: $E=\{x_2=0\}$ ; להישאר שם דורש $\dot x_2=-a\sin x_1=0$ , כלומר $\sin x_1=0$ ; בתחום $|x_1|<\pi$ זו רק הראשית. מסקנה:

(\theta,\dot\theta)=(0,0)\ \text{יציב אסימפטוטית (מקומית)}

למה רק מקומית?

כאן היציבות אינה גלובלית — ובצדק פיזיקלי: אם מסובבים את המטוטלת חזק מספיק, היא עוברת מעל הראש ומקיפה. מתמטית, למערכת יש נקודות שיווי משקל נוספות ( $x_1=\pm\pi$ ), ומערכת עם יותר משיווי משקל אחד לעולם אינה גלובלית־אסימפטוטית לאף אחד מהם. הניתוח שלנו תקף בעמק סביב הראשית, כלומר לתנודות שאינן עוברות מעל הראש. גם המגבלה הזו נקראת בשפת ליאפונוב: "אומדן אגן המשיכה" (region of attraction) — קווי הגובה הסגורים של $$V$$ נותנים אומדן מובטח שלו.

פרק 14

טעויות נפוצות ונקודות עדינות

טעות 1 — "לא מצאתי V, לכן המערכת לא יציבה"

המשפטים חד־כיווניים: פונקציית ליאפונוב מוכיחה יציבות, אבל כישלון של מועמדת ספציפית לא מוכיח כלום. אם $\dot V$ יצאה עם סימן לא עקבי — המועמדת גרועה, לא בהכרח המערכת. להוכחת אי־יציבות יש משפטים ייעודיים (משפט צ'טאייב, ומשפט אי־היציבות של ליאפונוב) עם מבנה משלהם.

טעות 2 — לשכוח לבדוק את המוגדרוּת של V עצמה

סטודנטים ממהרים לחשב $\dot V$ ושוכחים שהמשפט דורש שני תנאים: $$V$$ מוגדרת חיובית וגם $\dot V$ בסימן הנכון. פונקציה כמו $$V=x_1^2$$ במישור אינה מוגדרת חיובית (רק חצי), ומסקנות ממנה אינן תקפות כלשונן.

טעות 3 — $\dot V\leq0$ ≠ יציבות אסימפטוטית

חצי־מוגדרוּת נותנת רק יציבות (משפט 1). כדי לשדרג ליציבות אסימפטוטית צריך או $\dot V<0$ ממש (משפט 2), או טיעון לה־סאל מלא ומנומק (מציאת $$E$$ , בדיקת האינווריאנטיות, והסקה). לכתוב " $\dot V\leq0$ ולכן מתכנס" בלי לה־סאל — זו טעות שמורידה נקודות, ובצדק: מטוטלת ללא חיכוך מקיימת $\dot V=0\leq 0$ ואינה מתכנסת לעולם.

טעות 4 — לשכוח את אי־החסימות הרדיאלית במסקנות גלובליות

בלי $V(x)\to\infty$ כש־ $\|x\|\to\infty$ , כל המסקנות מקומיות בלבד — גם אם $\dot V<0$ בכל המרחב. עם מועמדות ריבועיות זה חינם; עם מועמדות אחרות חובה לבדוק.

נקודה עדינה — משפטי ה־converse

האם למערכת יציבה תמיד קיימת פונקציית ליאפונוב? כן — זהו תוכן משפטי ההיפוך (converse Lyapunov theorems): אם שיווי המשקל יציב אסימפטוטית, מובטח שקיימת פונקציית ליאפונוב (ולמערכות ליניאריות אף ריבועית, דרך משוואת ליאפונוב). הקושי הוא שההוכחות אינן קונסטרוקטיביות — הן בונות את $$V$$ מתוך הפתרון עצמו, שאותו אין לנו. המסר: הכלי תמיד קיים עקרונית; האתגר הוא רק למצוא אותו.

נקודה עדינה — מערכות תלויות־זמן

כל מה שהצגנו תקף למערכות אוטונומיות ( $$f$$ אינה תלויה מפורשות בזמן). למערכות $\dot x=f(t,x)$ יש תורה מקבילה עם פונקציות $$V(t,x)$$ ודרישות נוספות (decrescent), ועקרון לה־סאל מוחלף בלֶמת ברבלט (Barbalat). זה חומר לקורס בקרה לא־ליניארית — כאן חשוב רק לדעת שהגבול קיים.

פרק 15

סיכום

פונקציית ליאפונוב היא הכללה מתמטית של האנרגיה: פונקציה סקלרית מוגדרת חיובית שאינה עולה לאורך מסלולי המערכת. קיומה מוכיח יציבות בלי לפתור את המערכת — ולכן זהו הכלי המרכזי, ולעיתים היחיד, לניתוח מערכות לא־ליניאריות. ההיררכיה: $\dot V\leq0$ נותן יציבות; $\dot V<0$ נותן יציבות אסימפטוטית; בתוספת אי־חסימות רדיאלית — גלובלית; וכשהנגזרת רק חצי־מוגדרת, עקרון לה־סאל משלים את ההתכנסות אם הקבוצה האינווריאנטית ב־ $\{\dot V=0\}$ היא הראשית בלבד.

בעולם הליניארי, הצבת מועמדת ריבועית $V=x^{\mathsf T}Px$ הופכת את כל התורה למשוואה אלגברית אחת — משוואת ליאפונוב $A^{\mathsf T}P+PA=-Q$ — שהיא נושא העמוד הבא, והבסיס לגרמיאנים, LQR ומבחני LMI בבקרה המודרנית.

שמונה כללי אצבע לפונקציות ליאפונוב

1. יציבות במובן ליאפונוב = "מתחיל קרוב, נשאר קרוב"; אסימפטוטית = גם מתכנס.

2. פונקציית ליאפונוב:

\(V>0\)

(מלבד בראשית) וגם

\dot V\leq0

\dot V=\nabla V\cdot f(x)

— מחושבת בלי לפתור את המערכת.

\dot V<0

ממש

\Rightarrow

יציבות אסימפטוטית.

5. מסקנה גלובלית דורשת גם

V\to\infty

כש־

\|x\|\to\infty

\dot V\leq0

בלבד? לה־סאל: אם רק הראשית שורדת ב־

\{\dot V=0\}

— יש התכנסות.

7. כישלון מועמדת אינו מוכיח אי־יציבות — המשפטים תנאי מספיק בלבד.

8. מערכת ליניארית +

V=x^{\mathsf T}Px

\Rightarrow

משוואת ליאפונוב

A^{\mathsf T}P+PA=-Q

השלב הבא

העמוד הבא מוקדש כולו למשוואת ליאפונוב — הגרסה האלגברית של השיטה למערכות ליניאריות: משפט הקיום והיחידות, הפתרון האינטגרלי, שיטות פתרון ידניות ונומריות, והקשר לגרמיאנים של שליטוּת וצפיוּת ולכל מבני ה־LMI של הבקרה המודרנית.

חזרה לראש העמוד