הצגת המקור של הדף "הנחיית Q"

'''הנחיית Q''' היא שיטת הנחיה שנעשה בה שימוש ב[[טיל בליסטי|טילים בליסטיים]], שמסלולם מורכב משלב ממונע קצר יחסית, שבה מערכת ההנעה של הטיל פועלת, ולאחריו שלב בליסטי אשר במהלכו הטיל מתכוון אל המטרה בהשפעת ה[[כבידה]] ([[טיל שיוט|טילי שיוט]] משתמשים בשיטות הנחיה שונות). היעד של הנחיית Q הוא לפגוע במטרה ספציפית בזמן מסוים, שכן בגלל [[סיבוב כדור הארץ]] פגיעה במטרה בעלת [[קואורדינטות]] מסוימות פירושה פגיעה בה במיקום ובזמן מסוים.

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

== יישומים מוקדמים ==
בזמן שהנחיית Q פותחה השיטה המתחרה העיקרית נקראה הנחיית דלתא. לפי מקנזי, הטילים [[טיטן 2|טיטאן]], מספר גרסאות של [[אטלס (סדרת טילים)|אטלס]], ומיניוטמן 1 ו-2 השתמשו בהנחיית דלתא, בעוד שבהנחיית Q נעשה שימוש ב[[תור (משגר)|ת'ור]] ו[[פולריס]]. מניתוח סרטוני שיגור, עולה, כי מספר [[טיל בליסטי בין יבשתי|טילים בליסטיים בין יבשתיים]] סובייטים מוקדמים השתמשו בהנחיית דלתא.

==רקע: הנחיית דלתא==
הנחיית דלתא מבוססת על הצמדות למסלול מתוכנן מראש, שמפותח לפני השיגור באמצעות מחשבים קרקעיים ונשמר במערכת ההנחיה של הטיל. במהלך המעוף, מערכת ההנחיה של הטיל מנטרת סטיות מן מסלול הייחוס, והמסלול האמיתי ממודל מתמטית כפיתוח [[טור טיילור]] מסביב למסלול הייחוס. בניסוח מתמטי, מערכת ההנחיה דלתא עושה שימוש בקשר פונקציונלי מסוים מהצורה <math>F(r,v,r_T) = 0</math> (הקשר הפונקציונלי מתאפס על מסלול הייחוס). היעדר פרמטר זמן (זמן מוגדר לפגיעת הטיל) הוא מאפיין אינהרנטי של מערכת ההנחיה דלתא (קיימים מסלולי ייחוס אפשריים רבים). רגע כיבוי המנוע ( cut-off moment) מתוכנן להיות הרגע בו מסלול הטיל חותך את מסלול הייחוס. מערכת הבקרה של הטיל מנסה לאפס את סכום האיברים הליניאריים של הביטוי הזה, כלומר להביא את הטיל לקונפיגורצית מיקום-מהירות כזאת שמסלול הטיל החל מנקודת החיתוך יעבור דרך המטרה (האיברים מסדרים גבוהים יותר זניחים). מסיבה זו, להנחיית דלתא לעיתים מתייחסים כ-"לעוף לאורך הכבל" כאשר הכבל (הדמיוני) מתייחס למסלול הייחוס.

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

==המהירות הדרושה להשגה (Velocity to be gained)==
בזמן נתון t ובעבור מיקום נתון r של כלי הרכב, וקטור המהירות המתאים V<sub>c</sub> מוגדר כדלהלן: אם לכלי הרכב הייתה המהירות V<sub>c</sub> ומערכת ההנעה הייתה נכבית, אז הטיל יגיע למטרה הרצויה בזמן הרצוי תחת השפעת הכבידה.

המהירות האמיתית מסומנת ב-V<sub>m</sub> והטיל נתון לתאוצה בהשפעת הכבידה g ובהשפעת המנועים a<sub>t</sub>. המהירות הדרושה להשגה המתאימה מוגדרת כהפרש בין V<sub>c</sub> ו- V<sub>m</sub>:

<math>V_{TBG} = V_c - V_m</math>

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

אסטרטגיית הנחיה פשוטה היא להפעיל תאוצה (דחף מנוע) תמיד בכיוון של VTBG. לזה יהיה את האפקט של קירוב המהירות האמיתית ל- V<sub>c</sub>. כאשר שתי המהירויות נעשות שוות (ו- VTBG שווה לאפס), זה הזמן לכבות את המנועים, מכיוון שהטיל לפי ההגדרה מסוגל להגיע למטרה הרצויה בזמן הרצוי לבדו.

מה שנותר הוא איך לחשב את VTBG בקלות מהמידע על הטיל הנגיש למחשב ההנחיה (מיקומו, מהירותו ותאוצתו).

==מטריצת Q ==
במשוואה הדיפרנציאלית הבאה נעשה שימוש כדי לחשב את המהירות הדרושה להשגה המתאימה ברגע נתון:

:<math> \frac{dV_{TBG}}{dt} = -a_T - Q V_{TBG}</math>

כאשר מטריצת Q מוגדרת כ-: <math>Q = \left. \frac{\partial V_c}{\partial r}\right|_{r_T,t_f}</math>.

כאשר Q היא [[מטריצה]] מסדר 3 על 3, סימטרית, ומשתנה בזמן. (הקו האנכי מתייחס לעובדה שהנגזרות החלקיות חייבות להיות מחושבות ביחס למיקום מטרה נתון rT וזמן מעוף חופשי tf). חישוב המטריצה ואיבריה מתבצע באמצעות שיקולים אנליטיים (חישוב שטח גזרה אליפטית) וגאומטריים על חיתוך של [[אליפסה]] ומעגל בעלי מוקד משותף (מרכז כדור הארץ) ומשיקולים פיזיקליים כמו שימור [[תנע זוויתי]] ([[החוק השני של קפלר]]) ושימור אנרגיה.

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

הבעיה היא כיצד לחשב את ערכי מטריצת Q בכל נקודה במרחב בהינתן זמן המעוף <math>t_f</math> ומיקום המטרה <math>r_T</math>.

=== הוכחת המשוואה הדיפרנציאלית ===

ראשית נרשום <math>V_{{TBG}} = V_C - V</math>. תחילה נוכיח את המקרה בו <math>a_T = 0</math>, כלומר שהתאוצה הקווית המורגשת על ידי [[מד תאוצה|מדי התאוצה]] היא 0 ([[נפילה חופשית]] ב[[שדה כבידה]]). נפרק את הנגזרת <math>\frac{dV_{TBG}}{dt}</math>:

<math>\frac{dV_{TBG}}{dt} = \frac{{d(V_C - V)}}{{dt}} = \frac{dV_C}{dt} - \frac{dV}{dt} = V*Q - g</math>.

במעבר האחרון השתמשנו בהגדרה של מטריצת Q ובעובדה ש-<math>\frac{dV}{dt}</math> שווה לוקטור [[תאוצת הכובד]] הרגעית g. מן ההגדרה של מטריצת Q נובעת הזהות : <math>g = Q*V_C</math>. לכן:

<math>\frac{dV_{{TBG}}}{dt} = V*Q - V_C*Q = -V_{{TBG}}*Q</math>.

במקרה שהתאוצה הנמדדת שונה מאפס יש פשוט להוסיף רכיב ל-<math>\frac{dV}{dt}</math> נוסף ל-g וכך מתקבלת המשוואה הדיפרניצאלית לעיל.

== חישוב מטריצת Q עבור תנועה פרבולית דו־ממדית ==
לאחר סיום השלב הממונע המהירות בציר x קבועה לכן ברור ש- <math>X = V_x*t_f</math> כאשר X הוא המרחק האופקי למטרה ומכאן <math>dV_x/dx</math> שווה ל- <math>1/t_f</math>. התנועות בציר x וציר y בלתי תלויות ולפיכך האיברים לאורך האלכסון במקומות (1,2) ו- (2,1) הם 0. לצורך חישוב האיבר dVy/dy נבצע את הניסוי המחשבתי הבא: נניח כי משחררים ממנוחה כדור מגובה y אל הרצפה והוא נוחת על הרצפה כעבור זמן <math>t_f</math> . כעת נניח כי לוקחים כדור נוסף ומשחררים אותו בגובה נמוך יותר ב- ∆y ומשחררים אותו בו זמנית ביחד עם הכדור הראשון. ברור כי יש לשחרר אותו עם מעט מהירות אנכית מעלה (או מטה ∆y שלילי) אם על מנת ששניהם יפגעו בקרקע באותו זמן. מרגע שחרורם על שני הכדורים פועלת תאוצה אנכית זהה מטה g ולכן הפרש המהירויות או המהירות היחסית שלהם נותר קבוע. לפיכך על מנת ששניהם ינחתו באותו זמן צריך להתקיים ש- Vy או Vy ∆ יקיים 
∆y = tf Vy*. לפיכך האיבר (2,2) במטריצה גם הוא 1/tf. לפיכך מטריצת Q המתאימה היא: 

<math>\begin{vmatrix} 1/t_f & 0 \\ 0 & 1/t_f \end{vmatrix}</math> כאשר tf = T-t ׁׁ(T הוא הזמן בו הטיל יגיע אל המטרה ו-t הוא הזמן הנוכחי)
והמשוואה הדיפרנציאלית המתארת את השתנות המהירות הדרושה להשגה היא:

<math>\frac{dV_{TBG}}{dt} = -a_T - \begin{vmatrix} 1/t_f & 0 \\ 0 & 1/t_f \end{vmatrix}V_{TBG}</math>.

== חישוב מטריצת Q במקרה הכללי (של שדה היפוך ריבועי) ==
{{בבדיקה}}
=== פיתוח מתמטי ===
בחלק זה נציע [[אלגוריתם]] לחישוב איברי מטריצת Q במקרה הכללי של שדה היפוך ריבועי (Inverse square law).
נסמן את שתי שמורות התנועה (ה[[אנרגיה]] הכוללת וה[[תנע זוויתי|תנע הזוויתי]]) ב-E ו-L. מתקיים:
<math>E_k(R) = E + GMm/R</math>. כמו כן מתקיים בנקודות של המרחק המרבי והמינימלי מהמוקד: <math>p(R)^2/2m= L^2/2mR^2 = E_k(R)</math> . כך נקבל [[משוואה ריבועית]] שפתרונותיה הם המרחק המרבי <math>R_1</math> והמרחק הקטן ביותר <math>R_2</math> מהמוקד:
.<math>ER^2+ GMmR-L^2/2m=0</math> 
פתרונות המשוואה הם:
ה[[אקסצנטריות (מתמטיקה)|אקסצנטריות]] היא : <math>\epsilon = \frac {{R_1 - R_2}}{{R_1 + R_2}} = \sqrt {{1 + 4EL^2/2m(GMm)^2}}</math> .

והלטוס רקטום (Latus Rectum) הוא: <math>l = a(1 - \epsilon^2) = 1/2(R_1 + R_2)(1 - \epsilon^2) = \frac {{L^2}}{{GMm^2}}</math> 

'''חישוב שטח גזרה אליפטית''':

משוואת אליפסה בקואורדינאטות פולריות היא: <math>r(\theta) = \frac{{l}}{{1 - \epsilon\cos\theta}}</math>. לכן שטח גזרה אליפטית הוא (פתרון האינטגרל מסובך למדי ונציג כאן רק תשובה סופית):

<math>\frac{{l^2}}{{2}} \int \frac {{1}}{{(1 - \epsilon\cos\theta)^2}}d\theta </math> = <math>\frac{{l^2}}{{2}}</math> (<math>\frac {{\epsilon\sin\theta}}{{(\epsilon^2-1)(\epsilon\cos\theta - 1)}}</math> + <math>2\tanh^{{-1}}))\frac {{(\epsilon + 1)(\tan\theta/2)}}{{\sqrt {{\epsilon^2 - 1}}}})/(\epsilon^2 - 1)^{{3/2}}</math>.

