fig514.m

realize analysis and design for computer-controlled system

% Generates Fig 5.14
% Computer Controlled Systems (3rd ed)
% Author: B. Wittenmark
% Last edit: 1997-07-02
% 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'))
 
block510;

lw=0.9;
K=1;
a=2;
z=0.707;
an=0.01;
wn=2;
td=50;
tn=100;
l0=-0.05;
x10=0;
x20=0;
h=1;
w=0.2;
%
wh=w*h;
z1=sqrt(1-z*z);
a1m=-2*exp(-z*wh)*cos(z1*wh);
a2m=exp(-2*z*wh);
ams=1+a1m+a2m;

ae=exp(-a*h);
a1o=-2*ae;
a2o=ae*ae;
aos=1+a1o+a2o;
amp=2+a1m;
aop=2+a1o;

p1=ams*aos;
pp=amp*aos+ams*aop;
p2p=2*(ams+amp*aop+aos);
h2=h*h;
s0=(p1/4-pp/2+p2p/2)/h2;
s1=(-p1+2*pp-p2p)/h2;
s2p=(7*p1/4-3*pp/2+p2p/2);
s2=s2p/h2;
r=s2p/2-a2m*a2o;
r1=r-1;
r2=-r;
t0=ams/h2;
t1=t0*a1o;
t2=t0*a2o;
%
clf;
K=0.4;
%
[t,x,y]=sim('block510',[0, 100]);
subplot(3,2,1);
plot(t,y(:,3),'k--',t,y(:,1),'k','Linew',lw);
axis([0,100,0,1.3]);
set(gca,'FontSize',9,'Fontname','NewCenturySchlbk')
set(gca,'xtick',[0 50 100]);
set(gca,'ytick',[0 1]);
ylabel('Output','Fontname','NewCenturySchlbk','Fontsize',9)
text(-20,1.3,'(a)','Fontname','NewCenturySchlbk','Fontsize',9)
%title('Fig 5.14','fontname','NewCenturySchlbk','FontSize',9);
drawnow;
%
K=0.8;
%
[t,x,y]=sim('block510',[0, 100]);
subplot(3,2,2);
plot(t,y(:,3),'k--',t,y(:,1),'k','Linew',lw);
axis([0,100,0,1.3]);
set(gca,'FontSize',9,'Fontname','NewCenturySchlbk')
set(gca,'xtick',[0 50 100]);
set(gca,'ytick',[0 1]);
ylabel('Output','Fontname','NewCenturySchlbk','Fontsize',9)
text(-20,1.3,'(b)','Fontname','NewCenturySchlbk','Fontsize',9)
drawnow;
%
K=1.2;
%
[t,x,y]=sim('block510',[0, 100]);
subplot(3,2,3);
plot(t,y(:,3),'k--',t,y(:,1),'k','Linew',lw);
axis([0,100,0,1.3]);
set(gca,'FontSize',9,'Fontname','NewCenturySchlbk')
set(gca,'xtick',[0 50 100]);
set(gca,'ytick',[0 1]);
ylabel('Output','Fontname','NewCenturySchlbk','Fontsize',9)
text(-20,1.3,'(c)','Fontname','NewCenturySchlbk','Fontsize',9)
xlabel('Time','Fontname','NewCenturySchlbk','Fontsize',9)
drawnow;
%
K=2.5;
%
[t,x,y]=sim('block510',[0, 100]);
subplot(3,2,4);
plot(t,y(:,3),'k--',t,y(:,1),'k','Linew',lw);
axis([0,100,0,1.3]);
set(gca,'FontSize',9,'Fontname','NewCenturySchlbk')
set(gca,'xtick',[0 50 100]);
set(gca,'ytick',[0 1]);
ylabel('Output','Fontname','NewCenturySchlbk','Fontsize',9)
text(-20,1.3,'(d)','Fontname','NewCenturySchlbk','Fontsize',9)
xlabel('Time','Fontname','NewCenturySchlbk','Fontsize',9)
drawnow;