vibration61

vibration61 - x4...

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

View Full Document Right Arrow Icon
%% Problem 6.1 clc;clear m = 1; %kg w = 1.4; %rad/s wn = 1.87; %rad/s z = 0.2/(2*wn); wd = wn*sqrt(1-z^2); theta = acos(z); syms s t k = s^2+0.2*s+3.5; Xs = ((1.5*w)/(s^2+w^2))*(1/k)+(4*exp(-3*s))/k+(s*.1)/k+(.1*.2)/k; X x1 =0.1*(-wn/wd)*exp(-z*wn*t)*sin(wd*t-theta); x2 =0.1714*(-wn/wd)*exp(-z*wn*t)*sin(wd*t-theta); x3 =0.02*(1/wd)*exp(-z*wn*t)*sin(wd*t);
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: x4 =4*(1/wd)*exp(-z*wn*(t-3))*sin(wd*(t-3))*heaviside(t-3); x5 =-0.1714*cos(1.4*t); x6 =-1.2857*(1/wd)*exp(-z*wn*t)*sin(wd*t); x7 =(1.32/w)*sin(w*t); x Xt = x1 + x2 + x3 + x4 + x5 + x6 + x7; Xt2 = ilaplace(Xs) figure(1) ezplot(Xt, [0 70]) figure(2) ezplot(Xt2, [0 70])...
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