I haven't been posting here recently, since I've been busy with research, but the stuff I'm doing at the moment is a bit slow-moving, so I thought I'd stop in and give an update. In February, I mentioned work I was doing on approximations to a piece of solar system dynamics code for LIGO. I'm still working on it, and I wanted to talk about an interesting issue I just resolved.
The coordinates we use to specify a direction in the sky are called equatorial coordinates, given in right ascension and declination. These correspond to taking the longitude and latitude at the spring equinox, and projecting them into space. Since the coordinates are referenced to a specific position of the Earth, as the Earth moves the stars stay fixed in the coordinates, e.g. Procyon is always at 114.8° RA, 5.2° dec. However, since this system assumes a fixed Earth, the planets and even the Sun will move around. That means that as data is being collected from a source, the sun could cross the path and introduce the time delays I mentioned in the previous post.
Over a short time scale, say a week, the dominant frequency will be daily, due to Earth's rotation, but what kind of day? The ones we're used to are called solar days, which correspond to the time it takes for the sun to return to the same position in the sky. However, since the Earth is also moving around the sun, this is a bit longer than the time for the stars to return to their previous positions. That period is called a sidereal day, and is about 4 minutes shorter than a solar day.
I tried fitting the parameters I need using these two frequencies, but comparing the fitted values to the data still showed some structure:
When two sinusoidal waves with slightly different frequencies are added, they result in a wave with a frequency related to the sum and the difference of the two original frequencies: