Coordinate Transformations โ€” First Order Systems of ODEs

I just learned a new trick to transform systems of ordinary differential equations in Cartesian coordinates into polar coordinates. It is kind of basic, but leads to interesting physical interpretations, so I will record it here.

As everyone knows, ๐‘ฅ=๐‘Ÿcos(๐œƒ) and ๐‘ฆ=๐‘Ÿsin(๐œƒ). We will use these formula below.

By definition ๐‘Ÿ2=๐‘ฅ2+๐‘ฆ2. Taking the derivative with respect to time ๐‘ก,

2๐‘Ÿ๐‘Ÿห™=2๐‘ฅ๐‘ฅห™+2๐‘ฆ๐‘ฆห™โ‡’๐‘Ÿ๐‘Ÿห™=๐‘ฅ๐‘ฅห™+๐‘ฆ๐‘ฆห™.

A formula for ๐‘Ÿห™. We can derive a similar formula for ๐œƒห™. Differentiate the definitions of ๐‘ฅ and ๐‘ฆ with respect to ๐‘ก:

๐‘ฅห™=๐‘Ÿห™cos(๐œƒ)โˆ’๐‘Ÿ๐œƒห™sin(๐œƒ)๐‘ฆห™=๐‘Ÿห™sin(๐œƒ)+๐‘Ÿ๐œƒห™cos(๐œƒ),

which follows from application of the product rule and chain rule. To isolate ๐œƒห™, we use a technique. We multiply the equation for ๐‘ฅห™ by a term proportional to sin(๐œƒ), and the equation for ๐‘ฆห™ with a term proportional to cos(๐œƒ). We do this so we can apply the well known identity cos2(๐œƒ)+sin2(๐œƒ)=1.

๐‘ฅห™๐‘ฆ=๐‘Ÿห™cos(๐œƒ)๐‘Ÿsin(๐œƒ)โˆ’๐‘Ÿ2๐œƒห™sin2(๐œƒ)๐‘ฆห™๐‘ฅ=๐‘Ÿห™sin(๐œƒ)๐‘Ÿcos(๐œƒ)+๐‘Ÿ2๐œƒห™cos2(๐œƒ),

Now, subtracting the two expressions allows us to isolate ๐œƒห™ with the identity:

๐‘ฆห™๐‘ฅโˆ’๐‘ฅห™๐‘ฆ=๐‘Ÿห™sin(๐œƒ)๐‘Ÿcos(๐œƒ)+๐‘Ÿ2๐œƒห™cos2(๐œƒ)โˆ’(๐‘Ÿห™cos(๐œƒ)๐‘Ÿsin(๐œƒ)โˆ’๐‘Ÿ2๐œƒห™sin2(๐œƒ))=๐‘Ÿ2๐œƒห™cos2(๐œƒ)+๐‘Ÿ2๐œƒห™sin2(๐œƒ)=๐‘Ÿ2๐œƒห™(cos2(๐œƒ)+sin2(๐œƒ))=๐‘Ÿ2๐œƒห™.

Relatively simple combinations of the original expressions in terms of ๐‘ฅ, ๐‘ฆ, ๐‘ฅห™, and ๐‘ฆห™ yield simple expressions in polar coordinates. The term ๐‘Ÿ2๐œƒห™ bears similarity to the ๐‘ง-component of a cross product ๐‘Ÿโƒ—ร—๐‘ฃโƒ—, where ๐‘Ÿโƒ— is the position vector and ๐‘ฃโƒ— is the velocity vector. The term ๐‘Ÿ๐‘Ÿห™ bears similarity to a dot product ๐‘Ÿโƒ—โ‹…๐‘ฃโƒ—. These are fascinating similarities to me.