fig521.m

realize analysis and design for computer-controlled system

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

hnoise=0.01;
 
block521;

lw=0.9;
td=15;
tn=30;
l0=-0.5;
h=0.2;
% Calculation of Besselfilter polynomials
fc=1;
om=2*pi*fc;
z1=0.49;
om1=1.90*om;
z2=0.82;
om2=1.69*om;
z3=0.98;
om3=1.61*om;
bbes=om1*om1*om2*om2*om3*om3;
abes=conv([1 2*z1*om1 om1*om1],[1 2*z2*om2 om2*om2]);
abes=conv(abes,[1 2*z3*om3 om3*om3]);
%
r1=-1.735822;
r2=1.5876159;
r3=-0.8240501;
r4=-2.74110716e-03;
r5=-2.50026747e-02;
s0=8.7195875;
s1=-16.316408;
s2=7.7985257;
s3=0;
s4=0;
s5=0;
t0=1.8296419;
t1=-3.8639301;
t2=3.6366312;
t3=-1.8478182;
t4=0.50703278;
t5=-5.98520427e-02; 
%Areg=[-r1 1 0 0 0;-r2 0 1 0 0;-r3 0 0 1 0;-r4 0 0 0 1;-r5 0 0 0 0]
%
clf;
%
[t,x,y]=sim('block521',[0, 40]);
subplot(3,2,1);
plot(t,y(:,3),'k--',t,y(:,1),'k','Linew',lw);
axis([0,40,0,1.3]);
set(gca,'FontSize',9,'Fontname','NewCenturySchlbk')
set(gca,'xtick',[0 20 40]);
set(gca,'ytick',[0 1]);
ylabel('Output','Fontname','NewCenturySchlbk','Fontsize',9)
%title('Fig 5.21','fontname','NewCenturySchlbk','FontSize',9);
xlabel('Time','Fontname','NewCenturySchlbk','Fontsize',9)
drawnow;
%
subplot(3,2,2);
[ts,ys]=stairs(t,y(:,2));
plot(ts,ys,'m','Linew',lw);
axis([0,40,-0.5,3.5]);
set(gca,'FontSize',9,'Fontname','NewCenturySchlbk')
set(gca,'xtick',[0 20 40]);
set(gca,'ytick',[0 2]);
ylabel('Input','Fontname','NewCenturySchlbk','Fontsize',9)
xlabel('Time','Fontname','NewCenturySchlbk','Fontsize',9)
drawnow;