fig1210.m

realize analysis and design for computer-controlled system

% Generates Fig 12.10 (outputvar)
% Computer Controlled Systems (3rd ed)
% Author: B. Wittenmark
% Last edit: 1997-07-10
% Copyright (c) 1996 by K. J. 舠tr鰉, B. Wittenmark and
% Department of Automatic Control, Lund Institute of
% Technology, Lund, Sweden

newplot;
set(gcf,'PaperUnits','centimeters','PaperPosition',[2 2 13 8])
set(gcf,'Units','centimeters','position',get(gcf,'PaperPosition'))

lw=0.9;

hv=[0.2 0.5 1 2];
ve=1;
vv=1;
py=zeros(4,201);
t=zeros(1,201);
for j=1:4,
  h=hv(j);
  for i=1:201,
    t(i)=(i-1)*0.01;
    s1=h*vv/ve;
    c=-1-s1/2+sqrt(s1+s1*s1/4);
    s2=-ve/c;
    s=rem(t(i),h);
    p1=(1-s/h*(1+c))*(1-s/h*(1+c))*(s2-ve);
    p2=(s/h*(1+c))*(s/h*(1+c))*ve+vv*s+ve;
    py(j,i)=p1+p2;
  end;
end;
clf;
plot(t,py(1,:),'k-','Linew',lw);
hold on;
plot(t,py(2,:),'k--','Linew',lw);
plot(t,py(3,:),'k-.','Linew',lw);
plot(t,py(4,:),'k:','Linew',lw);
axis([0,2,0,4]);
%title(['Output variance, h=0.2, 0.5, 1, 2']);
set(gca,'Fontsize',9,'Fontname','NewCenturySchlbk');
set(gca,'xtick',[0 1 2]);
set(gca,'ytick',[0 2 4]);
ylabel('Output variance','Fontname','NewCenturySchlbk','Fontsize',9)
xlabel('Time','Fontname','NewCenturySchlbk','Fontsize',9)