המשוואה הכללית לשינוי בתדירות עקב העברת דופלר היא:
כאשר c היא מהירות האור, f O הוא תדירות הפעולה, ו- Δv הוא המהירות היחסית.
בעת החלת נוסחה זו על תצפית קרקעית של לווין, Δv מתואר טוב יותר כמהירות טווח (קצב של שינוי מרחק מהלוויין למתבונן). מכיוון שלווייני LEO רבים נמצאים במסלול אליפטי, המתמטיקה להשגת Δv כפונקציה של זמן או מיקום אינה נפתרת באמצעות חישובים זוויתיים פשוטים ולכן הופכת להיות מייגעת מאוד לביצוע גם כאשר מיישמים מתמטיקה מטריצתית לחלק מהחישובים.
ברצוני להציע, עם זאת, כי אין ערך רב לשחזר את כל החישובים למעט עניין אקדמי. במקום זאת, אני מעדיף להשתמש בספריית תוכנה מתוכננת היטב בשם PyEphem שעושה כמעט כל סוג של חישוב ארעי שניתן להזדקק לו (אני אפילו משתמש בספרייה זו במערכת האוטומציה הביתית שלי כדי לחשב את זמני ההפעלה / כיבוי של האורות החיצוניים שלי כדי להשיג עונתיות אוטומטית. התאמה). זוהי ספריית פייתון חינמית המופעלת בקלות כמעט בכל מערכת הפעלה ופלטפורמה, כולל Raspberry Pi בכל מקום.
יישום הספרייה לווייני רדיו חובבים הוא די פשוט. ראשית, רכש את ה- TLE (שני אלמנטים בקו) עבור הלוויין הרצוי. אלה זמינים מ אתר TLE Info. קבוצה טיפוסית של נתוני TLE מכילה את כל המשתנים הארעיים הדרושים ונראית כך:
OSCAR 71 7530U 74089B 17170.24378275 -.00000031 + 00000-0 + 84707-4 0 99922 7530 101.6303 138.8875 0011838 320.4872 153.4129 12.53627377948883
מידע זה משמש בספריית PyEmphem באופן הבא. ראשית קבע את מיקום הצופה על פני כדור הארץ (אשתמש כדוגמה בשיקגו, אילינוי, ארה"ב):
my_loc = ephem.Observer () my_loc.lon = '87 .6298'my_loc.lat = ' 41.8781'my_loc.elevation = 181
ואז החל את נתוני ה- TLE ליצירת אובייקט גוף של my_sat:
my_sat = ephem.readtle (name, line1, line2);
The כעת הספרייה יכולה לחשב מגוון נתונים הקשורים למעבר הבא של my_sat ב- my_loc:
info = my_loc.next_pass (my_sat)
אז למשל , אנו יכולים להדפיס את זמני העלייה וההגדרה, הגובה המרבי וזמן הגובה המרבי עבור המעבר הבא:
הדפסה ("AOS:% s LOS:% s גובה מרבי:% זמן גובה מקסימאלי:% s "% (מידע [0], מידע [4], מידע [3], מידע [2]))
שים לב שההזזה המקסימלית של הדופלר תתרחש ב- AOS ו- LOS וכי זמן הגובה המרבי הוא כאשר הסטת הדופלר תהיה אפס מכיוון שמהירות הטווח (Δv) תרד לאפס.
אנו יכולים גם לברר באשר למהירות הטווח בכל עת במהלך המעבר. לדוגמא, בתחילת המעבר:
my_sat.compute (info [0]) הדפס ("מהירות טווח:% s"% (my_sat.range_velocity))
על ידי יישום טכניקה זו כדי להשיג את מהירות הטווח לאורך המעבר, אנו יכולים לחשב בקלות את Δf של מעבר מסוים כפי שנשקף ממיקום ספציפי על פני כדור הארץ. וכל זה רק לכמה עשרות שורות קוד!
תוכל ללמוד עוד על ספריית PyEphem הרחבה כאן.
אם ברצונך לחקור יותר מהמתמטיקה הגולמית שמאחורי חישובי הלוויין של LEO, הקישורים הבאים עשויים להועיל:
חישוב מסלולים ותיקון דופלר
תיקון דופלר מסתגל
יסודות מסלול הלוויין