Understanding Motion in Space Calculations

AI Thread Summary
The discussion focuses on using Meeus' Astronomical Algorithms to transform star coordinates, specifically addressing the challenges of calculating motion in space. After deriving the Cartesian coordinates x', y', z' and their changes, the user seeks guidance on converting these back into spherical coordinates. The explanation outlines the necessary steps, including calculating the new radius r_1 and using inverse trigonometric functions to find the declination and right ascension. The user reports success in achieving accurate declination and is close to perfecting the right ascension. Overall, the thread emphasizes the importance of understanding the conversion process between coordinate systems in astronomical calculations.
Scott S
Messages
22
Reaction score
0
Hello All,
Noobie question. I am using Meeus' Astronomical Algorithims to creat a spreadsheet to transform coordinates. I had no problem with precession, nutation, aberration, etc. But, I don't have as good a math background for Motion in Space.
I can follow down to the final formulas for x', y' z'.
I researched further and found Meeus' source in the Sky Catalogue 2000.0 and that really didn't help me either.
I just don't have the math background.
Can anyone expound on what I do AFTER I solve for x', y' and z'?

Any help greatly appreciated.
Thanks, Scott

P.S. I added a zip of a bmp of the relevant page in the Sky Catalogue.
In any explanation, please keep my genetic deficit in mathematics in mind.
 

Attachments

Last edited:
Astronomy news on Phys.org
The text describes how you from the position of a star given by \alpha_0, \delta_0, r_0 for some time t_0 (like 2000.0) and its change per year given by \Delta \alpha, \Delta \delta, \Delta r can find the position of the star t years later for another time t_1 = t_0 + t, assuming that the star's proper motion can be approximated with a line in three-dimensional space over that time.

Imagine that you have a globe (the surface of which are representing the spherical coordinates \alpha, \delta) and you let a pencil represent the true motion of a star from the start to the end of the pencil, and you place it somewhere near the pole of your globe, then it is not hard to see that it is a complicated matter to describe the spherical coordinates of the end-point of the pencil given only the length of the pencil and its starting point in spherical coordinates. If the pencil was very short relative to the globe and close to equator you may just "add the spherical coordinates", but not when its near the pole.

So, the text describes how you can find the spherical coordinates of a star after t years by first transforming the spherical coordinate \alpha_0, \delta_0, r_0 and its rate of change to the Cartesian coordinates x_0, y_0, z_0 and rate of change \Delta x, \Delta y, \Delta z. From this you can find the new Cartesian position x_1, y_1, z_1 of the star after t years by assuming it follows a simple line (where I have taken the liberty to use subscript 1 instead of a accent).

Now, finally, you want to change this new Cartesian position back into a spherical coordinate by using the first set of equation "in reverse", but now relating the new position, that is, starting with

<br /> \begin{flalign*}<br /> x_1 &amp; = &amp; &amp; r_1 \cos \delta_1 \cos \alpha_1 \\<br /> y_1 &amp; = &amp; &amp; r_1 \cos \delta_1 \sin \alpha_1 \\<br /> z_1 &amp; = &amp; &amp; r_1 \sin \delta_1<br /> \end{flalign*}<br />

you can first find

<br /> r_1 = \sqrt{x_1^2 + y_1^2 + z_1^2}<br />

which you can insert in the equation for z_1 to get

<br /> \delta_1 = \arcsin\left( z_1/r_1 \right)<br />

where you use the inverse sine function. Next you can divide the equations for y_1 and x_1 getting

<br /> \frac{y_1}{x_1} = \frac{r_1 \cos \delta_1 \sin \alpha_1}{r_1 \cos\delta_1 \cos\alpha_1} = \frac{\sin\alpha_1}{\cos\alpha_1} = \tan\alpha_1<br />

which means you can find \alpha_1 by taking the inverse tangent:

<br /> \alpha_1 = \arctan\left( y_1 / x_1 \right)<br />

In a spreadsheet you should have the function ATAN2 available that takes the two arguments above, y_1 and x_1, and returns an angle in radians in the interval -\pi/2 to \pi/2 (or -180 to 180 if you convert the radians to degrees first). If you get a negative angle you can add 2\pi (or 360 degree) to it to get a positive angle.

Also, in case you are not already aware of it, note, that when \alpha is the equatorial rectascension coordinate it is often is given as an hour measure instead of degrees; in this case just multiply the value with 15 to get the angle in degrees.

I hope this helped a bit, otherwise feel free to ask.
 
Thank you Filip.

Excellent explanation and thanks for deriving the formulae.

I managed to get the declination perfect and have the ascension very close.
A little more tinkering and hopefully I should get it.

Thanks again, Scott
 
This thread is dedicated to the beauty and awesomeness of our Universe. If you feel like it, please share video clips and photos (or nice animations) of space and objects in space in this thread. Your posts, clips and photos may by all means include scientific information; that does not make it less beautiful to me (n.b. the posts must of course comply with the PF guidelines, i.e. regarding science, only mainstream science is allowed, fringe/pseudoscience is not allowed). n.b. I start this...
Asteroid, Data - 1.2% risk of an impact on December 22, 2032. The estimated diameter is 55 m and an impact would likely release an energy of 8 megatons of TNT equivalent, although these numbers have a large uncertainty - it could also be 1 or 100 megatons. Currently the object has level 3 on the Torino scale, the second-highest ever (after Apophis) and only the third object to exceed level 1. Most likely it will miss, and if it hits then most likely it'll hit an ocean and be harmless, but...
Back
Top