fig43.m

realize analysis and design for computer-controlled system

% Generates Fig 4.3
% 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=1;
clf;
%
N=5;
h=2*pi/(om*N*sqrt(1-zeta*zeta));
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,3.3]);
%title(['Fig 4.2']);
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,3.2,'(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,-2.5,1]);
%title(['N=5']);
set(gca,'Fontsize',9,'Fontname','NewCenturySchlbk');
set(gca,'xtick',[0 5 10]);
set(gca,'ytick',[-2 0]);
ylabel('Input','Fontname','NewCenturySchlbk','Fontsize',9)
text(-2.5,0.9,'(b)','fontname','NewCenturySchlbk','fontsize',9)
drawnow;
%
N=10;
h=2*pi/(om*N*sqrt(1-zeta*zeta));
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,-2.5,1]);
%title(['N=10']);
set(gca,'Fontsize',9,'Fontname','NewCenturySchlbk');
set(gca,'xtick',[0 5 10]);
set(gca,'ytick',[-2 0]);
xlabel('Time','Fontname','NewCenturySchlbk','Fontsize',9)
ylabel('Input','Fontname','NewCenturySchlbk','Fontsize',9)
text(-2.5,0.9,'(c)','fontname','NewCenturySchlbk','fontsize',9)
drawnow;
%
N=20;
h=2*pi/(om*N*sqrt(1-zeta*zeta));
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,-2.5,1]);
%title(['N=50']);
set(gca,'Fontsize',9,'Fontname','NewCenturySchlbk');
set(gca,'xtick',[0 5 10]);
set(gca,'ytick',[-2 0]);
xlabel('Time','Fontname','NewCenturySchlbk','Fontsize',9)
ylabel('Input','Fontname','NewCenturySchlbk','Fontsize',9)
text(-2.5,0.9,'(d)','fontname','NewCenturySchlbk','fontsize',9)
%
ylblalgn;
drawnow;