Retifique uma onda periódica 

Retifique uma onda senoidal.

In[1]:=

X sol = DSolve[{x''[t] + x[t] == 0, x[0] == 0, x'[0] == 1, WhenEvent[x[t] == 0, x'[t] -> -x'[t]]}, x, {t, 0, 4 \[Pi]}]

Out[1]=

Exiba o gráfico da onda retificada junto com a onda senoidal original.

In[2]:=

X {Plot[Sin[t], {t, 0, 4 \[Pi]}, Filling -> Axis, FillingStyle -> Orange], Plot[(x[t] /. sol), {t, 0, 4 \[Pi]}, Filling -> Axis, FillingStyle -> Green]}

Out[2]=

Retifique uma onda elíptica. Neste caso, os tempos dos eventos são obtidos através de um método numérico.

In[3]:=

X sol = DSolve[{x''[t] + Sin[x[t]] == 0, x[0] == 0., x'[0] == 1, WhenEvent[x[t] == 0, x'[t] -> -x'[t]]}, x[t] , {t, 0, 10}]

Out[3]=

In[4]:=

X {Plot[2 JacobiAmplitude[t/2, 4], {t, 0, 10.5}, Filling -> Axis, FillingStyle -> Orange], Plot[(x[t] /. sol), {t, 0, 10.5}, Filling -> Axis, FillingStyle -> Green]}

Out[4]=