Mathematica code:

P[t_] := (1 + 0.9 Cos[8 t]) (1 + 0.1 Cos[24 t]) (0.9 + 0.05 Cos[200 t]) (1 + Sin[t])

T[t_] := 
 If[0 < t < 17, .01 t, 
  If[16 < t < 25, .2 + .1 (t - 16)/2, 
   If[24 < t < 35, 0 + (t - 24)/2, 
    If[34 < t < 45, -3 + 20*(t - 34)/2, 0]]]]

ListAnimate[
 Table[
  Show[

   PolarPlot[
    Evaluate[
     Table[
      (T[Mod[s, 45]]*1/.85^n)*P[t],
      {n, 0, 20, 1}],
    {t, -Pi, Pi}], 
   PlotRange -> {{-7.5, 7.5}, {-4.3, 15.7}}, 
   Background -> Black, Axes -> False, ImageSize -> 500],

   PolarPlot[
    Evaluate[
     Table[
      (T[Mod[s + 5, 45]]*1/.85^n)*P[t],
      {n, 0, 20, 1}],
    {t, -Pi, Pi}], 
  PlotRange -> {{-7.5, 7.5}, {-4.3, 15.7}}, 
  Background -> Black]],

{s, 1, 45, 2}]]