vibration611

vibration611 - Xt=(1/(wn^2*m))*(x1-x2); Xt1 = ilaplace(Xs);...

Info iconThis preview shows page 1. Sign up to view the full content.

View Full Document Right Arrow Icon
%% Problem 6.11 clc; clear f = 1; t1 = 1; z = .1; m = 1; wn = 4; wd=wn*sqrt(1-z^2); theta= acos(z); syms s t Xs = (1/m)*((1-exp(-t1*s))/s)*(1/(s^2+2*z*wn*s+wn^2)); X x1 = 1-(wn/wd)*exp(-z*wn*t)*sin(wd*t+theta); x2 = (1-(wn/wd)*exp(-z*wn*(t-t1))*sin(wd*(t-t1)+theta))*heaviside(t-t1);
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: Xt=(1/(wn^2*m))*(x1-x2); Xt1 = ilaplace(Xs); X figure(1) ezplot(Xt, [0 20]) axis([0 20 -.1 .12]) figure(2) ezplot(Xt1, [0 20]) axis([0 20 -.1 .12])...
View Full Document

This note was uploaded on 02/13/2011 for the course ENME 361 taught by Professor Yoo during the Spring '11 term at Maryland.

Ask a homework question - tutors are online