Rectifique una onda periódica 

Rectifique una 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]=

Grafique la onda rectificada junto con la 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]=

Rectifique una onda elíptica. En este caso, los tiempos de eventos se obtienen usando un 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]=