fig114.m

realize analysis and design for computer-controlled system

% Generates Fig 11.4 (firstsl)
% 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;clf;
clear;
set(gcf,'PaperUnits','centimeters','PaperPosition',[2 2 13 9.5])
set(gcf,'Units','centimeters','position',get(gcf,'PaperPosition'))

Ms=9;
lw=0.9;
q1=1;
q2=10;
q0=10;
s=zeros(1,6);
t=zeros(1,6);
s(6)=q0;
t(6)=5;

for kk=1:5,
t(6-kk)=5-kk;
s(6-kk)=s(7-kk)+q1-s(7-kk)*s(7-kk)/(q2+s(7-kk));
l(6-kk)=s(7-kk)/(s(7-kk)+q2);
end;

x(1)=1;
for kk=1:5,
u1(kk)=-l(kk)*x(kk);
x(kk+1)=(1-l(kk))*x(kk); end;

subplot(3,1,1)
hpl=plot(t,s,'m.');
axis([0 5 0 10])
hold on
plot(t,s,'k--','Linew',lw)
set(gca,'Fontsize',9,'Fontname','NewCenturySchlbk');
set(gca,'ytick',[0 10]);
set(gca,'xtick',[0 5]);
set(hpl,'Markersize',Ms);

subplot(3,1,2)
plot(t(1:5),l,'k--','Linew',lw)
hold on
hpl=plot(t(1:5),l,'m.');
axis([0 5 0 0.5])
set(gca,'Fontsize',9,'Fontname','NewCenturySchlbk');
set(gca,'ytick',[0 0.5]);
set(gca,'xtick',[0 5]);
set(hpl,'Markersize',Ms);

subplot(3,1,3)
plot(t,x,'k--','Linew',lw)
hold on
hpl=plot(t,x,'m.');
axis([0 5 0 1])
set(gca,'Fontsize',9,'Fontname','NewCenturySchlbk');
set(gca,'ytick',[0 1]);
set(gca,'xtick',[0 5]);
set(hpl,'Markersize',Ms);


q0=0;
s(6)=q0;
t(6)=5;
for kk=1:5,
t(6-kk)=5-kk;
s(6-kk)=s(7-kk)+q1-s(7-kk)*s(7-kk)/(q2+s(7-kk));
l(6-kk)=s(7-kk)/(s(7-kk)+q2);
end;

x(1)=1;
for kk=1:5,
u2(kk)=-l(kk)*x(kk);
x(kk+1)=(1-l(kk))*x(kk); end;

subplot(3,1,1)
hpl=plot(t,s,'m.');
set(hpl,'Markersize',Ms);
plot(t,s,'b-.','Linew',lw)

subplot(3,1,2)
plot(t(1:5),l,'b-.','Linew',lw)
hpl=plot(t(1:5),l,'m.');
set(hpl,'Markersize',Ms);

subplot(3,1,3)
plot(t,x,'b-.','Linew',lw)
hpl=plot(t,x,'m.');
set(hpl,'Markersize',Ms);

q0=max(roots([1 -q1 -q1*q2]));
s(6)=q0;
t(6)=5;
for kk=1:5,
t(6-kk)=5-kk;
s(6-kk)=s(7-kk)+q1-s(7-kk)*s(7-kk)/(q2+s(7-kk));
l(6-kk)=s(7-kk)/(s(7-kk)+q2);
end;

x(1)=1;
for kk=1:5,
u3(kk)=-l(kk)*x(kk);
x(kk+1)=(1-l(kk))*x(kk); end;

subplot(3,1,1)
hpl=plot(t,s,'m.');
set(hpl,'Markersize',Ms);
plot(t,s,'g-','Linew',lw)
ylabel('s(k)','Fontname','NewCenturySchlbk','Fontang','italic','Fontsize',9)
subplot(3,1,2)
plot(t(1:5),l,'g-','Linew',lw)
hpl=plot(t(1:5),l,'m.');
set(hpl,'Markersize',Ms);
ylabel('l(k)','Fontname','NewCenturySchlbk','Fontang','italic','Fontsize',9)

subplot(3,1,3)
plot(t,x,'g-','Linew',lw)
hpl=plot(t,x,'m.');
set(hpl,'Markersize',Ms);
xlabel('Time','Fontname','NewCenturySchlbk','Fontsize',9)
ylabel('x(k)','Fontname','NewCenturySchlbk','Fontang','italic','Fontsize',9)
drawnow;
%
ylblalgn;