'''תיאור האלגוריתם''':

E ו-L במיקום הנתון של הטיל (עבורו יש לחשב את מטריצת Q) <math>r_0</math>, מגדירים באופן חד חד ערכי את וקטור מהירות הטיל המתאים (גודלו וכיוונו). שני התנאים שהמסלול האליפטי חייב לעבור דרך מיקום המטרה <math>r_T</math>, וזאת כעבור זמן <math>t_f</math>, מאפשרים לבנות מערכת של שתי משוואות עם שני נעלמים E ו-L, אותם אפשר למצוא ובכך להגדיר את מהירות הטיל המתאימה <math>v_C</math> בהתאם לכללים הבאים: <math>(v_r,v_\theta) = (\sqrt {{2(E/m + GM/R) - (L/mR)^2}},L/mR) </math>

ממשוואת האליפסה בקואורדינטות פולריות ניתן לחלץ את <math>cos\theta_0</math>, כאשר <math>\theta_0</math> הוא הזווית עם הציר הראשי של האליפסה של [[וקטור מיקום]] הטיל, כך (שכן המרחק למוקד שהוא מרכז כדור הארץ ידוע): <math>cos\theta_0 = 1/\epsilon(1 - l/R)</math>.
הזווית בין וקטור מיקום הטיל <math>r_0</math> לוקטור מיקום המטרה <math>r_T</math> נתונה ונסמנה <math>\vartriangle\theta</math>, ולכן ניתן להביע את <math>cos(\theta_0 + \vartriangle\theta)</math> בהתאם לזהות הטריגונומטרית: <math>cos(\theta_0 + \vartriangle\theta) = \cos\theta_0cos\vartriangle\theta - \sin\theta_0\sin\vartriangle\theta </math>. הצבת <math>\epsilon</math> ו-l בביטויים ל-<math>cos\theta_0</math> ו-<math>sin\theta_0</math> ושימוש בעובדה ש-<math>R_T = R_E</math> (מרחק המטרה מהמוקד הוא [[רדיוס]] כדור הארץ) ובביטוי לאליפסה בקואורדינטות פולריות מאפשר להרכיב משוואה אחת עם E ו-L, ולחלץ ממנה את L כתלות ב-E. לאחר מכן ניתן להיעזר בתנאי השני על זמן המעוף <math>t_f</math> כדי למצוא את E ולאחר מכן למצוא את L, ובכך למצוא את וקטור המהירות המתאים <math>v_C = (v_r,v_\theta)</math>.

