Loading [MathJax]/jax/element/mml/optable/LetterlikeSymbols.js
Skip to main content
LibreTexts - Ukrayinska

13.15: Обчислення елементів

Тепер ми можемо відразу обчислити semi latus прямої кишки з Рівняння 13.12.6a (нагадуючи2f3=v2v1, що, так що все,l крім Рівняння, вже відомо.) Насправді у нас є три можливості для обчислення semi latus прямої кишки за допомогою кожного з Рівняння 13.12.6a, b, c, і це служить перевіркою на арифметику. Для нашого числового прикладу я отримую

l=2.61779

однаково (принаймні до одинадцяти значущих цифр) для кожної з трьох перестановок.

Тепер, посилаючись на Рівняння 2.3.37, ми пам'ятаємо, що полярне рівняння до еліпса є

r=l1+ecosv.

Тому ми маємо, для першого і третього спостережень,

ecosv1=l/r11

і, визнаючи цеv3=v1+2f2,

ecos(v1+2f2)=l/r31.

Ми спостерігаємо, що в рівняннях\ ref {13.15.2} і\ ref {13.15.3} єдині величини, які ми ще не знаємо, єv1 іe - тому ми збираємося знайти наш перший орбітальний елемент, ексцентриситет!

Підказка для розв'язання рівнянь\ ref {13.15.2} і 3: Розгорнутиcos(v1+2f2). Візьмітьesinv в ліву сторону, і Equation\ ref {13.15.3} стане

esinv1=(l/r11).cos2f2(l/r31)sin2f2.

Після цього легко вирішувати Рівняння\ ref {13.15.2} і\ ref {13.15.4} дляe і дляv1. Інші справжні аномалії задаютьсяv2=v1+2f3 іv3=v1+2f2. Перевірка на арифметику може (і повинна) виконуватися шляхом проведення однакових обчислень для першого і другого спостережень і для другого і третього спостережень. За всі три я отримав

e=0.23875_

У нас є перший орбітальний елемент!

(MPCЗначення ексцентриситету для цієї епохи є0.22994 - але це ґрунтується на всіх наявних спостереженнях, і ми не можемо очікувати отриманняMPC значення лише з трьох гіпотетичних «спостережень».)

Справжні аномалії за часів трьох спостережень

v1=191.99814v2=192.68221v3=194.05377

Після цього напівпровідна вісь легко з Рівняння 2.3.10l=a(1e2), для прямої кишки прямої кишки еліпса. знаходимо

a=2.77602 au_

Період в сидеричних роках даєтьсяP2=a3, і тому є4.62524 сидеричними роками. Це не один із шести незалежних елементів, оскільки він завжди пов'язаний із напівосновною віссю третім законом Кеплера, тому він не заслуговує на додаткову гідність підкреслення. Однак, безумовно, варто перетворити його на середні сонячні дні, множивши на365.25636. Ми знаходимо тіP=1689.39944 дні.

Наступним елементом для поступки буде час проходження перигелію. Ми знаходимо ексцентричні аномалії для кожного з трьох спостережень з будь-якого з Рівнянь 2.3.16, 17a, 17b або 17c. Наприклад:

cosE=e+cosv1+ecosv.

Тоді час проходження перигелію прийде з Рівнянь 9.6.4 і 9.6.5:

T=tP2π(EesinE)+nP.

Зn=1 я роблю цеT=t1+756d.1319_

Наступним кроком є обчисленняP s іQ s, які визначені в Рівнянні 10.9.40. Вони являють собою косинуси напряму, що стосуються геліоцентричної основи площини орбіти, встановленого до геліоцентричного екваторіального базового набору.

Вправа. Застосуйте рівняння 10.9.50 до першого та третього спостережень, щоб показати, що

Px=ξ1r3sinv3ξ3r1sinv1r1r3sin2f2

іQx=ξ3r1cosvξ1r3cosv3r1r2sin2f2

З Рівняння 10.9.51 і 52 знайти подібні рівняння дляPy,Qy,Pz,Qz.

Чисельну роботу можна і потрібно перевіряти, обчислюючи ці напрямні косинуси також з першого і другого, а з другого і третього спостережень. Перевірте також, щоP2x+P2y+P2z=Q2x+Q2y+Q2z=1. Я отримую

Px=0.48044Py=+0.86568Pz=0.14059

Qx=0.87392Qy=0.45907Qz=+0.15978

(Пам'ятайте, що мій комп'ютер несе всі значні цифри з подвійною точністю, хоча я роздруковую тут лише обмежену кількість значущих цифр. Ви не отримаєте точно мої цифри, якщо ви теж не несете всі значущі цифри і не передчасно округляєте.)

Напрямки косинуси пов'язані з кутами Ейлера, звичайно, рівняннями 10.9.41- 46 (як ви могли б забути?!). Всі (!) Ви повинні зробити, то, щоб вирішити ці шість рівнянь для кутів Ейлера. (Вам потрібно шість рівнянь, щоб усунути неоднозначність квадрантів з кутів. Запам'ятайтеATAN2 функцію на вашому комп'ютері - це величезна допомога з квадрантами.)

Вправа. Показати, що (або перевірити в будь-якому випадку), що:

sinωsini=PzcosεPysinε

іcosωsini=QzcosεQysinε.

Тепер ви можете вирішити це для аргументу перигеліюω. Ще не намагайтеся вирішити це за схильність. (Чому б і ні?!) Використовуючиε=23.438 960 для конусності екліптики дати (обчислюється зі сторінкиB18 Астрономічного альманаху 2002 року), отримую

ω=304.81849_

Вправа. Показати, що (або перевірити в будь-якому випадку), що:

\sin Ω = (P_y \cos ω - Q_y \sin ω ) \sec ε \label{13.15.11} \tag{13.15.11}

і\cos Ω = P_x \cos ω - Q_x \sin ω. \label{13.15.12} \tag{13.15.12}

З них я знаходжу:

\underline{Ω = 172^\circ .64776}

Ще один, щоб піти!

Вправа. Показати, що (або перевірити в будь-якому випадку), що:

\cos i = -(P_x \sin ω + Q_x \cos ω) \csc Ω . \label{13.15.13} \tag{13.15.13}

Тепер ви можете вирішити це за допомогою Equation\ ref {13.15.9} або\ ref {13.15.10} (або обидва, як перевірка арифметики) для нахилу. Я отримую

\underline{i = 35^\circ .20872}

Ось вони, всі разом:

\ почати {масив} {c c}
a = 2.77602\ текст {АС} & я = 35^\ коло .20872\
е = 0,23875 & Ω = 172^\ коло .64776\\
T = t_1 + 756^\ текст {d} .1319 & ω = 304^\ цирк .81849\\
\ кінець {масив}

Чи допустили ми якісь помилки? Ну, імовірно, після того, як ви прочитали главу 10, ви написали програму для створення ефемериди. Отже, тепер використовуйте ці елементи, щоб побачити, чи будуть вони відтворювати оригінальні спостереження! До речі, для побудови ефемериди немає необхідності фактично використовувати елементи - замість них можна використовуватиP s іQ s.