Questions tagged [astronomy]

Questions on the astronomical applications of Mathematica.

Filter by
Sorted by
Tagged with
17
votes
2answers
5k views

What is the closest star to Earth other than the sun?

How to find the Sun's nearest neighbor star using Mathematica? I tried but it didn't work
3
votes
1answer
141 views

How to plot the hour of the sunsets

i have this code: ...
5
votes
2answers
308 views

How to find the Andromeda Galaxy

I would like to find the Andromeda Galaxy ($M31$) with the unaided eye. One way is to use the following image and graphic but I can’t find it just by looking in a dark sky. It is hard to detect ...
12
votes
1answer
695 views

Can I create a radar chart of sunrise/sunset times for a year

Can I create a radar chart (or just a simple curve) of sunrise/sunset times, for a given location thougout a year? It would be something like this, but with a smooth curve: (source: hassanhameed.com)...
8
votes
2answers
645 views

Extracting and displaying data from GW150914 hdf5 files

Has anyone looked at the LIGO data using Mathematica or know how to do this? I'm not familiar with the hdf5 format but I've been able to import the data, consisting primarily of 131,072 floating ...
10
votes
1answer
340 views

3D orbits of moons around their respective planets

I have seen a lot of solar system models (both animated and not), but I am trying to replicate the same thing with moons. Either just one planet and its lunar system or the our whole solar system with ...
3
votes
1answer
101 views

Lack of precision in solar declination leads to odd Fourier result

Here's the sun's declination for a 10 year period (roughly 2000-2010): ...
3
votes
1answer
136 views

Convert star name into entity value

