Ex5_15 - tf = 20.0; npts = 200; dt = (tf-t0)/npts; t =...

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

View Full Document Right Arrow Icon
% program for Prob. 5.15 % m = [1 0;0 2]; c = [4 -1;-1 2]; k = [5 -2;-2 3]; F = [1;2]; omega = 3; o for i = 1:2 for l = 1:2 Z(i,l) = -omega^2*m(i,l)+j*omega*c(i,l)+k(i,l); end end X = inv(Z)*F; X X1mag = abs(X(1,1)); X1angle = angle(X(1,1)); X2mag = abs(X(2,1)); X2angle = angle(X(2,1)); disp(' X magnitude X phase ') disp([X1mag X1angle]) disp([X2mag X2angle]) d t0 = 0;
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: tf = 20.0; npts = 200; dt = (tf-t0)/npts; t = [t0:dt:tf]'; for i = 1:npts+1 x(1,i) = X1mag*cos(omega*t(i)+X1angle); x(2,i) = X2mag*cos(omega*t(i)+X2angle); end e subplot(211) title('Rao Example 5.15'); plot(t, x(1, :)'); xlabel('t'); ylabel('x(1)'); subplot(212) plot(t, x(2, :)'); xlabel('t'); ylabel('x(2)'); y...
View Full Document

Ask a homework question - tutors are online