The Xiddlerian astronomers have identified three non-Keplerian planets that circularly orbit a neighboring star. Planet A is three astronomical units away from its star and completes its orbit in three years. Planet B is four astronomical units away from the star and completes its orbit in four years. Finally, Planet C is five astronomical units away from the star and completes its orbit in five years.
They report their findings to Xiddler’s Grand Minister, along with the auspicious news that all three planets are currently lined up with their star. However, the Grand Minister is far more interested in the three planets than the star and wants to know how long it will be until the planets are next aligned. How many years will it be until the three planets are again collinear (not necessarily including the star)?
Let's denote the positions of planets A, B and C at time t (in years) as A(t)=(3cos2πt3,3sin2πt3),B(t)=(4cosπt2,4sinπt2),C(t)=(5cos2πt5,5sin2πt5), respectively. In order to test for collinearity of A(t), B(t) and C(t) we will want to test whether the area of the triangle with ABC has zero area. Given two vectors →u,→v∈R2 the area of the parallelogram spanned by →u and →v is given by the norm of their cross product ‖→u×→v‖ therefore we can calculate the area of the triangle spanned by A(t), B(t) and C(t) is given by S(t)=12‖(B(t)−A(t))×(C(t)−A(t))‖=12abs(|4cosπt2−3cos2πt34sinπt2−3sin2πt35cos2πt5−3cos2πt35sin2πt5−3sin2πt3|)=12|(4cosπt2−3cos2πt3)(5sin2πt5−3sin2πt3)−(5cos2πt5−3cos2πt3)(4sinπt2−3sin2πt3)|=12|20cosπt2sin2πt5−12cosπt2sin2πt3−15cos2πt3sin2πt5+9cos2πt3sin2πt3−20cos2πt5sinπt2+12cos2πt3sinπt2+15cos2πt5sin2πt3−9cos2πt3sin2πt3|=12|−20sinπt10−12sinπt6+15sin4πt15|
Obviously in the context of area, we would want to keep the absolute value, though since we are trying to find zeros of this trigonometric function, it might be easier to just deal with the signed version ˜S(t)=12(15sin4πt15−20sinπt10−12sinπt6). The signed area is shown in the figure below for the first full period (until t=60 years).
We see from eyeballing it, as all good Xiddlerian astronomers are wont to do, that the first zero of ˜S is between t=5 and t=10 years, with ˜S(5)=−19.4952 and ˜S(10)=11.6913. From here we can resort to whatever friendly neighborhood root solving technique floats our Xiddlerian boats to deduce that the next alignments of Planets A, B and C will occur in ˆt=7.76678 years. For instance, the secant method gets with initial a−1=5,a0=10 and an=an−2˜S(an−1)−an−1˜S(an−2)˜S(an−1)−˜S(an−2) with stopping criterion |˜S(an)|<10−6 solution quoted above within 6 steps. The figure below shows the alignment of the planets with respect to their shared sun at the center of their concentric orbits.
No comments:
Post a Comment