ASE 370L - HW 4 Codes and Plots

ASE 370L - HW 4 Codes and Plots - Problem E 4.8 Plots and...

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

View Full Document Right Arrow Icon
Problem E 4.8 Plots and Code % ------------------------------------------------ % Flight Controls HW#4 Problem E 4.8 Solution % ------------------------------------------------ clear all close all clc %Inputs k=40:10:400; tmax=1; tstep=.001; for i=1:length(k) numopen = 46.24*k(i)*[1 50]; denopen = [1 (200+16.7) ((200*16.7)+72.9) (200*72.9)]; G=tf(numopen,denopen);
Background image of page 1

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

View Full DocumentRight Arrow Icon
numh=[425]; denh=[1 425]; H=tf(numh,denh); sys=feedback(G,H); sssys=ss(sys); t=[0:tstep:tmax]; [y,x]=step(sssys,t); ymax=max(y); eSS(i)=6.306/(6.306+k(i)); OS(i)=abs(ymax-y(end)); OSp(i)=abs(OS(i)/y(end))*100; end figure hold on plot(k,OSp, 'LineWidth' ,2) grid on xlabel( 'Gain "K"' ) ylabel( 'Percentage Overshoot' ) hold off figure hold on plot(k,OS, 'LineWidth' ,2) grid on xlabel( 'Gain "K"' ) ylabel( 'Overshoot' ) hold off figure hold on plot(k,eSS, 'LineWidth' ,2) grid on xlabel( 'Gain "K"' ) ylabel( 'Steady State Error' ) hold off
Background image of page 2
Problem E 4.10 Plots and Code % ------------------------------------------------ % Flight Controls HW#4 Problem E 4.10 Solution
Background image of page 3

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

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

Page1 / 8

ASE 370L - HW 4 Codes and Plots - Problem E 4.8 Plots and...

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

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