function tr = Trans(Sys,xo,uo,delt,tf,LLD,SBA,LON,DIM) % LTI Transient Response using 'lsim' % December 7, 2012 global tTran xTran tr = 1; if LLD >=1 % Reduced-Order Model if SBA >= 1 'Stability Axes' xo = [xo(4),xo(5),xo(6),xo(7),xo(8),xo(9),... xo(10),xo(11),xo(12),xo(13),xo(14),xo(15)]; if DIM < 6 'Fourth Order' if LON >= 1 'Longitudinal' xo = [xo(1),(xo(11) - xo(3)),xo(8),xo(3)]; uo = [uo(1),uo(4),uo(6),uo(7)]; else 'Lateral-Directional' xo = [xo(9),xo(2),xo(7),xo(10)]; uo = [uo(2),uo(3),uo(5)]; end else 'Sixth Order' if LON >= 1 'Longitudinal' xo = [xo(1),xo(3),xo(8),xo(11),xo(4),xo(6)]; uo = [uo(1),uo(4),uo(6),uo(7)]; else 'Lateral-Directional' xo = [xo(2),xo(7),xo(9),xo(10), xo(12), xo(5)]; uo = [uo(2),uo(3),uo(5)]; end end else 'Body Axes' xo = [xo(1),xo(2),xo(3),xo(7),xo(8),xo(9),... xo(10),xo(11),xo(12),xo(13),xo(14),xo(15)]; if DIM < 6 'Fourth Order' if LON >= 1 'Longitudinal' xo = [xo(1),xo(3),xo(8),xo(11)]; uo = [uo(1),uo(4),uo(6),uo(7)]; else 'Lateral-Directional' xo = [xo(2),xo(7),xo(9),xo(10)]; uo = [uo(2),uo(3),uo(5)]; end else 'Sixth Order' if LON >= 1 'Longitudinal' xo = [xo(1),xo(3),xo(8),xo(11),xo(4),xo(6)]; uo = [uo(1),uo(4),uo(6),uo(7)]; else 'Lateral-Directional' xo = [xo(2),xo(7),xo(9),xo(10), xo(12), xo(5)]; uo = [uo(2),uo(3),uo(5)]; end end end else 'Original Model' xo = [xo(1),xo(2),xo(3),xo(7),xo(8),xo(9),... xo(10),xo(11),xo(12),xo(13),xo(14),xo(15)]; end % Transient Response Simulation t = 0:delt:tf; u = []; for i = 1:length(t) u = [u;uo]; end [y,t,x] = lsim(Sys,u,t,xo); tTran = t; xTran = x; figure subplot(2,2,1) plot(t,x(:,1)) xlabel('Time, s'), ylabel('x(1)'), grid subplot(2,2,2) plot(t,x(:,2)) xlabel('Time, s'), ylabel('x(2)'), grid subplot(2,2,3) plot(t,x(:,3)) xlabel('Time, s'), ylabel('x(3)'), grid subplot(2,2,4) plot(t,x(:,4)) xlabel('Time, s'), ylabel('x(4)'), grid if (LLD <= 0) figure subplot(2,2,1) plot(t,x(:,5)) xlabel('Time, s'), ylabel('x(5)'), grid subplot(2,2,2) plot(t,x(:,6)) xlabel('Time, s'), ylabel('x(6)'), grid subplot(2,2,3) plot(t,x(:,7)) xlabel('Time, s'), ylabel('x(7)'), grid subplot(2,2,4) plot(t,x(:,8)) xlabel('Time, s'), ylabel('x(8)'), grid figure subplot(2,2,1) plot(t,x(:,9)) xlabel('Time, s'), ylabel('x(9)'), grid subplot(2,2,2) plot(t,x(:,10)) xlabel('Time, s'), ylabel('x(10)'), grid subplot(2,2,3) plot(t,x(:,11)) xlabel('Time, s'), ylabel('x(11)'), grid subplot(2,2,4) plot(t,x(:,12)) xlabel('Time, s'), ylabel('x(12)'), grid end if (LLD >= 1 & DIM >= 6) figure subplot(1,2,1) plot(t,x(:,5)) xlabel('Time, s'), ylabel('x(5)'), grid subplot(1,2,2) plot(t,x(:,6)) xlabel('Time, s'), ylabel('x(6)'), grid end