Continued

II. (a) Use Maple to solve the differential equation describing a damped harmonic oscillator with parameters the same as in I(a) except now there is a damping force characterized by -bv(t) where the damping constant, b = b/2m.

Remember in part Ia there was a mass, m, attached to a spring of force constant, k. The initial position of the mass was at x0 and its initial velocity v0 and the angular frequency was w0 = (k/m)^1/2.

Hint:You will solve the equation diff(x(t),t,t) + 2*beta*diff(x(t),t) + omega0^2*x(t) = 0 with x(0)=x0,D(x)(0)=v0 (like last week's lab);

Assign your solution

(b) Let x0 = 1, v0 = 0, and w₀ = 1. Plot the position x(t) for t from 0 to 4p for the case of b = 0.25 (underdamping), b = 1 (critical damping), b = 4 (overdamping).

Hints: In the case of beta = 1, you will need to plot the limit as beta goes to 1, you may as well do this for all of them, eg plot(limit(x(t),beta = 0.25),t = 0..4*Pi);.

You probably should restart before part c

(c) For the underdamping case of part (b), obtain a simple expression for x(t) in terms of cosines, sines and e^-bt by defining w₁² = w₀² - b². Obtain the simple expression for the case when omega1 = 0.9375, beta = 0.25, x0=1.0, v0=0.0

Hints Resolve the same equation for omega0 exactly as you did and then assign the solution and assume omega1, beta and t are > 0. Then use x2(t) := subs(omega0=sqrt(beta^2+omega1^2),x(t)): Define x3(t) :=simplify(evalc(subs({omega1 = 0.9375, beta = 0.25, x0=1.0, v0=0.0},%)));

(d) Make a phase plot for this showing the decay in the position and velocity.

Hint: Define X(t) and V(t): X:= unapply(x3(t),t); V:= unapply(diff(x3(t),t),t); and plot them like last time