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

סיכום

פרוטוקול זה יוצר פרוצדורת אבטיפוס מהירה סטנדרטית להערכת בקרת מעקב לאחור מותאמת על ידי Particle Swarm Optimization בסביבת סימולציה בזמן אמת בשלב קבוע.

תקציר

המטרה העיקרית של פרוטוקול זה היא לספק מסגרת סימולציה קבועה שניתן לשחזר להערכת כוונון רווח מבוסס Particle Swarm (PSO) במערכות בקרה לא ליניאריות. המימוש מתחיל בניסוח מודל מטוטלת מסוג פורוטה, ואחריו אינטגרציה של בקר צעד אחורה בתוך סביבת ביצוע צעד קבוע של 2 מילישניות. המתודולוגיה כוללת תהליך שיטתי בן ארבעה שלבים: אפיון מגבלות יישום לא אידיאליות, הגדרת מרחב חיפוש אופטימיזציה רב-מטרותית ל-Particle Swarm, ביצוע כוונון אוטומטי לא מקוון, והערכת הפרמטרים המתקבלים באמצעות מערך סטנדרטי של תרחישי מעקב מסלול ודחיית הפרעות. הגדרה זו, המשתמשת בתחנות עבודה תעשייתיות בעלות ביצועים גבוהים ובממשקי אות סטנדרטיים, תומכת בהשוואות עקביות בניסוי חוזר בתוך אותה ארכיטקטורת בקרה. העיצוב משווה בין בקר אחורה, מכוון ידנית, לבין גרסה מותאמת ל-PSO, שמשתפת את אותו מבנה בקרה בדיוק, ובכך מבודדת את השפעת בחירת הגיין. ביצועי הבקרה מוערכים בשלושה תרחישים תפעוליים מובחנים: מעקב מסלול צעד, מעקב סינוסואידלי בתדרים מעורבים, ודחיית הפרעות. ניתוח סטטיסטי של 10 ניסויים חוזרים הראה כי אופטימיזציה מבוססת PSO הפחיתה את קצב RMSE במעקב צעדים מ-0.065 ל-0.050 רדיאציה והפחית את שיא המטוטלת ב-33.1%. שיפורים אלו הושגו לצד הפחתה של 22.1% במאמצי בקרת ה-RMS, מה שמעיד שהפרמטרים האופטימליים אפשרו הפצת אנרגיה יעילה יותר במסגרת המבוססת על ליאפונוב. בסופו של דבר, מתודולוגיה זו מספקת מסגרת סימולציה מובנית להערכת אסטרטגיות בקרה לא ליניאריות לפני כל יישום חומרה פיזית עתידי.

מבוא

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

בספרות הרחבה יותר, מחקרים רבים חקרו שינויים מבניים וכיוונון פרמטרים לשיפור השליטה במטוטלות הפוכות סיבובית. עיצובים מבוססי אופטימיזציה נפוצים; לדוגמה, PSO שימש לבחירת פרמטרי בקרה ושולב בארכיטקטורות בקרה מטושטשת ^4,5. בהשוואה למטה-היוריסטיקות אחרות בהשראת ביולוגיה, PSO נבחר במיוחד למסגרת זו בשל ההתכנסות המהירה שלו במרחבי חיפוש רציפים בממדים נמוכים ודרישות הכוונון המינימלית של היפרפרמטרים. ספרות עדכנית מדגישה יותר ויותר את הצורך באלגוריתמים לאופטימיזציה חכמה בתרחישי בקרה מגוונים ומורכבים. לדוגמה, טכניקות אופטימיזציה מתקדמות שולבו ביעילות עם בקרה אדפטיבית (MRAC) ומסגרות סדר חלקי כדי לשפר את דיוק המעקב של תחנות סרבו לא-ליניאריות ^6,7. יתרה מזאת, כיוון מבוסס אופטימיזציה הוכיח את עצמו כיעיל מאוד בניהול הדינמיקה המקושרת והמגבלות הטבועות במערכות אלקטרומכניות מורכבות ^8,9.

מחקרים עדכניים מאשרים כי PSO והווריאציות ההיברידיות שלו משפרים משמעותית את יעילות מעקב נקודות ההספק המקסימלית במערכי פוטו-וולטאיים, ומדגימים זיהוי פרמטרים עמיד בתנאי הצללה חלקית¹⁰. ברובוטיקה, PSO שימש בהצלחה לאופטימיזציה של עיצובים ליניאריים מוגברים ולא ליניאריים של בקרה פרופורציונלית-נגזרת עבור מניפולטורים מקביליים, תוך מזעור שגיאות מעקב מסלול¹¹. בנוסף, שילוב PSO עם בקרת מצב החלקה אדפטיבית לאחור הוכיח את עצמו כקריטי לדיכוי רטט במסת שריר מלאכותית^{תלויה 12}. מעבר לבחירת פרמטרים בסיסיים, שילוב עיבוד אותות מודרני ואסטרטגיות אופטימיזציה עמידות הוא קריטי לשמירה על יציבות בלולאה סגורה בתנאים פיזיקליים ריאליסטיים ורועשים^13,14.

מאמצים אחרונים התמקדו גם במעקב סימולטני אחר זווית מפרקים וייצוב מטוטלת בתנאים לא ודאיים, תוך שימוש בהיפוך דינמי כללי חזק¹⁵, וכן בהערכה עצבית אדפטיבית¹⁶. יתרה מזאת, ארכיטקטורות בקרת צעד אחורה פותחו רבות כדי להתמודד עם הפרעות מורכבות במערכות מכניות שונות, כגון שילוב עיצובים של מצבי החלקה לדיכוי רעידות בניין, שימוש במשקיפים כמעט-מחליקים לשסתומי מצערת אלקטרוניים, ואינטגרציה של צופים לא ליניאריים לוויסות מהירות מנוע DC מדויק 17,18,19. יישומים מגוונים אלו מדגישים את הרב-שימושיות של עיצובים של צעד אחורה בשילוב עם אסטרטגיות הערכה או אופטימיזציה חזקות.

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

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

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

