Masalah Pemusatan yang Jelas Saat Menggunakan Sistem Perifokal untuk Menggambar Elips Orbital
Saya mencoba membuat elips orbital dalam perangkat lunak. Diberikan satu set elemen orbital, saya menggunakan sistem Perifocal untuk menentukan satu set titik di sepanjang kurva elips. Saya menggunakan persamaan:
r = r cos ( v ) P + r sin ( v ) Q ,
di mana v adalah sudutnya, P dan Q adalah vektor perifokal yang berkoresponden dengan I dan J , dan r adalah persamaan kutub dari kerucut:
r = p / 1 + e cos ( v ),
di mana p adalah rektum semi-latus, e adalah eksentrisitas, dan v , dalam kedua persamaan, adalah sudutnya.
Saya menentukan vektor P dan Q menggunakan persamaan berikut:
Pi = cos (lan) cos (aop) - sin (lan) cos (inc) sin (aop)
Pj = sin (lan) cos (aop) + cos (lan) cos (inc) sin (aop)
Pk = sin (inc) cos (aop)
Qi = - cos (lan) sin (aop) - sin (lan) cos (inc) cos (aop)
Qj = - sin (lan) sin (aop) + cos (lan) cos (inc) cos (aop)
Qk = sin (inc) cos (aop)
dimana lan adalah Bujur Node Naik, aop adalah Argumen Periapsis, dan inc adalah kemiringan.
Proses ini menghasilkan elips yang benar! Namun, elips tampaknya tidak berada di tengah dengan benar pada fokusnya; fokusnya tampak terlalu dekat dengan pusat elips. Berikut ini contohnya:
Anda dapat melihat bahwa dua elips yang sangat eksentrik ini memiliki pusat massa yang sama di pusatnya - bukan di fokusnya masing-masing.
Apakah ada beberapa langkah dalam proses saya yang saya lewatkan? Saya hanya bisa menerapkan offset ke elips, tetapi itu akan menjadi peretasan, saya ingin menyelesaikan masalah.
Bantuan apa pun dihargai! Terima kasih!
EDIT PERTANYAAN DAN KLARIFIKASI:
Gambar ini diambil dari titik pandang yang hampir normal ke bidang orbit; hampir tidak ada perspektif di sini.
Kedua, pengguna yang sangat membantu telah menunjukkan masalah dalam persamaan untuk Pk. Ekspresi yang benar adalah:
Pk = sin (inc) sin (aop)
Sayangnya, ini adalah benar dalam kode saya - maka sementara koreksi itu bagus, itu bukan sumber masalah saya.
Jawaban
Menurut definisi Bingkai Perifokal di wikipedia , persamaan Pk membutuhkan penyesuaian. Saya tidak yakin apakah itu sumber kesalahan yang Anda amati.
Dari pada:
Pk = sin (inc) cos (aop)
Wikipedia mengatakan:
Pk = sin (inc) sin (aop)
Saya memecahkan masalah tersebut.
Semua yang saya posting di OP benar; Saya membuat kesalahan tanda yang sederhana dan bodoh saat mentranskripsikan ekspresi menjadi kode.
Perlu juga dicatat bahwa dalam OP, ekspresi Pk salah. Ekspresi yang benar adalah
Pk = sin (inc) sin (aop)
Kalau-kalau ada yang mencoba melakukan ini sendiri!
Terima kasih atas idenya semuanya!
Orbit yang Benar:
