Matlab help
jeddawe
archive.zipspringmass(1)-3.m
function springmass(p,omega0,y0,v0,tf) %springmass(p,omega0,y0,v0,tf,omega) solves the damped harmonic problems % y'' + 2 p y' + omega0^2 y = 0 % with y(t=0)=y0, y'(t=0)=v0 from t=0 to tf using fourth-order Runge-Kutta t0 = 0; tspan = [t0 tf]; w0 = [y0; v0]; opts = odeset('RelTol',10^-9,'AbsTol',10^-12); [t,w] = ode45(@f,tspan,w0,opts); figure plot(t,w(:,1),'b-') xlabel('t','FontSize',16) ylabel('y','FontSize',16) xlim([t0 tf]) function wprime = f(t,w) wprime = [w(2); -omega0^2*w(1)-2*p*w(2)]; end end
__MACOSX/._springmass(1)-3.m
springmassdriven(1).m
function springmassdriven(p,omega0,y0,v0,tf,omega) %springmassdriven(p,omega0,y0,v0,tf,omega) solves the forced damped %harmonic problems % y'' + 2 p y' + omega0^2 y = cos(omega t) % with y(t=0)=y0, y'(t=0)=v0 from t=0 to tf using fourth-order Runge-Kutta t0 = 0; tspan = [t0 tf]; w0 = [y0; v0]; opts = odeset('RelTol',10^-9,'AbsTol',10^-12); [t,w] = ode45(@f,tspan,w0,opts); figure plot(t,w(:,1),'b-') xlabel('t','FontSize',16) ylabel('y','FontSize',16) xlim([t0 tf]) function wprime = f(t,w) wprime = [w(2); -omega0^2*w(1)-2*p*w(2)+cos(omega*t)]; end end