לא נפרט כאן את האלגברה, אך בתום תהליך אלגברי ארוך הביטוי המתקבל ל-E כתלות ב-L מהתנאי הראשון הוא:

<math>E =\frac {{(GM)^2m^3}}{{2L^2}}(\cos\vartriangle\theta - 1 + L^2/GMm^2(1/R_E - \cos\vartriangle\theta/R) + (1 - L^2/GMm^2R)^2)/((sin\vartriangle\theta)^2 - 1)</math>.

זהו למעשה התנאי לקשר בין E ו-L הדרוש כדי שמסלול פלנטרי בהשפעת כבידה (אליפטי) יעבור דרך שתי נקודות נתונות (כשהמוקד נתון). 

כדי להכניס את זמן המעוף <math>t_f</math> למשוואות יש לחלק את שטח הגזרה האליפטית בקצב כיסוי השטח של הרדיוס וקטור (שכן לפי [[החוק השני של קפלר]] הרדיוס וקטור מכסה שטחים בקצב קבוע). קצב כיסוי השטח הוא <math>v_\theta r/2 = L/2m</math>. לפיכך מתקיים: <math>2 S m/L = t_f</math>. 

במקרה של טילים בליסטיים ניתן להניח בקירוב טוב <math>R = R_E</math> (החלק הממונע של מעוף הטיל מתרחש ברובו כשהטיל עדיין קרוב מאוד לפני הקרקע), ובכך מתקבל <math>\cos(\theta_0 + \vartriangle\theta) = \cos\theta_0</math> לכן: <math>\theta_0 + \vartriangle\theta = -\theta_0 = \vartriangle\theta/2 </math> (שכן האליפסה סימטרית מסביב לציר הראשי שלה). לכן במקרה של טילים בליסטיים הביטוי לאינטגרל בין הגבולות <math>\theta_0</math> ו-<math>\theta_0 + \vartriangle\theta </math> מתפשט ל-:

