Hmwk_4_plots

Hmwk_4_plots - %Mike Iwanicki %AAE 340 % %Problem Set 4...

Info iconThis preview shows pages 1–2. Sign up to view the full content.

View Full Document Right Arrow Icon
%Mike Iwanicki %AAE 340 %Problem Set 4 % clear all close all clc c %Problem 1 t = linspace(0,.5,100); theta = 9.82617.*exp(-12.*t).*sin(2*sqrt(34).*t); theta_dot = -117.914.*exp(-12.*t).*sin(2*sqrt(34).*t)+114.592.*cos(2*sqrt(34).*t); dec_env = 9.82617.*exp(-12.*t); d plot(t,theta) hold on plot(0.0833,0,'rx') plot(t,dec_env,'g') plot(t,-dec_env,'g') title({'System Response';'Problem 1';'Mike Iwanicki'}) xlabel('Time (s)') ylabel('Angle (deg)') legend('Response','Damping Time Constant','Decay Envelope') grid on hold off h figure plot(theta,theta_dot) hold on plot(0,114.592,'rx') plot(2.985,28.8114,'kx') plot(1.2398,-56.5379,'mx') title({'Phase Plot';'Problem 1';'Mike Iwanicki'}) xlabel('Theta (deg)') ylabel('Theta dot (deg/s)') legend('Phase Curve','t=0','t=tau','t=2*tau') grid on hold off h %Problem 2 t = linspace(0,3*1.1808,100); phi = 0.1077.*sin(5.32.*t); phi_dot = 5.32*0.1077.*cos(5.32.*t); p figure plot(t,phi) title({'System Response';'Problem 2 w/o Dashpot';'Mike Iwanicki'}) xlabel('Time (s)')
Background image of page 1

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

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

Page1 / 2

Hmwk_4_plots - %Mike Iwanicki %AAE 340 % %Problem Set 4...

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