דטקטביליות

פרק 1

מהי דטקטביליות ולמה היא נחוצה?

במערכת בקרה אמיתית כמעט לעולם איננו מודדים את כל המשתנים הפנימיים. במנוע אנו אולי מודדים מהירות אך לא זרם; בלוויין אנו מודדים זווית אך לא קצב זוויתי; בתהליך כימי אנו מודדים טמפרטורה אך לא ריכוז. המצב הפנימי $x\in\mathbb{R}^n$ חי "מאחורי הקלעים", ואנו רואים רק חלון צר אליו — את המוצא הנמדד $y\in\mathbb{R}^p$ .

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

ההגדרה האינטואיטיבית

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

למה לא פשוט לדרוש צפיוּת מלאה?

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

האנלוגיה — צוללת בערפל

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

מטרת העמוד

נצא מהמושג של צפיוּת, נראה כיצד הפירוק הקנוני מפצל כל מערכת לחלק צפוי וחלק סמוי, ונגדיר דטקטביליות בדיוק כתנאי על החלק הסמוי. נציג שני כלי בדיקה מעשיים — מבחן PBH ותנאי ליאפונוב — ונקשור את הכול לתצפתן ולדואליות עם ייצוּב (stabilizability). בסוף יהיו בידיכם הן ההבנה האינטואיטיבית והן היכולת המעשית לבדוק דטקטביליות בכל מערכת LTI.

פרק 2

רקע הכרחי: צפיוּת (Observability)

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

\dot x(t)=A\,x(t)+B\,u(t),\qquad y(t)=C\,x(t)+D\,u(t)

כאשר $A\in\mathbb{R}^{n\times n}$ , $C\in\mathbb{R}^{p\times n}$ , ושאר המטריצות בממדים תואמים. לעניין הצפיוּת והדטקטביליות מה שחשוב הוא אך ורק הזוג $$(A,C)$$ — הדינמיקה הפנימית ומה שהמוצא "רואה" ממנה. הקלט $$u$$ אינו משנה תכונות אלה, ולכן נוכל לבחון בלי הגבלת הכלליות את המערכת הלא־מאולצת $\dot x=Ax,\; y=Cx$ .

ההגדרה

הזוג $$(A,C)$$ נקרא צפוי אם ניתן לקבוע באופן יחיד את המצב ההתחלתי $$x(0)$$ מתוך ידיעת המוצא $$y(t)$$ על פני קטע זמן סופי כלשהו $$[0,T]$$ . באופן שקול: אם שני מצבים התחלתיים שונים מייצרים בדיוק את אותו מוצא — המערכת אינה צפויה.

מטריצת הצפיוּת

המבחן האלגברי הקלאסי משתמש במטריצת הצפיוּת, הבנויה מערימת $$C$$ כפול חזקות עולות של $$A$$ :

\mathcal{O}= \begin{bmatrix} C\\ CA\\ CA^2\\ \vdots\\ CA^{\,n-1} \end{bmatrix} \in\mathbb{R}^{\,np\times n}

הזוג $$(A,C)$$ צפוי אם ורק אם מטריצת הצפיוּת היא בעלת דרגה מלאה בעמודות:

\operatorname{rank}(\mathcal{O})=n

למה דווקא עד החזקה $$n-1$$ ?

לפי משפט קיילי־המילטון, כל חזקה $$A^k$$ עבור $k\geq n$ היא צירוף ליניארי של $I,A,\dots,A^{n-1}$ . לכן הוספת שורות $CA^n, CA^{n+1},\dots$ אינה יכולה להגדיל את הדרגה — הן כבר תלויות ליניארית בקיימות. די לעצור ב־ $CA^{n-1}$ .

תת־מרחב הלא־צפיוּת — לב העניין

כאשר המערכת אינה צפויה, ה"חוסר" מתרכז בתת־מרחב מוגדר היטב הנקרא תת־מרחב הלא־צפיוּת (unobservable subspace), המסומן $\mathcal{N}$ . הוא בדיוק גרעין מטריצת הצפיוּת:

\mathcal{N}=\ker(\mathcal{O})=\{x_0:\;Ce^{At}x_0=0\;\;\forall t\geq 0\}

במילים: $\mathcal{N}$ הוא אוסף כל המצבים ההתחלתיים שמייצרים מוצא אפס לחלוטין — מצבים ש"בלתי נראים" מהמוצא. אם $\mathcal{N}=\{0\}$ , המערכת צפויה (רק המצב האפסי בלתי נראה). זהו בדיוק תת־המרחב שדטקטביליות תעסוק בו: השאלה תהיה לא "האם הוא קיים?" אלא "האם מה שחי בתוכו דועך לאפס?"

