vibration61

vibration61 - Xt2 = ilaplace(Xs) figure(1) ezplot(Xt, [0...

Info iconThis preview shows pages 1–2. 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; 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); 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); Xt = x1 + x2 + x3 + x4 + x5 + x6 + x7;
Background image of page 1

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
Background image of page 2
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: Xt2 = ilaplace(Xs) figure(1) ezplot(Xt, [0 70]) figure(2) ezplot(Xt2, [0 70]) Xt2 = (33*sin((7*t)/5))/35 - (6*cos((7*t)/5))/35 + (cos((349^(1/2)*t)/10) - (349^(1/2)*sin((349^(1/2)*t)/10))/349)/(10*exp(t/10)) + (6*(cos((349^(1/2)*t)/10) - (76*349^(1/2)*sin((349^(1/2)*t)/10))/349))/(35*exp(t/10)) + (349^(1/2)*sin((349^(1/2)*t)/10))/(1745*exp(t/10)) + (40*349^(1/2)*sin((349^(1/2)*(t - 3))/10)*heaviside(t - 3)*exp(3/10 - t/10))/349 Published with MATLAB 7.8...
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.

Page1 / 2

vibration61 - Xt2 = ilaplace(Xs) figure(1) ezplot(Xt, [0...

This preview shows document pages 1 - 2. Sign up to view the full document.

View Full Document Right Arrow Icon
Ask a homework question - tutors are online