fig123.m

realize analysis and design for computer-controlled system

% Generates Fig 12.3 (pred)
% 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
set(gcf,'PaperUnits','centimeters','PaperPosition',[2 2 13 13.5])
set(gcf,'Units','centimeters','position',get(gcf,'PaperPosition'))
 
lw=0.9;
N=100;

Apol=[1 -1.5 0.7];
Cpol=[1 -0.2 0.5];

%m=1
Gpol=[1.3 -0.2];
[t,x,y]=sim('blockpred',[0, 100]);
subplot(4,1,1)
plot(t,y(:,1),'k--',t,y(:,2),'m-','Linew',lw);
axis([0,N,-10  10]);
set(gca,'Fontsize',9,'Fontname','NewCenturySchlbk');
set(gca,'xtick',[0 50 100],'ytick',[-10 0 10]);
ylabel('{\it y} and {\it y}','Fontsize',9,'Fontname','NewCenturySchlbk')
text(-13,10.5,'(a)','Fontsize',9,'Fontname','NewCenturySchlbk');

subplot(4,1,4)
err=y(:,1)-y(:,2);
loss=cumsum(err.*err);
plot(t,loss,'k-.','Linew',lw);
axis([0,N,0 600]);
set(gca,'Fontsize',9,'Fontname','NewCenturySchlbk');
set(gca,'xtick',[0 50 100],'ytick',[0 500]);
ylabel('Loss','Fontsize',9,'Fontname','NewCenturySchlbk')
xlabel('Time','Fontsize',9,'Fontname','NewCenturySchlbk')
text(-13,600,'(d)','Fontsize',9,'Fontname','NewCenturySchlbk');
hold on

%m=3
Gpol=[1.715 -1.225];
[t,x,y]=sim('blockpred',[0, N]);
subplot(4,1,2)
plot(t,y(:,1),'k--',t,y(:,4),'m-','Linew',lw);
axis([0,N,-10  10]);
set(gca,'Fontsize',9,'Fontname','NewCenturySchlbk');
set(gca,'xtick',[0 50 100],'ytick',[-10 0 10]);
ylabel('{\it y} and {\it y}','Fontsize',9,'Fontname','NewCenturySchlbk')
text(-13,10.5,'(b)','Fontsize',9,'Fontname','NewCenturySchlbk');

subplot(4,1,4)
err=y(:,1)-y(:,4);
loss=cumsum(err.*err);
plot(t,loss,'k-','Linew',lw);

%m=5
Gpol=[0.8208 -0.9433];
[t,x,y]=sim('blockpred',[0, N]);
subplot(4,1,3)
plot(t,y(:,1),'k--',t,y(:,6),'m-','Linew',lw);
axis([0,N,-10  10]);
set(gca,'Fontsize',9,'Fontname','NewCenturySchlbk');
set(gca,'xtick',[0 50 100],'ytick',[-10 0 10]);
ylabel('{\it y} and {\it y}','Fontsize',9,'Fontname','NewCenturySchlbk')

text(-13,10.5,'(c)','Fontsize',9,'Fontname','NewCenturySchlbk');

subplot(4,1,4)
err=y(:,1)-y(:,6);
loss=cumsum(err.*err);
plot(t,loss,'k:','Linew',lw);

%m=2
Gpol=[1.75 -0.91];
[t,x,y]=sim('blockpred',[0, N]);
subplot(4,1,4)
err=y(:,1)-y(:,3);
loss=cumsum(err.*err);
plot(t,loss,'k--','Linew',lw);

ylblalgn
subplot(4,1,1)
a=hex2dec('d9');
text(-10.6,4.1,setstr(a),'fontname','symbol','fontsize',7,'rotation',90)
subplot(4,1,2)
text(-10.6,4.1,setstr(a),'fontname','symbol','fontsize',7,'rotation',90)
subplot(4,1,3)
text(-10.6,4.1,setstr(a),'fontname','symbol','fontsize',7,'rotation',90)