תכונה מבנית חשובה — $\mathcal{N}$ שמור תחת $$A$$

תת־מרחב הלא־צפיוּת הוא A־אינווריאנטי: אם $x\in\mathcal{N}$ אז גם $Ax\in\mathcal{N}$ . המשמעות עמוקה — דינמיקת המערכת "שומרת" מצב סמוי בתוך התחום הסמוי; הוא לעולם לא "ידלוף" החוצה אל המוצא. לכן יש היגיון לדבר על הדינמיקה הפנימית של החלק הסמוי בנפרד, וזה בדיוק מה שמאפשר את הפירוק הקנוני בפרק 4.

פרק 3

הרעיון המרכזי: לראות מול לדעוך

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

שאלה על המצב	מה היא בודקת	הכלי
האם הוא נראֶה מהמוצא?	צפיוּת — האם המצב משפיע על $$y(t)$$	מטריצת הצפיוּת $\mathcal{O}$
האם הוא דועך מעצמו?	יציבות — האם המצב הולך לאפס כש־ $$u=0$$	הערכים העצמיים של $$A$$

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

הסיסמה שכדאי לזכור

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

ארבע משבצות — מפת המצב המלאה

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

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

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

פרק 4

הפירוק הקנוני לחלק צפוי וחלק סמוי

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

הבנייה

מאחר שתת־מרחב הלא־צפיוּת $\mathcal{N}$ הוא A־אינווריאנטי, נוכל לבחור בסיס חדש שבו ה־ $n_{\bar o}$ הקואורדינטות הראשונות פורשׂות את $\mathcal{N}$ (החלק הסמוי), והשאר משלימות אותו. בקואורדינטות אלה, באמצעות מטריצת טרנספורמציה הפיכה $$T$$ ו־ $\tilde x=T^{-1}x$ , המערכת מקבלת צורה משולשת בלוקים:

\tilde A=T^{-1}AT= \begin{bmatrix} A_{\bar o} & 0\\[2pt] A_{21} & A_{o} \end{bmatrix}, \qquad \tilde C=CT= \begin{bmatrix} 0 & C_{o} \end{bmatrix}

כאן $A_{\bar o}$ מתאר את הדינמיקה של החלק הסמוי (unobservable), ו־ $A_{o}$ את הדינמיקה של החלק הצפוי (observable). שימו לב לשני דברים מכריעים: הבלוק הימני־עליון ב־ $\tilde A$ הוא אפס (החלק הצפוי אינו "מזין" את הסמוי), והבלוק השמאלי ב־ $\tilde C$ הוא אפס — כלומר המוצא אינו תלוי כלל בחלק הסמוי. זה בדיוק מה ש"סמוי" אומר.

איור 1: המבנה המשולש של הפירוק הקנוני. החלק הסמוי (A_ō) עשוי להזין את החלק הצפוי (A_o) דרך A₂₁, אך לא להפך, ורק החלק הצפוי מגיע אל המוצא y. דטקטביליות פירושה שהבלוק הסמוי יציב מעצמו.

מדוע המבנה המשולש כל כך חשוב?

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

\operatorname{eig}(A)=\operatorname{eig}(A_{\bar o})\;\cup\;\operatorname{eig}(A_{o})

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

פרק 5

ההגדרה הפורמלית של דטקטביליות

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

ניסוח א — דרך החלק הסמוי

הזוג $$(A,C)$$ נקרא דטקטבילי אם הדינמיקה הסמויה $A_{\bar o}$ שבפירוק הקנוני היא יציבה הורביץ — כלומר כל ערכיה העצמיים נמצאים בחצי המישור השמאלי הפתוח:

(A,C)\ \text{דטקטבילי}\iff \operatorname{Re}\bigl(\lambda\bigr)<0\ \ \forall\,\lambda\in\operatorname{eig}(A_{\bar o})

ניסוח ב — דרך האופנים הלא־צפויים

באופן שקול ואלגנטי: ערך עצמי $\lambda$ של $$A$$ נקרא אופן לא־צפוי אם הוקטור העצמי המתאים לו שייך לתת־מרחב הלא־צפיוּת. אז:

(A,C)\ \text{דטקטבילי}\iff \text{כל אופן לא־צפוי הוא יציב}\;(\operatorname{Re}\lambda<0)