I have a list of star names: stan = Select[#, Last[#] < Quantity[153.5, "LightYears"] &] &@ StarData[StarData[], {"Name", "DistanceFromEarth"}] ...
6
votes
2answers
598 views

Series solution of the Lane-Emden equation

I have the following problem: I'd like to show how the Lane-Emden equation would look like if someone would solve it as a series. Here's the code: ...
4
votes
1answer
135 views

Impact Crater Sim. with added criteria

I have been trying with out success to edit this impact crater simulation script into producing a variable which calculates crater density. In addition I have been trying to figure out how to ...
6
votes
2answers
547 views

Solving the path of Earth around Sun

This maybe isn't a universally helpful question. Maybe a little of a code dump. But here goes. I'm trying to solve the path Earth moves around Sun from Earths mass, Suns mass, Earths initial velocity ...
5
votes
1answer
175 views

Confused by result of MoonPosition when running on Linux [closed]

I've just switched from Windows to Linux and am running Mathematica 10.1 on Ubuntu 14.04. MoonPosition[] Gives me: ...
6
votes
0answers
106 views

Heliocentric coordinates in PlanetData

I tried to use planets' HeliocentricLongitude and I wonder if anyone has any insight on what coordinate system they use. I took vernal equinox times from http://ns1763.ca/equinox/eqindex.html and ...
11
votes
1answer
241 views

How can I get curated data to be expressed consistently in SI units?

AstronomicalData objects are consistently in SI units: ...
10
votes
1answer
394 views

Finding planetary conjunctions with Mathematica (project-level) [closed]

Summary: What's an efficient way of finding planetary conjuctions with Mathematica using the SPICE kernels, and without using AstronomicalData (which I believe is less efficient). NASA publishes ...
3
votes
1answer
193 views

Computing planet conjunctions with 2D circular orbits still hard?

All methods I've seen for computing planetary/stellar conjuctions (when two planets or a planet and a star appear close together in the sky because their angular separation, as measured from Earth, ...
14
votes
1answer
664 views

GeoProjection for astronomical data - wrong ticks

Here is an excellent answer to what I wanted to do, and it works like a charm. Unfortunately, when I want to change it a bit, it fails. For setting the stage, here is the datafile, and this is my ...
2
votes
0answers
61 views

PlanetData “VelocityAroundSun” giving inconsistent answers

My goal is to create a ListPlot showing the Earth's orbital speed over the course of a year, with one datapoint per day. I'm seeing problems long before I even get to the point of plotting the data! ...
5
votes
1answer
590 views

Using ParametricNDSolve and NMinimize to find accurate interplanetary trajectory arrival periapse radius

I have an interplanetary trajectory simulation that calculates the required velocity (using a Lambert solver) and escape angle for a spacecraft to travel from Earth to Mars. This model does not use a ...
8
votes
1answer
399 views

ParametricNDSolve[ ] and NMinimize[ ] in Combination

I have the following code that is used to simulate a spacecraft orbiting Mars, and uses ParametricNDSolve to vary the spacecraft's initial conditions using a parameter value ...
6
votes
0answers
134 views

Is there a problem with some of the curated Star Data?

I'm not really sure if I'm doing something wrong, or if there's a problem with some of the Mathematica curated data on stars. My problem is this: I'm trying to find the location (latitude and ...
2
votes
2answers
166 views

Finding Minimium Periapse Distance Between Spacecraft and Planet

I have the following simple model that simulates a spacecraft orbiting Mars. What I'm hoping to do is be able to find the periapse radius of the spacecraft as I did below, except instead of using ...
5
votes
1answer
218 views

Inaccurate Heliocentric coordinates of PlanetaryMoonData

I have noticed that the heliocentric coordinates obtained with PlanetaryMoonData are very inaccurate. Fore example, consider the Galilean moons of Jupiter: ...
3
votes
1answer
134 views

How to find µ (Gravitational constant mass product) with PlanetData?

I would like to use PlanetData to get µ, the product of the Gravitational constant G and a ...
8
votes
2answers
518 views

Generate a list of symbols for the planets

I would like to generate a list of the symbols for the planets to include in a nice table with other data. How can I programmatically generate (eg. using a pure function) the symbols from a list of ...
4
votes
1answer
93 views

Accessing Coordinates of Astronomical Data, e.g. Declination of all Galaxies

I am trying to access the coordinate information of all galaxies, stars, etc, available to the new astronomical data incorparated to Mathematica. However, the suggested (documented) approach fails. ...
21
votes
1answer
1k views

Right Ascension and Declination in astronomical functions

In the process of answering this question, I was forced to confront the various astronomical coordinate systems used by Mathematica. Background In astronomy, positions of celestial objects (stars, ...
7
votes
0answers
161 views

MoonPosition and TimeZone support

Bug introduced in 10.0 and fixed in 10.2 MoonPosition is new in 10.0 With this question, I noticed some some weird results of ...
13
votes
1answer
1k views

Geocentric, Earth-fixed, coordinates of Sun, Moon (and others objects eventually)

I know I can get "HelioCoordinates" of the Earth and the Moon with ...
7
votes
1answer
273 views

How to get star coordinates for a given time

Fixed in 10.3.1 In the documentation of StarData, there is an example that shows how to get the altitude of Sirius at a given time (in the horizontal coordinate ...
15
votes
2answers
970 views

Align stars in multiple images

I hope that you find this question interesting and that this will help others with similar problems as this configuration is not so particular. My problem I have multiple images of a particular ...
19
votes
2answers
1k views

Using map projections with astronomical data

I noticed that all important "Geoprojections" are available in projections for a spherical reference models: GeoProjectionData function; 1 - I am trying using the sinusoidal projection for ...
4
votes
0answers
392 views

HEALPix grid on a sphere

Are there any Mathematica tools for using the HEALPix ( Hierarchical Equal Area isoLatitude Pixelation) grid on a sphere? It is commonly used in astronomy, but I have not found any Mathematica ...
2
votes
1answer
968 views

Astronomy: Transform coordinates from horizon to equatorial

I'm using Mathematica 10. Giving a longitude & latitude on Earth, giving a local time, how can I transform Horizon Coordinte (az , alt) to Equatorial Coordinate (RA , DEC). Use function so I can ...
1
vote
1answer
170 views

Plotting Earth's position with units

I'm trying to plot the Earth position with the Sun in the origin. I have two lists like this: ...
2
votes
0answers
202 views

Creating an orrery

I am interested in creating an orrery, ideally with SystemModeler. I have been experimenting with the MultiBody library since that allows you to do animations. The most relevant example seems to be ...
5
votes
1answer
4k views

Combining Gravity Turn and Orbit Models

I have a mathematical model for the motion of an orbiting spacecraft about Earth: ...
20
votes
1answer
517 views

Segmentation of sunspots

I recently went for a walk with my camera and snapped a few pictures of the sun. The pictures didn't turn out too well, but this was to be expected without a telescope. I thought it might be a nice ...
5
votes
1answer
213 views

Obsolete AstronomicalData - how to convert to PlanetData

As of Mathematica version 10, AstronomicalData is obsolete, deprecated in favor of PlanetData, ...
2
votes
2answers
296 views

Problem evaluating a complicated integral

Apologies for stating my problem poorly in the first instance, thank you for the help in narrowing down the issue. Problem For those interested, it's for the exact unbiased inverse of the ...
7
votes
3answers
1k views

RK4 Gravity Simulator

I have the following RK4 solver which splits the two 2nd order ODEs, used to calculate x and y positions under the influence of a gravitating body where $$x''(t)=\frac{G m x(t)}{(x(t)^2+y(t)^2)^{3/2}}$...
7
votes
2answers
970 views

Combining three R, G, B matrices into a single matrix of 3-tuples?

Working with Mathematica 8 on OSX here. I have three $1024 \times 1024$ matrices representing the counts from consecutive exposures onto a CCD with $R$, $G$, and $B$ filters on it, and am looking to ...
10
votes
1answer
614 views

How can I get this Graphics3D to rotate around a particular point?

I made this Graphics3D of the Sun with the plotted orbits of a few astronomical objects. ...
1
vote
1answer
748 views

Combining solar system planets orbits and positions in Graphics3D and Show

I am building an animated model of the solar system, where I would like to display the OrbitPaths of the planets along with their ...
0
votes
1answer
118 views

AstronomicalData not including units?

I like using mathematica to do my physics homework, and AstronomicalData is useful for that, but it doesn't seem to work properly with units. [the question is, at ...
6
votes
1answer
216 views

Finding minimum fly-by radius between Mars and spacecraft from interpolating function

I've written an interplanetary trajectory solver/plotter that plots the path taken by a spacecraft on an Earth-Mars mission, but have run into a little trouble when the spacecraft actually reaches ...
2
votes
1answer
240 views

Using DeleteCases with multiple AstronomicalData properties

I am attempting to use AstronomicalData and DeleteCases together to remove invalid cases of two properties while listing the ...
1
vote
2answers
167 views

Computing lists using AstronomicalData

Edit: Initially in my question I asked how I could search the database of AstronomicalData["Exoplanets"] and return 2 filtered values that only contained valued ...
5
votes
1answer
1k views

p-iteration algorithm to solve Lambert's problem for interplanetary trajectories

I've been reading up on using the p-iteration method to solve Lambert's problem for choosing the correct interplanetary trajectory between two planets given the ...
10
votes
2answers
530 views

Impact Crater Simulation

I have been trying to model crater formation on a given planetary surface - $500\,\textrm{km}^2$. The locations of impacts are random, however, if an impact is within $30\,\textrm{km}$ of another, the ...