fig912.m

来自「realize analysis and design for computer」· M 代码 · 共 48 行

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

lw=0.9;
N=500;
tstep=5/N;
 
for k=1:N  
al=k*tstep;
t(k)=al;

r1=sqrt(max(0,1-1/(4*al*al)));   
r2=sqrt(max(0,1-9/(4*al*al))) ;  
r3=sqrt(max(0,1-25/(4*al*al)));  
r4=sqrt(max(0,1-49/(4*al*al)));  
r5=sqrt(max(0,1-81/(4*al*al))); 
 
if al<0.5 y=0; 
elseif al<1.5 y=r1;
elseif al<2.5 y=r1+r2;
elseif al<3.5 y=r1+r2+r3;
elseif al<4.5 y=r1+r2+r3+r4;
elseif al<5.5 y=r1+r2+r3+r4+r5; else y=0;
end; 
 
yc(k)=4*y/(pi*al);
end;
 
subplot(2,1,1)
plot(t,yc,'k','Linew',lw)
axis([0 5 0 1.4])
set(gca,'FontSize',9,'Fontname','NewCenturySchlbk')
set(gca,'xtick',[0 2 4]);
set(gca,'ytick',[0 1]);
xlabel('\it{A}/\delta','Fontname','NewCenturySchlbk','Fontsize',9)
ylabel('\it{Y_c(A)}','Fontname','NewCenturySchlbk','Fontsize',9)
%title(['Fig 9.12']);
drawnow;