ניסוח ג — דרך התכנסות אסימפטוטית

ההגדרה ה"התנהגותית", הקרובה ביותר לאינטואיציה: המערכת דטקטבילית אם מוצא אפס מאלץ מצב מתכנס לאפס. כלומר, אם $y(t)=Ce^{At}x_0=0$ לכל $t\ge 0$ , אזי בהכרח $e^{At}x_0\to 0$ כש־ $t\to\infty$ :

Ce^{At}x_0\equiv 0 \;\Longrightarrow\; \lim_{t\to\infty}e^{At}x_0=0

זהו ההבדל המדויק מול צפיוּת: בצפיוּת, מוצא אפס מאלץ מצב אפס מיד ( $$x_0=0$$ ). בדטקטביליות, מוצא אפס מאלץ רק התכנסות לאפס בסופו של דבר. ההרפיה היא בדיוק המעבר מ"אפס מיידי" ל"אפס בגבול".

שרשרת ההיכלות — שננו אותה

הקשר בין שני המושגים הוא חד־כיווני וברור:

\text{צפיוּת}\;\Longrightarrow\;\text{דטקטביליות}

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

פרק 6

היחס לצפיוּת ולשליטוּת

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

תכונה	הזוג	שאלה	גרסה חזקה / חלשה
צפיוּת	$$(A,C)$$	האם כל מצב נראה מהמוצא?	הגרסה החזקה
דטקטביליות	$$(A,C)$$	האם כל אופן לא־צפוי יציב?	ההרפיה של צפיוּת
שליטוּת	$$(A,B)$$	האם ניתן להזיז כל מצב ע"י הקלט?	הגרסה החזקה
ייצוּב	$$(A,B)$$	האם כל אופן לא־שליט יציב?	ההרפיה של שליטוּת

הדואליות המושלמת

הזוגות מסודרים בדואליות מדויקת: צפיוּת היא הדואל של שליטוּת, ודטקטביליות היא הדואל של ייצוּב. הדואליות אינה רק אנלוגיה — היא זהות אלגברית מדויקת:

(A,C)\ \text{דטקטבילי} \;\iff\; (A^{\mathsf T},\,C^{\mathsf T})\ \text{ניתן לייצוב}

כלומר, כל שאלה על דטקטביליות של $$(A,C)$$ ניתנת לתרגום מיידי לשאלה על ייצוּב של הזוג המשוחלף $(A^{\mathsf T},C^{\mathsf T})$ . זה שימושי מאוד: כל משפט שנוכיח לצד אחד, מתקבל "בחינם" בצד השני.

למה זה כל כך שימושי בפועל

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

פרק 7

מבחן PBH לדטקטביליות

הפירוק הקנוני יפה תיאורטית, אך לחישוב מעשי הוא מסורבל — יש למצוא בסיס, לבנות $$T$$ , ולחשב בלוקים. למרבה המזל קיים מבחן ישיר ואלגנטי הרבה יותר: מבחן Popov–Belevitch–Hautus, או בקיצור PBH. הוא בודק דטקטביליות ערך עצמי אחר ערך עצמי, בלי שום פירוק.

הנוסח

הזוג $$(A,C)$$ דטקטבילי אם ורק אם עבור כל ערך עצמי $\lambda$ של $$A$$ שנמצא בחצי המישור הימני הסגור ( $\operatorname{Re}\lambda\ge 0$ ), מתקיים:

\operatorname{rank} \begin{bmatrix} A-\lambda I\\[2pt] C \end{bmatrix}=n

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

מדוע המבחן עובד — ההיגיון

נניח שהדרגה נופלת עבור $\lambda$ מסוים. אז קיים וקטור $v\neq 0$ בגרעין של שתי המטריצות יחד, כלומר:

(A-\lambda I)v=0\quad\text{וגם}\quad Cv=0

התנאי הראשון אומר ש־ $$v$$ הוא וקטור עצמי של $$A$$ עם ערך עצמי $\lambda$ . התנאי השני אומר ש־ $$v$$ סמוי מהמוצא ( $$Cv=0$$ , ולמעשה $Ce^{At}v=e^{\lambda t}Cv=0$ לכל זמן). אם בנוסף $\operatorname{Re}\lambda\ge 0$ — מצאנו אופן גם סמוי וגם לא־יציב, וזה בדיוק השילוב האסור. לכן המבחן בודק רק את הערכים העצמיים בחצי המישור הימני: רק שם טמונה הסכנה.

ההבדל מ־PBH לצפיוּת

מבחן PBH לצפיוּת דורש שדרגת המטריצה תהיה מלאה עבור כל הערכים העצמיים של $$A$$ , ללא יוצא מן הכלל. מבחן PBH לדטקטביליות מקל: הוא דורש זאת רק עבור הערכים העצמיים עם $\operatorname{Re}\lambda\ge 0$ . את הערכים העצמיים היציבים (בחצי השמאלי) מותר "לפספס" — אם הם סמויים, אין בכך כל בעיה, שהרי הם דועכים בלאו הכי. זו בדיוק ההרפיה, מתורגמת למבחן.

יתרון מעשי גדול

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

פרק 8

תנאי ליאפונוב לדטקטביליות

הכלי השלישי מגיע מתורת ליאפונוב, והוא חשוב במיוחד בהקשר של תכנון תצפתנים ושל בקרה אופטימלית (LQR/LQG). במקום לבדוק ערכים עצמיים, אנו פותרים אי־שוויון מטריצי ושואלים אם קיים פתרון חיובי מתאים.

הנוסח (משוואת ליאפונוב לתצפיתן)

הזוג $$(A,C)$$ דטקטבילי אם ורק אם קיימת מטריצה סימטרית חיובית מוגדרת $P=P^{\mathsf T}\succ 0$ ומטריצה $$L$$ כך ש:

(A-LC)^{\mathsf T}P+P\,(A-LC)\prec 0

הסימון $\prec 0$ פירושו "שלילית מוגדרת". אי־שוויון זה הוא בדיוק התנאי שמטריצת התצפתן $$A-LC$$ תהיה יציבה הורביץ עבור בחירה כלשהי של $$L$$ (מטריצת הגֵבר של התצפתן). וזה, כפי שנראה בפרק 10, הוא בדיוק מה שמאפשר לבנות תצפתן מתכנס.

שימו לב לכיוון הלוגי

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

הקשר ל־LMI ולכלים מספריים

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

פרק 9

שפת האופנים: מה באמת קורה במערכת

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

x(t)=e^{At}x_0=\sum_{i=1}^{n} c_i\,e^{\lambda_i t}\,v_i

כל אופן מתפתח בזמן כמו $e^{\lambda_i t}$ : אם $\operatorname{Re}\lambda_i<0$ הוא דועך, ואם $\operatorname{Re}\lambda_i\ge 0$ הוא גָדֵל או נשאר. המוצא "רואה" כל אופן דרך $$Cv_i$$ :

y(t)=Cx(t)=\sum_{i=1}^{n} c_i\,e^{\lambda_i t}\,(Cv_i)

אם $$Cv_i=0$$ , האופן ה־ $$i$$ פשוט נעלם מהמוצא — הוא סמוי. אם $Cv_i\neq 0$ , הוא תורם למוצא — הוא צפוי. וכעת נוכל לנסח את דטקטביליות בצורה הפשוטה ביותר:

דטקטביליות בשפת האופנים

המערכת דטקטבילית אם ורק אם אין אף אופן שהוא בו־זמנית סמוי ולא־יציב. כלומר, לכל אופן $$i$$ מתקיים: אם $$Cv_i=0$$ (סמוי), אז בהכרח $\operatorname{Re}\lambda_i<0$ (יציב). אופן יכול להיות סמוי, או לא־יציב, אך לעולם לא שניהם יחד.

שלוש דוגמאות־מחשבה זריזות

נניח מערכת אלכסונית פשוטה עם שני אופנים בעלי ערכים עצמיים $\lambda_1,\lambda_2$ ומטריצת מוצא $C=[c_1\;\;c_2]$ :

ערכים עצמיים	$$C$$	האופן הסמוי	דטקטבילי?
$-2,\,-5$	$[1\;\;0]$	$\lambda_2=-5$ (יציב)	כן
$-2,\,+3$	$[1\;\;0]$	$\lambda_2=+3$ (לא־יציב)	לא!
$-2,\,+3$	$[1\;\;1]$	אין אופן סמוי	כן (צפוי לגמרי)

השורה האמצעית היא הלקח: ערך עצמי לא־יציב ( $$+3$$ ) ביחד עם $$C$$ שלא "רואה" אותו ( $$c_2=0$$ ) — זהו בדיוק הצירוף ההרסני. ברגע ש־ $$C$$ כן רואה אותו (שורה אחרונה), הכול בסדר.

פרק 10

הקשר לתצפתן (Observer) — הסיבה האמיתית שדטקטביליות חשובה

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

מהו תצפתן?

תצפתן (או "מעריך מצב") הוא מערכת דינמית מלאכותית שאנו מריצים במקביל למערכת האמיתית. הוא מקבל את אותו קלט $$u$$ ואת המוצא הנמדד $$y$$ , ומייצר אומדן $\hat x$ של המצב הפנימי. תצפתן לואנברגר נראה כך:

\dot{\hat x}=A\hat x+Bu+L\,(y-C\hat x)

האיבר $L(y-C\hat x)$ הוא "תיקון": הוא משווה את המוצא הנמדד $$y$$ למוצא החזוי $C\hat x$ , ומזריק את ההפרש בחזרה כדי לתקן את האומדן. $$L$$ היא מטריצת הגֵבר שאנו בוחרים.

דינמיקת השגיאה — כאן הכול מתחבר

נגדיר את שגיאת האומדן $e=x-\hat x$ . בחיסור משוואת התצפתן ממשוואת המערכת מקבלים דינמיקה נקייה ויפהפייה לשגיאה:

\dot e=(A-LC)\,e

זוהי משוואה הומוגנית. השגיאה תתכנס לאפס (האומדן יתכנס למצב האמיתי) אם ורק אם המטריצה $$A-LC$$ היא יציבה הורביץ. השאלה היחידה היא: האם קיימת בחירה של $$L$$ שמייצבת אותה?

המשפט המרכזי

קיים גֵבר תצפתן $$L$$ ההופך את $$A-LC$$ ליציבה אם ורק אם הזוג $$(A,C)$$ דטקטבילי. יתרה מזו, אם המערכת צפויה במלואה, ניתן למקם את הערכים העצמיים של $$A-LC$$ בכל מקום שנרצה (כמה שרוצים מהיר). אם המערכת רק דטקטבילית, נוכל למקם את הערכים הצפויים כרצוננו, אך הסמויים יישארו במקומם — וזה בסדר גמור, כי הם כבר יציבים.

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

מה קורה אם המערכת לא דטקטבילית?

אם יש אופן סמוי ולא־יציב, אז $$A-LC$$ יכילה אותו ערך עצמי לא־יציב לכל בחירת $$L$$ — כי $$LC$$ לא נוגע באופן הסמוי. השגיאה תתבדר, והאומדן יהיה חסר תועלת. אין תצפתן שיכול לעקוב אחר מצב שגם לא־יציב וגם בלתי נראה. זה לא עניין של תכנון חכם יותר — זה חסם עקרוני.

פרק 11

דואליות: דטקטביליות וייצוּב (Stabilizability)

הזכרנו את הדואליות בפרק 6; כעת נעמיק בה, כי היא מאחדת את כל התמונה. ייצוּב הוא ההרפיה של שליטוּת בדיוק כפי שדטקטביליות היא ההרפיה של צפיוּת. נעמיד את ההגדרות זו מול זו:

הזוג $$(A,C)$$ — כל אופן לא־צפוי יציב. תנאי לתצפתן: קיים $$L$$ כך ש־ $$A-LC$$ יציבה.

ייצוּב

הזוג $$(A,B)$$ — כל אופן לא־שליט יציב. תנאי לבקר: קיים $$K$$ כך ש־ $$A-BK$$ יציבה.

מבחן PBH (דטקט')

לכל $\operatorname{Re}\lambda\ge0$ : $\operatorname{rank}\!\begin{bmatrix}A-\lambda I\\C\end{bmatrix}=n$ . ערימה אנכית.

מבחן PBH (ייצוב)

לכל $\operatorname{Re}\lambda\ge0$ : $\operatorname{rank}\bigl[\,A-\lambda I\ \ B\,\bigr]=n$ . ערימה אופקית.

הזרקת המוצא

$$A-LC$$ — מזריקים את המוצא דרך $$L$$ כדי לתקן את האומדן ולייצב את השגיאה.

משוב המצב

$$A-BK$$ — מזינים את הקלט דרך $$K$$ כדי לייצב את המערכת עצמה.

עקרון ההפרדה — הפרס הגדול

הדואליות הזו מובילה לאחת התוצאות היפות בבקרה: עקרון ההפרדה. אם המערכת ניתנת לייצוב וגם דטקטבילית, ניתן לתכנן את הבקר ואת התצפתן בנפרד לחלוטין, ולחבר אותם יחד — והמערכת המשולבת תהיה יציבה. הערכים העצמיים של המערכת בלולאה סגורה הם פשוט איחוד הערכים העצמיים של הבקר ( $$A-BK$$ ) ושל התצפתן ( $$A-LC$$ ):

\operatorname{eig}(\text{closed loop})=\operatorname{eig}(A-BK)\;\cup\;\operatorname{eig}(A-LC)

מה זה אומר בפועל

ייצוּב מבטיח שנוכל לבנות בקר משוב־מצב שמייצב; דטקטביליות מבטיחה שנוכל לבנות תצפתן שמשחזר את המצב הדרוש לבקר זה. יחד, השניים מספקים פתרון בקרה מלא אפילו כשאיננו מודדים את כל המצבים ואיננו שולטים בכל הכיוונים. זוהי בדיוק הסיבה ששתי ההרפיות האלה — ייצוּב ודטקטביליות — הן הדרישות "הנכונות" בתורת הבקרה המודרנית, ולא שליטוּת וצפיוּת המלאות.

פרק 12

מתכון בדיקה שלב־אחר־שלב

נסכם את כל הכלים לכדי תהליך עבודה אחד שמכריע דטקטביליות עבור כל זוג $$(A,C)$$ :

מצאו את הערכים העצמיים של $$A$$

חשבו את $\det(A-\lambda I)=0$ וקבלו את $\lambda_1,\dots,\lambda_n$ . סמנו את אלה עם $\operatorname{Re}\lambda\ge0$ — רק הם "מסוכנים".

בדיקה מהירה — האם בכלל יש סכנה?

אם כל הערכים העצמיים יציבים ( $\operatorname{Re}\lambda<0$ ), המערכת דטקטבילית אוטומטית — אין מה לבדוק. סיימתם. אחרת, המשיכו לערכים הלא־יציבים.

הפעילו PBH על כל ערך עצמי לא־יציב

לכל $\lambda$ עם $\operatorname{Re}\lambda\ge0$ , בנו את $\begin{bmatrix}A-\lambda I\\ C\end{bmatrix}$ ובדקו אם דרגתה $$n$$ .

הכריעו

אם הדרגה מלאה ( $$=n$$ ) עבור כל הערכים הלא־יציבים — המערכת דטקטבילית. אם לפחות אחד נכשל (דרגה $\() — המערכת אינה דטקטבילית.$

(לפי הצורך) פירוק לאימות והבנה

אם תרצו לראות את החלק הסמוי במפורש, חשבו את $\mathcal{N}=\ker(\mathcal{O})$ , בצעו את הפירוק הקנוני, ואמתו ש־ $A_{\bar o}$ יציבה.

(אם דטקטבילית) תכננו תצפתן

בחרו $$L$$ כך ש־ $$A-LC$$ יציבה (מיקום קטבים על החלק הצפוי), והריצו $\dot{\hat x}=A\hat x+Bu+L(y-C\hat x)$ .

פרק 13

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

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

דוגמה 1 — דטקטבילית אך לא צפויה

נתונה המערכת:

A=\begin{bmatrix}-2 & 0\\ 0 & -7\end{bmatrix}, \qquad C=\begin{bmatrix}1 & 0\end{bmatrix}

שלב 1 — ערכים עצמיים. המטריצה אלכסונית, ולכן הערכים העצמיים גלויים: $\lambda_1=-2,\ \lambda_2=-7$ . שני הערכים בחצי המישור השמאלי — שניהם יציבים.

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

בדיקת צפיוּת. נבנה את מטריצת הצפיוּת ( $$n=2$$ ):

\mathcal{O}=\begin{bmatrix}C\\ CA\end{bmatrix} =\begin{bmatrix}1 & 0\\ -2 & 0\end{bmatrix}

העמודה השנייה כולה אפס, ולכן $\operatorname{rank}(\mathcal{O})=1<2$ — המערכת אינה צפויה. תת־מרחב הלא־צפיוּת הוא $\mathcal{N}=\operatorname{span}\{[0\;\;1]^{\mathsf T}\}$ , כלומר האופן השני ( $\lambda_2=-7$ ) סמוי.

שלב 3 — אימות PBH. אין ערכים עצמיים לא־יציבים, ולכן אין מה לבדוק — התנאי מתקיים ריקנית. אך לשם ההמחשה, נבדוק את הערך הסמוי $\lambda_2=-7$ :

\begin{bmatrix}A-(-7)I\\ C\end{bmatrix} =\begin{bmatrix}5 & 0\\ 0 & 0\\ 1 & 0\end{bmatrix} \;\Rightarrow\;\operatorname{rank}=1<2

הדרגה אכן נופלת — מה שמאשר ש־ $\lambda_2$ לא־צפוי. אבל מאחר ש־ $$-7<0$$ (יציב), זה לא פוגע בדטקטביליות. PBH לא דורש לבדוק ערך זה כלל.

המסקנה

המערכת דטקטבילית אך לא צפויה — בדיוק התרחיש המעניין. האופן הסמוי ( $e^{-7t}$ ) דועך מהר מעצמו, ולכן איננו צריכים לראותו. תצפתן יוכל לעקוב אחר המצב המלא בלי כל בעיה, אף שאינו "רואה" את הקואורדינטה השנייה ישירות.

דוגמה 2 — לא דטקטבילית

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

A=\begin{bmatrix}-2 & 0\\ 0 & +4\end{bmatrix}, \qquad C=\begin{bmatrix}1 & 0\end{bmatrix}

שלב 1 — ערכים עצמיים. $\lambda_1=-2$ (יציב) ו־ $\lambda_2=+4$ (לא־יציב — בחצי הימני). יש סכנה; חובה לבדוק.

שלב 3 — PBH על הערך הלא־יציב $\lambda_2=+4$ :

\begin{bmatrix}A-4I\\ C\end{bmatrix} =\begin{bmatrix}-6 & 0\\ 0 & 0\\ 1 & 0\end{bmatrix}

העמודה השנייה כולה אפס, ולכן הדרגה היא $$1<2=n$$ . המבחן נכשל.

שלב 4 — הכרעה. קיים אופן (הוקטור העצמי $v_2=[0\;\;1]^{\mathsf T}$ ) שהוא גם סמוי ( $$Cv_2=0$$ ) וגם לא־יציב ( $\lambda_2=+4>0$ ). זהו בדיוק הצירוף האסור. המערכת אינה דטקטבילית.

מה המשמעות הפיזיקלית

הקואורדינטה השנייה גָדֵלה כמו $e^{4t}$ — מתבדרת — אך המוצא $$y=x_1$$ אינו רואה אותה כלל. שום מדידה לא תתריע על הבדידה הזו, ושום תצפתן לא יוכל לעקוב אחריה. אם ננסה לבנות תצפתן, ערך עצמי $$+4$$ יישאר ב־ $$A-LC$$ לכל בחירת $$L$$ , ושגיאת האומדן תתפוצץ.

דוגמה 3 — PBH על מערכת בצורת ז'ורדן (ערך עצמי כפול)

המקרים העדינים ביותר הם ערכים עצמיים חוזרים. נבחן מערכת $3\times3$ עם בלוק ז'ורדן:

A=\begin{bmatrix} 3 & 1 & 0\\ 0 & 3 & 0\\ 0 & 0 & -1 \end{bmatrix}, \qquad C=\begin{bmatrix}0 & 1 & 0\end{bmatrix}

שלב 1 — ערכים עצמיים. המטריצה משולשת בלוקים: $\lambda=3$ (ריבוי אלגברי 2, בלוק ז'ורדן) ו־ $\lambda=-1$ (פשוט, יציב). הערך $\lambda=3>0$ לא־יציב — חובה לבדוק. את $$-1$$ (יציב) נתעלם.

שלב 3 — PBH על $\lambda=3$ :

\begin{bmatrix}A-3I\\ C\end{bmatrix} =\begin{bmatrix} 0 & 1 & 0\\ 0 & 0 & 0\\ 0 & 0 & -4\\ 0 & 1 & 0 \end{bmatrix}

נבדוק את דרגת העמודות. העמודה הראשונה כולה אפס, ולכן הדרגה לכל היותר $$2$$ . אכן: העמודה השנייה ( $[1,0,0,1]^{\mathsf T}$ ) והשלישית ( $[0,0,-4,0]^{\mathsf T}$ ) בלתי־תלויות, אך העמודה הראשונה אפסית. לכן:

\operatorname{rank}\begin{bmatrix}A-3I\\ C\end{bmatrix}=2<3=n

הכרעה. המבחן נכשל עבור הערך הלא־יציב $\lambda=3$ . המערכת אינה דטקטבילית. הסיבה: הוקטור העצמי של בלוק ז'ורדן הוא $v=[1,0,0]^{\mathsf T}$ , ומתקיים $$Cv=0$$ — הכיוון הזה לא־יציב וסמוי בו־זמנית.

הלקח מצורת ז'ורדן

ערך עצמי כפול עם בלוק ז'ורדן מסוכן במיוחד: גם אם $$C$$ "רואה" חלק מהבלוק (כאן את הקואורדינטה השנייה), הוקטור העצמי האמיתי עשוי להישאר סמוי. תמיד בדקו את PBH על הוקטור העצמי, לא על קואורדינטות שרירותיות. כאן, אילו $$C$$ היה $[1\;\;0\;\;0]$ (רואה את ראש הבלוק), המערכת הייתה דטקטבילית.

פרק 14

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

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

טעות 1 — בלבול בין דטקטביליות ליציבות

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

טעות 2 — לבדוק PBH על כל הערכים העצמיים

לדטקטביליות בודקים PBH רק על ערכים עצמיים עם $\operatorname{Re}\lambda\ge0$ . בדיקה על ערך עצמי יציב היא בזבוז — ואם תכשילו את המערכת בגללה, תטעו: כשלון PBH על ערך עצמי יציב פירושו רק שהוא לא־צפוי, וזה לגיטימי לחלוטין בדטקטביליות.

טעות 3 — שכחת הציר המדומה

חצי המישור הימני במבחן הוא סגור — כולל את הציר המדומה ( $\operatorname{Re}\lambda=0$ ). ערכים עצמיים על הציר המדומה (כמו $\pm j\omega$ טהורים, או $\lambda=0$ ) אינם יציבים אסימפטוטית — הם אינם דועכים. לכן הם חייבים להיות צפויים, ויש לכלול אותם בבדיקת PBH. השמטתם היא טעות שכיחה.

נקודה עדינה — דטקטביליות וזמן בדיד

כל מה שאמרנו תקף למערכות רציפות, שבהן "יציב" פירושו $\operatorname{Re}\lambda<0$ . במערכות בזמן בדיד ( $$x[k+1]=Ax[k]$$ ), הגבול משתנה: "יציב" פירושו $|\lambda|<1$ (בתוך מעגל היחידה). בהתאם, מבחן PBH לדטקטביליות בזמן בדיד בודק את הערכים העצמיים עם $|\lambda|\ge1$ . הרעיון זהה לחלוטין — רק ה"אזור המסוכן" משתנה ממחצית מישור למעגל.

נקודה עדינה — דטקטביליות חזקה (detectability vs. strong detectability)

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

פרק 15

סיכום

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

הפירוק הקנוני מפצל את המערכת לחלק צפוי וחלק סמוי, והערכים העצמיים מתפצלים בהתאם. דטקטביליות פירושה שהבלוק הסמוי $A_{\bar o}$ יציב הורביץ. שלושה כלים מעשיים בודקים זאת: הפירוק (תיאורטי), מבחן PBH (ישיר, ערך עצמי אחר ערך עצמי), ותנאי ליאפונוב/LMI (נומרי). וכל זה משרת מטרה אחת: דטקטביליות היא התנאי המדויק לקיומו של תצפתן מתכנס, הדואל של ייצוּב בעולם הבקרה.

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

1. דטקטבילי

\iff

כל אופן לא־צפוי יציב.

2. צפיוּת

\Rightarrow

דטקטביליות (לא להפך).

3. כל הערכים יציבים

\Rightarrow

דטקטבילי אוטומטית.

4. PBH:

\operatorname{rank}\!\begin{bmatrix}A-\lambda I\\C\end{bmatrix}=n

לכל

\operatorname{Re}\lambda\ge0

5. בודקים PBH רק על ערכים עצמיים לא־יציבים.

6. דטקטביליות

\ne

יציבות. מערכת לא־יציבה יכולה להיות דטקטבילית.

7. הדואל של ייצוּב:

\((A,C)\)

דטקט'

\iff (A^{\mathsf T},C^{\mathsf T})

ניתן לייצוב.

8. תנאי לתצפתן: קיים

\(L\)

כך ש־

\(A-LC\)

יציבה.

השלב הבא

לאחר שליטה בדטקטביליות ובדואל שלה — ייצוּב — ניתן להתקדם אל עקרון ההפרדה ובניית בקר־תצפתן מלא, ומשם אל בקרה אופטימלית (LQR/LQG ומסנן קלמן), שם דטקטביליות תהיה תנאי הכרחי להתכנסות המסנן ולקיום פתרון יציב.

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