Như tiêu đề cho thấy tôi đang cần một bộ dữ liệu của hệ mặt trời. Tương tự như http://bima.astro.umd.edu/nemo/archive/#iau25 , cụ thể là dữ liệu về cách thức và dữ liệu của sữa mẹ.
Tôi đang tạo một mô phỏng cơ thể n cho trường học và dường như không thể tìm thấy bất kỳ điều kiện hạt ban đầu nào mà tôi có thể sử dụng để mô phỏng hệ mặt trời của chúng ta trong phần mềm tôi đang phát triển. Tôi cần vị trí ban đầu, vận tốc và khối lượng.
Bất cứ ý tưởng nơi tôi có thể tìm thấy điều này?
Cảm ơn bạn.
Câu trả lời:
Trình tạo Ephemeris của Horizons có thể cung cấp cho bạn các vectơ vị trí và vận tốc của một hành tinh tại một thời điểm xác định.
Đây là một trong các tùy chọn có thể:
Nhấp vào Tạo Tạo thiên thạch trên trang đó sẽ cung cấp cho bạn các vectơ vị trí và vận tốc:
Trên các vectơ vị trí và vận tốc là ngày Julian cũng như ngày thông thường hơn.
Wikipedia có thể cung cấp cho quần chúng mặt trời và các hành tinh.
Tôi đã có thể có được các vectơ quỹ đạo của Cartesian cho tất cả các cơ quan chính từ HORIZON ở kỷ nguyên J2000. Tôi có thể mở rộng phạm vi bảo hiểm về phía trước thời gian thông qua. Thật dễ dàng để có được quá tải dữ liệu làm điều này. Mô phỏng của tôi được mô hình hóa bằng cách sử dụng Định luật hấp dẫn và chuyển động. Điều này mang lại kết quả gần đáng ngạc nhiên với những người được công bố. Vận hành hệ mặt trời ngược (bằng cách đảo ngược các vectơ vận tốc) đã cho tôi các vectơ ban đầu trở về năm 1900. Đây là tất cả những gì tôi cần và kết quả đã đủ gần với mục đích của tôi. Tôi vẫn có các tệp CSV.
Tôi cũng đã có tất cả các loại vấn đề với giao diện chân trời. Chẳng hạn, việc thay đổi ngày không ảnh hưởng đến giá trị của vectơ. tức là: tất cả các ngày bắt đầu được chỉ định có cùng giá trị. Gần đây, tôi đã không thể lặp lại kỳ tích này. Rõ ràng có một số vấn đề nghiêm trọng với giao diện này, đặc biệt là gần đây.
Tôi biết dữ liệu tôi nhận được là chính xác vì nó tương quan, hoàn hảo, với các sự kiện được công bố, ví dụ: quá cảnh gần đây của Sao Thủy.
Tôi cũng vẫn đang tìm kiếm loại dữ liệu này.
Ngôn ngữ nào bạn đang viết mô phỏng này? Nó là 2D hay 3D? Bạn chỉ cần vị trí và vận tốc cho các hành tinh hệ mặt trời của chúng tôi?
Tôi đã thực hiện điều này chính xác (mô phỏng hệ mặt trời ở Fortran) và tôi không cần vị trí ban đầu chính xác, tất cả những gì tôi cần là bán kính ban đầu (tính bằng AU từ Mặt trời / trung tâm khối lượng) và vận tốc ban đầu. Sử dụng một trình tạo số ngẫu nhiên để phân phối các hành tinh tại các vị trí ngẫu nhiên dọc theo quỹ đạo của chúng. Ở Fortran, nó trông giống như:
CALL RANDOM_NUMBER(randNum)
degrees = 2*3.141592653
theta(1:15) = degrees*randNum(1:15)
Và ở đó tôi có một mảng gồm 15 vị trí xuyên tâm ngẫu nhiên. Bạn có thể có được vận tốc ban đầu của các hành tinh trong hệ mặt trời của chúng ta từ bất kỳ nguồn tài nguyên có uy tín nào.
Có lẽ bạn đã chuyển từ lâu, nhưng, chỉ để tham khảo, các điều kiện ban đầu mà HORIZONS sử dụng được đề cập ("tiêu đề.431_572") trong ftp://ssd.jpl.nasa.gov/pub/eph/planets/Linux/ README.txt nhưng nơi duy nhất tôi có thể tìm thấy chúng trong "biểu mẫu bảng" là trong kho git của riêng tôi:
https://github.com/barrycarter/bcapps/blob/master/ASTRO/header.431_572
Các giá trị được giải thích trong http://ilrs.gsfc.nasa.gov/docs/2014/196C.pdf bắt đầu từ trang 39, "VI. Điều kiện và hằng số ban đầu", đặc biệt là trong các bảng bắt đầu với Bảng 4 trên trang 47 và kết thúc với Bảng 13 trên trang 74.
Tôi đã viết các tập lệnh để thiết lập các điều kiện ban đầu và giải quyết số lượng các phương trình vi phân bằng Mathicala, vì vậy những điều sau đây có thể hữu ích:
https://github.com/barrycarter/bcapps/blob/master/ASTRO/README
https://github.com/barrycarter/bcapps/blob/master/ASTRO/bc-header-values.pl
https://github.com/barrycarter/bcapps/blob/master/ASTRO/bc-integrate.m
Khi tôi giảm đủ kích thước bước (kích thước bước mặc định của Mathistica quá lớn), kết quả của tôi rất khớp với kết quả của HORIZONS:
https://github.com/barrycarter/bcapps/blob/master/ASTRO/bc-integrate-compare.m