clear
ad033init
sim('ad033s');

t = data.time;
u = data.signals.values(:,1);
y = data.signals.values(:,2);
w = data.signals.values(:,3);
[tt,ww]=stairs(t,w);
plot(t,u,'k--', t,y, 'k-',tt,ww,'k:')
title('(a)')
legend('u','y','w')
xlabel('t')
ylabel('u,y,w')
set(1, 'PaperUnits', 'centimeters')
set(1, 'PaperPosition', [3 10 11 8.25])
print -deps 'oad034a.eps'

t = param.time;
b1 = param.signals.values(:,1);
b0 = param.signals.values(:,2);
a2 = param.signals.values(:,3);
a1 = param.signals.values(:,4);
[tt,bb1]=stairs(t,b1);
[tt,bb0]=stairs(t,b0);
[tt,aa1]=stairs(t,a1);
[tt,aa2]=stairs(t,a2);
plot(tt,bb1,'k-',tt,bb0,'k-.',tt,aa1,'k--',tt,aa2,'k:')
title('(b)')
legend('b1','b0','a1','a2')
xlabel('t')
ylabel('Param')
set(1, 'PaperUnits', 'centimeters')
set(1, 'PaperPosition', [3 10 11 8.25])
print -deps 'oad034b.eps'