<math>\frac{l^2}{2} (\frac {{2\epsilon \sin\vartriangle\theta/2}}{{(\epsilon^2 - 1)(\epsilon \cos\vartriangle\theta/2 - 1)}} + 4\tanh^{{-1}}))\frac {{(\epsilon + 1)(tan\vartriangle\theta/4)}}{{\sqrt {{\epsilon^2 - 1}}}})/(\epsilon^2 - 1)^{{3/2}})</math>

הדרך למצוא את איברי מטריצת Q היא [[אנליזה נומרית|נומרית]]. הניסיון מראה שמטריצת Q משתנה לאט מאוד, כך שרק מספר ערכים שלה המתאימים לחלקים שונים של שלב המעוף הממונע צריכים להישמר במחשב הטיל. דרך אחת היא להציב ערכים של L בביטוי לעיל ולקבל כך זוגות (E,L) מהם מפתחים את האקצנטריות ואת שאר הקבועים שבתנאי לזמן המעוף, ולאחר מכן מציבים את הקבועים ומוצאים את זמן המעוף. L צריך להיבחר כך 
שיתקיים שהוא מסדר גודל של: <math>R_E \vartriangle\theta = v_\theta t_f = L/(mR_E) t_f, L = m(R_E)^2\vartriangle\theta/t_f </math> (כלומר שהטיל יעבור את המרחק בין שתי הנקודות על כדור הארץ בזמן <math>t_f</math>, בהנחה שהוא נע ב[[עקום גאודזי]] ביניהם). ניתן להציב כפולות רציונליות גדולות מ-1 של הגודל של L שמופיע לעיל (מונה [[מספר רציונלי|המספרים הרציונליים]] נבחר לפי דרגת הקרבה הרצויה של זמן המעוף ל-<math>t_f</math>) ולבחור את הכפולה עבורה מתקבל זמן מעוף קרוב ביותר ל-<math>t_f</math>. לאחר מכן יש לבצע אותו תהליך לגבי מיקום מעט שונה של הטיל, ולמצוא את 9 הנגזרות החלקיות המתאימות. הדרך המעניינת והמאתגרת יותר, מבוססת על פיתוח אלגוריתמים יעילים להערכת הביטוי לזמן המעוף. אלגוריתמים כאלה צריכים למצוא קיצורי דרך אלגבריים בביטויים המתקבלים באמצעות הצבות מתאימות, ולהיעזר בשיטות איטרטיביות להערכת הביטוי לזמן המעוף. פיתוח כזה דורש ידע רב ב[[מתמטיקה טהורה]]. בימים הראשונים של הטילים הבליסטיים, כאשר המחשבים היו [[מחשב אנלוגי|אנלוגיים]], האתגר שעמד בפני המדענים שפיתחו את מערכת ההנחיה Q היה בפיתוח תהליך מתמטי יעיל להערכת הביטויים המתקבלים, כאלה שלא יעמיסו יותר מידי על מחשבי הטיל ויצריכו כח מחשובי רב מידי. 