פרוטוקול

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

1. בניית מפעל והקמת כנסי אותות

1. לבנות מודל סימולציה באיכות גבוהה של מטוטלת סיבובית הפוכה מסוג פורוטה. הגדר שתי קואורדינטות כלליות למערכת: זווית הזרוע הסיבובית, θ, וזווית סטיית המטוטלת מהאנכי הישרי, α.
2. מגדירים את וקטור מצב המערכת כ-x = [θ, α, θ̇, α̇]T, ומגדירים את α = 0 כשיווי משקל זקוף לא יציב. ניסח את מודל הצמח המתמטי באמצעות משוואות אוילר-לגראנז' כדי להגדיר במפורש את מטריצת האינרציה הלא-ליניארית, כוחות קוריוליס וקטורים כבידתיים המחברים בין הזרוע הסיבובית הפעילה למטוטלת הפסיבית²².
3. הגדר θ חיובי וחיובי α בעקביות תוך שימוש באותה פולריות מקודד בין סקריפטים של המפעל, הבקר ועיבוד הפוסט-עיבוד.
4. הקם את פותר הסימולציה בזמן אמת בשלב קבוע, תוך שימוש מדויק בפרמטרים הנומינליים של מכאניקה ואלקטרומכנית המפורטים בטבלה 1. הגדר את מערכת הלולאה הסגורה בתוך סביבת הביצוע בזמן אמת המיועדת כפי שמצוין בטבלת החומרים. הגדר את תקופת העדכון של הבקר ל-2 מילישניות וודא שתקופת רישום הנתונים קבועה ב-20 מילישניות.
5. שלבו שלושה אפקטים לא אידיאליים ספציפיים במודל המפעל כדי לדמות מגבלות יישום ריאליסטיות. הזריק רעש מדידה לבן בממוצע אפס עם סטיות תקן של 0.003 רדיאט ו-0.004 רדיאיום לתעלות המשוב של הזרוע והמטוטלת, בהתאמה.
6. הכנס עיכוב העברה של 4 מילישניות לאות המטוטלת הנמדד. במקביל, הגדל את מקדמי הדיכוי של המפעל המדומים ב-8% ביחס לערכים הנומינליים בעיצוב הבקרה כדי לדמות אי-התאמה במודל.
7. הגבל את פלט פקודת המנוע ל-±10 וולט בלבד.
   אזהרה: הגדר עצירות בטיחות ברמת תוכנה כדי לסיים מיד את ביצוע הניסוי אם |θ| עובר 0.70 רדיאט או אם |α| עולה על 0.35 רדיאטור. בנוסף, יש ליישם טיימר שמירה כדי להפעיל איפוס מערכת אם פקודת הבקרה נשארת במגבלת רוויה של ±10 וולט ליותר מ-100 מילישניות.

2. הגדרה ראשונית למדינה ומיקום מראש

1. מאתחל את המערכת במצב מנוחה כלפי מטה לפני כל בדיקה.
2. בצע משימת מיקום מראש כדי להניף את המטוטלת למצב ההפוך ולהפעיל את בקר המייצב. החריגו את שלב ה-prepositioning הזה מהערכת הביצועים הפורמלית.
3. התחל באיסוף נתונים רשמי רק לאחר אימות שהסטייה של המטוטלת נשארה ברציפות בטווח של ±0.05 רדיאדה מהשיווי משקל הישרי למשך לפחות 1.0 שנייה, ובכך לשמור על תנאי ההתחלה הדינמיים ניתנים להשוואה מעשית לכל הניסויים שנערכו.
4. איפס את זמן הניסיון ל-0.0 שניות לאחר עמידה בתנאי זה והתחל מיד לבצע את פרופיל נקודת ההגדרה המיועד. ודא שכל הבקרים משתמשים באותו תהליך פריפוזיציה ועומדים בדיוק באותם סטנדרטים של קבלה.

3. מימוש בקר בסיסי לאחור

1. הגדר את שגיאת המעקב eθ(t) לפי הפניה הזווית-הזרועית θr(t) באמצעות משוואה
   eθ(t) = θr(t) − θ(t) (1)
2. בנה את משטח ההחלקה המורכב, s(t), לפי משוואה 2. בתוך ניסוח זה, הגדר λ כמקדם השיפוע של משטח השגיאה החיובי לחלוטין.
   s(t) = ėθ(t) + λeθ(t) (2)
3. ניסחו את בקר הצעד האחורי הייחוס על ידי אינטגרציה של מונח פיצוי מבוסס מודל נומינלי עם מונח פרופורציונלי מייצב שגיאות. הגדר מתמטית את משתני שגיאת המעקב כ-z1 = eθ(t) ו-z2 = s(t).
4. בצע ניתוח יציבות ליאפונוב להערכת התכנסות שגיאת מעקב. גזיר את הבקרה הנבחרת פונקציית ליאפונוב V = 1/2z1² + 1/2z2 ביחס לזמן כדי לקבל V̇ = z1ż1 + z2ż2. החליפו את חוק הבקרה המנוסח בנגזרת זו כדי לאמת ש-V̇ ≤ −k1z1² − k2z2² ≤ 0. תנאי שלילי-חצי-מוגדר זה מבטיח תיאורטית ששגיאות מערכת הלולאה הסגורה מתכנסות למקור.
5. להגביל את אלגוריתם ה-PSO הבא לחקור את מרחב הפרמטרים הרב-ממדי בתוך גבולות שנבחרו כדי לשמר תנאי יציבות מבוסס ליאפונוב זה.
6. החליפו את פונקציית הסימן הלא רציפה בפונקציית רוויה רציפה, sat(x/φ), כדי לצמצם רעש בתדר גבוה בקרבת שיווי המשקל, בהתאם לפרקטיקת צעד אחורי שכבת הגבול^{הסטנדרטית 23}. הגדר φ כרוחב שכבת הגבול המשמש להחלקת מעברי הבקרה בתוך משוואה 3.
   sat(x/φ) = 1 עבור x/φ > 1; sat(x/φ) = x/φ עבור |x/φ| ≤ 1; sat(x/φ) = −1 עבור x/φ < −1 (3)
