ראשי > כללי > העולם דרך עיניהם של מהנדסי חשמל (מטריצות, פרק הסיום) – על ייצוג במרחב המצב

העולם דרך עיניהם של מהנדסי חשמל (מטריצות, פרק הסיום) – על ייצוג במרחב המצב

זהו. הגיע רגע האמת.

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

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

זהירות מתמטיקה

***

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

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

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

עבור בעיות אלה (בעגה: אוסילטור מאולץ) נקבל משוואה אחרת שיש בה איבר מסוג חדש שנקרא לו איבר מקור.

הנה המשוואה המקורית:

Picture1

L הוא המרחק של הגוף מנקודת שיווי משקל, L עם שתי נקודות מעליו מסמל נגזרת שניה בזמן של המרחק מנקודת שיווי משקל.

והנה המשוואה המעודכנת:

Picture2

F הכוח החיצוני המופעל על הגוף.

הפתרון, אם כך, יהיה תלוי גם בתכונות הבעיה המקורית (בעגה: הבעיה ההומוגנית) וגם בתכונות הכוח החיצוני.

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

אבל,

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

ככה מהנדסי חשמל רואים את העולם.

***

בואו ונניח שניתן לתאר את המערכת על ידי שתי המשוואות הבאות:

Picture3

x הוא וקטור משתני המצב, u המקור, כלומר הכוח החיצוני, x עם נקודה למעלה מסמל נגזרת אחת בזמן של וקטור משתני המצב. y מסמל את אות היציאה של המערכת, למשל באוסילטור את המרחק מנקודת שיווי המשקל בכל רגע. A,B,C,D הם קבועים שאינם תלויים בזמן (בעגה מערכת כזאת נקראת LTI, כלומר linear-time-invariant).

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

המשוואה התחתונה מגדירה את אות היציאה שהחלטנו למדוד.

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

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

Picture4

סימנתי את איבר המקור F באות u מטעמי נוחות והרגל.

משתני העזר שלי הם:

Picture5

לכן שתי המשוואות שמייצגות אותן הן:

Picture6

נסגור את שתי המשוואות בכתב מטריצי:

Picture7

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

Picture8

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

Picture9

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

***

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

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

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

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

Picture10

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

קיבלנו שתי משוואות אלגבריות, כאשר אנחנו זוכרים שהמקדמים A,B,C,D הם מטריצות. נבודד את X מתוך המשוואה הראשונה באמצעות אלגברה של מטריצות ונציב אותו במשוואה השניה:

Picture11

I היא מטריצת היחידה.

קיבלנו ביטוי בתחום התדר עבור מוצא המערכת Y בהינתן המקור U. אם נחלק ביניהם נקבל ביטוי שנקרא פונקצית התמסורת (transfer function) של המערכת, כלומר מה יוצא ביחס למה שנכנס, הכל כתלות בתדר הנדנוד.

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

Picture12

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

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

***

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

זהו.

מטריצה אובר-אנד-אאוט.

מודעות פרסומת
  1. עדיין אין תגובות.
  1. No trackbacks yet.

כתיבת תגובה

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

הלוגו של WordPress.com

אתה מגיב באמצעות חשבון WordPress.com שלך. לצאת מהמערכת / לשנות )

תמונת Twitter

אתה מגיב באמצעות חשבון Twitter שלך. לצאת מהמערכת / לשנות )

תמונת Facebook

אתה מגיב באמצעות חשבון Facebook שלך. לצאת מהמערכת / לשנות )

תמונת גוגל פלוס

אתה מגיב באמצעות חשבון Google+ שלך. לצאת מהמערכת / לשנות )

מתחבר ל-%s

%d בלוגרים אהבו את זה: