fig45.m

realize analysis and design for computer-controlled system

% Generates Fig 4.5
% Computer Controlled Systems (3rd ed)
% Author: B. Wittenmark
% Last edit: 1997-06-30
% 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 10])
set(gcf,'Units','centimeters','position',get(gcf,'PaperPosition'))

lw=0.9;
t0=0;
m=0;
x10=1;
x20=1;
zeta=0.707;
om=1e3;
clf;
%
h=2;
p1=-2*exp(-zeta*om*h)*cos(om*h*sqrt(1-zeta*zeta));
p2=exp(-2*zeta*om*h);
l1=(1+p1+p2)/(h*h);
l2=(3+p1-p2)/(2*h);
[t,x,y]=sim('block42',[0, 10]);
subplot(3,2,1);
plot(t,y(:,1),'k-.','Linew',lw);
axis([0,10,-0.5,2.5]);
%title(['Fig 4.5']);
set(gca,'Fontsize',9,'Fontname','NewCenturySchlbk');
set(gca,'xtick',[0 5 10]);
set(gca,'ytick',[0 2]);
ylabel('Output','Fontname','NewCenturySchlbk','Fontsize',9)
text(-2.5,2.4,'(a)','fontname','NewCenturySchlbk','fontsize',9)
hold on
drawnow;
%
subplot(3,2,2);
[ts,ys]=stairs(t,y(:,3));
plot(ts,ys,'m-.','Linew',lw);
axis([0,10,-8,5.5]);
%title(['h=2']);
set(gca,'Fontsize',9,'Fontname','NewCenturySchlbk');
set(gca,'xtick',[0 5 10]);
set(gca,'ytick',[-5 0 5]);
ylabel('Input','Fontname','NewCenturySchlbk','Fontsize',9)
text(-2.5,5,'(b)','fontname','NewCenturySchlbk','fontsize',9)
drawnow;
%
h=1;
p1=-2*exp(-zeta*om*h)*cos(om*h*sqrt(1-zeta*zeta));
p2=exp(-2*zeta*om*h);
l1=(1+p1+p2)/(h*h);
l2=(3+p1-p2)/(2*h);
[t,x,y]=sim('block42',[0, 10]);
subplot(3,2,1);
plot(t,y(:,1),'k--','Linew',lw);
drawnow;
%
subplot(3,2,3);
[ts,ys]=stairs(t,y(:,3));
plot(ts,ys,'m--','Linew',lw);
axis([0,10,-8,5.5]);
%title(['h=1']);
set(gca,'Fontsize',9,'Fontname','NewCenturySchlbk');
set(gca,'xtick',[0 5 10]);
set(gca,'ytick',[-5 0 5]);
xlabel('Time','Fontname','NewCenturySchlbk','Fontsize',9)
ylabel('Input','Fontname','NewCenturySchlbk','Fontsize',9)
text(-2.5,5,'(c)','fontname','NewCenturySchlbk','fontsize',9)
drawnow;
%
h=0.5;
p1=-2*exp(-zeta*om*h)*cos(om*h*sqrt(1-zeta*zeta));
p2=exp(-2*zeta*om*h);
l1=(1+p1+p2)/(h*h);
l2=(3+p1-p2)/(2*h);
[t,x,y]=sim('block42',[0, 10]);
subplot(3,2,1);
plot(t,y(:,1),'k-','Linew',lw);
drawnow;
hold off
%
subplot(3,2,4);
[ts,ys]=stairs(t,y(:,3));
plot(ts,ys,'m-','Linew',lw);
axis([0,10,-8,5.5]);
%title(['h=0.5']);
set(gca,'Fontsize',9,'Fontname','NewCenturySchlbk');
set(gca,'xtick',[0 5 10]);
set(gca,'ytick',[-5 0 5]);
xlabel('Time','Fontname','NewCenturySchlbk','Fontsize',9)
ylabel('Input','Fontname','NewCenturySchlbk','Fontsize',9)
text(-2.5,5,'(d)','fontname','NewCenturySchlbk','fontsize',9)
drawnow;
ylblalgn;