7. הערכת המהירויות הזוויתיות ממדידות מיקום על ידי יישום אלגוריתם הפרש לאחור מסונן מסדר ראשון עם תדירות חיתוך של 25 הרץ.
8. הפעילו את הרווחים הקבועים הבסיסיים המפורטים בטבלה 2 (k1 = 3.60, k2 = 1.95, λ = 2.10, ו-φ = 0.12). אל תכוונו מחדש ערכים אלו לאחר תחילת הבדיקות הרשמיות.

4. אופטימיזציה של רווח צעד אחורה באמצעות אופטימיזציה של להקת חלקיקים (PSO)

1. בצע כיוון לא מקוון לפני ההערכה בזמן אמת. אופטימיזציה של ארבעת פרמטרי הבקרה (k1, k2, λ ו-φ) תוך שימוש ב-PSO המוגדר עם 20 חלקיקים והגבלת הביצוע ל-35 איטרציות²⁴. סיימו את האופטימיזציה בסף זה כדי למנוע עומס חישובי מיותר, שכן כוונון אמפירי לא מקוון הראה שהנחיל התכנס בעקביות לערך כושר גלובלי יציב בתוך 20 האיטרציות הראשונות.
2. אתחול החלקיקים באופן אחיד בתוך גבולות החיפוש שנקבעו: k1 ∈ [2.5, 5.5], k2 ∈ [1.2, 3.0], λ ∈ [1.2, 2.8], ו-φ ∈ [0.05, 1.20].
3. הערך כל חלקיק באמצעות תרחיש אימון לא מקוון של 20 שניות. ודא שסביבת ההערכה הזו משחזרת בקפדנות את מגבלות הרוויה, רעש החיישנים, העיכוב וחוסר ההתאמה של הדיכוי שנקבעו לניסויים הרשמיים בזמן אמת.
4. נסח את פונקציית המטרה למזעור כך J = 0.55RMSE + 0.25max|α| + 0.20RMSU.
5. נרמל את שלושת רכיבי המטרה הללו באמצעות ערכי ייחוס בסיסיים קבועים, המחושבים כממוצע של חמש ריצות בסיס ראשוניות שבוצעו לפני תחילת חיפוש הלהקה. שמור על קבועי נרמול אלו בקשיחות לאורך כל תהליך האופטימיזציה.
6. הפחתת משקל האינרציה באופן ליניארי מ-0.90 ל-0.40 במהלך ריצת האופטימיזציה. קבע את שני מקדמי ההאצה ב-1.50.
7. הקצה עלות קנס קבועה מראש לכל תצורת חלקיקים שמפעילה עצירת בטיחות במהלך שלב ההערכה.
8. להוציא את הפתרון הגלובלי-מיטבי ישירות עם השלמת האיטרציה הסופית ולהגדיר אותו כבקר האופטימלי.
9. מימוש הרווחים האופטימליים המפורטים בטבלה 2 לכל ניסויי הזמן אמת הפורמליים המותאמים ל-PSO (k1 = 4.19, k2 = 2.44, λ = 1.83, ו-φ = 0.92).

5. תצורת פלטפורמת ביצוע בזמן אמת

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

6. ביצוע בדיקות מעקב שלבים

1. יש להחיל רצף צעדים של 20 שניות בחלקים המוגדרים כ-θr = 0 רדיאטור מ-0.0 שניות עד 2.0 שניות, θr = 0.50 רדיאטור מ-2.0 שניות עד 10.0 שניות, ו-θr = 0.20 רדיאטור מ-10.0 שניות עד 20.0 שניות.
2. בצע 10 ניסויים תקפים עבור בקר הבסיס ו-10 ניסויים תקפים עבור בקר מותאם ל-PSO.
3. רשמו את הזמן שחלף, זווית ייחוס, זווית הזרוע הנמדדת, שגיאת מעקב, סטיית מטוטלת ומתח בקרה ברציפות במהלך כל ניסוי.
4. הערכו את ביצועי המעבר באופן ספציפי לאחר כניסת צעד רדיונל של 0.50 שמתרחשת ב-t = 2.0 שניות.
5. חשב את זמן העלייה, המוגדר כמשך הזמן הנדרש לתגובה הזווית-זרועית הנמדדת לעבור מ-10% ל-90% מנקודת ההגדרה היעדית.
6. חשב חריגה כאחוז שבו סטיית הזווית השיא עולה על יעד 0.50 רדיו.
7. חשב את זמן ההתייצבות, שנקבע כמקרה הראשון שאחריו התגובה נשארת מוגבלת בטווח של ±2% מהיעד של 0.50 רדיאטור, למשך מינימום של 1.0 שניות.
8. סמנו כל משפט שלא עומד בתנאי ההסדר הזה לפני סימון 10.0 שניות כלא פתור. החריגו אותו לחלוטין מחישובי ממוצע זמן ההתייצבות תוך שמירה על כל המדדים האחרים.
9. חשב את שגיאת המצב היציב על ידי ממוצע שגיאת המעקב על חלון הזמן המשתרע מ-9.0 שניות ל-10.0 שניות.

7. ביצוע מבחן מעקב סינוסואידלי

1. הזרקו מסלול ייחוס סינוסואידלי בתדר מעורב של 20 שניות, הנשלט על ידי משוואה 4.
   θr(t) = 0.26sin(2π·0.10T) + 0.12sin(2π·0.30T + 0.40) רדיאטור (4)
