Taking the equation for simple harmonic motion:
We can introduce the derivative operator:
This gives us:
This looks like a matrix, so let's see what happens in 2D:
Solving gives us either or .
If we assume that then this implies which cannot be solved in the reals, so
it must be that :
We still have . This only has 2 degrees of freedom (you can only pick 2 of and the
third will be fixed) so lets try rewriting it:
Let :
We can see that gives us the rotation matrix that we expect to see in the perfectly circular case
since we then have that the velocity is 90 degress out of phase with the position, which is itself 90 degress out of
phase with the acceleration:
This can be simplified using the identity :
And here we have the general solution for the instantaneous velocity and acceleration of a point undergoing simple
harmonic motion in 2D. Note that can be folded into since but
.
The demo bit
Since we have a closed form solution for the general case, here's an interactive demo. are as
above, and define the offset of the orbiting body (the black line). defines how many seconds to
simulate into the future (the blue line). The red and green lines show the instantaneous velocity and acceleration
respectively. You can click to place the point at an arbitrary position.