* '''הערה''' - טיל בליסטי מתחיל את מעופו (מכן השיגור) כשהוא נע אנכית לפני הקרקע. הזווית המדויקת באזימוט אליה יש לפנות לאחר מכן, כמו גם ה-<math>\vartriangle\theta</math> (הזווית המרכזית ב[[מעגל גדול|מעגל הגדול]] המתאים) שמופיע במשוואות לעיל, ניתנים לחישוב מ[[טריגונומטריה ספירית]] לגבי מיקום נקודות השיגור והמטרה (כשנתונים [[קו אורך|קווי האורך]] וה[[קו רוחב|רוחב]] המדויקים שלהם). אם נתונים מיקומיהם בקווי אורך ורוחב של [[כן שיגור]] הטיל הבליסטי והמטרה כ-<math>(\theta_1,\phi_1)</math> ו-<math>(\theta_2,\phi_2)</math> אז כדי לחשב את <math>\vartriangle\theta = \alpha</math> נעזרים בנוסחה הבאה: <math> cos\alpha =cos\phi_1cos\phi_2cos(\theta_2 - \theta_1) + sin\phi_1sin\phi_2</math>. כדי לחשב את הזווית (ה[[אזימוט]]) נעזרים בנוסחה הבאה: <math>cos\beta = \frac {{sin(\phi_2 - \phi_1)}}{{2tan(\alpha/2)}} + sin(\frac{{\alpha}}{{2}})sin(\frac {{\phi_2 - \phi_1}}{{2}})</math>. תוצאה זאת מתקבלת ממציאת ההיטלים של הקשתות במשולש הקשתי המתקבל על מישור קודקודי המשולש ולאחר מכן הפעלת [[משפט פיתגורס#משפט פיתגורס בגאומטריה לא אוקלידית|משפט פיתגורס]] בטריגונומטריה ספירית פעמיים.

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

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

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

*

[[קטגוריה:טכנולוגיה צבאית]]
[[קטגוריה:הנחיית טילים]]