2. השלים 10 ריצות ניסיון תקפות לכל קונפיגורציה של בקר.
3. הוצא את RMSE המעקב המלא של הניסוי, שגיאת המעקב המקסימלית המוחלטת, מתח השליטה ב-RMS, סטיית המטוטלת המקסימלית והשהיית פאזה עבור כל ניסוי בנפרד.
4. הערכת השהיית הפאזה על ידי ביצוע קורלציה צולבת בין מסלול זווית הזרוע הנמדד לבין ההתייחסות הפקודה בתוך חלון חיפוש מוגדר של ±0.50 שניות. רשום את ההשהיה כערך חיובי כאשר הפלט הנמדד עוקב זמנית אחרי אות הייחוס.

8. ביצוע מבחן דחיית הפרעה

1. שמור על התייחסות זווית הזרוע על 0 רדיאציה לאורך כל 20 השניות.
2. הזריק פולס מומנט נוסף בצד המנוע בדיוק ב-6.0 שניות, ומרח פולס זהה שני ב-12.2 שניות.
3. הגדר את משרעת ההפרעה ל-0.030 ניוטון-מטר והגדר את רוחב הפולס ל-0.12 שניות לשני אירועי ההזרקה. השלים 10 ניסויים תקפים לכל בקר.
4. חשב את שגיאת המעקב האבסולוטי השיא אך ורק בתוך חלון ה-0.50 שניות מיד לאחר כל התחלת פולס. זהה ושמור את הערך הגדול מבין השניים כשגיאת השיא המייצגת לאחר ההפרעה עבור אותו ניסוי.
5. חשב את זמן ההתאוששות, המוגדר כמשך הזמן שחלף מתחילת הפולס ועד לרגע הדגימה הראשון שבו |eθ| ≤ 0.010 רדיאטור. אין לבצע אינטרפולציה של זמני שחזור של תת-דגימות, שכן רישום הנתונים מתבצע אך ורק בפרקי זמן של 20 מילישניות.
6. אינסט את שגיאת המעקב המוחלטת וחשב את מתח ה-RMS של הבקרה באופן מצטבר לאורך כל 20 שנות הניסוי.

9. ייצוא נתונים וסיכום סטטיסטי

1. ייצא את מערכי הנתונים הגולמיים של תחום הזמן מבדיקות הצעד, הסינוסואידלי וההפרעה לקבצי נתונים נפרדים.
2. ארגן את מדדי הניסוי החוזר בפורמט מטריצה מובנה, תוך הקצאת שורה אחת לכל ניסוי ועמודה ייעודית אחת לכל מדד מוערך.
3. יש להוציא לחלוטין כל ניסוי שסומן כלא גמור או שהופסק על ידי עצירות בטיחות. החלף כל הרצה מוחרגת על ידי ביצוע ניסוי נוסף בתנאים זהים כדי להבטיח מערך נתונים סופי הכולל בדיוק 10 ניסויים תקפים לכל בקר לכל תרחיש בדיקה.
4. דווח על תוצאות הניסויים החוזרים המצטברות כסטיית תקן ממוצעת ±.
5. הערך את נורמליות ההתפלגות של המדדים שהופקו באמצעות מבחן שפירו-וילק²⁵.
6. לנתח את ההבדלים הסטטיסטיים בין הבקרים הבסיסיים והאופטימליים באמצעות מבחני t במדגם עצמאי דו-זנבי עבור נתונים שמפוזרים נורמלית. החלו את מבחן ה-U של מאן-וויטני לכל פרמטרים המציגים התפלגויות לא נורמליות.
7. קבע את סף המובהקות הסטטיסטית ב-p < 0.05.

תוצאות

אלגוריתם אופטימיזציית להקת חלקיקים (PSO) הדגים הפחתה ראשונית מהירה בפונקציית המטרה המורכבת, ואחריה תקופה של התכנסות הדרגתית. באופן ספציפי, ערך הכושר האופטימלי ירד ל-1.8757 במחזור החישוב הראשון, בעוד שהכשירות הממוצעת ללהקות ירדה מ-2.0573 ל-1.2518 בגרסה ה-35. רוב ההתכנסות הזו התרחשה במהלך 15 עד 20 האיטרציות הראשונות. מעבר לשלב זה, מסלול הפתרון הגלובלי-הטוב ביותר התכנס, מה שמעיד שהנחיל הגיע לאזור יציב בתוך מרחב החיפוש המוגדר מראש במקום לנדוד באקראי. קבוצת הרווח הסופית שנגזרה מ-PSO שמרה על שלמות מבנית של חוק הבסיס לאחור אך חילקה מחדש את המשקל היחסי בין המונחים המייצבים (טבלה 2). כפי שמודגם באיור 1, אלגוריתם האופטימיזציה אפשר דעיכה מהירה של פונקציית ההתאמה, כאשר עבר משלב חקירה ראשוני לפתרון גלובלי יציב ב-35 איטרציות. בהשוואה לקו הבסיס המכוון ידנית, הבקר האופטימלי אימץ ערכים גדולים יותר עבור k1 ו-k2, שיפוע שטח שגיאה מופחת λ, ושכבת גבול רחבה יותר φ.

ביצועי מעקב צעדים

ניתוח השוואתי חושף הבדלים ברורים בביצועים בין שני הבקרים במהלך ניסוי מעקב הצעדים, וממחיש את התוצאות החיוביות והלא אופטימליות הצפויות של פרוטוקול זה. בקרת הצעד הבסיסי לאחור מייצגת תוצאה טיפוסית תת-אופטימלית; בעוד ששמר על יציבות כוללת לאחר פקודת רדיאציה 0.50 שהוזרקה ב-2.0 שניות, הוא הציג קצב עלייה ראשוני מהיר שסטה מוקדם מהמסלול המטרה, מה שהוביל לתנודות שאריות מתמשכות ולצורך משך זמן ממושך להתכנסות. לעומת זאת, הבקר המותאם ל-PSO מציג תוצאה חיובית מייצגת, המאופיינת בגישה מהירה ומדוכאת לנקודת ההגדרה מבלי להפעיל רוויה פקודית. בעשרה ניסויים רצופים, הבקר האופטימלי הוריד את מעקב הצעדים של RMSE מ-0.065 ±-0.002 רדיאטור ל-0.050 ± 0.002 רדיאטור (בדיקת וולש t, עמוד < 0.001, טבלה 4). מדדים חולפים השתפרו גם הם: המעבר ירד מ-4.644 ± 0.604% ל-3.627 ± 0.620% (מאן–וויטני U, p = 0.009), וזמן העלייה קוצר מ-0.788 ± 0.025 שניות ל-0.548 ± 0.025 שניות (Mann–Whitney U, עמ' < 0.001). הבדלים אלו בביצועים מוצגים באיור 2 ובאיור 3, שם התגובה האופטימלית מציגה שיפור בדיכוך והתכנסות נקודת הגדרה מהירה יותר. דיוק משופר זה הושג לצד הפחתה של 22.1% במתח הבקרה של השורש-ממוצע-ריבוע (מ-0.669 ± 0.004 וולט ל-0.521 ± 0.003 וולט). במקביל, נדנוד המטוטלת המרבי הופחת ב-33.1%, וירד מ-5.812° ל-3.889° (טבלה 4). יתרה מזאת, אף שניסויי בסיס לא עמדו בקריטריון ההתייצבות של 1.0 שניות בתוך חלון ה-10.0 שניות, הבקר האופטימלי השיג התייצבות מוצלחת ב-70% מהניסויים, עם זמן התייצבות ממוצע של 1.057 ±-0.348 שניות.

ביצועי מעקב סינוסואידלי

שתי אסטרטגיות הבקרה שמרו על מעקב מוגבל תחת ההתייחסות הסינוסואידלית בתדר מעורב לאורך כל תקופת הבדיקה של 20 שניות. למרות שהפער בביצועים היה פחות בולט ויזואלית מאשר במבחני התגובה הצעדית, תוצאות סטטיסטיות בניסויים חוזרים הצביעו על שיפור בביצועים בגישה האופטימלית. אות הבקרה הבסיסי הציג טווח שגיאות מעקב רחב יותר והיה מפגר באופן ניכר מאחורי צורת הגל המטרה. לעומת זאת, הבקר האופטימלי עקב מקרוב יותר אחרי נקודת ההגדרה במהלך שינויים בתדרים נמוכים וגבוהים, מה שיצר התפלגות מדויקת יותר של שגיאות שיא לאורך מסלול הייחוס. כתוצאה מכך, ה-RMSE בניסוי מלא הופחת מ-0.061 ±0.003 רדאם ל-0.040 ± 0.003 רדיאם. כפי שמפורט בטבלה 5, הבקר האופטימלי גם סייע בהפחתת מתח ה-RMS של הבקרה (מ-0.557 ± 0.002 וולט ל-0.402 ± 0.002 וולט) והגביל את ההעתקה המקסימלית של המטוטלת מ-2.156 ± 0.079° ל-1.464 ± 0.061°, מה שמעיד על שיפור בסנכרון פאזה ויעילות אנרגטית (p < 0.001, טבלה 5; איור 4).

דחיית הפרעה ועמידות חוצת תרחישים

הזרקת פולסי מומנט מצטברים ב-6.0 שניות ו-12.2 שניות גרמה לסטיות זוויות זרוע גדולות יותר ולסטיות רחבות יותר של מטוטלת בתגובת בקר הבסיס. בעוד שהבקר האופטימלי גם זיהה את ההפרעות מיד, הנסיעות לאחר ההפרעה היו נמוכות יותר. תהליך ההתאוששות לעבר שיווי המשקל האפס הציג מגמה איכותנית מהירה יותר (מ-0.014 שניות ל-0.000 שניות); עם זאת, מדד זה לא הגיע למשמעות סטטיסטית (p = 0.078), בעיקר משום שרזולוציית רישום הנתונים של 20 מילישניות יצרה אפקט רצפה מדיד. במהלך ניסויים חוזרים, הפרמטריזציה האופטימלית הפחיתה את שגיאת השיא לאחר ההפרעה. אינטגרל השגיאה האבסולוטית (IAE) שחושב לאורך 20 שניות הופחת מ-0.377 ± 0.012 ל-0.248 ± 0.015, לצד ירידה מקבילה במתח ה-RMS של הבקרה מ-0.656 ± 0.003 וולט ל-0.508 ± 0.005 וולט. המשרעת המקסימלית של תנודת המטוטלת הוחלשה באופן דומה מ-11.456 ± 0.117° ל-7.561 ± 0.063° (טבלה 6). דירוג הבקרים נשאר עקבי בכל שלושת פרופילי הבדיקה. אסטרטגיית הצעד האחורי המותאמת ל-PSO הניבה שגיאות חולפות נמוכות יותר, הפחיתה אי-דיוקים במצב יציב ומזעור יציאות פסיביות של מטוטלת (p < 0.001, טבלה 6; איור 5) בלי לדרוש הגדלת מאמץ המפעיל הכולל. ההתפלגויות בין ניסויים של מדדי הביצועים הראשיים מסוכמים באיור 6.

זמינות נתונים:

מאגר הנתונים התומך ששימש ליצירת האיורים והטבלאות, כולל נתונים מייצגים, מדדים ברמת ניסוי, זרעים אקראיים מוגדרים מראש ורשומות כיוונון PSO, הופקד ב-Zenodo וזמין לציבור תחת DOI 10.5281/zenodo.20395966.

איור 1: התכנסות אופטימיזציה של עדר חלקיקים בכיוון בקר. (א) ערכי הכושר הטובים ביותר של הלהקה ב-35 איטרציות; (ב) ערכי כושר ממוצעים לכל נקודת זמן. המטרה המורכבת התפרקה במהירות בשלב החיפוש הראשוני אך האטה בהדרגה ככל שהתקרבה לאזור יציב של חוק בקרת הצעד האחורי המותאם ל-PSO. אנא לחצו כאן כדי לצפות בגרסה מוגדלת של הדמות הזו.

איור 2: תגובות מייצגות למעקב צעדים של בקר הבסיס ובקר החזרה המותאם ל-PSO. (א) הפניה בזווית זרוע פקודה ותגובת זווית זרוע נמדדת לרצף צעדים מלא של 20 שניות; (ב) תצוגה מורחבת של התגובה מיד לאחר צעד רדיאציה של 0.50 רדיאציה שהוחל ב-2.0 שניות; (ג) תצוגה מורחבת של התגובה לאחר מקטע הפיקוד הרדיאטיבי השני 0.20 מתחיל ב-10.0 שניות. לאחר אופטימיזציה, הבקר המותאם ל-PSO הראה זמן עלייה מהיר יותר, שגיאת שיא נמוכה יותר, ושיכוך טוב יותר בהשוואה לבקר הבסיסי. אנא לחצו כאן כדי לצפות בגרסה מוגדלת של הדמות הזו.

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

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

איור 5: תגובות ייצוגיות להפרעה-דחייה תחת רגולציה של אפס ייחוס. (א) שגיאת מעקב בזווית הזרוע במהלך מבחן דחיית ההפרעה, כולל הפרעות מומנט-פולס ב-6.0 שניות ו-12.2 שניות; (ב) זווית סטיית המטוטלת במהלך אותו ניסוי; (ג) אות מתח בקרה במהלך דחיית הפרעה. בהשוואה לבקר הבסיסי, בקרה המותאמת ל-PSO הראה סטייה מקסימלית נמוכה יותר לאחר ההפרעה והתאוששות מהירה יותר לכיוון נקודת ההגדרה. אנא לחצו כאן כדי לצפות בגרסה מוגדלת של הדמות הזו.

איור 6: התפלגות בין ניסויים של מדדי הביצועים הראשיים בכל שלושת מצבי הבדיקה. (א) RMSE במעקב צעדים; (ב) זמן עלייה במעקב צעדים; (ג) RMSE במעקב סינוסואידלי; (D) השהיית פאזה סינוסואידלית; (ה) שגיאת מעקב שיא לאחר הפרעה; (ו) הזזת המטוטלת המקסימלית בתנאי דחיית הפרעה. כל פאנל מציג את התפלגות 10 ניסויים תקפים תחת כל תנאי בקר. הנטייה והפיזור המרכזיים של מערכות בקרת צעד אחורה-בסיסי ומותאמות ל-PSO משווים ישירות. אנא לחצו כאן כדי לצפות בגרסה מוגדלת של הדמות הזו.

פרמטר	סמל	ערך	יחידה
אורך זרוע סיבובית	Lr	0.215	m
אורך מרכז המסה של המטוטלת	Lp	0.168	m
מסת זרוע סיבובית	מרר	0.254	ק"ג
מסת מטוטלת	MP	0.097	ק"ג
דממה צמיגית בזרוע	br	0.0031	N·m·s/rad
הפחתה צמיגית של מטוטלת	בפ	0.0018	N·m·s/rad
קבוע מומנט המנוע	Kt	0.053	N·m/A
רזולוציית המקודד	Nenc	4096	רוזנים/רף
תקופת עדכון הבקר	Ts	0.002	s

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

בקר	רווח	ערך	תיאור
צעד אחורה בקו הבסיס	K1	3.6	הגיין ייצוב ראשוני מכוון ידנית
צעד אחורה בקו הבסיס	K2	1.95	הגבר שליטה וירטואלי מכוון ידנית
צעד אחורה בקו הבסיס	λ	2.1	שיפוע שטח שגיאה
צעד אחורה בקו הבסיס	φ	0.12	רוחב שכבת הגבול
צעד אחורה מותאם ל-PSO	K1	4.19	החלקיק הטוב ביותר בעולם לאחר 35 איטרציות
צעד אחורה מותאם ל-PSO	K2	2.44	החלקיק הטוב ביותר בעולם לאחר 35 איטרציות
צעד אחורה מותאם ל-PSO	λ	1.83	החלקיק הטוב ביותר בעולם לאחר 35 איטרציות
צעד אחורה מותאם ל-PSO	φ	0.92	החלקיק הטוב ביותר בעולם לאחר 35 איטרציות

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

פריט	מפרט טכני
מחשב מטרה בזמן אמת	מחשב תעשייתי, מעבד אינטל i7, זיכרון RAM 16GB
מצב ביצוע	סימולציה בזמן אמת בשלב קבוע
תקופת עדכון הבקר	2 מילישניות
מרווח ייצוא אות	20 מילישניות
ממשק אותות	קלט ומפלט אנלוגי למקודד, ±10 וולט
תצורת המפעל	מודל מטוטלת סיבובית לא ליניארית מסוג פורוטה
השפעות לא אידיאליות מוזרקות	רעש מדידה, עיכוב ערוץ מטוטלת וחוסר התאמה בהחלקה
רעש מדידה	0.003 רדיאטור בזווית זרוע; 0.004 רדיאטור בזווית מטוטלת
עיכוב ערוץ מטוטלת	4 מילישניות
אי-התאמה בשיכוך הצמחים	+8% יחסית למודל הבקר הנומינלי
רוויית פקודה	±10 וולט
קריטריוני עצירות בטיחות	|θ| > 0.70 רדיטור, |α| > 0.35 רדיאטור, או רוויה > 100 מילישניות
קבלת מיקום מוקדם	|α| ≤ 0.05 רדיאטור לפחות ל-1.0 שניות לפני תחילת הניסיון
פרופיל ייחוס צעדים	0 רדיונינג (0.0-2.0 שניות), 0.50 רדיונר (2.0-10.0 שניות), 0.20 רדיונר (10.0-20.0 שניות)
פרופיל ייחוס סינוסואידלי	0.26 SIN (2π·0.20T) + 0.12 SIN (2π·0.30T + 0.40) רדיור
פרוטוקול הפרעה	שני פולסים של מומנט מצטבר ב-6.0 שניות ו-12.2 שניות
משרעת ורוחב הפרעה	0.030 ניוטון-מטר, 0.12 שנייה
ניסויים תקפים לכל בקר ותרחיש	10

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

מטריקה	צעד אחורה בקו בסיס n=10	צעד אחורה מותאם ל-PSO n=10	ערך p	מבחן סטטיסטי
מעקב RMSE (רדיאטור)	0.065 ± 0.002	0.050 ± 0.002	<0.001	מבחן t וולש
שגיאת מצב יציב, 9.0–10.0 שניות (רדיאט)	-0.007 ± 0.001	-0.004 ± 0.001	<0.001	מאן–וויטני U
אוברשוט (%)	4.644 ± 0.604	3.627 ± 0.620	0.009	מאן–וויטני U
זמן עלייה	0.788 ± 0.025	0.548 ± 0.025	<0.001	מאן–וויטני U
זמן ההתיישבות	אין משפטים תקפים שהוסדרו (0/10)	1.057 ± 0.348 (7 מתוך 10 משפטים תקפים)	NA	לא בהשוואה
Control RMS (V)	0.669 ± 0.004	0.521 ± 0.003	<0.001	מבחן t וולש
סטיית מטוטלת מרבית (deg)	5.812 ± 0.156	3.889 ± 0.109	<0.001	מבחן t וולש

טבלה 4: מדדי ביצוע של מעקב שלבים בין 10 ניסויים תקפים לכל בקר. טבלה זו מציגה את הסטטיסטיקות הסיכום של ניסוי מעקב הצעדים, הכוללות מעקב אחר RMSE, שגיאת מצב יציב, חריגה מהעבר, זמן עלייה, מהירות התכנסות נקודות הגדרה, מתח RMS בקרה, וסטיית מטוטלת מרבית. הערכים מדווחים כממוצעים ± סטיות תקן אלא אם צוין אחרת.

צעד אחורה מותאם ל-PSO n=10	ערך p	מבחן סטטיסטי
0.040 ± 0.003	<0.001	מבחן t וולש
0.088 ± 0.006	<0.001	מבחן t וולש
0.232 ± 0.019	<0.001	מבחן t וולש
0.402 ± 0.002	<0.001	מבחן t וולש
1.464 ± 0.061	<0.001	מבחן t וולש

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

ערך p	מבחן סטטיסטי
<0.001	מבחן t וולש
0.078	מאן–וויטני U
<0.001	מבחן t וולש
<0.001	מבחן t וולש
<0.001	מבחן t וולש

טבלה 6: מדדי ביצועי דחיית הפרעות ב-10 ניסויים תקפים לכל בקר. נתוני ניסוי חוזר מניסויי דחיית הפרעה כוללים שגיאת מעקב שיא לאחר ההפרעה, זמן התאוששות, אינטגרל השגיאה המוחלטת (IAE), מתח RMS בקרה, וסטיית מטוטלת מרבית. כפי שנקבע על ידי הגדרות ייצוא האות, זמן ההתאוששות המדיד מוגבל על ידי תקופת רישום של 20 מילישניות.

דיון

המטרה העיקרית של מחקר זה אינה רק להכריז על בקר אחד כעליון על אחר בתנאים מבודדים, אלא להראות כי מסגרת אופטימיזציית הרווח מבוססת PSO שנבדקה כאן יכולה לשפר את ביצועי בקר החזרה לאחור בתרחישים לא אידיאליים מוגדרים. המטוטלת ההפוכה הסיבובית משמשת כנקודת ייחוס מצוינת להערכה זו בשל מאפייניה הלא-ליניאריים מאוד, פאזה לא מינימלית וחסרת הפעלה²⁶. מעבר למדדי ביצועים, היושרה הטכנית של פרוטוקול זה תלויה באפיון מכוון של סביבת הסימולציה בזמן אמת. שלב יישום קריטי כולל הזרקה מבוקרת של חוסר התאמה בשיכוך של 8% והשהיית חיישן של 4 מילישניות, שכן השפעות לא אידיאליות אלו מאפשרות לאלגוריתם PSO למקסם רווחים תחת אותן הפרעות ששימשו בהערכה הפורמלית ולא רק בתנאים אידיאליים.

על ידי בדיקת אסטרטגיות הבקרה באמצעות פרוטוקולי מעקב צעדים, מעקב סינוסואידלי ודחיית הפרעות תחת הגדרת ביצוע אחידה של שלבים קבועים, תוצאות הסימולציה הראו מגמת ביצועים עקבית. ראוי לציין כי השיפורים שנצפו הן בתגובות צעדיות והן בתגובות סינוסואידליות—בפרט, הפחתת שגיאות מעקב, צמצום חריגה והפחתת השהיית פאזה—הושגו לצד ירידה מדידה במתח הבקרה של RMS. דבר זה מצביע על כך שתהליך האופטימיזציה לא שיפר רק את דיוק המעקב על ידי רוויה של יחידת ההנעה; במקום זאת, הוא חילק מחדש את מאמצי הבקרה בתוך מסגרת הצעד האחורי הקיימת המבוססת על ליאפונוב. העקביות של חלוקת פרמטרים זו מושפעת משמעותית מקריטריוני קבלת הפריפוזיציה המוגדרים בפרוטוקול שלב 2.4. הבטחה קפדנית שהמטוטלת נשארת בתוך חלון שיווי המשקל המיועד למשך זמן ממושך היא חיונית; אחרת, טרנזיינטים ראשוניים לא סדירים עלולים להוביל לקיפאון מוקדם בפונקציית עלות או להתנהגות עדר לא מתכנסת במהלך שלב הכוונון הלא מקוון הבא. במערכות מכניות עם פעילות חלשה, חלוקה מחדש זו חשובה כי היא שולטת כיצד מאמץ הבקרה מוקצה לניהול הקישור בין הקואורדינטה האקטיבית (זרוע סיבובית) לקואורדינטה הפסיבית (מטוטלת)²⁷.

יתרה מזאת, ההפחתה הסימולטנית הן בשגיאת המעקב בזווית הזרוע והן ביציאת המטוטלת הפנימית מדגישה יתרון של פרוטוקול הכיוון המוצע. בשליטה על מערכות עם פעילות נמוכה, השגת מעקב קואורדינטות מדויק על חשבון מצבים פנימיים תנודתיים נחשבת בדרך כלל לעיצוב תת-אופטימלי או אפילו כושל. חוק הבקרה האופטימלי הפחית את אמפליטודת התגובה של המטוטלת הלא מופעלת ב-33.1%, תוך שיפור איכות המעקב בו זמנית. שיפור משולב זה מרמז כי אלגוריתם PSO זיהה שילוב רווח שמאזן טוב יותר את הקישור הלא-ליניארי בין הקואורדינטות האקטיביות והפסיביות, במקום להסתמך על משוב בעל רווח גבוה שיכול להגביר תנודות טפיליות. ביצועים מאוזנים אלו הם אינדיקטור אינפורמטיבי יותר ליעילות הבקרה מאשר מעקב מסלול טהור בלבד, שכן הם מצביעים על כך שדינמיקת המצב הפנימית הייתה מוגבלת כראוי במהלך תהליכי המרת האנרגיה²⁸. באופן דומה, תוצאות ההפרעה-דחייה תומכות בעמידות המעשית של הרווחים האופטימליים בתוך תנאי הבדיקה המדומים. בעוד שמרווח הרישום של 20 מילישניות יצר אפקט רצפה שהסתיר זמני התאוששות מדויקים של תת-דגימה, הפחתות מקבילות בשגיאת השיא לאחר הפרעה, השגיאה המוחלטת המצטברת וההעתקה המקסימלית של המטוטלת מצביעות על כך שהפרמטרים האופטימליים הגבילו את מידת הזרקות הפולס החיצוניות שמתפשטות בדינמיקת המערכת²⁹.

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

למרות התוצאות המבטיחות הללו, המתודולוגיה נשארת מוגבלת על ידי ההפשטות הטבועות בסביבות סימולציה. הפרוטוקול המיושם לקח בחשבון עיכובים בדגימה, רוויה פקודתית, חוסר התאמה בדיכוך ורעש מדידה — גורמים שלעיתים קרובות מוזנחים באינטגרציה נומרית לא מקוונת טהורה. עם זאת, המחקר הנוכחי לא כלל ניסויים פיזיים של מטוטלת סיבובית הפוכה, אימות חומרה בתוך הלולאה, או בדיקות פריסה פיזיות. לכן, יש לפרש את הממצאים כהוכחה מבוססת סימולציה להשפעת הביצועים של כיוונון הגבר מבוסס PSO בתנאי סימולציה קבועים ולא אידיאליים מוגדרים. עבודות עתידיות עשויות להרחיב את מערך אי-הוודאות המודגם על ידי שילוב חיכוך, מפעיל, ציות וכימות נוספים בתוך הסימולציה, לפני שכל מחקר חומרה נפרד ייחשב 30,31,32.

מחקר זה קובע פרוטוקול סימולציה מהיר לסימולציית אב-טיפוס בקרה להערכת אופטימיזציה של רווח מבוסס PSO לבקרת אחורה. תוצאות כמותיות מ-10 ניסויים חוזרים מראות שהבקר האופטימלי מפחית את קצב ה-RMSE במעקב צעדים מ-0.065 ל-0.050 רדיאטור, ומשיג ירידה של 22.1% במאמץ בקרת ה-RMS בהשוואה לקו הבסיס המכוון ידנית. ההחלשה של 33.1% ביציאות מטוטלת מצביעה על שיפור ביציבות המצב הפנימי בתנאים הלא אידיאליים המדומים. בעוד שסביבת הסימולציה הזו משלבת רעש, עיכוב, חוסר התאמה בדיכוך ורוויה, היא אינה משחזרת במלואה את החיכוך המורכב של קולון, רצועות המתות של המנוע, עמידות מבנית או קוונטיזציה של חיישנים הנמצאים במערכות פיזיקליות. העבודה העתידית תתמקד בהרחבת מודל הסימולציה לכלול השפעות לא אידיאליות נוספות ובבדיקת האם אותו דפוס דירוג נשאר יציב בתנאי אי-ודאות מדומים רחבים יותר.

חומרים

רשימת החומרים שנעשה בהם שימוש במאמר זה

שם	חברה	מספר קטלוג	הערות
High-performance Workstation (Windows 11 Pro)	Various / Custom Build	N/A	מארח מכונה לסימולציה בזמן אמת עם צעד קבוע ופוסט-עיבוד.
MATLAB (Version R2024a)	MathWorks	https://www.mathworks.com/products/matlab.html	דוגמנות צמח לא לינארית, קידוד בקר, יצוא נתונים וניהול פרמטרים.
Simulink (Version R2024a)	MathWorks	https://www.mathworks.com/products/simulink.html	בניית מודל דיאגרמת בלוק עבור מטוטלת מסוג Furuta וביצוע בקר.
Simulink Desktop Real-Time (Version R2024a)	MathWorks	https://www.mathworks.com/products/simulink-desktop-real-time.html	גרעין בזמן אמת עם צעד קבוע עבור ביצוע דגם הבקרה בשולחן העבודה.
Python (Version 3.11)	Python Software Foundation	https://www.python.org/	עיבוד נתונים משני, טיפול סטטיסטי והכנת תרשים.
NumPy (Version 1.26)	NumPy Developers	https://numpy.org/	פעולות מערך מספרי עבור נתוני ניסוי מיוצאים.
pandas (Version 2.2)	pandas Developers	https://pandas.pydata.org/	ארגון נתונים מניסויים חוזרים והכנת טבלת סיכום.
SciPy (Version 1.13)	SciPy Developers	https://scipy.org/	בדיקות סטטיסטיות וכלי ניתוח אותות.
Matplotlib (Version 3.8)	Matplotlib Developers	https://matplotlib.org/	יצירת גרפים עבור התכנסות, מעקב ודמויות חלוקה.